Impurity Combination Effect on Oxygen Absorption in α 2 -Ti 3 Al

: The effect of substitutional impurities of the transition metals of VB–VIIB groups on the oxygen absorption in the doped α 2 -Ti 3 Al alloy was studied by the projector-augmented wave method within the density functional theory. It is established that all considered impurities prefer to substitute for a Ti atom rather than an Al atom. Changes in the absorption energy due to impurities being in the ﬁrst neighbors of the oxygen atom were estimated. It was demonstrated that the doping resulted in a decrease in the oxygen absorption energy, which is mainly caused by the chemical contribution to it. The interaction energy between impurity atoms was calculated in the dependence on the interatomic distance. It was shown that the conﬁguration with the impurity atoms being in the second neighbors of each other was stable in comparison with other possible conﬁgurations. The inﬂuence of two impurity atoms being in the ﬁrst neighbors of oxygen is additively enhanced. It was revealed that the effect of two impurity atoms on the oxygen absorption energy can be estimated as the sum of the effects of separate impurities with an accuracy of more than ~90%. atoms substitute for Ti with probability of ~100%. In this connection, further discussion is given only in case of impurities on the Ti sublattice. impurity and the effect becomes stronger with the ﬁlling of its d -shell; (iii) if oxygen locates in the region between impurity atoms, their effect is additively enhanced, which makes it possible to predict the joint effect of impurities with an accuracy of ~91%; (iv) regions near impurities can prevent the oxygen diffusion, since its penetration is energetically unfavorable. The latter is expected to be studied in more detail in our forthcoming papers.


