Two nearby states in the $X(3872)$ region: Resolving the radiative-decay ratio tension with $η_{c2}$

To resolve the significant tension between LHCb and BESIII measurements of the radiative-decay ratio in the $X(3872)$ region, this paper proposes a two-state scenario involving a shallow $D^{*0}\bar{D}^0$ bound state and a nearby $2^{-+}$ charmonium candidate ($\eta_{c2}$), which successfully describes the experimental data and predicts specific helicity-angle distributions for future verification.

The Gist

Imagine you are a detective trying to solve a mystery about a very strange, elusive particle called X(3872). For over 20 years, physicists have been studying this particle, but they've hit a wall: the clues don't add up.

Here is the story of the mystery, the conflicting clues, and the new theory proposed in this paper to solve it.

The Mystery: Two Different Stories

The X(3872) is a "heavy" particle made of charm quarks. It's like a tiny, unstable Lego structure that falls apart almost instantly. When it falls apart, it sometimes shoots out a flash of light (a photon) and turns into a different particle.

Physicists measure how often it turns into one specific type of light-flash versus another. They call this the "Radiative-Decay Ratio." Think of it like a baker counting how many chocolate chip cookies vs. oatmeal raisin cookies they sell.

• The LHCb Detective (at the Large Hadron Collider): They looked at the cookies and said, "Hey, we found 1.67 chocolate chip cookies for every oatmeal raisin one!"
• The BESIII Detective (at the Beijing Electron-Positron Collider): They looked at the same type of cookies and said, "No way. We found -0.04 chocolate chip cookies. Basically zero."

The Problem: These two numbers are wildly different. It's like one bakery saying they sell mostly chocolate, and the other saying they sell none. In the world of physics, a difference this big (called a "4.6 sigma" tension) means something is seriously wrong with our understanding. It suggests we aren't looking at just one baker; we might be looking at two different bakers working in the same kitchen, and we've been mixing their orders up.

The New Theory: The "Twin" Hypothesis

The author of this paper, Satoshi Nakamura, proposes a simple solution: There isn't just one X(3872). There are two distinct particles living right next to each other.

Imagine a crowded dance floor where two dancers are spinning so close together that they look like one giant blur.

1. Dancer A (The Known X(3872)): This is the famous one we've known since 2003. It's a "molecule" made of two other particles loosely stuck together (like a D* and a D). It loves to dance to the "oatmeal raisin" tune (decaying into J/ψγ).
2. Dancer B (The Hidden ηc2): This is the new suspect. It's a "charmonium" candidate (a tighter, more compact ball of quarks) with a slightly different spin. It sits just a tiny bit higher in energy than Dancer A. It loves to dance to the "chocolate chip" tune (decaying into ψ'γ).

Why the confusion?

• In the LHCb kitchen (B-decays): The production process favors Dancer B. So, when they count the cookies, they see a lot of chocolate chips (ψ'γ) and think, "Wow, the ratio is high!"
• In the BESIII kitchen (e+e- collisions): The production process favors Dancer A. So, they see almost no chocolate chips and mostly oatmeal raisin.

If you assume there is only one dancer, the math breaks. But if you assume there are two dancers, and they get produced in different amounts depending on the "kitchen" (the experiment), the math works perfectly.

The Evidence: The "Fingerprint" Test

The author didn't just guess; they built a complex computer simulation (a model) to test this "Twin Hypothesis."

1. The Fit: They fed the model all the data from both LHCb and BESIII.

   • Result: The "Two-Dancer" model fit the data like a glove. It explained why the ratios were different and why the shapes of the particle peaks looked the way they did.
   • The Fail: When they tried to fit the data with only one dancer (the old theory), the model failed miserably. It couldn't explain the chocolate chip vs. oatmeal raisin discrepancy.

2. The Missing Piece: The second dancer, the ηc2, has been predicted by theory for years but has never been seen directly. It's the "missing link" in the puzzle. The paper suggests this particle is hiding just above the energy threshold where the X(3872) lives.

The Future: How to Catch the Ghost

So, how do we prove there are two dancers and not just one? The author suggests looking at the dance moves (helicity angles).

• The Analogy: Imagine taking a photo of the dancers. If it's just one dancer, the photo will look the same no matter which angle you shoot from.
• The Prediction: If it's two dancers, the "chocolate chip" moves (ψ'γ) will look very different from the "oatmeal raisin" moves (J/ψγ) when you look at them from different angles.
• The Test: Future experiments can look at these specific angles. If the angles for the two types of light-flashes are different, it proves there are two different particles. If they are the same, the theory is wrong.

