Correlators of heavy-light quark currents in HQET: Perturbative contribution up to 4 loops and beyond

This paper calculates the perturbative contribution to heavy-light quark current correlators in HQET up to four loops and to all orders in the leading large-$\beta_0$ limit, revealing that renormalon poles appear in the Borel images and that naive nonabelianization performs poorly for these coefficient functions.

The Gist

Imagine the universe is a giant, bustling city made of tiny, invisible Lego bricks. In this city, there are two main types of residents: the Heavy Residents (like the "Heavy Quarks" found in particles called B-mesons) and the Light Residents (the "Light Quarks" that zip around them).

This paper is a massive, high-precision construction manual written by a physicist named Andrey Grozin. It describes how to calculate the "vibe" or "connection" between a Heavy Resident and a Light Resident when they hang out together.

Here is the breakdown of the paper's achievements, translated into everyday language:

1. The Goal: Measuring the "Handshake"

Physicists use something called a Correlator to measure how strongly these two types of particles interact. Think of it like measuring the strength of a handshake between a giant (the heavy quark) and a tiny person (the light quark).

To understand this handshake perfectly, you have to account for the fact that the tiny person isn't perfectly still; they have a little bit of mass (weight). The author calculated this interaction with extreme precision, going up to 4 loops.

• The Analogy: Imagine trying to predict the weather. A 1-loop calculation is like saying, "It might rain." A 2-loop calculation is, "It will rain tomorrow." This paper goes up to a 4-loop calculation, which is like predicting the exact minute it will start raining, how hard it will pour, and exactly how wet your shoes will get. It is the most detailed forecast possible right now.

2. The Toolkit: The "Mathematical Factory"

To do this, the author didn't just use a calculator; they used a massive digital factory.

• qgraf: A machine that draws millions of possible ways the particles could interact (like drawing every possible path a car could take through a city).
• FORM: A super-fast computer program that does the heavy lifting of the math, crunching numbers that would take a human a lifetime to solve.
• LiteRed: A tool that simplifies the messy math, turning a tangled ball of yarn into a neat, straight line.

3. The Big Surprise: The "Naive Guess" Failed

One of the most interesting parts of the paper is a test of a shortcut method physicists often use called "Naive Non-Abelianization."

• The Analogy: Imagine you are trying to guess the price of a house in a new neighborhood. You know the price of houses in a similar neighborhood down the street. So, you just take that price and add a little bit of "neighborhood tax" to guess the new price. This is what "Naive Non-Abelianization" is—it's a clever shortcut based on past experience.
• The Result: The author ran the numbers using the ultra-precise 4-loop method and compared them to the shortcut. The shortcut failed miserably. It was like guessing a house costs $200,000 when it actually costs $2 million. The author notes that for these specific particle interactions, you cannot rely on the shortcut; you have to do the hard, detailed math.

4. The "Renormalons": The Ghosts in the Machine

The paper also looks at something called Renormalons.

• The Analogy: Imagine you are trying to add up an infinite list of numbers to get a total sum. Usually, the numbers get smaller and smaller, and the sum settles down. But sometimes, the numbers start getting huge and wild, making the sum impossible to calculate. These "wild numbers" are the Renormalons.
• The Fix: The author found that these wild numbers (singularities) appear in the math. However, nature is fair. The "wildness" in the calculation of the handshake is exactly canceled out by "wildness" in the vacuum of space (the empty space between particles). It's like two people pushing a car in opposite directions with equal force; the car doesn't move, and the chaos cancels out. This ensures the final result is stable and makes sense.

5. Why Does This Matter?

Why spend years calculating these tiny numbers?

• Unlocking the Secrets of the Universe: These calculations help scientists understand B-mesons, which are particles that decay in very specific ways. By understanding them better, we can test the Standard Model of physics and look for cracks that might reveal "New Physics" (like Dark Matter).
• Checking the Computer Simulations: Scientists use giant supercomputers to simulate the universe (Lattice QCD). This paper provides a "gold standard" answer. If the supercomputer simulation matches this 4-loop calculation, we know the simulation is working correctly. If they don't match, we know something is wrong with the simulation or our understanding of the laws of physics.

Summary

In short, Andrey Grozin built the most detailed map possible of how heavy and light particles interact. He proved that old shortcuts don't work for this specific job and showed that the universe has a built-in "error correction" system (Renormalons) that keeps the math from falling apart. This map is now a vital tool for anyone trying to understand the fundamental building blocks of our reality.

Technical Summary

Here is a detailed technical summary of the paper "Correlators of heavy-light quark currents in HQET: Perturbative contribution up to 4 loops and beyond" by Andrey G. Grozin.

1. Problem Statement

The paper addresses the calculation of the perturbative contribution to the correlator of two Heavy Quark Effective Theory (HQET) heavy-light quark currents. Specifically, the author focuses on:

• Expansion in Light-Quark Masses: The correlator is expanded in powers of light-quark masses ($m$) up to quadratic terms ($m^0, m^1, m^2$).
• High-Order Perturbation Theory: The calculation is performed up to 4 loops (fourth order in the strong coupling constant $\alpha_s$).
• Large-$\beta_0$ Limit: The paper investigates the behavior of these correlators in the limit of a large number of fermion flavors ($n_f \to -\infty$, or large $\beta_0$), which allows for the extraction of terms with the highest powers of $n_f$ to all orders in perturbation theory.
• Renormalon Analysis: The study aims to identify renormalon singularities in the Borel plane and assess the validity of "naive nonabelianization" (a method to approximate full QCD results by replacing $C_F$ with $n_f$ terms) for these specific coefficient functions.

