Equilibrium Spin Density Matrix

• Equilibrium spin density matrix is a mathematical formulation that encodes the statistical and quantum characteristics of spin degrees of freedom at thermal equilibrium.
• It is constructed via expectation values of spin-resolved number operators and extends to electrons, relativistic fermions, and multi-spin systems under varying theoretical frameworks.
• Advanced computational methods, including DMRG and SRCAS, enable precise evaluations even in strongly correlated and spin-orbit coupled systems.

An equilibrium spin density matrix encodes the statistical and quantum structure of spin degrees of freedom in a system at thermal equilibrium. It provides the fundamental object for analyzing polarization, alignment, and coherence phenomena in quantum many-body theory, quantum chemistry, condensed matter, and relativistic quantum kinetic frameworks. The structure and computation of the equilibrium spin density matrix depend on the system class (e.g., electrons in molecules, relativistic fermions, spin-1 or higher systems) and the level of theoretical description (ab initio wave function, reduced system, kinetic theory, or statistical ensemble).

1. Fundamental Definitions and Formulations

In a basis of orthonormal spin-orbitals $\{\phi_{p}(\mathbf r)\otimes\sigma\}$, the equilibrium spin density matrix elements are defined via the expectation value of spin-resolved number operators. For electrons,

$$\rho_{pq}^{\sigma} = \langle\Psi| a_{p\sigma}^\dagger a_{q\sigma} | \Psi \rangle, \quad \sigma = \alpha, \beta,$$

with $a_{p\sigma}^\dagger$ ($a_{p\sigma}$) the creation (annihilation) operator for spin $\sigma$ in orbital $p$. The spin density matrix elements are then

$$T_{pq} = \langle\Psi|\hat T_{pq}|\Psi\rangle = \frac12(\rho_{pq}^{\alpha} - \rho_{pq}^{\beta}) = \frac12\langle\Psi|a_{p\alpha}^\dagger a_{q\alpha} - a_{p\beta}^\dagger a_{q\beta}| \Psi\rangle.$$

The spatially-resolved spin density operator reads

$$\hat{\delta}^{\rm spin}(\mathbf r) = \frac12\sum_{p,q} \phi_p^*(\mathbf r)\phi_q(\mathbf r)(a_{p\alpha}^\dagger a_{q\alpha} - a_{p\beta}^\dagger a_{q\beta}),$$

with its expectation value over $|\Psi\rangle$ yielding the real-space spin density profile

$$\rho^{\rm spin}(\mathbf r) = \sum_{p,q} \phi_p^*(\mathbf r)\phi_q(\mathbf r) T_{pq}.$$

This construction underpins both quantum chemical and condensed matter approaches (Boguslawski et al., 2012).

2. Equilibrium Spin Density in Quantum Statistical Mechanics

For finite quantum systems interacting with thermal baths, the canonical density matrix takes the form

$$\rho_{pq}^{\sigma} = \langle\Psi| a_{p\sigma}^\dagger a_{q\sigma} | \Psi \rangle, \quad \sigma = \alpha, \beta,$$0

where $$\rho_{pq}^{\sigma} = \langle\Psi| a_{p\sigma}^\dagger a_{q\sigma} | \Psi \rangle, \quad \sigma = \alpha, \beta,$$1. For spin-$$\rho_{pq}^{\sigma} = \langle\Psi| a_{p\sigma}^\dagger a_{q\sigma} | \Psi \rangle, \quad \sigma = \alpha, \beta,$$2 systems in magnetic fields, this specializes to

$$\rho_{pq}^{\sigma} = \langle\Psi| a_{p\sigma}^\dagger a_{q\sigma} | \Psi \rangle, \quad \sigma = \alpha, \beta,$$3

For $$\rho_{pq}^{\sigma} = \langle\Psi| a_{p\sigma}^\dagger a_{q\sigma} | \Psi \rangle, \quad \sigma = \alpha, \beta,$$4 and $$\rho_{pq}^{\sigma} = \langle\Psi| a_{p\sigma}^\dagger a_{q\sigma} | \Psi \rangle, \quad \sigma = \alpha, \beta,$$5,

$$\rho_{pq}^{\sigma} = \langle\Psi| a_{p\sigma}^\dagger a_{q\sigma} | \Psi \rangle, \quad \sigma = \alpha, \beta,$$6

Canonical typicality theorems guarantee that for weak coupling and large baths, the reduced and conditional density matrices closely concentrate around the above canonical form—specifically, the full statistical ensemble yields a deterministic equilibrium spin density matrix with exponentially suppressed fluctuations in system-bath scenarios (Pandya et al., 2013).

3. Relativistic and Hydrodynamic Equilibrium Spin Density Matrices

