Theoretical Study on Structural Stability and Elastic Properties of Fe25Cr25Ni25TixAl(25-x) Multi-Principal Element Alloys

Material genetic engineering studies the relationship between the composition, microstructure, and properties of materials. By adjusting the atomic composition, structure, or configuration of the material and combining different processes, new materials with target properties obtained. In this paper, the design, and properties of the ordered phases in Fe25Cr25Ni25TixAl(25-x) (subscript represents the atomic percentage) multi-principal element alloys are studied. By adjusting the percentages of Ti and Al atoms, the effect of the atomic percentage content on ordered phases’ structural stability in multi-principal element alloys are studied. Thermodynamic analysis predicted the composition phase and percentage of the alloy. Formation heat, binding energy, and elastic constants confirmed the structural stability and provide a theoretical basis for designing alloys with target properties. The results showed that the disordered BCC A2 phase and the ordered BCC B2 phase are the ductile phases, while the Laves phase is brittle. The research method in this paper is used to design multi-principal element alloys or other various complex materials that meet the target performance.


Introduction
Multi-principal element alloys, also known as high-entropy alloys (HEAs), have been proven to have superior properties such as ultra-high strength, high hardness, hightemperature oxidation resistance, wear resistance, corrosion resistance, and thermal stability, in materials science [1][2][3][4][5][6][7] and have received considerable attention the engineering field. Compared with conventional alloys, multi-principal element alloys are alloys composed of five or more elements in equimolar ratios. The concentration of each element is between 5% and 35%, and no element content exceeds 50% [8,9]. Therefore, the main element alloy is characterized by the synergy of various elements. The design of multi-principal element alloys represents a new strategy for developing unique engineered materials with targeted properties that cannot be achieved with traditional alloy designs that rely on one element and then add other elements to improve the properties of the alloy.
At present, studies on multi-principal element alloys show that the main factor affecting the entropy value is still the structure of the alloy. Secondly, due to the different atom types and sizes of the multi-principal element alloys, a severe lattice distortion is caused, and the shear modulus between the constituent atoms does not match, which may contribute to hardening. The structural classification of a multi-principal element alloys is classified into one or more microstructures composed of disordered solid solutions. Microstructures having one or more ordered phases (or intermetallic compounds) and a disordered solid solution A microstructure consisting of a mixture of intermetallic phases. The most common constituent structures, at present, are the disordered FCC phase and the BCC phase, followed by the ordered BCC B2 phase and the hexagonal Laves phase.
For the study of multi-principal element alloys, the empirical method is to calculate the formation of the disordered solid solution by the Hume-Rothery rule. The factors considered are the atomic size factor, crystal structure, electronegativity, electron concentration, and valence. The conditions for thermodynamic formation are mixed enthalpy, mixing entropy, and melting temperature [10]. The Gibbs free energy is obtained by establishing the equation by the number and concentration of alloying elements without considering the Gibbs free energy of the ordered phase or intermetallic compound. Another method can be used to form entropy and an intermetallic phase to create a crucible. The modeling can distinguish whether the alloy is a single-phase disordered solid solution or contains an ordered phase (intermetallic compound). The third method is the phase diagram calculation, obtaining the phase diagram of the multi-principal element alloys, analyzing the phase structure, and the phase volume fraction. The mechanical properties of the multi-principal element alloys depend, to a large extent, on the composition and microstructure of the alloy. In the case of fixed composition and phase content, the size, shape, and distribution of the phase can also be adjusted to improve the properties of the alloy. Second, factors affecting mechanical properties include defects such as vacancies, dislocations, and microscopic defects such as grain boundaries, as well as the porosity, segregation, cracks, and residual stress. For the study of the hardness of Al x CoCrCuFeNi multi-principal element alloys, the structure changed from single-phase fcc to BCC + FCC and then to single-phase BCC structure, and its hardness increased with the increase of Al content [11]. It has also shown that in the Al x CoCrCuFeNi multi-principal element alloys, since it contains a large electronegativity difference element such as Al/Ni, increasing the Al content leads to chemical order, which acts as a solid solution strengthening and precipitation strengthening, so the hardness increases [12][13][14]. In the study of CrCoFeNiMn multi-principal element alloys, the grain size and strain rate are the main factors affecting the tensile properties. The method of cooking material design has become an effective way to design new materials [15,16].
Based on the research ideas of material genetic engineering, the chemical composition, phase composition, structure, and properties of Fe 25 Cr 25 Ni 25 TixAl (25−x) (x = 0; 6.25; 12.5; 18.75; 25) (Subscript represents atomic percentage) multi-principal element alloys were studied through the thermodynamic phase diagram calculation, first-principles calculation and experimental verification methods. The relationship between the purposes is to design a new multi-principal element alloys, hoping to improve the ductility of the alloy. The thermodynamic phase diagram predicted the phase composition and content of the alloy and the main constituent elements of the phase. Then the stability and elastic properties of the phase structure were predicted by first-principles calculation. Finally, the reliability of the estimate was verified by experiments.

Thermodynamic Methods
For many types of alloys, the advantage of being able to calculate their thermo physical properties through JMatPro software is obvious. We can easily obtain information about the alloy and calculate where the alloy cannot be measured. Because of the lack of obvious and reliable measurement of solidification performance for many alloys, this is a precious advantage. The change of alloy phase content with alloy composition and allowable specification range can be obtained quickly and self-consistently through calculation. Beside, many calculation paths are accompanied by detailed information; otherwise, it is impossible or complicated to determine experimentally. For example, in FeCrNiTi x Al (1-x) multi-principal alloy. Due to the lack of data on thermo physics and physics, the properties of similar alloys can be given through calculations. Some alloys are very sensitive to small composition changes, while other alloys (for example, solid solution alloys) may not be sensitive. The composition specification may also be because the multi-principal alloy is very wide, and many researchers have taken advantage of this alloy.
The thermodynamic calculation software JMatPro used in this article is Public Release Version 7.0.0 provided by Sente Software Ltd. Surrey Technical Center [17][18][19]. Use the nickelbased alloy module of the software to calculate Fe 25 (25-x) (x = 0; 6.25; 12.5; 18.75; 25) (subscript represents atomic percentage). General physical properties of multi-principal element alloys, alloy phases and cooling curves, the setting temperature range are 0-1400 • C. The phase composition and content of the alloy during solidification were studied.

