First-principles investigations on structural stability, mechanical, and thermodynamic properties of LaT2Al20 (T = Ti, V, Cr, Nb, and Ta) intermetallic cage compounds
Quan Shanyu1, Zhang Xudong1, †, Liu Cong1, Jiang Wei1, 2
School of Science, Shenyang University of Technology, Shenyang 110870, China
School of Materials Science and Engineering, Shenyang University of Technology, Shenyang 110870, China

 

† Corresponding author. E-mail: zxdwfc@163.com

Project supported by the Program for Ph. D Start-up Fund of Liaoning Province of China (Grant No. 201601161).

Abstract

First principles calculations were used to explore the structural stability, mechanical properties, and thermodynamic properties of LaT2Al20 (T = Ti, V, Cr, Nb, and Ta) intermetallics. The calculated formation enthalpy and phonon frequencies indicate that LaT2Al20 intermetallics exhibit the structural stability. The elastic moduli (B, G, E, and Hv) indicate that these intermetallics possess the better elastic properties than pure Al. The values of Poisson’s ratio v and B/G demonstrate that LaT2Al20 intermetallics are all brittle materials. The anisotropy of elasticity and Young’s modulus (three- and two-dimensional figures) indicate that LaT2Al20 compounds are anisotropic. Importantly, the calculated thermal quantities demonstrate that LaT2Al20 intermetallics possess the better thermal physical properties than pure Al at high temperatures.

1. Introduction

Aluminum alloys have been playing an important role in many industrial fields due to their outstanding physical properties.[1] But the poor creep-resistance properties of aluminum alloys have hindered their wide applications at high temperatures. In the past few decades, many researchers found that the addition of appropriate transition metal elements (TM) and rare earth elements (RE) to aluminum alloys can form L12, D022, D023 or D019-ordered binary Al3TM and Al3RE intermetallics. These binary intermetallics exhibit some outstanding physical properties, including low mass densities, high melting points, and good oxidation resistances, and so on, which can effectively improve the poor creep-resistance of aluminum alloys.[215] In the meantime, some ternary AT2Al20 intermetallics systems (where A = actinide/lanthanide/calcium and T = transition metal) were synthesized with CeCr2Al20-type structure.[1632] Kangas and Thiede synthesized the new LaT2Al20 (T = Ti, V, Cr, Nb, and Ta) intermetallics by using a step-wise arc melting technique and by self-flux method.[31,32] For these new LaT2Al20 phases, although their intrinsic physical properties have been explored, their mechanical properties and thermodynamics properties are not clarified. It is necessary to investigate the related physical properties for further applications. In the present work, the first-principles calculations were used to comparatively investigate the structural stability, mechanical, and thermodynamics properties for Al and LaT2Al20 intermetallics.

2. Computational methods

The first-principles calculations were performed by using the VASP code based on the density functional theory (DFT). The electron–ion interaction and the exchange–correlation energy were described under the projector-augmented wave (PAW) method and the generalized gradient approximation (GGA) of Perdewe–Burkee–Ernzerhof (PBE) scheme, respectively.[33] The f electrons of La atom were frozen into the core states in PAW potentials. The energy cut off for the plane wave basis is 600 eV for all compounds. The 9 × 9 × 9 k-points for LaT2Al20 unit cell with the Monkhorst–Pack scheme were constructed for the Brillouin zone sampling. Spin polarization was used in the present calculations. The phonon dispersions, density of states, and thermodynamics quantities were calculated by using the supercell method and linear response method with PHONOPY code.[34] The stress–strain method was used to calculate the elastic constants of LaT2Al20 intermetallics.[35]

3. Results and discussion
3.1. Structural stability

