REVIEW 3 major objections 3 minor 2 cited by
FewBodyToolkit.jl claims to be a general, open-source Julia solver for quantum two- and three-body bound states and resonances.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
FewBodyToolkit.jl is a documented, open-source Julia solver for two- and three-body quantum bound and resonant states that matches reference results for Coulomb, low-dimensional universal, nuclear resonance, and positronium-ion benchmarks.
T0 review reviewed 2026-08-04 challenge →
load-bearing objection A well-executed software paper for a useful Julia few-body toolkit; the main caveat is that parameter sensitivity can silently produce wrong states, and that deserves a clearer look. the 3 major comments →
FewBodyToolkit.jl: a Julia package for solving quantum few-body problems
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
Core claim
The paper introduces FewBodyToolkit.jl as an open-source Julia implementation of the Gaussian expansion method for general two- and three-body quantum systems. The wavefunction is expanded in Gaussian basis functions whose widths form a geometric progression, and an optional complex range adds oscillatory structure to reach highly excited states and resonances. The Schrödinger equation becomes a generalized eigenvalue problem, H c = E S c, with mostly analytic matrix elements; resonances are obtained by complex scaling. The package is organized into three solver modules—a two-body solver, a one-dimensional three-body solver, and a three-dimensional three-body solver—and is validated against
What carries the argument
Gaussian expansion method: expand each relative-motion wavefunction as a superposition of Gaussians with widths in a geometric progression, using products of two such Gaussians in Jacobi coordinates for three-body Faddeev components. This converts the Schrödinger equation into a dense generalized eigenvalue problem with analytic overlap, kinetic, and (for several potentials) interaction matrix elements; complex-ranged Gaussians and complex scaling extend the same machinery to resonances.
Load-bearing premise
The solver's accuracy rests on the user supplying converging numerical parameters (the number and widths of the Gaussian basis functions), and the paper itself shows that unoptimized defaults can miss bound states and generate spurious eigenvalues.
What would settle it
Recompute the two-body Coulomb example with the paper's non-optimized parameters (nmax=10, r1=0.1, rnmax=30): the paper predicts exactly four bound states with energies near -0.499876, -0.124543, -0.054437, -0.028644 and two spurious positive eigenvalues near +0.007342 and +0.603176; any materially different pattern would contradict the reported behavior.
If this is right
- Researchers can compute bound and resonant states for two- and three-body systems with arbitrary central pair potentials, including potentials given only as numerical tables.
- Complex-ranged Gaussian bases make highly excited states and resonances tractable with modest basis sizes, as shown by the Coulomb spectrum through n=40.
- Automatic Faddeev decomposition and identical-particle symmetrization in the three-body modules remove a major manual step in setting up three-body calculations.
- The package offers a reproducible reference for teaching, benchmarking, and testing new few-body methods.
- If adopted, it could become a common baseline that different few-body approaches compare against, similar to standard packages in other numerical disciplines.
Where Pith is reading between the lines
- Editorial inference: the accessibility claim depends on users selecting good Gaussian width parameters; the paper's own Table 2 shows that defaults can miss bound states and produce spurious positive-energy eigenvalues, so automatic basis optimization would materially strengthen the package.
- Editorial inference: because the Gaussian expansion is variational for bound states, the same solver output could be extended to give rigorous upper-bound certificates and convergence estimates for each eigenvalue.
- Editorial inference: the interpolation scheme used for three-body interaction integrals could be adapted to nonlocal or momentum-dependent potentials, a natural extension the paper lists as future work.
- Editorial inference: the one-dimensional three-body module's reproduction of universal energy ratios suggests the package could serve as a cheap testbed for universal few-body predictions across mass ratios and interaction ranges.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces FewBodyToolkit.jl, an open-source Julia package implementing the Gaussian expansion method for two- and three-body quantum few-body problems. It describes three solver modules (GEM2B, GEM3B1D, ISGL), covering bound states and resonances via complex scaling, several spatial dimensions, symmetries, and general pairwise interaction shapes. The package is validated against exact Coulomb energies (up to n = 40, differences of order 10^-6 or smaller), resonance positions and widths for a nuclear-similar potential (agreement to ~0.5% or better), universal 1D three-body energy ratios (four digits), and Ps− bound-state properties (binding energy within 0.073%, radii within 0.45%). The paper claims that the package is an accessible, documented general few-body solver suitable for research and teaching.
Significance. The package addresses a genuine gap—there is no widely available, documented general few-body solver with this feature combination. The manuscript's positive evidence is strong and appropriately external: the Coulomb comparison uses the independent Antique.jl code, resonances are checked against Lazauskas [46], and Ps− against Frolov [47]. The inclusion of runnable examples, documentation, and the honest reporting of non-converged default parameters in Tables 2 and 3 are commendable. If the caveats discussed below are addressed, the paper would be a useful resource for benchmark studies, method development, and teaching. The central method itself is standard and the derivation is sound.
major comments (3)
- [§5.2, Tables 2–3] The accessibility claim is undercut by the behavior shown in Tables 2 and 3. With the non-optimized parameters in the worked example (nmax=10, r1=0.1, rnmax=30), the Coulomb solver returns only 4 of the first 6 bound states and two spurious positive-energy eigenvalues (+0.0073, +0.603) without any warning or diagnostic. The 2-body optimizer GEM_Optim_2B restores state 6 but degrades the ground state from -0.499876 to -0.489956 (~2% error). Section 4.2 lists input validation but no spectral-quality/convergence check, and no optimizer or automatic diagnostic is described for GEM3B1D or ISGL. The paper should either provide automated guidance/diagnostics for choosing Gaussian ranges and detecting spurious states in all modules, or explicitly frame the package as requiring expert parameter tuning and narrow the 'accessible' claim accordingly.
- [§3.2, §7, and abstract] The abstract's 'arbitrary pair-interactions' is stronger than what is implemented. Section 3.2 restricts interactions to those satisfying f(r→0)∼r^{−α} with α<d+l+l′, and Section 7 concedes that broad strongly repulsive interactions remain challenging. None of the benchmarks probes such a boundary case (all use Coulomb, exponential-screened, Gaussian, or contact potentials). Please replace 'arbitrary' with a precise scope statement, or add a benchmark/example that demonstrates the advertised flexibility at the stated boundary.
- [Table 1 and §5] Table 1 lists 'Resonances via CSM: Yes' for all three modules, and the abstract promises calculation of resonant states. Yet Section 5 demonstrates complex scaling only for GEM2B (Table 5); no three-body resonance calculation is shown for GEM3B1D or ISGL. Because three-body resonances are a core motivation of the package, please add at least one three-body resonance benchmark/example, or explicitly state in Table 1 and the text that three-body CSM is implemented but not yet validated.
minor comments (3)
- [§5.3, Table 6] The 'Reference' column in Table 6 is from Ref. [32], a previous publication by the same author. The agreement is valuable but is not a fully independent validation; please state this explicitly or add an independent reference.
- [General (codebase publication)] The manuscript does not specify a version number or archived DOI for the released package. For reproducibility, please cite a tagged release (e.g., Zenodo) in addition to the GitHub repository.
- [§5.2 code snippet] The variables simax and redind are defined in the example script but not used; removing them would make the pedagogical code clearer.
Circularity Check
No significant circularity: all central validations are against external benchmarks; the one self-benchmark (Ref. [32]) is independent substantive support.
full rationale
The paper is a software/methods paper, not a derivation of new physical predictions. Its load-bearing validations are external: the Coulomb spectrum is compared with exact values from Antique.jl (Tables 2-4), the complex-scaling resonances with results of Lazauskas [46] (Table 5), and the Ps- binding energy and radii with results of Frolov [47] (Table 7). These references do not depend on the package's inputs or fitted parameters. The only benchmark against work co-authored by the present author is Table 6, which uses universal three-body energy ratios from Ref. [32]; those ratios are published universal quantities, not fitted numbers, and the package reproduces them from independently constructed Gaussian and contact interactions. No equation in the paper defines a predicted quantity in terms of the fitted parameters, and the optimization example in Tables 2-3 is explicitly presented as a trade-off rather than as a prediction. The 'arbitrary pair-interactions' claim is explicitly scoped by the boundedness condition in Sec. 3.2 and by admitted challenges in Sec. 7, which is a limitation, not a circular step. The Gaussian expansion method itself is the standard method of Ref. [19] and is not derived from the package's outputs. Accordingly, no circular step is present.
Axiom & Free-Parameter Ledger
free parameters (3)
- Two-body binding energy scale (target_e2 = -1e-3) =
-1e-3
- Optimized Gaussian ranges (r1, rnmax) =
r1=0.871259, rnmax=45.664907
- Complex-scaling angle theta_csm =
40.0 degrees
axioms (4)
- domain assumption Center-of-mass motion separates from internal dynamics (no external forces, or harmonic confinement only)
- domain assumption Gaussian basis expansion converges for the target states with tractable basis sizes
- domain assumption Complex scaling (CSM) exposes resonance poles without distorting bound-state results
- domain assumption Integrability/boundedness condition f(r→0) ~ r^-alpha with alpha < d + l + l'
Cite this review
Pith. "Pith review of FewBodyToolkit.jl: a Julia package for solving quantum few-body problems." pith.science (2026). https://pith.science/paper/6IUDKRSF
@misc{pith2026251004447,
author = {Pith},
title = {Pith review of: FewBodyToolkit.jl: a Julia package for solving quantum few-body problems},
year = {2026},
howpublished = {\url{https://pith.science/paper/6IUDKRSF}},
note = {Machine review of arXiv:2510.04447}
}
read the original abstract
Few-body physics explores quantum systems of a small number of particles, bridging the gap between single-particle and many-body regimes. To provide an accessible tool for such studies, we present FewBodyToolkit.jl, a Julia package for quantum few-body simulations. The package supports general two- and three-body systems in various spatial dimensions with arbitrary pair-interactions, and allows to calculate bound and resonant states. The implementation is based on the well-established Gaussian expansion method and we illustrate the package's capabilities through benchmarks and research examples. The package comes with documentation and examples, making it useful for research, teaching, benchmarking, and method development.
Figures
Forward citations
Cited by 2 Pith papers
-
Stabilization of three-body resonances to bound states in a continuum
A two-channel model demonstrates stabilization of three-body resonances into bound states in the continuum via parameter tuning, shown in 1D mass-imbalanced and 3D Efimov systems with magnetic field control.
-
From three-body resonances to bound states in a continuum: pole trajectories
In a 1D three-body model, each of three parameter scans—interaction strength, range, and mass ratio—drives a resonance pole onto the real axis, forming at least one bound state in the continuum.
Reference graph
Works this paper leans on
-
[1]
Newton, Philosophiae naturalis principia mathematica , Jussu Societatis Regiae ac Typis Josephi Streater
I. Newton, Philosophiae naturalis principia mathematica , Jussu Societatis Regiae ac Typis Josephi Streater. Prostat apud plures bibliopolas, London (1687)
-
[2]
Euler, De motu rectilineo trium corporum se mutuo attrahentium, Novi Commentarii academiae scientiarum Petropolitanae pp
L. Euler, De motu rectilineo trium corporum se mutuo attrahentium, Novi Commentarii academiae scientiarum Petropolitanae pp. 144--151 (1767)
-
[3]
J. L. Lagrange, Essai sur le Probl \`e me des Trois Corps , Prix de l'Acad \'e mie Royale des Sciences de Paris IX(6), 229 (1772)
-
[4]
Burrau, Berechnung des Energiewertes des Wasserstoffmolekel-Ions (H2+) im Normalzustand , Naturwissenschaften 15(1), 16 (1927), doi:10.1007/BF01504875
-
[5]
B. D. Esry, C. H. Greene and J. P. Burke, Recombination of Three Atoms in the Ultracold Limit , Phys. Rev. Lett. 83(9), 1751 (1999), doi:10.1103/PhysRevLett.83.1751
-
[6]
Efimov, Energy levels arising from resonant two-body forces in a three-body system, Phys
V. Efimov, Energy levels arising from resonant two-body forces in a three-body system, Phys. Lett. B 33(8), 563 (1970), doi:10.1016/0370-2693(70)90349-7
-
[7]
P. Naidon and S. Endo, Efimov physics: a review, Rep. Prog. Phys. 80(5), 056001 (2017), doi:10.1088/1361-6633/aa50e8
-
[8]
J. R. Johansson, P. D. Nation and F. Nori, Qutip: An open-source python framework for the dynamics of open quantum systems, Comp. Phys. Comm. 183(8), 1760 (2012), doi:10.1016/j.cpc.2012.02.021
-
[9]
J. R. Johansson, P. D. Nation and F. Nori, Qutip 2: A python framework for the dynamics of open quantum systems, Comp. Phys. Comm. 184(4), 1234 (2013), doi:10.1016/j.cpc.2012.10.019
-
[10]
N. Lambert, E. Gigu \`e re, P. Menczel, B. Li, P. Hopf, G. Su \'a rez, M. Gali, J. Lishman, R. Gadhvi, R. Agarwal, A. Galicia, N. Shammah et al., QuTiP 5: The Quantum Toolbox in Python , doi:10.48550/arXiv.2412.04705 (2025), 2412.04705
-
[11]
J. M. Hutson and C. R. Le Sueur, MOLSCAT : A program for non-reactive quantum scattering calculations on atomic and molecular collisions , Comp. Phys. Comm. 241, 9 (2019), doi:10.1016/j.cpc.2019.02.014
-
[12]
B. Schmidt and U. Lorenz, WavePacket : A Matlab package for numerical quantum dynamics. I : Closed quantum systems and discrete variable representations , Comp. Phys. Comm. 213, 223 (2017), doi:10.1016/j.cpc.2016.12.007
-
[13]
B. Schmidt and C. Hartmann, WavePacket : A Matlab package for numerical quantum dynamics. II : Open quantum systems, optimal control, and model reduction , Comp. Phys. Comm. 228, 229 (2018), doi:10.1016/j.cpc.2018.02.022
-
[14]
R. De la Fuente and M. Lamorena, QMsolve: A module for solving and visualizing the Schrödinger equation , doi:10.5281/zenodo.11181977
-
[15]
S. B. S. Miller, A. Ekstr\"om and K. Hebeler, Neutron-deuteron scattering cross sections with chiral nn interactions using wave-packet continuum discretization, Phys. Rev. C 106, 024001 (2022), doi:10.1103/PhysRevC.106.024001
-
[16]
FewBodyECG.jl , https://github.com/JuliaFewBody/FewBodyECG.jl
-
[17]
Twobody.jl , https://github.com/ohno/TwoBody.jl
-
[18]
JPublicThreeBodySolver , https://github.com/roudnev/JPublicThreeBodySolver
-
[19]
E. Hiyama, Y. Kino and M. Kamimura, Gaussian expansion method for few-body sys., Prog. Part. Nucl. Phys. 51(1), 223 (2003), doi:10.1016/S0146-6410(03)90015-9
-
[20]
T. Yoshida, E. Hiyama, A. Hosaka, M. Oka and K. Sadato, Spectrum of heavy baryons in the quark model, Phys. Rev. D 92(11), 114029 (2015), doi:10.1103/PhysRevD.92.114029
-
[21]
H.-X. Zhu, L. Meng, Y. Ma, N. Li, W. Chen and S.-L. Zhu, Constraining the DDD ^ * three-body bound state via the Z _ c (3900) pole , Phys. Rev. D 111(9), 094022 (2025), doi:10.1103/PhysRevD.111.094022
-
[22]
L.-Z. Wen, Y. Ma, L. Meng and S.-L. Zhu, NN , J / NN , _ c NN systems based on HAL QCD interactions , Phys. Rev. D 111(11), 114004 (2025), doi:10.1103/mvqk-n377
-
[23]
A. J. Arifi, L. Happ, S. Ohno and M. Oka, Structure of heavy mesons in the light-front quark model, Phys. Rev. D 110(1), 014020 (2024), doi:10.1103/PhysRevD.110.014020
-
[24]
E. Hiyama and M. Kamimura, Three-body structure of 6He = 4He +n+n using realistic n-n potentials , Nucl. Phys. A 588(1), c35 (1995), doi:10.1016/0375-9474(95)00096-J
-
[25]
N. Yamanaka and E. Hiyama, Enhancement of the CP -odd effect in the nuclear electric dipole moment of ^ 6 Li , Phys. Rev. C 91(5), 054005 (2015), doi:10.1103/PhysRevC.91.054005
-
[26]
C. H. Schmickler, H.-W. Hammer and E. Hiyama, Tetramer bound states in heteronuclear systems, Phys. Rev. A 95(5), 052710 (2017), doi:10.1103/PhysRevA.95.052710
-
[27]
L. Happ, P. Naidon and E. Hiyama, Mass Ratio Dependence of Three-Body Resonance Lifetimes in 1D and 3D , Few-Body Sys. 65(2), 38 (2024), doi:10.1007/s00601-024-01900-w
-
[28]
L. Happ and P. Naidon, Stabilization of three-body resonances to bound states in a continuum, doi:10.48550/arXiv.2503.02037 (2025), 2503.02037
-
[29]
F. Serwane, G. Z \"u rn, T. Lompe, T. B. Ottenstein, A. N. Wenz and S. Jochim, Deterministic Preparation of a Tunable Few-Fermion System , Science 332(6027), 336 (2011), doi:10.1126/science.1201351
-
[30]
S. I. Mistakidis, A. G. Volosniev, R. E. Barfknecht, T. Fogarty, Th . Busch, A. Foerster, P. Schmelcher and N. T. Zinner, Few-body Bose gases in low dimensions--- A laboratory for quantum dynamics , Phys. Rep. 1042, 1 (2023), doi:10.1016/j.physrep.2023.10.004
-
[31]
E. Garrido and A. S. Jensen, Confinement of N-Body Systems and Non-integer Dimensions , Few-Body Sys. 65(2), 35 (2024), doi:10.1007/s00601-024-01906-4
-
[32]
L. Happ, M. Zimmermann, S. I. Betelu, W. P. Schleich and M. A. Efremov, Universality in a one-dimensional three-body system, Phys. Rev. A 100(1), 012709 (2019), doi:10.1103/PhysRevA.100.012709
-
[33]
L. Happ and M. A. Efremov, Proof of universality in one-dimensional few-body sys. including anisotropic interactions, J. Phys. B: At. Mol. Opt. Phys. 54(21), 21LT01 (2021), doi:10.1088/1361-6455/ac3b3f
-
[34]
L. Happ, M. Zimmermann and M. A. Efremov, Universality of excited three-body bound states in one dimension, J. Phys. B: At. Mol. Opt. Phys. 55(1), 015301 (2022), doi:10.1088/1361-6455/ac3cc8
-
[35]
T. Schnurrenberger, L. Happ and M. A. Efremov, Borromean states in a one-dimensional three-body system, Phys. Rev. Res. 7(1), 013090 (2025), doi:10.1103/PhysRevResearch.7.013090
-
[36]
Nishida, Semisuper Efimov effect induced by resonant pair exchange in mixed dimensions , Phys
Y. Nishida, Semisuper Efimov effect induced by resonant pair exchange in mixed dimensions , Phys. Rev. A 111(5), L051306 (2025), doi:10.1103/PhysRevA.111.L051306
-
[37]
P. A. Belov, Linewidths and energy shifts of electron-impurity resonant states in quantum wells with infinite barriers, Phys. Rev. B 105(15), 155417 (2022), doi:10.1103/PhysRevB.105.155417
-
[38]
P. A. Belov, F. Morawetz, S. O. Kr \"u ger, N. Scheuler, P. Rommel, J. Main, H. Giessen and S. Scheel, Energy states of Rydberg excitons in finite crystals: From weak to strong confinement , Phys. Rev. B 109(23), 235404 (2024), doi:10.1103/PhysRevB.109.235404
-
[39]
L. K \"u hner, P. Rommel, J. Main, S. Scheel and P. A. Belov, Wave functions and oscillator strengths of Rydberg excitons in cuprous oxide quantum wells , Phys. Rev. B 112(7), 075409 (2025), doi:10.1103/9pct-k79b
-
[40]
A. Aslanidis, J. Main, P. Rommel, S. Scheel and P. A. Belov, Bound states in the continuum in cuprous oxide quantum wells, Phys. Rev. B 111(12), L121103 (2025), doi:10.1103/PhysRevB.111.L121103
-
[41]
L. Happ, M. A. Efremov, H. Nha and W. P. Schleich, Sufficient condition for a quantum state to be genuinely quantum non- Gaussian , New J. Phys 20(2), 023046 (2018), doi:10.1088/1367-2630/aaac25
-
[42]
M. Walschaers, Non- Gaussian Quantum States and Where to Find Them , PRX Quantum 2(3), 030204 (2021), doi:10.1103/PRXQuantum.2.030204
-
[43]
W. P. Schleich, I. Tk \'a c ov \'a and L. Happ, Insights Into Complex Functions , In R. Citro, M. Lewenstein, A. Rubio, W. P. Schleich, J. D. Wells and G. P. Zank, eds., Sketches of Physics : The Celebration Collection , pp. 127--159. Springer International Publishing, Cham, doi:10.1007/978-3-031-32469-7_5 (2023)
-
[44]
N. Moiseyev, Quantum theory of resonances: calculating energies, widths and cross-sections by complex scaling, Phys. Rep. 302(5), 212 (1998), doi:10.1016/S0370-1573(98)00002-7
-
[45]
Antique.jl , https://github.com/ohno/Antique.jl
-
[46]
Lazauskas, Numerical Aspects of Resonant States in Quantum Mechanics , Few-Body Sys
R. Lazauskas, Numerical Aspects of Resonant States in Quantum Mechanics , Few-Body Sys. 64(2), 24 (2023), doi:10.1007/s00601-023-01808-x
-
[47]
A. M. Frolov, Bound-state properties of the positronium negative ion Ps ^- , Phys. Rev. A 60(4), 2834 (1999), doi:10.1103/PhysRevA.60.2834
-
[48]
A. J. Arifi and K. Suzuki, Structure of heavy quarkonia under strong magnetic field, doi:10.48550/arXiv.2507.18894 (2025), 2507.18894
-
[49]
N. Yamanaka, Light nuclei under magnetic field and the lithium problem, doi:10.48550/arXiv.2509.03684 (2025), 2509.03684
-
[50]
H. A. Bethe, Zur Theorie der Metalle. i. Eigenwerte und Eigenfunktionen der linearen Atomkette , Zeit. f \"u r Phys. 71 , 205 (1931), 10.1007\
1931
-
[51]
Ginsparg, It was twenty years ago today
P. Ginsparg, It was twenty years ago today... , http://arxiv.org/abs/1108.2700
-
[52]
, " * write output.state after.block =
ENTRY address archive author booktitle chapter doi edition editor eid eprint howpublished institution isbn journal key month note number organization pages publisher school series title type url volume year label INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence := #2 'af...
-
[53]
write newline
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize " " * FUNCT...
This paper was first reviewed by deepseek-v4-flash on August 4, 2026.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.