Spin-lattice and electron-phonon coupling in 3$d$/5$d$ hybrid Sr$_3$NiIrO$_6$

TL;DR

This study investigates the complex spin, lattice, and electronic interactions in Sr$_3$NiIrO$_6$, revealing strong spin-lattice and electron-phonon couplings that influence its magnetic properties and excitations.

Contribution

It combines experimental spectroscopy and first-principles calculations to uncover novel spin-lattice and vibronic couplings in a 3d/5d hybrid material.

Findings

Identification of spin-lattice coupling of three phonons affecting Ir environment

Discovery of inter-chain motion influencing magnetic interactions

Revelation of vibronic coupling and large crystal field effects in Ir ions

Abstract

While 3$d$-containing materials display strong electron correlations, narrow band widths, and robust magnetism, 5$d$ systems are recognized for strong spin-orbit coupling, increased hybridization, and more diffuse orbitals. Combining these properties leads to novel behavior. Sr$_3$NiIrO$_6$, for example, displays complex magnetism and ultra-high coercive fields - up to an incredible 55~T. Here, we combine infrared and optical spectroscopies with high-field magnetization and first principles calculations to explore the fundamental excitations of the lattice and related coupling processes including spin-lattice and electron-phonon mechanisms. Magneto-infrared spectroscopy reveals spin-lattice coupling of three phonons that modulate the Ir environment to reduce the energy required to modify the spin arrangement. While these modes primarily affect exchange within the chains, analysis also uncovers important inter-chain motion. This provides a mechanism by which inter-chain interactions can occur in the developing model for ultra-high coercivity. At the same time, analysis of the on-site Ir$^{4+}$ excitations reveals vibronic coupling and extremely large crystal field parameters that lead to a t$_{2g}$-derived low-spin state for Ir. These findings highlight the spin-charge-lattice entanglement in Sr$_3$NiIrO$_6$ and suggest that similar interactions may take place in other 3$d$/5$d$ hybrids.

Tables

Table 1. Table 1: Summary of 10 D ​ q 𝐷 𝑞 Dq and Racah parameters for the 3 d 𝑑 d /5 d 𝑑 d materials in this study compared to the crystal field parameters of other Ir-containing compounds along with information about the vibronically coupled phonons. No vibronic coupling analysis was found for Sr 2 IrO 4 or Li 2 IrO 3 , and the Ni 2+ excitations in Ni 3 TeO 6 are not vibronically activated, so no phonons are implicated.

Material	Element	Electronic	10$Dq$ (eV)	$B$ (eV)	Coupled phonon	Displacement	Refs.
		state			frequency (cm-1)
Sr3NiIrO6	Ir	5$d^5$	3.24	1.18	177	Ni in-plane motion	This work
Sr3CuIrO6	Ir	5$d^5$	2.33	0.86	273	O-Ir-O bend	This work
Sr2IrO4	Ir	5$d^5$	3.8	0.93	-	-	51; 52
Li2IrO3	Ir	5$d^5$	2.7	0.95	-	-	53
Ni3TeO6	Ni	3$d^8$	1.10	0.11	-	-	55
$\alpha$-Fe2O3	Fe	3$d^5$	1.59	0.09	525	In-plane, in-phase Fe-O stretch	54; ONeal2014a

Full text

Introduction

Interest in 5$d$ materials and 3$d$/5$d$ hybrids has blossomed in response to recent scientific advances and applications in hard magnets, topological insulators, multiferroics, superconductors, and thermoelectrics Rau2016 ; Witczak-Krempa2014 ; Feng2017 ; Feng2018 ; Liu2019 ; Isobe2019 . 5$d$ materials are unique for several reasons. First, strong spin-orbit coupling competes with magnetic, crystal-field, many-body Coulomb, and other interactions to drive new physical behaviors Cao2018 , such as the $J_{\rm eff}=1/2$ state in certain iridates Kim2009 ; Clancy2018 . Second, the bonding interactions associated with the larger 5$d$ orbitals promote inter-cation dimerization in pairwise, chain-like, and other complex orderings Yang2012a ; Hsu2013 . Third, the relativistic shifts in orbital energies, combined with spin-orbit and bandwidth effects, can drive band inversions leading to topological phases and enhanced Rashba splittings Zheng2015 ; Goldman2011 ; Bihlmayer2007 ; Sie2019a . In contrast, 3$d$ transition metal compounds typically display much narrower bandwidths, more robust magnetism, and stronger electron-electron interactions, and correlations Charnukha2017 ; Liu2018 . When these two sets of properties are brought together, as in Sr3NiIrO6, new and potentially useful behaviors can emerge.