LaT2Al20 intermetallics can crystallize in the face-centered cubic structure with Fd-3m space group. In the prototypic LaT2Al20 compounds, La atoms occupy the 8a Wyckoff sites, T atoms occupy the 16d Wyckoff sites, and Al atoms occupy the 16c, 48f, and 96g Wyckoff sites. The unit cell of structure is shown in Fig. 1.[25] The optimized atom coordinates are listed in Table 1, together with experimental ones. In the meantime, the ratios of volume change (V/V0) under different pressures were shown in Fig. 2. In Table 1, the calculated data are slightly larger than experimental ones[31,32] and the error is less than 5%, which indicates that the computational method is reliable. LaNb2Al20 has the largest bulk modulus B0 among five intermetallics, which indicates the LaNb2Al20 has the best ability to resist the compression. In Fig. 2, the curves of V/V0 are smooth and have no inflection points until 50 GPa, which indicates that these intermetallics can keep structural stability under high pressures. The curves of LaNb2Al20 and LaV2Al20 locate at upper part and lower part, which indicates that LaNb2Al20 and LaV2Al20 have the best/worst resistance ability to compression among five compounds.

Fig. 1. (color online) Crystal structure of LaTAl20 intermetallics with CeCr2Al20-type structure. The violet balls are La atoms, the red balls are T atoms, and the blue balls are Al atoms.
Fig. 2. (color online) Equation of state for LaT2Al20 intermetallics.

Generally speaking, the structural stability of the solid materials mainly depends on two conditions: one is thermodynamic condition, which depends on the formation enthalpy; the other is the dynamical condition, which depends on the phonon frequencies.[36,37] The formation enthalpy was calculated using the following expression[3841]

where ETotal represents the total energy of LaT2Al20; x, y, and z are the numbers of atoms in the unit cell; ELa, ET, and EAl represent the total energy of per atom in solids. The calculated formation enthalpies ΔH for LaT2Al20 are all negative values in Table 1, indicating that LaT2Al20 compounds have the stronger thermodynamic stability and is more inclined to form LaTa2Al20 due to the lower value of ΔH.

Table 1.

The lattice constant a0, bulk modulus B0 and formation enthalpy. ΔH of LaT2Al20 (T = Ti, V, Cr, Nb, and Ta) intermetallics.

.

In order to explore the dynamical stability, the phonon dispersion and the density of states are calculated in Figs. 3 and 4. The solid lines are calculated values and the scatter points are experimental ones of the inelastic neutron scattering measurement in Fig. 3.[42] The calculated values are close to the experimental ones, which indicate that the calculated results are reasonable. Figure 4 shows the phonon dispersion curves and the density of states for LaT2Al20 intermetallics. For five figures of phonon dispersion curves, the imaginary phonon frequencies are absent along the Brillouin zone path. It indicates that the five compounds possess the dynamical stability. Figure 5 displays the total phonon density of states for LaT2Al20 intermetallics. For the phonon density of states (PHDOS), the imaginary frequencies are also absent. It indicates that the five compounds possess the dynamical stability. In summary, according to the obtained formation enthalpies and phonon frequencies, the five intermetallics exhibit the structural stability.

Fig. 3. Phonon dispersion curves of Al and the circles correspond to the experimental data.
Fig. 4. (color online) Phonon dispersion curves and density of states for LaT2Al20 intermetallics.
Fig. 5. (color online) The total phonon density of states for LaT2Al20 intermetallics.
3.2. Elastic properties

The elastic constants and elastic moduli are important parameters for engineering applications. Table 2 listed the calculated elastic constants Cij and compliance matrices Sij, together with the experimental values of pure Al for comparison.[4345] The elastic constants of Al are consistent with experimental ones (as shown in Table 2). The calculated results indicate that our computational method is reliable. For the cubic structure intermetallics, there are three independent constants C11, C12, and C44. The elastic constants are satisfied with the mechanically stable conditions:[46] C11 ˃ 0, C44 ˃ 0, C11 ˃ |C12|, (C11 + 2C12)˃ 0. After obtaining the elastic constants Cij and Sij, the bulk modulus B, shear modulus G, Young’s modulus E, and hardness Hv were calculated by using the VRH method with the following equations[46]

