First-Principles Study of the Structural, Mechanical and Thermodynamic Properties of Al 11 RE 3 in Aluminum Alloys

: The stability and mechanical and thermodynamic properties of Al 11 RE 3 intermetallics (RE = Sc, Y and lanthanide La-Lu) have been investigated by combining ﬁrst-principles and Debye model calculations. It was found that the formation enthalpies of the Al 11 RE 3 intermetallics are all negative, indicating that they are stable; moreover, the experimental values of Al 11 La 3 and Al 11 Ce 3 are in good agreement with the predicted values, which are − 0.40 kJ/mol and − 0.38 kJ/mol, respectively. The calculated results of the mechanical properties reveal that the Young’s modulus E and shear modulus G of Al 11 RE 3 (RE = La, Ce, Pr, Nd and Sm) intermetallics are obviously greater than that of Al, implying that the stiffness, toughness, and tensile strength of them are signiﬁcantly greater than those of aluminum, and that they, as strengthen phases, can effectively improve the mechanical property of aluminum alloys. The Poisson’s ratio v of Al 11 Sc 3 (0.37) is the largest, and the heterogeneity is obvious. All the Al 11 RE 3 intermetallics can enhance the thermostability of the aluminum because of their lower Gibbs free energy F in the range of − 5.002~ − 4.137 eV/atom and thermal expansion coefﬁcient α of Al in the range of 2.34~2.89 × 10 − 5 /K at 300K, as well as higher entropy and constant volume-speciﬁc heat than aluminum at ﬁnite temperatures. With an increase in the atomic number, different change trends were observed for the formation enthalpy ∆ H f , bulk modulus B , Young’s modulus E , and shear modulus G . This paper can provide ideas and help for designing a high-performance, heat-resistant aluminum alloy.


Introduction
Aluminum (Al) alloys have become an important lightweight material for energy conservation and emission reduction in the automobile, rail transit, aerospace, and other industries due to their many excellent performances, including their light weight, high specific strength, good formability, excellent weldability, and corrosion resistance [1][2][3][4][5][6].As engineering materials, Al alloys have been widely used at room temperature; however, it is difficult to meet the requirements of high-temperature devices, such as aircraft and automobile engines, that have been used for a long time at more than 200 • C [7,8].The cause of this issue is the rapid coarsening and dissolution of the strengthen phases in the Al alloys at high temperature, which leads to a sharp decrease in the mechanical properties of the alloys [9,10].Therefore, improving the thermal stability of the strengthen phases at high temperatures is the key to expanding the application of Al alloys in hightemperature devices.
In recent years, alloying by adding rare-earth elements (RE) has become an important and effective method to improve the high-temperature performance of the Al alloys, mainly because of the formation of the high temperature-resistant strengthen phases and Al x RE y

Computation Detail
First-principles calculations based on density functional theory (DFT) [26] employed the plane-wave, total-energy methodology, which solved the Kohn-Sham equation to obtain the valence electron density and the wave function as implemented in the Vienna ab initio Simulation Package (VASP) [27].The interaction between the ionic core and valence electrons [28] was described by the Projector Augmented Wave (PAW) [29] method, and the exchange correlation function was calculated using the Perdew-Burke-Ernzerhof (PBE) version of the generalized gradient approximation (GGA) [30].To ensure that the results are accurate, the cutoff energy of the plane-wave basis was chosen as 350 eV, and Brillouin zone sampling was performed using the Gamma-centered Monkhorst-Pack method with 8 × 3 × 3 k-mesh grids [31].To obtain a better calculation and accurate results, the ground state geometry was relaxed using conjugate gradient (CG) minimization, and the cell shape and internal atomic coordinates were completely relaxed until the electron self-consistent convergence criterion of the electron iteration was 10 −6 eV/atom and the Hellmann-Feynman force [32] was less than 0.01 eV/Å.In the calculation process, the PAW_GGA pseudopotentials of Al, Sc, La, Y_sv, and (Ce-Lu)_3 are used to simulate the size of the real potential.The energy-strain method was used for the mechanical calculation, the VASPKIT code [33] was used to fit the equation of state (EOS), and the strain values were −9% ~9%.The thermodynamic properties were calculated by establishing the Debye-Grüneisen model [34,35] and using the Gibbs2 code [36].Although it is possible to calculate the thermodynamic properties of the material using the quasi harmonic approximation model (QHA), it is still a complex and expensive task to calculate the full vibration spectrum of the crystal on the volume grid.At the same time, compared with a fixed Poisson's ratio v in Debye-Grüneisen model, the exact Poisson's ratio v of each alloy is calculated using energy strain method in this experiment, which further improves the accuracy on the basis of the original.However, the Debye model only considers the contribution of the long wave part but not the short wave part, which makes the error of calculation at a high temperature larger.At the same time, the Al 11 RE 3 unit cell model required for the calculation was derived by the program software of VESTA to visualize the 3D structure, as shown in Figure 1a.It can be seen from the following figure that the structure belongs to the orthorhombic crystal system.
accuracy on the basis of the original.However, the Debye model only considers the contribution of the long wave part but not the short wave part, which makes the error of calculation at a high temperature larger.At the same time, the Al11RE3 unit cell model required for the calculation was derived by the program software of VESTA to visualize the 3D structure, as shown in Figure 1a.It can be seen from the following figure that the structure belongs to the orthorhombic crystal system.

