Orbital-Resolved Spin Fluctuations

• Orbital-resolved spin fluctuations are the decomposition of spin fluctuation spectra by orbital channel, enabling precise analysis of magnetic, orbital, and charge dynamics in correlated systems.
• Techniques such as multi-orbital RPA, DFT+DMFT, and polarized RIXS quantify orbital contributions, revealing dominant channels in materials like iron-based superconductors and ruthenates.
• Understanding orbital-resolved spin dynamics facilitates insights into magnetic instabilities, orbital-selective Mott phases, and the emergence of unconventional superconductivity and nematic order.

Orbital-resolved spin fluctuations refer to the analysis and theoretical calculation of spin fluctuation spectra with explicit decomposition into contributions from different orbital channels. This approach is essential for understanding the intertwined spin, orbital, and charge dynamics in correlated electron systems, especially in multi-orbital materials such as iron-based superconductors, ruthenates, heavy fermion compounds, and transition metal oxides. Orbital-resolved analysis provides fundamental insights into magnetic instabilities, mechanisms of unconventional superconductivity, and the emergence of complex ordered states.

1. Formulation in Multi-Orbital Models

The starting point for studying orbital-resolved spin fluctuations is a multi-orbital Hamiltonian, typically of Hubbard-type, combining hopping, local Coulomb interactions, and, where relevant, additional interactions such as quadrupolar electron-phonon couplings. In a general $N$-orbital system, the Hamiltonian takes the form: $H = H_0 + H_U + H_{\mathrm{orbital}}$ where $H_0$ describes kinetic energy and inter-orbital dispersion (from tight-binding or ab initio methods), $H_U$ contains the rotationally invariant multi-orbital on-site Coulomb interactions (parametrized by $U,~U',~J,~J'$), and $H_{\mathrm{orbital}}$ includes orbital-quadrupole or other relevant couplings (Saito et al., 2013).

The local spin operator for orbital $m$ is expressed as

$$S_m = \sum_{\alpha,\beta} c_{m\alpha}^\dagger \frac{\vec{\sigma}_{\alpha\beta}}{2} c_{m\beta}$$

where $c_{m\alpha}^\dagger$ creates an electron of spin $\alpha$ in orbital $m$.

The bare (irreducible) susceptibility in orbital space reads

$$[\chi^0(q)]_{l_1 l_2, l_3 l_4} = -\frac{T}{N} \sum_k G_{l_1 l_3}(k+q)G_{l_4 l_2}(k)$$

so that all subsequent physical (RPA-level or beyond) spin and charge/orbital susceptibilities carry explicit orbital indices (Saito et al., 2013, Boehnke et al., 2018).

2. Random Phase Approximation and Orbital Decomposition

The RPA formalism captures the enhancement or suppression of spin fluctuations via interaction vertices,

$$\chi^s(q) = [1 - \Gamma^s \chi^0(q)]^{-1} \chi^0(q)$$

where the spin interaction vertex $\Gamma^s$ encodes intra- and inter-orbital Hund's coupling effects. For the orbital (charge/quadrupole) channel,

$$\chi^c(q) = [1 - \Gamma^c \chi^0(q)]^{-1} \chi^0(q)$$

with $\Gamma^c$ obtained from the Coulomb multiplet structure and, where relevant, electron-phonon-quadrupole contributions (Saito et al., 2013, Singh et al., 2018).

Orbital-resolved spin-fluctuation spectra are then represented by components $\chi^s_{l l, l l}(q)$, allowing assignment of dominant fluctuation weight to particular orbitals at each $\mathbf{q}$ (Yoshida et al., 2013, Saito et al., 2013).

3. Orbital Selectivity: Key Systems and Spectral Features

Iron-Based Superconductors

In BaFe$_2$(As,P)$_2$, spin susceptibility is largest at stripe-type wave vectors $\mathbf{Q}=(\pi,0)$ in $d_{xz}$ and $d_{yz}$ intra-orbital channels, with $d_{xy}$ also contributing substantially but less strongly. Orbital selectivity is crucial for the nodal structure of the superconducting gap due to competition between repulsive spin-fluctuation mediated pairing (intra-orbital, sign-changing) and attractive orbital-fluctuation (inter-orbital, sign-preserving) channels (Saito et al., 2013, Yoshida et al., 2013):

• $d_{xz}$/ $d_{yz}$ dominate $\chi^s(\mathbf{Q})$
• $d_{xy}$ controls the appearance of loop-shaped gap nodes on electron pockets when spin/orbital fluctuation amplitudes are comparable

Ruthenates

For Sr$_2$RuO$_4$, DFT+DMFT calculations show the spin susceptibility peak at incommensurate $\mathbf{Q}_i=(0.3,0.3,0)$ is an equal-weight superposition of all three $t_{2g}$ orbitals: $$v_{d_{xz}}(Q_i) = v_{d_{yz}}(Q_i) = v_{d_{xy}}(Q_i) = 0.41$$ revealing a genuinely cooperative multi-orbital origin of dynamic spin fluctuations (Boehnke et al., 2018).

Heavy-Fermion and Strong Spin-Orbit Systems

