{"id":"b4eecb0d-5d26-4981-890e-e68dc384fe5e","arxiv_id":"2412.10775","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"Salpeter wave functions with S, P, and D components predict the D*0 width as 54.0 keV and show the radiative decay is dominated by relativistic corrections.","lead":"This paper computes the strong and electromagnetic decay widths of the D* meson using a relativistic Bethe-Salpeter wave function, predicting the unmeasured D*0 total width to be about 54 keV. A generalist might care because the result is a precise, testable prediction from quark models, and the paper claims that relativistic wave function components, not just the leading S-wave, control the radiative decay rate.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The predicted D*0 widths and the claimed dominance of relativistic corrections in D*→Dγ rest on small P/D-wave components of the Salpeter wave functions that are not constrained by the mass fit and that numerical cancellation amplifies; no sensitivity or error analysis is provided.","rationale":"The reader's CONDITIONAL verdict is appropriate. The strong-decay ratios Γ(D*+→D0π+)/Γ(D*+→D+π0) = 2.16 and Γ(D*+→D0π+)/Γ(D*0→D0π0) = 1.61 agree well with other models and experiment (Table I), suggesting the computation is not wildly off. However, the absolute widths, especially D*0→D0γ = 19.4 keV, are driven by the small P/D waves. Because the mass fit does not constrain these waves and the paper gives no error estimate, the central claim is not fully supported. The strongest check would be a sensitivity analysis or an independent extraction of the P/D admixtures. The abstract's claim of non-relativistic dominance in D*→Dπ is contradicted by Table II and should be corrected, but this is an interpretive flaw rather than a numerical one. Overall, the paper is plausible but requires a robustness check before the results can be accepted without qualification.","tokens_in":11289,"tokens_out":8760,"duration_ms":77694,"concrete_test":"Perform a sensitivity analysis by recomputing the P/D-wave components from the Salpeter equation with the same kernel but vary the vector-potential scale Λ_QCD between 0.15 and 0.21 GeV and the quark masses by ±5%, then recompute Γ(D*0→D0π0) and Γ(D*0→D0γ) using Eqs. (3) and (5). If the widths change by more than ~20% or if the P-wave admixture changes by more than a factor of 2, the central results are not robust to the unconstrained small components. As a sharper check, extract the D and D* P-wave admixtures from lattice QCD correlation functions (or from a different relativistic quark model) and compare them with the model's prediction; a mismatch would invalidate the destructive-interference mechanism that produces the 34.6 and 19.4 keV widths.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is that the small relativistic P- and D-wave components of the D and D* Salpeter wave functions are quantitatively reliable. These components enter the transition amplitudes via Eq. (3) (strong) and Eq. (5) (EM). The kernel parameters (λ, α_s, quark masses) are fit only to meson masses (Sec. III), so the P/D components are unconstrained predictions. Table II shows that for D*0→D0π0 the pure S-wave contribution (S×S′) is 91.0 keV, while the complete width is 34.6 keV; the reduction comes from destructive interference involving P-wave terms. Table V shows the same pattern for D*0→D0γ: S×S′ = 2.60 keV, complete = 19.4 keV, so relativistic corrections dominate. A modest change in the magnitude or phase of the P-wave amplitudes (B1,B2,B7,B8 in Eq. (10) and A3,A4 in Eq. (12)) could therefore shift the predicted widths by a factor of several. The paper provides no error bars, no sensitivity study, and no comparison of these small components with any independent determination. In addition, the abstract's statement that 'in a strong decay D*→Dπ, the non-relativistic contribution is dominant' is internally inconsistent with Table II, where the non-relativistic S×S′ term is 91.0 keV but the full width is 34.6 keV, i.e., the relativistic correction is large and destructive.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript solves the complete Salpeter equation for the D and D* mesons and uses the resulting relativistic wave functions, which contain S-, P-, and D-wave components, to compute the strong decays D* → Dπ and the radiative decays D* → Dγ. The authors report Γ(D*(2007)0 → D0π0) = 34.6 keV and Γ(D*(2007)0 → D0γ) = 19.4 keV, leading to an estimated D*0 full width of 54.0 keV, and they emphasize that the radiative decay is an M1+E2+M3+E4 transition in which relativistic corrections dominate, whereas they state that in the strong decay the non-relativistic contribution dominates.","tokens_in":11639,"tokens_out":2926,"duration_ms":27505,"significance":"If the predictions are reliable, the paper provides useful results for the unmeasured D*0 total width and for the branching ratios of D*0, and it highlights the possible importance of relativistic partial-wave admixtures in heavy-light meson decays. The strong decay widths for D*(2010)+ agree well with experiment (55.8 vs 56.5 keV and 25.8 vs 25.6 keV), and the predicted branching ratios for D*0 agree with the PDG values; the use of ratios between decay widths also helps reduce model dependence. The main limitations are the absence of uncertainty estimates and the reliance on small P- and D-wave components that are not constrained by the meson masses used to fix the kernel parameters.","major_comments":[{"comment":"The statement that in the strong decay D* → Dπ the non-relativistic contribution is dominant is contradicted by the authors' own Table II. For D*0 → D0π0 the pure S×S' term gives 91.0 keV, while the complete width is 34.6 keV; the corresponding numbers for D*+ → D0π+ in Table III are 143 keV and 55.8 keV. Thus the relativistic corrections are large and destructive rather than small, and the abstract and Sec. IV should be revised to reflect this.","section":"Abstract, Sec. III.A, Table II, Sec. IV"},{"comment":"The central prediction Γ(D*0 → D0γ) = 19.4 keV hinges on the small P- and D-wave components of the Salpeter wave functions, since the non-relativistic S×S' contribution is only 2.60 keV while the complete result is 19.4 keV. These components are not constrained by the meson masses used in the Sec. III kernel fit, and no sensitivity study or error estimate is given. Because the final width is dominated by interference among terms of different partial waves, a modest change in the magnitudes or phases of the B1, B2, B7, B8 and A3, A4 amplitudes could shift the predicted width by a factor of several. The authors should provide quantitative sensitivity information or independent constraints on these components.","section":"Sec. III.B, Eq. (10), Eq. (12), Tables V-VI"},{"comment":"The manuscript does not report the numerical values of two parameters that enter the calculation: the decay constant fπ used in Eq. (2) and the constant V0 in the linear potential of the kernel. Since these values are needed to reproduce the results, their omission makes the calculation incomplete as presented.","section":"Sec. II.A, Eq. (2), Sec. III"},{"comment":"No uncertainty estimates are provided for any computed decay width. Given that the theoretical inputs (quark masses, kernel parameters, and the small partial-wave amplitudes) carry uncertainties, the agreement claimed with experiment cannot be fully assessed. At minimum, the authors should quantify the sensitivity of the widths to reasonable variations of the input parameters.","section":"Throughout, Tables I, IV-VI"}],"minor_comments":[{"comment":"There is a typo in the sentence 'and the quark masses mu = 0.374 GeV, md = 0.38 GeV, and mc = 1.62 GeV are usde', where 'usde' should be 'used'.","section":"Sec. III, paragraph 1"},{"comment":"The heading 'DISSCUSSION' should be 'DISCUSSION'.","section":"Sec. IV heading"},{"comment":"The comparison with '68±17 keV in Ref.[ ? ]' contains an unresolved citation placeholder; the reference should be filled in.","section":"Sec. III.B, Eq. (18)"},{"comment":"The table captions use the notation 'D∗0' and 'D0π0' but do not define the prime on the final-state partial waves; a short definition would improve readability.","section":"Tables II-III"}],"recommendation":"major_revision","confidential_remarks":"The paper relies heavily on the authors' own prior Salpeter-equation framework, which is appropriate but means that the novelty is mainly in the application to D* decays. The main scientific issue is that the headline radiative width is dominated by contributions from wave-function components that are not pinned down by the mass fit, so the missing sensitivity analysis is essential rather than cosmetic."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Plainly: this is a workmanlike Salpeter-equation calculation of D* decays that lands close to experiment for the charged strong widths and predicts a D*0 total width of 54 keV with a branching ratio matching PDG. The genuinely new piece is the partial-wave decomposition: the radiative decay is shown to be M1+E2+M3+E4, with the relativistic P/D components supplying the higher multipoles, and in the strong decay a large destructive interference cuts the naive S-wave width from 91 keV down to 34.6 keV. If those small wave-function components are right, this is a useful, testable result.\n\nThat 'if' is the main soft spot. The P- and D-wave components are not constrained by the mass fit that sets the kernel parameters, and the paper gives no sensitivity study or error estimate. Tables II and V show the final widths emerge from a cancellation (strong) or an enhancement (EM) that depends on these small components; a modest change in their magnitude or phase could shift the widths by a factor of a few. The numeric values of V0 and f_pi are not given, and there is no code or data release, so reproducing the calculation is harder than it should be.\n\nThe abstract's claim that in D*→Dπ the non-relativistic contribution is dominant is contradicted by the paper's own numbers: Table II has S×S' = 91.0 keV versus a complete width of 34.6 keV, so the relativistic correction is large and destructive. The text in Sec. III.A actually says the relativistic effect is significant; the abstract and conclusion need correction.\n\nWhat the paper does well: the strong widths for D*(2010)+ agree well with experiment (55.8 vs 56.5 keV, 25.8 vs 25.6 keV), which is real evidence the machinery works, and the comparison with other models is fair. The citations are appropriate, including prior light-front results like Jaus that already found a large M1 width. The prediction for the unmeasured D*0 width is a legitimate contribution.\n\nMy take: this deserves a serious referee, not a desk reject. A referee should demand a sensitivity analysis of the P/D components, an uncertainty estimate, the missing parameter values, and a fix to the abstract. With those, the paper would be a solid addition to hadron spectroscopy. I would not cite it in my own work in the next year, but I would bring it to a reading group as a good example of where relativistic quark-model calculations need extra scrutiny.","headline":"Useful Salpeter calculation of D* decays with a testable D*0 width prediction, but the reliability of the small relativistic wave-function components that drive the results is unquantified, and the abstract contradicts the paper's own Table II.","tokens_in":12178,"tokens_out":2481,"would_cite":false,"duration_ms":22253,"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":"This paper claims that the $D^*$ meson's strong and radiative widths are controlled by relativistic partial-wave components in the Salpeter wave function, giving $\\Gamma(D^{*0})=54.0$ keV and turning $D^{*0}\\to D^0\\gamma$ into an…","keywords":["D* meson","Salpeter equation","relativistic wave function","strong decay","radiative decay","M1 transition","multipole expansion","heavy-light meson"],"falsifier":"A precise measurement of the $D^{*}(2007)^0$ full width, together with the photon angular distribution in $D^{*0}\\to D^0\\gamma$, would settle it: a width far from 54 keV or a purely dipole photon distribution would contradict the claimed relativistic dominance.","tokens_in":11090,"feed_emoji":"⚛️","tokens_out":9548,"duration_ms":80343,"temperature":0.7,"pith_summary":"This paper aims to show that a fully relativistic treatment of the charmed vector meson $D^*$ changes both its strong and radiative decay pattern. Solving the complete Salpeter equation, the authors predict $\\Gamma(D^{*0}\\to D^0\\pi^0)=34.6$ keV and $\\Gamma(D^{*0}\\to D^0\\gamma)=19.4$ keV, for a total width around 54 keV. Their central point is that the $D^*$ and $D$ wave functions contain small relativistic P- and D-wave components, so $D^*\\to D\\gamma$ is not a pure M1 transition but a superposition of M1, E2, M3, and E4 multipoles, with the relativistic corrections actually dominating the rate. If right, this means the $D^*$ width is far below the current 2.1 MeV upper limit and that radiative charmed-meson decays are a sensitive test of quark-model wave functions.","feed_headline":"Relativistic corrections, not M1, drive the D* gamma width","feed_subtitle":"Full Salpeter wave functions put D*0 at 54 keV and recast its radiative decay as M1+E2+M3+E4.","key_machinery":"The central object is the complete instantaneous Salpeter equation and its positive-energy wave-function solutions for a $1^-$ vector meson and a $0^-$ pseudoscalar meson. The $1^-$ wave function is expanded in $S$-, $P$-, and $D$-wave parts; the $0^-$ wave function in $S$- and $P$-wave parts. Transition amplitudes are formed by tracing these wave functions with $\\gamma_\\mu\\gamma_5$ for the strong decay (with a pion entering through PCAC) or with $\\gamma_\\mu$ for the electromagnetic decay, with quark charges $Q_1$ and $Q_2$. The numerical kernel is a linear confining potential plus a vector Coulomb potential whose parameters are fit to meson spectroscopy. This machinery lets the authors separate the naive $S\\times S'$ contribution from relativistic $P$- and $D$-wave contributions and expose the interference that sets the final widths.","core_discovery":"The authors solve the complete Salpeter equation for the charmed mesons and use the positive-energy wave functions to compute the $D^*\\to D\\pi$ and $D^*\\to D\\gamma$ transition amplitudes. In their solutions the $D^*$ is not a pure $S$-wave state: its wave function is a mixture of $S$-, $P$-, and $D$-waves, while the $D$ meson mixes $S$- and $P$-waves. Consequently $D^*\\to D\\gamma$ is calculated as an M1+E2+M3+E4 transition rather than a pure M1. They report $\\Gamma(D^{*0}\\to D^0\\pi^0)=34.6$ keV and $\\Gamma(D^{*0}\\to D^0\\gamma)=19.4$ keV, estimating $\\Gamma(D^{*0})\\simeq54.0$ keV, and note that the nonrelativistic M1 piece contributes only 2.60 keV of the radiative width, so the relativistic correction dominates.","pith_inferences":["If the small partial waves are this decisive, then using $D^*\\to D\\gamma$ as a clean magnetic-dipole probe of $D^*$ properties would be misleading; a photon angular-distribution measurement is a direct way to check the multipole content.","The interference pattern that reduces the strong width depends on the relative sign of the $S$- and $P$-wave amplitudes, so an independent lattice-QCD calculation of the $D^*\\to D\\pi$ matrix element could confirm or exclude the kernel's parameter choice.","The same relativistic components are likely to affect other heavy-light vector mesons such as $D_s^*$ and $B^*$; extending the calculation would predict radiative widths that deviate from pure-M1 estimates by comparable factors.","Because the $D^{*0}$ width is predicted so precisely, a future high-statistics measurement would turn the Salpeter-equation wave functions into a quantitative test of the instantaneous approximation itself."],"forward_implications":["The $D^{*}(2007)^0$ total width is predicted to be about 54 keV, far below the current experimental upper limit of 2.1 MeV, so a dedicated measurement is a sharp test of the calculation.","The predicted branching fractions, 64.1% for $D^{*0}\\to D^0\\pi^0$ and 35.9% for $D^{*0}\\to D^0\\gamma$, match the measured values within uncertainties.","The strong-decay width of $D^{*+}\\to D^0\\pi^+$ (55.8 keV) and the ratio to $D^{*+}\\to D^+\\pi^0$ (about 2.16) agree with experiment and with other model calculations, making these quantities a stable benchmark.","For $D^{*+}\\to D^+\\gamma$, the predicted 0.84 keV is close to the measured $1.33\\pm0.4$ keV, and the large neutral-versus-charged radiative width difference is traced to quark charges adding rather than cancelling in the two photon-emission graphs.","In the strong decay $D^{*0}\\to D^0\\pi^0$, the $S\\times S'$ term alone would give 91.0 keV, so relativistic $P$-wave interference is required to bring the width down to 34.6 keV."],"supporting_citations":[{"why":"Supplies the experimental masses, branching ratios, and width limits used to benchmark the calculated strong and radiative widths.","marker":"[13]"},{"why":"Original Bethe-Salpeter equation from which the instantaneous framework descends.","marker":"[26]"},{"why":"Introduces the instantaneous (Salpeter) approximation used to define the positive-energy wave functions.","marker":"[27]"},{"why":"Mandelstam formula connecting the BS wave functions to transition matrix elements.","marker":"[34]"},{"why":"Gives the strong-decay transition amplitude in terms of positive-energy Salpeter wave functions.","marker":"[35]"},{"why":"Provides the interaction kernel (linear plus vector potential) and the recipe for numerically solving the Salpeter equation.","marker":"[36]"},{"why":"Supplies the electromagnetic transition amplitude with photon emission from quark and antiquark.","marker":"[37]"},{"why":"Gives the general relativistic $1^-$ wave function, including the $S$-, $P$-, and $D$-wave decomposition used in the calculation.","marker":"[38]"}],"fun_headline_variants":["Relativistic waves recast D* radiative decay beyond M1","D* width 54 keV: relativistic correction rules the gamma decay","Not pure M1: D* to D gamma has E2, M3, E4 from Salpeter","Salpeter solutions mix D-wave, making relativistic gamma dominant","D*0 gamma width 19.4 keV: relativistic part beats M1"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything hangs on the small P- and D-wave components of the wave functions being quantitatively right, since those components provide the claimed relativistic dominance in $D^*\\to D\\gamma$ and the interference that lowers the strong width.","fun_headline_variants_meta":{"raw":{"variants":["Relativistic waves recast D* radiative decay beyond M1","D* width 54 keV: relativistic correction rules the gamma decay","Not pure M1: D* to D gamma has E2, M3, E4 from Salpeter","Salpeter solutions mix D-wave, making relativistic gamma dominant","D*0 gamma width 19.4 keV: relativistic part beats M1"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000239,"raw_usage":{"total_tokens":1566,"prompt_tokens":1048,"completion_tokens":518,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":664,"completion_tokens_details":{"reasoning_tokens":415}},"tokens_in":664,"tokens_out":518,"duration_ms":5069,"temperature":1.0,"reasoning_tokens":415,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T15:37:19.908196+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A precise measurement of the $D^{*}(2007)^0$ full width, together with the photon angular distribution in $D^{*0}\\to D^0\\gamma$, would settle it: a width far from 54 keV or a purely dipole photon distribution would contradict the claimed relativistic dominance.","supporting_citations":[{"cited_title":"Navas et al ., (Particle Data Group), Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental masses, branching ratios, and width limits used to benchmark the calculated strong and radiative widths."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the instantaneous (Salpeter) approximation used to define the positive-energy wave functions."},{"cited_title":"Wang, G.-L","cited_arxiv_id":null,"evidence_quote":"Mandelstam formula connecting the BS wave functions to transition matrix elements."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the strong-decay transition amplitude in terms of positive-energy Salpeter wave functions."},{"cited_title":"Mandelstam, Proc","cited_arxiv_id":null,"evidence_quote":"Provides the interaction kernel (linear plus vector potential) and the recipe for numerically solving the Salpeter equation."},{"cited_title":"Chang, C","cited_arxiv_id":null,"evidence_quote":"Supplies the electromagnetic transition amplitude with photon emission from quark and antiquark."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the general relativistic $1^-$ wave function, including the $S$-, $P$-, and $D$-wave decomposition used in the calculation."}],"review_version":1}