Digital Quantum Simulation of Nonequilibrium Dynamics in the Schwinger Model under a Strong External Electric Field

Technical Summary: Digital Quantum Simulation of Nonequilibrium Dynamics in the Schwinger Model under a Strong External Electric Field

Problem Statement
The paper addresses the challenge of simulating nonequilibrium dynamics in strong-field quantum electrodynamics (QED), specifically the Schwinger effect, where a strong electric field induces vacuum decay and particle-antiparticle pair production. Traditional perturbative methods fail to describe this nonperturbative, real-time phenomenon reliably. Furthermore, standard Euclidean path-integral Monte Carlo methods encounter the "sign problem" when dealing with real-time evolution or complex weights, leading to rapid growth in statistical errors. While the Hamiltonian formalism avoids the sign problem, it requires handling an exponentially large Hilbert space. The authors propose using digital quantum simulation to overcome these limitations by directly representing quantum many-body states and unitary time evolution on qubits. The specific focus is on the (1+1)-dimensional Schwinger model under a constant external electric field, investigating the finite-size nonequilibrium dynamics following an electric-field quench.

Methodology
The study employs a digital quantum simulation protocol based on the (1+1)-dimensional lattice Schwinger model with open boundary conditions. The methodology consists of three main stages:

1. Model Formulation and Qubit Mapping:

   • The continuum theory is discretized using the Kogut-Susskind staggered-fermion formalism.
   • Gauss's law is used to eliminate the gauge degrees of freedom, expressing the electric field on each link as a function of the matter charge distribution and the boundary field. This reduces the model to a purely fermionic Hamiltonian with long-range interactions.
   • The fermionic operators are mapped to Pauli operators via the Jordan-Wigner transformation, resulting in a qubit Hamiltonian composed of kinetic, mass, and electric-field energy terms.

2. State Preparation (VQE):

   • The initial state for the dynamics is the zero-field vacuum state ($\epsilon=0$).
   • This state is prepared using the Variational Quantum Eigensolver (VQE) with a hardware-efficient RealAmplitudes ansatz and the SLSQP classical optimizer.
   • To study the static properties of the system under varying external fields, the authors also use VQE to find the ground state and Variational Quantum Deflation (VQD) to estimate the first excited state energy.

3. Real-Time Evolution (Trotter-Suzuki):

   • The system is subjected to a sudden "quench" where an external electric field is switched on at $t=0$.
   • The real-time evolution is simulated using a second-order Trotter-Suzuki decomposition to approximate the time-evolution operator $e^{-iHt}$.
   • The simulation is performed on classical computers using the Qiskit framework, without including quantum noise, to benchmark the algorithm against Exact Diagonalization (ED).

Key Contributions and Results
The paper presents numerical results for lattice sizes $N=8$ and $N=12$ (with lattice units $a=m=g=1$), comparing the quantum simulation protocol against Exact Diagonalization (ED) benchmarks.

• Static Ground-State Properties:

  • The VQE-prepared zero-field ground state achieves high fidelity ($>99.99\%$) compared to ED results.
  • Scans over the external field strength $\epsilon$ reveal a step-like behavior in the chiral condensate and a vanishing energy gap between the ground and first excited states at specific critical field strengths. This indicates a field-induced reconstruction of the vacuum structure, corresponding to a first-order phase transition in finite volumes.
  • The critical field strength $\epsilon_c$ exhibits a clear finite-size dependence on $1/N$. Linear extrapolation yields an infinite-volume critical field of $\epsilon_c(\infty) \approx 0.469$, which the authors attribute to the influence of open boundary conditions reducing the energy gap.

• Nonequilibrium Dynamics:

  • Conservation Laws: The second-order Trotter evolution successfully preserves the total $U(1)$ charge and total energy (within numerical precision), validating the stability of the algorithm.
  • Charge Separation: Upon switching on the field, the spatial charge distribution shows periodic production and annihilation of charge-anticharge pairs near the open boundaries. The amplitude and period of these oscillations increase with the external field strength.
  • Energy Redistribution: The total electric-field energy $H_E(t)$ exhibits quasiperiodic behavior, reflecting energy exchange between the electric-field term and the fermionic sector. This is distinct from the simpler periodicity of the charge distribution, suggesting complex contributions from higher-energy modes.
  • Vacuum Fidelity: The vacuum-state fidelity $P_{vac}(t)$ decays as the system is driven away from the initial vacuum. In the early-time window ($0 \le t/a \le 1.0$), this decay follows an effective exponential law, $P_{vac}(t) \sim e^{-\gamma_{eff} t}$. The effective decay rate $\gamma_{eff}$ increases monotonically with the external field strength.
  • Finite-Size Effects: Comparisons between $N=8$ and $N=12$ systems confirm that the observed dynamical features (boundary charge separation, fidelity decay) are not numerical artifacts of small system sizes but persist in larger lattices, with larger systems showing faster deviation from the initial vacuum.

Significance and Claims
The authors claim that their work demonstrates the viability of combining VQE-based state preparation with digital real-time evolution (via Trotter-Suzuki decomposition) for studying nonequilibrium dynamics in strong-field lattice gauge theories.

• Validation: The protocol successfully reproduces the main features of field-induced phenomena—specifically boundary charge separation, vacuum-state fidelity decay, and quasiperiodic energy redistribution—when compared to exact diagonalization.
• Methodological Utility: The study serves as a benchmark for digital quantum simulation in high-energy physics, showing that explicit digital Trotter evolution can effectively capture complex strong-field dynamics in finite lattice systems.
• Future Outlook: The authors modestly suggest that this approach provides a foundation for future extensions to larger system sizes, improved variational ansatzes, implementation on noisy quantum hardware, and the study of time-dependent external fields. They do not claim immediate experimental realization on current noisy hardware but emphasize the algorithmic verification on classical simulators as a necessary step.