A Unified Formulation for $\langle \hat{S}^2 \rangle$ in Two-Component TDDFT

This paper presents a unified formalism for calculating the spin expectation value $\langle \hat{S}^2 \rangle$ in two-component time-dependent density functional theory (TDDFT), deriving specialized equations for collinear references to demonstrate that excited-state spin contamination arises from both the reference state and the excitation process itself.

The Gist

Imagine you are trying to describe the "spin" of an electron. In the quantum world, spin isn't just a tiny top spinning; it's a fundamental property that determines how electrons behave in magnetic fields and how they bond with each other. Usually, scientists treat electrons as if they are either "spin-up" or "spin-down," like coins that are either heads or tails. This works well for simple situations.

However, in complex molecules or when dealing with heavy atoms, electrons can do something trickier: they can exist in a messy mix of both heads and tails at the same time, or their spins can point in weird, diagonal directions. This is called a "noncollinear" state. To handle this, scientists use a sophisticated mathematical framework called Two-Component Time-Dependent Density Functional Theory (TDDFT). Think of this framework as a high-tech camera that can capture these messy, diagonal spins in 3D, rather than just a flat 2D picture.

The Problem: The "Spin Mess"
When scientists use this high-tech camera to look at excited states (electrons that have been kicked up to a higher energy level), they run into a problem. The math sometimes gets "contaminated." It's like trying to count the number of red and blue marbles in a jar, but the jar is slightly transparent, and you accidentally count some of the background light as marbles.

In quantum mechanics, we have a specific number we want to calculate called $\langle \hat{S}^2 \rangle$ (the expectation value of the total spin squared). This number tells us the "spin multiplicity"—essentially, whether the electrons are behaving like a calm, paired-up couple (a singlet) or a rowdy, unpaired group (a triplet). If the math is contaminated, this number comes out wrong, making it hard to know what kind of chemical reaction is actually happening.

The Solution: A Unified Recipe
Xiaoyu Zhang, the author of this paper, has written a new "recipe" (a unified formulation) to calculate this spin number correctly, no matter how messy the electron spins are.

Here is how the paper breaks it down, using simple analogies:

1. The Blueprint (Second Quantization):
   The author starts by rewriting the rules of spin using a language called "second quantization." Imagine the electrons as actors on a stage. Instead of describing the whole play at once, this method describes every single actor's entrance and exit. By doing this, the author shows that the math for calculating spin ($\hat{S}^2$) looks almost exactly like the math for calculating energy ($\hat{H}$). It's like realizing that the recipe for a cake is just a slightly modified version of the recipe for bread.

2. The Two Sources of Spin:
   The paper discovers that the total spin of an excited state comes from two distinct places:

   • The Base Spin ($\langle \hat{S}^2 \rangle_0$): This is the spin the molecule had before it got excited. It's the "foundation" of the building.
   • The Excitation Change ($\Delta \langle \hat{S}^2 \rangle$): This is the extra spin added or changed when the electron jumps to a higher energy level. It's the "renovation" done to the building.
     The paper provides a way to calculate both parts separately and then add them together to get the true total.

3. The "Casida" Machine:
   The author uses a mathematical machine known as the "Casida equation" (which is like a standard calculator for finding excited states in chemistry). Usually, this machine calculates energy. The author's big trick was to swap the "energy" settings on the machine for "spin" settings. Because the math is so similar, the machine can now spit out the spin numbers just as easily as it spits out energy numbers.

4. Testing the Recipe:
   To prove the recipe works, the author tested it on three different types of molecules:

   • Water ($H_2O$): A standard, stable molecule.
   • Water Ion ($H_2O^+$): A charged version of water.
   • Hydrogen Triplet ($H_3$): A tricky, unstable molecule where spins get very messy.

   The results showed that for simple molecules, the spin numbers were very clean. But for the messy $H_3$ molecule, the method correctly identified that the spins were "contaminated" (mixed up), which is a crucial piece of information for chemists trying to understand how these molecules react.

Why This Matters
Before this paper, if you wanted to know the spin of an excited electron in a complex, non-straightforward system, you might have had to use different, inconsistent methods depending on the situation. This paper provides one single, unified rulebook that works for all of them.

It's like having a universal translator that can speak every dialect of a language perfectly, whereas before, you needed a different translator for every village. This allows scientists to be much more confident when they are studying things like chemical reactions, how light interacts with matter, or how molecules behave in magnetic fields, ensuring they aren't being fooled by mathematical "noise."

In Summary
This paper gives scientists a new, reliable tool to measure the "spin" of excited electrons in complex systems. It breaks the measurement down into a "base" part and an "excitation" part, uses a clever mathematical swap to calculate it efficiently, and proves it works on a variety of test molecules. It doesn't promise to cure diseases or build new batteries directly, but it fixes a fundamental tool in the chemist's toolbox, ensuring that the theoretical maps they use to navigate the quantum world are accurate.

Technical Summary

1. Problem Statement

Linear-response Time-Dependent Density Functional Theory (TDDFT) is a standard tool for calculating excited states. However, in two-component formalisms (which handle noncollinear magnetism and relativistic effects), the spin operator $\hat{S}^2$ is not necessarily conserved, leading to spin contamination in excited states.

