arxiv: 2602.16325 · v2 · submitted 2026-02-18 · ⚛️ physics.chem-ph

A Unified Formulation for langle hat{S}² rangle in Two-Component TDDFT

Xiaoyu Zhang This is my paper

Pith reviewed 2026-05-15 21:23 UTC · model grok-4.3

classification ⚛️ physics.chem-ph

keywords two-component TDDFTspin contaminationexcited stateslinear responseS squared expectation valuespin-flip TDDFTnoncollinear excitationsreference state decomposition

0 comments

The pith

Two-component TDDFT decomposes excited-state ⟨Ŝ²⟩ into the reference determinant's value plus an extra term from the excitation.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a single set of equations to compute the spin squared expectation value for electronically excited states in two-component time-dependent density functional theory. This framework covers noncollinear excitations from noncollinear references as well as both spin-conserving and spin-flip excitations from collinear references. The central result is that the excited-state ⟨Ŝ²⟩ is the sum of the reference state's ⟨Ŝ²⟩₀ and a distinct increment Δ⟨Ŝ²⟩ generated by the linear-response excitation. The same equations are specialized for collinear cases and tested on calculations that start from two-component DFT, unrestricted Kohn-Sham, and restricted open-shell Kohn-Sham determinants. This separation makes it possible to trace the origin of spin contamination to either the ground-state reference or the excitation step itself.

Core claim

Within the two-component TDDFT linear-response formalism, the expectation value ⟨Ŝ²⟩ of an excited state is expressed as ⟨Ŝ²⟩₀ + Δ⟨Ŝ²⟩, where ⟨Ŝ²⟩₀ is evaluated on the reference determinant and Δ⟨Ŝ²⟩ is the contribution arising from the excitation amplitudes. For collinear reference states the general expression reduces to the conventional forms used in spin-conserving and spin-flip TDDFT, allowing direct comparison with earlier derivations. Numerical results are obtained by applying the formalism to two-component DFT, UKS, and ROKS reference states.

What carries the argument

The general working equations for ⟨Ŝ²⟩ derived from the two-component TDDFT linear-response density matrix.

If this is right

For collinear references the unified equations recover separate spin-conserving and spin-flip TDDFT expressions for ⟨Ŝ²⟩.
The same ⟨Ŝ²⟩ formula applies unchanged to noncollinear reference states and noncollinear excitations.
Spin contamination can be partitioned into reference and excitation contributions for any of the three reference types (two-component DFT, UKS, ROKS).
Previous ad-hoc expressions for ⟨Ŝ²⟩ in spin-flip TDDFT are recovered as special cases of the general result.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The decomposition could guide the choice of reference state to minimize unwanted spin mixing before the excitation step is applied.
Similar partitioning might be useful in other response methods such as TD-HF or equation-of-motion coupled cluster to locate the source of spin contamination.
Numerical benchmarks on molecules with known exact multiplicities would show whether the Δ⟨Ŝ²⟩ term is dominated by the exchange-correlation kernel or by the orbital relaxation.

Load-bearing premise

The two-component TDDFT linear-response equations and their collinear decomposition capture the exact spin properties of the excited states without large higher-order or basis-set corrections.

What would settle it

Compute ⟨Ŝ²⟩ for the lowest spin-flip excitation of a small open-shell molecule such as the oxygen atom or methylene radical using both this formalism and a high-accuracy wave-function method; systematic deviation larger than basis-set error would falsify the decomposition.

read the original abstract

Two-component linear-response time-dependent density functional theory (TDDFT) provides a unified framework that encompasses noncollinear excitations in noncollinear reference states, as well as both spin-conserving and spin-flip excitations in collinear reference states. In this work, we present a general formalism for evaluating the expectation value $\langle \hat{S}^2 \rangle$ of electronically excited states obtained within two-component TDDFT. We then derive and analyze specialized forms of the resulting equations for collinear reference determinants, for which the two-component formalism decomposes into conventional spin-conserving and spin-flip TDDFT. The resulting working equations are systematically compared with previously proposed theoretical approaches. On the basis of our analysis, $\langle \hat{S}^2 \rangle$ in the excited states is shown to arise from two distinct sources: (i) $\langle \hat{S}^2 \rangle_0$ in the reference state and (ii) additional $\Delta\langle \hat{S}^2 \rangle$ introduced by the excitation process itself. Finally, we evaluate the expectation value $\langle \hat{S}^2 \rangle$ by performing two-component TDDFT calculations based on two-component DFT, unrestricted Kohn-Sham (UKS), and restricted open-shell Kohn-Sham (ROKS) reference states, respectively.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper gives a clean unified expression for excited-state <S²> in two-component TDDFT that splits into the reference contribution plus an additive term from the response amplitudes.

