Effects of Short-Range Order on the Magnetic and Mechanical Properties of FeCoNi(AlSi)x High Entropy Alloys

The properties of a material are sensitive to chemically-ordered structure in multi-element alloys. Understanding the effects of chemical short-range order (SRO) on magnetic and mechanical properties is important. In this work, we use the Monte Carlo method in combination with density functional theory to investigate atomic nearest neighbor distribution, magnetic moment and elastic modulus in FeCoNi (AlSi)x alloys. It is found that the prominent feature of the FeCoNi (AlSi)x alloys is the change of SRO parameters: the SRO parameters are positive between Al-Al, Al-Si, Si-Si pairs and negative between Ni-Al, Co-Si, Fe-Co, Ni-Si and Fe-Si pairs. The Al and Si elements tend to bond with Fe, Co, Ni elements to form an SRO structure. The change of the atomic nearest neighbor environment leads to a reduction in the atomic magnetic moments of magnetic elements. The calculated saturation magnetizations by considering the effect of SRO are in good accord with the experimental values. We further show that SRO leads to an increase of the elastic modulus, by sacrificing ductility and isotropy. In the study of the structure and properties of high entropy alloys, the effect of SRO should not be ignored.


Introduction
In recent years, there have been extensive experimental investigations of multi-component alloys after the new alloy design concept, high entropy alloys (HEAs), was proposed [1,2]. HEAs consist of at least five principal elements, each element with a 5-35% concentration. These alloys have many exceptional properties compared with the traditional alloys due to the high mixing entropy, sluggish diffusion and lattice distortion effects, including good thermal stability, high mechanical strength and excellent corrosion resistance [3][4][5]. The properties of HEAs, such as hardness, ductility and magnetic properties can be tuned by changing the type and concentration of components [6][7][8]. These novel properties of HEAs have potential applications in refractory materials, wear-resistant materials and magnetic materials.
In general, HEAs are considered to be random solid solutions (RSS), and all constituent atoms randomly distribute on the available lattice sites. But is this really true? In fact, the diversity in atomic radii and the different attractive interactions between the constituent elements will result in a short-range ordered structure [9]. The SRO structures have been found in many engineering materials such as FeV, FeAl and FeCr alloys [10][11][12]. It will significantly affect the structural stability, magnetic and mechanical properties. The presence of SRO in these binary alloys suggests that some and BCC special quasi-random structures (SQS) of FeCoNi(AlSi) 0.8 alloy were generated by using the mcsqs code with the Alloy-Theoretic Automated Toolkit (ATAT) package [28].
The structure relaxation and energy calculation were performed using DFT implemented with the Vienna Ab-Initio Simulation Package (VASP, University of Vienna, version: 5.4.1.05, Vienna, Austria) [21]. The exchange-correlation functional was treated within the generalized gradient approximation as described by Perdew-Burke-Ernzerhof (PBE) [29]. Energy cutoffs for the plane wave in all calculations were set to 300 eV. The reciprocal space energy integration was performed by the Monkhorst-Pack technique with 2 × 2 × 2 mesh. The self-consistence convergence criterion for electron iterations was set to 10 −5 eV per atom. The valence electrons of Fe, Co, Ni, Al and Si atoms were eight, nine, 10, three and four, respectively. From Wikipedia, we know that the Curie temperature of Fe, Co and Ni are 1043, 1388 and 628 K. Körmann et al. predicted the Curie temperature of FeCoNi is 868 K [30]. The Curie temperature of Co 2 FeSi and Co 2 FeAl are all above 1000 K. Dobrzynski et al. reported that the Curie temperature of Fe 2.94 Al 0.38 Si 0.68 is 794 K [31]. There is not any Curie temperature of FeCoNi (AlSi) x alloys reported. However, Zhang et al. showed that the saturation magnetizations of FeCoNi (AlSi) x alloys range from 1.32 to 0.46 T (x = 0 to 0.8) [22]. We predict that the Curie temperature of FeCoNi (AlSi)x alloy is above 300 K. Thus, the magnetic state must be considered. All the calculations done were spin-polarized, with the magnetic moments of Fe, Co and Ni atoms initialized with 3, 2 -and 1 µ B , respectively. The calculations in this work did not include the orbital magnetic moment.
In order to probe SRO in FeCoNi (AlSi) x alloys, the Monte Carlo method was carried out in the simulation. The generation of a new atomic configuration is a two-step process. Firstly, the potential energy with the current atomic structure was calculated. In the next step, two different kinds of constituent atoms were selected and the potential energy after the swap of their positions was calculated. The position swap is based on the Metropolis-Hastings sampling [32] with the acceptance probability p = exp(−β∆E); ∆E is the energy difference between swapped and not-swapped structures, β = k B T, T = 300 K. The FeCoNi (FCC), FeCoNi(AlSi) 0.2 (FCC), FeCoNi(AlSi) 0.4 (BCC) and FeCoNi(AlSi) 0.8 (BCC) structures generated by the MaxEnt method, and FeCoNi(AlSi) 0.8 (BCC) generated by the mcsqs method, were selected as the initial starting point in the MC simulations. Each MC simulation ran for 3000 MC steps. The structure relaxations and energy calculations for the final structures with SRO and the initial structures were fully conducted with a higher number of k-points (3 × 3 × 3) to increase accuracy.
The Warren-Cowley SRO parameter α ij [33] was used to described the degree of short-range order in the FeCoNi(AlSi) x alloys. The nearest neighbor SRO parameter α ij can be calculated by Equation (1), where c j is the molar fraction of the type j element and p ij is the probability of finding the type j element around the type i element in the nearest neighbor shell. The SRO parameter vanishes if p ij = c j , meaning that there is no site preference between the type i and type j element. The negative α ij indicates the increase in the number of i and j pairs, while the positive value corresponds to the opposite.
Formation enthalpy E f is defined in Equation (2), where E tot is the zero-temperature total energy per atom, c i is the molar fraction of the type i element. E re f i is the ground state energy of type i element in its stable phase and is calculated by using the same DFT code and pseudopotential.

