Influence of Mg Doping on the Structure and Mechanical Properties of Al2Cu Precipitated Phase by First-Principles Calculations

Al-Cu-Mg high-strength alloys are widely used in industrial production because of their excellent mechanical performance and good machining properties. In this study, first-principles calculations based on density functional theory were carried out to investigate the influence of Mg doping on the structural stability and mechanical properties of the Al2Cu (θ) precipitated phase in Al-Cu-Mg alloys. The results show that the structural stability, electronic structure, bulk modulus, mechanical anisotropy, and thermodynamic properties of the precipitated Al2CuMgX phase change with the concentration of Mg doping (X = 2, 4, 6, and 8). The cohesive energy calculation and electronic structure analysis show that Al2CuMg6 has a high structural stability. The criterion based on elastic constants indicates that Al2CuMg2, Al2CuMg4, and Al2CuMg8 have a brittle tendency and show strong anisotropy of mechanical properties, while Al2CuMg6 shows better comprehensive mechanical properties. The thermodynamic analysis results based on the quasi-harmonic Debye model show that the Al2CuMg6 precipitated phase has good stability at high temperatures and pressure.


Introduction
As a kind of high-performance material, Al-Cu-Mg high-strength aluminum alloy has been widely used in aerospace, weapon equipment, new energy vehicles, and other fields due to its excellent mechanical properties and good machinability.After the supersaturated solid solution is formed by solution treatment or quenching, the supersaturated solute atoms are diffused and precipitated in the aging process and form the precipitated phase.For Al-Cu-Mg alloys, different types of strengthening phases can be precipitated in the aging process according to the alloy components, such as θ and θ ′ phases (Al 2 Cu), θ ′′ phase (Al 3 Cu 2 ), S phases (Al 2 CuMg), and Ω phases (Al 2 CuMgAg) [1][2][3][4][5].These precipitated phases have an important effect on the properties of the Al-Cu-Mg alloy.
Previous researchers have studied the evolution of doping and precipitated phases in Al-Cu-Mg alloys.For example, Cai et al. [6] investigated the effect of Ag on the precipitation behavior of the Ω phase in Al-Cu-Mg alloys.It was found that Ag could promote the precipitation of the Ω phase and increase the activation energy of the θ ′ phase, thus inhibiting the precipitation of the θ ′ phase.Xiao et al. [7] analyzed the effect of different concentrations of Cu and Mg doping on the aging behavior of Al-Cu-Mg alloys.The results indicate that different types of precipitated phases can be obtained by adjusting the Cu and Mg concentrations.Zhu et al. [8] investigated the formation process of the S precipitated phase in Al-Cu-Mg alloys using the first-principles method.The results show that substituting Zr atoms for Al atoms on the (021) surface can promote the precipitation of the S phase and enhance its bulk modulus and shear modulus.Therefore, elemental doping has an important effect on the structural stability and mechanical properties of Al-Cu-Mg alloys.
Al 2 Cu (θ phase) is one of the common precipitated phases in Al-Cu-Mg alloys [9,10], and its structure and properties were first described by Friaul et al. [11].They found that after solid solution treatment, the θ phase is desaturated and precipitated from the matrix during aging to form clusters of solute atoms, leading to the hardening of the alloy.Liu et al. [12] showed that Al 2 Cu could effectively improve the corrosion resistance of materials.It is also beneficial for the recrystallization of the microstructures and improves the fracture toughness of the material.Chen et al. [13] showed that doping elements would affect the nucleation, stability, and properties of the θ phase.They found that both co-and semico-grid interfaces can be formed between Al 2 Cu, and solute elements show different bias behavior at the Al 2 Cu interface.In addition, computational simulations were carried out to investigate the properties of the Al 2 Cu precipitated phase.Perakis et al. [14] investigated the interfacial structures of the θ phase in Al-Cu alloys at high temperatures after long aging using transmission electron microscopy and first-principles calculations.It was found that the lateral growth of the θ phase is influenced by several major mechanisms.Atomic diffusion of the θ phase transition can occur either through a gap mechanism based on additional Al atoms or directly through Al atom diffusion, both of which are quite different from previous theories based on direct diffusion of Cu atoms.Shao et al. [15] studied the effect of Nb, Ge, and Ni doping on the structural, stability, mechanical, and electronic properties of the Al 2 Cu precipitated phase by using first-principles calculations.The results show that doped elements improve the stiffness of the precipitated phase but reduce its ductility.In addition, doping elements can greatly reduce the anisotropy of the precipitated phase and significantly improve its binding strength and stability.
Due to the constraints of the variety of precipitated phases, complex structure, and difficulty of characterization, it is not easy to fully understand the relationship between the structure and properties of precipitated phases through experiments and characterization methods, which restricts the further development of high-performance alloys.The atomistic simulation method, such as first-principles calculations, provides a new idea for the composition design of alloys.The calculation method has been proven to be effective in studying the influence of doping elements on the precipitated phase [16][17][18][19] and can be used to optimize the structure and study the properties [20].In this paper, the effect of Mg doping on the structural stability, bulk modulus, mechanical properties, and electronic structure of the Al 2 Cu (θ phase) have been studied by first-principles calculations, which provides theoretical guidance for composition design and property optimization of the high-performance Al-Cu-Mg alloy.

