{"id":"4643d613-ab1d-4042-af59-38f3cea1fb62","arxiv_id":"1908.04861","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":2.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A review of the Faddeev-Yakubovsky equations in configuration space for N=3, 4, and 5 particles, including the numerical methods used to solve them.","lead":"This paper reviews the Faddeev-Yakubovsky equations, a mathematical framework for solving quantum few-body problems with three to five particles. It is a useful introduction for researchers entering few-body physics, though it contains no new results or data.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Self-admitted lack of formal derivation for N=4/N=5 leaves the 'rigorous' claim dependent on cited proofs; a closure check of eqs. (32) and (41) would settle it.","rationale":"The reader's weakest assumption identified exactly the same load-bearing point: the recursive FY components for N=4 and N=5 are claimed to be complete and closed, but the paper explicitly provides no formal demonstration and relies on citations to Yakubovsky and Sasakawa. My read agrees with that assessment. The concern is genuine but does not force a stronger verdict than the reader's CONDITIONAL: the paper is a review, not a new proof, and the cited literature is the standard source of the rigorous derivation. The concrete closure test would settle whether the concern is merely a missing exposition or an actual error in the displayed N=5 system. I did not find an independent internal inconsistency that would overturn the equations, but the definitional ambiguity in eq. (33) and the unproved completeness of eqs. (32) and (41) support keeping the acceptance conditional on verification.","tokens_in":14730,"tokens_out":13047,"duration_ms":120054,"concrete_test":"Perform the closure check that the paper omits. For N=4, substitute the definitions (31) into the right-hand sides of eq. (32), sum all 18 components, and verify that the total equals G_0 V Psi from the original Schrödinger equation. For N=5, substitute (39)-(40) into eq. (41), sum over all 180 components, and compare with the N=5 analogue of eq. (18). If each summation reproduces the original Faddeev equations for arbitrary short-range pairwise potentials, the completeness/closedness concern is resolved; if not, eqs. (41) are not equivalent to the stated N-body problem and the N=5 claim would need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that the configuration-space FY equations are a mathematically rigorous route to the N-body problem, including the cited N=5 solutions. The load-bearing condition is that the recursive decompositions in eqs. (31) and (38)-(40) are complete and closed: every partition chain appears exactly once, and summing the FY equations reproduces the original Schrödinger equation. The paper itself states in Section 4 that the N=4 derivation is presented 'without any formal demonstration', and Section 5 extends to N=5 by analogy. Moreover, eq. (41) is written using the \\psi aggregates of eqs. (39)-(40), so the reader cannot directly see whether the system is closed without repeating the derivation or consulting refs. [5,25]. A related ambiguity appears in eq. (33), where the H-type component is displayed with a V_{ik} kernel although eq. (29) and the preceding definitions imply V_{ij}; this needs correction or clarification. If a partition chain were omitted or double-counted, the numerical solutions in refs. [7,8] would not solve the stated scattering problem. This is a gap in self-containedness rather than a demonstrated error, but it is exactly the point on which the strongest claim rests.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":14965,"tokens_out":7256,"duration_ms":67174,"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":[{"comment":"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":"Section 4, eqs. (31) and (33)"},{"comment":"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].","section":"Section 4 and Section 5"}],"minor_comments":[{"comment":"The potential sum is written as V12+V23+V23; it should be V12+V13+V23.","section":"Section 2, eq. (1)"},{"comment":"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.","section":"Section 3, eq. (19)"},{"comment":"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":"Figure 2 caption"},{"comment":"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":"Section 5, eq. (44)"},{"comment":"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.","section":"Section 6, after eq. (48)"},{"comment":"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":"Table 1"},{"comment":"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}.","section":"Section 4, eq. (31)"}],"recommendation":"minor_revision","confidential_remarks":"The paper is a review article with standard scientific content. The main technical issue is the V_{ik} versus V_{ij} inconsistency in the H-type component definitions, which must be fixed; the other typos are cosmetic but should be cleaned up. The reliance on the authors' own earlier papers for the N=5 equations (refs. 6-8) is normal for a review, though it might be worth a brief editorial check that the presentation is not overly self-referential. Overall, the manuscript is suitable for publication after minor revisions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Take this as what it is: a short, clear review of the Faddeev-Yakubovsky equations for N=3, 4, and 5, written by people who have actually solved them. No new equations or numerical results appear; the N=3 and N=4 material is textbook, and the N=5 system is from Sasakawa and the authors' earlier papers. That is not a reason to dismiss it. The tree-diagram presentation and the scalability table give a useful picture that is hard to get from the primary papers, and the numerical-method summary is a readable entry point to practice.\n\nThe soft spots are mostly editorial, but a few are more than cosmetic. Eq. (1) lists V12+V23+V23 instead of V12+V13+V23; Eq. (19) has y3 where y2 belongs; Fig. 2's caption repeats xH three times. More substantially, Sec. 4 explicitly says the N=4 derivation is given 'without any formal demonstration,' and Sec. 5 extends to N=5 by analogy. The stress-test note is right that closure/completeness of the recursive decomposition is load-bearing. For a review, however, this is acceptable: the paper cites Yakubovsky and Sasakawa for the proofs, and the reader can chase those. It would be better if Sec. 7 said plainly that the rigor is inherited from the cited literature rather than demonstrated here. One actual apparent typo is in Eq. (33): the H-type FY component should be G_{ij} V_{ij} G0 V_{kl} to be consistent with Eqs. (28)-(29); the printed V_{ik} looks wrong and needs correction or clarification.\n\nThe citation pattern is fine. The self-citations are justified because refs. [6-8] are the papers that actually solved N=5. In short, this review deserves a serious referee, but only for light revision: fix the typos, correct Eq. (33), and soften the 'mathematically rigorous' claim so it points to the literature. I would not cite it in my own work, but I would hand it to a student entering few-body physics.","headline":"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.","tokens_in":15570,"tokens_out":2860,"would_cite":false,"duration_ms":30282,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["21.45.-v","03.65.Nk"],"model":"deepseek-v4-flash","headline":"The Faddeev-Yakubovsky decomposition turns the N-body Schrödinger equation into coupled integro-differential systems, now solved through five particles.","keywords":["Faddeev-Yakubovsky equations","few-body physics","N-body problem","configuration space","scattering theory","integro-differential equations","partial wave expansion","five-body problem"],"falsifier":"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.","tokens_in":14509,"feed_emoji":"⚛️","tokens_out":10401,"duration_ms":93275,"temperature":0.7,"pith_summary":"The paper is a guided tour of the Faddeev-Yakubovsky (FY) equations, a reformulation of the quantum-mechanical N-body problem designed to fix the failures of the bare Schrödinger equation in scattering situations. The authors show how the N-body wavefunction is split recursively into FY components, one for each way of dismantling the cluster down to a single interacting pair, and how the resulting systems for N=3, 4, and 5 become coupled integro-differential equations after a partial-wave expansion. The central claim is that these equations are mathematically rigorous and that the five-body case has now been solved numerically, so that exact few-body solutions are available with only the interparticle potential as input. A sympathetic reader would care because this separates errors coming from the interaction model from errors coming from the numerical approximation, which is exactly what the three-nucleon problem has needed since the beginning.","feed_headline":"Five-body scattering is now solvable, not just three-body","feed_subtitle":"The Faddeev-Yakubovsky decomposition turns the N-body equation into coupled equations with controlled numerical error.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Introduces the paired-component decomposition of the three-body wavefunction from which the whole recursion grows.","marker":"[1]"},{"why":"Supplies the induction scheme for arbitrary N that underlies the four- and five-body component systems.","marker":"[5]"},{"why":"Source of the 180 five-body equations reproduced as (41).","marker":"[25]"},{"why":"Provides the configuration-space differential form and boundary conditions for the four-body equations.","marker":"[27, 28, 29]"},{"why":"Recent numerical solutions of the five-body system, cited as evidence that the equations are solvable.","marker":"[6, 7, 8]"},{"why":"Defines the integral kernels entering the projected equations (22) and (35).","marker":"[30]"},{"why":"Introduces the tensor trick used as the preconditioner in the numerical method.","marker":"[36]"},{"why":"Documents the N=4 generalization of the numerical protocol and preconditioner.","marker":"[37]"}],"fun_headline_variants":["Five-body scattering now solvable via Faddeev-Yakubovsky","Faddeev-Yakubovsky equations now solve five-body scattering","From three to five particles: FY equations solved","N=5 scattering now solvable with Faddeev-Yakubovsky","Five-body challenge solved: Faddeev-Yakubovsky method"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Five-body scattering now solvable via Faddeev-Yakubovsky","Faddeev-Yakubovsky equations now solve five-body scattering","From three to five particles: FY equations solved","N=5 scattering now solvable with Faddeev-Yakubovsky","Five-body challenge solved: Faddeev-Yakubovsky method"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000764,"raw_usage":{"total_tokens":3308,"prompt_tokens":784,"completion_tokens":2524,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":400,"completion_tokens_details":{"reasoning_tokens":2430}},"tokens_in":400,"tokens_out":2524,"duration_ms":19023,"temperature":1.0,"reasoning_tokens":2430,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:30:35.056270+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the paired-component decomposition of the three-body wavefunction from which the whole recursion grows."},{"cited_title":"Yakubovsky, Sov","cited_arxiv_id":null,"evidence_quote":"Supplies the induction scheme for arbitrary N that underlies the four- and five-body component systems."},{"cited_title":"Sasakawa, Prog","cited_arxiv_id":null,"evidence_quote":"Source of the 180 five-body equations reproduced as (41)."},{"cited_title":"Ciesielski, J","cited_arxiv_id":null,"evidence_quote":"Defines the integral kernels entering the projected equations (22) and (35)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the tensor trick used as the preconditioner in the numerical method."},{"cited_title":"Lazauskas PhD, [http://tel.ccsd.cnrs.fr/document s/archives0/00/00/41/78/], Univ","cited_arxiv_id":null,"evidence_quote":"Documents the N=4 generalization of the numerical protocol and preconditioner."}],"review_version":1}