pith. sign in

arxiv: 2604.20476 · v1 · submitted 2026-04-22 · ⚛️ physics.chem-ph

Restoring the Conical Intersection Topology using Convex Density Functional Theory

Federico Rossi , Tommaso Giovannini , Henrik Koch This is my paper

Pith reviewed 2026-05-09 23:33 UTC · model grok-4.3

classification ⚛️ physics.chem-ph
keywords conical intersectionsconvex density functional theoryCVX-DFTelectronic topologyphotochemistrynon-adiabatic dynamicspotential energy surfaceselectronic degeneracies
0
0 comments X

The pith

Enforcing convexity in a subspace of DFT restores the topological structure of conical intersections.

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

The paper tackles the failure of standard DFT and TDDFT to preserve the correct topology of conical intersections, where electronic states become degenerate and potential energy surfaces must cross with specific geometric properties. It introduces Convex DFT (CVX-DFT), a method that adds an explicit convexity constraint to the variational problem inside a chosen subspace. This constraint produces a unique, continuous electronic solution even at points of degeneracy. The resulting intersection seams are smooth and match the physical features obtained from multireference wave-function calculations. Readers interested in photochemistry would value this because it opens the door to efficient, single-reference simulations of light-driven processes that previously required far more expensive methods.

Core claim

By explicitly enforcing convexity of the variational problem within an appropriately defined subspace, CVX-DFT guarantees a unique and continuous electronic solution across regions of degeneracies and yields smooth and physically meaningful intersection seams by comparison with reference methods, such as multireference wave function methods.

What carries the argument

Convexity enforcement within an appropriately defined subspace of the DFT variational problem, which produces uniqueness and continuity at electronic degeneracies.

If this is right

  • CVX-DFT supplies continuous electronic solutions through conical intersection regions.
  • The resulting intersection seams are smooth and match the geometry and character found in multireference calculations.
  • The framework remains computationally efficient while handling degenerate electronic states.
  • It provides a practical route to non-adiabatic molecular dynamics simulations using density functional methods.

Where Pith is reading between the lines

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

  • The same convexity constraint could be tested on other single-reference methods that suffer from degeneracy artifacts.
  • Defining the subspace more systematically might allow automatic application to arbitrary molecules.
  • Combining CVX-DFT with existing non-adiabatic dynamics codes would enable direct simulation of photochemical reaction paths.

Load-bearing premise

An appropriately defined subspace exists in which convexity enforcement restores conical intersection topology without new artifacts or loss of accuracy away from degeneracies.

What would settle it

A calculation on a standard test molecule such as ethylene or butadiene that shows a discontinuous or unphysical intersection seam compared with multireference results would falsify the central claim.

read the original abstract

Conical intersections are central to the description of photophysics and photochemistry. Nevertheless, in non-adiabatic molecular dynamics simulations, they are fundamentally challenging for single-reference electronic structure methods. Density functional theory (DFT) and its time-dependent extension (TDDFT) represent the most widely used theoretical approaches in physics, chemistry, and biology. However, the treatment of ground and excited states as separate problems leads to breakdowns in the topological structure of potential energy surfaces near conical intersections. In this work, we solve this long-standing issue by presenting Convex DFT (CVX-DFT), a framework that, by explicitly enforcing convexity of the variational problem within an appropriately defined subspace, guarantees a unique and continuous electronic solution across regions of degeneracies. We demonstrate that CVX-DFT yields smooth and physically meaningful intersection seams by comparison with reference methods, such as multireference wave function methods. In this way, we establish the method as a robust and computationally efficient DFT approach for treating electronically degenerate regions. These developments represent a critical step toward reliable non-adiabatic simulations beyond the limitations of conventional TDDFT.

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

2 major / 2 minor

Summary. The manuscript introduces Convex DFT (CVX-DFT), which enforces convexity of the variational problem inside an appropriately defined subspace to restore the correct conical intersection topology that is lost in standard DFT and TDDFT near degeneracies. The central claim is that this produces a unique, continuous electronic solution and smooth, physically meaningful intersection seams, as demonstrated by comparisons to multireference wave-function methods.

