Physics inspired quantum algorithm for QCD splitting functions

This paper introduces a modular quantum circuit primitive that models QCD parton splitting dynamics by encoding helicity entanglement and momentum-sharing fractions, successfully validating the approach against LHC data and demonstrating its feasibility on current superconducting quantum hardware.

The Gist

Imagine a high-energy particle collision at the Large Hadron Collider (LHC) as a chaotic game of "billiards," but instead of solid balls, we are dealing with tiny, invisible particles called gluons. When these gluons crash into each other, they don't just bounce off; they split apart, creating new gluons, which then split again, creating a cascading shower of particles. This process is called a parton shower.

For decades, scientists have simulated these showers using classical computers. They treat every split as a simple, random decision, like flipping a coin. But the authors of this paper argue that this misses a crucial piece of the puzzle: quantum entanglement. In the quantum world, when two particles are created from a split, they remain mysteriously linked, no matter how far apart they go. Classical computers ignore this link, but the universe doesn't.

Here is how the paper tackles this problem, explained through simple analogies:

1. The "Magic Split" (The Quantum Primitive)

The authors built a tiny, modular "building block" for a quantum computer. Think of this block as a magic splitter.

• The Goal: When a parent particle splits into two children, the magic splitter needs to do two things at once:
  1. Decide how much "momentum" (energy/motion) each child gets.
  2. Create the correct amount of "quantum entanglement" (the invisible link) between them, exactly as nature dictates.
• The Innovation: Instead of just guessing the split, they used the laws of physics (Quantum Chromodynamics, or QCD) to calculate exactly how much entanglement should exist. They found a mathematical formula for this "entanglement" based on how the momentum is shared.

2. The "Two-Qubit Circuit" (The Machine)

To mimic this magic splitter, they designed a simple circuit using just two qubits (the quantum equivalent of bits).

• Imagine the two qubits as two spinning coins.
• The authors programmed the circuit so that if you look at the coins, their behavior tells you exactly how the momentum was shared (e.g., 70% to one, 30% to the other).
• Crucially, the way the coins spin is also "entangled." If you measure one, it instantly affects the state of the other, perfectly matching the complex math of the real-world particle split.

3. Learning from the Real World (Calibration)

The team didn't just guess the settings for their quantum circuit. They went to the AspenOpenJets dataset, which contains real data from the LHC.

• They looked at real "jets" (sprays of particles) and measured how momentum was shared in the first split (the "two-prong" structure).
• They then adjusted the knobs (parameters) on their quantum circuit until its output matched the real-world data.
• The Result: The circuit learned to replicate the real-world momentum sharing while keeping the correct quantum entanglement.

4. Building a Tower (From Two to Many)

The real power of this approach is composition.

• Once they had a working "two-prong" splitter, they could stack them.
• Imagine taking the "heavier" child from the first split and feeding it into a second magic splitter. That child splits again, creating two more.
• By chaining these blocks together, they created circuits that could simulate three-prong and four-prong structures (three or four final particles).
• They tested this against real LHC data and found that their quantum-built towers matched the real-world particle sprays almost perfectly.

5. The Real-World Test (Running on Hardware)

Finally, they didn't just simulate this on a supercomputer; they actually ran the three-prong version on a real quantum computer (an IBM machine called ibm_Marrakesh).

• The Challenge: Real quantum computers are noisy and prone to errors.
• The Success: Despite the noise, the results were very close to the simulation and the real data. This worked because their circuit was so simple (only a few qubits and a shallow depth) that the errors didn't ruin the picture.

The Bottom Line

This paper introduces a new way to simulate particle physics. Instead of treating particle splits as simple, random events, they created a quantum-native tool that respects the "spooky" connections (entanglement) nature demands.

They proved that:

1. You can calculate exactly how much entanglement a particle split creates.
2. You can build a simple quantum circuit that mimics this split and the entanglement.
3. You can stack these circuits to simulate complex particle showers.
4. This works on real quantum hardware and matches real experimental data.

This is a foundational step toward a future where quantum computers don't just calculate numbers, but naturally "act out" the quantum dance of the universe's smallest building blocks.

Technical Summary

Technical Summary: Physics-Inspired Quantum Algorithm for QCD Splitting Functions

Problem Statement
Event generation in High-Energy Physics (HEP), particularly for the High-Luminosity LHC (HL-LHC), faces significant computational bottlenecks. Current state-of-the-art generators (e.g., Pythia, Herwig) rely on Markov Chain Monte Carlo algorithms to sample QCD splitting functions. These methods treat branching steps classically, discarding quantum coherence and entanglement between emissions. As Next-to-Next-to-Leading Order (NNLO) calculations become necessary, the computational cost of simulating high jet multiplicities is expected to grow significantly. There is a need for new algorithms that can naturally capture quantum correlations, specifically entanglement, inherent in parton showers to improve efficiency and physical consistency.