Single-Phase High Entropy Alloy (HEA) Stability
It was reported that non-equimolar FeCoNi(AlSi) x alloys form the FCC phase for 0 ≤ x ≤ 0.2 and the BCC phase for 0.4 ≤ x ≤ 0.8 [22]. In Figure 1, we present the MaxEnt structures of FeCoNi (AlSi) x for x = 0.2 (FCC), x = 0.4 (BCC) and x = 0.8 (BCC) and SQS structure with x = 0.8 (BCC) in their ideal, unrelaxed forms. We performed the structure relaxation and total energy calculation by using the settings described in Section 2. Table 1 shows the lattice parameters, formation enthalpies E f with BCC and FCC structures and the energy difference between the BCC and FCC phases The positive ∆E t indicates that the FCC structure is more stable, and the negative ∆E t shows that the BCC structure is more stable. As shown in Table 1, when the fraction of Al and Si elements is smaller (x ≤ 0.2), the energy difference ∆E t is positive, indicating that the FCC structure is more stable. When the fraction of Al and Si elements is larger (x ≥ 0.4), the energy difference ∆E t is negative, demonstrating that the BCC structure is more stable. The experimental lattice parameters are a = 3.571 Å for x = 0 [34] and a = 3.594 Å for x = 0.2 [14]. The calculated results are a = 3.545 Å (x = 0) and a = 3.554 Å (x = 0.2). When x = 0.8, there is no experimental lattice parameter available, however, the calculated lattice parameter a = 2.853 Å with the MaxEnt structure and a = 2.874 Å with the SQS structure were obtained. The consistency of the calculated and experimental values validates the MaxEnt models employed.

