Operator Renormalization using Emergent Supersymmetries

Imagine you are trying to solve a massive, incredibly complex jigsaw puzzle. This puzzle represents the laws of physics that govern how particles interact (specifically, a theory called Quantum Chromodynamics, or QCD, which describes the strong force holding atomic nuclei together).

For decades, physicists have been trying to finish this puzzle. But the pieces are so tiny and the connections so complicated that calculating even a small section can take a supercomputer months to run. It's like trying to count every grain of sand on a beach by hand.

This paper introduces a clever new trick: borrowing a shortcut from a parallel universe.

The Problem: The "Hard Mode" Universe

The universe we are studying (the Standard Model) is "Hard Mode." It's messy, non-symmetric, and every calculation requires brute force. There are no shortcuts, no hidden patterns to simplify the math.

The Solution: The "Easy Mode" Universe (Supersymmetry)

There is a theoretical "parallel universe" called Supersymmetry (SUSY). In this universe, nature is perfectly balanced. For every particle, there is a "super-partner." Because of this perfect balance, the math is much cleaner. There are "Ward identities"—think of these as magic rules or shortcuts that tell you: "If you know the answer for Particle A, you automatically know the answer for Particle B without doing any extra work."

The problem? We don't think our real universe is Supersymmetric. So, we can't just use these magic rules.

The Breakthrough: "Emergent" Supersymmetry

The authors, Mrigankamauli Chakraborty and Sven-Olaf Moch, came up with a brilliant idea. They built a "Universal Simulator" (called a Generalized Lagrangian).

Think of this simulator as a video game engine that can run two different games at once:

Game A: Our messy, real-world physics.
Game B: The perfect, balanced Supersymmetric physics.

They set up the engine so that if you tweak the settings (specifically, the number of "flavors" of particles), the game suddenly switches from "Hard Mode" to "Easy Mode."

Here is the magic trick:

They run the simulation in "Hard Mode" (our real world).
They tweak the settings to hit a specific point where the game accidentally becomes "Easy Mode" (Supersymmetry).
At that exact point, the magic rules (Ward identities) kick in.
They use those rules to simplify the math back in the Hard Mode version.

It's like trying to solve a difficult math problem. You realize that if you pretend the numbers were slightly different, the problem becomes a simple equation you can solve in your head. You solve the simple version, then use that answer to instantly figure out the answer to the difficult version, saving you hours of work.

The Real-World Result: A 25% Speed Boost

The authors tested this on a specific model called the Gross–Neveu–Yukawa (GNY) model.

Without the trick: Calculating certain particle behaviors took about 14 days of computer time.
With the trick: By using the "borrowed" magic rules from the Supersymmetric version, they reduced the work by 25%.

This might sound small, but in the world of high-energy physics, a 25% reduction is huge.

If a calculation takes 14 days, they now save 3.5 days.
If a calculation for our real universe (QCD) takes 6 months, this trick could save 1.5 months.

Why Does This Matter?

These calculations are crucial for understanding Parton Distribution Functions. In plain English, these are the "maps" that tell us how the tiny particles inside protons and neutrons are arranged. These maps are essential for:

Predicting what happens in the Large Hadron Collider (LHC).
Understanding the early universe.
Designing future particle accelerators.

The Bottom Line

The authors aren't trying to prove that Supersymmetry exists in nature. Instead, they are treating it as a computational tool. They are saying, "Even if the universe isn't supersymmetric, we can pretend it is for a split second to unlock a cheat code, use that cheat code to do the math faster, and then apply the result to our real, non-supersymmetric universe."

It's a bit like using a perfectly symmetrical snowflake to understand how to build a messy, asymmetrical snowman. The snowflake gives you the blueprint to build the snowman much faster than you could on your own.

Here is a detailed technical summary of the paper "Operator Renormalization using Emergent Supersymmetries" by Mrigankamauli Chakraborty and Sven-Olaf Moch.

1. Problem Statement

The primary challenge addressed in this work is the immense computational cost associated with high-order perturbative calculations in non-supersymmetric (non-SUSY) phenomenological theories, such as Quantum Chromodynamics (QCD) and the Standard Model.

