Spin–Orbit Coupling in Solids

Spin–orbit coupling in a solid is the relativistic coupling that entangles an electron’s spin with its orbital motion through the crystal potential. It begins with the same local physics that produces atomic fine structure, but a crystal adds orbital hybridization, crystal fields, Bloch momentum, surfaces, interfaces, and space-group symmetry. Those ingredients decide whether spin–orbit coupling merely shifts local levels, splits a band, changes magnetic anisotropy, generates a transverse spin response, or reorganizes the topology of the occupied states.

The phrase “strong spin–orbit coupling” is incomplete without a comparison scale. A matrix element $\lambda$ can be large relative to a magnetic exchange energy and small relative to a crystal-field gap; it can strongly alter a nearly degenerate band crossing while barely changing remote bands. The useful question is therefore not whether spin–orbit coupling is present—it is—but which projected Hilbert space it acts in, which symmetries constrain it, and which competing energy scales it exceeds.

Spin–Orbit Coupling is the canonical home for the angular-momentum-addition algebra of $\mathbf L\cdot\mathbf S$. Fine Structure develops its atomic spectroscopic role, and the Pauli Equation gives its nonrelativistic origin. This article starts from those results and develops their consequences in crystalline materials.

Required background. Spin–Orbit Coupling supplies angular-momentum algebra, while Symmetry of Bloch States supplies little-group, double-group, and degeneracy labels.

Helpful background. Tight-Binding Models supplies basis projections, Time Reversal for Spin-1/2 supplies Kramers constraints, and Berry Curvature supplies the geometric response object.

Convention and scale ledger

We use the following conventions throughout:

• $m_0>0$ is the free-electron mass, while $m^\ast$ is a signed band effective mass when a one-band approximation is valid.
• $V(\mathbf r)$ is the electron’s potential energy, not the electrostatic potential.
• $\boldsymbol{\sigma}$ denotes Pauli matrices in an explicitly stated two-state subspace. It represents physical spin only when that projection is justified.
• $\mathbf k$ is crystal momentum. A continuum expansion is local to a specified high-symmetry point and is not automatically periodic over the Brillouin zone.
• Time reversal for a physical spin-$1/2$ doublet is $\mathcal T=i\sigma_yK$, up to an overall phase, and $\mathcal T^2=-1$.
• Inversion sends $\mathbf k\mapsto-\mathbf k$ but leaves spin unchanged because spin is an axial vector.
• The signs of Rashba and Dresselhaus coefficients depend on axis, basis, interface-normal, and Pauli-matrix conventions. Spectra determine magnitudes more directly than signs.

Several energy scales must be kept separate:

$$\lambda_{\mathrm{at}}, \qquad \Delta_{\mathrm{CF}}, \qquad t, \qquad \Delta_{\mathrm{ex}}, \qquad \Delta_{\mathrm{SO}}(\mathbf k), \qquad \hbar/\tau .$$

Here $\lambda_{\mathrm{at}}$ is an atomic or local spin–orbit scale, $\Delta_{\mathrm{CF}}$ a crystal-field splitting, $t$ a hopping or bandwidth scale, $\Delta_{\mathrm{ex}}$ an exchange splitting, $\Delta_{\mathrm{SO}}(\mathbf k)$ a band splitting or avoided-crossing gap, and $\hbar/\tau$ a disorder broadening. There is no universal ordering among them.

Relativistic origin

For a static scalar potential, the Foldy–Wouthuysen expansion of the Dirac Hamiltonian contains the Pauli spin–orbit term

$$H_{\mathrm{SO}} = \frac{\hbar}{4m_0^2c^2} \boldsymbol{\sigma}\cdot \left( \boldsymbol{\nabla}V \times \mathbf p \right).$$

The coefficient includes the Thomas-precession factor. If one instead writes an electrostatic potential $\phi$ for an electron of charge $-e$, then $V=-e\phi$; silently replacing $V$ by $\phi$ reverses a sign.

For a central potential $V(r)$,

$$\boldsymbol{\nabla}V \times \mathbf p = \frac{1}{r} \frac{\mathrm dV}{\mathrm dr} \mathbf L ,$$

so that

$$H_{\mathrm{SO}} = \frac{1}{2m_0^2c^2r} \frac{\mathrm dV}{\mathrm dr} \mathbf L\cdot\mathbf S .$$

The large gradient near a heavy nucleus is why atomic spin–orbit scales generally increase strongly with atomic number. In a solid, however, a band state is a coherent mixture of atomic orbitals. The relevant matrix element depends on which atomic regions it samples, the orbital characters and covalency, and the projection used to construct a low-energy model. “Made from a heavy element” is useful chemical guidance, not a quantitative answer.

For a local shell one often writes

$$H_{\mathrm{loc}} = H_{\mathrm{CF}} + \lambda\,\mathbf L\cdot\mathbf S + H_{\mathrm{cov}} + H_{\mathrm{int}},$$

where $H_{\mathrm{CF}}$ is the crystal field, $H_{\mathrm{cov}}$ represents hybridization with ligands, and $H_{\mathrm{int}}$ contains local interactions. The fitted $\lambda$ in this Hamiltonian is an effective matrix element; it need not equal a free-ion constant.

