Effects of Different Contents of Each Component on the Structural Stability and Mechanical Properties of Co-Cr-Fe-Ni High-Entropy Alloys

: The structural stability of high-entropy alloys (HEAs) is closely related to their mechanical properties. The precise control of the component content is a key step toward understanding their structural stability and further determining their mechanical properties. In this study, ﬁrst-principle calculations were performed to investigate the effects of different contents of each component on the structural stability and mechanical properties of Co-Cr-Fe-Ni HEAs based on the supercell model. Co-Cr-Fe-Ni HEAs were constructed based on a single face-centered cubic (FCC) solid solution. Elemental components have a clear effect on their structure and performance; the Cr and Fe elements have an obvious effect on the structural stability and equilibrium lattice constant, respectively. The Ni elements have an obvious effect on stiffness. The Pugh ratios indicate that Cr and Ni addition may increase ductility, whereas Co and Fe addition may decrease it. With increasing Co and Fe contents or decreasing Cr and Ni contents, the structural stability and stiffness of Co-Cr-Fe-Ni HEAs are improved. The structural stability and mechanical properties may be related to the strength of the metallic bonding and covalent bonding inside Co-Cr-Fe-Ni HEAs, which, in turn, is determined by the change in element content. Our results provide the underlying insights needed to guide the optimization of Co-Cr-Fe-Ni HEAs with excellent mechanical properties. Author Contributions: Conceptualization, H.L., C.X., L.L., and C.Z.; Investigation, H.L., C.X., L.L., and C.Z.; Methodology, H.L., C.X., L.L., and C.Z.; Software, L.L.; Validation, H.L.; Visualization, H.L.; Writing—original draft, H.L.; Writing—review and editing, H.L., C.X., L.L., C.Z.


Introduction
High-entropy alloys (HEAs) have stimulated widespread interest due to their outstanding properties since they were independently proposed by Yeh and Cantor [1,2]. In general, HEAs are composed of four or more principal elements in equal or similar molar ratios, and the percentage range of each element is 5-35 at.%. HEAs have a high entropy of mixing. Lattice distortion [3], cocktail [4], and sluggish diffusion effects [5] in HEAs have also been reported. HEAs are a single-phase random solid solution. For example, CoCrFeNi HEAs [1] have a FCC structure. MoNbTaV refractory HEAs [6] have a body-centered cubic (BCC) structure. At present, some important works have shown that HEAs have many excellent mechanical properties, such as high strength, wear resistance, high hardness, and high-temperature oxidation resistance and corrosion resistance [6][7][8].

