{"id":"50b779f3-be8d-476a-bc91-531586c8ea7f","arxiv_id":"2506.23476","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Femtoscopic correlation functions for D0K+ and D0Dbar*0 pairs are predicted to be sensitive to the molecular versus bare-state composition of D_s0*(2317) and X(3872).","lead":"This paper computes momentum correlation functions for D0K+ and D0Dbar*0 pairs to see whether the internal structure of two exotic hadrons, D_s0*(2317) and X(3872), changes their femtoscopic signals. It finds that the shapes depend on how much of each state is a hadron molecule versus a bare quark state, which could turn femtoscopy into a compositeness probe.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed lineshape discrimination may be a regulator/source artifact: qmax and R are fixed to favorable values, and the paper's own uncertainty bands are broad exactly for P ~ 0.7; a regulator/source scan is needed to test whether the predicted separations survive.","rationale":"The reader's conditional verdict identifies the right weak point: the predicted lineshape differences are computed with a single Gaussian source and a sharp-cutoff contact EFT at favorable parameter values. My stress-test narrows this to a concrete, testable risk: the qmax-dependence is explicitly acknowledged to be large precisely in the physically relevant low-compositeness regime, so the central discrimination claim stands or falls on whether the compositeness curves remain separated when the regulator and source are varied. This is not an internal inconsistency but an external robustness question, and it is resolvable by a numerical scan. I also note that the inverse-problem demonstration in Sec. III.C is self-consistency rather than independent validation, since synthetic data are generated from and fit with the same potential form; this reinforces the need for a regulator/source scan. The paper is transparent about the Wigner-bound issue and about cutoff sensitivity, which makes the proposed check meaningful rather than adversarial. For these reasons, I recommend keeping the verdict as CONDITIONAL pending the robustness test, rather than upgrading to acceptance or moving to rejection.","tokens_in":12878,"tokens_out":5667,"duration_ms":68397,"concrete_test":"Recompute the D0K+ CFs in Scenarios II and IV for P = 0.4, 0.7, 1.0 and for the two bare-state masses, using qmax in {0.5, 0.75, 1.0, 1.25, 1.5} GeV and a smooth Gaussian regulator, with R in {0.8, 1.0, 1.5} fm. At each k below 50 MeV, compare the separation between compositeness curves with the quadrature sum of the regulator/source bands; require the separation to exceed the combined band width to retain the 'clearly distinguishable' claim. For X(3872), repeat the same scan with a finite-range interaction tuned to r0 = -1 fm (Wigner-consistent) and check whether the CF changes by more than the current Fig. 11 separation between compositeness curves.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is that the low-momentum lineshape differences attributed to compositeness and bare-state admixture survive when the short-range regulator and source parameters are varied. In Sec. II (Eqs. (2) and (4)) the CF depends on qmax, and Sec. III fixes qmax = 1 GeV and R = 1 fm because that choice maximizes the deviation from unity. The paper itself notes (Sec. III.B, text after Fig. 8) that for small compositeness the low-momentum CF has a broad qmax dependence. Since the Ds0* scenarios of interest use P1 ~ 0.7, the claim of 'clearly distinguishable' lineshapes requires the P=0.7 and P=1 curves to be separated by more than the qmax-induced band. If a qmax scan from 0.5 to 1.5 GeV makes the bands overlap, the predicted sensitivity is an artifact of the regulator choice rather than a robust observable signature. The X(3872) section has a separate acknowledged red flag: the single-channel zero-range fit in Sec. III.D gives r0 = 0.28 fm, violating the Wigner bound, so the quantitative X(3872) CF predictions rest on a model whose low-energy effective range is not physical.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies femtoscopic correlation functions as a probe of the internal structure of the D_s0*(2317) and X(3872) states. Using a contact-range effective field theory, the authors construct DK and Dbar*D potentials for four scenarios of D_s0*(2317) (pure molecule, molecule plus bare state, DK-D_s eta coupled channel, and the latter plus a bare state) by reproducing the known mass and an assumed compositeness. They then compute D0K+ correlation functions and find that the line shape depends on the molecular fraction and on the mass of the bare state. For X(3872), they use a similar zero-range model, with and without a bare charmonium state, and compute D0Dbar*0 correlation functions, claiming clearly distinguishable line shapes for different compositeness values. The paper also presents a synthetic-data inverse-problem exercise to extract compositeness from correlation functions.","tokens_in":13156,"tokens_out":8856,"duration_ms":83061,"significance":"If the predicted line-shape differences survive a broader variation of model inputs, the paper would establish femtoscopy as a genuinely useful discriminator of molecular versus compact content for D_s0*(2317) and X(3872), complementing spectroscopy and lattice QCD. The paper has concrete strengths: it works within a well-defined EFT framework, it compares the derived DK scattering lengths with lattice QCD results, and it explicitly tests the inversion of correlation functions into compositeness using synthetic data. These are useful proof-of-principle steps. The significance is currently limited, however, by the strong model dependence of the central claims: the potentials are fitted to the quantities whose sensitivity is then advertised, and the regulator and source-size dependence is not systematically quantified.","major_comments":[{"comment":"The central claim that the D0K+ correlation-function line shapes are 'clearly distinguishable' for different compositeness values is not supported by a robustness analysis. The authors choose q_max = 1 GeV and R = 1 fm, explicitly noting in Sec. III.A that R = 1 fm gives the largest deviation from unity, and they state in Sec. III.B that for smaller compositeness the low-momentum correlation function has a broad q_max dependence. Since the D_s0*(2317) scenarios of interest use P ~ 0.7, the separation between the P = 0.7 and P = 1 curves must be shown to exceed the bands obtained by scanning q_max (e.g., 0.5-1.5 GeV) and R (e.g., 0.5-2 fm). Without such a scan, the claimed discrimination could be an artifact of the selected regulator and source parameters rather than a robust observable signature.","section":"Sec. III.A, Fig. 4 and Sec. III.B, Fig. 7"},{"comment":"The inverse-problem demonstration does not support the Summary's claim of a 'bijective relationship between compositeness and CFs'. The synthetic data are generated from the same potential model (Scenario IV) and fitted with the same functional form, the same q_max, and the same R, so the exercise tests the numerical stability of the fit, not whether a measured correlation function uniquely determines compositeness in a model-independent way. To support the stronger claim, the authors would need to show that alternative interaction models (different V forms, cutoffs, or additional coupled channels) that reproduce the same correlation function do not lead to different extracted compositeness, or that the extracted value is stable under such variations.","section":"Sec. III.C"},{"comment":"The X(3872) correlation functions are computed with a single-channel zero-range model that yields r0 = 0.28 fm, which the authors acknowledge violates the Wigner bound for a zero-range interaction. Since X(3872) is a shallow bound state with a large scattering length, the low-momentum correlation function is sensitive to the effective range, and the acausal zero-range model may produce line shapes that differ from those of a causal model (e.g., the bare-state model, which gives r1 = -4.72 fm). The paper should either use a potential with r0 < 0 throughout the X(3872) section or explicitly demonstrate that the correlation-function predictions are insensitive to the effective range over the momentum range shown.","section":"Sec. III.D, Figs. 9 and 10"},{"comment":"The channel basis used in the coupled-channel correlation-function calculations is not specified consistently. Equation (12) defines a potential matrix for DK-D_s eta, but the text in Sec. III.B attributes the coupled-channel effect on the correlation functions to D0K+ and D+K0 rather than to D_s eta. This makes it impossible to determine whether the T-matrix in Eq. (10) is evaluated in the (DK, D_s eta) basis or in a (D0K+, D+K0) basis, and it undermines the stated claim that the correlation functions are sensitive to D+K0 admixture. The authors should define the channels explicitly and, if D+K0 is included as a coupled channel, provide the corresponding potential matrix and loop functions.","section":"Sec. III.B, Eqs. (10)-(13) and Figs. 6-8"}],"minor_comments":[{"comment":"The name 'Godfrey-Isgur' is misspelled as 'Goldfrey-Isgur' in the first paragraph of the Introduction.","section":"Introduction"},{"comment":"The sentence 'The corresponding CFs are shown in Fig. 7' appears to refer to Fig. 5, which displays the correlation functions for different bare-state masses; Fig. 7 shows the compositeness dependence in the coupled-channel case.","section":"Sec. III.A"},{"comment":"The phrase 'we take a point about every 2MeV at at momentum less than 40' contains a duplicated 'at' and the momentum value lacks units; it should presumably read 'every 2 MeV at momenta less than 40 MeV/c' or similar.","section":"Sec. III.C"},{"comment":"The text 'we obtain the scattering length a_D0Dbar*0 ≈ 19.58 fm and effective range r = -1 fm' uses 'r' without a subscript, inconsistent with Eq. (15); please use r0 or r1 consistently for the effective range.","section":"Sec. III.D"},{"comment":"The sentence 'we assume that the weights of the DK channel and the D_s eta channel are the same' is inconsistent with the immediately following ratio omega_DK/omega_Ds eta = 1/0.35 obtained from Eq. (21); please clarify which weights are assumed equal and which are obtained from the thermal-weight estimate.","section":"Sec. III.B"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the journal's scope well. The main technical concerns are the regulator/source robustness of the central discrimination claim, the circularity of the inverse-problem argument, the acausal X(3872) model, and the unclear channel basis in the coupled-channel calculation. All four are addressable within the manuscript's scope, so I do not recommend rejection, but they need to be fixed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I've read it. The genuinely new piece is the channel-by-channel predictions for D0K+ and D0Dbar*0 femtoscopic correlation functions with bare-state dressing, and the idea that the low-momentum CF could locate the bare-state pole. The formalism is standard Koonin-Pratt with a contact EFT, the paper is clearly written, and the scattering lengths are checked against lattice QCD, which is a useful sanity check. The X(3872) Wigner-bound violation is acknowledged, which is honest. The soft spots are real, though. The most load-bearing claim, that lineshapes are 'clearly distinguishable' as a function of compositeness, is not yet supported. The authors fix R=1 fm because it gives the largest deviation from unity, and qmax=1 GeV, and the paper itself notes that for P~0.7 the low-momentum CF has a broad qmax dependence. Without a scan over R and qmax, the claimed sensitivity could be a regulator/source artifact. That's a testable concern; the authors should run it. The inverse-problem section is more circular than it looks. Synthetic data generated from the same model and then refitted recovers the input; that's a self-consistency check, not a demonstration that compositeness can be extracted from real data. The word 'bijective' is too strong. The X(3872) section has an additional problem: the zero-range model gives r0=0.28 fm, which violates the Wigner bound, and the paper's explanation—'we set the cutoff to 1 GeV'—doesn't resolve it. So the quantitative X(3872) predictions rest on a model whose low-energy effective range is not physical. Bottom line: this is a serious paper worth refereeing, but the central claim needs a robustness scan and the inverse-problem language should be softened. I'd send it out, with a request for those changes.","headline":"A clearly written paper with genuinely new channel-specific predictions, but the central claim of compositeness-sensitive lineshapes needs a regulator/source robustness scan before I would trust it.","tokens_in":13718,"tokens_out":2664,"would_cite":true,"duration_ms":27902,"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":"The paper claims that femtoscopic correlation functions of D0K+ and D0Dbar*0 pairs can distinguish whether the D_s0*(2317) and X(3872) are pure molecules, mixtures with a bare quark-model state, or coupled-channel composites, and that the…","keywords":["femtoscopy","correlation functions","compositeness","hadronic molecules","exotic hadrons","effective field theory","D_s0(2317)","X(3872)"],"falsifier":"Measure the D0K+ correlation function in high-multiplicity collisions with a source radius near 1 fm and fine momentum resolution below k = 50 MeV: if no low-momentum peak or enhancement appears in a scenario where a bare state near 2.348 GeV is present, the claim that the correlation function can probe the bare-state position would be contradicted.","tokens_in":12637,"feed_emoji":"🔬","tokens_out":5788,"duration_ms":57527,"temperature":0.7,"pith_summary":"The paper argues that femtoscopic correlation functions measured in high-energy collisions can act as a structural probe for hadrons that sit near two-particle thresholds. By modeling the D_s0*(2317) as a DK molecule, with or without a DK-D_s eta coupled channel and with or without an underlying bare csbar state, the authors predict D0K+ correlation functions and find that their low-momentum lineshape is sensitive to the bare-state admixture and to the coupled D+K0 channel. They further predict that a bare state's mass, even if it lies below or above the DK threshold, leaves a distinctive low-momentum signature in the correlation function, so the function can probe a bare state that is otherwise invisible. For the X(3872), modeled as a shallow D0Dbar*0 molecule with a bare ccbar component, the D0Dbar*0 correlation functions change shape strongly with the molecular fraction, offering a way to extract compositeness. The paper also shows, through an inverse-problem fit to synthetic data, that the compositeness can be recovered from correlation functions.","feed_headline":"Femtoscopy lineshapes expose bare states inside exotic hadrons","feed_subtitle":"The D0K+ correlation function can reveal a hidden bare state's mass, and D0Dbar*0 lineshapes map the X(3872) molecular fraction.","key_machinery":"The calculations rest on the Koonin-Pratt formula, which writes the correlation function as an integral over the emission source (taken as a Gaussian of radius R) of the squared relative wave function. The wave function is built from the T-matrix, obtained by solving a Lippmann-Schwinger equation with a contact-range potential; a bare quark-model state is included as an energy-dependent pole term alpha/(sqrt(s)-m_bare), and coupled channels enter through a matrix potential with a fixed ratio between DK and D_s eta couplings. The compositeness is computed from the pole residue and the derivative of the loop function, and the scattering length and effective range are read off from the inverse T-matrix near threshold. The key mechanism is that the energy dependence introduced by the bare pole changes the low-momentum correlation function in a way that a purely constant contact interaction cannot.","core_discovery":"The central discovery is that the lineshape of a two-hadron momentum correlation function encodes more than the scattering length: it distinguishes a pure molecular state from a state mixed with a bare quark-model seed. Using contact effective-field-theory potentials fixed to the physical masses and assumed compositeness values, the authors show that adding a bare-state pole term alpha over (sqrt(s) minus m_bare) substantially changes the D0K+ correlation function, and that the position of the pole, whether m_bare sits below or above the DK threshold, produces a low-momentum peak or enhancement. In the coupled DK-D_s eta case, the D_s eta channel contributes little, while the D+K0 coupled channel matters. For X(3872), the D0Dbar*0 correlation function falls below unity in a compositeness-dependent way, and the bare-state dressing modifies both scattering length and effective range, with the effective range even becoming negative in some scenarios. The paper establishes a bijective (one-to-one) relationship between compositeness and correlation functions, at least within the model, by recovering the input compositeness from a fit to synthetic correlation-function points.","pith_inferences":["The same formalism could be applied to other threshold-bound exotics, such as the T_cc tetraquark or candidate pentaquarks, to map out bare-state admixtures wherever a quark-model seed is suspected.","If future data resolve the low-momentum correlation function with high precision, the predicted peak position for a bare state could be converted into a mass measurement with accuracy limited by the source size R.","The sensitivity to cutoff and source size at low momentum suggests that combining correlation functions measured at several source sizes (from different collision systems) might isolate short-range from long-range contributions."],"forward_implications":["If the D0K+ correlation function is measured in high-energy collisions, a low-momentum peak or enhancement would indicate the presence and position of a bare csbar state below or near the DK threshold.","The compositeness of the D_s0*(2317) could be extracted from the measured D0K+ correlation function, distinguishing a 100% DK molecule from a roughly 70% molecular mixture.","For the X(3872), measured D0Dbar*0 correlation functions would provide a compositeness diagnostic through the depth and sign of the correlation function below unity.","The D_s eta channel has negligible influence on the D0K+ correlation functions, so a single-channel DK analysis may suffice for this observable."],"supporting_citations":[{"why":"Supplies the foundational Koonin formula that relates the correlation function to the source size and the two-particle wave function.","marker":"[78]"},{"why":"Extends the formula to the relativistic source averaging used in heavy-ion femtoscopy.","marker":"[79]"},{"why":"Provides the model-independent effective field theory approach used to extract hadron-hadron interactions and compositeness.","marker":"[77]"},{"why":"Supplies the coupled-channel DK-D_s eta potential used in Scenarios III and IV.","marker":"[32]"},{"why":"Gives the lattice QCD determination of the DK scattering length used to validate the extracted potential.","marker":"[44]"},{"why":"Provides the bare-state pole potential alpha/(sqrt(s)-m_bare) used to dress the molecule.","marker":"[81]"},{"why":"States the Wigner bound/causality constraint on the effective range, cited for the X(3872) zero-range effective range discussion.","marker":"[88]"},{"why":"Provides the coupled D0Dbar*0 - D+D*- potential used for the X(3872) analysis.","marker":"[89]"}],"fun_headline_variants":["Correlation lineshapes reveal hidden bare states in exotic hadrons","Femtoscopy maps molecular and bare fractions in X(3872)","D0K+ correlation function probes bare-state pole position","Hadron femtoscopy distinguishes molecule from bare admixture"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predicted lineshape differences rely on the assumption that a single Gaussian source of radius R = 1 fm and a contact-range interaction with a sharp cutoff q_max = 1 GeV accurately describes the low-momentum correlation function; if real sources are non-Gaussian or the cutoff dependence is stronger than modeled, the distinguishing features could wash out.","fun_headline_variants_meta":{"raw":{"variants":["Correlation lineshapes reveal hidden bare states in exotic hadrons","Femtoscopy maps molecular and bare fractions in X(3872)","D0K+ correlation function probes bare-state pole position","Hadron femtoscopy distinguishes molecule from bare admixture"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000222,"raw_usage":{"total_tokens":1533,"prompt_tokens":1107,"completion_tokens":426,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":723,"completion_tokens_details":{"reasoning_tokens":355}},"tokens_in":723,"tokens_out":426,"duration_ms":4486,"temperature":1.0,"reasoning_tokens":355,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:41:30.747575+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the D0K+ correlation function in high-multiplicity collisions with a source radius near 1 fm and fine momentum resolution below k = 50 MeV: if no low-momentum peak or enhancement appears in a scenario where a bare state near 2.348 GeV is present, the claim that the correlation function can probe the bare-state position would be contradicted.","supporting_citations":[{"cited_title":"$P$-wave charmonium contribution to hidden-charm states from reanalysis of lattice QCD data","cited_arxiv_id":"2410.19563","evidence_quote":"Provides the bare-state pole potential alpha/(sqrt(s)-m_bare) used to dress the molecule."},{"cited_title":"Implication of a negative effective range on the $D\\bar{D}^*$ interaction and the nature of $X(3872)$","cited_arxiv_id":"2409.06409","evidence_quote":"Provides the coupled D0Dbar*0 - D+D*- potential used for the X(3872) analysis."}],"review_version":1}