Two-loop amplitude for the associated production of a top-anti-top pair and a $W$ boson at hadron colliders

High-energy physics is the study of the universe's most fundamental building blocks and the forces that bind them. At the heart of this inquiry lies the Standard Model, a theoretical framework that describes how particles like quarks and electrons interact. To test the limits of this model, scientists use massive machines called particle colliders, where they smash protons together at nearly the speed of light. These collisions recreate conditions similar to those just after the Big Bang, producing a shower of new particles. Among the most massive and interesting of these particles is the top quark, which is so heavy it decays almost instantly. When a top quark is created, it is often accompanied by its antimatter partner, the anti-top, and sometimes by a W boson, a particle that carries the weak nuclear force. Understanding exactly how often these specific combinations appear is crucial. If the observed rates differ from theoretical predictions, it could signal the presence of new, unknown physics. However, calculating these rates with extreme precision is a monumental task, requiring scientists to account for complex quantum effects that occur when particles interact multiple times within a single collision.

In a recent contribution to the field, Mattia Pozzoli and his collaborators have taken a significant step forward in this precise calculation. They focused on the process where a top-anti-top pair is produced alongside a W boson. While previous studies had provided good estimates using approximations, this team performed the first exact calculation of the two-loop quantum amplitude for this specific process. In the language of particle physics, a "loop" represents a virtual particle that briefly pops into existence and then disappears during an interaction. A "two-loop" calculation accounts for these complex, double-layered quantum fluctuations, which are essential for reaching the highest possible level of accuracy, known as next-to-next-to-leading order. The researchers worked within the "leading colour" approximation, a mathematical simplification that captures the dominant behavior of the strong nuclear force while ignoring smaller, less significant effects. This allowed them to tackle the immense algebraic complexity of the problem without losing the core physical picture.

The path to this result was fraught with mathematical obstacles. The process involves seven different variables describing the motion and energy of the particles, leading to incredibly intricate equations. Furthermore, the presence of heavy particles inside the calculation meant that standard mathematical tools were insufficient; the equations were tied to shapes known as elliptic curves, which are far more difficult to solve than the simpler shapes used in previous calculations. To overcome this, the team developed a new strategy to express the complex amplitude using a specific set of special functions. They treated these functions like a library of building blocks, calculating their values using a method that tracks how they change as the energy of the collision shifts. For the numerical coefficients that weight these functions, they employed a technique involving finite fields, which allowed them to reconstruct exact values from modular arithmetic rather than struggling with unwieldy algebraic expressions.

Because calculating the full amplitude for every possible collision scenario would take far too long, the team created a five-dimensional grid. They computed the finite remainder of the calculation at nearly 225,000 specific points across this grid and then used interpolation to estimate the values for the entire range of possibilities. This approach allowed them to determine the inclusive cross-section, which is essentially the probability of this specific event happening, for proton collisions at an energy of 13 tera-electronvolts. Their results showed that the exact calculation yields a probability of approximately 242 femtobarns for the negative W boson case and 489 femtobarns for the positive W boson case. These numbers are remarkably close to previous predictions that relied on a combination of approximations, differing by only a few percent. This agreement validates the earlier approximations while confirming that the new, exact method works as intended.

The study also shed light on the impact of the effects they chose to ignore. By comparing their leading-colour results with known one-loop calculations, they estimated that the missing sub-leading colour effects contribute about 22 percent to the virtual correction part of the calculation. Since the two-loop correction itself is about 10 percent of the total cross-section, the uncertainty introduced by ignoring these smaller effects is roughly 2.2 percent. This level of precision is vital for future experiments at the Large Hadron Collider. The researchers' work provides a robust theoretical benchmark, ensuring that when experimentalists measure these rare events, they have a reliable map to distinguish between standard physics and potential discoveries of new particles. The team's success in navigating the complex mathematics of elliptic curves and massive propagators demonstrates that even the most difficult quantum calculations can be tamed with the right combination of advanced algorithms and computational power.