REVIEW 1 major objections 1 minor 1 cited by
Analytic first-order non-adiabatic coupling matrix elements of spin-adapted open-shell time-dependent density functional theory
T0 review · 1 major / 1 minor · reviewed 2026-07-01 · grok-4.3
Pith's one-line read X-TDDFT reduces NACME errors of U-TDDFT by one-third to two-thirds for open-shell molecules
desk verdict X-TDDFT now has analytic NACMEs that reduce U-TDDFT errors by 1/3-2/3 versus multireference values and shift IC rates and pathways in open-shell cases, but the reference convergence is not demonstrated. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The analytic NACME formulas for X-TDDFT that incorporate spin adaptation through the restricted open-shell Kohn-Sham reference and the contribution of doubly excited determinants.
What would settle it
An experimental measurement of internal conversion rates and pathway preferences in copper(II) porphyrin that aligns with the X-TDDFT predictions rather than the U-TDDFT ones.
Extended reading notes
Core claim
X-TDDFT can supply analytic NACMEs for spin-conserving excitations in open-shell molecules. The implementation reuses U-TDDFT code with adjustments for the ROKS reference and implicit double excitations. Benchmarks show that the resulting NACMEs cut the error relative to multireference values by a factor of one-third to two-thirds, which in turn revises internal conversion rates by as much as two orders of magnitude and alters the ranking of excited-state relaxation pathways for copper(II) porphyrin.
Load-bearing premise
The multireference NACME values used as reference are accurate and complete enough to benchmark the two TDDFT variants across the molecules studied.
Editorial extensions
If this is right
- Internal conversion rates receive corrections of up to two orders of magnitude.
- The relative importance of excited state relaxation pathways changes qualitatively for copper(II) porphyrin.
- Substituent effects on the internal conversion rates are revised.
- X-TDDFT NACMEs are useful for photophysics and photochemistry studies of open-shell systems such as radicals and transition metal complexes.
Reading between the lines
- Unrestricted TDDFT errors in NACMEs are large and unsystematic, pointing to the need for spin adaptation in coupling calculations.
- The method opens the door to more reliable non-adiabatic dynamics simulations for open-shell transition metal complexes.
- Similar accuracy gains may apply when computing other derivative properties in spin-adapted TDDFT.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives and implements analytic first-order NACMEs for ground-to-excited and excited-to-excited transitions within the spin-adapted open-shell TDDFT framework (X-TDDFT). The implementation builds on an existing U-TDDFT NACME code by incorporating the ROKS reference and implicit doubly excited determinants. Benchmark results are presented showing that X-TDDFT NACMEs reduce errors relative to U-TDDFT by factors of 1/3–2/3 when compared against multireference reference values; these improvements are reported to produce corrections to internal conversion rates of up to two orders of magnitude and qualitative changes in pathway ordering and substituent effects for copper(II) porphyrin.
Significance. If the numerical claims hold, the work would be significant as the first reported analytic NACME implementation for any spin-adapted TDDFT method. The modest computational overhead and reuse of U-TDDFT infrastructure are practical strengths. The reported impact on IC rates for open-shell systems would be relevant to photophysics studies of radicals and transition-metal complexes.
major comments (1)
- [Benchmark calculations] Benchmark calculations: The central numerical claim—that X-TDDFT reduces U-TDDFT NACME error by 1/3–2/3 versus multireference values, producing order-of-magnitude rate changes and qualitative pathway revisions for Cu(II) porphyrin—rests on the assumption that the multireference NACME references are converged and more accurate than the DFT results. No explicit convergence tests (active-space enlargement, basis-set extrapolation, state-averaging sensitivity, or cross-validation against a second multireference method) are described for the reference values. This is load-bearing for the error-reduction percentages and downstream rate revisions.
minor comments (1)
- [Abstract] The abstract states the numerical improvements but does not reference the explicit NACME working equations or the benchmark data tables; readers must reach the main text to locate these.
Simulated Author's Rebuttal
We thank the referee for the careful reading of the manuscript and the constructive comment on the benchmark section. We address the point below.
read point-by-point responses
-
Referee: [Benchmark calculations] Benchmark calculations: The central numerical claim—that X-TDDFT reduces U-TDDFT NACME error by 1/3–2/3 versus multireference values, producing order-of-magnitude rate changes and qualitative pathway revisions for Cu(II) porphyrin—rests on the assumption that the multireference NACME references are converged and more accurate than the DFT results. No explicit convergence tests (active-space enlargement, basis-set extrapolation, state-averaging sensitivity, or cross-validation against a second multireference method) are described for the reference values. This is load-bearing for the error-reduction percentages and downstream rate revisions.
Authors: We agree that explicit convergence tests for the multireference NACME reference values are not described in the manuscript and that this information is important for supporting the reported error reductions. In the revised manuscript we will add a dedicated subsection presenting convergence checks with respect to active-space size, basis-set quality, state averaging, and cross-validation against a second multireference method for the key benchmark systems, including copper(II) porphyrin. These additional data will be used to confirm that the reference values are sufficiently converged for the comparisons presented. revision: yes
Circularity Check
No circularity in NACME derivation chain
full rationale
The paper presents an analytic derivation of ground-to-excited and excited-to-excited NACMEs for X-TDDFT, implemented as an extension of an existing U-TDDFT code that incorporates the ROKS reference and implicit doubly excited determinants. No load-bearing step reduces the NACME expressions to fitted parameters, self-citations, or prior ansatzes by construction; the benchmark comparisons to external multireference values serve as independent validation rather than definitional inputs. The derivation chain is therefore self-contained.
Assumptions & free parameters
Cite this review
Pith. "Pith review of Analytic first-order non-adiabatic coupling matrix elements of spin-adapted open-shell time-dependent density functional theory." pith.science (2026). https://pith.science/paper/JAGGSJLI
@misc{pith2026260526594,
author = {Pith},
title = {Pith review of: Analytic first-order non-adiabatic coupling matrix elements of spin-adapted open-shell time-dependent density functional theory},
year = {2026},
howpublished = {\url{https://pith.science/paper/JAGGSJLI}},
note = {Machine review of arXiv:2605.26594}
}
read the original abstract
While spin-adapted time-dependent density functional theory (TDDFT) approaches significantly improve the excitation energies and gradients of open-shell molecules, the effect of spin-adaptation on non-adiabatic coupling matrix elements (NACMEs) remains unknown for spin-conserving excitations. In this article, we report the derivation, implementation and benchmark studies of the ground state-excited state and excited state-excited state NACMEs of our spin-adapted TDDFT method, X-TDDFT; to our best knowledge, this represents the first implementation of the analytic NACMEs of a spin-adapted TDDFT method. Similar to the X-TDDFT analytic gradients, X-TDDFT NACMEs can be easily implemented on top of an existing U-TDDFT NACME implementation taking into account the restricted open-shell Kohn-Sham (ROKS) reference and the implicit involvement of doubly excited determinants, with acceptable computational overhead. Benchmark calculations reveal that X-TDDFT reduces the error of U-TDDFT NACMEs by 1/3-2/3 (referenced against high-level multireference NACMEs), which leads to large corrections of internal conversion rates (up to two orders of magnitude). In particular, for copper(II) porphyrin, X-TDDFT leads to qualitative revisions of the relative importance of the excited state relaxation pathways, as well as the substituent effects of the internal conversion (IC) rates, suggesting that the error of U-TDDFT NACMEs is not only large but also unsystematic. It is therefore expected that X-TDDFT NACMEs will prove useful in the photophysics/photochemistry studies of open-shell systems such as radicals and transition metal complexes.
Figures
Figures from the paper (4 more)
Forward citations
Cited by 1 Pith paper
-
Spin-Consistency Constraints in Noncollinear Tensor TDA
A two-step kernel reconstruction, derived from spin-consistency constraints, enforces spin purity and removes artifact states in noncollinear tensor TDA for open-shell references.
Reference graph
Works this paper leans on
-
[1]
J.; Gu, Q.; Ding, J.; Chen, Z.; Hele, T
(1) Li, F.; Gillett, A. J.; Gu, Q.; Ding, J.; Chen, Z.; Hele, T. J.; Myers, W. K.; Friend, R. H.; Evans, E. W. Singlet and triplet to doublet energy transfer: improving organic light- emitting diodes with radicals.Nat. Commun.2022,13,
work page 2022
-
[2]
(2) Pfund, B.; Wenger, O. S. Excited Organic Radicals in Photoredox Catalysis.JACS Au 2025,5, 426–447. (3) Huang, T.; Du, P.; Lin, Y.-M. Recent Advances in Photoactive First-Row Transition Metal Complexes for Organic Synthesis.Chin. J. Chem.2025,43, 2566–2587. (4) Runge, E.; Gross, E. K. Density-functional theory for time-dependent systems.Phys. Rev. Lett...
work page 2025
-
[3]
(5) Casida, M. E. InRecent Advances in Density Functional Methods: (Part I); Chong, D. P., Ed.; WORLD SCIENTIFIC, 1995; pp 155–192. (6) Chattopadhyay, S.; Chaudhuri, R. K.; Mahapatra, U. S.; Ghosh, A.; Ray, S. S. State- specific multireference perturbation theory: development and present status.WIREs Comput. Mol. Sci.2016,6, 266–291. (7) Song, Y.; Guo, Y....
work page 1995
-
[4]
BDF: A relativistic electronic structure program package.J
(33) Zhang, Y.; Suo, B.; Wang, Z.; Zhang, N.; Li, Z.; Lei, Y.; Zou, W.; Gao, J.; Peng, D.; Pu, Z.; Xiao, Y.; Sun, Q.; Wang, F.; Ma, Y.; Wang, X.; Guo, Y.; Liu, W. BDF: A relativistic electronic structure program package.J. Chem. Phys.2020,152, 064113. (34) Shiozaki, T. BAGEL: brilliantly advanced general electronic-structure library.WIREs Comput. Mol. Sci...
work page 2020
-
[5]
Toward reliable density functional methods without adjustable parameters: The PBE0 model.J
(42) Adamo, C.; Barone, V. Toward reliable density functional methods without adjustable parameters: The PBE0 model.J. Chem. Phys.1999,110, 6158–6170. (43) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu.J. Chem. Phys.2010,132, 1...
work page 1999
-
[6]
(58) Lei, Y.; Liu, W.; Hoffmann, M. R. Further development of SDSPT2 for strongly corre- lated electrons.Mol. Phys.2017,115, 2696–2707. (59) Lei,Y.; Suo,B.; Liu,W.iCAS:imposedautomaticselectionandlocalizationofcomplete active spaces.J. Chem. Theory Comput.2021,17, 4846–4859. 60 (60) Guo, Y.; Zhang, N.; Lei, Y.; Liu, W. iCISCF: an iterative configuration i...
work page 2017
-
[7]
(62) Zhou, Y.; Wang, Z.; Wu, C.; Li, J.; He, L.; Zhang, B.; Guo, J.; Zhang, J. Inhibition of Imidacloprid Photolysis by Particulate Matter: Experimental and Computational Evidence.Environ. Sci. Technol.2026,60, 8681–8691. (63) Humphrey, W.; Dalke, A.; Schulten, K. VMD: visual molecular dynamics.J. Mol. Graphics1996,14, 33–38. (64) Lu,T.; Chen,F.Multiwfn: ...
work page 2026
Reviewed July 1, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.