Strong-field Quantum Electrodynamics Research Articles

Strong-field QED experiments using the BELLA PW laser dual beamlines

The Petawatt (PW) laser facility of the Berkeley Lab Laser Accelerator (BELLA) Center has recently commissioned its second laser pulse transport line. This new beamline can be operated in parallel with the first beamline and enables strong-field quantum electrodynamics (SF-QED) experiments at BELLA. In this paper, we present an overview of the upgraded BELLA PW facility with a SF-QED experimental layout in which intense laser pulses collide with GeV-class laser-wakefield-accelerated electron beams. We present simulation results showing that experiments will allow the study of laser-particle interactions from the classical to the SF-QED regime with a nonlinear quantum parameter of up to $\chi\sim$2. In addition, we show that experiments will enable the study and production of GeV-class, mrad-divergence positron beams via the Breit-Wheeler process.

Charged particle motion and radiation in strong electromagnetic fields

The dynamics of charged particles in electromagnetic fields is an essential component of understanding the most extreme environments in our Universe. In electromagnetic fields of sufficient magnitude, radiation emission dominates the particle motion and effects of nonlinear quantum electrodynamics (QED) are crucial, which triggers electron-positron pair cascades and counterintuitive particle-trapping phenomena. As a result of recent progress in laser technology, high-power lasers provide a platform to create and probe such fields in the laboratory. With new large-scale laser facilities on the horizon and the prospect of investigating these hitherto unexplored regimes, we review the basic physical processes of radiation reaction and QED in strong fields, how they are treated theoretically and in simulation, the new collective dynamics they unlock, recent experimental progress and plans, as well as possible applications for high-flux particle and radiation sources.

Dense Polarized Positrons from Laser-Irradiated Foil Targets in the QED Regime.

Many works have shown that dense positrons can be effectively generated from laser-solid interactions in the strong-field quantum electrodynamics (QED) regime. Whether these positrons are polarized has not yet been reported, limiting their potential applications. Here, by polarized QED particle-in-cell simulations including electron-positron spin and photon polarization effects, we investigate a typical laser-solid setup that an ultraintense linearly polarized laser irradiates a foil target with micrometer-scale-length preplasmas. We find that once the positron yield becomes appreciable with the laser intensity exceeding 10^{24} W/cm^{2}, the positrons are obviously polarized. Around 30nC positrons can acquire >30% polarization degree with a flux of 10^{12} sr^{-1}. The angle-dependent polarization is attributed to the asymmetrical laser fields that positrons undergo near the skin layer of overdense plasmas, where radiative spin flip and radiation reaction play significant roles. The polarization mechanism is robust and could generally appear in future 100-PW-class laser-solid experiments.

Electron response to radiation under linear acceleration: Classical, QED, and accelerated frame predictions

A model detector undergoing constant, infinite-duration acceleration converges to an equilibrium state described by the Hawking-Unruh temperature ${T}_{a}=(a/2\ensuremath{\pi})(\ensuremath{\hbar}/c)$. To relate this prediction to experimental observables, a pointlike charged particle, such as an electron, is considered in place of the model detector. Instead of the detector's internal degree of freedom, the electron's low-momentum fluctuations in the plane transverse to the acceleration provide a degree of freedom and observables which are compatible with the symmetry and thermalize by interaction with the radiation field. General arguments in the accelerated frame suggest thermalization and a fluctuation-dissipation relation but leave underdetermined the magnitude of either the fluctuation or the dissipation. Lab frame analysis reproduces the radiation losses, described by the classical Lorentz-Abraham-Dirac equation, and reveals a classical stochastic force. We derive the fluctuation-dissipation relation between the radiation losses and stochastic force as well as equipartition $⟨{p}_{\ensuremath{\perp}}^{2}⟩=2m{T}_{a}$ from classical electrodynamics alone. The derivation uses only straightforward statistical definitions to obtain the dissipation and fluctuation dynamics. Since high accelerations are necessary for these dynamics to become important, we compare classical results for the relaxation and diffusion times to strong-field quantum electrodynamics results. We find that experimental realization will require development of more precise observables. Even wakefield accelerators, which offer the largest linear accelerations available in the lab, will require improvement over current technology as well as high statistics to distinguish an effect.