What makes Sr3NiIrO6 so remarkable is the extraordinary coercivity - up to 55 T depending on sample details Singleton2016 . By contrast, traditional hard magnets like Fe/Pt, Nd31-xDyxFebalCo2B1 ($x$ = 7 wt %), and LuFe2O4 have coercivities on the order of 1, 3, and 9 T, respectively Ristau1999 ; Bai2007 ; Wu2008 . The extraordinarily high coercive field is not due to ferromagnetic domains since the material is antiferromagnetic, though the exact mechanism remains to be understood Singleton2016 ; Ou2014 ; Zhang2010 ; Birol2018 .

The crystal structure consists of chains of alternating face-sharing NiO6 trigonal prisms and IrO6 octahedra stacked along the $c$ axis with Sr atoms separating the chains. These chains are arranged in a triangular configuration in the $ab$ plane Mikhailova2012 . Because Ir has a 4+ charge, it was originally thought to have a $J_{\rm eff}=1/2$ state Lefrancois2014 ; Toth2016 ; Zhang2010 . A trigonal distortion, however, breaks the symmetry and creates a slight deviation from the pure $J_{\rm eff}=1/2$ Lefrancois2016 . Sr3NiIrO6 is part of a chemical family of quasi-one-dimensional 3$d$/5$d$ hybrids with formula $A_{3}BB^{\prime}$O6. Notably, these materials display increasing coercive magnetic fields as the $B^{\prime}$ site ion evolves from a 3$d$ to a 4$d$ to a 5$d$ magnetic ion Mikhailova2012 ; Stitzer2002 ; Singleton2016 . This suggests that the unusual properties of Ir and similar 5$d$ ions must play an important role in creating the high coercivity. Recent calculations indeed suggest that it is the strongly anisotropic exchange interaction characteristic of the Ir4+ ion that plays a key role the high coercive magnetic field Birol2018 . The strongest magnetic exchange is antiferomagnetic between Ni and Ir ions within $c$-axis chains, which creates one-dimensional order below $\approx$200 K Toth2016 . Elastic neutron diffraction measurements identified two potential magnetic ground states at zero field: a partially disordered antiferromagnet or a more complex, spatially modulated state Lefrancois2014 . In light of the large Ni$\cdots$Ir intrachain exchange interactions Toth2016 , the partially disordered antiferromagnet model is most likely and has been suggested for several $A_{3}BB^{\prime}$O6 materials with high coercive field Niitaka2001 ; Mohapatra2007 ; Murthy2019 . In this model, $c$-axis chains act like giant spin units due to the difference between the Ni and Ir moments, with frustrated antiferromagnetic interactions between chains in the $ab$-plane. Consequently, the ground state consists of two oppositely aligned chains and one with random orientation. In this model, gradual flipping of the randomly oriented third chains accounts for the spin glass-like dynamics with magnetic hysteresis and coercivity that is seen below 15 K in Sr3NiIrO6, and for the initial evolution of the magnetization up to the coercive field. At the coercive field, the second chain flips, leading to a sudden jump and the high coercive field. The extraodrinary coercivity must result from the barrier to flipping of the second chain due to its many-body nature and interactions with the lattice including the anisotropic exchange interaction of Ir. Thus if we are to achieve a complete understanding of the extraordinary coercive magnetic field in Sr3NiIrO6 we need to quantify how the magnetic order and Ir electronic state interact with the lattice.

