$B_c$, $B_s$ and $D_s$ production at lepton-hadron colliders

Basic Information

• Paper ID: 2501.00550
• Title: $B_c$, $B_s$ and $D_s$ production at lepton-hadron colliders
• Authors: Qi-Wei Hu, Cong-Feng Qiao, Li-Ping Sun
• Classification: hep-ph (High Energy Physics - Phenomenology)
• Publication Date: December 31, 2024
• Paper Link: https://arxiv.org/abs/2501.00550

Abstract

This study investigates electroproduction processes at lepton-hadron colliders within the Non-Relativistic Quantum Chromodynamics (NRQCD) framework: $e+p\to e+B_c+\overline{c}+b$, $e+p\to e+B_s+\overline{s}+b$, and $e+p\to e+D_s+\overline{c}+s$. The paper presents differential cross sections in $\cos\theta$ and $p_T^2$ at the HERA collider. Results demonstrate that $B_c$ meson production is feasible at the EIC collider but not at HERA. The cross sections for $B_s$ and $D_s$ mesons at HERA are quite significant but exhibit sensitivity to variations in quark mass and renormalization scale.

Research Background and Motivation

Problem Background

1. Importance of Heavy Quarkonium Studies: Heavy quarkonium systems have been a paradigmatic system for QCD research since the discovery of $J/\psi$ in 1974. The low relative velocity of constituent quarks allows the scattering amplitude to be decomposed into hadronization and perturbative components.
2. Uniqueness of $B_c$ Mesons: The $B_c$ meson consists of two different heavy-flavor quarks (bottom and charm), and its decay is entirely controlled by weak interactions, providing a unique opportunity to probe weak quark interactions.
3. Advantages of Electron-Proton Colliders: Compared to pp colliders, ep colliders, despite lower energy, offer higher luminosity and different scientific objectives, particularly the forthcoming EIC and EicC colliders.

Research Motivation

1. Filling Theoretical Gaps: Current theoretical studies lack comprehensive treatment of $B_c$ electroproduction, necessitating valuable estimates for future collider experiments.
2. Validating NRQCD Framework: Assessing the applicability of the NRQCD framework to mesons containing single heavy quarks through calculations of $B_s$ and $D_s$ meson production.
3. Providing Experimental Guidance: Offering theoretical support for the design and data analysis of future high-luminosity collider experiments.

Core Contributions

1. First Systematic Calculation: First systematic calculation within the NRQCD framework of electroproduction of $B_c$, $B_s$, and $D_s$ mesons at ep colliders.
2. Multi-Collider Predictions: Providing production cross section predictions for three colliders: HERA, EIC, and EicC.
3. Differential Distribution Analysis: Detailed $\cos\theta$ and $p_T^2$ differential cross section distributions.
4. Observability Assessment: Evaluating the feasibility of observing these mesons at different colliders.
5. Theoretical Uncertainty Analysis: Detailed analysis of the impact of quark mass and renormalization scale variations on results.

Methodology Details

Task Definition

Investigation of heavy-flavor meson production in electron-proton collisions through virtual photon-gluon interactions:

• Input: Electron and proton beams
• Output: Electron + heavy-flavor meson + quark pair
• Constraints: NRQCD framework, leading-order calculation

Theoretical Framework

NRQCD Factorization