Materials and Methods
First-principles calculations were performed using the projector-augmented plane wave (PAW) method in the Vienna ab initio simulation package (VASP) code [21,22].The exchange-correlation term is treated with the Perdew-Burke-Ernzerhof (PBE) function [23].At this point, the cutoff energy of the plane wave is set to 700 eV to ensure the accuracy of the results.The supercell of 2 × 2 × 2 (Al 64 Cu 32 ) was built to achieve low dopant element concentration.The integration in the inverse space was performed using a Monkhorst-Pack grid centered at the Γ-point with a K-spacing of 0.1 Å −1 .The convergence criterion for the total energy is 10 −8 eV.All the related structures have been relaxed using the conjugate gradient algorithm.The structural optimization was performed until all residual forces acting on each atom's residual forces are less than 1 meV/Å.The electronic configurations of the potentials utilized in the analysis are Al (3s 2 3p 1 ), Cu (3d 10 4s 1 ), and Mg (2s 2 2p 6 ).
The experimentally determined θ phase has a tetragonal structure with space group I4/mcm (No. 140) and t/12 symmetry.Its unit cell contains 4 formula units, as shown in Figure 1a.The lattice constants are a = 6.063Å and c = 4.872 Å. Al and Cu atoms are located on θ (0.1541, 0.6541, 0) and 4a (0, 0, 0.25) sites, respectively, and Figure 1b is the supercell model.The lattice constants obtained after the crystal structure optimiza-tion are shown in Table 1, which agrees well with the previous experimental results [24].In order to study the preferential site occupation of Mg addition in the Al 2 Cu phase, the formation energies of the impurities-doped Al 2 Cu phases were calculated through [25].Lower formation energy corresponds to a more preferential site occupation tendency.When the doping Mg atom occupied the Al and Cu site, the formation energy was 4. on θ (0.1541, 0.6541, 0) and 4a (0, 0, 0.25) sites, respectively, and Figure 1b is the supercell model.The lattice constants obtained after the crystal structure optimization are shown in Table 1, which agrees well with the previous experimental results [24].In order to study the preferential site occupation of Mg addition in the Al2Cu phase, the formation energies of the impurities-doped Al2Cu phases were calculated through Δ E x E 1 x E [25].Lower formation energy corresponds to a more preferential site occupation tendency.When the doping Mg atom occupied the Al and Cu site, the formation energy was 4.

Structural Properties
The structural stability of the precipitated phase is closely associated with its cohesive energy (Ecoh), which is the energy released during the process of atoms combining from a free state into a crystalline compound [28].The negative value of cohesive energy indicates

