Spin-Orbital von Neumann Entropy

• Spin-orbital von Neumann entropy is a metric that quantifies the entanglement between spin and orbital states by evaluating the entropy of the reduced density matrix.
• It distinguishes classical product phases from entangled quantum states by displaying discontinuous jumps or continuous peak structures across quantum phase transitions.
• Spectral generalizations of the metric reveal universal scaling laws in excitations, linking theoretical models with experimental observables such as magnetic moments and RIXS signatures.

Spin-orbital von Neumann entropy (vNE) rigorously quantifies the entanglement between spin and orbital degrees of freedom in quantum systems where each elementary constituent carries both spin and orbital states. It is universally defined as the vNE of the reduced spin (or orbital) density matrix after tracing out the complementary subsystem from the pure state density operator. This metric, and its spectral generalizations, have revealed essential features of quantum phase transitions, elementary excitation structures, and the interplay of spin-orbit interactions in both theoretical models and real materials. Spin-orbital entanglement entropy underpins phase diagrams, dynamical signatures accessible to experiment (e.g., RIXS), and offers an invariant probe of genuine strong correlation distinct from spin-coupling effects.

1. Formal Definition and Construction

Consider a quantum system described by a pure state $\vert\Psi\rangle$ in a composite Hilbert space $$\mathcal{H} = \mathcal{H}_\mathrm{spin} \otimes \mathcal{H}_\mathrm{orbital}$$. The total density matrix is $\rho_\mathrm{tot} = |\Psi\rangle\langle\Psi|$. The reduced spin density matrix is obtained by tracing out the orbital degrees,

$\rho_S = \mathrm{Tr}_O[\rho_\mathrm{tot}]\,,$

and analogously, $\rho_O = \mathrm{Tr}_S[\rho_\mathrm{tot}]$.

The spin-orbital von Neumann entropy is

$$S^0_\mathrm{vN} = -\mathrm{Tr}_S\!\big[\rho_S \log_2 \rho_S\big] = -\mathrm{Tr}_O\!\big[\rho_O \log_2 \rho_O\big]\,,$$

which vanishes if and only if $|\Psi\rangle$ factorizes. Nonzero $S^0_\mathrm{vN}$ directly quantifies the degree of non-separability between spin and orbital sectors (You et al., 2015, You et al., 2015).

For excited states $|\Psi_n\rangle$, generalization yields the spectral entropy function

$$\mathcal{S}_\mathrm{vN}(\omega) = -\sum_n \mathrm{Tr}_S \big[\rho_S^{(n)} \log_2 \rho_S^{(n)}\big]\ \delta(\omega - \omega_n)\,,$$

where $$\mathcal{H} = \mathcal{H}_\mathrm{spin} \otimes \mathcal{H}_\mathrm{orbital}$$0 (You et al., 2015, You et al., 2012).

2. Entropy in Spin-Orbital Models: Ground States and Phases

In one-dimensional SU(2)$$\mathcal{H} = \mathcal{H}_\mathrm{spin} \otimes \mathcal{H}_\mathrm{orbital}$$1XXZ models with negative exchange coupling, the spin-orbital vNE differentiates classical product phases (I–IV), which have $$\mathcal{H} = \mathcal{H}_\mathrm{spin} \otimes \mathcal{H}_\mathrm{orbital}$$2, from entangled quantum phases. Notably:

• Majumdar-Ghosh-like dimer phase (V): At $$\mathcal{H} = \mathcal{H}_\mathrm{spin} \otimes \mathcal{H}_\mathrm{orbital}$$3, $$\mathcal{H} = \mathcal{H}_\mathrm{spin} \otimes \mathcal{H}_\mathrm{orbital}$$4, ground states manifest maximal spin-orbital entanglement ($$\mathcal{H} = \mathcal{H}_\mathrm{spin} \otimes \mathcal{H}_\mathrm{orbital}$$5), invariant with chain length. The transition into this phase is marked by a discontinuous jump in $$\mathcal{H} = \mathcal{H}_\mathrm{spin} \otimes \mathcal{H}_\mathrm{orbital}$$6 at critical $$\mathcal{H} = \mathcal{H}_\mathrm{spin} \otimes \mathcal{H}_\mathrm{orbital}$$7 values (You et al., 2015, You et al., 2015).
• Antiferromagnetic-spin/alternating-orbital phase (III): Finite vNE emerges when orbital fluctuations are present ($$\mathcal{H} = \mathcal{H}_\mathrm{spin} \otimes \mathcal{H}_\mathrm{orbital}$$8), showing area-law scaling. The III–V transition changes from first-order to continuous with increasing orbital exchange anisotropy; $$\mathcal{H} = \mathcal{H}_\mathrm{spin} \otimes \mathcal{H}_\mathrm{orbital}$$9 transitions smoothly, and $\rho_\mathrm{tot} = |\Psi\rangle\langle\Psi|$0 peaks rather than jumps.
• Orbital-dimer phase (VI): The $\rho_\mathrm{tot} = |\Psi\rangle\langle\Psi|$1 counterpart of the spin dimer phase, attaining a plateau $\rho_\mathrm{tot} = |\Psi\rangle\langle\Psi|$2 (You et al., 2015).

