Pith. Sign in

REVIEW

Gravitating semirelativistic N-boson systems

Mathematical Physics · Submitted 25 Feb 2006 · 37 references

Pith's verdict: The upper bound is solid and the lower bound is plausible, but the main derivation in Sec. 5 has a load-bearing gap: the constrained minimizer is not an eigenfunction, so Eq. (10) is unproven as written.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 28, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 4 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDmath-ph/0602059 (version 2)
Pith IDUYGJRBW5, Pith's permanent record of this paper and its review
CitationRichard L. Hall, Wolfgang Lucha. "Gravitating semirelativistic N-boson systems." arXiv:math-ph/0602059 (2006). https://pith.science/paper/UYGJRBW5
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

For N identical bosons with relativistic kinetic energy and gravitational -v/r interactions, the ground-state energy satisfies a closed-form lower bound, E ≥ Nm sqrt((1+sqrt(1-γv^2))/2), with γ = N(N-1)/2 and γv^2 < 1.

The reading

Gravitating semirelativistic N-boson systems—N identical spin-zero particles with kinetic energy $\sqrt$($p^{2}$+$m^{2}$) attracting each other through -v/r—are hard to solve directly. This paper shows that boson permutation symmetry, expressed in Jacobi relative coordinates, reduces the problem to an effective one-body operator, and from there yields a closed-form lower bound on the ground-state energy: E ≥ Nm $\sqrt$((1+$\sqrt$(1-$γv^{2}$))/2) with γ=N(N-1)/2. A scale-optimized Gaussian trial function gives an accompanying upper bound. When the coupling is scaled as v=c/N to avoid collapse as N grows, both bounds stay meaningful for arbitrarily many bosons and reproduce the known nonrelativistic parabolic bounds in the weak-coupling limit. The result matters because it turns a genuinely many-body relativistic problem into simple formulas useful for estimating bound-state energies.

What carries the argument

The load-bearing device is the Jacobi-coordinate reduction of the N-boson problem combined with the identity <(p+δP)^2>=2λ<$p^{2}$>, where λ=(N-1)/N and δ=$\sqrt$((N-2)/N). This identity follows from the vanishing of momentum cross-correlations between different Jacobi directions in the boson-symmetric ground state, and it is what lets the authors square the kinetic term N $\sqrt$((p+δP)^2+$m^{2}$) and re-express the expectation value as an effective one-body operator. Once squared, the potential produces a Kratzer-type operator—a Coulomb term -vE/r plus a 1/$r^{2}$ spike—whose known one-body lower bound supplies the N-body bound (10). The upper bound uses Gaussian trial wave functions in relative coordinates, whose width is scaled as a variational parameter and whose kinetic expectation is evaluated in momentum space in terms of the modified Bessel function K_1.

Extended reading notes

Core claim

Starting from the N-body Hamiltonian H = sum_i $\sqrt$($p_i^{2}$+$m^{2}$) - sum_{i<j} v/r_ij, the paper derives an improved lower bound on the ground-state energy E. Using Jacobi relative coordinates, the boson-symmetric ground state is shown to satisfy the N-representability relations <ρ_i·ρ_j>=δ_ij<$ρ_2^{2}$> and <π_i·π_j>=δ_ij<$π_2^{2}$>, which imply in particular that the squared kinetic term separates as <(p+δP)^2>=2λ<$p^{2}$> with λ=(N-1)/N. Squaring the eigenvalue equation for the reduced Hamiltonian and minimizing over states satisfying this identity converts the problem into a one-body Kratzer operator, whose known lower bound yields E ≥ Nm $\sqrt$((1+$\sqrt$(1-$γv^{2}$))/2) whenever $γv^{2}$<1 with γ=N(N-1)/2. A Gaussian trial function in relative coordinates, with its width optimized, supplies the upper bound E ≤ Nm $\sqrt$(2/π) min_{μ>0}[2g($μ^{2}$/4)-$\sqrt$(γ)v/μ]. In the scaling v=c/N the lower bound is valid for c<$\sqrt$(2) and the upper bound for c<2sqrt(2); at weak coupling both reduce to Nm(1-$λc^{2}$/16) ≤ E(c) ≤ Nm(1-$λc^{2}$/(6π)), matching earlier nonrelativistic results.

Load-bearing premise

The whole improved bound rests on the identity <(p+δP)^2>=2λ<$p^{2}$>, which relies on the ground state of identical bosons having vanishing momentum cross-correlations between different Jacobi directions; if that symmetry-based relation fails, the derivation of the lower bound collapses.