Crystal fields, orbital quenching, and pseudospin

A crystal field removes some or all of an atomic multiplet’s orbital degeneracy. If an orbital eigenstate is nondegenerate and can be chosen real, its first-order expectation value of $\mathbf L$ vanishes. This is the usual sense in which the orbital moment is quenched. It does not mean that $\mathbf L$ has disappeared as an operator. Virtual admixture of excited crystal-field levels produces corrections of order

$$\delta E_{\mathrm{SO}} \sim \frac{\lambda^2}{\Delta_{\mathrm{CF}}},$$

as well as anisotropic $g$ factors and magnetocrystalline anisotropy. Near an orbital degeneracy, spin–orbit coupling can instead act in first order.

Cubic d orbitals

In an ideal octahedral environment, the five $d$ orbitals split into a lower $t_{2g}$ triplet and an upper $e_g$ doublet for the usual ligand-field sign. Within the $t_{2g}$ manifold, projected orbital matrices can be represented by an effective angular momentum $\ell_{\mathrm{eff}}=1$ with the opposite sign:

$$P_{t_{2g}}\,\mathbf L\,P_{t_{2g}} = -\mathbf l_{\mathrm{eff}}.$$

Consequently,

$$P_{t_{2g}}H_{\mathrm{SO}}P_{t_{2g}} = -\lambda\, \mathbf l_{\mathrm{eff}}\cdot\mathbf S .$$

The statement is about matrices in a specified orbital basis, not a claim that the electron’s microscopic orbital angular momentum has literally changed sign. Electron versus hole descriptions, distortions, covalency, and the sign convention for $\lambda$ must all be declared before assigning level order.

In the ideal $e_g$ subspace, the direct projection $P_{e_g}\mathbf LP_{e_g}$ vanishes. Spin–orbit effects nevertheless survive through mixing with $t_{2g}$ states, ligand orbitals, and lower-symmetry fields. Thus “orbital quenching” and “negligible spin–orbit coupling” are not synonyms.

Three useful coupling regimes

The order in which one diagonalizes local terms should follow the scale hierarchy:

1. If $\Delta_{\mathrm{CF}}\gg|\lambda|$, first diagonalize the crystal field and project $\mathbf L\cdot\mathbf S$ into its low-energy manifold.
2. If $|\lambda|\gg\Delta_{\mathrm{CF}}$, first form atomic $j$ multiplets and then apply the crystal field within them.
3. If the scales are comparable, neither limiting label is exact; diagonalize the full local Hamiltonian and report wave-function composition rather than forcing a pure $LS$- or $jj$-coupling name.

The low-energy two-state object is often a Kramers pseudospin rather than a bare spin. Its magnetic moment has the general form

$$\boldsymbol{\mu} = -\mu_{\mathrm B}\, \mathbf g\cdot \frac{\boldsymbol{\sigma}}{2},$$

where $\mathbf g$ is a tensor. Neutron scattering, magnetic resonance, and Zeeman coupling probe matrix elements of the physical moment, not the abstract Pauli labels alone.

Three distinct uses of spin–orbit coupling. Left: a crystal field first organizes local orbitals, after which $\lambda$ acts within and between the surviving manifolds. Center: inversion asymmetry permits momentum-dependent Rashba or Dresselhaus fields. Right: spin–orbit coupling can gap a symmetry-allowed crossing and concentrate Berry curvature, but the existence of a gap alone does not determine a topological invariant.

Bloch Hamiltonians and symmetry

Near a chosen point, an isolated two-state band subspace can be written

$$H(\mathbf k) = \varepsilon_0(\mathbf k)I + \mathbf g(\mathbf k)\cdot\boldsymbol{\sigma}.$$

Its eigenvalues are

$$E_\pm(\mathbf k) = \varepsilon_0(\mathbf k) \pm |\mathbf g(\mathbf k)|.$$

For $\boldsymbol{\sigma}$ representing physical spin, the two eigenstates have

$$\langle\boldsymbol{\sigma}\rangle_\pm = \pm \frac{\mathbf g}{|\mathbf g|}.$$

For a spin–orbital pseudospin, this equation describes the pseudospin texture; the physical spin texture follows only after projecting the spin operator into the same subspace.

Time reversal and inversion

Time-reversal invariance requires

$$\mathcal T H(\mathbf k) \mathcal T^{-1} = H(-\mathbf k),$$

which gives

$$\varepsilon_0(\mathbf k) = \varepsilon_0(-\mathbf k), \qquad \mathbf g(\mathbf k) = -\mathbf g(-\mathbf k).$$

At a time-reversal-invariant momentum $\boldsymbol{\Lambda}$ satisfying

$$\boldsymbol{\Lambda} = -\boldsymbol{\Lambda} + \mathbf G,$$

with reciprocal vector $\mathbf G$, the odd field vanishes and the pair is Kramers degenerate.

If inversion is also present, then

$$\mathcal P H(\mathbf k) \mathcal P^{-1} = H(-\mathbf k).$$

Within a same-parity spin doublet this makes $\mathbf g$ even. Time reversal makes it odd, so $\mathbf g=0$. More generally, the antiunitary operation $\mathcal P\mathcal T$ leaves every $\mathbf k$ fixed and satisfies

