Electron-positron pair creation in a supercritical static asymmetric potential well

Basic Information

• Paper ID: 2510.11002
• Title: Electron-positron pair creation in a supercritical static asymmetric potential well
• Authors: Z. L. Li, A. R. Sun, J. H. Xia, J. X. Wu, Y. J. Li (China University of Mining and Technology, Beijing)
• Classification: quant-ph (Quantum Physics), hep-ph (High Energy Physics - Phenomenology)
• Publication Date: October 14, 2025
• Paper Link: https://arxiv.org/abs/2510.11002

Abstract

This paper investigates electron-positron pair creation in supercritical static asymmetric potential wells using computational quantum field theory. The potential well consists of a subcritical potential and a supercritical potential separated by a fixed distance. To explain the discrete peak structures observed in the positron energy spectrum, the authors derive analytical formulas determining the bound state positions in subcritical asymmetric potential wells and extend these formulas to supercritical asymmetric potential wells through two methods. One method not only predicts bound state positions but also provides particle pair production rates. The study reveals that the height of the subcritical potential can optimize the energy distribution of produced electrons, offering a new pathway for generating energy-concentrated high-energy electron beams in experiments.

Research Background and Motivation

Core Problem

This research addresses the mechanism of electron-positron pair creation in supercritical asymmetric potential wells, particularly the physical origin of discrete peak structures appearing in the positron energy spectrum.

Research Significance

1. Fundamental Physics: Schwinger pair creation is an important prediction of quantum electrodynamics (QED), with vacuum breakdown requiring a critical field strength of $E_{cr} = m^2c^3/(e\hbar) \sim 10^{16}$ V/cm
2. Experimental Challenges: Due to the extremely high field strength required, experimental verification remains challenging, making theoretical research crucial for understanding this phenomenon
3. Application Prospects: Provides new mechanisms for enhancing pair creation or lowering production thresholds, with significant implications for high-intensity laser experiments and astrophysical environments

Limitations of Existing Methods

1. For supercritical symmetric potential wells, transmission coefficient formulas can determine resonance state positions
2. However, supercritical asymmetric potential wells lack corresponding analytical expressions, making rigorous determination of resonance state positions difficult
3. Although complex scaling methods are powerful numerical tools for computing resonance states, analytical approaches remain insufficiently explored

Core Contributions

1. Derived analytical formulas for bound state energy levels in subcritical asymmetric potential wells, extended to supercritical cases through two methods
2. Established a fitting method: Based on analytical formulas, fitting the functional relationship between bound state energy and potential height, applicable to supercritical potential wells
3. Proposed complex energy plane extension method: Directly extending analytical formulas to the complex energy region to calculate resonance states
4. Discovered optimization effects of subcritical potential height: Can control the energy distribution of produced electrons, achieving energy concentration
5. Revealed synergistic enhancement mechanism: Pair production rates in supercritical asymmetric potential wells exceed the sum of individual potential contributions

Methodology Details

Theoretical Framework

The research is based on computational quantum field theory (CQFT), starting from the Dirac equation: $i\frac{\partial\hat{\psi}}{\partial t} = h\hat{\psi}$

where the simplified one-dimensional Hamiltonian is: $h(z) = c\sigma_1 p_z + \sigma_3 c^2 - V(z)$

Potential Well Model

An asymmetric potential well composed of two Sauter potentials is employed: $-V(z) = \frac{V_1}{2}[1+\tanh(z/w)] - \frac{V_2}{2}[1+\tanh(z+d/w)]$

When $w \to 0$, the potential well divides into three regions:

• Region I ($z < -d$): $-V(z) = 0$
• Region II ($-d \leq z \leq 0$): $-V(z) = -V_2$
• Region III ($z > 0$): $-V(z) = V_1 - V_2$

Derivation of Bound State Energy Level Formula

By solving the Dirac equation in each region and applying boundary continuity conditions, the analytical formula for bound state energy levels is obtained: $i\tan(p_2 d) = \frac{\gamma\tau + 1}{\gamma + \tau}$

where: $\gamma = \frac{p_1 c}{E + c^2} \frac{E + V_2 + c^2}{p_2 c}, \quad \tau = \frac{p_2 c}{E + V_2 + c^2} \frac{E + V_2 - V_1 + c^2}{p_3 c}$

Extension Methods

1. Fitting Method: Establishing linear fitting relationships between bound state energy and potential height for $V_1 > c^2$, extrapolating to the supercritical region
2. Complex Energy Extension: Directly extending the formula to the complex energy plane, with the real part giving resonance state energy and the imaginary part giving width

Experimental Setup

Computational Parameters

• Numerical box length: $L = 8-16$ a.u.
• Grid points: $N_z = 2048-8192$
• Potential width: $w = 0.075/c - 0.3/c$
• Potential well spacing: $d = 0.2$ a.u.
• Typical potential heights: $V_1 = 2.5c^2$, $V_2 = 0.25c^2$

Evaluation Metrics

