{"id":"7440ccbc-8e3c-4f9c-96f4-5d4eec8b92b9","arxiv_id":"hep-ph/0402204","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A heavy-quark dispersion relation links the neutral D meson mass difference to an integral of the width difference over heavy meson mass, and a phase-space model suggests the Standard Model can give a mass difference of the same order as the width difference.","lead":"Physicists derive a heavy-quark dispersion relation connecting the mass difference between the two neutral D meson states to an integral of their decay-width difference over meson mass. A phase-space model then suggests the Standard Model can produce a mass difference comparable to the width difference, which changes how D mixing experiments are read as new-physics signals.","discovery_kind":"unclear","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"The derivation is sound and the model dependence of the numerical estimate is real but explicitly acknowledged; the one under-quantified point is the low-E weight in the 4P dispersion integral, which determines whether the HQET relation is used within its regime.","rationale":"The derivation of Eq. (19) is the solid part of the paper: stripping the heavy-quark mass from the correlator via HQET gives a function with a physical interpretation at every E, and the unsubtracted dispersion relation is justified by the 1/E^2 falloff of ΔΓ(E) in the model. I checked the factor-of-two in Eq. (23) and the analytic structure of the correlator; both are consistent. The numerical claim is softer, as the reader says, and the authors are unusually candid about it. My only addition is that the claim 'the 4P integral is dominated by E near m_D' is asserted rather than shown, and it is the load-bearing reason the 4P result is treated as more trustworthy than the PP result. The proposed test would settle whether that reason holds. If it does, the reader's CONDITIONAL verdict stands as written; if it does not, the verdict should remain CONDITIONAL but with the headline range explicitly widened. Either way, the reader's identification of the modeled g_F(E) as the weakest assumption is correct, and I would not move the verdict.","tokens_in":13533,"tokens_out":20005,"duration_ms":362983,"concrete_test":"Compute the integrand of Eq. (23) for the (4P)_8 representation using g_4P(E) with m1 = 0.8 GeV and m2 = 3 GeV, and report the fraction of the principal-value integral contributed by E < 1 GeV. If that fraction exceeds roughly 30%, the 4P result is dominated by the region where the O(Λ_QCD/E) corrections are O(1), and the headline range should be widened or the claim downgraded accordingly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader's weakest_assumption is correct: the headline number r(4P) ~ -1 in Fig. 4 is set by the assumed form of g_4P(E) in Eq. (30), and the two models differ by a factor of two. The authors acknowledge this repeatedly (Sec. IV.B, Sec. V: 'we can trust it only at the order of magnitude level'). I do not present that as a new finding. What remains under-quantified is a specific technical point on which the reliability of the 4P result rests: the paper asserts that 'because the decay is close to threshold, the dispersion integral is dominated by E near m_D' (Sec. IV.B), which is what distinguishes the 4P case from the PP case where the authors themselves say 'we should not trust this result.' That assertion is not demonstrated. The dispersion relation Eq. (19) carries an explicit O(Λ_QCD/E) correction which the paper says is O(1) at small E; if the 4P integrand of Eq. (23) has non-negligible support below ~1 GeV, the 4P result inherits the same low-energy unreliability as the PP result, and the quoted range |x| ~ 10^-3 to 10^-2 is correspondingly less secure. This is a checkable gap, not a re-statement of the acknowledged model dependence.","agreement_with_reader":"agree"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"The manuscript derives a heavy-quark dispersion relation, Eq. (19), that expresses the neutral D mass difference Δm as a principal-value integral over E of the width difference ΔΓ(E) of a heavy meson of mass E. The derivation rests on the heavy-quark-mass independence of the HQET correlator Σ_v(q) defined in Eq. (15) and on the identification Σ_v = −2Δm(E) + iΔΓ(E) at q² = E², Eq. (17). The relation is then evaluated using the phase-space model of SU(3) breaking of Ref. [3] together with assumed energy dependences g_F(E) for the PP and 4P partial widths, Eqs. (27) and (30). The output is the ratio r_FR = x_FR/y_FR, of order −0.1 to −1 for the 4P multiplets, which translates into |x| roughly between 10⁻³ and 10⁻² if y is near 1%, with x and y of opposite sign.","tokens_in":27812,"tokens_out":11791,"duration_ms":217833,"significance":"The dispersion relation itself is the durable contribution: it is derived without free parameters and converts a model of ΔΓ(E) — an observable that can in principle be computed or measured at other heavy-meson masses — into a prediction for Δm, which is not directly calculable. It yields a falsifiable, if crude, prediction (x and y comparable in magnitude and of opposite sign) that bears directly on how the current and forthcoming D-mixing bounds should be interpreted, and it sharpens the argument that a percent-level x need not by itself signal new physics. The authors are explicit that the numerical estimate is order-of-magnitude only and that the PP channel is not under control; that candor is appropriate and should be preserved in revision.","major_comments":[{"comment":"The statement that 'because the decay is close to threshold, the dispersion integral is dominated by E near m_D' is the stated reason the 4P result is more trustworthy than the PP result, but it is asserted rather than shown. This matters because Eq. (19) carries an explicit O(Λ_QCD/E) correction that the text itself describes as O(1) at small E; if the 4P integrand of Eq. (23) has non-negligible support below about 1 GeV, the 4P number inherits the same low-energy unreliability that leads the authors to discount the PP channel. Please show the integrand d r_FR/dE, or the cumulative integral, for (4P)_8, (4P)_27 and (4P)_27′, and quote the fractional contribution to r_FR from E below 1 GeV and below 1.5 GeV. A companion cutoff study, varying the lower limit of the integral between roughly 0.8 and 1.5 GeV and quoting the change in r_FR, would settle the question directly.","section":"Sec. IV.B, Eq. (23) and Fig. 4","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The opposite-sign conclusion is presented as a prediction but its basis is not given. The sign follows if ΔΓ(E) is positive and decreasing, so that the E < m_D side of the principal value dominates the integral in Eq. (19); however, the paper's own model produces negative y_FR for (PP)_8 and (PP)_27, so ΔΓ(E) is not positive-definite multiplet by multiplet. Please state explicitly the assumption on ΔΓ(E) under which the sign prediction holds, and indicate whether it survives if individual multiplets enter with either sign or if ΔΓ(E) is non-monotonic over the relevant range. Since this is one of the few falsifiable statements in the paper, it should be stated with its conditions attached.","section":"Sec. V, conclusion that x and y are of opposite sign","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The working equation requires y_FR(E), the width-difference-to-width ratio at heavy-meson mass E, but the manuscript never states how this quantity is obtained for E ≠ m_D. Presumably the phase space factors Φ(n) entering Eq. (21) are recomputed at mass E, but this should be said in one explicit sentence, and the resulting asymptotic behaviour (the m_s²/E² falloff argued for in Sec. IV.A) should be shown to be what actually enters Eq. (23). As written, a reader cannot reproduce Fig. 4 from the text alone, which is a reproducibility gap in the one place where the paper produces numbers.","section":"Sec. IV, Eq. (23)","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The O(Λ_QCD/E) term inside the square bracket is dimensionally inconsistent as displayed: ΔΓ has mass dimension while Λ_QCD/E is dimensionless. It should be written as a relative correction multiplying the integrand, i.e. as ΔΓ(E)/(E − m_D) × [1 + O(Λ_QCD/E)].","section":"Eq. (19)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The assertion that Σ_v(q) 'only depends on q²' is stronger than what is needed and is not justified as stated. With an external state of fixed velocity v, the correlator can in general depend on q·v as well as q²; the dispersion relation only requires analyticity in q₀ at fixed v, along the path q = (E, 0) where q·v = E. Either justify the stronger claim or weaken it to the statement actually used.","section":"Sec. II, Eqs. (17)–(18)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The local ΔC = 2 contribution is dropped on the grounds that 'we are interested in the question whether x could be near the percent level.' The quoted range, however, extends down to |x| ~ 10⁻³, where the neglected term (estimated at or below 10⁻³ in Refs. [10–12]) is comparable. The lower end of the range should carry that caveat explicitly.","section":"Sec. II, following Eq. (13)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The vertical scales of the two figures differ by roughly an order of magnitude (0–20 versus 0.5–2.5), which makes the GMO sensitivity difficult to read off. A common scale or an inset would help, and the text should state the factor by which y_PP must change to keep x_PP approximately fixed when the GMO relation is imposed.","section":"Figs. 2 and 3","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Setting X = 0 is supported by Refs. [19]; it would strengthen the argument to quote the bound on |X| that those references actually give, so the reader can see how much the E-dependence of g_PP could shift.","section":"Sec. IV.A, Eq. (25)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The text notes that 'one can reconstruct x from x_FR using a relation analogous to Eq. (22),' but the total x is never actually formed; only the range of r_FR is quoted. Please state explicitly that the estimate of total x assumes r_FR to be roughly universal across the 4P multiplets that dominate y, since that assumption is what licenses quoting a range of r as a range of x/y.","section":"Sec. IV.B, reconstruction of total x","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"There is a typographical error in Sec. I: 'comes form double insertion' should read 'comes from double insertion.' Please also check the typesetting of Eq. (2), where x = Δm/Γ appears in the source as 'x = ∆mΓ'.","section":"Sec. I and Sec. III","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The paper distinguishes its relation from earlier dispersion treatments (Ref. [13]) by noting that in those works Im Σ(s) has no physical interpretation for s ≠ m_D². One concrete sentence explaining why the earlier constructions could not relate Δm to a physical ΔΓ(E) would make the novelty of Eq. (19) sharper than the current brief remark.","section":"Sec. II, comparison with Ref. [13]","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"major_revision"},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"The part worth your time is Sec. II, not the number in the abstract. Falk et al. show that by stripping the heavy quark mass out of the correlator — writing the charm field in HQET and keeping only the term that survives mc → ∞ — you get a function Σ_v(q) with a physical meaning at every q, not just at q = p_D. That is what the earlier dispersion treatments (their Ref. [13]) could not do, and it is what makes Eq. (19) possible: Δm as a principal-value integral of ΔΓ(E) over heavy meson mass. The analyticity argument is textbook and the derivation is explicit about what gets dropped, including the oscillating term replaced by local ΔC=2 operators, which they argue sits below 10^-3. That is a genuine construction and it will get cited on its own merits regardless of the numbers that follow it. The numerical claim is much softer, and they say so. r(4P) lands between roughly -0.1 and -1 depending on whether you use g_4P or g'_4P in Eq. (30), and for the PP case they tell you outright not to trust the result. The honest reading is that the sign and rough magnitude of x/y are what the model gives, not what the theory predicts. On the stress-test point: the paper asserts in Sec. IV.B that the 4P integral is dominated by E near m_D without arguing it. But the left panel of Fig. 4 is the argument — r(4P) is flat as m1 runs from 0.2 to 1 GeV, which is exactly what you expect if the low-E region contributes little. Contrast Fig. 2, where r_PP swings by a factor of several over the same range. So the claim is supported by their own numerics; they just did not connect the dots in the text. The residual worry is the uncomputed O(Λ_QCD/E) correction, which they flag as O(1) at small E. If the low-E weight really is small, that correction is weighted down too, but nobody has checked it. Minor, and fixable with a sentence. The interpretive payoff is what matters for the field: a percent-level x is not by itself evidence of new physics, and the only clean signal would be CP violation. Useful to have on record. Who should read it: anyone quoting D-mixing bounds or building a new-physics case on x ≫ y. Cite the derivation; cite the number with the caveat attached. Send it out.","headline":"The HQET dispersion relation in Sec. II is the real contribution and it is clean; the x/y number is a model estimate the authors themselves label order-of-magnitude, and they are right to.","tokens_in":14477,"tokens_out":2871,"duration_ms":53247},"pith_extraction":{"msc":[],"pacs":["12.39.Hg","13.25.Ft","14.40.Lb"],"model":"selfhosted/deepseek-v4.1-flash","headline":"The authors derive a heavy-quark dispersion relation that expresses the neutral D meson mass difference as a principal-value integral of the width difference over a range of heavy-meson masses, then use a phase-space model of SU(3) breaking","keywords":["D0-D0bar mixing","dispersion relation","heavy quark effective theory","SU(3) breaking","final state phase space","charm meson decays","mass and width differences"],"falsifier":"A first-principles determination of how the heavy-meson width difference varies with heavy-meson mass, or improved measurements of both $x$ and $y$: a measured $x$ far above $y$, or a percent-level $x$ with the same sign as $y$, would contradict the estimate.","tokens_in":14966,"feed_emoji":"⚛️","tokens_out":2777,"duration_ms":50149,"temperature":0.7,"pith_summary":"Neutral D mesons can turn into their antiparticles, but in the Standard Model this mixing vanishes in the limit of exact SU(3) flavor symmetry, so its size depends entirely on how badly that symmetry is broken. The authors ask whether the mass difference $x = 2M_{12}/2m_D$-driven quantity $x = 2M_{12}/2m_D$ can be as large as the width difference $y = 2M_{12}/2m_D$... ","feed_headline":"Standard Model may give D mixing near one percent","feed_subtitle":"A heavy-quark dispersion relation links the D mass difference to the width difference, suggesting the two are comparable.","key_machinery":"The machinery is a correlator of two weak Hamiltonians sandwiched between heavy-meson states, $g_F(E)$, the energy dependence of the partial width to a given final state.","core_discovery":"The paper's central claim is a new relation plus a numerical consequence drawn from it. The relation, Eq. (19), states that the neutral D mass difference is a principal-value integral over the width difference of a heavy meson taken at varying heavy-meson mass, $x$ is comparable in magnitude to $y$ and of opposite sign, with $|x|$ between roughly $10^{-3}$ and $10^{-2}$.","pith_inferences":["Because the relation is an integral transform, a lattice or dispersive determination of how $\\Delta\\Gamma$ falls with heavy-meson mass would convert this estimate into a real prediction without any new theoretical idea.","The same transform could be applied to the $B_s$ system as a consistency check: there $\\Delta\\Gamma$ is measured and $\\Delta m$ is measured, so the relation's convergence and its neglected $O(\\Lambda_{\\rm QCD}/E)$ corrections could be tested where both sides are known.","If the sign correlation between $x$ and $y$ survives better modeling, it is a cheap discriminant: a same-sign measurement at the percent level would be hard to accommodate without new physics or a very different SU(3)-breaking pattern.","The paper's logic suggests that any charm mixing prediction built purely on short-distance operator products will underestimate $x$, since the effect lives at thresholds that the expansion cannot resolve."],"forward_implications":["A percent-level measurement of $x$ in D mixing would no longer be, by itself, a signal of physics beyond the Standard Model.","The sign of $x$ relative to $y$ becomes a diagnostic: the phase-space mechanism predicts opposite signs, so a same-sign percent-level result would point elsewhere.","The result gives experimental searches a concrete target region, $|x| \\sim 10^{-3}$ to $10^{-2}$, rather than a generic 'very small' expectation.","It identifies threshold and multi-body final states as the place where the Standard Model's D mixing is generated, which is precisely where operator-product-expansion treatments are blind."],"supporting_citations":[{"why":"Supplies the phase-space model of SU(3) breaking that generates the input width difference $y$ and its estimates for the $PP$ and $4P$ representations.","marker":"[3]"},{"why":"Originates the technique of deriving a dispersion relation in the heavy quark mass within heavy quark effective theory, which this paper adapts.","marker":"[8]"},{"why":"Provides the heavy quark effective theory conventions and the mass-independent state normalization used to define the correlator $Sigma_v$.","marker":"[9]"},{"why":"Gives the large-mass scaling of the heavy-to-light form factor, used to fix the high-energy behavior of $g_{PP}(E)$.","marker":"[18]"},{"why":"Estimate the local $Delta C = 2$ contributions at or below the $10^{-3}$ level, justifying their neglect in the derivation.","marker":"[10-12]"},{"why":"Earlier dispersion-relation treatments of the D mass difference, whose spectral density lacked a physical interpretation for $s eq m_D^2$; the contrast defines the novelty here.","marker":"[13]"},{"why":"Current experimental bounds on $x$ and $y$ at the few-percent level, which set the scale against which the estimate is judged.","marker":"[4]"}],"fun_headline_variants":["Dispersion relation ties D mass gap to width difference","Heavy-quark integral links D mixing mass and width","D meson mass difference may rival its width difference","SU(3) phase-space breaking shifts D mixing estimate","D mass gap may be comparable to width, opposite sign"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The numerical answer rests entirely on the hand-modeled energy dependence of the partial widths $g_F(E)$ in Eqs. (27) and (30), whose free parameters and two alternative functional forms move the result by about a factor of two, and the authors say plainly that for the two-pseudoscalar case 'we should not trust this result.'","fun_headline_variants_meta":{"raw":{"variants":["Dispersion relation ties D mass gap to width difference","Heavy-quark integral links D mixing mass and width","D meson mass difference may rival its width difference","SU(3) phase-space breaking shifts D mixing estimate","D mass gap may be comparable to width, opposite sign"]},"model":"selfhosted/deepseek-v4.1-flash","effort":"low","attempts":1,"cost_usd":0.0,"raw_usage":{"total_tokens":2224,"prompt_tokens":720,"completion_tokens":1504,"prompt_tokens_details":null,"completion_tokens_details":{"reasoning_tokens":1437}},"tokens_in":720,"tokens_out":1504,"duration_ms":26864,"image_pages":[],"temperature":1.0,"reasoning_budget":50,"reasoning_tokens":1437,"transport_provider":"selfhosted","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-11-13T08:52:35.718582+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"pipeline_version":"pith-pipeline@v1.5.0-archive-selfhost","input_kind":null,"input_chars":null,"input_total_chars":null,"passes_missing":null,"published_standing":null,"published_confidence":null,"objection_audit":null,"run_record":null,"basis":{"kind":"not_recorded","source":null,"chars_given":null,"chars_total":null},"receipt":null,"falsifier":"A first-principles determination of how the heavy-meson width difference varies with heavy-meson mass, or improved measurements of both $x$ and $y$: a measured $x$ far above $y$, or a percent-level $x$ with the same sign as $y$, would contradict the estimate.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the phase-space model of SU(3) breaking that generates the input width difference $y$ and its estimates for the $PP$ and $4P$ representations."},{"cited_title":"Global duality in heavy flavor hadronic decays","cited_arxiv_id":"hep-ph/0112323","evidence_quote":"Originates the technique of deriving a dispersion relation in the heavy quark mass within heavy quark effective theory, which this paper adapts."},{"cited_title":"Cambridge Monogr. Part. Phys . Nucl. Phys. Cosmol. 10, 1","cited_arxiv_id":null,"evidence_quote":"Provides the heavy quark effective theory conventions and the mass-independent state normalization used to define the correlator $Sigma_v$."},{"cited_title":"Heavy-to-Light Form Factors in the Final Hadron Large Energy Limit of QCD","cited_arxiv_id":"hep-ph/9812358","evidence_quote":"Gives the large-mass scaling of the heavy-to-light form factor, used to fix the high-energy behavior of $g_{PP}(E)$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier dispersion-relation treatments of the D mass difference, whose spectral density lacked a physical interpretation for $s eq m_D^2$; the contrast defines the novelty here."},{"cited_title":"Limits on D0-D0bar Mixing and CP Violation from the Ratio of Lifetimes for Decay to K-pi+, K-K+, and pi-pi+","cited_arxiv_id":"hep-ex/0306003","evidence_quote":"Current experimental bounds on $x$ and $y$ at the few-percent level, which set the scale against which the estimate is judged."}],"lane_accuracy":{"status":"not yet measured","lane":{"model":"selfhosted/deepseek-v4.1-flash","input_kind":null},"url":"https://pith.science/accuracy"},"review_version":1}