$$(\mathcal P\mathcal T)^2=-1$$

for spinful electrons when inversion and time reversal commute. Every Bloch level is then at least twofold degenerate at every $\mathbf k$, not only at the time-reversal-invariant momenta. Kramers Degeneracy develops the theorem and its hypotheses.

Breaking inversion while preserving time reversal allows spin-split bands away from the invariant momenta. Breaking time reversal allows a Zeeman or exchange term even at those momenta. A measured band splitting therefore does not identify Rashba coupling until the relevant crystal, surface, magnetic, and little-group symmetries have been established.

Rashba coupling

For an approximately isotropic two-dimensional band with a polar axis $\hat{\mathbf z}$, the leading structural-inversion-asymmetry term is

$$H_{\mathrm R} = \alpha_{\mathrm R} \left( k_y\sigma_x-k_x\sigma_y \right) = \alpha_{\mathrm R} \left( \mathbf k\times\hat{\mathbf z} \right)\cdot \boldsymbol{\sigma}.$$

This convention gives

$$\mathbf g_{\mathrm R} = \alpha_{\mathrm R} \left( k_y,-k_x,0 \right).$$

The sign changes if the surface normal, coordinate orientation, or Pauli basis is reversed. For a parabolic band,

$$H = \frac{\hbar^2k^2}{2m^\ast}I + H_{\mathrm R},$$

and

$$E_\pm(k) = \frac{\hbar^2k^2}{2m^\ast} \pm |\alpha_{\mathrm R}|k .$$

For $m^\ast>0$, useful scales are

$$k_{\mathrm R} = \frac{m^\ast|\alpha_{\mathrm R}|}{\hbar^2}, \qquad E_{\mathrm R} = \frac{m^\ast\alpha_{\mathrm R}^2}{2\hbar^2}.$$

The lower branch has a ring minimum at $k=k_{\mathrm R}$ lying $E_{\mathrm R}$ below the unsplit crossing energy in this ideal model. The two pseudospins wind tangentially in opposite directions around the crossing.

It is tempting to identify $\alpha_{\mathrm R}$ with a bare electric field in the vacuum Pauli term. In actual solids, $\alpha_{\mathrm R}$ is a band parameter produced by inversion-asymmetric orbital hybridization together with atomic spin–orbit matrix elements. Wave-function weight near heavy atoms, avoided crossings, confinement, electrostatic gating, and structural relaxation can all matter. A large potential gradient without suitable orbital mixing need not produce a large band splitting.

Crystal symmetry also permits higher harmonics. Hexagonal warping, for example, can generate an out-of-plane spin component and a noncircular contour. The linear isotropic Hamiltonian is a controlled low-$k$ model, not a universal definition of every surface splitting.

Dresselhaus coupling

Bulk inversion asymmetry in a zinc-blende crystal permits the cubic Dresselhaus invariant

$$\begin{aligned} H_{\mathrm D}^{(3\mathrm D)} =\gamma\big[ &k_x(k_y^2-k_z^2)\sigma_x \\ &+k_y(k_z^2-k_x^2)\sigma_y \\ &+k_z(k_x^2-k_y^2)\sigma_z \big]. \end{aligned}$$

The angular dependence follows the crystal axes and is qualitatively different from the rotationally symmetric Rashba field. Along selected high-symmetry directions, the leading splitting vanishes.

Confinement in a [001] quantum well replaces $k_z^2$ by a subband expectation value and produces a linear in-plane term. Absorbing projection signs into $\beta$, we write

$$H_{\mathrm D}^{(1)} = \beta \left( k_x\sigma_x-k_y\sigma_y \right),$$

while residual cubic terms become important as the in-plane Fermi momentum grows. The fitted $\beta$ depends on well width, strain, composition, subband occupation, and convention; it is not simply a universal bulk $\gamma$ multiplied by one geometric number.

Rashba and Dresselhaus together

For a [001] two-dimensional system with both linear terms,

$$\begin{aligned} H_{\mathrm{RD}} = & \alpha \left( k_y\sigma_x-k_x\sigma_y \right) \\ &+ \beta \left( k_x\sigma_x-k_y\sigma_y \right), \end{aligned}$$

so

$$\mathbf g_{\mathrm{RD}} = \left( \alpha k_y+\beta k_x,\, -\alpha k_x-\beta k_y,\, 0 \right).$$

The splitting is anisotropic:

$$|\mathbf g_{\mathrm{RD}}|^2 = (\alpha^2+\beta^2) (k_x^2+k_y^2) + 4\alpha\beta k_xk_y .$$

At $\alpha=\beta$, the field points along the fixed spin axis $[1\bar 10]$:

$$H_{\mathrm{RD}} = 2\alpha k_+\sigma_-,$$

where

$$k_+ = \frac{k_x+k_y}{\sqrt 2}, \qquad \sigma_- = \frac{\sigma_x-\sigma_y}{\sqrt 2}.$$

The component $\sigma_-$ is conserved in the ideal parabolic, linear-coupling model. The associated exact SU(2) structure supports a persistent spin helix with wavevector magnitude

$$Q = \frac{4m^\ast|\alpha|}{\hbar^2}.$$

