{"id":"7b404bde-93f7-4c2e-837a-d2fc0bddfabf","arxiv_id":"2505.24569","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The physical neutron-neutron scattering length is too far from unitarity for 3n and 4n formation to show the predicted unitarity-modified threshold exponents.","lead":"This paper asks whether the neutron-neutron interaction is strong enough to make the energy dependence of 3- and 4-neutron reactions look like the unitarity limit. Using hyperspherical and JWKB methods, it finds that the physical scattering length is about 3 to 5 times too small, so the ordinary threshold law should be observed instead.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Concern: The 4n JWKB exponent is only shown for physical a_nn, never for tuned 3-5x larger scattering lengths, so the 3-5x estimate for 4n is unsupported; the physical-a_nn non-observation is supported but rests on a single 20% deviation in Fig. 6(b).","rationale":"I read the paper in good faith and identify the 4n analysis as the least secure link in the central claim. The 3n analysis is cross-validated against Faddeev calculations (Golak et al., Fig. 2), and that agreement is the paper's strongest evidence that JWKB tunneling in the lowest adiabatic potential captures the near-threshold energy dependence at physical a_nn. The 4n analysis, by contrast, is explicitly a gedanken experiment with no independent check. The headline conclusion is that a_nn is not sufficiently large for the unitarity exponent to be observed in 3n or 4n; for 4n, this rests on the single 20%-above-unitarity deviation in Fig. 6(b). The paper states the 3-5x estimate for both systems, but the 4n section does not show the tuned-scattering-length analogs of Fig. 3. If the adiabatic single-channel description or the JWKB treatment breaks down at tuned unitary parameters, the quantitative 3-5x estimate for 4n is unsupported. That is a load-bearing concern for the strength of the conclusion, although not for the qualitative direction. My proposed test is narrow and computational: compute tuned-4n γ(E) curves and check their sensitivity to the fitted coefficients and ρsmall. The reader already partially identified this by noting the 4n section lacks an independent cross-check. I agree with that assessment; my concrete test asks whether the 4n conclusion itself is robust within the same JWKB framework, not whether the method is universally validated.","tokens_in":13880,"tokens_out":1610,"duration_ms":16155,"concrete_test":"Recompute the 4n JWKB effective exponent γ(E) for s-wave scattering lengths tuned to 5a_nn, 10a_nn, and 50a_nn, using the same AV8' potential and hyperradial fitting (Eq. 14); if those curves do not approach γ=2.517 by E≈2 MeV, or if they depend sensitively on the fitted C_2-C_6 coefficients or on ρsmall in the 1.25-5 fm range, then the 3-5x estimate is not uniformly supported and the 4n section should be explicitly re-scoped.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The physical-nucleus conclusion for 4n rests on Fig. 6 alone, where γ(E) for the physical a_nn stays above the unitarity value (2.517) by about 20% at the upper end of Region II. The abstract's stronger claim that a_nn is not yet sufficiently large for the 3n or 4n systems to demonstrate the unitarity threshold exponent is supported for 4n only by that single panel, which does not include tuned scattering lengths. Consequently, the paper does not demonstrate that the parameter window is insensitive to the adiabatic-JWKB approximation in the case where the 4n system is actually close to unitarity. The absence of tuned-4n curves is a gap between the stated 3-5x estimate and the evidence shown, not an internal inconsistency.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper addresses the question of how close to the S-wave unitarity limit a two-body scattering length must be for the unitarity-modified Wigner threshold exponent to be observable in low-energy reactive processes producing three or four neutrons. Using the adiabatic hyperspherical representation, the authors compute the lowest hyperradial potential curves for the 3n (Jπ=3/2−) and 4n (Jπ=0+) systems, extract the JWKB tunneling amplitude, and define an energy-dependent effective exponent γ(E) via Eq. (11). They find that at the physical neutron-neutron scattering length ann=-18.5 fm, the low-energy exponent is the ordinary non-unitary value (γ=3 for 3n, γ=11/2 for 4n) in Region I, and although it decreases in Region II, it does not reach the unitarity values (1.7727 for 3n, 2.517 for 4n); only for scattering lengths roughly 3–5 times more negative does the unitarity exponent emerge over an appreciable energy range. The 3n result is validated by comparison with Faddeev calculations of the capture reactions 3H(π−,γ)3n and related processes; the 4n result is presented as a gedanken experiment based solely on the JWKB analysis of the lowest 0+ adiabatic curve.","tokens_in":14096,"tokens_out":6361,"duration_ms":72700,"significance":"If the analysis holds, the paper provides a quantitative answer to an open question and corrects an over-reading of earlier work: the apparent agreement of the 3n capture rate with the unitarity exponent in a limited energy window is fortuitous, and the physical nn scattering length is too far from unitarity for the modified threshold law to be observable. The paper's strengths include the direct comparison of the JWKB scale factor with the Faddeev capture-rate data for 3n, the use of well-established adiabatic hyperspherical methods, and the consistency of the unitarity-limit exponents with the conformal field theory results of Hammer and Son. The three-region decomposition of the threshold energy dependence is a useful conceptual framework. The main weakness is that the 4n conclusions rest on a single JWKB channel without an independent cross-check or tuned-scattering-length calculations, so the quantitative 3–5 times estimate for 4n is not directly demonstrated.","major_comments":[{"comment":"The 4n analysis displays only the physical scattering-length case: Fig. 6(a) and (b) show the JWKB density and effective exponent for ann only, with no curves for 5×, 10×, or 50× ann as in Fig. 3. The Introduction and Conclusion nevertheless state that for the 3n or 4n system the unitarity exponent would be observable only if the scattering length were approximately 3–5 times larger. For 4n this quantitative statement is therefore an extrapolation by analogy, not a result of the calculation shown. Please either add tuned 4n curves (e.g., a panel showing γ(E) for several multiples of ann) or explicitly restrict the 3–5× estimate to 3n and state the 4n statement as a conjecture.","section":"3.2, Fig. 6"},{"comment":"The conclusion that physical ann is insufficient for the 4n system rests on a single JWKB calculation in the lowest 0+ adiabatic curve, and the deviation from the unitarity exponent is only about 20% at the upper end of Region II. Unlike the 3n case, there is no independent Faddeev or other few-body calculation to validate the adiabatic single-channel and JWKB approximations for 4n. Please provide a robustness check—for example, a numerical solution of Eq. (12) compared with the JWKB result, a sensitivity study with respect to ρ_small, or a discussion of the effect of the fitted coefficients C2–C6 on γ(E)—so that the 4n negative conclusion is supported by more than one figure.","section":"3.2, Fig. 6(b)"}],"minor_comments":[{"comment":"The word \"quantities\" in the abstract should be \"quantifies.\"","section":"Abstract"},{"comment":"The spelling of \"Faddeev\" is inconsistent; both \"Fadeev\" and \"Faddeev\" appear (e.g., §3.1 and the Conclusion).","section":"Throughout"},{"comment":"The notation \"4He(µ−µ−,νν)4n\" is unclear: specify the neutrino flavors and use a more standard notation such as \"νμνμ\" or \"2νμ\" to avoid ambiguity.","section":"3.2"},{"comment":"The statement that the scattering length \"would have to be at least 3–5 times more negative\" is based on visual inspection of Fig. 3 without an explicit criterion; defining \"appreciable energy range\" (e.g., the range where γ(E) is within 20% of the unitarity value) would make the claim more quantitative.","section":"3.1, Fig. 3"},{"comment":"The caption refers to \"JWKB density\" but the text does not define whether this is τ_JWKB from Eq. (10) or the squared wave function |Ψ|²; please clarify.","section":"Fig. 2"}],"recommendation":"major_revision","confidential_remarks":"The paper is at the boundary between nuclear few-body physics and atomic/molecular physics; for a hep-ph journal the fit is marginal, though the connection to Hammer and Son's 'unnuclear physics' makes it topical. The main quantitative claim for the 4n system needs to be backed by tuned-scattering-length curves and a robustness check; if the authors supply those, the paper would be a solid contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline: this paper asks the right question—how close to unitarity do you actually have to be before the modified threshold exponents show up—and gives a clear answer for 3n and 4n: the physical neutron-neutron scattering length is not close enough. The negative conclusion looks right for 3n, and the 4n argument is plausible but thinner.\n\nWhat's genuinely new is the observability estimate. Previous work (Hammer-Son, the Greene group's own papers) established what the unitarity-limit exponents are; nobody had quantified how far you are from seeing them at finite scattering length. The paper's main quantitative result is the 3-5x estimate: you'd need a_nn somewhere around -55 to -90 fm for the unitarity exponent to hold over an appreciable energy range. That's a useful, concrete claim.\n\nThe method is straightforward but appropriate. They compute JWKB tunneling amplitudes in the lowest adiabatic hyperradial potential, and for the 3n case they validate this against the Faddeev capture-rate calculations of Golak et al. The agreement in Region (I) and (II) is reassuring, and they make a good point that the apparent agreement with the unitarity exponent in the 1-3 MeV range is fortuitous—the effective exponent gamma(k) never gets within 20% of 1.7727 at the physical a_nn. That's the strongest part of the paper.\n\nThe soft spots are real but not fatal. First, the \"3-5 times\" number is qualitative. There's no explicit criterion for what counts as \"observable\" (how broad an energy range, what statistical precision), so the estimate is more a useful rule of thumb than a sharp prediction. Second, the 4n analysis is a gedanken experiment, as the authors admit, with no Faddeev cross-check. The stress-test note is on target: Fig. 6(b) only shows the physical a_nn curve, not tuned values, so the 3-5x extrapolation for 4n is an analogy from the 3n case rather than a demonstrated result. That said, the actual claim being made for 4n—that gamma stays ~20% above the unitarity value—is directly visible in the figure, so the main negative conclusion for 4n is supported by the evidence shown. The tuned-4n gap is a limitation, not a fatal flaw. Third, they rely on potential curves from their own prior work; the expansion coefficients are given but the fitted forms and the Gaussian-interaction parameters are only partially specified, which will make reproduction annoying but not impossible.\n\nBottom line: this is a solid, honest paper that settles the observability question in the negative for the physical n-n scattering length within the adiabatic-JWKB framework. It deserves a serious referee. I'd send it to review, and I'd probably cite it.","headline":"This paper gives a credible negative answer to a previously unquantified question: at the physical n-n scattering length, the unitarity-modified threshold exponents are not observable in 3n or 4n formation, and you'd need roughly 3-5 times larger |a_nn| to see them.","tokens_in":14619,"tokens_out":2208,"would_cite":true,"duration_ms":23526,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"At the physical neutron-neutron scattering length, the unitarity-modified threshold law is not observable for 3n or 4n; only with |a_nn| enlarged by roughly a factor of 3-5 would the modified exponent be seen over an appreciable energy…","keywords":["neutron-neutron scattering length","Wigner threshold law","unitarity","three-neutron system","four-neutron system","adiabatic hyperspherical representation","JWKB tunneling amplitude","effective threshold exponent"],"falsifier":"Compute the 4n effective exponent gamma(k) with physical a_nn using a fully converged four-body continuum method that does not rely on the single-channel adiabatic potential curve; if gamma approaches 2.517 over an appreciable energy range before node structure sets in, the paper's 4n conclusion is wrong. For 3n, measure the pion-capture 3n spectrum below about 1 MeV and extract the exponent directly; a clear approach to 1.7727 there would contradict the claimed ordinary Wigner behavior.","tokens_in":13595,"feed_emoji":"⚛️","tokens_out":9871,"duration_ms":107760,"temperature":0.7,"pith_summary":"The paper asks how close a few-neutron system must be to s-wave unitarity before the predicted unitarity-modified threshold law actually shows up in reaction rates. It answers that the physical neutron-neutron scattering length is not close enough: for both 3n and 4n formation, the low-energy exponent is the ordinary Wigner-law value, and the unitarity exponent appears only when the scattering length is made roughly 3-5 times more negative. This matters because recent theoretical work has identified apparent near-threshold behavior in 3n capture spectra with the unitarity limit; the paper argues that apparent agreement is accidental and that existing calculations are consistent with the ordinary threshold law once the energy is low enough.","feed_headline":"Unitarity exponent stays hidden for 3n and 4n reactions","feed_subtitle":"At -18.5 fm the ordinary Wigner law wins; |a| must be 3-5 times larger to reveal the unitarity exponent.","key_machinery":"The object carrying the argument is the long-range coefficient of the lowest adiabatic hyperradial potential curve W0(rho), written as l_eff(l_eff+1) $hbar^{2}$/(2 mu $rho^{2}$). The value of l_eff controls the Wigner threshold exponent |T|^2 ~ $k^{{2l_eff+1}}$: the non-interacting value applies for any finite scattering length sufficiently close to threshold (Region I), while the unitarity-limit value would apply only at exact resonance. The paper's diagnostic is the JWKB tunneling amplitude in this potential, evaluated at small hyperradius, and the effective exponent gamma(k) = k d ln|psi|/dk, which tracks the transition between the two exponent limits as k rises past ~1/|a|.","core_discovery":"The paper's central claim is that the unitarity-modified Wigner threshold law predicted for short-range reactive processes producing three or four neutrons cannot be seen at the physical neutron-neutron scattering length a_nn = -18.5 fm. Using the lowest adiabatic hyperradial potential curves for the 3n J^pi = 3/2^- and 4n 0+ symmetries, the authors compute the JWKB tunneling amplitude and the effective threshold exponent gamma(k) = k d ln|psi|/dk. They find that at sufficiently low energy the exponent is the ordinary non-unitarity value (gamma = 3.0 for 3n, gamma = 11/2 for 4n), and that the approach toward the unitarity values (gamma = 1.7727 for 3n, gamma = 2.517 for 4n) begins only near k ~ 1/|a|, well inside the energy window where a short-range formation process could be observed. Only when the scattering length is increased in magnitude by a factor of about 3-5 does the unitarity exponent become visible over an appreciable energy range.","pith_inferences":["Inference: the same three-region picture implies that any experimental claim of a unitarity-modified threshold law in nuclear reactions should be accompanied by a check that the energy window lies beyond k ~ 1/|a| and that the extracted exponent is stable over an appreciable range.","Inference: ultracold atoms with a Feshbach resonance could test the predicted 3-5 tuning requirement directly, since the analysis should transfer to any short-range s-wave system with the same effective-centrifugal-constant structure.","Inference: if a future exact four-body continuum calculation at physical a_nn finds significant deviations from the single-channel JWKB tunneling exponent, the paper's 4n conclusion would need revision, while its 3n conclusion would stand."],"forward_implications":["Existing 3n capture data in the near-threshold region should be read as following the ordinary Wigner law E^3 with l_eff = 5/2, not as evidence for the unitarity-modified exponent.","At the physical scattering length, the effective exponent for 4n formation stays well above the unitarity value gamma = 2.517 throughout the threshold energy regions where a simple power law applies.","The unitarity-modified threshold exponent becomes observable only when the neutron-neutron scattering length is tuned to about 3-5 times the physical magnitude; for any finite a there is a low-energy region where the ordinary Wigner law always wins.","There is a universal transition in the effective exponent near k ~ 1/|a|: below it the non-interacting exponent holds, above it the exponent moves toward the unitarity value."],"supporting_citations":[{"why":"Defines the original Wigner threshold law whose finite-scattering-length form is the baseline exponent.","marker":"[1]"},{"why":"Conformal-field-theory prediction of unitarity-modified threshold exponents that the present paper tests against the physical scattering length.","marker":"[2]"},{"why":"Provides the adiabatic-hyperspherical effective centrifugal constants for three- and four-fermion symmetries at unitarity that set the modified exponents.","marker":"[4]"},{"why":"Exact three-body calculation of muon capture on 3H whose near-threshold final-state spectrum is compared with the JWKB tunneling density.","marker":"[9]"},{"why":"Exact three-body pion-capture calculation giving the differential capture rate dGamma/dE used to identify Region (I) and Region (II) behavior.","marker":"[10]"},{"why":"Supplies the low-energy three- and four-neutron adiabatic potential curves and the unitarity-limit effective centrifugal constants used in the JWKB analysis.","marker":"[12]"},{"why":"Establishes the nonresonant density-of-states enhancement and potential curves for 3n and 4n at low energies that underlie the threshold-region analysis.","marker":"[14]"},{"why":"Realistic nucleon-nucleon interaction that determines the physical neutron-neutron scattering length and the hyperradial potential curves.","marker":"[16]"}],"fun_headline_variants":["Unitarity threshold exponent requires 3-5x larger neutron length","Wigner law prevails: unitarity exponent invisible at physical a_nn","3-5x larger neutron scattering length needed to see unitarity exponent","Neutron-neutron length far from unitarity threshold for 3n and 4n"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The three-neutron part is cross-checked against exact three-body calculations, but the four-neutron part is a thought experiment that assumes the energy dependence of any four-neutron formation process is set by tunneling through a single lowest potential curve; if that single-channel tunneling description fails, the four-neutron conclusion is unsupported.","fun_headline_variants_meta":{"raw":{"variants":["Unitarity threshold exponent requires 3-5x larger neutron length","Wigner law prevails: unitarity exponent invisible at physical a_nn","3-5x larger neutron scattering length needed to see unitarity exponent","Neutron-neutron length far from unitarity threshold for 3n and 4n"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000799,"raw_usage":{"total_tokens":3488,"prompt_tokens":894,"completion_tokens":2594,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":510,"completion_tokens_details":{"reasoning_tokens":2509}},"tokens_in":510,"tokens_out":2594,"duration_ms":24051,"temperature":1.0,"reasoning_tokens":2509,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:17:35.665644+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the 4n effective exponent gamma(k) with physical a_nn using a fully converged four-body continuum method that does not rely on the single-channel adiabatic potential curve; if gamma approaches 2.517 over an appreciable energy range before node structure sets in, the paper's 4n conclusion is wrong. For 3n, measure the pion-capture 3n spectrum below about 1 MeV and extract the exponent directly; a clear approach to 1.7727 there would contradict the claimed ordinary Wigner behavior.","supporting_citations":[{"cited_title":"Physical Review C 94, 034002 (2016) https://doi.org/10.1103/physrevc.94.034002 17","cited_arxiv_id":null,"evidence_quote":"Exact three-body calculation of muon capture on 3H whose near-threshold final-state spectrum is compared with the JWKB tunneling density."},{"cited_title":"Physical Review C 98, 054001 (2018) https://doi.org/10.1103/ PhysRevC.98.054001","cited_arxiv_id":null,"evidence_quote":"Exact three-body pion-capture calculation giving the differential capture rate dGamma/dE used to identify Region (I) and Region (II) behavior."}],"review_version":1}