Tests of the Envelope Theory in One Dimension

Claude Semay; Lorenzo Cimino

Tests of the Envelope Theory in One Dimension

T0 review · 0 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read The envelope theory yields simple analytical bounds for one-dimensional many-body ground states, with accuracy set by interaction range and particle number.

desk verdict A clean, honest test of the envelope theory in 1D; the bounds hold up, and the only real weakness is a missing explicit verification of the envelope inequalities, which are actually satisfied. read the letter →

arxiv 1908.07385 v1 pith:UL533TWK submitted 2019-08-19 quant-ph

Claude Semay , Lorenzo Cimino This is my paper · ORCID

classification quant-ph

keywords envelopetheoryone-dimensionalmany-bodysystemsCalogeromodelGaussianpotentialground-stateenergyboundsidenticalparticlesLagrange-meshmethodLambertWfunction

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 tests whether the envelope theory—replace each kinetic term and pairwise potential by a tangent quadratic, solve the resulting harmonic-oscillator problem, then optimize—can provide trustworthy ground-state energies for one-dimensional systems of identical particles. For the Calogero model of N fermions with linear plus inverse-cube forces, it establishes a rigorous lower bound from a closed formula. For up to 100 bosons interacting through a Gaussian potential, it obtains an analytical upper bound whose relative error decreases as the particle number or the potential range increases. The practical payoff is that a few algebraic equations can bracket exact or accurately computed energies in regimes where direct numerical work is heavy.

What carries the argument

The load-bearing object is the envelope pair: a quadratic $\tilde{V}(x)=c_1 x^2+c_2$ tangent to the true potential $V$ at at least one point, and a quadratic $\tilde{T}(p)=d_1 p^2+d_2$ tangent to the kinetic term $T$. Replacing $H$ by $\tilde{H}$ makes the $N$-body problem exactly solvable, and optimizing the tangency parameters yields equations (5)–(7) together with $x_0 p_0 = Q$, where $Q$ takes the bosonic or fermionic ground-state value. The resulting energies are rigorous bounds whenever the two quadratics stay on one side of the true curves for all $x$ and $p$; this pointwise comparison is what turns a variational estimate into a theorem. For the Gaussian potential the optimized formula is closed in terms of the Lambert $W$ function.

What would settle it

Take the Gaussian-boson system at a parameter set where (14) returns a real finite value, compute the ground-state energy independently on a fine grid or with the Lagrange-mesh method of [12], and test whether $E_{\mathrm{ET}} < E_{\mathrm{LM}}$; one such violation would disprove the claimed upper-bound character, while a direct check of $\tilde{V}(x)-V(x)\ge 0$ and $\tilde{T}(p)-T(p)\ge 0$ at the optimal envelope parameters would show whether the theorem's precondition actually holds.

Watch

Extended reading notes

Core claim

The central claim is that in one dimension the envelope equations (5)–(7) produce quantitative bounds, not just estimates. For fermionic Calogero ground states, the envelope energy $E_{\mathrm{ET}}$ is always below the exact energy $E_{\mathrm{ex}}$ ($E_{\mathrm{ET}}^2 < E_{\mathrm{ex}}^2$), coincides with $E_{\mathrm{ex}}$ at zero inverse-cube coupling, and its relative error saturates at the value (12) as $N$ grows. For bosons in a Gaussian well, the envelope energy (14), expressed with the principal Lambert $W$ branch, lies above the accurate Lagrange-mesh energies when the potential is not too singular, and the comparison in Table I shows the error falling from about 0.41 at $(N,a)=(20,1.0)$ to about 0.054 at $(100,1.0)$. The paper also notes that for small $N$ or small range the method can return complex or positive energies, in which case no bound exists.

Load-bearing premise

The load-bearing premise is that the bounding theorem of the envelope method—pointwise inequalities $\tilde{V}(x)\ge V(x)$ and $\tilde{T}(p)\ge T(p)$ give an upper bound, and reversed inequalities a lower bound—holds for these one-dimensional systems; for the Gaussian potential the paper does not verify those inequalities explicitly at the adopted parameter values, and the appearance of complex or positive energies shows the preconditions are not automatic.

Editorial extensions

If this is right

For the Calogero model, the envelope lower bound is exact at $g'=0$ and remains accurate for $g'<1$, with the limiting relative error given by (12).
For Gaussian bosons, the upper bound's relative error falls with both $N$ and the range $a$, from 0.41 at $(N,a)=(20,1.0)$ to 0.054 at $(100,1.0)$.
The closed-form envelope energies can serve as a quick independent check on accurate numerical methods for one-dimensional many-body systems.
In parameter regions where complex or positive energies appear, the method simply supplies no relevant bound, so the variational character is not universal.

