Skip to content

Exchange Interactions

A magnetic exchange interaction is the dependence of a system’s low-energy electronic energy on the relative quantum state of magnetic moments. It is not an additional fundamental force placed beside electromagnetism. Rather, fermionic antisymmetry changes the allowed many-electron states, while Coulomb matrix elements, orbital overlap, hopping, charge excitation energies, itinerant response, and spin–orbit coupling determine whether those states have different energies.

The result is often summarized by a spin Hamiltonian,

Hex=∑i<jJij si⋅sj,H_{\mathrm{ex}} = \sum_{i<j} J_{ij}\, \mathbf s_i\cdot\mathbf s_j,

but the number JijJ_{ij} is the end of a reduction, not the beginning of an explanation. Different microscopic pathways can produce the same apparent bilinear coupling. Several pathways can coexist on one bond, and some mechanisms do not reduce accurately to a field-independent pairwise JijJ_{ij} at all.

The central question is which electronic process couples the moments. Direct exchange, superexchange, double exchange, and RKKY coupling arise in different limits, with different signs and energy scales. Exchange and Correlation derives atomic Coulomb integrals; the Heisenberg Model develops the spin model that often results.

Required background. Spin and Spatial Wavefunctions explains the required fermionic symmetry, while the Hubbard Model supplies hopping and charge-excitation scales used in virtual exchange.

Helpful background. Singlet and Triplet States fixes coupled-spin basis conventions, while the Schrieffer–Wolff Transformation supplies a controlled low-energy projection.

For two identical electrons, the total state must be antisymmetric under particle exchange. A spin singlet is antisymmetric in spin and must therefore be paired with a symmetric spatial state. A spin triplet is symmetric in spin and must be paired with an antisymmetric spatial state.

Let a(r)a(\mathbf r) and b(r)b(\mathbf r) be orthonormal spatial orbitals. The normalized exchange eigenstates are

Φ±(1,2)=12[a(1)b(2)±b(1)a(2)].\Phi_{\pm}(1,2) = \frac{1}{\sqrt2} \left[ a(1)b(2) \pm b(1)a(2) \right].

For a spin-independent two-body interaction V(1,2)V(1,2), define the direct and exchange matrix elements

Cab=∬d1 d2 ∣a(1)∣2V(1,2)∣b(2)∣2,C_{ab} = \iint d1\,d2\, |a(1)|^2 V(1,2) |b(2)|^2, Kab=∬d1 d2 a∗(1)b∗(2)V(1,2)b(1)a(2).K_{ab} = \iint d1\,d2\, a^*(1)b^*(2) V(1,2) b(1)a(2).

Their interaction energies are

⟨Φ±∣V∣Φ±⟩=Cab±Kab.\langle\Phi_{\pm}|V|\Phi_{\pm}\rangle = C_{ab}\pm K_{ab}.

For electrons, Φ+\Phi_+ accompanies the spin singlet and Φ−\Phi_- the triplet. The interaction contribution to the singlet–triplet splitting is therefore

ET−ES=−2Kab.E_T-E_S = -2K_{ab}.

For the positive Coulomb kernel and fixed orbitals, Kab≥0K_{ab}\ge0, so this contribution favors the triplet or leaves it degenerate. That conclusion is narrower than the slogan “Pauli exclusion causes ferromagnetism.” The full energy also contains kinetic, one-body, orthogonalization, hybridization, and correlation contributions. Those can reverse the net splitting.

Most importantly, antisymmetry by itself is kinematic. If the relevant orbitals have disjoint support and no process connects them, KabK_{ab} and the hopping matrix element vanish. The spin states can then remain degenerate even though the total wavefunction is still exactly antisymmetric. An exchange splitting requires state-dependent Hamiltonian matrix elements.

Exchange signs, spin units, and bond counting

Section titled “Exchange signs, spin units, and bond counting”

Throughout this page, si\mathbf s_i is a dimensionless spin operator:

Si=ℏsi,si2=si(si+1).\mathbf S_i = \hbar\mathbf s_i, \qquad \mathbf s_i^2 = s_i(s_i+1).

The isotropic pair Hamiltonian is

Hij=Jij si⋅sj.H_{ij} = J_{ij}\, \mathbf s_i\cdot\mathbf s_j.

With this convention:

  • Jij>0J_{ij}>0 favors antiparallel correlations and is called antiferromagnetic;
  • Jij<0J_{ij}<0 favors parallel correlations and is called ferromagnetic;
  • JijJ_{ij} has units of energy.

For two spin-1/21/2 moments,

s1⋅s2=12[stot2−s12−s22].\mathbf s_1\cdot\mathbf s_2 = \frac12 \left[ \mathbf s_{\mathrm{tot}}^2 - \mathbf s_1^2 - \mathbf s_2^2 \right].

The singlet and triplet energies are

ES=−3J4,ET=J4,E_S = -\frac{3J}{4}, \qquad E_T = \frac{J}{4},

so

J=ET−ES.J = E_T-E_S.

This dimer identity is the simplest operational definition of the sign.

The permutation operator on two spin-1/21/2 states is

P12=2s1⋅s2+12.P_{12} = 2\mathbf s_1\cdot\mathbf s_2 + \frac12.

It has eigenvalue −1-1 on the singlet and +1+1 on the triplet. This identity explains why a two-state exchange splitting can be rewritten as a Heisenberg coupling, but it does not identify the microscopic mechanism that generated the splitting.

Other common conventions include

H=−∑i<jJijsi⋅sjH = - \sum_{i<j} \mathcal J_{ij} \mathbf s_i\cdot\mathbf s_j

and Hamiltonians written using unit vectors ei\mathbf e_i rather than spin operators. Some sums count every bond twice. Consequently, a quoted positive “exchange constant” may mean ferromagnetic coupling in one paper and antiferromagnetic coupling in another.

For example, changing H=Jsi⋅sjH=J\mathbf s_i\cdot\mathbf s_j to H=−Jei⋅ejH=-\mathcal J\mathbf e_i\cdot\mathbf e_j, with unit vectors ei=si/s\mathbf e_i=\mathbf s_i/s, gives J=−Js2\mathcal J=-Js^2 in a classical fixed-length description. An ordered-pair sum instead of one unordered-bond sum adds another factor of two. Comparing the Hamiltonians resolves the sign and normalization before any numerical values are compared.

Direct, superexchange, double-exchange, and RKKY pathways between magnetic moments

Four routes by which electronic motion and interaction couple magnetic moments. The superexchange arrows indicate virtual charge transfer; their direction depends on the occupied orbitals. The double-exchange carrier remains mobile. In RKKY coupling the intervening electron gas responds to a local spin. The diagram identifies processes, not universal signs for every material.

These mechanisms are not mutually exclusive labels for materials. A transition-metal oxide can have antiferromagnetic superexchange, ferromagnetic Hund coupling, weaker direct overlap, and longer-range itinerant exchange simultaneously. The effective JijJ_{ij} is the low-energy sum after the active orbitals and charge fluctuations have been specified.

Direct exchange couples magnetic orbitals that overlap without requiring a distinct ligand or conduction sea as mediator. At the two-orbital level, its characteristic Coulomb contribution is the exchange integral

Kij=∬dr dr′ ϕi∗(r)ϕj(r)e24πϵ∣r−r′∣ϕj∗(r′)ϕi(r′).K_{ij} = \iint d\mathbf r\,d\mathbf r'\, \phi_i^*(\mathbf r) \phi_j(\mathbf r) \frac{e^2}{ 4\pi\epsilon |\mathbf r-\mathbf r'| } \phi_j^*(\mathbf r') \phi_i(\mathbf r').

It depends on the transition density ϕi∗ϕj\phi_i^*\phi_j, so it decreases rapidly as localized orbitals separate. Orbital orientation, radial extent, screening, covalency, and orthogonalization all matter.

