$N$-Jettiness Soft Functions Made Simple

This paper introduces a simplified method for computing $N$-Jettiness soft functions at arbitrary $N$ and high perturbative orders by decomposing the dipole contribution into an analytically calculable inclusive part and a finite remainder, enabling efficient NNLO calculations for up to five jets and outlining prospects for N$^3$LO applications.

The Gist

Imagine you are trying to predict exactly how a massive, chaotic crowd of people will move through a stadium after a game ends. You want to know not just the total number of people leaving, but exactly how they are grouped, where the bottlenecks are, and how they interact with each other.

In the world of particle physics, the "stadium" is the Large Hadron Collider (LHC), the "people" are subatomic particles called quarks and gluons, and the "movement" is the creation of jets of particles after a collision. Physicists call this QCD (Quantum Chromodynamics).

The paper you provided, "N-Jettiness Soft Functions Made Simple," is about a new, much smarter way to calculate the behavior of these particles, specifically when they are moving very slowly (the "soft" part) and interacting in complex groups.

Here is the breakdown using everyday analogies:

1. The Problem: The "Traffic Jam" of Math

For decades, physicists have been trying to predict the results of particle collisions with extreme precision. To do this, they use a mathematical tool called N-Jettiness. Think of N-Jettiness as a way to measure how "messy" the traffic is.

• N stands for the number of main "lanes" or jets of particles flying out.
• Jettiness measures how much the particles are "jittering" or spreading out.

The problem is that calculating the interactions of these particles gets incredibly hard as you add more particles (more jets) or try to calculate them with higher precision (like going from a rough sketch to a 4K movie). The math involves "infrared singularities," which is a fancy way of saying the equations blow up and become infinite when particles get too close to each other or move too slowly.

Previously, calculating this for complex scenarios (like 5 jets flying out) was like trying to solve a Sudoku puzzle while someone is screaming in your ear and shaking the table. It took years of supercomputer time and brilliant mathematicians to get even one step forward.

2. The Solution: The "Smart Split"

The authors of this paper realized they were trying to solve the whole messy puzzle at once. Their new method is like realizing you can split the problem into two much easier parts:

Part A: The "Inclusive" Crowd (The Easy Part)
Imagine you only care about the total number of people leaving the stadium, regardless of which specific exit they use. This is a simpler calculation that physicists have already solved perfectly. The authors call this the Inclusive Soft Function. It's like knowing the total traffic volume on the highway.

Part B: The "Specific" Crowd (The Hard Part)
Now, imagine you need to know exactly how many people are stuck in a specific bottleneck between two specific exits. This is the "remainder."

• The Magic Trick: The authors realized that the "hard" part (the specific bottleneck) is actually finite. It doesn't blow up to infinity.
• The Analogy: If the "Inclusive" part is the ocean (huge, deep, known), the "Remainder" is just a small puddle on the sidewalk. You don't need a boat to measure the puddle; you just need a cup.

3. How They Did It: The "Subtraction" Technique

The core of their method is a clever subtraction:

1. Calculate the Ocean: They use the known, easy math for the total traffic (Inclusive).
2. Calculate the Puddle: They calculate the difference between the real, messy traffic and the smooth "ocean" traffic.
3. The Result: Because the "ocean" part handles all the dangerous infinities, the "puddle" (the remainder) is safe, finite, and easy to measure with a computer.

They also discovered that for the next level of complexity (called N3LO, or "Next-to-Next-to-Next-to-Leading Order"), this "puddle" will be even easier to measure—roughly as easy as a standard traffic count, rather than a complex simulation.

4. The "Tripole" Surprise

In the paper, they also mention a "Tripole" contribution.

• Dipole: Two people interacting (like a conversation between two friends).
• Tripole: Three people interacting (like a heated argument between three friends).
• The Breakthrough: They derived a very simple formula for this three-way interaction. Before, calculating a three-way argument in this context was a nightmare. Now, they have a "shortcut" formula that allows computers to crunch the numbers in seconds instead of days.

5. Why This Matters

The Large Hadron Collider is about to get even more powerful (the High-Luminosity LHC). It will produce millions of collisions. To find new physics (like dark matter or new particles), scientists need to know the "Standard Model" predictions with percent-level precision.

If you are looking for a needle in a haystack, you need to know exactly what the hay looks like. If your calculation of the hay is off by even 1%, you might mistake a piece of hay for a needle.

This paper provides a scalable, fast, and simple toolkit to calculate the "hay" (the background noise) for complex jet events.

• Before: Calculating 5 jets took a massive effort and was prone to errors.
• Now: They can calculate 5 jets (and potentially more) quickly and accurately.

Summary

Think of this paper as the invention of a new GPS algorithm for particle physics.

• Old GPS: Tries to calculate every single car's path, every pothole, and every traffic light simultaneously. It crashes when the map gets too big.
• New GPS: Calculates the main highway flow (which is easy) and then only calculates the specific detours where cars get stuck. It's faster, simpler, and works for any size of city (any number of jets).

This allows physicists to finally tackle the most complex collisions at the LHC with the precision needed to discover the secrets of the universe.

Technical Summary

1. Problem Statement

The calculation of high-precision QCD predictions for processes involving jets at the Large Hadron Collider (LHC) is a critical challenge for the High-Luminosity LHC (HL-LHC) era. To achieve percent-level accuracy, theoretical predictions must reach Next-to-Next-to-Leading Order (NNLO) and eventually Next-to-Next-to-Next-to-Leading Order (N$^3$LO).