Short-Range Order (SRO) in FeCoNi(AlSi) x Alloys
The evolution of the relative potential energy and SRO parameters of FeCoNi (AlSi) 0.8 alloy are shown in Figure 2. The relative potential energy is defined as the energy difference between the potential energy of the ith MC step and the initial potential energy. In the case that the 92-atom SQS structure is selected as the initial starting structure, all the constituent atoms are randomly distributed in the lattice positions, α ii = 0 and α ij = 0 (Figure 2a,b). Figure 2b shows that there are major deviations of SRO parameters between the finial structure and the SQS structure. The average numbers of Al-Al and Si-Si pairs are reduced by~80%, the nearest neighbor number changes from 2.1 for the random case to about 0.5 for the SRO case, and the average number of Al-Si pairs is reduced by~50%, indicating Al (or Si) atoms separate from each other. Meanwhile, it can be seen that there is rapid growth of the average number of Ni-Al pairs, followed by Co-Si, Co-Al and Ni-Si pairs. The rapid growth of Ni-Al pairs shows that the Ni-Al pair first formed during the solidification from the liquid phase. The Ni-Al pairs will form a B2-orderd parent crystal structure, which is an ordered structure based on BCC with Pearson symbol of cP2. Si element also have negative mixing enthalpies with Ni, Co and Fe elements, the binary mixing enthalpies of Ni-Si, Co-Si and Fe-Si pairs are −40, −38 and −35 kJ/mol, respectively. Si atoms prefer to separate with Al and Si atoms, and bond with Ni, Co and Fe atoms. Raja et al. studied the structural properties of Fe 3-x Co x Si alloys (0 ≤ x ≤ 1) using X-ray powder diffraction [35]. The results revealed that Si atoms occupy the body-centered positions, Fe and Co atoms occupy the body corners and form B2 and L2 1 ordered phases. The L2 1 phase is an ordered structure based on BCC with Pearson symbol of cF16. Just like Al, Si also tends to stabilize BCC structure. In the FeCoNi (AlSi) 0.4 alloy, although the fraction of Al element (x = 0.4) is smaller as compared to the case in Al x CoCrFeNi alloy (BCC phase in x ≥ 1.25 ), there will be a large amount of Ni-Al and Co-Si pairs first formed and other constituents will dissolve in the parent crystal structure due to the mixing entropy effect. Thus, the system forms a BCC structure. In the case where the MaxEnt structure is selected as the initial starting structure (Figure 2c,d), since each constituent atoms in the system strives for the maximum free space, α ii = 0.92 and α ij = −0.22 for FeCoNi(AlSi) 0.8 alloy. We obtained the lower potential energy at the 2200 MC step. After 800 MC steps searches, the MC method did not find the lower energy structure. The emergence of the SRO reduced the potential energy, which made the potential energy change from −325 to −359 meV per atom. In Figure 2d, the SRO parameters of α SiSi and α Al Al decrease slightly. α SiSi changes from 0.92 to 0.86, and α Al Al changes from 0.92 to 0.84. We started from two different initial starting points to probe the change of the atomic nearest neighbor environment in FeCoNi (AlSi) 0.8 alloy, and obtained a similar trend to the SRO. It can be concluded that the potential energy curve is relatively converged after 3000 MC steps simulation. From the above comparison, it can be seen that in the process of searching for the SRO structure, using the MaxEnt model as the initial starting structure is more efficient. This is based on the following two points: on the one hand, the potential energy of the MaxEnt structure is lower than that of the SQS structure, which can be validated from the difference in formation energy listed in Table 1. On the other hand, the SRO parameters α ii and α ij of the MaxEnt structure are closer to the SRO parameters of the finial SRO structure. Therefore, in the following probing of SRO, and the calculation of the magnetic and mechanical properties, the MaxEnt structure is selected as the initial structure.
The partial pair distribution function can be used to describe the relative positional preferences of different constituent elements. The shape and position of the peaks provide detailed information about the atomic local environment [36]. We used pair PDF to show the SRO structure of the FeCoNi(AlSi) 0.2 alloy, and the corresponding data are plotted in Figure 3. We observe that there is no obvious lattice distortion in the alloy. The average bond lengths of Co-Si and Ni-Si are r~2.45 Å, which is shorter than the metal-metal bonds because the Si atom has a smaller radius than the metal atoms (Fe, Co, Ni, and Al). The Ni-Al pair has a higher peak intensity, followed by the Co-Si and Fe-Co, Ni-Si pairs. The peaks of the Al-Al and Al-Si pairs have moved to the second nearest neighbor shell, and the peak of the Si-Si pair has moved to the third nearest neighbor shell, indicating that Al and Si prefer to bond with other elements rather than themselves. The results suggest the existence of the preferred Ni-Al, Co-Si, Fe-Co and Ni-Si pairs in the FeCoNi(AlSi) 0.2 alloy. The partial PDF of the FeCoNi(AlSi) 0.8 alloy after the 3000 MC steps is shown in Figure 4. It can be seen that the most preferred pair is Ni-Al, followed by the Co-Si, Co-Al, Fe-Si and Ni-Si pairs. The least favored pairs are the Si-Si, Al-Al pairs. Al and Si atoms prefer to bond with other elements to lower the potential energy. The results suggest the existence of a preferred short-range order of Ni-Al, Co-Si, Co-Al, Fe-Si and Ni-Si pairs in the alloy. In addition, there is obvious lattice distortion when the fraction of Al and Si atoms is larger. The lattice distortion makes the second intensity peak smear out and become less distinct, indicating that the atoms in the alloy deviate from the ideal lattice positions. The lattice distortion can significantly scatter free electrons, shortened electrons mean free paths, which will reduce the thermal and electrical conductivity of the alloy. In Table 2, the SRO parameters of the FeCoNi(AlSi) x alloys as a function of Al and Si fraction x averaged after 3000 MC steps are shown. For FeCoNi alloy, we observed that the SRO parameter of the Fe-Ni pair is negative, and the Fe-Fe, Ni-Ni and Co-Co pairs are positive. Fe atoms prefer to bond with Ni atoms and form an SRO structure. Tamm et al. investigated the SRO behavior of NiCrCoFe alloys and found a negative Ni-Fe pair and positive Fe-Fe, Ni-Ni and Co-Co pairs [16]. When Al and Si atoms are added to the FeCoNi-based alloy, the Fe-Ni pair become positive. Al and Si elements have relative stronger attractive interaction with other constituent elements. They will separate from each other and bond with Fe, Co and Ni elements to form an SRO structure. The strength of the attractive interaction can be seen from the binary mixing enthalpies [37], that is, the more negative the binary mixing enthalpy, the stronger the attractive interaction between the binary pair. This is applicable to other HEAs as well. The SRO depends on the type of constituent elements present in the HEAs. If a constituent element has a stronger attractive interaction with the other constituent elements, the atoms of this element will separate from each other and bond with other constituent elements to form an SRO structure. Also, the concentration of the constituent elements can change the degree of SRO. We observed that the SRO parameter of Co-Al changes from 0.34 (x = 0.2) to −0.56 (x = 0.8) and Co-Co changes from 0.15 (x = 0.2) to 0.84 (x = 0.8). It indicates that the element concentrations also have a certain influence on the degree of SRO. We think that a possibility for the formation of SRO in the HEAs is the enrichment of the preferred element pairs during solidification. At high temperature, the atomistic structure of liquid alloy is generally thought to be totally random due to the high mixing entropy effect. As the temperature decreases, the enthalpy effect becomes more important, and the non-random configuration shows a tendency toward phase separation or chemical short-range order, and some ordered phases may form during the solidification. Santodonato et al. studied the structural evolution of Al 1.3 CoCrCuFeNi alloy from the high temperature liquid phase to the room temperature phase [38]. The results demonstrated that the alloy is a liquid above 1315 K with Al-Ni, Cr-Fe and Cu-Cu the preferred nearest-neighbor pairs. The results demonstrated that the alloy is a liquid above 1315 K with Al-Ni, Cr-Fe and Cu-Cu preferred nearest-neighbor pairs. During the cooling of the melt, the ordered phases will be preserved. Due to the preferred Ni-Al, Co-Al and Fe-Si pairs in the BCC phase of FeCoNi (AlSi) 0.4 and FeCoNi (AlSi) 0.8 alloys, a B2-ordered phase structure may be formed. In the multi-component ordered BCC alloys, the crystal structure can be divided in two interpenetrating sublattices, designed by an α sublattice, and a β sublattice. x iα and x iβ denote the molar fraction of the ith element on α and β sublattices. The B2-ordered parameter η i was used to quantify the degree of ordering of the ith element in the sublattices [38]. η i is described using: The total order parameter η is calculated from η i with the equation: where x i is the overall molar fraction. η = 0 indicates the random solid solution, and η = 1 shows the fully B2-ordered structure. The mixing entropy for the case of a B2 structure can be calculated: In Table 3, we presented the ordering parameter η i , η and the mixing entropy for the ideal, and taking into account the B2-ordered structure. For the FeCoNi (AlSi) 0.4 alloy, η Al = 0.692, η Si = −0.231 and η Ni = −0.294, indicating that Al prefers to occupy the α sublattice and Ni and Si prefer the β sublattice to form the B2 phase. The total B2-ordered parameter was η = 0.294 when x = 0.8, η Al = 0.652, η Co = −0.428 and η Ni = −0.481, indicating that Al prefers to occupy the α sublattice and Ni and Co prefer the β sublattice to form the B2 phase. With the increase of Al and Si content, the total B2-ordered parameter η reaches 0.427. The emergence of the B2 phase in the alloys reduces the configuration entropy from 1.521 to 1.475 R for the FeCoNi (AlSi) 0.4 alloy, forming 1.606 R to 1.509 R for the FeCoNi(AlSi) 0.8 alloy.

