
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39232054
71058
10.1038/s41598-024-71058-2
Article
Doping position estimation for FeRh-based alloys
Rumiantsev Egor egor.rumiantsev@epfl.ch

12
Khrabrov Kuzma khrabrov@airi.net

1
Tsypin Artem tsypin@airi.net

1
Peresypkin Nikita D. 3
Gimaev Radel R. 3
Zverev Vladimir 3
Eremin Roman 1
Kadurin Artur kadurin@airi.net

14
1 https://ror.org/014a87f14 AIRI, Kutuzovskiy Prospect 32 Bld. 1, Moscow, Russia 121170
2 https://ror.org/02s376052 grid.5333.6 0000 0001 2183 9049 École Polytechnique Fédérale de Lausanne (EPFL), Station Z, 1015 Lausanne, Switzerland
3 https://ror.org/010pmpe69 grid.14476.30 0000 0001 2342 9668 Faculty of Physics, Moscow State University, Moscow, Russia 119991
4 https://ror.org/01yqewm58 grid.26083.3f 0000 0000 9000 3133 Kuban State University, Stavropolskaya Street, 149, Krasnodar, Russia 350040
4 9 2024
4 9 2024
2024
14 2061220 5 2024
23 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
FeRh-based alloys have attracted significant attention due to their magnetic phase transition and significant magnetocaloric effects. These properties position them as promising candidates for fundamental research and practical applications, including magnetic cooling and targeted drug delivery. The study of FeRh alloys, particularly those where Rhodium or Iron atoms are substituted with other transition metals, is crucial as certain substitutions preserve the alloy’s magnetocaloric properties. However, even within a specific structural type and without considering competing phases, determining which atom (Fe or Rh) is replaced upon introducing a third element remains unclear. This paper addresses this ambiguity through ab initio calculations. We propose an approach to predict whether a dopant will replace Fe or Rh, offering insights into the electronic and structural factors influencing the substitution. Additionally, we present a dataset of ab initio calculations on doped FeRh alloys, which will support future data-driven modeling efforts. Our findings not only advance the understanding of FeRh-based alloys but also contribute to the design of novel materials for experimental and industrial applications.

Subject terms

Materials chemistry
Magnetic materials
Computational chemistry
Density functional theory
Electronic properties and materials
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Currently, FeRh-based alloys are generating considerable interest within the scientific community1–10. These alloys have recently attracted the attention of researchers because of their first-order magnetic phase transition occurring at temperatures close to those of the human body, along with their magnetic and crystalline structure properties and significant magnetocaloric effect. These properties mark FeRh-based alloys as promising materials for practical applications in fields such as magnetic cooling, magnetic recording technology, and targeted drug delivery11–16.

A number of recent works17–29 have investigated doped FeRh alloys. They demonstrate that doping FeRh alloys with Cu, Pd, and Ni in low concentrations (from approximately 1–5%) increases the refrigerant capacity of the resulting ternary alloy. However, doping these alloys may introduce specific defects where the doping occurs alongside vacancy defect. In this setting, two vacancy defects (one among Fe atoms and one among Rh atoms) occur in the crystal structure. The doping atom then substitutes one of two vacancies. Naturally, the doping atom favors one of two possible substitutions, and it is crucial to understand the underlying energetics of this process as it directly influences the properties of doped solutions30. From the computational point of view, comparing the per-atom energies of these two cases is straightforward, as the atomic composition stays constant no matter which substitution is favored by the doping atom. However, this issue has been only partially addressed in existing literature20–27, with recent studies indicating the absence of a comprehensive theory to predict which atom in the alloy is substituted by the doping element31.

From the computational perspective, the thermodynamic properties of defects and favorability of specific substitutions can be efficiently evaluated for binary compounds32,33. However, implementation of such convex-hull-based approaches still needs to be improved for ternary and more chemically complex systems due to the absence of complete lists of competing phases34. It was recently proposed28,29 to study the magnetocaloric properties of FeRh alloys as well as its doped variant FeRh1−xPdx using ab initio methods. Another work35 focuses on the magnetic and electronic properties of FePd alloy doped with Rh. However, these studies either consider pure FeRh alloys without doping elements or consider a single doping element such as Pd. Moreover, all these papers mainly focus on the magnetic properties of alloys.

In this work, we propose a unique computational scheme based on ab initio methods to study the energetics and thermodynamic properties of doped FeRh alloys where doping occurs together with a vacancy defect. We provide calculations for the full set of possible configurations for 16 different dopants and dopant concentration levels from 1.5 to 5%. Additionally, this work proposes a novel method for determining whether a dopant is more likely to substitute Rh or Fe. Specifically, we hypothesize that the favorability of substitution does not depend on dopant concentration or its position in the lattice and can be determined by considering specific dopant atom features such as ionic radius, electronegativity, and maximum valence. We demonstrate that a decision tree that operates in a three-dimensional feature space induced by the features mentioned above trained on the collected ab initio data perfectly separates all the dopants (see Fig. 4).

Finally, we compose a dataset comprising geometry optimization trajectories for doped FeRh alloys with a vacancy defect. This dataset covers various dopant concentration levels and different dopant and vacancy positions. The dataset can be utilized in future work to train data-driven models such as Graph Neural Networks36–42.

Our contributions are twofold:We introduce a method to determine the favorability of dopant substituting Rh or Fe atom upon introduction into the alloy, verified through extensive ab initio calculations.