Structural Properties
The structural stability of the precipitated phase is closely associated with its cohesive energy (E coh ), which is the energy released during the process of atoms combining from a free state into a crystalline compound [28].The negative value of cohesive energy indicates that the structure can exist stably, thus indicating the structural stability of the Al 2 Cu precipitated phase in this study.The larger the absolute value of cohesive energy, the more stable the structure.E coh can be obtained by Equation ( 1).
where E coh is the cohesive energy, E AB total is the energy per number of atoms for each intermetallic compound, E A atom and E B atom are the energy of free atoms A and B, and X A represents the atomic fraction of A in compounds.From the calculation results shown in Table 1, it can be seen that the cohesive energy of the Al 2 Cu phase is −3.82 eV/atom, which agrees with the previous calculation results [27].In addition, the value of cohesive energy of In order to investigate the dynamical stability of the crystal of Al 2 CuMg X , the phonon spectra are calculated by first-principles calculations.The model established in this work is a doping model, and the symmetry corresponding to different concentration gradients is not the same, and the symmetry is related to the acoustic branch and the optical branch of the phonon spectrum.The abscissa of the phonon spectrum is taken as the point of high symmetry, and the ordinate is the frequency of the phonon spectrum.The acoustic branch represents the overall vibration of the primitive cell, while the optical branch represents the relative vibration between atoms within the primitive cell.Based on the determined lattice constants, the phonon dispersion for the 2 × 2 × 2 supercell of Al 2 CuMg X is shown in Figure 2. Note that the phonon-related properties in this work are calculated based on the density functional perturbation theory.Figure 2 shows that no negative vibrational modes appear in the phonon spectrum, and there are no imaginary frequencies, indicating the dynamical stability of the cubic Al 2 CuMg X .
Materials 2023, 16, x FOR PEER REVIEW 4 of 13 that the structure can exist stably, thus indicating the structural stability of the Al2Cu precipitated phase in this study.The larger the absolute value of cohesive energy, the more stable the structure.Ecoh can be obtained by Equation (1).
where Ecoh is the cohesive energy, E is the energy per number of atoms for each intermetallic compound, E and E are the energy of free atoms A and B, and XA represents the atomic fraction of A in compounds.From the calculation results shown in Table 1, it can be seen that the cohesive energy of the Al2Cu phase is −3.82 eV/atom, which agrees with the previous calculation results [27].In addition, the value of cohesive energy of Al2Cu doped with different Mg concentrations is negative and lower than that of the Al2Cu phase, which means that the doped Al2Cu phases can exist stably, and Mg can improve the stability of the Al2Cu phase.In addition, the absolute values of Ecoh of Al2CuMg6 and Al2CuMg8 are greater than that of Al2CuMg2 and Al2CuMg4, indicating the configurations of Al2CuMg6 and Al2CuMg8 are more stable than the other two configurations.
In order to investigate the dynamical stability of the crystal of Al2CuMgX, the phonon spectra are calculated by first-principles calculations.The model established in this work is a doping model, and the symmetry corresponding to different concentration gradients is not the same, and the symmetry is related to the acoustic branch and the optical branch of the phonon spectrum.The abscissa of the phonon spectrum is taken as the point of high symmetry, and the ordinate is the frequency of the phonon spectrum.The acoustic branch represents the overall vibration of the primitive cell, while the optical branch represents the relative vibration between atoms within the primitive cell.Based on the determined lattice constants, the phonon dispersion for the 2 × 2 × 2 supercell of Al2CuMgX is shown in Figure 2. Note that the phonon-related properties in this work are calculated based on the density functional perturbation theory.Figure 2 shows that no negative vibrational modes appear in the phonon spectrum, and there are no imaginary frequencies, indicating the dynamical stability of the cubic Al2CuMgX.

Electronic Structure
In order to reveal the bonding characteristics of Al2CuMgX and further understand the structural stability differences of the three phases, their partial density of states (PDOS)

