First-Principles Calculations of High-Pressure Physical Properties of Ti0.5Ta0.5 Alloy

In this paper, an in-depth theoretical study on some physical properties of Ti0.5Ta0.5 alloy with systematic symmetry under high pressure is conducted via first-principles calculations, and relevant physical parameters are calculated. The results demonstrate that the calculated parameters, including lattice parameter, elastic constants, and elastic moduli, fit well with available theoretical and experimental data when the Ti0.5Ta0.5 alloy is under T = 0 and P = 0, indicating that the theoretical analysis method can effectively predict the physical properties of the Ti0.5Ta0.5 alloy. The microstructure and macroscopic physical properties of the alloy cannot be destroyed as the applied pressure ranges from 0 to 50GPa, but the phase transition of crystal structure may occur in the Ti0.5Ta0.5 alloy if the applied pressure continues to increase according to the TDOS curves and charge density diagram. The value of Young’s and shear modulus is maximized at P = 25 GPa. The anisotropy factors A(100)[001] and A(110)[001] are equal to 1, suggesting the Ti0.5Ta0.5 alloy is an isotropic material at 28 GPa, and the metallic bond is strengthened under high pressure. The present results provide helpful insights into the physical properties of Ti0.5Ta0.5 alloy.


Introduction
Ti alloys have many outstanding properties, including high strength, low density, corrosion resistance, thermostability, and superior biocompatibility, and thus are extensively applied into aerospace, marine, chemical, and biological fields [1][2][3][4][5]. Among them, it has been found that TiTa alloy is one of the best biomaterials since it is cheaper than pure Ta and more corrosion-resistant than pure Ti; hence, it has attracted more and more attention [6][7][8]. Zhou et al. [9] found that TiTa alloys showed different microstructures with respect to various Ta contents, such as the α phase of lamellar HCP structure under 20% Ta contents, the α phase of needle-like orthorhombic structure between 30%-50% Ta contents, the β + α phase of metastable structure for Ta contents of 60%, along with a single metastable β phase beyond Ta contents of 60%. Yin et al. [10] found that in Ti-28 Ta at % alloy, the cooling rate had a remarkable impact on the percentage of martensitic phase α and metastable phase β, and α phases tended to be induced at grain boundaries. Via the powder metallurgy method, Dercz et al. [11] evaluated the possibilities of synthesis of porous Ti-50Ta alloys; the microstructure and porosity of Ti-50 Ta alloys were greatly affected by the sintering time, and the corrosion resistance of sintered Ti-50 Ta alloys was better than that of pure Ti or Ta. Sing et al. [12] reported that TiTa alloy was composed of β phase with the help of its rapid solidification and stabilizing effect, and the strength was higher with respect to pure Ti and Ti6Al4V. Behera et al. [13] measured the enthalpy change of Ti-xTa alloys at the temperatures ranging 463 to 1257 K, and then clearly revealed that the enthalpy came mainly from the contributions of two parts: (i) one part from untransformed α and coexisting β phases, (ii) another part from the diffusional phase transformation from α phase to β phase. Furthermore, plenty of studies on TiTa alloys were also carried out by others [14][15][16][17].
On the basis of the above studies, plenty of physical properties of Ti 0.5 Ta 0.5 alloy were revealed through theory and experiment, as far as we know, but high-pressure physical properties in Ti 0.5 Ta 0. 5 alloy have not been studied yet, such as the structural, mechanical, and electronic properties. Therefore, the aim of this paper is to apply the first-principles calculations to investigate the high-pressure physical properties in Ti 0.5 Ta 0.5 alloy in the frame of density functional theory (DFT), and related physical parameters are computed to reveal the relationships between these parameters and applied pressure acting on Ti 0.5 Ta 0.5 alloy. In this paper, the main framework is listed below: In Section 2, the theoretical methodology and design parameters are introduced in detail. In Section 3, the calculated parameters of Ti 0.5 Ta 0.5 alloy are presented and discussed. Finally, a brief summary is provided in Section 4.