We create a large dataset comprising geometry optimization trajectories for FeRh alloys with various dopants in different spatial configurations to facilitate the training of data-driven models.

Methods and results

This section describes how we compiled the dataset using ab initio calculations. Then, we propose a method to determine whether the dopant atom is more likely to substitute Fe or Rh from the ab initio calculations. Lastly, we show that the favorability of substitution is independent of the dopant concentration or its position in the lattice and propose a model that relies on the dopant’s atomic radius, electronegativity, and valency.

Supercells and dopants

In this part, we describe the part of the configuration space that will be studied. We start with a standard Fe–Rh cubic lattice structure. The primitive cell of this alloy is defined with orthogonal lattice vectors equal to l0≈2.99A˚ in which one Fe atom and one Rh atom are placed with coordinates (0, 0, 0) and (0.5, 0.5, 0.5), respectively, in the lattice basis. We study doped FeRh alloys where doping occurs alongside a vacancy defect. To achieve that, we remove one atom of Fe and one atom of Rh from the lattice structure and then substitute one of the vacancies with a dopant atom. This procedure ensures that the solution’s atomic composition does not depend on the substitution type (whether the atom of Rh or the atom of Fe is substituted), which allows for the straightforward comparison of the per-atom energy.

To calculate the concentration of the dopant, we use c=1N, where N is the total amount of atoms in a lattice structure excluding the vacancies. This work considers dopant concentration levels in the alloy of 1% to 5%. To obtain the desired concentration percentages, we expanded the original (l0,l0,l0) configuration by periodically repeating our initial primitive cell in space. We assume that there is only a single dopant atom and a single vacancy in each resulting configuration. For example, to obtain a dopant concentration of approximately 5% percent, we repeat the original lattice three times along one dimension and two times along two other dimensions. We call the resulting lattice (3l0,2l0,2l0) Small Supercell (SS). To get lower percentages of dopant concentration, we perform an analogous procedure and get Medium Supercell (MS) and Large Supercell (LS). The details on these configurations can be found in Table 1. Our work investigates the following M=16 dopants: Ir, Co, Cr, Ru, Pd, Pt, Mn, Ni, Te, Ti, Ta, Re, Nb, Os, Sn, Sb.Table 1 Atomic composition of considered supercells.

Supercell type	Total atoms (N)	Fe atoms	Rh atoms	Dopant atoms	Dopant concentration (c) (%)	Lattice type	
Small supercell (SS)	23	11	11	1	4.35	(3l0,2l0,2l0)	
Medium supercell (MS)	35	17	17	1	2.86	(3l0,3l0,2l0)	
Large supercell (LS)	53	26	26	1	1.89	(3l0,3l0,3l0)	

Note that the formation of interstitial solid solutions is unlikely for these transition metals since the atomic radii of these elements differ by not more than 5% to the atomic radius of Fe and Rh. Thus, according to the assumptions of Hume-Rothery rules43, all of them form substitutional solid solutions and occupy sites in the nodes of the crystalline lattice of FeRh. As indicated by previously published experimental studies, the first eight chosen components are elements that maximally preserve the significant magneto-thermal properties of binary alloys22,27,31.Fig. 1 3D representations of all possible non-equivalent double substitutions. Red atoms correspond to Fe; blue atoms correspond to Rh. The dotted bounding box represents the volume spanned by lattice vectors. Yellow atoms indicate pairs of atoms that are first removed and then successively replaced by the dopant atom.

Non-equivalent configurations

Generally, there are three different scenarios in double substitutions: we can either substitute two Fe atoms (Fe–Fe type), two Rh atoms (Rh–Rh type), or both Fe and Rh atoms (Fe–Rh type). For the reasons described in Section “Determining the replacement type”, we only consider Fe–Rh type substitutions in this work. Our goal is to study the thermodynamic properties of each symmetrically inequivalent realization of the structure, which can correspond to fixed atomic occupancies and, consequently, hinder applications of the virtual crystal44 and coherent potential approximations45.

Additionally, we obtain different resulting conformations by varying the locations from which Fe and Rh atoms are removed. Note that many double substitutions are equivalent to each other when translations and reflections are considered. To find all possible unique double substitutions, we employ Supercell software46 to obtain all non-equivalent pairs of atoms. We have two options for placing the dopant for any such pair, so the total number of double substitutions K=2∗[total number of non-equivalent pairs]. We illustrate non-equivalent double substitutions in Fig. 1. There are KSS=4 unique double substitutions for small supercell, KMS=6 unique double substitutions for medium supercell, and KLS=8 unique double substitutions for large supercell. The total number of double substitutions is thus (KSS+KMS+KLS)∗M=288.

Overall database

To complete the database, we extend it with additional structures. First, we include the baseline supercells without any substitutions, thus adding three additional conformations to the database. Second, we include the so-called “single substitutions,” in which a single Fe (Rh) atom is replaced with a dopant. The total number of single substitutions for all considered supercells is 2∗3(number of supercells)∗M=96. Additionally, it is interesting to consider the case when some of the Fe and Rh atoms are interchanged with each other. Since this procedure could be done for each combinatorially non-equivalent case, thus we increase our database by 12(KSS+KMS+KLS)=9 entities. Our dataset consists of 288+3+96+9=396 initial configurations. To get the final dataset, we obtain geometry optimization trajectories for all initial configurations (see Section “Determining the replacement type”).Fig. 2 Absolute energy differences |ΔED,SC| per atom in relaxed states for different dopants and supercell types. All the values are in milielectronvolts (meV). The red and blue colors indicate double substitutions where ED,SC<0 and ED,SC>0 respectively.

