Theoretical Analysis of Stacking Fault Energy, Elastic Properties, Electronic Properties, and Work Function of MnxCoCrFeNi High-Entropy Alloy

The effects of different Mn concentrations on the generalized stacking fault energies (GSFE) and elastic properties of MnxCoCrFeNi high-entropy alloys (HEAs) have been studied via first-principles, which are based on density functional theory. The relationship of different Mn concentrations with the chemical bond and surface activity of MnxCoCrFeNi HEAs are discussed from the perspectives of electronic structure and work function. The results show that the plastic deformation of MnxCoCrFeNi HEAs can be controlled via dislocation-mediated slip. But with the increase in Mn concentration, mechanical micro twinning can still be formed. The deformation resistance, shear resistance, and stiffness of MnxCoCrFeNi HEAs increase with the enhancement of Mn content. Accordingly, in the case of increased Mn concentration, the weakening of atomic bonds in MnxCoCrFeNi HEAs leads to the increase in alloy instability, which improves the possibility of dislocation.


Introduction
High-entropy alloys (HEAs), also known as multicomponent alloys, are alloys consisting of five or more equal or approximately equal amounts of metals.The concentration of each element in the alloy is between 5 and 35% [1,2].It has a high entropy effect, extreme lattice distortion effect, cocktail effect, and sluggish diffusion effect [3][4][5][6].In recent years, more and more people have paid attention to high-entropy alloys because of their unique compositions and excellent properties [7][8][9][10].In particular, the high-entropy alloy Fe 20 Cr 20 Mn 20 Ni 20 Co 20 , with equal proportions, also known as the Cantor alloy, was first introduced in 2004 [11].Therefore, this kind of CrMnFeCoNi alloy has long been a focus and hotspot of research.
So far, many researchers have extensively studied CoCrFeNi-based high-entropy alloys.Li et al. [12] designed a kind of metastable high-entropy dual-phase alloy, Fe 80-x Mn x Co 10 Cr 10 .They concluded that the single face-centered cubic (FCC) phase structure could be satisfied when the Mn content is 40 and 45% (Fe 35 Mn 45 Co 10 Cr 10 and Fe 40 Mn 40 Co 10 Cr 10 , respectively).The Mn content plays an important role in phase composition, regulating phase stability and improving phase transition mechanism.Sun et al. [13] compared the lattice stability of Al x CrMnFeCoNi and Al x CrFeCoNi high-entropy alloys.It was found that Mn decreases the stable field of FCC phase and widens the width of the two-phase region.Moreover, via the exact muffin-tin orbitals (EMTO) method, Zhang et al. [14] calculated the elastic properties of body-centered cubic (BCC) and face-centered cubic (FCC) Al x CrMnFeCoNi (0 ≤ x ≤ 5) HEAs, demonstrating that there is a complex dependence between the elastic parameters and the composition, and the elastic anisotropy of both phases is extremely high.Shi Y.Z.et al. [15] studied the homogenization effect of 1250 • C heat treatment on Al x CoCrFeNi HEAs.After heat treatment, the homogenization effect of microstructure leads to the decrease in work function and the improvement of corrosion resistance.Zhang et al. [16] prepared CoCrFeNi-Nb x (x = 0, 1, 3, 5, 7, 9 wt%) high-entropy alloys through high-energy ball milling and discharge plasma sintering.The effects of niobium on the microstructure and properties of cobalt-nickel alloys were studied systematically.Nb atoms cause lattice distortion of the alloy, and the microstructure of CoCrFeNi HEAs changes from a single-phase structure of FCC to a bi-phase structure of FCC and Laves, which increases the tensile strength, yield strength, and hardness of the HEA.
Moreover, Kivy M.B. et al. [17] investigated the effects of Cu, Mn, Al, Ti, and Mo on the generalized stacking fault energies (GSFE), Rice-criterion ductilities, and twin ability of CoCrFeNi-based HEAs with FCC structure via density functional theory.The results presented that the addition of Ti and Mo increases the tendency of dislocation slip and deformation twinning.Furthermore, the addition of Mn, Cu, or Al with high content promoted dislocation slip and martensitic phase transition, while a low amount of Al led to dislocation slip.Achmad et al. [18] have calculated pure cobalt and Co-9 at.%X solid-solution alloys (X = Cr, W, Mo, Ni, Mn, Al, Fe) using first-principles density-functional theory.They found that the alloying of Cr, W, and Mo increases the interlayer distance distortion of the stacking fault plane and reduces the GSFE value of pure Co due to the increase in charge accumulation.Also, the stacking fault energy (SFE) values of several typical FCC high-entropy alloys (HEAs) were measured by Liu et al. [19] using experimental means.The experimental results indicated that the lower SFE is helpful for the formation of deformation twins under the loading condition, and the smaller the thickness is, the better the mechanical properties are at low temperature.
It is obvious from the above research results that the Mn element plays an important role in the CoCrFeNi-based high-entropy alloy.Although a great deal of research has been conducted on the effects of alloying elements, especially Mn, on high-entropy alloys, there is still a lack of theoretical research on the GSFE and work function of CoCrFeNi-based high-entropy alloys with different Mn.In general, the GSFE has been widely studied in characterization of the mechanical properties and brittle-to-ductile transition of alloys.The GSFE is considered to be a measure of energy penalty between two adjacent planes during shear deformation along a specific slip direction on a given slip plane.It represents the nature of the slip and involves stable and unstable stacking and twin fault energy.In particular, the intrinsic stacking fault energy (ISF) of the alloy can be calculated via transmission electron microscopy (TEM) [20,21] and X-ray powder diffraction (XRD) [22].However, the unstable stacking fault energy (USF) cannot be measured by experiment, and the first principle calculation is an effective method by which to measure the USFE [23].On the other hand, virtual crystal approximation (VCA) is considered a method by which to study the properties of solid solutions.The method uses "virtual" atoms inserted between atomic behaviors in the original compound to study crystals in primitive periodicity.This method makes the calculation simpler and the cost lower.Of course, previous work [24,25] has proved that VCA has good accuracy in some HEAs.
Therefore, based on the first-principle calculation, the theoretical models of Mn x CoCrFeNi HEAs are established via the VCA method in the present work.In order to study the effect of different Mn concentrations on slip and twin in Mn x CoCrFeNi HEAs, 13 close-packed (111) atomic layers are used to calculate GSFE, as well as USF, ISF, unstable twin-fault energy (UTF), and two-layers twinning stacking fault energy (TSF) parameters.More importantly, the elastic properties and electron density difference were systematically studied with the aim of further uncovering the nature of the strength and ductility of Mn x CoCrFeNi HEAs.Finally, the work functions of (111), (110), and (100) planes in Mn x CoCrFeNi HEAs are estimated, which reveals the analysis results in more details.Therefore, theoretical analysis was conducted on HEAs containing different proportions of Mn elements in this study.The influence of Mn element on the GSFE and other properties of CoCrFeNi-based HEAs can be obtained, which can provide useful guidance for the design of high-performance CoCrFeNi-based HEAs.

