## Theoretical prediction of a highly responsive material: spin fluctuations and superconductivity in FeNiB 2 system

Renhai Wang, 1, 2 Yang Sun, 2, ∗ Vladimir Antropov, 2, † Zijing Lin, 1 Cai-Zhuang Wang, 2, 3 and Kai-Ming Ho 2, 3, 4

1

Ames Laboratory, US DOE, Ames, Iowa 50011, USA

Department of Physics, University of Science and Technology of China, Hefei 230026, China 2

3 Department of Physics, Iowa State University, Iowa 50011, USA

4 International Center for Quantum Design of Functional Materials (ICQD),

Hefei National Laboratory for Physics Sciences at the Microscale,

University of

Science and Technology of China,

Hefei

230026,

China

(Dated: Oct. 20, 2019)

By analyzing Fe-Ni-B compositional diagram we predict an energetically and dynamically stable FeNiB 2 compound. This system belongs to the class of highly responsive state of material, as it is very sensitive to the external perturbations. This state is also characterized by a high level of spin fluctuations which strongly influence possible magnetic long- and short-range orders. Furthermore, we demonstrate that these antiferromagnetically dominating fluctuations could lead to the appearance of spin mediated superconductivity. The obtained results suggest a promising avenue for the search of strong spin fluctuation systems and related superconductors.

Introduction Spin fluctuations (SF) in itinerant electron magnets play an important role in many metallic systems including weak itinerant magnets, heavy fermion compounds, actinides, Invar alloys and many magnetoresistive materials. SF also have been identified in high-temperature superconductors [1] and recently discovered Fe-based superconductors [2] with a clear suggestion about their crucial influence on the mechanisms of unconventional superconductivity. In the majority of cases superconductivity appears near magnetic quantum critical point, where magnetism has a pure itinerant character with no local moments involved. This is exactly a case when quantum SF (including spin zero-point motion) are strong and provide a dominating contribution to many observable physical properties. However, as a typical example, in a case of superconductivity, the direct calculations and prediction of critical temperature from first principle is still not possible. In this situation one can focus on searching for new materials with a strong level of SF, leading to possible development of promising physical properties.

Here we will introduce qualitatively the concept of a highly responsive state (HRS) (see Fig. 1). This term is used here to indicate that a physical system is in a state characterized by an order parameter which is very sensitive to the external perturbations (temperature, pressure, magnetic fields, etc.). The HRS state is naturally related with a closeness to the instability and, as usual in such situations, one can expect that fluctuations are very important.

To search for new material of such type we start our study with a less known system Fe-Ni-B in the concentration regime where magnetism is not supposed to be very strong (equal concentration of magnetic 3d atoms and nonmagnetic B atoms). Further we proceed with a structural search for new stable compounds and then study the stability of predicted antiferromagnetic (AFM)

FIG. 1. Schematics of HRS from a viewpoint of mean field description. Regions I and IV represent nonmagnetic and saturated regimes; Region II corresponds to quantum critical point and a system in the area III is in HRS regime with possible large gradients of magnetization and strong spin fluctuations. M is the order parameter while ϕ is the external perturbation (pressure, magnetic field, etc).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01025v2-figure-3ee45e1c-e2f9-4fda-9dad-6868d040d0eb.jpg)

order by adding possible SF in the ground state (quantum zero-point motion). At the end we estimate the influence of quantum SF on possible appearance of SF induced superconductivity.

Structure prediction -Crystal structures of Fe-Ni-B were searched by an adaptive genetic algorithm (AGA), which integrates auxiliary interatomic potentials and first-principles calculations together in an adaptive manner to ensure the high efficiency and accuracy (see Supplementary Materials for more details) [3, 4]. The structure searches were only constrained by the chemical composition, without any assumption on the Bravais lattice type, symmetry, atom basis or unit cell dimensions. A range of different composition around the FeNiB (111) compositional ratio (i.e., 211, 121, 112, 311, 131, 113, 221, 122, 212, 133, 313, 331, 321, 312, 123, 132, 231, 213, 223, 232, 322, 233, 323 and 332) were selected with 2 and 4 formula units to perform the AGA search. First-

