The Machines Are Learning to See the Intractable

Richa Sharma, February 16, 2026

Tags:physicsautomated-reasoningllms
The Machines Are Learning to See the Intractable

In December 1959, Richard Feynman stood before the American Physical Society and pointed toward a territory almost no one was exploring. “There's plenty of room at the bottom,” he said. The laws of physics did not forbid engineering at the atomic scale—we simply hadn’t built the tools to get there. He sketched a world where the entire Encyclopedia Britannica could be written on the head of a pin and machines could be built to build smaller machines. Biology had already proven that information storage and mechanical action were possible at scales far smaller humans ever attempted. The talk was visionary not because it violated physics, but because it took physics literally: if the laws permit something, then the remaining hill to climb is merely an engineering problem. Sixty-six years later, we stand at a different kind of threshold. We are no longer only pushing against the limits of matter. We are pushing inward, against the limits of human symbolic reasoning.

The laws of physics are encoded in scattering amplitudes. When particles collide—at CERN, in the cores of stars, in the first microseconds after the Big Bang—the probabilities of what emerges are governed by laws that obey nature's deepest symmetries. To compute them from first principles, physicists rely on the Feynman diagram expansion, perturbatively summing over every conceivable quantum process. But this approach carries a devastating computational cost: for an event involving nn gluons, the number of required diagrams grows faster than exponential, involving order n!n! terms. It is a combinatorial explosion that rapidly exceeds the capacity of human working memory.

Yet, a profound mystery lies hidden within this chaos. When theorists actually grind through the thousands of diagrams, the final answers often miraculously collapse into simple, elegant formulas—a glaring signal that our present understanding of quantum physics is incomplete and obscures a deeper, more efficient architecture. Recently, a collaboration of physicists discovered a new pillar of this architecture: they proved that "single-minus" gluon interactions, long presumed by textbooks to vanish entirely, actually occur under restricted "half-collinear" kinematics. However, calculating these non-vanishing amplitudes up to just six particles (n=6n=6) resulted in a chaotic algebraic knot of 32 terms. The human researchers knew a beautiful simplification had to exist, but the pattern remained entirely intractable.

To penetrate this mathematical noise, the researchers brought in GPT-5.2 Pro. The language model traversed the human-calculated data, identified hidden patterns, and conjectured a brilliantly simple, piecewise-constant closed-form expression for an arbitrary number of particles. Shortly after, a new internal OpenAI scaffolded-reasoning model independently derived and rigorously proved the formula, satisfying strict physical consistency conditions like Weinberg's soft theorem. The key intuition was entirely human. Physicists recognized that the usual textbook-level no-go argument fails on the half-collinear support. But the episode still marks something important: the point at which frontier models begin to act like collaborators in the discovery pipeline. As these systems improve, what does the future hold for theoretical physics? Should we expect a better mechanical tool than Mathematica, or a new kind of collaborator that guides us through the symbolic maze that obscures our cognition?

1. Background: The Combinatorial Explosion of Scattering Amplitudes

To understand the breakthrough, we must establish the baseline mechanics of how particle physicists calculate the probability of quantum events.

1.1 Gluons and Feynman Diagrams

GluonsInside the proton, gluons bind quarks together through the exchange of color charge. Because gluons themselves carry color, they can also interact with each other. are the exchange particles that mediate the strong nuclear force. Because gluons themselves carry "color charge," they interact with one another, creating incredibly complex webs of quantum events. To calculate the probability of nn gluons scattering (interacting) in a particle collider, physicists compute the scattering amplitude (An\mathcal{A}_n).

Historically, these amplitudes are derived using the Feynman diagramsFeynman Diagram Gluon RadiationA tree-level Feynman diagram illustrating gluon radiation. The complexity of these diagrams grows factorially with the number of particles, necessitating new mathematical frameworks., expansion, which perturbatively sums over all possible quantum processes. But as the number of interacting gluons increases, the number of required diagrams grows faster than exponential, reaching order n!n! terms. To survive this combinatorial explosion, physicists needed a better mathematical language.

1.2 The Spinor-Helicity Formalism

Because gluons are massless, their momentum vectors travel at the speed of light, meaning their invariant mass (the "length" of their 4-vector) is exactly zero. The spinor-helicity formalism takes advantage of this zero-mass property to radically simplify the math. First, the standard 4-momentum vector of a gluon is translated into a 2×22 \times 2 matrix. Because the particle is massless, the determinant of this matrix is zero. In linear algebra, a matrix with a zero determinant can be perfectly factored into two simpler pieces. We decompose this momentum matrix into two 2-component Weyl spinors: pαα˙=λαλ~α˙p_{\alpha\dot{\alpha}} = \lambda_{\alpha}\tilde{\lambda}_{\dot{\alpha}}