Fitting XMM-Newton observations of the AXP 1RXS J170849.0−400910 with four magnetar surface emission models, and predictions for X-ray polarization observations with IXPE

Context. Phase-resolved spectral and spectropolarimetric X-ray observations of magnetars present us with the opportunity to test models of the origin of the X-ray emission from these objects, and to constrain the properties of the neutron star surface and atmosphere. Aims. Our first aim is to use archival XMM-Newton observations of the magnetar 1RXS J170849.0−400910 to ascertain how well four emission models describe the phase-resolved XMM-Newton energy spectra. Our second aim is to evaluate the scientific potential of future spectropolarimetric observations of 1RXS J170849.0−400910 with the Imaging X-ray Polarimetry Explorer (IXPE) scheduled for launch in late 2021. The most salient questions are whether IXPE is able to distinguish between the different emission models, and whether IXPE can unambiguously detect the signatures of quantum electrodynamics (QED) effects in strong magnetic fields. Methods. We used numerical radiation transport calculations for a large number of different system parameters to predict the X-ray flux and polarization energy spectra of the source 1RXS J170849.0−400910. Based on the numerical results, we developed a new model to fit phase-resolved and phase-averaged X-ray spectral (i.e., XMM-Newton and IXPE) and spectropolarimetric (IXPE) data. In order to test the sensitivity of IXPE to strong-field QED effects, we fit a simulated IXPE observation with two versions of the model, i.e., with and without QED effects accounted for. Results. The fixed-ions condensed surface model gives the best description of the phase-resolved XMM-Newton spectra, followed by the blackbody and free-ions condensed surface models. The magnetized atmosphere model gives a poor description of the data and seems to be largely excluded. Simulations show that the IXPE observations of sources such as 1RXS J170849.0−400910 will allow us to cleanly distinguish between high-polarization (blackbody, magnetized atmosphere) and low-polarization (condensed surface) models. If the blackbody or magnetized atmosphere models apply, IXPE can easily prove QED effects based on ∼200 ksec observations as studied here; longer IXPE observation times will be needed for a clear detection in the case of the condensed surface models. Conclusions. The XMM-Newton data have such a good signal-to-noise ratio that they reveal some limitations of the theoretical models. Notwithstanding this caveat, the fits clearly favor the fixed-ions condensed surface and blackbody models over the free-ions condensed surface and magnetized atmosphere models. The IXPE polarization information will greatly help us to figure out how to improve the models. The first detection of strong-field QED effects in the signal from astrophysical sources seems possible if an adequate amount of time is dedicated to the observations.

PICSAR-QED: a Monte Carlo module to simulate strong-field quantum electrodynamics in particle-in-cell codes for exascale architectures

Physical scenarios where the electromagnetic fields are so strong that quantum electrodynamics (QED) plays a substantial role are one of the frontiers of contemporary plasma physics research. Investigating those scenarios requires state-of-the-art particle-in-cell (PIC) codes able to run on top high-performance computing (HPC) machines and, at the same time, able to simulate strong-field QED processes. This work presents the PICSAR-QED library, an open-source, portable implementation of a Monte Carlo module designed to provide modern PIC codes with the capability to simulate such processes, and optimized for HPC. Detailed tests and benchmarks are carried out to validate the physical models in PICSAR-QED, to study how numerical parameters affect such models, and to demonstrate its capability to run on different architectures (CPUs and GPUs). Its integration with WarpX, a state-of-the-art PIC code designed to deliver scalable performance on upcoming exascale supercomputers, is also discussed and validated against results from the existing literature.

Single particle detection system for strong-field QED experiments

Measuring signatures of strong-field quantum electrodynamics (SF-QED) processes in an intense laser field is an experimental challenge: it requires detectors to be highly sensitive to single electrons and positrons in the presence of the typically very strong x-ray and γ-photon background levels. In this paper, we describe a particle detector capable of diagnosing single leptons from SF-QED interactions and discuss the background level simulations for the upcoming Experiment-320 at FACET-II (SLAC National Accelerator Laboratory). The single particle detection system described here combines pixelated scintillation LYSO screens and a Cherenkov calorimeter. We detail the performance of the system using simulations and a calibration of the Cherenkov detector at the ELBE accelerator. Single 3 GeV leptons are expected to produce approximately 537 detectable photons in a single calorimeter channel. This signal is compared to Monte-Carlo simulations of the experiment. A signal-to-noise ratio of 18 in a single Cherenkov calorimeter detector is expected and a spectral resolution of 2% is achieved using the pixelated LYSO screens.