Enthalpy of Formation and Stability
For the purpose of evaluating the relative stabilities of the Al11RE3 intermetallics in Al alloys, we calculated the formation enthalpy ΔHf according to the following equation [37] where  is the total energy of different compounds;  and  are the calculated total energy of the elements of Al and RE in a steady state, respectively.Parameters x and y are the number of Al and RE atoms in the unit cell, respectively.Based on the calculated formation enthalpy ΔHf of Al11RE3 intermetallics and the previous theoretical data [20] at 0 K and 0 GPa, as shown in Figure 1b, we can obtain the following conclusions: (1) The calculated values were in good agreement with the previous theoretical predictions, which indicates that the present calculations are accurate and reliable; (2) the formation enthalpies of the Al11RE3 intermetallics are all negative, indicating that they are stable compounds; (3) with an increase in the atomic number, the formation enthalpy decreases first and then increases.The maximum and minimum values of the ΔHf appear at the Sc and La elements, respectively, which indicates that the thermodynamic stability of Al11La3 is the strongest, followed by Al11Ce3, while the thermodynamic stability of Al11Sc3 is the worst.These results are in good agreement with the experimental findings that Al11La3 and Al11Ce3 particles are the most found compounds in the experiment [38].
In this experiment, the formation enthalpy  of the aluminum alloy Al11RE3 doped with 11 rare earth elements was calculated according to the formula.Compared with the existing experimental values [20], as shown in Figure 2, it can be clearly seen that the existing experimental values of Al11La3 and Al11Ce3 are in good agreement with the existing theoretical predicted values, which are −0.40 kJ/mol and −0.38 kJ/mol, respectively.We can observe that the formation enthalpies of the binary aluminum alloys are all negative, indicating that there are exothermic reactions when they form compounds, and that they are all stable compounds [20].At the same time, Figure 1b also helps us better understand the relationship between the formation enthalpy Δ and the atomic number of rare earth elements.We can clearly observe that the overall curve shows a trend that first decreases and then increases.Moreover, the formation enthalpy Δ reaches the