Computational Details
To study mechanical properties, it is necessary to evaluate structural stability. A solid-solution model was used to study the phase structure based on empirical parameters. The valence electron concentration (VEC) was used to determine the phase structure (BCC or FCC) [40,41]. For HEAs, the VEC is expressed as [40]: where c i represents the atomic percentage of the ith element, and (VEC) i is the VEC of the current element. For 3d transition elements, VEC includes all electrons together with d orbitals. The valence electrons of elements are Cr3d 5 4s 1 , Fe3d 6 4s 2 , Co3d 7 4s 2 , and Ni3d 8 4s 2 . The atomic size difference (δ) is an important parameter related to the solid-solution phase. A larger δ corresponds to severe lattice distortion and sluggish diffusion effects, resulting in a decrease in the stability and transformation rate of the solid-solution phase and precipitate formation in the HEAs [42]. The hardness changes monotonously with δ, which means that the lattice distortion hinders the dislocation motion of plastic deformation, resulting in solid-solution strengthening [40]. Accordingly, δ is defined as [42]: where r i represents the atomic radius of the ith element. We use the atomic radius measured by X-ray diffraction in the elemental state [43]. The enthalpy of mixing (∆H mix ) affects the solid-solution phase or compound formation in HEAs [44]. A negative ∆H mix means that different elements form intermetallic compounds, while positive a ∆H mix leads to phase separation in alloys. Only when the absolute value of ∆H mix is close to zero is the solid-solution phase stable in the solid phase. ∆H mix is defined as [45]: where ∆H mix AB is the enthalpy of mixing of the binary liquid alloy. The value of ∆H mix AB can be found from Ref. [46], which is calculated by the Miedema macroscopic model [47].
High-entropy alloys mean higher entropy of mixing (∆S mix ) and the value of ∆S mix is positive in the alloy. A higher ∆S mix corresponds to an increase in the degree of confusion in the multicomponent alloy system, and a decrease in the ordering and separation of different elements in the alloy, resulting in easier and more stable random solid-solution formation. ∆S mix is expressed as [42]: where R is the gas constant. HEAs of equal molar ratios have the largest entropy of mixing. According to thermodynamics, the energy T∆S mix increases with increasing temperature, and T∆S mix is related to the formation of the solid-solution phase. The phase formation of the alloy usually occurs near the melting point (T m ) if the solid-state phase transformation is not considered. In terms of the Gibbs free energy of multicomponent alloy systems, Yang and Zhang [42] defined a parameter Ω for predicting the solid-solution formation for HEAs as: c i (T m ) i is the melting point of HEAs following the mixing rules of multicomponent HEAs [42]. The T m of the atom can be found from Ref. [48]. Equations (1)- (5) or their combination can be used to evaluate the phase stability of HEAs, thereby laying a theoretical foundation for the study of mechanical properties.
The elastic constant calculation of Co-Cr-Fe-Ni HEAs can be implemented in the CASTEP package [49] by the stress-strain method based on Hooke's law. Hooke's law describes the linear relationship between stress component σ i and strain component ε j , and it is defined as [50]: where C ij represents the elastic constant and is described by a 6 × 6 tensor. The elastic constants in all directions can be calculated by the finite displacement method. According to the Voigt-Reuss-Hill approximation method [51], the polycrystal elastic moduli of the alloys was calculated from the elastic constants, including bulk modulus (B), shear modulus (G), Young's modulus (E), Pugh ratio (B/G), and Poisson's ratio (ν). In this work, a series of supercell solid-solution alloy models containing 20 and 32 atoms were established through the SQS method, and they were expanded to 1 × 1 × 5 and 2 × 2 × 2 on the basis of a single cell. The Co-Cr-Fe-Ni HEAs model and the percentage of each element are shown in Table 1 and Figure 1. The first-principles calculation uses the CASTEP package [49] based on plane-wave pseudopotential. This module is a quantum mechanics program developed by the University of Cambridge and has been widely used in the study of properties of complex multicomponent alloys [34,52]. We choose spin polarization and used formal spin as initial conditions because the ferromagnetic state is more stable than the paramagnetic state of room temperature at 0 K [29,53,54]. The Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) [55] was used to describe the exchange-correlation function. The ultrasoft pseudopotential was used to treat the electron-ion interaction. The energy cutoff for the plane-wave basis set was 400 eV. The self-consistent field tolerance is 1 × 10 −6 eV/atom, and the Brillouin zone was integrated using a Monkhorst-Pack grid. When x is 0.2 and 0.56, 1 and 1.57, and 2, the k-point set in the reciprocal space is 4 × 4 × 3, 4 × 3 × 4, and 1 × 7 × 7, respectively. The convergence tolerance is that the maximum force is less than 0.03 eV/Å, the energy is less than 1 × 10 −5 eV/atom, the maximum stress is less than 0.05 GPa, and the maximum displacement is less than 0.001 Å. Based on the above parameters and alloy models, the structure was first optimized, and then the elastic properties were calculated.

