---
title: Superexchange Interactions in Magnetic Systems
url: https://www.emergentmind.com/topics/superexchange-interactions
type: topic
---

# Superexchange Interactions in Magnetic Systems

Superexchange interactions are indirect magnetic couplings produced by virtual charge motion through high-energy intermediate states. In the canonical insulating case, localized moments on magnetic sites couple through nonmagnetic ligands, so that charge fluctuations are frozen in the ground state but still generate an effective spin Hamiltonian after downfolding. Contemporary work shows that this mechanism is broader than the textbook ligand-mediated antiferromagnetic exchange of simple Mott insulators: superexchange can be ferromagnetic, bond-directional, multi-anion, controlled by nominally nonmagnetic cations, reshaped by spin-orbit coupling, and engineered in quantum dots, driven optical lattices, and trapped-ion polariton systems [2102.04607] [1912.13274] [1910.07118] [1405.6071].

## 1. Microscopic formulation and effective-spin reduction

At the formal level, superexchange is a low-energy consequence of eliminating virtual charge excitations. In insulating \(d\)-\(p\)-\(d\) models this frequently appears at fourth order in hopping, whereas in already localized manifolds—such as polaritonic Mott states or detuned charge states of quantum dots—it can appear as a second-order process in the residual hopping. A compact example is the single-ligand antiferromagnetic model of a Mott insulator in a DC electric field, where the effective Hamiltonian takes the form \(\mathcal H_{\rm A}^{\rm eff}=J_{\rm A}(\mathbf E)\,\mathbf S_0\cdot\mathbf S_1-h(S_0^z+S_1^z)\), with \(J_{\rm A}(\mathbf 0)>0\) under the conditions \(\Delta_j(\mathbf 0)\ge 0\) and \(U_d>U_p>0\) [2102.04607].

In multi-orbital charge-transfer systems, the exchange constant depends explicitly on the virtual hopping geometry. For one-dimensional corner-shared cuprates, the effective antiferromagnetic scale is written as
\[
J=\frac{t_{pd,1}^2 t_{pd,2}^2}{\Delta_{dp}^2}\left(\frac{1}{\Delta_{dp}+U_p/2}+\frac{1}{U_d}\right),
\]
so that both Cu–O hoppings and the charge-transfer gap enter directly [2205.12678]. This structure already makes clear that superexchange is not a purely geometric label; it is a perturbative amplitude built from orbital-resolved transfer integrals and interaction energies.

The resulting spin model is not unique across the literature, and even the sign convention is model dependent. In graphene trimers the effective Heisenberg model is written as \({\cal H}=-(J/2)\sum \vec S_i\cdot\vec S_j\), so \(J<0\) denotes antiferromagnetic coupling [1808.09343]. In EuZn\(_2\)P\(_2\), by contrast, the effective model is \(\mathcal H_{\rm eff}=-\sum J_{ij}\,{\bm e}_i{\bm e}_j\), so positive \(J\) favors parallel alignment [2307.11924]. For this reason, the physical content of a reported “positive” or “negative” exchange constant must always be read together with the Hamiltonian convention.

## 2. Geometry, bridge chemistry, and the sign of exchange

The sign and magnitude of superexchange are controlled by the geometry and electronic activity of the bridge. In EuZn\(_2\)P\(_2\), the intralayer Eu–P–Eu path has bond angle \(85.52^\circ\) and is assigned a ferromagnetic sign, whereas the interlayer Eu–P–P–Eu path has Eu–P–P angle \(105.90^\circ\) and is assigned an antiferromagnetic sign. DFT mapping gives \(J_{\parallel}=1.23\) meV and \(J_{\perp}=-0.35\) meV, and Monte Carlo with these two couplings yields \(T_{\rm N}^{\rm MC}\approx 21\) K, close to the experimental \(T_{\rm N}=23.5\) K [2307.11924]. This is an “extended” superexchange in the paper’s terminology, because the interlayer path involves two phosphorus atoms rather than a single anion.

Bridge identity can be as important as bond angle. In the double perovskites Sr\(_2\)CuTeO\(_6\) and Sr\(_2\)CuWO\(_6\), formally nonmagnetic Te\(^{6+}\) and W\(^{6+}\) produce qualitatively different super-superexchange networks. Quantum-chemistry calculations give \(J_1=7.38\) meV and \(J_2=0.05\) meV for Sr\(_2\)CuTeO\(_6\), but \(J_1=0.68\) meV and \(J_2=8.33\) meV for Sr\(_2\)CuWO\(_6\). The reason is that low-lying empty W \(5d\) states strongly enhance the \(180^\circ\) Cu–O–W–O–Cu path, whereas in the Te compound the dominant \(J_1\) is mainly a Cu–O–O–Cu process and the Te-centered diagonal path remains weak [1902.09376]. This is a direct demonstration that nominally nonmagnetic cations can be active elements of the exchange bridge rather than passive spacers.