• The Gap: While methods exist to calculate $\langle \hat{S}^2 \rangle$ for spin-conserving TDDFT or specific spin-flip cases, there is no unified, general formalism for evaluating the expectation value $\langle \hat{S}^2 \rangle$ within the broader two-component TDDFT framework.
• The Need: Accurate calculation of $\langle \hat{S}^2 \rangle$ is crucial for determining spin multiplicity, analyzing spin contamination, and quantifying spin mixing (e.g., singlet-triplet admixture) in nonadiabatic molecular dynamics. Existing approaches often lack a systematic derivation that covers both collinear and noncollinear reference states within a single theoretical framework.

2. Methodology

The authors develop a rigorous theoretical framework based on second-quantized Equation-of-Motion (EOM) theory.

• Second-Quantized Representation:
  • The total spin-squared operator $\hat{S}^2$ is expressed in second quantization. The authors demonstrate that its algebraic structure parallels that of the electronic Hamiltonian $\hat{H}$.
  • $\hat{S}^2$ is decomposed into one-body and two-body terms involving spin operators $\hat{s}_u$ and creation/annihilation operators.
• General Formalism Derivation:
  • Reference State: The expectation value $\langle \hat{S}^2 \rangle_0$ for a general Slater determinant (reference state) is derived using Slater-Condon rules, expressed in terms of the one-particle reduced density matrix ($D$) and the overlap matrix ($S$) in a two-component basis.
  • Excited States: The change in spin expectation value, $\Delta \langle \hat{S}^2 \rangle$, is derived by substituting the Hamiltonian $\hat{H}$ with the operator $\hat{S}^2$ in the standard TDDFT linear-response (Casida) equations.
  • Casida-like Eigenvalue Equation: The authors derive a specific eigenvalue equation:
    $$\mathbf{M}^s \mathbf{X} = \Delta \langle \hat{S}^2 \rangle \mathbf{N} \mathbf{X}$$
    where $\mathbf{M}^s$ is the "spin-Fock" and "spin-kernel" matrix constructed from the second derivatives of the spin operator, and $\mathbf{X}$ represents the excitation vectors.
• Specialization to Collinear References:
  • The general formalism is specialized for collinear reference states (derived from Unrestricted Kohn-Sham (UKS) or Restricted Open-Shell Kohn-Sham (ROKS)).
  • In this limit, the two-component formalism decomposes into spin-conserving and spin-flip channels. The authors derive simplified working equations for these specific blocks, allowing direct comparison with previous literature.
• Noncollinear Kernels: The framework incorporates noncollinear exchange-correlation kernels derived via the multi-collinear approach, ensuring the theory preserves zero-energy excitations for open-shell systems (crucial for maintaining degeneracy between spin states).

3. Key Contributions

1. Unified Theory: The paper provides the first unified derivation for $\langle \hat{S}^2 \rangle$ in two-component TDDFT, applicable to both noncollinear reference states and collinear references (covering both spin-conserving and spin-flip excitations).
2. Decomposition of Spin Contamination: The analysis reveals that the total $\langle \hat{S}^2 \rangle$ in an excited state arises from two distinct sources:
   • $\langle \hat{S}^2 \rangle_0$: The intrinsic spin contamination of the reference state.
   • $\Delta \langle \hat{S}^2 \rangle$: The additional spin contamination introduced specifically by the excitation process.
3. Resolution of Ambiguities: The derived equations for collinear spin-flip TDDFT correct and clarify previous formulations (specifically addressing discrepancies in the work of Li et al. and Ipatov et al. regarding the definition of the metric matrix $\mathbf{N}$ and the handling of anticommutators).
4. Handling Zero-Energy Modes: The authors provide a rigorous treatment for zero-eigenvalue modes (improper modes arising from symmetry breaking), enforcing $\Delta \langle \hat{S}^2 \rangle = 0$ for these physically trivial modes to ensure numerical stability.

4. Results and Numerical Benchmarks

The methodology was implemented in an in-house JAX-based quantum chemistry code (IQC) and tested on three systems: H$_2$O, H$_2$O$^+$, and H$_3$.

• Systems: Calculations were performed using various functionals (HF, SVWN, PBE, B3LYP) and reference formalisms (GKS, UKS, ROKS).
• Observations:
  • H$_2$O / H$_2$O$^+$: These systems exhibited negligible spin contamination, with $\langle \hat{S}^2 \rangle$ values close to the theoretical limits for singlets and doublets.
  • H$_3$: This noncollinear molecule showed severe spin contamination in excited states, attributed to the incompleteness of the spin space in the basis set.
  • ROKS Anomaly: In H$_2$O$^+$ calculations using a ROKS reference, a negative excitation energy appeared. The authors explain this as a consequence of the ROKS reference not rigorously maintaining the zero-energy excitation (spin-flip transition) that is preserved in UKS/GKS due to the lack of full Fock matrix diagonalization.
• Consistency: The results for spin-conserving transitions matched established ab initio formulas (Myneni and Casida), validating the general derivation.

5. Significance

• Spectroscopy and Dynamics: The ability to accurately compute $\langle \hat{S}^2 \rangle$ is vital for interpreting spectroscopic data and for nonadiabatic molecular dynamics, where spin mixing drives electronic hopping probabilities.
• Relativistic Compatibility: While presented in the non-relativistic limit, the formalism is designed to be transferable to relativistic two-component methods (e.g., X2C) and potentially four-component TDDFT (for $\langle \hat{J}^2 \rangle$).
• Theoretical Foundation: By establishing a unified second-quantized framework, this work allows for the systematic calculation of expectation values for various observables beyond just spin, enhancing the predictive power of TDDFT for complex magnetic and relativistic systems.