Quantum state tomography, entanglement detection and Bell violation prospects in weak decays of massive particles

This paper presents a general method for reconstructing the spin density matrix of multi-particle systems from angular decay data using Bloch parameterization and Wigner-Weyl transforms, and applies it to Monte Carlo simulations of massive particle decays to propose measurements for detecting entanglement and Bell inequality violations.

The Gist

Imagine you are trying to figure out the exact shape and orientation of a spinning top, but you can't touch it, look at it directly, or stop it from spinning. All you can do is watch the tiny pebbles it throws off as it spins.

This is essentially what the physicists in this paper are doing, but instead of a toy top, they are studying the most fundamental particles in the universe: W and Z bosons (heavy particles that carry the weak nuclear force) and top quarks.

Here is a breakdown of their work using simple analogies.

1. The Problem: The "Black Box" Mystery

In the quantum world, particles have a property called spin. Think of spin like a compass needle pointing in a specific direction, but it can point in many directions at once until it is measured.

When these particles are created in a collider (like the Large Hadron Collider at CERN), they are in a "quantum state." This state is described by a Density Matrix.

• The Analogy: Imagine the Density Matrix is a secret recipe for a cake. You can't see the recipe, but you can taste the cake (the decay products). If you only taste a few crumbs, you might guess the flavor, but you can't be sure of the exact ingredients.
• The Goal: The authors want to reconstruct the entire secret recipe (the full Density Matrix) just by tasting the crumbs (the angles at which particles fly out).

2. The Old Way vs. The New Way

Previously, scientists mostly studied particles with "spin 1/2" (like electrons or top quarks). These are simple, like a coin that can only be Heads or Tails.

• The New Challenge: The authors are looking at particles with Spin 1 (like W and Z bosons). These are more complex, like a die that can land on 1, 2, 3, 4, 5, or 6. In quantum computing terms, electrons are "qubits" (2 states), but these bosons are "qutrits" (3 states).
• The Difficulty: You can't just look at a coin to see if it's heads or tails; you need a more complex map to understand a spinning die. The old methods didn't work well for these heavier, more complex particles.

3. The Solution: The "Magic Map" (Wigner-Weyl Transform)

The authors developed a new mathematical toolkit to solve this. They used a method called Quantum State Tomography.

• The Analogy: Imagine you are in a dark room with a spinning, glowing object. You can't see the object, but you have a special camera that takes pictures of the light it casts on the walls (the angular decay distributions).
• The Innovation: The authors created a "Magic Map" (based on something called the Generalised Gell-Mann parameterisation). This map translates the shadows on the wall (the data from the detectors) directly back into the 3D shape of the object (the Density Matrix).
• The "Projective" vs. "Non-Projective" Trick:
  • Sometimes, the particle decays in a way that gives a very clear, sharp shadow (like a coin landing flat). This is easy to read.
  • Sometimes, the decay is "fuzzy" (like a spinning top wobbling). The authors figured out how to read even these fuzzy shadows to get the exact recipe.

4. Why Does This Matter? (The "Spooky" Stuff)

Once you have the full recipe (the Density Matrix), you can check for two very strange quantum phenomena:

A. Entanglement (The "Twin Telepathy")

• The Concept: Two particles can be "entangled," meaning they share a single existence. If you measure one, the other instantly knows what to do, no matter how far apart they are. Einstein called this "spooky action at a distance."
• The Paper's Finding: The authors simulated collisions where pairs of W and Z bosons are created. They used their new map to prove that these pairs are highly entangled. It's as if two dice rolled in different rooms always land on the same number, not by chance, but because they are linked.

B. Bell Inequalities (The "Reality Check")

