Discovering the $D_0^\ast(2100)$ in $B$ semileptonic decays

Basic Information

• Paper ID: 2509.12133
• Title: Discovering the $D_0^\ast(2100)$ in $B$ semileptonic decays
• Authors: M.-L. Du, F.-K. Guo, C. Hanhart, F. Herren, B. Kubis, R. van Tonder
• Classification: hep-ph (High Energy Physics - Phenomenology), hep-ex (High Energy Physics - Experiment), hep-lat (High Energy Physics - Lattice)
• Publication Date: November 23, 2025 (v2)
• Paper Link: https://arxiv.org/abs/2509.12133

Abstract

The mass and width of the lightest scalar open-charm state $D_0^\ast(2300)$ listed in the Particle Data Group (PDG) exhibit puzzling tension with predictions from unitarized chiral perturbation theory (UChPT) and lattice QCD, which favor a lighter state around 2100 MeV. However, to date, direct experimental evidence for this lighter state remains absent. To facilitate direct observation, this paper introduces angular asymmetries in $B\rightarrow D\pi\ell\nu$ decays, enabling direct extraction of $D\pi$ S-wave phase shifts and discussing new measurement strategies for the Belle II experiment. Sensitivity studies reveal that Belle II can determine the pole position with sufficient precision using existing datasets, thereby establishing the existence of $D_0^\ast(2100)$. Additionally, the statistical requirements for measuring the $D\pi$ isospin-1/2 scattering length are investigated, with precision sufficient to distinguish between UChPT and lattice QCD predictions as well as ALICE measurements using femtoscopy methods.

Research Background and Motivation

Core Problem

This study addresses a long-standing controversy in charm meson spectroscopy: the mass of the lightest scalar open-charm state. The Particle Data Handbook (PDG) lists $D_0^\ast(2300)$ at 2300 MeV, while theoretical predictions from both UChPT and lattice QCD point toward a lighter state around 2100 MeV.

Problem Significance

1. Fundamental Structure Question: Determining this state's mass is crucial for understanding its internal structure—whether it is a quark-antiquark bound state, a compact tetraquark, or a hadronic molecule.
2. Symmetry Tests: Its strange partner $D_{s0}^\ast(2317)$ has nearly identical mass (~2317 MeV), with the mass proximity interpreted as evidence for a compact tetraquark structure, yet contradicting theoretical predictions.
3. Theoretical Consistency: Both UChPT and lattice QCD, two independent methods, predict a lighter mass, suggesting a molecular nature that requires experimental verification.

Limitations of Existing Methods

1. Indirect Evidence: Although $D\pi$ phase shifts extracted from angular momentum analysis of $B^- \rightarrow D^+\pi^-\pi^-$ appear consistent with a 2100 MeV pole, the complexity of the three-body final state (dominated by triangle diagrams) makes conclusions ambiguous.
2. Femtoscopy Controversy: The scattering length measured by the ALICE collaboration using femtoscopy ($a_0^{(1/2)} = 0.02(3)(1)$ fm) is far smaller than theoretical predictions (0.4 fm), exhibiting significant discrepancy.
3. Lack of Direct Measurement: No direct experimental determination of low-energy $D\pi$ scattering parameters and corresponding pole positions exists to date.

Research Motivation

1. Provide independent, direct measurements of $D\pi$ scattering length to verify or challenge ALICE results.
2. Directly extract phase shifts using the Watson final-state interaction theorem in semileptonic decays without model dependence.
3. Design feasible measurement strategies for the Belle II experiment, leveraging existing and future data.

Core Contributions

1. Novel Observable: Introduces four angular asymmetry observables enabling direct extraction of $D\pi$ S-wave phase shifts without relying on model assumptions.
2. Theoretical Framework: Establishes a complete form factor factorization for $B\rightarrow D\pi\ell\nu$ decays, including S-wave and P-wave contributions, correctly handling chiral structure (Adler zero).
3. Measurement Strategy: Designs a practical measurement scheme for Belle II based on hadronic tagging.
4. Sensitivity Studies: Through pseudodata analysis, demonstrates:
   • Existing Belle II data can exclude the PDG average at 2.2σ level
   • Combined Belle+Belle II data achieves 3σ exclusion
   • 2 ab$^{-1}$ data enables 5σ confirmation of $D_0^\ast(2100)$
   • 5-10 ab$^{-1}$ data allows precise scattering length measurement, challenging ALICE results
