⚛️ phenomenology

Integration-by-parts identities in the presence of measurement delta functions

This paper derives a systematic method for generating integration-by-parts identities for dimensionally regulated loop integrals with non-linear measurement constraints by utilizing vectors tangent to the constraint hypersurface, thereby avoiding derivatives of distributions and enabling the calculation of differential distributions like transverse-momentum spectra in processes such as $t\bar tH$ and $t\bar t$ production.

Original authors: Saimeng Zhou

Published 2026-08-25

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Saimeng Zhou

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ 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

In the high-stakes world of particle physics, scientists smash atoms together at nearly the speed of light to uncover the fundamental building blocks of the universe. To make sense of the debris from these collisions, they rely on complex mathematical calculations that predict how often specific particles appear and with what energy. These predictions are essential for comparing theory against the data collected by massive detectors at facilities like the Large Hadron Collider. However, calculating these probabilities is notoriously difficult when physicists want to know not just the total number of particles produced, but how they are distributed across different energies and directions. This requires tracking the momentum of particles as they fly out from the collision point, a task that turns the standard mathematical tools into a tangled web of equations that are often impossible to solve with existing methods.

A researcher has developed a new mathematical strategy to untangle this web, specifically for cases where the measurement of a particle's motion involves a curved, non-linear relationship. In their work, they derived a set of rules that allow these difficult calculations to be simplified without losing the crucial information about the particle's path. By focusing on the geometry of the measurement itself, they found a way to generate the necessary equations directly, avoiding the mathematical dead ends that have previously blocked progress. They tested this new approach on two specific scenarios involving the production of top quarks and Higgs bosons, and the results matched perfectly with independent, established calculations. This success confirms that their method works, opening a clear path for future studies of complex particle distributions that were previously out of reach.

The core challenge in these calculations lies in how physicists handle the "measurement" of a particle's properties. When a particle is created in a collision, its behavior is described by a set of variables, such as its transverse momentum, which is a measure of how hard it kicks sideways relative to the beam. In simple cases, this sideways kick relates to other variables in a straight, linear way, making it easy to plug into standard calculation tools. However, for many important measurements, the relationship is curved and non-linear. Imagine trying to navigate a city where the streets twist and turn in complex patterns; a map designed for a grid of straight streets would fail to guide you. Similarly, the standard mathematical tools used by physicists, known as integration-by-parts identities, are designed for straight-line relationships. When faced with a curved measurement constraint, these tools break down because they cannot easily handle the non-linear geometry without introducing errors or becoming unsolvable.

To solve this, the researcher realized they did not need to force the curved measurement into a straight-line mold. Instead, they looked at the shape of the constraint itself. They discovered that if they chose their mathematical "vectors"—the directions in which they move through the calculation—to run parallel to the surface defined by the measurement, they could bypass the difficulties entirely. By moving along the curve rather than trying to cut across it, they could generate the necessary equations while keeping the measurement constraint intact as a single, unbroken factor. This approach meant they never had to differentiate the measurement function, a step that usually causes the equations to explode in complexity. The result was a clean, systematic set of rules that worked exactly like the standard tools but were adapted to handle the curved geometry of the problem.

The researcher put this new method to the test in two demanding real-world scenarios. First, they calculated the distribution of the Higgs boson's transverse momentum in collisions where top quarks and Higgs bosons are produced together. This is a leading-order calculation, meaning it is a fundamental prediction of the theory. Second, they tackled the distribution of a single top quark's transverse momentum in collisions where top-antitop pairs are produced, a calculation that is significantly more complex because it includes the next level of quantum corrections. In both cases, they implemented their new rules into a computer program designed to handle these massive algebraic reductions. The program successfully reduced thousands of complex terms down to a manageable set of master integrals, which could then be solved to produce the final predictions.

The validation of this work was rigorous. The researcher compared their results against two other independent methods that are already trusted by the physics community. For the Higgs boson calculation, their results agreed with a well-known software package called MCFM to within a fraction of a percent. For the more complex top quark calculation, they compared their findings against a different framework known as Stripper. The agreement was so precise that the differences were consistent with the tiny statistical uncertainties inherent in the computer simulations. This level of agreement confirms that their new mathematical construction is not just a theoretical curiosity but a robust, working tool. It proves that by respecting the geometry of the measurement, one can navigate the complex landscape of particle physics calculations without getting lost.

This breakthrough is significant because it removes a major bottleneck in the field. Previously, calculating distributions for non-linear measurements required either approximations that might miss subtle effects or entirely different, often cumbersome, techniques that did not generalize well. The new method provides a unified, systematic route to these calculations. It allows physicists to keep the measurement constraint explicit throughout the process, ensuring that the final results are exact and reliable. The author demonstrated that this approach works for single particles, but their framework is general enough to apply to more complex systems, such as the collective recoil of multiple particles or differences in rapidity between them. By solving the problem of how to handle these curved constraints, the researcher has equipped the community with a powerful new lens through which to view the intricate details of particle collisions.

The implications of this work extend beyond the specific calculations performed. It offers a new way of thinking about how to handle constraints in high-energy physics. By showing that the measurement delta function can be treated as a constant factor when the right mathematical directions are chosen, the researcher has simplified a process that was previously fraught with technical obstacles. This clarity allows for the exploration of more complex observables that were previously too difficult to compute with high precision. As experiments at the Large Hadron Collider continue to collect data with increasing accuracy, the ability to make equally precise theoretical predictions becomes ever more critical. This new method ensures that the theoretical tools are ready to meet the demands of the next generation of discoveries, allowing scientists to probe the deepest secrets of the universe with greater confidence and clarity.