• The Concept: Albert Einstein believed the universe was "local" (things only affect their immediate surroundings) and "real" (particles have definite properties before you look at them). Quantum mechanics says no.
• The Test: There is a mathematical rule (Bell's Inequality) that "local reality" must obey. If you break this rule, you prove the universe is quantum.
• The Paper's Finding: The authors showed that in specific scenarios (like when a Higgs boson decays into W or Z pairs), the particles break this rule. They violate the "local reality" limit, proving that the universe is indeed "spooky" and quantum.

5. The Simulation (The "Virtual Lab")

Since we can't perfectly isolate these particles in a real lab yet (due to noise and background interference), the authors ran massive computer simulations (Monte Carlo).

• They simulated millions of collisions.
• They applied their "Magic Map" to the fake data.
• The Result: The map worked perfectly. It successfully reconstructed the quantum states and confirmed that entanglement and Bell violations are detectable in these heavy particle collisions.

Summary

This paper is a user manual for reading the quantum mind of heavy particles.

The authors built a new mathematical lens that allows physicists to look at the debris left behind by decaying W and Z bosons and say, "Aha! I know exactly how these particles were spinning, and I can prove they were 'telepathically' linked to their partners."

This is a major step forward because it moves quantum experiments from simple "coins" (electrons) to complex "dice" (bosons), opening the door to testing the deepest foundations of reality using the heaviest particles we can create.

Technical Summary

Here is a detailed technical summary of the paper "Quantum state tomography, entanglement detection and Bell violation prospects in weak decays of massive particles" by Ashby-Pickering, Barr, and Wierzchucka.

1. Problem Statement

The paper addresses the challenge of performing Quantum State Tomography (QST) on systems of massive particles with spin $j \geq 1/2$ (specifically $j=1$ vector bosons like $W^\pm$ and $Z^0$) produced in high-energy collisions.

• Limitation of Existing Methods: While QST is well-established for spin-1/2 systems (qubits) via Bloch vector measurements, it is insufficient for higher-spin systems (qudits, where dimension $d = 2j+1$). For $j=1$ ($d=3$), the spin polarization vector alone does not fully characterize the density matrix $\rho$; $d^2-1 = 8$ real parameters are required.
• Experimental Gap: There is a lack of a general, analytical framework to reconstruct the full spin density matrix of multipartite systems (e.g., $W^+W^-$, $ZZ$, $WZ$) directly from angular decay distributions, particularly for non-projective decays where the daughter particles do not act as perfect spin analyzers.
• Objective: To develop a general method to reconstruct the spin density matrix of arbitrary spin particles from decay data, enabling the calculation of quantum information metrics such as entanglement and Bell inequality violations in collider physics.

2. Methodology

The authors propose a unified framework combining Generalised Gell-Mann (GGM) parameterization with the Wigner-Weyl formalism.

A. Generalised Gell-Mann Parameterization

• The density matrix $\rho$ for a system of dimension $d$ is expanded in a basis of traceless, Hermitian Gell-Mann operators $\lambda^{(d)}_i$ (generalizing Pauli matrices for $d=2$ and Gell-Mann matrices for $d=3$).
• $\rho = \frac{1}{d}I_d + \sum_{i=1}^{d^2-1} a_i \lambda^{(d)}_i$.
• The goal is to determine the real coefficients $a_i$ (the generalized Bloch vector) from experimental data.

B. Wigner-Weyl Formalism for Spin

The paper utilizes the invertible mapping between operators in Hilbert space and functions on the sphere ($S^2$):

1. Forward Mapping (Q-Symbols): The probability density function (PDF) of decay products in direction $\hat{n}$ is related to the density matrix via the Wigner Q symbol (or generalized Q symbol for non-projective cases):
   $$p(\hat{n}; \rho) \propto \text{tr}(\rho \Pi_{\hat{n}}) = \Phi^Q_\rho(\hat{n})$$
   where $\Pi_{\hat{n}}$ represents the measurement operator associated with the decay direction.
2. Inverse Mapping (P-Symbols): To reconstruct $\rho$, the authors derive Wigner P symbols ($\Phi^P_i$) which are orthogonal to the Q symbols.
   $$a_i = \frac{1}{2} \int d\Omega \, p(\hat{n}; \rho) \, \Phi^P_i(\hat{n})$$
   This allows the extraction of each density matrix parameter $a_i$ as a simple angular average over the experimental decay distribution.

C. Handling Projective vs. Non-Projective Decays

• Projective Decays: For $W^\pm \to \ell \nu$ (with massless leptons), the decay acts as a von Neumann projective measurement. The authors derive explicit P and Q symbols for $d=3$ (spin-1).
• Non-Projective Decays: For cases where the daughter mass is non-negligible (e.g., heavy leptons) or for $Z \to \ell\ell$ (where chiral couplings are mixed), the measurement is a Positive Operator-Valued Measure (POVM). The authors generalize the formalism using Kraus operators to define generalized Q and P symbols that account for the specific decay dynamics (spin analyzing power $\kappa$).

D. Multipartite Systems

The method is extended to bipartite systems (e.g., $W^+W^-$). The joint density matrix is parameterized using tensor products of GGM matrices. The correlation parameters ($c_{ij}$) are extracted via double angular integrals over the decay products of both parents.

3. Key Contributions

1. General Analytical Framework: Developed a rigorous, analytical procedure to reconstruct the full spin density matrix for arbitrary spin $j$ and arbitrary numbers of parents, moving beyond the limitations of spin-1/2 tomography.
2. Derivation of Wigner Symbols: Explicitly calculated the Wigner P and Q symbols for:
   • Spin-1 ($d=3$, qutrits) for both projective ($W^\pm$) and non-projective ($Z$, massive lepton) decays.
   • Spin-3/2 ($d=4$) and Spin-1/2 ($d=2$) cases.
3. Integration with Quantum Information: Linked particle physics decay kinematics directly to quantum information concepts (Kraus operators, POVMs, concurrence, Bell operators), allowing the direct application of entanglement criteria to collider data.
4. Monte Carlo Validation: Performed extensive simulations of LHC processes ($pp \to W^+W^-, ZZ, WZ, H \to WW^*, H \to ZZ^*$) to demonstrate the feasibility of the method.

4. Results

The authors applied their method to Monte Carlo simulations of various diboson systems at $\sqrt{s} = 13$ TeV:

• Tomography Reconstruction: Successfully reconstructed the full $3 \times 3$(qutrit-qutrit) and$2 \times 2$ (qubit-qubit) density matrices. The reconstructed parameters matched theoretical expectations for idealized singlet states and Higgs decay scenarios.
• Entanglement Detection:
  • Used a lower bound on the square of the concurrence ($c^2_{MB}$) to detect entanglement.
  • Higgs Decays: Found strong evidence of entanglement in $H \to WW^*$ and $H \to ZZ^*$ decays ($c^2_{MB} > 0$).
  • Continuum Production: Inclusive $pp \to W^+W^-$ showed negative bounds (inconclusive for this specific measure), but entanglement was detected in specific kinematic regions (high invariant mass, specific angles).
  • Signal Purity: Demonstrated that entanglement can be detected in mixed signal/background states provided the signal fraction $\alpha \gtrsim 0.55$ for $WW^*$ and $\alpha \gtrsim 0.5$ for $ZZ^*$.
• Bell Inequality Violation:
  • Tested the CGLMP inequality (generalized for qutrits) using the Bell operator $B_{CGLMP}$.
  • Higgs Decays: Simulations of $H \to WW^*$ and $H \to ZZ^*$ showed violations of the local realist bound ($I_3 > 2$), with expectation values reaching up to $\approx 2.87$ (close to the quantum maximum).
  • Continuum: Standard $pp \to W^+W^-$ and $WZ$ continuum did not show violation with the standard operator, but $ZZ$ continuum showed potential violation in high-mass, forward-angle regions.
  • Note: The authors note that unphysically large values in some $ZZ^*$ simulations are likely due to higher-order perturbative corrections and kinematic selection effects not fully modeled in leading-order simulations.

5. Significance

• Foundational Physics: This work provides a concrete pathway to test the foundations of quantum mechanics (non-locality, entanglement) in the high-energy regime of massive vector bosons, a domain previously unexplored experimentally.
• Experimental Feasibility: It demonstrates that with current LHC luminosity and future colliders ($e^+e^-$, $\mu^+\mu^-$), it is possible to perform full quantum state tomography on diboson systems.
• New Observables: The paper establishes a set of new observables (GGM parameters, concurrence bounds, Bell operator expectations) that can be measured to probe anomalous couplings and the nature of the Higgs boson.
• Methodological Bridge: It successfully bridges the gap between the theoretical formalism of quantum information theory and the practical analysis of particle physics decay data, offering a template for future studies of quantum effects in hadron colliders.

In conclusion, the paper establishes that massive vector bosons produced in weak decays are highly entangled and capable of violating Bell inequalities, provided that full quantum state tomography is performed using the proposed Wigner-Weyl/GGM methodology.