2. Methodology

The author employs a rigorous computational framework combining symbolic algebra and advanced integration techniques:

• Diagram Generation and Calculation:
  • Feynman diagrams are generated using qgraf.
  • Dirac traces and tensor contractions are handled by FORM.
  • Color factors are computed using the color package.
• Integral Reduction:
  • Integration-by-Parts (IBP) reduction is performed using LiteRed (version 2) to reduce the large number of loop integrals to a set of master integrals.
  • Master integrals are expanded in the dimensional regularization parameter $\epsilon$.
• Gauge Invariance Check: The calculation is performed exactly in the covariant gauge parameter $\xi$ (up to 3 loops) and up to linear order in $\xi$ at 4 loops. The cancellation of gauge-dependent terms serves as a primary validation of the results.
• Renormalization Group (RG) Structure:
  • The Wilson coefficients $C_{mn}(\tau)$ and spectral densities $R_{mn}(\omega)$ are analyzed using RG equations involving the anomalous dimensions of the heavy-light current ($\gamma_j$) and the quark mass ($\gamma_m$).
  • The author utilizes known 3-loop results for $\gamma_j$ and $\gamma_m$, and 2-loop results for the QCD $\beta$-function.
• Large-$\beta_0$ Approach:
  • The bare coefficient functions are expressed as a series in $1/\beta_0$.
  • The author derives functions $F_n(\epsilon, u)$ and $S_n(u)$ (Borel images) which encode the all-order behavior in the large-$\beta_0$ limit.
  • Renormalon poles are identified by analyzing the singularities of these Borel functions.

3. Key Contributions and Results

A. Four-Loop Perturbative Results

The paper provides the explicit 4-loop coefficients for the Wilson coefficients of the Operator Product Expansion (OPE) in coordinate space ($\tau$) and momentum space ($\omega$).

• Operators Considered: Dimension 0 (unit operator), Dimension 1 (linear mass $m$), and Dimension 2 (quadratic mass $m^2$ and sum of squared masses $\sum m_i^2$).
• New Coefficients: The paper presents the 3-loop coefficients ($c^{(mn)}_3$ and $r^{(mn)}_3$) which were previously unknown.
  • For the unit operator ($C_1$), the 4-loop expansion in QCD ($n_f=4$) shows a rapidly growing series: $1 + 7.05 \alpha_s/\pi + 10.15 (\alpha_s/\pi)^2 + 125.94 (\alpha_s/\pi)^3$.
  • Similar high-order expansions are provided for $C_m$, $C_{m^2}$, and $C_{\sum m_i^2}$.
• Spectral Densities: The corresponding spectral densities $R(\omega)$ are derived, showing similar high-order coefficients (e.g., $R_1$ has a 3-loop coefficient of $\approx 939$).

B. Large-$\beta_0$ Limit and All-Order Terms

The author calculates the terms proportional to $n_f^l$ to all orders in $\alpha_s$ (up to 8 loops in the expansion of the Borel integral).

• Borel Images: Explicit expressions for the Borel functions $S_n(u)$ are derived for $n=0, 1, 2$.
• Renormalon Poles:
  • UV Renormalons: Found at $u = 1/2$, related to the ambiguity of the HQET residual mass parameter $\bar{\Lambda}$.
  • IR Renormalons:
    • For the unit operator ($n=0$), the leading IR pole is at $u=3$ (dimension-6 operators). The poles at $u=1$ and $u=2$ are absent.
    • For mass-dependent operators ($n=1, 2$), IR poles exist at integer $u \ge 1$.
• Ambiguity Compensation: The paper confirms that IR renormalon ambiguities in the perturbative series are compensated by UV renormalon ambiguities in the vacuum condensates of higher-dimensional operators (e.g., dimension-6 four-quark condensates $\langle (\bar{q}q)^2 \rangle$).

C. Assessment of Naive Nonabelianization

A critical finding of the paper is the failure of naive nonabelianization (setting $\beta_0 = \frac{11}{3}C_A - \frac{4}{3}T_F n_f$) for these specific correlators.

• When the large-$\beta_0$ results are "nonabelianized" to approximate full QCD, the resulting numerical coefficients deviate significantly from the exact 4-loop QCD results.
• Example: For the unit operator $C_1$, the exact 4-loop coefficient is $\approx 125.9$, whereas the naive nonabelianization estimate yields $\approx 290.5$.
• Contrast: The method works reasonably well for the $\sum m_i^2$ contribution but fails for the unit and linear mass terms.

4. Significance

• Precision for QCD Sum Rules: The 4-loop results significantly improve the precision of QCD sum rules used to determine the properties of heavy-light mesons (such as $B$ and $B_s$ mesons).
• Flavor Symmetry Breaking: The inclusion of light-quark mass terms up to quadratic order is essential for studying $SU(3)$ flavor symmetry breaking effects, particularly the mass differences between $B_s$ and $B$ mesons.
• Lattice QCD Comparison: These high-order perturbative results provide a necessary benchmark for comparing with non-perturbative lattice QCD simulations.
• Renormalon Physics: The detailed analysis of renormalon poles clarifies the structure of the OPE in HQET, specifically identifying which condensates are required to cancel perturbative ambiguities.
• Methodological Insight: The paper demonstrates that naive nonabelianization is not a reliable proxy for full QCD calculations in the context of heavy-light current correlators, cautioning against its use for high-precision phenomenology in this sector.

In summary, this work represents a major advancement in the perturbative understanding of HQET correlators, providing the first 4-loop results and a comprehensive analysis of their asymptotic behavior and renormalon structure.