Real radiative decays of heavy pseudoscalar mesons

This paper presents an ongoing lattice QCD study of radiative leptonic decays for charged pseudoscalar mesons ($D$, $D_s$, $B$, and $B_c$) using a single JLQCD ensemble, aiming to reduce theoretical uncertainties in CKM matrix element extraction and provide first-principles form factor estimates.

The Gist

Imagine the universe is built from tiny, invisible Lego bricks called quarks. When these bricks snap together in specific ways, they form larger structures called mesons (like the $D$, $B$, and $B_c$ particles mentioned in the paper).

This paper is about a team of scientists trying to understand exactly how these heavy Lego structures fall apart and shoot out a flash of light (a photon) while doing so.

Here is the breakdown of their work, translated into everyday language:

1. The Big Picture: The "Sparkle" Decay

Usually, when a heavy particle like a $B$ meson decays, it just splits into a lepton (like an electron or muon) and a neutrino. It's a quiet, invisible breakup.

But sometimes, the breakup is messy. The particle doesn't just split; it gets excited and shoots out a photon (a particle of light) right at the moment of the explosion. This is called a radiative decay.

• The Analogy: Imagine a firework rocket (the heavy meson) flying through the sky. Usually, it just explodes into sparks (the lepton and neutrino). But sometimes, as it explodes, it also shoots out a bright, blinding flare (the photon).
• Why it matters: That flare isn't random. The way the light is emitted tells us about the internal shape and structure of the rocket before it exploded. By studying this light, scientists can measure fundamental constants of the universe with extreme precision.

2. The Problem: We Need a Better Map

Scientists want to use these decays to measure something called the CKM matrix. Think of this as the "ID card" of the universe's particles. It tells us how likely one type of particle is to turn into another.

• Currently, our measurements of these "ID cards" have a little bit of fuzziness (uncertainty).
• The scientists in this paper want to remove that fuzziness. To do that, they need to calculate exactly how the photon is emitted from first principles, using the laws of physics (Quantum Chromodynamics, or QCD) rather than just guessing.

3. The Tool: The "Digital Universe" (Lattice QCD)

You can't build a real $B$ meson in a lab and watch it decay in slow motion with a camera. It happens too fast and is too small.

• The Solution: They build a supercomputer simulation of the universe.
• The Analogy: Imagine a giant 3D grid (like a chessboard, but in 4D space and time). They place their virtual Lego bricks (quarks) on this grid. They then run a simulation to see how the bricks interact, how they move, and how they eventually break apart.
• The Challenge: This grid is huge, and the math is incredibly complex. It requires supercomputers (like the Fugaku and DiRAC systems mentioned) to crunch the numbers.

4. What They Did in This Paper

This paper is a "progress report." They haven't finished the whole puzzle yet, but they have laid the foundation.

• The Setup: They used a specific, high-quality grid (called the JLQCD ensemble) with a very fine resolution (small squares) to make sure their simulation is accurate.
• The Method: They calculated the "three-point function."
  • Analogy: Imagine taking a photo of the meson at the start, a photo of the photon being emitted in the middle, and a photo of the leftovers at the end. They are trying to stitch these three moments together mathematically to see the whole story.
• The Twist: To get enough data points, they used a clever trick called "twisted boundary conditions."
  • Analogy: Imagine trying to measure the wind speed in a room. If you only stand in the corners, you miss the wind in the middle. By "twisting" the rules of the room (the simulation), they can effectively "stand" in the middle and catch more wind data points.

5. The Results So Far

They have successfully run the simulation for the lighter heavy-mesons (the $D$ and $D_s$) and are starting to tackle the heavier ones ($B$ and $B_c$).

• They have generated the raw data (the "photos" of the decay).
• They have verified that their computer code works correctly by checking the signals.
• They are currently working on the next step: filtering out the "noise" (unwanted background signals) to get a clear picture of the photon emission.

6. Why Should You Care?

You might think, "Who cares about a virtual Lego rocket?"

• The Big Reveal: If we can calculate these decays perfectly, we can check if the Standard Model of physics is actually correct.
• The Mystery: If the real-world experiments (like those at CERN or LHCb) disagree with this supercomputer calculation, it means there is New Physics hiding in the shadows. It could mean there are particles we haven't discovered yet!

Summary

This paper is about a team of scientists building a digital microscope to watch heavy particles explode and shoot out light. By doing this with extreme precision, they hope to either confirm our current understanding of the universe or discover a crack in the foundation that leads to new, unknown physics. They are currently calibrating the microscope and taking the first few test shots.

Technical Summary

1. Problem Statement and Motivation

The paper addresses the theoretical calculation of radiative leptonic decays of heavy charged pseudoscalar mesons ($P \to \ell \nu_\ell \gamma$), specifically for the $D$, $D_s$, $B$, and $B_c$ families.

• Physics Goal: These decays are crucial for extracting the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements $|V_{cd}|$ and $|V_{cs}|$ with reduced theoretical uncertainty. Currently, experimental measurements of purely leptonic decays ($P \to \ell \nu$) are "radiation inclusive," meaning they include soft photons below detector thresholds. To match theory with experiment, one must account for radiative corrections and isospin-breaking effects, which requires precise knowledge of the form factors governing the real radiative process.
• Theoretical Gap: While the RM123 collaboration previously calculated form factors for $D$ and $D_s$ mesons, their results were in the electroquenched approximation (neglecting quark-disconnected diagrams). For heavier $B$ and $B_c$ mesons, theoretical knowledge is sparse, relying heavily on Heavy Quark Effective Theory (HQET) where hadronic structure coefficients are poorly constrained.
• Objective: To perform a complete QCD+QED lattice calculation of these decays from first principles, including both quark-connected and (future) quark-disconnected contributions, to provide accurate form factors ($F_V$ and $F_A$) across the physical kinematic range.

