TY - JOUR

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

AU - Shi, Yuan

AU - Xiao, Jianyuan

AU - Qin, Hong

AU - Fisch, Nathaniel J.

N1 - Funding Information:
The authors are grateful to Qun Wang and Sebastian Meuren for helpful discussions. This research is supported by NNSA Grant No. DE-NA0002948 and DOE Research Grant No. DEAC02-09CH11466. J.X. is supported by the National Magnetic Confinement Fusion Energy Research Project (2015GB111003, 2014GB124005), National Natural Science Foundation of China (NSFC-11575185, 11575186, 11305171), JSPS-NRF-NSFC A3 Foresight Program (NSFC-11261140328), Chinese Scholar Council (201506340103), Key Research Program of Frontier Sciences CAS (QYZDB-SSW-SYS004), and GeoAlgorithmic Plasma Simulator (GAPS) Project. APPENDIX A:
Funding Information:
This research is supported by NNSA Grant No. DE-NA0002948 and DOE Research Grant No. DEAC02-09CH11466. J.X. is supported by the National Magnetic Confinement Fusion Energy Research Project (2015GB111003, 2014GB124005), National Natural Science Foundation of China (NSFC-11575185, 11575186, 11305171), JSPS-NRF-NSFC A3 Foresight Program (NSFC-11261140328), Chinese Scholar Council (201506340103), Key Research Program of Frontier Sciences CAS (QYZDB-SSW-SYS004), and GeoAlgorithmic Plasma Simulator (GAPS) Project.
Publisher Copyright:
© 2018 American Physical Society.

PY - 2018/5/9

Y1 - 2018/5/9

N2 - Real-time lattice quantum electrodynamics (QED) provides a unique tool for simulating plasmas in the strong-field regime, where collective plasma scales are not well separated from relativistic-quantum scales. As a toy model, we study scalar QED, which describes self-consistent interactions between charged bosons and electromagnetic fields. To solve this model on a computer, we first discretize the scalar-QED action on a lattice, in a way that respects geometric structures of exterior calculus and U(1)-gauge symmetry. The lattice scalar QED can then be solved, in the classical-statistics regime, by advancing an ensemble of statistically equivalent initial conditions in time, using classical field equations obtained by extremizing the discrete action. To demonstrate the capability of our numerical scheme, we apply it to two example problems. The first example is the propagation of linear waves, where we recover analytic wave dispersion relations using numerical spectrum. The second example is an intense laser interacting with a one-dimensional plasma slab, where we demonstrate natural transition from wakefield acceleration to pair production when the wave amplitude exceeds the Schwinger threshold. Our real-time lattice scheme is fully explicit and respects local conservation laws, making it reliable for long-time dynamics. The algorithm is readily parallelized using domain decomposition, and the ensemble may be computed using quantum parallelism in the future.

AB - Real-time lattice quantum electrodynamics (QED) provides a unique tool for simulating plasmas in the strong-field regime, where collective plasma scales are not well separated from relativistic-quantum scales. As a toy model, we study scalar QED, which describes self-consistent interactions between charged bosons and electromagnetic fields. To solve this model on a computer, we first discretize the scalar-QED action on a lattice, in a way that respects geometric structures of exterior calculus and U(1)-gauge symmetry. The lattice scalar QED can then be solved, in the classical-statistics regime, by advancing an ensemble of statistically equivalent initial conditions in time, using classical field equations obtained by extremizing the discrete action. To demonstrate the capability of our numerical scheme, we apply it to two example problems. The first example is the propagation of linear waves, where we recover analytic wave dispersion relations using numerical spectrum. The second example is an intense laser interacting with a one-dimensional plasma slab, where we demonstrate natural transition from wakefield acceleration to pair production when the wave amplitude exceeds the Schwinger threshold. Our real-time lattice scheme is fully explicit and respects local conservation laws, making it reliable for long-time dynamics. The algorithm is readily parallelized using domain decomposition, and the ensemble may be computed using quantum parallelism in the future.

UR - http://www.scopus.com/inward/record.url?scp=85047014339&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85047014339&partnerID=8YFLogxK

U2 - 10.1103/PhysRevE.97.053206

DO - 10.1103/PhysRevE.97.053206

M3 - Article

C2 - 29906856

AN - SCOPUS:85047014339

VL - 97

JO - Physical Review E

JF - Physical Review E

SN - 2470-0045

IS - 5

M1 - 053206

ER -