5. Methodological Innovation: Proposes a modified effective range expansion (ERE) to correctly handle the pole in $q_1\cot\delta_0$ arising from the Adler zero.

Detailed Methodology

Task Definition

Input: Five-fold differential cross section of $B\rightarrow D\pi\ell\nu$ semileptonic decay, depending on:

• $q^2$: invariant mass squared of the lepton pair
• $M_{D\pi}$: invariant mass of the $D\pi$ system
• $\theta$: helicity angle of the D meson
• $\theta_\ell$: helicity angle of the lepton
• $\chi$: azimuthal angle between lepton and hadronic decay planes

Output:

1. $D\pi$ S-wave phase shift $\delta_0(M_{D\pi})$
2. $D_0^\ast$ resonance pole position $M_{\text{pole}} + i\Gamma/2$
3. Isospin-1/2 $D\pi$ scattering length $a_0^{(1/2)}$

Form Factor Factorization Framework

S-Wave Form Factor

Matrix element decomposition (simplified form): $\langle D\pi|A^\mu|B\rangle = P^\mu f_+(q^2, M_{D\pi}^2) + \text{P-wave terms} + \ldots$

Key Innovation: Using unitarity relations $\text{Im} f_{+,i} = \sum_k T_{ik}^* \rho_k f_{+,k}$

Solved as: $f_{+,i} = T_{ij}P_j(q^2)$

Through SU(3) symmetry, relating production amplitudes in three coupled channels ($D\pi$, $D\eta$, $D_sK$): $A(B^- \rightarrow D\pi^{(I=1/2)}) : A(B^- \rightarrow D\eta^{(I=1/2)}) : A(B^- \rightarrow D_s\bar{K}^{(I=1/2)}) = \sqrt{\frac{3}{2}} : \frac{1}{\sqrt{6}} : 1$

Final S-wave form factor: $f_{D^+\pi^-}(q^2, s) = P(q^2)\left[\sqrt{\frac{3}{2}}T_{D\pi\rightarrow D\pi}^{I=1/2} + \frac{1}{\sqrt{6}}T_{D\eta\rightarrow D\pi}^{I=1/2} + T_{D_sK\rightarrow D\pi}^{I=1/2}\right]$

Chiral Behavior: Correctly includes the Adler zero at $s = M_D^2 - M_\pi^2$, a necessary consequence of QCD chiral symmetry.

P-Wave Form Factor

Due to the narrow width of the $D^\ast$ resonance, a factorization approximation is adopted: $F_{1,1}(q^2, M_{D\pi}^2) \approx \hat{F}_{1,1}(q^2)h(M_{D\pi}^2)$

The line shape uses a relativistic Breit-Wigner form, with energy-dependent width expressed through the Chew-Mandelstam function: $h(M_{D\pi}^2) = \frac{g_1 F^{(1)}(M_{D\pi}^2, q_0)}{(M_{D\pi}^2 - M_{R,1}^2) + iM_{R,1}\Gamma_R(M_{D\pi}^2)}$

Dispersion integral: $\Sigma_l(s+i\epsilon) = \frac{s-s_+}{\pi}\int_{s_+}^\infty \frac{\rho(s')n_l^2(s', q_0)}{(s'-s_+)(s'-s-i\epsilon)}ds'$

Phase Shift Extraction Method

Angular Asymmetry Observables

Interference term (S-wave–P-wave) in the differential decay rate: $\frac{1}{K}\frac{d\Gamma^{(5)}}{dq^2 dM_{D\pi}^2 d\chi d\cos\theta d\cos\theta_\ell}\bigg|_{\text{int}} = \sin\theta\sin\theta_\ell\left[\text{4 terms containing}\ f_+f_1^*, f_+g_1^*\right] + \ldots$

Through specific angular integration projection operators, four independent combinations are extracted: $\text{Re}(f_+f_1^*),\ \text{Im}(f_+f_1^*),\ \text{Re}(f_+g_1^*),\ \text{Im}(f_+g_1^*)$