Density Functional Theory
The plane wave ultra-soft pseudopotential method based on density functional theory has been widely applied to the structural design and performance prediction of materials. Based on this method, the first principle modeling of the constituent phases in the alloy is carried out according to the calculation results of the phase diagram. Using the CASTEP module in Materials Studio 8.0 [20], select the utilized density functional PBE in generalized gradient approximation (GGA) [21], select 350 eV for the energy cutoff value, and select k point for 4 × 4 × 4 in the first Brillouin region of the reciprocal space. The volume and atomic position of the supercell structure optimized until the maximum force, maximum displacement, maximum stress, and energy change of each atom are below 0.01 eV/Å, 5.0 × 10 −4 Å, 0.02 GPa, and 5.0 × 10 −6 eV/atom, respectively. The stability and elastic properties of the alloy structure, such as second-order elastic constant, bulk modulus, shear modulus, Young's modulus, and Poisson's ratio were comprehensively calculated.

Thermodynamic Analysis
To study the types and relative contents of various constituent phases of Fe 25 Cr 25 Ni 25 Ti x Al (25−x) (x = 0; 6.25; 12.5; 18.75; 25) (Subscript represents atomic percentage) multi-principal element alloys at different temperatures, and to determine the composition of each equilibrium phase, thermodynamic analysis The software JMat-Pro performs phase diagram calculation on it. To combine the thermodynamic analysis with the order of the structure in the multi-principal element alloys, try to combine the thermodynamics with the first-principles calculations to explore a new multi-principal alloy that meets the target performance.
Among the constituent elements of the Fe 25 Cr 25 Ni 25 Ti x Al (25−x) (subscript represents atomic percentage) system multi-principal element alloys, the Fe, Ni, and Ti elements are the main constituent elements of the conventional metal-based superalloy. The Ni, Ti, and Cr elements contribute to the high-temperature properties of traditional alloys. Due to the large atomic radius of the Al element, it has a great influence on the crystal structure of the alloy system. The crystal structure of the alloy can be adjusted by changing the content of the Al element to achieve the purpose of improving the properties of the alloy. In the Fe 25 Cr 25 Ni 25 Ti x Al (25−x) system multi-principal element alloys, it is desirable to adjust the content of the Ti and Al elements to improve the toughness. Figure 1 shows the thermodynamic calculation results for the alloy. When x = 0, the Fe 25 Cr 25 Ni 25 Al 25 alloy phase diagram shows that the disordered BCC A2 phase and the ordered BCC B2 phase content are 50%, respectively ( Figure 1a). As the Ti content increases (x = 6.25), the Al content reduces, the content of the ordered BCC B2 phase and the disordered BCC A2 phase increased simultaneously, and the Laves phase decreased (Figure 1b). When x = 12.5, the content of the disordered BCC A2 phase changes little, the content of the ordered BCC B2 phase continues to increase, and the Laves phase content continues to decrease ( Figure 1c). As the Ti content continued to increase (x = 18.75), the Al content continued to reduce, the content of the disordered BCC A2 phase and the ordered BCC B2 phase increases, and the Laves phase content decreases, less than 50% (Figure 1d). When x = 25, the phase diagram of the Fe 25 Cr 25 Ni 25 Ti 25 alloy shows that the alloy consists of a disordered BCC A2 phase, an ordered BCC B2 phase, and a Laves phase, a small BCC B2 phase content, and a Laves phase content greater than 75% (Figure 1e). Further analysis showed that with the increase of Ti content and the decrease of Al content, the content of Al phase and B2 phase decreased gradually, the content of A2 phase decreased rapidly, and the content of Laves phase increased gradually. Table 1 shows the mole percent of each phase of the Fe 25 Cr 25 Ni 25 Ti x Al (25-x) alloy at 800 • C. It can be seen from the data in the table that as the Al content in the alloy increases, the contents of the BCC A2 phase and the BCC B2 phase gradually increase. The Laves phase content is gradually reduced. content of the disordered BCC A2 phase and the ordered BCC B2 phase increases, and the Laves phase content decreases, less than 50% (Figure 1d). When x = 25, the phase diagram of the Fe25Cr25Ni25Ti25 alloy shows that the alloy consists of a disordered BCC A2 phase, an ordered BCC B2 phase, and a Laves phase, a small BCC B2 phase content, and a Laves phase content greater than 75% (Figure 1e). Further analysis showed that with the increase of Ti content and the decrease of Al content, the content of Al phase and B2 phase decreased gradually, the content of A2 phase decreased rapidly, and the content of Laves phase increased gradually. Table 1 shows the mole percent of each phase of the Fe25Cr25Ni25TixAl(25-x) alloy at 800 °C. It can be seen from the data in the table that as the Al content in the alloy increases, the contents of the BCC A2 phase and the BCC B2 phase gradually increase. The Laves phase content is gradually reduced.        Figure 2b is the fraction liquid (at)-temperature ( • C) curve of Fe 25 Cr 25 Ni 25 Ti 18.75 Al 6.25 alloy (subscript is atomic percentage). The figure shows that as the temperature increases, the volume fraction of the liquid phase in the alloy changes accordingly. The multi-principal element alloy has a fixed melting temperature, the melting temperature is 1284.85 • C, currently fraction liquid (at): 0.0. It absorbs heat during the melting process. When the temperature rises to 1356.1 • C, fraction Liquid (at): 0.9. Then as the temperature increases, the alloy gradually melts. When the temperature increases to 1382.36 • C, fraction liquid (at): 1.0. Figure 2c shows the element in liquid (wt.%)-temperature curve of the alloy. As the temperature increases, the percentage content of various elements in the alloy liquid phase will change. Figure 2d shows the cooling curve of the alloy. The cooling rate obtained for the cooling curve is 1.0 • C/S. When the cooling time is 423.31 s and the temperature is 1288.53 • C, the cooling rate can be adjusted to obtain cooling curves with different cooling rates. Figure 3a is the fraction solid (at)-temperature ( • C) curve of Fe 25 Cr 25 Ni 25 Ti 6.25 Al 18.75 alloy (subscript is atomic percentage). The figure shows that as the temperature rises, the alloy gradually melts, and the melting point is 1140.41 • C. Fraction solid (at): 1.0. Then as the temperature increases, the solid phase volume fraction gradually decreases. When the temperature increases to 1334.85 • C, fraction solid (at): 0.0. Figure 3b is the fraction liquid (at)-temperature ( • C) curve of Fe 25 Cr 25 Ni 25 Ti 18.75 Al 6.25 alloy (subscript is atomic percentage). The figure shows that as the temperature increases, the volume fraction of the liquid phase in the alloy changes accordingly. The multi-principal element alloy has a fixed melting temperature, which is 1140.41 • C, currently fraction liquid (at): 0.0. It absorbs heat during the melting process. When the temperature rises to 1146.89 • C, fraction liquid (at): 0.16. When the temperature rises to 1289.85 • C, fraction liquid (at): 0.88. Then as the temperature rises, the alloy gradually melts. When the temperature rises to 1334.85 • C, fraction liquid (at): 1.0. As the temperature rises, fraction liquid (at) remains unchanged. Figure 3c shows the element in liquid (wt.%)-temperature curve of the alloy. As the temperature increases, the percentage content of various elements in the alloy liquid phase will change. Figure 3d shows the cooling curve of the alloy. The cooling rate obtained for the cooling curve is 1.0 • C/s. when the cooling time is 516.5 s and the temperature is 1140.41 • C, the cooling rate can be adjusted to obtain cooling curves with different cooling rates.    Figure 3b is the fraction liquid (at)-temperature (°C) curve of Fe25Cr25Ni25Ti18.75Al6.25 alloy (subscript is atomic percentage). The figure shows that as the temperature increases, the volume fraction of the liquid phase in the alloy changes accordingly. The multi-principal element alloy has a fixed melting temperature, which is 1140.41°C, currently fraction liquid (at): 0.0. It absorbs heat during the melting process. When the temperature rises to 1146.89 °C, fraction liquid (at): 0.16. When the temperature rises to 1289.85 °C, fraction liquid (at): 0.88. Then as the temperature rises, the alloy gradually melts. When the temperature rises to 1334.85 °C, fraction liquid (at): 1.0. As the temperature rises, fraction liquid (at) remains unchanged. Figure 3c shows the element in liquid (wt.%)-temperature curve of the alloy. As the temperature increases, the percentage content of various elements in the alloy liquid phase will change. Figure 3d shows the cooling curve of the alloy. The cooling rate obtained for the cooling curve is 1.0 °C/s. when the cooling time is 516.5 s and the temperature is 1140.41 °C, the cooling rate can be adjusted to obtain cooling curves with different cooling rates.    Figure 4b is the fraction liquid (at)-temperature (°C) curve of Fe25Cr25Ni25Ti12.5Al12.5 alloy (subscript is atomic percentage). The figure shows that as the temperature increases, the volume fraction of the liquid phase in the alloy changes accordingly Variety. The multi-principal element alloys has a fixed melting temperature, the melting temperature is 1129.91 °C, at this time fraction liquid (at): 0.0. It absorbs heat during the melting process. When the temperature rises to 1199.85 °C, fraction liquid (at): 0.8. Then as the temperature rises, the alloy gradually melts. When the temperature rises to 1314.85 °C, fraction liquid (at): 1.0. As the temperature rises, the fraction liquid (at) remains unchanged. Figure 4c shows the element in liquid (wt.%)-temperature curve of the alloy. As the temperature increases, the percentage content of various elements in the alloy liquid phase will change. Figure 4d shows the cooling curve of the alloy. The cooling rate obtained for the cooling curve is 1.0 °C/s.   Figure 4b is the fraction liquid (at)-temperature ( • C) curve of Fe 25 Cr 25 Ni 25 Ti 12.5 Al 12.5 alloy (subscript is atomic percentage). The figure shows that as the temperature increases, the volume fraction of the liquid phase in the alloy changes accordingly Variety. The multi-principal element alloys has a fixed melting temperature, the melting temperature is 1129.91 • C, at this time fraction liquid (at): 0.0. It absorbs heat during the melting process. When the temperature rises to 1199.85 • C, fraction liquid (at): 0.8. Then as the temperature rises, the alloy gradually melts. When the temperature rises to 1314.85 • C, fraction liquid (at): 1.0. As the temperature rises, the fraction liquid (at) remains unchanged. Figure 4c shows the element in liquid (wt.%)-temperature curve of the alloy. As the temperature increases, the percentage content of various elements in the alloy liquid phase will change. Figure 4d shows the cooling curve of the alloy. The cooling rate obtained for the cooling curve is    Figure 5b is the fraction liquid (at)-temperature (°C) curve of Fe25Cr25Ni25Ti18.75Al6.25 (subscript is atomic percentage), which shows that as the temperature increases, the volume fraction of the liquid phase in the alloy changes accordingly. The multi-principal element alloy has a fixed melting temperature, which is 1117.17 °C, at this time fraction liquid (at): 0.0. As the temperature rises, the alloy gradually melts. When the temperature rises to 1280.37 °C, fraction liquid (at): 1.0. As the temperature rises, fraction liquid (at) remains unchanged. Figure 5c shows the element in liquid (wt.%)-temperature curve of the alloy. As the temperature increases, the percentage content of various elements in the alloy liquid phase will change. Figure 5d shows the cooling curve of the alloy. The cooling rate of the obtained cooling curve is 1.0 °C/S. When the cooling time is 589.38 s, the temperature is 1117.92 °C. The cooling rate can be adjusted to obtain cooling curves with different cooling rates.   Figure 5b is the fraction liquid (at)-temperature ( • C) curve of Fe 25 Cr 25 Ni 25 Ti 18.75 Al 6.25 (subscript is atomic percentage), which shows that as the temperature increases, the volume fraction of the liquid phase in the alloy changes accordingly. The multi-principal element alloy has a fixed melting temperature, which is 1117.17 • C, at this time fraction liquid (at): 0.0. As the temperature rises, the alloy gradually melts. When the temperature rises to 1280.37 • C, fraction liquid (at): 1.0. As the temperature rises, fraction liquid (at) remains unchanged. Figure 5c shows the element in liquid (wt.%)-temperature curve of the alloy. As the temperature increases, the percentage content of various elements in the alloy liquid phase will change. Figure 5d shows the cooling curve of the alloy. The cooling rate of the obtained cooling curve is 1.0 • C/S. When the cooling time is 589.38 s, the temperature is 1117.92 • C. The cooling rate can be adjusted to obtain cooling curves with different cooling rates. shows that as the temperature increases, the volume fraction of the liquid phase in the alloy changes accordingly. The multi-principal element alloys has a fixed melting temperature, the melting temperature is 1159.85 • C, at this time fraction liquid (at): 0.0. As the temperature rises, the alloy gradually melts. When the temperature rises to 1159.85 • C, fraction liquid (at): 0.28; when the temperature rises to 1344.44 • C, fraction liquid (at): 1.0. As the temperature rises, fraction liquid (at) remains unchanged. Figure 6c shows the element in liquid (wt.%)-temperature curve of the alloy. As the temperature increases, the percentage content of various elements in the alloy liquid phase will change. Figure 6d shows the cooling curve of the alloy. The cooling rate of the obtained cooling curve is 1.0 • C/s. When the cooling time is 561.6 s and the temperature is 1150.53 • C, the cooling rate can be adjusted to obtain cooling curves with different cooling rates.   Figure 6b is the fraction liquid (at)-temperature (°C) curve of Fe25Cr25Ni25Ti25 (subscript as atomic percentage). The figure shows that as the temperature increases, the volume fraction of the liquid phase in the alloy changes accordingly. The multi-principal element alloys has a fixed melting temperature, the melting temperature is 1159.85 °C, at this time fraction liquid (at): 0.0. As the temperature rises, the alloy gradually melts. When the temperature rises to 1159.85 °C, fraction liquid (at): 0.28; when the temperature rises to 1344.44 °C, fraction liquid (at): 1.0. As the temperature rises, fraction liquid (at) remains unchanged. Figure 6c shows the element in liquid (wt.%)-temperature curve of the alloy. As the temperature increases, the percentage content of various elements in the alloy liquid phase will change. Figure 6d shows the cooling curve of the alloy. The cooling rate of the obtained cooling curve is 1.0 °C/s. When the cooling time is 561.6 s and the temperature is 1150.53 °C, the cooling rate can be adjusted to obtain cooling curves with different cooling rates.

