Pith. sign in

REVIEW 2 major objections 6 minor 24 references

Preparing Small Gaussian Basis for Highly Accurate Ab Initio Description of Lithium Rydberg States

T0 review · 2 major / 6 minor · reviewed 2026-07-30 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A compact tailored Gaussian basis recovers lithium Rydberg states to better than 0.01 eV.

desk verdict Solid methods paper: a smaller ExTG-contracted Li Rydberg basis that hits few-meV EOM-CCSD accuracy through n=7 S–G and beats a Kaufmann-style universal set at lower cost. read the letter →

arxiv 2607.26728 v1 pith:OG37RFCK submitted 2026-07-29 physics.chem-ph physics.atom-ph

Jan Šmydke , Benjamín Andreides , Simona Dubcová This is my paper · ORCID
classification physics.chem-phphysics.atom-ph
keywords RydbergstatesGaussianbasissetslithiumEOM-CCSDexponentiallytemperedGaussiansKoopmanstheorembasis-setoptimizationdiffuseorbitals
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper builds a small Gaussian basis that accurately describes the Rydberg excited states of atomic lithium up to principal quantum number n = 7 and angular momenta S through G. The authors variationally optimize exponentially tempered diffuse Gaussians on the virtual orbitals of Li+, contract them, and add the resulting few functions to a standard correlation basis. At the EOM-CCSD level the excitation energies and ionization potential then agree with the best available non-relativistic references to a few meV, while the basis stays smaller than both an earlier even-tempered set and a comparable universal Rydberg Gaussian set. One-dimensional orbital plots show a clean, regular nodal pattern on a logarithmic radial scale that extends tens of angstroms from the nucleus. The protocol is presented as a practical step toward routine ab initio Rydberg work on larger atoms and molecules without enormous diffuse bases.

What carries the argument

Exponentially tempered Gaussian (ExTG) primitives, variationally optimized as Hartree–Fock virtual orbitals of closed-shell Li+ under Koopmans’ theorem and then contracted into a handful of Rydberg basis functions that are simply added to aug-ANO-pVQZ.

What would settle it

Recompute the same ROHF-EOM-CCSD excitation energies for the higher-n F and G states against a large explicitly correlated or STO Full-CI reference; if the few-meV systematic offset grows or reverses while the universal basis does not, the claimed accuracy and superiority of the tailored ExTG contraction fail.

Watch

Extended reading notes

Core claim

At the EOM-CCSD level the new AAQZ-ExTG-Ry7G basis, and its still smaller state-selective form AAQZ-ExTG-Ry3P plus one extra nL function, yields lithium Rydberg excitation energies up to n = 7 (S–G) and the ionization potential accurate to better than 10^{-2} eV, outperforming an analogous universal Kaufmann-style Gaussian Rydberg basis while using fewer functions.

Load-bearing premise

That Hartree–Fock virtual orbitals optimized on the Li+ cation remain accurate enough contracted functions for the correlated Rydberg states of neutral lithium once few-meV accuracy is required.

Editorial extensions

If this is right

  • State-selective high-n calculations need only two extra S, two extra P, and one nL contracted Rydberg function on top of a standard basis.
  • The same optimization protocol can generate compact Rydberg bases for larger atoms.
  • Properly tempered and contracted Gaussians can describe orbitals that extend tens of angstroms (and, with the earlier larger set, even to ~10 µm).
  • The ExTG form reduces the near-linear-dependence problems that appear with ordinary even-tempered diffuse sets.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the HF-to-EOM-CCSD transfer holds for open-shell molecular cores, the same few-function augmentation could make Rydberg PESs routine in polyatomics without custom continuum bases.
  • The observed logarithmic regularity of nodes suggests a simple radial scaling that could supply starting ExTG parameters for other atoms without full re-optimization.
  • Breakdown of the map under strong core polarization would likely appear first as n-dependent error growth in the ionization-potential limit rather than in the lowest valence states.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 6 minor