Relativistic flying laser focus by a laser-produced parabolic plasma mirror

The question of electromagnetic field intensification towards the values typical for strong field Quantum Electrodynamics is of fundamental importance. One of the most promising intensification schemes is based on the relativistic-flying mirror concept, which shows that the electromagnetic radiation reflected by the mirror will be frequency up-shifted by a factor of 4 gamma^2 (gamma: the Lorentz factor of the mirror). In laser-plasma interactions, such a mirror travels with relativistic velocities and typically has a parabolic form, which is advantageous for light intensification. Thus, a relativistic-flying parabolic mirror reflects the counter-propagating radiation in a form of focused and flying electromagnetic wave with a high frequency. The relativistic-flying motion of the laser focus makes the electric and magnetic field distributions of the focus complicated, and the mathematical expressions describing the field distributions of the focus is important. We present analytical expressions describing the field distribution formed by an ideal flying mirror having a perfect reflectance over the entire surface and wavelength range. The peak field strength of an incident laser pulse with a center wavelength of lambda_0 and an effective beam radius of w_e is enhanced by a factor proportional to gamma^3 (w_e/lambda_0) in the relativistic limit. Electron-positron pair production is investigated in the context of invariant fields based on the enhanced electromagnetic field. The pair production rate under the relativistic-flying laser focus is modified by the Lorentz gamma-factor and the beam radius-wavelength ratio (w_e/lambda_0). We show that the electron-positron pairs can be created by colliding two counter-propagating relativistic-flying laser focuses in vacuum, each of which is formed when a 180 TW laser pulse is reflected by a relativistic-flying parabolic mirror with a gamma = 12.2.

Compressing magnetic field into a high-intensity electromagnetic field with a relativistic flying mirror

Ultrahigh electromagnetic fields (≥~1023 W cm−2) are necessary for the study of strong-field quantum electrodynamics (QED). In this study, for the first time, we propose the compression of a pre-seeding static magnetic field with a relativistic flying mirror to generate a high electromagnetic field. The produced field intensity can be further amplified to be 5 × 1023 W cm−2 owing to the multiple reflections between the flying mirror and a stationary solid target; this produced field intensity is approximately four orders of magnitude larger than that of the seeding field and far exceeds that of the driver laser field (9.6 × 1022 W cm−2). Therefore, the ultrahigh electromagnetic field can significantly facilitate strong-field QED effects such as high-energy gamma photon emission. An analytical theory is developed to self-consistently describe the motion of the flying mirror and the field amplification. The predications from the theory are well demonstrated by numerical simulations. The scheme of producing high-intensity electromagnetic fields proposed in this letter provides a new, powerful means to study strong-field QED with a relatively low laser intensity.

Vacuum instability in time-dependent electric fields: New example of an exactly solvable case

A new exactly solvable case in strong-field quantum electrodynamics with a time-dependent external electric field is presented. The corresponding field is given by an analytic function, which is asymmetric (in contrast to Sauter-like electric field) with respect to the time instant, where it reaches its maximum value, that is why we call it the analytic asymmetric electric field. We managed to exactly solve the Dirac equation with such a field, which made it possible to calculate characteristics of the corresponding vacuum instability nonperturbatively. We construct the so-called in- and out-solutions and with their help calculate mean differential and total numbers of created charged particles, probability of the vacuum to remain a vacuum, vacuum mean values of current density and energy-momentum tensor of the particles. We study the vacuum instability in regimes of rapidly and slowly changing analytic asymmetric electric field, and compare the obtained results with corresponding ones obtained earlier for the case of the symmetric Sauter-like electric field. We also compare exact results in the regime of slowly changing field with corresponding results obtained within the slowly varying field approximation recently proposed by two of the authors, thus demonstrating the effectiveness of such an approximation.

Probing Strong-Field QED with Doppler-Boosted Petawatt-Class Lasers

We propose a scheme to explore regimes of strong-field quantum electrodynamics (SF QED) otherwise unattainable with the currently available laser technology. The scheme relies on relativistic plasma mirrors curved by radiation pressure to boost the intensity of petawatt-class laser pulses by Doppler effect and focus them to extreme field intensities. We show that very clear SF QED signatures could be observed by placing a secondary target where the boosted beam is focused.