FIG. 2. (a) The configuration of FeNiB 2 from AGA search. Atom color: Fe atoms colored with brown; Ni atoms colored with silver; B atoms colored with green. (b) Phonon density of States.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01025v2-figure-b1653ecf-34ba-4dc6-95af-2b25ede03f16.jpg)

principles calculations were carried out using the projector augmented wave (PAW) method [5] within density functional theory (DFT) as implemented in VASP code [6, 7]. The exchange and correlation energy is treated within the spin-polarized generalized gradient approximation (GGA) and parameterized by PerdewBurke-Ernzerhof formula (PBE) [8]. A plane-wave basis was used with a kinetic energy cutoff of 520 eV. During the AGA search, Monkhorst-Pack's sampling scheme [9] was adopted for Brillouin zone sampling with a kpoint grid of 2 π × 0 033 ˚ . A -1 , and the ionic relaxations stopped when the forces on every atom became smaller than 0 01 eV . / ˚. A The energy convergence criterion is 10 -4 eV. The phonon calculations were performed with the finite difference method via Phonopy code [10]. The magnetic properties were calculated in a higher-quality calculation using a finer k-point sampling grid of spacing 2 π × 0 02 ˚ . A -1 and a stricter energy convergence criterion of 10 -6 eV.

Among all AGA searched ternary phases, a new FeNiB 2 phase, with 2 formula units and a space group of P 2 1 /m (space group number 11), is found to show the energy even smaller than the Gibbs triangle formed by the existing stable phases on the ternary phase diagram. The structure of this new phase is shown in Fig. 2(a). There is one Wyckoff site for Fe, Ni and two Wyckoff sites for B, respectively (the crystallographic details are given in Supplementary Materials). The Fe, Ni and B atoms formed separated 1D chain. Phonon calculation further confirmed the phase is dynamically stable without any soft phonon mode at T = 0 K, which is shown in Fig. 2(b).

With the newly found FeNiB 2 structure, the existing Fe-Ni-B phase diagram should be updated. In the current DFT-calculated database (i.e., the Materials Project [11]), there are eight reported binary stable phases however there is no ternary stable phase reported. It should be noted that while the experimental database (i.e. Inorganic Crystal Structure Database [12]) shows that a few ternary phases such as FeNiB, Fe NiB, 2 Fe Ni B 3 3 2 and Fe Ni 3 20 B 6 previously synthesized by experiments, they

FIG. 3. DFT-calculated convex hull of the Fe-Ni-B system based on (a) previously reported phases and (b) after including the newly discovered FeNiB 2 phase. All DFT calculations were carried out at 0 K. Red solid triangles represent stable compounds and the green solid star represent newly found stable phase, while blue and orange open triangles indicate the metastable phases from AGA search and experiments, respectively. The black lines separate the compositional space to Gibbs triangles.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01025v2-figure-f6bacce6-41e0-4e1d-bc94-4958418e199d.jpg)

are all metastable phases according to DFT calculations, indicating they may decompose to the surrounding stable phases on the Gibbs triangle in the equilibrium state [13, 14]. The existing ternary convex hull computed from DFT calculation are shown in Fig. 3(a) (the detailed energy values are provided in Supplementary Materials). At the composition of FeNiB , the thermodynamic sta2 bility is determined by the Gibbs triangle formed by the nearby Ni B , FeB, and FeB 4 3 2 phases. Considering the following reaction FeNiB 2 → (Ni B 4 3 +3FeB + FeB ) 4, 2 / FeNiB 2 is a stable phase if the formation energy of the left hand side is lower than that of the right hand side. The DFT calculations indeed show that the formation energy of FeNiB 2 in the P 2 1 /m structure is equal to -340 9 meV/atom, while at the right-hand side, i.e. . the energy on the convex hull at the composition of FeNiB , is calculated to be 2 -326 57 meV/atom. . Therefore, FeNiB 2 is thermodynamically stable at the temperature of 0 K and the convex hull of the Fe-Ni-B system should be updated with the inclusion of the new FeNiB 2 phase. The new phase diagram of Fe-Ni-B system at 0 K are show in Fig. 3(b).