A major bottleneck in extending N-Jettiness slicing (a method used to isolate infrared-sensitive regions of phase space) to N$^3$LO is the computation of soft functions involving an arbitrary number of light-like directions ($N$).

• Complexity: Traditional methods for computing soft functions at three loops (N$^3$LO) or even two loops (NNLO) for arbitrary $N$ involve extremely complex multi-loop integrals and intricate infrared (IR) subtraction schemes.
• Current Limitations: While NNLO results for arbitrary $N$ have been achieved recently, extending these to N$^3$LO has proven prohibitively difficult due to the growth in complexity and the need for sophisticated subtraction techniques.

2. Methodology

The authors propose a novel computational strategy that drastically simplifies the calculation of N-Jettiness soft functions by decomposing the problem into an analytically solvable part and a numerically tractable finite remainder.

Core Concept: Decomposition of the Dipole Contribution

The soft function is expanded in terms of color structures. The dominant contribution comes from dipole correlators (interactions between two hard emitters). The authors decompose the $n$-th order dipole contribution $\tilde{s}^{(n)}_{T, \{i,j\}}$ as:
$$\tilde{s}^{(n)}_{T, \{i,j\}} = \tilde{s}^{(n)}_{T \text{ inc.}, \{i,j\}} + \Delta\tilde{s}^{(n)}_{T, \{i,j\}}$$

1. Inclusive Soft Function ($\tilde{s}^{(n)}_{T \text{ inc.}}$):

   • This term is calculated using an inclusive approximation of the N-Jettiness observable, where the observable depends only on the total momentum of the real radiation ($k_{tot}$) rather than individual momenta.
   • Key Insight: The UV divergences of the full N-Jettiness soft function and this inclusive approximation are identical.
   • Benefit: The inclusive function is analytically known up to three loops (derived from fully differential soft functions and threshold expansions). This allows the authors to subtract the known singularities analytically.

2. Finite Remainder ($\Delta\tilde{s}^{(n)}_{T, \{i,j\}}$):

   • This term represents the difference between the full observable and the inclusive approximation.
   • Finiteness: Because the UV poles cancel exactly, $\Delta\tilde{s}^{(n)}$ is finite as $\epsilon \to 0$ (in dimensional regularization).
   • Reduced Complexity: The IR structure of this remainder is significantly simpler than the full function.
     • At NNLO ($n=2$), the remainder requires no IR subtraction; it can be evaluated directly in 4 dimensions.
     • At N$^3$LO ($n=3$), the structure resembles an NLO cross-section calculation, requiring only standard NLO-type subtractions.

Tripole Correlator

For processes with $N \ge 2$ jets, tripole correlators (involving three hard emitters) appear starting at NNLO. The authors derive a new, compact analytical formula for the tripole contribution. By parametrizing the phase space using Sudakov variables and specific sector decompositions, they isolate the singularities and provide a representation that is ready for fast numerical integration.

3. Key Contributions

• Generalized Method for Arbitrary $N$: The method is applicable to any number of hard emitters, overcoming the limitation of previous approaches that were often restricted to specific low-$N$ cases.
• Analytic-Numeric Hybrid: The strategy separates the calculation into a fully analytic part (handling all divergences) and a finite numerical part. This avoids the need for complex, order-by-order IR subtraction codes for the full soft function.
• Simplified Tripole Formula: The derivation of a compact, numerically efficient formula for the tripole contribution at NNLO, which was previously a significant computational hurdle.
• Scalability to N$^3$LO: The framework is explicitly designed to be scalable. The authors argue that the complexity of the remainder term at N$^3$LO is comparable to standard NLO calculations, making the N$^3$LO soft function computationally feasible.

4. Results

The authors applied their method to compute the N-Jettiness soft function at NNLO for general $N$.

• Validation:
  • 0-Jettiness: The results were compared against known analytic expressions for the thrust soft function, showing exquisite agreement.
  • 1-, 2-, and 3-Jettiness: Numerical results were benchmarked against existing literature (Refs [40, 41]), confirming excellent agreement for both gluonic and fermionic channels.
• New Predictions:
  • The paper provides the first numerical results for 4-Jettiness and 5-Jettiness soft functions at NNLO for specific benchmark phase-space points.
• Performance:
  • The numerical evaluation is extremely fast. On a single thread (MacBook Pro M3), achieving 0.1% precision takes:
    • $\sim 0.1$–$1$ second for 1-Jettiness.
    • $\sim 1$–$10$ seconds for 5-Jettiness.
  • The speed is attributed to the fact that the remainder term is finite and requires no subtraction, allowing for direct Monte Carlo integration.

5. Significance

This work represents a major step forward in precision QCD phenomenology for the LHC:

1. Enabling N$^3$LO Jet Physics: By simplifying the soft function calculation, this method removes a primary bottleneck preventing N$^3$LO predictions for processes with jets (e.g., Higgs + jet, Drell-Yan + jet).
2. Factorization Robustness: It provides a systematic way to handle color correlations and factorization breaking effects that become relevant at higher orders.
3. Efficiency: The ability to compute soft functions for up to 5 jets in seconds makes it practical to include these corrections in global fits and experimental analyses.
4. Future Outlook: The authors outline that this approach paves the way for fully differential N$^3$LO predictions for generic jet processes, which are essential for maximizing the physics potential of the HL-LHC.

In summary, the paper transforms the calculation of N-Jettiness soft functions from a prohibitively complex multi-loop problem into a manageable task by leveraging the universality of inclusive soft functions and isolating the process-dependent finite remainders.