Enthalpy of Formation and Stability
For the purpose of evaluating the relative stabilities of the Al 11 RE 3 intermetallics in Al alloys, we calculated the formation enthalpy ∆H f according to the following equation [37] ∆H where E tot is the total energy of different compounds; E Al and E RE are the calculated total energy of the elements of Al and RE in a steady state, respectively.Parameters x and y are the number of Al and RE atoms in the unit cell, respectively.Based on the calculated formation enthalpy ∆H f of Al 11 RE 3 intermetallics and the previous theoretical data [20] at 0 K and 0 GPa, as shown in Figure 1b, we can obtain the following conclusions: (1) The calculated values were in good agreement with the previous theoretical predictions, which indicates that the present calculations are accurate and reliable; (2) the formation enthalpies of the Al 11 RE 3 intermetallics are all negative, indicating that they are stable compounds; (3) with an increase in the atomic number, the formation enthalpy decreases first and then increases.The maximum and minimum values of the ∆H f appear at the Sc and La elements, respectively, which indicates that the thermodynamic stability of Al 11 La 3 is the strongest, followed by Al 11 Ce 3 , while the thermodynamic stability of Al 11 Sc 3 is the worst.These results are in good agreement with the experimental findings that Al 11 La 3 and Al 11 Ce 3 particles are the most found compounds in the experiment [38].
In this experiment, the formation enthalpy ∆H f of the aluminum alloy Al 11 RE 3 doped with 11 rare earth elements was calculated according to the formula.Compared with the existing experimental values [20], as shown in Figure 2, it can be clearly seen that the existing experimental values of Al 11 La 3 and Al 11 Ce 3 are in good agreement with the existing theoretical predicted values, which are −0.40 kJ/mol and −0.38 kJ/mol, respectively.We can observe that the formation enthalpies of the binary aluminum alloys are all negative, indicating that there are exothermic reactions when they form compounds, and that they are all stable compounds [20].At the same time, Figure 1b also helps us better understand the relationship between the formation enthalpy ∆H f and the atomic number of rare earth elements.We can clearly observe that the overall curve shows a trend that first decreases and then increases.Moreover, the formation enthalpy ∆H f reaches the maximum value at the Sc element and then shows a downward trend; the lanthanide rare earth element La reaches a minimum value of −0.40 eV/atom, and the atomic number increases until the Lu element shows a generally upward trend; however, the formation enthalpy ∆H f of the Sc element is not exceeded.This shows that the thermodynamic stability of Al 11 La 3 is the strongest [38], while the thermodynamic stability of Al 11 Sc 3 is the weakest.maximum value at the Sc element and then shows a downward trend; the lanthanide rare earth element La reaches a minimum value of −0.40 eV/atom, and the atomic number increases until the Lu element shows a generally upward trend; however, the formation enthalpy Δ of the Sc element is not exceeded.This shows that the thermodynamic stability of Al11La3 is the strongest [38], while the thermodynamic stability of Al11Sc3 is the weakest.

Elastic Properties
Elastic constants and elastic moduli are important components in understanding the mechanical properties of materials [39].The elastic constant can be used to determine the mechanical stability of the material, reflect the ability of the crystal to resist external stress [20], and obtain important information on the stability of the aluminum alloy unit cell structure and the anisotropic characteristics of the crystal.By calculating the elastic constant of the alloy [40], the elastic modulus such as bulk modulus B, shear modulus G, Young's modulus E, and Poisson's ratio ν can be further determined to determine the mechanical properties of the material.The elastic constant is first calculated from a small strain of a uniform unit cell and is then determined by using the corresponding change in the total energy [41]; finally, the elastic modulus is derived using the Voigt-Reuss-Hill approximation [42].
Cij is the crystal's response to an applied external stress, while the number of elastic constants lies in the symmetry of the structure.Al11RE3 is an orthorhombic structure with nine independent elastic constants ( ,  ,  ,  ,  ,  ,  ,  ,   ), whose elastic constant matrix is in the form of 6 × 6.On the basis of the Born criterion [43], according to the necessary and sufficient Born criterion under the simplified orthogonal system [44], the following criteria can be obtained: Table 1 and Figure 2 present the relevant elastic constants and elastic moduli, respectively, which are used to determine the mechanical stability of aluminum alloys and the mechanical properties of materials.
Hill proved that the equations obtained using the Voigt-Reuss description of the bulk modulus B and the shear modulus G are the upper and lower limits of the true crystal constant, respectively, which is exactly the arithmetic mean provided by Voigt and Reuss; moreover, the following obtains the orthorhombic crystal formulas for calculating the elastic modulus and shear modulus [42]

