Energy absorption in the laser-QED regime.

A theoretical and numerical investigation of non-ponderomotive absorption at laser intensities relevant to quantum electrodynamics is presented. It is predicted that there is a regime change in the dependence of fast electron energy on incident laser energy that coincides with the onset of pair production via the Breit-Wheeler process. This prediction is numerically verified via an extensive campaign of QED-inclusive particle-in-cell simulations. The dramatic nature of the power law shift leads to the conclusion that this process is a candidate for an unambiguous signature that future experiments on multi-petawatt laser facilities have truly entered the QED regime.

Introduction

With the impending completion and commissioning of the Extreme Light Infrastructure - in particular ELI Beamlines[1], and the Apollon laser facility[2] amongst other multi-petawatt laser systems, it will soon be possible to investigate an entirely new regime of plasma physics in the laboratory. At this frontier of high-energy density physics it will be possible to conduct experimental investigations into a range of topics including non-linear quantum electrodynamic (QED) processes that spawn electron-positron pair production[3] and laser wakefield acceleration of electron bunches up to multi-GeV energies[4]. There will also be the opportunity to make advances in several other fields within laser-plasma interactions such as coherent harmonic generation and focusing[5-7], attosecond science[8] ion beam characterisation and acceleration[9-13], electron beam generation via laser-channeling and hole-boring[14-16], and laboratory astrophysics[17].

In order for future experiments to be able to claim with categorical certainty that results are due to QED effects, it is necessary that there be an easily verifiable signature that the laser-plasma interactions have entered the QED regime. While there has been significant work exploring electron-positron cascade production[3,18-21] and laser energy absorption[22-25] far into the QED regime, in this paper, we consider the modifications to laser energy absorption on the cusp of entering the QED regime. In this paper, we suggest that such a signature could be found by considering the changes to absorption processes at extremely high laser intensities. It is already known that the processes by which plasma absorbs incident laser energy moves through several regimes as the incident laser intensity is increased[26]. These regimes are often identified by how the “hot” electron energy scales with the electron density, n, and/or the normalised amplitude of the laser, a0:where A is the vector potential, I the intensity, and λ the wavelength of the incident laser pulse, c is the speed of light in vacuum, and e and me are the electronic charge and rest mass. Here, as in previous work[27], and throughout the rest of this manuscript, the vector potential is defined in the Lorenz gauge.

At relatively modest laser intensities (), the dominant absorption processes are inverse bremsstrahlung[28] and resonant absorption[29]. As the laser power is increased to relativistic intensities - identified by - a number of other absorption mechanisms have been proposed theoretically and verified experimentally. Notable, and dominant, examples include the explicitly ponderomotive mechanism[30] and Brunel (or vacuum) heating[31,32]. Theoretical and numerical work in recent years has predicted that when moving into the regime of “ultra-relativistic” lasers incident on relativistically over-dense targets, absorption processes are dominated by non-ponderomotive mechanisms such as the zero-vector-potential (ZVP) absorption mechanism[27,33]. This regime is identified by , and where nc is the critical density of the plasma, above which the plasma is opaque to light of wavelengths less than or equal to λ0:

Here, is the permittivity of free space.

In the ZVP mechanism, the plasma frequency of the electrons is so high (due to the high density) that the electrons are able to respond adiabatically to the ponderomotive force, thus setting up a potential well that balances the radiation pressure of the laser with the electrostatic pressure of a plasma with spatially separated electron and ion fluids[33]. Thus any energy absorption can be considered to be “non-ponderomotive” in nature and occurs over multi-attosecond timescales. In the case of the ZVP mechanism, the radiation pressure from the laser field displaces the electron fluid from the ion background, setting up a pseudo-capacitor system. When the vector potential of the laser field passes through zero, the radiation pressure instantaneously vanishes and the electric field of the pseudo-capacitor causes the electrons to shuttle across the pseudo-capacitor, and be accelerated to energies on the order of several keV up to MeV[27]. These high-momentum bunches of electrons, such as those shown in Fig. 1 are seen to co-propagate with the zeroes in the vector potential, hence the nomenclature ZVP absorption[27,33].

