The standard Kohn–Sham DFT developed in Chapter 3 treats electrons as spinless particles: the density $\rho(\mathbf{r})$ carries no spin label, and the KS equations yield doubly degenerate orbital energies. This is adequate for non-magnetic systems, but fails entirely when the ground state breaks time-reversal symmetry through a net electron spin polarisation.
To describe magnetic materials — from ferromagnetic iron to the half-metallic Heusler alloys and antiferromagnetic insulators — we must extend DFT to spin-polarised DFT (SDFT), and, for materials with strong spin–orbit coupling, to the full non-collinear formalism. This chapter develops the theoretical framework and connects it to the key observables: magnetic moments, exchange coupling constants $J$, and magnetic anisotropy energies (MAE).
The Spin-Density Functional
Collinear SDFT
In collinear SDFT, spins are quantised along a fixed axis (conventionally $z$). The fundamental variable is replaced by the spin-resolved density: \begin{equation} \rho_\sigma(\mathbf{r}) = \sum_{i} f_{i\sigma},|\phi_{i\sigma}(\mathbf{r})|^2, \qquad \sigma \in {\uparrow, \downarrow}, \end{equation} where $f_{i\sigma}$ are the orbital occupation numbers and the sum runs over all KS orbitals of spin $\sigma$. The total and magnetisation densities are: \begin{equation} \rho(\mathbf{r}) = \rho_\uparrow(\mathbf{r}) + \rho_\downarrow(\mathbf{r}), \qquad m(\mathbf{r}) = \rho_\uparrow(\mathbf{r}) - \rho_\downarrow(\mathbf{r}). \end{equation}
The HK theorem is generalised to show that the ground state is uniquely determined by the pair $(\rho(\mathbf{r}), m(\mathbf{r}))$, or equivalently by $(\rho_\uparrow, \rho_\downarrow)$. The total energy functional becomes: \begin{equation} E[\rho_\uparrow, \rho_\downarrow] = T_s[\rho_\uparrow, \rho_\downarrow] + E_{\rm H}[\rho] + E_{\rm xc}[\rho_\uparrow, \rho_\downarrow] + \int V_{\rm ext}(\mathbf{r}),\rho(\mathbf{r}),d\mathbf{r}. \end{equation}
The two sets of spin-channel KS equations decouple: \begin{equation} \left[-\frac{1}{2}\nabla^2 + V_{\rm eff}^\sigma(\mathbf{r})\right]\phi_{i\sigma}(\mathbf{r}) = \epsilon_{i\sigma},\phi_{i\sigma}(\mathbf{r}), \label{eq:SDFT-KS} \end{equation} with spin-dependent effective potentials: \begin{equation} V_{\rm eff}^\sigma(\mathbf{r}) = V_{\rm ext}(\mathbf{r}) + V_{\rm H}(\mathbf{r}) + V_{\rm xc}^\sigma(\mathbf{r}), \qquad V_{\rm xc}^\sigma = \frac{\delta E_{\rm xc}}{\delta \rho_\sigma(\mathbf{r})}. \end{equation}
The Hartree potential $V_{\rm H}$ is spin-independent (it depends on the total $\rho$), while $V_{\rm xc}^\sigma$ differs between spin channels, producing a spin-splitting of the eigenvalues $\epsilon_{i\uparrow} \neq \epsilon_{i\downarrow}$. This is the quantum-mechanical origin of the exchange splitting that drives ferromagnetism.
The Spin-Polarised XC Functional
The XC functional must be generalised to take both $\rho_\uparrow$ and $\rho_\downarrow$ as arguments. For the LDA, the local spin-density approximation (LSDA) uses the XC energy of a spin-polarised uniform electron gas: \begin{equation} E_{\rm xc}^{\rm LSDA}[\rho_\uparrow, \rho_\downarrow] = \int \varepsilon_{\rm xc}^{\rm UEG}(\rho_\uparrow(\mathbf{r}), \rho_\downarrow(\mathbf{r})),\rho(\mathbf{r}),d\mathbf{r}, \end{equation} with the spin-scaling relation (von Barth–Hedin): \begin{equation} \varepsilon_{\rm xc}(\rho_\uparrow, \rho_\downarrow) = \varepsilon_{\rm xc}(\rho, 0) + [\varepsilon_{\rm xc}(\rho/2, \rho/2) - \varepsilon_{\rm xc}(\rho, 0)],f(\zeta), \end{equation} where $\zeta = (\rho_\uparrow - \rho_\downarrow)/\rho$ is the spin polarisation and $f(\zeta)$ interpolates between the paramagnetic ($\zeta=0$) and fully polarised ($\zeta=1$) limits. GGA and meta-GGA functionals are extended analogously.
Magnetic Moments
The local magnetic moment on site $A$ is: \begin{equation} \mu_A = \int_{\Omega_A} m(\mathbf{r}),d\mathbf{r} = \int_{\Omega_A} [\rho_\uparrow(\mathbf{r}) - \rho_\downarrow(\mathbf{r})],d\mathbf{r}, \end{equation} where $\Omega_A$ is the Wigner–Seitz (or PAW augmentation) sphere of atom $A$. The total magnetisation is: \begin{equation} M = \int m(\mathbf{r}),d\mathbf{r} = N_\uparrow - N_\downarrow, \end{equation} where $N_\sigma = \int \rho_\sigma,d\mathbf{r}$ is the total number of electrons of spin $\sigma$.
For elemental ferromagnets, SDFT recovers the experimental moments well: Fe ($2.2,\mu_B$ exp., $2.1$–$2.3,\mu_B$ DFT), Co ($1.7,\mu_B$), Ni ($0.6,\mu_B$). Orbital contributions (from spin–orbit coupling) are not included at the collinear SDFT level.
Magnetic Ordering: Ferro-, Antiferro-, and Ferrimagnetic States
By setting up different initial spin configurations in the SCF cycle, collinear SDFT can target different magnetic orderings:
- Ferromagnetic (FM): all local moments aligned parallel, $M \neq 0$.
- Antiferromagnetic (AFM): moments on sublattices aligned antiparallel, $M = 0$ globally. Requires a supercell large enough to contain both sublattices.
- Ferrimagnetic (FiM): unequal antiparallel moments on distinct sublattices, $M \neq 0$. Common in half-Heusler and spinel compounds.
To determine which ordering is the ground state, compute the total energy for each configuration and compare: $\Delta E = E_{\rm AFM} - E_{\rm FM}$.
Exchange Coupling: The Heisenberg Model
The energy difference between magnetic configurations can be mapped onto the Heisenberg spin Hamiltonian: \begin{equation} \hat{\mathcal{H}}{\rm Heis} = -\sum{\langle i,j\rangle} J_{ij},\hat{\mathbf{S}}_i \cdot \hat{\mathbf{S}}_j, \end{equation} where $\hat{\mathbf{S}}i$ is the spin operator on site $i$ and $J{ij}$ is the exchange coupling constant between sites $i$ and $j$. The sign convention: $J > 0$ favours FM alignment; $J \leq 0$ favours AFM.
Extraction of $J$ from Total Energies
For a two-sublattice system with spin $S$ per site, the energy difference between FM and AFM configurations (with $z$ equivalent nearest-neighbour pairs per formula unit) is: \begin{equation} \Delta E = E_{\rm AFM} - E_{\rm FM} = 2J,z,S^2. \label{eq:J-extraction} \end{equation}
Hence: \begin{equation} J = \frac{E_{\rm AFM} - E_{\rm FM}}{2,z,S^2}. \end{equation}
For more complex magnetic structures with multiple shells of neighbours, the four-state method (Xiang et al.) or the spin-spiral approach (Section below) provides a systematic way to disentangle $J_1, J_2, J_3, \ldots$
Liechtenstein–Katsnelson–Antropov–Gubanov (LKAG) Method
An alternative to total energy differences is the LKAG (linear response) method, which computes $J_{ij}$ directly from the Green’s function: \begin{equation} J_{ij} = \frac{1}{4\pi}\int_{-\infty}^{E_F} {\rm Im},{\rm Tr}L\left[\Delta_i G{ij}^\uparrow(E),\Delta_j G_{ji}^\downarrow(E)\right],dE, \end{equation} where $\Delta_i = V_{\rm xc,i}^\uparrow - V_{\rm xc,i}^\downarrow$ is the local exchange splitting on site $i$ and $G_{ij}^\sigma$ are the spin-resolved intersite Green’s functions. This method gives the full $J_{ij}$ tensor in a single DFT calculation and is available in codes such as FLEUR, SPR-KKR, and Questaal.
Practical note: Total energy methods (equation \eqref{eq:J-extraction}) are simpler and widely used. They require DFT+U or hybrid functionals for correlated $d$/$f$-electron systems where GGA overdelocalises the magnetic orbitals and underestimates $J$.
Non-Collinear SDFT
In the collinear formalism, the magnetisation is constrained to lie along $z$. This is insufficient for:
- Spin spirals and incommensurate magnetic structures.
- Spin frustration in triangular or kagome lattices.
- Spin–orbit coupling (SOC) effects: the magnetisation direction matters when SOC is present, and it may vary in space.
- Dzyaloshinskii–Moriya interaction (DMI) and skyrmions.
In non-collinear SDFT, the density variable is replaced by the full $2\times 2$ density matrix in spin space: \begin{equation} n_{\alpha\beta}(\mathbf{r}) = \sum_i f_i,\phi_i^\alpha(\mathbf{r})\left(\phi_i^\beta(\mathbf{r})\right)^*, \qquad \alpha,\beta \in {\uparrow,\downarrow}, \end{equation} which is related to the magnetisation density vector $\mathbf{m}(\mathbf{r})$ via: \begin{equation} n_{\alpha\beta} = \frac{1}{2}\left[\rho,\delta_{\alpha\beta} + \mathbf{m}\cdot\boldsymbol{\sigma}_{\alpha\beta}\right], \end{equation} where $\boldsymbol{\sigma} = (\sigma_x, \sigma_y, \sigma_z)$ are the Pauli matrices.
The KS orbitals become two-component spinors $\boldsymbol{\phi}_i = (\phi_i^\uparrow, \phi_i^\downarrow)^T$, and the KS equations take the matrix form: \begin{equation} \left[-\frac{1}{2}\nabla^2,\mathbf{I} + \mathbf{V}_{\rm eff}(\mathbf{r})\right]\boldsymbol{\phi}_i(\mathbf{r}) = \epsilon_i,\boldsymbol{\phi}_i(\mathbf{r}), \end{equation} with the $2\times 2$ effective potential: \begin{equation} \mathbf{V}_{\rm eff}(\mathbf{r}) = \left[V_{\rm ext} + V_{\rm H} + V_{\rm xc}^{(0)}\right]\mathbf{I} - \mathbf{B}_{\rm xc}(\mathbf{r})\cdot\boldsymbol{\sigma}, \end{equation} where $V_{\rm xc}^{(0)} = \frac{1}{2}(V_{\rm xc}^\uparrow + V_{\rm xc}^\downarrow)$ is the spin-averaged XC potential and $\mathbf{B}_{\rm xc} = \frac{1}{2}(V_{\rm xc}^\uparrow - V_{\rm xc}^\downarrow)\hat{\mathbf{m}}$ is the XC magnetic field aligned along the local magnetisation direction $\hat{\mathbf{m}}(\mathbf{r})$.
Spin–Orbit Coupling (SOC)
Spin–orbit coupling is a relativistic effect arising from the interaction of an electron’s spin magnetic moment with the magnetic field it experiences due to its orbital motion around the nucleus. In the Pauli approximation it takes the form: \begin{equation} \hat{H}_{\rm SOC} = \frac{\hbar^2}{4m_e^2c^2}\frac{1}{r}\frac{dV}{dr},\hat{\mathbf{L}}\cdot\hat{\mathbf{S}}, \end{equation} where $\hat{\mathbf{L}}$ and $\hat{\mathbf{S}}$ are the orbital and spin angular momentum operators, and $dV/dr$ is the radial derivative of the ionic potential. SOC is strong near heavy nuclei (scales as $Z^4$), making it important for $4d$, $5d$, $4f$, and $5f$ elements.
In SDFT, SOC is most commonly included via the second-variational method:
- Solve the scalar-relativistic (collinear) SDFT problem first.
- Construct the SOC matrix elements in the basis of scalar-relativistic eigenstates.
- Diagonalise the full Hamiltonian $\hat{H}{\rm SR} + \hat{H}{\rm SOC}$.
Full self-consistent inclusion of SOC requires the non-collinear spinor formalism.
Magnetic Anisotropy Energy (MAE)
The magnetic anisotropy energy is the difference in total energy between magnetisation oriented along the easy axis and a hard axis: \begin{equation} {\rm MAE} = E(\hat{\mathbf{n}}{\rm hard}) - E(\hat{\mathbf{n}}{\rm easy}), \end{equation} typically measured in $\mu$eV/atom for $3d$ metals or meV/atom for heavy-element compounds. MAE determines the magnetic hardness of a material: a large positive MAE (perpendicular easy axis) is required for permanent magnets and perpendicular magnetic recording media.
MAE arises entirely from spin–orbit coupling: without SOC, all directions are degenerate. It vanishes in spin-orbit-free collinear SDFT and requires either:
- Self-consistent non-collinear calculation with SOC for large MAE (e.g. $L1_0$-FePt, $> 1$ meV/atom).
- Force theorem (perturbative SOC): For small MAE, compute the non-self-consistent energy with SOC applied to the converged scalar-relativistic density. Computationally cheaper but breaks down when SOC significantly modifies the density.
The MAE is dominated by the band energy contribution (summed KS eigenvalue difference): \begin{equation} {\rm MAE} \approx \sum_i (f_i^{\rm hard},\epsilon_i^{\rm hard} - f_i^{\rm easy},\epsilon_i^{\rm easy}), \end{equation} which requires very dense $k$-meshes to converge (fine Brillouin zone features near the Fermi level). Typical convergence requires $40\times 40\times 40$ or finer grids for transition metal systems.
Contribution from Orbital Moments
By the second sum rule of XMCD (Bruno, 1989): \begin{equation} {\rm MAE} \approx -\frac{\xi}{4\mu_B}(\langle L_z^{\rm easy}\rangle - \langle L_z^{\rm hard}\rangle), \end{equation} where $\xi$ is the SOC constant and $\langle L_z\rangle$ is the orbital moment along the magnetisation axis. Larger orbital moments correlate with larger MAE — a design principle exploited in rare-earth permanent magnets (Nd$2$Fe${14}$B) and $L1_0$ alloys.
Application to Heusler Alloys
Heusler alloys ($X_2YZ$ full Heusler or $XYZ$ half-Heusler) are a rich playground for SDFT because they exhibit a wide range of magnetic phenomena:
Half-metallic ferromagnetism: In compounds like Co$_2$MnSi, one spin channel ($\uparrow$) is metallic while the other ($\downarrow$) has a band gap at $E_F$, giving a spin polarisation of 100% in principle. SDFT with GGA predicts this correctly for many Heuslers; a scissor correction or hybrid functional is needed to accurately reproduce the minority-spin gap.
Slater–Pauling rule: The total magnetic moment per formula unit follows $M = N_v - 24$ (full Heusler) or $M = N_v - 18$ (half-Heusler), where $N_v$ is the total valence electron count. This rule, reproduced by SDFT, provides a rapid screening criterion for half-metallicity.
Exchange coupling and $T_C$: The Curie temperature can be estimated from $J_{ij}$ values via mean-field theory ($k_B T_C^{\rm MF} = \frac{2}{3}J_0$) or random phase approximation (RPA), with RPA giving more accurate results for itinerant magnets. SDFT+LKAG is the standard approach.
MAE in Heuslers: Most cubic Heuslers have small MAE ($\sim \mu$eV/atom) due to cubic symmetry. Tetragonally distorted Heuslers (e.g. Mn$_3$Ga) or inverse Heuslers can have significantly larger MAE, relevant for spin-transfer torque applications.
Practical Checklist for Magnetic DFT Calculations
- Initialise spin correctly: Supply a non-zero initial magnetic moment to each magnetic site; the SCF loop will not spontaneously break spin symmetry from a non-magnetic starting point.
- Check for multiple magnetic minima: Try several initial configurations (FM, AFM, ferrimagnetic) and compare final energies. Collinear SDFT can get trapped in local minima.
- Use appropriate XC: PBE is standard; consider PBE+U or hybrid functionals for correlated $d/f$ systems where GGA under-localises the magnetic orbitals.
- For MAE: Use a non-collinear calculation with SOC; converge the $k$-mesh carefully (MAE is extremely sensitive to $k$-point density).
- For $J$: Cross-check total-energy and Green’s function methods; the two should agree within $\sim 10$–$20%$.
- Report orbital moments: When SOC is included, report both spin and orbital magnetic moments for comparison with XMCD experiments.
Outlook
Collinear SDFT provides an accurate and efficient framework for most magnetic materials, with non-collinear SDFT required for frustrated systems, SOC-driven anisotropy, and DMI. A systematic limitation shared by all SDFT methods is the treatment of strongly correlated $d$ and $f$ electrons, where the mean-field XC potential fails to capture Mott–Hubbard physics. The DFT+U method, which introduces an on-site Hubbard correction for localised orbitals, is the subject of the next chapter.