Few-particle lepton bound states in variational approach

This paper calculates the ground-state energy levels and hyperfine structure of three- and four-particle lepton bound states in quantum electrodynamics using a variational method with Gaussian basis functions that accounts for pairwise spin-spin interactions.

The Gist

Imagine the universe as a giant, cosmic dance floor. Usually, we think of particles like electrons and protons as solo dancers or simple pairs holding hands. But what happens when four of them try to dance together in a tight circle? Do they stick together, or do they fly apart?

This paper is a detailed mathematical study of exactly that: four-particle "dances" made of light particles called leptons.

Here is the breakdown of the research in simple terms:

1. The Cast of Characters

The scientists are studying "exotic" atoms. You know normal atoms (like Hydrogen) have a nucleus and one electron. These new systems are made entirely of leptons (lightweight particles like electrons, positrons, and muons).

• Positronium: An electron and its anti-matter twin (a positron) holding hands.
• Muonium: A heavy electron (muon) and a positron.
• The Goal: They are looking at systems where three or four of these particles are bound together, like a tiny, fragile molecule made of pure light and anti-light.

2. The Problem: It's Hard to Predict the Dance

In physics, when you have two particles, it's relatively easy to calculate how they move and how much energy they have. But when you add a third or fourth particle, the math gets messy. They push and pull on each other in complex ways.

• The Challenge: How do you predict exactly how tightly these four particles will hug each other, and how much energy is needed to keep them together?

3. The Solution: The "Gaussian" Guessing Game

The authors used a method called the Variational Method. Think of this like trying to find the lowest point in a foggy valley.

• The Analogy: Imagine you are blindfolded in a valley and want to find the very bottom (the most stable energy state). You can't see the whole map. Instead, you take a guess at where the bottom is, check the height, and then adjust your guess to go lower. You keep doing this until you are sure you can't go any lower.
• The Tool: They used "Gaussian functions" as their guesses. In math, these are bell-shaped curves. By stacking hundreds of these curves together, they built a super-accurate model of the particles' positions.
• The Supercomputer: They wrote a computer program (using MATLAB) to run this "guessing game" thousands of times, refining the model until they found the most stable energy levels.

4. The Spin-Spin Interaction (The Magnetic Handshake)

Particles have a property called "spin," which acts like a tiny internal magnet.

• The Metaphor: Imagine the particles are little compass needles. If they point in the same direction, they repel; if they point opposite, they attract.
• The paper calculates how these tiny magnets affect the energy of the system. This is called Hyperfine Structure. It's a tiny tweak to the energy, but it's crucial for understanding exactly how the system behaves.

5. Why Compare Them to "Tetraquarks"?

The authors make a fascinating comparison.

• Tetraquarks: These are heavy particles made of four quarks (the building blocks of protons) held together by the "strong" nuclear force. They are like heavy, dense bricks glued together.
• Tetra-leptons: These are the light particles in this paper, held together by the "electromagnetic" force (like magnets).
• The Connection: Even though the forces are different (one is strong glue, the other is magnetic), the math of how four things stick together is surprisingly similar. By studying the light, easy-to-calculate lepton dances, the scientists hope to learn secrets about the heavy, difficult-to-study quark dances.

6. What Did They Find?

They calculated the exact "binding energy" (how tightly they hold hands) for several exotic molecules, including:

• Positronium Hydride (HPs): A proton, two electrons, and a positron.
• Muonium Hydride (HMu): A proton, a muon, and two electrons.
• True Muonium Molecules: Four muons dancing together.

Their results match what other scientists have found, but they did it with a new, highly efficient method. They also calculated the average distance between the particles, showing that these systems are indeed "molecules" where two neutral pairs are gently attracted to each other, like two magnets floating near each other.

The Big Picture

Why does this matter?

1. Testing the Rules: It tests our understanding of Quantum Electrodynamics (QED), the rulebook for how light and matter interact.
2. Future Tech: Understanding these states helps us imagine new materials or even new ways to store energy.
3. The "What If": It opens the door to studying even stranger things, like "resonance states" (particles that almost stick together but then fly apart), which could reveal new physics beyond our current theories.

In short: These scientists built a super-precise mathematical model to watch four tiny particles dance together. They found that even in the chaotic world of subatomic particles, there is a beautiful, predictable order to how they stick together, and this order might help us understand the heaviest, most mysterious particles in the universe.

Technical Summary

Here is a detailed technical summary of the paper "Few-particle lepton bound states in variational approach" by A. V. Eskin et al.

1. Problem Statement

The paper addresses the theoretical calculation of energy levels and structural properties for few-body lepton systems within Quantum Electrodynamics (QED). While two-particle systems (like positronium and muonium) and three-particle systems (like the positronium ion, $Ps^-$) have been extensively studied, four-particle bound states (such as the positronium molecule $Ps_2$, positronium hydride $HPs$, and muonium hydride $HMu$) present greater computational complexity.

The specific challenges addressed include:

• Calculating the ground state binding energies of four-particle systems with high precision.
• Accounting for the hyperfine structure arising from pairwise spin-spin interactions.
• Determining the internal spatial structure (mean square distances) of these exotic systems.
• Extending variational methods previously used for three-particle systems to the more complex four-particle Coulomb systems.

2. Methodology

The authors employ the Stochastic Variational Method (SVM) using Gaussian basis functions in the coordinate representation.