Elastic Properties
Elastic constants and elastic moduli are important components in understanding the mechanical properties of materials [39].The elastic constant can be used to determine the mechanical stability of the material, reflect the ability of the crystal to resist external stress [20], and obtain important information on the stability of the aluminum alloy unit cell structure and the anisotropic characteristics of the crystal.By calculating the elastic constant of the alloy [40], the elastic modulus such as bulk modulus B, shear modulus G, Young's modulus E, and Poisson's ratio ν can be further determined to determine the mechanical properties of the material.The elastic constant is first calculated from a small strain of a uniform unit cell and is then determined by using the corresponding change in the total energy [41]; finally, the elastic modulus is derived using the Voigt-Reuss-Hill approximation [42].
C ij is the crystal's response to an applied external stress, while the number of elastic constants lies in the symmetry of the structure.Al 11 RE 3 is an orthorhombic structure with nine independent elastic constants (C 11 , C 12 , C 13 , C 22 , C 23 , C 33 , C 44 , C 55 , and C 66 ), whose elastic constant matrix is in the form of 6 × 6.On the basis of the Born criterion [43], according to the necessary and sufficient Born criterion under the simplified orthogonal system [44], the following criteria can be obtained: (2) Table 1 and Figure 2 present the relevant elastic constants and elastic moduli, respectively, which are used to determine the mechanical stability of aluminum alloys and the mechanical properties of materials.
Hill proved that the equations obtained using the Voigt-Reuss description of the bulk modulus B and the shear modulus G are the upper and lower limits of the true crystal constant, respectively, which is exactly the arithmetic mean provided by Voigt and Reuss; moreover, the following obtains the orthorhombic crystal formulas for calculating the elastic modulus and shear modulus [42] where Using the elastic constants, Young's modulus E and Poisson's ratio ν can be expressed as: According to the criterion of mechanical stability and the calculated value of the elastic constant, we can observe from Figure 2a that their values are positive, and that all the rare earth-aluminum alloys meet the conditions of mechanical stability, indicating that they are stable.We can also see from Figure 2a that there are two trends in relation to the La element, as the node shows a trend of first increasing and then decreasing, or first decreasing and then increasing.
Bulk modulus B can be used to measure the average atomic bond strength of a material and is inversely proportional to compressibility; the greater the bulk modulus, the greater the resistance of the material to compression [45].Observing the bulk modulus B of Al 11 RE 3 and aluminum in Figure 2b, it is found that the bulk modulus of all aluminum rare earth alloys except Al 11 Sc3 (83.21 GPa) is smaller than that of Al (81.71 GPa), of which Al 11 La 3 has the smallest bulk modulus B, and that the minimum is 73.15GPa.This shows that Al 11 Sc 3 has the strongest compressive ability, the most stable structure, is not easily deformed, has the strongest rigidity and the highest hardness.On the contrary, Al 11 La 3 is easy to deform, weak in rigidity, and is the lowest in hardness.
Young's modulus E and shear modulus G are the most important parameters for evaluating the strength of materials.Young's modulus E can represent the ability to resist deformation; the larger the value, the harder the material.In addition, the shear modulus G can represent the ability to resist shear stress, and the larger the shear modulus, the more obvious the directional bonds between atoms [45].It can be seen from Figure 2b that the shear modulus G and Young's modulus E of Al 11 RE 3 have the same trend, with La atoms as the nodes, of first rising and then decreasing with an increase in the atomic number.With Al 11 Sc 3 as the starting point, Young's modulus E increases from 64.00 GPa to 115.81 GPa for Al 11 La 3 and then slowly decreases to Al 11 Lu 3 (81.91GPa) with an increase in the atomic number.Similarly, the shear modulus G also showed the same trend, rising from 23.33 GPa at Al 11 Sc 3 to 46.88 GPa at Al 11 La 3 and then slowly decreasing to 30.78 GPa.The Young's modulus E and shear modulus G of aluminum are 76.11GPa and 28.30GPa, respectively.It is obvious that the addition of RE elements makes the Young's modulus E and shear modulus G of the aluminum alloy higher than that of aluminum, except for Sc.Among them, Al 11 La 3 has the highest Young's modulus E value and shear modulus G value, which are greatly improved compared with aluminum, which are increased by about 52% and 66%, respectively.After adding rare earth elements (except Sc), the stiffness of the aluminum alloy material increases, the toughness is good, the material is not easily deformed, and the tensile strength is significantly improved.In contrast, the shear modulus G of Al 11 Sc 3 is about 18% lower than that of aluminum, and the Young's modulus E is about 16% lower than that of aluminum, indicating that the element Sc does not enhance its elastic properties, and that the formed aluminum-based alloy is volatile.Moreover, stiffness and strength decrease, and elastic deformation increases.
Poisson's ratio ν can be used to predict the plasticity of the material, thereby quantifying the stability of the crystal against shearing [20,46].It also represents the ratio of the lateral strain to the longitudinal strain of the solid medium; therefore, the range of Poisson's ratio is greater than 0, the Poisson's ratio is in the range of 0.1~0.3, and the material is hard; when Poisson's ratio is 0.4~0.5, the material is generally fluid or soft.Figure 2c shows that Al 11 La 3 has the lowest Poisson's ratio of 0.24, which is also in the range of 0.1 to 0.3, which belongs to an alloy with high hardness.Generally speaking, aluminum has low hardness, high plasticity, and is easy to process.It can be seen from Figure 2c  B/G can also be used to characterize the plasticity of the material: with a B/G of 1.75 as the boundary [47], below 1.75, the material is ductile, and, above 1.75, the material is brittle [20,45].As shown in Figure 2c, the Poisson's ratio v and B/G of the aluminum alloy both decrease first and then increase with an increase in the rare earth atomic number, taking the La element as the node.Among them, the B/G of Al 11 La 3 , Al 11 Ce 3 , and Al 11 Pr 3 are all lower than 1.75, with values of 1.56, 1.67, and 1.71, respectively, showing certain ductility.The value of Al 11 La 3 is the smallest, which indicates that the ductility is good and that the rolling degree is improved, which is beneficial to the production of precision materials.The B/G value of Al 11 Sc 3 is up to 3.567, which is not conducive to rolling.