Electronic Structure
In order to reveal the bonding characteristics of Al 2 CuMg X and further understand the structural stability differences of the three phases, their partial density of states (PDOS) and total density of states (DOS) were calculated, respectively, as shown in Figures 3 and 4. As shown in Figure 3, there is a peak in the total density of states, located in the range of −6 to −2 eV, which is mainly caused by the 3d orbitals of Cu, and the hybridization between electrons is not very obvious.The free electrons at the Fermi level mainly come from the p-orbitals of Al 3p.
Materials 2023, 16, x FOR PEER REVIEW 5 of 13 and total density of states (DOS) were calculated, respectively, as shown in Figures 3 and  4. As shown in Figure 3, there is a peak in the total density of states, located in the range of −6 to −2 eV, which is mainly caused by the 3d orbitals of Cu, and the hybridization between electrons is not very obvious.The free electrons at the Fermi level mainly come from the p-orbitals of Al 3p.The spin-resolved total and density of states (DOS) of the Mg-doped Al2Cu phases were calculated by the tetrahedron method with Blöchl corrections.The DOS for all the cases were calculated by their equilibrium lattice constants and optimized internal parameters.The Fermi level (Ef) was set as zero and was used as the reference.The valence electrons of the orbital electrons contributing to the calculations are Al (3s 2 3p 1 ), Cu (3d 10 4p 1 ), and Mg (2s 2 2p 6 ).In addition, structural stability is closely related to the electrons' orbital contribution.To further interpret the contributions of different orbitals to the DOS of the Al2Cu precipitated phases, the PDOS at atomic resolution is presented in Figure 4.The common feature and total density of states (DOS) were calculated, respectively, as shown in Figures 3 and  4. As shown in Figure 3, there is a peak in the total density of states, located in the range of −6 to −2 eV, which is mainly caused by the 3d orbitals of Cu, and the hybridization between electrons is not very obvious.The free electrons at the Fermi level mainly come from the p-orbitals of Al 3p.The spin-resolved total and density of states (DOS) of the Mg-doped Al2Cu phases were calculated by the tetrahedron method with Blöchl corrections.The DOS for all the cases were calculated by their equilibrium lattice constants and optimized internal parameters.The Fermi level (Ef) was set as zero and was used as the reference.The valence electrons of the orbital electrons contributing to the calculations are Al (3s 2 3p 1 ), Cu (3d 10 4p 1 ), and Mg (2s 2 2p 6 ).In addition, structural stability is closely related to the electrons' orbital contribution.To further interpret the contributions of different orbitals to the DOS of the Al2Cu precipitated phases, the PDOS at atomic resolution is presented in Figure 4.The common feature The spin-resolved total and density of states (DOS) of the Mg-doped Al 2 Cu phases were calculated by the tetrahedron method with Blöchl corrections.The DOS for all the cases were calculated by their equilibrium lattice constants and optimized internal parameters.The Fermi level (E f ) was set as zero and was used as the reference.The valence electrons of the orbital electrons contributing to the calculations are Al (3s 2 3p 1 ), Cu (3d 10 4p 1 ), and Mg (2s 2 2p 6 ).
In addition, structural stability is closely related to the electrons' orbital contribution.To further interpret the contributions of different orbitals to the DOS of the Al 2 Cu precipitated phases, the PDOS at atomic resolution is presented in Figure 4.The common feature of the four configurations is that there are several characteristic peaks in the DOS.
The width of the peak is mainly between −6 and −2 eV, and the peak is relatively narrow.The peak values are 403.26,372.65, 364.42, and 387.54 states/eV, respectively.This is a typical characteristic of transition metals caused by the d electron orbits of Cu, indicating that the d electrons are relatively localized.The binding peaks between −6 and −2 eV principally originate from the contribution of Al (3s and 3p) and Mg (2s and 2p) orbits.The hybridization between the electrons in all precipitated phases is not apparent.The free electrons near the Fermi level mainly come from the contributions of the state electrons of Al (3p) and Mg (2p), and as the Mg concentration increases, more free electrons are generated, leading to a certain degree of hybridization and increasing the stability of the precipitate phase, which is consistent with the results of cohesive energy.

