Analysis of the strong decay $X(4140)\rightarrow J/\psi \phi$ via the light-cone QCD sum rules

1. Problem Statement

The nature of the X(4140) resonance, observed in the $J/\psi\phi$ invariant mass spectrum by the CDF, CMS, D0, and LHCb collaborations, remains a subject of intense debate. While its quantum numbers have been determined as $J^{PC} = 1^{++}$, its internal structure is uncertain. Proposed assignments include:

• A $D_s^*\bar{D}_s^*$ molecular state (ruled out for $1^{++}$ by the $D_s^*\bar{D}_s^*$ threshold and quantum number constraints, though $0^{++}$ or $2^{++}$ molecules remain possible).
• A hybrid state.
• A rescattering effect.
• A tetraquark state (specifically a diquark-antidiquark configuration).

Previous QCD sum rule studies successfully predicted the mass of the $X(4140)$ as a $1S$ $[sc]_S[\bar{s}\bar{c}]_A + [sc]_A[\bar{s}\bar{c}]_S$ axialvector tetraquark state. However, a critical test for this assignment is whether the predicted strong decay width ($\Gamma(X(4140) \to J/\psi\phi)$) matches the experimental value. This paper aims to calculate this decay width using Light-Cone QCD Sum Rules (LCQSR) to validate the tetraquark hypothesis.

2. Methodology

The authors employ the Light-Cone QCD Sum Rules (LCQSR) framework to study the two-body strong decay $X(4140) \to J/\psi\phi$.

A. Correlation Function and Interpolating Currents

• Interpolating Currents:
  • $J/\psi$: Standard vector current $J_{J/\psi}^\alpha = \bar{c}\gamma^\alpha c$.
  • $X(4140)$: A specific axialvector tetraquark current representing the $[sc]_S[\bar{s}\bar{c}]_A + [sc]_A[\bar{s}\bar{c}]_S$ structure:
    $$J_X^\beta = \frac{\epsilon_{ijk}\epsilon_{imn}}{\sqrt{2}} \left[ s_j^T C\gamma_5 c_k \bar{s}_m \gamma^\beta C \bar{c}_n^T + s_j^T C\gamma^\beta c_k \bar{s}_m \gamma_5 C \bar{c}_n^T \right]$$
• Two-Point Correlation Function: The authors construct a correlation function involving the $X(4140)$ current and the $J/\psi$ current in the presence of a $\phi$ meson background:
  $$\Pi_{\alpha\beta}(p, q) = i \int d^4x e^{iq\cdot x} \langle \phi(p) | T \{ J_{J/\psi}^\alpha(x) J_X^{\beta\dagger}(0) \} | 0 \rangle$$

B. Phenomenological Side (Hadronic Representation)

• The correlation function is saturated by intermediate hadronic states. The ground state contribution involves the hadronic coupling constant $g_{XJ/\psi\phi}$.
• Novelty in Treatment of Continuum: Unlike traditional LCQSR which often uses a simple continuum subtraction, the authors introduce free parameters ($C_{\psi'}$ and $C_{X'}$) to parameterize the net effects of higher resonances and continuum states. This allows for a more flexible matching between the hadron and QCD sides.
• Rigorous Quark-Hadron Duality: The authors apply a rigorous duality relation suggested in their previous works. They match the hadronic spectral density with the QCD spectral density by integrating over one variable ($s'$) to eliminate the double dispersion relation, leading to a single dispersion relation that is more stable.

C. QCD Side (Operator Product Expansion)

• The correlation function is calculated in the deep Euclidean region using the Operator Product Expansion (OPE).
• The calculation involves contracting quark fields using Wick's theorem and expanding the non-local matrix elements of the $\phi$ meson in terms of Light-Cone Distribution Amplitudes (LCDAs) up to twist-4.
• Contributions include perturbative terms and non-perturbative vacuum condensates (gluon condensate $\langle \frac{\alpha_s}{\pi} G^2 \rangle$).
• The result is expressed as a spectral density $\rho_{QCD}$ integrated over the Borel parameter $T^2$.

D. Sum Rule Derivation

• By equating the phenomenological and QCD sides and applying the Borel transformation, the authors derive the sum rule for the coupling constant $g_{XJ/\psi\phi}$.
• The free parameter $C_{X'}$ is tuned to obtain a stable Borel platform (a region where the result is insensitive to variations in the Borel parameter $T^2$).

3. Key Contributions

1. Decay Width Calculation: This is one of the first rigorous calculations of the strong decay width of the $X(4140)$ specifically assuming the $[sc]_S[\bar{s}\bar{c}]_A + [sc]_A[\bar{s}\bar{c}]_S$ tetraquark structure using LCQSR.
2. Improved Sum Rule Technique: The introduction of parameters $C$ to handle higher resonance/continuum contaminations and the application of rigorous quark-hadron duality for three-point correlation functions represent a methodological refinement over traditional LCQSR approaches, leading to more stable sum rules.
3. Validation of Quantum Numbers: The study provides a dynamical test (decay width) to complement mass predictions, strengthening the case for the $1^{++}$ tetraquark assignment.

4. Results

• Input Parameters: The calculation uses standard PDG values for masses ($m_{J/\psi}, m_\phi, m_X$) and decay constants, along with the running charm quark mass evolved to $\mu = 1.275$ GeV.
• Borel Window: A stable platform is found in the Borel window $T^2 = (2.4 - 3.4) \text{ GeV}^2$.
• Hadronic Coupling Constant:
  $$g_{XJ/\psi\phi} = 2.88 \pm 0.21$$
• Decay Width: Using the calculated coupling constant, the partial decay width is determined to be:
  $$\Gamma(X(4140) \to J/\psi\phi) = 145 \pm 21 \text{ MeV}$$
• Comparison with Experiment: The LHCb collaboration (2021) measured the width as $162 \pm 21^{+24}_{-49}$ MeV. The theoretical result ($145 \pm 21$ MeV) is in excellent agreement with the experimental data within uncertainties.

5. Significance

• Confirmation of Tetraquark Nature: The close agreement between the calculated width and the experimental measurement provides strong support for the hypothesis that the $X(4140)$ is a compact axialvector tetraquark state with the specific diquark-antidiquark structure $[sc]_S[\bar{s}\bar{c}]_A + [sc]_A[\bar{s}\bar{c}]_S$.
• Methodological Benchmark: The successful application of the rigorous quark-hadron duality and the parameterization of continuum effects offers a robust framework for analyzing other exotic hadron candidates (such as $X(4274)$, $X(4500)$, and $X(4700)$) in future studies.
• Resolution of Ambiguity: Since mass calculations alone cannot unambiguously assign the nature of a hadron (as different structures can yield similar masses), the decay width serves as a crucial discriminator. This work demonstrates that the tetraquark model is consistent with both the mass and the decay dynamics of the $X(4140)$.