Enthalpy of Mixing ΔHmix, Entropy of Mixing ΔSmix and Gibbs Free Energy Calculation of Alloy
In thermodynamics, the entropy value represents the degree of chaos. The greater the degree of disorder in a system, the greater the entropy value, as shown in Formula (1) [22]: R: Represents the gas constant; : Represents the percentage of member i.
In multi-principal alloys, the mixing entropy plays a significant role in the formation of a solid solution, and the mixing enthalpy factor cannot be ignored. The formula for mixing enthalpy is shown in Formulas (2) and (3) [23]:

Enthalpy of Mixing ∆Hmix, Entropy of Mixing ∆Smix and Gibbs Free Energy Calculation of Alloy
In thermodynamics, the entropy value represents the degree of chaos. The greater the degree of disorder in a system, the greater the entropy value, as shown in Formula (1) [22]: R: Represents the gas constant; C i : Represents the percentage of member i.
In multi-principal alloys, the mixing entropy plays a significant role in the formation of a solid solution, and the mixing enthalpy factor cannot be ignored. The formula for mixing enthalpy is shown in Formulas (2) and (3) [23]: C i : The atomic percentage of the i component; C j : The atomic percentage of the j component; ∆H mix AB binary system mixing enthalpy. Enthalpy of mixing, also known as the heat of mixing, is one of the mixing functions, which refers to the linear relationship between the alloy's transition temperature and the absolute value of the alloy's mixing enthalpy. The mixing entropy and the mixing enthalpy are in a competitive position. Regardless of whether the mixing enthalpy is positive or negative, the increase of the mixing entropy will inevitably lead to a decrease in the system's free energy, thereby improving the stability of the system. This article uses JMat-Pro software to predict the mixing enthalpy and mixing entropy of Fe 25 Cr 25 Ni 25 Ti x Al (25−x) multi-principal element alloys, as shown in Figure 7. It can be observed from the figure that as the temperature increases, the mixing entropy and mixing enthalpy of the alloy increase, indicating that the system is more stable. The greater the absolute value of the mixing enthalpy, the higher the alloy transition temperature. The absolute value of the mixing enthalpy reflects the excess entropy value at the transition point of the alloy. By adding alloy elements or adjusting the alloy composition based on of known alloys, beneficial alloy transition temperature values can be obtained. It can be observed from Figure 8 that as the Ti content increases, the mixing entropy of the alloy gradually increases, but the overall stability tends to be stable. The mixing enthalpy gradually decreases and tends to stabilize, as shown in Figure 9. observed from Figure 8 that as the Ti content increases, the mixing entropy of the alloy gradually increases, but the overall stability tends to be stable. The mixing enthalpy gradually decreases and tends to stabilize, as shown in Figure 9.   Gibbs free energy, also known as the Gibbs function, also known as free enthalpy, is an important thermodynamics parameter. It is often represented by G, and its definition is G = H − TS, Where H is enthalpy, T is temperature (absolute temperature, K), and S is entropy. Define ΔG = ΔH − TΔS (KJ/mol), G is called Gibbs free energy. Because H, T, and S are all state functions, G is the state function. The amount of Gibbs free energy change -ΔG = −(G2 − G1) ≥ W is not. It shows that the state function G is the ability of the system to do non-volume work under isothermal pressure and pressure. The reduction of G during the reaction-ΔG is the maximum non-volume work done by the system. This maximum is achieved in a reversible way. ΔG can judge the direction and mode of the reaction: −ΔG > W non-reaction proceeds spontaneously in an irreversible manner; −ΔG = W nonreaction proceeds in a reversible manner; −ΔG < W is not possible; if the reaction is carried out under isothermal pressure and pressure, no non-volume work is done, that is, W is not = 0 then ΔG < 0 proceed spontaneously; ΔG = 0 cannot proceed; ΔG > 0, the reverse reaction proceeds spontaneously. It can be seen that the direction in which the Gibbs free energy of the system decreases under isothermal pressure and pressure is the direction in which the chemical reaction that does not do non-volume work proceeds. The Gibbs free energy of any spontaneous process that does not do non-volume work under isothermal pressure and pressure will decrease. This paper uses JMatPro software to calculate the Gibbs free energy of Fe25Cr25Ni25TixAl(25−x) (x = 0; 6.25; 12.5; 18.75; 25) multi-principal element alloy, as shown in Figure 10. It can be seen from the figure that as the temperature rises high, the Gibbs free energy of the system gradually decreases, indicating that the reaction proceeds spontaneously. In the isothermal and isotactic state, as the Ti content increases, the Gibbs free energy decreases and the reaction proceeds spontaneously. Gibbs free energy, also known as the Gibbs function, also known as free enthalpy, is an important thermodynamics parameter. It is often represented by G, and its definition is G = H − TS, Where H is enthalpy, T is temperature (absolute temperature, K), and S is entropy. Define ∆G = ∆H − T∆S (KJ/mol), G is called Gibbs free energy. Because H, T, and S are all state functions, G is the state function. The amount of Gibbs free energy change −∆G = −(G2 − G1) ≥ W is not. It shows that the state function G is the ability of the system to do non-volume work under isothermal pressure and pressure. The reduction of G during the reaction-∆G is the maximum non-volume work done by the system. This maximum is achieved in a reversible way. ∆G can judge the direction and mode of the reaction: −∆G > W non-reaction proceeds spontaneously in an irreversible manner; −∆G = W nonreaction proceeds in a reversible manner; −∆G < W is not possible; if the reaction is carried out under isothermal pressure and pressure, no non-volume work is done, that is, W is not = 0 then ∆G < 0 proceed spontaneously; ∆G = 0 cannot proceed; ∆G > 0, the reverse reaction proceeds spontaneously. It can be seen that the direction in which the Gibbs free energy of the system decreases under isothermal pressure and pressure is the direction in which the chemical reaction that does not do non-volume work proceeds. The Gibbs free energy of any spontaneous process that does not do non-volume work under isothermal pressure and pressure will decrease. This paper uses JMatPro software to calculate the Gibbs free energy of Fe 25 Cr 25 Ni 25 Ti x Al (25−x) (x = 0; 6.25; 12.5; 18.75; 25) multi-principal element alloy, as shown in Figure 10. It can be seen from the figure that as the temperature rises high, the Gibbs free energy of the system gradually decreases, indicating that the reaction proceeds spontaneously. In the isothermal and isotactic state, as the Ti content increases, the Gibbs free energy decreases and the reaction proceeds spontaneously. Gibbs free energy, also known as the Gibbs function, also known as free enthalpy, is an important thermodynamics parameter. It is often represented by G, and its definition is G = H − TS, Where H is enthalpy, T is temperature (absolute temperature, K), and S is entropy. Define ΔG = ΔH − TΔS (KJ/mol), G is called Gibbs free energy. Because H, T, and S are all state functions, G is the state function. The amount of Gibbs free energy change -ΔG = −(G2 − G1) ≥ W is not. It shows that the state function G is the ability of the system to do non-volume work under isothermal pressure and pressure. The reduction of G during the reaction-ΔG is the maximum non-volume work done by the system. This maximum is achieved in a reversible way. ΔG can judge the direction and mode of the reaction: −ΔG > W non-reaction proceeds spontaneously in an irreversible manner; −ΔG = W nonreaction proceeds in a reversible manner; −ΔG < W is not possible; if the reaction is carried out under isothermal pressure and pressure, no non-volume work is done, that is, W is not = 0 then ΔG < 0 proceed spontaneously; ΔG = 0 cannot proceed; ΔG > 0, the reverse reaction proceeds spontaneously. It can be seen that the direction in which the Gibbs free energy of the system decreases under isothermal pressure and pressure is the direction in which the chemical reaction that does not do non-volume work proceeds. The Gibbs free energy of any spontaneous process that does not do non-volume work under isothermal pressure and pressure will decrease. This paper uses JMatPro software to calculate the Gibbs free energy of Fe25Cr25Ni25TixAl(25−x) (x = 0; 6.25; 12.5; 18.75; 25) multi-principal element alloy, as shown in Figure 10. It can be seen from the figure that as the temperature rises high, the Gibbs free energy of the system gradually decreases, indicating that the reaction proceeds spontaneously. In the isothermal and isotactic state, as the Ti content increases, the Gibbs free energy decreases and the reaction proceeds spontaneously.