The Thermodynamic Properties
In order to understand the thermodynamic properties of the second-phase particles of the material, this calculation uses the Debye theory (DT) [35] to obtain a model of crystal vibration [48][49][50], namely, the Debye-Grüneisen model, and the key is to calculate the Debye function Θ D (V) through the Debye cutoff frequency ω D [51][52][53].At high temperature [54][55][56], the Al 11 RE 3 phase will be decomposed into Al and Al 2 RE; in this work, we adopt P = 0 at a temperature of 0-1000 K [57] for the thermodynamic calculation of the aluminum alloy, according to the following equation [51][52][53]: Crystals 2023, 13, 347 where the h and κ B are Plank's and Boltzmann's constants, respectively.The ω D can be obtained using the Debye sound velocity v D : where the V is the volume of the solid.There are three different types of sound velocity in solids, all of which are anisotropic.For isotropic media, Anderson expresses the low temperature average for the sound velocity for an isotropic crystal as: where v L and v S are longitudinal and transverse sound velocities and can be expressed as v L = (L/ρ) and v S = (S/ρ), respectively.With Equations ( 14)-( 17), a simpler method can be obtained to calculate v D and Θ D ; then, they are used to calculate the longitudinal and transverse moduli as follows: Then, the v D can be obtained: It has been proved that the method can obtain reasonable Debye sound velocity.Therefore, the Θ D can be calculated from the bulk modulus and density using the following equation: where M is the atomic weight and r0 is the equilibrium Wigner-Seitz radius.
After inserting the Θ D into quasi harmonic formula, the Helmholtz free energy F, entropy S, heat capacity CV, and thermal expansion α can be, respectively, described by [34,53]: At a temperature of 0-1000 K, the Gibbs free energy F, entropy S, thermal melt C V , and thermal expansion coefficient α are shown in Figure 3a-d.
At a temperature of 0-1000 K, the Gibbs free energy F, entropy S, thermal melt CV, and thermal expansion coefficient  are shown in Figure 3a-d.As shown in Figure 3a-c, it is found that the thermal expansion coefficient α, constant volume specific heat CV, and entropy S of the aluminum and rare earth aluminum alloys all increase with an increase in temperature of ~1000 K; however, the increase range is different.
The thermal expansion coefficient α can be used to characterize the ability of an object to expand and contract due to temperature changes.In Figure 3a, the initial and final values of thermal expansion coefficients and growth rates of 17 Al11RE3 alloys are basically the same, i.e., in the order of ~2.89 × 10 −5 /K.In detail, at 0-200 K, there is a relatively large increase, and, after 200 K, the growth rate is slow.The change of thermal expansion coefficient α of an Al of ~5.00 × 10 −5 /K with T is significantly larger than that of Al11RE3, indicating that the thermal stability of the aluminum element is poor, and that the addition of rare earth elements greatly improves the thermal stability of the aluminum element.Aluminum alloy can maintain a good stable state even at high temperature, suitable for high temperature resistant aluminum alloy materials.
As shown in Figure 3b, the constant volume specific heat CV of Al and Al11RE3 both showed a significant increase, then a slow increase, and, finally, a gradual increase.The constant volume specific heat CV of RE aluminum alloys grows the most at 0-200 K, increasing exponentially, and then the growth rate slows down significantly at 200-600 K, and basically stabilizes after 600 K.The specific heat CV value Al11RE3 of ~24.86 J/mol/K and the growth rate are both higher than those of the Al of ~24.81 J/mol/K.The growth rates of rare earth aluminum alloys are basically the same, and the constant volume specific heat CV of Al11RE3 tends to be consistent.Since the constant volume specific heat CV As shown in Figure 3a-c, it is found that the thermal expansion coefficient α, constant volume specific heat C V , and entropy S of the aluminum and rare earth aluminum alloys all increase with an increase in temperature of ~1000 K; however, the increase range is different.
The thermal expansion coefficient α can be used to characterize the ability of an object to expand and contract due to temperature changes.In Figure 3a, the initial and final values of thermal expansion coefficients and growth rates of 17 Al 11 RE 3 alloys are basically the same, i.e., in the order of ~2.89 × 10 −5 /K.In detail, at 0-200 K, there is a relatively large increase, and, after 200 K, the growth rate is slow.The change of thermal expansion coefficient α of an Al of ~5.00 × 10 −5 /K with T is significantly larger than that of Al 11 RE 3 , indicating that the thermal stability of the aluminum element is poor, and that the addition of rare earth elements greatly improves the thermal stability of the aluminum element.Aluminum alloy can maintain a good stable state even at high temperature, suitable for high temperature resistant aluminum alloy materials.
As shown in Figure 3b, the constant volume specific heat C V of Al and Al 11 RE 3 both showed a significant increase, then a slow increase, and, finally, a gradual increase.The constant volume specific heat C V of RE aluminum alloys grows the most at 0-200 K, increasing exponentially, and then the growth rate slows down significantly at 200-600 K, and basically stabilizes after 600 K.The specific heat C V value Al 11 RE 3 of ~24.86 J/mol/K and the growth rate are both higher than those of the Al of ~24.81 J/mol/K.The growth rates of rare earth aluminum alloys are basically the same, and the constant volume specific heat C V of Al 11 RE 3 tends to be consistent.Since the constant volume specific heat C V of an RE aluminum alloy is greater than that of aluminum, it means that, at a certain temperature, the heat absorption or heat dissipation ability of an aluminum alloy after adding RE elements is stronger than that of aluminum.
The entropy S value indicates that the system changes from an ordered state to a disordered state.The entropy S value increases with temperature because molecular motion becomes more violent with an increasing temperature and will deviate more from its equilibrium position.Looking at Figure 3c, the entropy S curves of aluminum of Crystals 2023, 13, 347 9 of 12 ~61.27J/mol/K and aluminum alloys of ~66.86 J/mol/K have been increasing with an increase in temperature; moreover, although the growth rate has decreased, the growth rate and growth value are higher than those of aluminum.The entropy growth and growth trend among the 17 kinds of aluminum alloys are very close, indicating that the addition of RE elements does not mean that the aluminum alloys are greatly affected by temperature and have good consistency.
As shown in Figure 3d, it is found that the Gibbs free energies F of both aluminum and Al 11 RE 3 alloys are less than 0, which have an Al value of ~−3.712 eV/atom and Al 11 RE 3 alloys of ~4.137 eV/atom, indicating that they exist stably.At 0 and 1000 K, the Gibbs free energies of lanthanides are also very close, while the initial and final values of Sc and Y are smaller than those of lanthanides.It is also clearly observed that the variation curves of the Gibbs free energy F of the aluminum alloy with the temperature T all decrease, and that the decreasing range is consistent.At the same time, the Gibbs free energy of 17 RE aluminum alloys is lower than that of aluminum, which means that the addition of rare earth elements improves the stability of aluminum.
The variation of the thermal expansion coefficient a of the aluminum element with the temperature T is significantly larger than that of Al 11 RE 3 , indicating that the thermal stability of the aluminum is poor, and that the addition of rare earth elements greatly improves the thermal stability of the aluminum alloy, making the aluminum alloy exist even at high temperatures.Nevertheless, it can still maintain good stability in the state and is suitable for high temperature-resistant aluminum alloy materials.