To this end, we measured the infrared and optical properties of Sr3NiIrO6 and compared the response to high field magnetization, first principles lattice dynamics calculations, and simulations of spin-phonon coupling. While the phonons are nearly rigid across the magnetic ordering and freezing transitions, magneto-infrared work reveals three modes that track the magnetic response. What distinguishes these modes is how they modulate superexchange pathways around the Ir centers. The O-Ir-O bending mode near 310 cm*-1* is the dominant feature with an exceptionally large spin-lattice coupling constant - on the order of 10 cm*-1*. That very specific local lattice distortions are involved in the approach to the coercive field provides further confirmation that the high field transition in Sr3NiIrO6 is more than a simple domain reorientation process Singleton2016 . At the same time, we reached beyond spin-lattice interactions to explore electron-phonon coupling. Analysis of the Ir4+ on-site $d$-to-$d$ excitations reveals that they are vibronically activated. An oscillator strength analysis uncovers electron-phonon coupling involving the 177 cm*-1* Ni in-plane mode which modulates the Ni-O-Ir bond angles (and thus the local Ir4+ environments). We discuss these findings in terms of inter- and intra-chain coupling as well as how the diffuse character of the Ir 5$d$ orbitals enlarges the crystal field parameters to give a spin configuration that derives from the $t_{2g}$-derived energy levels.

Results and discussion

Strong spin-lattice coupling and intra-chain interactions in Sr3NiIrO6

Figure 1 (a) displays the infrared spectrum of Sr3NiIrO6. A symmetry analysis of the $R\bar{3}c$ space group yields seven $A_{2u}$ and twelve $E_{u}$ infrared-active phonons - consistent with our spectra. Our mode assignments are based upon lattice dynamics calculations which also provide symmetries and displacement patterns. These assignments are summarized in Table S1 (Supplementary Information) and allow the microscopic interactions indicated by the spectra to be fully understood. Spin-lattice coupling across the 75 K magnetic ordering and 15 K freezing temperatures is extremely weak (Supplementary Information), likely because these energy scales do not correspond to collective transitions Mikhailova2012 ; Lefrancois2014 ; Singleton2016 . The signature in other thermodynamic probes is not strong either. Magnetoelastic coupling in applied field, however, is robust.

In order to search for elastic contributions to the record high coercivity, we measured the magneto-infrared response of Sr3NiIrO6 [Fig. 1]. Here, magnetic field was increased step-wise from 0 to 35 T, and spectra were acquired at each step to avoid complications from magnetic hysteresis. The measurements therefore follow the line from the origin in Fig. 1 (b), shown more closely in Fig. 1 (c). The absorption difference ($\Delta\alpha=\alpha(30$ T$)-\alpha(0$ T)) in the upper portion of Fig. 1 (a) highlights spectral differences. Three phonons are sensitive to magnetic field. The Eu symmetry mode at 310 cm*-1* displays the largest magneto-infrared response [Fig. 1 (f)]. The A2u symmetry modes at 133 and 534 cm*-1* are also sensitive to changes in the microscopic spin arrangement [Fig. 1 (d,h)]. Interestingly, the two lower frequency phonons harden with applied field, whereas the higher frequency mode softens.

