Accurate B meson and Bottomonium masses and decay constants from the tadpole improved clover ensembles

Using anisotropic clover fermion discretization on 16 $2+1$ flavor QCD ensembles, the authors present a non-perturbative renormalization framework that accurately determines the bottom quark mass, the full S-wave bottomonium spectrum, and decay constants with uncertainties of 0.1% or less, even on lattices with relatively coarse spacing.

The Gist

Imagine the universe is built out of tiny, invisible Lego bricks called quarks. Most of the time, these bricks are light and easy to handle, like standard Lego pieces. But there are two special, super-heavy bricks: the Charm quark and the Bottom quark.

This paper is about the Bottom quark. It's so heavy that trying to study it on a computer is like trying to film a hummingbird's wings with a camera that only takes one picture per second. You miss all the action, and your picture looks blurry.

Here is a simple breakdown of what the scientists in this paper did, using some everyday analogies:

1. The Problem: The "Heavy" Quark is Too Fast

To understand the universe, scientists need to know the exact weight and behavior of these heavy Bottom quarks. They use a method called Lattice QCD, which is like drawing a giant grid (a lattice) over space and time to simulate how particles interact.

• The Issue: The Bottom quark is so heavy and moves so fast that on a standard grid, it "jumps" over the lines. It's like trying to measure the speed of a bullet with a ruler that has inch-long markings. The measurement is too coarse, leading to big errors.
• The Old Way: Previously, scientists had to guess the answer. They would simulate a lighter version of the quark and then use math to "extrapolate" (guess) what the heavy one would do. It's like trying to guess the weight of an elephant by weighing a mouse and multiplying by a huge number. It's risky and prone to error.

2. The Solution: A Specialized "High-Speed Camera"

The team (CLQCD Collaboration) built a new kind of grid specifically for heavy particles.

• The Analogy: Imagine you are trying to film a race.
  • Standard Grid: You take a photo every second. You see the runners, but you miss the details of their footwork.
  • Their New Grid: They used an anisotropic grid. This means they kept the grid square on the ground (space) but made the "time" slices much thinner.
  • The Result: It's like switching to a high-speed camera that takes 1,000 photos per second. Even though the Bottom quark is zooming by, this "time-lens" catches every tiny movement. This allows them to simulate the actual heavy Bottom quark directly, without needing to guess or extrapolate.

3. The Calibration: Tuning the Engine

Even with a high-speed camera, you need to make sure the settings are perfect.

• The Tuning: They adjusted their simulation until the "Bottomonium" (a particle made of a Bottom quark and its anti-particle) weighed exactly what we see in real life (specifically, the mass of the $\Upsilon$ particle).
• The "Renormalization" (Cleaning the Lens): When you look through a new lens, the colors might look slightly off. In physics, the raw numbers from the computer need to be "cleaned up" to match reality. They developed a new, non-perturbative way to clean these numbers, ensuring their results aren't distorted by the grid itself.

4. The Results: A Crystal Clear Picture

Because they could simulate the real thing directly, their results are incredibly precise.

• The Mass: They calculated the mass of the Bottom quark with an uncertainty of less than 1%. It's like weighing a gold bar and being off by less than the weight of a single grain of sand.
• The Decay Constants: They measured how quickly these particles break down (decay). This is crucial because the way particles decay is linked to why the universe has more matter than antimatter (which is why we exist).
• The Comparison: Previous studies were like looking at a blurry photo and trying to guess the details. This paper is like looking at a 4K HD photo. They confirmed previous theories but with much sharper detail and less guesswork.

Why Does This Matter?

The universe is a mystery. One of the biggest mysteries is: Why is there more matter than antimatter? If they were equal, they would have annihilated each other, and we wouldn't be here.

The behavior of the Bottom quark is a key clue in solving this puzzle. By measuring its properties with such high precision, this paper gives scientists a better "rulebook" to test against the Standard Model of physics. If their precise measurements don't match the predictions of our current theories, it could mean there is New Physics waiting to be discovered—perhaps a hidden force or a new particle that explains the universe's existence.

In short: They built a super-fine microscope to look at the heaviest building blocks of the universe, allowing them to measure them with unprecedented accuracy and helping us understand why we are here at all.

Technical Summary

1. Problem Statement

The precise determination of bottom quark physics (masses, decay constants, and mixing parameters) is crucial for testing the Standard Model (SM), particularly in understanding CP violation and the matter-antimatter asymmetry of the universe. However, direct Lattice QCD (LQCD) simulations of the physical bottom quark ($b$) are numerically challenging due to the large mass of the $b$-quark ($m_b \approx 4.2$ GeV).

• Discretization Errors: To keep discretization errors of order $(m_b a)^2$ small, standard isotropic lattices require extremely fine spacings ($a \lesssim 0.02$ fm), which are computationally prohibitive.
• Current Limitations: Most high-precision studies rely on effective field theories (HQET, NRQCD) or extrapolations from lighter quark masses. These approaches introduce additional non-perturbative low-energy constants and systematic uncertainties associated with the extrapolation or expansion.
• Renormalization: Non-perturbative renormalization of heavy-light currents in anisotropic setups has historically relied on lattice perturbation theory, introducing further uncertainties.

2. Methodology

The authors perform a first-principles calculation using 2+1 flavor QCD ensembles generated by the CLQCD collaboration. The key methodological innovations include:

• Anisotropic Clover Fermion Action:

  • They employ an anisotropic clover fermion action for the heavy quark (bottom) while using isotropic gauge ensembles.
  • The action includes a bare anisotropy parameter $\nu$ and tadpole-improved clover coefficients ($c_E, c_B$).
  • This approach utilizes a finer temporal lattice spacing relative to the spatial one, allowing for better resolution of the heavy quark's energy scales without requiring an ultra-fine spatial lattice.
  • The simulations are performed on 16 ensembles spanning 6 lattice spacings ($a \approx 0.037$ fm to $0.105$ fm) with pion masses ranging from 135 to 350 MeV.