Determining the replacement type

In this section, we propose a method to determine whether the dopant atom is more likely to substitute Fe or Rh in case of the double vacancy defect through ab initio calculations (see Section “Additional details” for the details on calculations). For a given dopant type D and supercell type SC, we divide the set of all non-equivalent double substitutions into two subsets SCFe(D) and SCRh(D), where an atom of Iron and atom of Rh is substituted by a dopant atom, correspondingly. Next, for each configuration c, we calculate the per-atom energy E(c) in the relaxed state. To get the relaxed structure and the potential energy for a substitution, we perform geometry relaxation using Vienna Ab initio Simulation Package(VASP)47–49 (see Section “Additional details”). We simultaneously perform volume relaxation to account for different sizes of dopant atoms. Finally, we compare the minimum energies of two types: minc∈SCFe(D)E(c) and minc∈SCRh(D)E(c). The configuration with lower energy indicates a more desirable substitution type for the dopant. Note that the direct comparison of energies for two configurations is only possible for structures with the same atomic composition. This is the main reason we only consider double substitutions with one atom of each Fe and Rh removed. Such alloys can be denoted as Fex-1Rhx-1D1, where x is the total number of Fe(Rh) atoms in a supercell, and D denotes a dopant.

To illustrate the preferred replacement types for various dopants, we calculate the difference in minimal total energies after the relaxation between the cases where the dopant atom replaces either Fe or Rh: ΔED,SC=minc∈SCRh(D)E(c)-minc∈SCFe(D)E(c), where D stands for a dopant type, and SC stands for the supercell type. The results are summarized in Fig. 2. We color entries in Fig. 2 depending on the ΔED,SC sign. If ΔED,SC<0 (replacing the atom of Fe is more energetically preferable for the dopant atom), the color is set to purple; if ΔED,SC>0 (replacing the atom of Rh is more energetically preferable for the dopant atom), the color is set to blue. Qualitatively, we can observe that Cr, Mn, Te, Ti, Ta, Re, Nb, Sn, and Sb tend to replace Fe, while Ir, Co, Ru, Pd, Pt, Ni, and Os are more likely to substitute Rh. Moreover, for each dopant, we observe that the concentration determined by the supercell type does not influence (with precision up to 1.25 meV) the preferability of the substitution.

Rules for prediction

Fig. 3 Barplots of atomic radii, electronegativities, and maximum valencies of dopants. Dopants replacing Fe, according to the proposed computational method, are colored red, while elements replacing Rh are colored blue. Dashed lines indicate the atomic radii, electronegativities, and maximum valencies of Fe and Rh.

In this work, we are concerned with substitutional solid solutions of three metals. There is a basic set of rules50 that indicates when an element could dissolve in a metal, forming a substitutional solid solution. This set of rules operates with three numerical features of atoms: empirically measured covalent radius (also known as Slater radius51), electronegativity by Pauling scale52, and maximum valency. However, this set of rules was initially derived for two-element alloys, whereas our work studies three-element alloys. Moreover, these rules do not provide an intuition on this manuscript’s central question: whether the dopant atom is more likely to substitute Fe or Rh in case of the double vacancy defect in the alloy.

We hypothesize that the favorability of the substitution does not depend on the concentration of the dopant or its position in the lattice. We propose using the aforementioned numerical features to infer rules to determine the favorability directly from ab initio data using classification decision trees53. To train a classification decision tree, we must first map dopant atoms to vectors in the feature space defined by atomic radius, electronegativity, and valency of the atom. This procedure results in a feature matrix X∈RM×3, which is visualized in Fig. 3.

To define the targets for classification, we use ab initio calculations described in Section “Determining the replacement type” and presented in Fig. 2. We say that an element belongs to class “0” (replaces Fe) if it replaces Fe in the majority of concentrations (see Fig. 2). Analogously, we say that an element belongs to class “1” (replaces Rh) if it replaces Rh in the majority of concentrations. This results in a target vector y∈{0,1}M.

To infer the classification rules, we use the DecisionTreeClassifier from the sklean package54 and fit it using X as the feature matrix and y as the target vector. We specifically limit the maximum depth of the tree to 3 to avoid complex rules. The resulting tree is visualized in Fig. 4. Every leaf node is a terminal node with a certain class (“replaces Fe” or “replaces Rh”) assigned to it. The class of the leaf node is determined during the tree fitting as the class that the majority of the elements in the leaf have. The leaf nodes also contain the Gini impurity53 value that indicates how homogeneous the objects in the leaf are in terms of their target value. The Gini impurity of 0 means that all objects in the leaf belong to the same class. Note that the Gini impurity equals 0 in all leaf nodes, so the decision tree perfectly separates all dopants.

Every non-leaf node contains a rule consisting of a feature and a threshold. For example, the root node contains the following rule: If the dopant’s electronegativity is larger than 2.15, the dopant replaces Rh (the right child of the root node is the leaf node with class “1” assigned to it); otherwise, continue with the rule in the left child. Note that all non-leaf nodes contain rules with unique features, which indicates that all three features are essential to classify all the dopants perfectly.Fig. 4 A visualization of the classification decision tree. A leaf node contains the Gini impurity value, the total amount of samples from data X that ended up in the leaf, class values for these samples (i.e., values=[0,5] stands for 0 samples from class “0” and 5 samples from class “1”), and the class name of the majority of samples in the leaf. A non-leaf node additionally contains a rule by which to partition the samples.