Instead of calculating collisions using clunky 4-vectors, physicists can now compute interactions using the Lorentz-invariant inner products of these fundamental spinors. We denote them using a simple bracket notation:

  • Angle bracket: ij=ϵαβλiαλjβ\langle ij \rangle = \epsilon_{\alpha\beta}\lambda_i^{\alpha}\lambda_j^{\beta}
  • Square bracket: [ij]=ϵα˙β˙λ~iα˙λ~jβ˙[ij] = \epsilon_{\dot{\alpha}\dot{\beta}}\tilde{\lambda}_i^{\dot{\alpha}}\tilde{\lambda}_j^{\dot{\beta}}

If we were calculating this in standard, everyday spacetime (Minkowski space), the angle spinor λ\lambda and the square spinor λ~\tilde{\lambda} would be mathematically tethered together as complex conjugates. But to expose the hidden structure of these amplitudes, the researchers performed their calculations in a theoretical geometry known as (2, 2) Klein signature spacetime. In Klein space, the λ\lambda and λ~\tilde{\lambda} spinors are real and entirely independent of one another. As we will see, this geometric freedom—the ability to manipulate the angle bracket without affecting the square bracket—is exactly what allowed physicists to break the "Zero" assumption.

1.3 The MHV Baseline: The n2n-2 Limit

In 1986, Parke and Taylor collapsed the n!n! complexity of Feynman diagrams into a single elegant formula for Maximally Helicity Violating (MHV) tree amplitudes. Gluons carry a helicity (spin state) of either +1+1 (plus) or 1-1 (minus). By definition, an nn-gluon MHV amplitude contains exactly 2 minus-helicity gluons and n2n-2 plus-helicity gluons.

Crucially, for generic tree-level collisions, physicists accepted that n2n-2 was the absolute maximum number of plus-gluons allowed in a non-zero event. While the MHV formula revealed deep symmetries in the strong force, it also erected a strict boundary. It raised the obvious question: if n2n-2 is the maximum, what happens if you attempt to calculate an amplitude with n1n-1 plus-gluons—a "single-minus" interaction?

2. The "Zero" Assumption and The Loophole

For decades, the standard textbook position was that the nn-gluon "single-minus" tree amplitude (one 1-1 helicity, n1n-1 +1+1 helicity) was mathematically equal to exactly zero.

2.1 The Power-Counting Argument (Why n2n-2 was the limit)

To calculate a scattering amplitude, you have to mathematically contract (multiply) the polarization vectors of all the interacting gluons together. For generic collisions, physicists use a mathematical trick: they choose a specific reference spinor that makes all of the incoming and outgoing gluon polarization vectors completely orthogonal to one another. Because they are orthogonal, multiplying them together yields exactly zero. The only way the amplitude can survive and yield a real probability is if these polarization vectors contract with the momentum variables generated by the internal interaction vertices of the Feynman diagramsFeynman diagrams with 4 and 6 particlesFeynman diagrams with 4 particles have 2 vertices; 6-particle tree-level diagrams have 4 vertices.. And here lies the fatal bottleneck: A standard tree-level diagram has, at most, n2n-2 internal vertices. Therefore, there are only enough momentum powers to successfully contract with n2n-2 polarization vectors.

By the pigeonhole principle, the math fails. If you have n1n-1 plus-gluons, you have too many polarization vectors and not enough vertices to "soak them up." For generic kinematics, the single-minus amplitudes strictly vanish.

2.2 The Half-Collinear Loophole (How they proved n1n-1 is possible)

This "Zero" assumption stood unquestioned until the researchers of this paper found a crack in the foundation. The entire power-counting argument hinges on the assumption that you can safely choose a generic reference spinor to make all the polarizations orthogonal. The physicists realized this assumption breaks down under highly restricted kinematics. Specifically, they looked at the half-collinear regime. In this physical setup, the momenta of the particles align so perfectly that all of their spinor angle brackets vanish (ij=0\langle ij \rangle = 0 for all i,ji, j). If you try to use the standard reference spinor trick in this aligned regime, the mathematics shatters—the polarization vectors become singular because you end up dividing by zero.

Because the textbook argument fails in this exact configuration, you cannot conclude the amplitude is zero. The researchers demonstrated that in this narrow, half-collinear locus, n1n-1 (single-minus) amplitudes are not only allowed, but they are supported as non-vanishing entities.The physics permitted the intractable. The only remaining problem was the math—calculating these newly discovered non-zero amplitudes for a large number of particles resulted in a chaotic algebraic knot. They had proven the territory existed, but they needed an AI to map it.

3. The Human Roadblock: n=6n=6

Having proved the single-minus amplitude is non-zero in the half-collinear regime, the team needed to calculate it. They stripped the core kinematics (the stripped amplitude, A1nA_{1\cdots n}) from the delta functions enforcing the half-collinear state, and applied the Berends-Giele recursion relation.

The Berends-Giele recursion computes amplitudes by iteratively building "off-shell" currents (F\mathcal{F}) and merging them at vertices. While powerful, the iterative nature of the recursion creates an explosion of terms based on sign functions (sg(x)sg(x)).

The researchers calculated the stripped amplitude by hand up to n=6n=6. The results were brutal. For 6 gluons, the equation exploded into a 32-term polynomial of nested sign functions:

A123456=18[sg1,23sg12,3sg123,4sg56+ (31 more terms)]A_{123456} = \frac{1}{8} \Big[ -sg_{1,23}sg_{12,3}sg_{123,4}sg_{56} + \dots \text{ (31 more terms)} \Big]

Human intuition had hit a wall. Much like the MHV amplitudes before Parke-Taylor, it was obvious a simpler generalization existed, but the algebraic noise was too dense to parse.

4. Enter AI: Conjecturing and Proving the All-nn Formula

To penetrate the mathematical noise, the researchers utilized GPT-5.2 and an internal scaffolded reasoning model. The AI contribution occurred in two distinct phases: Pattern Recognition and Rigorous Proof Generation.

4.1 Simplifying the Search Space (Region R1\mathcal{R}_1)

Through iterative prompting with GPT-5.2 Pro, the team isolated a restricted kinematic sub-region called R1\mathcal{R}_1, where particle 1 has a negative frequency (ω1<0\omega_1 < 0) and all others are positive (ωa>0\omega_a > 0). In this region, the AI identified dramatic factorizations of the human-derived equations:

A1234R1=14(sg12+sg23)(sg34+sg41)A_{1234}|_{\mathcal{R}_1} = \frac{1}{4}(sg_{12} + sg_{23})(sg_{34} + sg_{41})

A12345R1=18(sg12+sg23)(sg34+sg1,23)(sg45+sg51)A_{12345}|_{\mathcal{R}_1} = \frac{1}{8}(sg_{12} + sg_{23})(sg_{34} + sg_{1,23})(sg_{45} + sg_{51})

4.2 The GPT-5.2 Conjecture

By analyzing the factorized outputs for n{3,4,5,6}n \in \{3,4,5,6\}, GPT-5.2 successfully extrapolated the underlying mathematical pattern, conjecturing a closed-form solution for any arbitrary number of nn gluons:

A1nR1=12n2m=2n1(sgm,m+1+sg1,2m)A_{1\cdots n}|_{\mathcal{R}_1} = \frac{1}{2^{n-2}}\prod_{m=2}^{n-1}(sg_{m,m+1} + sg_{1,2\cdots m})

This is a beautiful result. Because each paired sign function evaluates to (±1±1)( \pm 1 \pm 1 ), every factor is simply 0,2,0, 2, or 2-2. The entire, superexponentially complex Feynman sum collapses into a piecewise-constant integer.

4.3 Autonomous Proof Generation

A conjecture is not a proof. To validate Eq. 39, OpenAI deployed a scaffolded internal model tasked with formally proving the conjecture from first principles.

Over the course of roughly 12 hours of automated reasoning, the AI proved the formula by showing that within Region R1\mathcal{R}_1, a massive interaction vertex term (VV) identically vanishes:

Vλ~2λ~nR1=0V_{\tilde{\lambda}_2\cdots\tilde{\lambda}_n}|_{\mathcal{R}_1} = 0

Because V=0V=0, the entire Berends-Giele recursion relation collapses. The AI rigorously documented this collapse, outputting a formal proof that matched the conjectured formula perfectly, while independently verifying that the formula satisfied cyclic symmetry and Weinberg's soft theorem.

import sympy
from sympy import prod, sign

def verify_amplitude(n, frequencies):
    """
    Programmatically verifies the vanishing of the interaction 
    vertex for the single-minus gluon amplitude.
    """
    # R1 Region: omega_1 < 0, all other omega > 0
    if not (frequencies[0] < 0 and all(f > 0 for f in frequencies[1:])):
        return "Invalid Kinematic Region"
        
    # The AI discovered that the vertex V collapses to zero 
    # identically when these conditions are met.
    vertex_value = sum(frequencies) # Simplification for demonstration
    return vertex_value == 0

5. Implications: The Verifier Era of Physics

When physicists find simple formulas (like Parke-Taylor or Eq. 39) hiding inside complex QFT calculations, it usually indicates that our foundational understanding of the physics is incomplete—that there is a more efficient, unified geometry underlying the universe that we haven't mapped yet.

Extensions to Gravity: The mathematical framework used for gluon scattering (Yang-Mills theory) maps closely to graviton scattering (the theoretical particle of gravity) via "double-copy" relations. The methods pioneered here provide a direct template for AI to explore the unmapped regions of quantum gravity.

AI as a Peer Collaborator: The AI did not just retrieve information; it navigated an exponential search space to find a generalized formula, and then engaged in multi-step logical deduction to prove it. As physicist Nathaniel Craig noted, this establishes a template where humans frame the physical constraints, and AI handles the combinatorial mathematics.

Feynman noted that the barrier to "the bottom" was merely an issue of lacking the right tools—better electron microscopes, tinier lathes. Today, the barrier to "the frontier" of theoretical physics is human cognitive bandwidth for algebraic complexity. By effectively solving the single-minus gluon amplitude, GPT-5.2 has proven that LLMs are the new electron microscopes for the mathematics of the universe.