• Parameter Tuning:

  • The bare heavy quark mass ($m_Q$) and anisotropy parameter ($\nu$) are tuned non-perturbatively on each ensemble.
  • Constraints: The vector meson mass is matched to the physical $\Upsilon$ mass ($m_\Upsilon \approx 9.46$ GeV), and the effective speed of light is tuned to unity.
  • The study demonstrates that the effective anisotropy parameter approaches unity in the continuum limit with controllable $O(a^2)$ corrections.

• Non-Perturbative Renormalization (NPR):

  • A fully non-perturbative renormalization procedure is developed using the RI/SMOM scheme.
  • They calculate renormalization constants ($Z_V, Z_A, Z_P, Z_S$) for local vector, axial-vector, and pseudoscalar currents.
  • Mixed-Action Strategy: Since the bottom quark uses unsmeared gauge links while light/strange/charm quarks use stout-smeared links, the authors develop a specific strategy to handle mixed-action renormalization constants ($Z_{ns,s}$) in the chiral limit.
  • Ratios of renormalization constants are computed to suppress mass-dependent discretization errors.

• Extrapolation:

  • A combined chiral and continuum extrapolation is performed using a fit ansatz that accounts for valence and sea quark mass deviations and discretization effects ($O(a^2)$ and $O(a^4)$).
  • Calculations are performed directly at the physical $b$-quark mass, avoiding heavy-quark effective theory extrapolations.

3. Key Contributions

• Direct Simulation at Physical Mass: This is one of the first high-precision LQCD calculations to simulate the bottom quark directly at its physical mass on lattices with $a \sim 0.1$ fm (where $m_b a \sim 2.5$), successfully controlling discretization errors to the $\sim 1\%$ level.
• Validation of Anisotropy: The paper provides strong evidence that the effective anisotropy parameter for the heavy quark approaches unity with $O(a^2)$ corrections, validating the use of anisotropic actions on isotropic gauge ensembles.
• Rigorous Renormalization: The development and validation of a fully non-perturbative renormalization framework for heavy-light and heavy-heavy currents in this specific anisotropic setup eliminates reliance on perturbative matching for these quantities.
• Comprehensive Spectrum: The study provides a complete spectrum of S-wave bottom mesons ($B, B_s, B_c, \Upsilon, \eta_b$, etc.) and their decay constants with unprecedented precision.

4. Key Results

The authors report the following results with high precision (uncertainties generally $\le 0.1\%$ for masses and $\sim 1-2\%$ for decay constants):

• Bottom Quark Mass:

  • Determined in the $\overline{\text{MS}}$ scheme: $m_b(m_b) = 4.185(37)$ GeV.
  • This represents the most precise lattice determination to date using a relativistic heavy-quark action without heavy-quark expansions.

• Meson Masses:

  • Full spectrum of S-wave bottom mesons is provided with uncertainties $\le 0.1\%$.
  • Hyperfine Splittings:
    • $\Delta m_{\Upsilon-\eta_b} = 57.8(1.2)(1.0)$ MeV (consistent with PDG and previous high-precision results).
    • $\Delta m_{B^*-B} = 41.9(4.4)$ MeV.
    • $\Delta m_{B_c^*-B_c} = 58.1(1.0)$ MeV (a new precise prediction).

• Decay Constants:

  • Pseudoscalar ($f_P$) and Vector ($f_V$) decay constants for $B, B_s, B_c, \eta_b, \Upsilon$ are provided.
  • Key Values:
    • $f_B = 189.8(3.5)$ MeV.
    • $f_{B_s} = 229.1(2.4)$ MeV.
    • $f_{B_c} = 451.0(5.2)$ MeV.
    • $f_{\eta_b} = 726.4(10.5)$ MeV.
    • $f_{\Upsilon} = 701.4(11.3)$ MeV.
  • The ratio $f_{B_s}/f_B = 1.207(17)$ is slightly larger than previous determinations.
  • Using the experimental value $f_B |V_{ub}| = 0.77 \pm 0.12$ MeV, the authors extract $|V_{ub}| = 4.06(64) \times 10^{-3}$.

• Discretization Control:

  • On the coarsest lattice ($a \approx 0.1$ fm), discretization errors for hadron masses are suppressed to the 1% level.
  • On the finest lattice ($a \approx 0.04$ fm), renormalization constant discretization errors are reduced to $\sim 2\%$, significantly better than previous HISQ calculations at similar spacings.

5. Significance

• Precision Flavor Physics: The results provide critical inputs for interpreting experimental data from LHCb and Belle II, particularly for constraining the CKM unitarity triangle and searching for New Physics via CP violation.
• Methodological Benchmark: The successful application of anisotropic clover fermions with non-perturbative renormalization at the physical $b$-mass sets a new benchmark for heavy-quark lattice simulations. It demonstrates that direct relativistic simulations are feasible and superior to effective theory extrapolations when computational resources allow.
• Reduced Systematics: By working directly at the physical mass and avoiding HQET/NRQCD expansions, the study eliminates a major source of systematic uncertainty present in many previous lattice determinations.
• Future Applications: The framework established here enables the calculation of semileptonic $B$-decay form factors directly in the recoil frame, extending the accessible kinematic range and paving the way for percent-level precision predictions for $B$-physics observables.

In summary, this work represents a major advancement in Lattice QCD, delivering the most precise relativistic determination of bottom quark observables to date and validating a robust methodology for future high-precision heavy-flavor physics.