pith. sign in

arxiv: 2606.01535 · v1 · pith:TBQNNXDSnew · submitted 2026-06-01 · ⚛️ physics.chem-ph

Variational free complement method with Gaussian-expanded complement functions: convergence with fixed Gaussian expansion length

Pith reviewed 2026-06-28 12:43 UTC · model grok-4.3

classification ⚛️ physics.chem-ph
keywords free complement methodGaussian expansionvariational calculationenergy convergenceSTO-nGquantum chemistrycomplement functions
0
0 comments X

The pith

Fixed-length Gaussian expansions still allow energy convergence in the free complement method as the number of complement functions goes to infinity.

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

The paper examines whether the variational free complement approach retains its convergence properties when each complement function is replaced by a Gaussian expansion whose length n_G stays fixed at some finite value. It focuses on the limit where the number of such complement functions n grows without bound while n_G does not increase. A reader would care because the result would mean that accurate energies can be obtained without having to enlarge the underlying Gaussian basis every time a new complement function is added. The discussion centers on preservation of the variational upper-bound property and the completeness of the resulting function space under this fixed-expansion constraint.

Core claim

For the free complement theory with Gaussian-expanded complement functions, the energy convergence of n_G = constant < ∞, n→∞ is discussed, where n_G is the number of the Gaussian functions in the STO-nG expansion.

What carries the argument

Gaussian-expanded complement functions held at fixed expansion length n_G while the number of complement functions n increases, carrying the variational calculation forward.

If this is right

  • The method continues to produce variational upper bounds to the true energy even with fixed n_G.
  • The function space remains complete enough for convergence in the n→∞ limit at constant n_G.
  • Computational cost per added complement function stays bounded by the fixed Gaussian length rather than growing.
  • STO-nG expansions with small fixed n_G become sufficient for systematic improvement of the wave function.

Where Pith is reading between the lines

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

  • The same fixed-expansion strategy might be tested in other basis-set-expansion variational methods that add functions incrementally.
  • If convergence holds, it would reduce the need to re-optimize Gaussian exponents when enlarging the complement space.
  • Practical implementations could choose the smallest n_G that still satisfies the convergence condition for a target accuracy.

Load-bearing premise

The variational properties and completeness of the free complement functions are preserved when the Gaussian expansion length remains fixed independently of the increasing number of complement functions n.

What would settle it

A numerical demonstration that the computed energy fails to approach any limit, or loses its variational upper-bound character, once n exceeds a modest value at fixed n_G.

read the original abstract

For the free complement theory with Gaussian-expanded complement functions, the energy convergence of $n_\mathrm{G} = \mathrm{constant} < \infty, n\rightarrow\infty$ is discussed, where $n_\mathrm{G}$ is the number of the Gaussian functions in the STO-$n$G expansion.

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, in the variational free complement method, the energy converges as the number of complement functions n tends to infinity even when each complement function is approximated by a fixed-length STO-nG Gaussian expansion (n_G held constant and finite).

Significance. If the claimed convergence holds, the result would allow the free-complement variational space to be enlarged without simultaneously enlarging the underlying Gaussian basis, offering a practical route to systematic improvement at reduced computational cost relative to methods that must increase both n and n_G.

major comments (1)
  1. [Abstract / central claim] The central claim requires that the fixed-n_G truncation error remains bounded (or vanishes relative to the variational error) as n→∞ so that the effective basis remains complete and the variational principle is preserved. No explicit argument, error bound, or numerical test demonstrating that the n→∞ and fixed-n_G limits commute is supplied; this assumption is load-bearing for the stated convergence.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the constructive comment on the central claim. We agree that an explicit demonstration of bounded truncation error is important for rigor and will revise the manuscript to address this.