First-Principles Calculations
All the first-principles calculations are performed using the Cambridge sequential total energy package (CASTEP) (MS 7.0) [26,27] based on density functional theory (DFT) [28,29].The exchange correlation functional uses the Pardew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) [30,31].Clearly, the premise for simulation is geometry optimization (minimum energy of atoms at different volumes).Therefore, according to the termination of structural relaxation, the convergence parameters are selected.The convergence parameters were set as follows: total energy tolerance is 10 −5 eV/atom; force tolerance is 0.03 eV/Å; maximum stress is 0.05 GPa; and maximum displacement is 0.001 Å.After the convergence test, the plane wave energy cutoff is 600 eV and the Monkhorst-Pack scheme [32,33] k-points set is 10 × 10 × 10 in the Brillouin zone.According to the selected setting, the error reaches less than 1% through the convergence of the total energy of the calculated model.

VCA Models
The construction of VCA unit cells is achieved by replacing real alloy atoms in the cell structure with "virtual" atoms, which are obtained from the weighted average technique of different alloy elements.The atomic percentages of each constituent atom within the cell of the Mn x CoCrFeNi high-entropy alloy are presented in Table 1.Accurately determining the lattice constant of the unit cell is an essential prerequisite for predicting the properties of compounds in theoretical calculations.In this paper, the comparison between the optimized lattice constant and the lattice constant obtained from other calculation methods and experimental results is detailed in Table 2.In Table 2, SQS and CPA represent the estimation of lattice constants of high-entropy alloys using a special quasi random structure (SQS) and coherent-potential approximation (CPA), respectively.As shown in Table 2, the calculated lattice constant is consistent with other calculation methods and experimental results.

Stacking Fault Model and GSFE
The structure diagram of the FCC structural unit along [111], [110], and [112] directions is shown in Figure 1a.Based on the fact that dislocations in the FCC structure mainly occur on the (111) surface of dense packing, a supercell model with 13 closed-stacking (111) is constructed, in which the stacking sequence from bottom to top is "ABCABCABCABCA".In Figure 1b, the GSFE value was obtained vis shear slip along [112] in the plane (111).In the first process, the upper 7-13 layers of atoms are sheared along the [112] direction.With the step size of 0.1b p , the unstable stacking fault (USF) is generated at the shear displacement of 0.0 to 0.5 b p , in which the Burgers vector b p = a 0 / √ 6.When the shear displacement increases to 1.0 b p , the structure becomes the ISF of ABCABCBCABCAB sequence.The second process is to move the upper 8-13 layers of atoms at the same distance along the [112] direction.When the shear displacement value increases to 1.0 b p , the structural sequence becomes ABCABCBABCABC, which is the TSF.The GSFE is calculated as follows [18]: where E u is the total energy of the supercell after shear displacement, E 0 is the total energy of the supercell without defect, and A is the area of the fracture plane.To avoid periodic atomic interactions, a 15 Å vacuum layer is added to the structure.After adequate convergence testing for the GSFE calculation, and considering the calculation efficiency, the plane wave energy cutoff and the k-points are set as 400 eV and 10 × 10 × 1, respectively.

