Combescot (Phys. Rev. X 7, 041035, 2017) reformulated the exact three-dimensional three-body Coulomb bound-state problem as a small set of coupled integral equations in which the interaction enters only through the known two-body Coulomb T-matrix, with all internal frequency integrations eliminated so that only on-the-shell quantities appear; the method reproduced the helium ground-state energy and was subsequently applied by the same author to compute the trion binding energy over the full mass-ratio range (Phys. Rev. B 100, 245201, 2019). The paper explicitly left open whether this strategy can be generalized 'along the same lines' to four-particle Coulomb bound states, a step it described as extremely useful, notably for treating exciton-exciton interactions in semiconductors. This four-body extension has not been derived: four-body Coulomb complexes are still treated by variational, quantum Monte Carlo, hyperspherical, or model-potential momentum-space methods rather than by an exact scheme built on the analytic two-body Coulomb T-matrix. The open problem is to derive the four-body analog of this integral-equation scheme — a closed equation system in which the pair interactions enter only through the two-body Coulomb T-matrices, with the analogous on-shell frequency elimination — and to validate it numerically on a benchmark four-body Coulomb bound state such as the three-dimensional biexciton.

Background

A number of physically central systems are few charged particles bound by pairwise Coulomb (inverse-distance) forces: the helium atom (a nucleus plus two electrons), the negative hydrogen ion, the trion in a semiconductor, and the biexciton. The trion is a 'charged exciton' — two electrons and one hole (or two holes and one electron), where a hole is the quasiparticle left behind when an electron is removed from a full valence band; it is usually modeled in the effective-mass approximation, in which band electrons and holes behave as free particles with effective masses interacting via a dielectrically screened Coulomb potential. The biexciton is the corresponding four-particle complex of two electrons and two holes; an exciton is one electron-hole pair. Few-body Coulomb problems have no closed-form solutions, and high accuracy has historically come from coordinate-space variational calculations, correlated-Gaussian stochastic-variational expansions, hyperspherical-harmonics expansions, and quantum Monte Carlo methods — all of which work directly with the wavefunction and the explicit Coulomb potential.

A different, scattering-based route exploits the fact that the two-body problem is exactly solvable. The two-body T-matrix, defined as the solution of the Lippmann-Schwinger equation T2=V+VgT2T_2=V+VgT_2 (with VV the potential and gg the free two-particle propagator), contains the complete two-body scattering and bound-state information. For short-range interactions, the scattering of a third particle off a bound pair obeys the Skorniakov-Ter-Martirosian equation — a low-dimensional integral equation in which the potential never appears explicitly, only through the two-body propagator; the systematic three-body generalizations are Faddeev's equations, which decompose the three-body amplitude into components associated with each interacting pair.

Combescot's 2017 paper carried this program over to the Coulomb interaction, which is long-ranged and therefore a notorious source of difficulties in momentum-space integral equations. Working diagrammatically, he derived the exact three-body integral equations for a general interaction in which the potential enters only through the two-body T-matrix, and then eliminated all internal frequency (energy) variables by contour integration, so that only 'on-shell' quantities — evaluated at the physical energy-momentum relation — remain. For the Coulomb case the required two-body T-matrix has a closed analytic form that follows from Schwinger's 1964 Coulomb Green's function. The result for the general three-body Coulomb bound-state problem is two coupled three-dimensional linear integral equations, whose eigenvalue condition yields the bound-state energies; for the helium ground state (infinitely heavy nucleus) the equations reduce to two-dimensional ones, and the calculation reproduces the known result 2K2=0.725932K^{2}=0.72593 with precision improving from 10−210^{-2} to 10−410^{-4} under mesh refinement, with the wavefunction obtained from the solution. In a 2019 follow-up, the same author applied the method to the trion in a three-dimensional semiconductor with bare Coulomb interaction in the effective-mass approximation, covering the complete electron-to-hole mass-ratio range in agreement with existing variational calculations.