Summary. The manuscript develops a compact Gaussian basis (AAQZ-ExTG-Ry7G and the state-selective AAQZ-ExTG-Ry3P+nL) for 2Li Rydberg states up to n=7 and angular momenta S–G. Building on a Koopmans/HF virtual-orbital optimization of Li+ with exponentially tempered Gaussian (ExTG) primitives added to aug-ANO-pVQZ, the authors report ROHF-EOM-CCSD excitation energies and a CCSD ionization potential accurate to a few meV (better than 10^{-2} eV) relative to high-level non-relativistic references and experiment. The new bases are compared favorably to the authors’ earlier ETG-based 2023Sm set and to a Kaufmann-style universal primitive set of comparable scope, while using substantially fewer primitives and contracted functions. One-dimensional orbital cuts are used to illustrate a regular logarithmic nodal structure of the Rydberg orbitals out to tens of Å (and, with the prior basis, to ~10 µm for 25S).

Significance. If the tabulated accuracy and compactness hold, the work is a practical methodological contribution to ab initio Rydberg chemistry: it shows that carefully contracted, system-tailored Gaussians can describe high-n Rydberg states of Li at few-meV accuracy with basis sizes compatible with standard quantum-chemistry codes, including a minimal state-selective construction that adds only a handful of functions. The external benchmarks against ECG/Hy-CI-type references, experiment, and a Kaufmann-style universal set, together with the explicit size comparison in Table I, are genuine strengths. The ExTG parameterization and the documented reduction of near-linear-dependence issues relative to ETG are useful technical results. The molecular-extension framing is appropriately prospective and does not underwrite the Li claim.

major comments (2)
  1. [§III.B, Tables II–IV] §III.B and Tables II–III: The residual errors for S–F states systematically approach ~−3.5 meV and track the ~4 meV IP underestimation (Table IV), while G states show an opposite-sign ~+6 meV residual. The manuscript attributes the overall few-meV accuracy to the basis at EOM-CCSD, but does not cleanly separate basis incompleteness from EOM-CCSD/correlation (and missing relativity) error. A short, explicit discussion—e.g., that the common offset matches the IP error and is therefore largely method- rather than basis-limited, and why G differs—would make the central “better than 10^{-2} eV” claim more precise and less easy to over-read as pure basis quality.
  2. [§II.A–C; Abstract; §IV] §II.A–C and the transferability claim: The contracted functions are optimized exclusively on closed-shell RHF virtual energies of Li+ under Koopmans’ theorem, then used for ROHF-EOM-CCSD of neutral 2Li. Tables II–IV already show that this map works to a few meV for Li, so the central atomic claim stands. However, the abstract and §IV present the protocol as “an important step towards” polyatomic Rydberg work without stating what would have to be re-validated for open-shell cores or multi-center systems (orbital relaxation, core–Rydberg correlation, choice of ion-core reference). A brief, concrete caveat on the scope of the demonstrated transfer would keep the prospective claim proportionate to the evidence.
minor comments (6)
  1. [Abstract; §IV] Abstract and §IV: the phrase “higher than 10^{-2} eV accuracy” is ambiguous (larger accuracy vs. accuracy better than 10^{-2} eV). Prefer “better than 10^{-2} eV” or “errors below 10^{-2} eV,” consistent with the body text.
  2. [§I] §I: “but my need to obtain just a few” → “may need”.
  3. [§III] §III: “thoroughfully” → “thoroughly”.
  4. [Table I; §III.A] Table I footnote and basis naming: AAQZ-ExTG-Ry3P is described as adding two S and two P contracted Rydberg functions; ensure the table’s contracted counts and the “+ nL” usage in the text/tables are fully consistent so a reader can reconstruct the state-selective sets without the GENBAS file.
  5. [§III.C; Figs. 1–3] Figures 1–3 and SI: the logarithmic radial cuts are central to the nodal-structure claim; adding a brief note on normalization/scale of the plotted amplitude (and that only the m-aligned component is shown for ℓ>0) would aid reproducibility.
  6. [§III.B] References: several high-precision Li Rydberg benchmarks are cited; a single sentence placing the present EOM-CCSD/basis errors relative to residual relativistic/QED effects (already partly in the IP comparison) would help non-specialist readers.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: HF-Koopmans basis fit is validated by independent EOM-CCSD benchmarks, not by the fit objective