Theoretical Methodology
In this work, all DFT energy calculations are carried out via the Cambridge Serial Total Energy Package Program [18][19][20]. Herein, we adopt the generalized-gradient approximation (GGA) within the Perdew-Burke-Ernzerhof (PBE) as the exchange-correlation functional for calculations [21]. The pseudopotential is described by Vanderbilt-type ultrasoft pseudopotentials (USPP) to decide the interatomic interactions, which is determined by the shell electrons, namely, the valence electrons [22]. The valence electrons of Ti and Ta elements are 4s 2 3d 2 and 6s 2 5d 3 , respectively. To ensure the convergence accuracy of the electronic calculations of the Ti 0.5 Ta 0.5 alloy, the k-mesh in the Brillouin zone is set as 13 × 13 × 13, and the cutoff energy is optimized as 400 eV. The energy convergence criterion in the self-consistent calculation is optimized for 1.0 × 10 -6 eV/atom. For Ti 0.5 Ta 0.5 alloy, it has the space group of Im-3m, and Figure 1 displays the crystal structure. The algorithm of Broyden-Fletcher-Goldfarb-Shanno (BFGS) [23] is selected for optimizing the geometric configuration of the alloy in the pressure range of -10 to 50 GPa. Meanwhile, the Hellmann-Feyman force of each atom is accurate to 0.01 eV/Å. clearly revealed that the enthalpy came mainly from the contributions of two parts: (i) one part from untransformed  and coexisting  phases, (ii) another part from the diffusional phase transformation from  phase to  phase. Furthermore, plenty of studies on TiTa alloys were also carried out by others [14][15][16][17]. On the basis of the above studies, plenty of physical properties of Ti0.5Ta0.5 alloy were revealed through theory and experiment, as far as we know, but high-pressure physical properties in Ti0.5Ta0.5 alloy have not been studied yet, such as the structural, mechanical, and electronic properties. Therefore, the aim of this paper is to apply the first-principles calculations to investigate the highpressure physical properties in Ti0.5Ta0.5 alloy in the frame of density functional theory (DFT), and related physical parameters are computed to reveal the relationships between these parameters and applied pressure acting on Ti0.5Ta0.5 alloy. In this paper, the main framework is listed below: In Section 2, the theoretical methodology and design parameters are introduced in detail. In Section 3, the calculated parameters of Ti0.5Ta0.5 alloy are presented and discussed. Finally, a brief summary is provided in Section 4.

Theoretical Methodology
In this work, all DFT energy calculations are carried out via the Cambridge Serial Total Energy Package Program [18][19][20]. Herein, we adopt the generalized-gradient approximation (GGA) within the Perdew-Burke-Ernzerhof (PBE) as the exchange-correlation functional for calculations [21]. The pseudopotential is described by Vanderbilt-type ultrasoft pseudopotentials (USPP) to decide the interatomic interactions, which is determined by the shell electrons, namely, the valence electrons [22]. The valence electrons of Ti and Ta elements are 4s 2 3d 2 and 6s 2 5d 3 , respectively. To ensure the convergence accuracy of the electronic calculations of the Ti0.5Ta0.5 alloy, the k-mesh in the Brillouin zone is set as 13 × 13 × 13, and the cutoff energy is optimized as 400 eV. The energy convergence criterion in the self-consistent calculation is optimized for 1.0 × 10 -6 eV/atom. For Ti0.5Ta0.5 alloy, it has the space group of Im-3m, and Figure 1 displays the crystal structure. The algorithm of Broyden-Fletcher-Goldfarb-Shanno (BFGS) [23] is selected for optimizing the geometric configuration of the alloy in the pressure range of -10 to 50 GPa. Meanwhile, the Hellmann-Feyman force of each atom is accurate to 0.01 eV/Å.