Cubic Dresselhaus terms, unequal coefficients, spatial disorder, additional subbands, and other spin-relaxation channels limit the lifetime in real samples. The symmetry point is robust physics, not a promise of infinite experimental lifetime.

Magnetic consequences

Spin–orbit coupling lets the lattice distinguish spin directions. After high-energy orbital states are eliminated, it can generate:

• anisotropic $g$ tensors and single-ion anisotropy;
• symmetric anisotropic exchange and antisymmetric Dzyaloshinskii–Moriya exchange;
• magnetocrystalline anisotropy and easy axes or easy planes;
• spin–momentum locking and current-induced nonequilibrium spin polarization;
• gaps at magnon or electronic band crossings;
• orbital moments and mixed spin–orbital collective modes.

Exchange Interactions gives the canonical exchange-tensor decomposition. Magnetic Anisotropy develops the resulting easy directions, shape and interface competition, texture length scales, and material inference. The important power counting is that a local spin–orbit term can enter an effective spin Hamiltonian only after it is combined with hopping, hybridization, exchange, or crystal-field processes. A Dzyaloshinskii–Moriya vector is therefore constrained by bond symmetry and microscopic virtual paths; it is not obtained by attaching an arbitrary vector to every bond.

In an itinerant ferromagnet, exchange already breaks time reversal and separates spin-mixed bands. Spin–orbit coupling then changes their avoided crossings, orbital moments, Berry curvature, and orientation-dependent total energy. Describing such states as pure “spin up” and “spin down” is generally an approximation.

Spin Hall response

The spin Hall effect is a transverse spin response to a longitudinal electric field. A commonly used conventional spin-current operator is

$$J_i^a = \frac{1}{2} \left\{ v_i,s_a \right\},$$

where $i$ is the flow direction and $a$ is the spin-polarization axis. Because spin–orbit coupling produces a torque,

$$\frac{\mathrm ds_a}{\mathrm dt} = \frac{i}{\hbar} \left[ H,s_a \right] \ne 0,$$

this current is generally not associated with a strictly conserved bulk density. Alternative conserved-current definitions shift torque-dipole terms between “current” and “source.” A quoted spin Hall conductivity is meaningful only with the operator, axes, units, frequency, and boundary protocol stated.

Intrinsic contribution

For a clean multiband crystal, an intrinsic dc contribution can be written in Kubo form as

$$\sigma_{ij}^{a,\mathrm{int}} = e\hbar \sum_n \int_{\mathrm{BZ}} \frac{\mathrm d^dk}{(2\pi)^d} f_{n\mathbf k}\, \Omega_{n,ij}^{a}(\mathbf k),$$

with the spin Berry curvature

$$\Omega_{n,ij}^{a} = -2\,\operatorname{Im} \sum_{m\ne n} \frac{ \langle n|J_i^a|m\rangle \langle m|v_j|n\rangle }{ (\varepsilon_n-\varepsilon_m)^2 }.$$

All states and energies here are evaluated at the same $\mathbf k$. Near spin–orbit-induced avoided crossings, the small denominator can produce large hot spots. The total response still requires the occupation, all relevant bands, and symmetry-related cancellations.

Time reversal forbids an ordinary intrinsic charge Hall conductivity in a nonmagnetic equilibrium crystal because charge Berry curvature is odd in $\mathbf k$. It does not generally forbid spin Hall response: both velocity and spin reverse under time reversal, so their symmetrized product is even.

Extrinsic contribution

Disorder can generate two distinct leading mechanisms:

• skew scattering, an asymmetric scattering probability whose conductivity often scales with the transport lifetime in a dilute-disorder regime;
• side jump, a coordinate displacement during scattering whose leading conductivity is often independent of the lifetime, even though scattering events are essential.

The decomposition into intrinsic, side-jump, and skew terms can depend on model and representation, especially in multiband systems. Scaling with longitudinal conductivity, impurity control, thickness dependence, and comparison with microscopic calculations are more informative than assigning a mechanism from one transverse voltage.

A widely used conversion is

$$\theta_{\mathrm{SH}} = \frac{2e}{\hbar} \frac{j_s}{j_c},$$

but some authors absorb $\hbar/2e$ into the definition of $j_s$ or reverse an axis. A spin Hall angle without its convention cannot be compared safely across papers.

The ideal linear Rashba two-dimensional electron gas is a useful warning. Its clean-band calculation suggests a simple intrinsic dc value, but vertex corrections cancel that result for a common model of spin-independent short-range disorder in the uniform, zero-frequency limit. The cancellation is not a theorem for real multiorbital crystals, finite frequency, finite size, nonlinear coupling, or arbitrary disorder.

Hall Effect owns the general conductivity-tensor, sign, multiband, and experimental-inference framework. Spintronics owns spin injection, inverse conversion, torques, and device architectures rather than the definition of spin–orbit coupling itself.

Bridge to topological materials

Spin–orbit coupling can gap a band crossing, invert orbital characters, or change the symmetry representations of occupied states. None of those observations alone proves topology. A topological phase is diagnosed by the global occupied-state bundle and the protecting symmetries.

In a two-band avoided crossing,

