Automated computation of spin-density matrices and quantum observables for collider physics
Pith reviewed 2026-05-18 05:55 UTC · model grok-4.3
The pith
A framework automates the computation of production spin-density matrices for generic collider processes at tree level.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper establishes a systematic procedure that takes helicity amplitudes for generic processes and constructs event-by-event production spin-density matrices. These matrices are output in the LHE format with run information, supporting various final state dimensions and initial state configurations, so that quantum-information quantities can be evaluated directly in post-processing.
What carries the argument
Event-by-event production spin-density matrices constructed from helicity amplitudes, which preserve the full quantum information of the multi-particle final state for subsequent analysis.
If this is right
- Spin-density matrices for top-quark pair production and vector boson pair production can be computed and compared to existing results.
- Quantum measures including entanglement witnesses become accessible for processes with three or four top quarks.
- Configurable reference frames and initial state polarizations allow flexible studies of spin correlations.
- A library of routines computes D-coefficients, correlation matrices, and stabiliser-based magic measures for the final states.
Where Pith is reading between the lines
- This automation opens the door to scanning quantum properties across large samples of simulated events without manual amplitude calculations.
- Connections to quantum information theory could lead to new observables sensitive to new physics effects in high-energy collisions.
- The approach might be extended to loop-level calculations to study how quantum correlations change with higher-order corrections.
Load-bearing premise
The helicity amplitudes for the processes can be assembled into compact production spin-density matrices that retain all quantum information without errors in the automation.
What would settle it
Computing the spin-density matrix for top pair production in proton-proton collisions and finding it differs from the known analytical result in a specific matrix element would falsify the correctness of the assembly method.
read the original abstract
We present a fully automated framework to compute production spin-density matrices for generic collider processes at tree level within \textsc{MadGraph5\_aMC@NLO}. The method assembles helicity amplitudes into event-by-event production matrices. These are written to the LHE file in a compact form, together with run metadata, enabling direct post-processing of quantum observables. The implementation supports bi- and multipartite qubit and qutrit final states, configurable reference frames, and both polarised and unpolarised initial states. A companion, easy-to-extend library provides analysis routines to determine key quantum-information measures and witnesses. These include purity, concurrence, and entanglement of formation for qubits; Peres--Horodecki tests and negativity; spin-polarisation vectors and correlation matrices; $D$-coefficients; and stabiliser-based ``magic'' measures. As a result, multi-particle quantum correlations can be quantified systematically. We validate the implementation against known results for $t\bar t$ and $VV$ ($V=W^\pm,Z$) production in $pp$ and $e^+e^-$ collisions and in heavy-resonance decays. We then consider new applications and study quantum correlations in several LHC final states: $t\bar t W^\pm$, $tW^-$ vs.\ $t(\bar t\to W^- \bar b)$, and $t\bar t t$ vs. $t\bar t t\bar t$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a fully automated framework integrated into MadGraph5_aMC@NLO to compute production spin-density matrices at tree level for generic collider processes. Helicity amplitudes are assembled into event-by-event matrices that are written to LHE files together with metadata, enabling post-processing via a companion library to extract quantum observables and measures such as purity, concurrence, negativity, spin-polarisation vectors, correlation matrices, D-coefficients, and stabiliser-based magic measures. The implementation supports bi- and multipartite qubit/qutrit states with configurable reference frames and both polarised and unpolarised beams. Validation is reported against known ttbar and VV results, followed by applications to ttW, tW, ttt, and tttt final states.
Significance. If the matrix construction preserves quantum information without phase or frame errors, the work provides a practical tool for systematic studies of multi-particle entanglement and spin correlations at colliders. Integration with a standard generator and the provision of an extendable analysis library are clear strengths that could enable routine exploration of quantum-information quantities in complex LHC final states.
major comments (2)
- [Section 3] Section 3 (matrix assembly): the central step that forms ρ_{λ,λ′} ∝ ∑ A_λ A^*_{λ′} from MadGraph helicity amplitudes is described at a high level, but no explicit account is given of global-phase conventions or consistent boost/reference-frame definitions across three or more legs. This is load-bearing for the multipartite claim, as relative phases directly affect off-diagonal elements used in negativity or concurrence for ttt and tttt.
- [Section 5] Section 5 (applications): the new results for ttt and tttt are presented without an independent cross-check (e.g., comparison to a manually computed toy case or symbolic verification of a small subset of matrix elements). Validation is confined to two-particle systems (ttbar, VV) where external benchmarks exist; extending the same code path to three- and four-particle states therefore rests on untested assumptions about the automation layer.
minor comments (1)
- [Abstract] Abstract: the term 'D-coefficients' appears without a short definition or reference; a parenthetical clarification would help readers outside the spin-correlation community.
Simulated Author's Rebuttal
We thank the referee for the positive evaluation of our work and for the constructive comments, which have helped us improve the clarity and robustness of the manuscript. We respond to each major comment below.
read point-by-point responses
-
Referee: [Section 3] Section 3 (matrix assembly): the central step that forms ρ_{λ,λ′} ∝ ∑ A_λ A^*_{λ′} from MadGraph helicity amplitudes is described at a high level, but no explicit account is given of global-phase conventions or consistent boost/reference-frame definitions across three or more legs. This is load-bearing for the multipartite claim, as relative phases directly affect off-diagonal elements used in negativity or concurrence for ttt and tttt.
Authors: We agree that an explicit discussion of phase and frame conventions is necessary to support the multipartite results. MadGraph5_aMC@NLO computes helicity amplitudes with a fixed global-phase convention based on standard spinor and polarization-vector choices (consistent with the conventions in the MadGraph literature and the HELAS library). For processes with three or more final-state legs, all momenta are defined in a common lab frame, and boosts to individual particle rest frames are performed using the same Lorentz-transformation routines, preserving relative phases between helicity configurations. In the revised manuscript we have added a dedicated paragraph in Section 3 that spells out these conventions, includes a brief derivation of how the overall phase cancels in the density-matrix construction, and provides a worked example for a three-particle final state. revision: yes
-
Referee: [Section 5] Section 5 (applications): the new results for ttt and tttt are presented without an independent cross-check (e.g., comparison to a manually computed toy case or symbolic verification of a small subset of matrix elements). Validation is confined to two-particle systems (ttbar, VV) where external benchmarks exist; extending the same code path to three- and four-particle states therefore rests on untested assumptions about the automation layer.
Authors: We acknowledge that external analytic benchmarks for three- and four-particle spin-density matrices do not exist in the literature. The matrix-assembly algorithm itself is number-of-particles agnostic and was validated on two-body processes against known results. For the new applications we have added internal consistency tests (trace normalisation, hermiticity, and positive-semidefiniteness of the extracted matrices) that are reported in the revised Section 5. In addition, we now include a comparison against a manually computed toy model with a restricted set of helicity states, where the density-matrix elements agree to within numerical precision. While a complete symbolic verification of the full ttt and tttt processes is impractical, these checks address the core concern about the automation layer. revision: partial
Circularity Check
No circularity: automated assembly of helicity amplitudes into density matrices validated on external benchmarks
full rationale
The paper presents a software implementation that takes existing MadGraph helicity amplitudes as input and assembles them into standard production spin-density matrices using the conventional definition ρ ∝ ∑ A_λ A^*_λ'. This assembly is the direct computational realization of the quantum-mechanical definition rather than a derived prediction. Validation is performed against independently known results for ttbar and VV processes, with no fitted parameters, self-referential definitions, or load-bearing self-citations that reduce the central claim to its own inputs. The contribution is the automation layer and analysis library, which remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Tree-level perturbation theory suffices for the processes under study.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The method assembles helicity amplitudes into event-by-event production matrices... Rh1h2 ≡ 1/N ∑ M∗h1 Mh2
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
validation against known results for ttbar and VV production
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 3 Pith papers
-
Radiation effects on the entanglement of fermion pairs at colliders
Energetic radiation induces decoherence that significantly reduces entanglement in fermion pairs at colliders, with statistically significant signals observable in ttbar(g) at the LHC and tau pairs at Belle II.
-
Disentangling new physics with quantum entanglement in $t\bar{t}$ production at future lepton colliders
Entanglement in ttbar production at lepton colliders is typically reduced by scalar mediators but shows sizable deviations in U(1)B-L and Randall-Sundrum models, positioning quantum-information observables as probes f...
-
Quantum Tomography and Entanglement in Semi-Leptonic $h\to VV^*$ Decays at Higher Orders
Semi-leptonic h to VV* decays retain an effective two-qutrit description for quantum tomography and entanglement after including finite fermion masses and NLO corrections.
Reference graph
Works this paper leans on
-
[1]
Y. Afik and J. R. M. n. de Nova,Entanglement and quantum tomography with top quarks at the LHC,Eur. Phys. J. Plus136(2021) 907, [2003.02280]
-
[2]
M. Fabbrichesi, R. Floreanini and G. Panizzo,Testing Bell Inequalities at the LHC with Top-Quark Pairs,Phys. Rev. Lett.127(2021) 161801, [2102.11883]
-
[3]
C. Severi, C. D. E. Boschi, F. Maltoni and M. Sioli,Quantum tops at the LHC: from entanglement to Bell inequalities,Eur. Phys. J. C82(2022) 285, [2110.10112]. – 62 –
-
[4]
C. Severi and E. Vryonidou,Quantum entanglement and top spin correlations in SMEFT at higher orders,JHEP01(2023) 148, [2210.09330]
- [5]
-
[6]
Y. Afik and J. R. M. n. de Nova,Quantum information with top quarks in QCD,Quantum6 (2022) 820, [2203.05582]
- [7]
-
[8]
Y. Afik and J. R. M. n. de Nova,Quantum Discord and Steering in Top Quarks at the LHC, Phys. Rev. Lett.130(2023) 221801, [2209.03969]
- [9]
- [10]
- [11]
- [12]
- [13]
- [14]
- [15]
- [16]
- [17]
- [18]
-
[19]
F. Maltoni, C. Severi, S. Tentori and E. Vryonidou,Quantum detection of new physics in top-quark pair production at the LHC,JHEP03(2024) 099, [2401.08751]
-
[20]
F. Maltoni, C. Severi, S. Tentori and E. Vryonidou,Quantum tops at circular lepton colliders, JHEP09(2024) 001, [2404.08049]
-
[21]
Probing new physics in the top sector using quantum information,
R. Aoude, H. Banks, C. D. White and M. J. White,Probing new physics in the top sector using quantum information,2505.12522
-
[22]
M. Fabbrichesi, R. Floreanini and L. Marzola,Local vs. nonlocal entanglement in top-quark pairs at the LHC,2505.02902. – 63 –
- [23]
-
[24]
M. Fabbrichesi, R. Floreanini and E. Gabrielli,Constraining new physics in entangled two-qubit systems: top-quark, tau-lepton and photon pairs,Eur. Phys. J. C83(2023) 162, [2208.11723]
-
[25]
K. Ehatäht, M. Fabbrichesi, L. Marzola and C. Veelken,Probing entanglement and testing Bell inequality violation with e+e-→τ+τ- at Belle II,Phys. Rev. D109(2024) 032005, [2311.17555]
-
[26]
M. Fabbrichesi and L. Marzola,Dipole momenta and compositeness of theτlepton at Belle II, 2401.04449
-
[27]
M. Fabbrichesi and L. Marzola,Quantum tomography withτleptons at the FCC-ee: Entanglement, Bell inequality violation, sinθW, and anomalous couplings,Phys. Rev. D110 (2024) 076004, [2405.09201]
-
[28]
Y. Zhang, B.-H. Zhou, Q.-B. Liu, S. Li, S.-C. Hsu, T. Han et al.,Entanglement and Bell Nonlocality inτ +τ − at the LHC using Machine Learning for Neutrino Reconstruction, 2504.01496
work page internal anchor Pith review Pith/arXiv arXiv
- [29]
-
[30]
N. A. Tornqvist,Suggestion for Einstein-podolsky-rosen Experiments Using Reactions Like e+e− →Λ ¯Λ→π −pπ+¯p,Found. Phys.11(1981) 171–177
work page 1981
-
[31]
R. A. Bertlmann and W. Grimus,A Model for decoherence of entangled beauty,Phys. Rev. D 64(2001) 056004, [hep-ph/0101160]
work page internal anchor Pith review Pith/arXiv arXiv 2001
- [32]
- [33]
-
[34]
M. Fabbrichesi, R. Floreanini, E. Gabrielli and L. Marzola,Bell inequality is violated in B0→J/ψK*(892)0 decays,Phys. Rev. D109(2024) L031104, [2305.04982]
-
[36]
E. Gabrielli and L. Marzola,Entanglement and Bell Inequality Violation in B→ϕϕDecays, Symmetry16(2024) 1036, [2408.05010]
-
[37]
M. Fabbrichesi, R. Floreanini, E. Gabrielli and L. Marzola,Bell inequality is violated in charmonium decays,Phys. Rev. D110(2024) 053008, [2406.17772]
-
[38]
K. Cheng and B. Yan,Bell Inequality Violation of Light Quarks in Back-to-Back Dihadron Pair Production at Lepton Colliders,2501.03321
-
[39]
K. Chen, Z.-P. Xing and R. Zhu,Test of Bell Locality Violation in Flavor Entangled Neutral Meson Pair,2407.19242. – 64 –
-
[40]
M. Fabbrichesi, R. Floreanini, E. Gabrielli and L. Marzola,Quantum pions,2506.05464
- [41]
- [42]
- [43]
- [44]
-
[45]
R. Ashby-Pickering, A. J. Barr and A. Wierzchucka,Quantum state tomography, entanglement detection and Bell violation prospects in weak decays of massive particles,JHEP05(2023) 020, [2209.13990]
- [46]
- [47]
- [48]
- [49]
-
[50]
M. Fabbrichesi, R. Floreanini, E. Gabrielli and L. Marzola,Bell inequalities and quantum entanglement in weak gauge boson production at the LHC and future colliders,Eur. Phys. J. C 83(2023) 823, [2302.00683]
-
[51]
M. Fabbrichesi, R. Floreanini, E. Gabrielli and L. Marzola,Stringent bounds on HWW and HZZ anomalous couplings with quantum tomography at the LHC,JHEP09(2023) 195, [2304.02403]
- [52]
- [53]
-
[54]
Constraining New Physics withh→V VTomography,
M. Sullivan,Constraining New Physics withh→V VTomography,2410.10980
- [55]
- [56]
- [57]
- [58]
- [59]
- [60]
- [61]
- [62]
-
[63]
M. Fabbrichesi, R. Floreanini, E. Gabrielli and L. Marzola,Quantum contextuality of spin-1 massive particles,Phys. Rev. A112(2025) L020402, [2503.14587]
-
[64]
Quantum properties of H→VV ∗: precise predictions in the SM and sensitivity to new physics,
M. Del Gratta, F. Fabbri, P. Lamba, F. Maltoni and D. Pagani,Quantum properties of H→V V ∗: precise predictions in the SM and sensitivity to new physics,2504.03841
-
[65]
Decoherence effects in entangled fermion pairs at colliders
R. Aoude, A. J. Barr, F. Maltoni and L. Satrioni,Decoherence effects in entangled fermion pairs at colliders,2504.07030
work page internal anchor Pith review Pith/arXiv arXiv
- [66]
-
[67]
D. Gonçalves, A. Kaladharan, F. Krauss and A. Navarro,Quantum Entanglement is Quantum: ZZ Production at the LHC,2505.12125
-
[68]
Z-boson quantum tomography at next-to-leading order,
M. Del Gratta, F. Fabbri, M. Grossi, F. Maltoni, D. Pagani, G. Pelliccioli et al.,Z-boson quantum tomography at next-to-leading order,2509.20456
- [69]
-
[70]
Entanglement measures and their properties in quantum field theory
S. Hollands and K. Sanders,Entanglement measures and their properties in quantum field theory,1702.04924
-
[71]
Jozsa,Fidelity for Mixed Quantum States,J
R. Jozsa,Fidelity for Mixed Quantum States,J. Mod. Opt.41(1994) 2315–2323
work page 1994
-
[72]
M. B. Plenio,Logarithmic negativity: A full entanglement monotone that is not convex,Phys. Rev. Lett.95(Aug, 2005) 090503
work page 2005
-
[73]
H. Ollivier and W. H. Zurek,Introducing Quantum Discord,Phys. Rev. Lett.88(2001) 017901, [quant-ph/0105072]
work page internal anchor Pith review Pith/arXiv arXiv 2001
-
[74]
F. Mintert, M. KuÅŻ and A. Buchleitner,Concurrence of mixed multipartite quantum states, Physical Review Letters95(Dec., 2005)
work page 2005
-
[75]
F. Mintert and A. Buchleitner,Observable entanglement measure for mixed quantum states, Phys. Rev. Lett.98(Apr, 2007) 140505
work page 2007
-
[76]
F. Mintert and A. Buchleitner,Observable entanglement measure for mixed quantum states, Physical Review Letters98(Apr., 2007)
work page 2007
-
[77]
A computable measure of entanglement
G. Vidal and R. F. Werner,Computable measure of entanglement,Phys. Rev. A65(2002) 032314, [quant-ph/0102117]
work page internal anchor Pith review Pith/arXiv arXiv 2002
- [78]
-
[79]
M. Fabbrichesi, R. Floreanini and L. Marzola,About testing Bell locality at colliders, 2503.18535. – 66 –
-
[80]
M. Fabbrichesi, M. Low and L. Marzola,Trace distance between density matrices: A nifty tool in new-physics searches,Phys. Rev. D112(2025) 013003, [2501.03311]
-
[81]
M. Fabbrichesi, R. Floreanini, E. Gabrielli and L. Marzola,Tests of quantum contextuality in particle physics,Phys. Rev. D112(2025) 033005, [2504.12382]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.