Research Article  Open Access
The Effect of Indium Concentration on the Structure and Properties of Zirconium Based Intermetallics: FirstPrinciples Calculations
Abstract
The phase stability, mechanical, electronic, and thermodynamic properties of InZr 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 Inp and Zrd 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 bulk modulus and Young’s modulus. The calculated theoretical hardness of αIn_{3}Zr is higher than the other InZr compounds.
1. Introduction
Saitovitch et al. made a study on the nuclear probe Cd on In sites of intermetallic compounds of the InZr system [1]. However, there is still short of experimental information available on the zirconiumindium system. Experiments on the phase diagram of InZr system present the phase compositions of InZr 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 metalindium compounds by hightemperature direct synthesis calorimetry and provided some new thermochemical data for InZr binary systems [4].
To the best of our knowledge, there is no systematically relevant research on the structural stability and mechanical properties of the InZr 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 indiumrelated materials is still a challenge.
Recently, advances in computeraided design and computeraided manufacturing (CAD/CAM) technology are used in novel manufacturing routes [5]. Firstprinciples calculation is a reliable and precise approach for studying groundstate 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 InZr compounds have been analyzed and discussed in detail by performing the firstprinciples calculations. It may improve the knowledge of the InZr alloys and provide available theoretical basis for further optimization and design for both indium and zirconium alloys.
2. 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 interaction between ions and electrons was described by using the ultrasoft pseudopotential [7]. The generalized gradient approximation (GGA) with the PerdewBurkeErnzerhof (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 MonkhorstPack 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.

3. Results and Discussions
3.1. Structural Stability
In this section, we established the initial crystal structures of four InZr compounds based on the experimental crystallographic data. Then the lattice parameters and internal coordinates of these compounds were optimized with the measure of firstprinciple calculation. Figure 2 presents the optimized structural models of four InZr 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 InZr compounds. Enthalpy of formation of InZr compounds is defined as follows:where is the total energy of and and 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 InZr 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 ZrIn compounds.
3.2. Mechanical Properties
Before considering the mechanical properties, we prefer to study the mechanical stability of these InZr 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, , , and . Counterpart criterion of mechanical stability is based on the following formula [14]:
(a) InZr3
(b) InZr
(c) In2Zr
(d) αIn3Zr
(a) InZr3
(b) InZr
(c) In2Zr
(d) αIn3Zr
For tetragonal system, there are six independent elastic constants, that is, , , , , , and . The qualification for stability can be measured in the following expression [14]:
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 InZr compounds conform to the corresponding standard of stability. That is, the InZr compounds considered in this paper are all mechanically stable. The mechanical prospects for all compounds show that four ZrIn intermetallic compounds have stable structure [16].
For researching the mechanical properties of the InZr 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 [18]
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 InZr 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 InZr system, indicating that shear modulus limits the mechanical stability of the InZr 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.509 GPa, which is higher than that of the other InZr 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 and are much lower than 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 brittleductile 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 InZr intermetallic compounds are ductile. Poisson’s ratio can be used to forecast the ductility of crystals, as well [25]. From Table 4, it is shown that B/G and have similar variation trend with increasing Zr content. It is found that when B/G is higher than 1.75, the values of 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 .
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 InZr binary compounds.
3.3. Electronic Properties
In order to study the electronic properties of the InZr intermetallic compounds, we calculate the electronic energy band structures and the results are depicted in Figure 2, where the dot dash line zeropoint 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 InZr intermetallics are conductive phases.
To further understand the electronic structure and mechanical properties of four InZr compounds, density of states was calculated and analyzed in detail. Figure 3 represents the total and partial density of states of four InZr 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 Inp and Zrd states. It indicates that the 5p states of In hybridize strongly with the Zr 4d states and the InZr atomic bonds formed along the pd 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].
(a) InZr3
(b) InZr
(c) In2Zr
(d) αIn3Zr
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 for InZr among these binary InZr 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.
3.4. 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, , , and 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.9 K 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.
4. Conclusions
In summary, we have performed firstprinciple calculations to investigate the structural, mechanical, electronic, and thermodynamic properties of InZr 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.509 GPa, respectively, which are higher than the other compounds. All InZr 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 InZr compounds.



Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This research work is supported by the National Natural Science Foundation of China (51361002 and 51161002), the Program for New Century Excellent Talents in University of China (NCET120650), and the Training Plan of HighLevel Talents of Guangxi University (2015).
References
 H. Saitovitch, P. R. J. Silva, J. T. Cavalcante, and M. Forker, “Zirconium–indium intermetallic compounds investigated by measurements of nuclear electric quadrupole interactions,” Journal of Alloys and Compounds, vol. 505, no. 1, pp. 157–162, 2010. View at: Publisher Site  Google Scholar
 H. Okamoto, “The InZr (IndiumZirconium) system,” Bulletin of Alloy Phase Diagrams, vol. 11, no. 2, pp. 150–152, 1990. View at: Publisher Site  Google Scholar
 M. F. Zumdick, G. A. Landrum, R. Dronskowski, R.D. Hoffmann, and R. Pöttgen, “Structure, chemical bonding, and properties of ZrIn2, IrIn2, and Ti3Rh2In3,” Journal of Solid State Chemistry, vol. 150, no. 1, pp. 19–30, 2000. View at: Publisher Site  Google Scholar
 S. V. Meschel and O. J. Kleppa, “Standard enthalpies of formation of some transition metal indium compounds by high temperature direct synthesis calorimetry,” Journal of Alloys and Compounds, vol. 333, no. 12, pp. 91–98, 2002. View at: Publisher Site  Google Scholar
 K. Yamanaka, M. Mori, and A. Chiba, “Developing high strength and ductility in biomedical CoCr cast alloys by simultaneous doping with nitrogen and carbon,” Acta Biomaterialia, vol. 31, pp. 435–447, 2016. View at: Publisher Site  Google Scholar
 M. D. Segall, P. J. D. Lindan, M. J. Probert et al., “Firstprinciples simulation: ideas, illustrations and the CASTEP code,” Journal of Physics: Condensed Matter, vol. 14, no. 11, pp. 2717–2744, 2002. View at: Publisher Site  Google Scholar
 D. Vanderbilt, “Soft selfconsistent pseudopotentials in a generalized eigenvalue formalism,” Physical Review B, vol. 41, no. 11, pp. 7892–7895, 1990. View at: Publisher Site  Google Scholar
 J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple,” Physical Review Letters, vol. 77, no. 18, pp. 3865–3868, 1996. View at: Publisher Site  Google Scholar
 K. Anderko, “Beitrag zu den binären Systemen des Titans mit Gallium, Indium und Germanium und des Zirkons mit Gallium und Indium,” Zeitschrift für Metallkunde, vol. 49, no. 4, pp. 165–172, 1958. View at: Google Scholar
 J. C. Uy and A. A. Burr, “The solute metallic valence as an index of phase stabilization in zirconiumbased alloys,” Transactions of the Metallurgical Society of AIME, vol. 224, no. 2, p. 204, 1962. View at: Google Scholar
 E. C. Moshopoulou, R. M. Ibberson, J. L. Sarrao, J. D. Thompson, and Z. Fisk, “Structure of Ce_{2}RhIn_{8}: an example of complementary use of highresolution neutron powder diffraction and reciprocalspace mapping to study complex materials,” Acta Crystallographica Section B: Structural Science, vol. 62, no. 2, pp. 173–189, 2006. View at: Publisher Site  Google Scholar
 B. D. Lichter, “Precision lattice parameter determination of zirconiumoxygen solid solution,” Transactions of the Metallurgical Society of AIME, vol. 218, 1960. View at: Google Scholar
 Q.J. Liu, N.C. Zhang, F.S. Liu, and Z.T. Liu, “Structural, mechanical and electronic properties of OsTM and TMOs_{2} (TM = Ti, Zr and Hf): firstprinciples calculations,” Journal of Alloys and Compounds, vol. 589, pp. 278–282, 2014. View at: Publisher Site  Google Scholar
 Z.J. Wu, E.J. Zhao, H.P. Xiang, X.F. Hao, X.J. Liu, and J. Meng, “Crystal structures and elastic properties of superhard IrN_{2} and IrN_{3} from first principles,” Physical Review B, vol. 76, no. 5, Article ID 054115, 15 pages, 2007. View at: Publisher Site  Google Scholar
 Q.J. Liu, H. Tian, and Z.T. Liu, “Mechanical properties and electronic structures of the HfSi system: firstprinciples calculations,” Solid State Communications, vol. 205, pp. 39–45, 2015. View at: Publisher Site  Google Scholar
 Z. Zhou, X. Zhou, and K. Zhang, “Phase stability, electronic structure and mechanical properties of IrBx (x = 0.9,1.1): firstprinciples calculations,” Computational Materials Science, vol. 113, pp. 98–103, 2016. View at: Publisher Site  Google Scholar
 L.H. Li, W.L. Wang, L. Hu, and B.B. Wei, “Firstprinciple calculations of structural, elastic and thermodynamic properties of FeB compounds,” Intermetallics, vol. 46, pp. 211–221, 2014. View at: Publisher Site  Google Scholar
 A. Kumar, A. Chernatynskiy, M. Hong, S. R. Phillpot, and S. B. Sinnott, “An ab initio investigation of the effect of alloying elements on the elastic properties and magnetic behavior of Ni_{3}Al,” Computational Materials Science, vol. 101, pp. 39–46, 2015. View at: Publisher Site  Google Scholar
 Y. Pan, M. Wen, L. Wang, X. Wang, Y. H. Lin, and W. M. Guan, “Iridium concentration driving the mechanical properties of iridiumaluminum compounds,” Journal of Alloys and Compounds, vol. 648, Article ID 34748, pp. 771–777, 2015. View at: Publisher Site  Google Scholar
 B. Huang, Y.H. Duan, W.C. Hu, Y. Sun, and S. Chen, “Structural, anisotropic elastic and thermal properties of MB (M=Ti, Zr and Hf) monoborides,” Ceramics International, vol. 41, no. 5, pp. 6831–6843, 2015. View at: Publisher Site  Google Scholar
 V. V. Bannikov, I. R. Shein, and A. L. Ivanovskii, “Structural, elastic, electronic properties and stability trends of 1111like silicide arsenides and germanide arsenides MCuXAs (M = Ti, Zr, Hf; X = Si, Ge) from first principles,” Journal of Alloys and Compounds, vol. 533, pp. 71–78, 2012. View at: Publisher Site  Google Scholar
 H. Hu, X. Wu, R. Wang, W. Li, and Q. Liu, “Phase stability, mechanical properties and electronic structure of TiAl alloying with W, Mo, Sc and Yb: firstprinciples study,” Journal of Alloys and Compounds, vol. 658, pp. 689–696, 2016. View at: Publisher Site  Google Scholar
 G. Yi, X. Zhang, J. Qin et al., “Mechanical, electronic and thermal properties of Cu_{5}Zr and Cu_{5}Hf by firstprinciples calculations,” Journal of Alloys and Compounds, vol. 640, pp. 455–461, 2015. View at: Publisher Site  Google Scholar
 E. Jain, G. Pagare, S. S. Chouhan, and S. P. Sanyal, “Electronic structure, phase stability and elastic properties of ruthenium based four intermetallic compounds: Abinitio study,” Intermetallics, vol. 54, pp. 79–85, 2014. View at: Publisher Site  Google Scholar
 S. Liu, Y. Zhan, J. Wu, and X. Wei, “Insight into structural, mechanical, electronic and thermodynamic properties of intermetallic phases in Zr–Sn system from firstprinciples calculations,” Journal of Physics and Chemistry of Solids, vol. 86, pp. 177–185, 2015. View at: Publisher Site  Google Scholar
 Y. Zhao, H. Hou, Y. Zhao, and P. Han, “Firstprinciples study of the nickelsilicon binary compounds under pressure,” Journal of Alloys and Compounds, vol. 640, pp. 233–239, 2015. View at: Publisher Site  Google Scholar
 Y. Pan, Y. Lin, X. Wang et al., “Structural stability and mechanical properties of Pt–Zr alloys from firstprinciples,” Journal of Alloys and Compounds, vol. 643, pp. 49–55, 2015. View at: Publisher Site  Google Scholar
 A. Lee, G. Etherington, and C. N. J. Wagner, “Partial structure functions of amorphous Ni_{35}Zr_{65},” Journal of NonCrystalline Solids, vol. 6162, no. 1, pp. 349–354, 1984. View at: Publisher Site  Google Scholar
 A. W. Weeber and H. Bakker, “Extension of the glassforming range of Ni–Zr by mechanical alloying,” Journal of Physics F: Metal Physics, vol. 18, no. 7, pp. 1359–1369, 1988. View at: Publisher Site  Google Scholar
 J. Du, B. Wen, R. Melnik, and Y. Kawazoe, “Firstprinciples studies on structural, mechanical, thermodynamic and electronic properties of Ni–Zr intermetallic compounds,” Intermetallics, vol. 54, pp. 110–119, 2014. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Fuda Guo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.