• Coordinate System: The problem is formulated using Jacobi coordinates ($\rho, \lambda, \sigma$) to describe the relative distances between particles in the center-of-mass frame. This reduces the four-body problem to three relative coordinates, simplifying the kinetic energy operator.
• Hamiltonian: The non-relativistic Hamiltonian includes the kinetic energy of the four particles and pairwise Coulomb interaction potentials.
• Wave Function: The trial wave function for the ground state is constructed as a linear combination of Gaussian functions:
  $$\Psi(\rho, \lambda, \sigma) = \mathcal{A} \sum_{I=1}^{K} C_I \exp\left[ -\frac{1}{2} \left( A_{11}\rho^2 + 2A_{12}\rho\lambda + A_{22}\lambda^2 + \dots \right) \right]$$
  where $C_I$ are linear parameters and $A_{ij}$ are non-linear variational parameters optimized stochastically. $\mathcal{A}$ is an operator ensuring the correct permutation symmetry (antisymmetry for identical fermions).
• Analytical Integration: A key feature of the work is the derivation of analytical expressions for matrix elements of the kinetic energy, potential energy (Coulomb terms), and the Dirac $\delta$-function (for hyperfine splitting). This avoids numerical integration errors and allows for high precision.
• Hyperfine Structure: The hyperfine splitting is calculated using a spin-spin interaction Hamiltonian. The authors explicitly construct spin wave functions for systems with identical fermions (e.g., two electrons in $HPs$) to satisfy the Pauli exclusion principle and calculate the expectation values of the contact terms $\langle \delta(r_{ij}) \rangle$.
• Numerical Implementation: The generalized eigenvalue problem ($HC = E\lambda BC$) is solved numerically using a custom MATLAB program, based on algorithms from previous works (Varga-Suzuki). The basis size is systematically increased (from ~200 to ~800 functions) to ensure convergence.

3. Key Contributions

• Extension to Four-Body Systems: The authors successfully extended their variational framework, previously limited to three particles, to four-particle Coulomb systems.
• Analytical Formulas: They provided explicit analytical formulas for the normalization constants, kinetic/potential energy matrix elements, and mean square distances in terms of the variational parameters and particle masses.
• Hyperfine Calculations: Detailed calculation of hyperfine splitting for specific systems like Positronium Hydride ($HPs$) and Muonium-Positronium ($MuPs$), including the derivation of spin wave functions for mixed lepton systems.
• Structural Analysis: Calculation of mean square inter-particle distances, offering insight into the "cluster" nature of these molecules (e.g., viewing $MuPs$ as a muonium atom interacting with a positronium atom).

4. Results

The study presents binding energies for several lepton bound states in electronic atomic units (e.a.u.), comparing them with existing literature.

Binding Energies (Ground State):

• $Ps_2$ (Positronium Molecule): $-0.515982$ e.a.u. (Agrees with $-0.515989$ from Ref [30]).
• $HPs$ (Positronium Hydride): $-0.788819$ e.a.u. (Agrees with $-0.788867$ from Ref [4]).
• $HMu$ (Muonium Hydride): $-1.148355$ e.a.u.
• $MuPs$ (Muonium-Positronium): $-0.786262$ e.a.u.
• $Mu_2$ (Muonium Molecule): $-1.140230$ e.a.u.
• $Ps^-$ (Positronium Ion): $-0.261954$ e.a.u.
• $Mu^-$ (Muonium Ion): $-0.525054$ e.a.u.
• $TMu_2$ (True Muonium Molecule): $-106.6603$ e.a.u. (Calculated for the first time in this context, consisting of $\mu^+\mu^-\mu^+\mu^-$).

Hyperfine Splitting:

• For $HPs$ (with total electron spin $S_{12}=0$), the hyperfine splitting is calculated as 2.692 MHz.
• For $MuPs$, the splitting is 22.458 MHz.

Structural Data:

• Mean square distances were calculated for $MuPs$. For example, the distance between the electron and the muon ($\delta_{12}$) is $2.80$a.u., while the distance between the two electrons ($\delta_{13}$) is$4.00$ a.u., confirming the molecule's structure as two weakly interacting neutral clusters.

5. Significance

• Validation of QED: The results show excellent agreement (within 0.04%) with other high-precision calculations, validating the stochastic variational method with Gaussian bases for four-body QED problems.
• Analogy to QCD: The paper draws a significant theoretical parallel between four-lepton bound states and tetraquarks (four-quark states) in Quantum Chromodynamics (QCD). Both systems involve short-range interactions (Coulomb vs. Color-Electric) and can be studied using similar variational techniques, potentially aiding the understanding of exotic hadrons.
• Experimental Guidance: The calculated energy levels and decay widths provide benchmarks for experimentalists. The paper highlights the potential for creating and studying $MuPs$ and $TMu_2$ using high-intensity muon beams, which could test the Standard Model in the lepton sector.
• Foundation for Corrections: The high-precision non-relativistic results serve as a necessary baseline for future inclusion of relativistic ($O(\alpha^2)$) and radiative corrections, which are essential for ultimate precision tests of QED.

In conclusion, this work provides a robust theoretical framework and precise numerical data for four-particle lepton systems, bridging the gap between simple two-body QED and complex multi-body interactions, with implications for both atomic physics and high-energy particle physics.