The calculated elastic moduli of LaT2Al20 are listed in Tables 3 and 4, together with the experimental ones of Al.[43,45] The bigger values of bulk modulus B and shear modulus G correspond to the better resistance of solid materials against volume change under hydrostatic pressure and shape deformation by the shearing force. In Table 3, the values of elastic modulus B follow the order: LaNb2Al20 ˃ LaTa2Al20 ˃ LaCr2Al20 ˃ Al ˃ LaTi2Al20 ˃ LaV2Al20, this indicates that LaNb2Al20 has the best ability of resistance to the volume deformation. The values of elastic modulus G follow the order: LaNb2Al20 ˃ LaTa2Al20 ˃ LaCr2Al20 ˃ Al, indicating that LaNb2Al20 has the best ability of resistance to shape change under external conditions. Five LaT2Al20 intermetallics have the stronger shear deformation resistance in comparison with the pure Al. The larger elastic constants C44, elastic modulus G, and E correspond to the higher hardness of solid materials. The values of C44, G, and E of five intermetallics are remarkably larger than the ones of pure Al, indicating that these ternary intermetallics possess the higher hardness than the pure Al. In the meantime, for LaTi2Al20, LaV2Al20, LaCr2Al20, LaNb2Al20, LaTa2Al20, and Al, the values of hardness are also calculated by using the semi-empirical equations with Eq. (4): 12.249 GPa, 12.191 GPa, 7.747 GPa, 12.590 GPa, 9.359 GPa, and 2.868 GPa. The calculated results indicate the LaT2Al20 compounds are harder than Al.

Table 2.

Calculated elastic constants Cij (in unit GPa) and compliance matrices Sij.

.
Table 3.

Calculated elastic moduli (in unit GPa) B, G, E, B/G, Hv, and Poisson’s ratio v.

.
Table 4.

Calculated bulk modulus BV, BR (in unit GPa), shear modulus GV, GR (in unit GPa), and elastic anisotropic factors.

.

The ductile or brittle behavior influences the high temperature application of solid materials. The values of B/G ratio and Poisson’s ratio v can determine the brittle and ductile behaviors of solid materials. When the values of B/G and v are less than 1.75 and 0.26, the solid materials possess brittle behavior. Otherwise, the solid materials possess ductile behavior. The calculated values of B/G and v are shown in Table 3. The values of B/G and v of LaT2Al20 intermetallics are all smaller than the critical value 1.75 and 0.26, which indicates that LaT2Al20 intermetallics are all brittle materials. But Al possesses the ductile behavior because of the larger B/G and v values (2.774 and 0.339). In the meantime, LaV2Al20 has the smallest B/G and v, this indicates that it possesses the largest brittleness among these compounds. The values of elastic constants (C11C12) and Young’s modulus are the criterion to estimate the plasticity of solid materials. Generally speaking, if the solid materials have the smaller values of elastic constants and Young’s modulus, they will possess better plasticity.[47] In order to intuitively describe the plastic behavior, the results of the LaT2Al20 compounds and the experimental ones for pure Al[45] are shown in Fig. 6. In Fig. 6, the pure Al has the best plasticity and LaT2Al20 compounds have the worse plasticity than the pure Al. In order to present the brittleness of LaT2Al20 compounds, figure 7 shows the renormalized hyperbolic correlation derived by (C12C44)/E and G/B, including the calculated values of LaT2Al20 with some typical materials, Au, Pd, Al, Si, Ge, BN, diamond, and so on.[48] According to the Pugh criterion and Pettifor criterion, the LaT2Al20 compounds possess the brittle nature and high hardness, which is close to that of Si, Ge, and BN.

Fig. 6. (color online) The values of C11C12 and Young’s modulus for pure Al and LaT2Al20 intermetallics.
Fig. 7. (color online) The renormalized hyperbolic correlation derived by (C12C44)/E and G/B.
3.3. Anisotropic properties

The anisotropic properties are closely related to the micro-cracks behavior of solid materials. The anisotropic properties were described by the universal and shear anisotropic index AU and A[46]

The calculated anisotropic indexes AU and A are shown in Table 4. When the anisotropic values of AU and A approach to and 1 respectively, the crystalline phases possess the isotropic properties. The anisotropic data listed in Table 4 confirm that these compounds are anisotropic. The indexes AU and A correspond to the same order: LaTa2Al20 ˃ LaV2Al20 ˃ LaTi2Al20 ˃ Al ˃> LaCr2Al20 ˃ LaNb2Al20, which indicates that LaTa2Al20 has the largest anisotropy, but LaNb2Al20 possesses the smallest anisotropy. But the degree of anisotropy of LaT2Al20 compounds is small and LaNb2Al20 approaches to the isotropy.

