pith. sign in

arxiv: 2605.22584 · v1 · pith:WINWRSWOnew · submitted 2026-05-21 · 🧮 math.NA · cs.NA· physics.chem-ph

On the Regularity and Interpolation of Coupled Cluster Amplitudes in Canonical Orbital Basis

Pith reviewed 2026-05-22 03:13 UTC · model grok-4.3

classification 🧮 math.NA cs.NAphysics.chem-ph
keywords coupled clusteramplitudesregularityinterpolationnuclear coordinatesanalyticityHartree-Fockcanonical orbitals
0
0 comments X

The pith

Coupled cluster amplitudes are real analytic functions of nuclear coordinates under non-degeneracy assumptions.

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

The paper establishes that single-reference coupled cluster amplitudes behave as real analytic functions of nuclear displacements when certain non-degeneracy conditions hold at the Hartree-Fock and coupled cluster levels. This regularity would allow the use of interpolation techniques to compute amplitudes at many geometries using only a few reference calculations, reducing the overall computational burden for studying chemical phenomena over ranges of nuclear coordinates. The authors also examine practical issues arising from canonical orbitals in computations and propose ways to address them, with validation through numerical interpolation tests.

Core claim

Under certain non-degeneracy assumptions on the Hartree-Fock level of theory and the coupled cluster level of theory, the coupled cluster amplitudes are real analytic functions of the nuclear coordinates.

What carries the argument

Real analyticity of the coupled cluster amplitudes with respect to nuclear displacements under non-degeneracy assumptions, which supports interpolation from limited reference geometries.

Load-bearing premise

The non-degeneracy assumptions on the Hartree-Fock and coupled cluster levels of theory hold true.

What would settle it

Compute amplitudes at multiple nuclear geometries, then test whether a high-order polynomial or spline interpolant reproduces exact amplitudes with error decreasing faster than any algebraic rate; failure to do so would contradict the analyticity claim.

read the original abstract

Arguably the most widely used approaches for obtaining highly accurate molecular ground-state energies are coupled cluster methods. Despite introducing two layers of approximation, a linear and a nonlinear one, coupled cluster methods remain computationally intensive, with the complexity scaling as $O(poly(N))$, where $N$ is the number of electrons. Moreover, this method must be applied over a large set of different nuclear coordinates in order to study certain chemical phenomena. Therefore, in this work, we investigate the regularity of single-reference coupled cluster amplitudes with respect to nuclear coordinate displacements, with the aim of enabling interpolation or extrapolation approaches that rely on only a limited number of reference geometries. We show that, in theory, under certain non-degeneracy assumptions on the Hartree-Fock level of theory, and the coupled cluster level of theory the amplitudes behave real analytic. Furthermore, we analyze the artifacts that arise in practical calculations that use canonical orbitals, which hinder this high degree of regularity, and suggest strategies to mitigate these issues. Finally, we validate our findings through numerical experiments by interpolating the amplitudes and comparing the performance of the interpolants with that of the exact amplitudes.

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

1 major / 0 minor

Summary. The manuscript claims that under non-degeneracy assumptions on the Hartree-Fock orbital equations and the coupled-cluster amplitude equations, the CC amplitudes are real-analytic functions of nuclear coordinates in the canonical orbital basis. It analyzes practical artifacts that disrupt this regularity in computations, proposes mitigation strategies, and validates the approach via numerical experiments that compare amplitude interpolants against exact values at additional geometries.

Significance. If the global regularity result holds, the work would support reduced-cost interpolation or extrapolation of CC amplitudes across nuclear geometries, lowering the expense of constructing potential energy surfaces. The numerical validation provides concrete evidence that interpolation can be competitive when artifacts are controlled, and the explicit invocation of the analytic implicit-function theorem supplies a clear theoretical route.

major comments (1)
  1. The central claim of real-analytic dependence on nuclear coordinates (Abstract) rests on the analytic implicit-function theorem applied to the HF and CC equations. This requires the relevant Jacobians to remain invertible throughout the connected domain of geometries. The manuscript states the non-degeneracy assumptions but supplies neither a proof that eigenvalue gaps and Jacobian singular values stay bounded away from zero under continuous nuclear motion nor numerical checks of minimal singular values along representative paths. This gap is load-bearing for the global analyticity needed to justify reliable interpolation over extended domains.