$$H(\mathbf q) = v_xq_x\sigma_x + v_yq_y\sigma_y + m_{\mathrm{SO}}\sigma_z,$$

the gap is $2|m_{\mathrm{SO}}|$, and the Berry curvature is concentrated near $\mathbf q=0$. Whether the integrated contribution is topologically nontrivial depends on how this valley is completed elsewhere in the Brillouin zone and on the signs of all symmetry-related masses.

For a time-reversal-invariant two-dimensional insulator, Kramers structure permits a $\mathbb Z_2$ classification even when no spin component is conserved. In the spin-conserving limit of the Kane–Mele model, one can picture two opposite-Chern sectors related by time reversal. Rashba coupling mixes those sectors, yet the quantum spin Hall phase can survive as long as time reversal and the bulk gap remain intact. The spin Chern-number picture is then no longer the fundamental definition.

Berry Curvature develops the local geometric response, Chern Numbers develops the integer invariant, and Topological Quantum Numbers compares invariant families. Topological Insulators owns the full class-AII $\mathbb Z_2$ construction and material evidence ledger. A material claim additionally needs a bulk invariant, a gap or mobility-gap statement, symmetry control, and boundary evidence consistent with the same phase.

How to infer spin–orbit physics in a material

No single measurement returns “the spin–orbit coupling.” A defensible analysis separates model identification, parameter extraction, and mechanism attribution.

Symmetry first

1. Determine the bulk space group and whether the relevant surface, interface, distortion, magnetic order, or applied field breaks inversion or time reversal.
2. Identify the little group along the measured momentum path and the degeneracies it protects.
3. Decide whether a two-band Pauli model is legitimate or whether orbital and sublattice degrees of freedom must remain explicit.
4. Enumerate the lowest-order invariants before fitting their coefficients.

Symmetry of Bloch States owns the little-group, double-group, and nonsymmorphic labels and the degeneracy or crossing taxonomy used in the first two steps. This page owns the microscopic spin–orbit mechanisms and the resulting Rashba, Dresselhaus, spin-Hall, and material-facing consequences.

Bulk inversion asymmetry, structural inversion asymmetry, and exchange splitting can produce superficially similar band separations but obey different angular and symmetry constraints.

Spectroscopy and transport

Angle-resolved photoemission can measure dispersions and, with spin resolution, spin expectation values. Matrix elements, overlapping bands, final-state effects, surface band bending, and finite resolution can mimic or obscure a splitting. Quantum oscillations resolve extremal orbit areas but require care: two frequencies can arise from warping, multiple pockets, domains, or magnetic breakdown rather than spin splitting.

Weak antilocalization, spin resonance, optical orientation, and spin-relaxation anisotropy constrain spin precession and decoherence. Their extracted “spin–orbit field” is model dependent. Fits should state the diffusive or ballistic regime, dimensionality, intervalley and magnetic scattering assumptions, and the hierarchy among $\tau$, the spin-precession time, and the phase-coherence time.

Spin Hall measurements often infer a bulk conversion through edge spin accumulation, inverse spin Hall voltages, spin-torque ferromagnetic resonance, nonlocal transport, or terahertz emission. Interface transparency, current shunting, spin diffusion, thermoelectric backgrounds, and spin-memory loss must be calibrated before assigning a bulk angle.

Computation

A useful first-principles workflow is:

1. compare scalar-relativistic and fully relativistic bands using the same converged structure;
2. project wave functions onto atomic and symmetry-adapted orbitals;
3. verify symmetry-enforced degeneracies and transformation properties;
4. construct a spinor Wannier or other low-energy model only after checking the disentanglement window;
5. converge avoided-crossing gaps and spin Berry-curvature hot spots on dense momentum meshes;
6. test the stability of fitted Rashba or Dresselhaus coefficients over the claimed momentum range.

Turning spin–orbit coupling on and off in a code is a diagnostic counterfactual, not an experimentally tunable operation. It can reveal which feature is SOC-enabled, but it does not by itself identify the atomic path or establish a topological invariant.

Common mistakes

Mistake	Why it fails	Better check
Calling every nonmagnetic band splitting “Rashba”	Bulk inversion asymmetry, nonsymmorphic symmetry, hidden sector polarization, or unresolved bands may be responsible	Establish bulk and local symmetry, momentum dependence, and spin texture
Treating Pauli matrices as physical spin	A low-energy doublet may be strongly spin–orbital entangled	Project the physical spin and magnetic-moment operators
Equating a heavy atom with a large target-band splitting	Orbital weight, hybridization, and energy denominators control the projection	Inspect orbital-resolved wave functions and avoided crossings
Inferring topology from an SOC-opened gap	A local mass does not determine a global invariant	Compute the appropriate bulk invariant and check its protecting symmetry
Comparing spin Hall angles without conventions	Current normalization and axis signs differ	Report $J_i^a$, geometry, units, and the $\hbar/2e$ conversion
Assuming inversion plus time reversal forbids SOC	It forbids spin splitting of each Bloch level, not spin–orbital entanglement	Distinguish degeneracy from absence of SOC matrix elements
Dropping cubic Dresselhaus terms at all densities	Their relative importance grows with in-plane momentum	Fit angular dependence over the experimental momentum range
Reading a fitted $\alpha_{\mathrm R}$ as a microscopic electric field	The coefficient is a multiband material parameter	Track structural, orbital, and gate dependence in a microscopic model