Pulse shape effects in high-field Bethe–Heitler pair production

Electron–positron pair production via the nonlinear Bethe–Heitler effect in the combined fields of a bare nucleus and a high-intensity laser pulse is studied theoretically. The calculations are performed within the framework of strong-field quantum electrodynamics using a flat-top laser profile with raising and falling edges. This way, the dependence of the pair production process on the precise shape of the laser field is analyzed. Our approach allows us, in particular, to follow the evolution of the created particles’ energy spectra from ultra-short few-cycle pulses to the monochromatic infinite pulse-train limit. We show how the various portions of the pulse influence these spectra and determine conditions for which the outcome from a laser pulse closely resembles the predictions from monochromatic theory.

Collisional strong-field QED kinetic equations from first principles

Starting from nonequilibrium quantum field theory on a closed time path, we derive kinetic equations for the strong-field regime of quantum electrodynamics (QED) using a systematic expansion in the gauge coupling $e$. The strong field regime is characterized by a large photon field of order $\mathcal{O}(1/e)$, which is relevant for the description of, e.g., intense laser fields, the initial stages of off-central heavy ion collisions, and condensed matter systems with net fermion number. The strong field enters the dynamical equations via both quantum Vlasov and collision terms, which we derive to order $\mathcal{O}(e^2)$. The kinetic equations feature generalized scattering amplitudes that have their own equation of motion in terms of the fermion spectral function. The description includes single photon emission, electron-positron pair photoproduction, vacuum (Schwinger) pair production, their inverse processes, medium effects and contributions from the field, which are not restricted to the so-called locally-constant crossed field approximation. This extends known kinetic equations commonly used in strong-field QED of intense laser fields. In particular, we derive an expression for the asymptotic fermion pair number that includes leading-order collisions and remains valid for strongly inhomogeneous fields. For the purpose of analytically highlighting limiting cases, we also consider plane-wave fields for which it is shown how to recover Furry-picture scattering amplitudes by further assuming negligible occupations. Known on-shell descriptions are recovered in the case of simply peaked ultrarelativistic fermion occupations. Collisional strong-field equations are necessary to describe the dynamics to thermal equilibrium starting from strong-field initial conditions.

Multistep pulse compressor for 10s to 100s PW lasers.

High-energy tens (10s) to hundreds (100s) petawatt (PW) lasers are key tools for exploring frontier fundamental researches such as strong-field quantum electrodynamics (QED), and the generation of positron-electron pair from vacuum. Recently, pulse compressor became the main obstacle on achieving higher peak power due to the limitation of damage threshold and size of diffraction gratings. Here, we propose a feasible multistep pulse compressor (MPC) to increase the maximum bearable input and output pulse energies through modifying their spatiotemporal properties. Typically, the new MPC including a prism pair for pre-compression, a four-grating compressor (FGC) for main compression, and a spatiotemporal focusing based self-compressor for post-compression. The prism pair can induce spatial dispersion to smooth and enlarge the laser beam, which increase the maximum input and output pulse energies. As a result, as high as 100 PW laser with single beam or more than 150 PW through combining two beams can be obtained by using MPC and current available optics. This new optical design will simplify the compressor, improve the stability, and save expensive gratings/optics simultaneously. Theoretically, the output pulse energy can be increased by about 4 times using the MPC method in comparison to a typical FGC. Together with the multi-beam tiled-aperture combining method, the proposed tiled-grating based tiled-aperture method, larger gratings, or negative chirp pulse based self-compression method, several 100s PW laser beam is expected to be obtained by using this MPC method in the future, which will further extend the ultra-intense laser physics research fields.

Realization of laser intensity over 1023 W/cm2

High-intensity lasers are critical for the exploration of strong field quantum electrodynamics. We report here a demonstration of laser intensity exceeding 1 0 23 W / c m 2 with the CoReLS petawatt (PW) laser. After wavefront correction and tight focusing with a two-stage adaptive optical system and an f/1.1 ( f = 300 m m ) off-axis parabolic mirror, we obtained near diffraction-limited focusing with a spot size of 1.1 µm (FWHM). From the measurement of 80 consecutive laser shots at 0.1 Hz, we achieved a peak intensity of ( 1.1 ± 0.2 ) × 1 0 23 W / c m 2 , verifying the applicability of the ultrahigh intensity PW laser for ultrahigh intensity laser–matter interactions. From the statistical analysis of the PW laser shots, we identified that the intensity fluctuation originated from air turbulence in the laser beam path and beam pointing. Our achievement could accelerate the study of strong field quantum electrodynamics by enabling exploration of nonlinear Compton scattering and Breit–Wheeler pair production.