Phase Shift Difference Tangent

The form factor ratio eliminates form factor magnitudes, directly yielding the phase shift difference: $\tan(\delta_0 - \delta_1) = \frac{\text{Im}(f_+f_1^*)}{\text{Re}(f_+f_1^*)} = \frac{\text{Im}(f_+g_1^*)}{\text{Re}(f_+g_1^*)}$

Since the P-wave phase shift $\delta_1$ is completely determined by the known $D^\ast$ (equal to $\pi$ above threshold), this provides a model-independent extraction of the S-wave phase shift $\delta_0$.

Experimental Measurement Scheme

Define one-dimensional asymmetries (practical for Belle II): $A_x(M_{D\pi}) = \frac{N_x^+(M_{D\pi}) - N_x^-(M_{D\pi})}{N_x^+(M_{D\pi}) + N_x^-(M_{D\pi})}$ where $x = \sin\chi, \cos\theta_\ell\cos\chi, \cos\chi, \cos\theta_\ell\sin\chi$

Advantages:

1. Experimental systematic errors cancel in ratios
2. Improved reconstruction efficiency
3. Reduced statistical fluctuations

Parameter Extraction

Pole Position Extraction

Adopting a modified K-matrix form: $T(s) = \frac{K(s)}{1 - \Sigma_0(s+i\epsilon)K(s)}$$K(s) = \frac{E_\pi(g_0 + g_1 s)^2}{s - M_R^2} + E_\pi g_2$

The energy factor $E_\pi$ ensures correct chiral behavior (Adler zero).

Fitting procedure:

1. Minimize likelihood function on pseudodata using Migrad algorithm
2. Calculate confidence intervals using Minos algorithm
3. Determine pole position at confidence interval boundaries

Scattering Length Extraction

Problem: Traditional effective range expansion (ERE) $q_1\cot\delta_0 = \frac{1}{a_0} + \frac{1}{2}r_0 q_1^2 + O(q_1^4)$ has convergence radius limited by the Adler zero (producing a pole in $q_1\cot\delta_0$ just below threshold).

Solution: Modified ERE $q_1\cot\delta_0^{(1/2)} = \frac{c_{-1}}{E_\pi} + c_0 + c_1 E_\pi + O(E_\pi^2)$

Scattering length extracted from threshold value: $a_0^{(1/2)} = \frac{M_\pi}{c_{-1} + c_0 M_\pi + c_1 M_\pi^2}$

This parameterization significantly extends the convergence radius.

Experimental Setup

Data Sources

Benchmark Data: Belle experiment measurement of $B\rightarrow D\pi\ell\nu$ Ref. 48

• Integrated luminosity: 772 fb$^{-1}$ ($\Upsilon(4S)$ resonance)
• Using hadronic tagging method
• $M_{D\pi} \in [2.05, 2.3]$ GeV region: 400±23 events

Extrapolation: Extended to $M_{D\pi} \in [(M_D + M_\pi), 2.3]$ GeV

Study Scenarios

Five luminosity scenarios:

1. Existing Belle II: ~500 fb$^{-1}$ → ~260 events
2. Belle + Belle II: ~1.3 ab$^{-1}$ → ~670 events
3. 2 ab$^{-1}$: Expected 2027 → ~1040 events
4. 5 ab$^{-1}$: Expected 2030 → ~2590 events
5. 10 ab$^{-1}$: Expected 2032 → ~5180 events

Pseudodata Generation

1. Sample events from five-fold differential rate (based on central form factor values)
2. Bin in $M_{D\pi} \in [2.0, 2.3]$ GeV range with 6 bins
3. Calculate $N_x^+$ and $N_x^-$ for each bin
4. Determine statistical uncertainties via bootstrap
5. Calculate asymmetries $A_x$ and take ratios
6. Subtract P-wave phase shift $\delta_1$ to obtain S-wave phase shift data

Fitting Strategy

Minimal Scenario: $g_1 = g_2 = 0$ (retaining only $g_0$)

• 68% confidence region (1σ)
• Using Migrad and Minos algorithms implemented in iminuit

Extended Scenario: Allow $g_1$ or $g_2$ non-zero to assess truncation uncertainties

Experimental Results