Simulated Author's Rebuttal

1 responses · 1 unresolved

We thank the referee for the careful reading and for identifying a key point regarding the scope of our global regularity claim. We respond to the major comment below and outline the revisions we will undertake.

read point-by-point responses
  1. Referee: The central claim of real-analytic dependence on nuclear coordinates (Abstract) rests on the analytic implicit-function theorem applied to the HF and CC equations. This requires the relevant Jacobians to remain invertible throughout the connected domain of geometries. The manuscript states the non-degeneracy assumptions but supplies neither a proof that eigenvalue gaps and Jacobian singular values stay bounded away from zero under continuous nuclear motion nor numerical checks of minimal singular values along representative paths. This gap is load-bearing for the global analyticity needed to justify reliable interpolation over extended domains.

    Authors: We agree that the analytic implicit-function theorem yields local real-analytic dependence wherever the Jacobians are nonsingular. Global analyticity over a connected domain therefore presupposes that the stated non-degeneracy conditions (invertibility of the HF and CC Jacobians) hold uniformly along all paths within that domain. The manuscript explicitly invokes these assumptions but does not contain a general proof that the eigenvalue gaps remain bounded away from zero for arbitrary continuous nuclear displacements; such a proof would be system-dependent and lies outside the present scope. To strengthen the practical justification for interpolation, we will add, in the revised numerical section, explicit checks of the minimal singular values of the HF and CC Jacobians evaluated at the geometries sampled along the representative paths. These diagnostics will confirm that the assumptions remain satisfied in the tested cases. We will also insert a clarifying paragraph in the theoretical development that distinguishes the local character of the implicit-function theorem from the global assumption required for the full domain. This constitutes a partial revision focused on numerical support and textual clarification. revision: partial

standing simulated objections not resolved
  • A general proof that eigenvalue gaps and Jacobian singular values remain bounded away from zero under arbitrary continuous nuclear motions for all molecular systems.

Circularity Check

0 steps flagged

No significant circularity: derivation applies standard analytic implicit-function theorem to HF and CC equations

full rationale

The central claim establishes real-analytic dependence of coupled-cluster amplitudes on nuclear coordinates by invoking the analytic implicit-function theorem on the Hartree-Fock orbital equations and the CC amplitude equations, conditioned on explicit non-degeneracy assumptions (invertibility of the Fock-matrix eigenvalue gaps and the CC Jacobian). This is a direct application of a classical external theorem rather than a self-referential construction; the assumptions are stated as hypotheses and not derived from the target result. Numerical interpolation experiments are presented as independent empirical checks against exact amplitudes, not as fitted predictions that tautologically reproduce the inputs. No load-bearing step reduces by construction to a parameter fit, a self-citation chain, or a renaming of known patterns. The derivation therefore remains self-contained against external mathematical benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim depends on non-degeneracy assumptions at the Hartree-Fock and coupled cluster levels; these are domain assumptions rather than free parameters or new entities.

axioms (2)
  • domain assumption non-degeneracy assumptions on the Hartree-Fock level of theory
    Invoked to guarantee real-analytic behavior of the amplitudes as stated in the abstract.
  • domain assumption non-degeneracy assumptions on the coupled cluster level of theory
    Invoked to guarantee real-analytic behavior of the amplitudes as stated in the abstract.

pith-pipeline@v0.9.0 · 5735 in / 1211 out tokens · 65238 ms · 2026-05-22T03:13:55.539856+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