SRO on Magnetic Properties
As a reference, a series of calculations were carried out to determine the magnetic moments of Fe with BCC, Ni with FCC and Co with HEX structure, respectively. The calculated magnetic moments per Fe, Co and Ni atom are 2.18 µ B , 1.67 µ B , and 0.62 µ B , which are in good agreement with the experimental values [39]. The deviations between the calculated and the experimental lattice parameters (a) are within 1%. Chandran et al. [40] investigated Fe 1−x Co x alloys in BCC structures and found that with the increase of Co contents, the magnetic moment per Fe atom increased from 2.22 µ B (x = 0) to 2.76 µ B (x = 0.5). Apiñaniz et al. [41] studied the magnetic properties of ordered Fe x Al 1-x alloys and found that the magnetic moment per Fe atom decreases from 2.22 to 0.64 µ B with the increase in Al contents. The FeM, CoM and NiM binary ordered alloys (M = Fe, Co, Ni, Al and Si) with BCC primitive cells were constructed. Taking the FeM alloy as an example, the Fe atom was assigned to the body-centered position and the M atom to the body corners. Table 4 listed the calculated magnetic moments per Fe, Co and Ni atoms with the M atom in the nearest neighbor shell. It can be seen that Co and Ni atoms in the nearest neighbor shell of Fe atoms greatly increase the magnetic moment of Fe atoms, while Al and Si atoms in the nearest neighbor shell of magnetic atoms will drastically decrease the magnetic moments of magnetic atoms. Therefore, the atomic nearest neighbor environment has considerable influence on the atomic magnetic moments in alloys. There is obvious SRO behavior in the FeCoNi (AlSi) x alloys. The SRO will significantly change the atomic nearest neighbor environment, which has an impact on the magnetic properties of the alloys. The saturation magnetizations and average atomic magnetic moments of the FeCoNi (AlSi) x alloys with MaxEnt and SRO structures are summarized in Table 5. For the FeCoNi alloy, the average atomic magnetic moments per Fe, Co and Ni atoms are 2.65, 1.64 µ B and 0.62 µ B with the initial structure. The Co and Ni atoms in the nearest neighbor of Fe atoms increase the magnetic moment of the Fe atoms. The effect of SRO on the magnetic moments of the FeCoNi alloy is relatively small. The calculated magnetic moments (M s = 1.65 T) are consistent with the previous calculated values (M s = 1.62 T) [22]. For the FeCoNi(AlSi) 0.2 alloy, the SRO parameter of the Ni-Al pair reduces by −100% and the nearest neighbor number of the Ni to Al atom changes from 3.5 to 7 for the SRO structure. Similarly, the SRO parameters of the Co-Si, Ni-Si and Fe-Co pairs are also negative. According to the atomic coordinate analysis, Al atoms occupy the body corners, while Ni atoms occupy the face-centered positions and form an Ni-Al rich region. The negative SRO parameters of Ni-Al, Ni-Si and Fe-Co pairs present a similar trend with the experimental work [14] that the dendritic area is rich in Fe and Co atoms, while the inter-dendritic area is rich in Al, Ni, and Si atoms in FeCoNi (AlSi) 0.2 alloy. In Figure 5, the average atomic magnetic moment µ i,j of the FeCoNi (AlSi) 0.2 and FeCoNi (AlSi) 0.8 alloys are presented. µ i,j is the average magnetic moment of the type i element, which is calculated that there is type j element in its nearest neighbor shell. The average magnetic moments µ Fe,j , µ Co,j and µ Ni,j vary with different nearest neighbor elements. The results are different from the case of random solid solution. If each lattice is randomly occupied and the probability is proportional to the concentration of the element, the magnetic moment will be homogeneous. The SRO behavior decreases the atomic magnetic moments and the saturation magnetization of the system. The average magnetic moment of the Fe atom changes from 2.50 to 2.32 µ B , the Co atom from 1.43 to 1.19 µ B and the Ni atom from 0.51 to 0.32 µ B . Compared to the previously-calculated saturation magnetization value M s = 1.31 T with the random structure model [22], the saturation magnetization with SRO is M s = 1.17 T, which is closer to the experimental value M s = 1.15 T [22]. Table 5. The saturation magnetizations M s , the errors of calculated saturation magnetizations to experimental values and average magnetic moments µ i for FeCoNi(AlSi) x alloys with the MaxEnt structures and the finial structures with SRO, the calculated and experimental data from previous publications are also shown.