Exercises

1. Recover the central-potential form

Starting from

$$H_{\mathrm{SO}} = \frac{\hbar}{4m_0^2c^2} \boldsymbol{\sigma}\cdot \left( \boldsymbol{\nabla}V\times\mathbf p \right),$$

derive the $\mathbf L\cdot\mathbf S$ form for $V=V(r)$. State where the factor of $1/2$ relative to a naive rest-frame argument appears.

Solution

For a central potential,

$$\boldsymbol{\nabla}V = \frac{\mathrm dV}{\mathrm dr}\hat{\mathbf r}, \qquad \hat{\mathbf r}\times\mathbf p = \frac{\mathbf r\times\mathbf p}{r} = \frac{\mathbf L}{r}.$$

Using $\mathbf S=\hbar\boldsymbol{\sigma}/2$ gives

$$\begin{aligned} H_{\mathrm{SO}} &= \frac{\hbar}{4m_0^2c^2} \frac{1}{r} \frac{\mathrm dV}{\mathrm dr} \boldsymbol{\sigma}\cdot\mathbf L \\ &= \frac{1}{2m_0^2c^2r} \frac{\mathrm dV}{\mathrm dr} \mathbf L\cdot\mathbf S. \end{aligned}$$

The Pauli coefficient already contains the Thomas-precession correction. A naive transformation to the electron’s instantaneous rest frame misses this kinematic factor and gives twice the correct spin–orbit coupling.

2. Degeneracy with inversion and time reversal

For a spinful, nonmagnetic centrosymmetric crystal, show that every Bloch energy is at least doubly degenerate at each $\mathbf k$. Why does this not imply that the eigenstates are pure spin states?

Solution

Inversion maps $\mathbf k$ to $-\mathbf k$, and time reversal does the same. Their product $\mathcal P\mathcal T$ therefore leaves $\mathbf k$ fixed. For ordinary inversion, $\mathcal P^2=1$ and $[\mathcal P,\mathcal T]=0$, while $\mathcal T^2=-1$. Hence

$$(\mathcal P\mathcal T)^2 = \mathcal P^2\mathcal T^2 = -1.$$

If $|u_{n\mathbf k}\rangle$ is an eigenstate, $\mathcal P\mathcal T|u_{n\mathbf k}\rangle$ is an orthogonal eigenstate at the same energy and momentum. Spin–orbit coupling may still entangle spin with orbital and sublattice amplitudes within this doublet. Degeneracy constrains the energy; it does not provide a globally unique spin quantization axis.

3. Rashba momentum and energy scales

For $m^\ast>0$ and

$$E_-(k) = \frac{\hbar^2k^2}{2m^\ast} -|\alpha_{\mathrm R}|k,$$

find the ring minimum and its energy relative to $E_-(0)$. Give the spin direction for $\alpha_{\mathrm R}>0$ on the $+k_x$ axis.

Solution

Differentiation gives

$$\frac{\mathrm dE_-}{\mathrm dk} = \frac{\hbar^2k}{m^\ast} -|\alpha_{\mathrm R}| = 0,$$

so

$$k_{\mathrm R} = \frac{m^\ast|\alpha_{\mathrm R}|}{\hbar^2}.$$

Substitution yields

$$E_-(k_{\mathrm R})-E_-(0) = -\frac{m^\ast\alpha_{\mathrm R}^2}{2\hbar^2} = -E_{\mathrm R}.$$

At $\mathbf k=(k,0)$,

$$\mathbf g_{\mathrm R} = (0,-\alpha_{\mathrm R}k,0).$$

The $+$ helicity has spin along $-\hat{\mathbf y}$ and the $-$ helicity along $+\hat{\mathbf y}$. Which energy branch is called “inner” or “outer” depends on the chosen Fermi energy.

4. Persistent-spin-helix symmetry

Set $\alpha=\beta$ in the linear Rashba–Dresselhaus Hamiltonian. Show that one spin component is conserved and derive the helix wavevector for a parabolic band.

Solution

At equal couplings,

$$\begin{aligned} H_{\mathrm{RD}} &= \alpha(k_x+k_y)(\sigma_x-\sigma_y) \\ &= 2\alpha k_+\sigma_-. \end{aligned}$$

Therefore $[H,\sigma_-]=0$ for

$$\sigma_- = \frac{\sigma_x-\sigma_y}{\sqrt 2}.$$

The kinetic and SOC terms along $k_+$ complete the square:

$$\frac{\hbar^2k_+^2}{2m^\ast} + 2\alpha k_+\sigma_- = \frac{\hbar^2}{2m^\ast} \left( k_+ + \frac{2m^\ast\alpha}{\hbar^2}\sigma_- \right)^2 - \frac{2m^\ast\alpha^2}{\hbar^2}.$$

The two conserved-spin branches are shifted by

$$Q = \frac{4m^\ast|\alpha|}{\hbar^2},$$

which is the spatial helix wavevector. Terms that fail to commute with $\sigma_-$ make the lifetime finite.

5. Dresselhaus nodes