We quantify the magnetic field dependence of the 133, 310, and 534 cm*-1* phonons by integrating the absolute value of the absorption difference over a narrow frequency window ($\int_{\omega_{1}}^{\omega_{2}}|\Delta\alpha|d\omega$) at each field. This quantity is proportional to the field-induced frequency shift, although error bars on ${\Delta}{\alpha}$ are much smaller Brinzari2012 . Comparison with the bulk magnetization [Fig. 1 (c)] reveals that these changes grow as magnetization squared Granado1999 . This implies that the field-driven transition is not just a spin reorientation process; it also involves specific local lattice distortions. The magnetic field drives the spin rearrangement, while cooperative lattice distortions reduce the exchange interactions and the energy required to modify the microscopic spin arrangement. Similar mechanisms are active in other materials Granado1999 ; Poirier2007 ; Poirier2011 . While the 48 T spin-reorientation field is beyond the current reach of resistive magnets, the observed spin-phonon coupling is quite strong. We can estimate the spin-phonon coupling constants in Sr3NiIrO6 as $\omega=\omega_{0}+\lambda<S_{i}\cdot S_{j}>$. Here, $\omega_{0}$ is the unperturbed (zero field) mode frequency, $\omega$ is the perturbed (high field) mode frequency, and $<S_{i}\cdot S_{j}>$ is the spin-spin correlation function Casto2015 . Taking the limiting low temperature value of the spin-spin correlation function as $<S_{i}\cdot S_{j}>=<{1\cdot\frac{1}{2}}>=\frac{1}{2}$, we find $\lambda\approx 2$, 10, and 5 cm*-1* for the 133, 310, and 534 cm*-1* modes, respectively. By comparison, the highest value of $\lambda$ in the 3$d$/4$p$ hybrid Ni3TeO6 at the same applied magnetic field (i.e. across the spin-flop transition) is 4 cm*-1* Yokosuk2015 . This suggests that higher fields may drive even stronger spin-lattice interactions, especially across the coercive field Kim2015 ; Jaime2012 .

Examination of the calculated displacement patterns uncovers important similarities in the spin-lattice coupled modes that provide mechanistic insight into the field-driven process. The Eu symmetry O-Ir-O bending mode at 310 cm*-1* is the most sensitive to field [Fig. 1 (g)]. This displacement pattern - with its strong O component - unambiguously modulates the superexchange interactions between Ir and its neighbors by changing the Ir environment. This motion also reduces the symmetry of the crystal and (as discussed below) introduces new terms into the magnetic Hamiltonian. Local lattice distortions of this type are thus important to the development of various spin rearrangements - such as that across the coercive field. The other two magneto-infrared-active features modulate the exchange interactions around the Ir centers as well - although by less effective routes. For instance, the A2u symmetry mode at 133 cm*-1* consists of Sr out-of-plane displacement and in-phase Ni displacement along the $c$ axis [Fig. 1 (e)] that slightly varies the Ni-O-Ir angle as a second-order effect. The A2u symmetry mode at 534 cm*-1*, on the other hand, consists mainly of O stretching around the Ir center [Fig. 1 (i)] which also impacts the superexchange angles. The aforementioned displacement patterns primarily affect exchange within the chains. A view along the chains (see images and animations of the displacement patterns in Supplementary Information) reveals that these modes have inter-chain motion as well. This supports the role of inter-chain interactions in the developing magnetic model Lefrancois2014 ; Singleton2016 and provides a mechanism by which such effects can occur.

The form of the spin-phonon Hamiltonian for Sr3NiIrO6 is very complicated - even when considering only intra-chain interactions - because various phonons break lattice symmetries and allow novel magnetic interaction terms (such as Dzyaloshinskii-Moriya) to emerge. (See the supplementary information for details.) This makes the first principles calculation of every spin-phonon coupling parameter practically impossible. We therefore developed a simplified approach to predict which phonon modes have a stronger effect on the magnetic structure (and vice versa). We began by considering the displacement pattern of the 310 cm*-1* mode and calculated the effect of this type of modulation on the ground state spin arrangement. Similar calculations using the 177 cm*-1* mode pattern are included for comparison. As a reminder, the Eu symmetry O-Ir-O bending mode at 310 cm*-1* engages in spin-lattice coupling, whereas (as we shall see below) the Eu symmetry mode at 177 cm*-1* is the phonon that vibronically activates on-site $d$-to-$d$ excitations of Ir4+. When the atoms are not displaced, the Ni and Ir moments are predicted to be collinear and parallel to the $c$ axis (within error bars), consistent with the observed magnetism. Displacing atoms according to the pattern of the 310 cm*-1* mode leads to very significant tilting of the Ni spin moment as well as weaker tilting of the Ir spin [Fig. 1 (j)]. As a result, the system becomes non-collinear. The 177 cm*-1* mode, on the other hand, leads to a much smaller effect on the magnetic structure, in line with a lack of field-induced changes of this mode. This reveals the complex nature of the spin-lattice coupling and suggests that the induced tilting may lower the magnetic switching barrier.

