Simulations of relativistic quantum plasmas using real-time lattice scalar QED

Applicability of the Klein-Gordon equation for pair production in vacuum and plasma

In this paper, a phase-space description of electron-positron pair creation will be applied, based on a Wigner transformation of the Klein-Gordon equation. The resulting theory is similar in many respects to the equations from the Dirac-Heisenberg-Wigner formalism. However, in the former case, all physics related to particle spin is neglected. In the present paper we compare the pair-production rate in vacuum and plasmas, with and without spin effects, in order to evaluate the accuracy and applicability of the spinless approximation. It is found that for modest frequencies of the electromagnetic field the pair production rate of the Klein-Gordon theory is a good approximation to the Dirac theory, provided the matter density is small enough for Pauli blocking to be neglected, and a factor of 2 related to the difference in the vacuum energy density is compensated for.

Plasma dynamics and vacuum pair creation using the Dirac-Heisenberg-Wigner formalism

We derive a system of coupled partial differential equations for the equal-time Wigner function in an arbitrary strong electromagnetic field using the Dirac-Heisenberg-Wigner formalism. In the electrostatic limit, we present a system of four coupled partial differential equations, which are completed by Ampères law. This electrostatic system is further studied for two different cases. In the first case, we consider linearized wave propagation in a plasma accounting for the nonzero vacuum expectation values. We then derive the dispersion relation and compare it with well-known limiting cases. In the second case, we consider Schwinger pair production using the local density approximation to allow for analytical treatment. The dependence of the pair production rate on the perpendicular momentum is investigated and it turns out that the spread of the produced pairs along with perpendicular momentum depends on the strength of the applied electric field.

Variational schemes and geometric simulations for a hydrodynamic-electrodynamic model of surface plasmon polaritons

A class of variational schemes for the hydrodynamic-electrodynamic model of lossless free-electron gas in a quasineutral background is developed for high-quality simulations of surface plasmon polaritons. The Lagrangian density of lossless free-electron gas with a self-consistent electromagnetic field is established, and the dynamical equations with the associated constraints are obtained via a variational principle. Based on discrete exterior calculus, the action functional of this system is discretized and minimized to obtain the discrete dynamics. Newton-Raphson iteration and the biconjugate gradient stabilized method are equipped as a hybrid nonlinear-linear algebraic solver. Instead of discretizing the partial differential equations, the variational schemes have better numerical properties in secular simulations, as they preserve the discrete Lagrangian symplectic structure, gauge symmetry, and general energy-momentum density. Two numerical experiments were performed. The numerical results reproduce characteristic dispersion relations of bulk plasmons and surface plasmon polaritons, and the numerical errors of conserved quantities in all experiments are bounded by a small value after long-term simulations.