Gauge invariant canonical symplectic algorithms for real-time lattice strong-field quantum electrodynamics

A class of high-order canonical symplectic structure-preserving geometric algorithms are developed for high-quality simulations of the quantized Dirac-Maxwell theory based strong-field quantum electrodynamics (SFQED) and relativistic quantum plasmas (RQP) phenomena. With minimal coupling, the Lagrangian density of an interacting bispinor-gauge fields theory is constructed in a conjugate real fields form. The canonical symplectic form and canonical equations of this field theory are obtained by the general Hamilton’s principle on cotangent bundle. Based on discrete exterior calculus, the gauge field components are discreted to form a cochain complex, and the bispinor components are naturally discreted on a staggered dual lattice as combinations of differential forms. With pull-back and push-forward gauge covariant derivatives, the discrete action is gauge invariant. A well-defined discrete canonical Poisson bracket generates a semi-discrete lattice canonical field theory (LCFT), which admits the canonical symplectic form, unitary property, gauge symmetry and discrete Poincaré subgroup, which are good approximations of the original continuous geometric structures. The Hamiltonian splitting method, Cayley transformation and symmetric composition technique are introduced to construct a class of high-order numerical schemes for the semi-discrete LCFT. These schemes involve two degenerate fermion flavors and are locally unconditional stable, which also preserve the geometric structures. Admitting Nielsen-Ninomiya theorem, the continuous chiral symmetry is partially broken on the lattice. As an extension, a pair of discrete chiral operators are introduced to reconstruct the lattice chirality. Equipped with statistically quantization-equivalent ensemble models of the Dirac vacuum and non-trivial plasma backgrounds, the schemes are expected to have excellent performance in secular simulations of relativistic quantum effects, where the numerical errors of conserved quantities are well bounded by very small values without coherent accumulation. The algorithms are verified in detail by numerical energy spectra. Real-time LCFT simulations are successfully implemented for the nonlinear Schwinger mechanism induced e-e+ pairs creation and vacuum Kerr effect, where the nonlinear and non-perturbative features captured by the solutions provide a complete strong-field physical picture in a very wide range, which open a new door toward high-quality simulations in SFQED and RQP fields.

Disentangling x-ray dichroism and birefringence via high-purity polarimetry

High-brilliance synchrotron radiation sources have opened new avenues for x-ray polarization analysis that go far beyond conventional polarimetry in the optical domain. With linear x-ray polarizers in a crossed setting, polarization extinction ratios down to 10 − 10 can be achieved. This renders the method sensitive to probe the tiniest optical anisotropies that would occur, for example, in strong-field quantum electrodynamics due to vacuum birefringence and dichroism. Here we show that high-purity polarimetry can be employed to reveal electronic anisotropies in condensed matter systems with utmost sensitivity and spectral resolution. Taking CuO and L a 2 C u O 4 as benchmark systems, we present a full characterization of the polarization changes across the Cu K-absorption edge and their separation into dichroic and birefringent contributions. At diffraction-limited synchrotron radiation sources and x-ray lasers, where polarization extinction ratios of 10 − 12 can be achieved, our method has the potential to assess birefringence and dichroism of the quantum vacuum in extreme electromagnetic fields.

Resummation of background-collinear corrections in strong field QED

Scattering processes in laser backgrounds are degenerate to the emission of unobservable photons collinear with the laser. We identify processes and observables for which such degeneracies factorize and exponentiate, obtaining the leading-order intensity dependence of these inclusive observables at high laser intensity, correct to all orders in the fine-structure constant (all loops, all emissions). The results show an exponential intensity dependence distinct from that predicted by the Narozhny-Ritus conjecture on the high-intensity behavior of quantum electrodynamics in strong fields.

Laser-driven radiation-reaction effect and polarized particle acceleration