Technical Summary: Integration-by-parts identities in the presence of measurement delta functions

Problem Statement
In perturbative QCD, computing differential cross sections for infrared-safe observables often requires integrating over phase space with specific measurement constraints. The "reverse unitarity" method facilitates this by rewriting on-shell delta functions as cut propagators, allowing phase-space integrals to be treated as cut loop integrals. This enables the use of Integration-by-Parts (IBP) identities to reduce integrals to a basis of master integrals, which are then solved via differential equations.

However, a significant obstruction arises when the measurement constraint $\delta(\phi(\ell))$ involves a non-linear function $\phi$ of the loop momentum $\ell$ . While linear constraints (e.g., rapidity distributions) can be incorporated as additional cut propagators within a standard Laporta family, non-linear constraints (e.g., transverse momentum $p_T$ ) introduce quadratic relations among scalar products. Promoting such constraints to propagators leads to non-linear maps between propagators and scalar products, often resulting in square-root branches that prevent a closed IBP reduction. Previous approaches to non-linear constraints, such as those for energy-energy correlations, required supplementing the IBP system with ad-hoc relations to lower propagator powers, while other methods (e.g., Higgs-differential approaches) relied on factorizing the phase space, which is not universally applicable.

Methodology
The author proposes a systematic method to derive IBP identities for dimensionally regulated integrals containing a general non-linear measurement constraint $\delta(\phi(\ell))$ without promoting the delta function to a propagator. The core of the methodology relies on the geometric properties of the hypersurface defined by $\phi(\ell) = 0$ .

Tangent Vector Construction:
Standard IBP identities are derived from the vanishing of total derivatives $\int d^d\ell \, \partial_\mu [v^\mu f(\ell) \delta(\phi)]$ . When expanding the derivative, a term involving $\partial_\mu \delta(\phi) = \delta'(\phi) \partial_\mu \phi$ appears. The author demonstrates that if the generating vector $v^\mu$ is chosen to be tangent to the hypersurface $\phi=0$ (i.e., $v^\mu \partial_\mu \phi = 0$ ), the term involving the derivative of the delta function vanishes. Consequently, the measurement delta function $\delta(\phi)$ remains as an overall factor, and the resulting identity takes the form of an ordinary IBP relation:
$0 = \int d^d\ell \, \delta(\phi) \, \partial_\mu (t^\mu f)$
where $t^\mu$ is a tangent vector. This allows the use of conventional Laporta elimination algorithms without differentiating distributions.
Polynomial Formulation (Weak Tangency):
To avoid introducing inverse powers of scalar products (which complicate the reduction), the author introduces a "weak" tangency condition. They construct vectors $\tilde{t}^\mu$ such that $\tilde{t}^\mu \partial_\mu \phi = h \phi$ for some regular function $h$ . Using the distributional identity $x \delta'(x) = -\delta(x)$ , these vectors generate closed IBP relations with polynomial coefficients. Specifically, for the $p_T$ constraint, they derive a "weak inverse normal" that yields a polynomial differential operator, avoiding the need for the exact inverse normal which would introduce denominators like $1/(p_1 \cdot \ell)$ .
Differential Equations:
The same tangent vectors provide the derivative operators required to evolve master integrals with respect to kinematic variables (such as $s$ or $p_T^2$ ). By differentiating the integral with respect to a variable $x$ and using the weak inverse normal to integrate by parts, the author derives differential equations where the right-hand side consists of ordinary integrands multiplied by $\delta(\phi)$ , reducible via the tangent-vector IBP identities.
Boundary Conditions:
To solve the differential equations, boundary conditions are obtained by "unfolding" the measurement constraint. This involves restoring the cut associated with the measured particle to its on-shell phase-space form and integrating over the remaining continuous degrees of freedom (e.g., rapidity). This reduces the problem to lower-loop integrals. The author addresses endpoint singularities in these boundary integrals by subtracting the singular behavior at fixed $\epsilon$ before performing the numerical integration.

Key Contributions and Results

Systematic IBP Generation: The paper provides a general construction for generating IBP identities for non-linear measurement constraints, keeping the delta function as a spectator factor. This avoids the complications of non-linear propagator maps.
Implementation in Kira: The method is implemented as a source-level extension of the IBP reduction code Kira.
Validation: The framework is validated against independent calculations for two specific processes:
1. $pp \to t\bar{t}H$ at Leading Order (LO): The Higgs transverse-momentum distribution is computed and compared with MCFM, showing agreement within 0.02% ( $q\bar{q}$ channel) and 0.06% ($gg$ channel).
2. $pp \to t\bar{t}$ at Next-to-Leading Order (NLO): The single-inclusive top-quark transverse-momentum distribution is computed. The real and virtual contributions reproduce the expected infrared pole structure with high precision. The finite parts agree with FONLL kernels at the level of $10^{-4}$ , and the full hadronic spectrum matches independent calculations using the Stripper framework within Monte Carlo uncertainties.
Generalizability: The appendix demonstrates that the construction extends to other non-linear observables, including the transverse recoil of multiparticle systems, rapidity differences, and energy-energy correlations.

Significance
The paper claims that this construction provides a systematic route to differential distributions within reverse unitarity for non-linear measurement functions, specifically addressing the long-standing difficulty of $p_T$ distributions. By keeping the measurement delta function explicit and utilizing tangent vectors, the method allows for the application of standard multi-loop techniques (IBP reduction and differential equations) to differential cross sections without requiring ad-hoc modifications to the integral families or the promotion of constraints to non-linear propagators. The successful validation in $t\bar{t}H$ and $t\bar{t}$ production demonstrates the practical viability of the approach for high-precision QCD calculations.

More like this