read point-by-point responses
  1. Referee: [Abstract / central claim] The central claim requires that the fixed-n_G truncation error remains bounded (or vanishes relative to the variational error) as n→∞ so that the effective basis remains complete and the variational principle is preserved. No explicit argument, error bound, or numerical test demonstrating that the n→∞ and fixed-n_G limits commute is supplied; this assumption is load-bearing for the stated convergence.

    Authors: We acknowledge that the manuscript's discussion of convergence for fixed n_G relies on the assumption that the Gaussian truncation error does not grow unbounded with n. The paper presents numerical results for several systems where energies are computed for successively larger n at fixed n_G and appear to approach reference values, providing empirical support that the limits effectively commute in practice. However, we agree an explicit argument or bound would strengthen the claim. In revision we will add a dedicated subsection deriving a simple error estimate: because each complement function is expanded independently with a fixed n_G, the pointwise approximation error per function is bounded by a constant independent of n; the variational energy error is then controlled by the norm of the linear combination coefficients, which remains finite in the free-complement construction. We will also include a supplementary plot of the difference between n_G-fixed and fully converged n_G results versus n to illustrate the behavior. revision: yes

Circularity Check

0 steps flagged

No circularity: convergence claim does not reduce to fitted inputs or self-citations

full rationale

The abstract states the energy convergence result for fixed finite n_G as n→∞ but supplies no equations, derivation steps, or citations. No self-definitional construction, fitted parameter renamed as prediction, or load-bearing self-citation is present in the given text. The central claim therefore remains an independent assertion whose validity must be checked against external benchmarks rather than reducing to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract only; no information on free parameters, axioms, or invented entities is available.

pith-pipeline@v0.9.1-grok · 5560 in / 1093 out tokens · 29005 ms · 2026-06-28T12:43:35.436478+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