Electron-phonon coupling and strength of the crystal field interactions in Sr3NiIrO6

Equipped with a microscopic understanding of the elastic distortions in response to the magnetic field, we sought to determine whether similar phonons contribute to the electronic properties - specifically the vibronically-activated crystal field excitations of Ir4+. The Ir4+ electronic configuration is key to creating the anisotropic exchange interactions that have been identified as an important contributor to the high coercive field, suggesting that additional insight would prove useful. We therefore measured the optical properties of Sr3NiIrO6 as well as the Cu analog, Sr3CuIrO6, to track the Ir4+ on-site excitations near 0.7 eV [Fig. 2 (a) and Supplementary Information]. While the position and general behavior of these excitations agree with prior resonant inelastic x-ray scattering work Lefrancois2016 ; Liu2012a , our measurements offer significantly higher spectral resolution. The presence of intra-$t_{2g}$ on-site excitations indicates that the Ir4+ site symmetry is slightly distorted from octahedral, creating a distortion away from a pure $J_{\rm eff}=1/2$ state Lefrancois2016 . Moreover, analysis reveals that these excitations are vibronically activated by a phonon that is distinct from those that contribute to the coercivity via spin-lattice pathways.

Inter-band excitations, like those in Sr3NiIrO6, are responsible for the colors of transition metal-containing materials and have been extensively studied Sell1967 ; Lohr1972 ; Ballhausen1962 . Vibronic coupling, in which an odd-parity phonon mixes with a $d$-to-$d$ excitation, is a common activation mechanism Sell1967 ; Lohr1972 ; Ballhausen1962 . Here, both spin and parity selection rules are broken due to coupling with a phonon. In this scenario, the temperature dependence of the oscillator strength is modeled as $f=f_{x}+f_{0}\coth(h\nu/2k_{\rm B}T)$, where $\nu$ is the frequency of the activating phonon, $f_{0}$ is the oscillator strength at base temperature, $f_{x}$ represents oscillator strength from other mixing processes, and $h$, $k_{\rm B}$, and $T$ have their usual meanings ONeal2017 . This model can be used to determine which phonon activates the electronic transition. While phonon-assisted $d$-to-$d$ excitations have been reported in Sr3Ir2O7 Park2014 , quantitative vibronic coupling analyses such as we do here are rare in 4- and 5$d$- systems and particularly so in 3$d$/5$d$ hybrids.

Figure 2(b) displays the oscillator strength analysis of the Ir4+ on-site excitations in Sr3NiIrO6. There is a small deviation from the overall trend near 70 K that may be due to the spin ordering transition Flahaut2003 , but the size of our error bars precludes a detailed analysis. Examination reveals that the Ir-related $d$-to-$d$ excitations are vibronically coupled with the 177 cm*-1* Eu phonon - which is present in the infrared spectrum [Fig. 2 (c)]. According to our calculations, this mode consists of Ni in-plane motion against Sr counter-motion that indirectly modulates the Ir environment [Fig. 2 (d)]. Thus, the vibronically coupled mode is separate and distinct from those involved in magnetoelastic coupling (133, 310, and 534 cm*-1*). This separation offers the possibility of selective property control. We also carried out the same analysis for the more distorted Cu analog and unveiled coupling to a different displacement in which O motion more directly affects the Ir environment (Supplementary Information). This contrast likely emanates from the dissimilar chain configurations (linear for the Ni system and zigzag for the Cu analog), highlighting the importance of local symmetry in coupling processes.