For fixed orthonormal orbitals, the bare Coulomb exchange term lowers the triplet relative to the singlet and maps to a ferromagnetic contribution,

Jdirect(V)=−2Kij.J_{\mathrm{direct}}^{(V)} = -2K_{ij}.

This equation should not be promoted into a universal sign theorem for an interatomic bond. Once orbital relaxation and hopping are allowed, virtual charge fluctuations often add an antiferromagnetic term. What a fitted model calls “nearest-neighbor exchange” may be the sum of both.

Different orbitals on the same ion also have direct Coulomb exchange. A common local convention writes

HHund=−JH∑m<m′sim⋅sim′,JH>0.H_{\mathrm{Hund}} = - J_{\mathrm H} \sum_{m<m'} \mathbf s_{im}\cdot\mathbf s_{im'}, \qquad J_{\mathrm H}>0.

This favors a high-spin atomic multiplet. Hund exchange is crucial in both superexchange and double exchange, but it is not identical to an intersite Heisenberg coupling. Crystal-field splitting, pair hopping, and the full multiorbital Coulomb tensor are needed for quantitative work.

Superexchange is an indirect coupling generated by virtual hopping through high-energy charge configurations. It is especially important in magnetic insulators, where real dc charge transport is suppressed but virtual charge transfer remains allowed.

The cleanest example is the half-filled, one-orbital two-site Hubbard model. Choose the phase of the orbitals so that its single hopping amplitude tt is real:

H=−t∑σ(c1σ†c2σ+c2σ†c1σ)+U∑i=12ni↑ni↓,U>0.H ={} - t \sum_{\sigma} \left( c_{1\sigma}^{\dagger}c_{2\sigma} + c_{2\sigma}^{\dagger}c_{1\sigma} \right) + U \sum_{i=1}^{2} n_{i\uparrow}n_{i\downarrow}, \qquad U>0.

At half filling and ∣t∣≪U|t|\ll U, the low-energy subspace has one electron on each site. A hop temporarily creates one empty site and one doubly occupied site, costing energy UU. A second hop returns to the low-energy sector.

To second order,

Heff=4∣t∣2U(s1⋅s2−14)H_{\mathrm{eff}} = \frac{4|t|^2}{U} \left( \mathbf s_1\cdot\mathbf s_2 - \frac14 \right)

within the singly occupied subspace. The constant is chosen so that the triplet energy shift is zero. The singlet is lowered by 4∣t∣2/U4|t|^2/U, and therefore

JSE=4∣t∣2U>0.J_{\mathrm{SE}} = \frac{4|t|^2}{U} >0.

This antiferromagnetic sign follows from Pauli blocking and virtual kinetic energy: the singlet can access the relevant doubly occupied intermediate state, whereas the triplet cannot in a one-orbital model.

The result is controlled only when charge excitations remain well separated:

∣t∣U≪1.\frac{|t|}{U} \ll 1.

Near a metal–insulator transition, higher orders, longer-range hopping, ring exchange, and residual charge fluctuations can become important.

In a transition-metal–ligand–transition-metal path, virtual charge transfer between metal dd and ligand pp states connects the two magnetic sites. The allowed direction depends on occupancy: for a filled ligand shell, a first electron hop can leave a ligand hole, rather than add an electron to an already occupied spin orbital. If tpdt_{pd} is a metal–ligand hopping and Δ\Delta a charge-transfer energy, a schematic downfolding gives

teff∼tpd2Δ,t_{\mathrm{eff}} \sim \frac{t_{pd}^2}{\Delta},

and an antiferromagnetic scale of order

JAF∼4tpd4Δ2Ueff.J_{\mathrm{AF}} \sim \frac{4t_{pd}^4}{ \Delta^2 U_{\mathrm{eff}} }.

This scaling identifies the small parameters, not a complete material formula. Real denominators include metal and ligand Coulomb repulsions, crystal-field levels, Hund couplings, and multiple orbital paths.