Near-\(90^\circ\) geometries do not force a unique outcome. In CaMnCrSbO\(_6\), the nearest-neighbor Mn–O–Mn and Mn–O–Cr channels are all antiferromagnetic but weak, with fitted values such as \(J'_{1a}\approx -0.41\) meV and \(J'_{1b}\approx -0.06\) meV for Mn–Cr exchange, while the next-nearest-neighbor Cr–O–O–Cr super-superexchange is ferromagnetic with \(J'_2\approx +0.32\) meV. The paper argues that the correct sign of this Cr–O–O–Cr term requires a four-site model appropriate to doubly occupied oxygen \(p\) orbitals rather than the usual hole-model formula, which would give the wrong sign [2205.02452].

The strongest departure from textbook Anderson superexchange is the explicit ferromagnetic kinetic channel identified for diamagnetic metal bridges. In a three-orbital model with two localized magnetic orbitals and one bridge orbital, the exchange contains a third-order term
\[
J_{K3}=-\frac{4t_{MD}t_{DM}t_{MM}}{\Delta^2},
\]
which is ferromagnetic when \(t_{MD}t_{DM}t_{MM}>0\). In the Fe\(^{3+}\)–Co\(^{3+}\)–Fe\(^{3+}\) complex analyzed in detail, the decomposition yields \(K1=27.1\) meV, \(K2=23.7\) meV, \(K3=-75.1\) meV, \(K4=-6.2\) meV, \(PE=-5.1\) meV, and total \(J=-10.4\) meV, so the large ferromagnetic coupling is driven predominantly by the cyclic kinetic term rather than by conventional potential exchange [1912.13274]. A common misconception is therefore that kinetic exchange through diamagnetic bridges is generically antiferromagnetic; this paper shows that the conclusion depends on how explicitly the bridge orbitals and magnetic-center transfer \(t_{MM}\) are retained.

## 3. Spin-orbit-entangled and anisotropic superexchange

When spin-orbit coupling is large, superexchange is not exhausted by isotropic Heisenberg exchange. In the \(5d^5\) \(t_{2g}\) problem treated in the \(J\)-\(J\) coupling scheme, the local basis is reorganized into \(\Gamma_7\) and \(\Gamma_8\) spin-orbit-entangled orbitals. Two virtual channels then appear naturally: \(\Gamma_7\)-to-\(\Gamma_8\) processes generate a ferromagnetic Ising interaction identified with a quantum-compass term, whereas \(\Gamma_7\)-to-\(\Gamma_7\) processes generate antiferromagnetic Heisenberg exchange. The paper further argues that increasing the spin-orbit coupling \(\zeta\) suppresses the ferromagnetic Ising contribution, so the dominant exchange changes from ferromagnetic Ising type to antiferromagnetic Heisenberg type as \(\zeta\) increases [1407.0811].

A second development is the explicit role of spin-orbit coupling on the nonmagnetic bridge itself. In a multi-orbital Hubbard model with SOC on both magnetic cations and nonmagnetic anions, local SOC rotates the hopping into the form \(b+\mathbf C\cdot\boldsymbol\sigma\), so anisotropic exchange emerges already at the level of the effective hopping matrices. The resulting fourth-order superexchange contains the isotropic Heisenberg term, a Dzyaloshinskii–Moriya vector \(\mathbf D_{ij}\), a symmetric anisotropic tensor \(\mathbf\Gamma_{ij}\), and single-ion anisotropy; in the weak-SOC limit, \(\mathbf D_{ij}\) is linear in SOC and \(\mathbf\Gamma_{ij}\) quadratic. The central claim is that SOC on nonmagnetic anions can induce these anisotropies on an equal footing with SOC on magnetic ions [1910.07118].

This shift in viewpoint is significant for heavy-ligand compounds. It implies that ligand SOC is not merely a correction to an otherwise cation-controlled exchange problem, but can be a primary source of bond-directional and antisymmetric exchange. A plausible implication is that the conventional separation between “exchange geometry” and “single-ion spin-orbit physics” becomes inadequate once the bridge orbitals themselves are spin-orbit active.

## 4. Exchange topology, dimensionality, and materials realization

Superexchange determines not only the sign of a pairwise coupling but also the dimensionality of the magnetic network. In tetragonal multiferroic \(0.50\mathrm{BiFeO_3}-0.50\mathrm{PbTiO_3}\), reducing the particle size weakens the ferroelectric distortion and shortens the long apical Fe–O bond from \(2.571\) Å at \(120\) nm to \(2.348\) Å at \(18\) nm, while tetragonality decreases from \(14.16\%\) to \(5.75\%\). The authors interpret this as restoring the Fe–O–Fe exchange path along the \(c\) axis, producing a crossover from essentially \(2d\) antiferromagnetic interactions in bulk-like particles to \(3d\) antiferromagnetic interactions in nanoparticles. Experimentally, \(T_{\rm N}\) rises from \(\sim 120\) K to \(\sim 350\) K, opposite to ordinary finite-size suppression in non-multiferroic antiferromagnets [1407.6485].

In iron pnictides, the dominant exchange path can be shifted away from nearest neighbors. For LaFePO, first-principles calculations give \(J_1\sim 0.5\) meV/\(S^2\) and \(J_2\sim 7.5\) meV/\(S^2\), with the next-nearest-neighbor \(J_2\) attributed to P \(3p\)-bridged superexchange on the Fe square lattice. The small \(J_1\) and dominant \(J_2\) stabilize the collinear stripe state and are tied to the same pnicogen-mediated exchange picture previously advanced for LaFeAsO [1009.0726]. In low-symmetry two-dimensional magnets, the same logic generalizes from “which sign?” to “which Hamiltonian?”: in monolayer Cr\(_3\)Te\(_6\), DFT plus Wannier path counting yields six distinct nearest-neighbor exchange constants rather than a single \(J_1\), with \(J_{23}S^2=-37.4\) meV the strongest term and a calculated \(T_c=328\) K close to the reported \(344\) K [2301.01923].

Carbon-based systems illustrate how interface symmetry reshapes superexchange. In zigzag graphene nanoribbons embedded in graphane, symmetric C–CH/C–CH interfaces generate strong, long-range interedge antiferromagnetic superexchange for all \(N\ge 2\), with \(\Delta E_{\mathrm{NM-AFM}}\) rising from \(8.0\) meV/unit cell at \(N=2\) to \(72.7\) meV at \(N=8\). In zigzag ribbons embedded in h-BN, by contrast, asymmetric C–B and C–N edges make the interedge coupling short-ranged and weak, so the ground state is nonmagnetic except at \(N=5\) and \(6\), where only a weakly stabilized half-semimetallic state appears with \(\Delta E_{\mathrm{NM-HS}}=1.7\) and \(1.5\) meV/unit cell [1504.02180]. In graphene trimers, the interpretation is even less canonical: Cr and Mn trimers are antiferromagnetic-like at \(N=3\) but ferromagnetic at \(N=5\) and \(7\), and the authors describe the mechanism as “RKKY-like super-exchange” because it combines substitutional local bonding with host-mediated, metallicity-dependent, distance-sensitive coupling rather than fitting a purely local superexchange picture [1808.09343].

Fluoride polymorphs show a related control of exchange topology by structural reorganization. In AgMF\(_4\) (\(M=\) Ni, Cu), quasi-linear Ag–F–M bridges become the dominant antiferromagnetic pathways once the lattice connectivity departs from stacked AgF\(_2\)-like layers. Low-pressure AgCuF\(_4\) remains mainly layer-dominated, with strong Ag–Ag and Cu–Cu superexchange, but compression yields \(J_1^{\rm mix}=-45.9\) meV in AgCuF\(_4\)-HP. AgNiF\(_4\) likewise develops dominant mixed-chain couplings, with \(J_1^{\rm mix}=-32.8\) meV in the low-pressure phase and \(-33.3\) meV in the high-pressure phase [2108.12241].

## 5. External fields, lattice dynamics, and nonequilibrium control

Because superexchange depends on virtual excitation energies, it is electrically tunable even in insulators. A general DC-field theory shows that static electric fields modify the site-energy differences \(\Delta_j(\mathbf E)\) entering the virtual processes, thereby changing both antiferromagnetic and ferromagnetic superexchange. On the square lattice, this permits control of \(J_1/J_2\) and can even split the two diagonal couplings into \(J_2\pm \delta J_2\), effectively deforming a square-lattice \(J_1\)-\(J_2\) magnet into a triangular one. On a spin chain, a field-induced bond alternation \(\delta J_1(\mathbf E)\) produces singlet-dimer or Haldane-dimer order; for the Heisenberg antiferromagnetic chain the induced gap scales as \(\Delta(\mathbf E)\propto |E|^{2/3}\). The paper estimates observable field scales of order \(1~\mathrm{MV/cm}\) for inorganic Mott insulators and \(0.1~\mathrm{MV/cm}\) for organic ones, and gives \(\Delta\sim 22\) K for \(\mathrm{KCuMoO_4(OH)}\) at \(1~\mathrm{MV/cm}\), about \(8.6\%\) of \(J_1(\mathbf 0)=238\) K [2102.04607].

Lattice fluctuations can renormalize the same virtual processes in the opposite direction. In a four-orbital Hubbard–Su–Schrieffer–Heeger model for one-dimensional cuprates, bond-stretching oxygen phonons modulate the two Cu–O hoppings oppositely, so the exchange is renormalized as \(J\to J(1-g^2x^2)^2\). Exact diagonalization and determinant quantum Monte Carlo show that increasing the electron-phonon coupling suppresses the effective superexchange, softens the spin spectrum, and enhances the uniform susceptibility. A narrow dimerized regime exists only just below a critical coupling \(g_c\sim 0.8\); beyond that, the linear SSH model becomes unstable because the effective hopping changes sign [2205.12678]. This directly contradicts the common expectation that phonons necessarily enhance exchange by reducing an effective \(U\); in the charge-transfer geometry studied here, the dominant effect is instead the reduction of the hopping product entering superexchange.

Time-periodic driving provides a nonequilibrium control knob. In AC-driven optical double wells, a magnetic-field gradient \(G\) suppresses static superexchange, but lattice modulation restores it resonantly at \(\hbar\omega=2G\). In the symmetric case \(\Delta=0\), the photon-assisted transverse exchange is
\[
J_\mathrm{ex}^\perp=8JK\frac{U}{U^2-G^2},
\]
while the longitudinal term is
\[
J_\mathrm{ex}^z=\frac{4J^2}{U}\frac{U^2}{U^2-G^2}.
\]
This produces an effective XXZ Hamiltonian with independently tunable transverse and Ising couplings, and the experiment reports a one-photon superexchange resonance at \(\omega/2\pi=2.6(1)\,\mathrm{kHz}\) with \(J_\mathrm{ex}^\perp/h=560(20)\,\mathrm{Hz}\) [1104.1833].

## 6. Mesoscopic and synthetic realizations

In semiconductor triple quantum dots, superexchange becomes a directly gate-defined long-range resource. For a two-electron linear triple dot with an empty mediator, full configuration interaction finds that the outer-spin superexchange \(J_\mathrm{se}=E_{|T(101)\rangle'}-E_{|S(101)\rangle'}\) can be non-monotonic and can switch sign as a function of middle-dot detuning when the outer-dot detunings are leveled. The same study shows that increasing left-right detuning strengthens an originally positive \(J_\mathrm{se}\) and weakens an originally negative one, and that simplified Hubbard descriptions miss this behavior unless long-range Coulomb exchange \(J^e_{\mathrm{L,R}}\) is included [2203.15521].

Mediator occupancy adds a second control axis. Configuration-interaction calculations for a larger central dot show that an empty mediator produces at most about a \(20\%\) enhancement of the remote coupling, whereas loading two electrons into the mediator can enhance the outer-dot exchange by more than two orders of magnitude at small angle. In the linear geometry, the coupling remains above \(1\,\mu\mathrm{eV}\) out to roughly \(130\) nm, compared with about \(56\) nm for the direct double-dot reference [2003.03416]. With a four-electron mediator, the sign itself becomes occupancy and field dependent: a two-electron singlet mediator gives positive exchange, a four-electron singlet mediator negative exchange, and a four-electron triplet mediator under a non-uniform magnetic field can be tuned through zero as \(B_{L/R}\) is increased relative to \(B_M\); leakage from the logical subspace is estimated to fall below \(10^{-3}\) when \(B_{L/R}-B_M\gtrsim 0.42\) T [2204.02615].

Quantum-dot superexchange can also display interference. In a triple dot with two distinct virtual intermediate states, \((1,0,1)\) and \((0,2,0)\), the effective singlet amplitude is
\[
\tau_\text{eff}^\text{S}=\tau_\text{LC}\tau_\text{CR}\left(\delta_{101}^{-1}+2\delta_{020}^{-1}\right),
\]
whereas the triplet amplitude lacks the \((0,2,0)\) path. Destructive interference at \(\delta_{101}=-\delta_{020}/2\) gives \(\tau_\text{eff}^\text{S}=0\), producing a dark state and a “superexchange blockade” of current [1312.2512]. This is a mesoscopic analogue of path-interference control over an otherwise conventional virtual-exchange mechanism.

Synthetic many-body platforms realize the same logic with different microscopic carriers. In a linear trapped-ion crystal, a V-type Jaynes–Cummings–Hubbard model creates Mott-localized ion-phonon polaritons, and second-order phonon hopping generates effective spin models. At one excitation per site the low-energy theory is an XXZ spin-\(\tfrac12\) model; at two excitations per site it becomes a spin-1 Heisenberg-type model with additional anisotropies. The exchange strengths are tunable through trap frequencies, laser intensities, and detuning [1405.6071].

Taken together, these results suggest that “superexchange interaction” is best understood not as a single fixed antiferromagnetic mechanism, but as a family of virtual-exchange processes whose sign, anisotropy, range, and dimensionality are set by the detailed structure of the intermediate manifold—ligand identity, bridge multiplicity, orbital symmetry, spin-orbit entanglement, lattice dynamics, and external control parameters.

Source: https://www.emergentmind.com/topics/superexchange-interactions