In SmB$_6$, orbital-resolved spectral weights $S^{\alpha\beta}(\mathbf{q},\omega)$ reveal that certain spin excitations are dominated by intra-orbital channels (trace $\sim$ sum), while others have large inter-orbital (spin–orbital entangled) character (sum $\gg$ trace), highlighting nontrivial mixing and entanglement between $f$-orbital multiplets due to strong Coulomb and relativistic couplings (Singh et al., 2018). In iridates, spin–orbital fluctuations between $J_{\mathrm{eff}}=1/2$ and $J_{\mathrm{eff}}=3/2$ sectors are suppressed by Mott gap formation and spin–orbit coupling, freezing orbital-resolved spin dynamics in the correlated insulator (Martins et al., 2011).

4. Experimental Probes and Selection Rules

Resonant Inelastic X-ray Scattering (RIXS)

Polarization- and momentum-resolved RIXS at the $L_3$ edge provides a powerful means for orbital-resolved measurement of spin fluctuations. The RIXS cross section in the spin-flip channel can be written as

$$I^{\text{spin}}_{\text{RIXS}}(\mathbf{q},\omega) \propto -\Im \left\{ \sum_{\mu\nu,\mu'\nu',\alpha\beta} [s^\alpha_{\mu\nu}]^* s^\beta_{\mu'\nu'} \chi^{\alpha\beta}_{\mu\nu,\mu'\nu'}(\mathbf{q},\omega) \right\}$$

where $s^\alpha_{\mu\nu}$ is a polarization- and geometry-dependent form factor that enables selective probing of, e.g., $d_{xz}$, $d_{yz}$, or $d_{xy}$ spin-fluctuations by adjusting the outgoing photon direction and scattering geometry. In iron-based systems, $\sigma\to\pi'$ geometry with $k_{\mathrm{out}}$ aligned along $x$, $y$, or $z$ enables strict orbital resolution (Yao et al., 2016).

In cuprates, full polarization control in RIXS on NdBa$_2$Cu$_3$O$_{7-\delta}$ unequivocally disentangles $dd$- and spin-excitations, permitting identification of softening and broadening of $d_{xy}$ features upon hole doping (Fumagalli et al., 2019).

Electron Spin Resonance (ESR) and Other Techniques

Multi-frequency ESR in hexagonal Ba$_3$CuSb$_2$O$_9$ accesses orbital dynamics via $g$-factor anisotropy and motional narrowing, with the orbital fluctuation time $\tau_c$ ($\sim$100 ps) directly extracted from frequency-dependent crossovers in lineshape, indicating the presence of an orbital-liquid regime dynamically modulating local spin exchange (Han et al., 2015).

5. Impact on Magnetism, Nematicity, and Pairing

Orbital-resolved spin fluctuations underlie:

• the stabilization of exotic magnetic states, including four-sublattice noncollinear phases in Kugel-Khomskii models (driven by spin–orbital entanglement and higher-order superexchange couplings) (Brzezicki et al., 2012, Brzezicki et al., 2014)
• the emergence of orbital-selective Mott phases, nematic order, and directionally dependent magnetic response in multi-band correlated metals and insulators
• the formation of unconventional superconducting states whose gap anisotropy, presence or absence of nodes, and orbital selectivity depend sensitively on the detailed structure of orbital-resolved $\chi^s(\mathbf{q},\omega)$ (Saito et al., 2013, Yoshida et al., 2013, Boehnke et al., 2018, Yao et al., 2016)

In one-dimensional systems, orbital fluctuations produce satellite spin–orbital branches and Kohn anomalies in the spin-wave spectrum due to the interplay with fermionized orbital excitations (Herzog et al., 2011).

6. Theoretical and Computational Approaches

Comprehensive theoretical analysis employs:

• Multi-orbital RPA, DFT+DMFT, and Bethe–Salpeter techniques for quantitative computation of $\chi^s_{l_1 l_2,l_3 l_4}(\mathbf{q},\omega)$ (Saito et al., 2013, Boehnke et al., 2018, Singh et al., 2018)
• Perturbative expansions and cluster mean-field/entanglement renormalization to capture entangled spin–orbital excitations and emergent longer-range couplings in orbitally degenerate Mott insulators (Brzezicki et al., 2012, Brzezicki et al., 2014)
• Analytical selection rules and matrix-element engineering for experimental orbital selectivity in photon or neutron probes (Yao et al., 2016, Fumagalli et al., 2019)

Table: Representative Orbital Contributions to Spin Susceptibility Peaks

Material (orbital set)	Peak $\mathbf{q}$	Dominant orbital channel(s)	Reference
BaFe$_2$(As,P)$_2$ (Fe 3d)	$(\pi,0)$	$d_{xz}$, $d_{yz}$ (large); $d_{xy}$ (medium)	(Saito et al., 2013)
Sr$_2$RuO$_4$ (Ru $t_{2g}$)	$(0.3,0.3,0)$	$d_{xz}$, $d_{yz}$, $d_{xy}$ (equal)	(Boehnke et al., 2018)
SmB$_6$ (Sm 4f, 5d)	various	$J=5/2$ intra- vs. inter-orbital mixing	(Singh et al., 2018)

7. Outlook and Significance

Orbital-resolved spin fluctuations are a central organizing principle in correlated electron physics. The ability to selectively compute, measure, and manipulate the detailed orbital content of spin dynamics enables precision control over emergent phenomena such as high-temperature superconductivity, quantum spin liquids, orbital nematics, and multipolar order. Recent advances in both theory (DFT+DMFT, RPA, tensor network states) and experiment (polarized RIXS, multi-frequency ESR) ensure this field remains at the forefront of quantum materials research.