Figure 1

A momentum-position phase space plot in the laser-forward direction demonstrating the high-momentum peaks characteristic of the ZVP mechanism. The laser amplitude, was set to and the electron density was 50 times critical.

The non-relativistic, ponderomotive, and non-ponderomotive regimes are distinguished experimentally by the scaling of the “hot electron” energy with a0. Specifically, , where: for the non-relativistic case[29], in the ponderomotive regime[26], and in the non-ponderomotive regime[27,33]. In this paper, we propose an alteration to the ZVP model devised by Baeva et al.[33] that accounts for the possibility of generating electron-positron pairs and model how such pair production would affect the scaling of “hot” electron energy with a0. By considering this energy scaling, it is suggested that laser-energy absorption by the plasma electrons offers a possible experimentally verifiable signature of entrance into the QED regime.

The structure of the paper is as follows: firstly, a theoretical justification for the ZVP’s continued effect at QED-relevant intensities is considered. Next, a simple physical model capturing the dynamics of energy absorption is presented, including the effect of electron-positron pair production. Finally, the results of a series of particle-in-cell simulations are then presented which numerically validate the predictions of our model and pictorially demonstrate the signature we expect will indicate the entrance of future experimental campaigns into the QED regime.

Results

ZVP at QED-relevant intensities

The existence of the ZVP mechanism relies on a shift in perception when changing between reference frames. If one considers a laser pulse incident on a plasma with an ablating front surface, there are two principal frames of interest. The first is the rest frame of the bulk plasma, or “lab frame”; and the second is the rest frame of the expanding ablation front which will be referred to here as the “rest frame”. In the rest frame, the incident laser radiation changes form at the critical density surface from a propagating oscillation to an evanescent decay[33]. Previous work has established that for non-relativistic ablation velocities, for any angle of incidence, transforming into the lab frame yields a short skin depth, δ, over which the oscillations in the laser pulse’s vector potential propagate beyond the critical density surface[27]. As δ is independent of the decay length-scale of the ablation front, λS, it is possible to enter the ZVP regime provided that λS is sufficiently small.

As laser intensities are increased, the velocity of the ablation front will also rise[34,35]. As such, it is prudent to verify that the ZVP mechanism is still a viable absorption process at relativistic velocities. If one considers an obliquely incident laser pulse, its vector potential beyond the critical density surface in the rest frame, will be an evanescent wave because propagating waves cannot be sustained beyond the critical density surface. As only the temporal component of the vector potential will retain its oscillatory nature can be described as:where A0 is the amplitude and the frequency of the laser pulse, t is the time of propagation, r′ is the distance of propagation into the plasma, is the unit vector in the direction of polarisation of the laser pulse, and primed variables indicate that the quantity is defined and measured in the rest frame (consequentially, non-primed variables are defined and measured in the lab frame).

By constraining the problem to a p-polarised pulse propagating in the x–y plane, with being the axis of plasma ablation, the unit vector , in the lab frame can be written as:where θ is the angle of incidence of the laser pulse onto the expanding plasma. By considering the headlight effect and the Lorentz transformation of the electromagnetic 4-potential[36] (cf. Methods Section), a series of deductions can be made about the vector potential in the lab frame:

By applying a Lorentz transformation to the position 4-vector and converting variables to a pseudo-equivalent in the lab frame, the final form for A can be deduced:where , , , , and all coordinate variables are defined in their usual way as outlined in the Methods section.

It is clear from Equation 8 that zeroes in the vector potential can continue to propagate into the over-dense plasma, even for relativistic ablation velocities, a simple schematic of this is shown in Fig. 2. With this principle established, it is possible now to consider how QED processes will impact the electron energy scaling with laser intensity.

Figure 2

A simple schematic of how transforming from the rest frame of the expanding plasma into the laboratory frame leads to a continued propagation of the vector potential wave (red) beyond the plasma’s critical surface - indicated by the boundary between dark (vacuum) and white (plasma).