What the 2017 paper explicitly did not do, and posed as an open path, is the four-body extension: deriving the analogous exact integral equations for four-particle Coulomb bound states, for which it expected the scheme to be 'extremely useful for many problems, in particular to treat appropriately exciton-exciton interaction in semiconductors which is of importance for the Bose-Einstein condensation of excitons'. The existing four-body Coulomb literature does not contain such a scheme: biexcitons and other four-body Coulomb complexes are computed by the variational and Monte Carlo methods mentioned above; momentum-space four-body Coulomb scattering has been treated only within static-exchange or model-potential approximations; and the Faddeev-type four-body integral equations (the Yakubovsky and Alt–Grassberger–Sandhas (AGS) equations) that are standard in nuclear few-body physics are formulated for short-range potentials, not for the Coulomb interaction entering through its analytic two-body Coulomb T-matrix with the frequency-eliminated, on-shell, dimension-reduced structure that characterizes Combescot's three-body scheme.

Problem Statement

The problem, as left open by Combescot (2017), is how the exact two-body-T-matrix integral-equation strategy of the three-body Coulomb problem can be generalized to four-particle Coulomb bound states along the same lines. A complete answer should:

  1. Derive the closed system of integral equations for the four-body Coulomb bound-state problem within this scattering/diagrammatic framework, in which all pair interactions (electron-electron, electron-hole, and hole-hole, or their atomic analogues with the appropriate reduced masses) enter only through the corresponding known two-body Coulomb T-matrices, with the correct treatment of particle-exchange symmetries for identical fermions.

  2. Carry out the analog of the three-body frequency elimination — the contour-integration step that removes all internal energy variables — so that only on-the-shell quantities appear, and characterize the resulting system: how many coupled equations, on how many independent momentum variables, the kernels depend.

  3. Demonstrate the scheme numerically by computing a benchmark four-body Coulomb bound-state energy — for example the ground-state binding energy of the three-dimensional biexciton (two electrons and two holes) as a function of the electron-to-hole mass ratio, or another four-body Coulomb system with a well-established reference value — in agreement with the known high-precision results for the chosen benchmark to within a stated tolerance, and with controlled numerical convergence.

These three components jointly serve one objective: an exact, compact, and verified four-body Coulomb scheme of the same type as the three-body one. The derivation alone, without the numerical check, would not establish that the generalization is workable 'along the same lines'; a four-body numerical result obtained by other means (variational, Monte Carlo) would not answer the methodological question.

References

  1. R. Combescot, Three-Body Coulomb Problem, Phys. Rev. X 7, 041035 (2017). DOI: 10.1103/PhysRevX.7.041035; arXiv:1705.06177, https://arxiv.org/abs/1705.06177
  2. R. Combescot, Trion ground-state energy: Simple results, Phys. Rev. B 100, 245201 (2019). DOI: 10.1103/PhysRevB.100.245201; arXiv:1909.13815, https://arxiv.org/abs/1909.13815
  3. Roman Ya. Kezerashvili, Few-Body Systems in Condensed Matter Physics, Few-Body Systems 60, 52 (2019). DOI: 10.1007/s00601-019-1520-1; arXiv:1902.07843, https://arxiv.org/abs/1902.07843
  4. Roman Ya. Kezerashvili, Bound Trions in Two-Dimensional Monolayers: A Review (2026). arXiv:2603.08346, https://arxiv.org/abs/2603.08346
  5. Julian Schwinger, Coulomb Green's Function, Journal of Mathematical Physics 5, 1606 (1964). DOI: 10.1063/1.1931195, https://doi.org/10.1063/1.1931195
  6. Hasi Ray, The collision between two ortho-Positronium (Ps) atoms: a four-body Coulomb problem (2014). arXiv:1409.3371, https://arxiv.org/abs/1409.3371
  7. K. Mohseni, M. R. Hadizadeh, T. Frederico, D. R. da Costa, A. J. Chaves, Trion clustering structure and binding energy in 2D semiconductor materials: Faddeev equations approach, Phys. Rev. B 107, 165427 (2023). DOI: 10.1103/PhysRevB.107.165427; arXiv:2212.09600, https://arxiv.org/abs/2212.09600
  8. Yasufumi Nakano, Meera M. Parish, Jesper Levinsen, Coupled-channel approach to scattering of hybrid excitons (2026). arXiv:2607.19634, https://arxiv.org/abs/2607.19634