The observation of Ir4+ excitations offers an opportunity to compare crystal field parameters of a 5$d$ center to the more commonly studied 3$d$-containing oxides. Based on the position and shape of the Ir-related intra-band excitations of Sr3NiIrO6 Lefrancois2016 and the $d^{5}$ Tanabe-Sugano diagram Tanabe1954 ; Tanabe1954b ; Tanabe1956 , we estimate 10$D$q = 3.24 eV and the Racah parameter $B$ = 1.18 eV. As a reminder, 10$Dq$ and $B$ describe the strength of crystal field interactions, and because 5$d$ orbitals are highly diffuse, the crystal field parameters are large. We find 10$Dq$ = 2.33 eV and $B$ = 0.86 eV for Sr3CuIrO6 Liu2012a . As summarized in Table 1, these values are similar to other Ir4+-containing materials including Sr2IrO4 and Li2IrO3 Sala2014 ; Ishii2011 ; Gretarsson2013 but much higher than prototypical transition metal oxides like $\alpha$-Fe2O3 Marusak1980 and even the 3$d$/4$p$ hybrid Ni3TeO6 Yokosuk2016 . We attribute this difference to the heavy mass of the Ir center, which is predicted to increase spin-orbit coupling and Racah parameters in free ions Ma2014 . Importantly, the large 10$Dq$ values separate the $t_{2g}$ and $e_{g}$ levels such that the $e_{g}$-derived bands in Sr3NiIrO6 play no role in determining the Ir spin configuration.

In summary, we combined infrared and optical spectroscopy with high field magnetization and first principles calculations to explore coupling processes involving the fundamental excitations of the lattice in Sr3NiIrO6 - a material with significant spin-orbit interactions. These include both spin-lattice and electron-phonon processes. Magneto-infrared spectroscopy reveals that three phonons - all of which modulate the magnetic pathways around and the symmetry of the Ir centers - display strong spin-lattice interactions, demonstrating that the approach to the coercive field takes place with very specific local lattice distortions - different from expectations for simple domain reorientation in a ferromagnet. Examination of the mode displacement patterns also provides a specific mechanism for inter-chain interactions, a finding that is crucial to the development of the working magnetic model in Sr3NiIrO6 and related materials. At the same time, analysis of the on-site Ir4+ excitations unveils vibronic coupling and extremely large crystal field parameters. For instance, 10$Dq$ is a factor of two larger than that in traditional transition metal oxides, and the Racah parameter $B$ is a factor of 10 higher. Moreover, the phonon that activates the vibronic coupling has a completely different displacement pattern than those that are sensitive to magnetic field. We therefore find that certain phonons in Sr3NiIrO6 are strongly entangled with the spin and charge channels. In addition, the diffuse character of the Ir 5$d$ orbitals determines the ground state spin structure of Sr3NiIrO6 whereas spin-lattice interactions reduce the energy required to modify microscopic spin arrangements. Returning to the developing magnetic model of Sr3NiIrO6, we qualitatively expect that the ultra-high coercive field results from the high barrier to flipping a chain of strongly-coupled Ni2+ and Ir4+ magnetic moments, which is effectively a many-body magnetic system. This barrier to flipping results from the interaction of the magnetic order with the lattice. This interaction has two parts: the interaction of the magnetic exchange interactions with lattice distortions, and the mixing of the local spin-orbit coupled state of the Ir4+ ion with the lattice. We quantified both in this work, thereby providing necessary insight to develop a complete model.

Methods

High quality single crystals were grown as described previously Singleton2016 ; Nguyen1995 and either polished or combined with a transparent matrix to control optical density due to strong phonon absorption. Absorption was calculated as $\alpha=\frac{-1}{hd}\ln(\cal{T}(\omega))$, where $h$ is concentration, $d$ is thickness, and $\cal{T}(\omega)$ is the measured transmittance. Magneto-infrared and magnetization measurements were carried out at the National High Magnetic Field Laboratory using the 35 T resistive and 65 T short-pulse magnets, respectively. Absorption was obtained at zero field, and magneto-infrared measurements tracked changes. Absorption differences are calculated as $\Delta\alpha=\alpha($B$)-\alpha(0$ T). We integrate $|\Delta\alpha|$ over small energy windows (127-150, 290-335, and 530-545 cm*-1*) to quantify changes, and renormalize to the 35 T magnetization squared to match energy scales and provide a proper comparison.