New energy scaling

In the original model for the ZVP absorption mechanism, the fast electron energy was deduced by assuming that the laser pressure displaces the electron fluid from the ionic background to set up a pseudo-capacitor system with the plates replaced by two charge fluids of charge Q, separated by a small displacement, Δr[27,33]. This model yielded a dependence of the fast electron energy on the laser amplitude of , which was successfully numerically validated up to a value of [27].

Beyond , QED effects – such as radiation reaction and pair production – have the potential to contribute to laser-plasma interactions. As the intensity of the laser is further increased, the extent of pair production is predicted to increase logarithmically[3] and radiation reaction is expected to inhibit the onset of relativistic transparency, thus maintaining the opaqueness of the plasma[23]. Ultra-high intensity laser pulses are expected to act as a background electromagnetic field that can stimulate the Breit-Wheeler pair production process[37]. When applied to the preexisting ZVP model, this leads to the assumption that at the point of the laser’s peak intensity, the rate of pair production will be maximal. This position of peak intensity coincides with the position to which the electron fluid is displaced by the laser pulse, a representation of which can be seen in Fig. 3.

Figure 3

A pictorial representation of the proposed modification to the ZVP mechanism to include QED effects. The original pseudo-capacitor, characterised by a displaced electron fluid (blue) and a consequent region of net positive space charge (red) is augmented by pair production at the peak intensity of the incident laser pulse (faded purple) on the negative “plate”. The positrons (yellow), not feeling a restoring force, propagate into the bulk plasma and leave the region of interest. In contrast, some of the pair-produced electrons (green), are subject to a restoring force which keeps a fraction of them within the region of interest, thus increasing the net charge in the region of the original negative “plate”.

With electron-positron pairs being produced on one of the pseudo-capacitor’s “plates”, it is necessary to modify the ZVP model to account for these additional particles. As the e−-e+ pairs are produced at the negative plate, one would expect that the pairs will be subject to the electric field of the pseudo-capacitor. Thus, as the pairs will be created with high energy and momentum predominantly in the laser-forward direction, the positrons will propagate unhindered into the quasi-neutral bulk plasma.

In contrast, some number, f, of the newly created electrons could be captured onto the negative “plate” by the pseudo-capacitor’s electric field, as shown in Fig. 3. This problem is best considered by comparing the kinetic energy, Tesc, required to escape from the pseudo-capacitor’s electric field to the average energy, , of the newly created electronic energy distribution. The escape energy will be of the same order of the energy that an electron would gain by crossing the capacitor, which was found – in prior work[27,33] – to be proportional to , while numerical investigations into the Bethe-Heitler process[38] have found that the average energy of pair-produced positrons (and by inference, pair-produced electrons) is proportional to 1/a0[39].

If one considers the pair-produced electronic energy distribution function , where is the ratio of energy to , then the number of captured electrons can be thought of as:

In the limit that the escape energy is much less than the average pair-produced electron energy it is possible, in principle, to Taylor expand the distribution function and by preserving only first-order terms, Equation 9 becomes:

The final consideration, therefore, is to examine how the additional electrons, f, on the negative “plate” of the pseudo-capacitor impacts the energy scaling of the fast electrons. If the pseudo-capacitor has a cross-sectional area, σ, then the associated electric field in the non-QED case is given by . However, at QED-relevant intensities, the capacitor adds some number of electrons, f, to the negative “plate” such that while the charge of the positive “plate” remains as Q, the negative “plate” is disrupted and the region acquires additional negative charge such that the total charge of the highly localised negative region can be expressed as:

By applying Gauss’ law to this system of two regions of charge, Q and , it can be deduced that the electrons comprising the region are subject to an electric field given by:

The kinetic energy, T, of the fast electrons is calculated by considering the energy gained as an electron crosses the quasi-neutral region, under the influence of the electric field given in Equation 12, from the area of uncompensated negative charge to the uncompensated ion charge:

By using the scaling relation derived in Equation 10 this yields:provided that . With a cubic increase in the dependence of T on a0 it is expected that on entering the QED regime, the change in electron energy scaling with incident laser intensity will become swiftly apparent.

While the scaling relation for the energy of the most energetic electrons is predicted to move from Equations 15 to 16 upon entry into the QED regime, this does not necessarily mean that the scaling of laser-energy absorption will follow the same relation. Indeed, it has been shown previously[33], that while in the non-ponderomotive regime the fastest electrons’ energies scale with , it has been predicted that the total amount of laser-energy absorbed will scale with . This scaling arises logically by considering that the total amount of energy absorbed will be proportional to the number of electrons participating in the absorption mechanism. In the case of the non-QED regime, this number can be approximated to be the number of electrons required to set-up the pseudo-capacitor which has previously been shown to be proportional to a0[33]. The case for total absorption after entering the QED regime becomes even more difficult to predict, as the energy could be partitioned amongst the electrons, ions, re-radiated as coherent x-rays[40], distributed in surface waves as plasmons[41], or in pair-production processes. For this reason, it is expected that the scaling of the highest energy in the electron energy spectrum with laser intensity will provide the clearest experimental signature of entry into the QED regime.

Numerical validation

In lieu of experimental measurements to verify the above model, fully relativistic and kinetic particle-in-cell codes were used to run QED-inclusive simulations of intense laser pulses interacting with critically over-dense plasma.

The particle-in-cell codes OSIRIS[42] and EPOCH[43] were both used to run a parameter scan of the fast electron energy as the value of a0 was increased from 100 up to 560. The simulations were run in two dimensions, as the dynamics of the ZVP interaction are known to be confined to the plane of propagation of the laser pulse[27]. The simulations modelled an intense p-polarised laser pulse of varying a0, with a wavelength of 1 μm and a Gaussian temporal profile with a pulse duration of 85 fs, incident upon an aluminum plasma 1 μm thick with an exponentially decaying density profile of decay-scale length 0.2 μm on the laser-incident side. The density of the plasma was set to be 5.55 × 1022 cm−3, corresponding to 50 times the critical density, this is a density typical of templated foams - a type of target finding increasing usefulness in laser-plasma experiments[44]. A package comprising a QED Monte Carlo integrator[45] was included in the simulation to model pair production. The energy of the fast electrons, and the number of positrons at the end of the simulation were extracted from each simulation run to determine the trend of fast electron energy and positron number against incident laser intensity.