Orbital symmetry provides useful qualitative rules:

  • a nearly 180∘180^\circ path between half-filled orbitals that hybridize with the same ligand orbital commonly gives strong antiferromagnetic exchange;
  • a near-90∘90^\circ geometry can suppress that shared-orbital path, allowing ligand or metal Hund coupling to favor ferromagnetic exchange;
  • half-filled–empty and half-filled–full combinations can have signs different from the half-filled–half-filled case;
  • distortions can change both hopping magnitude and the active orbital channel.

These are chemistry-aware limiting rules, not a lookup table based on bond angle alone. The valence, orbital occupancy, local axes, ligand states, and competing paths must all be specified.

Double exchange couples mixed-valence local moments through the kinetic energy of a mobile carrier. A minimal form contains hopping and a strong ferromagnetic Hund coupling between the carrier spin and a core moment:

H=−∑⟨i,j⟩,σ(tijciσ†cjσ+tij∗cjσ†ciσ)−JH∑iℓi⋅si,JH>0.H ={} - \sum_{\langle i,j\rangle,\sigma} \left( t_{ij} c_{i\sigma}^{\dagger}c_{j\sigma} +t_{ij}^{\ast} c_{j\sigma}^{\dagger}c_{i\sigma} \right) - J_{\mathrm H} \sum_i \boldsymbol\ell_i\cdot\mathbf s_i, \qquad J_{\mathrm H}>0.

Here ℓi\boldsymbol\ell_i is a dimensionless core spin of length ℓ\ell, and si=ci†σci/2\mathbf s_i=c_i^\dagger\boldsymbol\sigma c_i/2 is the dimensionless carrier spin. Each unordered bond is counted once. For the following semiclassical result, treat the cores as fixed classical directions and take JHℓ≫∣t∣J_{\mathrm H}\ell\gg |t|: the carrier spin then follows the local core. Anderson and Hasegawa also analyze quantum cores; their finite-spin problem is not exactly the same as two fixed directions. If neighboring core moments subtend an angle θij\theta_{ij}, the overlap of their aligned spinors has magnitude

∣⟨χi|χj⟩∣=cos⁡θij2.\left| \left\langle \chi_i \middle| \chi_j \right\rangle \right| = \cos\frac{\theta_{ij}}{2}.

The projected amplitude is tijeff=tij⟨χi∣χj⟩t_{ij}^{\mathrm{eff}}=t_{ij}\langle\chi_i|\chi_j\rangle. Its magnitude is therefore

∣tijeff∣=∣tij∣cos⁡θij2.|t_{ij}^{\mathrm{eff}}| = |t_{ij}| \cos\frac{\theta_{ij}}{2}.

The complex phase also matters around a loop; replacing every amplitude by its magnitude would discard the associated geometric phase. For a single bond that phase can be removed by a change of orbital phases.

Parallel moments maximize each projected hopping magnitude and can lower the occupied-band kinetic energy. Antiparallel moments suppress hopping in the ideal limit. This produces a ferromagnetic tendency tied to carrier motion.

The comparison highlights the different physical processes:

Property Superexchange Double exchange
Charge process Virtual excursion and return Real carrier delocalization
Typical regime Integer filling, charge localized Mixed valence or partial filling
Energy scale Often proportional to t2/Ut^2/U Often proportional to occupied-band kinetic energy
Angular form Bilinear at leading order Can scale as cos⁡(θ/2)\cos(\theta/2) and be non-Heisenberg

The ferromagnetic tendency is not unconditional. At an empty or completely filled active band there may be no kinetic-energy gain, and antiferromagnetic superexchange between core moments can compete strongly. Lattice distortions, orbital order, polarons, and phase separation can also reshape the simple picture.

RKKY interaction through a conducting host

Section titled “RKKY interaction through a conducting host”

