Pith

open record

sign in

arxiv: 2605.22584 · v2 · pith:WINWRSWO · submitted 2026-05-21 · math.NA · cs.NA· physics.chem-ph

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

Jonas Beck , Benjamin Stamm This is my paper

Reviewed by Pith T0 review T1 audit T2 compute T3 formal T4 kernel 2026-07-02 23:34 UTCgrok-4.3pith:WINWRSWOrecord.jsonopen to challenge →

classification math.NA cs.NAphysics.chem-ph
keywords coupled clusteramplitudesregularityreal analyticitynuclear coordinatesinterpolationHartree-Fockcanonical orbitals
0
0 comments X

The pith

Under non-degeneracy assumptions, coupled cluster amplitudes are real analytic functions of nuclear coordinates.

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 vary in a real analytic way with nuclear coordinates when certain non-degeneracy conditions hold at the Hartree-Fock and coupled cluster levels. This regularity would allow accurate interpolation of amplitudes from a small number of computed geometries instead of recalculating at every point. The authors also examine how canonical orbital choices introduce artifacts that reduce this smoothness in practice and propose ways to lessen those effects. Numerical tests confirm that interpolants based on the amplitudes perform close to exact calculations when the assumptions are met.

Core claim

The paper shows that, under non-degeneracy assumptions on the Hartree-Fock level of theory and the coupled cluster level of theory, the coupled cluster amplitudes behave as real analytic functions of the nuclear coordinates. It further identifies artifacts arising from the use of canonical orbitals in practical calculations that disrupt this analyticity and suggests mitigation strategies. The findings are validated by comparing interpolated amplitudes to exact ones in numerical experiments.

What carries the argument

Real analyticity of the amplitudes with respect to nuclear displacements, proven under non-degeneracy assumptions on Hartree-Fock and coupled cluster theory, which enables interpolation from limited reference geometries.

If this is right

  • Amplitudes can be interpolated from a limited set of reference geometries while retaining high accuracy.
  • Canonical orbital artifacts can be mitigated to restore closer to analytic behavior in practical calculations.
  • Computational cost for studies over many nuclear coordinates drops because full coupled cluster solves are needed at fewer points.

Where Pith is reading between the lines

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

  • Similar regularity arguments might apply to other single-reference methods if their defining equations satisfy comparable non-degeneracy conditions.
  • The approach could reduce the number of geometry samples required in molecular dynamics or potential energy surface scans.
  • Testing whether the non-degeneracy assumptions hold for typical molecular systems would determine how often the analyticity result applies in practice.

Load-bearing premise

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

What would settle it

A numerical example in which the non-degeneracy conditions are violated at either level and the amplitudes exhibit a non-analytic point or discontinuity as nuclear coordinates vary.

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 / 3 minor

Summary. The manuscript claims that, under explicit non-degeneracy assumptions on both the Hartree-Fock and coupled-cluster levels of theory, single-reference coupled-cluster amplitudes are real analytic functions of nuclear displacements. It separates this theoretical result from an analysis of canonical-orbital artifacts that reduce observed regularity in practice, proposes mitigation strategies, and reports numerical experiments that compare amplitude interpolants against exact amplitudes computed at additional geometries.

Significance. If the conditional analyticity result holds, the work supplies a parameter-free theoretical basis for amplitude interpolation over nuclear coordinate spaces, which could reduce the number of full CC calculations required for potential-energy surfaces. The manuscript credits the conditional nature of the claim and supplies both the analytic argument and numerical validation; these elements strengthen the contribution within numerical analysis of quantum-chemical methods.

major comments (1)
  1. [§5] §5 (numerical experiments): the reported interpolation comparisons do not verify that the non-degeneracy assumptions of §3 hold for the chosen molecular systems and displacements; without such verification the numerical results cannot directly confirm that the observed performance stems from the proven analyticity rather than from other factors.
minor comments (3)
  1. [Abstract, §2] Abstract and §2: the phrase 'poly(N)' scaling is used without specifying the degree; a brief parenthetical (e.g., O(N^7) for CCSD(T)) would improve precision.
  2. Figure captions: several figures lack explicit labels for the x-axis (nuclear displacement coordinate) and y-axis (amplitude component or error); this reduces immediate readability.
  3. [§4] §4 (mitigation strategies): the description of orbital-rotation corrections would benefit from a short pseudocode or explicit formula for the transformation applied to the amplitudes.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and the constructive comment on the connection between the theoretical assumptions and the numerical experiments. We agree that the point merits attention and will revise the manuscript to address it.

read point-by-point responses

  1. Referee: [§5] §5 (numerical experiments): the reported interpolation comparisons do not verify that the non-degeneracy assumptions of §3 hold for the chosen molecular systems and displacements; without such verification the numerical results cannot directly confirm that the observed performance stems from the proven analyticity rather than from other factors.

    Authors: We agree with the referee that an explicit verification of the non-degeneracy assumptions would strengthen the link between the analyticity result of §3 and the observed interpolation performance in §5. In the revised manuscript we will add a short verification subsection (or paragraph) to §5. For each molecular system and displacement range considered, we will report the minimal distance to degeneracy at the Hartree–Fock level (via the smallest singular value of the orbital Hessian or the HOMO–LUMO gap, as appropriate) and at the coupled-cluster level (via the smallest singular value of the CC Jacobian). These quantities are already defined in §3; we will simply evaluate and tabulate them for the geometries used in the numerical tests. If any near-degeneracies are detected we will note them and discuss their possible influence on the interpolation results. This addition will make the numerical evidence more directly supportive of the conditional analyticity claim without altering the existing figures or tables. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper derives real-analyticity of CC amplitudes w.r.t. nuclear displacements directly from stated non-degeneracy assumptions on the HF and CC levels; the abstract and described structure present this as a conditional mathematical result, not a fit, self-definition, or self-citation chain. No equations reduce the claimed regularity to input data or prior author work by construction. Numerical experiments are presented only as validation after the theoretical statement, and canonical-orbital artifacts are analyzed separately. The derivation is therefore self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The theoretical result rests on domain assumptions from quantum chemistry; no free parameters or invented entities appear in the abstract.

axioms (1)
  • domain assumption non-degeneracy assumptions on the Hartree-Fock level of theory and the coupled cluster level of theory
    These assumptions are invoked to guarantee that the amplitudes behave real analytic with respect to nuclear coordinate displacements.

pith-pipeline@v0.9.1-grok · 5735 in / 1114 out tokens · 31675 ms · 2026-07-02T23:34:18.924737+00:00 · methodology

Review history (2 revisions) →

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

    work page 1993

  2. [2]

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

    work page 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

    work page 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

    work page 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

    work page 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

    work page 2000

  7. [7]

    Coester and H

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

    work page 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

    work page 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

    work page 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

    work page 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

    work page 1910

  12. [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. [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

    work page 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

    work page 1928

  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

    work page 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

    work page 2023

  17. [17]

    Helgaker, P

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

    work page 2000

  18. [18]

    Hohenberg and W

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

    work page 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

    work page 1957

  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

    work page 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

    work page 1951

  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

    work page 2019

  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

    work page 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

    work page 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

    work page 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

    work page 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

    work page 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

    work page 1969

  29. [29]

    Rohwedder

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

    work page 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

    work page 2013

  31. [31]

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

    work page 1951

  32. [32]

    Schneider

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

    work page 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

    work page 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

    work page 1962

  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

    work page 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

    work page 1996

  37. [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. [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

    work page 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...

    work page 1966