First-Principles Study of Elastic Properties and Electronic Properties of Al-Ni-Ce Ternary Intermetallic Compounds

: The materials of the Al-Ni-Ce ternary system have promising application potential in additive manufacturing, and the systematic study of Al-Ni-Ce IMCs is an important part of the design and performance development of Al-Ni-Ce alloys. In this paper, the mechanical properties and electronic properties of seven major Al-Ni-Ce IMCs are studied using ﬁrst-principles calculations, and the differences in their mechanical properties are discussed. The enthalpy of formation and cohesion energy of the Al-Ni-Ce IMCs showed that Ce 4 Al 23 Ni 6 had the highest enthalpy of formation and cohesion energy, which were − 49.57 kJ/mol and − 4.47 kJ/mol, respectively. This demonstrated that all Al-Ni-Ce IMCs had excellent thermodynamic stability. The elastic modulus calculation of the Al-Ni-Ce IMCs showed that the Young’s modulus of CeAl 3 Ni 2 was the largest at 178.15 GPa, the bulk modulus of CeAlNi 4 was the largest at 125.78 GPa, and the shear modulus of CeAlNi 2 was the largest at 62.53 GPa. This proved that among the Al-Ni-Ce IMCs, CeAl 3 Ni 2 had the greatest stiffness, CeAlNi 4 had the strongest resistance to uniform deformation, and CeAl 3 Ni 2 had the strongest resistance to plastic deformation. The electronic density of states of Al-Ni-Ce IMCs was analyzed, and the results showed that the main bonding types of Al-Ni-Ce IMCs were metallic and covalent bonds. This work can provide a theoretical basis for the material design and development of the Al-Ni-Ce ternary system.