read the letter

The main thing to know is that the work derives a single set of equations for in two-component TDDFT excited states and shows it comes from two sources: the reference state's value plus a correction generated by the excitation amplitudes. This covers noncollinear cases directly and reduces to the usual spin-conserving and spin-flip formulas when the reference is collinear. They also line up the new equations against earlier specialized versions, which clarifies how the pieces fit together. Numerical runs with two-component DFT, UKS, and ROKS references illustrate the approach on concrete systems. The derivation itself is straightforward: insert the first-order response density into the exact spin operator and simplify. It recovers the known limits without extra assumptions or contradictions, so the additive split holds up on its own terms. The main soft spot is the validation side. The paper includes basic illustrations but does not show extensive benchmarks against exact results or other methods on standard test molecules, nor does it quantify how much the new expressions change practical outcomes in cases with strong spin contamination. That keeps the contribution more formal than applied. This is for computational chemists who run TDDFT on open-shell or magnetic systems and need consistent spin analysis for excited states. It fills a small but real gap in the two-component toolkit. I would send it to peer review; the central formalism is logically sound and the comparisons are useful, so referees can focus on tightening the numerical checks.

Referee Report

0 major / 2 minor

Summary. The manuscript develops a unified formalism for computing the expectation value ⟨Ŝ²⟩ of electronically excited states in two-component linear-response TDDFT. Starting from the exact spin operator expression and inserting the first-order response density, it derives general working equations that encompass noncollinear excitations in noncollinear references as well as spin-conserving and spin-flip excitations in collinear references. For collinear cases the equations reduce to the known specialized forms. The central result is the decomposition ⟨Ŝ²⟩ = ⟨Ŝ²⟩₀ + Δ⟨Ŝ²⟩, where ⟨Ŝ²⟩₀ is the reference-state contribution and Δ⟨Ŝ²⟩ is the additive correction from the TDDFT amplitudes. The formalism is compared to prior approaches and illustrated numerically with two-component DFT, UKS, and ROKS references.

Significance. If the derivation and implementation are correct, the work supplies a consistent, reference-independent route to quantify spin contamination in TDDFT excited states. The additive decomposition clarifies the physical origin of spin impurity and should improve the interpretability of results for open-shell systems and spin-flip processes. The systematic comparison to earlier methods also helps place the two-component framework in context.

minor comments (2)

[Comparison paragraph] The abstract states that the working equations are 'systematically compared with previously proposed theoretical approaches,' but the main text should explicitly name the specific prior formulations (with citations) at the point of comparison to make the differences transparent.
[Numerical results] The numerical illustrations would be strengthened by a compact table that reports ⟨Ŝ²⟩ values for the same test systems across the three reference types (two-component, UKS, ROKS) together with any available reference values or error metrics.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment and recommendation of minor revision. The provided summary accurately captures the scope and contributions of our work on the unified formalism for ⟨Ŝ²⟩ in two-component TDDFT.

Circularity Check

0 steps flagged

Derivation self-contained; no circularity in ⟨Ŝ²⟩ decomposition

full rationale

The paper derives ⟨Ŝ²⟩ for TDDFT excited states by direct substitution of the first-order response density matrix (obtained from standard two-component linear-response equations) into the exact definition of the spin operator Ŝ². This produces an additive split ⟨Ŝ²⟩ = ⟨Ŝ²⟩₀ + Δ⟨Ŝ²⟩ whose working equations are shown to recover the known spin-conserving and spin-flip limits for collinear references. No parameter is fitted to the target quantity, no self-citation supplies a uniqueness theorem or ansatz that forces the result, and the decomposition follows immediately from the operator definition plus the linear-response ansatz already present in TDDFT. The numerical examples with UKS/ROKS references serve only as illustration, not as input to the formalism. Consequently the central claim remains independent of its own outputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work rests on standard linear-response TDDFT assumptions without introducing new free parameters, invented entities, or ad-hoc axioms visible in the abstract.

axioms (1)

domain assumption Linear-response TDDFT equations apply to two-component systems and decompose for collinear references into spin-conserving and spin-flip forms.
The formalism extends established TDDFT theory to compute ⟨Ŝ²⟩.

pith-pipeline@v0.9.0 · 5532 in / 1190 out tokens · 32652 ms · 2026-05-15T21:23:56.704855+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

⟨Ŝ²⟩ in the excited states is shown to arise from two distinct sources: (i) ⟨Ŝ²⟩₀ in the reference state and (ii) additional Δ⟨Ŝ²⟩ introduced by the excitation process itself.
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

M s X s = Δ⟨Ŝ²⟩ N s X s (eq. 58) with Casida-like kernel K s derived from W_PQRS

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.