$K \pi$ scattering as a step towards $B \to K^* \ell^+ \ell^-$ from Lattice QCD

This paper presents the status of an exploratory lattice QCD calculation utilizing a dual heavy-quark strategy and the $1+J\to2$ finite-volume formalism to determine the hadronic matrix elements for the resonant $K\pi$ final state, a crucial step toward achieving precise Standard Model predictions for the $B\to K^*\ell^+\ell^-$ decay.

The Gist

Imagine the universe as a giant, incredibly complex LEGO set. Physicists are the master builders trying to figure out exactly how the pieces snap together. For decades, they've had a perfect blueprint called the "Standard Model," but sometimes, when they look at how certain heavy LEGO bricks (particles called B-mesons) break apart, the pieces don't fit quite right according to the blueprint. This suggests there might be hidden, invisible forces or new types of LEGO pieces we haven't discovered yet.

This paper is a progress report from a team of scientists (led by Felix Erben) who are trying to build a super-accurate simulation of one of these tricky break-ups on a computer. Their goal? To solve a specific puzzle: How does a heavy B-meson decay into a K-star particle and two leptons?

Here is the breakdown of their journey, explained with everyday analogies:

1. The Problem: The "Unstable" LEGO Brick

In the real world, when a B-meson decays, it often turns into a K-star ($K^*$) particle. The problem is that the K-star is like a wobbly, unstable LEGO tower. It doesn't sit still; it immediately falls apart into a Kaon and a Pion ($K\pi$).

• The Old Way: Previous computer simulations treated the K-star like a solid, stable brick. They ignored the fact that it's actually a wobbly tower falling apart. This worked fine for very stable particles, but for the wobbly K-star, it's like trying to measure the weight of a collapsing sandcastle by pretending it's a solid rock. The results are too rough to spot the tiny cracks in the Standard Model.
• The New Goal: These scientists want to simulate the falling apart process itself. They need to understand the "dance" between the Kaon and the Pion as they fly apart, which requires a much more sophisticated computer model.

2. The Challenge: The "Goldilocks" Dilemma

Simulating this on a computer is incredibly hard because of two conflicting requirements:

• The Room Size: To see the wobbly K-star dance properly, the computer needs a huge "room" (a large simulation volume) so the particles don't bump into the walls.
• The Resolution: To simulate the heavy B-meson (which is very heavy and moves fast), the computer needs a very high-resolution grid (tiny pixels) to catch all the details.

Usually, you can't have a huge room and tiny pixels at the same time without the computer crashing or taking a million years to run. It's like trying to film a high-speed race car in a massive stadium using a camera that takes 1000 photos per second, but your hard drive is too small to store the video.

3. The Solution: A "Dual-Engine" Strategy

To solve this, the team is using a clever two-pronged approach, like driving a car with two different engines:

1. The Heavy Engine: They use a special method to simulate the actual heavy B-meson directly.
2. The Light Engine: They also simulate lighter versions of the heavy quark (like the charm quark) and then mathematically "stretch" the results to guess what the heavy B-meson would do.

By comparing the results from both engines, they can check their math and ensure they aren't making mistakes. It's like weighing a heavy box by weighing it directly, and also by weighing a lighter version of the box and using a formula to calculate the heavy one's weight. If both methods agree, you know the answer is solid.

4. The Tool: "Distillation" (The Magic Filter)

To handle the massive amount of data, they use a technique called Distillation.

• The Analogy: Imagine you are trying to hear a specific conversation in a noisy crowded room. Instead of recording everyone talking, you put on noise-canceling headphones that only let through the voices of the two people you care about.
• In the Paper: They use a mathematical filter to ignore the "noise" of the computer simulation and focus only on the specific particles involved in the decay. This saves massive amounts of computer power.

5. Current Status: The "First Draft"

This paper is a "status report." They haven't finished the movie yet; they've just finished the storyboard and the first few scenes.

• What they did: They successfully set up the simulation and checked that the "wobbly tower" (the K-star) behaves correctly in their virtual room. They showed that their computer code can see the energy levels of the particles correctly.
• What's next: They need to run the simulation on many more computer configurations (like taking more photos to make a clear picture) to reduce the "fuzziness" (statistical errors). Once they have enough data, they will calculate the final numbers that physicists can compare with real-world experiments.

Why Does This Matter?

If their final numbers match the real-world experiments perfectly, it confirms our current understanding of the universe. But if there is a mismatch—even a tiny one—it could be the smoking gun that proves New Physics exists. It could mean there are invisible particles or forces interacting with the B-meson that we haven't discovered yet.

In short: These scientists are building a super-precise, high-definition movie of a particle breaking apart. They are fixing the camera angles, tuning the focus, and checking the lighting. Once the movie is finished, it will tell us if the universe is exactly as we thought, or if there's a secret plot twist waiting to be discovered.

Technical Summary

1. Problem Statement and Motivation

The rare decay $B \to K^* \ell^+ \ell^-$ is a critical observable for testing the Standard Model (SM) and searching for New Physics. However, precise theoretical predictions are hindered by the complexity of the hadronic final state:

• Resonant Nature: The $K^*$ meson is a broad resonance ($K\pi$ system), not a stable particle. Treating it as a stable particle (narrow-width approximation) is insufficient for the precision required to distinguish New Physics signals.
• Dual Challenges: The calculation requires handling:
  1. Multi-hadron dynamics: A rigorous treatment of the $K\pi$ scattering and the $K^*$ resonance in a finite volume.
  2. Heavy-quark dynamics: Precise control over the $b$-quark, which requires small lattice spacings to suppress discretization errors, while simultaneously requiring large physical volumes to suppress finite-volume effects for the light $K\pi$ system.