Pole Position Determination

Main Results (Figure 4)

One-sigma confidence ellipses in the mass-width plane for each scenario:

Key Findings:

1. Existing Data Already Meaningful: Using only existing Belle II data, one can challenge the PDG's $D_0^\ast(2300)$ at 2.2σ level.
2. Near-term Confirmation: 2 ab$^{-1}$ (2027) enables 5σ confirmation of $D_0^\ast(2100)$.
3. Theory Comparison: Even with 10 ab$^{-1}$, experimental precision remains slightly below theoretical calculations Ref. 28, but any efficiency improvements would make experimental extraction competitive.

Systematic Uncertainty Assessment

Allowing $g_1$ or $g_2$ non-zero:

• Belle II data still excludes PDG value at 2σ level
• High statistics require assessment of K-matrix parameterization truncation uncertainties
• Simultaneously releasing $g_1$ and $g_2$ leads to large correlations but conclusions remain unchanged

Scattering Length Determination

Extraction Results (Table 1)

Theoretical Predictions Comparison:

• UChPT: $a_0^{(1/2)} \approx 0.4$ fm (after unitarization)
• ALICE Femtoscopy: $a_0^{(1/2)} = 0.02(3)(1)$ fm

Key Observations (Figure 5)

Fitting example for 5 ab$^{-1}$ scenario:

• Pseudodata agrees well with fitted curve
• ALICE prediction curve (dashed line) shows significantly different slope near threshold
• Nearly all scenarios can set lower limits at 90% CL, challenging ALICE results

Measurement Challenges

1. High Statistical Requirements: Due to large correlations among ERE coefficients, precise measurement requires large statistics.
2. Low-Statistics Limitations: Existing Belle II data can only set lower limits.
3. Precision Requirements: Distinguishing 0.4 fm (theory) from 0.02 fm (ALICE) requires 5-10 ab$^{-1}$.
4. Comparison with $K\rightarrow\pi\pi\ell\nu$: This classical measurement used two orders of magnitude more events.

Method Verification

Dual-Check Consistency

Two independent phase shift difference ratios: $\tan(\delta_0 - \delta_1) = \frac{\text{Im}(f_+f_1^*)}{\text{Re}(f_+f_1^*)} = \frac{\text{Im}(f_+g_1^*)}{\text{Re}(f_+g_1^*)}$ provide internal consistency checks, enhancing result reliability.

Angular Projection Equivalence

Figures 2 and 3 demonstrate:

• Equivalence between two-dimensional asymmetry categories and simplified one-dimensional categories
• Verification of experimental implementation feasibility

Related Work

Theoretical Predictions

1. UChPT Studies:
   • Kolomeitsev & Lutz (2004): First prediction of molecular nature
   • Guo et al. (2006, 2018): Systematic study of $D_0^\ast$ as molecular state
   • Du et al. (2018): Two-pole structure interpretation, light pole ~2100 MeV, heavy pole ~2300 MeV
   • Albaladejo et al. (2017): UChPT necessarily produces two poles
2. Lattice QCD:
   • Hadron Spectrum Collaboration (2016, 2021): Extrapolation to physical pion mass, supporting ~2100 MeV
   • Guo et al. (2018, 2019): Lattice calculations consistent with UChPT
   • Similarity of scalar and axial-vector charm meson spectra excludes compact diquark structure
3. Flavor Symmetry Analysis:
   • Yeo et al. (2024): Heavy pole belongs to flavor 6 representation, incompatible with $c\bar{q}$ core
   • Supports molecular or diquark-antidiquark structure

Experimental Status

1. Indirect Evidence:
   • LHCb $B^- \rightarrow D^+\pi^-\pi^-$ angular momentum analysis (Du et al. 2021)
   • Phase shifts appear consistent with 2100 MeV, but three-body complexity limits conclusions
2. Femtoscopy Measurement:
   • ALICE (2024): $a_0^{(1/2)} = 0.02(3)(1)$ fm, $a_0^{(3/2)} = 0.01(2)(1)$ fm
   • Severe disagreement with theoretical predictions
   • Methodological controversy (Epelbaum et al. 2025 criticism)
