pith. sign in

arxiv: 2603.29948 · v2 · submitted 2026-03-31 · ⚛️ physics.chem-ph

TDDeltaSCF: Time-Dependent Density Functional Theory with a Non-Aufbau Reference for near-degenerate states

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

classification ⚛️ physics.chem-ph
keywords TDDFTDeltaSCFnon-Aufbau referencenear-degenerate statesdiradicalssinglet-triplet gapsbond dissociation
0
0 comments X

The pith

TDΔSCF uses a non-Aufbau ΔSCF reference for TDDFT to treat near-degenerate electronic states while keeping standard response contributions.

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

The paper proposes TDΔSCF as a linear-response method that starts from a non-Aufbau ΔSCF determinant instead of the usual ground-state reference. This lets the calculation target states with near-degeneracy from an already excited configuration, preserving the ordinary Coulomb and exchange-correlation kernel terms that collinear spin-flip TDDFT drops. Tests on ethylene torsion, diradical singlet-triplet gaps, benzyne isomer geometries, and bond-dissociation curves show reduced sensitivity to functional choice and fewer spurious low-lying states than spin-flip TDDFT. The approach works best when the chosen ΔSCF reference closely matches the character of the final target state. A reader would care because it offers a relatively cheap route to balanced singlet energies in systems where conventional single-reference DFT breaks down.

Core claim

TDΔSCF is a linear-response scheme in which a non-Aufbau ΔSCF determinant serves as the reference for a subsequent TDDFT calculation; in contrast to collinear spin-flip TDDFT, this formulation preserves the usual Coulomb and exchange-correlation response contributions while describing the target states from an electronically excited reference.

What carries the argument

TDΔSCF method: non-Aufbau ΔSCF reference combined with standard TDDFT linear response to generate target states from an excited configuration.

If this is right

  • Torsional potentials become smooth without artificial barriers or discontinuities.
  • Singlet-triplet gaps improve and become more balanced for representative diradicals.
  • Singlet m-benzyne consistently optimizes to the monocyclic geometry rather than ring-opening.
  • Bond-dissociation curves for HF and F2 lack the spurious low-lying states seen in SF-TDDFT.
  • Functional dependence is markedly weaker across the tested systems.

Where Pith is reading between the lines

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

  • The method could be combined with orbital-optimized or range-separated functionals to further reduce singlet overestimation when the reference quality is marginal.
  • Numerical instabilities traced to exchange-correlation potentials near nodal regions suggest that grid or basis-set refinements may be needed for routine use on larger molecules.
  • Because the reference is chosen by hand, the approach invites automated selection protocols that scan multiple non-Aufbau occupations before the TDDFT step.

Load-bearing premise

The non-Aufbau ΔSCF reference must stay well suited to the final target state; accuracy drops when this match fails.

What would settle it

A test case in which the ΔSCF reference produces a qualitatively wrong state ordering or energy errors larger than those of high-level reference methods would show the central claim does not hold.

Figures