Localized moments can interact even when their orbitals scarcely overlap. Each polarizes the conducting electrons, whose response is felt by the other. Use a lattice-normalized coupling

HK=JK∑aℓa⋅sa,sa=12ca†σca,H_K=J_K\sum_a\boldsymbol\ell_a\cdot\mathbf s_a, \qquad \mathbf s_a=\frac12c_a^\dagger\boldsymbol\sigma c_a,

where both spins are dimensionless and JKJ_K has units of energy. Define the static host response to Hh=−∑aha⋅saH_h=-\sum_a\mathbf h_a\cdot\mathbf s_a by δ⟨sa⟩=∑bχs,abhb\delta\langle\mathbf s_a\rangle=\sum_b\chi_{s,ab}\mathbf h_b. For a spin-isotropic reciprocal host, its quadratic free-energy change is −12∑abχs,abha⋅hb-\tfrac12\sum_{ab}\chi_{s,ab}\mathbf h_a\cdot\mathbf h_b. Substitution of ha=−JKℓa\mathbf h_a=-J_K\boldsymbol\ell_a gives, for distinct sites,

HRKKY=∑a<bJabRKKYℓa⋅ℓb,JabRKKY=−JK2χs,ab.H_{\mathrm{RKKY}} =\sum_{a<b}J_{ab}^{\mathrm{RKKY}}\boldsymbol\ell_a\cdot\boldsymbol\ell_b, \qquad J_{ab}^{\mathrm{RKKY}}=-J_K^2\chi_{s,ab}.

This minus sign follows from the stated physical response convention. It is not an independent choice once HhH_h and the plus-sign pair Hamiltonian are fixed. Continuum contact couplings have different units because the electron spin is then a density; RKKY Interaction gives that normalization and the full range-function derivation.

For a clean three-dimensional isotropic free-electron host, the large-distance coupling oscillates as cos⁡(2kFR)/R3\cos(2k_{\mathrm F}R)/R^3, so different separations can favor parallel or antiparallel moments. In a crystal, the full Fermi surface determines periods and directions. Disorder, temperature, dimensionality, and spin–orbit coupling change the range or tensor structure.

This is a weak-coupling result. Kondo Effect describes the competing screening of a local moment; Kondo Lattices develops the competition in a dense array of moments.

Spin–orbit coupling and low bond symmetry allow a general bilinear tensor,

Hij=siTJijsj.H_{ij} = \mathbf s_i^{\mathsf T} \mathsf J_{ij} \mathbf s_j.

Decomposing Jij\mathsf J_{ij} into scalar, antisymmetric, and symmetric-traceless parts gives

Hij=Jijisosi⋅sj+Dij⋅(si×sj)+si⋅Γijsj.H_{ij} ={} J_{ij}^{\mathrm{iso}} \mathbf s_i\cdot\mathbf s_j + \mathbf D_{ij}\cdot \left( \mathbf s_i\times\mathbf s_j \right) + \mathbf s_i\cdot \mathsf\Gamma_{ij} \mathbf s_j.

Here:

  • Jijiso=Tr⁡Jij/3J_{ij}^{\mathrm{iso}}=\operatorname{Tr}\mathsf J_{ij}/3 is isotropic exchange;
  • Dij\mathbf D_{ij} is the Dzyaloshinskii–Moriya vector;
  • Γij\mathsf\Gamma_{ij} is symmetric and traceless anisotropic exchange.

Reversing the bond orientation gives Dji=−Dij\mathbf D_{ji}=-\mathbf D_{ij}; this compensates the sign change of si×sj\mathbf s_i\times\mathbf s_j.

An inversion center at the midpoint of a bond forbids Dij\mathbf D_{ij} for that bond. Broken inversion is necessary but not sufficient; the remaining point-group operations constrain its direction. Strong spin–orbit-coupled materials can also generate bond-dependent Ising, compass, or Kitaev-like terms. Spin–Orbit Coupling in Solids explains how local orbital physics and crystal symmetry feed these effective interactions.