In relativistic quantum kinetic theory, the equilibrium spin density matrix for a spin-$$\rho_{pq}^{\sigma} = \langle\Psi| a_{p\sigma}^\dagger a_{q\sigma} | \Psi \rangle, \quad \sigma = \alpha, \beta,$$7 system is obtained via the local equilibrium Wigner function formalism. The Wigner function is expressed in terms of matrix-valued (spinor) phase-space distributions: $$\rho_{pq}^{\sigma} = \langle\Psi| a_{p\sigma}^\dagger a_{q\sigma} | \Psi \rangle, \quad \sigma = \alpha, \beta,$$8 The equilibrium (Becattini–Ferroni) spin density matrix is

$$\rho_{pq}^{\sigma} = \langle\Psi| a_{p\sigma}^\dagger a_{q\sigma} | \Psi \rangle, \quad \sigma = \alpha, \beta,$$9

with $a_{p\sigma}^\dagger$0 the "spin potential" (thermal vorticity in global equilibrium) and $a_{p\sigma}^\dagger$1 (Peng et al., 2021). The Dirac-matrix representation,

$a_{p\sigma}^\dagger$2

ensures direct access to scalar and axial (spin) components; the Pauli–Lubanski vector encodes physical polarization.

Advanced treatments correct known deficiencies in normalization for large spin polarization by expressing the spin density matrix in exponential form involving the axial vector constructed from the dual spin potential, ensuring $a_{p\sigma}^\dagger$3 for arbitrary $a_{p\sigma}^\dagger$4 (Bhadury et al., 5 May 2025).

4. Many-Body and Strongly Correlated Regimes

For quantum many-body systems (e.g., correlated Fermi liquids, open quantum spin systems), the equilibrium spin density matrix can be derived via cluster expansions or via trace over the environment:

• In correlated Fermi liquids, the one-body spin-resolved density matrix at finite temperature is given by (Serhan, 2010)

$a_{p\sigma}^\dagger$5

with explicit inclusion of Jastrow correlations and FHNC cluster resummation producing spin- and temperature-dependent momentum distributions.

• In open quantum systems, the equilibrium reduced density matrix (mean-force density) of a spin system $a_{p\sigma}^\dagger$6 coupled to environment $a_{p\sigma}^\dagger$7 is

$a_{p\sigma}^\dagger$8

which can be realized numerically via stochastic trace estimators and Krylov subspace projection, yielding the effective mean-force Hamiltonian $a_{p\sigma}^\dagger$9 and thus the equilibrium spin density matrix even for strong coupling (Chen et al., 2022).

For spin–boson models, the exact equilibrium reduced density matrix is obtained via imaginary-time path integrals or ordered exponentials, incorporating an influence functional determined by the bath spectral density, and encoding full system–bath entanglement (1712.06397).

5. Equilibrium Spin Density Matrix in Systems with Spin-Orbit Coupling and Disorder

In systems with strong spin-orbit interaction or disorder, the equilibrium spin density matrix contains both diagonal (population) and off-diagonal (coherence) elements in band indices. For a general multi-band system, the equilibrium spin density matrix in the crystal momentum ($a_{p\sigma}$0) domain is (Xiao et al., 2018)

$a_{p\sigma}$1

The first term encodes pure spin populations, while the second captures disorder-induced interband coherences essential for phenomena such as side-jump contributions to spin and anomalous Hall effects.

6. Spin-1 and Higher Spin: Alignment and Tensor Polarization

For spin-1 systems, the equilibrium spin density matrix naturally decomposes into vector (polarization) and tensor (alignment) components in the adjoint spin representation. The equilibrium spin density matrix is

$a_{p\sigma}$2

with normalization $a_{p\sigma}$3, $a_{p\sigma}$4, $a_{p\sigma}$5, and $a_{p\sigma}$6 (Florkowski et al., 31 Jan 2026). This structure leads to a density matrix characterized in the laboratory frame by both mean spin polarization $a_{p\sigma}$7 and alignment tensor $a_{p\sigma}$8, with direct implications for observables like the $a_{p\sigma}$9 spin alignment parameter probed in vector meson polarization experiments.

Gradient expansions (in the hydrodynamic or nonequilibrium context) reveal that, to first order, both thermal vorticity and shear do not shift the $\sigma$0 component in vector mesons away from $\sigma$1, with leading deviations appearing only at second order in gradients (Yang et al., 2024).

7. Computational Approaches and Convergence Criteria

The density-matrix renormalization group (DMRG) and sampling-reconstruction complete active space (SRCAS) schemes provide ab initio methods for constructing equilibrium spin densities in quantum chemistry. In DMRG, the spin density matrix elements are computed as explicit expectation values in matrix product states, while SRCAS reconstructs a CASCI wave function from sampled determinants, allowing direct evaluation and convergence checks based on the completeness measure

$\sigma$2

Both methods yield spin density matrices that converge systematically to the exact result as the active space and sampling threshold are extended, as demonstrated in challenging iron-nitrosyl molecular systems (Boguslawski et al., 2012).

---

Equilibrium spin density matrices thus constitute a central construct bridging quantum statistical mechanics, many-body theory, relativistic hydrodynamics, and quantum chemistry—governing the macroscopic signatures of spin polarization, alignment, coherence, and emergent dynamical phenomena across diverse physical settings.