Methodology
The authors propose a modular, quantum-native approach to modeling QCD parton splitting, specifically for the pure-gluon channel ($g \to gg$). The methodology proceeds in three stages:

1. Analytic Derivation of Entanglement:
   Using quantum information tools, the authors analyze the $g \to gg$ scattering process. By defining the $R$-matrix from scattering amplitudes and tracing out color degrees of freedom, they derive the spin density matrix for the two outgoing gluons. They quantify the entanglement using concurrence ($C$), obtaining an analytic expression dependent on the momentum sharing fraction $z$:
   $$C_{QCD}(z) = \left( \frac{z(1-z)}{1-z(1-z)} \right)^2$$
   This function peaks at $z=0.5$ (equal momentum sharing) and vanishes as $z \to 0$ or $1$.

2. Quantum Circuit Construction:
   A two-qubit quantum circuit is designed to reproduce both the momentum fractions and the entanglement structure derived above.

   • Encoding: The momentum fractions $z$ and $1-z$ are encoded as the expectation values of the Pauli $\sigma_3$ operator for two qubits ($\langle \sigma_3 \rangle_A = z$, $\langle \sigma_3 \rangle_B = 1-z$).
   • Parameters: The circuit utilizes free parameters $\{\gamma_1, \gamma_2, \gamma_3\}$ constrained such that the concurrence of the circuit output, $C_{circuit}(\gamma_1, \gamma_3)$, matches $C_{QCD}(z)$.
   • Modularity: To model multi-prong structures (parton showers), the primitive circuit ($U_S$) is iterated. A CNOT gate connects the qubit carrying the highest momentum fraction from a previous splitting to the input of the next splitting block. This ensures kinematic consistency (momentum conservation) and allows the circuit to generate $n$-prong distributions by successive application.

3. Data-Driven Calibration:
   Instead of relying solely on perturbative theory, the circuit parameters are calibrated against experimental data. The authors use the AspenOpenJets dataset (derived from 2016 CMS Open Data).

   • Jets are reconstructed and declustered to extract momentum fractions for two-prong, three-prong, and four-prong topologies.
   • For each observed momentum fraction $z_i$, the corresponding circuit parameters $(\gamma_{1,i}, \gamma_{3,i})$ are solved numerically to satisfy both the momentum constraint and the entanglement matching condition.

Key Contributions

• Analytic Entanglement Formula: Derivation of the analytic concurrence $C_{QCD}(z)$ for the $g \to gg$ splitting, linking QCD dynamics directly to a quantum information metric.
• Composable Quantum Primitive: Construction of a two-qubit circuit that encodes both momentum sharing and QCD-predicted entanglement. The circuit is designed to be composable, allowing the simulation of complex shower evolution by iterating the primitive on the highest-momentum qubit.
• Experimental Validation: Successful calibration of the circuit parameters using LHC jet substructure data. The framework demonstrates that a physics-informed ansatz can reproduce empirical momentum distributions.
• Hardware Execution: Execution of the three-prong circuit on superconducting quantum hardware (ibm_Marrakesh). The results show consistency with noiseless simulations and experimental data after standard quality cuts and minor post-processing.

Results

• Two-Prong Validation: The calibrated circuit successfully reproduces the empirical momentum fraction distribution of two-prong jets while maintaining concurrence values consistent with the QCD prediction.
• Multi-Prong Generalization: By iterating the two-prong primitive, the authors generated three- and four-prong momentum fraction distributions. These distributions showed excellent agreement with the AspenOpenJets dataset without re-calibrating parameters for higher prong counts.
• Hardware Performance: On real quantum hardware, the three-prong results matched the ideal simulation and experimental data after excluding unphysical outputs (e.g., negative momentum fractions) and applying a shift to compensate for noise-induced deviations in values near $z=1$. The low qubit count and shallow circuit depth were cited as key factors enabling this success.
• Dominant Configuration: The study confirmed that the dominant physical configuration involves successive splittings of the highest-momentum particle. Configurations where lower-momentum particles split again contributed negligibly to the final distribution in this specific setup.

Significance and Claims
The paper positions this work as a proof-of-concept and a foundational step toward a complete quantum algorithm for parton showers.

• Physics-Informed Ansatz: The authors claim that by basing the quantum circuit on the physical splitting function and entanglement structure, they provide a structured, physics-informed starting point for Quantum Machine Learning (QML) applications. This approach is argued to be superior to generic variational circuits.
• Quantum-Native Simulation: The work demonstrates a concrete framework for "quantum-native" parton-shower modules that encode quantum correlations at the level of splitting dynamics, a feature absent in classical Markov Chain Monte Carlo generators.
• Future Outlook: The authors modestly note that while the current circuit handles perturbative splittings, a complete formulation requires integrating Sudakov factors (to handle virtuality evolution) and modeling non-perturbative hadronization. They suggest that future work could train circuit parameters directly on jet substructure observables to inherently capture non-perturbative corrections. The current work does not claim to replace existing generators but offers a new, entanglement-aware building block for future quantum algorithms in HEP.