{"id":"6eb7c3b4-5205-4361-813a-1342083b0d9b","arxiv_id":"2411.17098","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A simulation-based feasibility study shows that near-threshold Dbar Kbar mass distributions from Lambda_b decays could determine the D*_s0(2317) mass to about 4 MeV and its DK molecular probability to about 11%.","lead":"This paper proposes a way to learn about the D*_s0(2317) particle from the shapes of Dbar Kbar mass distributions in Lambda_b decays. It simulates how precisely future LHCb measurements could pin down the particle's binding energy and its Dbar Kbar molecular nature.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed 4 MeV pole and 11% DK-probability precisions are internal to the authors' production model and unitarization scheme; an equally plausible production amplitude or loop regulator could shift these beyond the quoted errors.","rationale":"The paper is an honest feasibility study: it clearly states that pseudodata are taken from the local hidden gauge approach, and it acknowledges that the method's usefulness depends on the data being measured. The reader's CONDITIONAL verdict already captures the central risk that real LHCb data may differ from the assumed model. My concern sharpens this: even granting the pseudodata are measured with 5% errors and no backgrounds, the quoted precision for the pole and, especially, for the compositeness is not a property of the data alone but of the specific unitarized coupled-channel scheme and of the assumption that the production vertex is a constant. Compositeness is not a direct observable; Eq. (14) gives a scheme-dependent quantity, and the 42 MeV binding means the pole is far from the data region, so extrapolation depends on the regulator and the energy dependence of V. The proposed test would quantify this by using a different production shape, which is exactly the kind of variation that can be expected in a real weak decay. I therefore do not change the reader's verdict, because the conditionality is already present, but I would make the condition explicit: the 4 MeV and 11% figures require the production mechanism and unitarization scheme assumed in the paper to hold, and should not be quoted as standalone experimental projections without a systematic-model study.","tokens_in":13058,"tokens_out":6594,"duration_ms":67972,"concrete_test":"Generate new pseudodata in which the tree-level amplitude A in Eqs. (18)-(19) is multiplied by a smooth factor (1 + c (M_inv - m_th)/Lambda) with c = +/-1 and Lambda = 1 GeV, keeping the same underlying T-matrix, and rerun the resampling extraction. If the reconstructed pole mass or DK probability shifts by more than the quoted 4 MeV or 0.08, the stated precisions are not robust to production-model systematic uncertainty; if the shifts are smaller, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical claims (pole mass 2319.3 ± 3.9 MeV, DK probability 0.70 ± 0.08) are obtained by fitting roughly 60 pseudodata points generated from the authors' local hidden gauge model, with only 50 resampling fits and 5% assumed errors. The quoted uncertainties therefore propagate only the statistical resampling spread within one assumed model class. They do not include the leading systematic uncertainty for the observables in question. The production amplitude in Eqs. (18)-(19) is taken as a single energy-independent constant A with equal weights for the two channels, and the compositeness extracted via Eq. (14) is known to depend on the loop-function regulator and on the energy-dependent terms added in Eqs. (24)-(26). The data window lies entirely above threshold while the pole sits 42 MeV below it; the mapping from threshold data to the pole and to the probability rests on the assumed analytic form of V(s) and on the regulator. A different smooth production form factor or a different unitarization scheme could change the extracted pole by more than 4 MeV and P1+P2 by more than 0.08, without degrading the fit quality in the observed window. The paper does not test this, so the headline precision is conditional on the production and interaction models being the true ones.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a method to extract properties of the D*_{s0}(2317) from future measurements of the \\bar D^0 K^- and D^- \\bar K^0 invariant mass distributions in \\Lambda_b \\to \\Lambda_c (\\bar D \\bar K)^- decays. The authors first compute these distributions using a unitarized coupled-channel formalism (channels \\bar D^0 K^-, D^- \\bar K^0, D_s^- \\eta) based on the local hidden gauge approach, finding a strong threshold enhancement due to the subthreshold pole. They then generate pseudodata from this model with 5% relative Gaussian errors, fit them with a more general energy-dependent potential within a 50 MeV window above threshold, and use the resampling method to estimate uncertainties. They report that this procedure recovers the scattering lengths, an I=0 bound state at 2319.3 \\pm 3.9 MeV, and a \\bar D \\bar K molecular probability P1+P2 = 0.70 \\pm 0.08, and conclude that LHCb data could determine the nature of the D*_{s0}(2317) with better precision than correlation-function analyses.","tokens_in":13362,"tokens_out":7100,"duration_ms":68691,"significance":"If the claimed precision survived contact with real data, the paper would provide a valuable and timely method for extracting hadronic observables from decay mass distributions, improving on the much larger uncertainties quoted for femtoscopic correlation functions. The internal consistency of the inverse problem is a genuine strength: the fit potential in Eqs. (23)-(26) is more general than the local hidden gauge potential that generated the pseudodata, and the pole position is an output rather than a fitted input. The resampling procedure is also appropriate for the correlated parameter sets. However, the numerical claims are established only for one generating model, and the paper does not quantify the leading systematics (production amplitude, regulator scheme, possible nonmolecular component). Since the headline result is a precision claim, this model dependence is the central issue.","major_comments":[{"comment":"The quoted uncertainties of \\pm 3.9 MeV on the pole and \\pm 0.08 on P1+P2 are resampling dispersions for pseudodata generated by the local hidden gauge model. The abstract and conclusions state these as the precision obtainable from the data, but the data window lies entirely above threshold and the extrapolation to a pole 42 MeV below rests on the assumed unitarized form T=[1-VG]^{-1}V and the loop function in Eq. (5). A different but equally plausible unitarization or regulator could shift the pole and probability by more than the quoted errors while reproducing the in-window distributions. The authors should perform a closure test using pseudodata generated from an alternative model (e.g., a different loop regulator or a potential with an explicit genuine-state term) and report the resulting shifts in the extracted observables. The paper itself in Sec. II C acknowledges that the fit freedom could in principle render the method useless, but it does not address this particular source of systematic error.","section":"Section III, Tables IV and VI"},{"comment":"The production amplitude is taken as a single energy-independent constant A with equal weights for \\bar D^0 K^- and D^- \\bar K^0. Real LHCb data will contain momentum-dependent production vertices and possibly unequal weights; because the unitarized amplitude is linear in A, such effects can be absorbed into the fitted potential parameters and bias the extracted pole and compositeness. The authors should test the inversion with A_1 \\neq A_2 or with a smooth form factor in the production vertex, and show how the extracted pole position and P1+P2 change.","section":"Section II B, Eqs. (18)-(19)"},{"comment":"The compositeness extracted from Eq. (14) depends on the derivative of the loop function G with respect to s, not only on the on-shell amplitude. Although qmax is a fitted parameter, the functional form of G in Eq. (5) is fixed; a different regularization (for instance dimensional regularization, as commonly used in chiral unitary approaches) would give a different dG/ds and hence different P_i for the same in-window data. The paper should quantify this sensitivity before presenting 11% uncertainty on P1+P2.","section":"Section II C, Eq. (14)"},{"comment":"The paper does not report any goodness-of-fit measure or the number of converged fits in the resampling procedure. With eight parameters and sixty data points, the stability of the quoted dispersions should be demonstrated, for example by showing the distribution of \\chi^2 and checking that the results are stable when the number of resampled fits is increased beyond 50.","section":"Section III, resampling procedure"}],"minor_comments":[{"comment":"The phase-space comparison mentioned in the text is not identified in the caption; please specify how the phase-space curve is normalized and whether it includes the same kinematic prefactors as the full distributions.","section":"Section III, Fig. 3"},{"comment":"The symbols \\sum\\sum in Eq. (20) are not defined; please state the spin sums and whether a spin average is included.","section":"Section II B, Eq. (20)"},{"comment":"The entry r_{0,3} is effectively undetermined; please mark it explicitly in the table (e.g., as unconstrained) rather than only noting this in the text.","section":"Section III, Table V"},{"comment":"The phrase \"model independent analysis\" is too strong; the analysis still fixes the unitarization form and the isospin structure of the potential. A more precise term would be \"minimal model\" analysis.","section":"Section II C"}],"recommendation":"major_revision","confidential_remarks":"The paper is a feasibility study based on pseudodata, and the authors are transparent about this in the body. The main issue is that the headline precision is presented without the caveat that it is conditional on the generating model. I think the paper can be made publishable if the authors add closure tests with alternative production amplitudes and loop regulators, and if the abstract is reworded to state the conditional nature of the precision claims. There are no concerns about citation patterns or novelty disclosure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a feasibility study for extracting D*_s0(2317) properties from Dbar Kbar mass distributions in Lambda_b -> Lambda_c Dbar Kbar decays. What's new: the forward calculation of these specific distributions using the hidden gauge approach, and the inverse test showing that a flexible energy-dependent potential fit to threshold pseudodata can recover the pole to about 4 MeV and the DK molecular probability to about 11% via resampling. That is a useful demonstration, and the authors are honest that the method follows Ref. [34] for correlation functions.\n\nThe paper does several things well. The inverse problem is not circular: the generating model is the local hidden gauge potential, while the fit uses a more general potential with free parameters including energy dependence, and the pole is an output, not a fit parameter. The resampling approach suits correlated parameters, and the authors correctly note that parameter uncertainties can be large while observables are more stable. Including the D_s eta channel and energy-dependent terms allows for non-molecular components.\n\nThe soft spots are real. The precision claims are within-model: the 'data' are pseudodata generated by the authors' own model, with an assumed 5% error and a 50 MeV window. The production amplitude is a single constant A with equal weights for the two channels. Real LHCb data will include production form factors, backgrounds, and efficiencies that could shift the distributions in ways not captured by A. The stress-test concern is on point: a different production form factor or a different unitarization/regularization scheme could move the pole and probability beyond the quoted errors. Fitting qmax as a free parameter helps for the cutoff, but the functional form of the loop function is still fixed, and the energy-dependent terms in V(s) are a specific ansatz. The paper does not test alternative assumptions, and it should. Generating pseudodata with a different production amplitude or loop regulator and seeing if the extraction still works would make the claims much more robust.\n\nAlso, the abstract overstates the effective-range results: the paper's own numbers give r0,1 with a 36% error, r0,2 with huge uncertainty, and r0,3 undetermined. 'Precise values of the scattering lengths and effective ranges' is not supported for the effective ranges. Calling the analysis 'model independent' is a stretch even with the I=0 caveat.\n\nThat said, the central conclusion - that threshold data contain enough information to pin the bound-state pole and molecular probability within the model - appears to hold. The 4 MeV and 11% numbers are honest within the assumptions, but conditional on them. As a feasibility study, it is a reasonable contribution.\n\nWho is this for? People working on D_s0*(2317) compositeness and on extracting scattering information from production data. It deserves a serious referee, but the authors should correct the abstract, add a model-variation test, and ideally release code or data. My recommendation: send to review with a request for major revision; I would not cite it in my own work until the robustness checks are done.","headline":"Useful within-model feasibility study showing threshold DK mass distributions could pin the D*_s0(2317) pole to ~4 MeV, but the claimed precision is conditional on the production model and the abstract overstates the effective-range results.","tokens_in":13867,"tokens_out":5094,"would_cite":false,"duration_ms":41131,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper argues that the near-threshold $\\bar{D}^0K^-$ and $D^-\\bar{K}^0$ mass distributions in $\\Lambda_b \\to \\Lambda_c \\bar{D}^0K^-$ and $\\Lambda_b \\to \\Lambda_c D^-\\bar{K}^0$ decays can determine the $D^*_{s0}(2317)$ pole mass to…","keywords":["D*_s0(2317)","DK molecular state","Lambda_b decays","invariant mass distribution","coupled-channel unitarization","local hidden gauge approach","resampling method","compositeness"],"falsifier":"Measure the $\\bar{D}^0K^-$ and $D^-\\bar{K}^0$ invariant mass distributions in the observed $\\Lambda_b \\to \\Lambda_c \\bar{D}^0K^-$ decay and the companion $\\Lambda_b \\to \\Lambda_c D^- \\bar{K}^0$ decay: if the sharp near-threshold enhancement over phase space is absent, or if fitting the inverse procedure to data generated by a different production model shifts the recovered pole by much more than 4 MeV, the central precision claim is falsified.","tokens_in":12851,"feed_emoji":"⚛️","tokens_out":11223,"duration_ms":85817,"temperature":0.7,"pith_summary":"The paper sets out to show that the invariant mass distributions of $\\bar{D}^0K^-$ and $D^-\\bar{K}^0$ close to threshold in the decays $\\Lambda_b \\to \\Lambda_c \\bar{D}^0 K^-$ and $\\Lambda_b \\to \\Lambda_c D^- \\bar{K}^0$ carry enough information to pin down the properties of the $D^*_{s0}(2317)$ resonance. The authors compute these distributions from a coupled-channel model in which the resonance is a $\\bar{D}\\bar{K}$ bound state, finding a strong enhancement relative to phase space at threshold. They then solve the inverse problem: assuming the distributions are measured with 5% relative errors, they fit a flexible potential to the pseudodata using a resampling method. The recovered $I=0$ bound-state pole lies at the right mass with an uncertainty of about 4 MeV, and the summed $\\bar{D}\\bar{K}$ molecular probability comes out with an uncertainty of about 11%. If the actual experiment delivers data of this quality, the method would give a tighter, experiment-driven determination of the binding and composition of the $D^*_{s0}(2317)$ than current correlation-function analyses.","feed_headline":"D*_s0(2317) mass pinned to 4 MeV from decay shapes","feed_subtitle":"Near-threshold Lambda_b decay data could fix the state's binding and its DK molecule probability to roughly 11 percent.","key_machinery":"The machinery is a coupled-channel unitarized amplitude $T = [1 - VG]^{-1}V$, with a vector-meson-exchange potential $V$ from the local hidden gauge approach, diagonal loop functions $G_i$ regulated by a cutoff $q_{\\rm max}$ tuned to reproduce the $D^*_{s0}(2317)$ mass, and decay amplitudes built from a tree-level hadronization vertex followed by rescattering through the coupled channels. For the inverse problem, the paper replaces the model potential by a general isospin-symmetric matrix with constant and linearly energy-dependent entries, where the $\\alpha,\\beta,\\gamma$ terms are meant to absorb possible genuine nonmolecular components. The resampling method—generating many Gaussian-perturbed copies of the pseudodata and refitting each—propagates data uncertainties into derived observables: scattering lengths and effective ranges from the effective-range expansion, the pole position of the amplitude, and channel probabilities $P_i = -g_i^2\\, \\partial G_i/\\partial s$ evaluated at the pole.","core_discovery":"On the paper's own terms, the central claim is that the $D^*_{s0}(2317)$ is coupled strongly enough to the $\\bar{D}\\bar{K}$ channels to leave a sharp, measurable fingerprint in the $\\bar{D}^0K^-$ and $D^-\\bar{K}^0$ mass distributions, even though the resonance sits about 42 MeV below the $\\bar{D}\\bar{K}$ threshold and cannot decay into it. Using the local hidden gauge approach with three coupled channels ($\\bar{D}^0K^-$, $D^-\\bar{K}^0$, $D^-_s\\eta$), the authors build the decay amplitudes from a tree-level hadronization vertex followed by rescattering, and obtain distributions that rise steeply at threshold. Treating these model distributions as pseudodata with 5% relative errors, they fit a general energy-dependent potential with free parameters; despite strong correlations among the parameters, the physical observables come out stable. Averaging over resampled data sets yields a pole at $2319.3 \\pm 3.9$ MeV, scattering lengths for the two $\\bar{D}\\bar{K}$ channels with uncertainties around 14–27%, and a summed $\\bar{D}\\bar{K}$ probability $P_1+P_2 \\approx 0.70$ with about 11% error.","pith_inferences":["A direct experimental test follows: the $\\bar{D}^0K^-$ and $D^-\\bar{K}^0$ spectra from the observed decay should show a pronounced threshold peak; if it is absent or far weaker than the model prediction, the predominantly molecular picture of the $D^*_{s0}(2317)$ would be in doubt.","The same inverse-problem strategy could be transferred to other putatively molecular states near heavy-flavor thresholds, such as the $T_{cc}(3875)$ or $X(3872)$, where production mass distributions might give comparable compositeness precision.","The authors find that the energy-dependent parameters $\\alpha,\\beta,\\gamma$ are poorly determined, which suggests the data alone cannot cleanly separate a genuine nonmolecular component from a purely dynamical bound state; the claimed model independence applies to the extracted observables, not to the underlying production mechanism.","A natural extension would be to fold a more realistic error model—backgrounds, normalization, detector resolution—into the resampling, since the quoted 4 MeV and 11% uncertainties assume 5% point-to-point statistical errors."],"forward_implications":["If the $\\bar{D}^0K^-$ and $D^-\\bar{K}^0$ mass distributions are measured with roughly 5% differential-width errors over about 50 MeV above threshold, the $D^*_{s0}(2317)$ pole mass is recovered to about 4 MeV even though the pole lies about 42 MeV below threshold.","The same fit determines the summed $\\bar{D}\\bar{K}$ molecular probability of the state to about 11%, a substantial improvement over the roughly 60% uncertainty quoted for correlation-function analyses.","The data fix the $\\bar{D}^0K^-$ and $D^-\\bar{K}^0$ scattering lengths with 14–27% uncertainties; the effective ranges come out with larger errors, and the distant $D^-_s\\eta$ channel is only poorly constrained.","Because the inverse-analysis potential includes energy-dependent terms, a pure molecular state and a state with a genuine nonmolecular component would be accommodated differently; the size of the threshold enhancement in the data is what tells whether the resonance is coupled to $\\bar{D}\\bar{K}$.","The decay $\\Lambda_b \\to \\Lambda_c \\bar{D}^0K^-$ has already been observed, so only the measurement of the mass distribution, rather than the discovery of a new decay mode, is needed to apply the method."],"supporting_citations":[{"why":"The experimental observation of the decay, which motivates the study and would supply the mass distributions to which the inverse method is applied.","marker":"[26]"},{"why":"The local hidden gauge approach that provides the vector-meson-exchange potential used to generate the D*_s0(2317) bound state and the pseudodata.","marker":"[27–30]"},{"why":"The source of the coupled-channel potential matrix, loop function, effective-range relations, and the correlation-function baseline with 20 MeV and 60% uncertainties that this paper improves upon.","marker":"[34]"},{"why":"The lattice QCD reanalysis that gives roughly 72% for the DK probability, the number the extracted ~70% result is compared with.","marker":"[25]"},{"why":"The resampling/bootstrap procedure used to propagate the assumed 5% data errors into uncertainties on the fitted observables.","marker":"[31–33]"},{"why":"The derivation of channel probabilities from couplings and loop-function derivatives, used to compute the molecular probability of the state.","marker":"[35, 37, 38]"}],"fun_headline_variants":["D*_s0(2317) mass pinned within 4 MeV by decay shapes","Lambda_b decays reveal D*_s0(2317) binding to 4 MeV","Near-threshold shapes fix D*_s0(2317) mass and DK probability","DK molecule probability of D*_s0(2317) from Lambda_b decays","Precise binding of D*_s0(2317) from near-threshold spectra"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quoted 4 MeV and 11% uncertainties assume that the pseudodata generated by the authors' local hidden gauge model, with 5% relative errors and a fitting window of 50 MeV above threshold, faithfully represent what the actual experiment will measure for these decay distributions.","fun_headline_variants_meta":{"raw":{"variants":["D*_s0(2317) mass pinned within 4 MeV by decay shapes","Lambda_b decays reveal D*_s0(2317) binding to 4 MeV","Near-threshold shapes fix D*_s0(2317) mass and DK probability","DK molecule probability of D*_s0(2317) from Lambda_b decays","Precise binding of D*_s0(2317) from near-threshold spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000685,"raw_usage":{"total_tokens":3235,"prompt_tokens":1200,"completion_tokens":2035,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":816,"completion_tokens_details":{"reasoning_tokens":1919}},"tokens_in":816,"tokens_out":2035,"duration_ms":15513,"temperature":1.0,"reasoning_tokens":1919,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:32:59.812365+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the $\\bar{D}^0K^-$ and $D^-\\bar{K}^0$ invariant mass distributions in the observed $\\Lambda_b \\to \\Lambda_c \\bar{D}^0K^-$ decay and the companion $\\Lambda_b \\to \\Lambda_c D^- \\bar{K}^0$ decay: if the sharp near-threshold enhancement over phase space is absent, or if fitting the inverse procedure to data generated by a different production model shifts the recovered pole by much more than 4 MeV, the central precision claim is falsified.","supporting_citations":[],"review_version":1}