QEDtool: A Python package for numerical quantum information in quantum electrodynamics

QEDtool is an object-oriented Python package designed to perform numerical quantum electrodynamics calculations by evaluating Feynman amplitudes in the momentum-helicity basis to reconstruct and quantify full state correlations, polarization, and entanglement for arbitrary initial scattering states across different inertial frames.

The Gist

Imagine you are a director filming a high-speed action movie. The actors are subatomic particles (like electrons and photons), and the scene is a violent collision at nearly the speed of light. In the old days, if you wanted to know how these actors would behave after the crash, you had to do the math on a chalkboard for weeks, and you could only get the "average" result.

QEDtool is a new, open-source "special effects engine" (a Python software package) that lets scientists simulate these collisions instantly and, more importantly, see the secret relationships between the actors.

Here is a breakdown of what this paper is about, using simple analogies:

1. The Problem: The "Quantum Entanglement" Mystery

In the quantum world, particles can be "entangled." Think of this like a pair of magical dice. No matter how far apart you roll them, if one lands on a 6, the other instantly knows and lands on a 1. They are linked in a way that classical physics can't explain.

Scientists are now trying to use these entangled particles for new technologies, like super-precise medical scanners (Quantum PET) or ultra-sharp microscopes (Quantum Lithography). But there's a catch: these particles are moving at relativistic speeds (close to the speed of light).

When things move that fast, the rules of space and time change (thanks to Einstein). A "left-handed" particle in one frame of reference might look "right-handed" to someone moving past it. Calculating how these "magical links" (entanglement) survive or change when particles zoom around at near-light speed is incredibly hard. It's like trying to predict how a spinning top behaves while the room itself is twisting and stretching.

2. The Solution: QEDtool (The Digital Lab)

The authors built QEDtool, a software toolkit that acts as a digital laboratory. Instead of doing the math by hand, you tell the computer:

• "Here are two particles entering the room."
• "Here is how they are spinning (polarization)."
• "Here is the speed of the room."

The software then:

1. Calculates the Crash: It uses the rules of Quantum Electrodynamics (QED)—the physics of light and matter—to figure out exactly how they bounce off each other.
2. Reconstructs the Aftermath: It doesn't just tell you where they went; it tells you their quantum state. It answers: "Are they still entangled? How strong is the link? Did the spin get mixed up?"
3. Changes the Camera Angle: It can instantly "boost" the view to a different speed or angle (Lorentz transformation) to see how the entanglement looks to a different observer.

3. How It Works (The "Lego" Analogy)

Think of the software as a set of Lego bricks:

• The Bricks: The basic building blocks are the particles (electrons, photons) and their "handedness" (spin).
• The Instructions: The software uses "Feynman diagrams" (which are like flowcharts for particle collisions) to snap the bricks together.
• The Output: Once the collision is simulated, the software builds a "report card" for the particles. It calculates:
  • Concurrence: A score from 0 to 1. If it's 1, the particles are perfectly entangled (maximum magic link). If it's 0, they are just regular, unlinked particles.
  • Stokes Parameters: A way to describe the "color" or "direction" of their spin, similar to how polarized sunglasses filter light.

4. Why This Matters

The paper shows that this tool is already working. The authors ran simulations of electron-positron annihilation (where matter and antimatter crash and turn into pure light).

• They found that at certain speeds, the entanglement is very strong.
• They found that if you speed up the observer (change the reference frame), the entanglement doesn't disappear, but it changes shape and looks different.
• They proved the software is accurate by comparing its results to known mathematical formulas, finding they matched perfectly (down to the 16th decimal place!).

The Bottom Line

QEDtool is a translator. It translates the incredibly complex, abstract math of high-energy physics into a user-friendly program. It allows scientists to ask, "What happens to the quantum link between these particles if we zoom past them at 90% the speed of light?" and get an answer in seconds rather than months.

This is a crucial step for building the next generation of quantum technologies that rely on particles moving at extreme speeds, ensuring we can harness the power of the quantum world even when it's moving at relativistic speeds.

Technical Summary

1. Problem Statement

The paper addresses the growing intersection between Quantum Technologies (e.g., quantum lithography, entangled positron emission tomography) and Relativistic Quantum Field Theory (QFT).

• The Challenge: Many emerging quantum technologies operate at relativistic energy scales where electromagnetic interactions dominate. Describing these systems requires Quantum Electrodynamics (QED). However, analytically determining tree-level cross-sections, spin correlations, and entanglement for arbitrary initial multi-particle states (pure or mixed) and their subsequent Lorentz transformations is mathematically complex, time-consuming, and prone to error.
• The Gap: While standard QFT tools exist for calculating scattering cross-sections, there is a lack of accessible, automated software specifically designed to reconstruct the full quantum state in the helicity basis, quantify entanglement, and analyze polarization correlations under relativistic conditions.

2. Methodology

The authors developed QEDtool, an object-oriented Python package designed to perform numerical QED calculations. The methodology is built upon the following theoretical and computational pillars:

• Framework: The package operates within a relativistic framework using natural units ($\hbar = c = \epsilon_0 = 1$). It utilizes the momentum-helicity basis to evaluate Feynman amplitudes.
• State Representation:
  • Initial States: Users can define both pure and mixed initial scattering states in polarization space. Mixed states are represented by density operators ($\rho$) constructed from convex combinations of pure states.
  • Scattering Formalism: The package uses the perturbative S-matrix formalism. It calculates transition amplitudes from pre-scattering asymptotes to post-scattering states using Feynman rules (propagators, vertices, and external line spinors/polarizations).
  • Final State Reconstruction: By applying the scattering operator to the initial density matrix and projecting onto specific momentum-helicity eigenstates, the package reconstructs the final density matrix $\rho^{(out)}$.
• Lorentz Transformations:
  • The package implements active Lorentz transformations for 4-vectors, Dirac spinors, and quantum states.
  • Crucially, it handles the mixing of helicity eigenstates for massive fermions (electrons/positrons) under boosts. While photon helicity is Lorentz invariant, massive spin-1/2 particles undergo Wigner rotations, causing their helicity states to mix. The package calculates these mixing coefficients (derived in Appendix A) to accurately transform quantum states between reference frames.
• Quantum Information Metrics:
  • Entanglement: Calculates the Concurrence ($C$) for two-qubit systems (e.g., photon pairs or electron-positron pairs) using the spin-flipped conjugate method (Wootters formula).
  • Polarization: Computes Stokes parameters ($S_{\mu}$) and the Degree of Polarization ($P^{(1)}$ for single particles, $P^{(2)}$ for pairs) to characterize the polarization state and correlations.

3. Key Contributions

• Novel Software Package: Introduction of QEDtool, the first Python package dedicated to numerical QED calculations focused on quantum information properties (entanglement and correlations) rather than just cross-sections.
• Hierarchical Architecture: The package is designed hierarchically:
  • Low-level: Customizable primitives for 4-vectors, Dirac spinors, and Green's functions.
  • High-level: Pre-packaged functions for standard scattering processes (e.g., Bhabha scattering, $e^+e^- \to 2\gamma$) and automatic calculation of quantum information metrics.
• Mixed State Handling: Unlike many analytical approaches that assume pure states, QEDtool natively supports mixed initial states (density matrices), allowing for the simulation of realistic, noisy, or thermal initial conditions.
• Relativistic Entanglement Analysis: The package explicitly demonstrates how Lorentz boosts affect entanglement and polarization correlations, a critical factor for high-energy quantum technologies where detectors may be in different inertial frames.

4. Results and Demonstrations

The paper validates QEDtool through several examples and comparisons:

• Validation against Analytical Theory:
  • The package's calculation of the differential scattering probability for unpolarized $e^+e^- \to 2\gamma$ annihilation was compared against standard QFT literature (Eq. 34). The results matched with an absolute difference of $\sim 10^{-16}$, confirming numerical accuracy.
• Electron-Positron Annihilation ($e^+e^- \to 2\gamma$):
  • Simulations showed that as collision momentum increases, the concurrence (entanglement) generally decreases, and the differential cross-section shifts toward forward/backward scattering.
  • The package successfully mapped the Stokes parameters and degree of polarization across solid angles, revealing that forward and backward scattered pairs remain circularly polarized.
• Bhabha Scattering ($e^+e^- \to e^+e^-$):
  • Using the standard_scattering function, the authors analyzed a mixed initial state. They observed that for specific kinematic conditions (backscattering at $\approx 400$ keV), the system approaches a maximally entangled Bell state ($|\Psi^-\rangle$), with concurrence $C \approx 1$.
• Lorentz Boost Effects:
  • A boost was applied to the annihilation process. The results demonstrated that while Feynman amplitudes are Lorentz invariant, the quantum correlations (entanglement and Stokes parameters) change significantly when viewed from a boosted frame due to the Wigner rotation of the helicity basis.

5. Significance and Future Outlook

• Technological Impact: QEDtool provides a necessary tool for the design and analysis of relativistic quantum technologies, such as Quantum Entangled PET (QE-PET) and Quantum Lithography, where understanding how entanglement survives scattering and frame transformations is vital.
• Educational Value: The package serves as an educational resource for learning the connection between QED scattering theory and quantum information concepts.
• Limitations and Future Work:
  • Currently, the package handles tree-level (leading order) processes only; loop corrections and renormalization are not included.
  • It is limited to electromagnetic interactions; the weak interaction is not yet implemented.
  • Initial states are currently restricted to momentum eigenstates; future versions aim to support general momentum wave functions.
  • The authors plan to extend the tensor classes to handle more general Lorentzian tensors and spin tensors.

In conclusion, QEDtool bridges the gap between high-energy particle physics and quantum information science, offering a robust, open-source (MIT license) platform for numerically exploring relativistic quantum correlations.