Editorial extensions

If this is right

  • For v=c/N, after dividing by Nm the bounds become functions of c and λ alone, so they remain valid for arbitrarily large N as long as c<sqrt(2).
  • At weak coupling the bounds reduce to Nm(1-λc^2/16) ≤ E(c) ≤ Nm(1-λc^2/(6π)), reproducing the known nonrelativistic parabolic energy estimates.
  • The improved lower bound lies strictly above the earlier N/2 bound (5), so the two-sided bracket around the true energy is narrowed.
  • The overlap of the lower and upper bounds is restricted to γv^2<1, giving a concrete range of couplings in which the ground-state energy is controlled from both sides.

Reading between the lines

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

  • The same squaring-plus-boson-symmetry scheme should extend to other attractive pair potentials with a one-body Kratzer-type bound, such as Yukawa or truncated Coulomb potentials; the paper only treats V=-v/r.
  • Because the bound is exact for N=2, since it reduces to the one-body relativistic Coulomb bound, the N-dependence of the gap between bound and true energy is a natural target for numerical tests on N=3,4 systems.
  • The condition γv^2<1 offers a quantitative estimate for the maximum particle number N at fixed v before the semirelativistic Hamiltonian loses its lower bound, which is relevant to models of boson-star collapse.
  • An external trapping potential would break translation invariance and complicate the Jacobi reduction; extending the bound to trapped bosons would require a new ingredient.
Original abstract

Analytic energy bounds for N-boson systems governed by semirelativistic Hamiltonians of the form H=\sum_{i=1}^N(p_i^2 + m^2)^{1/2} - sum_{1=i<j}^N v/r_{ij}, with v>0, are derived by use of Jacobi relative coordinates. For gravity v=c/N, these bounds are substantially tighter than earlier bounds and they are shown to coincide with known results in the nonrelativistic limit.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 2,410 input tokens.

What the abstract itself states, as probabilities: claims first 0.05 · states limitation 0.12 · human participants 0.01 · controlled comparison 0.13 · randomized assignment 0.01 · states quantitative result 0.53 · compares with prior results 0.96 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

Analytic energy bounds for N-boson systems governed by semirelativistic Hamiltonians of the form H=\sum_{i=1}^N(p_i^2 + m^2)^{1/2} - sum_{1=i<j}^N v/r_{ij}, with v>0, are derived by use of Jacobi relative coordinates. For gravity v=c/N, these bounds are substantially tighter than earlier bounds and they are shown to coincide with known results in the nonrelativistic limit.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Analytic energy bounds for N-boson systems governed by semirelativistic Hamiltonians of the form H=\sum_{i=1}^N(p_i^2 + m^2)^{1/2} - sum_{1=i<j}^N v/r_{ij}, with v>0, are derived by use of Jacobi relative coordinates. contribution 0.84 0.93 0.40 asserted 1.00
s1 For gravity v=c/N, these bounds are substantially tighter than earlier bounds and they are shown to coincide with known results in the nonrelativistic limit. result 0.88 0.86 0.48 asserted 1.00

What kind of paper the classifier reports this is: theoretical result 1.00 · evidence from proof or derivation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Energy bounds E(v) (in dimensionless units) for the semirelativistic N-body problem with N = 5 and mass m = 1, as functions of the dimensionless coupling parameter v. The figure shows the simple lower bound (SL), the improved lower bound (L), and a scale-optimized Gaussian upper bound (U). where φ(r) = (α/π) 3/4 exp(−1 2 αr2 ). The kinetic-energy expectation value may be expressed in terms of the modified Bessel function of the second kind K1(x) [26]. Evaluating the integral, for convenience, in momentum space, yields  φ, p λp2 + m2φ  = 2m µ r 2 π g  µ 2 4  ,
Figure 2. Upper and lower bounds to the ground-state energy E(c) (in dimensionless units) of the gravitational N-particle problem, N = 2, 3, . . . , 6, with mass m = 1, as functions of the dimensionless coupling parameter c. 7. Results and conclusion Any consideration of an arbitrarily large number, N, of self-gravitating bosons, that is, the inspection of the limit N → ∞, obliges us to force the gravitational coupling v to decrease for increasing N in some appropriate manner, in order to avoid the (“relativistic”) collapse E → −∞ for N → ∞. Keeping v fixed, inevitably implies collapse [27].…

Figures reproduced from arXiv:math-ph/0602059 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.