Elastic Properties
Elastic constants (C ij ) are investigated to understand the mechanical and physical properties of the Al 2 Cu precipitated phase with doping Mg atoms.Elastic constants reflect the bond strength between atoms and are dependent on the effect of external forces on the crystal, playing an important role in determining the mechanical properties of materials.There is a quadratic linear relationship between energy and lattice strain in the crystal structure.Based on the energy-strain relationship, a certain elastic constant or a combination of crystal elastic constants can be obtained by applying strains in different directions and sizes to crystals, and finally obtaining the constant flexible matrix of the crystal.Some criteria are utilized to judge the mechanical stability of the precipitated phases based on the elastic constants.For tetragonal structures, the criteria are determined by Equation (2) [29]: The calculated elastic constants of Al 2 CuMg X are listed in Table 2.According to the elastic constants and mechanical stability criteria, all of them satisfy Born's criterion [30], indicating that they can exist stably, which is consistent with the calculation results of cohesive energy.The elastic properties of the precipitated phases were obtained based on the calculations of the elastic constants of the precipitated phases and using the Voigt-Reuss-Hill approximation [31], as shown in Table 3.For the tetragonal crystal system, the bulk modulus (B) and the shear modulus (G) can be calculated by Equations ( 3)-( 10) [32]: The theoretical value of the modulus of elasticity can be approximated by Equations ( 11) and ( 12) [33]: Poisson's ratio (ν) and elastic modulus (E) can be calculated according to Equations ( 13) and ( 14) [33]: High bulk modulus (B) indicates that the material has strong resistance to volume change, while shear modulus (G) represents the crystal's resistance to plastic deformation, both of which are closely related to the hardness of materials.The high values of B and G reflect the strong bond between atoms.In Table 3, it can be seen that with the addition of the Mg element, the bulk modulus, shear modulus, and Young's modulus of Al 2 Cu are all improved.The values of the bulk modulus and shear modulus of Al 2 CuMg 8 are 58.87 and 60.42 GPa, and the values of Al 2 CuMg 6 are 60.42 and 37.59 GPa, respectively, which are higher than the values of the other two configurations of precipitated phases.The results indicate that Al 2 CuMg 6 and Al 2 CuMg 8 could show good compressive and shear resistance properties.In addition, the hardness of the material (Hv) is closely related to the elasticity of the material, which can be predicted by a semi-empirical formula Hv = 2(K 2 •G H ) 0.585 − 3.0, where K = G H /B H and G H is the shear modulus [34].The hardness values of Al 2 CuMg X with different Mg concentrations are 20.35,21.69, 23.58, and 24.09GPa, respectively.
According to Pugh's theory, the ductility of the materials is closely related to the B/G value, where a value higher than 1.75 shows a favorable ductility [35].The results displayed in Table 3 reveal that with the increase in Mg concentration, a gradual transition to a ductile phase becomes apparent.However, only Al 2 CuMg 6 is the ductile phase, and Al 2 CuMg 2 , Al 2 CuMg 4 , and Al 2 CuMg 8 are the brittle phases, which indirectly indicates that concentration changes affect the mechanical properties of the precipitated phase.Furthermore, the value of C 11 -C 12 can also be used to determine the ductility of materials [36].The smaller value means the better ductility of the materials.From Table 3, the value again proves that Al 2 CuMg 6 has the best ductility.In addition, Poisson's ratio (ν) is also an important index to measure the ductility or brittleness of materials.The high ductility is associated with a high value of Poisson's ratio.According to the above criteria, it is believed that Al 2 CuMg 6 is the most likely precipitated phase component to exhibit good ductility.

Elastic Anisotropy
The elastic anisotropy of materials is closely related to the initiation and propagation of micro-cracks in the process of deformation.The universal anisotropy index (A U ) was introduced to describe the elastic anisotropy behavior of the Al 2 CuMg X precipitated phases.A U is calculated by Equations ( 15)-( 17) [37]: where B V and B R are the Voigt and Reuss bulk modulus, respectively.G V and G R are the Voigt and Reuss shear modulus [38].Note that an isotropic material has A U = 0, while the magnitude of any deviation from zero indicates the degree of elastic anisotropy.If the value of A U is equal to zero, the system is isotropic.On the contrary, the more significant the difference from zero, the greater the anisotropy.It can be seen from  To further discuss the mechanical anisotropy of the Al 2 CuMg X precipitated phases in different orientations, we used the spherical coordinate method to give the threedimensional direction of Young's modulus and the two-dimensional projections in the (001), (100), and (010) planes, as shown in Figure 5.The color bar represents the magnitudes of Young's modulus.Based on the directional dependence of the three-dimensional closed surface on the elastic modulus, the degree of elastic anisotropy can be evaluated, and the relationship between S ij and C ij is an inverse matrix.The directional dependence of Young's modulus for a cubic crystal can be obtained by Equations ( 18) and ( 19) [39]: where S ij are the elastic compliance constants and li (i = 1, 2, 3) are the direction cosines [39]: In an isotropic system, the 3D closed surface representing the relation between the elastic modulus and the crystallographic direction should be spherical.The degree of deviation from the sphere reflects the degree of anisotropy.Aniso represents the anisotropy of materials, and crystal anisotropy refers to the periodicity and density of the atomic arrangement along different directions of the lattice, resulting in different physical and chemical properties of crystals in different directions, and aniso uses value 1 as a reference.It can be combined and compared with A U values to evaluate the anisotropy of Al2CuMgX.As shown in Figure 5, Young's modulus of Al2CuMg6 deviates slightly from the sphere with the 3D characterization of crystallographic direction, and the value of aniso is 1.0713, which is the smallest among the four cases, indicating that the configuration has weak elastic anisotropy.As shown in Figure 5a,b,d, Al2CuMg2, Al2CuMg4, and Al2CuMg8 deviate noticeably from the spherical shape.The result indicates that these configurations have strong elastic anisotropy, which is negative to the ductility or toughness of the materials.

