REVIEW 2 major objections 7 minor 49 references
The Faddeev-Yakubovsky symphony
T0 review · 2 major / 7 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The Faddeev-Yakubovsky decomposition turns the N-body Schrödinger equation into coupled integro-differential systems, now solved through five particles.
desk verdict A clear, useful review of known FY formalism with no new results; worth refereeing after proofreading and a more honest framing of where the rigor lives. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
At the centre is the FY component: a part of the total wavefunction tied to one complete partition chain, i.e. one tree diagram showing how the N-particle cluster is broken, step by step, until only one interacting pair remains. The recursion starts with the Faddeev partition $\Psi = \sum_{i<j} \Phi_{ij}$ and continues with relations of the form $\Phi^l_{ij,k} = G_{ij}V_{ij}(\Phi_{ik}+\Phi_{jk})$ and $\Phi_{ij,kl} = G_{ij}V_{ik}\Phi_{kl}$, where $G_{ij}$ is the interacting two-body Green function embedded in the N-body space. Each component is written in its own Jacobi coordinates, so the pair potential acting on it is a function of a single coordinate, and the permutation symmetry of identical particles collapses the full set of components to a handful. After a partial-wave expansion, the equations become the coupled integro-differential equations of the paper, whose right-hand sides carry smooth (N$-$2)-dimensional integral kernels. The numerical engine is a spline or Lagrange local expansion combined with a tensor-trick preconditioner that makes the large linear systems solvable by iterative methods.
What would settle it
Take one fixed short-range potential, solve the same four-body bound state or scattering process with the 18 FY equations and with an independent method such as direct diagonalisation or hyperspherical harmonics, and check that energies or phase shifts agree within numerical tolerance; any systematic discrepancy would falsify the completeness claim. The same check for the n-4He five-body system would probe the imported N=5 equations.
Extended reading notes
Core claim
The paper's central claim, stated outright in the summary, is that the Faddeev-Yakubovsky equations provide a mathematically rigorous approach for the full solution of the N-body problem. Concretely, the claim is that the recursive decomposition of the wavefunction, starting from $\Psi = \sum_{i<j} \Phi_{ij}$ for interacting pairs and splitting each $\Phi_{ij}$ by repeated insertion of the pair Green function, produces closed systems of 3, 18, and 180 coupled equations for N=3, 4, and 5, and that after projection onto partial waves these become the solvable integro-differential systems (22), (35), and (44). The paper presents the N=3 case as rigorous, presents the N=4 construction explicitly while noting it is given without any formal demonstration, and imports the N=5 equations from earlier derivations. Its evidence that the machinery works is that the N=5 equations have recently been solved for n-4He scattering and for the 5H resonance. In the paper's view this makes FY the reference method for disentangling genuine interaction physics from numerical artefacts in few-body scattering.
Load-bearing premise
The load-bearing assumption is that the recursive splitting of the wavefunction into FY components is complete and closed, meaning that the 18 four-body and 180 five-body equations capture every physical solution of the Schrödinger equation, since the paper gives the four-body construction without a formal proof.
Editorial extensions
If this is right
- For any short-range pairwise interaction, the N=5 FY equations give numerically controlled scattering observables and resonance parameters; the reported n-4He and 5H calculations are the first instances.
- Experimental mismatches in few-body observables can now be interpreted as failures of the interaction model, since the numerical error is controlled by grid density and partial-wave number.
- The same recursion formally extends to N=6, but the paper's Table 1, with 2,700 equations and roughly 15 independent components for identical particles, shows the obstacle is computational scale, not principle.
- The configuration-space boundary conditions make the FY form especially useful when many scattering channels are open, including break-up, where a single Schrödinger function cannot easily carry all the required asymptotics.
Reading between the lines
- Inference: the same Merkuriev-style splitting of long-range Coulomb forces that the paper describes for N=3 should carry over to N=4 and N=5, extending the method to charged-particle scattering; the paper only discusses Coulomb regularisation for the three-body case.
- Inference: because the component structure is purely combinatorial, one could generate the FY equations for arbitrary N algorithmically from partition chains; the practical frontier would then be the linear-algebra cost, not the derivation.
- Inference: a direct comparison of the 18 four-body FY equations with an independent solution method on a benchmark potential would settle the completeness question the paper leaves open for N=4.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript is a review of the Faddeev-Yakubovsky (FY) equations for nonrelativistic N-body quantum scattering, covering N=3, 4, and 5 in configuration space. The authors present the recursive decomposition of the wave function into FY components associated with partition chains, the resulting coupled integro-differential equations after partial-wave expansion, and a unified numerical method based on spline collocation with a tensor-trick preconditioner. They also summarize recent applications to 5-body systems (n-4He scattering and 5H resonances) and discuss the combinatorial scalability limits of the approach, including a table of the growth of the number of equations and amplitudes with N.
Significance. The paper fills a useful niche as a compact pedagogical account of the FY formalism for N up to 5, including explicitly the N=5 equations (eqs. 41-44) that are otherwise scattered in the literature. Its strengths are the clear tree diagrams for partition chains, the unified presentation of the numerical scheme, and the honest discussion of the combinatorial growth in Table 1. The manuscript is a review rather than a new derivation; the ''mathematically rigorous'' claim therefore rests on the cited literature, and the authors are transparent about this in Sections 4 and 5. If the identified technical typos are corrected, the paper will serve as a reliable reference for practitioners entering the field.
major comments (2)
- [Section 4, eqs. (31) and (33)] The H-type component is defined with a V_{ik} kernel in eqs. (31) and (33), but the explicit N=4 derivation in eqs. (27)-(29) uses V_{ij} (e.g., (E-H0-V12)Phi_{12,34}=V12 Phi_{34}). This inconsistency is load-bearing because substituting eq. (31) into eq. (32) does not reproduce the displayed right-hand sides; the system would not close. Please correct the potential index in (31) and (33) to V_{ij}, or justify a different convention explicitly.
- [Section 4 and Section 5] The N=4 system is introduced 'without any formal demonstration' and the N=5 system is extended by analogy, yet Section 7 claims these equations provide a 'mathematically rigorous approach for the full solution of the N-body problem.' The completeness and closure of the recursive decompositions (31) and (38)-(40)—that every partition chain appears exactly once and that summing the FY components recovers the Schrödinger solution—are the load-bearing properties for this claim. To make the review more self-contained, I recommend adding a short paragraph stating this property explicitly for N=4 and N=5, either with a brief counting argument (e.g., the number of components in (32) and (41)) or with a precise pointer to the proofs in Refs. [5] and [27-29].
minor comments (7)
- [Section 2, eq. (1)] The potential sum is written as V12+V23+V23; it should be V12+V13+V23.
- [Section 3, eq. (19)] The second component is written as Phi_2(x2, y3); the variable y3 should be y2, consistent with the coordinate notation for the (13) pair.
- [Figure 2 caption] The caption states the Jacobi coordinates for the H component as (xH, xH, xH); this should presumably be (xH, yH, zH) or another triple of distinct coordinates.
- [Section 5, eq. (44)] The effective potential ve_alpha is written with minus signs before the angular momentum barrier terms, whereas eqs. (22) and (36) use plus signs in the same convention. The minus signs would place the barriers inside the operator incorrectly; they should likely be plus signs.
- [Section 6, after eq. (48)] The sentence 'starting with a trial value lambda_0 and an initial guess x_0. One can show...' is followed by the fragment 'k = 0, 1,....nite.' which appears to be a typo for 'k = 0,1,... finite' or similar; please reformat.
- [Table 1] The last row for general N is garbled in the typesetting ('N N!(N−1)!/2N−1 Int(2(N−1)!/(π/2)N)'); please reformat the formula and the heading so that the entries for N=6 and the general expressions are legible.
- [Section 4, eq. (31)] The superscript 'l' on the H-type component Phi^l_{ij,kl} is inconsistent with the notation in eq. (27) and eq. (32), where the H-type component is written without a superscript. This is presumably a typo and should be corrected to Phi_{ij,kl}.
Circularity Check
No circularity: review paper with no fitted parameters or predictions; self-citations document prior work and are not load-bearing.
full rationale
This is a review article rather than a derivation of new predictions. The central construction is an exact rewriting: the Faddeev components are defined by Phi_ij = G0 V_ij Psi (eq. 17), from which the coupled system (18) follows algebraically; the N=4 and N=5 decompositions (eqs. 31, 39-41) recursively define FY components in terms of cluster partitions, so there is no fitted input or predicted quantity that reduces to an input by construction. The paper's citations to the authors' own earlier work (refs. 6, 7, 8) report prior numerical solutions of the 5-body equations, but these are not used as load-bearing premises to derive the equations; they are the reported state of the art. The explicit admissions that the N=4 derivation is presented 'without any formal demonstration' (Section 4) and that the N=5 case is built by analogy (Section 5) are genuine rigor/self-containedness gaps: the equivalence of the recursively defined FY system to the Schroedinger equation is cited to Yakubovsky [5], Sasakawa [25], and Merkuriev et al. rather than proved here. That is a missing proof, not circularity under the definitions used here, because the cited results are external (not self-citations of the present authors) and the review does not define any quantity in terms of the claim it supports. Accordingly, the circularity score is 0.
Assumptions & free parameters
assumptions (4)
- domain assumption The N-body interaction is a sum of pairwise potentials V = Σ_{i<j} V_{ij}.
- domain assumption Short-range pair interactions make the FY components decouple asymptotically, enabling boundary conditions; Coulomb forces require Merkuriev's artificial splitting into short-range and long-range parts.
- domain assumption The recursive decomposition into FY components (eqs 31 and 40) produces a complete, closed system equivalent to the Schrödinger equation for N=4 and N=5.
- domain assumption The partial wave expansion (eqs 21, 42) and the local polynomial basis (45) converge, so truncation gives a controlled approximation.
Cite this review
Pith. "Pith review of The Faddeev-Yakubovsky symphony." pith.science (2026). https://pith.science/paper/HAGVXYKR
@misc{pith2026190804861,
author = {Pith},
title = {Pith review of: The Faddeev-Yakubovsky symphony},
year = {2026},
howpublished = {\url{https://pith.science/paper/HAGVXYKR}},
note = {Machine review of arXiv:1908.04861}
}
read the original abstract
We briefly summarize the main steps leading to the Faddeev-Yakubovsky equations in configuration space for N=3, 4 and 5 interacting particles.
Reference graph
Works this paper leans on
- [5]
-
[1]
L. D. Faddeev, JETP 39, 1459 (1960)
work page 1960
- [2]
-
[3]
L.D. Faddeev, Mathematical aspects of the three-body pro blem in the quantum scattering theory Israel Program for Scientic Translations, Jerusalem 1965
work page 1965
-
[4]
L.D. Faddeev and S.P. Merkuriev, Quantum Scattering Theo ry for several particle sys- tems, Kluwer Academic Publishers (1993)
work page 1993
- [6]
- [7]
- [8]
Show all 49 references
-
[9]
Merkuriev, Theoretical and Mathematical Physics 8, 798 (1971)
S.P. Merkuriev, Theoretical and Mathematical Physics 8, 798 (1971). The Faddeev-Yakubovsky symphony 19 T able 1 Scalability of the Faddeev-Yakubovsky scheme as a function of the number of interacting particles N.: ne is the number of equations,(eq ual to the number of FY compo...
1971
-
[10]
Merkuriev, Sov
S.P. Merkuriev, Sov. J. Nucl. Phys. 19, 22 (1974)
1974
-
[11]
Gignoux, A
C. Gignoux, A. Laverne, Phys. Rev. Lett 29, 436 (1972)
1972
-
[12]
Gignoux, A
C. Gignoux, A. Laverne, S.P. Merkuriev, Phys. Rev. Lett 3 3, 1350 (1974)
1974
-
[13]
Merkurev, C
S.P. Merkurev, C. Gignoux, A. Laverne, Annals Phys. 99, 3 0 (1976)
1976
-
[14]
Gl¨ ockle, The Quantum Mechanical Few-Body Problems, Springer-Verlag (1983)
W. Gl¨ ockle, The Quantum Mechanical Few-Body Problems, Springer-Verlag (1983)
1983
-
[15]
Merkuriev, Ann
S.P. Merkuriev, Ann. Phys. (N.Y.) 130, 395 (1980)
1980
-
[16]
Merkuriev, Acta Physica Austriaca, Supplementum X XIII 65 (1981)
S.P. Merkuriev, Acta Physica Austriaca, Supplementum X XIII 65 (1981)
1981
-
[17]
Doleschall and Z
P. Doleschall and Z. Papp, Phys. Rev. C, 72(4), 044003 (20 05)
-
[18]
Kvitsninsky, J
A. Kvitsninsky, J. Carbonell, C. Gignoux, Phys. Rev. A46 , 1310 (1992)
1992
-
[19]
M Valdes et al., Phys. Rev. A Phys. Rev. A 97, 012709 (2018)
2018
-
[20]
V. A. Gradusov et al., Journal of Physics B: Atomic, Molec ular and Optical Physics, 52 , 055202 (2019)
2019
-
[21]
Noyes, H
P. Noyes, H. Fideldey, Proc. of Texas A.M. Conference, Ed . by J. Gillespie and J. Nuttall, 195-293, New-York Benjamin (1968)
1968
-
[22]
Faddeev, Proceedings of the First International Co nference of the three body problems in Nuclear and Particle Physics, Birmingham, 8-10 july (1969), Ed
L.D. Faddeev, Proceedings of the First International Co nference of the three body problems in Nuclear and Particle Physics, Birmingham, 8-10 july (1969), Ed. by J.S.C. Mc Kee and P.M. Rolph, North-Holland, London (1970), pag 154
1969
-
[23]
W einberg, Phys
S. W einberg, Phys. Rev. 130, 776 (1963); S. W einberg, Phys. Rev. 131, 440 (1963)
1963
-
[24]
Narodetsky, O.A.Yakubovksy, ITEP-1-1980 preprin t (1980)
I.M. Narodetsky, O.A.Yakubovksy, ITEP-1-1980 preprin t (1980)
1980
-
[25]
Sasakawa, Prog
T. Sasakawa, Prog. Theor. Phys. Suppl 61, 149 (1977)
1977
-
[26]
Lazauskas, Application of the complex scaling method in quantum scattering theory, Habilitation Thesis (2019), arXiv:1904.04675 [nucl-th]
R. Lazauskas, Application of the complex scaling method in quantum scattering theory, Habilitation Thesis (2019), arXiv:1904.04675 [nucl-th]
2019 arXiv
-
[27]
Merkuriev and S.L
S.P. Merkuriev and S.L. Yakovlev, Akad. Nauk (SSSR) 261, 591 (1982)
1982
-
[28]
Merkuriev and S.L
S.P. Merkuriev and S.L. Yakovlev, , Theor. Math. Phys. 56 , 673 (1983)
1983
-
[29]
Merkuriev, S.L
S.P. Merkuriev, S.L. Yakovlev and C. Gignoux, Nucl. Phys . A431, 125 (1984)
1984
-
[30]
Ciesielski, J
F. Ciesielski, J. Carbonell, Phys. Rev. C58, 58 (1998); R . Lazauskas, PhD Univ. J. Fourier Grenoble, (2003); R. Lazauskas, J. Carbonell, Phys . Rev. C70, 044002 (2007)
1998
-
[31]
Nogga et al., Phys.Rev
A. Nogga et al., Phys.Rev. C65, 054003 (2002)
2002
-
[32]
Uzu et al., Phys
E. Uzu et al., Phys. Rev. C 68, 061001 (2003)
2003
-
[33]
Deltuva and A
A. Deltuva and A. C. Fonseca, Phys. Rev. C 86, 011001, (201 2)
-
[34]
Lazauskas, Phys
R. Lazauskas, Phys. Rev. C 91, 041001 (2015)
2015
-
[35]
Saad, Iterative Methods for Sparse Linear Systems
Y. Saad, Iterative Methods for Sparse Linear Systems. Ph iladelphia SIAM (2003)
2003
-
[36]
N. W. Schellingerhout, L. P. Kok, and G. D. Bosveld Phys. R ev. A 40, 5568 (1989)
1989
-
[37]
Lazauskas PhD, [http://tel.ccsd.cnrs.fr/document s/archives0/00/00/41/78/], Univ
R. Lazauskas PhD, [http://tel.ccsd.cnrs.fr/document s/archives0/00/00/41/78/], Univ. J. Fourier, Grenoble, 2003
2003
-
[38]
Alt E.O., Grassberger P., Sandhas W., Nucl. Phys. B2, 167 (1967)
1967
-
[39]
Deltuva and A
A. Deltuva and A. C. Fonseca Phys. Rev. C 76, 021001(R) (20 07); A. Deltuva and A. C. Fonseca Phys. Rev. C 95, 024003 (2017)
2017
-
[40]
Glockle, H
W. Glockle, H. Witala, D. Huber, H. Kamada, J. Golak, Phys .Rept. 274, 107 (1996); W. Glockle, H. Kamada, Phys. Rev. Lett. 71, 971 (1993); W. Glo ckle, H. Kamada, Nucl. A 560, 541 (1993). A. Nogga, D. Huber, H. Kamada, W. Glockle, P LB 409, 19 (1997). 20 Rimantas Lazauskas*, J...
1996
-
[41]
B. S. Pudliner et al., Phys. Rev. C 56, 1720 (1997); K.M. No llett et al., Phys. Rev. Lett. 99, 022502 (2007)
1997
-
[42]
Gattobigio, A
M. Gattobigio, A. Kievsky and M. Viviani, Few-Body Syst. 54, 657 (2013)
2013
-
[43]
Kievsky et al., J.Phys
A. Kievsky et al., J.Phys. G35, 063101 (2008)
2008
-
[44]
Hiyama, Y
E. Hiyama, Y. Kino and M. Kamimura, Prog. Theor. Phys. 51, 223 (2003)
2003
-
[45]
Navr´ atil et al., Phys
P. Navr´ atil et al., Phys. Rev. C 61, 044001 (2000); P. Nav r´ atil et al., Physica Scripta 5, 053002 (2016)
2000
-
[46]
Hupin, S
G. Hupin, S. Quaglioni, P. Navrtil, Nature Commun. 10, 35 1 (2019)
2019
-
[47]
Efros, W
V.D. Efros, W. Leidemann, and G. Orlandini, Phys. Lett. B 338, 130 (1994)
1994
-
[48]
Giuseppina Orlandini, in An Advanced Course in Computat ional Nuclear Physics, Lec- ture Notes in Physics 936 (2017) pag. 263
2017
-
[49]
Kievsky et al., Nucl.Phys
A. Kievsky et al., Nucl.Phys. A577, 511 (1994)
1994
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.