Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSimulations of relativistic-quantum plasmas using real-time lattice scalar QED
52
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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 1D 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.
rate research
Read More
We present a scheme for simulating relativistic quantum physics in circuit quantum electrodynamics. By using three classical microwave drives, we show that a superconducting qubit strongly-coupled to a resonator field mode can be used to simulate the dynamics of the Dirac equation and Klein paradox in all regimes. Using the same setup we also propose the implementation of the Foldy-Wouthuysen canonical transformation, after which the time derivative of the position operator becomes a constant of the motion.
Collisionless plasmas, mostly present in astrophysical and space environments, often require a kinetic treatment as given by the Vlasov equation. Unfortunately, the six-dimensional Vlasov equation can only be solved on very small parts of the considered spatial domain. However, in some cases, e.g. magnetic reconnection, it is sufficient to solve the Vlasov equation in a localized domain and solve the remaining domain by appropriate fluid models. In this paper, we describe a hierarchical treatment of collisionless plasmas in the following way. On the finest level of description, the Vlasov equation is solved both for ions and electrons. The next courser description treats electrons with a 10-moment fluid model incorporating a simplified treatment of Landau damping. At the boundary between the electron kinetic and fluid region, the central question is how the fluid moments influence the electron distribution function. On the next coarser level of description the ions are treated by an 10-moment fluid model as well. It may turn out that in some spatial regions far away from the reconnection zone the temperature tensor in the 10-moment description is nearly isotopic. In this case it is even possible to switch to a 5-moment description. This change can be done separately for ions and electrons. To test this multiphysics approach, we apply this full physics-adaptive simulations to the Geospace Environmental Modeling (GEM) challenge of magnetic reconnection.
Magnetic confinement fusion reactors suffer severely from heat and particle losses through turbulent transport, which has inspired the construction of ever larger and more expensive reactors. Numerical simulations are vital to their design and operation, but particle collisions are too infrequent for fluid descriptions to be valid. Instead, strongly magnetised fusion plasmas are described by the gyrokinetic equations, a nonlinear integro-differential system for evolving the particle distribution functions in a five-dimensional position and velocity space, and the consequent electromagnetic field. Due to the high dimensionality, simulations of small reactor sections require hundreds of thousands of CPU hours on High Performance Computing platforms. We develop a Hankel-Hermite spectral representation for velocity space that exploits structural features of the gyrokinetic system. The representation exactly conserves discrete free energy in the absence of explicit dissipation, while our Hermite hypercollision operator captures Landau damping with few variables. Calculation of the electromagnetic fields becomes purely local, eliminating inter-processor communication in, and vastly accelerating, searches for linear instabilities. We implement these ideas in SpectroGK, an efficient parallel code. Turbulent fusion plasmas may dissipate free energy through linear phase mixing to fine scales in velocity space, as in Landau damping, or through a nonlinear cascade to fine scales in physical space, as in hydrodynamic turbulence. Using SpectroGK to study saturated electrostatic drift-kinetic turbulence, we find that the nonlinear cascade suppresses linear phase mixing at energetically-dominant scales, so the turbulence is fluid-like. We use this observation to derive Fourier-Hermite spectra for the electrostatic potential and distribution function, and confirm these spectra with simulations.
We show how to implement a Rydberg-atom quantum simulator to study the non-equilibrium dynamics of an Abelian (1+1)-D lattice gauge theory. The implementation locally codifies the degrees of freedom of a $mathbf{Z}_3$ gauge field, once the matter field is integrated out by means of the Gauss local symmetries. The quantum simulator scheme is based on current available technology and scalable to considerable lattice sizes. It allows, within experimentally reachable regimes, to explore different string dynamics and to infer information about the Schwinger $U(1)$ model.
We consider the reflection of relativistically strong radiation from plasma and identify the physical origin of the electrons tendency to form a thin sheet, which maintains its localisation throughout its motion. Thereby we justify the principle of the relativistic electronic spring (RES) proposed in [A. Gonoskov et al. PRE 84, 046403 (2011)]. Using the RES principle we derive a closed set of differential equations that describe the reflection of radiation with arbitrary variation of polarization and intensity from plasma with arbitrary density profile for arbitrary angle of incidence. PIC simulations show that the theory captures the essence of the plasma dynamics. In particular, it can be applied for the studies of plasma heating and surface high-harmonic generation with intense lasers.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