3. Form Factor Measurements:
   • Belle (2023): $B\rightarrow D\pi\ell\nu$ mass spectrum
   • Belle/Belle II (2023, 2024): $B\rightarrow D^\ast\ell\nu$ angular asymmetries
   • Provide benchmark data and methodological precedent for this work

Novel Aspects of This Work

1. vs. Indirect Methods: Provides direct, model-independent phase shift extraction
2. vs. Femtoscopy: Independent verification pathway, not relying on complex correlation function analysis
3. vs. Spectrum Fitting: Angular information provides richer physical constraints
4. Methodological: First systematic proposal to extract scattering parameters from semileptonic decay angular asymmetries

Conclusions and Discussion

Main Conclusions

1. Method Feasibility:
   • Successfully establishes theoretical framework for direct extraction of $D\pi$ S-wave phase shifts from $B\rightarrow D\pi\ell\nu$ angular asymmetries
   • Proposed observables are experimentally measurable with controllable systematic errors
2. Short-term Prospects (2-3 years):
   • Existing Belle II data (~500 fb$^{-1}$) can already challenge PDG's $D_0^\ast(2300)$ at 2.2σ level
   • Combined Belle+Belle II analysis achieves 3σ exclusion
   • 2 ab$^{-1}$ (expected 2027) enables 5σ confirmation of $D_0^\ast(2100)$
3. Long-term Prospects (5-10 years):
   • 5-10 ab$^{-1}$ data enables precise measurement of scattering length $a_0^{(1/2)}$
   • Precision sufficient to distinguish UChPT prediction (~0.4 fm) from ALICE measurement (~0.02 fm)
   • Provides crucial experimental input for understanding low-energy $D\pi$ scattering and charm meson spectrum
4. Theoretical Significance:
   • Confirmation of 2100 MeV pole will support molecular state picture
   • Excludes compact tetraquark structure
   • Validates predictions of UChPT and lattice QCD

Limitations

1. Statistical Requirements:
   • Precise scattering length measurement requires 5-10 years of data accumulation
   • Low statistics only allows setting lower limits, limiting short-term physics conclusions
   • Event samples are two orders of magnitude smaller than classical $K\rightarrow\pi\pi\ell\nu$ measurements
2. Systematic Uncertainties Insufficiently Assessed:
   • Angular resolution and acceptance effects not detailed in simulations
   • Background contributions assumed controllable but not quantified
   • Uncertainty in extrapolating Belle efficiency to Belle II not discussed
   • Form factor $q^2$ dependence uncertainty impact not fully evaluated
3. Simplified Assumptions:
   • Only S-wave and P-wave considered, D-wave neglected (may be important at high $M_{D\pi}$)
   • P-wave form factor factorization approximation (narrow width approximation)
   • SU(3) symmetry breaking corrections not included
   • $q^2$ dependence of production amplitude $P(q^2)$ treated somewhat crudely
4. Insufficient Experimental Details:
   • Complete simulation studies not provided (only based on event number extrapolation)
   • Detector effects, trigger efficiency not discussed
   • Quantitative comparison between semileptonic tagging and non-tagging methods lacking
   • Justification for excluding $B^0$ channel needs more detailed argument
5. Limited Comparison with Other Methods:
   • Detailed comparison with LHCb three-body analysis insufficient
   • Criticism of ALICE femtoscopy method mainly cites others' work
   • Other possible experimental approaches (e.g., production experiments) not discussed

Impact Assessment

Contributions to the Field

1. Short-term (1-3 years):
   • Provides clear analysis targets for Belle II
   • May resolve $D_0^\ast(2300)$ vs. $D_0^\ast(2100)$ controversy
   • Stimulates renewed research interest in charm meson spectroscopy
2. Medium-term (3-5 years):
   • If 2100 MeV pole confirmed, will reshape understanding of positive-parity charm mesons
   • Promotes further development of UChPT and lattice QCD methods
   • May inspire similar analyses in other flavor sectors
3. Long-term (5-10 years):
   • Precise scattering length measurement becomes benchmark data for low-energy QCD
   • Provides crucial experimental support for understanding hadronic molecules
   • Methodology generalizable to other resonance studies

Practical Value