Structural Stability
To understand the structural properties of Co-Cr-Fe-Ni HEAs with different contents of each component, the VEC, δ, ∆ , ∆ , Ω, and tetragonal shear modulus (c′) were calculated and are shown in Figure 2. The VEC, δ, ∆ , ∆ and Ω for the equimolar ratio of the CoCrFeNi alloy are consistent with previous works [40,42]. According to previous results [40], HEAs prefer to form a single BCC phase when the VEC is 4.33-7.55, whereas the structure is FCC for VEC between 7.80 and 9.50. As shown in Figure 2a, all calculated VEC values are higher than 7.8, indicating that Co-Cr-Fe-Ni HEAs may maintain a single-phase FCC structure. This result is consistent with previous works [9][10][11]38]. The VEC of Cr and Fe is lower than the average valence electron concentration of the four elements, which causes the VEC of the alloys to decrease with increasing Cr and Fe contents. On the contrary, the VEC of the alloys increases with Co and Ni contents because the VEC of Co and Ni are higher. The δ is related to the atomic radius. Co has the largest atomic radius, and Fe has the smallest. This indicates that δ increases with Co and Fe contents, and decreases with Cr and Ni contents are shown in Figure 2b. HEAs form solid solutions with −15 ≤ ∆ ≤5 kJ mol −1 and 11 ≤ ∆ ≤ 19.5 J K −1 mol −1 [40]. As shown in Figure 2c,d, the values of ∆ and ∆ are in the above range, indicating that Co-Cr-Fe-Ni HEAs may be a solid-solution phase. This was confirmed by previous works [25,[36][37][38]. Only when the absolute value of ∆ is close to zero can the solid solution be stable in the solid phase. The results show that the solid-solution phase stability of the alloys decreases with increasing Cr and Ni contents and increases with increasing Co and Fe contents. This is consistent with previous work [38]. However, the entropy of mixing result shows that only the alloys with an equal molar ratio have the largest entropy of mixing, which is not affected by the component. For HEAs, Ω > 1 means that the contribution of ∆ will be higher than ∆ , corresponding to the formation of the solid-solution phase in the alloy [42]. As shown in Figure 2e, all alloys are solid solutions and become more stable with decreasing Cr and Ni content and increasing Co and Fe content. According to a previous report [42], a random solid solution is stable when δ ≤ 6.6% and Ω ≥ 1.1. Our results indicate that the Co-Cr-Fe-Ni HEAs studied in this work may be a random solid solution. This result is consistent with previous works [25,[36][37][38]. In addition, the parameter Ω

Structural Stability
To understand the structural properties of Co-Cr-Fe-Ni HEAs with different contents of each component, the VEC, δ, ∆H mix , ∆S mix , Ω, and tetragonal shear modulus (c ) were calculated and are shown in Figure 2. The VEC, δ, ∆H mix , ∆S mix and Ω for the equimolar ratio of the CoCrFeNi alloy are consistent with previous works [40,42]. According to previous results [40], HEAs prefer to form a single BCC phase when the VEC is 4.33-7.55, whereas the structure is FCC for VEC between 7.80 and 9.50. As shown in Figure 2a, all calculated VEC values are higher than 7.8, indicating that Co-Cr-Fe-Ni HEAs may maintain a single-phase FCC structure. This result is consistent with previous works [9][10][11]38]. The VEC of Cr and Fe is lower than the average valence electron concentration of the four elements, which causes the VEC of the alloys to decrease with increasing Cr and Fe contents. On the contrary, the VEC of the alloys increases with Co and Ni contents because the VEC of Co and Ni are higher. The δ is related to the atomic radius. Co has the largest atomic radius, and Fe has the smallest. This indicates that δ increases with Co and Fe contents, and decreases with Cr and Ni contents are shown in Figure 2b. HEAs form solid solutions with −15 ≤ ∆H mix ≤5 kJ mol −1 and 11 ≤ ∆S mix ≤ 19.5 J K −1 mol −1 [40]. As shown in Figure 2c,d, the values of ∆H mix and ∆S mix are in the above range, indicating that Co-Cr-Fe-Ni HEAs may be a solid-solution phase. This was confirmed by previous works [25,[36][37][38]. Only when the absolute value of ∆H mix is close to zero can the solid solution be stable in the solid phase. The results show that the solid-solution phase stability of the alloys decreases with increasing Cr and Ni contents and increases with increasing Co and Fe contents. This is consistent with previous work [38]. However, the entropy of mixing result shows that only the alloys with an equal molar ratio have the largest entropy of mixing, which is not affected by the component. For HEAs, Ω > 1 means that the contribution of T∆S mix will be higher than ∆H mix , corresponding to the formation of the solid-solution phase in the alloy [42]. As shown in Figure 2e, all alloys are solid solutions and become more stable with decreasing Cr and Ni content and increasing Co and Fe content. According to a previous report [42], a random solid solution is stable when δ ≤ 6.6% and Ω ≥ 1.1. Our results indicate that the Co-Cr-Fe-Ni HEAs studied in this work may be a random solid solution. This result is consistent with previous works [25,[36][37][38]. In addition, the parameter Ω shows a similar trend to that of the δ. The tetragonal shear modulus (c = (C 11 −C 12 )/2) is related to the energy difference of different crystal structures and also can reflect the structural stability [56]. As shown in Figure 2f, c decreases with increasing Cr and Ni contents and increases as the contents of Co and Fe increase. For HEAs composed of 3d transition elements, c is determined by the energy difference between the FCC and BCC structures of the multicomponent alloy system, which, in turn, is determined by the d orbital electrons [56]. Therefore, the energy difference of the alloys increases with Co and Fe contents and decreases with Cr and Ni contents. This result is consistent with the enthalpy of mixing, indicating that the structure is more stable as Co and Fe increase and Cr and Ni decrease. Overall, all parameters indicate that the Co-Cr-Fe-Ni HEAs studied in this work may be a single FCC solid-solution phase, and the structural stability decreases with increasing Cr and Ni contents and increases with increasing Co and Fe contents. It is worth noting that the range of VEC, ∆H mix , Ω, and c of the Cr element is the largest, which is consistent with previous work [38], indicating that Cr has a significant impact on structural stability. In addition, the Fe element has a significant effect on δ, and the influence on c is almost the same as that of Cr.
shows a similar trend to that of the δ. The tetragonal shear modulus (c′ = (C11−C12)/2) is related to the energy difference of different crystal structures and also can reflect the structural stability [56]. As shown in Figure 2f, c′ decreases with increasing Cr and Ni contents and increases as the contents of Co and Fe increase. For HEAs composed of 3d transition elements, c′ is determined by the energy difference between the FCC and BCC structures of the multicomponent alloy system, which, in turn, is determined by the d orbital electrons [56]. Therefore, the energy difference of the alloys increases with Co and Fe contents and decreases with Cr and Ni contents. This result is consistent with the enthalpy of mixing, indicating that the structure is more stable as Co and Fe increase and Cr and Ni decrease. Overall, all parameters indicate that the Co-Cr-Fe-Ni HEAs studied in this work may be a single FCC solid-solution phase, and the structural stability decreases with increasing Cr and Ni contents and increases with increasing Co and Fe contents. It is worth noting that the range of VEC, ∆ , Ω, and c′ of the Cr element is the largest, which is consistent with previous work [38], indicating that Cr has a significant impact on structural stability. In addition, the Fe element has a significant effect on δ, and the influence on c′ is almost the same as that of Cr. The equilibrium lattice constants and mass density of the Co-Cr-Fe-Ni HEAs based on the FCC supercell structure are presented in Figure 3. The arithmetic mean cell length of the supercell is used as the lattice constants of the alloy for comparison with experimental data.
is calculated as follows [57]: The equilibrium lattice constants and mass density of the Co-Cr-Fe-Ni HEAs based on the FCC supercell structure are presented in Figure 3. The arithmetic mean cell length of the supercell a avg is used as the lattice constants of the alloy for comparison with experimental data. a avg is calculated as follows [57]: where n represents the number of unit cells in the supercell, and a i , b i , and c i are lattice constants. The lattice constant for the equimolar ratio is 0.35 nm, consistent with previous works, which found it as 0.3498-0.3587 nm [27][28][29]32,33,36,58]. In Figure 3a, the lattice constant decreases with increasing Co and Fe content and increases with increasing Cr and Ni content. The lattice constant of the Fe element has the largest variation range, which is consistent with that of δ, indicating that Fe has the most significant influence on the equilibrium lattice constant. The lattice constants and atomic size difference of Co-Cr-Fe-Ni HEAs vary with the content of each element in the opposite direction. Fe has the smallest atomic radius, which can explain the change in lattice constants and atomic size difference with Fe content. However, the lattice constants decrease with Co contents even though Co has the largest atomic radius. This shows that the atomic radius does not determine the lattice constant of the alloy. This result is consistent with previous work [54], which demonstrated that lattice distortion cannot be accurately described by empirical parameter δ based on the atomic radius, especially for 3d FCC alloys. The variation in VEC can cause charge transfer resulting in distortions even for atoms of FCC Co-Ce-Fe-Ni HEAs with similar atomic sizes [59,60]. Chemical interactions between the elements, electronic spin-flip, and magnetovolume effects in HEAs can maximize the effect of lattice distortion caused by Cr [33,34]. For Co and Fe atoms, the (Fe ↔ Cr) and (Co ↔ Cr) bonds are stronger than (Fe ↔ Ni) and (Co ↔ Ni) bonds (A ↔ B means the atomic bond of A and B elements), hence (Fe ↔ Cr) and (Co ↔ Cr) bonds will be preferred over (Fe ↔ Ni) and (Co ↔ Ni) bonds [33]. For Co-Cr-Fe-Ni HEAs, the increase in Cr-related bonds will enhance the lattice distortion effect caused by Cr. This may be the reason why the influence of the Cr element on the equilibrium lattice constant is almost the same as that of Ni. Therefore, the lattice distortion caused by the increase in Cr and Ni content or decrease in Co and Fe content may be the cause of the increase in the lattice constant.
Fe-Ni HEAs vary with the content of each element in the opposite direction. Fe has the smallest atomic radius, which can explain the change in lattice constants and atomic size difference with Fe content. However, the lattice constants decrease with Co contents even though Co has the largest atomic radius. This shows that the atomic radius does not determine the lattice constant of the alloy. This result is consistent with previous work [54], which demonstrated that lattice distortion cannot be accurately described by empirical parameter δ based on the atomic radius, especially for 3d FCC alloys. The variation in VEC can cause charge transfer resulting in distortions even for atoms of FCC Co-Ce-Fe-Ni HEAs with similar atomic sizes [59,60]. Chemical interactions between the elements, electronic spin-flip, and magnetovolume effects in HEAs can maximize the effect of lattice distortion caused by Cr [33,34]. For Co and Fe atoms, the (Fe ↔ Cr) and (Co ↔ Cr) bonds are stronger than (Fe ↔ Ni) and (Co ↔ Ni) bonds (A ↔ B means the atomic bond of A and B elements), hence (Fe ↔ Cr) and (Co ↔ Cr) bonds will be preferred over (Fe ↔ Ni) and (Co ↔ Ni) bonds [33]. For Co-Cr-Fe-Ni HEAs, the increase in Cr-related bonds will enhance the lattice distortion effect caused by Cr. This may be the reason why the influence of the Cr element on the equilibrium lattice constant is almost the same as that of Ni. Therefore, the lattice distortion caused by the increase in Cr and Ni content or decrease in Co and Fe content may be the cause of the increase in the lattice constant. As shown in Figure 3b, the mass density of Co-Cr-Fe-Ni HEAs increases with Co, Fe, and Ni contents and decreases with Cr contents. The Cr element has the most significant influence on the density. The decrease in mass density with Cr content is due to the smallest density of Cr and the increase in lattice constant. The densities of Co and Fe are larger and the lattice constants decrease, causing the mass density to increase with the content of Co and Fe. Ni has the highest density, but the lattice constants increase with Ni content, and therefore the density of the alloy increases first (x ≤ 1.57) and then decreases (x > 1.57). That is, the relative effect of Ni content on the mass and volume of the alloy is different with different content. When x > 1.57 of Ni, the mass increase of the alloy is less than the volume increase, and the density appears to decrease.  As shown in Figure 3b, the mass density of Co-Cr-Fe-Ni HEAs increases with Co, Fe, and Ni contents and decreases with Cr contents. The Cr element has the most significant influence on the density. The decrease in mass density with Cr content is due to the smallest density of Cr and the increase in lattice constant. The densities of Co and Fe are larger and the lattice constants decrease, causing the mass density to increase with the content of Co and Fe. Ni has the highest density, but the lattice constants increase with Ni content, and therefore the density of the alloy increases first (x ≤ 1.57) and then decreases (x > 1.57). That is, the relative effect of Ni content on the mass and volume of the alloy is different with different content. When x > 1.57 of Ni, the mass increase of the alloy is less than the volume increase, and the density appears to decrease.

Mechanical Properties
The mechanical stability of crystal systems is the basis for studying HEAs' hardness and stiffness. FCC crystals usually have three independent elastic constants, but in this work, there are nine independent elastic constants belonging to the orthogonal structure Appl. Sci. 2021, 11, 2832 8 of 13 due to the use of the supercell model. The mechanical stability of the orthogonal structures can be determined from the following equations [52]: The calculated elastic constants of Co-Cr-Fe-Ni HEAs are listed in Table 2, lacking experimental data comparison. Our results show that the elastic constants of Co-Cr-Fe-Ni HEAs satisfy the stability criteria, indicating that the alloys are mechanically stable in the ground state. Table 2. Elastic constant C ij of Co-Cr-Fe-Ni HEAs. The unit is GPa.
The elastic constants and polycrystal elastic moduli of CoCrFeNi alloy are shown in Table 3. For the other alloys, only limited data are available. The values for the CoCrFeNi alloy are in agreement with the calculations of Ref. [58]. The error of the B may be affected by the alloy model. B/G and ν of CoCrFeNi alloy also agree well with the values of Refs. [37,50]. However, the calculated C ij , B, G, and E are much higher than the previous results [37,50]. This may be caused by the alloy model, calculation package, and magnetism, which have a great influence on the calculation results [36,48,51]. The comparison of Refs. [50,58] shows that the alloy model may be the main influence factor of the elasticity results. Table 3. Elastic constants (C ij ) and polycrystal elastic moduli (B, G, E, ν, and B/G) for CoCrFeNi alloy. The unit is GPa, except for the dimensionless ν and B/G. The brittle-ductile behavior is an important mechanical performance indicator of the energy absorption and plastic deformation ability of materials under external forces. The Pugh ratio is usually used to describe the brittle-ductile behavior of HEAs [61]. Ductile materials typically show B/G > 1.75, while brittle materials show B/G < 1.75. Poisson's ratio can also reflect the brittle-ductile behavior of HEAs. For isotropic materials, B/G > 1.75 implies ν > 0.26 [37], which has been confirmed by Gu et al. [62]. As shown in Figure 4a,b, according to B/G and ν, we find that all Co-Cr-Fe-Ni HEAs are predicted to be ductile, except for Fe-1.57. For CoCrFeNi alloy, this is in line with experiments [15,18]. Unfortunately, no available experimental information on the brittle-ductile behavior was reported for the other HEAs that can be compared with our results. ν has the same changing trend as B/G. Alloys have high values of ν and B/G as the contents of Cr and Ni increase and the contents of Co and Fe decrease, indicating that the ductility of the alloys enhances. Fe element has the largest variation range of B/G and ν, indicating that Fe may have the most significant effect on ductility. However, the effect of the Cr element on ductility is almost the same as Fe. The Cauchy pressure (C 12 -C 44 ) reflects the bonding characteristics of the alloy. A more positive C 12 -C 44 means more metallic bonding, and otherwise more covalent bonding [63]. As shown in Figure 4c, alloys have lower negative C 12 -C 44 as the content of Cr and Ni decreases and the content of Co and Fe increases, which indicates that the covalent nature of the metallic bond is enhanced. Covalent bonds are more stable than metal bonds, which provides an explanation for the structural stability of the Co-Cr-Fe-Ni HEAs. The stronger the alloy's metallicity, the better the ductility of the materials. The trend of B/G and ν is consistent with that of C 12 -C 44 . In addition, the severe lattice distortion corresponding to the higher δ hinders the dislocation slip of HEAs in plastic deformation and increases yield strength. Therefore, our results indicate that the structural stability of the Co-Cr-Fe-Ni HEAs decreases but the ductility enhances with increasing Cr and Ni content and decreasing Co and Fe content.
The bulk modulus (B) represents the ability of the alloys to resist volume deformation under external pressure, the shear modulus (G) represents the ability to resist reversible deformation under shear stress, and Young's modulus (E) represents the ability to resist compression or tension within the elastic limit and can describe the stiffness of the materials. As shown in Figure 4d, the bulk modulus increases with increasing Fe content and decreases with increasing Ni content. Ni element has the most significant effect on the bulk modulus. The variation of Co and Cr has a slight effect on the bulk modulus, and the calculated bulk modulus is 255-262 GPa, indicating that the incompressibility of all Co and Cr series alloys is similar. This is consistent with previous work [38]. The variation of the shear modulus with composition is consistent with the variation of Young's modulus, as shown in Figure 4e,f. They increase with increasing Co and Fe contents and decrease with increasing Cr and Ni contents. The G and E variation ranges of Fe and Ni elements are almost the same. With regard to Cr, the shear modulus and Young's modulus increase firstly and then decrease. This change is consistent with the change in the electronic charge gain of the atoms in FCC-Cr material [34]. Among the four elements, only the 4s orbital unsaturation of Cr may change the chemical interactions with other elements [33,34], as well as the unique magnetic characteristics of Cr [53], which affects the structure and elastic properties. Stiffer materials generally have a larger Young's modulus. Hence, increasing Co and Fe contents or decreasing Cr and Ni contents would improve the stiffness of Co-Cr-Fe-Ni HEAs, and the Ni element has the most significant effect on the stiffness. Combining structural stability and brittle-ductile behavior, we find that the variation of stiffness with composition is consistent with that of structural stability but opposite to the brittle-ductile behavior. The relationship between stiffness and ductility is consistent with previous work on Al x FeCuCrNi HEAs [64]. This may be because the covalent nature of the metal bond improves the structural stability, and stiffness, but reduces the ductility. For FCC Co-Cr-Fe-Ni HEAs, the atomic-level pressure variation originating from the charge transfer between neighboring elements in the alloy will also improve the stiffness of the system [59]. Overall, the incompressibility, resistance to shear deformation, and stiffness of Co-Cr-Fe-Ni HEAs improve with the increasing Co and Fe contents or decreasing Cr and Ni contents, corresponding to a more stable solid-solution structure and weakened ductility. stiffness with composition is consistent with that of structural stability but opposite to the brittle-ductile behavior. The relationship between stiffness and ductility is consistent with previous work on AlxFeCuCrNi HEAs [64]. This may be because the covalent nature of the metal bond improves the structural stability, and stiffness, but reduces the ductility. For FCC Co-Cr-Fe-Ni HEAs, the atomic-level pressure variation originating from the charge transfer between neighboring elements in the alloy will also improve the stiffness of the system [59]. Overall, the incompressibility, resistance to shear deformation, and stiffness of Co-Cr-Fe-Ni HEAs improve with the increasing Co and Fe contents or decreasing Cr and Ni contents, corresponding to a more stable solid-solution structure and weakened ductility.

Conclusions
In summary, we used a supercell model based on empirical parameters and firstprinciple calculations to study the effects of different contents of each component on the structural stability and mechanical properties of Co-Cr-Fe-Ni HEAs. The valence electron

Conclusions
In summary, we used a supercell model based on empirical parameters and firstprinciple calculations to study the effects of different contents of each component on the structural stability and mechanical properties of Co-Cr-Fe-Ni HEAs. The valence electron concentration, enthalpy of mixing, entropy of mixing, atomic size difference, and Ω indicate that the Co-Cr-Fe-Ni HEAs studied in this work may be a single FCC solid-solution phase, and all alloys are mechanically stable at the ground state. The structural stability increases with increasing Co/Fe content and with decreasing Cr/Ni content. The equilibrium lattice constants increase with increasing Cr/Ni content and with decreasing Co/Fe content. The mass density increases with Co, Fe, and Ni content and decreases with Cr content. The ductility may enhance with increasing Cr/Ni content and decreasing Co/Fe content. The stiffness improves with increasing Co/Fe contents or decreasing Cr/Ni contents. Certain components of the alloy significantly affect its performance; specifically, Cr affects the structural stability, Fe impacts the equilibrium lattice constant, Cr affects the mass density, Fe affects the ductility, and Ni impacts the stiffness. The results show that the stiffness and structural stability of Co-Cr-Fe-Ni HEAs are improved, while the ductility may be weakened with increasing Co/Fe content or decreasing Cr/Ni content. The structural stability and mechanical properties may be related to the strength of the metallic bonding and covalent bonding inside the Co-Cr-Fe-Ni HEAs, which, in turn, is determined by the change in element content.