Single-ion anisotropy, magnetic dipole interactions, and Zeeman coupling are not exchange, even though they coexist in the same spin Hamiltonian. Magnetic Anisotropy explains their orientation-dependent material consequences and their competition with anisotropic exchange.

An exchange parameter is controlled by several linked scales:

Direct overlap and hopping usually decrease rapidly with distance. Because superexchange is proportional to several hopping amplitudes, modest bond-length or bond-angle changes can strongly alter JJ. Pressure, epitaxial strain, and structural transitions can therefore tune exchange.

Virtual exchange denominators depend on onsite repulsion, charge-transfer energy, crystal-field splitting, and Hund coupling. Changing valence or gating a narrow band can move the system between superexchange-dominated and carrier-mediated regimes.

RKKY and itinerant exchange inherit wavevectors and anisotropy from the susceptibility χ(q)\chi(\mathbf q). Nesting, multiple pockets, and spin–orbit splitting can select noncollinear or long-period order.

Two bonds of equal length need not have equal exchange if their ligand chemistry differs. Several paths can cancel, making a small net JJ highly sensitive to structural details. A sign change is often evidence for competition, not for one mechanism changing its intrinsic nature.

Downfolding into different orbital sets redistributes physics among hopping, onsite interactions, and exchange terms. A static pairwise JijJ_{ij} extracted near one reference state can differ from a coupling fitted over large spin rotations or at another temperature. Exchange parameters are model coordinates, not basis-independent observables.

No single measurement reads out every JijJ_{ij} directly. Reliable inference combines complementary information.

For an isolated spin-1/21/2 dimer in this convention,

J=ET−ES.J = E_T-E_S.

Inelastic neutron scattering, electron-spin resonance, optical spectroscopy, or tunneling spectroscopy can resolve this gap. Extra anisotropy splits the triplet and must be included rather than absorbed blindly into one JJ.

In an ordered magnet, exchange controls the momentum dependence of spin-wave energies. Inelastic neutron scattering measures the dynamic structure factor S(q,ω)S(\mathbf q,\omega), allowing multiple neighbor couplings to be fitted across the Brillouin zone. Magnons develops the quantization and spectral-weight framework.

A zone-center stiffness alone usually constrains only a weighted moment of the JijJ_{ij}. Several exchange networks can share the same long-wavelength slope but differ at the zone boundary.

Susceptibility, specific heat, and ordering temperature constrain exchange scales, but they are inverse problems. Frustration, dimensionality, anisotropy, disorder, and quantum fluctuations can suppress the ordering temperature far below the largest bond energy. The estimate kBTord∼Jk_{\mathrm B}T_{\mathrm{ord}}\sim J is dimensional guidance, not an extraction formula.

One can calculate total energies for several magnetic configurations and fit them to a spin Hamiltonian. For one isolated classical pair of equal spin length ss,

E∥−Eantiparallel=2Js2.E_{\parallel} -E_{\mathrm{antiparallel}} = 2Js^2.

In a lattice, bond multiplicities and all simultaneously changed pairs must be counted. The fitted result depends on the spin normalization, reference state, electronic-structure approximation, orbital projection, and chosen interaction terms.

Infinitesimal-rotation or magnetic-force-theorem methods extract the curvature of electronic energy around a reference state. Correlated methods can add frequency-dependent self-energies and vertex information. Agreement with measured magnon dispersions is a stronger test than agreement with one collinear energy difference.

Competing exchange paths and the limits of a spin model

Section titled “Competing exchange paths and the limits of a spin model”

Consider a bond with t=0.10 eVt=0.10\,\mathrm{eV} and U=4.0 eVU=4.0\,\mathrm{eV}. Its one-orbital superexchange estimate is

JSE=4t2U=0.010 eV=10 meV.J_{\mathrm{SE}}=\frac{4t^2}{U}=0.010\,\mathrm{eV}=10\,\mathrm{meV}.

