The Effect of Indium Concentration on the Structure and Properties of Zirconium Based Intermetallics : First-Principles Calculations

The phase stability, mechanical, electronic, and thermodynamic properties of In-Zr compounds have been explored using the firstprinciples calculation based on density functional theory (DFT). The calculated formation enthalpies show that these compounds are all thermodynamically stable. Information on electronic structure indicates that they possess metallic characteristics and there is a common hybridization between In-p and Zr-d states near the Fermi level. Elastic properties have been taken into consideration. The calculated results on the ratio of the bulk to shear modulus (B/G) validate that InZr 3 has the strongest deformation resistance. The increase of indium content results in the breakout of a linear decrease of the bulkmodulus and Young’s modulus.The calculated theoretical hardness of α-In 3 Zr is higher than the other In-Zr compounds.


Saitovitch et al. made a study on the nuclear probe Cd on
In sites of intermetallic compounds of the In-Zr system [1].However, there is still short of experimental information available on the zirconium-indium system.Experiments on the phase diagram of In-Zr system present the phase compositions of In-Zr compounds [2].Zumdick et al. examined the crystal structures and properties of ZrIn 2 and discussed the chemical bonding of ZrIn 2 on the basis of density functional calculations [3].Meschel and Kleppa studied the standard enthalpies of formation of some transition metal-indium compounds by high-temperature direct synthesis calorimetry and provided some new thermochemical data for In-Zr binary systems [4].
To the best of our knowledge, there is no systematically relevant research on the structural stability and mechanical properties of the In-Zr compounds.Indium possesses unique physical and chemical properties and has been widely applied in the fields of medicine and health, solar battery, national defense and military, aerospace, and modern information industry.The applications of indium compounds and alloys are restricted due to technical level and high cost.Actually, fundamental investigation on indium-related materials is still a challenge.
Recently, advances in computer-aided design and computer-aided manufacturing (CAD/CAM) technology are used in novel manufacturing routes [5].First-principles calculation is a reliable and precise approach for studying ground-state properties such as the structure, elastic properties, and phase stability.In this paper, we made a systemic investigation on the structure and properties of zirconium based intermetallics.Structural, mechanical, electronic, and thermodynamic properties of four In-Zr compounds have been analyzed and discussed in detail by performing the firstprinciples calculations.It may improve the knowledge of the In-Zr alloys and provide available theoretical basis for further optimization and design for both indium and zirconium alloys.

Computational Methods
All calculations were performed by first principle based on density functional theory (DFT) implemented in Cambridge sequential total energy package (CASTEP) code [6].The Advances in Condensed Matter Physics interaction between ions and electrons was described by using the ultrasoft pseudopotential [7].The generalized gradient approximation (GGA) with the Perdew-Burke-Ernzerhof (PBE) [8] was used for the exchange correlation functional.The electronic configurations in this study were In 4d 10 5s 2 5p 1 and Zr 4s 2 4p 6 4d 2 5s 2 , respectively.We choose various  points here mainly because the different phases possess diverse structure.All parameters including  point convergence test have been carefully tested.To ensure the precision of the calculations, we select firstly the  point and cutoff energy at the fine precision as the initial value and then improve continually precision by change the  point and cutoff energy until the calculated crystal parameters have a great agreement with previous data.The corresponding results were displayed in Table 1.According to the calculation at the different  points and cutoff energy, the cutoff energy and  point were selected for different phase.In this work, the cutoff energy of 450 eV was set to InZr and -In 3 Zr and 350 eV was set to In 2 Zr.Meanwhile, the cutoff energy of 550 eV was selected for all other phases.Brillouin zone integrations were performed by using a Monkhorst-Pack  point mesh and grids of 12 × 12 × 12, 7 × 7 × 7, 6 × 6 × 1, 6 × 6 × 2, 15 × 15 × 10, and 15 × 15 × 8  point were set to be the sampling of InZr 3 , InZr, In 2 Zr, -In 3 Zr, In, and -Zr, respectively.During the structural optimization, all atoms were fully relaxed using conjugate gradient method until the total energy, force, and atomic displacement were less than 1 × 10 −5 eV/ Å, 0.03 eV/ Å, and 0.001 Å, respectively.