Considering the stability of FeNiB 2 phase, we search the literature to see if there is any hint of its existence. First, we find the stable CoB phase (tetragonal, space group Pnma ) shows a same motif as the current FeNiB 2 phase. From the chemical viewpoint, Co atoms should easily be replaced by Fe and Ni atoms. Moreover, while the ground state of FeB is the tetragonal phase ( I 4 1 /amd ), experiments always result the metastable FeB analogous to the CoB structure, due to the vibrational properties as suggested in Ref. [15]. Very interestingly, Gianoglio and Badini [14] reported that Ni can substitute Fe in the FeB-Pnma phase up to 70%. They showed the lattice of 50% Ni substituted FeB were 0.40,

FIG. 4. (a) The energy difference and (b) the Fe magnetic moment at AFM and FM states as a function of the volume for the FeNiB 2 structure. The insert in (a) shows the total energies for both configurations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01025v2-figure-24ba448b-b645-44de-9dca-2365bf6f319a.jpg)

0.30 and 0.54 nm, which is consistent with our current results of 0.395, 0.30 and 0.532 nm except that the cell angle β becomes 91 35 . ◦ in our DFT calculations. With these evidences, we believe the FeNiB 2 exists however it is not explicitly reported by experiments.

With the verified stability of the new FeNiB 2 structure we are now in position to investigate its magnetic ground state. The DFT calculations are performed for the ferromagnetic (FM) state and the simplest AFM state, in which two Fe atoms in the primitive cell show opposite spin directions. It shows a strong competition of energetic stability between the FM and AFM states in the FeNiB 2 structure at 0 K. As shown in Fig. 4(a), the two states exhibit very similar energies, while the AFM state has a slightly lower energy ∼ 0 45 meV/atom at the equi-. librium volume 7 8 ˚ . A /atom. 3 As the volume increases, the energy difference becomes larger and the AFM state is always slightly more stable than FM state. More interestingly, in Fig. 4(b), the equilibrium magnetic moments demonstrate a rapid increase (nearly by 10 times). While the current results were based on a fixed geometry of the lattice, the calculation with relaxed structures shows almost same results (see Supplementary Material). By using the mean field approximation of Heisenberg model, the Neel temperature is estimated to be nearly 50-60K (the Heisenberg model parameters were calculated using Ref. [16]). It suggests the FeNiB 2 structure is locating in the highly responsive state (region III in Fig. 1). Coexistence of different magnetic long-range orders even for small values of atomic magnetic moments indicates that energy profile should have many local minima without well-defined global minimum. Correspondingly, the system is very sensitive to any external perturbation and could easily change its magnetic state (magnetic chameleon) including magnetic tunneling. We note that this HRS can be a characterization of FeNiB 2 phase compared to the analogous FeB-Pnma phase which has a ferromagnetic ground state.

The effects of zero-point spin fluctuations on the magnetic stability of the ground state of FeNiB 2 -Let us first

TABLE I. The partial density of state at the Fermi level and the Stoner parameter at each atom for FeNiB . 2

|               |   DOS @ E f (State/eV) | I (eV)   |
|---------------|------------------------|----------|
| Total ( N t ) |                 1.962  |          |
| Fe            |                 1.477  | 0.84     |
| Ni            |                 0.312  | 0.87     |
| B             |                 0.0895 | 2.20     |

address the issue of magnetic instability in FeNiB 2 . Our nonmagnetic calculations of this system revealed that the density of states (DOS) at the Fermi level is dominated by Fe electrons. So the Stoner criterion is fulfilled

$$1 - I \chi _ { 0 } > 0.$$

Here χ 0 is non-enhanced magnetic susceptibility, while I is the effective Stoner parameter of static LDA. For multicomponent system the corresponding Stoner criterion was obtained in Ref. [18] and is written as