Stacking Fault Model and GSFE
The structure diagram of the FCC structural unit along [111], [1 10], and [112 ] tions is shown in Figure 1a.Based on the fact that dislocations in the FCC structure m occur on the (111) surface of dense packing, a supercell model with 13 closed-sta (111) is constructed, in which the stacking sequence from bottom to t "ABCABCABCABCA".In Figure 1b, the GSFE value was obtained vis shear slip [112 ] in the plane (111).In the first process, the upper 7-13 layers of atoms are sh along the [112 ] direction.With the step size of 0.1bp, the unstable stacking fault (U generated at the shear displacement of 0.0 to 0.5 bp, in which the Burgers vector  0 √6 ⁄ .When the shear displacement increases to 1.0 bp, the structure becomes the ABCABCBCABCAB sequence.The second process is to move the upper 8-13 lay atoms at the same distance along the [112 ] direction.When the shear displacement increases to 1.0 bp, the structural sequence becomes ABCABCBABCABC, which is th The GSFE is calculated as follows [18]: where Eu is the total energy of the supercell after shear displacement, E0 is the total e of the supercell without defect, and A is the area of the fracture plane.To avoid pe atomic interactions, a 15 Å vacuum layer is added to the structure.After adequate co gence testing for the GSFE calculation, and considering the calculation efficienc plane wave energy cutoff and the k-points are set as 400 eV and 10 × 10 × 1, respecti

Surface Energies
In order to better understand the surface properties of Mn x CoCrFeNi HEAs, the surface energies (γ s ) of ( 100), (110), and (111) planes were calculated.Figure 2 illustrates the (100), (110), and (111) surface models of Mn x CoCrFeNi HEAs, which contain 10 atomic layers.A vacuum layer with a thickness of 5 Å is set on the bottom and top of the surface Materials 2024, 17, 4378 5 of 16 model to eliminate the influence of the interaction between periodic structures, as shown in Figure 2a,c,e.With the aim of avoiding the interaction between atoms in the slab, an 8 Å vacuum layer is set between the atomic layers of the slab model (Figure 2b,d,f).The surface energy can be calculated by the following formula [25,37]: here, E slab is the total energy of the surface model.E Bulk is the total energy of the perfect supercell, and A is the surface area.

Surface Energies
In order to better understand the surface properties of MnxCoCrFeNi HEAs, th face energies (γs) of ( 100), (110), and (111) planes were calculated.Figure 2 illustrate (100), (110), and (111) surface models of MnxCoCrFeNi HEAs, which contain 10 a layers.A vacuum layer with a thickness of 5 Å is set on the bottom and top of the su model to eliminate the influence of the interaction between periodic structures, as sh in Figure 2a,c,e.With the aim of avoiding the interaction between atoms in the slab Å vacuum layer is set between the atomic layers of the slab model (Figure 2b,d,f) surface energy can be calculated by the following formula [25,37]: here,  is the total energy of the surface model.
is the total energy of the fect supercell, and A is the surface area.

Figure 3 and Table 1 exhibit the GSFE curves and related values of MnxCoCr
HEAs with different Mn concentrations, which were calculated based on first-princ According to the Figure 3, with the increase in shear displacement along the directi [112 ], the structure exhibits an energy barrier at the first highest energy point on the curve, which is called the unstable stacking fault energy (γusf).This energy barrier is sidered as the minimum energy or critical stress required for local dislocation nucle

Figure 3 and Table 1 exhibit the GSFE curves and related values of Mn x CoCrFeNi
HEAs with different Mn concentrations, which were calculated based on first-principles.According to the Figure 3, with the increase in shear displacement along the direction of [112], the structure exhibits an energy barrier at the first highest energy point on the GSFE curve, which is called the unstable stacking fault energy (γ usf ).This energy barrier is considered as the minimum energy or critical stress required for local dislocation nucleation.As the structure is further sheared, some dislocations begin to spread, and a stacking fault defect is created.The first minimum energy point on the GSFE curve is formed, and this energy is called intrinsic or stable stacking fault energy (γ isf ).The second maximum energy point on the GSFE curve is called the unstable twin fault energy (γ utf ).This energy is considered to be the minimum energy barrier to produce an extrinsic or twinning stacking fault.The second minimum energy point is defined as the extrinsic stacking fault energy or two-layer twinning stacking fault energy (2γ tsf ).The GSFE calculation based on the firstprinciple separates the formation energy from any other related effects on total energy, so the calculated value may be more accurate than the experimental measurement value [18].
Van Swygenhoven et al. [38] found that only using SFE parameters as a criterion for judging the twinning and partial dislocation deformation mechanism of materials is not sufficient.curves of MnxCoCrFeNi HEAs gradually decreased.The USF energy decreases 808.84 mJ/m 2 to 276.17 mJ/m 2 with increasing Mn content from 0.0 to 1.0.The re demonstrate that the alloy with higher Mn content easily forms dislocation due to its nuclear dislocations resistance.The two-layer twinning SFE of HEAs was lower than of the stable SFE.And the decrease degree of the two-layer twinning SFE was greater that of stable SFE.When the Mn content was 0, the two-layer twinning SFE of HEAs ( = 574.71mJ/m 2 ) was lower than the stable SFE (γisf = 578.93mJ/m 2 ).When the Mn con was 1, the value of the two-layer twinning SFE is 60.36 mJ/m 2 , and the value of the s SFE is 110.04 mJ/m 2 .The ratio of ISF energy to USF energy (γisf/γusf) is often used to indicate the tend of complete dislocation dissociation into partial dislocation.The lower the ratio greater the tendency of full dislocation dissociation.When the leading partial disloca nucleate by overcoming the energy barrier γusf, the tailing part needs to exceed this en As can be seen from the Figure 3, with the increase in Mn concentration, the GSFE curves of Mn x CoCrFeNi HEAs gradually decreased.The USF energy decreases from 808.84 mJ/m 2 to 276.17 mJ/m 2 with increasing Mn content from 0.0 to 1.0.The results demonstrate that the alloy with higher Mn content easily forms dislocation due to its low nuclear dislocations resistance.The two-layer twinning SFE of HEAs was lower than that of the stable SFE.And the decrease degree of the two-layer twinning SFE was greater than that of stable SFE.When the Mn content was 0, the two-layer twinning SFE of HEAs (2γ tsf = 574.71mJ/m 2 ) was lower than the stable SFE (γ isf = 578.93mJ/m 2 ).When the Mn content was 1, the value of the two-layer twinning SFE is 60.36 mJ/m 2 , and the value of the stable SFE is 110.04 mJ/m 2 .
The ratio of ISF energy to USF energy (γ isf /γ usf ) is often used to indicate the tendency of complete dislocation dissociation into partial dislocation.The lower the ratio, the greater the tendency of full dislocation dissociation.When the leading partial dislocations nucleate by overcoming the energy barrier γ usf , the tailing part needs to exceed this energy for nucleation.The critical stress for the nucleation of the trailing part is a function of (γ usf − γ isf ).The increased value of γ usf − γ isf is beneficial to the partial dislocation, and the stacking fault more likely to form.The difference between γ usf and γ isf in Mn x CoCrFeNi HEAs was large, implying that the deformation mechanism was only extended partial dislocation.Moreover, with the increase in Mn content, the difference value decreased gradually from the original 229.91 to 166.13 mJ/m 2 .The result displays that the tendency of stacking faults decreases with the increase in Mn content.As indicated in Figure 4a, the ratio of γ isf and γ usf in Mn x CoCrFeNi HEAs was larger, demonstrating that the energy barrier needed to form the trailing part was smaller, and it was difficult to produce stacking faults.In addition, Tadmor et al. [39] gave a standard by which the main deformation mechanism can be judged when mechanical twinning becomes an ideal crack tip.By this method, the tendency of partial dislocations to form total dislocations can be determined, resulting in dislocation-mediated slip or mechanical twins.Its expression is as follows.
ing fault more likely to form.The difference between γusf and γisf in MnxCoCrFeNi was large, implying that the deformation mechanism was only extended partial d tion.Moreover, with the increase in Mn content, the difference value decreased gra from the original 229.91 to 166.13 mJ/m 2 .The result displays that the tendency of st faults decreases with the increase in Mn content.As indicated in Figure 4a, the ratio and γusf in MnxCoCrFeNi HEAs was larger, demonstrating that the energy barrier n to form the trailing part was smaller, and it was difficult to produce stacking fau addition, Tadmor et al. [39] gave a standard by which the main deformation mech can be judged when mechanical twinning becomes an ideal crack tip.By this meth tendency of partial dislocations to form total dislocations can be determined, resul dislocation-mediated slip or mechanical twins.Its expression is as follows.3 that the  values of MnxCoCrFeNi HEAs are positiv gesting that the energy barrier for unstable twin formation (γutf) is greater than the barrier for partial dislocation propagation (γusf).It can be considered that the defor mechanism of MnxCoCrFeNi HEAs is not conducive to the transformation from d tion to twins.This conclusion suggests that plastic deformation can be dominated location-mediated slip.However, with the increase in Mn concentration, its val creases.It is hinted that mechanical twin can still be formed with the increase in M centration.The existence of twin propagation starts from local dislocation or twin cr which is related to the γutf/γusf ratio.If the ratio is low, twin deformation is more li occur.It is obvious in Table 1 that the ratio of γutf/γusf has little correlation with the c in Mn concentration.For further study, Tadmor et al.'s [39] criterion was used to the tendency of FCC metal forming mechanical twins.The formula is as follows: Table 3.The calculated unstable γusf (mJ/m 2 ), stable γisf (mJ/m 2 ), unstable twinning γutf (mJ/m ning 2γtsf (mJ/m 2 ), energy barrier height  , γisf/γusf, γutf/γusf parameter, Rice-criterion D a face energies of (100), (110) and ( 111) planes (J/m 2 ) in MnxCoCrFeNi HEAs.It is clear from Table 3 that the δ ut f us f values of Mn x CoCrFeNi HEAs are positive, suggesting that the energy barrier for unstable twin formation (γ utf ) is greater than the energy barrier for partial dislocation propagation (γ usf ).It can be considered that the deformation mechanism of Mn x CoCrFeNi HEAs is not conducive to the transformation from dislocation to twins.This conclusion suggests that plastic deformation can be dominated by dislocationmediated slip.However, with the increase in Mn concentration, its value decreases.It is hinted that mechanical twin can still be formed with the increase in Mn concentration.The existence of twin propagation starts from local dislocation or twin crystals, which is related to the γ utf /γ usf ratio.If the ratio is low, twin deformation is more likely to occur.It is obvious in Table 1 that the ratio of γ utf /γ usf has little correlation with the change in Mn concentration.For further study, Tadmor et al.'s [39] criterion was used to study the tendency of FCC metal forming mechanical twins.The formula is as follows: Among them, 1.136 and 0.151 are the general coefficients of the FCC lattice.A higher τ a value can achieve higher twinning tendency.According to Figure 4b, the τ a values increase with the increase in the concentration of Mn in Mn x CoCrFeNi HEAs, but the change in τ a was little.When the Mn concentration increases from 0.0 to 1.0, the τ a value increases from 0.956 to 0.969, and the difference is only 0.013.It follows that the increase in Mn concentration can lead to twin formation, but the effect is not significant.
An alternative parameter τ, given by Asaro et al. [40], can reflect the competition between dislocation propagation and grain boundary source mechanical twins, and the expression is expressed as follows: where β = 1 − γ is f /γ us f .When τ > 1, the twin phenomenon is considered to be the dominant mechanism, which is more conducive to complete dislocation.It can be observed from Figure 4c that all the τ values of HEAs are greater than 1, and the value gradually increases from 1.164 to 1.337 with the increase in Mn concentration.This shows that the increase in Mn concentration is more favorable for twin generation.Kibey et al. [41] proposed a continuous, multi-scale method by which to predict twin stresses.The method uses the dislocation-based twin nucleation model to calculate the critical twin stress of GSFE in FCC metal, and the expression is as follows: Twin phenomena in FCC structures are caused by shear along the direction of [112] on the (111) surface, so b twin is defined as a/ √ 6.According to Equation ( 6), τ crit is determined by the four typical GSFE values (γ isf , γ usf , γ utf , and 2γ tsf ), instead of only considering the correlation of γ isf .This method has been used to predict twin stresses of some FCC metals.When the τ crit value decreases, the tendency to produce twins increases.Figure 4d demonstrates that the τ crit value decreases gradually with the increase in Mn content, from 77.83 to 23.20.It is obvious that alloys with higher Mn content require less energy for critical twinning stress than alloys with lower Mn content, and twins are more likely to be produced.
The influence of different Mn contents on the ductility of Mn x CoCrFeNi HEAs can be analyzed using the Rice-criterion [42].This analysis explains the competitive relationship between crack tip dislocation formation and crack cleavage.Its expression can be expressed as follows: where D is the ductility parameter and γ s is the surface energy along the direction of [111].When D > 1, the dissociation energy of the crack is greater than that of the dislocation nuclear energy.As a result, the alloy will exhibit ductile behavior.However, in the event of D < 0.3, the failure was caused by crack cleavage rather than dislocation slip.The calculated results listed in Figure 4e and Table 3 indicate that the D values of all alloys are greater than 1 and increase continuously with the increase in Mn concentration.The results show the Mn x CoCrFeNi HEAs' ductility, and there was dislocation formation of the crack tip.Kivy M.B. et al. [17] studied the Rice-criterion ductilities for the alloy (D > 1), indicating that the display ductility of the alloy is consistent with the calculation results in this paper.
The calculated surface energies of Mn x CoCrFeNi HEAs (111), (110), and (100) surfaces are summarized in Figure 4f.It is obvious that in HEAs with different Mn concentrations, the close-packed (111) plane has the lowest surface energy.With the increase in Mn content, the surface energy of the three planes tends to decrease.When the x value changes from 0.1 to 1.0, in Mn x CoCrFeNi HEAs, the surface energy differences of (111), (110), and (100) is 1.447 J/m 2 , 1.102 J/m 2 , and 1.071 J/m 2 , respectively.Therefore, the surface energy of surface ( 111) is most susceptible to the increase in Mn concentration.

Elastic Properties
With the goal of studying the elastic properties of Mn x CoCrFeNi HEAs, the elastic constants of the alloys were also calculated.Furthermore, the elastic modulus, Poisson's ratio, and anisotropy of the alloys are analyzed.Mn x CoCrFeNi HEAs are a cubic crystal system, which has three independent elastic constants: C 11 , C 12 , and C 44 .The change in elasticity constant of Mn x CoCrFeNi HEAs with Mn concentration is expressed in Figure 5a.As displayed in the Figure 5a, the elastic constant of the alloy conforms to the mechanical criterion of the cubic crystal system [43]: C 11 > 0; C 44 > 0; C 11 −C 12 > 0; C 11 + 2C 12 > 0. This is evidence that these alloys are mechanically stable.The three elastic constants increase with the increase in Mn content.The change in C 11 was more significant.As the value of x increases from 0 to 1, the value of C 11 increases from 237.37 GPa to 357.48 GPa.Compared with C 11 , the variation trend of C 12 and C 44 is relatively gentle.The bulk modulus (B), shear modulus (G), Young's modulus (E), Vickers har (Hv), and Poisson's ratio (σ) of the alloy can be calculated according to the followin mulae [44,45]: The bulk modulus (B), shear modulus (G), Young's modulus (E), Vickers hardness (H v ), and Poisson's ratio (σ) of the alloy can be calculated according to the following formulae [44,45]: where subscripts V and R are the approximate values obtained by Voigt and Reuss, respectively.Subscript H represents the average values of these two approximations.The bulk modulus reflects the ability of the alloy to resist bulk deformation.Shear modulus can be used to measure the shear deformation resistance of the alloy.Young's modulus can reflect the stiffness of the alloy.As can be seen from Figure 5b, with the increase in Mn concentration, the bulk modulus of HEAs increases gradually.The bulk modulus of the HEAs gradually increases from 188.23 GPa to 267.280 GPa.The change trend of Mn x CoCrFeNi HEAs' shear modulus and Young's modulus was consistent with the bulk modulus.When x is 0.7, the shear modulus and Young's moduli of HEAs reach 109.951GPa and 288.592GPa, respectively.As the concentration of Mn continues to increase, the change trend tended to be gentle.The results demonstrate that along with the raise in Mn content, the volume deformation resistance, shear deformation resistance, and alloy stiffness of HEAs increases continuously.The Vickers hardness can be macroscopic reaction alloy hardness.It can be seen from Equation (11) that Vickers hardness is affected by bulk modulus and shear modulus.It can be observed from the Figure 5b that the Vickers hardness of the alloy increases with the increase in Mn concentration.When the value of x is 0.7, the Vickers hardness of the alloy achieves 8.615 GPa.After that, with the increase in Mn concentration, the Vickers hardness of the HEAs is almost unchanged.It is obvious that when Mn content is greater than 0.7, the Vickers hardness is no longer affected by Mn content.The critical value of the Poisson's ratio, which is 0.26, is often used to judge the ductility and brittleness of alloys.When the Poisson's ratio is greater than 0.26, the compound is ductility.On the contrary, when the Poisson's ratio is less than 0.26, the compound is brittle.It can be visualized in Figure 5c that the Poisson's ratios of Mn x CoCrFeNi HEAs were all greater than 0.26, hinting that these alloys are ductile.B/G is another criterion for judging the brittleness and ductility of alloys.According to Pugh criterion [46], if the B/G value of an alloy is greater than 1.75, the material has ductility.In it is not, the material is brittle.It is presented in the Figure 5c that the values of B/G are all greater than 1.75, and their variation trend is similar to that of Poisson's ratio.This is evidence that the alloys are ductile, which is consistent with the Poisson's ratio analysis.
The anisotropy of crystal elasticity plays an important role in the study of macroscopic mechanical properties.Therefore, it is of great significance to study the influence of anisotropy of Mn x CoCrFeNi HEAs on the mechanical behavior of the alloy.The anisotropy index (A B and A G ) and the general anisotropy index (A U ) can be calculated using the following formulae [47]: When the values of A B , A G , and A U are zero, this indicates that the crystal appears isotropic.On the contrary, when the anisotropy index is not equal to zero, the crystal shows anisotropy.The greater the deviation from 0, the greater the anisotropy.Since the crystal structure of Mn x CoCrFeNi HEAs had FCC cubic structure, the values of B V and B R were equal.The calculated bulk modulus anisotropy index is zero, implying that the bulk modulus of the HEAs is isotropic.The variation trend of A G and A U values with Mn content is revealed in Figure 5d.The trend of A G and A U is consistent.The A G and A U of the CoCrFeNi HEA are the largest, which are 15.589 and 1.847, respectively.With the increase in Mn content, their values decreased slightly.When the value of x reaches 0.7, the values of A G and A U decrease to the lowest values of 7.96 and 0.853.As the x value continues to increase, the values of A G and A U increase slightly, but the change is not significant.
In order to more intuitively display the anisotropy of the alloy's elastic modulus, the bulk modulus, shear modulus, and Young's modulus of Mn x CoCrFeNi HEAs are drawn in different directions using the spherical coordinate method.The directional relation of each cubic crystal modulus can be determined via the following equations [48]: where S ij denotes the elastic compliance constants; l 1 , l 2 , and l 3 denote the directional cosines; and l 1 = sin θcos ϕ; l 2 = sin θsin ϕ; l 3 = cos θ.When the 3D diagram is spherical, this means that the alloy is isotropic.Conversely, when the 3D diagram deviates from the sphere, the alloy is anisotropic.The greater the deviation, the greater the anisotropy.As illustrated in Figure 6, the anisotropic 3D diagrams of the bulk modulus, shear modulus and Young's modulus of CoCrFeNi, Mn 0.2 CoCrFeNi, Mn 0.5 CoCrFeNi, Mn 0.7 CoCrFeNi, and MnCoCrFeNi are drawn, respectively.In Figure 6, the 3D figure of the bulk modulus is single in color and spherical, demonstrating that the bulk moduli of these HEAs are isotropic.The shear modulus and Young's modulus are all deviated from the sphere, meaning that the shear modulus and Young's modulus are anisotropic.The 3D diagram of the shear modulus of CoCrFeNi deviates most from the sphere, which clearly indicates that the shear modulus anisotropy of CoCrFeNi is the most severe.On the contrary, the deviation of the 3D diagram of the shear modulus of Mn 0.7 CoCrFeNi from the sphere is relatively small, so the shear modulus anisotropy of Mn 0.7 CoCrFeNi is small.The 3D diagram of the Young's modulus indicates a similar trend to that of the shear modulus.These results are consistent with the previous analysis results of A B and A G .where Sij denotes the elastic compliance constants; l1, l2, and l3 denote the directional sines; and l = sin θ cos ϕ ; l = sin θ sin ϕ ; l = cos θ.When the 3D diagram is spher this means that the alloy is isotropic.Conversely, when the 3D diagram deviates from sphere, the alloy is anisotropic.The greater the deviation, the greater the anisotropy illustrated in Figure 6, the anisotropic 3D diagrams of the bulk modulus, shear modu and Young's modulus of CoCrFeNi, Mn0.2CoCrFeNi, Mn0.5CoCrFeNi, Mn0.7CoCrFeNi, MnCoCrFeNi are drawn, respectively.In Figure 6, the 3D figure of the bulk modulu single in color and spherical, demonstrating that the bulk moduli of these HEAs are tropic.The shear modulus and Young's modulus are all deviated from the sphere, me ing that the shear modulus and Young's modulus are anisotropic.The 3D diagram of shear modulus of CoCrFeNi deviates most from the sphere, which clearly indicates the shear modulus anisotropy of CoCrFeNi is the most severe.On the contrary, the d ation of the 3D diagram of the shear modulus of Mn0.7CoCrFeNi from the sphere is r tively small, so the shear modulus anisotropy of Mn0.7CoCrFeNi is small.The 3D diagr of the Young's modulus indicates a similar trend to that of the shear modulus.These sults are consistent with the previous analysis results of AB and AG.

Electronic Properties
In order to study the electronic properties of MnxCoCrFeNi high-entropy alloys w different Mn concentrations, the density of states and charge density difference CoCrFeNi, Mn0.5CoCrFeNi, and MnCoCrFeNi are also calculated.The total and par density of states is plotted from −15 to 40 eV.The Fermi energy level (Ef) at 0 eV is mar

Electronic Properties
In order to study the electronic properties of Mn x CoCrFeNi high-entropy alloys with different Mn concentrations, the density of states and charge density difference of CoCrFeNi, Mn 0.5 CoCrFeNi, and MnCoCrFeNi are also calculated.The total and partial density of states is plotted from −15 to 40 eV.The Fermi energy level (E f ) at 0 eV is marked with a dashed line.As can be seen from Figure 7, the d-orbital electrons have a strong value compared with the s-orbital and p-orbital electrons in the whole range.The TDOS is mainly contributed by the d-orbital, and the peak of the TDOS is mainly distributed in the range of −5-2 eV near the Fermi level.In the case of x = 0, 0.5, and 1.0, the total densities of states at Fermi level are 45.05, 120.44, and 124.69, respectively.Obviously, the total density of states increases as the value of x increases.It can be inferred that the stability of the HEAs decrease with the increase in Mn concentration.
Materials 2024, 17, x FOR PEER REVIEW 13 of states at Fermi level are 45.05, 120.44, and 124.69, respectively.Obviously, the tota sity of states increases as the value of x increases.It can be inferred that the stability HEAs decrease with the increase in Mn concentration.In the VCA model, all atoms in a cell are the same, so the calculation of charge de difference is related to the atomic density difference in FCC cell and the VCA mod the end, the differences in electron density of MnxCoCrFeNi HEAs are simply comp Figure 8 shows a two-dimensional slice of the charge density difference of (110) pla MnxCoCrFeNi (x = 0, 0.5, and 1) HEAs.The plotting range of charge density differe from −0.1 to 0.1 e/Å 3 .Electron loss is represented in blue and electron accumulation in In Figure 8, is evident that with the increase in Mn content, the degree of charge tra gradually decreases.The weakening of atomic bonding leads to the instability of CoCrFeNi, which can increase the possibility of dislocation.In the VCA model, all atoms in a cell are the same, so the calculation of charge density difference is related to the atomic density difference in FCC cell and the VCA model.In the end, the differences in electron density of Mn x CoCrFeNi HEAs are simply compared.Figure 8 shows a two-dimensional slice of the charge density difference of (110) plane in Mn x CoCrFeNi (x = 0, 0.5, and 1) HEAs.The plotting range of charge density difference is from −0.1 to 0.1 e/Å 3 .Electron loss is represented in blue and electron accumulation in red.In Figure 8, is evident that with the increase in Mn content, the degree of charge transfer gradually decreases.The weakening of atomic bonding leads to the instability of Mn x CoCrFeNi, which can increase the possibility of dislocation.

Work Function
Generally, the electron work function (EWF) is the minimum amount of ener quired for an electron in a metal to move to its surface at the Fermi level.It is one basic electronic properties of metals, depending on their composition and surface tions.It fundamentally reflects the interaction between atoms and is directly related physical properties.In order to further understand the properties of MnxCoCrFeNi H the degrees of EWF of low exponent surfaces (100), (110), and (111) were calculated calculated trend of the EWF of MnxCoCrFeNi HEAs with a change in Mn concentra drawn in Figure 9.It can be seen that the EWF of surface (110) is the largest, follow surface (100) and surface (111).With the increase in Mn concentration, the EWF of

Work Function
Generally, the electron work function (EWF) is the minimum amount of energy required for an electron in a metal to move to its surface at the Fermi level.It is one of the basic electronic properties of metals, depending on their composition and surface conditions.It fundamentally reflects the interaction between atoms and is directly related to the physical properties.In order to further understand the properties of Mn x CoCrFeNi HEAs, the degrees of EWF of low exponent surfaces (100), (110), and (111) were calculated.The calculated trend of the EWF of Mn x CoCrFeNi HEAs with a change in Mn concentration is drawn in Figure 9.It can be seen that the EWF of surface (110) is the largest, followed by surface (100) and surface (111).With the increase in Mn concentration, the EWF of (111), (110), and (100) decreases continuously.As the EWF decreases, the atomic bonds weaken.Therefore, with the raise of Mn concentration, dislocations are more likely to occur, which is in accordance with the results of previous GSFE analysis.The research results of Shi Y.Z.et al. [15] proved that with the decrease in the work function, the alloy system containing high passivation elements (Cr, Ni, Mo, Ti, etc.) has higher corrosion resistance.There were passivation elements Cr and Ni in Mn x CoCrFeNi HEAs, and the EWF gradually decreased.It can be inferred that the corrosion resistance of Mn x CoCrFeNi HEAs heightens when accompanied by Mn concentration.The study of Lu H. et al. [49] demonstrated that the Young's modulus and hardness would increase with the increase in EWF and the strengthening of atomic bonds.It is worth noting that the work functions of Mn x CoCrFeNi HEAs have different trends with respect to the Young's modulus and hardness.

Conclusions
In this paper, the stacking fault energy, elastic properties, electronic propertie work function of the MnxCoCrFeNi high-entropy alloy are analyzed using density tional theory.The calculated results can be summarized as follows.
The results of the GSFE curves demonstrate that the alloy with higher Mn c has lower anti-dislocation nucleation ability and easily forms dislocation.The relati rier heights ( ) of the HEAs are all positive, which indicates that the deform mechanism of the HEAs is not conducive to the transformation from dislocation to The results of the GSFE and surface energies suggest that plastic deformation is con by dislocation-mediated slip.However, with the increase in Mn concentration, the of  decreases, which suggests that mechanical twinning can still be formed alloy.Based on the Rice-criterion, it is inferred that the MnxCoCrFeNi HEAs are d and have the formation of crack tip dislocation.
With the increase in Mn concentration, the bulk modulus, shear modulu Young's modulus of MnxCoCrFeNi HEAs tend to increase.It can be concluded tha the enhancement of Mn concentration, the deformation resistance, shear resistanc stiffness of the HEAs are improved.The density of states of MnxCoCrFeNi (x = 0, 0 1) HEAs at the Fermi level decreases with the increase in x, which indicates that t bility of the alloy decreases.The calculation of charge density difference demonstrat

Conclusions
In this paper, the stacking fault energy, elastic properties, electronic properties, and work function of the Mn x CoCrFeNi high-entropy alloy are analyzed using density functional theory.The calculated results can be summarized as follows.
The results of the GSFE curves demonstrate that the alloy with higher Mn content has lower anti-dislocation nucleation ability and easily forms dislocation.The relative barrier heights (δ ut f us f ) of the HEAs are all positive, which indicates that the deformation mechanism of the HEAs is not conducive to the transformation from dislocation to twin.The results of the GSFE and surface energies suggest that plastic deformation is controlled by dislocation-mediated slip.However, with the increase in Mn concentration, the value of δ ut f us f decreases, which suggests that mechanical twinning can still be formed in the alloy.Based on the Rice-criterion, it is inferred that the Mn x CoCrFeNi HEAs are ductile and have the formation of crack tip dislocation.
With the increase in Mn concentration, the bulk modulus, shear modulus, and Young's modulus of Mn x CoCrFeNi HEAs tend to increase.It can be concluded that with the enhancement of Mn concentration, the deformation resistance, shear resistance, and stiffness

Figure 1 .Figure 1 .
Figure 1.The computational cell used in the present work: (a) the atomic configuration of th cell with stacking sequence A, B, and C; (b) the supercell consisting of 13 atomic layers of (111)Figure 1.The computational cell used in the present work: (a) the atomic configuration of the FCC cell with stacking sequence A, B, and C; (b) the supercell consisting of 13 atomic layers of (111) plane.

Figure 5 .
Figure 5.The effect of Mn concentration on Mn x CoCrFeNi HEAs elastic properties: (a) elastic constants; (b) bulk modulus, shear modulus, Young's modulus, and Vickers hardness; (c) Poisson's ratio and G/B; (d) anisotropy indexes A G and A U .

Figure 6 .
Figure 6.The anisotropic characteristics of the bulk modulus, shear modulus, and Young's mod of MnxCoCrFeNi HEAs.

Figure 6 .
Figure 6.The anisotropic characteristics of the bulk modulus, shear modulus, and Young's modulus of Mn x CoCrFeNi HEAs.

Table 2 .
Calculated lattice parameter (in Å) of HEAs in the VCA model compared with previous experiment (Expt.).