Structure and Stability
In the light of the set parameters, the adequate structural optimization of the unit cell is first performed to get the theoretical E-V data, herein, E denotes total energy and V indicates the volume of the unit cell, and we fit the E-V data to determine the equilibrium crystal structure of Ti0.5Ta0.5 alloy according to the Birch-Murnaghan equation of state, as depicted in Figure 2, in which the volume range is set to 0.9-1.1 0 V , and each unit cell is optimized sufficiently. From the variation curve, the total energy gets the minimum value (

Structure and Stability
In the light of the set parameters, the adequate structural optimization of the unit cell is first performed to get the theoretical E-V data, herein, E denotes total energy and V indicates the volume of the unit cell, and we fit the E-V data to determine the equilibrium crystal structure of Ti 0.5 Ta 0.5 alloy according to the Birch-Murnaghan equation of state, as depicted in Figure 2, in which the volume range is set to 0.9-1.1 V 0 , and each unit cell is optimized sufficiently. From the variation curve, the total energy gets the minimum value (E t = −19.401 eV) when the volume is V 0 = 34.660 Å 3 , which shows that the crystal structure reaches the most stable state when these physical parameters are taken, and then the theoretical prediction of corresponding lattice constant is obtained by derivation, namely, a 0 = 3.260 Å, where V 0 and a 0 represents the primitive cell volume and equilibrium lattice constant under T = 0 and P = 0, respectively. Meanwhile, Table 1 lists the comparative results of the calculated lattice constant a 0 versus other data for Ti 0.5 Ta 0.5 alloy, and it fits well with these available data, indicating the validity of the theoretical analysis method.   A series of structural optimizations are performed on the crystal structure of the Ti0.5Ta0.5 alloy under different applied pressures to get the corresponding volumes and equilibrium lattice constants, which are used to investigate the impacts of applied pressure on the unit cell volume and lattice constant. Subsequently, the dependencies of dimensionless ratios 0 V V and 0 a a on the applied pressure are obtained in the range of −10 to 50 GPa, where the negative pressure denotes the tension, as shown in Figure 3. From the variation curves, it can be found that the two ratios decrease when applied pressure increases and the descent speed for 0 V V is significantly faster compared to 0 a a , suggesting that applying pressure can greatly reduce interatomic distance, and then strengthen electron interactions between adjacent atoms in the Ti0.5Ta0.5 alloy.  A series of structural optimizations are performed on the crystal structure of the Ti 0.5 Ta 0.5 alloy under different applied pressures to get the corresponding volumes and equilibrium lattice constants, which are used to investigate the impacts of applied pressure on the unit cell volume and lattice constant. Subsequently, the dependencies of dimensionless ratios V/V 0 and a/a 0 on the applied pressure are obtained in the range of −10 to 50 GPa, where the negative pressure denotes the tension, as shown in Figure 3. From the variation curves, it can be found that the two ratios decrease when applied pressure increases and the descent speed for V/V 0 is significantly faster compared to a/a 0 , suggesting that applying pressure can greatly reduce interatomic distance, and then strengthen electron interactions between adjacent atoms in the Ti 0.5 Ta 0.5 alloy.
As we all know, elastic constants are the important physical parameters in influencing mechanical properties for anisotropic materials, and they are often used for assessing structural stability of crystalline materials and obtain relevant physical quantities, such as elastic moduli, anisotropy factors, Poisson's ratio, and so on, thereby determining the deformation resistance of materials. For cubic structures, these three elastic constants C 11 , C 12 , and C 44 are usually used for estimating their structural stability, and they must satisfy the corresponding stability criterion, namely, C 11 > 0, C 44 > 0, C 11 − C 12 > 0, and C 11 + 2C 12 > 0 [26]. Through calculation, the theoretical predictions for them are presented, and Table 2 lists the theoretical calculation results for comparisons with other theoretical data at T = 0 and P = 0, indicating that they agree well with other theoretical results [24,25].  As we all know, elastic constants are the important physical parameters in influencing mechanical properties for anisotropic materials, and they are often used for assessing structural stability of crystalline materials and obtain relevant physical quantities, such as elastic moduli, anisotropy factors, Poisson's ratio, and so on, thereby determining the deformation resistance of materials. For cubic structures, these three elastic constants 11 C , 12 C , and 44 C are usually used for estimating their structural stability, and they must satisfy the corresponding stability criterion, namely,  [26]. Through calculation, the theoretical predictions for them are presented, and Table 2 lists the theoretical calculation results for comparisons with other theoretical data at 0 T  and 0 P  , indicating that they agree well with other theoretical results [24,25]. Table 2. Comparisons of three elastic constants with available theoretical data in Ti0.5Ta0.5 alloy.

Ti0.5Ta0.5 Alloy Present
Theoretical Data 180.34 163.40 [24], 181.80 [25] 12 C (GPa) 142.16 132.80 [24], 138.86 [25] 44 C (GPa) 30.14 39.00 [24], 45.18 [25] At the same time, Figure 4 demonstrates the variation curves of elastic constants ij C in Ti0.5Ta0.5 alloy under the applied pressure ranging 0 to 50 GPa. Results display that 11 C and 12 C have the same trend that the values of theoretical predictions increase monotonously with the increase of applied pressure, and the elastic constant 44 C increases slowly and then decreases gradually, where the increase of elastic constants indicates that the deformation resistance of Ti0.5Ta0.5 alloy becomes stronger. Meanwhile, the calculated values of ij C always satisfy the stability criterion of crystal structure in the range of applied pressure, suggesting that the high-pressure mechanical stability of Ti0.5Ta0.5 alloy cannot be destroyed.  30.14 39.00 [24], 45.18 [25] At the same time, Figure 4 demonstrates the variation curves of elastic constants C ij in Ti 0.5 Ta 0.5 alloy under the applied pressure ranging 0 to 50 GPa. Results display that C 11 and C 12 have the same trend that the values of theoretical predictions increase monotonously with the increase of applied pressure, and the elastic constant C 44 increases slowly and then decreases gradually, where the increase of elastic constants indicates that the deformation resistance of Ti 0.5 Ta 0.5 alloy becomes stronger. Meanwhile, the calculated values of C ij always satisfy the stability criterion of crystal structure in the range of applied pressure, suggesting that the high-pressure mechanical stability of Ti 0.5 Ta 0.5 alloy cannot be destroyed.

Mechanical Properties
Ti alloys have outstanding mechanical performances, such as high strength, wear resistance, and nontoxicity, and they are extensively used in modern industries, which are closely connected with the material modulus against deformation, such as Young's modulus E, shear modulus G, and bulk modulus B [27]. Herein, the elastic moduli of materials with cubic crystal can be obtained through

Mechanical Properties
Ti alloys have outstanding mechanical performances, such as high strength, wear resistance, and nontoxicity, and they are extensively used in modern industries, which are closely connected with the material modulus against deformation, such as Young's modulus E, shear modulus G, and bulk modulus B [27]. Herein, the elastic moduli of materials with cubic crystal can be obtained through elastic constants C ij according to the research work of Iotova et al. [28], and the specific expressions are represents the Voigt modulus of shear resistance, and denotes the Reuss modulus of shear resistance.
Through calculation, the theoretical predictions of elastic moduli are obtained, and Figure 5 depicts the corresponding variation curves for the Ti 0.5 Ta 0.5 alloy under the applied pressure in the range of 0 to 50 GPa. Obviously, bulk modulus B increases gradually as applied pressure increases, which indicates that high pressure for Ti 0.5 Ta 0.5 alloy can enhance its resistance to volume change. E and G first increase and then decrease gradually, and both of them get the maximum values at P = 25 GPa, which indicates that increasing pressure can improve the resistances to elastic and shear deformations between 0 and 25 GPa, but decline their resistances beyond 25 GPa. Furthermore, we compare the theoretical predictions of elastic moduli with available theoretical and experimental data for Ti 0.5 Ta 0.5 alloy under T = 0 and P = 0, and Table 3 shows that the calculated data of this work are well consistent with the research results of others [7,12,25].  Researches show that putting pressure on materials will affect the ductile-brittle transition of crystal structure, thereby influencing the mechanical properties of materials [29,30]. Pugh proposed a significant conclusion that the brittleness or ductility of materials depends on the modulus ratio B G [27], and 1.75 B G  denotes the ductility, and suggesting that the Ti0.5Ta0.5 alloy has good ductility. Additionally, there is a critical pressure P = 25 GPa; namely, the ductility of Ti0.5Ta0.5 alloy is almost invariable at the pressure between 0 and 25 GPa, whereas the ductility of the alloy significantly increases when the applied pressure further increases.  Researches show that putting pressure on materials will affect the ductile-brittle transition of crystal structure, thereby influencing the mechanical properties of materials [29,30]. Pugh proposed a significant conclusion that the brittleness or ductility of materials depends on the modulus ratio B/G [27], and B/G > 1.75 denotes the ductility, and B/G < 1.75 indicates the brittleness. In light of Symmetry 2020, 12, 796 6 of 12 the calculated elastic moduli, Figure 6 depicts the changing curve of the modulus ratio B/G under pressure changing from 0 to 50 GPa. The result shows B/G > 1.75 at any pressure, suggesting that the Ti 0.5 Ta 0.5 alloy has good ductility. Additionally, there is a critical pressure P = 25 GPa; namely, the ductility of Ti 0.5 Ta 0.5 alloy is almost invariable at the pressure between 0 and 25 GPa, whereas the ductility of the alloy significantly increases when the applied pressure further increases.
Young's modulus E 71.44 93.71 [25] 69.00 [7], 71.73 [12] Shear modulus G 25.10 33.51 [25] Researches show that putting pressure on materials will affect the ductile-brittle transition of crystal structure, thereby influencing the mechanical properties of materials [29,30]. Pugh proposed a significant conclusion that the brittleness or ductility of materials depends on the modulus ratio B G [27], and 1.75 B G  denotes the ductility, and suggesting that the Ti0.5Ta0.5 alloy has good ductility. Additionally, there is a critical pressure P = 25 GPa; namely, the ductility of Ti0.5Ta0.5 alloy is almost invariable at the pressure between 0 and 25 GPa, whereas the ductility of the alloy significantly increases when the applied pressure further increases.

Anisotropy
In crystalline materials, the elastic anisotropy of crystal structure influences the mechanical properties of materials. At present, researchers usually use the anisotropy factor A to investigate this physical characterize, and materials with A = 1 are isotropic, otherwise anisotropic [31,32]. In 1986, the impacts of elastic anisotropy on mechanical properties of materials had been thoroughly studied by Yoo [33], he found an interesting phenomenon that a large A can enhance the tangential force to drive the cross slip of screw dislocations, and then proposed a significant model, namely, the cross-slip and pinning model. A (110)[001] = C 44 (C 11 + 2C 12 + C )/ C 11 C − C 2 12 , in which C = C 44 + (C 11 + C 12 )/2 [34,35]. Through calculation, Figure 7 plots the variation curves of anisotropy factors of Ti 0.5 Ta 0.5 alloy under applied pressure. The values of two anisotropy factors are equal to 1 when the applied pressure is about 28 GPa, indicating that the Ti 0.5 Ta 0.5 alloy is an isotropic material at 28 GPa. In the light of the trend of variation curves, the elastic anisotropy of Ti 0.5 Ta 0.5 alloy decreases with increasing pressure between 0 and 28 GPa, but increases under applied pressure ranging from 28 to 50 GPa, thereby the tangential force decreases first and then increases in promoting the cross slip of screw dislocations.
Poisson's ratio σ is a significant physical parameter that can be used to determine the plasticity of materials; thus, we investigate the change of Poisson's ratio σ of Ti 0.5 Ta 0.5 alloy under the applied pressure in this paper. For Poisson's ratio σ, the value changes from −1 to 0.5, and a large value indicates good plasticity. Reed et al. [36] put forward an important conclusion that σ = 0.25 and σ = 0.5 were the smallest and largest bounds for central-force solids, respectively. Additionally, Fu et al. [35] found that the magnitude of σ was used for defining the sort of interatomic bonding. In this paper, Poisson's ratios σ [001] and σ [111] in the [001] and [111] crystallographic directions are calculated according to the theoretical predictions of elastic constants C ij , and the calculation formulas are σ [001] = C 12 /(C 11 + C 12 ) and σ [111] = (C 11 + 2C 12 − 2C 44 )/2(C 11 + 2C 12 + C 44 ) [37]. by Yoo [33], he found an interesting phenomenon that a large A can enhance the tangential force to drive the cross slip of screw dislocations, and then proposed a significant model, namely, the crossslip and pinning model. 2 C C C C     [34,35].
Through calculation, Figure 7 plots the variation curves of anisotropy factors of Ti0.5Ta0.5 alloy under applied pressure. The values of two anisotropy factors are equal to 1 when the applied pressure is about 28 GPa, indicating that the Ti0.5Ta0.5 alloy is an isotropic material at 28 GPa. In the light of the trend of variation curves, the elastic anisotropy of Ti0.5Ta0.5 alloy decreases with increasing pressure between 0 and 28 GPa, but increases under applied pressure ranging from 28 to 50 GPa, thereby the tangential force decreases first and then increases in promoting the cross slip of screw dislocations. Poisson's ratio  is a significant physical parameter that can be used to determine the plasticity of materials; thus, we investigate the change of Poisson's ratio  of Ti0.5Ta0.5 alloy under the applied pressure in this paper. For Poisson's ratio  , the value changes from −1 to 0.5, and a large value indicates good plasticity. Reed et al. [36] put forward an important conclusion that 0.25   and 0.5

 
were the smallest and largest bounds for central-force solids, respectively. Additionally, Fu et al. [35] found that the magnitude of  was used for defining the sort of interatomic bonding. In      [32,38], and Figure 9a draws the variation curves of material moduli under applied pressure ranging from 0 to 50 GPa. It shows that the total varying trends of  Furthermore, we calculate the other moduli of Ti 0.5 Ta 0.5 alloy, including G (100) [010] , G (110) [110] , and E <100> , to further study the related mechanical performances, where G (100)[010] represents the shear modulus in the (100)[010] direction, and G (110) [110] indicates the one in the (110) 110 direction, respectively, and E <100> denotes Young's modulus along the < 100 > direction. These parameters can also be obtained by elastic constants C ij , and the calculation formulas are G (100)[010] = C 44 , G (110)[110] = (C 11 − C 12 )/2, and E <100> = C 11 − C 12 + C 11 C 12 − C 2 12 /(C 11 + C 12 ) [32,38], and Figure 9a draws the variation curves of material moduli under applied pressure ranging from 0 to 50 GPa. It shows that the total varying trends of G (110) [110] and E <100> increase gradually as the applied pressure increases, although E <100> decreases slightly at 50 GPa, indicating that high pressure strengthens the resistances to shear and elastic deformations of Ti 0.5 Ta 0.5 alloy in the (110) 110 and < 100 > crystallographic directions. Different from G (110) [110] and E <100> , shear modulus G (100)[010] initially increases and then decreases when applied pressure increases, and reaches the maximum value at 25 GPa, which shows that the Ti 0.5 Ta 0.5 alloy has the strongest shear resistance along the (100)[010] direction under P = 25 GPa, while the shear deformation resistance of the alloy will decline if the pressure continues to increase. Meanwhile, another important physical parameter, Cauchy pressure 12

C C 
, is calculated for evaluating the bonding characteristics from the atomic level under high pressure. Herein, the positive value of Cauchy pressure stands for the metallic characteristics of atomic bonding, indicating that there is a uniform electron gas near the spherical atom, which is composed of electrons without directionality and regionality. However, the atomic bonding exhibits the directional characteristic as the Cauchy pressure is negative, and a large value denotes the strong directional characteristic of atomic bonding [32,39]. According to elastic constants ij C , Figure 9b shows the changing curve of 12 44 C C  versus different pressure in Ti0.5Ta0.5 alloy. It is noted from the curve that the value of 12 44 C C  is positive at any pressures, and quickly increases for the Ti0.5Ta0.5 alloy between 0 and 50 GPa, suggesting that the metallic bond is the main bonding form of Ti0.5Ta0.5 alloy, and the metallic bond can be strengthened in the case of high pressure.

Electronic Properties
The valence electrons determine the bonding form between atoms, thereby influencing the microstructure and macroscopic physical properties of the crystal. Therefore, the electronic structures are given and discussed to study the structural stability of Ti0.5Ta0.5 alloy under applied pressure in this paper. Figure 10 exhibits the changing curves of the partial density of states (PDOS) and the total density of states (TDOS) at P = 0 GPa, where the red dash line stands for the Fermi level ( F 0 eV E  ). For the TDOS curve, it can be found that the value is not zero at F 0 eV E  , suggesting that the Ti0.5Ta0.5 alloy exhibits the metallicity, which verifies the conclusion in Figure 9b. For PDOS curves, they show that the Ti-3d and Ta-5d states have the greatest effect on DOS at 0 eV, while the effects of the Ti-4s and Ta-6s states on DOS are almost negligible. Meanwhile, another important physical parameter, Cauchy pressure C 12 − C 44 , is calculated for evaluating the bonding characteristics from the atomic level under high pressure. Herein, the positive value of Cauchy pressure stands for the metallic characteristics of atomic bonding, indicating that there is a uniform electron gas near the spherical atom, which is composed of electrons without directionality and regionality. However, the atomic bonding exhibits the directional characteristic as the Cauchy pressure is negative, and a large value denotes the strong directional characteristic of atomic bonding [32,39]. According to elastic constants C ij , Figure 9b shows the changing curve of C 12 − C 44 versus different pressure in Ti 0.5 Ta 0.5 alloy. It is noted from the curve that the value of C 12 − C 44 is positive at any pressures, and quickly increases for the Ti 0.5 Ta 0.5 alloy between 0 and 50 GPa, suggesting that the metallic bond is the main bonding form of Ti 0.5 Ta 0.5 alloy, and the metallic bond can be strengthened in the case of high pressure.

Electronic Properties
The valence electrons determine the bonding form between atoms, thereby influencing the microstructure and macroscopic physical properties of the crystal. Therefore, the electronic structures are given and discussed to study the structural stability of Ti 0.5 Ta 0.5 alloy under applied pressure in this paper. Figure 10 exhibits the changing curves of the partial density of states (PDOS) and the total density of states (TDOS) at P = 0 GPa, where the red dash line stands for the Fermi level (E F = 0 eV). For the TDOS curve, it can be found that the value is not zero at E F = 0 eV, suggesting that the Ti 0.5 Ta 0.5 alloy exhibits the metallicity, which verifies the conclusion in Figure 9b. For PDOS curves, they show that the Ti-3d and Ta-5d states have the greatest effect on DOS at 0 eV, while the effects of the Ti-4s and Ta-6s states on DOS are almost negligible. Meanwhile, Figure 11 describes the changing curves of TDOS versus pressures of 0, 15, 25, 35, and 50 GPa. From these changing curves, the results show that energy of the valence-band maximum increases as applied pressure increases, while the one of the conduction-band bottoms was reversed, indicating that high pressure decreases the band gap between the conduction band and valence band, and the outer electrons tend to easily transfer from the valence-band maximum to the conductionband bottom. Hence, high pressure may destroy the structural stability of Ti0.5Ta0.5 alloy, which affects the microstructure and macroscopic physical properties of Ti0.5Ta0.5 alloy. At the same time, according to the results of Figure 4, it can be found that the value of elastic constant 44 C may become negative as the applied pressure increases further, which indicates that the Ti0.5Ta0.5 alloy may produce the structural phase transition under higher pressure, and the analysis result is consistent with the previous conclusion in Figure 11.  Meanwhile, Figure 11 describes the changing curves of TDOS versus pressures of 0, 15, 25, 35, and 50 GPa. From these changing curves, the results show that energy of the valence-band maximum increases as applied pressure increases, while the one of the conduction-band bottoms was reversed, indicating that high pressure decreases the band gap between the conduction band and valence band, and the outer electrons tend to easily transfer from the valence-band maximum to the conduction-band bottom. Hence, high pressure may destroy the structural stability of Ti 0.5 Ta 0.5 alloy, which affects the microstructure and macroscopic physical properties of Ti 0.5 Ta 0.5 alloy. At the same time, according to the results of Figure 4, it can be found that the value of elastic constant C 44 may become negative as the applied pressure increases further, which indicates that the Ti 0.5 Ta 0.5 alloy may produce the structural phase transition under higher pressure, and the analysis result is consistent with the previous conclusion in Figure 11. Meanwhile, Figure 11 describes the changing curves of TDOS versus pressures of 0, 15, 25, 35, and 50 GPa. From these changing curves, the results show that energy of the valence-band maximum increases as applied pressure increases, while the one of the conduction-band bottoms was reversed, indicating that high pressure decreases the band gap between the conduction band and valence band, and the outer electrons tend to easily transfer from the valence-band maximum to the conductionband bottom. Hence, high pressure may destroy the structural stability of Ti0.5Ta0.5 alloy, which affects the microstructure and macroscopic physical properties of Ti0.5Ta0.5 alloy. At the same time, according to the results of Figure 4, it can be found that the value of elastic constant 44 C may become negative as the applied pressure increases further, which indicates that the Ti0.5Ta0.5 alloy may produce the structural phase transition under higher pressure, and the analysis result is consistent with the previous conclusion in Figure 11.  Finally, the charge density diagrams of isosurface contours versus various pressures are presented to understand the local chemical bonding in the Ti 0.5 Ta 0.5 alloy, as shown in Figure 12, where the isosurface levels are set to 0.0425 r −3 0 (r 0 is the Bohr radius). With increasing applied pressure, the local chemical bonding between Ti and Ta atoms gradually increases, implying that the high pressure can enhance the electron interaction between the two atoms, thereby increasing the local chemical bonding between them, and the change increases as applied pressure increases. However, the strong interactions between two atoms may destroy the stability of the crystal structure. Therefore, the structural phase transition may occur in the Ti 0.5 Ta 0.5 alloy under higher pressure beyond 50 GPa, which coincides with the result in Figure 11. Finally, the charge density diagrams of isosurface contours versus various pressures are presented to understand the local chemical bonding in the Ti0.5Ta0.5 alloy, as shown in Figure 12, where the isosurface levels are set to 0.0425 3 0 r  ( 0 r is the Bohr radius). With increasing applied pressure, the local chemical bonding between Ti and Ta atoms gradually increases, implying that the high pressure can enhance the electron interaction between the two atoms, thereby increasing the local chemical bonding between them, and the change increases as applied pressure increases. However, the strong interactions between two atoms may destroy the stability of the crystal structure. Therefore, the structural phase transition may occur in the Ti0.5Ta0.5 alloy under higher pressure beyond 50 GPa, which coincides with the result in Figure 11.

Conclusions
In conclusion, DFT calculations are implemented to study the influences of applied pressure on the relevant physical parameters of Ti0.5Ta0.5 alloy, including structural, mechanical, and electronic properties. Results reveal that the calculated data for equilibrium lattice parameter, elastic constants, elastic moduli are consistent with the research results of others for Ti0.5Ta0.5 alloy under 0 T  and 0 P  . 11 C and 12 C increase monotonously as applied pressure increases, but 44 C increases slowly and then decreases gradually. The applied pressure cannot destroy the mechanical stability of Ti0.5Ta0.5 alloy because the stability criterion is always satisfied in the range of 0-50GPa. Bulk modulus B gradually increases with increasing applied pressure, but E and G reach the maximum values at P = 25 GPa, suggesting that high pressure increases the resistance to volume deformation, and the resistances to elastic and shear deformation are maximized under P = 25 GPa in Ti0.5Ta0.5 alloy. B G indicates that the ductility of the alloy significantly increases as applied pressure increases further when P > 25 GPa. The Ti0.5Ta0.5 alloy is an isotropic material when the applied pressure is 28 GPa because the values of anisotropy factors   

001
A and   

001
A are equal to 1. The major atomic bonding in Ti0.5Ta0.5 alloy is characterized by a metallic bond, and high pressure can enhance the metallic bond. The TDOS curves and the isosurface contours of charge density indicate that the

Conclusions
In conclusion, DFT calculations are implemented to study the influences of applied pressure on the relevant physical parameters of Ti 0.5 Ta 0.5 alloy, including structural, mechanical, and electronic properties. Results reveal that the calculated data for equilibrium lattice parameter, elastic constants, elastic moduli are consistent with the research results of others for Ti 0.5 Ta 0.5 alloy under T = 0 and P = 0. C 11 and C 12 increase monotonously as applied pressure increases, but C 44 increases slowly and then decreases gradually. The applied pressure cannot destroy the mechanical stability of Ti 0.5 Ta 0.5 alloy because the stability criterion is always satisfied in the range of 0-50GPa. Bulk modulus B gradually increases with increasing applied pressure, but E and G reach the maximum values at P = 25 GPa, suggesting that high pressure increases the resistance to volume deformation, and the resistances to elastic and shear deformation are maximized under P = 25 GPa in Ti 0.5 Ta 0.5 alloy. B/G indicates that the ductility of the alloy significantly increases as applied pressure increases further when P > 25 GPa. The Ti 0.5 Ta 0.5 alloy is an isotropic material when the applied pressure is 28 GPa because the values of anisotropy factors A (100) [001] and A (110)[001] are equal to 1. The major atomic bonding in Ti 0.5 Ta 0.5 alloy is characterized by a metallic bond, and high pressure can enhance the metallic bond. The TDOS curves and the isosurface contours of charge density indicate that the structural phase transition may occur in the Ti 0.5 Ta 0.5 alloy if the applied pressure increases further. The present results are valuable for the application of Ti 0.5 Ta 0.5 alloy under high pressure in the future, such as the wide applications in aerospace, marine, chemical, and biological fields under high pressure.