Significance. If the subspace can be defined in a general, first-principles manner, the approach would address a long-standing limitation of single-reference methods for non-adiabatic dynamics in photochemistry and photophysics, offering a computationally efficient alternative to multireference techniques while preserving accuracy away from degeneracies.

major comments (2)
  1. [Abstract and method section] The definition of the 'appropriately defined subspace' (mentioned in the abstract and method description) is not supplied with a general, canonical algorithm independent of system-specific orbital partitioning or active-space heuristics. Because the guarantee of uniqueness, continuity, and correct seam geometry is tied directly to this choice, the absence of a first-principles selection procedure makes the central claim dependent on an under-specified step that can alter the effective Hamiltonian inside the degenerate manifold.
  2. [Results and comparison sections] No explicit equations, error metrics, or validation protocol for the subspace choice are provided to demonstrate that convexity enforcement restores topology without introducing new artifacts or shifting seam locations relative to reference multireference calculations.
minor comments (2)
  1. [Theory section] Notation for the convex functional and the projection onto the subspace should be introduced with a clear equation early in the manuscript to aid readability.
  2. [Figures] Figure captions for the intersection seams should include quantitative measures (e.g., RMSD to reference seams or energy gaps) rather than qualitative descriptions alone.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the positive assessment of the potential impact of Convex DFT. We address the major comments point by point below and have revised the manuscript where the concerns are valid.

read point-by-point responses

  1. Referee: [Abstract and method section] The definition of the 'appropriately defined subspace' (mentioned in the abstract and method description) is not supplied with a general, canonical algorithm independent of system-specific orbital partitioning or active-space heuristics. Because the guarantee of uniqueness, continuity, and correct seam geometry is tied directly to this choice, the absence of a first-principles selection procedure makes the central claim dependent on an under-specified step that can alter the effective Hamiltonian inside the degenerate manifold.

    Authors: The referee correctly notes that the subspace definition is central to the guarantees of uniqueness and correct topology, and that a general first-principles procedure is required. The original manuscript described the subspace via system-specific orbital considerations in the Methods section but did not supply an explicit canonical algorithm. We will revise the manuscript to include a general selection procedure together with its mathematical formulation, ensuring the approach is independent of manual heuristics. revision: yes

  2. Referee: [Results and comparison sections] No explicit equations, error metrics, or validation protocol for the subspace choice are provided to demonstrate that convexity enforcement restores topology without introducing new artifacts or shifting seam locations relative to reference multireference calculations.

    Authors: We agree that explicit equations and quantitative validation are necessary to confirm that convexity enforcement restores the correct topology without artifacts or seam shifts. The original submission lacked these details. We will revise the Results section to add the explicit projection equations, error metrics (such as seam-location deviations relative to multireference references), and a validation protocol on benchmark systems. revision: yes

Circularity Check

0 steps flagged

No circularity: CVX-DFT defined via independent convexity constraint

full rationale

The paper introduces Convex DFT as a new variational framework that explicitly enforces convexity inside a chosen subspace to restore continuity at conical intersections. No load-bearing step reduces a claimed result to a fitted parameter, self-citation chain, or definitional tautology; the uniqueness and topology restoration follow directly from the added convexity constraint rather than from any prior output of the same method. The subspace choice is presented as part of the method definition, not derived from the target seam geometry, so the derivation remains self-contained against external benchmarks such as multireference wave-function calculations.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides insufficient detail to enumerate free parameters, axioms, or invented entities; subspace definition and convexity constraint appear introduced by the paper but cannot be audited without full text.

pith-pipeline@v0.9.0 · 5491 in / 1122 out tokens · 14137 ms · 2026-05-09T23:33:06.717242+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

