pith. sign in

arxiv: 2605.19361 · v2 · pith:MD3WAEUOnew · submitted 2026-05-19 · ⚛️ physics.chem-ph

Spin-Adapted TDDFT

Pith reviewed 2026-06-30 18:01 UTC · model grok-4.3

classification ⚛️ physics.chem-ph
keywords TDDFTspin-adaptedRPAexcited statesspin contaminationhybrid functionalsconical intersections
0
0 comments X

The pith

Spin-adapted TDDFT extends RPA by casting its Fock matrix and kernels as energy derivatives to remove spin contamination.

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

The paper develops a spin-adapted form of time-dependent density functional theory for calculating excited states. It begins with an existing spin-adapted random phase approximation that uses tensor decoupling and the Wigner-Eckart theorem. Recasting the RPA Fock matrix and kernels as derivatives of the total energy produces the corresponding TDDFT equations. Hybrid functionals that combine Hartree-Fock exchange with spin-unpolarized exchange-correlation parts are introduced to enforce degeneracy among spin components. Benchmarks on chromium dimer dissociation and phenol O-H conical intersections illustrate the resulting behavior for both spin-conserving and spin-flip transitions.

Core claim

Casting the RPA Fock matrix and kernels as energy derivatives gives a TDDFT extension. To restore spin-component degeneracy, hybrids combining HF exchange with spin-unpolarized pure XC parts are used. This removes spin contamination that otherwise appears in linear-response TDDFT for spin-conserving and spin-flip channels.

What carries the argument

Energy-derivative formulation of the RPA Fock matrix and kernels, extended to TDDFT with hybrid functionals that mix HF exchange and spin-unpolarized XC.

If this is right

  • Spin-conserving and spin-flip excitation energies become free of contamination artifacts.
  • Bond-dissociation curves for open-shell systems such as Cr2 become more reliable.
  • Conical intersections involving different spin states, as in phenol O-H, are described without spurious splitting.
  • The same hybrid construction can be used for both pure and hybrid exchange-correlation functionals in the TDDFT kernel.

Where Pith is reading between the lines

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

  • The energy-derivative route may allow similar spin adaptation to be applied to other response properties beyond linear TDDFT.
  • Photochemical simulations that rely on accurate spin-state ordering could adopt the method with minimal change to existing TDDFT codes.
  • The approach might be tested on systems with strong spin-orbit coupling to check whether degeneracy restoration persists.

Load-bearing premise

The tensor decoupling and Wigner-Eckart theorem used in spin-adapted RPA transfer directly to the TDDFT equations obtained from energy derivatives without creating new spin-contamination artifacts.

What would settle it

A calculation on a closed-shell molecule in which the TDDFT response energies for equivalent spin components differ after applying the hybrid construction would show that the transfer of the RPA spin adaptation has failed.

Figures

Figures reproduced from arXiv: 2605.19361 by Tai Wang, Xiaoyu Zhang.