$$N _ { t } \left ( E _ { F } \right ) I > 1 \text{ or } \sum _ { i } I _ { i } N _ { i } ^ { 2 } > N _ { t }, \text{ \quad \ } \text{(2)}$$

where the effective Stoner parameter I for multicomponent system is defined as I = ∑ i I i ( N i N t ) 2 , where I i and N i are partial Stoner parameter and partial DOS at the Fermi level, correspondingly, and N t is the total DOS at the Fermi level. The obtained numbers are shown in Table I and the resulting criteria of possible long-range AFM order is fulfilled. Such criteria are also fulfilled at many other q -vectors, including FM one. We also found that criteria of a local magnetic moment stability (Anderson criteria [19, 21]) is satisfied as well.

However, LDA traditionally overestimates the magnetism for weakly magnetic systems due to neglect of quantum SF at T = 0 K (spin zero-point motion). To study the stability of predicted above weak magnetism in FeNiB 2 we used a modified version of SF theory from Ref.[20]. First, we obtain the ab initio SF energy (absent in LDA) using RPA technique [21]. Then we calculate a corresponding SF contribution to the spin susceptibility and finally obtain renormalized Stoner criterion for magnetic instability.

The following RPA expression for the SF contribution to the total energy has been used

$$E _ { S F } = \lim _ { T \rightarrow 0 } F _ { S F } = 2 \hbar { \lim } _ { T \rightarrow 0 } \sum _ { v } \int _ { 0 } ^ { I } d I \int _ { \Omega _ { B Z } } d \mathfrak { q } \int _ { - \infty } ^ { \infty } d \omega \left ( \text{Im} \chi _ { v } \left ( \mathfrak { q }, \omega \right ) - \text{Im} \chi _ { v 0 } \left ( \mathfrak { q }, \omega \right ) \right ) \coth \left ( \beta \omega / 2 \right ).$$

This expression can also be rewritten as

