{"id":"b4106899-3a90-40fb-9981-123c2a260c13","arxiv_id":"2411.10660","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A QCD sum rule calculation with dimension-six condensates gives f_D+ = 203.0 ± 1.5 MeV and, with BESIII data, |Vcd| = 0.227 +0.002 -0.001.","lead":"The authors use QCD sum rules with background field theory to compute the D+ meson decay constant and the resulting leptonic branching fractions. A generalist should care because the same calculation is combined with BESIII data to extract the CKM matrix element |Vcd|, a key Standard Model parameter.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Quoted f_D+ = 203.0 ± 1.5 MeV excludes the dominant sum-rule systematics: the continuum threshold s0 and the dimension-six truncation are fixed without assigned uncertainty, so the error budget is underestimated and the |V_cd| extraction inherits it.","rationale":"The reader's weakest assumption identifies exactly the load-bearing issue: the quoted precision of f_D+ and hence of |V_cd| depends on a single continuum threshold and on the truncation of the OPE at dimension-six, with no systematic uncertainty attached to either. My reading agrees that this is the central soft spot. The paper does provide concrete, falsifiable branching-fraction predictions, which is a genuine strength, but those predictions inherit the same understated theory error. The suggested sensitivity scan would settle whether the concern is quantitative: if f_D+ is stable under threshold variation and a reasonable dimension-eight estimate, then the concern is disarmed; if not, the central claim should be rephrased with a larger theory uncertainty. Because the reader's verdict is already CONDITIONAL and my concern does not move the assessment to a different category, the appropriate verdict remains UNCHANGED.","tokens_in":12814,"tokens_out":8845,"duration_ms":91582,"concrete_test":"Perform a controlled sensitivity scan of Eq. (9): keep all Eq. (13) inputs fixed, vary s0 over [5.0, 5.8] GeV^2, and for each s0 reselect the Borel window using the same three criteria from Ref. [39]; record f_D+ and the window. Then add a representative dimension-eight term (e.g., a vacuum-saturation estimate proportional to <q q>^3 / Lambda_QCD^4 with an O(1) coefficient) and recompute in the same window. If f_D+ moves by more than 1.5 MeV under either variation, the Eq. (15) error is understated and the |V_cd| uncertainty must be inflated; if both shifts remain below 1.5 MeV, the quoted error is robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is Eq. (9) evaluated with s0 = 5.4 GeV^2 and a Borel window selected by the three criteria of Ref. [39] (continuum contribution <35%, dimension-six condensate contribution <5%, and flatness of f_D+). The error quoted in Eq. (15), ±1.5 MeV, is propagated only from the input parameters in Eq. (13); no uncertainty is assigned to s0, to the Borel-window choice, or to the truncation of the OPE at dimension-six. This matters because the continuum-subtracted spectral integral is sizeable and because the dimension-six terms are suppressed only by powers of 1/M^2. A modest change of s0 (the excited-state mass motivating 5.4 GeV^2 is not known to better than tens of MeV) or an O(1) vacuum-saturation estimate of the first omitted dimension-eight condensates can plausibly shift f_D+ by more than 1.5 MeV. The final result is directly affected: |V_cd| from Eq. (1) scales as 1/f_D+, so the claimed few-per-mille uncertainty on |V_cd| is not justified until threshold and truncation systematics are quantified. This is a correctness risk distinct from the tension with FLAG's f_D = 209.0(2.4) MeV and PDG's 212.0(7) MeV; that tension is a symptom of the same understated error budget rather than an independent refutation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses QCD sum rules in the background-field framework to compute the D+ meson decay constant f_D+, including condensates up to dimension-six. The authors obtain f_D+ = 203.0 ± 1.5 MeV (Eq. (15)), use it to predict the branching fractions for D+ → ℓ+νℓ (ℓ = e, μ, τ) in Eq. (16) and Table I, and extract |V_cd| = 0.227^{+0.002}_{-0.001} from the BESIII measurement of B(D+ → μ+νμ). The paper claims agreement with BESIII and older lattice results; the disagreement with newer averages is noted only in passing.","tokens_in":13204,"tokens_out":6904,"duration_ms":61981,"significance":"The paper presents an explicit QCD sum rule expression for f_D+ including dimension-six condensates, which is a nontrivial calculation, and gives predictions for decay widths and branching fractions that are testable at BESIII and future facilities. If the error budget were complete, the result would be a competitive non-lattice determination of f_D+ and a clean CKM extraction. The main strengths are the transparent presentation of Eq. (9), the use of a systematic OPE framework, and the broad comparison with experimental and lattice results. However, the quoted precision relies on an incomplete treatment of sum-rule systematics, and the derivation of the central formula is delegated to earlier work, which reduces the self-containedness of the manuscript.","major_comments":[{"comment":"The uncertainty in the quoted f_D+ = 203.0 ± 1.5 MeV accounts only for the input parameters in Eq. (13). No uncertainty is propagated from the continuum threshold s0 = 5.4 GeV^2, from the Borel window selection via the three criteria of Ref. [39], or from the truncation of the OPE at dimension-six. The threshold is motivated by the mass of D0(2550)0, which is not known to high precision, and the dimension-six terms in Eq. (9) are suppressed only by powers of 1/M^2, so a modest variation of s0 or an estimate of dimension-eight contributions could shift f_D+ by more than 1.5 MeV. Since the final |V_cd| scales as 1/f_D+, the claimed few-per-mille uncertainty on |V_cd| is not justified until these systematics are quantified. Please specify the Borel window used, vary s0 over a range consistent with the excited-state mass uncertainty, and estimate the effect of higher-dimensional condensates (e.g., by adding a dimension-eight term and studying its contribution).","section":"Numerical Analysis, Eqs. (9) and (15)"},{"comment":"Equation (9) is the core sum rule that determines f_D+, but its derivation is not presented in this manuscript; the text only refers to Ref. [35]. Since the paper's central claim is that the calculation includes dimension-six condensates and improves the accuracy of f_D+, the derivation (or a substantial outline, including the treatment of infrared divergences and the matching of the OPE to the hadronic dispersion relation) should be given, at least in an appendix. Without this, the reader cannot verify the claimed accuracy or identify what is new relative to Ref. [35].","section":"Theoretical Framework, Eq. (9)"},{"comment":"The final value f_D+ = 203.0 ± 1.5 MeV is several standard deviations below the modern averages quoted by PDG 2024 (212.0 ± 0.7 MeV) and FLAG 2019 (209.0 ± 2.4 MeV). The paper notes agreement with BESIII'14 and LQCD'13, but does not address the two-to-four-sigma tension with the more recent averages. This discrepancy is not a refutation by itself, but it reinforces the need for the systematic error analysis requested in the first comment. The discussion should be updated with a direct comparison to the current PDG and FLAG averages and a comment on whether the sum rule is expected to reproduce them.","section":"Fig. 4 and Table II"}],"minor_comments":[{"comment":"There are several typos and grammatical errors: 'mestioned' should be 'mentioned', 'Moveover' should be 'Moreover', 'the six-dimensional condensates' should be 'dimension-six condensates', 'Fork state' should be 'Fock state', and 'BEIII' in Table II should be 'BESIII'.","section":"Throughout"},{"comment":"The renormalization scale at which Eq. (15) is evaluated is not stated explicitly. The RGE expressions in Eq. (14) contain µ and µ0, but the final numerical value should specify the scale used (e.g., µ = 2 GeV or µ = m_c), otherwise the result is not fully reproducible.","section":"Numerical Analysis, Eq. (14)"},{"comment":"The BESIII'19 entry reads '0.237 ± 0.024 ± 0.12 ± 0.001', but the Introduction quotes '0.237 ± 0.024 ± 0.012 ± 0.001'; the '0.12' appears to be a typo for '0.012'.","section":"Table II"},{"comment":"Reference [38] (Colangelo and Khodjamirian) lacks publication details; please provide the journal, volume, and year, or indicate that it is a preprint.","section":"References"},{"comment":"The Borel window boundaries are not given numerically. The figures show shaded regions, but the paper should state the M^2 range used to obtain Eq. (15) so that readers can reproduce the quoted central value and uncertainty.","section":"Numerical Analysis, Figs. 2 and 3"}],"recommendation":"major_revision","confidential_remarks":"The analysis relies heavily on input parameters and criteria from the same group's previous papers (Refs. [35–39]). An independent cross-check of the condensate values and the Borel window criteria would strengthen the credibility of the result, but this is not a blocking issue in itself. The main concern remains the incomplete systematic error budget, which directly affects the CKM extraction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague — This is a workmanlike QCD sum rule calculation of f_D+ with dimension-six condensates, and a useful cross-check, but the quoted precision is not supported. The paper does what it says: it uses background field theory to include dimension-six effects, gives the sum rule expression explicitly, and computes branching ratios and |V_cd| from BESIII's B(D+→μν). The central value f_D+ = 203.0 ± 1.5 MeV is within errors of older CLEO/BESIII determinations, so the phenomenology is not off by a mile. The derivation of Eq. (9) is delegated to Ref. [35], but the expression is stated, so a referee can check the numerical work.\n\nThe soft spot is the error budget. The ±1.5 MeV error is propagated only from the input parameters in Eq. (13). The continuum threshold s0 = 5.4 GeV^2 is fixed without an assigned uncertainty, and the Borel window is chosen by three criteria from the same group's Ref. [39]. No uncertainty is assigned to the OPE truncation at dimension-six. Since the continuum-subtracted spectral integral is sizeable and the dimension-six terms are suppressed by powers of 1/M^2, a modest change in s0 or an O(1) estimate of dimension-eight condensates can easily shift f_D+ by more than 1.5 MeV. The stress-test note is correct: the |V_cd| result inherits this because it scales as 1/f_D+. Also, the triple-gluon condensate in Eq. (13) has no error at all, which is another unquantified input. The tension with modern lattice values (FLAG 209.0(2.4), HPQCD/PDG ~212) is a symptom of the same understated error budget, not an independent refutation. The paper should either assign a systematic error to s0 and the Borel window, show the sensitivity to those choices, or discuss why the truncation uncertainty is small.\n\nNone of this makes the paper incoherent. The calculation is standard, the presentation is honest about what was computed, and the comparison tables are useful. It just does not support the few-per-mille error claim. For a QCD sum rule audience, this is a reasonable data point; for the CKM/charm community, lattice QCD remains the precision tool. With a revision that quantifies the missing systematics, it could be a solid contribution. I'd send it to a referee, with the request to focus on the threshold and truncation uncertainties.","headline":"A solid but incremental QCD sum rule calculation for f_D+, whose ±1.5 MeV error ignores the dominant threshold and truncation systematics.","tokens_in":13779,"tokens_out":3447,"would_cite":false,"duration_ms":33958,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["13.25.Hw","11.55.Hx","12.38.Aw","14.40.Be"],"model":"deepseek-v4-flash","headline":"The paper claims that a QCD sum rule extended to dimension-six condensates determines the D+ meson decay constant as $f_{D^+}=203.0\\pm1.5$ MeV, and that combining this with the measured $D^+\\to\\mu^+\\nu_\\mu$ branching fraction yields…","keywords":["QCD sum rules","D+ meson decay constant","dimension-six condensates","background field theory","CKM matrix element |Vcd|","leptonic branching fractions","charm quark","quark-hadron duality"],"falsifier":"A lattice QCD calculation of $f_{D^+}$ with total uncertainty below 1 MeV that returns a central value more than 2 MeV above 203 MeV, or a measurement of $\\mathcal{B}(D^+\\to\\mu^+\\nu_\\mu)$ that, combined with an independently known $|V_{cd}|$, implies $f_{D^+}$ above 207 MeV, would falsify the central claim; agreement within errors would support it.","tokens_in":12584,"feed_emoji":"⚛️","tokens_out":8465,"duration_ms":71516,"temperature":0.7,"pith_summary":"The paper aims to fix the D+ meson decay constant $f_{D^+}$ from QCD sum rules by including condensates up to dimension six, and then to use that constant to extract the Cabibbo-Kobayashi-Maskawa matrix element $|V_{cd}|$ from the measured leptonic branching fraction. The authors argue that their sum rule, built in the background field theory framework, yields $f_{D^+}=203.0\\pm1.5$ MeV, which is below the current lattice-QCD world average and closer to some collider determinations. Using the experimental branching fraction for $D^+\\to\\mu^+\\nu_\\mu$, they obtain $|V_{cd}|=0.227^{+0.002}_{-0.001}$, consistent with recent determinations. If correct, the predicted branching fractions for the electron, muon, and tau channels are $8.260\\times10^{-9}$, $3.508\\times10^{-4}$, and $0.935\\times10^{-3}$, respectively, with the muon and tau values slightly lower than lattice-based expectations. A sympathetic reader would care because an independent, purely QCD-based $f_{D^+}$ would help resolve the current spread among experimental and lattice results and sharpen tests of lepton flavor universality.","feed_headline":"QCD sum rules set D+ decay constant at 203 MeV","feed_subtitle":"Dimension-six terms pin the D+ decay constant at 203.0 MeV; with data they yield |Vcd|=0.227.","key_machinery":"The load-bearing object is the two-point QCD sum rule for $f_{D^+}$, built from the axial-vector current correlation function and evaluated in the background field theory approach, where quantum fluctuations carry the perturbative part and background fields encode vacuum condensates. The sum rule links the Borel-transformed correlation function to the hadronic side through quark-hadron duality, with a single continuum threshold $s_0=5.4~\\mathrm{GeV}^2$ and a Borel window fixed by three criteria: continuum contribution below 35%, dimension-six condensate contribution below 5%, and stability of $f_{D^+}$ within the window. Together, the dimension-six truncation and these criteria determine the quoted central value and error.","core_discovery":"The central claim is that extending the QCD sum rule for the D+ meson decay constant to include dimension-six condensates – the double-quark condensate, the quark-gluon condensate, the double-gluon condensate, and the triple-gluon condensate – produces $f_{D^+}=203.0\\pm1.5$ MeV. From this constant and the standard model formula for the leptonic branching fraction, the paper derives the integrated decay widths and branching fractions for all three charged-lepton channels, and then uses the experimentally measured $\\mathcal{B}(D^+\\to\\mu^+\\nu_\\mu)$ to extract $|V_{cd}|=0.227^{+0.002}_{-0.001}$. The author would state that the predicted branching fractions are consistent with existing measurements but are systematically below the values obtained when $f_{D^+}$ is taken from lattice QCD.","pith_inferences":["The quoted 1.5 MeV uncertainty rests on the single-threshold quark-hadron duality ansatz; a more realistic excited-state model could shift $f_{D^+}$ by more than the stated error.","A precise measurement of the ratio $\\mathcal{B}(D^+\\to\\tau^+\\nu_\\tau)/\\mathcal{B}(D^+\\to\\mu^+\\nu_\\mu)$ would test lepton flavor universality without needing $f_{D^+}$ at all.","The same background-field sum-rule machinery could be carried over to $D_s^+$ and $B^+$ decays, where the dimension-six terms would allow a comparative check of the truncation scheme.","If lattice QCD eventually reaches sub-percent precision on $f_{D^+}$ and confirms a value near 210 MeV, the gap with the sum-rule result would indicate missing higher-dimension contributions rather than experimental error."],"forward_implications":["The leptonic branching fractions of the D+ meson become sharp predictions: $\\mathcal{B}(D^+\\to e^+\\nu_e)=8.260\\times10^{-9}$, $\\mathcal{B}(D^+\\to\\mu^+\\nu_\\mu)=3.508\\times10^{-4}$, and $\\mathcal{B}(D^+\\to\\tau^+\\nu_\\tau)=0.935\\times10^{-3}$.","The extracted $|V_{cd}|=0.227^{+0.002}_{-0.001}$ provides an independent cross-check of CKM unitarity and of semileptonic determinations of the same matrix element.","Because the predicted branching fractions are lower than lattice-QCD-based values, a precise measurement of $D^+\\to\\mu^+\\nu_\\mu$ or $D^+\\to\\tau^+\\nu_\\tau$ can distinguish the two theoretical inputs.","Improved experimental precision on the leptonic channels will directly reduce the uncertainty of the $|V_{cd}|$ extraction, since $f_{D^+}$ is already fixed at the sub-percent level in this approach."],"supporting_citations":[{"why":"supplies the background-field propagator expressions and dimension-six sum-rule terms used to build the correlation function.","marker":"[35]"},{"why":"defines the three Borel-window criteria and the continuum threshold parametrization adopted in the numerical analysis.","marker":"[39]"},{"why":"provides the measured $D^+\\to\\mu^+\\nu_\\mu$ branching fraction from which the paper extracts $|V_{cd}|$.","marker":"[13]"},{"why":"provides the masses, lifetime, and world-average decay constant values used as inputs and comparisons.","marker":"[17]"},{"why":"supplies the values of the quark condensate and higher condensates up to dimension six.","marker":"[36]"},{"why":"supplies updated values for the double-quark and quark-gluon condensates used in the sum rule.","marker":"[37]"},{"why":"establishes the QCD sum rules approach with vacuum condensates that the present calculation builds on.","marker":"[33]"}],"fun_headline_variants":["D+ decay constant hits 203 MeV with dimension-six QCD","QCD sum rules pin f_D+ at 203 MeV, yield |V_cd|=0.227","Dimension-six condensates set D+ decay constant to 203 MeV","D+ meson: decay constant and CKM from leptonic decays","Leptonic D+ decays fix f_D+ and |V_cd| precisely"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central premise is that all higher-energy contributions and condensates beyond dimension six can be absorbed into a single continuum threshold and the known lower-dimension condensates, so that neglecting them changes $f_{D^+}$ by less than the quoted 1.5 MeV.","fun_headline_variants_meta":{"raw":{"variants":["D+ decay constant hits 203 MeV with dimension-six QCD","QCD sum rules pin f_D+ at 203 MeV, yield |V_cd|=0.227","Dimension-six condensates set D+ decay constant to 203 MeV","D+ meson: decay constant and CKM from leptonic decays","Leptonic D+ decays fix f_D+ and |V_cd| precisely"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000262,"raw_usage":{"total_tokens":1756,"prompt_tokens":1262,"completion_tokens":494,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":878,"completion_tokens_details":{"reasoning_tokens":389}},"tokens_in":878,"tokens_out":494,"duration_ms":5149,"temperature":1.0,"reasoning_tokens":389,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T19:27:35.793700+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A lattice QCD calculation of $f_{D^+}$ with total uncertainty below 1 MeV that returns a central value more than 2 MeV above 203 MeV, or a measurement of $\\mathcal{B}(D^+\\to\\mu^+\\nu_\\mu)$ that, combined with an independently known $|V_{cd}|$, implies $f_{D^+}$ above 207 MeV, would falsify the central claim; agreement within errors would support it.","supporting_citations":[{"cited_title":"Charmed and strange pseudoscalar meson decay constants from HISQ simulations","cited_arxiv_id":"1312.0149","evidence_quote":"defines the three Borel-window criteria and the continuum threshold parametrization adopted in the numerical analysis."},{"cited_title":"B Meson Decay Constants Using Relativistic Heavy Quarks","cited_arxiv_id":"2212.10125","evidence_quote":"provides the measured $D^+\\to\\mu^+\\nu_\\mu$ branching fraction from which the paper extracts $|V_{cd}|$."}],"review_version":1}