First-Principles Prediction of Alloy Stability
To quantitatively evaluate the stability and elastic properties of the phase structure of Fe 25 Cr 25 Ni 25 TixAl (25−x) multi-principal element alloys, the first-principle used to calculate the formation heat, bonding energy, bulk modulus, shear modulus, and Youngs' modulus of the alloy. The modulus, Poisson's ratio, and the nature of the anisotropy predicted. According to the phase diagram analysis described above, the composition phase of the Fe 25 Cr 25 Ni 25 Ti x Al (25−x) multi-principal element alloys includes a disordered BCC1 phase, an ordered BCC2 phase, and a Laves phase. The crystal structure model of the bodycentered cubic structure and the close-packed hexagonal structure is established based on the Fe supercell and the Ti supercell. The Fe crystal space group of the body-centered cubic structure is IM-3M, the lattice constant is a = b = c = 2.866 Å, and the unit cell contains two atoms, and the positions are (0 0 0) and (0.5 0.5 0.5), respectively. The disordered BCC A2 phase and the ordered BCC B2 phase established on this basis contain 16 atoms. The Ti crystal space group of the close-packed hexagonal structure is P63/mmc, and the lattice constant is a = b = 2.9506 Å, c = 4.6788 Å, α = β = 90 • , γ = 120 • . The atomic positions in the unit cell are (0.333, 0.667, 0.25), (0.667, 0.333, 0.75), and the established Laves phase supercell contains 12 atoms. Multi-principal element alloys have higher entropy values and greater chaos, and the parent phase is disordered. Common modeling methods include Coherent Potential Approximation (CPA) [24] and Virtual Crystal Approximation (VCA) [25]. Methods such as Special Quasirandom Structure (SQS), the modeling methods in this paper are all modeled by SQS [26]. For the Fe 25 Cr 25 Ni 25 Ti 12.5 Al 12.5 alloy studied, a 2 × 2 × 2 supercell structure containing 16 atoms was used. In the modeling stage, an alloy of a certain composition was first established, and based on this, the atomic ratio adjusted to form an alloy of different compositions. For example, a Fe 25 Cr 25 Ni 25 Ti 12.5 Al 12.5 multi-principal element alloys was first established to establish a disordered BCC A2 phase, an ordered BCC B2 phase, and a Laves phase, respectively. Figure 11 shows the supercell structure of the disordered BCC A2 phase, the ordered BCC B2 phase, and the fully chemically disordered Laves phase in the Fe 25 Cr 25 Ni 25 Ti 12.5 Al 12.5 alloy. Table 2 shows the lattice parameters of the supercells calculated by the direct optimization method for the