1. Experimental Guidance:
   • Provides concrete measurement strategy for Belle II
   • Sensitivity studies help optimize data collection and analysis priorities
   • Can be directly applied to actual data analysis
2. Theory Testing:
   • Provides experimental criteria for distinguishing different theoretical models
   • Validates or challenges existing theoretical calculations
   • Constrains future theoretical development directions
3. Methodological Value:
   • Demonstrates powerful capabilities of angular analysis in semileptonic decays
   • Provides analysis template for other decay channels
   • Extends application scope of Watson theorem

Reproducibility

1. Theoretical Calculations:
   • Form factor framework clear and reproducible
   • UChPT amplitudes from published work
   • K-matrix parameterization explicit
2. Sensitivity Studies:
   • Pseudodata generation method clearly described
   • Statistical analysis tools specified (iminuit)
   • Specific code not publicly available
3. Experimental Implementation:
   • Measurement strategy described in detail
   • Based on verified techniques
   • Requires actual implementation by Belle II collaboration

Applicable Scenarios

Direct Applications

1. Belle II Experiment:
   • $B\rightarrow D\pi\ell\nu$ analysis (primary target)
   • Similar analysis for $B\rightarrow D^\ast\ell\nu$ (precedent exists)
   • Other semileptonic $B$ decays
2. Other B Factories:
   • Method directly applicable if sufficient statistics available

Generalization Applications

1. Other Flavor Sectors:
   • $\Lambda_b \rightarrow \Lambda_c\pi\ell\nu$ (baryonic)
   • $B_s \rightarrow D_sK\ell\nu$ (strange)
   • $D \rightarrow K\pi\ell\nu$ (charm decay)
2. Other Interactions:
   • Any two-body system producible via semileptonic decay
   • Watson theorem guarantees final-state interaction phase shift extraction
3. Methodological Borrowing:
   • Angular asymmetry construction methods
   • Form factor unitarity constraints
   • Modified ERE handling threshold singularities

Limitation Scenarios

1. Insufficient Statistics: Method requires reasonable event samples (hundreds)
2. Multi-body Final States: Increased complexity requires more refined theoretical treatment
3. Broad Resonances: Narrow-width approximation fails, requiring more complete form factor parameterization
4. Strong Backgrounds: Requires excellent particle identification and signal extraction

Selected References

Key Theoretical Works

• 28 Du et al., PRD 98, 094018 (2018): Two-pole structure interpretation of $D_0^\ast(2300)$
• 37 Liu et al., PRD 87, 014508 (2013): UChPT $D\pi$ scattering amplitude
• 22,23 Hadron Spectrum, JHEP (2016, 2021): Lattice QCD charm meson spectrum
• 36 Guo et al., EPJA 40, 171 (2009): $D\pi$ scattering length LO and NLO calculations

Key Experimental Works

• 48 Belle, PRD 107, 092003 (2023): $B\rightarrow D\pi\ell\nu$ mass spectrum measurement
• 46,47 Belle/Belle II, PRL (2023, 2024): $B\rightarrow D^\ast\ell\nu$ angular asymmetries
• 38 ALICE, PRD 110, 032004 (2024): $D\pi$ femtoscopy scattering length
• 4 NA48/2, EPJC 70, 635 (2010): $K\rightarrow\pi\pi\ell\nu$ phase shifts (methodological precedent)

Methodological Development

• 40 Gustafson et al., PRD 110, L091502 (2024): $B\rightarrow D\pi\ell\nu$ form factor framework
• 41 Herren et al., PRD 112, 014037 (2025): Non-factorizable contribution parameterization
• 44 Hanhart et al., EPJC 84, 483 (2024): Chew-Mandelstam function for unequal masses

Overall Assessment: This is a high-quality theory-experiment combined work addressing an important problem in charm meson spectroscopy with innovative and feasible solutions. The methodology is rigorous, physics motivation clear, and experimental feasibility strong. While precise measurements require extended data accumulation, physically meaningful results can be obtained in the short term (2-3 years). This work promises to advance charm meson spectroscopy research significantly and provides a methodological example for studying strong hadron interactions using semileptonic decays. Main shortcomings include insufficient assessment of systematic uncertainties and incomplete quantification of certain simplified assumptions' impacts, though these do not affect core conclusion validity.