2. Methodology

The authors employ Lattice QCD simulations to compute the necessary hadronic matrix elements.

A. Theoretical Framework

The decay amplitude is governed by the time-ordered hadronic tensor $H^{\mu\nu}_{V-A}$, defined as the vacuum expectation value of the product of the electromagnetic current ($J^\mu_{em}$) and the weak $V-A$ current ($J^\nu_{V-A}$).

• Decomposition: The tensor is decomposed into Lorentz structures involving form factors $F_V$ (vector), $F_A$ (axial), and subleading terms $H_1, H_2$.
• On-shell Simplification: For real photons ($k^2=0$), the contributions from $H_1$ and $H_2$ vanish, leaving only $F_V$ and $F_A$ as the non-perturbative inputs required.
• Correlator: The form factors are extracted from the Euclidean three-point correlation function:
  $$C^{\mu\nu}_{3, V-A}(p, k, t_\gamma, t_W) = \int d^3x e^{-ik\cdot x} \langle J^\mu_{em}(x, t_\gamma+t_W) J^\nu_{V-A}(0, t_W) P^\dagger(p, 0) \rangle$$
  where $t_W$ is the source-sink separation.

B. Lattice Setup

• Ensemble: Simulations are performed on a single JLQCD ensemble with lattice spacing $a \approx 0.044$ fm.
• Action:
  • Gauge: Tree-level improved Symanzik action with stout-smeared links.
  • Fermions: Möbius domain-wall fermions (preserving chiral symmetry) with a Wilson-Dirac kernel.
• Quark Masses:
  • Light ($u,d$) and strange ($s$) quarks are tuned to physical values.
  • Heavy Quarks: Due to the coarse lattice spacing ($a \approx 0.044$ fm), direct simulation at physical $B$ and $B_c$ masses is impossible (residual mass effects break chiral symmetry for $am_h \gtrsim 0.7$). Instead, the authors simulate six different heavy-quark masses ($am_h$) and plan to extrapolate to the physical point later.
• Boundary Conditions: To access a dense set of photon momenta ($k$) and cover the physical kinematic region ($0 < x_\gamma < 1 - m_\ell^2/m_P^2$), the authors use Twisted Boundary Conditions (TBC) in addition to standard Fourier modes. This allows them to fill gaps in the $x_\gamma$ spectrum (specifically around $x_\gamma \sim 0.1-0.2$) that pure Fourier modes miss.

C. Computational Strategy

• Diagrams: The study focuses on quark-connected diagrams (Fig. 1a), where the photon is emitted by either valence quark. Quark-disconnected diagrams (Fig. 1b) are currently neglected but planned for future inclusion.
• Sources:
  • Meson interpolators use Jacobi smearing with stochastic $Z_2 \otimes Z_2$ time-wall sources.
  • Weak currents use point sources.
• Data Volume: The initial dataset consists of 18 gauge configurations with 32 time translations per configuration. They compute all non-zero $(\mu, \nu)$ components of the correlator for various source-sink separations ($t_W$) to control excited state contamination.

3. Key Contributions and Results

• Initial Dataset Generation: The authors successfully generated the first lattice data for the three-point correlators $C^{\mu\nu}_{3}$ for $D$ and $D_s$ mesons.
• Signal Demonstration: Figure 4 presents the signal for the vector ($C^V$) and axial ($C^A$) components of the correlators. The data shows clear plateaus for different source-sink separations, validating the methodology.
• Momentum Coverage: By combining Fourier modes ($|n| \le 3$) with a specific twist angle ($\theta \approx 0.2762$), the authors achieved coverage of 369 distinct photon momenta. This significantly improves the resolution of the form factors $F_V(x_\gamma)$ and $F_A(x_\gamma)$ compared to previous studies.
• Heavy Quark Extrapolation Strategy: The paper establishes the framework for simulating heavy mesons at unphysical masses to later perform a controlled extrapolation to the physical $B$ and $B_c$ masses, overcoming the limitations of the current lattice spacing.

4. Significance and Future Outlook

• Precision CKM Physics: This work is a critical step toward reducing theoretical uncertainties in the determination of $|V_{cd}|$ and $|V_{cs}|$. By providing first-principles estimates of radiative corrections, it enables a more rigorous comparison between experimental inclusive rates and theoretical predictions.
• Exploration of $B_c$ and $B$ sectors: Unlike previous lattice studies limited to charm mesons, this framework is designed to extend to the $B$ and $B_c$ sectors, where experimental data is scarce (e.g., $B_c$ radiative decays are currently unexplored).
• Future Steps:
  1. Excited State Control: Implementation of the Laplace filter technique to better isolate the ground state signal.
  2. Disconnected Diagrams: Inclusion of quark-disconnected contributions (Fig. 1b) to achieve a complete QCD+QED calculation.
  3. Physical Point: Performing calculations on multiple ensembles with varying lattice spacings and quark masses to extrapolate results to the physical $B$ and $B_c$ masses.

In summary, this paper presents a robust methodological framework and the first results for a high-precision lattice QCD calculation of radiative leptonic decays, aiming to bridge the gap between theoretical predictions and experimental measurements in heavy flavor physics.