A direct Coulomb exchange integral K=3 meVK=3\,\mathrm{meV} adds the leading contribution −2K=−6 meV-2K=-6\,\mathrm{meV}, giving J≃4 meVJ\simeq4\,\mathrm{meV}. This additive estimate assumes small charge admixture and neglects mixed higher-order corrections. The small net value reflects cancellation between larger paths. Modest structural changes can therefore reverse it even when neither microscopic process changes its own sign.

A fixed-length Heisenberg model works when charge and orbital excitations are well separated and the angular energy is approximately bilinear. Compare several spin orientations or observables: one collinear energy difference cannot reveal a biquadratic interaction or changing moment length. The review by Szilva and collaborators explains how exchange extraction depends on the reference electronic state and the retained variables.

If moment amplitudes change during rotation, orbital states remain active, or charge fluctuations are strong, include those degrees of freedom. Multispin, anisotropic, retarded, or state-dependent interactions can then matter. A useful effective model specifies the energy and temperature range in which its parameters predict the observations.

Let

∣χ(θ,ϕ)⟩=(cos⁡(θ/2)eiϕsin⁡(θ/2))|\chi(\theta,\phi)\rangle = \begin{pmatrix} \cos(\theta/2)\\ e^{i\phi}\sin(\theta/2) \end{pmatrix}

represent a spin aligned with a unit vector. Show that the magnitude of the overlap of two aligned spinors separated by angle Θ\Theta is cos⁡(Θ/2)\cos(\Theta/2).

Solution

For pure spin-1/21/2 states with Bloch vectors ni\mathbf n_i and nj\mathbf n_j,

∣⟨χi∣χj⟩∣2=1+ni⋅nj2.\left| \langle\chi_i|\chi_j\rangle \right|^2 = \frac{ 1+\mathbf n_i\cdot\mathbf n_j }{2}.

Since ni⋅nj=cos⁡Θ\mathbf n_i\cdot\mathbf n_j=\cos\Theta,

∣⟨χi∣χj⟩∣2=1+cos⁡Θ2=cos⁡2Θ2.\left| \langle\chi_i|\chi_j\rangle \right|^2 = \frac{1+\cos\Theta}{2} = \cos^2\frac{\Theta}{2}.

For 0≤Θ≤π0\le\Theta\le\pi, the magnitude is cos⁡(Θ/2)\cos(\Theta/2). Parallel core spins maximize hopping; antiparallel spins suppress it in the ideal strong-Hund limit.

2. The angular energy of one mobile electron

Section titled “2. The angular energy of one mobile electron”

For two classical core directions separated by θ\theta and one strongly Hund-aligned electron, the two-site hopping amplitude has magnitude tcos⁡(θ/2)t\cos(\theta/2), with t>0t>0. Find the electronic ground-state energy and its small-angle Heisenberg approximation Jfitℓ2cos⁡θJ_{\mathrm{fit}}\ell^2\cos\theta. Can this approximation reproduce the whole angular curve?

Solution

The two eigenvalues of the projected hopping matrix are ±tcos⁡(θ/2)\pm t\cos(\theta/2) for 0≤θ≤π0\le\theta\le\pi. The one-electron energy is

E(θ)=−tcos⁡θ2=−t+tθ28+O(θ4).E(\theta)=-t\cos\frac\theta2 =-t+\frac{t\theta^2}{8}+O(\theta^4).

Since Jfitℓ2cos⁡θJ_{\mathrm{fit}}\ell^2\cos\theta has quadratic term −Jfitℓ2θ2/2-J_{\mathrm{fit}}\ell^2\theta^2/2, matching the curvature gives Jfit=−t/(4ℓ2)J_{\mathrm{fit}}=-t/(4\ell^2). Adjusting the constant to match E(0)=−tE(0)=-t then predicts E(π)=−t/2E(\pi)=-t/2, whereas the exact projected value is zero. A good small-angle fit does not establish a bilinear interaction over large rotations.