Reading between the lines

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

A dimensionless measure of short-distance steepness, such as $a|V'(0)|/|V(0)|$, could predict when the Gaussian bound will be reliable; the paper only documents the trend without proposing such a predictor.
In one dimension the absence of angular momentum blocks the quantum-number improvement available in $D\ge 2$, so these results expose the baseline envelope construction rather than an optimized variant.
Extracting the large-$N$ asymptotic form of the Lambert-$W$ expression (14) could give closed-form estimates for many-boson systems beyond the $N=100$ benchmark considered here.

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.

Desk Editor's Note

read the letter

Hey —

The thing to know: this is a modest, honest test of an existing approximation method, not a breakthrough. The authors take the envelope theory, which they helped develop, and check it against two exactly/numerically solvable one-dimensional N-body systems. The new bits are the Calogero lower-bound formula (Eq. 10) and the application of their earlier Gaussian upper bound (Eq. 14) up to N=100. The comparisons are clean.

What's good: the math is straightforward and the bound property for Calogero is shown algebraically (E_ex^2 > E_ET^2). For the Gaussian case, they compare against independent Lagrange-mesh results from Timofeyuk and Baye, and the error table is honest. The paper even flags parameter regions where the method returns complex or positive 'irrelevant' energies—good scientific hygiene. No fitted parameters, no circularity. The self-citations are to the method's foundations, which is appropriate.

Where it's soft: the one thing a skeptical reader might trip on is that the paper asserts the pointwise envelope inequalities (Section I) justify the bounds, but never verifies them explicitly for the Gaussian potential at the parameter values used. I checked, and the inequalities do hold—the tangent envelope lies above V(x) because e^{-t} is convex, so Eq. (14) is a legit upper bound wherever the Lambert branch is real. The Calogero side is even easier. So this is a presentation gap, not a flaw. The 'irrelevant values' are a real-branch issue of the Lambert equation, and the authors acknowledge it.

Who it's for: envelope-theory practitioners and people doing few-body numerics who want a sanity check. It's not going to change the field, but it's a clean, useful calibration of a method that is indeed very simple to implement.

On peer review: yes, send it out. It deserves a serious referee—short, verifiable, and honest. I wouldn't cite it unless I was working on this exact method, but it's a legitimate contribution.

Referee Report

0 major / 5 minor

Summary. The manuscript reports tests of the envelope theory (ET) for identical particles in one spatial dimension. It treats two systems: the Calogero model with linear plus inverse-cube pair potentials, for which an exact fermionic ground state is known, and a Gaussian pair potential for up to N=100 bosons, for which accurate Lagrange-mesh results from Ref. [11] are available. The authors derive the ET approximations (Eqs. 10 and 14), compare them with the reference results (Figs. 1-2 and Table I), and conclude that good bounds can be obtained for a parameter-dependent range, while acknowledging the limitations of the method.

Significance. If the bounds are valid, the paper provides a useful extension of ET to D=1, where previous applications were mostly in D>=3. Its strengths are the use of external benchmarks (exact Calogero solution and independent Lagrange-mesh calculations), fully analytic formulas, and an honest discussion of the parameter regimes in which the bounds are reliable. The central claim is moderate but sound: ET gives quick, certified approximate energies for many-body systems in simple cases. The paper is reproducible and self-contained modulo the cited theorems.

minor comments (5)

[Section I / Section III (Eq. 14)] The upper-bound property of Eq. (14) relies on the pointwise envelope condition tilde V(x) >= V(x) for the Gaussian potential, but the manuscript does not verify this condition. I checked it: with u = x0^2/a^2, the tangent envelope gives tilde V(x) - V(x) = Vg[exp(-x^2/a^2) - exp(-u)(1+u-x^2/a^2)] >= 0 by convexity of exp(-t). Please include this short derivation or an explicit reference so that the bound is certified within the paper.
[Section II (Eqs. 9-10)] The statement that E_ex^2 > E_ET^2 is only described as 'easily checked'. Since the lower-bound character is a central result, please include the algebraic comparison or a footnote with the explicit expansion, so that the reader does not have to reproduce the algebra.
[Section III / Table I] The 'irrelevant values' (complex or positive energies) that occur for small N or small a are acknowledged but not characterized precisely. Please specify the parameter range in which the Lambert-W branch W0 gives a real, negative energy, and state explicitly when Eq. (14) ceases to provide a valid upper bound.
[Section I] The introduction says the method was previously developed for D > 2 dimensions. Please state explicitly that the D=1 formulas follow from Refs. [4,6] and summarize any changes in the definitions of x, p, and Q, since this is the first D=1 application.
[Eqs. (9) and (10)] The typeset formulas are ambiguous due to missing parentheses. Please rewrite them with clear brackets, e.g., using explicit fraction layout, so that the reader can unambiguously parse the expressions.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the envelope theory predictions are benchmarked against the exact Calogero solution and independent Lagrange-mesh calculations, with no parameters fitted to the reference data.