Fig. 8. (color online) (a)–(f) The surface contours of the Young’s modulus for LaTi2Al20 (a), LaV2Al20 (b), LaCr2Al20 (c), LaNb2Al20 (d), LaTa2Al20 (e), and Al (f).

In order to visually describe anisotropic properties of LaT2Al20 intermetallics, the three-dimensional surfaces of Young’s modulus E was adopted to further explore the anisotropy of crystal phases. The three-dimensional (3D) figures were displayed using the following equation[46]

If the 3D surfaces of Young’s modulus are spherical, they will be isotropic. Otherwise, they are anisotropic. The greater the degree of deviation from the sphere is, the stronger the anisotropic properties are. The 3D figures of Young’s modulus are shown in Fig. 8. According to the degree of deviation from the sphere, it can be concluded that LaT2Al20 intermetallics are all anisotropic. The degree of deviation from the sphere for LaNb2Al20 is smallest among six figures, suggesting that it is basically isotropic. On the contrary, LaTa2Al20 has the largest anisotropy. The figures of other compounds display the degree of deviation, but the degree of deviation is small, indicating that these intermetallics are anisotropic and the anisotropy is small. In order to furthermore describe the anisotropic properties of the LaT2Al20 intermetallics, the two-dimensional figures of Young’s modulus on different planes are plot in Fig. 9. If the contours of the closed curves are circles, the crystal materials possess the isotropic properties. If the contours of the closed curves deviate from circles, the crystal materials possess the elastic anisotropy. In Fig. 9, the contours of the LaNb2Al20 are closed to circles, but LaTa2Al20 and Al possess the larger degree of deviation, indicating that LaNb2Al20 is more closed to the elastic isotropy.

Fig. 9. (color online) The projections on the (001) and (110) planes of Al and LaT2Al20 intermetallics.
3.4. Electronic and magnetic properties

The total density of states (TDOS) and partial density of states (PDOS) of LaT2Al20 from −12 eV to 4 eV are shown in Fig. 10. The dash line at zero energy represents the Fermi level EF in Fig. 10. There are no energy gaps near the Fermi level EF, indicating that five intermetallics exhibit metallic property. The region of energy binding states from about −12 eV to −4 eV in TDOS is determined mainly by Al 3s and 3p states. The energy binding states near the Fermi level EF are hybridized mainly by La 5d, T (Ti/V/Cr) 3d, and Al 3p states. The strong hybridization creates the pseudogap near the Fermi level in TDOS for five LaT2Al20 intermetallics. The strong hybridization can guarantee the electronic interaction between the atoms and improve the bond strength. LaT2Al20 intermetallics exhibit the better elastic properties than pure Al, which is attributed to the bond strength of La–Al, La–T, and T–Al metallic bonds of LaT2Al20 intermetallics. Comparing the up states with the down states of TDOS in Fig. 11, we found that the spin-down and spin-up states are not symmetric for LaCr2Al20 intermetallics, but the spin states of other four LaT2Al20 (T = Ti, V, Nb, and Ta) compounds are symmetric. It indicates that LaCr2Al20 intermetallics exhibit magnetic properties and other four LaT2Al20 compounds are non-magnetic compounds. The results are consistent with the experiments ones.[17,29,31]

Fig. 10. (color online) Density of states of LaT2Al20 intermetallics.
Fig. 11. (color online) Total density of states of spin-up and sipn-down for LaT2Al20 intermetallics.
Fig. 12. (color online) Relationships between Gibbs free energy and temperature of pure Al and LaT2Al20 intermetallics.
3.5. Thermodynamic properties