First-Principles Prediction of Alloy Stability
To quantitatively evaluate the stability and elastic properties of the phase structure of Fe25Cr25Ni25TixAl(25−x) multi-principal element alloys, the first-principle used to calculate the formation heat, bonding energy, bulk modulus, shear modulus, and Youngs' modulus of the alloy. The modulus, Poisson's ratio, and the nature of the anisotropy predicted. According to the phase diagram analysis described above, the composition phase of the Fe25Cr25Ni25TixAl(25−x) multi-principal element alloys includes a disordered BCC1 phase, an ordered BCC2 phase, and a Laves phase. The crystal structure model of the body-centered cubic structure and the close-packed hexagonal structure is established based on the Fe supercell and the Ti supercell. The Fe crystal space group of the body-centered cubic structure is IM-3M, the lattice constant is a = b = c = 2.866 Å, and the unit cell contains two atoms, and the positions are (0 0 0) and (0.5 0.5 0.5), respectively. The disordered BCC A2 phase and the ordered BCC B2 phase established on this basis contain 16 atoms. The Ti crystal space group of the close-packed hexagonal structure is P63/mmc, and the lattice constant is a = b = 2.9506 Å, c = 4.6788 Å, α = β = 90°, γ = 120°. The atomic positions in the unit cell are (0.333, 0.667, 0.25), (0.667, 0.333, 0.75), and the established Laves phase supercell contains 12 atoms. Multi-principal element alloys have higher entropy values and greater chaos, and the parent phase is disordered. Common modeling methods include Coherent Potential Approximation (CPA) [24] and Virtual Crystal Approximation (VCA) [25]. Methods such as Special Quasirandom Structure (SQS), the modeling methods in this paper are all modeled by SQS [26]. For the Fe25Cr25Ni25Ti12.5Al12.5 alloy studied, a 2 × 2 × 2 supercell structure containing 16 atoms was used. In the modeling stage, an alloy of a certain composition was first established, and based on this, the atomic ratio adjusted to form an alloy of different compositions. For example, a Fe25Cr25Ni25Ti12.5Al12.5 multi-principal element alloys was first established to establish a disordered BCC A2 phase, an ordered BCC B2 phase, and a Laves phase, respectively. Figure 11 shows the supercell structure of the disordered BCC A2 phase, the ordered BCC B2 phase, and the fully chemically disordered Laves phase in the Fe25Cr25Ni25Ti12.5Al12.5 alloy. Table 2 shows the lattice parameters of the supercells calculated by the direct optimization method for the Fe25Cr25Ni25Ti12.5Al12.5 alloy. The theoretical calculations shown in the table are in good agreement with the experimental results. Figure 11. The supercells of the disordered BCC A2 structure (a), the ordered BCC B2 structure (b), and the Laves phase structure (c) employed in the first-principles calculations for Fe25Cr25Ni25Ti12.5Al12.5 multi-principle element alloys.   The strength and structural stability of crystals are closely related to their formation of heat and bonding energy. Heat is formed to determine which configuration is easier to form, and bonding energy is used to determine which configuration is more stable. The formation heat is defined as the difference between the energy of the compound and the corresponding composition. Usually, a negative value indicates that the substance can be stably present, and a positive value indicates that the substance is difficult to be formed. A smaller value indicates that the alloying ability is stronger. For Fe 25 Figure 12 is the calculated total ground state energy of Fe 25 Cr 25 Ni 25 Ti x Al (25-x) alloy atom assembly, forming heat and binding energy, and the result is negative, indicating that the alloy of the system is usually easy to form. From the energy point of view, the system of the alloy is stable. Further analysis, x = 0; 6.25; 12.5; 18.75; 25. With the increase of Ti content, the calculated total energy of bcc A2 phase and bcc B2 phase is getting smaller and smaller, indicating that the formation of the system is relatively stable. However, the calculated total energy for forming the Laves phase is getting larger and larger, indicating that the more difficult the Laves phase is to form, the system is not easy to stabilize. The formation heat of the bcc A2 phase gradually decreases and tends to be stable, and both are less than 0, indicating that it is easier to form. The binding energy gradually increases and tends to be stable and all are less than 0, indicating that the configuration with less Ti content is more stable. With the increase of Ti content, the formation heat and binding energy of Laves phase gradually decrease, indicating that Laves phase is easy to form in alloys with high Ti content and the configuration is more stable.

