Pith. Sign in

REVIEW

Kinetic energy operator approach to the quantum three-body problem with Coulomb interactions

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Xuguang Chi, Wuyi Hsiang, Ping Sheng

Computational Physics · Submitted 21 Feb 2006 · 96 references

Pith's verdict: A clever non-variational reduction with impressive-looking eigenvalues, but the printed l=1 odd-parity equations carry a sign error in the B1 term, so the central derivation is not reproducible 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 →

· 2 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDphysics/0602139 (version 1)
Pith IDC2QM34SX, Pith's permanent record of this paper and its review
CitationPing Sheng, Wuyi Hsiang, Xuguang Chi. "Kinetic energy operator approach to the quantum three-body problem with Coulomb interactions." arXiv:physics/0602139 (2006). https://pith.science/paper/C2QM34SX
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

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

A systematic kinetic-energy-operator reduction, using Jacobi-spherical coordinates, solves the three-body Coulomb Schrödinger problem as a one-dimensional algebraic eigenvalue problem.

The reading

The paper presents a way to solve the Schrödinger equation for three particles that interact through electric forces, such as two electrons and a nucleus. Instead of guessing a trial wavefunction and adjusting parameters, as variational methods do, the authors build the solution from the symmetries of the kinetic energy operator, the part of the equation that describes how the particles move. They first remove the overall motion of the center of mass, then project out the rotational motion, and then introduce special coordinates called Jacobi-spherical coordinates. In these coordinates, the remaining equation separates into a radial part and an angular part, and the angular part can be solved exactly with Jacobi polynomials. The result is a set of coupled one-dimensional equations in the radial variable, which is then expanded in Laguerre polynomials and solved as a large matrix eigenvalue problem.

The numerical scheme uses a special integration rule to compute the Coulomb interaction matrix elements accurately, and a truncation path with extrapolation to estimate the infinite-matrix limit. The authors apply the method to four systems: the negative hydrogen ion, helium and helium-like ions, the hydrogen molecular ion, and the positronium negative ion. For the ground states of these systems, their energies agree with very accurate variational calculations to many digits. For several excited states, they report energies and wavefunction properties that they say are more accurate or even new.

The main advantage claimed is that the method is non-variational, so it does not rely on choosing a good trial function, and it yields wavefunctions at the same time as eigenvalues.

Extended reading notes

Core claim

The load-bearing assertion is that the four-step reduction outlined in Section 3 reduces the nine-dimensional Coulomb three-body problem to a one-dimensional coupled system of ODEs that can be solved accurately as an infinite-dimensional algebraic eigenvalue problem, achieving up to thirteen-significant-figure accuracy for ground states and new results for several excited states.

Load-bearing premise

The accuracy of the final extrapolated energies depends on the assumed power-law convergence form E(x)=E_inf-c0 x^a(1+c1 x+...) with a state-dependent exponent a estimated from the data; if that functional form is wrong, the extrapolated eigenvalues could be biased. This assumption enters in Section 4.2.2 and is not independently justified.

Original abstract

We present a non-variational, kinetic energy operator approach to the solution of quantum three-body problem with Coulomb interactions, based on the utilization of symmetries intrinsic to the kinetic energy operator, i.e., the three-body Laplacian operator with the respective masses. Through a four-step reduction process, the nine dimensional problem is reduced to a one dimensional coupled system of ordinary differential equations, amenable to accurate numerical solution as an infinite-dimensional algebraic eigenvalue problem. A key observation in this reduction process is that in the functional subspace of the kinetic energy operator where all the rotational degrees of freedom have been projected out, there is an intrinsic symmetry which can be made explicit through the introduction of Jacobi-spherical coordinates. A numerical scheme is presented whereby the Coulomb matrix elements are calculated to a high degree of accuracy with minimal effort, and the truncation of the linear equations is carried out through a systematic procedure

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, 3,805 input tokens.

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

The abstract, sentence by sentence

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

We present a non-variational, kinetic energy operator approach to the solution of quantum three-body problem with Coulomb interactions, based on the utilization of symmetries intrinsic to the kinetic energy operator, i.e., the three-body Laplacian operator with the respective masses. Through a four-step reduction process, the nine dimensional problem is reduced to a one dimensional coupled system of ordinary differential equations, amenable to accurate numerical solution as an infinite-dimensional algebraic eigenvalue problem. A key observation in this reduction process is that in the functional subspace of the kinetic energy operator where all the rotational degrees of freedom have been projected out, there is an intrinsic symmetry which can be made explicit through the introduction of Jacobi-spherical coordinates. A numerical scheme is presented whereby the Coulomb matrix elements are calculated to a high degree of accuracy with minimal effort, and the truncation of the linear equations is carried out through a systematic procedure

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 We present a non-variational, kinetic energy operator approach to the solution of quantum three-body problem with Coulomb interactions, based on the utilization of symmetries intrinsic to the kinetic energy operator, i.e., the three-body Laplacian operator with the respective masses. contribution 0.95 0.94 0.25 asserted 1.00
s1 Through a four-step reduction process, the nine dimensional problem is reduced to a one dimensional coupled system of ordinary differential equations, amenable to accurate numerical solution as an infinite-dimensional algebraic eigenvalue problem. method 0.98 0.90 0.20 asserted 1.00
s2 A key observation in this reduction process is that in the functional subspace of the kinetic energy operator where all the rotational degrees of freedom have been projected out, there is an intrinsic symmetry which can be made explicit through the introduction of Jacobi-spherical coordinates. method 0.82 0.83 0.10 asserted 1.00
s3 A numerical scheme is presented whereby the Coulomb matrix elements are calculated to a high degree of accuracy with minimal effort, and the truncation of the linear equations is carried out through a systematic procedure method 0.98 0.90 0.09 asserted 1.00