3. Why midpoint inversion forbids Dzyaloshinskii–Moriya exchange

Section titled “3. Why midpoint inversion forbids Dzyaloshinskii–Moriya exchange”

Consider

HDM=Dij⋅(si×sj).H_{\mathrm{DM}} = \mathbf D_{ij}\cdot \left( \mathbf s_i\times\mathbf s_j \right).

Show why an inversion center at the bond midpoint requires Dij=0\mathbf D_{ij}=0.

Solution

Inversion about the midpoint exchanges sites ii and jj. Spin is an axial vector, so each spin direction is unchanged by spatial inversion, but the site labels swap:

si×sj⟶sj×si=−si×sj.\mathbf s_i\times\mathbf s_j \longrightarrow \mathbf s_j\times\mathbf s_i = - \mathbf s_i\times\mathbf s_j.

The bond is mapped onto itself. Invariance would therefore require the same coefficient multiplying both a term and its negative. The only possibility is Dij=0\mathbf D_{ij}=0.

A calculation uses two classical spins of length s=3/2s=3/2 and the convention H=Js1⋅s2H=J\mathbf s_1\cdot\mathbf s_2. It finds E∥−Eantiparallel=−18 meVE_{\parallel}-E_{\mathrm{antiparallel}}=-18\,\mathrm{meV} for an isolated pair. Determine JJ and its magnetic sign.

Solution

For parallel and antiparallel classical vectors,

E∥=Js2,Eantiparallel=−Js2.E_{\parallel} = Js^2, \qquad E_{\mathrm{antiparallel}} = -Js^2.

Thus

2Js2=−18 meV.2Js^2 = -18\,\mathrm{meV}.

With s2=9/4s^2=9/4,

J=−18 meV2(9/4)=−4 meV.J = \frac{-18\,\mathrm{meV}}{2(9/4)} = -4\,\mathrm{meV}.

The negative sign is ferromagnetic in this convention. In a lattice calculation, every changed bond and any moment-length change would also need to be counted.

The next step depends on the question. Ferromagnetism and Antiferromagnetism connect exchange to order; Spin Waves and Magnons in Materials connects an exchange network to measured excitation spectra. For two controlled spins, Silicon Spin Qubits uses exchange to generate gates and discusses charge-noise and leakage limits.

  • P. W. Anderson, “New Approach to the Theory of Superexchange Interactions,” Physical Review 115, 2–13 (1959), doi:10.1103/PhysRev.115.2.

  • P. W. Anderson and H. Hasegawa, “Considerations on Double Exchange,” Physical Review 100, 675–681 (1955), doi:10.1103/PhysRev.100.675.

  • J. B. Goodenough, “Theory of the Role of Covalence in the Perovskite-Type Manganites,” Physical Review 100, 564–573 (1955), doi:10.1103/PhysRev.100.564.

  • J. Kanamori, “Superexchange interaction and symmetry properties of electron orbitals,” Journal of Physics and Chemistry of Solids 10, 87–98 (1959), doi:10.1016/0022-3697(59)90061-7.

  • D. I. Khomskii, Transition Metal Compounds (Cambridge University Press, 2014), doi:10.1017/CBO9781139096782.

  • T. Moriya, “Anisotropic Superexchange Interaction and Weak Ferromagnetism,” Physical Review 120, 91–98 (1960), doi:10.1103/PhysRev.120.91.

  • M. A. Ruderman and C. Kittel, “Indirect Exchange Coupling of Nuclear Magnetic Moments by Conduction Electrons,” Physical Review 96, 99–102 (1954), doi:10.1103/PhysRev.96.99.

  • A. Szilva, Y. Kvashnin, E. A. Stepanov, L. Nordström, O. Eriksson, A. I. Lichtenstein, and M. I. Katsnelson, “Quantitative theory of magnetic interactions in solids,” Reviews of Modern Physics 95, 035004 (2023), doi:10.1103/RevModPhys.95.035004.