Introduction
Al alloys are widely used in industrial manufacturing fields such as aerospace and vehicle engineering due to their ease of processing, excellent strength-to-mass ratio, and unique mechanical properties [1,2].A large number of long-range ordered intermetallic compounds (IMCs) are present in Al alloys, and their characteristics clearly distinguish them from the Al matrix [3].Although the brittleness exhibited by IMCs can limit the application of IMCs, understanding the nature of IMCs in Al alloying systems is nonetheless extremely important for the strengthening mechanisms and performance enhancement of Al alloys.
In recent years, researchers have found that Al-Ni-Ce alloys have a narrow freezing range similar to that of Al-Si10-Mg alloys, which has immeasurable potential for additive manufacturing [4].There are also some correlations between the IMCs of the Al-Ni-Ce system and the creep properties of Al-Ni-Ce alloys.Michi et al. [5] found that Al-Ni-Ce alloys to which a small amount of Mn was added exhibited significant creep resistance due to the high-volume fraction of IMCs in the eutectic microstructure.The results of Wu et al. [6] showed that creep strengthening of Al-Ni-Ce alloys arises from load transfer from the Al matrix to IMCs.Therefore, the study of Al-Ni-Ce IMCs is valuable for the application development and strengthening mechanism of Al-Ni-Ce alloys.
to refine the thermodynamic parameters of the ternary equilibrium pha Tang et al. [8] studied the phase equilibria of the Al-Ce-Ni system in the at.% Ce at 800 °C and systematically investigated the ternary phases and rameters as a function of the solid solution composition, and based on rected the modified partial isothermal cross-section of the Al-Ce-Ni sy Wang et al. [9] investigated the effect of elemental Ce on the structural, e tronic properties of NiAl alloy structures and showed that Ce led to latt NiAl and increased the hardness and ductility of NiAl.Delsante et al. [ peculiarities of the formation of aluminum-rich phases in the Ce-Ni-Al and combined with microstructure and X-ray diffraction, the existence of CeNiAl4 was confirmed.These studies have continuously improved the parameters of Al-Ce-Ni IMCs.However, the load transfer between Al-C the Al matrix and the mechanical properties of Al-Ce-Ni IMCs are closely fore, Al-Ce-Ni IMCs' phase stability, mechanical properties, and other asp further studied.
First-principles calculations based on density-functional theory hav obtain information at the atomic and electronic microscopic level, wh unique advantage in studying the crystal structure, mechanical properties properties of IMCs [11][12][13].Therefore, first-principles calculations have b ingly popular among researchers in the study of IMCs [14 -16].In this pap IMCs (CeAl4Ni, Ce4Al23Ni6, CeAl5Ni2, CeAl3Ni2, CeAlNi4, CeAl2Ni, and been investigated using first-principles calculations in the Al-Ni-Ce system tributions in the Al-Ni-Ce system are presented in Figure 1 [17,18].Their cr enthalpy of formation and cohesion energy, elastic properties, and electr are investigated, and these parameters can reflect the stability of the pha stability and mechanical properties under stress conditions, and the bond tics of Al-Ni-Ce IMCs.This study is expected to provide a theoretical bas opment of creep-strengthening mechanisms and properties of Al-Ni-Ce a

Calculation Method
All calculations in this paper were performed using the open-sourc sion of the CASTEP program [19].The electron exchange correlation was ing the Generalized Gradient Approximation (GGA) with PBE generali ultrasoft pseudopotential was adopted because it can be calculated at a l ergy.The Kohn-Sham equation was solved using the self-concordant field the ground state through an electron minimization process [21], sett

Calculation Method
All calculations in this paper were performed using the open-source academic version of the CASTEP program [19].The electron exchange correlation was investigated using the Generalized Gradient Approximation (GGA) with PBE generalization [20].The ultrasoft pseudopotential was adopted because it can be calculated at a lower cutoff energy.The Kohn-Sham equation was solved using the self-concordant field (SCF) to reach the ground Metals 2023, 13,1921 3 of 15 state through an electron minimization process [21], setting the energy tolerance of the SCF to 5.0 × 10 −6 eV/atom.During the calculation process, the electronic configurations of Ce-4f5s5p5d6s, Ni-3d84s2, and Al-3s2s3p were taken as the valence electrons.The structure optimization of the equilibrium crystals was carried out through the BFGS method, with the total energy change per atom and the force convergence tolerance set to 5.0 × 10 −5 eV/atom and 0.01 eV/Å, respectively.The crystal structure of Al-Ni-Ce IMCs is shown in Figure 2

Crystal Structure
The Al-Ni-Ce IMCs studied have a total of three different crystal sys them, CeNi2Al3 and CeNiAl belong to the hexagonal crystalline system, w and P̅ 62m space groups, respectively.CeNiAl4, CeNi2Al5, CeNi4Al, and CeN to the orthorhombic crystal system, CeNiAl4, CeNi4Al, and CeNiAl2 all Cmmm space group, and CeNi2Al5 belongs to the Immm space group.C monoclinic crystal system with a C2/m space group.The crystal constan IMCs were calculated, and the results of the calculations and other papers references are presented in Table 1.The calculated results are derived fro rium lattice parameters obtained after each cell geometry optimization.Th stants of the Al-Ni-Ce IMCs are compared with the experimental values an the experimental values with the maximum and minimum error ranges of 2. respectively.

Crystal Structure
The Al-Ni-Ce IMCs studied have a total of three different crystal systems.Among them, CeNi 2 Al 3 and CeNiAl belong to the hexagonal crystalline system, with P6/mmm and P62m space groups, respectively.CeNiAl 4 , CeNi 2 Al 5 , CeNi 4 Al, and CeNiAl 2 all belong to the orthorhombic crystal system, CeNiAl 4 , CeNi 4 Al, and CeNiAl 2 all belong to the Cmmm space group, and CeNi 2 Al 5 belongs to the Immm space group.Ce 4 Ni 6 Al 23 is a monoclinic crystal system with a C2/m space group.The crystal constants of Al-Ni-Ce IMCs were calculated, and the results of the calculations and other papers' comparative references are presented in Table 1.The calculated results are derived from the equilibrium lattice parameters obtained after each cell geometry optimization.The crystal constants of the Al-Ni-Ce IMCs are compared with the experimental values and are close to the experimental values with the maximum and minimum error ranges of 2.6% and 0.49%, respectively.
In addition, the volume and density of the Al-Ni-Ce IMCs were taken separately and analyzed.Figure 3   Exp.Experiment.
In addition, the volume and density of the Al-Ni-Ce IMCs were taken separately analyzed.Figure 3

Enthalpy of Formation and Cohesion Energy
The calculation of the enthalpy of formation and cohesion energy allows the ph stability of Al-Ni-Ce IMCs to be analyzed from an energetic point of view.The entha of formation (∆    ) is equal to the energy of formation (    ) for Ni-Ce IMCs in the computational environment without taking into account temperat and air pressure [26], specifically expressed as The calculation of the enthalpy of formation and the cohesion energy can be pressed as [27]

Enthalpy of Formation and Cohesion Energy
The calculation of the enthalpy of formation and cohesion energy allows the phase stability of Al-Ni-Ce IMCs to be analyzed from an energetic point of view.The enthalpy of formation (∆H f (Al m Ni n Ce l )) is equal to the energy of formation (E f (Al m Ni n Ce l )) for Al-Ni-Ce IMCs in the computational environment without taking into account temperature and air pressure [26], specifically expressed as The calculation of the enthalpy of formation and the cohesion energy can be expressed as [ The enthalpy of formation is defined as the difference between the total energy of a compound and the energy of its constituent elements in a steady state.The cohesion energy is defined as the work required when crystals break down into free atoms.The enthalpy of formation and cohesion energy of Al-Ni-Ce IMCs are shown in Table 2.The calculation results of the enthalpy of formation showed that the absolute enthalpy of formation of CeAl 3 Ni 2 is the maximum, indicating that it has the strongest alloying ability, while the opposite is true for Ce 4 Al 23 Ni 6 .The alloying ability of the seven Al-Ni-Ce IMCs was in the order of CeAl 3 Ni 2 > CeAl 5 Ni 2 > CeAlNi > CeAl 4 Ni > CeAl 2 Ni > CeAlNi 4 > Ce 4 Al 23 Ni 6 .The calculated results are in good agreement with the experimental results, which confirms the rationality and reliability of the calculation method.The calculated enthalpy of formation and cohesion energy of Al-Ni-Ce IMCs are less than 0, which indicates that the process of combining these atoms to form the compound is exothermic.The thermodynamic stability of IMCs is closely related to the enthalpy of formation and cohesion energy, and the lower the enthalpy of formation and cohesion energy, the better the thermodynamic stability.Therefore, Al-Ni-Ce IMCs have great thermal stability.Due to the particularity of rare earth atoms, the enthalpy of formation and cohesion energy of Al-Ni-Ce IMCs will be affected by the content ratio of Ce atoms.To further illustrate the relationship between the enthalpy of formation and cohesion energy of Al-Ni-Ce IMCs and the Ce atom under vacuum conditions (0 K, 0 Pa), the relationship between the enthalpy of formation and the cohesion energy and the content of Ce atoms is shown in Figure 4.The relationship between the enthalpy of formation of Al-Ni-Ce IMCs and the Ce content is shown in Figure 4a.Clearly, the enthalpy of formation is higher on the side with less Ce content.According to other works [28], the solubilization of Ce in Al-Ni IMCs is difficult.Therefore, the content of Ce is less, which is more conducive to the alloying of Al-Ni-Ce IMCs.The relationship between the cohesion energy and Ce content of Al-Ni-Ce IMCs is shown in Figure 4b, and the absolute value of cohesion energy gradually increases with the increase in Ce content.This indicates that the presence of Ce atoms strengthens the association of Al-Ni-Ce IMCs, which are difficult to decompose into free atoms.

Elastic Constants
The elastic stiffness constant of the crystal,  , reflects the deformation capacity of the crystal under stress conditions, where  and  denote transverse deformation and  denotes longitudinal deformation, and the larger the value, the better the mechanical properties.The elastic stiffness constant  can be expressed as [29]  = (4) where  is the energy of the crystal,  is its equilibrium volume, and  denotes a strain.The  reflects the mechanical stability of the crystal, and the stability criterion is determined differently for different crystal systems.The calculated elastic stiffness constants of Al-Ni-Ce IMCs are presented in Table 3. From the calculation results, it can be seen that the elastic stiffness constants of the Al-Ni-Ce IMCs depend on the stability determination criterion of the crystalline systems they correspond to, which indicates that the Al-Ni-Ce IMCs can be mechanically stabilized.
For the orthorhombic crystals [30] (CeAl2Ni, CeAlNi4, CeAl5Ni2, CeAl4Ni): For the hexagonal crystals [31] (CeAl3Ni2, CeAlNi): For the monoclinic crystal [32] (Ce4Al23Ni6): where E is the energy of the crystal, V 0 is its equilibrium volume, and ε denotes a strain.The C ij reflects the mechanical stability of the crystal, and the stability criterion is determined differently for different crystal systems.The calculated elastic stiffness constants of Al-Ni-Ce IMCs are presented in Table 3. From the calculation results, it can be seen that the elastic stiffness constants of the Al-Ni-Ce IMCs depend on the stability determination criterion of the crystalline systems they correspond to, which indicates that the Al-Ni-Ce IMCs can be mechanically stabilized.

Elastic Modulus
The Voigt-Reuss-Hill method [33] was used to calculate shear modulus (G) and bulk modulus (B) for Al-Ni-Ce IMCs.As with the elastic constants, different crystal systems correspond to different equations of the calculation.
For the orthorhombic crystals (CeAl 2 Ni, CeAlNi 4 , CeAl 5 Ni 2 , CeAl 4 Ni): For the hexagonal crystals (CeAl 3 Ni 2 , CeAlNi): The average shear modulus (G) and bulk modulus (B), Young's modulus (E), and Poisson's ratio (ν) of Al-Ni-Ce IMCs were finally calculated as follows [34]: The calculated results of the Al-Ni-Ce IMCs elastic modulus are presented in Table 4. Pugh [35] proposed defining the brittleness of a material by the ratio of bulk modulus B to shear modulus, B/G.If the B/G ratio is greater than 1.75, the material is considered to be brittle.If the B/G ratio is greater than 1.75, the material is considered to be ductile, and if the B/G ratio is less than 1.75, the material is considered to be brittle.Based on the calculation of B/G ratios of Al-Ni-Ce intermetallic compounds, CeAl 2 Ni, CeAlNi, and CeAlNi 4 were defined as ductile materials, while CeAl 3 Ni 2 , CeAl 4 Ni, CeAl 5 Ni 2 , and Ce 4 Al 23 Ni 6 were defined as brittle materials.The shear modulus was used to evaluate the plastic deformation resistance and fracture behavior of the materials, with larger values indicating greater resistance to plastic deformation behavior.The calculated results in Table 4 show that the size order of the plastic deformation resistance of Al-Ni-Ce IMCs was CeAl  [36] concluded that the bulk modulus of a crystal is determined by the ability of the chemical bonding to resist compressive forces, and the atom's electronegativity strength determines the size of the bulk modulus of the crystal.The electronegativity strength of Al-Ni-Ce IMCs is influenced by the rare earth element Ce atoms, so there is a correlation between the bulk modulus of Al-Ni-Ce IMCs and the composition of Ce atoms.Young's modulus is a physical measure of the stiffness of a material; the higher the Young's modulus the stiffer the material.The order of magnitude of the Young's modulus of Al-Ni-Ce IMCs is CeAl The elastic modulus calculations of Al-Ni-Ce IMCs were compared with those of other materials, and the comparative results showed that the elastic modulus of Al-Ni-Ce IMCs was generally higher than that of Al-La IMCs [38], which proved that Al-Ni-Ce IMCs have more superior mechanical properties than some Al-RE binary IMCs.The shear modulus of Al-Ni-Ce IMCs was found to be generally higher than that of Al-Cu IMCs [39], which indicates that the Al-Ni-Ce system has a better plastic deformation capacity than the Al-Cu system.In comparison with the elastic modulus of the Mg-Al-Si IMCs [40], the elastic modulus of the Al-Ni-Ce IMCs was superior to that of the Mg-Al-Si system across the board.From these comparative results, it is clear that Al-Ni-Ce alloys have superior mechanical properties.

Anisotropic Elastic Properties
The spatial anisotropy coefficients of the modulus of elasticity can reflect the differences in the deformation of crystals subjected to forces in space.The contribution of the elastic modulus of the Al-Ni-Ce IMCs in all directions of space is shown in Figure 5.The elastic modulus can be presented using the visualization tool developed by Liao et al. [41].
When the three-dimensional image deviates from the spherical surface, it indicates that the elastic modulus behaves as anisotropies, and the larger the deviation, the more significant the anisotropies are.Therefore, as shown in Figure 5, the bulk modulus, shear modulus, Young's modulus, and Poisson's ratio of Al-Ni-Ce IMCs are all anisotropic.Among them, the most obvious anisotropy of bulk modulus is that of AlNiCe: the minimum value of bulk modulus is 49.0 GPa and the maximum value of bulk modulus is 31,315.5GPa along the axial directions <−0.99, −0.16, 0> and <−0.47, 0, 0.88>, respectively.The largest anisotropy of Young's modulus is that of CeAl 2 Ni: the minimum value of Young's modulus is 69.9 GPa and the maximum value is 164.1 GPa along the axial directions <0, 0, 1> and <−0.75, 0, −0.66>, respectively.The largest anisotropy of shear modulus is that of CeAl 2 Ni: the minimum value of shear modulus is 29.1 GPa and the maximum value is 58.6 GPa along the axial directions <−0.75, 0, 0.66> and <0, 0, 1>, respectively.The largest Poisson's ratio is that of CeAl 5 Ni 2 : the minimum value of Poisson's ratio is 0.2 and the maximum value is 0.4 along the axial directions <0.84, 0.54, 0> and <0, −1, 0>, respectively.However, the elastic anisotropy of the bulk modulus of CeAl 3 Ni 2 and CeAlNi 4 is very similar, which indicates that there will be no distortion, crack, or other defects between CeAl 3 Ni 2 and CeAlNi 4 under uniform stress conditions [42].
31,315.5 GPa along the axial directions <−0.99, −0.16, 0> and <−0.47, 0, 0.88>, respectively.The largest anisotropy of Young's modulus is that of CeAl2Ni: the minimum value of Young's modulus is 69.9 GPa and the maximum value is 164.1 GPa along the axial directions <0, 0, 1> and <−0.75, 0, −0.66>, respectively.The largest anisotropy of shear modulus is that of CeAl2Ni: the minimum value of shear modulus is 29.1 GPa and the maximum value is 58.6 GPa along the axial directions <−0.75, 0, 0.66> and <0, 0, 1>, respectively.The largest Poisson's ratio is that of CeAl5Ni2: the minimum value of Poisson's ratio is 0.2 and the maximum value is 0.4 along the axial directions <0.84, 0.54, 0> and <0, −1, 0>, respectively.However, the elastic anisotropy of the bulk modulus of CeAl3Ni2 and CeAlNi4 is very similar, which indicates that there will be no distortion, crack, or other defects between CeAl3Ni2 and CeAlNi4 under uniform stress conditions [42].The elasticity of each anisotropy can be expressed by the following equation [43]: When  = 1, this shows each isotropy, and the larger the value of  , the more The elasticity of each anisotropy can be expressed by the following equation [43]: When A X = 1, this shows each isotropy, and the larger the value of A X , the more drastic the crystal deformation.In addition, the elastic anisotropies of the Al-Ni-Ce IMCs were calculated separately and are presented in Table 5.The results of the elastic modulus anisotropy calculations are consistent with those in the 3D image.In addition, it can be seen from Table 5 that the elastic anisotropy of the Young's modulus of CeAl 4 Ni and CeAlNi 4 is the same, and the elastic anisotropy of the bulk modulus of CeAl 4 Ni and Ce 4 Al 23 Ni 6 is the same as that of the Poisson's ratio.

Partial Density of States
The PDOS is an excellent auxiliary tool for analyzing the electronic structure and chemical properties of materials.It can effectively reflect the bonding properties of IMCs.The results of the PDOS calculations for Al-Ni-Ce IMCs are presented in Figure 6.The results of the PDOS show that all Al-Ni-Ce IMCs have distinct metallic properties.In addition, the existence of resonance peaks near −17 eV for Ce and Al atoms with the same center of gravity proves the existence of covalent bonding properties between Al and Ce atoms.On the left side of the Fermi energy level, there is hybridization between the d-orbital electrons of the Ni atoms and the p-orbital electrons of the Al atoms, and the hybridization is more pronounced for CeAlNi 4 and CeAl 3 Ni 2 .There is hybridization of the f-orbital electrons of the Ce atoms and the p-orbital electrons of the Al atoms, in addition to the obvious hybridization of the s-orbital electrons of the Al atoms and the f-orbital electrons of the Ce atoms of Ce 4 Al 23 Ni 6 .The waveform curves of these electronic orbital hybridizations have the same position, but different centers of gravity.This suggests that the Al-Ni-Ce IMCs also have a metallic bond.The presence of covalent bonds can enhance the strength of Al-Ni-Ce IMCs compared to pure metal bonds [44].Therefore, the main bonding types of Al-Ni-Ce IMCs are covalent and metallic bonds with strong structural stability.
to the obvious hybridization of the s-orbital electrons of the Al atoms and the f-or electrons of the Ce atoms of Ce4Al23Ni6.The waveform curves of these electronic or hybridizations have the same position, but different centers of gravity.This suggests the Al-Ni-Ce IMCs also have a metallic bond.The presence of covalent bonds can enh the strength of Al-Ni-Ce IMCs compared to pure metal bonds [44].Therefore, the bonding types of Al-Ni-Ce IMCs are covalent and metallic bonds with strong struc stability.

Conclusions
The structure, thermodynamic properties, mechanical properties, and electronic structure of Al-Ni-Ce IMCs are systematically investigated in this work through first-principles calculations, and the conclusion is as follows: 1.
The calculated results of lattice parameters were verified to be close to the experimental values.

2.
The calculated enthalpy of formation and cohesion energy shows that the Al-Ni-Ce IMCs have good thermodynamic stability.The analysis of the fractional density of states reveals the prevalence of covalent and metallic bonds in the Al-NI-Ce IMCs.

Figure 1 .
Figure 1.Distribution of Al-Ni-Ce ternary intermetallic compounds in phase diag

Figure 1 .
Figure 1.Distribution of Al-Ni-Ce ternary intermetallic compounds in phase diagrams.
presents the relationship between the volume and density of Al-Ni-Ce IMCs and the Al, Ni, and Ce contents.The volume contour distribution of the Al-Ni-Ce system is shown in Figure 3a.It can be seen that the larger region of Al-Ni-Ce is mainly concentrated on the side of Al and Ce, and CeAl 2 Ni, CeAl 3 Ni 2 , CeAl 4 Ni, CeAl 5 Ni 2 , CeAlNi, and Ce 4 Al 23 Ni 6 are located in this region.However, there are small protruding areas at the Ni end, and CeAlNi 4 is present around it.The density contours of the Al-Ni-Ce system are shown in Figure 3b.It can be seen that the density of the Al-Ni-Ce ternary system is significantly affected by Al, and the lower the Al content, the greater the density.Combining the two parameters of volume and density, Al atoms have an important contribution to the lightening of Al-Ni-Ce IMCs.
presents the relationship between the volume and density of Al-Ni IMCs and the Al, Ni, and Ce contents.The volume contour distribution of the Al-Ni system is shown in Figure 3a.It can be seen that the larger region of Al-Ni-Ce is ma concentrated on the side of Al and Ce, and CeAl2Ni, CeAl3Ni2, CeAl4Ni, CeAl5Ni2, CeA and Ce4Al23Ni6 are located in this region.However, there are small protruding areas at Ni end, and CeAlNi4 is present around it.The density contours of the Al-Ni-Ce system shown in Figure 3b.It can be seen that the density of the Al-Ni-Ce ternary system is nificantly affected by Al, and the lower the Al content, the greater the density.Combin the two parameters of volume and density, Al atoms have an important contribution the lightening of Al-Ni-Ce IMCs.
where  denotes the total energy;    ,    , and    denote the energie individual atoms in the face-centered cubic (FCC) structure of Al, Ni, and Ce monom

Figure 4 .
Figure 4.The calculated enthalpies of formation and cohesion energy plotted as a function of mole fraction of Ce for Al-Ni-Ce IMCs.(a) Enthalpy of formation; (b) cohesion energy.

Figure 4 .
Figure 4.The calculated enthalpies of formation and cohesion energy plotted as a function of mole fraction of Ce for Al-Ni-Ce IMCs.(a) Enthalpy of formation; (b) cohesion energy.

3. 3 .
Elastic Constants and Elastic Modulus 3.3.1.Elastic ConstantsThe elastic stiffness constant of the crystal, C ij , reflects the deformation capacity of the crystal under stress conditions, where C 11 and C 12 denote transverse deformation and C 44 denotes longitudinal deformation, and the larger the value, the better the mechanical properties.The elastic stiffness constant C ij can be expressed as[29]

3
Ni 2 > Ce 4 Al 23 Ni 6 > CeAl 4 Ni > CeAlNi 4 > CeAl 5 Ni 2 > CeAlNi > CeAl 2 Ni.This indicates that CeAl 3 Ni 2 has the strongest resistance to plastic deformation of Al-Ni-Ce IMCs, while CeAl 2 Ni has the weakest resistance to plastic deformation.The bulk modulus reflects the ability of the material to resist uniform deformation.Calculations show that the order of magnitude of the ability of Al-Ni-Ce IMCs to resist uniform deformation is CeAlNi 4 > CeAl 3 Ni 2 > CeAl 5 Ni 2 > Ce 4 Al 23 Ni 6 > CeAlNi > CeAl 4 Ni > CeAl 2 Ni.Li et al.

3
Ni 2 > Ce 4 Al 23 Ni 6 > CeAl 4 Ni > CeAlNi 4 > CeAl 5 Ni 2 > CeAlNi > CeAl 2 Ni.This shows that CeAl 3 Ni 2 is the hardest of the Al-Ni-Ce intermetallic compounds, while CeAl 2 Ni has the lowest hardness.The magnitude of Poisson's ratio reflects the bonding force, which ranges from 0.25 to 0.45 for metallic materials [37].The Poisson's ratios of the Al-Ni-Ce IMCs CeAl 3 Ni 2 , CeAl 4 Ni, and Ce 4 Al 23 Ni 6 are all less than 0.25, which suggests that ionic bonding plays a dominant role in their bonding.The Poisson's ratios of CeAl 2 Ni, CeAl 5 Ni 2 , CeAlNi, and CeAlNi 4 are all greater than 0.25, indicating that IMCs play a major role in their bonding.
Al m Ni n Ce l total denotes the total energy; E Al solid , E Ni solid , and E Ce solid denote the energies of individual atoms in the face-centered cubic (FCC) structure of Al, Ni, and Ce monomers; E Al atom , E Ni atom , and E Ce atom denote the energies of individual atoms in the free states of Al, Ni, and Ce; and m, n, and l denote the number of atoms in the Al-Ni-Ce IMCs that correspond to Al, Ni, and Ce, respectively.

Table 2 .
Enthalpy of formation and cohesion energy of Al-Ni-Ce intermetallic compounds.

Table 3 .
The elastic constants of Al-Ni-Ce intermetallic compounds.

Table 3 .
The elastic constants of Al-Ni-Ce intermetallic compounds.
* Ce 4 Al 23 Ni 6 belongs to the monoclinic crystal system, and it is still necessary to add several elastic constants, which are C 15 = 2.84; C 25 = 9.53; C 35 = −6.09;and C 45 = 9.97.

Table 5 .
The elastic modulus anisotropy A X of Al-Ni-Ce intermetallic compounds.
3.The elastic modulus results show that CeAl 3 Ni 2 has the highest Young's modulus, CeAlNi 4 has the highest bulk modulus, and the Poisson's ratios of CeAlNi and CeAl 2 Ni are the highest.4. The analysis of the Al-Ni-Ce IMCs modulus of elasticity showed that CeAl 3 Ni 2 , CeAl 4 Ni, CeAl 5 Ni 2 , and Ce 4 Al 23 Ni 6 exhibited brittleness, while CeAlNi, CeAlNi 4 , and CeAl 2 Ni exhibited toughness.5.The results of elastic anisotropy analysis showed that the Young's modulus of CeAl 4 Ni and CeAlNi 4 was the same, the bulk modulus of Ce 4 Al 23 Ni 6 and CeAl 4 Ni was the same, and the Poisson's ratio of CeAlNi 4 and CeAl 2 Ni was the same.6.