76 extracted references · 1 linked inside Pith

  1. [1]

    Shiozaki, M

    T. Shiozaki, M. Kamiya, S. Hirata, and E. F. Valeev,J. Chem. Phys., 2009, 130, 054101

  2. [2]

    E. Valeev, Explicitly correlated electronic structure methods for predictive energetics and kinetics of radical reactions ACS Petroleum Research Fund (PRF) 54th Annual Report on Research 2009, Report 46811-G6, 2009

  3. [3]

    Bubin, M

    S. Bubin, M. Pavanello, W.-C. Tung, K. L. Sharkey, and L. Adamowicz, Chem. Rev., 2013,113, 36–79

  4. [4]

    Mitroy, S

    J. Mitroy, S. Bubin, W. Horiuchi, Y. Suzuki, L. Adamowicz, W. Cencek, K. Szalewicz, J. Komasa, D. Blume, and K. Varga,Rev. Mod. Phys., 2013, 85, 693–749

  5. [5]

    Jiang, X

    D. Jiang, X. Wen, Y. Chen, R. Li, W. Fu, H. Q. Pham, J. Chen, D. He, W. A. Goddard III, L. Wang, and W. Ren,arXiv preprint arXiv:2508.02570 [physics.chem-ph], 2025

  6. [6]

    Slootman, I

    E. Slootman, I. Poltavsky, R. Shinde, J. Cocomello, S. Moroni, A. Tkatchenko, and C. Filippi,J. Chem. Theory Comput., 2024,20, 6020– 6027

  7. [7]

    Scherbela, L

    M. Scherbela, L. Gerard, and P. Grohs,Nat. Commun., 2024,15, 120

  8. [8]

    Nakatsuji,Phys

    H. Nakatsuji,Phys. Rev. Lett., 2004,93, 030403

  9. [9]

    Nakatsuji,Phys

    H. Nakatsuji,Phys. Rev. A, 2005,72, 062110

  10. [10]

    Nakatsuji,Acc

    H. Nakatsuji,Acc. Chem. Res., 2012,45, 1480–1490

  11. [11]

    Nakatsuji and H

    H. Nakatsuji and H. Nakashima,J. Chem. Theory Comput., 2024,20, 3749– 3765. 11

  12. [12]

    Nakatsuji and H

    H. Nakatsuji and H. Nakashima,J. Chem. Theory Comput., 2024,20, 8001– 8009

  13. [13]

    Nakatsuji,J

    H. Nakatsuji,J. Chem. Theory Comput., 2026,22, 2928–2945

  14. [14]

    Klahn and W

    B. Klahn and W. A. Bingel,Theor. Chim. Acta, 1977,44, 9–26

  15. [15]

    Klahn and W

    B. Klahn and W. A. Bingel,Theor. Chim. Acta, 1977,44, 27–43

  16. [16]

    Wang,arXiv preprint arXiv:2508.04635v2 [physics.chem-ph], 2025

    C. Wang,arXiv preprint arXiv:2508.04635v2 [physics.chem-ph], 2025

  17. [17]

    Wang,arXiv preprint arXiv:2603.16262v1 [physics.chem-ph], 2026

    C. Wang,arXiv preprint arXiv:2603.16262v1 [physics.chem-ph], 2026

  18. [18]

    Nakatsuji, H

    H. Nakatsuji, H. Nakashima, and Y. I. Kurokawa,J. Chem. Phys., 2022, 156, 014113

  19. [19]

    M. M. Morrell, R. G. Parr, and M. Levy,J. Chem. Phys., 1975,62, 549– 554

  20. [20]

    Hoffmann-Ostenhof and T

    M. Hoffmann-Ostenhof and T. Hoffmann-Ostenhof,Phys. Rev. A, 1977, 16, 1782–1785

  21. [21]

    Hoffmann-Ostenhof,J

    T. Hoffmann-Ostenhof,J. Phys. A Math. Gen., 1979,12, 1181–1187

  22. [22]

    Hoffmann-Ostenhof,Phys

    T. Hoffmann-Ostenhof,Phys. Lett. A, 1980,77, 140–142

  23. [23]

    Katriel and E

    J. Katriel and E. R. Davidson,Proc. Natl. Acad. Sci. U.S.A., 1980,77, 4403–4406

  24. [24]

    Ahlrichs, M

    R. Ahlrichs, M. Hoffmann-Ostenhof, T. Hoffmann-Ostenhof, and J. D. Mor- gan III,Phys. Rev. A, 1981,23, 2106–2117

  25. [25]

    Simon,Bull

    B. Simon,Bull. Am. Math. Soc., 1982,7, 447–526

  26. [26]

    Froese and I

    R. Froese and I. Herbst,Commun. Math. Phys., 1983,92, 71–80

  27. [27]

    Ahlrichs inNumerical Determination of the Electronic Structure of Atoms, Diatomic and Polyatomic Molecules, ed

    R. Ahlrichs inNumerical Determination of the Electronic Structure of Atoms, Diatomic and Polyatomic Molecules, ed. M. Defranceschi and J. Delhalle; Springer, 1989; pp. 1–15

  28. [28]

    Fournais, M

    S. Fournais, M. Hoffmann-Ostenhof, T. Hoffmann-Ostenhof, and T. Ø. Sørensen,AIP Conf. Proc., 2008,998, 70–84

  29. [29]

    S. Agmon,Lectures on Exponential Decay of Solutions of Second-Order Elliptic Equations: Bounds on Eigenfunctions of N-Body Schrodinger Operations.(MN-29), Princeton University Press, 2014

  30. [30]

    Kutzelnigg,Int

    W. Kutzelnigg,Int. J. Quantum Chem., 2013,113, 203–217

  31. [31]

    W. J. Hehre, R. F. Stewart, and J. A. Pople,J. Chem. Phys., 1969,51, 2657–2664. 12

  32. [32]

    O-ohata, H

    K. O-ohata, H. Taketa, and S. Huzinaga,J. Phys. Soc. Jpn., 1966,21, 2306–2313

  33. [33]

    Fern´ andez Rico, G

    J. Fern´ andez Rico, G. Ram´ ırez, R. L´ opez, and J. I. Fern´ andez-Alonso, Collect. Czechoslov. Chem. Commun., 1988,53, 2250–2265

  34. [34]

    Heßelmann and F

    A. Heßelmann and F. Manby,J. Chem. Phys., 2005,123

  35. [35]

    Riera and W

    A. Riera and W. Meath,Int. J. Quantum Chem., 1972,6, 501–508

  36. [36]

    B. Klahn,J. Chem. Phys., 1985,83, 5749–5753

  37. [37]

    B. Klahn,J. Chem. Phys., 1985,83, 5754–5759

  38. [38]

    R. N. Hill,Int. J. Quantum Chem., 1998,68, 357–384

  39. [39]

    Kato,Commun

    T. Kato,Commun. Pure Appl. Math., 1957,10, 151–177

  40. [40]

    R. T. Pack and W. B. Brown,J. Chem. Phys., 1966,45, 556–559

  41. [41]

    von StecherTrapped ultracold atoms with tunable interactionsPhD thesis, University of Colorado at Boulder, 2008

    J. von StecherTrapped ultracold atoms with tunable interactionsPhD thesis, University of Colorado at Boulder, 2008

  42. [42]

    von Stecher and C

    J. von Stecher and C. H. Greene,Phys. Rev. A, 2009,80, 022504

  43. [43]

    Rakshit and D

    D. Rakshit and D. Blume,Phys. Rev. A, 2012,86, 062513

  44. [44]

    A. A. S. Kalaee, Deuteron photodisintegration using quasicontinuum of p-waves in correlated gaussian basis BSc thesis, Aarhus University, 2014

  45. [45]

    P. H. Mosegaard, Deuteron photodisintegration in a shifted correlated gaus- sian basis BSc thesis, Aarhus University, 2018

  46. [46]

    Moriya, W

    H. Moriya, W. Horiuchi, and B. Zhou,Eur. Phys. J. A, 2023,59, 197

  47. [47]

    Coomar, K

    A. Coomar, K. Jones, and L. Adamowicz,Chem. Phys. Lett., 2022,790, 139358

  48. [48]

    Szabo and N

    A. Szabo and N. S. Ostlund,Modern quantum chemistry: introduction to advanced electronic structure theory, Dover Publications, 1996

  49. [49]

    Bromley and J

    M. Bromley and J. Mitroy,Int. J. Quantum Chem., 2007,107, 1150–1161

  50. [50]

    Klopper and W

    W. Klopper and W. Kutzelnigg,J. Mol. Struct.: THEOCHEM, 1986,135, 339–356

  51. [51]

    Kutzelnigg,Int

    W. Kutzelnigg,Int. J. Quantum Chem., 1994,51, 447–463

  52. [52]

    Kutzelnigg inStrategies and Applications in Quantum Chemistry: From Molecular Astrophysics to Molecular Engineering, ed

    W. Kutzelnigg inStrategies and Applications in Quantum Chemistry: From Molecular Astrophysics to Molecular Engineering, ed. Y. Ellinger and M. Defranceschi; Springer, 1996; pp. 79–101

  53. [53]

    Kutzelnigg,Oberwolfach Reports, 2011,8, 1775–1784

    W. Kutzelnigg,Oberwolfach Reports, 2011,8, 1775–1784. 13

  54. [54]

    Kutzelnigg InAIP Conf

    W. Kutzelnigg InAIP Conf. Proc., Vol. 1504, pp. 15–30. American Institute of Physics, 2012

  55. [55]

    Bakken and T

    V. Bakken and T. Helgaker,Theor. Chem. Acc., 2004,112, 124–134

  56. [56]

    L. K. McKemmish and P. M. Gill,J. Chem. Theory Comput., 2012,8, 4891–4898

  57. [57]

    Bachmayr, H

    M. Bachmayr, H. Chen, and R. Schneider,Numer. Math., 2014,128, 137– 165

  58. [58]

    R. A. Shaw,Int. J. Quantum Chem., 2020,120, e26264

  59. [59]

    Wang,Phys

    C. Wang,Phys. Rev. A, 2013,88, 032511

  60. [60]

    Jensen,J

    F. Jensen,J. Chem. Phys., 1999,110, 6601–6605

  61. [61]

    Jensen,Theor

    F. Jensen,Theor. Chem. Acc., 2005,113, 267–273

  62. [62]

    Karton and J

    A. Karton and J. M. Martin,Theor. Chem. Acc., 2006,115, 330–333

  63. [63]

    Nakatsuji and H

    H. Nakatsuji and H. Nakashima,Int. J. Quantum Chem., 2009,109, 2248– 2262. [64]https://docs.python.org/3/library/string.html# format-specification-mini-language, Accessed: 2026-03-12. [65]https://docs.python.org/3/library/functions.html#round, Ac- cessed: 2026-03-12

  64. [64]

    Goldberg,ACM Comput

    D. Goldberg,ACM Comput. Surv., 1991,23, 5–48

  65. [65]

    Kikuchi,J

    R. Kikuchi,J. Chem. Phys., 1954,22, 148–148

  66. [66]

    Shavitt and M

    I. Shavitt and M. Karplus,J. Chem. Phys., 1962,36, 550–551. [69]NIST Digital Library of Mathematical Functionshttps://dlmf.nist. gov/, Release 1.2.6 of 2026-03-15, 2026

  67. [67]

    Ammar, A

    A. Ammar, A. Leclerc, and L. U. Ancarani inAdv. Quantum Chem., Vol. 88; Elsevier, 2023; pp. 133–149

  68. [68]

    Van Rossum and J

    G. Van Rossum and J. De Boer,CWI quarterly, 1991,4, 283–303. [72]https://github.com/tkem/cachetools, Accessed: 2026-05-21

  69. [69]

    mpmath: a Python library for arbitrary-precision floating-point arithmetic (version 1.3.0). T. mpmath development team 2023

  70. [70]

    Meurer, C

    A. Meurer, C. P. Smith, M. Paprocki, O. ˇCert´ ık, S. B. Kirpichev, M. Rock- lin, A. Kumar, S. Ivanov, J. K. Moore, S. Singh, T. Rathnayake, S. Vig, B. E. Granger, R. P. Muller, F. Bonazzi, H. Gupta, S. Vats, F. Johansson, F. Pedregosa, M. J. Curry, A. R. Terrel, v. Rouˇ cka, A. Saboo, I. Fernando, S. Kulal, R. Cimrman, and A. Scopatz,PeerJ Comput. Sci., ...

  71. [71]

    Virtanen, R

    P. Virtanen, R. Gommers, T. E. Oliphant, M. Haberland, T. Reddy, D. Cournapeau, E. Burovski, P. Peterson, W. Weckesser, J. Bright, S. J. van der Walt, M. Brett, J. Wilson, K. J. Millman, N. Mayorov, A. R. J. Nelson, E. Jones, R. Kern, E. Larson, C. J. Carey, ˙I. Polat, Y. Feng, E. W. Moore, J. VanderPlas, D. Laxalde, J. Perktold, R. Cimrman, I. Henriksen,...

  72. [72]

    J. D. Hunter,Computing in Science & Engineering, 2007,9, 90–95. [77]https://docs.scipy.org/doc/scipy/reference/generated/scipy. optimize.curve_fit.html, Accessed: 2026-05-22

  73. [73]

    Bezanson, A

    J. Bezanson, A. Edelman, S. Karpinski, and V. B. Shah,SIAM Rev., 2017, 59, 65–98. [79]https://github.com/JuliaLinearAlgebra/GenericLinearAlgebra.jl, Accessed: 2026-05-21. [80]https://github.com/JuliaIO/JSON.jl, Accessed: 2026-05-21. [81]https://www.cursor.com/, Accessed: 2026-03-12

  74. [74]

    Y. I. Kurokawa and H. Nakatsuji,J. Chem. Phys., 2023,159, 024103

  75. [75]

    Nakashima and H

    H. Nakashima and H. Nakatsuji,Phys. Rev. A, 2020,102, 052835

  76. [76]

    Nakatsuji, H

    H. Nakatsuji, H. Nakashima, and Y. I. Kurokawa,Phys. Rev. A, 2020,101, 062508. 15