Conclusions
In summary, the stability, mechanical properties, and thermodynamic properties of eleven Al 11 RE 3 intermetallics at 0-1000 K were investigated by combining the firstprinciples calculations and the Debye model.The results of the study are as follows: (1) The enthalpy of formation ∆H f of Al 11 RE 3 intermetallics are all negative values, all belong to exothermic reactions, and are all stable compounds.The thermodynamic stability of Al 11 La 3 is the strongest, while the thermodynamic stability of Al 11 Sc 3 is the weakest; (2) Comparing the mechanical constants of the eleven Al 11 RE 3 intermetallics, it is found that the stiffness, toughness, and tensile strength of the Al 11 RE 3 intermetallics (RE = La, Ce, Pr, Nd, Sm) are significantly greater than those of aluminum, which indicates that these Al 11 RE 3 intermetallics act as strengthen phases that can effectively improve the mechanical property of aluminum; (3) The Debye-Grüneisen model was used to analyze the thermodynamic properties (Gibbs free energy, constant volume specific heat, entropy, and thermal expansion coefficient) of the eleven Al 11 RE 3 intermetallics at a range of 0-1000 K, and the results consistently showed that the Al 11 RE 3 intermetallics act as strengthen phases that can improve the thermal stability of the aluminum.In the process of temperature change, Al 11 Pm 3 and Al 11 Sm 3 are most affected by temperature and easily become unstable; their heat absorption and heat dissipation capabilities are also stronger than other RE aluminum alloys; (4) With an increase in the atomic number, different change trends were observed for the formation enthalpy ∆H f , bulk modulus B, Young's modulus E, and shear modulus G, which decreases first.

