{"id":"0ebcc08a-a4b4-4f19-a81a-5127e379f04f","arxiv_id":"2504.14608","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"The Lambda_c^+ -> n pi+ pi0 decay is predicted to show peaks near 1500 MeV in the pi0 n and pi+ n invariant mass spectra, tracing the dynamically generated N(1535) resonance.","lead":"The paper calculates the decay of a charmed baryon into a neutron and two pions, predicting a clear bump in the pion-neutron mass spectrum caused by the elusive N(1535) nucleon excitation. It suggests this decay, measurable at current colliders, could distinguish between a molecule-like and a quark-core picture of N(1535).","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Line-shape claim rests on dropping Re G in the production loop; at sqrt(s) ~ 1500 MeV the stated reason (piN threshold is far away) is false, and no sensitivity test is provided.","rationale":"The paper is a standard chiral-unitary calculation, and the appearance of a peak near 1500 MeV is not surprising, since the N(1535) is put into the final-state interaction by construction. The genuinely testable part is the line-shape difference between the molecular picture and a Breit-Wigner description. That difference is most directly vulnerable at the production-loop level, not in the well-tested two-body T matrix. The reader's weakest assumption already points to the imaginary-part-only loop and the cutoff; I agree with that identification. The additional arbitrary A = B in Model II is real but secondary: it changes the comparison baseline, whereas the loop approximation changes the Model-I prediction itself. The concrete test above would settle whether the concern lands. A conditional verdict is appropriate; no rejection is warranted because the calculation is in-principle falsifiable and the missing test is well defined.","tokens_in":11671,"tokens_out":17922,"duration_ms":173102,"concrete_test":"Recompute Model I with the full loop function of Eq. (11) in the production amplitudes (Eq. 9), using qmax = 1000, 1150, and 1300 MeV, and compare peak positions and the Model-I vs Model-II line-shape difference in Fig. 10. If the peak shifts by more than about 20 MeV, or if the full-G calculation brings the Model-I line shape into agreement with the Breit-Wigner curve, the paper's central claim is not yet established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section II (after Eq. 15) states that only the imaginary part of the loop function is kept, on the ground that the piN-channel threshold is far from the N(1535) region. That ground is not correct: at sqrt(s) ~ 1500 MeV the piN channel (threshold about 1078 MeV) is open, and the real part of G(s) from Eq. (11) is not negligible. This matters because Eq. (9) is not a detached phenomenological insertion: h + h G T equals h[1 - V G]^{-1} only if the same full G appears in the production loop and in T. Replacing G by i Im(G^II) changes the phase and magnitude of the rescattering contribution, so it can shift the peak and modify the Model-I line shape in Fig. 10. The paper gives no check of this approximation or of the qmax = 1150 MeV cutoff dependence. The claimed distinction between molecular and Breit-Wigner N(1535) is therefore conditional on a test that is not reported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the three-body weak decay Λ_c^+ → n π^+ π^0 in a model where the N(1535) resonance is dynamically generated from S-wave coupled-channel meson-baryon interactions (πN, ηN, KΛ, KΣ) within the chiral unitary approach. The authors derive the production amplitudes for the two charge states, compute the π^0 n and π^+ n invariant mass distributions, and find a distinct peak near 1500 MeV in both. They then compare the π^0 n line shape with a conventional Breit-Wigner description of the N(1535), obtain a narrower shape in the molecular picture, and recommend experimental measurements to distinguish the two scenarios.","tokens_in":11918,"tokens_out":8831,"duration_ms":78129,"significance":"If the line-shape prediction survives closer scrutiny, the paper would provide a new, experimentally accessible observable for discriminating between molecular and quark-core interpretations of the N(1535). The technical framework is standard and the hadronization procedure from the quark-level weak decay is clearly presented. The paper also makes a serious effort to examine the dependence of the results on the free parameters α and C. However, the central claim rests on two approximations that are not quantitatively controlled in the current manuscript: the truncation of the meson-baryon loop function to its imaginary part, and an ad hoc high-energy suppression factor that is applied asymmetrically to the two models being compared.","major_comments":[{"comment":"The approximation of keeping only the imaginary part of the loop function is not justified at the energies considered. The text states that the πN-channel threshold is 'far away' from the N(1535) region, but at √s around 1500 MeV the πN channel (threshold about 1078 MeV) is open, and the ηN threshold (about 1486 MeV) lies inside the plotted range. The real part of G(s) from Eq. (11) is therefore not negligible. This matters because the production amplitude in Eq. (9) is formally equivalent to h[1−VG]^{-1} only when the same full G is used in the production loop and in the scattering equation; replacing G by i Im(G) changes both the phase and magnitude of the rescattering contribution, which can shift the peak position and modify the line shape. No numerical check of this approximation, nor any sensitivity study with respect to the cutoff qmax = 1150 MeV, is reported. This is a load-bearing step for the claimed line-shape distinction, and it must be addressed with a quantitative comparison.","section":"Section II, after Eq. (15)"},{"comment":"The comparison between Model I (chiral unitary) and Model II (Breit-Wigner) is not made on equal footing. The smooth suppression factor of Eq. (28) with Mcut = 1650 MeV is applied to the chiral-unitary amplitude, but no equivalent suppression is applied to the Breit-Wigner amplitude. Because this factor cuts off the high-invariant-mass tail of Model I, the statement that the Model II resonance shape is broader than Model I may be partly an artifact of the asymmetric treatment. The authors should either apply the same smooth suppression to Model II or show the comparison with both models computed without the suppression factor, so that the difference in line shape is a genuine model prediction.","section":"Section III, Fig. 10 and Eq. (28)"},{"comment":"The three-momentum cutoff qmax = 1150 MeV is adopted from previous fits to N(1535) properties, but the dependence of the predicted invariant mass distributions on this parameter is not tested. Since the loop function enters both the scattering T-matrix in Eq. (16) and the production amplitudes in Eq. (9), a moderate change in qmax could affect the peak position and especially the line-shape comparison with the Breit-Wigner model. A numerical sensitivity study (or a comparison with a dimensional-regularization scheme for G) is needed to establish that the central claim is robust against the choice of regulator.","section":"Section II, Eq. (11); Section III, Figs. 6-10"}],"minor_comments":[{"comment":"The notation 'N+(1535)' should be 'N(1535)+' for consistency with the use of 'N(1535)^0' and 'N(1535)^+' elsewhere.","section":"Section III, text near Fig. 6"},{"comment":"Please verify the phase-space normalization in the expression for d^2Γ/dM_{π+n}^2 dM_{π0n}^2; the factor '4 M_Λ_c+ M_n' in the numerator looks unusual and could be a typographical artifact of the rendering.","section":"Eq. (25)"},{"comment":"The assumption A = B for the weak vertex Λ_c^+ → π^+ N(1535) is ad hoc and directly affects the shape of the Breit-Wigner distribution in Model II. The authors should explicitly identify this as a source of uncertainty and, ideally, show how the line shape changes when A/B is varied.","section":"Section II, Eq. (19)"},{"comment":"The choice Mcut = 1650 MeV is close to the peak region, and the suppression factor modifies the high-energy tail of the distribution. The paper tests the sensitivity to α but not to Mcut; a brief discussion or a scan over Mcut would make the treatment of this ad hoc factor more transparent.","section":"Section II, Eq. (28)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a standard application of the chiral unitary approach to a weak three-body decay, and it proposes a testable line-shape distinction. The main weakness is that the line-shape claim currently depends on two uncontrolled approximations: the truncation of the loop function to its imaginary part and an asymmetric smooth-suppression factor in the comparison with a Breit-Wigner model. I believe these can be fixed within a revision by repeating the calculation with the full loop function, performing a cutoff-sensitivity study, and applying the same suppression to both models. The novelty is moderate but sufficient for a specialist journal if the central claim is made robust."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe one thing to know: this is a sensible, standard chiral-unitary calculation of Lambda_c -> n pi+ pi0, with a genuinely new comparison between a dynamically generated and a Breit-Wigner N(1535). The idea is good and the channel is a plausible discriminator. But the line-shape claim in Fig. 10 relies on dropping the real part of the meson-baryon loop in the production amplitude, and the stated justification for that is wrong. Treat the model comparison as promising but not yet established.\n\nWhat is new: the specific channel and the Model I vs Model II comparison appear not to have been computed before. The paper also does well in keeping the machinery transparent, checking dependence on alpha and C, and being explicit about arbitrary normalization. The peak near 1500 MeV is robust because the T-matrix is engineered to produce the N(1535) pole; that is not a weakness by itself, but it means the existence of the peak is not the test.\n\nThe main soft spot is after Eq. (15). The authors say the piN threshold is far from the N(1535) region, so only Im G is kept. At sqrt(s) ~ 1500 MeV the piN threshold (about 1078 MeV) is open, and the real part of G is not negligible. Eq. (9) has h + h G T; if G in the production loop is replaced by i Im G while T is built from the full G in Eq. (16), the resummation identity h[1 - V G]^{-1} no longer holds, and the phase and magnitude of the rescattering term change. There is no sensitivity test to qmax = 1150 MeV or to including Re G, so the shape difference in Fig. 10 is conditional. This is fixable but important. Minor caveats: the BW comparison uses A = B and a phenomenological smooth cutoff at 1650 MeV; these are acceptable as modeling choices but add uncertainty.\n\nWho it is for: hadron phenomenologists and experimentalists studying N* from charmed-baryon decays. I would send it to a serious referee; it is not a desk reject. But I would not cite the line-shape prediction until the loop approximation is addressed.\n\nRecommendation: peer review, with the request that the authors repeat the calculation keeping the full G and scanning qmax, and report whether the Model I vs Model II line-shape difference survives.","headline":"A useful new channel for N(1535) spectroscopy, but the Model I vs Breit-Wigner line-shape comparison rests on an unjustified 'keep only Im G' approximation; peer review should ask for sensitivity tests.","tokens_in":12418,"tokens_out":3707,"would_cite":false,"duration_ms":35994,"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":"In the decay $\\Lambda_c^+ \\to n\\pi^+\\pi^0$, the $N(1535)$ appears as a clear peak near 1500 MeV in the $\\pi^0n$ and $\\pi^+n$ spectra, with a line shape distinct from Breit-Wigner.","keywords":["N(1535) resonance","chiral unitary approach","dynamically generated resonance","meson-baryon molecular state","Lambda_c+ nonleptonic decay","invariant mass distribution","Bethe-Salpeter equation","coupled channels"],"falsifier":"Measure the $\\pi^0 n$ and $\\pi^+ n$ invariant mass distributions of $\\Lambda_c^+\\to n\\pi^+\\pi^0$ with high statistics. A clear peak near 1500 MeV with the predicted narrower chiral-unitary shape supports the claim; the absence of a peak, or a peak as broad as the Breit-Wigner model, would refute it.","tokens_in":11408,"feed_emoji":"⚛️","tokens_out":12231,"duration_ms":98775,"temperature":0.7,"pith_summary":"The paper predicts that the charmed-baryon decay $\\Lambda_c^+ \\to n\\pi^+\\pi^0$ carries a visible imprint of the $N(1535)$ resonance in both final-state invariant mass distributions. It describes $N(1535)$ not as an elementary quark-model state but as dynamically generated by $S$-wave interactions among the $\\pi N$, $\\eta N$, $K\\Lambda$, and $K\\Sigma$ channels in the $I=1/2$ sector. The calculation produces a clear peak near 1500 MeV in the $\\pi^0 n$ and $\\pi^+ n$ mass spectra, close to the experimental pole position of $N(1535)$, and shows that the chiral-unitary line shape is narrower than the Breit-Wigner shape obtained with an explicit resonance field. This matters because the internal structure of $N(1535)$, including the mass-inversion puzzle with $N(1440)$, remains unresolved, and a single decay could help distinguish the molecular and quark-core pictures.","feed_headline":"A charmed-baryon decay exposes the N(1535) near 1500 MeV","feed_subtitle":"The predicted π0n and π+n shapes would separate the molecular and quark-core N(1535) pictures.","key_machinery":"The machinery is the chiral unitary coupled-channel amplitude $T=[1-VG]^{-1}V$ in the $I=1/2$ sector, together with the SU(3) hadronization prescriptions that convert the weak $\\Lambda_c^+$ decay into a $\\pi$ plus a meson-baryon pair. The transition potential $V$ uses SU(3) coefficients $C_{ij}$ and is supplemented by the meson-baryon loop function $G$, computed with a three-momentum cutoff $q_{\\max}=1150\\ \\mathrm{MeV}$ and approximated by its imaginary part in the $N(1535)$ region. This object turns the weak-production vertex into resonant line shapes: the pole of $[1-VG]^{-1}$ in the $\\pi N$, $\\eta N$, $K\\Lambda$, $K\\Sigma$ coupled system is the $N(1535)$, and the same $G$ appears in the rescattering diagrams that build the final-state peaks.","core_discovery":"The central claim is that in $\\Lambda_c^+ \\to n\\pi^+\\pi^0$ the $\\pi^0 n$ and $\\pi^+ n$ invariant mass distributions each show a peak near 1500 MeV, produced by the dynamically generated $N(1535)^0$ and $N(1535)^+$, respectively. After weak hadronization of the charm quark, the final meson-baryon pairs scatter through the coupled channels $\\pi N$, $\\eta N$, $K\\Lambda$, and $K\\Sigma$, and the Bethe-Salpeter amplitude $T=[1-VG]^{-1}V$ generates the resonance pole without inserting an explicit $N(1535)$ field. When the same observable is computed with a Breit-Wigner amplitude for an elementary $N(1535)$, the line shape is broader; the difference is the paper's testable signature. The conclusion is that $\\Lambda_c^+ \\to n\\pi^+\\pi^0$ can serve as a probe of $N(1535)$'s nature, and the paper calls for precise measurements of this decay.","pith_inferences":["Extending the calculation to the full complex loop function, rather than only its imaginary part, would test whether the predicted 1500 MeV peak position is stable; a large shift would make the imaginary-part approximation the dominant systematic.","The same hadronization-plus-rescattering construction could be applied to sibling decays, for example $\\Lambda_c^+\\to p\\pi^-\\pi^+$, to see whether the $N(1535)$ line shape is universal or channel dependent.","If future data show no peak near 1500 MeV, the molecular production mechanism assumed here would be disfavored for these weak decays, unless the weak vertex is more complicated than the constant strength $V_p$ used in the paper."],"forward_implications":["The $\\pi^0 n$ and $\\pi^+ n$ invariant mass distributions in $\\Lambda_c^+\\to n\\pi^+\\pi^0$ will each show a peak near 1500 MeV, marking the $N(1535)^0$ and $N(1535)^+$ poles.","The chiral-unitary line shape is narrower than the Breit-Wigner shape from an explicit $N(1535)$ field, so the same observable can separate the molecular and quark-core descriptions.","Varying the suppression parameter $\\alpha$ or the color factor $C$ changes the signal strength but leaves the peak position essentially unchanged.","A single Dalitz-plot analysis of this decay can access both charged and neutral $N(1535)$ states and constrain the relative weight $C$ of the two weak-emission mechanisms."],"supporting_citations":[{"why":"Establishes the charmed-baryon decay with final-state rescattering and supplies the Breit-Wigner parameters used for the comparison model.","marker":"[10]"},{"why":"Provides the meson-baryon loop function G and the cutoff regularization employed for its analytic expression.","marker":"[32]"},{"why":"Supplies the three-momentum cutoff qmax = 1150 MeV that puts the N(1535) pole in the right place.","marker":"[46, 49]"},{"why":"Justifies keeping only the imaginary part of the loop function in the N(1535) energy region.","marker":"[55, 56]"},{"why":"Gives the transition potential V used inside the Bethe-Salpeter equation.","marker":"[57, 58]"},{"why":"Supplies the SU(3) coefficients C_ij for the coupled-channel scattering matrix.","marker":"[59]"},{"why":"Motivates the smooth high-energy suppression factor applied above the 1650 MeV cutoff.","marker":"[64, 65]"}],"fun_headline_variants":["Probing N(1535) via Lambda_c+ decay","N(1535) revealed in charm-baryon decay","Charmed decay sharpens N(1535) picture","Molecular N(1535) leaves distinct decay signature","New decay channel to test N(1535) nature"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that in the $N(1535)$ region the meson-baryon loop function can be reduced to its imaginary part while the three-momentum cutoff is fixed at $q_{\\max}=1150\\ \\mathrm{MeV}$; if the real part matters or the cutoff is wrong, the predicted peak and line shape shift.","fun_headline_variants_meta":{"raw":{"variants":["Probing N(1535) via Lambda_c+ decay","N(1535) revealed in charm-baryon decay","Charmed decay sharpens N(1535) picture","Molecular N(1535) leaves distinct decay signature","New decay channel to test N(1535) nature"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00082,"raw_usage":{"total_tokens":3581,"prompt_tokens":930,"completion_tokens":2651,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":546,"completion_tokens_details":{"reasoning_tokens":2569}},"tokens_in":546,"tokens_out":2651,"duration_ms":17929,"temperature":1.0,"reasoning_tokens":2569,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:44:36.773441+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the $\\pi^0 n$ and $\\pi^+ n$ invariant mass distributions of $\\Lambda_c^+\\to n\\pi^+\\pi^0$ with high statistics. A clear peak near 1500 MeV with the predicted narrower chiral-unitary shape supports the claim; the absence of a peak, or a peak as broad as the Breit-Wigner model, would refute it.","supporting_citations":[],"review_version":1}