Alloy
Ordering  The FeCoNi(AlSi) 0.8 phase contains more Al and Si components. The SRO parameters of Ni-Al, Co-Si, Co-Al, Fe-Si and Ni-Si pairs are more negative. It is obvious that the Fe, Co and Ni atoms prefer to gather around the Si atom, while the Ni and Co atoms prefer to gather around the Al atom. The SRO greatly alters the local environment of the magnetic atoms, which will further reduce the average magnetic moments of the magnetic atoms. The average magnetic moment of the Fe atoms changes from 2.10 to 2.03 µ B , Co atoms from 0.94 to 0.76 µ B . Due to the enrichment of Al and Si atoms in the nearest neighboring shell of Ni atoms, the magnetic moment of Ni atoms vanishes. The saturation magnetization changes from 0.67 T with the random structure to 0.6 T with the SRO structure. The saturation magnetization with SRO is closer to the experimental value = 0.46 T [22].

SRO on Elastic Properties
Mechanical properties are a vital aspect in material selection. It is important to understand the relationship of component, structure and mechanical properties for material applications. FeCoNi(AlSi) x alloys form FCC or BCC phase structures with different Al and Si fractions. There are three independent elastic constants c 11 , c 12 , and c 44 for the present cubic lattice. c 11 and c 12 can be determined from the bulk modulus (B) and tetragonal shear modulus (c ) with B = (c 11 + 2c 12 )/3 , and c = (c 11 − c 12 )/2. The bulk modulus can be obtained by fitting energy-volume data with the three-order Birch-Murnaghan equation of state [42]. Tetragonal shear modulus can be extracted from ∆E(δ 0 ) = 2V c δ 2 0 by applying orthorhombic strain (δ 0 ) to the cubic lattice. The elastic constant c 44 can be obtained from fitting the energy-strain equation: ∆E(δ m ) = 2Vc 44 δ 2 m by applying the monoclinic strain (δ m ) to the base lattice [43]. The shear strain, strain matrix and energy-strain equations are shown in Table 6. In order to keep the elastic behavior of crystals, the applied strains should be relatively small, so the strains adopted are δ = −0.009, −0.006, −0.003, 0, 0.003, 0.006, 0.009. The MaxEnt structure is a cubic lattice structure. We modified the basis vector matrix according to the strain matrix to apply different deformations to the supercell structures. The calculations of the total energy were conducted with 3 × 3 × 3 k-points to increase the accuracy. The energy curves of different deformations for the FeCoNi(AlSi) x alloys are shown in Figure 6. Table 6. The shear strain, strain matrix and the corresponding energy-strain equation to calculate c 11 , and c 44 for the cubic structure. The shear modulus G is calculated by the Hill average G = (G V + G R )/2, and the Voigt and Reuss bounds shear modulus can be obtained by using Equations (6) and (7):

