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
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.
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.
Referee Report
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)
- 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
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
-
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
- 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
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
axioms (2)
- domain assumption non-degeneracy assumptions on the Hartree-Fock level of theory
- domain assumption non-degeneracy assumptions on the coupled cluster level of theory
Reference graph
Works this paper leans on
-
[1]
R. Agarwal and P. Wong.Error Inequalities in Polynomial Interpolation and Their Applications. Mathematics and its applications. Kluwer Academic Publishers, 1993
work page 1993
-
[2]
J. Beck. Interpolation methods for post-hartree-fock calculations in electronic structure theory. Master’s thesis, 2025
work page 2025
-
[3]
G. Beylkin and S. Sharma. A fast algorithm for computing the boys function.The Journal of Chemical Physics, 155(17), November 2021
work page 2021
-
[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
work page 1950
- [5]
-
[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
work page 2000
-
[7]
F. Coester and H. Kümmel. Short-range correlations in nuclear wave functions.Nuclear Physics, 17:477–485, 1960
work page 1960
-
[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
work page 2000
-
[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
work page 2023
-
[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
work page 2023
-
[11]
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
work page 1910
-
[12]
V. Fock. Näherungsmethode zur Lösung des quantenmechanischen Mehrkörperproblems.Zeitschrift fur Physik, 61(1-2):126–148, January 1930
work page 1930
-
[13]
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
work page 2026
-
[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
work page 1928
-
[15]
M. Hassan and Y. Maday. Analysis of the single reference coupled cluster method for electronic structure calculations: The discrete coupled cluster equations, 2025
work page 2025
- [16]
-
[17]
T. Helgaker, P. Jørgensen, and J. Olsen.Molecular Electronic Structure Theory. John Wiley & Sons, LTD, Chichester, 2000
work page 2000
-
[18]
P. Hohenberg and W. Kohn. Inhomogeneous Electron Gas.Physical Review, 136(3B):864–871, November 1964
work page 1964
-
[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
work page 1957
-
[20]
W. Hunziker and I. M. Sigal. The quantum n-body problem.Journal of Mathematical Physics, 41(6):3448–3510, 06 2000
work page 2000
-
[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
work page 1951
-
[22]
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
work page 2019
-
[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
work page 1978
-
[24]
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
work page 1985
-
[25]
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
work page 1988
-
[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
work page 1977
-
[27]
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
work page 1989
-
[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
work page 1969
- [29]
-
[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
work page 2013
-
[31]
C. C. Roothaan. New developments in molecular orbital theory.Reviews of Modern Physics, 23(2):69–89, April 1951
work page 1951
- [32]
-
[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
work page 2023
- [34]
-
[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
work page 2018
-
[36]
A. Szabo and N. Ostlund.Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory. Dover Books on Chemistry. Dover Publications, 1996
work page 1996
-
[37]
E. Tadmor. The exponential accuracy of fourier and chebyshev differencing methods.SIAM Journal on Numerical Analysis, 23(1):1–10, 1986
work page 1986
-
[38]
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
work page 2016
-
[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...
work page 1966
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.