Figure 4 plots the fast electron energy against a0. It is clear that around , there is a distinct change in the dependence of T on a0. This corresponds to the detection of pair-produced positrons (also plotted in Fig. 4. By dividing the acquired data into two sets: the first corresponding to when pair production is not observed, and the second corresponding to when pair production is prevalent, it is possible to obtain qualitative information about the power law linking T and a0. In the non-QED regime:which is broadly in agreement with previous work investigating the ZVP mechanism[27,33]. In contrast, upon entering the QED regime, the dependence of the fast electron energy on the laser intensity changes dramatically to:

Figure 4

Plot of extracted fast electron energy against incident laser a0. Black triangles indicate EPOCH simulations for which no positrons were measured while the black crosses indicate equivalent OSIRIS simulations, the red line is a fit to both sets of data corresponding to . The blue asterisks indicate results from EPOCH simulations for which positrons were detected. The green line is a fit corresponding to . Filled-in black circles indicate the measured number density of positrons in the EPOCH simulations for which positrons were detected. There is a clear correlation between the onset of pair production and the shift in energy scaling. The region of interest for calculating the positron number density began at the point of laser focus (the boundary between the bulk plasma and skin layer), and extended 0.2 μm into the bulk plasma. Inset is a plot copmaring the fast electron energy in simulations where the EPOCH QED package was turned on (blue asterisks and green line), and off (black triangles and red line).

This is in excellent agreement with the prediction of Equation 14. It is noteworthy that if the QED calculations in the presented simulations are suppressed, the kinetic energy of the electrons does not deviate from the scaling of Equation 15, lending confidence to the hypothesis that entrance into the QED regime cause the scaling change seen in Equation 16.

Discussion

We have predicted and numerically verified that there is a distinct regime change for the electronic absorption of incident laser energy that is coincident with the onset of quantum electrodynamics effects. The prediction that the power law linking fast electron energy and incident laser amplitude would change from quadratic to quintic was borne out by the results of an extensive simulation campaign. Throughout the simulation campaign, absorption efficiencies ranging from 30% to 60% were observed. While it could be expected that the rapidity of a quintic scaling’s growth will cause the electron energy to exceed the laser energy beyond some threshold value of a0, this has not yet been observed in our simulations, thus leading us to reason that the model proposed in this paper does not violate energy conservation up to intensities as extreme as 5 × 1023 W cm−2. This leads us to the conclusion that QED effects such as pair production augment non-ponderomotive absorption mechanisms, such as the ZVP mechanism in a dramatic fashion.

Such a significant change in energy scaling, makes measurement of the ZVP mechanism a candidate for being a clear and unambiguous indicator that future high-power laser systems, such as the Extreme Light Infrastructure have truly breached the energy barrier into the QED regime. This clarity could allow experimental investigations into a wealth of unexplored physics to proceed with the confidence that quantum electrodynamics is truly playing a role in extreme high energy density physics.

Methods

Relativistic dynamics

The Lorentz transformation and headlight effect[36] were used when transforming variables between the laboratory and rest frames detailed in the main text:where β is the ratio of the ablation velocity to the speed of light in vacuum, is the well-known Lorentz factor, is the electromagnetic 4-potential, is the space-time 4-vector, and is the Lorentz transform given by:

Particle-in-cell simulations

Both OSIRIS[42] and EPOCH[43] operate according to the same principles of particle-in-cell (PIC) codes. Particles are modelled as “macro-particles”. The negatively charged macro-particles can be considered to be electrons as they bear the same charge:mass ratio as electrons. The positively charged macro-particles were modelled to have the same charge:mass ratio as deuterons. PIC codes update the electromagnetic fields’ spatial distributions and macro-particles’ positions and trajectories on alternate time-steps, i.e. on every odd time-step, the fields are updated by solving Maxwell’s equations, given the distribution of macro-particles. On each subsequent even time-step, the macro-particles are updated according to the Lorentz force law given the updated electromagnetic fields. QED packages include the addition of a Monte Carlo integrator[45] to calculate the probabilities and rates of pair-production, amongst other QED effects.

One key difference between EPOCH and OSIRIS is that while EPOCH takes Système International (SI) units, OSIRIS uses a dimensionless variation where all quantities are normalised to some user-defined reference time- and length-scales. For the benefits of discussion here, the quantities are listed in SI units for both codes, a simple numerical correction was computed for the OSIRIS parameters.

The total size of the simulation box was 20 μm × 20 μm. 100 cells per micron were used to resolve the -direction (corresponding to the axis of laser propagation), and 25 cells per micron were used to resolve the -direction. 20 macro-ions and 500 macro-electrons were modelled within each cell. The laser was polarised in the -direction. The time-step increment of the simulations was set to 2.5 as for a total duration of 200 fs. When extracting the fast electron energy plotted in Fig. 4, the energy of the most energetic 5% of electrons was considered. This ensured that no false positive result for the fast electron energy, that could potentially arise from numerical error on a single macro-particle, was recorded.

OSIRIS was used in prior work[27,33] and in this paper to obtain simulation data points from . EPOCH was used to obtain data from . Repeating the OSIRIS simulations using EPOCH verified that the unit conversion between the two simulations was correct and ensured that the results of this work were consistent with previous work[27,33]. EPOCH was used as the main code for the QED-relevant section of the results due to the established nature of its QED package[45], and the fact that EPOCH is an open-access code.