Structural Stability.
In this section, we established the initial crystal structures of four In-Zr compounds based on the experimental crystallographic data.Then the lattice parameters and internal coordinates of these compounds were optimized with the measure of first-principle calculation.Figure 2 presents the optimized structural models of four In-Zr compounds and Table 2 lists the lattice parameter compared with the available experimental data.The reliability of this work and the reasonability of the computational methodology can be justified from the good agreement between the available experimental data and present work in optimized lattice parameters and atomic coordinates.
The enthalpy of formation was calculated in this part so as to estimate the thermodynamic stability of In-Zr compounds.Enthalpy of formation of In-Zr compounds is defined as follows: where  total is the total energy of In  Zr  and  bulk In and  bulk Zr are the total energy of an indium atom and a zirconium atom in the bulk state, respectively.The calculated results are listed in Table 2.The obtained formation enthalpies of the In-Zr compounds are all negative, meaning that they are all stable at 0 K and 0 GPa [13].We can see that the formation enthalpy of InZr is the lowest contrast with -In 3 Zr, In 2 Zr, and InZr 3 , which implies that InZr is the most stable among four Zr-In compounds.

Mechanical Properties.
Before considering the mechanical properties, we prefer to study the mechanical stability of these In-Zr compounds.According to the structural feature (see Figure 1), the compounds discussed in this paper mainly consist of two crystal structures such as cubic structure (i.e., InZr 3 and InZr) and tetragonal structure (i.e., In 2 Zr and -In 3 Zr).For the cubic crystal, there are three independent elastic constants, that is,  11 ,  12 , and  44 .Counterpart criterion of mechanical stability is based on the following formula [14]: (2) For tetragonal system, there are six independent elastic constants, that is,  11 ,  12 ,  13 ,  33 ,  44 , and  66 .The  ( Elastic constants reflect the resistance of a crystal to an applied external stress, whose values will provide valuable information about structural stability [15].Based on the calculation results, all the In-Zr compounds conform to the corresponding standard of stability.That is, the In-Zr compounds considered in this paper are all mechanically stable.The mechanical prospects for all compounds show that four Zr-In intermetallic compounds have stable structure [16].
For researching the mechanical properties of the In-Zr compounds, the following computational formulas are used [17]: For cubic system [17], Particularly, the InZr compound breaks the symmetry of the cubic system, and the general elastic constants are [ For tetragonal system [14], From Table 4, it can be seen that InZr 3 , with the bulk modulus of 97.112 GPa, has the strongest resistance to volume change by applied pressure [19] and the strongest average bond strength of atoms for the given crystal [20].On the contrary, -In 3 Zr is the weakest one compared with the other In-Zr compounds.As shown in Table 4 the relation between bulk modulus and indium content is monotonically decreasing, which is equivalent to the resistance to volume deformation of the four compounds and the abrupt degradation of average bond strength of atoms with increasing indium content.The shape change capacity of material is usually reflected in the value of shear modulus [15].The difference of this value between InZr 3 and -In 3 Zr is quite small.The lowest shear modulus is 47.007 GPa for InZr and the highest shear modulus is only 51.6 GPa for InZr 3 .There is no obvious rule about shear modulus with increasing indium content.What is more, it can be seen that  >  exists in the In-Zr system, indicating that shear modulus limits the mechanical stability of the In-Zr compounds [21].Materials with high Young's modulus depict strong resistance to uniaxial tension [22] and indicate the elastic stiffness as well.It can be seen from Table 4 that calculated Young' modulus of InZr 3 is 131.509GPa, which is higher than that of the other In-Zr compounds.Young's modulus of -In 3 Zr is only about 118.350 GPa, which is lower than all the others.We can conclude from above analysis that  has similar change compared with .Young's modulus decreases linearly with increasing indium content.In addition, it should be noted that  12 and  44 are much lower than  11 for all compounds, which implies that all the crystals are more difficult to deform in the axial direction, compared with shear direction [23].
B/G ratio, in general, was applied to describe the brittleness and toughness of materials [24].The brittle-ductile transition of materials can be judged by the B/G ratio.If the B/G ratio exceeds 1.75, the material is ductile.Otherwise, it is brittle.From Table 4, the B/G value of -In 3 Zr is less than 1.75, which implies the brittleness of this compound.The B/G values of InZr 3 , InZr, and In 2 Zr are all higher than 1.75, demonstrating that these In-Zr intermetallic compounds are ductile.Poisson's ratio V can be used to forecast the ductility of crystals, as well [25].From Table 4, it is shown that B/G and V have similar variation trend with increasing Zr content.It is found that when B/G is higher than 1.75, the values of V are more than 0.26 for InZr 3 , InZr, and In 2 Zr.This is consistent with the conclusion that ductile compounds usually have high V.
The hardness parameter  can be calculated [26]: The calculated theoretical hardness was shown in Table 4.We can know that there is a transition at 0.67 indium content.When indium content is less than 0.67, the hardness may decrease with increasing indium content.While it is higher than 0.67, hardness exhibits a big jump and increases with the indium increasing.The calculated theoretical hardness of -In 3 Zr is 15.809, which indicates that -In 3 Zr is expected to be the hardest among all the In-Zr binary compounds.

Electronic Properties.
In order to study the electronic properties of the In-Zr intermetallic compounds, we calculate the electronic energy band structures and the results are depicted in Figure 2, where the dot dash line zero-point energy signifies the highest energy level occupied by the valence electrons at 0 K. From Figure 2, we can know that all compounds show metallic nature due to the fact that these bands cross   .Meanwhile, the density of states corresponding to the Fermi surface was all positive, which also implies that all these In-Zr intermetallics are conductive phases.
To further understand the electronic structure and mechanical properties of four In-Zr compounds, density of states was calculated and analyzed in detail.Figure 3 represents the total and partial density of states of four In-Zr compounds, where the black vertical dashed of DOS indicates the Fermi level (  ).As is vividly demonstrated in Figure 3, the valence electron near   mainly derived from In-p and Zr-d states.It indicates that the 5p states of In hybridize strongly with the Zr 4d states and the In-Zr atomic bonds formed along the p-d directions near the region Fermi level.Particularly, the DOS profile of InZr 3 presents an obvious pseudogap.The pseudogap indicates that the hybridization between In and Zr induces a separation of the bonding states and antibonding states [13].That is why InZr 3 own the greatest deformation resistance [27].
The distance between the zero energy symbolized bands and the nearest peaks may manifest the bonds transformation of compounds [28,29].In general, the longer the distance to the zero energy symbolized bands, the stronger the covalent bonds.Moreover, with the falling in metallicity of bond, alloy ductility usually decreases [30].As the DOS shown in Figure 3, the distance adjacent zero energy is the lowest Advances in Condensed Matter Physics 7 for InZr among these binary In-Zr phases, implying that it possesses the strongest metallicity and hence it has the best ductility.This result has a well consistency with the ductility estimated from the / and Poisson's ratio values.
Furthermore, we can estimate the structural stability of intermetallic compounds through different density of states at the Fermi levels [16].Usually, the lower the value at Fermi levels, the more strong the stability of the structure.From Figure 3, InZr has the minimum value of density of states and in consequence has the most stable structure, which has a great consistency to our results of formation enthalpy.

Thermodynamic Properties.
Debye temperature was calculated to provide insight into the thermodynamics behaviors.The values of Young's modulus, bulk modulus, and shear modulus were obtained from the above section and then put into the following formula [16]: where Θ  is Debye temperature, ℎ is Planck's constant,  is Boltzmann's constant,  is the number of atoms per formula,  is the molecular weight, V  , V  , and V  are the average wave velocity, the longitudinal, and shear velocity, respectively, and  is the theoretical density of the compound. and  are the bulk modulus and shear modulus, respectively.The values of Debye temperature, average wave velocity, longitudinal, and shear wave velocity for each compound are listed in Table 5.
As we all know the strength of chemical bonding can be characterized by Debye temperature which is inverse to the molecular weight [16,23].Stronger chemical bonding favors higher Debye temperature.As shown in Table 5, we can find that the highest Θ  value is 314.9K for InZr 3 , and the lowest one is ascribed to -In 3 Zr.Therefore, we can conclude that the covalent bonds in InZr 3 are stronger than the others.

Conclusions
In summary, we have performed first-principle calculations to investigate the structural, mechanical, electronic, and thermodynamic properties of In-Zr compounds with different indium concentrations.The calculated results show that the lattice parameters of all compounds are in good agreement with the obtained experimental data and theoretical results.The calculated formation enthalpy indicates that InZr is more stable than the others.The calculated elastic constants show that all the compounds are mechanically stable (Table 3).The bulk modulus and Young's modulus decrease linearly with increasing indium content.The calculated shear modulus and Young's modulus of InZr 3 are 51.6 GPa and 131.509GPa, respectively, which are higher than the other compounds.All In-Zr compounds exhibit ductile behavior but -In 3 Zr exhibits brittle behavior.We conclude that the variation of mechanical properties is related not only to the indium concentration, but also to the bonding state in the In-Zr compounds.

Figure 3 :
Figure 3: Total and partial electronic densities of states near Fermi level of In-Zr intermetallics.The dotted lines denote the Fermi level   .

Table 1 :
The lattice constants ( Å) of In-Zr compounds at different  points and cutoff energy.

Table 2 :
Calculated lattice parameters ( Å) and formation enthalpies (eV/atom) for In, -Zr, and In-Zr intermetallic compounds in comparison with the available experimental data.

Table 5 :
Sound velocity, density, and Debye temperature of In-Zr compounds.