39 extracted references · 39 canonical work pages

  1. [1]

    Agarwal and P

    R. Agarwal and P. Wong.Error Inequalities in Polynomial Interpolation and Their Applications. Mathematics and its applications. Kluwer Academic Publishers, 1993

  2. [2]

    J. Beck. Interpolation methods for post-hartree-fock calculations in electronic structure theory. Master’s thesis, 2025

  3. [3]

    Beylkin and S

    G. Beylkin and S. Sharma. A fast algorithm for computing the boys function.The Journal of Chemical Physics, 155(17), November 2021

  4. [4]

    S. F. Boys. Electronic wave functions - i. a general method of calculation for the stationary states of any molecular system.Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences, 200(1063):542–554, 02 1950

  5. [5]

    Cancès, G

    E. Cancès, G. Kemlin, and A. Levitt. Convergence analysis of direct minimization and self-consistent iterations.SIAM Journal on Matrix Analysis and Applications, 42(1):243–274, 2021

  6. [6]

    On the convergence of scf algorithms for the hartree-fock equations

    Cancès, Eric and Le Bris, Claude. On the convergence of scf algorithms for the hartree-fock equations. ESAIM: M2AN, 34(4):749–774, 2000. [7]F. Coester. Bound states of a many-particle system.Nuclear Physics, 7:421–424, 1958

  7. [7]

    Coester and H

    F. Coester and H. Kümmel. Short-range correlations in nuclear wave functions.Nuclear Physics, 17:477–485, 1960

  8. [8]

    T. D. Crawford and H. F. Schaefer III.An Introduction to Coupled Cluster Theory for Computational Chemists, pages 33–136. John Wiley & Sons, Ltd, 2000

  9. [9]

    M. A. Csirik and A. Laestadius. Coupled-cluster theory revisited: Part i: Discretization.ESAIM: Mathematical Modelling and Numerical Analysis, 57(2):645–670, March 2023

  10. [10]

    M. A. Csirik and A. Laestadius. Coupled-cluster theory revisited: Part ii: Analysis of the single-reference coupled-cluster equations.ESAIM: Mathematical Modelling and Numerical Analysis, 57(2):545–583, March 2023

  11. [11]

    Dieci, A

    L. Dieci, A. Papini, and A. Pugliese. Decompositions and coalescing eigenvalues of symmetric definite pencils depending on parameters.Numerical Algorithms, 91(4):1879–1910, June 2022. [13]F. F aulstich, B. Sturmfels, and S. Sverrisdóttir. Algebraic varieties in quantum chemistry, 2024

  12. [12]

    V. Fock. Näherungsmethode zur Lösung des quantenmechanischen Mehrkörperproblems.Zeitschrift fur Physik, 61(1-2):126–148, January 1930

  13. [13]

    Grazioli, Y

    L. Grazioli, Y. Hu, and E. Cancès. Critical point search and linear response theory for computing electronic excitation energies of molecular systems. i. general framework, application to hartree–fock and dft.The Journal of Chemical Physics, 164(6):064101, 02 2026

  14. [14]

    D. R. Hartree. The wave mechanics of an atom with a non-coulomb central field. part i. theory and methods. Mathematical Proceedings of the Cambridge Philosophical Society, 24(1):89–110, 1928. [17]D. Hartree.The Calculation of Atomic Structures. Structure of matter series. J. Wiley, 1957

  15. [15]

    Hassan and Y

    M. Hassan and Y. Maday. Analysis of the single reference coupled cluster method for electronic structure calculations: The discrete coupled cluster equations, 2025

  16. [16]

    Hassan, Y

    M. Hassan, Y. Maday, and Y. W ang. Analysis of the single reference coupled cluster method for electronic structure calculations: The full coupled cluster equations, 2023

  17. [17]

    Helgaker, P

    T. Helgaker, P. Jørgensen, and J. Olsen.Molecular Electronic Structure Theory. John Wiley & Sons, LTD, Chichester, 2000

  18. [18]

    Hohenberg and W

    P. Hohenberg and W. Kohn. Inhomogeneous Electron Gas.Physical Review, 136(3B):864–871, November 1964

  19. [19]

    J. Hubbard. The description of collective motions in terms of many-body perturbation theory.Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences, 240(1223):539–560, 1957. [23]N. Hugenholtz. Perturbation theory of large quantum systems.Physica, 23(1):481–532, 1957. ON THE REGULARITY OF COUPLED CLUSTER AMPLITUDES 29

  20. [20]

    Hunziker and I

    W. Hunziker and I. M. Sigal. The quantum n-body problem.Journal of Mathematical Physics, 41(6):3448–3510, 06 2000

  21. [21]

    T. Kato. Fundamental properties of hamiltonian operators of schrödinger type.Transactions of the American Mathematical Society, 70:195–211, 1951. [26]T. Kato.Perturbation Theory for Linear Operators. 1966. [27]T. Kolda and B. Bader. Tensor decompositions and applications.SIAM Review, 51:455–500, 08 2009

  22. [22]

    Laestadius and F

    A. Laestadius and F. M. F aulstich. The coupled-cluster formalism – a mathematical perspective.Molecular Physics, 117(17):2362–2373, January 2019. [29]J. Mason and D. Handscomb.Chebyshev Polynomials. CRC Press, 2002

  23. [23]

    L. E. McMurchie and E. R. Davidson. One- and two-electron integrals over cartesian gaussian functions. Journal of Computational Physics, 26(2):218–231, 1978

  24. [24]

    Mills and A

    I. Mills and A. G. Robiette. On the relationship of normal modes to local modes in molecular vibrations. Molecular Physics, 56(4):743–765, 1985

  25. [25]

    Obara and A

    S. Obara and A. Saika. General recurrence formulas for molecular integrals over cartesian gaussian functions. The Journal of Chemical Physics, 89(3):1540–1559, 08 1988

  26. [26]

    J. Paldus. Correlation problems in atomic and molecular systems. v. spin-adapted coupled cluster many-electron theory.The Journal of Chemical Physics, 67(1):303–318, 07 1977

  27. [27]

    Raghavachari, G

    K. Raghavachari, G. W. Trucks, J. A. Pople, and M. Head-Gordon. A fifth-order perturbation comparison of electron correlation theories.Chemical Physics Letters, 157(6):479–483, 1989

  28. [28]

    Rellich.Perturbation Theory of Eigenvalue Problems

    F. Rellich.Perturbation Theory of Eigenvalue Problems. New York University. Institute of Mathematical Sciences. Gordon and Breach, 1969

  29. [29]

    Rohwedder

    T. Rohwedder. The continuous coupled cluster formulation for the electronic schrödinger equation.preprint, submitted, 01 2011

  30. [30]

    Error estimates for the coupled cluster method.ESAIM: M2AN, 47(6):1553–1582, 2013

    Rohwedder, Thorsten and Schneider, Reinhold. Error estimates for the coupled cluster method.ESAIM: M2AN, 47(6):1553–1582, 2013

  31. [31]

    C. C. Roothaan. New developments in molecular orbital theory.Reviews of Modern Physics, 23(2):69–89, April 1951

  32. [32]

    Schneider

    R. Schneider. Analysis of the projected coupled cluster method in electronic structure calculation.Numerische Mathematik, 113:433–471, 08 2009

  33. [33]

    S. E. Schrader and S. Kvaal. Accelerated coupled cluster calculations with procrustes orbital interpolation. The Journal of Chemical Physics, 158(11):114116, 03 2023

  34. [34]

    Sinanoğlu

    O. Sinanoğlu. Many-electron theory of atoms and molecules. i. shells, electron pairs vs many-electron correlations. The Journal of Chemical Physics, 36(3):706–717, 02 1962. [42]J. C. Slater. Note on hartree’s method.Phys. Rev., 35:210–211, Jan 1930

  35. [35]

    Q. Sun, T. C. Berkelbach, N. S. Blunt, G. H. Booth, S. Guo, Z. Li, J. Liu, J. D. McClain, E. R. Sayfutyarova, S. Sharma, S. Wouters, and G. K.-L. Chan. Pyscf: the python-based simulations of chemistry framework.WIREs Computational Molecular Science, 8(1):e1340, 2018

  36. [36]

    Szabo and N

    A. Szabo and N. Ostlund.Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory. Dover Books on Chemistry. Dover Publications, 1996

  37. [37]

    E. Tadmor. The exponential accuracy of fourier and chebyshev differencing methods.SIAM Journal on Numerical Analysis, 23(1):1–10, 1986

  38. [38]

    W ang and C

    L.-P. W ang and C. Song. Geometry optimization made simple with translation and rotation coordinates.The Journal of Chemical Physics, 144(21):214108, 06 2016

  39. [39]

    J. Čížek. On the correlation problem in atomic and molecular systems. calculation of wavefunction components in ursell-type expansion using quantum-field theoretical methods.The Journal of Chemical Physics, 45(11):4256– 4266, 12 1966. ⋆ Jonas Beck, Institute of Applied Analysis and Numerical Simulation, University of Stuttgart, Pfaffenwaldring 57, 70569 S...