Use the cubic bulk Dresselhaus Hamiltonian to show that the spin splitting vanishes for momenta along $[100]$ and $[111]$. Is this enough to conclude that the entire line remains degenerate in a real material?

Solution

For $[100]$, take $(k_x,k_y,k_z)=(k,0,0)$. Every component contains either a vanishing prefactor or a vanishing difference, so $\mathbf g_{\mathrm D}=0$.

For $[111]$, take $k_x=k_y=k_z=k/\sqrt 3$. Each difference such as $k_y^2-k_z^2$ vanishes, again giving $\mathbf g_{\mathrm D}=0$.

This conclusion applies to the displayed lowest-order invariant in the chosen band and symmetry setting. Higher-order invariants, coupling to other bands, strain, surfaces, magnetic order, or lower actual symmetry can lift the degeneracy. The little-group representations of the full crystal decide what is exact.

6. Why time reversal permits a spin Hall response

Determine the time-reversal parity of the charge current $j_i=-ev_i$ and of the conventional spin current $J_i^a=\{v_i,s_a\}/2$. Use it to explain why a nonmagnetic crystal can have zero charge Hall conductivity and nonzero spin Hall conductivity.

Solution

Time reversal changes both velocity and spin:

$$\mathcal T v_i\mathcal T^{-1} = -v_i, \qquad \mathcal T s_a\mathcal T^{-1} = -s_a.$$

Thus charge current is odd, whereas

$$\mathcal T J_i^a\mathcal T^{-1} = \frac{1}{2} \left\{ -v_i,-s_a \right\} = J_i^a$$

is even. In a time-reversal-invariant equilibrium crystal, symmetry-related momenta cancel the transverse charge response. The extra time-reversal-odd spin insertion reverses that cancellation pattern, so a transverse spin current is symmetry allowed. Its actual value can still vanish because of additional crystal symmetries or dynamical cancellations.

7. A gap is not yet a topological invariant

Two time-reversal-related massive Dirac valleys contribute Berry-curvature hot spots. Explain why observing the same gap magnitude at both valleys does not determine whether their topological contributions add or cancel.

Solution

For a local Dirac Hamiltonian, the sign of the Berry-curvature contribution depends not only on $|m|$ but on the sign of the mass and the orientation encoded by the velocity Jacobian. Schematically,

$$C_{\mathrm{valley}} \propto \operatorname{sgn} \left( m\,v_xv_y \right).$$

Spectroscopy of the gap gives $2|m|$, not this signed product. Symmetry can relate the masses and valley orientations so that contributions cancel, add within one spin sector, or reorganize into a $\mathbb Z_2$ invariant. One must determine the global occupied wave functions or an equivalent symmetry indicator; equal local gap magnitudes are insufficient.

Connections

• Unconventional Superconductivity takes a spin–orbit-projected pseudospin or orbital basis as input for pairing classification; this page retains the material Hamiltonian, band spin texture, and physical-moment projection.
• Magnetic Moments in Matter takes the crystal-field and spin–orbit-projected manifold as input, interprets its physical $g$ tensor, and tests whether a fixed pseudospin moment is controlled.
• Spin–Orbit Coupling owns the coupled-basis spectrum of $\lambda\mathbf L\cdot\mathbf S$.
• Pauli Equation derives the nonrelativistic relativistic corrections and fixes electromagnetic conventions.
• Tight-Binding Models develops orbital Bloch Hamiltonians into which local and bond spin–orbit matrices are projected.
• Time Reversal for Spin One-Half develops $\mathcal T=i\sigma_yK$ and $\mathcal T^2=-1$.
• Fermi Surface explains how band splitting changes contours, extremal orbits, and momentum-resolved observables.
• Exchange Interactions develops symmetric anisotropy and Dzyaloshinskii–Moriya exchange.
• Hall Effect supplies the charge-transport baseline needed to interpret spin Hall conversion.
• Berry Curvature connects interband matrix elements and avoided crossings to geometric response.
• Topological Quantum Numbers distinguishes local curvature from global invariants.
• Weyl and Dirac Semimetals shows how spin–orbit coupling, inversion breaking, magnetism, and crystalline symmetry organize topological point crossings in solids.
• Quantum Wells supplies the confinement, heterointerface, subband, and electric-field ledger for Rashba- and Dresselhaus-coupled two-dimensional carriers.
• Two-Dimensional Electron Gases connects those couplings to gate-dependent occupation, weak-antilocalization and oscillation diagnostics, disorder lifetimes, and oxide interfaces.
• Transition-Metal Dichalcogenides applies atomic spin–orbit coupling and broken inversion to a material-specific spin–valley and exciton hierarchy.
• 2D Magnets and Ferroelectrics connects spin–orbit-induced anisotropy and valley locking to finite-temperature monolayer order, polarization, and magnetoelectric coupling.
• van der Waals Heterostructures develops the interface self-energy and evidence ledger for proximity-induced spin–orbit terms.

References

