{"id":"acd4173b-0080-4fd7-8d67-1bdd81a2e42d","arxiv_id":"2411.18134","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A quark-level hadronization plus meson rescattering model reproduces why D_s -> pi+ pi+ pi- contains only the f0(980) resonance while D+ -> pi+ pi+ pi- contains both f0(500) and f0(980).","lead":"This paper explains why the D_s and D meson decays into three pions show different sets of light scalar resonances. It models each decay as a weak production of meson pairs followed by rescattering, and fits the resulting amplitudes to BaBar, BESIII, and LHCb data.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The D_s+ no-f0(500) prediction rests on the exact ηη cancellation in Eq. (5), which follows only from the unstated equal-weight assumption for s-sbar pair creation in Eq. (1); a modest strangeness suppression would introduce a π+ηη source that can rescatter into ππ.","rationale":"The reader's weakest_assumption is correct, and I do not have a different primary objection. The reason this is the load-bearing point rather than the fits is that the D_s vs D+ distinction is made before any fitting: Eq. (7) is the only input that selects KKbar sources for D_s and five channels for D+. Normalizations and µ only adjust line shapes, not the presence or absence of f0(500). The cancellation in Eq. (5) is algebraic and depends on the SU(3)-symmetric pair-creation ansatz; the paper cites earlier work using the same ansatz but gives no independent justification or robustness study. The proposed λ scan directly tests whether the central qualitative claim survives a plausible breaking of that symmetry. If it survives, the paper's conclusion is strengthened; if not, the conclusion is an artifact of an unexamined assumption. I therefore keep the reader's CONDITIONAL verdict rather than moving to accept or reject. Agreement is 'agree' because the reader named the same assumption as weakest; my contribution is to make the test quantitative and to note the low-energy phase problem as the place a spurious f0(500) would appear.","tokens_in":18493,"tokens_out":20076,"duration_ms":194042,"concrete_test":"Modify the hadronization operator in Eq. (1) to (bar uu + bar dd + λ bar ss) and recompute the D_s+ amplitude: Eq. (7) becomes (1+β)π+(K+K− + K0Kbar0) + (2/3)(λ−1)π+ηη, so add δ Gηη Tηη→π+π− to Eq. (9) with δ = (2/3)(λ−1)/(1+β). Refit C (and optionally µ) to the BABAR/BESIII/LHCb D_s+ data below 1.2 GeV for λ = 0.5, 0.75, 1.25, 1.5 and compare the resulting magnitude and phase of the π+π− S-wave with the λ=1 baseline. If the curves shift by more than the experimental uncertainties in the √s<0.8 GeV region, or if a broad f0(500)-like structure appears, the no-f0(500) claim is not robust to the hadronization assumption; if they are stable, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central qualitative result—that D_s+ → π+π+π− contains f0(980) but no f0(500)—is derived from Eq. (7), in which only K+K− and K0Kbar0 survive the primary hadronization. That survival is not generic: the +(2/3)π+ηη term from (M·M)_33 in Eq. (5) is cancelled by the −(2/3)π+ηη term from η(M·M)_12 only because the vacuum insertion in Eq. (1) creates ubar-u, dbar-d, and sbar-s with exactly equal weights. If the insertion is replaced by (bar uu + bar dd + λ bar ss), the K+K− and K0Kbar0 terms are unchanged, but the net π+ηη source becomes (2/3)(λ−1)π+ηη from the external-W piece, and Eq. (9) must include a Gηη Tηη→π+π− contribution with relative strength (2/3)(λ−1)/(1+β). Within the same ChUA, Tηη→π+π− is not zero, so the question is quantitative: for moderate λ≠1, can this term produce a visible low-mass enhancement that would look like f0(500)? The paper neither justifies λ=1 independently nor checks the sensitivity. This matters because the Sec. III fit admits the model already underestimates the low-energy D_s data and describes the phase poorly below 0.8 GeV, exactly the region where f0(500) would appear. The equal-weight ansatz is the step that creates the qualitative D_s vs D+ difference, so it is load-bearing.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the three-pion decays D_s+ -> pi+ pi+ pi- and D+ -> pi+ pi+ pi- in a model that combines quark-level weak hadronization with final-state rescattering in the chiral unitary approach (ChUA). Starting from external and internal W-emission mechanisms, the authors derive the primary two-meson channels produced in each decay. For D_s+ they find in Eq. (7) that only K+K- and K0Kbar0 channels survive in S-wave, because an exact cancellation removes the pi+ eta eta term; consequently only f0(980), which couples to K Kbar, is generated. For D+ the five coupled channels pi+ pi-, pi0 pi0, K+K-, K0Kbar0, and eta eta are produced, so both f0(500) and f0(980) are generated. The model is fitted to BaBar/BESIII/LHCb measurements of the pi+ pi- S-wave magnitude and phase, with two renormalization schemes and an alternative fit without the eta eta channel. The authors conclude that the observed pattern supports the molecular/dynamically-generated interpretation of these scalar mesons.","tokens_in":18950,"tokens_out":16909,"duration_ms":156726,"significance":"If correct, the paper gives a very simple selection mechanism for a striking experimental pattern: the same three-pion final state is f0(980)-dominated when produced from D_s+, but contains both f0(500) and f0(980) when produced from D+. The derivation leading to Eq. (7) is transparent, and the phase prediction in Fit I for D_s+ is parameter-free. The paper is also commendably explicit about the low-energy discrepancies and about the alternative eta-eta-excluded fit. The significance is moderate rather than high: the conclusion supports a picture that is already widely explored in ChUA, and the robustness of the central rule depends on the equal-weight vacuum-pair assumption discussed below, which needs to be quantified before the selection rule can be regarded as a firm prediction.","major_comments":[{"comment":"The central prediction that only K+K- and K0Kbar0 are produced in the D_s+ decay follows from the exact cancellation of the pi+ eta eta term in Eq. (5). This cancellation uses the specific vacuum insertion written as ubar-u + dbar-d + sbar-s in Eq. (1) with equal weights for all three flavors. If the strange-pair insertion is weighted by lambda, the residual source is (2/3)(lambda-1) pi+ eta eta, which, through the nonzero ChUA amplitude T_{eta eta -> pi+ pi-}, feeds the pi+ pi- S-wave at low energies and can mimic f0(500). The paper neither justifies lambda = 1 independently nor studies the sensitivity to lambda. Since the D_s+ fit already underestimates the data below 0.8 GeV and describes the phase poorly in the 0.4-0.8 GeV range (Sec. III, Fig. 6), exactly the f0(500) region, this is a load-bearing gap for the qualitative D_s+ versus D+ difference.","section":"Sec. II A, Eqs. (1), (5), (7)"},{"comment":"The D_s+ fit combines data from BaBar, BESIII, and LHCb with a single global normalization C. These experiments use different amplitude conventions, with LHCb amplitude analyses typically reporting arbitrary-normalized amplitudes, so the authors must specify how the data sets are normalized to a common scale. Without this information, the combined fit quality and the claimed agreement are not well defined. In addition, no chi-squared per degree of freedom is quoted for any of the fits, which makes it difficult to assess the quantitative support for the model.","section":"Sec. III, Fig. 6, Tables I and II"},{"comment":"For D+, the quantitative agreement is obtained with two fitted normalizations D1 and D2, and in scheme II also with a fitted mu. In particular, the relative weight D2/D1 between the K Kbar contribution and the direct/eta eta contribution, which controls the relative strength of f0(980) and f0(500), is a free parameter rather than a prediction from the CKM and color input. The authors should report the fitted ratio D2/D1 and compare it with the naive expectation from V'_P V_cs V_us / (V_P V_cd V_ud), or explicitly state that this ratio is not predicted by the model.","section":"Sec. III, Tables I and II, Figs. 6 and 7"},{"comment":"The improved description of the D_s+ magnitude is obtained only after removing the eta eta channel from the coupled-channel calculation. This removal is not derived within the formalism and changes the unitarity structure of the T-matrix; the paper justifies it only by citing Ref. [56] for the resulting larger f0(980) width. As presented, this is an ad hoc model alteration rather than a controlled check. The central Eq. (7) argument is unaffected, but the quantitative support for the claim that the model naturally explains the D_s+ data rests in part on this altered model; the authors should either justify the removal as a separate scheme or soften the corresponding claim.","section":"Sec. III, Tables III and IV, Fig. 9"}],"minor_comments":[{"comment":"The vertex factor V'_P is introduced without definition; please clarify whether it denotes the same weak vertex factor as V_P or a different one, since the two enter the expressions for D1 and D2.","section":"Sec. II B, Eqs. (11), (13), and (18)"},{"comment":"There is a typo in the sentence 'In the present work. we use the dimensional regularization method'; the manuscript would benefit from a careful proofread of punctuation and spacing throughout.","section":"Sec. II C"},{"comment":"The restriction of the comparison to sqrt(s) < 1.2 GeV is stated only in the numerical section; since f0(500) is a broad low-energy structure, the authors should state this restriction earlier and discuss the possible impact of the excluded higher-energy region on the conclusions.","section":"Sec. III"},{"comment":"The sentence regarding the factor of 2 cancelled by the 1/2 factor in the pi0 pi0 and eta eta propagators is hard to follow; please rephrase to state explicitly which identical-particle symmetry factor is used.","section":"Eq. (20)"},{"comment":"The assumption eta = eta_8 ignores eta-eta' mixing; because the cancellation in Eq. (5) depends on the exact flavor coefficient of eta, the authors should at least acknowledge this approximation. Also, the figures and tables in the submitted version contain apparent encoding artifacts (e.g., repeated 'uni00000013' sequences) that must be cleaned in the final version.","section":"Eq. (4) and Figs. 6-9"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely question and the mechanism behind Eq. (7) is genuinely attractive. My main concern is not internal inconsistency but robustness: the exact cancellation that produces the no-f0(500) prediction for D_s+ depends on an unquantified equal-weight assumption for strange-pair creation, and the paper's own fit shows the model is weakest precisely in the f0(500) region. The alternative fit without the eta-eta channel also needs to be presented more carefully. These issues are fixable without changing the basic approach, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real contribution is the structural hadronization argument: in D_s+ -> pi+pi+pi-, the external and internal W-emission mechanisms leave only K+K- and K0Kbar0 as primary S-wave sources (Eq. 7), so f0(500) is absent and f0(980) comes from KKbar rescattering; in D+ -> pi+pi+pi-, the same formalism gives five coupled channels (Eq. 20) and generates both f0(500) and f0(980). The derivation from the quark-level diagrams to Eq. (7) is short, transparent, and parameter-free once you accept the standard vacuum-insertion ansatz u ubar + d dbar + s sbar with equal weights. That is genuinely new, and it connects the molecular interpretation of the two scalars directly to recent LHCb, BaBar, and BESIII data. The D_s part overlaps with Ref. [52], but the D+ comparison and the contrast in hadronization structure are new.\n\nThe fits are decent, not brilliant. The D+ side describes the LHCb amplitudes well. The D_s side reproduces the f0(980) peak and the phase near 1 GeV, but the text agrees that the magnitude is too small below 0.8 GeV and the phase is poor in 0.4-0.8 GeV—the region where a f0(500)-like signal would show up. So the abstract's claim that magnitudes and phases are in agreement is an overstatement.\n\nThe soft spot the stress-test identifies is real. The disappearance of the eta eta source in Eq. (5) comes from an exact cancellation between the pi+(M.M)_33 and eta(M.M)_12 terms. With a vacuum insertion (u ubar + d dbar + lambda s sbar), a net primary pi+ eta eta term proportional to (lambda-1) survives and can rescatter into pi+ pi- via T_{eta eta -> pi+ pi-}. The paper neither justifies lambda=1 nor tests the sensitivity. The practical impact may be moderate because the eta eta threshold sits near 1.1 GeV, so this extra source mostly affects the region above the f0(500) peak; but the sign matters and the model already has a low-mass shortfall. This is a referee question, not a fatal flaw. The later switch to a ChUA without the eta eta channel improves the D_s magnitude but is not deeply motivated here, and the fits leave several normalization parameters and sometimes mu free.\n\nBottom line: the structural argument is worth publishing and deserves a serious referee. I would recommend sending it to review and asking for a lambda sensitivity check, a softened abstract, and an explicit discussion of the low-mass D_s phase. Colleagues working on light scalars and final-state interactions in charm decays will want to read it.","headline":"The structural hadronization argument distinguishing D_s and D+ is clean and new, but the exact cancellation that removes f0(500) from D_s rests on an untested equal-weight vacuum ansatz.","tokens_in":19463,"tokens_out":8055,"would_cite":true,"duration_ms":69019,"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 $D_s^+\\to\\pi^+\\pi^+\\pi^-$ decays the $f_0(500)$ is absent because the weak decay hadronizes only into $K\\bar K$ pairs, while in $D^+\\to\\pi^+\\pi^+\\pi^-$ all five coupled channels are produced and both $f_0(500)$ and $f_0(980)$ appear.","keywords":["f0(500) meson","f0(980) meson","chiral unitary approach","final-state interactions","hadronization","D_s meson three-body decays","scalar meson molecular states","Bethe-Salpeter equation"],"falsifier":"A decisive test would be a precision measurement of the $D_s^+\\to\\pi^+\\pi^+\\pi^-$ $\\pi^+\\pi^-$ S-wave amplitude in the $0.4$–$0.8$ GeV region, where the model undershoots the data: if the phase shows the rapid motion characteristic of the $f_0(500)$ pole, the cancellation in Eq. (7) is incomplete and the central explanation fails. A lattice or QCD-sum-rule evaluation of the hadronization matrix elements that finds the strange-pair weight differing from the up/down weight would also break the derivation.","tokens_in":18278,"feed_emoji":"🔬","tokens_out":7621,"duration_ms":60840,"temperature":0.7,"pith_summary":"This paper sets out to explain a sharp experimental puzzle: in the Cabibbo-favored decay $D_s^+\\to\\pi^+\\pi^+\\pi^-$ only the $f_0(980)$ resonance shows up in the $\\pi^+\\pi^-$ S-wave amplitude, while the $D^+\\to\\pi^+\\pi^+\\pi^-$ decay shows both $f_0(500)$ and $f_0(980)$. The authors show that the difference follows from the dominant external and internal $W$-emission mechanisms at the quark level together with the final-state rescattering of the hadronized meson pairs. Because the $D_s^+$ weak decay hadronizes only into $K^+K^-$ and $K^0\\bar K^0$ pairs, the $\\pi\\pi$ channel that would seed $f_0(500)$ is never produced; the $D^+$ decay, in contrast, produces all five coupled channels $(\\pi^+\\pi^-,\\pi^0\\pi^0,K^+K^-,K^0\\bar K^0,\\eta\\eta)$, so both resonances are generated. If the calculation is right, it supports the picture of $f_0(500)$ as a $\\pi\\pi$ resonance and $f_0(980)$ as a $K\\bar K$ molecular state, both dynamically generated by the chiral unitary approach.","feed_headline":"Why f0(500) hides in D_s^+ decays but appears in D^+ decays","feed_subtitle":"A quark-pair creation rule plus meson rescattering predicts exactly which scalar resonances show up in each three-pion decay.","key_machinery":"The load-bearing objects are the hadronization identities of Eqs. (1)–(7) and the coupled-channel amplitudes of the chiral unitary approach. The hadronization step replaces each quark–antiquark pair with $(M\\cdot M)_{ij}$ using the pseudoscalar matrix $P$ of Eq. (4), with the vacuum creating $\\bar u u + \\bar d d + \\bar s s$ pairs; this is what collapses the $D_s^+$ final state to the $K\\bar K$ channels. The scattering of these pairs is then resummed through the Bethe–Salpeter equation $T=[1-VG]^{-1}V$, with the $S$-wave potentials $V_{ij}$ from the lowest-order chiral Lagrangian and the two-meson loop function $G$, producing the dynamically generated resonances that are identified with $f_0(500)$ and $f_0(980)$.","core_discovery":"The central claim is that the absence of $f_0(500)$ in $D_s^+\\to\\pi^+\\pi^+\\pi^-$ is not accidental but follows from the hadronization of the weak decay. Using the quark-pair-creation ansatz $\\bar u u + \\bar d d + \\bar s s$, the authors derive the hadronic amplitudes of the external and internal $W$-emission diagrams; after the $\\eta\\eta$ pieces cancel, the $D_s^+$ transition leaves only $\\pi^+ K^+ K^-$ and $\\pi^+ K^0 \\bar K^0$ in S-wave (Eq. (7)). The final $\\pi^+\\pi^-$ pair then arises purely from rescattering $K\\bar K\\to\\pi^+\\pi^-$, which the chiral unitary approach produces through the $f_0(980)$. For $D^+$, the corresponding derivation keeps all five coupled channels and a tree-level $\\pi^+\\pi^+\\pi^-$ term (Eq. (20)), so both $f_0(500)$ and $f_0(980)$ appear. The authors fit the measured magnitudes and phases and find the $D^+$ data are described well, with the $D_s^+$ data described reasonably above the $\\pi\\pi$ threshold region.","pith_inferences":["If the equal-weight quark-pair creation is replaced by a suppression of strange pairs, the $D_s^+$ decay would acquire a direct $\\pi\\pi$ component; the size of the $f_0(500)$ signal would then be a measure of the $SU(3)$ breaking in the hadronization vacuum, a quantity that could be extracted from data or from lattice calculations.","The same mechanism should govern other $D_s^+$ decays into three pseudoscalars, predicting, for example, which scalar resonances appear in $D_s^+\\to\\pi^+K^+K^-$ versus $D^+\\to\\pi^+K^+K^-$; those channels could be compared with existing experimental results.","The absence of $f_0(500)$ in $D_s^+$ decays and its presence in $D^+$ decays could serve as a production-side filter to separate $\\pi\\pi$-driven from $K\\bar K$-driven scalar mesons in other final states."],"forward_implications":["The $D_s^+\\to\\pi^+\\pi^+\\pi^-$ amplitude should show no $f_0(500)$ signal, exactly as observed in the current measurements.","The $D^+\\to\\pi^+\\pi^+\\pi^-$ amplitude contains both the broad $f_0(500)$ bump below 1 GeV and the narrow $f_0(980)$ cusp at the $K\\bar K$ threshold.","The phase of the $\\pi^+\\pi^-$ amplitude in the $D_s^+$ decay is parameter-free once the magnitude is fitted, so an improved phase measurement in the low-mass region provides a direct check of the calculation.","Dropping the $\\eta\\eta$ channel widens the $f_0(980)$ and improves the description of the $D_s^+$ magnitude, which the model ascribes to the small width of $f_0(980)$ in the five-channel calculation.","The framework identifies the produced isoscalar channels, so $a_0(980)$ and other $I=1$ contributions are predicted to be strongly suppressed in these decays."],"supporting_citations":[{"why":"Supplies the $D^+\\to\\pi^+\\pi^+\\pi^-$ $\\pi^+\\pi^-$ S-wave magnitude and phase data that the model fits.","marker":"[47]"},{"why":"Provides the first $D_s^+\\to\\pi^+\\pi^+\\pi^-$ amplitude data used in the fit.","marker":"[48]"},{"why":"Provides a second $D_s^+$ dataset for the same amplitude.","marker":"[49]"},{"why":"Completes the $D_s^+$ comparison with a third measured dataset.","marker":"[50]"},{"why":"Earlier work interpreting $f_0(980)$ as a molecular state in $D_s^+\\to\\pi^+\\pi^+\\pi^-$, which this paper extends to both decays.","marker":"[52]"},{"why":"Defines the external and internal $W$-emission quark-level topology that underlies the hadronization amplitudes.","marker":"[54]"},{"why":"Fixes the $\\eta=\\eta_8$ flavor assignment used in the hadronization matrices, which produces the cancellation in Eq. (5).","marker":"[58]"},{"why":"Provides the coupled-channel Bethe-Salpeter formalism $T=[1-VG]^{-1}V$ that dynamically generates $f_0(500)$ and $f_0(980)$.","marker":"[60]"}],"fun_headline_variants":["Why f0(500) vanishes in D_s^+ but shows in D^+ decays","Missing f0(500) in D_s^+ decay traced to weak hadronization","D_s^+ vs D^+: why different scalar mesons appear","Quark-pair creation explains f0(500) absence in D_s^+ decays","Rescattering model reveals distinct f0 states in D decays"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument depends on the assumption that the vacuum creates up, down, and strange quark–antiquark pairs with equal weight during hadronization, since that precise equality is what cancels the direct $\\pi\\pi$ content in the $D_s^+$ decay and leaves only $K\\bar K$ channels.","fun_headline_variants_meta":{"raw":{"variants":["Why f0(500) vanishes in D_s^+ but shows in D^+ decays","Missing f0(500) in D_s^+ decay traced to weak hadronization","D_s^+ vs D^+: why different scalar mesons appear","Quark-pair creation explains f0(500) absence in D_s^+ decays","Rescattering model reveals distinct f0 states in D decays"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001259,"raw_usage":{"total_tokens":5215,"prompt_tokens":1059,"completion_tokens":4156,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":675,"completion_tokens_details":{"reasoning_tokens":4049}},"tokens_in":675,"tokens_out":4156,"duration_ms":25090,"temperature":1.0,"reasoning_tokens":4049,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:28:38.274173+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be a precision measurement of the $D_s^+\\to\\pi^+\\pi^+\\pi^-$ $\\pi^+\\pi^-$ S-wave amplitude in the $0.4$–$0.8$ GeV region, where the model undershoots the data: if the phase shows the rapid motion characteristic of the $f_0(500)$ pole, the cancellation in Eq. (7) is incomplete and the central explanation fails. A lattice or QCD-sum-rule evaluation of the hadronization matrix elements that finds the strange-pair weight differing from the up/down weight would also break the derivation.","supporting_citations":[{"cited_title":"Dalitz Plot Analysis of D_s+ -> pi+ pi- pi+","cited_arxiv_id":"0808.0971","evidence_quote":"Provides the first $D_s^+\\to\\pi^+\\pi^+\\pi^-$ amplitude data used in the fit."},{"cited_title":"Ablikim et al","cited_arxiv_id":null,"evidence_quote":"Provides a second $D_s^+$ dataset for the same amplitude."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the external and internal $W$-emission quark-level topology that underlies the hadronization amplitudes."}],"review_version":1}