## 1 Introduction  
In recent decades, a major direction from materials science focuses on understanding and applying the properties of matter to develop new multi-component alloys, with improved physical and chemical properties, based on advanced technological methods. In the literature, there are various studies on multi-component alloys (such as grain growth [1], single-phase high-entropy alloys [2], non-equiatomic high-entropy alloys [3], phase stability of equiatomic high-entropy alloys [4], tensile properties of dual-phase high-entropy alloys [5], and so on). On the other hand, it is known that multi-component alloys are considered good candidates to improve the entropy coefficient. In 2004, for the first time, nanostructured high-entropy alloys with multiple principal elements were proposed by [PERSON] et al. [6] and [PERSON] et al. [7], and for the first time, scientists made extensive investigations on the maximum value of a high entropy alloy [8]. A number of studies have suggested that high-entropy alloys have unique properties in a variety of applications [9, 10, 11]. These new findings prompted researchers to explore medium-entropy alloys (MEAs) such as CoCrNi (three-component alloys), that exhibit superior strength to CoCrFeMnNi alloy [12, 13, 14, 15]. In the literature, there are different studies aimed at the mechanical properties of low entropy alloys [6, 8, 10, 12, 14, 16]. However, so far, only a few studies have used transition metals, such as Fe [17, 18], Ni [19, 20], Al [21], and Cu [22], in magnetism,thermal, shape memory, soft elastomers, and electrical materials (such as CoCrFeMnNi that exhibits highly complex low-temperature magnetic properties [23]). Recently, scientists have identified the effects of high entropy of mixing and low mixing temperatures, associated with alloys that can form solid solutions randomly in a structural phase, such as Face-centered cubic (FCC) or Body-centered cubic (BCC). In other studies, extensive experimental investigations on high entropy alloys were conducted [24], (e.g., an integrated electromagnetic inductor with FeNiCo [25], FeNiCoCrMn, and its subset [26; 27]). Besides the results of the experimental method, there is also a simulation method that plays a great role in the success of studying the properties of new materials. Among the simulation methods, there is Density Functional Theory (DFT) methods, Monte-Carlo (MC) methods, and Molecular Dynamics (MD) simulation methods. With the specified methods, the MD method plays a key role, because of its simplicity. This method, which has been used by researchers since the 1950s [28], was only developed in 1980 thanks to the first IPM computer system. To this day, this method is strongly developed thanks to supercomputer systems combined with theoretical models of solids [29] and has been proven in various studies [30]. In studies based on the simulation method, researchers are always interested in issues such as the stability of materials, and whether the selection of parameters that will be put into the model is consistent with the experimental results. In addition, scientists are also interested in issues such as the sensitive dependence of the initial conditions leading to chaos, and the accuracy of the integration method for the system of large nonlinear partial differential equations [31]. Scientists have published various studies concerning structural characteristics, phase transition, crystallization process, and mechanical and magnetic properties of 2-component alloys (such as NiFe nanoparticles [32], NiCu [33; 34; 35], AlNi nanoparticles [36], AgAu alloy [37], AuCu bulk [38; 39], NiAu alloy [40] and polymer [41], and Fe\({}_{2}\)O\({}_{3}\)[42; 43]). In the obtained results, the authors have successfully determined the lengths of links of the Fe-Fe, and Fe-Ni.