Thermodynamic Analysis
The above calculation results show that Al2CuMg6 has the best comprehensive mechanical properties among the four models of the Al2Cu precipitated phase.The thermodynamic properties of Al2CuMg6 were further investigated using the quasi-harmonic Debye model for Gibbs simulation [40,41].The relationship between energy and volume, as well as Poisson's ratio, were obtained through the matrix calculation of Al2CuMg6.The Brich-Murnaghan state equation was used for fitting to obtain the E-V curve [42]: The simulation temperature was set as 0 to 2000 K, and the system pressure was set as 0 to 20 GPa for the first-principles calculations.Although there are some errors in the results of the quasi-harmonic Debye model used in the thermodynamic calculations, they are in the same order of magnitude as the experimental values [43], so the fluctuations In an isotropic system, the 3D closed surface representing the relation between the elastic modulus and the crystallographic direction should be spherical.The degree of deviation from the sphere reflects the degree of anisotropy.Aniso represents the anisotropy of materials, and crystal anisotropy refers to the periodicity and density of the atomic arrangement along different directions of the lattice, resulting in different physical and chemical properties of crystals in different directions, and aniso uses value 1 as a reference.It can be combined and compared with A U values to evaluate the anisotropy of Al 2 CuMg X .As shown in Figure 5, Young's modulus of Al 2 CuMg 6 deviates slightly from the sphere with the 3D characterization of crystallographic direction, and the value of aniso is 1.0713, which is the smallest among the four cases, indicating that the configuration has weak elastic anisotropy.As shown in Figure 5a,b,d, Al 2 CuMg 2 , Al 2 CuMg 4 , and Al 2 CuMg 8 deviate noticeably from the spherical shape.The result indicates that these configurations have strong elastic anisotropy, which is negative to the ductility or toughness of the materials.