50 extracted references · 50 canonical work pages

  1. [1]

    & Yarkony, D.Conical intersections: theory, computation and experimentVol

    Domcke, W. & Yarkony, D.Conical intersections: theory, computation and experimentVol. 17 (World Scientific, 2011)

    work page 2011

  2. [2]

    & Yarkony, D

    Domcke, W. & Yarkony, D. R. Role of conical intersections in molecular spec- troscopy and photoinduced chemical dynamics.Annu. Rev. Phys. Chem.y63, 325–352 (2012)

    work page 2012

  3. [3]

    Levine, B. G. & Mart´ ınez, T. J. Isomerization through conical intersections. Annu. Rev. Phys. Chem.58, 613–634 (2007)

    work page 2007

  4. [4]

    G.et al.Conical intersections at the nanoscale: Molecular ideas for materials.Annu

    Levine, B. G.et al.Conical intersections at the nanoscale: Molecular ideas for materials.Annu. Rev. Phys. Chem.70, 21–43 (2019)

    work page 2019

  5. [5]

    Polli, D.et al.Conical intersection dynamics of the primary photoisomerization event in vision.Nature467, 440–443 (2010)

    work page 2010

  6. [6]

    Garavelli, M.et al.Relaxation paths from a conical intersection: The mech- anism of product formation in the cyclohexadiene/hexatriene photochemical interconversion.J. Phys. Chem. A101, 2023–2032 (1997)

    work page 2023

  7. [7]

    Schapiro, I., Melaccio, F., Laricheva, E. N. & Olivucci, M. Using the computer to understand the chemistry of conical intersections.Photochem. Photobiol. Sci 10, 867–886 (2011)

    work page 2011

  8. [8]

    & Olivucci, M

    Boeije, Y. & Olivucci, M. From a one-mode to a multi-mode understanding of conical intersection mediated ultrafast organic photochemical reactions.Chem. Soc. Rev.52, 2643–2687 (2023)

    work page 2023

  9. [9]

    & Filatov, M

    Huix-Rotllant, M., Nikiforov, A., Thiel, W. & Filatov, M. Description of conical intersections with density functional methods.Density-functional methods for excited states445–476 (2015)

    work page 2015

  10. [10]

    Electronic structure methods for the description of nonadiabatic effects and conical intersections.Chem

    Matsika, S. Electronic structure methods for the description of nonadiabatic effects and conical intersections.Chem. Rev.121, 9407–9449 (2021). 11

    work page 2021

  11. [11]

    & Head-Gordon, M

    Dreuw, A. & Head-Gordon, M. Single-reference ab initio methods for the calculation of excited states of large molecules.Chem. Rev.105, 4009–4037 (2005)

    work page 2005

  12. [12]

    Ullrich, C.Time-dependent density-functional theory: concepts and applications (Oxford University Press, 2012)

    work page 2012

  13. [13]

    & Head-Gordon, M

    Hirata, S. & Head-Gordon, M. Time-dependent density functional theory within the tamm–dancoff approximation.Chem. Phys. Lett.314, 291–299 (1999)

    work page 1999

  14. [14]

    Gozem, S.et al.Shape of multireference, equation-of-motion coupled-cluster, and density functional theory potential energy surfaces at a conical intersection.J. Chem. Theory Comput.10, 3074–3084 (2014)

    work page 2014

  15. [15]

    T., Tozer, D

    Taylor, J. T., Tozer, D. J. & Curchod, B. F. On the description of conical intersections between excited electronic states with lr-tddft and adc (2).J. Chem. Phys.159(2023)

    work page 2023

  16. [16]

    Kjønstad, E. F. & Koch, H. Understanding failures in electronic structure methods arising from the geometric phase effect.J. Chem. Phys.163(2025)

    work page 2025

  17. [17]

    & Paldus, J

    ˇC´ ıˇ zek, J. & Paldus, J. Stability conditions for the solutions of the hartree-fock equations for atomic and molecular systems. v.the nonanalytic behavior of the broken-symmetry solutions at the branching point.Phys. Rev. A3, 525–527 (1971)

    work page 1971

  18. [18]

    W., Jim´ enez-Hoyos, C

    Cui, Y., Bulik, I. W., Jim´ enez-Hoyos, C. A., Henderson, T. M. & Scuseria, G. E. Proper and improper zero energy modes in hartree-fock theory and their relevance for symmetry breaking and restoration.J. Chem. Phys.139, 154107 (2013)

    work page 2013

  19. [19]

    Capelle, K., Ullrich, C. A. & Vignale, G. Degenerate ground states and nonunique potentials: Breakdown and restoration of density functionals.Phys. Rev. A76, 012508 (2007)

    work page 2007

  20. [20]

    & Korsell, K

    Alml¨ of, J., Fægri Jr, K. & Korsell, K. Principles for a direct scf approach to licao–moab-initio calculations.J. Comput. Chem.3, 385–399 (1982)

    work page 1982

  21. [21]

    Assessment of density functional methods for obtaining geometries at conical intersections in organic molecules.J

    Filatov, M. Assessment of density functional methods for obtaining geometries at conical intersections in organic molecules.J. Chem. Theory Comput.9, 4526– 4541 (2013)

    work page 2013

  22. [22]

    & Casida, M

    Tapavicza, E., Tavernelli, I., Rothlisberger, U., Filippi, C. & Casida, M. E. Mixed time-dependent density-functional theory/classical trajectory surface hopping study of oxirane photochemistry.J. Chem. Phys.129(2008)

    work page 2008

  23. [23]

    L., Marenich, A

    Li, S. L., Marenich, A. V., Xu, X. & Truhlar, D. G. Configuration interaction- corrected tamm–dancoff approximation: A time-dependent density functional 12 method with the correct dimensionality of conical intersections.J. Phys. Chem. Lett.5, 322–328 (2014)

    work page 2014

  24. [24]

    Shu, Y., Parker, K. A. & Truhlar, D. G. Dual-functional tamm–dancoff approx- imation: a convenient density functional method that correctly describes s1/s0 conical intersections.J. Phys. Chem. Lett.8, 2107–2112 (2017)

    work page 2017

  25. [25]

    & Subotnik, J

    Teh, H.-H. & Subotnik, J. E. The simplest possible approach for simulating s 0–s 1 conical intersections with dft/tddft: Adding one doubly excited configuration. J. Phys. Chem. Lett.10, 3426–3432 (2019)

    work page 2019

  26. [26]

    & Krylov, A

    Shao, Y., Head-Gordon, M. & Krylov, A. I. The spin–flip approach within time- dependent density functional theory: Theory and applications to diradicals.J. Chem. Phys.118, 4807–4818 (2003)

    work page 2003

  27. [27]

    Spin-restricted ensemble-referenced kohn–sham method: basic princi- ples and application to strongly correlated ground and excited states of molecules

    Filatov, M. Spin-restricted ensemble-referenced kohn–sham method: basic princi- ples and application to strongly correlated ground and excited states of molecules. Wiley Interdiscip. Rev. Comput. Mol. Sci.5, 146–167 (2015)

    work page 2015

  28. [28]

    A., Shushkov, P

    Evangelista, F. A., Shushkov, P. & Tully, J. C. Orthogonality constrained density functional theory for electronic excited states.J. Phys. Chem. A117, 7378–7392 (2013)

    work page 2013

  29. [29]

    Schmerwitz, Y. L. A., Urgell Oll´ e, N., Levi, G. & J´ onsson, H. Saddle point search algorithms for variational density functional calculations of excited electronic states with self-interaction correction.Proceedings of the Platform for Advanced Scientific Computing Conference11 (2024)

    work page 2024

  30. [30]

    C., Tao, Z

    Duston, T., Bradbury, N. C., Tao, Z. & Subotnik, J. E. Conical intersections and electronic momentum as viewed from phase space electronic structure theory.J. Phys. Chem. Lett.16, 8994–9003 (2025)

    work page 2025

  31. [33]

    Tuna, D.et al.Assessment of approximate coupled-cluster and algebraic- diagrammatic-construction methods for ground-and excited-state reaction paths and the conical-intersection seam of a retinal-chromophore model.J. Chem. Theory Comput.11, 5758–5781 (2015)

    work page 2015

  32. [34]

    O., Taylor, P

    Roos, B. O., Taylor, P. R. & Sigbahn, P. E. A complete active space scf method (casscf) using a density matrix formulated super-ci approach.Chem. Phys.48, 157–173 (1980). 13

    work page 1980

  33. [35]

    & Werner, H.-J

    Shiozaki, T., Gy˝ orffy, W., Celani, P. & Werner, H.-J. Communication: Extended multi-state complete active space second-order perturbation theory: Energy and nuclear gradients.J. Chem. Phys.135(2011)

    work page 2011

  34. [36]

    Becke, A. D. A new mixing of hartree-fock and local density-functional theories. J. Chem. Phys.98, 1372–1377 (1993)

    work page 1993

  35. [39]

    Marrazzini, G.et al.Multilevel density functional theory.J. Chem. Theory Comput.17, 791–803 (2021)

    work page 2021

  36. [40]

    & Koch, H

    Giovannini, T., Scavino, M. & Koch, H. Time-dependent multilevel density functional theory.J. Chem. Theory Comput.20, 3601–3612 (2024)

    work page 2024

  37. [42]

    & Barone, V

    Adamo, C. & Barone, V. Toward reliable density functional methods without adjustable parameters: The pbe0 model.J. Chem. Phys.110, 6158–6170 (1999)

    work page 1999

  38. [43]

    Barbatti, M., Aquino, A. J. & Lischka, H. Ultrafast two-step process in the non-adiabatic relaxation of the ch2 molecule.Mol. Phys.104, 1053–1060 (2006)

    work page 2006

  39. [44]

    Van Lenthe, J., Zwaans, R., Van Dam, H. J. & Guest, M. Starting scf calculations by superposition of atomic densities.J. Comput. Chem.27, 926–932 (2006)

    work page 2006

  40. [45]

    Coriani, S.et al.An atomic-orbital-based lagrangian approach for calculating geometric gradients of linear response properties.J. Chem. Theory Comput.6, 1028–1047 (2010)

    work page 2010

  41. [46]

    Restoring the Conical Intersection Topology using Convex Density Functional Theory

    Larsen, H., Joørgensen, P., Olsen, J. & Helgaker, T. Hartree–fock and kohn– sham atomic-orbital based time-dependent response theory.J. Chem. Phys.113, 8908–8917 (2000). Acknowledgments We thank Jack T. Taylor and Basile F. Curchod for helpful correspondence that ensured the correct reproduction of the formaldimine results. This work was supported by the ...

    work page 2000

  42. [47]

    & Koch, H

    Rossi, F. & Koch, H. Convex hartree–fock theory for modeling ground state conical intersections.Commun. Chem.(2026)

    work page 2026

  43. [48]

    A., Oliveira, M

    Marques, M. A., Oliveira, M. J. & Burnus, T. Libxc: A library of exchange and correlation functionals for density functional theory.Comput. Phys. Commun. 183, 2272–2281 (2012)

    work page 2012

  44. [49]

    Lehtola, S., Steigemann, C., Oliveira, M. J. & Marques, M. A. Recent devel- opments in libxc—a comprehensive library of functionals for density functional theory.SoftwareX7, 1–5 (2018)

    work page 2018

  45. [50]

    D.et al.eT 1.0: An open source electronic structure program with emphasis on coupled cluster and multilevel methods.J

    Folkestad, S. D.et al.eT 1.0: An open source electronic structure program with emphasis on coupled cluster and multilevel methods.J. Chem. Phys.152, 184103 (2020)

    work page 2020

  46. [51]

    D.et al.et 2.0: An efficient open-source molecular electronic structure program.J

    Folkestad, S. D.et al.et 2.0: An efficient open-source molecular electronic structure program.J. Chem. Phys.164(2026)

    work page 2026

  47. [52]

    T., Tozer, D

    Taylor, J. T., Tozer, D. J. & Curchod, B. F. On the description of conical intersec- tions between excited electronic states with lr-tddft and adc (2).J. Chem. Phys. 159(2023)

    work page 2023

  48. [53]

    Bagel: brilliantly advanced general electronic-structure library.Wiley Interdiscip

    Shiozaki, T. Bagel: brilliantly advanced general electronic-structure library.Wiley Interdiscip. Rev. Comput. Mol. Sci.8, e1331 (2018)

    work page 2018

  49. [54]

    & Zhu, C

    Yu, L., Xu, C. & Zhu, C. Probing theπ→π* photoisomerization mechanism of cis-azobenzene by multi-state ab initio on-the-fly trajectory dynamics simulation. Phys. Chem. Chem. Phys.17, 17646–17660 (2015)

    work page 2015

  50. [55]

    Tuna, D.et al.Assessment of approximate coupled-cluster and algebraic- diagrammatic-construction methods for ground-and excited-state reaction paths and the conical-intersection seam of a retinal-chromophore model.J. Chem. Theory Comput.11, 5758–5781 (2015). S15

    work page 2015