Within the NRQCD framework, the cross section can be factorized as: $d\sigma(ep\to e +M + q + \bar{q}') = \int dx f_{g/p}(x, \mu^2) \times d\sigma(eg \to e+ q\bar{q}'[n] + q + \bar{q}')\langle O_M(n)\rangle$

where:

• $M$ is the final meson
• $\langle O_M(n)\rangle$ is the long-distance matrix element (LDME)
• $n = {}^{2S+1}L_J$ represents quantum numbers
• $f_{g/p}(x, \mu^2)$ is the gluon parton distribution function

Electroproduction Formula

For the electroproduction process, the lepton tensor formula is: $d\sigma(eg \to e+ q\bar{q}'[n] + q + \bar{q}') = \frac{1}{4xP \cdot k}\frac{1}{N_cN_s}\frac{1}{(Q^2)2}L^{\mu\nu}H_{\mu\nu}[n]dPS_4$

Calculation Methods

Phase Space Decomposition

To improve efficiency, the 4-body phase space is decomposed as: $\int dPS_4 = (2\pi)^2\int_{4m^2}^{E^2} dP^2_{M q\bar{q}}dPS_2(e+ g \to e + [M q\bar{q}]) \times \int_{m^2}^{M^2} dP^2_{q\bar{q}}dPS_2([M q\bar{q}] \to M + [q\bar{q}])dPS_2([q\bar{q}] \to q+ \bar{q})$

Projection Operators

For pseudoscalar mesons, the projection operator is: $\frac{1_c}{\sqrt{N_c}} \otimes \frac{1}{2\sqrt{M}}\gamma_5(\not{P} +m_1 +m_2)$

For vector mesons, $\gamma_5$ is replaced by $\not{\epsilon}$.

Feynman Diagram Calculation

• 24 leading-order Feynman diagrams generated using FeynArts
• Analytical amplitude calculations and SU(3) color algebra performed using FeynCalc
• Numerical phase space integration performed using CUBA

Experimental Setup

Parameter Settings

Quark Masses:

• $m_c = 1.5$ GeV
• $m_b = 4.9$ GeV
• $m_s = 0.5$ GeV

Coupling Constants:

• Charm mesons: $n_f = 4$, $\Lambda_{QCD} = 297$ MeV
• Bottom mesons: $n_f = 5$, $\Lambda_{QCD} = 214$ MeV
• Renormalization scale: $\mu = m_T = \sqrt{M^2 + p_T^2}$

Long-Distance Matrix Elements: Determined through wave function values at the origin:

• $|R_{B_c}(0)|^2 = 1.642$ GeV³
• $|R_{B_s}(0)|^2 = 0.299$ GeV³
• $|R_{D_s}(0)|^2 = 0.124$ GeV³

Collider Parameters

HERA:

• Electron energy: 27.5 GeV
• Proton energy: 920 GeV
• Experimental cuts: $p_T > 1$ GeV, $0.3 < z < 0.9$, $2 < Q^2 < 100$ GeV²

EIC:

• Electron energy: 10 GeV
• Proton energy: 275 GeV
• Experimental cuts: $Q^2 > 1$ GeV², $p_T > 1$ GeV, $20 < W < 80$ GeV

Experimental Results

$B_c$ Meson Production

HERA Results

Total cross sections:

• $\sigma_{B_c} = 42.13^{+38.14}_{-14.77}$ fb
• $\sigma_{B_c^*} = 164.17^{+151.34}_{-58.65}$ fb

Results indicate insufficient $B_c$ event numbers detectable at HERA, primarily due to:

1. Bottom meson production cross section being 2-3 orders of magnitude lower than charm mesons
2. Small reconstruction efficiency
3. $B_c^± \to J/\psi\pi^±$ branching ratio of only ~0.5%

EIC Results

Total cross sections:

• $\sigma^{EIC}_{B_c} = 5.44^{+6.36}_{-1.91}$ fb
• $\sigma^{EIC}_{B_c^*} = 22.12^{+25.6}_{-7.8}$ fb

At luminosity $10^{34}$ cm⁻²s⁻¹, approximately 3922-19209 $B_c$ mesons can be produced annually. Considering reconstruction efficiency, this yields approximately 5-28 events/year, indicating that observing $B_c$ at EIC is promising.

$B_s$ and $D_s$ Meson Production

HERA Total Cross Sections

$B_s$ Mesons:

• $\sigma_{B_s} = 0.29^{+0.25}_{-0.09}$ pb
• $\sigma_{B_s^*} = 0.73^{+0.63}_{-0.23}$ pb

$D_s$ Mesons:

• $\sigma_{D_s} = 9.35^{+5.13}_{-3.03}$ pb
• $\sigma_{D_s^*} = 20.78^{+11.74}_{-2.8}$ pb

Event Number Estimates

At integrated luminosity $L = 315$ pb⁻¹:

• $B_s^0$: 48-189 events
• $B_s^*$: 119-476 events
• $D_s^+$: 1575-5100 events
• $D_s^{*+}$: 3465-11440 events

Differential Distribution Characteristics

1. Angular Distribution: $\cos\theta$ distribution concentrated in the direction of the incident proton ($\cos\theta = -1$)
2. Transverse Momentum Distribution: Main contributions from the region $p_T^2 \lesssim 30$ GeV²
3. Vector Mesons: Differential cross section for $B_c^*$ decreases more rapidly than for $B_c$

Related Work

Historical Development

1. Hadronic Production Studies: Beginning with leading-order calculations in the 1990s, primarily focused on S-wave $B_c$ production at Tevatron and LHC
2. P-wave Excited States: Including P-wave $B_c$ hadronic production and color-octet mechanism studies
3. Other Production Mechanisms: Electron-positron collisions, heavy-ion collisions, etc.

Innovations in This Work

1. First Systematic Study: First systematic theoretical investigation of $B_c$ electroproduction at ep colliders
2. Multi-Meson Comparison: Simultaneous study of $B_c$, $B_s$, and $D_s$ to assess NRQCD applicability
3. Future Collider Predictions: Providing predictions for next-generation colliders such as EIC

Conclusions and Discussion

Main Conclusions

1. $B_c$ Observability: Observing $B_c$ at EIC is feasible, but not at HERA
2. Abundant $B_s$/$D_s$ Production: Sufficient $B_s$ and $D_s$ events at HERA for major decay channel analysis
3. Theoretical Uncertainties: Mesons containing light quarks are more sensitive to parameter variations, presenting both challenges and opportunities

Limitations

1. Leading-Order Calculation: LO uncertainties limit the precision of theoretical predictions
2. Parameter Sensitivity: Particularly large sensitivity to light quark masses
3. Rare Decay Channels: More events needed for studying rare decay channels

Future Directions

1. NLO Calculations: Performing next-to-leading-order calculations to reduce theoretical uncertainties
2. Polarization Studies: Investigating polarization distributions at high-luminosity colliders
3. New Physics Exploration: Utilizing clean data to search for new physics signals

In-Depth Evaluation

Strengths

1. Theoretical Completeness: Provides comprehensive theoretical calculations within the NRQCD framework
2. Practical Value: Offers important reference for future experimental design
3. Methodological Innovation: Employs efficient phase space decomposition and numerical calculation methods
4. Result Reliability: Detailed analysis of theoretical uncertainty sources

Limitations

1. Precision Constraints: Only considers leading-order contributions; higher-order corrections may be significant
2. Model Dependence: Results are relatively sensitive to LDME values
3. Experimental Details: Insufficient consideration of actual detector effects in experiments

Impact

1. Academic Contribution: Fills the theoretical gap in electroproduction of heavy-flavor mesons at ep colliders
2. Experimental Guidance: Provides important input for the physics program of new colliders such as EIC
3. Method Generalization: The employed methods can be extended to other heavy-flavor hadron production studies

Applicable Scenarios

1. Collider Design: Provides theoretical support for design optimization of ep colliders
2. Data Analysis: Furnishes theoretical expectations for experimental data analysis
3. New Physics Search: Provides background estimates for utilizing heavy-flavor mesons to search for new physics

References

The paper cites 36 important references covering classical and recent research in heavy quarkonium theory, NRQCD methods, collider physics, and other key fields, providing a solid theoretical foundation for this research.