$$E _ { S F } = \frac { 1 } { 2 } \int _ { 0 } ^ { I } d I \left ( \left \langle M ^ { 2 } \right \rangle - \left \langle M _ { 0 } ^ { 2 } \right \rangle \right ), \quad \ \ ( 4 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where 〈 M 2 〉 is a square mean of magnetic moment obtained with enhanced and nonenhanced susceptibility, correspondingly (see Fig.8 in Ref.[22] for detailed definition and calculations of these quantities in pure 3d metals). Computational details of our linear response method can be found in Ref.[22].

The transverse part of SF correction to the inverse magnetic susceptibility is obtained directly as

Our results indicated that the SF renormalization of the effective Stoner parameter in FeNiB 2 is nearly 7% with dominating suppression coming from Fe sites. Many different SF with q-vectors around AFM instability contribute to this suppression in the frequency range up to 0 3 . eV, but overall transitions inside electronic d-band (4 -5 eV) contribute to this renormalization. It suggests that the existence of predicted SF can be measured by neutron scattering experiments. Both longitudinal and transversal SF on Ni sites are also strong but due to small original static moment on Ni atoms, the Ni contribution to the total magnetization suppression is small.

$$\chi _ { S F } ^ { - 1 } = \frac { 1 } { 2 M } \frac { \partial } { \partial M } \left ( \int _ { 0 } ^ { I } d I \left ( \left \langle M ^ { 2 } \right \rangle - \left \langle M _ { 0 } ^ { 2 } \right \rangle \right ) \right ). \quad ( 5 ) \quad \text{retic} \\ \text{not } \text{i}$$

This expression is different from the one used in Refs.[20, 23] as we do not employ low frequency and long wavelength approximations [24]. The SF renormalized Stoner criterion then is written as

$$\left ( \chi _ { S F } ^ { - 1 } + \chi ^ { - 1 } \right ) \chi _ { 0 } = 1 - I ^ { * } \chi _ { 0 } > 0, \quad \ \left ( 6 \right ) \quad t o \ \Sigma$$

The appearance of AFM like magnetic instability can also be favorable for SF induced superconductivity. While s-wave and triplet superconductivity could theoretically appear near FM instability, such systems are not really known. Superconductivity reports in weakly FM system ZnZr 2 under pressure have not been confirmed. But more popular singlet superconductivity can naturally appear near AFM instability (cuprates [1], iron pnictides [2]). So, the appearance of superconductivity looks possible in this system and now we will estimate a strength of electron-spin interaction and its contribution to SF induced superconductivity.

where I ∗ is a renormalized Stoner parameter, which includes quantum spin zero points SF beyond usual static local density approximation. This renormalization was first noticed by Moriya [20] and studied recently using electronic structure calculations in Ref. [23] for several simple 3d metals. All these studies have been performed using ad-hoc low frequency and long wave approximations. Our approach is truly ab-initio and takes into account SF at all q vectors and frequencies.

The effects of zero-point spin fluctuations on superconductivity in FeNiB 2 -From practical point of view, e.g., for the search of new superconducting material, the qualitative estimation of λ is of primary interest. To obtain the required strength of spin-electron interaction we will use well known s-d exchange model [25]. To estimate λ we do not need to consider the superconducting phase since it enters just the mass renormalization in the normal state and one needs just to deal with the corresponding expression for the electron self-energy [23],

$$\Sigma _ { k } \left ( E \right ) = I ^ { 2 } \sum _ { q } \int _ { 0 } ^ { \infty } \frac { d \omega } { 4 \pi } I m \chi \left ( { \mathbf k } - { \mathbf q }, \omega \right ) \left [ \frac { 1 - f _ { q } } { E - \varepsilon _ { q } - \omega } + \frac { f _ { q } } { E - \varepsilon _ { q } + \omega } \right ],$$

where f q is the Fermi distribution with ε q being one electron spectrum. The expression for the coupling constant at T = 0 K is then written in standard way as

$$\lambda = - \, \frac { \partial \left \langle \Sigma _ { k } \left ( E \right ) \right \rangle } { \partial E } \Big | _ { E = 0 } \,. \quad \quad \left ( 8 \right ) \, \frac { \text{pr} } { \text{fec} } \, \\ \text{we}$$

To make the simplest possible estimation of electron magnon interaction one can replace the average over the

Fermi surface above by the average over the Brillouin zone. To further obtain the temperature of the superconducting transition using the approach advocated above we use usual MacMillan formula [26]. Under these approximations our calculated value of λ = 0 6 and the ef-. fective SF frequency is ω = 0 2 eV. With these numbers . we obtain a value of 20 -30 K for the critical temperature. Thus, our qualitative estimation predicts not only suppression of magnetic long-range order in the normal

state but also suggest a very promising scenario for the SF induced superconductivity near AFM quantum critical point.

In summary, we performed a structure search for a stable system near magnetic instability. Our AGA approach which integrates auxiliary interatomic potentials and first-principles calculations found a stable compound FeNiB , which appears to be a magnetically highly re2 sponsive system close to antiferromagnetic instability. The structure of FeNiB 2 is analogous to the FeB-Pnma phase and previous literature provide hints of its possible existence. Compared to the ferromagnetic long-range order of FeB, this system is characterized by strong antiferromagnetic-like spin fluctuations at T = 0K (quantum spin zero-point motion) and possible spin fluctuation induced superconductivity. The experimental verifications of our prediction are highly desirable. Our magnetostructural combined studies thus open a promising route for the search of strong SF state and superconductors.

Acknowledgements Work at Ames Laboratory was supported by the U.S. Department of Energy, Basic Energy Sciences, Materials Science and Engineering Division, under Contract No. DEAC02-07CH11358, including a grant of computer time at the National Energy Research Scientific Computing Center (NERSC) in Berkeley, CA. Work at University of Science and Technology of China was supported by the National Natural Science Foundation of China (11574284 &amp; 11774324) and the Supercomputing Center of the University of Science and Technology of China.

∗ Email: yangsun@ameslab.gov

- † Email: antropov@ameslab.gov
- [1] T. Moriya and K. Ueda, Adv. Phys. 49, 555 (2000).
- [2] D. C. Johnston, Adv. Phys. 59, 803 (2010).
- [3] S. Q. Wu, M. Ji, C. Z. Wang, M. C. Nguyen, X. Zhao, K. Umemoto, R. M. Wentzcovitch, and K. M. Ho, J. Phys.
- Condens. Matter 26, 035402 (2014).
- [4] X. Zhao, M. C. Nguyen, W. Y. Zhang, C. Z. Wang, M. J. Kramer, D. J. Sellmyer, X. Z. Li, F. Zhang, L. Q. Ke, V. P. Antropov, and K. M. Ho, Phys. Rev. Lett. 112, 045502 (2014).
- [5] P. E. Blchl, Phys. Rev. B 50, 17953 (1994).
- [6] G. Kresse and J. Furthmller, Comput. Mater. Sci. 6, 15 (1996).
- [7] G. Kresse and J. Furthmller, Phys. Rev. B 54, 11169 (1996).
- [8] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996).
- [9] H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1976).
- [10] A. Togo and I. Tanaka, Scr. Mater. 108, 1 (2015).
- [11] A. Jain, S. P. Ong, G. Hautier, W. Chen, W. D. Richards, S. Dacek, S. Cholia, D. Gunter, D. Skinner, G. Ceder, and K. A. Persson, APL Mater. 1, 011002 (2013).
- [12] M. Hellenbrandt, Crystallogr. Rev. 10, 17 (2004).
- [13] Y. Khan, Z. Met. 74, 385 (1983).
- [14] C. Gianoglio and C. Badini, J. Mater. Sci. 21, 4331 (1986).
- [15] A.N. Kolmogorov, S. Shah, E.R. Margine, A.F. Bialon, T. Hammerschmidt, and R. Drautz, Phys. Rev. Lett. 105, 217003 (2010)
- [16] A.I. Liechtenstein, M.I. Katsnelson, V.P. Antropov, and V.A. Gubanov, J. Magn. Magn. Mater. 67, 65 (1987).
- [17] G.D. Samolyuk and V.P. Antropov, Phys. Rev. B 79, 052505 (2009).
- [18] T. Jarlborg and A. J. Freeman, Phys. Rev. B 22, 2332 (1980).
- [19] P.W. Anderson, Phys. Rev. 124, 41 (1961).
- [20] T. Moriya and A. Kawabata, J. Phys. Soc. Japan 34, 639 (1973).
- [21] F. Englert and R. Brout, Phys. Rev. 120, 1085 (1960).
- [22] A. L. Wysocki, V. N. Valmispild, A. Kutepov, S. Sharma, J. K. Dewhurst, E. K. U. Gross, A. I. Lichtenstein, and V. P. Antropov, Phys. Rev. B 96, 184418 (2017).
- [23] V. P. Antropov and A. Solontsov, J. Appl. Phys. 109, 07E116 (2011).
- [24] V. P. Antropov, J. Magn. Magn. Mater. 262, L192 (2003).
- [25] V. Y. Irkhin and M. I. Katsnelson, Z. Phys. B 75, 67 (1989).
- [26] W. L. McMillan, Phys. Rev. 167, 331 (1968).
Renhai Wang Yang Sun Vladimir Antropov Zijing Lin Cai-Zhuang Wang Kai-Ming Ho
Theoretical prediction of a highly responsive material: Spin fluctuations and superconductivity in FeNiB2 system
By analyzing Fe-Ni-B compositional diagram we predict an energetically and
dynamically stable FeNiB2 compound. This system belongs to the class of highly
responsive state of material, as it is very sensitive to the external
perturbations. This state is also characterized by a high level of spin
fluctuations which strongly influence possible magnetic long- and short-range
orders. Furthermore, we demonstrate that these antiferromagnetically dominating
fluctuations could lead to the appearance of spin mediated superconductivity.
The obtained results suggest a promising avenue for the search of strong spin
fluctuation systems and related superconductors.