Summary

• The Problem: Two major experiments got completely different results about how the X(3872) particle decays.
• The Solution: There are actually two particles in that region, not one. One is the known X(3872), and the other is a hidden, slightly heavier particle called ηc2.
• The Proof: A mathematical model using two particles explains all the confusing data perfectly, while a model with one particle fails.
• The Next Step: Scientists need to look at the "angles" of the decay to confirm the existence of this missing ηc2 particle.

In short, the X(3872) isn't a solo act; it's a duet that we've been hearing as a single voice. Once we separate the voices, the music makes perfect sense.

Technical Summary

1. Problem Statement

The paper addresses a significant discrepancy (tension) in the experimental data regarding the radiative decay ratios of the $X(3872)$ particle.

• The Tension: The LHCb collaboration reported a radiative decay ratio $R_{\psi\gamma} \equiv \mathcal{B}[X(3872) \to \psi'\gamma] / \mathcal{B}[X(3872) \to J/\psi\gamma] = 1.67 \pm 0.25$ from $B^+$ decays. In contrast, the BESIII collaboration reported $R_{\psi\gamma} = -0.04 \pm 0.28$ from $e^+e^-$ annihilation.
• Significance: This represents a $\sim 4.6\sigma$ difference. If $X(3872)$ were a single isolated resonance, these ratios should be consistent across different production mechanisms ($B$-decay vs. $e^+e^-$), assuming experimental systematics are under control. The large discrepancy suggests that the assumption of a single dominant state in the $X(3872)$ region is flawed.
• Secondary Issues: Previous analyses hinted at mass shifts in $X(3872)$ when reconstructed via different channels (e.g., $J/\psi\omega$ and $D^{*0}\bar{D}^0$ vs. $J/\psi\pi^+\pi^-$), and the spin-parity assignment ($J^{PC}=1^{++}$) has been challenged by some $J/\psi\omega$ data.

2. Methodology

The author proposes a two-state scenario and develops a comprehensive coupled-channel model to test this hypothesis against a broad set of experimental data.

• The Two-State Hypothesis:

  1. $X_A$ (The established $X(3872)$): A shallow $D^{*0}\bar{D}^0$ bound state with $J^{PC}=1^{++}$, interpreted primarily as a hadronic molecule with a small $c\bar{c}$ admixture.
  2. $X_B$ (The $\eta_{c2}$ candidate): A $2^{-+}$ charmonium state ($\eta_{c2}$) located slightly above the $D^{*0}\bar{D}^0$ threshold. Theoretical models (quark models, Lattice QCD) predict such a state, but it has not been directly observed.
  • Key Mechanism: The hypothesis posits that $\psi'\gamma$ events arise predominantly from $X_A$ decays, while $J/\psi\gamma$ events receive significant contributions from $X_B$ ($\eta_{c2}$). Furthermore, the production rate of $\eta_{c2}$ relative to $X(3872)$ is significantly enhanced in $e^+e^-$ annihilations compared to $B$-meson decays.

• Theoretical Framework:

  • Coupled-Channel Dynamics: The model treats the system as a coupled-channel scattering problem involving $\{D^{*0}\bar{D}^0\}$, $\{D^{*+}D^-\}$, $J/\psi\omega$, $J/\psi\rho$, $J/\psi\gamma$, and $\psi'\gamma$.
  • Amplitudes:
    • $X(3872)$ is generated dynamically as a pole in the coupled-channel system via short-range contact interactions and rescattering.
    • $\eta_{c2}$ and $\chi_{c0}(3915)$ are included as explicit Breit-Wigner excitations.
    • Weak decays ($B \to K X$) and electromagnetic decays ($e^+e^- \to \gamma X$) are modeled using effective vertices with form factors to account for finite size effects.
  • Isospin Breaking: The model explicitly accounts for isospin violation arising from the mass difference between neutral ($D^{*0}\bar{D}^0$) and charged ($D^{*+}D^-$) thresholds, which is crucial for explaining the $X(3872)$ lineshape.
  • Fit Strategy: The author performs a global fit to invariant mass distributions, branching ratios, and decay ratios from LHCb, BESIII, Belle, and BABAR. Two models are compared:
    1. Model A/B (Default): Includes both $X(3872)$ and $\eta_{c2}$.
    2. Model /$\eta_{c2}$: Excludes the $\eta_{c2}$ component (single-state scenario).

3. Key Contributions

• Resolution of the $R_{\psi\gamma}$ Tension: The paper demonstrates that the single-state model fails to reproduce the observed difference between LHCb and BESIII radiative ratios. The two-state model naturally explains this by attributing the different ratios to the interference and varying production weights of $X(3872)$ and $\eta_{c2}$ in different production environments.
• Unified Description of Lineshapes: The model successfully describes invariant mass distributions for multiple final states ($J/\psi\pi^+\pi^-$, $J/\psi\omega$, $D^{*0}\bar{D}^0$) and explains the slight mass shifts observed in different channels as a result of the interference between the two nearby states.
• Prediction of Helicity-Angle Distributions: The paper provides specific predictions for angular distributions (helicity angles) that differ significantly between the single-state and two-state hypotheses, offering a clear experimental test.

4. Results

• Fit Quality:
  • The two-state models (A and B) provide a consistent description of all datasets, including the conflicting radiative ratios.
  • The single-state model (/$\eta_{c2}$) fails significantly. It predicts $R_{\psi\gamma}^B \approx R_{\psi\gamma}^{e^+e^-}$, which contradicts the experimental data ($1.67$ vs. $-0.04$).
• Pole Positions:
  • $X(3872)$ pole: $3871.4 - 0.022i$ MeV (slightly below threshold).
  • $\eta_{c2}$ pole: $3874. - 0.25i$ MeV (slightly above threshold).
• Branching Ratios: The model reproduces the observed branching fractions for $J/\psi\gamma$ and $\psi'\gamma$ by assigning different decay strengths to the two states. Specifically, $\eta_{c2}$ is predicted to favor $J/\psi\gamma$ over $\psi'\gamma$, while $X(3872)$ favors $\psi'\gamma$ in the context of the fit.
• Helicity Angles (Fig. 4):
  • The model predicts distinct angular distributions for $J/\psi\gamma$ (dominated by $\eta_{c2}$) and $\psi'\gamma$ (dominated by $X(3872)$).
  • For the $D^{*0}\bar{D}^0$ channel, selecting mass regions below and above the threshold ($X_L$ vs. $X_H$) isolates the contributions of $X(3872)$ and $\eta_{c2}$ respectively, leading to distinct angular shapes. A single state would produce flat distributions in both regions.

5. Significance

• Nature of Exotic Hadrons: The work challenges the "single resonance" paradigm for the $X(3872)$ region, suggesting that what is observed as a single peak may be a superposition of a molecular state and a nearby conventional charmonium state.
• Discovery Potential: The paper provides a concrete roadmap for discovering the missing $\eta_{c2}(1D)$ state. It suggests that future experiments (LHCb, BESIII, Belle II) should analyze helicity-angle distributions and specific mass windows to separate the two components.
• Theoretical Consistency: By reconciling the LHCb and BESIII data, the paper resolves a long-standing puzzle in heavy quark spectroscopy and validates the utility of coupled-channel models that include both molecular and compact quarkonium degrees of freedom.

In conclusion, Nakamura's work argues that the $X(3872)$ region is not a single state but a complex interplay between a $D^*\bar{D}$ molecule and a $2^{-+}$ charmonium state ($\eta_{c2}$), a hypothesis that resolves experimental tensions and offers clear, testable predictions for future high-precision measurements.