full rationale

The paper's central claims are the envelope-theory energies for two one-dimensional systems: the fermionic Calogero ground state and the bosonic Gaussian ground state. These predictions are compared against external benchmarks: the exact Calogero solution from Ref. [10] and the Lagrange-mesh results of Ref. [11]. No parameter appearing in the envelope-theory formulas is fitted to those benchmarks; the only inputs are the Hamiltonian parameters and the previously derived envelope-theory equations (5)-(7) from Refs. [4,6] and the Gaussian bound formula (14) from Ref. [7]. These self-citations are not load-bearing in a circular sense, because the cited prior work derives general formulas for the envelope theory rather than tuning them to the present systems. The bound character of the results is asserted via the envelope theorem of Hall (Refs. [1,2]) and, although the paper does not explicitly verify the pointwise inequalities, an independent check shows they hold: for the Calogero potential, V(x) minus its envelope is nonnegative for the lower bound, and for the Gaussian potential the convexity of exp(-x^2/a^2) guarantees the envelope lies above the potential for the upper bound wherever a real solution exists. The paper honestly acknowledges that for small N or small a the method returns irrelevant complex or positive energies, and Table I reports the relative errors relative to external data. Therefore the derivation chain does not reduce to its inputs by construction, and no circular step is present.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The paper introduces no new free parameters or invented entities. It relies on the envelope theory from prior work (assuming its validity in 1D), the exact Calogero solution, and the Lagrange-mesh numerical results as external benchmarks.

assumptions (4)

domain assumption The envelope theory equations (5)-(7), originally derived for D>2, remain valid in D=1 with scalar momenta and positions and with ground-state quantum numbers Q_B0 and Q_F0 from Eqs. (2)-(3).
Stated in Section I with reference to Refs. [4,6]; no proof is given for the 1D extension, but it is a direct transcription.
domain assumption The exact Calogero ground-state energy E_ex in Eq. (9) is the correct reference for the fermionic ground state.
Quoted from Calogero (Ref. [10]) and used as the external benchmark in Section II.
domain assumption The Lagrange-mesh results of Ref. [11] provide accurate reference ground-state energies for the Gaussian bosonic systems.
Adopted as ground truth in Section III without independent verification in this paper.
domain assumption For the computed energies to be rigorous bounds, the envelope functions must satisfy the pointwise inequalities described in Section I: tilde V(x) >= V(x) (for upper bounds) or <= (for lower bounds), and similarly for T(p).
The paper relies on the theorem from Refs. [2,4] but does not verify the inequalities pointwise for the Gaussian parameters; this is the main fragility.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Tests of the Envelope Theory in One Dimension." pith.science (2026). https://pith.science/paper/UL533TWK

@misc{pith2026190807385,
  author       = {Pith},
  title        = {Pith review of: Tests of the Envelope Theory in One Dimension},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/UL533TWK}},
  note         = {Machine review of arXiv:1908.07385}
}

read the original abstract

The envelope theory is a simple technique to obtain approximate, but reliable, solutions of many-body systems with identical particles. The accuracy of this method is tested here for two systems in one dimension with pairwise forces. The first one is the fermionic ground state of the analytical Calogero model with linear forces supplemented by inverse-cube forces. The second one is the ground state of up to 100 bosons interacting via a Gaussian potential. Good bounds can be obtained depending on values of the model parameters.

Figures

Figures reproduced from arXiv: 1908.07385 by the authors.