Computational Bottleneck: Calculating higher-loop corrections (e.g., three-loop or higher) for operator renormalization involves complex numerator structures. This leads to lengthy Integration-by-Parts (IBP) reduction times, often spanning days or weeks for specific models and months for QCD.
Lack of Symmetry Constraints: Unlike Supersymmetric (SUSY) theories, non-SUSY theories lack powerful Ward identities (relations between Green's functions generated by super-Poincaré symmetry) that could reduce the number of independent quantities to be computed.
Goal: The authors aim to utilize SUSY not as a phenomenological model (since experimental evidence is lacking), but as a computational optimization tool to streamline calculations in experimentally relevant non-SUSY theories.

2. Methodology

The authors propose a formalism based on a Generalized Lagrangian ( $\mathcal{L}_{gen}$ ) that unifies various scalar-fermion theories, allowing SUSY Ward identities to "emerge" at specific parameter points and be applied to non-SUSY limits.

A. The Generalized Lagrangian

They define a theory containing real scalars ( $\phi_I$ ) and Weyl fermions ( $\psi_A$ ) with a global $SO(n_s)$ flavor symmetry for scalars and $SU(n_f)$ for fermions. The Lagrangian includes:

Kinetic terms for scalars and fermions.
A ternary Yukawa interaction with flavor tensors $Z_I$ and $\bar{Z}_I$ .
A quartic scalar coupling with flavor tensor $M_{IJKL}$ .

B. Theorem on Non-Evanescent Flavour Terms

To ensure the generalized Lagrangian yields valid results for standard models (like Gross-Neveu-Yukawa, NJLY, and Wess-Zumino) without introducing "evanescent" terms (terms that vanish in $d=4$ but appear in dimensional regularization), the authors prove Theorem 1:

The Lagrangian contains no evanescent flavor terms if and only if the flavor tensors satisfy the relation:
$Z_I \bar{Z}_J + Z_J \bar{Z}_I = 2\delta_{IJ}$
This is the same algebraic relation satisfied by Pauli matrices.
This condition, combined with proportionality constraints on primitive Feynman graphs, completely determines the flavor factors ( $Z, \bar{Z}, M$ ) as functions of the number of scalars ( $n_s$ ) and fermions ( $n_f$ ).

C. Emergent Supersymmetry

The authors calculate the anomalous dimensions ( $\gamma$ ) and $\beta$ -functions for $\mathcal{L}_{gen}$ up to four loops using the FORCER package (via Qgraf and FORM).

They identify "emergent SUSY" points in the $(n_f, n_s)$ plane where the SUSY Ward identity $\gamma_f = \gamma_s$ (equality of fermion and scalar anomalous dimensions) holds.
Key Finding: Curves for $\gamma_f = \gamma_s$ $γ_{f} = γ_{s}$ at every loop order intersect at two specific points:
1. $(n_s, n_f) = (2, 1)$ : Corresponds to the 4D Wess-Zumino model (2 scalar, 2 fermionic degrees of freedom).
2. $(n_s, n_f) = (1, 1/2)$ : Corresponds to the 3D Wess-Zumino model (1 scalar, 1 fermionic degree of freedom).
Additionally, the curve for the highest transcendental number at any loop order coincides with the one-loop condition $n_s - 2n_f = 0$ (equality of bosonic and fermionic degrees of freedom).

3. Key Contributions

Unified Formalism: Development of a single generalized Lagrangian that encompasses both SUSY and non-SUSY theories (GNY, NJLY, Wess-Zumino) by simply substituting values for $n_s$ and $n_f$ .
Ward Identity Transfer: Demonstration that Ward identities valid at the emergent SUSY points can be used to derive relations among integrals in the generalized theory, which subsequently apply to the non-SUSY limit (e.g., GNY model where $n_s=1$ ).
Optimization Mechanism: A method to eliminate specific Feynman integrals using the derived Ward identities before performing the actual loop integration.
- For example, at two loops, the Ward identities allow the elimination of two integrals.
- At three loops, independent relations for non-planar topologies allow further reductions.

4. Results

Validation: The calculated quantities at the emergent SUSY points $(2,1)$ and $(1, 1/2)$ perfectly match known literature values for the 4D and 3D Wess-Zumino models, validating the formalism.
Operator Renormalization: The method was applied to the renormalization of twist-two flavor-singlet operators in the Gross-Neveu-Yukawa (GNY) model.
- Ward Identity Verification: The identity $2\gamma_{\psi\psi} + 4\gamma_{\phi\psi} = 2\gamma_{\phi\phi} + \gamma_{\psi\phi}$ was verified to hold at the emergent SUSY points up to three loops.
Computational Efficiency:
- In the GNY model, calculating all even moments ( $N$ ) up to $N=30$ at three loops normally takes ~14 days.
- By applying the emergent SUSY optimization (eliminating non-planar integrals), the computational time was reduced by 25% across all values of $N$ .
- This reduction is achieved with negligible extra overhead, as the constraints are derived algebraically from the generalized Lagrangian.

5. Significance and Future Outlook

Phenomenological Impact: The anomalous dimensions of twist-two operators are directly related to Parton Distribution Functions (PDFs) via the Wilson Operator Product Expansion (OPE). Reducing the computational time for these calculations is crucial for precision QCD phenomenology.
Scalability: While the 25% speedup in the GNY model saves days, the authors project that applying this to QCD (where calculations currently take months) could save weeks to months of computation time.
Future Work: The authors intend to extend this formalism to QCD. They have already constructed a generalized Lagrangian unifying QCD with $N=1, 2, 4$ Super Yang-Mills (SYM) theories. The next step involves handling the subtleties introduced by gauge symmetry to fully implement the emergent SUSY optimization in the QCD context.

In summary, the paper presents a novel "computational bridge" where supersymmetry is treated as a mathematical tool to generate constraints in non-supersymmetric theories, significantly accelerating high-order perturbative calculations essential for modern particle physics.