Shear Strain Strain Matrix Energy-Strain Equation
Young's modulus (E) and Poisson's ratio (ν) are calculated from the bulk modulus (B) and the shear modulus (G): The calculated lattice parameters, three elastic moduli B, G, E and derived elastic moduli of the FeCoNi(AlSi) x alloys with the MaxEnt and SRO structures are presented in Table 7. From Table 7, all the calculated elastic constants (c 11 , c 12 and c 44 ) fulfill the mechanical stability criteria: c 44 > 0, c 11 > |c 12 | and (c 11 + c 12 ) > 0, which demonstrates that the FeCoNi (AlSi) x alloys are mechanically stable. The experimental information on the elastic modulus is very limited. The calculated lattice parameter for the FeCoNi alloy is a = 3.545 Å, and the experimental value is 3.599 Å [22]. The calculated lattice parameter is closer to the value a = 3.563 Å calculated with EMTO [44]. The calculated bulk modulus is B = 184.7 GPa, shear modulus G = 84.5 GPa, and Young's modulus E = 219.8 GPa.
The corresponding results calculated with EMTO are B = 186.1 GPa, G = 90.2 GPa and E = 233.0 GPa. The bulk modulus is a measure of the resistance to compressibility of a material. From the calculations, it can be seen that the fractions of Al and Si have significant influence on bulk modulus of FeCoNi (AlSi) x alloys. It is observed that the elastic constants (c 11 , c 12 and c 44 ) and the three elastic moduli (B, G and E) decrease with increase of Si and Al fractions. The bulk modulus changes from 184.7 GPa (x = 0) to 157.1 GPa (x = 0.8). Note that Al and Si have a smaller bulk modulus compared to Fe, Co and Ni components. Also the addition of Al and Si make the alloy structure change from the FCC to BCC phase (the packing factor changes from 0.74 to 0.68), which will result in a decrease of the bulk modulus. Taking into account the effect of SRO, SRO behavior prefers to make the affinity atoms locate together, resulting in an increase of the average bonding strength and a decrease of alloy volume. The FeCoNi(AlSi) x alloys with SRO have a higher bulk modulus as compared to that of random solid structures. We also observed that SRO has similar effects on shear modulus and Young's modulus.
The Pugh's ratio (B/G) [45] and Poisson's ratio (ν) [46] can be used to qualify the ductile properties of a material. It was reported that an alloy material is ductile when B/G > 1.75 and ν > 0.31, otherwise it is brittle. For FeCoNi(AlSi) x alloys, the Cauchy pressure ( c 12 − c 44 ), Pugh's ratio and Poisson's ratio increase with the increase of Al and Si fractions. The trends of Cauchy pressure, Pugh's ratio B/G and Poisson's ratio νindicate the addition of Al and Si to FeCoNi alloys improves the ductility of alloys. However, it is observed that the SRO behavior reduces the Pugh's ratio, Poisson's ratio and Cauchy pressure, showing that SRO results in a reduction in the ductility of the material. The Zener ratio A z and the ratio A VR are used to describe the isotropy property of FeCoNi (AlSi) x alloys. A z = 2c 44 /(c 11 − c 12 ), and A VR = (G V − G R )/(G V + G R ). For an isotropic cubic material, A z = 1 and A VR = 0. The values of A z and A VR indicate the relative degree of the elastic anisotropy. For FeCoNi(AlSi) x alloys, the Zener ratio A z changes from 2.94 to 4.29 and A VR changes from 0.13 to 0.23 with the increase of Al and Si fractions. The Al and Si enhance the anisotropy of the FeCoNi (AlSi) x alloys. When taking into account the effect of SRO, the A z changes from 2.99 to 4.36 and the A VR changes from 0.14 to 0.24. The SRO behavior further enhances the anisotropy of the FeCoNi (AlSi) x alloys. From the above results, it can be seen that although SRO has a positive effect on B, G and E, it has a negative effect on the ductility and isotropic properties of the materials.

Conclusions
Using Monte Carlo simulation in combination with density functional theory, we investigated the SRO behavior of FeCoNi-based high entropy alloys. We took the SQS and MaxEnt structures as the initial starting point to probe the SRO structure in the system and found that using the MaxEnt structure as the initial starting point was more efficient. The calculated results show that there are obvious SRO structures in the FeCoNi (AlSi) x alloys. The binary mixing enthalpies between constituent elements were found to be the key factor in controlling the formation of SRO in a material. Al and Si elements have more negative binary mixing enthalpies with Fe, Co and Ni elements. The average numbers of Al-Al, Al-Si and Si-Si pairs decrease significantly, while those of Ni-Al, Co-Si, Fe-Si, Ni-Si and Fe-Co pairs correspondingly increase. The changes in the SRO parameters indicate that Al and Si tend to bond with Fe, Co, and Ni atoms to lower the potential energy. The element concentrations also have a certain influence on the degree of SRO. The emergence of SRO will change the properties of high entropy alloys, such as reducing the formation energy, altering the equilibrium volume of the system and decreasing the average atomic magnetic moments. The SRO further reduces the saturation magnetization and damages the ductility and isotropy. The present work shows that SRO behavior has an important influence on the magnetic and mechanical properties of high entropy alloys.