1. L. L. Foldy and S. A. Wouthuysen, “On the Dirac Theory of Spin 1/2 Particles and Its Non-Relativistic Limit,” Physical Review 78, 29–36 (1950), doi:10.1103/PhysRev.78.29.
2. J. H. Van Vleck, The Theory of Electric and Magnetic Susceptibilities, Oxford University Press (1932).
3. A. Abragam and B. Bleaney, Electron Paramagnetic Resonance of Transition Ions, Oxford University Press (1970).
4. G. Khaliullin, “Orbital Order and Fluctuations in Mott Insulators,” Progress of Theoretical Physics Supplement 160, 155–202 (2005), doi:10.1143/PTPS.160.155.
5. E. I. Rashba, “Properties of Semiconductors with an Extremum Loop. 1. Cyclotron and Combinational Resonance in a Magnetic Field Perpendicular to the Plane of the Loop,” Soviet Physics, Solid State 2, 1109–1122 (1960).
6. Y. A. Bychkov and E. I. Rashba, “Properties of a 2D Electron Gas with Lifted Spectral Degeneracy,” JETP Letters 39, 78–81 (1984).
7. G. Dresselhaus, “Spin-Orbit Coupling Effects in Zinc Blende Structures,” Physical Review 100, 580–586 (1955), doi:10.1103/PhysRev.100.580.
8. R. Winkler, Spin–Orbit Coupling Effects in Two-Dimensional Electron and Hole Systems, Springer (2003), doi:10.1007/b13586.
9. A. Manchon, H. C. Koo, J. Nitta, S. M. Frolov, and R. A. Duine, “New Perspectives for Rashba Spin–Orbit Coupling,” Nature Materials 14, 871–882 (2015), doi:10.1038/nmat4360.
10. B. A. Bernevig, J. Orenstein, and S.-C. Zhang, “Exact SU(2) Symmetry and Persistent Spin Helix in a Spin-Orbit Coupled System,” Physical Review Letters 97, 236601 (2006), doi:10.1103/PhysRevLett.97.236601.
11. J. D. Koralek, C. P. Weber, J. Orenstein, B. A. Bernevig, S.-C. Zhang, S. Mack, and D. D. Awschalom, “Emergence of the Persistent Spin Helix in Semiconductor Quantum Wells,” Nature 458, 610–613 (2009), doi:10.1038/nature07871.
12. M. I. Dyakonov and V. I. Perel, “Possibility of Orienting Electron Spins with Current,” JETP Letters 13, 467–469 (1971).
13. J. E. Hirsch, “Spin Hall Effect,” Physical Review Letters 83, 1834–1837 (1999), doi:10.1103/PhysRevLett.83.1834.
14. S. Murakami, N. Nagaosa, and S.-C. Zhang, “Dissipationless Quantum Spin Current at Room Temperature,” Science 301, 1348–1351 (2003), doi:10.1126/science.1087128.
15. J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jungwirth, and A. H. MacDonald, “Universal Intrinsic Spin Hall Effect,” Physical Review Letters 92, 126603 (2004), doi:10.1103/PhysRevLett.92.126603.
16. J. Sinova, S. O. Valenzuela, J. Wunderlich, C. H. Back, and T. Jungwirth, “Spin Hall Effects,” Reviews of Modern Physics 87, 1213–1260 (2015), doi:10.1103/RevModPhys.87.1213.
17. Y. K. Kato, R. C. Myers, A. C. Gossard, and D. D. Awschalom, “Observation of the Spin Hall Effect in Semiconductors,” Science 306, 1910–1913 (2004), doi:10.1126/science.1105514.
18. J. Wunderlich, B. Kaestner, J. Sinova, and T. Jungwirth, “Experimental Observation of the Spin-Hall Effect in a Two-Dimensional Spin-Orbit Coupled Semiconductor System,” Physical Review Letters 94, 047204 (2005), doi:10.1103/PhysRevLett.94.047204.
19. E. G. Mishchenko, A. V. Shytov, and B. I. Halperin, “Spin Current and Polarization in Impure Two-Dimensional Electron Systems with Spin-Orbit Coupling,” Physical Review Letters 93, 226602 (2004), doi:10.1103/PhysRevLett.93.226602.
20. C. L. Kane and E. J. Mele, “Quantum Spin Hall Effect in Graphene,” Physical Review Letters 95, 226801 (2005), doi:10.1103/PhysRevLett.95.226801.
21. C. L. Kane and E. J. Mele, “$Z_2$ Topological Order and the Quantum Spin Hall Effect,” Physical Review Letters 95, 146802 (2005), doi:10.1103/PhysRevLett.95.146802.
22. L. Fu and C. L. Kane, “Topological Insulators with Inversion Symmetry,” Physical Review B 76, 045302 (2007), doi:10.1103/PhysRevB.76.045302.
23. M. Z. Hasan and C. L. Kane, “Colloquium: Topological Insulators,” Reviews of Modern Physics 82, 3045–3067 (2010), doi:10.1103/RevModPhys.82.3045.
24. X.-L. Qi and S.-C. Zhang, “Topological Insulators and Superconductors,” Reviews of Modern Physics 83, 1057–1110 (2011), doi:10.1103/RevModPhys.83.1057.
25. I. Žutić, J. Fabian, and S. Das Sarma, “Spintronics: Fundamentals and Applications,” Reviews of Modern Physics 76, 323–410 (2004), doi:10.1103/RevModPhys.76.323.