Additional details

Spin-polarized total energy DFT calculations and total energy relaxations were executed using the Vienna Ab initio Simulation Package(VASP)47–49. We used a fully automatic generation scheme of the Bloch vectors with Rk=20A˚ where Rk defines subdivisions Ni along the reciprocal lattice vectors bi in the following way: Ni=int(max(1,Rk|bi|+0.5)). The energy cutoff for the plane-wave basis set was set at 600 eV. For self-consistent field conditions, the criteria for breaking the cycle were that the total energy change and the band structure energy change between two steps are smaller than 10-6 eV. Gaussian smearing was used for electronic occupancies with a smearing of 0.04 eV. All calculations were considered spin-polarized ones with initial magnetic moments equal to 2μb for each atom. Ionic relaxation was executed via the conjugate gradient algorithm until the norms of all the forces acting on the atoms were smaller than 10-2 eV/Å. We used the default blocked-Davidson-iteration scheme. Finally, the Projector Augmented-Wave (PAW)55 pseudopotentials were used to reproduce the atomic core effects in the electronic density of the valence electrons. In the standard mode, VASP performs a fully relativistic calculation for the core electrons and treats valence electrons in a scalar relativistic approximation.

Table 2 presents the averages and standard deviations of energies and maximum forces for each system, along with the number of structures analyzed. These properties characterize the systems at the final stage of ionic optimization. Notably, the average maximum force meets the criterion for halting the ionic relaxation, as the forces are smaller than 10-2 eV/Å.Table 2 Averages and standard deviations of energies per atom (eV) and maximum forces (eV/Å) for each system, alongside the number of analyzed structures at the final stage of ionic optimization.

System	Avg. Energy	Std. Dev. Energy	Avg. Max Force	Std. Dev. Force	Num. of Struc.	
SS baseline	- 7.7799	0.0000	0.000387	0.000000	1	
SS Fe–Rh exchange	- 7.7641	0.0008	0.007956	0.000455	2	
SS single substitution	- 7.7708	0.1210	0.007397	0.002083	32	
SS double substitution (1)	- 7.7350	0.1219	0.008203	0.001394	32	
SS double substitution (2)	- 7.7286	0.1246	0.008292	0.001186	32	
MS baseline	- 7.7799	0.0000	0.000566	0.000000	1	
MS Fe–Rh exchange	- 7.7723	0.0012	0.007211	0.001703	3	
MS single substitution	- 7.7712	0.0792	0.007937	0.001402	32	
MS double substitution (1)	- 7.7498	0.0804	0.008385	0.001317	32	
MS double substitution (2)	- 7.7460	0.0836	0.008477	0.001366	32	
MS double substitution (3)	- 7.7468	0.0838	0.008522	0.001234	32	
LS baseline	- 7.7799	0.0000	0.000455	0.000000	1	
LS Fe–Rh exchange	- 7.7750	0.0009	0.008095	0.000913	4	
LS single substitution	- 7.7727	0.0537	0.007196	0.002112	32	
LS double substitution (1)	- 7.7610	0.0532	0.008620	0.001066	32	
LS double substitution (2)	- 7.7550	0.0552	0.008697	0.000999	32	
LS double substitution (3)	- 7.7564	0.0552	0.008616	0.001097	32	
LS double substitution (4)	- 7.7552	0.0547	0.007517	0.001771	32	

Conclusions

In this study, we focused on FeRh-based alloys and their doped variants. We explored the issue of determining whether a dopant replaces a Rh or Fe atom. To address this, we collected a substantial dataset of ab initio calculations and developed a decision-tree-based algorithm that accurately categorizes the data. We hypothesize that the favorability of substitution depends only on the dopant’s atom type, which agrees with the computation data up to 1.25 meV. Furthermore, we compiled a comprehensive dataset of geometry optimization trajectories for FeRh-based alloys to aid in training data-driven models. We consider the exploration of the configurational space and properties of FeRh-based alloys using neural networks to be a promising direction for future research.

Author contributions

Conceptualization: E.R., K.K., R.E., V.Z., A.K. Methodology: E.R., K.K., R.E., V.Z., A.K. Software: E.R., K.K., A.T. Validation: K.K., V.Z., N.P., R.G. Investigation: E.R., K.K., V.Z., N.P., R.G. Data Curation: E.R., K.K. Writing—Original Draft: E.R., K.K., V.Z., N.P., R.G., A.K. Writing - Review and Editing: E.R, K.K., A.T., R.E., A.K. Visualization: E.R. Supervision: A.K. Project administration: A.K.

Data availability