full rationale

The paper variationally optimizes ExTG Gaussian exponents to minimize Li+ RHF virtual orbital energies (Koopmans map), contracts those functions, then reports ROHF-EOM-CCSD excitation energies and the CCSD ionization potential. Those correlated observables are compared to external non-relativistic references (ECG, Hylleraas-CI, STO Full-CI) and experiment (Tables II–IV), not to the HF fit target. Superiority versus the Kaufmann universal set and versus the authors’ prior ETG basis is likewise an external size/accuracy comparison. Adoption of the ExTG form from the authors’ earlier helium work and comparison to their 2023Sm basis are ordinary methodological self-citation; neither forces the few-meV accuracy claim by construction. The derivation chain is therefore self-contained against outside benchmarks, with no self-definitional loop, no fitted-input-as-prediction, and no load-bearing uniqueness theorem imported from the authors.

Assumptions & free parameters 3 free parameters · 5 assumptions · 1 invented entities

The central accuracy claim rests on standard electronic-structure machinery plus the modeling choice that optimized Li⁺ HF virtuals are good contracted Rydberg functions for neutral Li EOM-CCSD, the ExTG exponent family, and the fixed aug-ANO-pVQZ correlation backbone. Free parameters are the ExTG coefficients optimized under orbital-energy constraints; no new physical entities are postulated.

free parameters (3)
  • ExTG series parameters {ζ0,i, αi, βi} per angular momentum (and primitive count N) = Not numerically tabulated in main text; implied by contracted GENBAS SI
    Variationally tuned to minimize successive Rydberg virtual energies of Li⁺ while holding E_HF and lower Rydberg energies fixed; they define the new basis and are not predicted a priori.
  • Choice of n_max=7 and ℓ≤G optimization target = n≤7, S–G
    Hand-chosen scope that sets which contracted Rydberg functions exist and what ‘full’ vs state-selective means.
  • Orbital-energy convergence threshold 10^{-9} a.u. = 10^{-9} a.u.
    Optimizer stopping criterion that affects how tightly exponents are determined.
assumptions (5)
  • domain assumption Koopmans’ theorem: negative canonical RHF virtual energies of Li⁺ approximate electron affinities / Rydberg orbital energies of the ion+electron system.
    Stated as the construction principle in §II.A; load-bearing for using HF virtual LCAOs as contracted Rydberg basis functions.
  • domain assumption Rydberg–core correlation is weak enough that HF-shaped Rydberg orbitals remain accurate when used in ROHF-EOM-CCSD on neutral ²Li.
    Explicitly argued in §II.A; required to equate basis quality for HF virtuals with correlated excitation accuracy.
  • domain assumption aug-ANO-pVQZ adequately describes valence/core correlation so residual excitation errors are dominated by the diffuse Rydberg set.
    Protocol step 1 (§II.C); all production energies keep this fixed backbone.
  • ad hoc to paper Exponentially tempered Gaussian (ExTG) exponent series (Eqs. 3–5) is a flexible, near-linear-dependence-resistant primitive set for Rydberg orbitals.
    Adopted from authors’ prior He work and chosen over ETG; two-term expansion used throughout results.
  • domain assumption Standard non-relativistic EOM-CCSD / CCSD energy differences are valid comparators to the cited ‘exact non-relativistic estimates’ for assessing basis error at the ~meV level.
    Tables II–IV rest on this method and reference hierarchy; relativity/QED are outside scope.
invented entities (1)
  • AAQZ-ExTG-Ry7G / AAQZ-ExTG-Ry3P+nL contracted Rydberg basis independent evidence
    purpose: Provide small, high-accuracy Gaussian functions for ²Li Rydberg states n≤7, S–G, including state-selective use.
    The paper’s primary deliverable: named contracted sets built from optimized ExTG primitives on aug-ANO-pVQZ.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Preparing Small Gaussian Basis for Highly Accurate Ab Initio Description of Lithium Rydberg States." pith.science (2026). https://pith.science/paper/OG37RFCK