Thermodynamic Analysis
The above calculation results show that Al 2 CuMg 6 has the best comprehensive mechanical properties among the four models of the Al 2 Cu precipitated phase.The thermodynamic properties of Al 2 CuMg 6 were further investigated using the quasi-harmonic Debye model for Gibbs simulation [40,41].The relationship between energy and volume, as well as Poisson's ratio, were obtained through the matrix calculation of Al 2 CuMg 6 .The Brich-Murnaghan state equation was used for fitting to obtain the E-V curve [42]: The simulation temperature was set as 0 to 2000 K, and the system pressure was set as 0 to 20 GPa for the first-principles calculations.Although there are some errors in the results of the quasi-harmonic Debye model used in the thermodynamic calculations, they are in the same order of magnitude as the experimental values [43], so the fluctuations and trends of the thermodynamic properties of the precipitated phase can be predicted approximately.The E-V curve of the Al 2 CuMg 6 precipitated phase was plotted in Figure 6.The lowest energy point was determined according to the fitting of the E-V curve.The results show that the fitting degree of the E-V curve is high, and the midpoint of volume is the lowest point of energy, indicating that the calculation method and parameters are accurate, which is the premise of thermodynamic analysis.
Materials 2023, 16, x FOR PEER REVIEW 10 of 1 and trends of the thermodynamic properties of the precipitated phase can be predicted approximately.The E-V curve of the Al2CuMg6 precipitated phase was plotted in Figure 6.The lowest energy point was determined according to the fitting of the E-V curve.The results show that the fitting degree of the E-V curve is high, and the midpoint of volume is the lowest point of energy, indicating that the calculation method and parameters are accurate, which is the premise of thermodynamic analysis.Figure 7 shows the variation of the thermal expansion coefficient of Al2CuMg6 with temperature under different pressures.It is found that temperature and pressure have important effects on the thermal expansion coefficient of the precipitated phase structure Before 500 K, the coefficient of thermal expansion changes rapidly with temperature There is an obvious slope inflection point in the curve, and it basically occurs in the same temperature range.When the temperature rises above 500 K, the change rate of therma expansion coefficient decreases continuously.When the pressure is at a relatively low level, the coefficient of thermal expansion is sensitive to the trend of temperature.How ever, when the pressure increases to more than 6 GPa, the thermal expansion coefficien tends to be flat with the increase in temperature, and the greater the pressure, the smalle the thermal expansion coefficient.Figure 7 shows the variation of the thermal expansion coefficient of Al 2 CuMg 6 with temperature under different pressures.It is found that temperature and pressure have important effects on the thermal expansion coefficient of the precipitated phase structure.Before 500 K, the coefficient of thermal expansion changes rapidly with temperature.There is an obvious slope inflection point in the curve, and it basically occurs in the same temperature range.When the temperature rises above 500 K, the change rate of thermal expansion coefficient decreases continuously.When the pressure is at a relatively low level, the coefficient of thermal expansion is sensitive to the trend of temperature.However, when the pressure increases to more than 6 GPa, the thermal expansion coefficient tends to be flat with the increase in temperature, and the greater the pressure, the smaller the thermal expansion coefficient.
Figure 8 shows the isobaric heat capacity (C p ) and isothermal heat capacity (C v ) of Al 2 CuMg 6 under different pressures as a function of temperature.As can be seen from the result, the values of C v and C p of the precipitated phase are basically the same in the low-temperature region below 300 K, and their values increase very fast due to the quantum effect at low temperatures.However, as the temperature continues to rise, the two curves show different trends, and the value of C p increases continuously with the increase in temperature, while C v gradually increases slowly, and finally tends to a saturation value.In the high-temperature region, the value of C p is higher than C v , but the increasing speed also slows down with the temperature, which conforms to the Dulong-Petit law, indicating that the Al 2 CuMg 6 precipitated phase has good stability at high temperature and high pressure.
temperature range.When the temperature rises above 500 K, the change rate of thermal expansion coefficient decreases continuously.When the pressure is at a relatively low level, the coefficient of thermal expansion is sensitive to the trend of temperature.However, when the pressure increases to more than 6 GPa, the thermal expansion coefficient tends to be flat with the increase in temperature, and the greater the pressure, the smaller the thermal expansion coefficient.Figure 8 shows the isobaric heat capacity (Cp) and isothermal heat capacity (Cv) of Al2CuMg6 under different pressures as a function of temperature.As can be seen from the result, the values of Cv and Cp of the precipitated phase are basically the same in the lowtemperature region below 300 K, and their values increase very fast due to the quantum effect at low temperatures.However, as the temperature continues to rise, the two curves show different trends, and the value of Cp increases continuously with the increase in temperature, while Cv gradually increases slowly, and finally tends to a saturation value.In the high-temperature region, the value of Cp is higher than Cv, but the increasing speed also slows down with the temperature, which conforms to the Dulong-Petit law, indicating that the Al2CuMg6 precipitated phase has good stability at high temperature and high pressure.

Conclusions
First-principles calculations were carried out to investigate the structural stability, mechanical properties, and thermodynamic properties of the Al 2 Cu precipitated phase with different concentrations of Mg doping.The main conclusions are summarized as follows: (4) The thermodynamic analysis results based on the quasi-harmonic Debye model show that the volume thermal expansion coefficient and heat capacity of Al 2 CuMg 6 decrease with an increase in pressure and increase with an increase in temperature, which accords with the Dulong-Petit law, indicating that the Al 2 CuMg 6 precipitated phase has good stability at high temperature and high pressure.
43 eV and 8.62 eV per atom, respectively.It indicates that the doping Mg atom prefers to occupy the Al site in the Al 2 Cu structure.The change in Mg concentration in the Al 2 Cu precipitated phase is realized by changing the number of Mg atoms doped in the model.In this study, the number of Mg atoms doped in Al 2 CuMg X is set as X = 2, 4, 6, 8.The different concentration models are shown in Figure 1c-f, where two Mg atoms are defined as a set to obtain different doping concentrations from 2 to 8 Mg atoms, and the doped Mg atoms are distributed discretely between each other.Materials 2023, 16, x FOR PEER REVIEW 3 of 13 43 eV and 8.62 eV per atom, respectively.It indicates that the doping Mg atom prefers to occupy the Al site in the Al2Cu structure.The change in Mg concentration in the Al2Cu precipitated phase is realized by changing the number of Mg atoms doped in the model.In this study, the number of Mg atoms doped in Al2CuMgX is set as X = 2, 4, 6, 8.The different concentration models are shown in Figure 1c-f, where two Mg atoms are defined as a set to obtain different doping concentrations from 2 to 8 Mg atoms, and the doped Mg atoms are distributed discretely between each other.
Al 2 Cu doped with different Mg concentrations is negative and lower than that of the Al 2 Cu phase, which means that the doped Al 2 Cu phases can exist stably, and Mg can improve the stability of the Al 2 Cu phase.In addition, the absolute values of E coh of Al 2 CuMg 6 and Al 2 CuMg 8 are greater than that of Al 2 CuMg 2 and Al 2 CuMg 4 , indicating the configurations of Al 2 CuMg 6 and Al 2 CuMg 8 are more stable than the other two configurations.