• Current Limitations: Previous lattice calculations for vector final states (e.g., $K^*$) are over a decade old and relied on the narrow-width approximation. Recent advances in finite-volume formalism ($1+J \to 2$) and $K\pi$ scattering exist separately but have not yet been combined for heavy-meson decays into resonances.

2. Methodology

The authors present a new exploratory lattice QCD calculation designed to bridge these gaps using a unified framework.

A. Theoretical Framework

• Finite-Volume Formalism: The calculation utilizes the $1+J \to 2$ finite-volume formalism (an extension of Lüscher's method). This relates finite-volume matrix elements to infinite-volume physical amplitudes via a generalized Lellouch-Lüscher factor, which depends on the $K\pi$ scattering amplitude.
• Variational Method: To isolate the $K^*$ resonance, a variational basis of interpolating operators is constructed, including both quark-antiquark ($\bar{s}\gamma_\mu b$) and two-hadron ($K\pi$) operators. Solving the Generalized Eigenvalue Problem (GEVP) yields optimized interpolators for finite-volume energy eigenstates.
• Form Factors: The hadronic matrix elements for vector, axial, and tensor currents are parameterized by seven independent form factors ($V, A_{0,1,2}, T_{1,2,3}$) as functions of momentum transfer $q^2$.

B. Computational Strategy

• Ensemble: The study uses the RBC/UKQCD domain-wall fermion ensemble F1M with:
  • Inverse lattice spacing: $a^{-1} \approx 2.7$ GeV.
  • Light/strange masses: $M_\pi = 232$ MeV, $M_K = 510$ MeV.
  • Lattice size: $48^3 \times 96$.
• Dual Heavy-Quark Strategy: To address the heavy-quark challenge, two approaches are employed on the same ensemble:
  1. Relativistic Heavy Quark (RHQ): Tuned directly to the physical $b$-quark mass.
  2. Domain-Wall Heavy Masses: Three masses in the range $m_c \le m_h \lesssim 0.5 m_b$ to provide a lever arm for interpolation/extrapolation from charm to bottom.
• Distillation Technique: All correlation functions are computed using (stochastic) distillation.
  • Hybrid Setup: Exact distillation is used for source-to-sink lines, while stochastic distillation (with dilution) is used for sink-to-sink lines to manage computational cost.
  • Local Currents: Electroweak currents are inserted locally (unsmeared) to ensure correct operator matching, requiring significant I/O buffering (approx. 500 TB) but enabling a versatile dataset for various transitions ($b \to s, d, u$ and $c \to s, d, u$).

3. Key Contributions

• First Principles Resonance Treatment: This is a pioneering step toward calculating $B \to K^* \ell^+ \ell^-$ form factors where the $K^*$ is treated as a dynamical $K\pi$ resonance rather than a stable particle.
• Unified Heavy-Quark Approach: The simultaneous use of RHQ (physical $b$) and domain-wall heavy masses allows for a controlled study of heavy-quark discretization effects and interpolation.
• Versatile Dataset: The use of distillation creates a unified dataset capable of studying a broad range of heavy-to-light transitions into resonant final states, not limited to just $B \to K^*$.
• Kinematic Focus: The study focuses on the high-$q^2$ region. Due to the lattice size and spacing, reaching low $q^2$ (large recoil) is currently infeasible without encountering inelastic thresholds. High $q^2$ allows for better control of heavy-quark discretization errors.

4. Results and Status

• Two-Point Functions: The authors present preliminary two-point correlation functions for the $K^*$ channel.
  • A $5 \times 5$ correlation matrix was constructed using the operator basis $\{V, K_1\pi_1, K_2\pi_2, K_3\pi_3, K_4\pi_4\}$.
  • Data is available from 8 gauge-field configurations (with 24 source timeslices each, grouped into 3 bins).
  • Observation: The data shows good signal quality, and the resonance region ($M_{K^*} \approx 960$ MeV) is well-separated from the $K\pi$ threshold ($\approx 742$ MeV) on this ensemble, allowing for good resolution.
• Current Limitations: The statistics (8 configurations) are insufficient for a full GEVP analysis and stable spectrum extraction. The uncertainties in the preliminary plots are underestimated.
• Future Work: The target is to increase statistics to 30–60 configurations to stabilize the GEVP extraction and extract the finite-volume spectrum.

5. Significance and Outlook

• Phenomenological Impact: This work lays the groundwork for the first systematically improvable, first-principles calculation of $B \to K^* \ell^+ \ell^-$ form factors. This is essential for reducing theoretical uncertainties in flavor physics anomalies.
• Methodological Advancement: It demonstrates the feasibility of combining multi-hadron scattering formalism with heavy-quark physics on a single lattice ensemble.
• Next Steps:
  1. Complete the GEVP analysis to extract finite-volume energy levels and matrix elements.
  2. Incorporate $K\pi$ scattering input to map finite-volume matrix elements to physical amplitudes.
  3. Extend the calculation to lower $q^2$ (potentially requiring coarser lattices for larger momenta) and address charmonium resonance effects in the intermediate $q^2$ region, which are known to impact phenomenology significantly.

In summary, this paper represents a critical "proof-of-concept" stage in a major lattice QCD program, successfully establishing the technical infrastructure (ensemble, action, and distillation setup) required to tackle one of the most complex problems in heavy-flavor physics: the decay of a $B$ meson into a resonant $K\pi$ final state.