1. Energy Distribution Function: $N(E,t) = \frac{L}{2\pi(dE/dp)}\sum_n |U_{pn}(t)|^2$
2. Particle Number Density: $\rho(z,t) = \langle vac|\hat{\psi}_e^\dagger(z,t)\hat{\psi}_e(z,t)|vac\rangle$
3. Total Particle Number: $N(t) = \int dz \rho(z,t)$

Verification Comparison

• Comparison of analytical formula results with CQFT numerical results
• Comparison of pair production effects under different potential height configurations
• Analysis of temporal evolution behavior

Experimental Results

Main Results

Bound State Prediction Accuracy

Comparison of resonance state energies calculated by two extension methods with CQFT results:

• Fitting Method: Maximum relative error < 0.1%
• Complex Energy Method: Maximum relative error < 0.13%

Energy Spectrum Characteristics

1. Positron Energy Spectrum: Exhibits distinct discrete peak structures corresponding to bound state-continuum spectrum overlap
2. Electron Energy Spectrum: Shows oscillatory structures with increasing $V_2$, with most concentrated energy at $V_2 = 0.25c^2$

Temporal Evolution Patterns

• Early stage: Both bound-continuum and continuum-continuum spectrum overlap regions contribute to pair production
• Long time: Due to Pauli exclusion principle, only continuum-continuum spectrum overlap continues producing particles

Synergistic Enhancement Effect

For potential well with $V_1 = 3.54c^2$, $V_2 = 0.95c^2$:

• Pair production rate in asymmetric potential well: 558
• Sum of production rates from symmetric potential well + step potential: 454
• Enhancement effect: Approximately 23% improvement

Optimization Findings

When subcritical potential height $V_2 = 0.25c^2$, electron energy distribution is most concentrated, with energy range from $c^2$ to $1.25c^2$.

Related Work

Main Research Directions

1. Schwinger Pair Creation Enhancement Mechanisms: Dynamic assistance, multi-slit interference effects, frequency chirping, etc.
2. Non-perturbative Phenomena Research: Stokes phenomena, particle self-bunching, Ramsey fringes, etc.
3. Vacuum Control Research: Manipulating vacuum decay using control fields, vacuum information transmission, etc.

Innovations of This Work

Compared to existing work, this paper is the first to:

• Provide analytical calculation methods for bound states in asymmetric potential wells
• Discover mechanisms for subcritical potential optimization of electron energy distribution
• Reveal synergistic enhancement effects in asymmetric potential wells

Conclusions and Discussion

Main Conclusions

1. Validity of Analytical Formulas: The derived bound state formulas and their two extension methods accurately predict resonance states in supercritical asymmetric potential wells
2. Physical Mechanism Confirmation: Discrete peaks in the positron energy spectrum indeed originate from resonance states formed by bound state and negative energy continuum spectrum overlap
3. Optimization Strategy: Adjusting subcritical potential height can optimize electron energy distribution, providing new control methods for experiments
4. Enhancement Mechanism: Synergistic effects in asymmetric potential wells can significantly increase pair production rates

Limitations

1. Small Peak Structures: The origin of small peak structures between discrete peaks in the positron energy spectrum requires further investigation
2. Analytical Method Limitations: Supercritical asymmetric potential wells lack transmission coefficient formulas similar to symmetric potential wells
3. Parameter Range: Research primarily focuses on specific potential well parameter configurations

Future Directions

1. Explore the influence of bound state interference on small peak structures
2. Develop transmission coefficient analytical methods for supercritical asymmetric potential wells
3. Apply complex scaling methods to calculate Hamiltonian eigenvalues
4. Extend to more complex potential field configurations

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Starting from the Dirac equation, the derivation is mathematically rigorous with clear physical interpretation
2. Methodological Innovation: Proposes two methods for extending bound state formulas, providing new insights for complex potential field problems
3. Practical Value: Discovered optimization and enhancement mechanisms provide guidance for experiments
4. Comprehensive Verification: Multiple parameter configurations and temporal evolution analyses verify theoretical predictions

Weaknesses

1. Applicable Scope: Primarily focuses on specific Sauter potential models; universality for other potential forms requires verification
2. Experimental Feasibility: While providing theoretical guidance, practical experimental implementation still faces technical challenges
3. Detailed Mechanisms: Microscopic mechanisms of certain phenomena (such as small peak structures) require deeper investigation

Impact

1. Academic Contribution: Provides new analytical tools for non-perturbative phenomena research in quantum field theory
2. Application Prospects: Provides theoretical support for high-intensity laser experiments and extreme condition physics research
3. Method Generalization: Analytical methods may be applicable to research on other resonance phenomena

Applicable Scenarios

1. Basic Research: Non-perturbative QED effects research
2. Laser Physics: Ultra-intense laser-matter interactions
3. Astrophysics: Particle production processes in extreme astrophysical environments
4. Accelerator Physics: Generation and control of high-energy particle beams

References

The paper cites 37 important references covering theoretical foundations of Schwinger pair creation, enhancement mechanism research, and computational method development, demonstrating comprehensive understanding and deep consideration of related research.

This research makes important contributions at the intersection of theoretical and computational physics, deepening understanding of quantum phenomena under extreme conditions while providing valuable theoretical guidance for related experimental research.