Three-loop functions for a quartic model with a cutoff

✨

This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to bake the perfect cake (a physical model of the universe) using a specific recipe (a mathematical theory called a "quartic model"). You want to know exactly how the flavor changes as you tweak the ingredients. In physics, this "flavor" is how particles interact, and the "ingredients" are things like mass and force.

However, when you try to calculate the flavor of this cake using the laws of quantum mechanics, you run into a problem: the math explodes. It's like trying to measure the height of a mountain that goes up to infinity. To fix this, physicists use a "sieve" or a "cutoff" (a mathematical tool called regularization) to ignore the infinitely high, nonsensical parts of the mountain and focus on the realistic ones.

This paper is about checking the recipe for a very complex, three-layer cake (a three-loop calculation). Here is the breakdown in simple terms:

1. The Problem: Too Many Layers

For a long time, physicists could only perfectly calculate the first two layers of this cake. They knew exactly how the flavor changed at the "one-loop" and "two-loop" stages. But the third layer was a mystery, especially when using a specific type of sieve (called a cutoff regularization) that keeps the math physically realistic (ensuring energy doesn't go negative).

Previous attempts to calculate this third layer used a sieve that was mathematically convenient but physically "broken" (it allowed for impossible negative energies). The authors of this paper wanted to use a better, sturdier sieve (a specific function called $f_4$ ) that guarantees the physics stays sane.

2. The Solution: A New Numerical Map

The authors didn't just guess the answer; they built a massive, high-precision digital map of the "auxiliary integrals."

The Analogy: Imagine you need to measure the volume of a weirdly shaped, invisible jelly. You can't just pour water into it. Instead, you have to break the jelly down into millions of tiny, measurable cubes, calculate the volume of each, and add them up.
What they did: They took the complex math equations for these "cubes" (integrals) and used powerful computers (running Python code) to crunch the numbers. They calculated 13 different "volume measurements" (labeled $\alpha_1$ to $\alpha_{13}$ ) for their new, sturdy sieve.

3. The Results: A New Flavor Profile

Once they had these numbers, they plugged them into the master recipe to find the Beta Function and Anomalous Dimensions.

The Beta Function: Think of this as the "thermometer" for the cake. It tells you how the strength of the interaction changes as you change the temperature (or energy scale).
The Finding: They found that with their new, physically correct sieve, the "thermometer" reads slightly differently than it did with the old, broken sieve or the standard "dimensional" sieve used by most other physicists.
- Old broken sieve result: ~45.91
- New sturdy sieve result: ~45.75
- Standard method result: ~32.55

This proves that the "flavor" of the theory does depend on which sieve you use, but the authors' result is the one that respects the rule of "no negative energy."

4. How They Did It (The Kitchen Tools)

The paper describes their cooking process:

Simplification: They turned 8-dimensional math problems (which are impossible to visualize) into 4-dimensional problems by switching to spherical coordinates (like measuring a ball by its radius instead of x, y, z coordinates).
Fourier Transform: They used a mathematical "magic trick" (Fourier transform) to turn complicated multiplication problems into simpler addition problems.
Numerical Integration: They used a method called the "trapezoid rule" (basically, slicing the area under a curve into trapezoids and adding their areas) to get the final numbers. They did this with extreme precision, ensuring their error margin was tiny (less than one part in 100 million).

Why Does This Matter?

This paper is like a quality control check for theoretical physics.

It confirms that when we use a "physically safe" cutoff (one that doesn't break the laws of energy), we get a specific, calculable result for the three-loop interaction.
It provides a reference point. Now, other scientists can compare their own theories against these numbers to see if their models are working correctly.
It opens the door for even more complex calculations (like four-loop or five-loop cakes) in different types of models, ensuring that the math we use to describe the universe stays grounded in physical reality.

In short: The authors built a better mathematical "sieve," used a supercomputer to measure the tiny details of a complex quantum interaction, and provided a new, more accurate "flavor profile" for how particles interact at very high energies.

1. Problem Statement

The paper addresses the challenge of performing high-order perturbative calculations in Quantum Field Theory (QFT), specifically for a four-dimensional scalar field model with a quartic interaction ( $\phi^4$ theory).

Context: While dimensional regularization (DR) is the standard method for calculating multi-loop coefficients (up to 4-loops and beyond), cutoff regularization (using a momentum cutoff $\Lambda$ ) is crucial for studying the Functional Renormalization Group (FRG) and ensuring the non-negativity of the spectrum of the deformed free Laplace operator.
The Gap: Historically, cutoff regularization calculations were limited to two-loop coefficients. A previous study [22] attempted to calculate three-loop coefficients using a specific cutoff function ( $f=0$ ) but failed to satisfy the physical condition of non-negative spectrum.
Objective: The authors aim to perform a rigorous numerical analysis of auxiliary integrals required for the three-loop approximation of the $\beta$ -function and anomalous dimensions ( $\gamma_\phi, \gamma_\tau$ ) using a physically valid cutoff regularization function ( $f_4$ ) that satisfies the spectral non-negativity condition.

2. Methodology

The authors employ a hybrid analytical-numerical approach to evaluate complex multi-dimensional integrals arising in the three-loop approximation.

Regularization Scheme:
- They utilize a coordinate-space cutoff regularization defined by a continuous function $f(s)$ with support in $[0, 1]$ .
- The specific function chosen is $f_4(s) = 3 - 4s + \frac{2}{\pi}\left(\frac{1-s}{s}\right)^{1/2} + \frac{2}{\pi}\left(\frac{1}{s} - 4\right)\arcsin(\sqrt{s})$ . This function ensures the deformed free Laplace operator has a non-negative spectrum, a condition not met by the previously used $f=0$ case.
Integral Reduction:
- The core of the calculation involves 13 auxiliary functionals ( $\alpha_1(f)$ through $\alpha_{13}(f)$ ), which are high-dimensional integrals (up to 8 dimensions).
- Fourier Transform & Spherical Coordinates: To reduce computational complexity, the authors switch to spherical coordinates and apply the Fourier transform in 4D space. This utilizes the convolution theorem and Bessel functions ( $J_1, J_0$ ) to reduce the dimensionality of the integrals from 8D to 4D (or lower).
- Explicit Transforms: They utilize known analytical Fourier transforms for the specific regularization functions $R^1_0(x)$ and $f^1_2(x)$ under both $f=0$ and $f=f_4$ cases.
Numerical Integration:
- The reduced integrals are computed numerically.
- Algorithm: The infinite integration limits are truncated to finite intervals. The composite trapezoid rule (Newton–Cotes formula) is used for integration.
- Implementation: Calculations were performed using Python, leveraging the NumPy library for array operations and Numba for Just-In-Time (JIT) compilation to accelerate the numerical integration.
- Error Control: The error term is estimated based on the second derivative of the integrand, ensuring high precision (errors reported as $\pm 10^{-8}$ ).

3. Key Contributions

Validation of Cutoff Regularization at Three Loops: The paper provides the first successful numerical determination of three-loop coefficients for the $\phi^4$ model using a cutoff regularization that satisfies physical spectral constraints.
Numerical Values of Auxiliary Integrals: The authors computed precise numerical values for 13 auxiliary integrals ( $\alpha_1$ to $\alpha_{13}$ ) for both the previously studied $f=0$ case and the new physically valid $f=f_4$ case (Table 1).
Derivation of Three-Loop Coefficients: Using the computed integrals, they derived the third-order coefficients for the $\beta$ -function and anomalous dimensions in the Minimal Subtraction (MS) scheme adapted for cutoff regularization.
Comparison and Verification: The results for the $f=0$ case were recalculated to verify the computational algorithm against previous literature, confirming consistency.

4. Key Results

The study yielded specific numerical values for the renormalization group functions.

Three-Loop Coefficients ( $f=f_4$ case):
- $\beta$ -function coefficient ( $\beta_3$ ): $\approx 45.75371 \pm 10^{-5}$ .
- Anomalous dimension of the field ( $\gamma_{\phi,3}$ ): $\approx -0.0638262 \pm 10^{-7}$ .
- Anomalous dimension of the mass ( $\gamma_{\tau,3}$ ): $\approx -3.4982318 \pm 10^{-7}$ .
Comparative Hierarchy:
The authors established the following inequality for the $\beta_3$ coefficients across different schemes:
$\beta_3^{(f=0)} > \beta_3^{(f=f_4)} > \beta_3^{(DR)} > 0$
Where $\beta_3^{(DR)} \approx 32.55$ (Dimensional Regularization) and $\beta_3^{(f=0)} \approx 45.91$ . This confirms that while the specific values depend on the regularization scheme, the physical consistency is maintained.
Power Singularities:
The paper also analyzed the part of the renormalization constant proportional to $\Lambda^2$ (power singularities). They determined the first two coefficients for the associated $\gamma_\Lambda$ function:
- $\gamma_{\Lambda,1} = -3$ (scheme independent).
- $\gamma_{\Lambda,2} = 35/6$ (scheme dependent, MS-scheme).

5. Significance

Bridging Theory and Numerics: The work demonstrates that complex multi-loop calculations in QFT, previously restricted to dimensional regularization due to analytical intractability, can be successfully tackled using advanced numerical methods with physical cutoffs.
Functional Renormalization Group (FRG): By providing valid three-loop coefficients for a cutoff scheme, the paper supports more accurate FRG studies, which are essential for understanding critical phenomena and non-perturbative effects.
Future Directions: The authors suggest that this methodology could be adapted for the three-dimensional sextic model (where four-loop coefficients are now known) and potentially for adapting the standard Bogoliubov–Parasyuk procedure in cutoff regularization.
Robustness: The successful verification of the $f=0$ case against known results and the derivation of new results for $f=f_4$ establishes a robust computational framework for future higher-order calculations in regularized QFT models.