First-Principles Prediction of Alloy Elastic Properties
The elastic constant is the stiffness that characterizes the response of the crystal to the external strain ε. When a small strain is applied, the internal energy of the system has a quadratic linear relationship with the magnitude of the strain (Hooke's law), and the elastic constant is the coefficient of the quadratic term. The expression is: The strain tensor ε is defined as:

First-Principles Prediction of Alloy Elastic Properties
The elastic constant is the stiffness that characterizes the response of the crystal to the external strain ε. When a small strain is applied, the internal energy of the system has a quadratic linear relationship with the magnitude of the strain (Hooke's law), and the elastic constant is the coefficient of the quadratic term. The expression is: The strain tensor ε is defined as: where ∆E represents the energy difference before and after deformation of the unit cell, V 0 is the volume of the unit cell, C ij is the component of the elastic tensor, and e i and e j are the applied minor strains. For a stable structure, the Born mechanical stability standard for orthorhombic crystals at zero pressure is [27]: According to the calculated single-crystal elastic constant C ij , the polycrystalline elastic properties of the material such as the bulk modulus B, the shear modulus G, the elastic modulus E, and the Poisson can be obtained by the Voigt-Reuss-Hill (VRH) approximation method [28]. The ratios ν, Voigt, and Reuss approximate the maximum and minimum limits of polycrystalline elastic properties. The expressions are: B R = 1 (S 11 + S 22 + S 33 ) + 2(S 12 + S 13 + S 23 ) G R = 1 4(S 11 + S 22 + S 33 ) − 4(S 12 + S 13 + S 23 ) + 3(S 44 + S 55 + S 66 ) (10) where S ij is the elastic compliance of the material single crystal and S ij is the inverse matrix of C ij . The V and R in the subscript represent the calculation results of the Voigt and Reuss approximation models, respectively. For a hexagonal crystal, the standard for mechanical stability is [29]: From the calculated single-crystal elastic constants, the polycrystalline elastic property can be approximated by Voigt-Reuss-Hill (VRH), and the expression is [30]: [(C 11 + C 12 )C 33 − 2C 2 13 ]C 44 C 66 3B V C 44 C 66 + [(C 11 + C 12 )C 33 − 2C 2 13 ](C 44 + C 66 ) In the Fe 25 Cr 25 Ni 25 Ti x Al (25-x) system multi-principal element alloys, the single-crystal elastic constants C ij of the alloy calculated by the above formula are shown in Figure 13.
Through verification, the single-crystal elastic constants of the five alloys all meet the mechanical stability criterion, indicating that the five multi-principal element alloys structures established are mechanically stable. According to the calculated theoretical elastic constant of the single alloy crystal, the elastic properties of the alloy polycrystal obtained by the Voigt-Reuss-Hill approximation shown in Figure 12. The bulk modulus of a crystal is the resistance of a material to uniform compression of the outside under an elastic system. Shear modulus is the ratio of shear stress to strain, which characterizes the material's ability to resist shear strain. Young's modulus is one of the most essential and characteristic mechanical properties of elastic materials, reflecting the softness of materials to some extent hard. Poisson's ratio ν is commonly used to characterize the ability of a material to resist shear stress. The larger the Poisson's ratio, the better the shape of the material. In addition to Poisson's ratio, the ductile and brittleness of multi-principal element alloys can also be determined according to the Pugh empirical criteria [31], specifically defined as the ratio of bulk modulus to shear modulus (B/G). When B/G < 1.75, the material exhibits brittleness; when B/G > 1.75, the material exhibits toughness, and the larger the ratio, the stronger the toughness and the better the ductility of the material. Figure 14 shows the polycrystalline elastic properties of the five alloys. The values of shear modulus G and Young's modulus E in the alloy are found to be BCC B2 phase < BCC A2 phase < Laves phase; illustrated in the Fe 25 Cr 25 Ni 25 Ti x Al (25−x) system. Among the alloys, the Laves phase is the hardest and has the strongest resistance to shear deformation, followed by the BCC A2 phase, followed by the BCC B2 phase. When the value of bulk modulus B is x < 0.5, BCC B2 phase < BCC A2 phase < Laves phase; when x ≥ 0.5, BCC A2 phase < BCC B2 phase < Laves phase, indicating that Laves phase has the strongest volume deformation Resistance. By analyzing the Poisson's ratio and the value of B/G, it is found that the BCC B2 phase > BCC A2 phase > Laves phase indicates that the BCC B2 phase has excellent ductility, followed by the BCC A2 phase, and the Laves phase has the worst ductility. resistance of a material to uniform compression of the outside under an elastic system. Shear modulus is the ratio of shear stress to strain, which characterizes the material's ability to resist shear strain. Young's modulus is one of the most essential and characteristic mechanical properties of elastic materials, reflecting the softness of materials to some extent hard. Poisson's ratio ν is commonly used to characterize the ability of a material to resist shear stress. The larger the Poisson's ratio, the better the shape of the material. In addition to Poisson's ratio, the ductile and brittleness of multi-principal element alloys can also be determined according to the Pugh empirical criteria [31], specifically defined as the ratio of bulk modulus to shear modulus (B/G). When B/G < 1.75, the material exhibits brittleness; when B/G > 1.75, the material exhibits toughness, and the larger the ratio, the stronger the toughness and the better the ductility of the material. Figure 14 shows the polycrystalline elastic properties of the five alloys. The values of shear modulus G and Young's modulus E in the alloy are found to be BCC B2 phase < BCC A2 phase <Laves phase; illustrated in the Fe25Cr25Ni25TixAl(25−x) system. Among the alloys, the Laves phase is the hardest and has the strongest resistance to shear deformation, followed by the BCC A2 phase, followed by the BCC B2 phase. When the value of bulk modulus B is x < 0.5, BCC B2 phase < BCC A2 phase < Laves phase; when x ≥ 0.5, BCC A2 phase < BCC B2 phase < Laves phase, indicating that Laves phase has the strongest volume deformation Resistance. By analyzing the Poisson's ratio and the value of B/G, it is found that the BCC B2 phase > BCC A2 phase > Laves phase indicates that the BCC B2 phase has excellent ductility, followed by the BCC A2 phase, and the Laves phase has the worst ductility.   Figure 13. The calculated elastic constants (Cij) for single-phase Fe25Cr25Ni25TixAl(25−x). Figure 13. The calculated elastic constants (C ij ) for single-phase Fe 25

Conclusions
In this paper, a new type of Fe25Cr25Ni25TixAl(25−x) multi-principal element alloys is designed by using the research idea of material genetic engineering, which can improve the ductility of the alloy. It is realized by a combination of first-principles calculation and thermodynamic analysis. The main conclusions are as follows: (1) The phase composition, content and main constituent elements of Fe25Cr25Ni25TixAl(25−x) (x = 0; 6.25; 12.5; 18.75; 25) alloys were successfully predicted by thermodynamic phase diagram analysis. When x = 0, the alloy consists of a disordered BCC A2 phase and an ordered BCC B2 phase. When x = 6.25; 12.5; 18.75; 25, the phase consists of a disordered BCC A2 phase, an ordered BCC B2 phase, and a Laves phase. As the content of the Ti element continues to increase, the content of the disordered BCC A2 phase

Conclusions
In this paper, a new type of Fe 25 Cr 25 Ni 25 Ti x Al (25−x) multi-principal element alloys is designed by using the research idea of material genetic engineering, which can improve the ductility of the alloy. It is realized by a combination of first-principles calculation and thermodynamic analysis. The main conclusions are as follows: (1) The phase composition, content and main constituent elements of Fe 25 Cr 25 Ni 25 Ti x Al (25−x) (x = 0; 6.25; 12.5; 18.75; 25) alloys were successfully predicted by thermodynamic phase diagram analysis. When x = 0, the alloy consists of a disordered BCC A2 phase and an ordered BCC B2 phase. When x = 6.25; 12.5; 18.75; 25, the phase consists of a disordered BCC A2 phase, an ordered BCC B2 phase, and a Laves phase. As the content of the Ti element continues to increase, the content of the disordered BCC A2 phase and ordered BCC B2 phase in the alloy gradually decreases. The content of the Laves phase increases gradually.
(2) The structural stability and elastic properties of Fe 25 Cr 25 Ni 25 Ti x Al (25−x) (x = 0; 6.25; 12.5; 18.75; 25) alloys were predicted by first-principles calculation based on density functional theory. The calculation results of the formation heat and bonding energy of the alloy indicate that the alloy of the system is usually formed easily, and the alloy of the system has structural stability from the viewpoint of energy. By calculating the elastic properties of the alloy, such as bulk modulus B, shear modulus G, Young's modulus E, Poisson's ratio ν, and B/G, the results show that the disordered BCC A2 phase and the ordered BCC B2 phase are the ductile phases, the Laves phase is brittle. The ordered BCC B2 phase has the best ductility.
(3) We predicted the mixing entropy, mixing enthalpy and Gibbs free energy of Fe 25 Cr 25 Ni 25 Ti x Al (25−x) multi-principal element alloys. The results show that as the temperature increases, the mixing entropy and mixing enthalpy of the system are equally increase, indicating that the system tends to stabilize. With the increase of Ti content, the mixing entropy of the alloy gradually increases; the mixing enthalpy gradually decreases, but the overall stability tends to be stable. By calculating the Gibbs free energy of the system, the results show that as the temperature increases, the Gibbs free energy of the system gradually decreases, indicating that the reaction proceeds spontaneously. In the isothermal and isostatic state, as the Ti content increases, the Gibbs free energy decreases and the reaction proceeds spontaneously.