What kind of paper the classifier reports this is: method or tool 1.00 · evidence from simulation or computation 0.55.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
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. The triangle formed by three particles. 5. Results on three-body systems In this section we present results on some three-body systems, and discuss the properties of their relevant wave functions. In order to display the data in physical space, we first give a description of the relevant angles, as well as the procedure by which the wave functions are re-constituted from numerical data. 5.1. Euler angles and wave functions 5.1.1. Euler angles and the distribution function The three-body wave function depends upon the shape of the triangle formed by three particles and the three Euler
Figure 2. Distributions of r1, r2 ( or r3) and cos θ of the ground state of ∞H −. The solid line is the distribution of r2 or r3 (same), and the dashed line is the distribution of r1. Inset is the distribution of cos θ. main source of observed opacity in the solar atmosphere. The absorption coefficient of H − has been studied extensively and used in the theory of the solar atmosphere. In fact, discrepancies between early calculations and observational evidence on the sun’s radiation have pointed out the inaccuracies of the calculated H− wave functions available then. Calculation of the wave f
Figure 3. Radial density probability distribution P(r) for ∞H −( 3S), H(1s) and H(2s). bound in a short-range attractive potential well. In modern variational calculations, many r j 1 terms are necessary in the trial wave function in order to fully account for the electron-electron correlation. While many-parameter variational calculations can give great accuracy for the ground state energy, the best experimental values come from a high resolution (0.03cm−1 ) laboratory photodetachment laser experiment. The binding energy has been determined to be 6082.99 ± 0.15cm−1 for H− and 6086.2 ±…
Figure 4. Distributions of r1, r2 (or r3) and cos θ of the 3P e state of ∞H −. The solid line is the distribution of r2 or r3 (same), and the dashed line is the distribution of r1. Inset is the distribution of cos θ. figure 3, the electron distribution of ∞H − is seen to be very different from the neutral 1s or 2s state. Thus the electron-electron correlation plays an important role in H−, and independent-electron model can not give an accurate picture. 5.2.2. The 3P e state The even-parity 3P e state of ∞H − and H− is quasi-stable, and its energy is just below the n=2 threshold (-0.12
Figure 5. Radial density probability distribution P(r) for ∞H −( 3P), H(2p), H(3p) and H(4p). We compare our results with other calculations in table 2. The convergence is just as good in this case as for the ground state, and the leading term of all the sequences is estimated to be x 5 . Our energy eigenvalue is 0.125 355 451 242 a.u. for ∞H−, better than most of previous calculations. As a negative ion, the extra electron is significantly affected by the other electron
Figure 6. Distributions of r1, r2 (or r3) of the five low-lying 3P e states of ∞He (ground state plus four excited states). The solid line is the distribution of r2 or r3 (same), and the dashed line is the distribution of r1. For comparison, we also plot the radial density distribution of the neutral He(2p), H(2p) and H(3p) states in (5) (normalized to 1/2)
Figure 7. The coordinate system of H+ 2 . of r1, r2, r3 and cos θ are shown in figure 6 for ∞He. Some properties of these states are given in tables 7 and 8. It should be mentioned that only in the case of E2 of the 1P e state is our result higher in energy than the variational approach, but has more significant figures. We think that in this case our result is correct, because our excited states’ eigen-energies and their related wavefunctions were all obtained simultaneously with the ground state properties. Hence the accuracy of our excited states’ properties are insured by the demonstrate
Figure 8. Distributions of R for the five low-lying S states of H+ 2 . Similarity to oscillator is noted
Figure 9. The function |Ψ(X, Y, R = R0)| for the ground state of H+ 2 , where R0 = 0.5, 1.0. The electron density is noted to peak directly above the positions of the protons, with a sharp cusp. is used, and the wave function has the form: Ψ = Ψ(X, Y, R). (152) Here R is the distance between the two protons, and X, Y denote the coordinate of the electron relative to the two protons
Figure 10. Same as in figure 9, but R0 = 1.5, 2.0. stable against dissociation into a free electron and a positronium atom. The formation of Ps− has been discussed in the e + −He [85] and Ps−H [86] scattering calculations. The binding energy has been calculated variationally by several authors. Recent interests include the calculations of autoionizing doubly-excited states and the investigations of the possible existence of the so-called second bound state 3P e . These positronium negative ions have been observed in the laboratory by Mills [87, 88]. A large number of the doubly-excited states…
Figure 11. Distributions of r1, r2 (or r3) and cos θ of the ground state of Ps−. The solid line is the distribution of r2 or r3 (same), and the dashed line is the distribution of r1. Inset is the distribution of cos θ. the present case, the mass of the positive charge is the same as that of the electron, so the positron is envisioned to be much more mobile, whereas in the case of H− the proton is more localized. This comparison is perhaps helpful in answering the question [90]: Is the positronium system a molecule or an atom? 6. Further developments The kinetic energy operator approach is…

Figures reproduced from arXiv:physics/0602139 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.