Figure 1 .
Figure 1.(a) The unit cell model of Al11RE3; (b) the comparison between the calculated and experimental values of the formation enthalpy of Al11RE3.Ref. [a] represents [20].

Figure 1 .
Figure 1.(a) The unit cell model of Al 11 RE 3 ; (b) the comparison between the calculated and experimental values of the formation enthalpy of Al 11 RE 3 .Ref. [a] represents [20].

Figure 2 .
Figure 2. (a) Calculated elastic constant of Al 11 RE 3 ; (b) Young's modulus E, shear modulus G, bulk modulus B values; (c) Poisson's ratio v and B/G values.
that the Poisson's ratios of Al, Al 11 Sc 3 , and Al 11 Er 3 ~Al 11 Lu 3 are all between 0.3 and 0.4, with certain plasticity and low hardness.At the same time, Poisson's ratio v can reflect the heterogeneity of materials.Taking the Poisson's ratio v = 0.26 as the boundary, except for Al 11 Ce 3 and Al 11 La 3 , all aluminum-based alloys exhibit uneven deformation, and the Poisson's ratio v of Al 11 Sc 3 (0.37) is the largest, and the heterogeneity is obvious.

Figure 3 .
Figure 3. Aluminum and Al11RE3 at T = 0-1000 K, (a) thermal expansion coefficient  as a function of temperature T; (b) constant volume specific heat CV; (c) entropy S; (d) Gibbs free energy F.

Figure 3 .
Figure 3. Aluminum and Al 11 RE 3 at T = 0-1000 K, (a) thermal expansion coefficient α as a function of temperature T; (b) constant volume specific heat C V ; (c) entropy S; (d) Gibbs free energy F.

Table 1 .
The elastic constant C ij of Al 11 RE 3 intermetallics (units: GPa).