The full data obtained and used within the manuscript is available free of charge at https://github.com/AI4DD/DOPED.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Fallot M Les alliages du fer avec les métaux de la famille du platine Ann. Phys. 1938 11 291 332 10.1051/anphys/193811100291
Fallot, M. Les alliages du fer avec les métaux de la famille du platine. Ann. Phys. 11, 291–332. 10.1051/anphys/193811100291 (1938).10.1051/anphys/193811100291
2. Mariager SO Structural and magnetic dynamics of a laser induced phase transition in FeRh Phys. Rev. Lett. 2012 108 087201 10.1103/physrevlett.108.087201 22463562
Mariager, S. O. et al. Structural and magnetic dynamics of a laser induced phase transition in FeRh. Phys. Rev. Lett. 108, 087201. 10.1103/physrevlett.108.087201 (2012).22463562 10.1103/physrevlett.108.087201
3. Gray AX Electronic structure changes across the metamagnetic transition in FeRh via hard x-ray photoemission Phys. Rev. Lett. 2012 108 257208 10.1103/physrevlett.108.257208 23004654
Gray, A. X. et al. Electronic structure changes across the metamagnetic transition in FeRh via hard x-ray photoemission. Phys. Rev. Lett. 108, 257208. 10.1103/physrevlett.108.257208 (2012).23004654 10.1103/physrevlett.108.257208
4. Cooke DW Thermodynamic measurements of Fe–Rh alloys Phys. Rev. Lett. 2012 109 255901 10.1103/physrevlett.109.255901 23368481
Cooke, D. W. et al. Thermodynamic measurements of Fe–Rh alloys. Phys. Rev. Lett. 109, 255901. 10.1103/physrevlett.109.255901 (2012).23368481 10.1103/physrevlett.109.255901
5. Derlet PM Landau–heisenberg hamiltonian model for FeRh Phys. Rev. B 2012 85 174431 10.1103/physrevb.85.174431
Derlet, P. M. Landau–heisenberg hamiltonian model for FeRh. Phys. Rev. B 85, 174431. 10.1103/physrevb.85.174431 (2012).10.1103/physrevb.85.174431
6. de Vries MA Hall-effect characterization of the metamagnetic transition in FeRh New J. Phys. 2013 15 013008 10.1088/1367-2630/15/1/013008
de Vries, M. A. et al. Hall-effect characterization of the metamagnetic transition in FeRh. New J. Phys. 15, 013008. 10.1088/1367-2630/15/1/013008 (2013).10.1088/1367-2630/15/1/013008
7. Staunton JB Banerjee R Dias MdS Deak A Szunyogh L Fluctuating local moments, itinerant electrons, and the magnetocaloric effect: Compositional hypersensitivity of FeRh Phys. Rev. B 2014 89 054427 10.1103/physrevb.89.054427
Staunton, J. B., Banerjee, R., Dias, Md. S., Deak, A. & Szunyogh, L. Fluctuating local moments, itinerant electrons, and the magnetocaloric effect: Compositional hypersensitivity of FeRh. Phys. Rev. B 89, 054427. 10.1103/physrevb.89.054427 (2014).10.1103/physrevb.89.054427
8. Song S Vertical inhomogeneous magnetic order in FeRh film Appl. Surf. Sci. 2023 607 154870 10.1016/j.apsusc.2022.154870
Song, S. et al. Vertical inhomogeneous magnetic order in FeRh film. Appl. Surf. Sci. 607, 154870 (2023).10.1016/j.apsusc.2022.154870
9. Xie Y Control of coexistent phase by rotation of magnetic field in a metamagnetic FeRh thin film J. Magn. Magn. Mater. 2023 573 170674 10.1016/j.jmmm.2023.170674
Xie, Y. et al. Control of coexistent phase by rotation of magnetic field in a metamagnetic FeRh thin film. J. Magn. Magn. Mater. 573, 170674 (2023).10.1016/j.jmmm.2023.170674
10. Qin P Antiferromagnetic spintronics: Towards high-density and ultrafast information technology Sci. Bull. 2023 68 972 974 10.1016/j.scib.2023.04.024
Qin, P. et al. Antiferromagnetic spintronics: Towards high-density and ultrafast information technology. Sci. Bull. 68, 972–974 (2023).10.1016/j.scib.2023.04.024
11. Annaorazov M Alloys of the Fe–Rh system as a new class of working material for magnetic refrigerators Cryogenics 1992 32 867 872 10.1016/0011-2275(92)90352-b
Annaorazov, M. et al. Alloys of the Fe–Rh system as a new class of working material for magnetic refrigerators. Cryogenics 32, 867–872. 10.1016/0011-2275(92)90352-b (1992).10.1016/0011-2275(92)90352-b
12. Annaorazov M Ünal M Nikitin S Tyurin A Asatryan K Magnetocaloric heat-pump cycles based on the AF–F transition in Fe–Rh alloys J. Magn. Magn. Mater. 2002 251 61 73 10.1016/s0304-8853(02)00477-8
Annaorazov, M., Ünal, M., Nikitin, S., Tyurin, A. & Asatryan, K. Magnetocaloric heat-pump cycles based on the AF–F transition in Fe–Rh alloys. J. Magn. Magn. Mater. 251, 61–73. 10.1016/s0304-8853(02)00477-8 (2002).10.1016/s0304-8853(02)00477-8
13. Thiele J-U Maat S Fullerton EE FeRh/FePt exchange spring films for thermally assisted magnetic recording media Appl. Phys. Lett. 2003 82 2859 2861 10.1063/1.1571232
Thiele, J.-U., Maat, S. & Fullerton, E. E. FeRh/FePt exchange spring films for thermally assisted magnetic recording media. Appl. Phys. Lett. 82, 2859–2861. 10.1063/1.1571232 (2003).10.1063/1.1571232
14. Annaorazov M Güven H Bärner K Cop of cooling cycles around the AF–F transition in FeRh based on experimental data J. Alloys Compd. 2005 397 26 30 10.1016/j.jallcom.2005.01.016
Annaorazov, M., Güven, H. & Bärner, K. Cop of cooling cycles around the AF–F transition in FeRh based on experimental data. J. Alloys Compd. 397, 26–30. 10.1016/j.jallcom.2005.01.016 (2005).10.1016/j.jallcom.2005.01.016
15. Tishin A Spichkin Y Zverev V Egolf P A review and new perspectives for the magnetocaloric effect: New materials and local heating and cooling inside the human body Int. J. Refrig. 2016 68 177 186 10.1016/j.ijrefrig.2016.04.020
Tishin, A., Spichkin, Y., Zverev, V. & Egolf, P. A review and new perspectives for the magnetocaloric effect: New materials and local heating and cooling inside the human body. Int. J. Refrig. 68, 177–186. 10.1016/j.ijrefrig.2016.04.020 (2016).10.1016/j.ijrefrig.2016.04.020
16. Komlev AS Gimaev RR Zverev VI Smart magnetocaloric coatings for implants: Controlled drug release for targeted delivery Phys. Open 2021 7 100063 10.1016/j.physo.2021.100063
Komlev, A. S., Gimaev, R. R. & Zverev, V. I. Smart magnetocaloric coatings for implants: Controlled drug release for targeted delivery. Phys. Open 7, 100063. 10.1016/j.physo.2021.100063 (2021).10.1016/j.physo.2021.100063
17. Myagkov V Solid-state synthesis, magnetic and structural properties of epitaxial d03-fe3rh (001) thin films Intermetallics 2023 157 107871 10.1016/j.intermet.2023.107871
Myagkov, V. et al. Solid-state synthesis, magnetic and structural properties of epitaxial d03-fe3rh (001) thin films. Intermetallics 157, 107871 (2023).10.1016/j.intermet.2023.107871
18. Merkel D Iron self-diffusion in b2-FeRh thin film Vacuum 2023 218 112617 10.1016/j.vacuum.2023.112617
Merkel, D. et al. Iron self-diffusion in b2-FeRh thin film. Vacuum 218, 112617 (2023).10.1016/j.vacuum.2023.112617
19. Negi P Formation of b2-ordered FeRh alloy thin films on annealing of pure and nitrogen doped Fe/Rh multilayers J. Magn. Magn. Mater. 2023 581 170941 10.1016/j.jmmm.2023.170941
Negi, P. et al. Formation of b2-ordered FeRh alloy thin films on annealing of pure and nitrogen doped Fe/Rh multilayers. J. Magn. Magn. Mater. 581, 170941 (2023).10.1016/j.jmmm.2023.170941
20. Uebayashi K Shimizu H Yamada H Structure and magnetism of Fe(Rh, Pd) alloys Mater. Trans. 2006 47 456 459 10.2320/matertrans.47.456
Uebayashi, K., Shimizu, H. & Yamada, H. Structure and magnetism of Fe(Rh, Pd) alloys. Mater. Trans. 47, 456–459. 10.2320/matertrans.47.456 (2006).10.2320/matertrans.47.456
21. Yuasa S Akiyama T Miyajima H Otani Y Change in the resistivity of bcc and bct FeRh alloys at first-order magnetic phase transitions J. Phys. Soc. Jpn. 1995 64 3978 3985 10.1143/JPSJ.64.3978
Yuasa, S., Akiyama, T., Miyajima, H. & Otani, Y. Change in the resistivity of bcc and bct FeRh alloys at first-order magnetic phase transitions. J. Phys. Soc. Jpn. 64, 3978–3985 (1995).10.1143/JPSJ.64.3978
22. Baranov N Barabanova E Electrical resistivity and magnetic phase transitions in modified FeRh compounds J. Alloys Compd. 1995 219 139 148 10.1016/0925-8388(94)01375-6
Baranov, N. & Barabanova, E. Electrical resistivity and magnetic phase transitions in modified FeRh compounds. J. Alloys Compd. 219, 139–148. 10.1016/0925-8388(94)01375-6 (1995).10.1016/0925-8388(94)01375-6
23. Amirov A Direct measurements of the magnetocaloric effect of Fe49Rh51 using the mirage effect J. Appl. Phys. 2020 127 233905 10.1063/5.0006355
Amirov, A. et al. Direct measurements of the magnetocaloric effect of Fe49Rh51 using the mirage effect. J. Appl. Phys. 127, 233905 (2020).10.1063/5.0006355
24. Polesya S Mankovsky S Ködderitzsch D Minár J Ebert H Finite-temperature magnetism of FeRh compounds Phys. Rev. B 2016 93 024423 10.1103/physrevb.93.024423
Polesya, S., Mankovsky, S., Ködderitzsch, D., Minár, J. & Ebert, H. Finite-temperature magnetism of FeRh compounds. Phys. Rev. B 93, 024423. 10.1103/physrevb.93.024423 (2016).10.1103/physrevb.93.024423
25. Belov MP Syzdykova AB Abrikosov IA Temperature-dependent lattice dynamics of antiferromagnetic and ferromagnetic phases of FeRh Phys. Rev. B 2020 101 134303 10.1103/physrevb.101.134303
Belov, M. P., Syzdykova, A. B. & Abrikosov, I. A. Temperature-dependent lattice dynamics of antiferromagnetic and ferromagnetic phases of FeRh. Phys. Rev. B 101, 134303. 10.1103/physrevb.101.134303 (2020).10.1103/physrevb.101.134303
26. Nishimura K Nakazawa Y Li L Mori K Magnetocaloric effect of Fe (Rh1-xPdx) alloys Mater. Trans. 2008 49 1753 1756 10.2320/matertrans.mra2008080
Nishimura, K., Nakazawa, Y., Li, L. & Mori, K. Magnetocaloric effect of Fe (Rh1-xPdx) alloys. Mater. Trans. 49, 1753–1756. 10.2320/matertrans.mra2008080 (2008).10.2320/matertrans.mra2008080
27. Komlev AS Influence of structural disorder on the magnetic order in FeRhCr alloys Metals 2023 13 1650 10.3390/met13101650
Komlev, A. S. et al. Influence of structural disorder on the magnetic order in FeRhCr alloys. Metals 13, 1650. 10.3390/met13101650 (2023).10.3390/met13101650
28. Jiménez MJ Schvval AB Cabeza GF Ab initio study of FeRh alloy properties Comput. Mater. Sci 2020 172 109385 10.1016/j.commatsci.2019.109385
Jiménez, M. J., Schvval, A. B. & Cabeza, G. F. Ab initio study of FeRh alloy properties. Comput. Mater. Sci 172, 109385 (2020).10.1016/j.commatsci.2019.109385
29. Jiménez MJ Komlev AS Gimaev RR Zverev VI Cabeza GF Electronic and thermoelectric properties of FeRh Pd-doped alloys: Ab initio study J. Magn. Magn. Mater. 2021 538 168258 10.1016/j.jmmm.2021.168258
Jiménez, M. J., Komlev, A. S., Gimaev, R. R., Zverev, V. I. & Cabeza, G. F. Electronic and thermoelectric properties of FeRh Pd-doped alloys: Ab initio study. J. Magn. Magn. Mater. 538, 168258 (2021).10.1016/j.jmmm.2021.168258
30. Komlev AS Gimaev RR Davydov AS Zverev VI The influence of chemical impurities on the properties of heavy rare-earth metals (Tb, Dy, Ho): Experimental and theoretical approaches Materialia 2021 18 101166 10.1016/j.mtla.2021.101166
Komlev, A. S., Gimaev, R. R., Davydov, A. S. & Zverev, V. I. The influence of chemical impurities on the properties of heavy rare-earth metals (Tb, Dy, Ho): Experimental and theoretical approaches. Materialia 18, 101166 (2021).10.1016/j.mtla.2021.101166
31. Komlev AS Tuning magnetocaloric effect in ternary FeRh-based alloys by slight doping Metall. Mater. Trans. A 2023 54 3683 3690 10.1007/s11661-023-07138-1
Komlev, A. S. et al. Tuning magnetocaloric effect in ternary FeRh-based alloys by slight doping. Metall. Mater. Trans. A 54, 3683–3690 (2023).10.1007/s11661-023-07138-1
32. Solokha P New quasicrystal approximant in the Sc–Pd system: From topological data mining to the bench Chem. Mater. 2020 32 1064 1079 10.1021/acs.chemmater.9b03767
Solokha, P. et al. New quasicrystal approximant in the Sc–Pd system: From topological data mining to the bench. Chem. Mater. 32, 1064–1079. 10.1021/acs.chemmater.9b03767 (2020).10.1021/acs.chemmater.9b03767
33. Eremin RA Hybrid dft/data-driven approach for searching for new quasicrystal approximants in sc-x (x= rh, pd, ir, pt) systems Cryst. Growth Design 2022 22 4570 4581 10.1021/acs.cgd.2c00463
Eremin, R. A. et al. Hybrid dft/data-driven approach for searching for new quasicrystal approximants in sc-x (x= rh, pd, ir, pt) systems. Cryst. Growth Design 22, 4570–4581. 10.1021/acs.cgd.2c00463 (2022).10.1021/acs.cgd.2c00463
34. Bartel CJ Review of computational approaches to predict the thermodynamic stability of inorganic solids J. Mater. Sci. 2022 57 10475 10498 10.1007/s10853-022-06915-4
Bartel, C. J. Review of computational approaches to predict the thermodynamic stability of inorganic solids. J. Mater. Sci. 57, 10475–10498. 10.1007/s10853-022-06915-4 (2022).10.1007/s10853-022-06915-4
35. Pavlukhina O Sokolovskiy V Buchelnikov V Zagrebin M Structural, magnetic and electronic properties of ferhxpd1-x compounds: Ab initio study Phys. B: Condens. Matter 2020 578 411882 10.1016/j.physb.2019.411882
Pavlukhina, O., Sokolovskiy, V., Buchelnikov, V. & Zagrebin, M. Structural, magnetic and electronic properties of ferhxpd1-x compounds: Ab initio study. Phys. B: Condens. Matter 578, 411882 (2020).10.1016/j.physb.2019.411882
36. Gilmer, J., Schoenholz, S. S., Riley, P. F., Vinyals, O. & Dahl, G. E. Neural message passing for quantum chemistry. In International conference on machine learning 1263–1272, PMLR (2017).
37. Chanussot L Open catalyst 2020 (oc20) dataset and community challenges ACS Catal 2021 10.1021/acscatal.0c04525
Chanussot, L. et al. Open catalyst 2020 (oc20) dataset and community challenges. ACS Catal10.1021/acscatal.0c04525 (2021).10.1021/acscatal.0c04525
38. Merchant A Scaling deep learning for materials discovery Nature 2023 624 80 85 10.1038/s41586-023-06735-9 38030720
Merchant, A. et al. Scaling deep learning for materials discovery. Nature 624, 80–85 (2023).38030720 10.1038/s41586-023-06735-9
39. Tran R The open catalyst 2022 (oc22) dataset and challenges for oxide electrocatalysts ACS Catal. 2023 13 3066 3084 10.1021/acscatal.2c05426
Tran, R. et al. The open catalyst 2022 (oc22) dataset and challenges for oxide electrocatalysts. ACS Catal. 13, 3066–3084 (2023).10.1021/acscatal.2c05426
40. Pandey BK Pandey D Parametric optimization and prediction of enhanced thermoelectric performance in co-doped camno3 using response surface methodology and neural network J. Mater. Sci. Mater. Electron. 2023 34 1589 10.1007/s10854-023-10954-1
Pandey, B. K. & Pandey, D. Parametric optimization and prediction of enhanced thermoelectric performance in co-doped camno3 using response surface methodology and neural network. J. Mater. Sci. Mater. Electron. 34, 1589 (2023).10.1007/s10854-023-10954-1
41. Na GS Jang S Chang H Predicting thermoelectric properties from chemical formula with explicitly identifying dopant effects npj Comput. Mater. 2021 7 106 10.1038/s41524-021-00564-y
Na, G. S., Jang, S. & Chang, H. Predicting thermoelectric properties from chemical formula with explicitly identifying dopant effects. npj Comput. Mater. 7, 106 (2021).10.1038/s41524-021-00564-y
42. Wang J Kwon H Kim HS Lee B-J A neural network model for high entropy alloy design npj Comput. Mater. 2023 9 60 10.1038/s41524-023-01010-x
Wang, J., Kwon, H., Kim, H. S. & Lee, B.-J. A neural network model for high entropy alloy design. npj Comput. Mater. 9, 60 (2023).10.1038/s41524-023-01010-x
43. Callister WD Materials Science and Engineering: An Introduction 2007 New York Wiley
Callister, W. D. et al. Materials Science and Engineering: An Introduction (Wiley, New York, 2007).
44. Bellaiche L Vanderbilt D Virtual crystal approximation revisited: Application to dielectric and piezoelectric properties of perovskites Phys. Rev. B 2000 61 7877 10.1103/PhysRevB.61.7877
Bellaiche, L. & Vanderbilt, D. Virtual crystal approximation revisited: Application to dielectric and piezoelectric properties of perovskites. Phys. Rev. B 61, 7877 (2000).10.1103/PhysRevB.61.7877
45. Yonezawa F Morigaki K Coherent potential approximation. Basic concepts and applications Prog. Theor. Phys. Suppl. 1973 53 1 76 10.1143/PTPS.53.1
Yonezawa, F. & Morigaki, K. Coherent potential approximation. Basic concepts and applications. Prog. Theor. Phys. Suppl. 53, 1–76 (1973).10.1143/PTPS.53.1
46. Okhotnikov K Charpentier T Cadars S Supercell program: A combinatorial structure-generation approach for the local-level modeling of atomic substitutions and partial occupancies in crystals J. Cheminform. 2016 8 1 15 10.1186/s13321-016-0129-3 26807156
Okhotnikov, K., Charpentier, T. & Cadars, S. Supercell program: A combinatorial structure-generation approach for the local-level modeling of atomic substitutions and partial occupancies in crystals. J. Cheminform. 8, 1–15 (2016).26807156 10.1186/s13321-016-0129-3
47. Kresse G Furthmüller J Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set Phys. Rev. B 1996 54 11169 10.1103/PhysRevB.54.11169
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169 (1996).10.1103/PhysRevB.54.11169
48. Kresse G Furthmüller J Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set Comput. Mater. Sci. 1996 6 15 50 10.1016/0927-0256(96)00008-0
Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 15–50 (1996).10.1016/0927-0256(96)00008-0
49. Kresse G Hafner J Ab initio molecular dynamics for liquid metals Phys. Rev. B 1993 47 558 10.1103/PhysRevB.47.558
Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 47, 558 (1993).10.1103/PhysRevB.47.558
50. Hume-Rothery W Electrons, Atoms, Metals And Alloys 1963 New York Dover Publications
Hume-Rothery, W. Electrons, Atoms, Metals And Alloys (Dover Publications, New York, 1963).
51. Slater JC Atomic radii in crystals J. Chem. Phys. 1964 41 3199 3204 10.1063/1.1725697
Slater, J. C. Atomic radii in crystals. J. Chem. Phys. 41, 3199–3204. 10.1063/1.1725697 (1964).10.1063/1.1725697
52. Pauling L The nature of the chemical bond. iv. The energy of single bonds and the relative electronegativity of atoms J. Am. Chem. Soc. 1932 54 3570 3582 10.1021/ja01348a011
Pauling, L. The nature of the chemical bond. iv. The energy of single bonds and the relative electronegativity of atoms. J. Am. Chem. Soc. 54, 3570–3582. 10.1021/ja01348a011 (1932).10.1021/ja01348a011
53. Breiman L Classification and Regression Trees 2017 Routledge
Breiman, L. Classification and Regression Trees (Routledge, 2017).
54. Pedregosa F Scikit-learn: Machine learning in Python J. Mach. Learn. Res. 2011 12 2825 2830
Pedregosa, F. et al. Scikit-learn: Machine learning in Python. J. Mach. Learn. Res. 12, 2825–2830 (2011).
55. Blöchl PE Projector augmented-wave method Phys. Rev. B 1994 50 17953 10.1103/PhysRevB.50.17953
Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953 (1994).10.1103/PhysRevB.50.17953