@misc{pith2026260726728,
  author       = {Pith},
  title        = {Pith review of: Preparing Small Gaussian Basis for Highly Accurate Ab Initio Description of Lithium Rydberg States},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OG37RFCK}},
  note         = {Machine review of arXiv:2607.26728}
}
abstract

A new small and highly accurate Gaussian basis set has been developed for the ab initio description of $^2$Li Rydberg excited states up to $n = 7$ and S, P, D, F, and G angular momentum symmetry, and the appropriate optimization protocol is presented. The obtained Rydberg excitation energies are compared with results from other highly accurate approaches. At the EOM-CCSD level of theory the new basis exhibits higher than $10^{-2}$ eV accuracy of the excitation energies and the ionization potential and provides superior results than an analogous universal Gaussian basis set, while utilizing even smaller number of basis functions. Plots of the 1-dimensional Rydberg-orbital cuts reveal a regular nodal structure along the logarithmic scale of the atomic radius reaching tens of angstrom far from the nucleus. The presented Rydberg basis set generation methodology is an important step towards routine ab initio Rydberg-state related investigations of more complex systems, such as large atoms and polyatomic molecules.

Figures

Figures reproduced from arXiv: 2607.26728 by the authors.

Figure 1
Figure 1. FIG. 1. Rydberg S orbitals of the [PITH_FULL_IMAGE:figures/full_fig_p019_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Red: Rydberg P [PITH_FULL_IMAGE:figures/full_fig_p020_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Rydberg 25S orbital of the [PITH_FULL_IMAGE:figures/full_fig_p021_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 8 canonical work pages

  1. [1]

    Using Koopmans’ theorem for constructing basis sets: approaching high Rydberg excited states of lithium with a compact Gaussian basis

    Šmydke, Jan. Using Koopmans’ theorem for constructing basis sets: approaching high Rydberg excited states of lithium with a compact Gaussian basis. Phys. Chem. Chem. Phys. 2023. doi:10.1039/D2CP04633D

    work page doi:10.1039/d2cp04633d 2023

  2. [2]

    The 2S Rydberg series of the lithium atom

    Amir Bralin and Sergiy Bubin and Monika Stanke and Ludwik Adamowicz , abstract =. The 2S Rydberg series of the lithium atom. Calculations with all-electron explicitly correlated Gaussian functions , journal =. 2019 , issn =. doi:https://doi.org/10.1016/j.cplett.2019.06.051 , url =

    work page doi:10.1016/j.cplett.2019.06.051 2019

  3. [3]

    , title =

    Plakhutin,Boris N. , title =. The Journal of Chemical Physics , volume =. 2018 , doi =

    2018

  4. [4]

    and Davidson,Ernest R

    Plakhutin,Boris N. and Davidson,Ernest R. , title =. The Journal of Chemical Physics , volume =. 2014 , doi =

    2014

  5. [5]

    Gaussian basis sets for highly excited and resonance states of helium , Author =. J. Chem. Phys. , Year =

  6. [6]

    Excitation of helium Rydberg states and doubly excited resonances in strong extreme ultraviolet fields: Full-dimensional quantum dynamics using exponentially tempered Gaussian basis sets , Author =. J. Chem. Phys. , Year =

  7. [7]

    Hylleraas-configuration-interaction analysis of the low-lying states in the three-electron Li atom and Be

    Ruiz, Mar\'. Hylleraas-configuration-interaction analysis of the low-lying states in the three-electron Li atom and Be. Phys. Rev. A , volume =. 2013 , month =. doi:10.1103/PhysRevA.88.012505 , url =

    work page doi:10.1103/physreva.88.012505 2013

  8. [8]

    Variational energies and the Fermi contact term for the low-lying states of lithium: Basis-set completeness , author =. Phys. Rev. A , volume =. 2012 , month =. doi:10.1103/PhysRevA.85.052513 , url =

    work page doi:10.1103/physreva.85.052513 2012

Show all 24 references
  1. [9]

    Bubin,Sergiy and Adamowicz,Ludwik , title =. J. Chem. Phys. , volume =. 2012 , doi =

    2012

  2. [10]

    and Bubin, Sergiy and Adamowicz, Ludwik , journal =

    Sharkey, Keeper L. and Bubin, Sergiy and Adamowicz, Ludwik , journal =. Lower Rydberg. 2011 , month =. doi:10.1103/PhysRevA.83.012506 , url =

    2011 doi

  3. [11]

    , title =

    Neese, Frank and Valeev, Edward F. , title =. J. Chem. Theory Comput. , volume =. 2011 , doi =

    2011

  4. [12]

    and Plakhutin,Boris N

    Davidson,Ernest R. and Plakhutin,Boris N. , title =. The Journal of Chemical Physics , volume =. 2010 , doi =

    2010

  5. [13]

    Krylov , title =

    Hanna Reisler and Anna I. Krylov , title =. Int. Rev. Phys. Chem. , volume =. 2009 , publisher =

    2009

  6. [14]

    and Hagstrom, Stanley A

    Sims, James S. and Hagstrom, Stanley A. , journal =. Hylleraas-configuration-interaction study of the 2. 2009 , month =. doi:10.1103/PhysRevA.80.052507 , url =

    2009 doi

  7. [15]

    Relativistic, QED, and finite nuclear mass corrections for low-lying states of Li and

    Puchalski, Mariusz and Pachucki, Krzysztof , journal =. Relativistic, QED, and finite nuclear mass corrections for low-lying states of Li and. 2008 , month =. doi:10.1103/PhysRevA.78.052511 , url =

    2008 doi

  8. [16]

    Plakhutin,B. N. and Gorelik,E. V. and Breslavskaya,N. N. , title =. J. Chem. Phys. , volume =

  9. [17]

    King , editor =

    Frederick W. King , editor =. High-Precision Calculations for the Ground and Excited States of The Lithium Atom , series =. 1999 , issn =. doi:https://doi.org/10.1016/S1049-250X(08)60111-1 , url =

    1999 doi

  10. [18]

    doi:10.1088/0953-4075/22/14/007 , year =

    K Kaufmann and W Baumeister and M Jungen , title =. doi:10.1088/0953-4075/22/14/007 , year =

    doi

  11. [19]

    Attila Szabo and Neil S. Ostlund. Modern Quantum Chemistry. 1996

    1996

  12. [20]

    Martin, W. C. and Zalubas, R. , editor =. CRC Handbook of Chemistry and Physics , chapter =. 1992 , edition =

    1992

  13. [21]

    Multidimensional Optimizer for Command Line. See http://mdoptcli.sf.net

    mdoptcli , note = "Multidimensional Optimizer for Command Line. See http://mdoptcli.sf.net", url =

  14. [22]

    Galassi and Jim Davies and James Theiler and Brian Cough and Gerard Jungman and Patrick Alken and Michael Booth and Fabrice Rossi

    M. Galassi and Jim Davies and James Theiler and Brian Cough and Gerard Jungman and Patrick Alken and Michael Booth and Fabrice Rossi. GNU Scientific Library Reference Manual, 3rd Ed

  15. [23]

    Kállay and P

    M. Kállay and P. R. Nagy and D. Mester and Z. Rolik and G. Samu and J. Csontos and J. Csóka and P. B. Szabó and L. Gyevi-Nagy and B. Hégely and I. Ladjánszki and L. Szegedy and B. Ladóczki and K. Petrov and M. Farkas and P. D. Mezei and A. Ganyecz , title =. J. Chem. Phys. , v...

    2020

  16. [24]

    Barca,Giuseppe M. J. and Bertoni,Colleen and Carrington,Laura and Datta,Dipayan and De Silva,Nuwan and Deustua,J. Emiliano and Fedorov,Dmitri G. and Gour,Jeffrey R. and Gunina,Anastasia O. and Guidez,Emilie and Harville,Taylor and Irle,Stephan and Ivanic,Joe and Kowalski,Karol...

    2020

Pith tools

Reviewed July 30, 2026 · model on record in the stance chip above.