**Figure 2.** FIG. 2: Energy in K for the Gaussian potential as a function of [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

13 extracted references · 13 canonical work pages

[11]

Timofeyuk and D

N.K. Timofeyuk and D. Baye, Hyperspherical Harmonics E xpansion on Lagrange Meshes for Bosonic Systems in One Dimension. Few-Body Syst. 58, 157 (2017)

work page 2017
[1]

Hall, Energy trajectories for the N -boson problem by the method of potential envelopes

R.L. Hall, Energy trajectories for the N -boson problem by the method of potential envelopes. Phys. R ev. D 22, 2062 (1980)

work page 1980
[2]

Hall, A geometrical theory of energy trajectories i n quantum mechanics

R.L. Hall, A geometrical theory of energy trajectories i n quantum mechanics. J. Math. Phys. 24, 324 (1983)

work page 1983
[3]

R.L. Hall, W. Lucha, and F.F. Sch¨ oberl, Relativistic N -boson systems bound by pair potentials V (rij) = g(r2 ij). J. Math. Phys. 45, 3086 (2004)

work page 2004
[4]

Semay and C

C. Semay and C. Roland, Approximate solutions for N -body Hamiltonians with identical particles in D dimensions. Res. Phys. 3, 231 (2013)

work page 2013
[5]

Semay and F

C. Semay and F. Buisseret, Bound Cyclic Systems with the E nvelope Theory. Few-Body Syst. 58, 151 (2017)

work page 2017
[6]

Semay and G

C. Semay and G. Sicorello, Many-Body Forces with the Enve lope Theory. Few-Body Syst. 59, 119 (2018)

work page 2018
[7]

Semay, Numerical Tests of the Envelope Theory for Few- Boson Systems

C. Semay, Numerical Tests of the Envelope Theory for Few- Boson Systems. Few-Body Syst. 56, 149 (2015)

work page 2015

Show all 13 references

[8]

Semay, Improvement of the envelope theory with the dom inantly orbital state method

C. Semay, Improvement of the envelope theory with the dom inantly orbital state method. Eur. Phys. J. Plus 130, 156 (2015)

2015
[9]

L´ evy-Leblond, Generalized uncertainty relations for many-fermion system

J.M. L´ evy-Leblond, Generalized uncertainty relations for many-fermion system. Phys. Lett. A 26, 540 (1968)

1968
[10]

Calogero, Ground State of a One-Dimensional N -Body System

F. Calogero, Ground State of a One-Dimensional N -Body System. J. Math. Phys. 10, 2197 (1969)

1969
[12]

Baye, The Lagrange-mesh method

D. Baye, The Lagrange-mesh method. Phys. Rep. 565, 1 (2015)

2015
[13]

Corless, G.H

R.M. Corless, G.H. Gonnet, D.E.G. Hare, D.J. Jeﬀrey, D. E. Knuth, On the Lambert W Function. Adv. Comput. Math. 5, 329 (1996)

1996

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [11]

Timofeyuk and D

N.K. Timofeyuk and D. Baye, Hyperspherical Harmonics E xpansion on Lagrange Meshes for Bosonic Systems in One Dimension. Few-Body Syst. 58, 157 (2017)

work page 2017

[2] [1]

Hall, Energy trajectories for the N -boson problem by the method of potential envelopes

R.L. Hall, Energy trajectories for the N -boson problem by the method of potential envelopes. Phys. R ev. D 22, 2062 (1980)

work page 1980

[3] [2]

Hall, A geometrical theory of energy trajectories i n quantum mechanics

R.L. Hall, A geometrical theory of energy trajectories i n quantum mechanics. J. Math. Phys. 24, 324 (1983)

work page 1983

[4] [3]

R.L. Hall, W. Lucha, and F.F. Sch¨ oberl, Relativistic N -boson systems bound by pair potentials V (rij) = g(r2 ij). J. Math. Phys. 45, 3086 (2004)

work page 2004

[5] [4]

Semay and C

C. Semay and C. Roland, Approximate solutions for N -body Hamiltonians with identical particles in D dimensions. Res. Phys. 3, 231 (2013)

work page 2013

[6] [5]

Semay and F

C. Semay and F. Buisseret, Bound Cyclic Systems with the E nvelope Theory. Few-Body Syst. 58, 151 (2017)

work page 2017

[7] [6]

Semay and G

C. Semay and G. Sicorello, Many-Body Forces with the Enve lope Theory. Few-Body Syst. 59, 119 (2018)

work page 2018

[8] [7]

Semay, Numerical Tests of the Envelope Theory for Few- Boson Systems

C. Semay, Numerical Tests of the Envelope Theory for Few- Boson Systems. Few-Body Syst. 56, 149 (2015)

work page 2015