Figure 1
Figure 1. Figure 1: FIG. 1. Near-equilibrium potential-energy curves of Cr [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Relative energies of the three lowest singlet states [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
read the original abstract

Linear-response TDDFT is widely used for excited states but suffers from spin contamination in spin-conserving and spin-flip channels. Spin-adapted RPA was developed using tensor equation-of-motion and applying the Wigner-Eckart theorem with tensor decoupling. Casting the RPA Fock matrix and kernels as energy derivatives gives a TDDFT extension. To restore spin-component degeneracy, we use hybrids combining HF exchange with spin-unpolarized pure XC parts. Benchmarks, Cr$_2$ dissociation, and phenol O-H conical intersections demonstrate the method.

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.

Referee Report

3 major / 2 minor

Summary. The manuscript develops a spin-adapted linear-response TDDFT by casting the RPA Fock matrix and exchange-correlation kernels as energy derivatives, thereby extending the tensor equation-of-motion and Wigner-Eckart theorem framework previously used for spin-adapted RPA. Hybrids that combine Hartree-Fock exchange with spin-unpolarized pure XC functionals are introduced to restore spin-component degeneracy. The resulting method is applied to standard benchmarks, the Cr₂ dissociation curve, and the phenol O-H conical intersection.

Significance. If the derivative construction preserves spin purity without new contamination artifacts, the approach would supply a formally grounded route to spin-adapted TDDFT excitations, addressing a known limitation of conventional TDDFT for open-shell and photochemical applications. The explicit link to RPA via energy derivatives and the concrete demonstrations on bond dissociation and conical intersections constitute practical strengths.

major comments (3)
  1. [TDDFT extension section] TDDFT extension section: the central step of representing RPA quantities as energy derivatives is asserted to inherit the tensor decoupling and Wigner-Eckart spin adaptation without additional mixing, yet no explicit verification (analytic or numerical) is supplied that the derivative operation commutes with the spin-projection operators; this assumption is load-bearing for the claim of restored spin purity.
  2. [Hybrid construction paragraph] Hybrid construction paragraph: the use of HF exchange plus spin-unpolarized XC is presented as restoring degeneracy, but the manuscript does not quantify the residual spin contamination (e.g., via <S²> expectation values or degeneracy splitting) before versus after the hybrid correction, leaving the effectiveness of this fix unverified for the TDDFT linear-response equations.
  3. [Cr₂ dissociation benchmark] Cr₂ dissociation benchmark: while the curve is shown, the paper does not report the spin contamination metric or the degeneracy of the spin components along the dissociation path, which is required to substantiate that the spin-adaptation survives the energy-derivative extension in a strongly correlated regime.
minor comments (2)
  1. [Abstract and introduction] The abstract and introduction would benefit from a concise statement of the precise difference between the present energy-derivative route and earlier spin-adapted TDDFT formulations in the literature.
  2. [Notation section] Notation for the spin-adapted kernels should be introduced with an explicit reference back to the corresponding RPA definitions to aid readability.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their detailed and insightful comments on our manuscript. We have carefully considered each point and provide point-by-point responses below. Where appropriate, we will revise the manuscript to address the concerns raised.

read point-by-point responses
  1. Referee: [TDDFT extension section] the central step of representing RPA quantities as energy derivatives is asserted to inherit the tensor decoupling and Wigner-Eckart spin adaptation without additional mixing, yet no explicit verification (analytic or numerical) is supplied that the derivative operation commutes with the spin-projection operators; this assumption is load-bearing for the claim of restored spin purity.

    Authors: We agree that an explicit verification would strengthen the presentation. The energy derivatives are taken with respect to the density matrix elements that are themselves constructed as spin-adapted tensors from the tensor EOM-RPA framework. Because the derivative is a linear operation and the spin-projection operators act on the same tensor space, the commutation holds by construction. In the revised manuscript, we will add a short analytic paragraph demonstrating this commutation property in the TDDFT extension section. revision: yes

  2. Referee: [Hybrid construction paragraph] the use of HF exchange plus spin-unpolarized XC is presented as restoring degeneracy, but the manuscript does not quantify the residual spin contamination (e.g., via <S²> expectation values or degeneracy splitting) before versus after the hybrid correction, leaving the effectiveness of this fix unverified for the TDDFT linear-response equations.

    Authors: The referee is correct that quantitative measures of the improvement are valuable. We will add a table or supplementary data reporting <S²> values and the splitting between spin components for representative states in the benchmarks, both with and without the hybrid correction, to demonstrate the restoration of degeneracy. revision: yes

  3. Referee: [Cr₂ dissociation benchmark] while the curve is shown, the paper does not report the spin contamination metric or the degeneracy of the spin components along the dissociation path, which is required to substantiate that the spin-adaptation survives the energy-derivative extension in a strongly correlated regime.

    Authors: We acknowledge this omission. In the revised version, we will include additional panels or a table showing the <S²> expectation values and the energy degeneracy of the spin multiplet components as a function of the Cr-Cr distance for the spin-adapted TDDFT calculations. revision: yes

Circularity Check

0 steps flagged

No circularity detected; derivation is self-contained

full rationale

The paper's central step of casting the RPA Fock matrix and kernels as energy derivatives to obtain a TDDFT extension is presented as a direct methodological mapping from an established spin-adapted RPA framework (using tensor EOM and Wigner-Eckart theorem) without any quoted reduction of the TDDFT linear-response equations to fitted parameters, self-defined quantities, or load-bearing self-citations. The hybrid construction for restoring spin degeneracy is described as an additional practical choice rather than a renaming or ansatz smuggled via prior work. No equations in the abstract or description exhibit a prediction that equals its input by construction, and the derivation chain remains independent of the target result.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review provides no identifiable free parameters, axioms, or invented entities; no specific numbers fitted to data or unproved background assumptions can be extracted.

pith-pipeline@v0.9.1-grok · 5597 in / 1093 out tokens · 40991 ms · 2026-06-30T18:01:58.495950+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

45 extracted references · 6 canonical work pages · 1 internal anchor

  1. [1]

    Spin-Adapted TDDFT

    In addition, MRSF-TDDFT employs an empirical parameter for the pairing strength,c pq. In short, spin contamination in TDDFT remains a puzzle. We make a comparison among previously proposed methods in table I. TABLE I. Comparison of representative TDDFT-based ap- proaches for spin purification. Method spin purity ∆Sempirical parameters X-TDA/X-SF-TDA [12, ...

  2. [2]

    The development of more accurate hybrid ex- change–correlation functionals, in which the pure FIG. 2. Relative energies of the three lowest singlet states of phenol along the O–H dissociation coordinate, computed with SA-TDDFT/BHHLYP. Top: full scan ofS 0,S 1, andS 2 as functions ofR OH. Bottom left: enlarged view of theS 1/S2 near-degeneracy region aroun...

  3. [3]

    The integration of the present spin-adapted TDDFT framework with spin–orbit coupling (SOC). [36]

  4. [4]

    [37] Such gradients would enable the subsequent development of nonadiabatic coupling schemes

    The derivation and implementation of analytic en- ergy gradients within the spin-adapted TDDFT formalism. [37] Such gradients would enable the subsequent development of nonadiabatic coupling schemes. [38] We gratefully acknowledge the financial support and intellectually stimulating scientific discussions provided by Yunlong Xiao. This work was supported ...

  5. [5]

    Time-dependent density functional re- sponse theory for molecules,

    Mark E. Casida, “Time-dependent density functional re- sponse theory for molecules,” (1995) pp. 155–192. 6

  6. [6]

    Excited-state spin-contamination in time-dependent density-functional theory for molecules with open-shell ground states,

    Andrei Ipatov, Felipe Cordova, Lo¨ ıc Joubert Doriol, and Mark E. Casida, “Excited-state spin-contamination in time-dependent density-functional theory for molecules with open-shell ground states,” Journal of Molecu- lar Structure: THEOCHEM914, 60–73 (2009), time- dependent density-functional theory for molecules and molecular solids

  7. [7]

    Noncollinear spin-flip tddft for potential energy surface crossings: Conical intersections and spin crossings,

    Xiaoyu Zhang, Tai Wang, Yi Qin Gao, and Yun- long Xiao, “Noncollinear spin-flip tddft for potential energy surface crossings: Conical intersections and spin crossings,” Journal of Chemical Theory and Com- putation21, 11550–11561 (2025), pMID: 41208132, https://doi.org/10.1021/acs.jctc.5c01272

  8. [8]

    Cluster expansion of the wavefunction. symmetry-adapted-cluster expansion, its variational determination, and extension of open-shell or- bital theory,

    H. Nakatsuji and K. Hirao, “Cluster expansion of the wavefunction. symmetry-adapted-cluster expansion, its variational determination, and extension of open-shell or- bital theory,” The Journal of Chemical Physics68, 2053– 2065 (1978)

  9. [9]

    Cluster expansion of the wavefunction. pseudo-orbital theory based on the sac ex- pansion and its application to the spin density of open- shell systems,

    H. Nakatsuji and K. Hirao, “Cluster expansion of the wavefunction. pseudo-orbital theory based on the sac ex- pansion and its application to the spin density of open- shell systems,” The Journal of Chemical Physics68, 4279–4291 (1978)

  10. [10]

    Generalizations of the multiconfigurational time- dependent hartree–fock approach,

    Danny L. Yeager, Jeppe Olsen, and Poul Jørgensen, “Generalizations of the multiconfigurational time- dependent hartree–fock approach,” Faraday Symp. Chem. Soc.19, 85–95 (1984)

  11. [11]

    Polarization propagator methods in atomic and molec- ular calculations,

    Jens Oddershede, Poul Jørgensen, and Danny L. Yeager, “Polarization propagator methods in atomic and molec- ular calculations,” Computer Physics Reports2, 33–92 (1984)

  12. [12]

    General excita- tions in time-dependent density functional theory,

    Olav Vahtras and Zilvinas Rinkevicius, “General excita- tions in time-dependent density functional theory,” The Journal of Chemical Physics126, 114101 (2007)

  13. [13]

    A spin-complete version of the spin-flip approach to bond breaking: What is the impact of obtaining spin eigen- functions?

    John S. Sears, C. David Sherrill, and Anna I. Krylov, “A spin-complete version of the spin-flip approach to bond breaking: What is the impact of obtaining spin eigen- functions?” The Journal of Chemical Physics118, 9084– 9094 (2003)

  14. [14]

    The spin- flip extended single excitation configuration interaction method,

    David Casanova and Martin Head-Gordon, “The spin- flip extended single excitation configuration interaction method,” The Journal of Chemical Physics129, 064104 (2008)

  15. [15]

    Spin-adapted open-shell random phase approximation and time-dependent den- sity functional theory. i. theory,

    Zhendong Li and Wenjian Liu, “Spin-adapted open-shell random phase approximation and time-dependent den- sity functional theory. i. theory,” The Journal of Chemi- cal Physics133, 064106 (2010)

  16. [16]

    Spin-adapted open-shell time-dependent density functional theory. iii. an even better and simpler formulation,

    Zhendong Li and Wenjian Liu, “Spin-adapted open-shell time-dependent density functional theory. iii. an even better and simpler formulation,” The Journal of Chemi- cal Physics135, 194106 (2011)

  17. [17]

    Spin-flip, ten- sor equation-of-motion configuration interaction with a density-functional correction: A spin-complete method for exploring excited-state potential energy surfaces,

    Xing Zhang and John M. Herbert, “Spin-flip, ten- sor equation-of-motion configuration interaction with a density-functional correction: A spin-complete method for exploring excited-state potential energy surfaces,” The Journal of Chemical Physics143, 234107 (2015)

  18. [18]

    Spin-adapted spin-flip-down time-dependent density functional the- ory,

    Chima S. Chibueze and Lucas Visscher, “Spin-adapted spin-flip-down time-dependent density functional the- ory,” The Journal of Chemical Physics163, 094111 (2025)

  19. [19]

    Spin-adapted open- shell time-dependent density functional theory: towards a simple and accurate method for spin-flip-down excita- tions,

    Hewang Zhao and Zhendong Li, “Spin-adapted open- shell time-dependent density functional theory: towards a simple and accurate method for spin-flip-down excita- tions,” Molecular Physics0, e2631735 (2026)

  20. [20]

    Seunghoon Lee, Michael Filatov, Sangyoub Lee, and Cheol Ho Choi, “Eliminating spin-contamination of spin- flip time dependent density functional theory within lin- ear response formalism by the use of zeroth-order mixed- reference (mr) reduced density matrix,” The Journal of Chemical Physics149, 104101 (2018)

  21. [21]

    The spin–flip approach within time-dependent den- sity functional theory: Theory and applications to diradicals,

    Yihan Shao, Martin Head-Gordon, and Anna I. Krylov, “The spin–flip approach within time-dependent den- sity functional theory: Theory and applications to diradicals,” The Journal of Chemical Physics118, 4807–4818 (2003), https://pubs.aip.org/aip/jcp/article- pdf/118/11/4807/19210740/4807 1 online.pdf

  22. [22]

    Zero excitation energy theorem and the spin- flip kernel,

    Tai Wang, Hao Li, Yi Qin Gao, and Yunlong Xiao, “Zero excitation energy theorem and the spin- flip kernel,” Journal of Chemical Theory and Com- putation21, 6905–6921 (2025), pMID: 40638888, https://doi.org/10.1021/acs.jctc.5c00714

  23. [23]

    Higher random- phase approximation as an approximation to the equa- tions of motion,

    Tai-Ichi Shibuya and Vincent McKoy, “Higher random- phase approximation as an approximation to the equa- tions of motion,” Phys. Rev. A2, 2208–2218 (1970)

  24. [24]

    Equations-of-motion method and the extended shell model,

    D. J. ROWE, “Equations-of-motion method and the extended shell model,” Rev. Mod. Phys.40, 153–166 (1968)

  25. [25]

    Tensor equations of motion for the excitations of rotationally invariant or charge-independent systems,

    D. J. Rowe and C. Ngo-Trong, “Tensor equations of motion for the excitations of rotationally invariant or charge-independent systems,” Rev. Mod. Phys.47, 471– 485 (1975)

  26. [26]

    Noncollinear density functional theory,

    Zhichen Pu, Hao Li, Ning Zhang, Hong Jiang, Yiqin Gao, Yunlong Xiao, Qiming Sun, Yong Zhang, and Sihong Shao, “Noncollinear density functional theory,” Phys. Rev. Res.5, 013036 (2023)

  27. [27]

    Operator formal- ism for noncollinear functionals in the multicollinear approach,

    Xiaoyu Zhang and Taoni Bao, “Operator formal- ism for noncollinear functionals in the multicollinear approach,” Journal of Chemical Theory and Com- putation21, 9620–9630 (2025), pMID: 41039659, https://doi.org/10.1021/acs.jctc.5c01305

  28. [28]

    Noncollinear and spin-flip tddft in mul- ticollinear approach,

    Hao Li, Zhichen Pu, Qiming Sun, Yi Qin Gao, and Yunlong Xiao, “Noncollinear and spin-flip tddft in mul- ticollinear approach,” Journal of Chemical Theory and Computation19, 2270–2281 (2023)

  29. [29]

    A unified formulation for〈 ˆS2〉in two- component tddft,

    Xiaoyu Zhang, “A unified formulation for〈 ˆS2〉in two- component tddft,” Journal of Chemical Theory and Com- putation22, 4429–4438 (2026)

  30. [30]

    A mountaineering strategy to excited states: Highly accurate reference energies and benchmarks,

    Pierre-Fran¸ cois Loos, Anthony Scemama, Aymeric Blon- del, Yann Garniron, Michel Caffarel, and Denis Jacquemin, “A mountaineering strategy to excited states: Highly accurate reference energies and benchmarks,” Journal of Chemical Theory and Computation14, 4360– 4379 (2018)

  31. [31]

    XTDDFT: Spin- adapted td-dft for open-shell systems,

    Quantum-Chemistry-Group-BNU, “XTDDFT: Spin- adapted td-dft for open-shell systems,” (2026), commit 02efd408ba9d8c87be7a92f19e12a4c162438ed5, accessed 2026-05-22

  32. [32]

    OpenQP: Open quantum platform,

    Open-Quantum-Platform, “OpenQP: Open quantum platform,”https://github.com/ Open-Quantum-Platform/openqp, gitHub repository, accessed 2026-05-22

  33. [33]

    The chromium dimer: Closing a chapter of quantum chemistry,

    Henrik R. Larsson, Huanchen Zhai, C. J. Umrigar, and Garnet Kin-Lic Chan, “The chromium dimer: Closing a chapter of quantum chemistry,” Journal of the American Chemical Society144, 15932–15937 (2022)

  34. [34]

    Quasirel- ativistic theory equivalent to fully relativistic 7 theory,

    Werner Kutzelnigg and Wenjian Liu, “Quasirel- ativistic theory equivalent to fully relativistic 7 theory,” The Journal of Chemical Physics123, 241102 (2005), https://pubs.aip.org/aip/jcp/article- pdf/doi/10.1063/1.2137315/13513393/241102 1 online.pdf

  35. [35]

    Infinite-order two- component theory for relativistic quantum chemistry,

    Maria Barysz and Andrzej J. Sadlej, “Infinite-order two- component theory for relativistic quantum chemistry,” The Journal of Chemical Physics116, 2696–2704 (2002)

  36. [36]

    On the spin separation of algebraic two-component relativistic hamiltonians,

    Zhendong Li, Yunlong Xiao, and Wenjian Liu, “On the spin separation of algebraic two-component relativistic hamiltonians,” The Journal of Chemical Physics137, 154114 (2012), https://pubs.aip.org/aip/jcp/article- pdf/doi/10.1063/1.4758987/13274957/154114 1 online.pdf

  37. [37]

    Negative ion photoelectron spectroscopy of chromium dimer,

    Sean M. Casey and Doreen G. Leopold, “Negative ion photoelectron spectroscopy of chromium dimer,” The Journal of Physical Chemistry97, 816–830 (1993)

  38. [38]

    Full-dimensional potentials and state couplings and multidimensional tunneling calculations for the photodissociation of phenol,

    Ke R. Yang, Xuefei Xu, Jingjing Zheng, and Don- ald G. Truhlar, “Full-dimensional potentials and state couplings and multidimensional tunneling calculations for the photodissociation of phenol,” Chem. Sci.5, 4661– 4680 (2014)

  39. [39]

    Photodissociation dynamics of phenol: Multi- state trajectory simulations including tunneling,

    Xuefei Xu, Jingjing Zheng, Ke R. Yang, and Donald G. Truhlar, “Photodissociation dynamics of phenol: Multi- state trajectory simulations including tunneling,” Jour- nal of the American Chemical Society136, 16378–16386 (2014)

  40. [40]

    Combining spin-adapted open-shell td-dft with spin–orbit coupling,

    Zhendong Li, Bingbing Suo, Yong Zhang, Yunlong Xiao, and Wenjian Liu, “Combining spin-adapted open-shell td-dft with spin–orbit coupling,” Molecular Physics111, 3741–3755 (2013)

  41. [41]

    Analytic energy gradients of spin-adapted open- shell time-dependent density functional theory,

    Zikuan Wang, Zhendong Li, Yong Zhang, and Wenjian Liu, “Analytic energy gradients of spin-adapted open- shell time-dependent density functional theory,” The Journal of Chemical Physics153, 164109 (2020)

  42. [42]

    Nac- tddft: Time-dependent density functional theory for nonadiabatic couplings,

    Zikuan Wang, Chenyu Wu, and Wenjian Liu, “Nac- tddft: Time-dependent density functional theory for nonadiabatic couplings,” Accounts of Chemical Research 54, 3288–3297 (2021)

  43. [43]

    NIST Atomic Spectra Database (ver. 5.11),

    A. Kramida, Y. Ralchenko, J. Reader, and NIST ASD Team, “NIST Atomic Spectra Database (ver. 5.11),” (2023)

  44. [44]

    NIST Diatomic Spectral Database (ver. 2.1),

    F. J. Lovas, E. Tiemann, J. S. Coursey, S. A. Ko- tochigova, J. Chang, K. Olsen, and R. A. Dragoset, “NIST Diatomic Spectral Database (ver. 2.1),” (2001)

  45. [45]

    Spec- troscopic identification of the low-lying electronic states of s2 molecule,

    Zhengbo Qin, Lin Wang, Ran Cong, Chengxiang Jiao, Xianfeng Zheng, Zhifeng Cui, and Zichao Tang, “Spec- troscopic identification of the low-lying electronic states of s2 molecule,” The Journal of Chemical Physics150, 044302 (2019). 8 SUPPLEMENT AL MA TERIALS: SPIN-ADAPTED TDDFT I. W orking Equations Due to the non-zero condition ofW(Γ ′ΓSiSi, KSf), we have...