<sec>Laser-plasma interaction at intensities beyond 10<sup>22</sup> W/cm<sup>2</sup> enters a new regime where gamma-photon emission and the induced radiation-reaction effect dominate. In extreme laser fields, high energy electrons emit gamma-photons efficiently, which take considerable portion of energy away and impose strong reaction forces on radiating electrons. When the radiation power is comparable to the electron energy gained in a certain period of time, the radiation-reaction (RR) effect becomes significant, which fundamentally changes the picture of laser-plasma interaction. In this review article, we introduce the physics of radiation-reaction force, including both classical description and quantum description. The effects of stochastic emission and particle spins in the quantum-electrodynamics (QED) RR process are discussed. We summarize the RR-induced phenomena in laser-plasma interaction and some proposed measurements of RR. As a supplement, we also introduce the latest progress of producing spin polarized particles based on laser-plasma accelerations, which provides polarized beam sources for verifying the QED-RR effects.</sec><sec>In the classical picture, the RR force can be described by the Landau-Lifshitz (LL) equation, which eliminates the non-physical run-away solution from the Lorentz-Abraham-Dirac (LAD) equation. The damping force could induce the electron trajectories to instantaneously reverse, electrons to cool and even high energy electrons to be reflected by laser pulses. The latter leads to a “potential barrier” at a certain threshold that prevents the electrons of arbitrarily high energy from penetrating the laser field. In general, classical LL equation overestimates the RR effect, thus calling for more accurate quantum description.</sec><sec>When the emitted photon energy is close to the electron energy, radiation becomes discrete. Quantum effects arise such that the process, also known as nonlinear multi-photon Compton Scattering, must be considered in the strong-field QED picture. This is resolved in the Furry picture by using the laser-dressed Volkov state in the local constant cross-field approximation (LCFA). The QED model is applied to particle dynamics via Monte-Carlo (MC) sampling. We introduce the prominent feature of quantum RR-stochastic photon emission. It allows the processes forbidden in classical picture to emerge, such as quantum ‘quenching’, quantum ‘reflection’, etc. These observables validate the strong-field QED theory. Recently, there has been a rising interest in identifying the spin effect in the QED-RR force. We summarize the latest progress of this topic, showing that when spins are coupled with photon emission the electrons of different spin states undergo distinctive RR force. </sec><sec>The RR force has a significant effect on laser-plasma interaction. The review paper introduces recent QED-MC based PIC simulation results. Some key features include electron cooling in laser-driven radiation pressure acceleration and the radiation-reaction trapping (RRT) mechanism. In the RRT regime the laser pulse conveys over 10% of its energy to gamma-photons, facilitating the creation of a highly efficient gamma-ray source and electron-positron pair. In addition, the paper mentions the major efforts to measure the RR effect in recent years. It relies on high energy electrons either colliding with ultra-intense laser pulses or traversing crystals. Primitive observations indicate that existing theories do not match experimental results. Further investigation is required in both SF-QED theory and experiment.</sec><sec>Finally, the review paper discusses the idea of laser-driven polarized particle acceleration as a supplement. The all-optical approach integrates pre-polarized gas target into laser wakefield acceleration, offering a compact all-optical polarized particle source, which is highly favorable for strong-field QED studies, high-energy colliders and material science.</sec>

ML-Based Analysis of Particle Distributions in High-Intensity Laser Experiments: Role of Binning Strategy.

When entering the phase of big data processing and statistical inferences in experimental physics, the efficient use of machine learning methods may require optimal data preprocessing methods and, in particular, optimal balance between details and noise. In experimental studies of strong-field quantum electrodynamics with intense lasers, this balance concerns data binning for the observed distributions of particles and photons. Here we analyze the aspect of binning with respect to different machine learning methods (Support Vector Machine (SVM), Gradient Boosting Trees (GBT), Fully-Connected Neural Network (FCNN), Convolutional Neural Network (CNN)) using numerical simulations that mimic expected properties of upcoming experiments. We see that binning can crucially affect the performance of SVM and GBT, and, to a less extent, FCNN and CNN. This can be interpreted as the latter methods being able to effectively learn the optimal binning, discarding unnecessary information. Nevertheless, given limited training sets, the results indicate that the efficiency can be increased by optimizing the binning scale along with other hyperparameters. We present specific measurements of accuracy that can be useful for planning of experiments in the specified research area.