Introduction
Intermetallic alloys based on titanium and aluminum are known to be characterized by their low densities, high specific strength and stiffness, and high creep resistance under moderately elevated temperatures [1][2][3][4][5], which makes them ideal materials for various high-temperature structural applications in automobile, aerospace, and gas turbine industries [6][7][8]. Two of the most promising intermetallic phases of the binary system Ti-Al for technical applications are α 2 -Ti 3 Al and γ-TiAl. These materials are considered to be excellent for manufacturing static components or rotating parts [3,4,6,9]. At the same time, the Ti 3 Al alloy has a specific modulus and stress rupture resistance comparable to that of the superalloys; however, the complete absence of room temperature plasticity posed the primary challenge in using it as a structural material. Moreover, the single γ-phase TiAl alloy is brittle with practically no deformability at temperatures up to~1000 K [1,6]. At high temperature, mixed oxide scales grow by the competitive oxidation of Ti and Al, which prevents the formation of a continuous and dense α-alumina [10][11][12][13][14].
Two-phase TiAl/Ti 3 Al lamellar alloys possess a number of properties that distinguish them from the TiAl and Ti 3 Al constituent single-phase alloys. For example, toughness, fatigue and creep resistance, and other properties are appreciably improved [15][16][17][18]. It is demonstrated that the minimum creep rate of two-phase lamellar alloys is up to an order of magnitude lower than the creep rates of the constituent single-phase alloys [19,20]. In fully lamellar structures, a very low dislocation activity and almost no twinning structures are observed [21]. At the same time, the high-temperature corrosion resistance of the two-phase lamellar alloys is still insufficient.
The addition of some alloying elements can enhance both the oxidation resistance and mechanical properties of single-phase TiAl and Ti 3 Al and two-phase lamellar TiAl/Ti 3 Al alloys. For example, an addition of V, Cr, Nb, Mo, Ta, and W produces solid solution strengthening (see [1,6,[22][23][24][25][26] and references therein). The latter is also advantageous for the creep resistance. The Nb, Mo, Ta, and W elements are very effective in improving oxidation resistance [1,6,[27][28][29][30][31], while the addition of Cr or V does not improve it [29]. It is a reason why almost all multicomponent alloys contain at least 2 at.% Nb, which is found to be the most appropriate.
In several papers [32][33][34][35], the energetics of substitutional defect formation in the α 2 -Ti 3 Al and γ-TiAl alloys was investigated by experimental and theoretical methods. Summing up the obtained results, it is possible to conclude that the energy-preferable situation is if V, Cr, Nb, Mo, Ta, W, and Re atoms occupy the Ti sublattice in the Ti 3 Al alloy. In our earlier papers [36][37][38], the oxygen diffusion properties in both γ-TiAl and α 2 -Ti 3 Al alloys were investigated. It was established that the highest absorption energy corresponded to the positions with the highest number of Ti atoms in the nearest neighbors. The key migration barrier is connected with a jump from these sites to the nearest tetrahedral sites. In the paper [38], the impurity influence on the temperature-dependent diffusion coefficient in the Ti 3 Al alloy was investigated. It was concluded that the dominant factor determining the change in the diffusion coefficient upon doping is a change in the oxygen absorption energy. Note that in [38], only the influence of single impurity atom on the absorption and migration energies was considered. At the same time, real alloys are multicomponent ones and the influence of several dopants on the absorption energy is desirable to be known.
Thus, the main goal of the present paper is to study the influence of an increase of impurity concentration or the addition of other impurities on the oxygen absorption in the α 2 -Ti 3 Al alloy. The increase in the concentration of impurity atoms leads to oxygen interaction with two dopant atoms simultaneously and makes the absorption process more complicated. The transition metals of VB-VIIB groups, such as V, Cr, Nb, Mo, Ta, W, and Re, are considered because of their favorable effect on the corrosion resistance or the mechanical properties of the alloy.

Computational Details
First-principles calculations of the alloy electronic structure were performed by the plane-wave basis projector-augmented wave (PAW) method [39,40] implemented by VASP (The Vienna Ab initio Simulation Package) code [41][42][43]. The generalized gradient approximation for an exchange-correlation functional in the PBE form [44] was used. It should be noted that the α 2 -Ti 3 Al alloy has a hexagonal close-packed D0 19 structure and is characterized by space group 194 (P6 3 /mmc), where aluminum and titanium atoms occupy 2c (1/3,2/3,1/4) and 6h (x,2x,1/4) sites with x = −0.170, according to the Wyckoff classification. The theoretical lattice parameters of the alloy are a = 5.736 Å and c = 4.639 Å. Based on convergence tests, the plane-wave cutoff energy was set at 550 eV and a Γ-centered 13 × 13 × 17 k-point mesh was used in the primitive unit cell of the α 2 -Ti 3 Al alloy. The (2 × 2 × 2) supercell (i.e., 64 atoms) with mesh of 5 × 5 × 7 k-points was used to study the substitutional defect formation and the oxygen absorption. The full structural optimization scheme involving change in the atomic positions as well as in the cell shape and volume was applied. Convergence was considered to be achieved if the difference in total energies for the next two iterations did not exceed 10 -5 eV. Atomic positions were relaxed using the conjugate gradient algorithm until the forces at the atoms were at most 10 -3 eV/Å.
The oxygen absorption energy was calculated by the formula: where E(Ti 3 Al + O) and E(Ti 3 Al) are the total energies of the alloy with and without oxygen, respectively, and E(O 2 ) is the total energy of the oxygen molecule calculated in an empty cell with sides of 12 × 13 × 14 Å. Factor 1/2 corresponds to the absorption of one oxygen atom in the computational cell. The formation energy of a substitutional defect was estimated by following equation: where E(Ti 3 Al + Imp) and E(Ti 3 Al) are the total energies of the alloy with and without impurity, respectively, and E(Ti/Al) and E(Imp) is the total energy of the titanium/aluminum (in dependence on occupied sublattice) and the impurity metal in the ground state per atom, respectively. In accordance with the approach suggested by Lozovoi with co-authors [45] for embrittlement of grain boundaries, the effect of a substitutional (interstitial) defect on any physical quantity can be decomposed into three (two) contributions. In Figure 1a, several steps of the change in some physical quantity Z are shown. Following Lozovoi's notation, we will call the first step as "host removal" (HR), since it results in change in Z due to removal of a matrix atom, i.e., vacancy formation without optimization (system B in Figure 1a). The second step, "substitutional structure" (SS), is connected with a deformation of the undoped system with the unrelaxed vacancy so that the new structure corresponds to the doped structure with the unrelaxed vacancy (system C in Figure 1a). Finally, the value of Z is changed due to the impurity incorporation. The corresponding contribution we will call "chemical + compressed impurity" (CC). In Figure 1a, the undoped and doped systems are marked A and D, respectively. where E(Ti3Al + O) and E(Ti3Al) are the total energies of the alloy with and without oxygen, respectively, and E(O2) is the total energy of the oxygen molecule calculated in an empty cell with sides of 12 × 13 × 14 Å. Factor 1/2 corresponds to the absorption of one oxygen atom in the computational cell. The formation energy of a substitutional defect was estimated by following equation: where E(Ti3Al + Imp) and E(Ti3Al) are the total energies of the alloy with and without impurity, respectively, and E(Ti/Al) and E(Imp) is the total energy of the titanium/aluminum (in dependence on occupied sublattice) and the impurity metal in the ground state per atom, respectively. In accordance with the approach suggested by Lozovoi with co-authors [45] for embrittlement of grain boundaries, the effect of a substitutional (interstitial) defect on any physical quantity can be decomposed into three (two) contributions. In Figure 1a, several steps of the change in some physical quantity Z are shown. Following Lozovoi's notation, we will call the first step as "host removal" (HR), since it results in change in Z due to removal of a matrix atom, i.e., vacancy formation without optimization (system B in Figure 1a). The second step, "substitutional structure" (SS), is connected with a deformation of the undoped system with the unrelaxed vacancy so that the new structure corresponds to the doped structure with the unrelaxed vacancy (system C in Figure 1a). Finally, the value of Z is changed due to the impurity incorporation. The corresponding contribution we will call "chemical + compressed impurity" (CC). In Figure 1a, the undoped and doped systems are marked A and D, respectively. Summing up the above considerations, the contributions to ΔZ in the case of a substitutional impurity are: Summing up the above considerations, the contributions to ∆Z in the case of a substitutional impurity are: As a result, ∆Z can be expressed using Equations In the case of an interstitial impurity, the first (HR) and second (SS) stages should be replaced by a single "interstitial structure" (IS) (Figure 1b): and The interaction energy between two atoms X and Y was calculated using following expression: where E(Ti 3 Al + X + Y), E(Ti 3 Al + X), and E(Ti 3 Al + Y) are the total energy of the Ti 3 Al including both X and Y atoms, the X atom only, and the Y atom only, respectively, and E(Ti 3 Al) is the total energy of the perfect host crystal. The positive value of E int means that the interaction is repulsive, and negative means it is attractive. The larger the absolute value of E int , the stronger is the interaction.

Results and Discussion
3.1. Substitutional Defect Formation 3.1.1. Single Impurity Atom Figure 2 shows the calculated values of the substitutional defect formation energy for the Ti and Al sublattices in the case of elements of the VB-VIIB groups. It can be seen that all impurities prefer to occupy the Ti sublattice. This trend is in agreement with the results of earlier papers [33,35]. It should be noted that significant differences in specific values are due to the following factors: in [33], the chemical potential of Ti and Al was estimated based on thermodynamic equilibrium between the α 2 -Ti 3 Al and γ-TiAl alloys, while in [35], another approximation for the exchange-correlation potential was applied. The difference in the defect formation energy for the two sublattices reaches 0.7-2.0 eV (0.7-1.6 eV in [33,35]). This suggests that at moderate temperatures, all impurity atoms substitute for Ti with probability of~100%. In this connection, further discussion is given only in case of impurities on the Ti sublattice.
As a result, ΔZ can be expressed using Equations (3)-(5) as In the case of an interstitial impurity, the first (HR) and second (SS) stages should b replaced by a single "interstitial structure" (IS) (Figure 1b): The interaction energy between two atoms X and Y was calculated using following expression: , and E(Ti3Al + Y) are the total energy of the Ti3Al in cluding both X and Y atoms, the X atom only, and the Y atom only, respectively, and E(Ti3Al) is the total energy of the perfect host crystal. The positive value of Eint means tha the interaction is repulsive, and negative means it is attractive. The larger the absolut value of Eint, the stronger is the interaction.

Substitutional Defect Formation
3.1.1. Single Impurity Atom Figure 2 shows the calculated values of the substitutional defect formation energy for the Ti and Al sublattices in the case of elements of the VB-VIIB groups. It can be seen that all impurities prefer to occupy the Ti sublattice. This trend is in agreement with th results of earlier papers [33,35]. It should be noted that significant differences in specifi values are due to the following factors: in [33], the chemical potential of Ti and Al wa estimated based on thermodynamic equilibrium between the α2-Ti3Al and γ-TiAl alloys while in [35], another approximation for the exchange-correlation potential was applied The difference in the defect formation energy for the two sublattices reaches 0.7-2.0 eV (0.7-1.6 eV in [33,35]). This suggests that at moderate temperatures, all impurity atom substitute for Ti with probability of ~100%. In this connection, further discussion is given only in case of impurities on the Ti sublattice.  The calculated values of the HR, SS, and CC contributions to the substitutional defect formation energy in accordance with Equations (3)-(5) are presented in Table 1. Note that the SS contribution can be considered as a mechanical one (µ), since it is due to the difference in the size of the impurity and Ti atoms. It can be seen that this contribution is two orders of magnitude smaller than the contributions of the HR and CC mechanisms. Recall that the former one is associated with the breaking of chemical bonds when a titanium atom is removed, while the latter is associated with the formation of new impurity-matrix bonds. Therefore, the sum of HR and CC can be considered as a chemical contribution (χ) to the substitutional defect formation energy. Since HR does not depend on the impurity and represents the formation energy of the Ti vacancy without relaxation, it is convenient to use the chemical contribution χ alongside with the mechanical µ one in the further discussion. It is interesting that, taking into account the results of our earlier work [34], we can conclude that the relaxation energy of the α 2 -Ti 3 Al alloy to the energy of a vacancy is 0.25 eV, which is an order of magnitude greater than the mechanical contribution µ during the formation of a substitutional defect. Thus, E f is mainly determined by the difference in the binding energy between impurity-matrix and Ti-matrix, i.e., the chemical contribution ( Figure 3). The cases of exception are Nb and Re impurities, since for them the χ and µ contributions are comparable. fect formation energy in accordance with Equations (3)-(5) are presented in Table 1. that the SS contribution can be considered as a mechanical one (μ), since it is due t difference in the size of the impurity and Ti atoms. It can be seen that this contributi two orders of magnitude smaller than the contributions of the HR and CC mechani Recall that the former one is associated with the breaking of chemical bonds when a nium atom is removed, while the latter is associated with the formation of new impu matrix bonds. Therefore, the sum of HR and CC can be considered as a chemical co bution (χ) to the substitutional defect formation energy. Since HR does not depend o impurity and represents the formation energy of the Ti vacancy without relaxation convenient to use the chemical contribution χ alongside with the mechanical μ one i further discussion. It is interesting that, taking into account the results of our earlier w [34], we can conclude that the relaxation energy of the α2-Ti3Al alloy to the energy vacancy is ~0.25 eV, which is an order of magnitude greater than the mechanical co bution μ during the formation of a substitutional defect. Thus, Ef is mainly determine the difference in the binding energy between impurity-matrix and Ti-matrix, i.e. chemical contribution (Figure 3). The cases of exception are Nb and Re impurities, for them the χ and μ contributions are comparable.   Figure 4 shows the values of the contribution μ versus the difference (ΔRcov) i valent radii of impurity and titanium atoms [46]. It can be seen that μ is in good ag ment with ΔRcov, the correlation coefficient R is 0.88, and the standard deviation σ is 0 eV. At the same time, the chemical contribution is not only determined by the hybri tion contribution, which can be estimated by calculating the Crystal Orbital Ham Populations (COHP) [47,48]. Note that the COHP curve integrated up to the Fermi and, taken with the opposite sign (−ICOHP), can be considered as a measure of the bridization contribution. Thus, as seen from Figure 5, the correlation between χ and values of −ICOHP is weak enough: R equals −0.61 and σ = 0.216 eV. Unfortunately  Figure 4 shows the values of the contribution µ versus the difference (∆R cov ) in covalent radii of impurity and titanium atoms [46]. It can be seen that µ is in good agreement with ∆R cov , the correlation coefficient R is 0.88, and the standard deviation σ is 0.013 eV. At the same time, the chemical contribution is not only determined by the hybridization contribution, which can be estimated by calculating the Crystal Orbital Hamilton Populations (COHP) [47,48]. Note that the COHP curve integrated up to the Fermi level and, taken with the opposite sign (−ICOHP), can be considered as a measure of the hybridization contribution. Thus, as seen from Figure 5, the correlation between χ and the values of −ICOHP is weak enough: R equals −0.61 and σ = 0.216 eV. Unfortunately, it is quite difficult to estimate the metallic contribution to the chemical bonding, while the ionic contribution should be insignificant, since the electronegativity of impurities differs slightly from that of titanium. For example, the maximum difference reaches 0.62 eV 1/2 in the case of Mo and for other metals it does not exceed 0.36 eV 1/2 [46].
Metals 2022, 12, x FOR PEER REVIEW 6 quite difficult to estimate the metallic contribution to the chemical bonding, while ionic contribution should be insignificant, since the electronegativity of impurities di slightly from that of titanium. For example, the maximum difference reaches 0.62 eV the case of Mo and for other metals it does not exceed 0.36 eV 1/2 [46].

Two Impurity Atoms
In the general case for a lattice with the hexagonal symmetry, the coordinatio lipsoids should be considered instead of the coordination spheres [49]. In this regard nearest neighboring atoms (1NN) can be divided into two types: (i) atoms lying in same (0001) plane as the selected atom (1in) and (ii) atoms lying in neighboring pl (1out, see Figure 6a). In both groups there are six neighboring atoms (Table 2). In the of the second coordination ellipsoid (2NN), there are only six titanium atoms that l adjacent (0001) planes ( Figure 6b). As seen from Figure 6c, two atoms located in the jacent (0001) planes exactly above and below the selected atom are the third neigh (3NN). Finally, in the case of the fourth coordination ellipsoid, the atomic distribu manner is similar to that in the first coordination ellipsoid: there are six atoms in the s (0001) plane (4in) and twelve atoms in the adjacent (0001) planes (4out) (Figure 6d) data, including interatomic distances and multiplicities of atoms, are given in Tab Frankly speaking, since the value of 2c/a = 1.617 differed insignificantly from the valu   [46].

Two Impurity Atoms
In the general case for a lattice with the hexagonal symmetry, the coordinatio lipsoids should be considered instead of the coordination spheres [49]. In this regard nearest neighboring atoms (1NN) can be divided into two types: (i) atoms lying in same (0001) plane as the selected atom (1in) and (ii) atoms lying in neighboring pl (1out, see Figure 6a). In both groups there are six neighboring atoms (Table 2). In the of the second coordination ellipsoid (2NN), there are only six titanium atoms that l adjacent (0001) planes ( Figure 6b). As seen from Figure 6c, two atoms located in th jacent (0001) planes exactly above and below the selected atom are the third neigh (3NN). Finally, in the case of the fourth coordination ellipsoid, the atomic distribu manner is similar to that in the first coordination ellipsoid: there are six atoms in the s (0001) plane (4in) and twelve atoms in the adjacent (0001) planes (4out) (Figure 6d) data, including interatomic distances and multiplicities of atoms, are given in Tab Frankly speaking, since the value of 2c/a = 1.617 differed insignificantly from the val

Two Impurity Atoms
In the general case for a lattice with the hexagonal symmetry, the coordination ellipsoids should be considered instead of the coordination spheres [49]. In this regard, the nearest neighboring atoms (1NN) can be divided into two types: (i) atoms lying in the same (0001) plane as the selected atom (1 in ) and (ii) atoms lying in neighboring planes (1 out , see Figure 6a). In both groups there are six neighboring atoms ( Table 2). In the case of the second coordination ellipsoid (2NN), there are only six titanium atoms that lie in adjacent (0001) planes ( Figure 6b). As seen from Figure 6c, two atoms located in the adjacent (0001) planes exactly above and below the selected atom are the third neighbors (3NN). Finally, in the case of the fourth coordination ellipsoid, the atomic distribution manner is similar to that in the first coordination ellipsoid: there are six atoms in the same (0001) plane (4 in ) and twelve atoms in the adjacent (0001) planes (4 out ) (Figure 6d). All data, including interatomic distances and multiplicities of atoms, are given in Table 2. Frankly speaking, since the value of 2c/a = 1.617 differed insignificantly from the value of 1.633 for an ideal HCP structure, the coordination ellipsoid is quite close to a sphere, and there is an insignificant difference in the distance to atoms 1 in and 1 out as well as 4 in and 4 out .  The interaction energy between impurity atoms Eint was estimated as a function of interatomic distance. Calculations showed that in the case of pairs of identical atoms, regardless of their combination, the lowest interaction energy corresponds to the situation when the impurity atoms are in the second neighbors relative to each other ( Figure  7a). Moreover, only in the case of Mo, W, and Re is the interaction energy negative, which indicates an energetically preferable configuration of impurity atoms. The Cr, Nb, and Ta atoms prefer to be as far as possible from each other. A similar situation takes place in the case of pairs of different atoms (Figure 7b). The only exception is the Nb-Re pair, for which the interaction energy is negative. It is known that as the distance between atoms increases, the interaction energy should tend to zero. It can be seen from Figure 7 that, in the case when the impurity atoms are the sixth neighbors of each other, the interaction energy becomes less than 0.02 eV.  The interaction energy between impurity atoms E int was estimated as a function of interatomic distance. Calculations showed that in the case of pairs of identical atoms, regardless of their combination, the lowest interaction energy corresponds to the situation when the impurity atoms are in the second neighbors relative to each other (Figure 7a). Moreover, only in the case of Mo, W, and Re is the interaction energy negative, which indicates an energetically preferable configuration of impurity atoms. The Cr, Nb, and Ta atoms prefer to be as far as possible from each other. A similar situation takes place in the case of pairs of different atoms (Figure 7b). The only exception is the Nb-Re pair, for which the interaction energy is negative. It is known that as the distance between atoms increases, the interaction energy should tend to zero. It can be seen from Figure 7 that, in the case when the impurity atoms are the sixth neighbors of each other, the interaction energy becomes less than 0.02 eV.  Table 3 shows the calculated contributions to the interaction energy between impurity atoms at second neighbors' position. It can be seen that the mechanical and chemical contributions are comparable in magnitude. The negative values of Eint are found for V-V, Mo-Mo, W-W, Re-Re, and Nb-Re combinations. In the case of the Mo-Mo pair, both contributions are negative, the mechanical contribution is dominant for V-V, W-W, and Nb-Re, and the chemical one for Re-Re. In general, it can be seen from Table 3 that with an increase in the occupation of the impurity valence band the chemical contribution to the interaction energy increases in magnitude. The values of the correlation coefficient and the standard deviation are 0.78 and 0.037 eV, respectively. Thus, the dominance of the chemical contribution in the case of Re-Re can be explained by its high valence.   Table 3 shows the calculated contributions to the interaction energy between impurity atoms at second neighbors' position. It can be seen that the mechanical and chemical contributions are comparable in magnitude. The negative values of E int are found for V-V, Mo-Mo, W-W, Re-Re, and Nb-Re combinations. In the case of the Mo-Mo pair, both contributions are negative, the mechanical contribution is dominant for V-V, W-W, and Nb-Re, and the chemical one for Re-Re. In general, it can be seen from Table 3 that with an increase in the occupation of the impurity valence band the chemical contribution to the interaction energy increases in magnitude. The values of the correlation coefficient and the standard deviation are 0.78 and 0.037 eV, respectively. Thus, the dominance of the chemical contribution in the case of Re-Re can be explained by its high valence.

Impurity Influence on the Oxygen Absorption
The highest energy of oxygen absorption was found to be in the octahedral Ti-rich site (2a Wyckoff position) that is in agreement with the results of earlier experimental [50] and theoretical [51,52] papers. The oxygen absorption energy in other positions (4f, 6g, 6h, for Metals 2022, 12, 650 9 of 14 details see [52]) is by 1.5-2.3 eV lower. The probability (p i ) of oxygen to be absorbed in a specific site i can be estimated according to the Boltzmann distribution as where E abs (i) is the oxygen absorption energy in the site i, n i is the multiplicity of position i within the primitive cell, k B is the Boltzmann constant, and T is temperature; the summation is carried out over all absorption positions. Thus, at the moderate temperature all oxygen atoms should occupy 2a positions with 100% probability. Since there are only two 2a sites in the primitive cell of the alloy, the absorption in these sites is possible for an oxygen atom concentration of less than 25 at.%. Note that, in reality, this percentage is much lower due to the alloy oxidation accompanied by the formation of the titanium and aluminum oxides. In addition, the migration barrier from this position controls the oxygen diffusivity in the α 2 -Ti 3 Al alloy [38]. Since titanium has a high affinity for oxygen, the fact that all the considered impurities lead to a decrease in oxygen absorption energy is expected (Figure 8). In fact, the more distant from titanium the impurity element is in the periodic table, the stronger it lowers the oxygen absorption energy. It can be seen from Figure 8 that it is the chemical contribution that has the main effect. Even if some substituting atoms are smaller than titanium, the mechanical contribution is negative for all considered impurities. It results from the fact that the equilibrium length of the Me-O bond is larger than that of Ti-O. This explains the increase of µ (in magnitude) with the filling of the impurity d-shell. where Eabs(i) is the oxygen absorption energy in the site i, ni is the multiplicity of position i within the primitive cell, kB is the Boltzmann constant, and T is temperature; the summation is carried out over all absorption positions. Thus, at the moderate temperature all oxygen atoms should occupy 2a positions with 100% probability. Since there are only two 2a sites in the primitive cell of the alloy, the absorption in these sites is possible for an oxygen atom concentration of less than 25 at.%. Note that, in reality, this percentage is much lower due to the alloy oxidation accompanied by the formation of the titanium and aluminum oxides. In addition, the migration barrier from this position controls the oxygen diffusivity in the α2-Ti3Al alloy [38]. Since titanium has a high affinity for oxygen, the fact that all the considered impurities lead to a decrease in oxygen absorption energy is expected (Figure 8). In fact, the more distant from titanium the impurity element is in the periodic table, the stronger it lowers the oxygen absorption energy. It can be seen from Figure 8 that it is the chemical contribution that has the main effect. Even if some substituting atoms are smaller than titanium, the mechanical contribution is negative for all considered impurities. It results from the fact that the equilibrium length of the Me-O bond is larger than that of Ti-O. This explains the increase of μ (in magnitude) with the filling of the impurity d-shell. Taking into account the results from Section 3.1.2, the joint effect of two atoms will be considered in the second neighbor configuration. Obviously, the appearance of the second impurity atom near oxygen enhances the change in the absorption energy ( Figure  9). The dominant factor responsible for a decrease in the oxygen absorption energy is the chemical contribution. The sum of values of ΔEabs calculated for each atom separately (red and blue bars in Figure 9), shown by green diamonds in Figure 9, is quite close to the results of calculating ΔEabs in the presence of two impurity atoms (yellow bars in Figure  9). The maximum difference between estimated ΔEabs(Imp1) + ΔEabs(Imp2) and ab initio calculated ΔEabs(Imp1 + Imp2) values is obtained for the Re-Re pair and it is equal to 0.11 eV. Denoting this difference ε, ( Δ abs abs abs Taking into account the results from Section 3.1.2, the joint effect of two atoms will be considered in the second neighbor configuration. Obviously, the appearance of the second impurity atom near oxygen enhances the change in the absorption energy ( Figure 9). The dominant factor responsible for a decrease in the oxygen absorption energy is the chemical contribution. The sum of values of ∆E abs calculated for each atom separately (red and blue bars in Figure 9), shown by green diamonds in Figure 9, is quite close to the results of calculating ∆E abs in the presence of two impurity atoms (yellow bars in Figure 9). The maximum difference between estimated ∆E abs (Imp1) + ∆E abs (Imp2) and ab initio calculated ∆E abs (Imp1 + Imp2) values is obtained for the Re-Re pair and it is equal to 0.11 eV. Denoting this difference ε, ε = ∆E abs (Imp1) + ∆E abs (Imp2) − ∆E abs (Imp1 + Imp2).
(11) Metals 2022, 12, x FOR PEER REVIEW 10 of 14 Figure 9. Change in the oxygen absorption energy due to alloying by one and two impurity atoms.
We suppose that the value of ε should be correlated with the interaction energy between impurity atoms. Due to the presence of oxygen, these interaction energies differ from those earlier calculated and given in Table 3. A sufficiently high electronegativity of oxygen leads to O-Me bonds being highly ionic. Therefore, impurity atoms get a charge of the same sign that leads to a repulsion interaction. As a result, the interaction energy becomes positive for all considered pairs of atoms, which leads to an additional decrease in the oxygen absorption energy. Indeed, this can be seen from Figure 10. In this case, the correlation coefficient between the interaction energy and the energy difference ε is −0.82, with the standard deviation of 0.022 eV. Considering just the effects of each impurity atom separately, the effect of two impurity atoms can be predicted as the sum of two corresponding contributions (Figure 11). In the case of such an estimation, the maximal error is 0.11 eV, and the relative error is less than 9%. Thus, in the case of an α2-Ti3Al alloy doped with elements of the VB-VIIB groups, with an increase in the impurity concentration or the addition of a second alloying element, the following conclusions can be drawn: (i) the configuration with the impurity atoms substituting for Ti, and being the second neighbors relative to each other, is the energy preferable one; (ii) the oxygen absorption energy is decreased for the positions near the impurity and the effect becomes stronger with the filling of its d-shell; (iii) if oxygen locates in the region between impurity atoms, their effect is additively enhanced, which makes it possible to predict the joint effect of impurities with an accuracy of ~91%; Figure 9. Change in the oxygen absorption energy due to alloying by one and two impurity atoms.
We suppose that the value of ε should be correlated with the interaction energy between impurity atoms. Due to the presence of oxygen, these interaction energies differ from those earlier calculated and given in Table 3. A sufficiently high electronegativity of oxygen leads to O-Me bonds being highly ionic. Therefore, impurity atoms get a charge of the same sign that leads to a repulsion interaction. As a result, the interaction energy becomes positive for all considered pairs of atoms, which leads to an additional decrease in the oxygen absorption energy. Indeed, this can be seen from Figure 10. In this case, the correlation coefficient between the interaction energy and the energy difference ε is −0.82, with the standard deviation of 0.022 eV. Considering just the effects of each impurity atom separately, the effect of two impurity atoms can be predicted as the sum of two corresponding contributions (Figure 11). In the case of such an estimation, the maximal error is 0.11 eV, and the relative error is less than 9%. Figure 9. Change in the oxygen absorption energy due to alloying by one and two impuri We suppose that the value of ε should be correlated with the interaction between impurity atoms. Due to the presence of oxygen, these interaction energi from those earlier calculated and given in Table 3. A sufficiently high electronega oxygen leads to O-Me bonds being highly ionic. Therefore, impurity atoms get of the same sign that leads to a repulsion interaction. As a result, the interaction becomes positive for all considered pairs of atoms, which leads to an additional in the oxygen absorption energy. Indeed, this can be seen from Figure 10. In this correlation coefficient between the interaction energy and the energy difference ε with the standard deviation of 0.022 eV. Considering just the effects of each i atom separately, the effect of two impurity atoms can be predicted as the sum corresponding contributions (Figure 11). In the case of such an estimation, the m error is 0.11 eV, and the relative error is less than 9%. Thus, in the case of an α2-Ti3Al alloy doped with elements of the VB-VIIB with an increase in the impurity concentration or the addition of a second alloy ment, the following conclusions can be drawn: (i) the configuration with the i atoms substituting for Ti, and being the second neighbors relative to each othe energy preferable one; (ii) the oxygen absorption energy is decreased for the p near the impurity and the effect becomes stronger with the filling of its d-she oxygen locates in the region between impurity atoms, their effect is additively en Thus, in the case of an α 2 -Ti 3 Al alloy doped with elements of the VB-VIIB groups, with an increase in the impurity concentration or the addition of a second alloying element, the following conclusions can be drawn: (i) the configuration with the impurity atoms substituting for Ti, and being the second neighbors relative to each other, is the energy preferable one; (ii) the oxygen absorption energy is decreased for the positions near the impurity and the effect becomes stronger with the filling of its d-shell; (iii) if oxygen locates in the region between impurity atoms, their effect is additively enhanced, which makes it possible to predict the joint effect of impurities with an accuracy of~91%; (iv) regions near impurities can prevent the oxygen diffusion, since its penetration is energetically unfavorable. The latter is expected to be studied in more detail in our forthcoming papers.
Metals 2022, 12, x FOR PEER REVIEW (iv) regions near impurities can prevent the oxygen diffusion, since its penetrati ergetically unfavorable. The latter is expected to be studied in more detail in o coming papers. Figure 11. Comparison of ΔEabs(Imp1) + ΔEabs(Imp2) and ΔEabs(Imp1 + Imp2) values.
Finally, several remarks should be made. First, the above-mentioned res valid if impurity atoms can be freely placed in the first four coordination ellips cording to energy preference. This means that the total concentration of the im should not exceed ~5.3 at.% (2/38, where 38 is the total number of atoms in the f coordination ellipsoids). Moreover, treatment can also have an influence on the tion of the impurity atoms. Second, it is known that the α2-Ti3Al alloy can be dis transforming to α-structure. In this case, all obtained results remain valid, since amount of aluminum cannot block a coordination ellipsoid. Thirdly, a direct mental confirmation of the established trends is quite difficult. Nevertheless, i done indirectly. For example, a more or less uniform distribution of Nb atoms in substrate was observed in [27], whereas Y atoms showed a tendency to form Such distribution of Nb agrees with our conclusion that Nb atoms prefer to be possible from each other.

Conclusions
This manuscript is a theoretical study of the influence of the substitutional ties of elements from the VB-VIIB groups, such as V, Cr, Nb, Mo, Ta, W, and R oxygen absorption in the intermetallic α2-Ti3Al alloy. It is shown that the substitu a Ti atom by these impurities is more preferential in energy than substituting A former case, the defect formation energy spans the range from −0.11 to 0.76 eV, the latter case it is 1.46-2.15 eV. For both sublattices, the highest values of th formation energy correspond to chromium. It was found that this effect is mainl the chemical contribution to the substitutional defect formation energy. All con the impurities result in a decrease in the oxygen absorption energy. Moreover, chemical contribution that is responsible for this effect.
The interaction energy between two impurity atoms was calculated as a fun their relative positions. It was shown that the structure with impurities as the neighbors to each other is preferred one. At the same time, only for Mo-Mo, W Figure 11. Comparison of ∆E abs (Imp1) + ∆E abs (Imp2) and ∆E abs (Imp1 + Imp2) values.
Finally, several remarks should be made. First, the above-mentioned results are valid if impurity atoms can be freely placed in the first four coordination ellipsoids according to energy preference. This means that the total concentration of the impurities should not exceed~5.3 at.% (2/38, where 38 is the total number of atoms in the first four coordination ellipsoids). Moreover, treatment can also have an influence on the distribution of the impurity atoms. Second, it is known that the α 2 -Ti 3 Al alloy can be disordered, transforming to α-structure. In this case, all obtained results remain valid, since a small amount of aluminum cannot block a coordination ellipsoid. Thirdly, a direct experimental confirmation of the established trends is quite difficult. Nevertheless, it can be done indirectly. For example, a more or less uniform distribution of Nb atoms in the alloy substrate was observed in [27], whereas Y atoms showed a tendency to form clusters. Such distribution of Nb agrees with our conclusion that Nb atoms prefer to be as far as possible from each other.

Conclusions
This manuscript is a theoretical study of the influence of the substitutional impurities of elements from the VB-VIIB groups, such as V, Cr, Nb, Mo, Ta, W, and Re, on the oxygen absorption in the intermetallic α 2 -Ti 3 Al alloy. It is shown that the substitution for a Ti atom by these impurities is more preferential in energy than substituting Al. In the former case, the defect formation energy spans the range from −0.11 to 0.76 eV, while in the latter case it is 1.46-2.15 eV. For both sublattices, the highest values of the defect formation energy correspond to chromium. It was found that this effect is mainly due to the chemical contribution to the substitutional defect formation energy. All considered, the impurities result in a decrease in the oxygen absorption energy. Moreover, it is the chemical contribution that is responsible for this effect.
The interaction energy between two impurity atoms was calculated as a function of their relative positions. It was shown that the structure with impurities as the second neighbors to each other is preferred one. At the same time, only for Mo-Mo, W-W, Re-Re, and Nb-Re pairs is this energy is negative. In other words, most of the impurity atoms considered prefer to be as far as possible from each other. The energy preference of the second neighbor configuration is conditioned by the mechanical contribution in the case of W-W and Nb-Re, by the chemical contribution for Re-Re, and both contributions in the case of the Mo-Mo interaction. It was established that the effect of two impurity atoms on the oxygen absorption energy can be estimated as the sum of the effects of the separate impurities. Error in such estimation does not exceed 9%. The latter can be useful in modeling of effect of impurity combinations on the oxygen temperature-dependent diffusion coefficient.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.