First principles calculations were performed using projector augmented waves as implemented in VASP Kresse1996 ; Kresse1996a . The PBEsol functional was used to approximate the exchange correlation energy Perdew1996 ; Perdew2008 , and DFT+$U$ Liechtenstein1995 was utilized for the transition metals. For calculations of Sr3NiIrO6, $U$=5 eV for Ni and $U$=1 eV for Ir were used. For the Cu analog, $U$=5 eV for Cu and $U$=1.5 for Ir were used. Changing the values of $U$ and $J$ changes the obtained phonon frequencies but does not affect the physical picture. Phonons were calculated using both the direct method and density functional perturbation theory, and the results were identical within numerical noise. An $8\times 8\times 8$ Monkhorst-Pack grid Monkhorst1976 and a plane wave energy cutoff of 500 eV gave good convergence. The crystal structure used in the phonon calculations was obtained by relaxing the experimentally reported structure in the magnetically ordered phase.

The complexity of the crystal and magnetic structure of Sr3NiIrO6 makes it impractical to perform a detailed first-principles calculation of the spin-phonon or spin-lattice coupling in this material. In order to gain insight about the spin-lattice coupling strength, we therefore performed noncollinear magnetic calculations in a series of crystal structures that are obtained by displacing the atoms according to the displacement pattern of particular phonon modes. (Results are presented in Fig. 1 (j).) We considered displacements up to 0.5 Å in total for the 22 atom unit cell, which corresponds to an average of $\sim 0.1$ Å per atom. The magnetic ground state does not change, and is still predominantly ferrimagnetic along the $c$ axis for all of these structures. However, unlike the collinear ferrimagnetic state observed in the experiment and reproduced by first principles calculations, the magnetic ground state for these displaced structures are tilted ferrimagnetic. The 310 cm*-1* displacement gives rise to large tilting of Ni moments - as much as $\approx$10*∘* for the larger displacements we considered. For reference, we present also the results of the same calculation for the 177 cm*-1* mode of same symmetry, and find a much smaller tilting. This observation is not surprising, because the character of this mode does not change the exchange pathways significantly. (It is of 49% Sr, 13% Ni, 8% Ir, and only 29% O character.) The finding is also in line with the experimental result that this mode is not engaged in spin-lattice coupling.

Data Availability

Data is available from the corresponding author upon reasonable request.

Acknowledgments

Research at the University of Tennessee, Rutgers University, and University of Minnesota is supported by the National Science Foundation DMREF program (DMR-1629079, DMR-1629059, and DMR-1629260). The crystal growth was partially supported by the National Research Foundation of Korea (NRF) funded by the Ministry of Science and ICT (No. 2016K1A4A4A01922028). We also appreciate funding from the U.S. Department of Energy, Basic Energy Sciences, contract DE-FG02-01ER45885 (Tennessee), “Science at 100 Tesla” (LANL), and “Topological phases of quantum matter and decoherence” (LANL). The NHMFL facility is supported by the U.S. National Science Foundation through Cooperative Grant DMR-1644779, the State of Florida, and the U.S. Department of Energy.

Competing Interests

The authors declare no competing interests.

Author Contributions

J.L.M., T.B. and S.-W.C. devised the project. X.L. and S.-W.C. grew the samples. K.R.O., A.A., K.D.H., A.L.B., M.O., and J.L.M. carried out the spectroscopic measurements. V.Z., C.V.T., and J.S. measured the magnetic properties. A.P. and T.B. performed the theoretical calculations. K.R.O., J.L.M., V.Z., and T.B. wrote the manuscript, and all authors contributed to it.