Figure 3 .
Figure 3.Total and partial density of state of Al2Cu.The dashed line is the Fermi level.

Figure 3 .
Figure 3.Total and partial density of state of Al 2 Cu.The dashed line is the Fermi level.

Figure 3 .
Figure 3.Total and partial density of state of Al2Cu.The dashed line is the Fermi level.

Figure 6 .
Figure 6.The relationship between energy and volume (E-V) curve of Al2CuMg6 phase.

Figure 7 .
Figure 7. Thermal expansion coefficient of Al2CuMg6 phase as a function of temperature at differen

Figure 6 .
Figure 6.The relationship between energy and volume (E-V) curve of Al 2 CuMg 6 phase.

Figure 7 .
Figure 7. Thermal expansion coefficient of Al2CuMg6 phase as a function of temperature at different pressures.

Figure 7 .
Figure 7. Thermal expansion coefficient of Al 2 CuMg 6 phase as a function of temperature at different pressures.

Figure 8 .
Figure 8.Heat capacity of Al2CuMg6 phase as a function of temperature at different pressures.
calculations were carried out to investigate the structural stability, mechanical properties, and thermodynamic properties of the Al2Cu precipitated phase with different concentrations of Mg doping.The main conclusions are summarized as follows: (1) The calculations of cohesive energy and phonons show that the Al2CuMgX precipitated phase with different Mg doping can exist stably, and Al2CuMg6 shows the highest structural stability.(2) The results of DOS analysis show that the changes in concentration have similar contributions to the electronic structure orbitals and hybridization effects of the precipitates, and the peaks are caused by Cu(3d); the electron hybridization of the precipitated phase is not obvious.(3) Based on the Voigt-Reuss-Hill approximation criterion, the elastic constants of Al2CuMgX were calculated.C12-C44, B/G, and other criteria show that Al2CuMg6 has better comprehensive mechanical properties compared to the other three precipitated phases.Mechanical anisotropy analysis shows that the anisotropy of Al2CuMg6 is the least significant among the four Mg doping precipitated configurations.(4) The thermodynamic analysis results based on the quasi-harmonic Debye model show that the volume thermal expansion coefficient and heat capacity of Al2CuMg6 decrease with an increase in pressure and increase with an increase in temperature, which accords with the Dulong-Petit law, indicating that the Al2CuMg6 precipitated

Figure 8 .
Figure 8.Heat capacity of Al 2 CuMg 6 phase as a function of temperature at different pressures.

( 1 )
The calculations of cohesive energy and phonons show that the Al 2 CuMg X precipitated phase with different Mg doping can exist stably, and Al 2 CuMg 6 shows the highest structural stability.(2) The results of DOS analysis show that the changes in concentration have similar contributions to the electronic structure orbitals and hybridization effects of the precipitates, and the peaks are caused by Cu(3d); the electron hybridization of the precipitated phase is not obvious.(3) Based on the Voigt-Reuss-Hill approximation criterion, the elastic constants of Al 2 CuMg X were calculated.C 12 -C 44 , B/G, and other criteria show that Al 2 CuMg 6 has better comprehensive mechanical properties compared to the other three precipitated phases.Mechanical anisotropy analysis shows that the anisotropy of Al 2 CuMg 6 is the least significant among the four Mg doping precipitated configurations.

Table 2 .
The calculated elastic constants (C ij ) of the precipitated phases (GPa).

Table 3 .
The calculated mechanical properties of Al 2 CuMg X .
Table 4 that the calculated A U of Al 2 CuMg 6 is 0.24, which is the lowest among the four configurations, implying the best elastic isotropy.The A U value of Al 2 CuMg 8 is the highest (0.45), indicating that the Al 2 CuMg 8 exhibits strong anisotropy and is prone to generating micro-cracks.

Table 4 .
The calculated global anisotropy index (A U ) and partial anisotropy index (A G and A B ).