Figures reproduced from arXiv: 2603.29948 by Fumiya Mohri, Shuto Shibasaki, Takashi Tsuchimochi.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8 [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9 [PITH_FULL_IMAGE:figures/full_fig_p014_9.png] view at source ↗
read the original abstract

Near-degenerate electronic structures remain a major challenge for conventional single-reference density functional theory (DFT). To address this problem, we propose time-dependent $\Delta$SCF (TD$\Delta$SCF), a novel linear-response scheme in which a non-Aufbau $\Delta$SCF determinant serves as the reference for a subsequent TDDFT calculation. In contrast to collinear spin-flip (SF)-TDDFT, this formulation preserves the usual Coulomb and exchange-correlation response contributions while describing the target states from an electronically excited reference. We examine the performance of TD$\Delta$SCF for several prototypical problems involving near-degeneracy, including the torsional potential of ethylene, singlet--triplet gaps of representative diradicals, geometry optimizations of benzyne isomers, and bond-dissociation curves of hydrogen fluoride and F$_2$. Across these tests, TD$\Delta$SCF shows markedly weaker functional dependence than SF-TDDFT and often yields a more balanced description of challenging singlet states. In particular, it provides smooth torsional potentials, improved singlet--triplet gaps, a consistent monocyclic structure for singlet $m$-benzyne, and a more satisfactory description of bond dissociation without the spurious low-lying states found in SF-TDDFT. At the same time, the method exhibits a systematic tendency to overestimate singlet energies and can lose accuracy when the underlying $\Delta$SCF reference is not well suited to the final state. We also identify a numerical instability that can arise in non-Aufbau calculations and trace its origin to the exchange-correlation potential near uncompensated nodal regions. These results highlight both the promise and the practical limitations of TD$\Delta$SCF as a low-cost method for singlet states with near-degenerate electronic structures.

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 introduces TDΔSCF, a linear-response TDDFT formulation that employs a non-Aufbau ΔSCF determinant as the reference state for near-degenerate electronic structures. Unlike collinear spin-flip TDDFT, the approach retains the standard Coulomb and exchange-correlation response kernels. Performance is assessed on the ethylene torsional potential, singlet-triplet gaps in diradicals, benzyne isomer geometries, and HF/F2 bond-dissociation curves, with claims of reduced functional dependence, smoother potentials, improved gaps, and avoidance of spurious states relative to SF-TDDFT, while noting a tendency to overestimate singlet energies and loss of accuracy when the reference determinant is unsuitable.

Significance. If the central claims hold after detailed verification, TDΔSCF would constitute a practical, low-cost extension of standard TDDFT for singlet states in near-degenerate regimes, offering an alternative to spin-flip methods with less functional sensitivity. The preservation of unmodified response kernels is a conceptually attractive feature that could facilitate broader adoption in computational studies of diradicals, bond breaking, and photochemical processes.

major comments (3)
  1. [Abstract] Abstract: the assertion that TDΔSCF 'preserves the usual Coulomb and exchange-correlation response contributions' is central to the method's novelty yet is presented only qualitatively; no explicit equations, kernel definitions, or derivation showing equivalence to standard TDDFT response are supplied, leaving the claim unverified from the given text.
  2. [Abstract] Abstract: numerical performance is reported on multiple test systems but without error bars, convergence thresholds, basis-set details, or tabulated quantitative metrics (e.g., MAE or RMSD versus reference values), rendering the statements of 'markedly weaker functional dependence' and 'more balanced description' difficult to evaluate rigorously.
  3. [Abstract] Abstract: the critical assumption that the non-Aufbau ΔSCF reference remains 'well suited to the final target state' is acknowledged as a limitation, but no diagnostic criterion, orbital-occupation check, or quantitative measure of reference-target deviation is provided to determine when the linear-response correction can be trusted.
minor comments (2)
  1. The abstract lists specific test cases but omits the exchange-correlation functionals and basis sets employed in the comparisons; inclusion of these details would improve reproducibility.
  2. A summary table comparing TDΔSCF, SF-TDDFT, and reference values across all reported quantities would enhance clarity of the performance claims.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their thoughtful and constructive review. We address each major comment below and have revised the manuscript to improve clarity, quantitative rigor, and practical guidance.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the assertion that TDΔSCF 'preserves the usual Coulomb and exchange-correlation response contributions' is central to the method's novelty yet is presented only qualitatively; no explicit equations, kernel definitions, or derivation showing equivalence to standard TDDFT response are supplied, leaving the claim unverified from the given text.

    Authors: We agree the abstract is too terse on this central point. The full manuscript (Section II) derives the TDΔSCF working equations from the standard TDDFT linear-response formalism applied to a non-Aufbau reference determinant; because the reference enters only through the ground-state density and orbitals, the Coulomb and XC kernels remain exactly the usual expressions. We have revised the abstract to include a concise statement of this equivalence and added a cross-reference to the explicit derivation. revision: yes

  2. Referee: [Abstract] Abstract: numerical performance is reported on multiple test systems but without error bars, convergence thresholds, basis-set details, or tabulated quantitative metrics (e.g., MAE or RMSD versus reference values), rendering the statements of 'markedly weaker functional dependence' and 'more balanced description' difficult to evaluate rigorously.

    Authors: We acknowledge that the abstract lacks the requested quantitative summary. The main text already specifies basis sets, convergence criteria, and reports raw energies and geometries for each system. To strengthen the abstract, we have added a compact table of MAE/RMSD values across the test set and inserted error bars (where replicate calculations exist) into the revised figures. These additions make the claims of reduced functional dependence directly verifiable. revision: partial

  3. Referee: [Abstract] Abstract: the critical assumption that the non-Aufbau ΔSCF reference remains 'well suited to the final target state' is acknowledged as a limitation, but no diagnostic criterion, orbital-occupation check, or quantitative measure of reference-target deviation is provided to determine when the linear-response correction can be trusted.

    Authors: This is a fair and important practical concern. We have expanded the discussion section to include a short set of recommended checks: (i) inspection of the final ΔSCF orbital occupations relative to the target state, (ii) monitoring of the magnitude of the TDΔSCF energy correction, and (iii) comparison of the reference density to a known benchmark when available. These diagnostics are now summarized in the abstract as well. revision: yes

Circularity Check

0 steps flagged

No circularity in TDΔSCF formulation or derivation chain

full rationale

The paper introduces TDΔSCF as a novel linear-response TDDFT scheme that uses a non-Aufbau ΔSCF determinant as reference while retaining standard Coulomb and XC kernels. This is presented as a direct methodological extension, not a reduction of outputs to inputs by construction. Performance claims rest on independent test cases (ethylene torsion, diradical gaps, benzyne geometries, HF/F2 dissociation) rather than fitted parameters or self-citation chains. No self-definitional steps, fitted-input predictions, load-bearing self-citations, or smuggled ansatzes appear in the derivation; the acknowledged limitation (accuracy loss when the reference is unsuitable) is an explicit practical caveat, not a hidden circularity. The method is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central construction rests on the domain assumption that a non-Aufbau determinant can serve as a stable reference for standard TDDFT response kernels; no free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption A non-Aufbau ΔSCF determinant can be used as a valid reference for conventional TDDFT linear response without altering the Coulomb and XC kernels.
    Stated directly in the abstract as the key difference from SF-TDDFT.

pith-pipeline@v0.9.0 · 5633 in / 1309 out tokens · 24850 ms · 2026-05-13T23:15:59.935223+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

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

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.

Reference graph

Works this paper leans on

49 extracted references · 49 canonical work pages · 1 internal anchor

  1. [1]

    Kohn and L

    W. Kohn and L. J. Sham, Physical Review 140, A1133 (1965)

  2. [2]

    Runge and E

    E. Runge and E. K. U. Gross, Phys. Rev. Lett. 52, 997 (1984). 15

  3. [3]

    Casida and M

    M. Casida and M. Huix-Rotllant, Annual Review of Physical Chemistry 63, 287 (2012)

  4. [4]

    Y. Shao, M. Head-Gordon, and A. I. Krylov, The Journal of Chemical Physics 118, 4807 (2003)

  5. [5]

    A. I. Krylov, Accounts of Chemical Research 39, 83 (2006)

  6. [6]

    Minezawa and M

    N. Minezawa and M. S. Gordon, The Journal of Phys- ical Chemistry A 113, 12749 (2009) , pMID: 19905013, https://doi.org/10.1021/jp908032x

  7. [7]

    Harabuchi, K

    Y. Harabuchi, K. Keipert, F. Zahariev, T. Taket- sugu, and M. S. Gordon, The Journal of Physi- cal Chemistry A 118, 11987 (2014) , pMID: 25181251, https://doi.org/10.1021/jp5072428

  8. [8]

    L. Yue, Y. Liu, and C. Zhu, Phys. Chem. Chem. Phys. 20, 24123 (2018)

  9. [9]

    Li and W

    Z. Li and W. Liu, The Journal of Chemical Physics 136, 024107 (2012)

  10. [10]

    Y. A. Bernard, Y. Shao, and A. I. Krylov, The Journal of Chemical Physics 136, 204103 (2012)

  11. [11]

    Huix-Rotllant, B

    M. Huix-Rotllant, B. Natarajan, A. Ipatov, C. Muhavini Wawire, T. Deutsch, and M. E. Casida, Phys. Chem. Chem. Phys. 12, 12811 (2010)

  12. [12]

    A. T. Gilbert, N. A. Besley, and P. M. Gill, J. Phys. Chem. A 112, 13164 (2008)

  13. [13]

    Filatov and S

    M. Filatov and S. Shaik, The Journal of Chemical Physics 110, 116 (1999)

  14. [14]

    Kowalczyk, T

    T. Kowalczyk, T. Tsuchimochi, P.-T. Chen, L. Top, and T. Van Voorhis, J. Chem. Phys. 138, 164101 (2013)

  15. [15]

    Hait and M

    D. Hait and M. Head-Gordon, The Journal of Physi- cal Chemistry Letters 11, 775 (2020) , pMID: 31917579, https://doi.org/10.1021/acs.jpclett.9b03661

  16. [16]

    Hait and M

    D. Hait and M. Head-Gordon, J. Phys. Chem. Lett. 12, 4517 (2021)

  17. [17]

    Selenius, A

    E. Selenius, A. E. Sigurdarson, Y. L. A. Schmerwitz, and G. Levi, 20, 3809 (2024)

  18. [18]

    Paetow and J

    L. Paetow and J. Neugebauer, Phys. Chem. Chem. Phys. 27, 16354 (2025)

  19. [19]

    Carter-Fenk and J

    K. Carter-Fenk and J. M. Herbert, Journal of Chemi- cal Theory and Computation 16, 5067 (2020) , pMID: 32644792, https://doi.org/10.1021/acs.jctc.0c00502

  20. [20]

    Hait and M

    D. Hait and M. Head-Gordon, Journal of Chemical The- ory and Computation 16, 1699 (2020) , pMID: 32017554, https://doi.org/10.1021/acs.jctc.9b01127

  21. [21]

    G. Levi, A. V. Ivanov, and H. Jónsson, Faraday Discuss. 224, 448 (2020)

  22. [22]

    G. Levi, A. V. Ivanov, and H. Jnsson, J. Chem. Theory Comput. 16, 6968 (2020)

  23. [23]

    Berera, R

    R. Berera, R. van Grondelle, and J. T. M. Kennis, Pho- tosynthesis Research 101, 105 (2009)

  24. [24]

    Z. J. Knepp, D. R. Fertal, G. B. Masso, R. C. Hamburger, C. A. Guzman, E. R. Young, and L. A. Fredin, Journal of Chemical Theory and Computation 21, 9736 (2025) , pMID: 41031514, https://doi.org/10.1021/acs.jctc.5c00591

  25. [25]

    and Arias-Martinez, Juan E

    K. Carter-Fenk, L. A. Cunha, J. E. Arias-Martinez, and M. Head-Gordon, The Journal of Physical Chem- istry Letters 13, 9664 (2022) , pMID: 36215404, https://doi.org/10.1021/acs.jpclett.2c02564

  26. [26]

    Tsuchimochi, The Journal of Chemical Physics 161, 241102 (2024)

    T. Tsuchimochi, The Journal of Chemical Physics 161, 241102 (2024)

  27. [27]

    Tsuchimochi and B

    T. Tsuchimochi and B. Mokhtar, Journal of Chemical Theory and Computation 0, null (0) , https://doi.org/10.1021/acs.jctc.6c00182

  28. [28]

    Analytical Nuclear Gradients for State-Averaged Configuration Interaction Singles Variants: Application to Conical Intersections

    T. Tsuchimochi, “Analytical nuclear gradients for state-averaged configuration interaction singles vari- ants: Application to conical intersections,” (2026), arXiv:2602.15700 [physics.chem-ph]

  29. [29]

    R. G. Parr and Y. Weitao, Density-Functional Theory of Atoms and Molecules (Oxford University Press, 1995)

  30. [30]

    C. A. Ullrich, Time-Dependent Density-Functional Theory: Concepts and Applications (Oxford University Press, 2011)

  31. [31]

    Petersilka, U

    M. Petersilka, U. J. Gossmann, and E. K. U. Gross, Phys. Rev. Lett. 76, 1212 (1996)

  32. [32]

    Bauernschmitt and R

    R. Bauernschmitt and R. Ahlrichs, Chemical Physics Let- ters 256, 454 (1996)

  33. [33]

    Hirata and M

    S. Hirata and M. Head-Gordon, Chemical Physics Letters 314, 291 (1999)

  34. [34]

    Wang and T

    F. Wang and T. Ziegler, The Journal of Chemical Physics 121, 12191 (2004)

  35. [35]

    Wang and T

    F. Wang and T. Ziegler, The Journal of Chemical Physics 122, 074109 (2005)

  36. [36]

    A. D. Becke, The Journal of Chemical Physics 98, 5648 (1993)

  37. [37]

    Epifanovsky, A

    E. Epifanovsky, A. T. B. Gilbert, X. Feng, J. Lee, Y. Mao, N. Mardirossian, P. Pokhilko, A. F. White, M. P. Coons, A. L. Dempwolff, Z. Gan, D. Hait, P. R. Horn, L. D. Jacobson, I. Kaliman, J. Kussmann, A. W. Lange, K. U. Lao, D. S. Levine, J. Liu, S. C. McKen- zie, A. F. Morrison, K. D. Nanda, F. Plasser, D. R. Rehn, M. L. Vidal, Z.-Q. You, Y. Zhu, B. Ala...

  38. [38]

    Q. Sun, X. Zhang, S. Banerjee, P. Bao, M. Barbry, N. S. Blunt, N. A. Bogdanov, G. H. Booth, J. Chen, Z.-H. Cui, J. J. Eriksen, Y. Gao, S. Guo, J. Hermann, M. R. Hermes, K. Koh, P. Koval, S. Lehtola, Z. Li, J. Liu, N. Mardirossian, J. D. McClain, M. Motta, B. Mus- sard, H. Q. Pham, A. Pulkin, W. Purwanto, P. J. Robin- son, E. Ronca, E. R. Sayfutyarova, M. ...

  39. [39]

    Saade and H

    S. Saade and H. G. A. Burton, Journal of Chemical The- ory and Computation 20, 5105 (2024) , pMID: 38847452

  40. [40]

    A. I. Krylov, Chemical Physics Letters 338, 375 (2001)

  41. [41]

    Nikiforov, J

    A. Nikiforov, J. A. Gamez, W. Thiel, M. Huix-Rotllant, and M. Filatov, The Journal of Chemical Physics 141, 124122 (2014)

  42. [42]

    K. P. Huber and G. Herzberg, Constants of Diatomic Molecules (Van Nostrand Rein- hold, New York, 1979)

  43. [43]

    P. G. Wenthold, R. R. Squires, and W. C. Lineberger, Journal of the American Chemical Society 120, 5279 (1998)

  44. [44]

    Furche and R

    F. Furche and R. Ahlrichs, The Journal of Chemical Physics 117, 7433 (2002)

  45. [45]

    Kraka, J

    E. Kraka, J. Anglada, A. Hjerpe, M. Filatov, and D. Cre- mer, Chemical Physics Letters 348, 115 (2001)

  46. [46]

    Winkler and W

    M. Winkler and W. Sander, The Journal of Physical Chemistry A 105, 10422 (2001)

  47. [47]

    W. A. Al-Saidi and C. J. Umrigar, The Journal of Chem- ical Physics 128, 154324 (2008)

  48. [48]

    W. T. Zemke, W. C. Stwalley, J. A. Coxon, and P. G. Hajigeorgiou, Chemical Physics Letters 177, 412 (1991)

  49. [49]

    P. Wang, S. Gong, and Y. Mo, The Journal of Physical Chemistry Letters 15, 12594 (2024)