For cubic structural crystals, they have only one-independent lattice parameter. Thus, we performed the lattice dynamic calculations for the variable lattice parameter from 0.94a to 1.06a in a step of 0.01a within the PHONOPY code. By using F (V, T)–V data, we fit them to Vinet EOS and obtained the related thermodynamic parameters. The relationships between thermodynamic quantities and temperature for pure Al and LaT2Al20 intermetallics are shown in Figs. 1215. Because the melt point temperature of Al is about 960 K, so we gave out the thermodynamic quantities between 0 K and 900 K for pure Al. Figure 12 shows the relationships between Gibbs free energy and the elevated temperature. The Gibbs free energy is negative and decreases gradually with increasing temperature, indicating that these compounds are thermodynamically stable in the whole temperature range. The smaller the Gibbs free energy is, the better the thermal stability of compound is. Five LaT2Al20 intermetallics exhibit the better thermal stability than pure Al in the high temperature section. Figure 13 shows the relationships between bulk modulus and the elevated temperature. It is clear that when the temperature is less than 100 K, the bulk modulus nearly remains constant, but it drops remarkably in the high temperature section. The change of bulk modulus of LaT2Al20 intermetallics is smaller than Al, indicating that the ability of creep resistance of LaT2Al20 intermetallics are better than pure Al. Figure 14 shows the relationships between thermal expansion coefficients and temperature. It is noted that the expansion coefficients increase exponentially in the low temperature section and increase slowly (approaches to the linear increases) in the high temperature regime. But in the high temperature section, the values of pure Al are remarkably larger than those of five LaT2Al20 compounds at a given temperature, indicating that the pure Al becomes easier to expand at high temperature. Figure 15 shows the relationships between entropy and temperature. The values of entropy of pure Al and RET2Al20 compounds rapidly increase almost linearly with the elevated temperature. In order to prove the correctness of our calculation method, we compared the theoretical values and the experimental ones[4951] of bulk modulus and entropy of pure Al in Figs. 13 and 15. The calculated results are very close to the experimental ones, which indicate that the calculated method and results are reliable.

Fig. 13. (color online) Relationships between bulk modulus and temperature of LaT2Al20 intermetallics.
Fig. 14. (color online) Relationships between thermal expansion and temperature of pure Al and LaT2Al20 intermetallics.
Fig. 15. (color online) Relationships between entropy and temperature of pure Al and LaT2Al20.
4. Conclusions

In summary, the first-principles method is used to explore the structural stability, mechanical properties, and thermodynamic properties of pure Al and LaT2Al20 compounds. The calculated formation enthalpy and phonon frequencies indicate that LaT2Al20 intermetallics exhibit the structural stability. The elastic constants Cij and moduli (B, G, and E) indicate that LaT2Al20 intermetallics possess the better elastic properties than pure Al. The anisotropy in elastic constants and Young’s modulus indicates that LaT2Al20 compounds are anisotropic. The calculated thermal quantities indicate that LaT2Al20 intermetallics possess the better thermal physics properties at high temperatures.