In the SU(2)$\rho_\mathrm{tot} = |\Psi\rangle\langle\Psi|$3SU(2) limit ($\rho_\mathrm{tot} = |\Psi\rangle\langle\Psi|$4), the entanglement entropy at the I–III transition line scales extensively with system size: $\rho_\mathrm{tot} = |\Psi\rangle\langle\Psi|$5, indicating breakdown of mean-field factorization.

3. Behavior Across Quantum Phase Transitions

The evolution of $\rho_\mathrm{tot} = |\Psi\rangle\langle\Psi|$6 and its spectral function across phase transitions provides a fingerprint of underlying entanglement changes:

• First-order transitions: $\rho_\mathrm{tot} = |\Psi\rangle\langle\Psi|$7 jumps discontinuously, e.g., II$\rho_\mathrm{tot} = |\Psi\rangle\langle\Psi|$8V boundary at $\rho_\mathrm{tot} = |\Psi\rangle\langle\Psi|$9.
• Continuous transitions: $\rho_S = \mathrm{Tr}_O[\rho_\mathrm{tot}]\,,$0 demonstrates peak structure in its derivative, with restructuring of low-lying entangled excitations (You et al., 2015).
• In the spectral entropy $\rho_S = \mathrm{Tr}_O[\rho_\mathrm{tot}]\,,$1, entangled ground states correspond to level crossings, while disentangled regimes exhibit zero-entropy excitations (pure magnons, orbitons). In dimer phases, even elementary spin excitations manifest finite spin-orbital entanglement.

4. Entanglement of Excitations and Universal Scaling

Beyond ground states, the vNE of excited states reveals universal properties:

• Bound states (BS) and on-site spin–orbit excitons (SOEX) exhibit logarithmic system-size scaling: $\rho_S = \mathrm{Tr}_O[\rho_\mathrm{tot}]\,,$2 with $\rho_S = \mathrm{Tr}_O[\rho_\mathrm{tot}]\,,$3 dependent on the mode (e.g., $\rho_S = \mathrm{Tr}_O[\rho_\mathrm{tot}]\,,$4 for SOEX, $\rho_S = \mathrm{Tr}_O[\rho_\mathrm{tot}]\,,$5 for BS branches) (You et al., 2015, You et al., 2012).
• Continuum states possess vNE saturating to O(1), highlighting the critical distinction between localized bound entanglement and extended excitations.
• This log-scaling is robust under variation of orbital anisotropy and exchange parameters and is experimentally testable via spectroscopic probes such as RIXS (You et al., 2012).

5. Spin-Orbital von Neumann Entropy in Quantum Materials

In single-ion systems such as $\rho_S = \mathrm{Tr}_O[\rho_\mathrm{tot}]\,,$6 compounds, the spin-orbital vNE is extracted from relativistic crystal field theory:

• The entropy is defined as

$\rho_S = \mathrm{Tr}_O[\rho_\mathrm{tot}]\,,$7

where $\rho_S = \mathrm{Tr}_O[\rho_\mathrm{tot}]\,,$8 is the total system vNE, and the difference isolates spin–orbital entanglement (García-Rojas et al., 22 Nov 2025).

• Empirical analysis shows that $\rho_S = \mathrm{Tr}_O[\rho_\mathrm{tot}]\,,$9 correlates directly with the measured effective magnetic moment $\rho_O = \mathrm{Tr}_S[\rho_\mathrm{tot}]$0, not with the spin–orbit coupling strength $\rho_O = \mathrm{Tr}_S[\rho_\mathrm{tot}]$1 or crystal field splitting $\rho_O = \mathrm{Tr}_S[\rho_\mathrm{tot}]$2 alone.
• Table: Representative values in select perovskites and molecular systems:

Compound	$\rho_O = \mathrm{Tr}_S[\rho_\mathrm{tot}]$3 (eV)	$\rho_O = \mathrm{Tr}_S[\rho_\mathrm{tot}]$4 (eV)	$\rho_O = \mathrm{Tr}_S[\rho_\mathrm{tot}]$5	$\rho_O = \mathrm{Tr}_S[\rho_\mathrm{tot}]$6	$\rho_O = \mathrm{Tr}_S[\rho_\mathrm{tot}]$7	$\rho_O = \mathrm{Tr}_S[\rho_\mathrm{tot}]$8
K$\rho_O = \mathrm{Tr}_S[\rho_\mathrm{tot}]$9TaCl$$S^0_\mathrm{vN} = -\mathrm{Tr}_S\!\big[\rho_S \log_2 \rho_S\big] = -\mathrm{Tr}_O\!\big[\rho_O \log_2 \rho_O\big]\,,$$0	0.211	0.293	1.046	–0.0020	0.0205	0.0522
Ba$$S^0_\mathrm{vN} = -\mathrm{Tr}_S\!\big[\rho_S \log_2 \rho_S\big] = -\mathrm{Tr}_O\!\big[\rho_O \log_2 \rho_O\big]\,,$$1MgReO$$S^0_\mathrm{vN} = -\mathrm{Tr}_S\!\big[\rho_S \log_2 \rho_S\big] = -\mathrm{Tr}_O\!\big[\rho_O \log_2 \rho_O\big]\,,$$2	0.584	0.498	0.787	–0.0115	0.0882	0.2112
Ba$$S^0_\mathrm{vN} = -\mathrm{Tr}_S\!\big[\rho_S \log_2 \rho_S\big] = -\mathrm{Tr}_O\!\big[\rho_O \log_2 \rho_O\big]\,,$$3NaOsO$$S^0_\mathrm{vN} = -\mathrm{Tr}_S\!\big[\rho_S \log_2 \rho_S\big] = -\mathrm{Tr}_O\!\big[\rho_O \log_2 \rho_O\big]\,,$$4	0.614	0.538	0.778	–0.0114	0.0874	0.2094

A plausible implication is that spin-orbital entanglement, as rigorously quantified by $$S^0_\mathrm{vN} = -\mathrm{Tr}_S\!\big[\rho_S \log_2 \rho_S\big] = -\mathrm{Tr}_O\!\big[\rho_O \log_2 \rho_O\big]\,,$$5, is tied to observable magnetic properties even for single-ion systems, and cannot be inferred solely from SOC strength (García-Rojas et al., 22 Nov 2025).

6. Spin-Free Orbital Entropy: Invariance and Correlation Diagnostics

In quantum chemistry, particularly in density matrix renormalization group (DMRG) frameworks, the raw spin-orbital entropy can conflate static spin-coupling (multiplet) effects with genuine spatial multireference correlation. Spin-free orbital entropy, defined by collapsing spin-resolved microstates into aggregated occupation states, yields quantities invariant with respect to spin projection $$S^0_\mathrm{vN} = -\mathrm{Tr}_S\!\big[\rho_S \log_2 \rho_S\big] = -\mathrm{Tr}_O\!\big[\rho_O \log_2 \rho_O\big]\,,$$6. This enables:

• Extraction of a three-state local basis entropy per orbital (occupied, singly-occupied, doubly-occupied).
• $$S^0_\mathrm{vN} = -\mathrm{Tr}_S\!\big[\rho_S \log_2 \rho_S\big] = -\mathrm{Tr}_O\!\big[\rho_O \log_2 \rho_O\big]\,,$$7-independent mapping of spatial correlation, discriminating true strong correlation from mere spin-coupling.
• Mututal information and pair entropy measures that robustly identify entangled active spaces (Pittner, 7 Feb 2025).

7. Experimental Access and Broader Context

Spin-orbital vNE and its spectral analogs underpin experimental identification of entangled excitations. Resonant inelastic x-ray scattering protocols directly probe dynamical spin-orbital correlation functions, which couple uniquely to entangled bound states with large vNE (You et al., 2012). Material-specific studies leveraging optical and magnetic measurements combined with crystal-field theory extract vNE fingerprints in transition-metal systems (García-Rojas et al., 22 Nov 2025).

This metric, whether in condensed-matter chains, quantum chemistry active spaces, or single-ion environments, serves as a universal quantifier of spin-orbital entanglement, illuminating underlying quantum mechanisms and emergent physical observables.