Reference
[1] Schlapbach L Züttel A 2001 Nature 414 353
[2] Taend J Orthacker A Amenitsch H 2016 Acta Mater. 117 43
[3] Yu L B Wang J Qu F S Wang M Wang W 2018 J. Alloys Compd. 737 655
[4] Suwanpreecha C Pandee P Patakham U Limmaneevichitr C 2018 Mater. Sci. Eng. 709 46
[5] Dorin T Ramajayam M Lamb J Langanc T 2017 Mater. Sci. Eng. 707 58
[6] Tian T Wang X F Li W 2013 Solid State Commun. 156 69
[7] Marquis E A Seidman D N Dun D C 2002 Acta Mater. 50 4021
[8] Iwamura S Miura Y 2004 Acta Mater. 52 591
[9] Mikhaylovskay A V Mochugovskiy A G Levchenko V S Mufalo W 2018 Mater. Charact. 139 30
[10] Chen Z Zhang P Chen D Wu Y Wang M 2015 J. Appl. Phys. 117 085904
[11] Clouet E Barbu A La L Martin G 2005 Acta Mater. 53 2313
[12] Vo N Q Dun D C Seidman D N 2014 Acta Mater. 63 73
[13] Spierings A B Dawson K 2017 Mater. Design 115 52
[14] Lin J D Okle P Dun D C 2017 Mater. Sci. Eng. 680 64
[15] Wen S P Wang W Zhao W H 2016 J. Alloys Compd. 687 143
[16] Halevy I Sterer E Aizenshtein M Kimmel G Regev D Yahel E 2001 J. Alloys Compd. 319 19
[17] Yamada A Higashinaka R Matsuda T D Aoki Y 2018 J. Phys. Soc. Jpn. 87 033707
[18] Bram A I Venkert A Meshi L 2015 J. Alloys Compd. 641 1
[19] Wakiya K Onimaru T Matsumoto K T 2017 J. Phys. Soc. Jpn. 86 3
[20] Winiarski M J Griveau J C Colineau E Wochowski K 2017 J. Alloys Compd. 696 1113
[21] Liu P C Xian Y J Wang X Zhang Y T Zhang P C 2017 J. Nucl. Mater. 493 147
[22] Winiarski M J Klimczuk T 2017 J. Solid State Chem. 245 10
[23] Moussa C Berche A Pasturel M Barbosa J Stepnik B Dubois Tougait O 2017 J. Alloys Compd. 691 893
[24] Winiarski M J Klimczuk T 2017 Intermetallics 85 103
[25] Winiarski M J Wiendlocha B Sternik M Winiewski P O’Brien J R Kaczorowski D Klimczuk T 2016 Phys. Rev. 93 134507
[26] Sakai A Nakatsuji S 2011 Phys. Rev. 84 201106
[27] Swatek P Kaczorowski D 2016 J. Magn. Magn. Mater. 416 348
[28] Tsujimoto M Matsumoto Y Tomita T Sakai A Nakatsuji S 2014 Phys. Rev. Lett. 113 267001
[29] Kumar K R Nair H S Christian R Thamizhavel A Strydom A M 2016 J. Phys.: Condens. Matter 28 436002
[30] Verbovytsky Y Latka K Tomala K 2007 J. Alloys Compd. 442 334
[31] Kangas M J Schmitt D C Sakai A Nakatsuji S Chan J Y 2012 J. Solid State Chem. 196 274
[32] Thiede V M T Jeitschko W Niemann S Ebel T 1998 J. Alloys Compd. 267 23
[33] Kresse G Joubert D 1999 Phys. Rev. 59 1758
[34] Togo A Tanaka I 2015 Scr. Mater. 108 1
[35] Shang S Wang Y Liu Z K 2007 Appl. Phys. Lett. 90 101909
[36] Pan Y Guan W M 2018 Ceram. Int. 44 9893
[37] Pan Y Guan W M Li Y Q 2018 Phys. Chem. Chem. Phys. 20 15863
[38] Song Y Zuo X Yue Y L 2018 Chin. Phys. 27 037102
[39] Li Y Dai X F Liu G D Wei Z Y Liu E K Han X L Du Z W Xi X K Wang W H Wu G H 2018 Chin. Phys. 27 026101
[40] Liu D R Han D Huang M L Zhang X Zhang T Dai C M Chen S Y 2018 Chin. Phys. 27 018806
[41] He W Q Huang H B Liu Z H Ma X Q 2018 Chin. Phys. 27 016201
[42] Stedman R Nilson G 1966 Phys. Rev. 145 492
[43] Li Z Tse J S 2000 Phys. Rev. 61 14531
[44] Mao Z Chen W Seidman D N Wolverton C 2011 Acta Mater. 59 3012
[45] Grimvall G 1999 Thermo Physical Properties of Materials Amsterdam North Holland 256
[46] Chen S Sun Y Duan Y H Huang B Peng M J 2015 J. Alloys Compd. 630 202
[47] Hu W C Liu Y Li D J Zeng X Q Xu C S 2013 Physica 427 85
[48] Niu H Y Chen X Q Liu P T Xing W W Cheng X Y Li D Z 2012 Sci. Rep. 2 718
[49] Vinet P Rose J Ferrante J Smith J 1989 J. Phys.: Condens. Matter 1 1941
[50] Wallace D C 1972 Thermodynamics of Crystals New York Dover 241
[51] Barin I 1995 Thermochemical Data of Pure Substance Beijing Science Press 528