Repositório ISCTE-IUL

Scattering relativistic electrons with optical lasers can result in a signiﬁcant frequency upshift for the photons, potentially producing γ -rays. This is what linear Compton scattering taught us. Ultra-intense lasers oﬀer nowadays a new paradigm where multi-photon absorption eﬀects come into play. These eﬀects can result in higher harmonics, higher yields and also electron-positron pairs. This article intends to discriminate the diﬀerent laser scenarios that have been proposed over the past years as well as to give scaling laws for future experiments. The energy conversion from laser or particles to high-frequency photons is addressed for both the well-known counter propagating electron beam-laser interaction and for Quantum-electrodynamics cascades triggered by various lasers. Constructing bright and energetic gamma-ray sources in controlled conditions is within an ace of seeing the light of day.


I. INTRODUCTION
Ultra-intense laser-matter interactions can generate γrays. Several recent experiments have obtained hard Xrays in the laboratory using counter-propagating electron beams and laser pulses 1-5 . A review on laser-wakefield acceleration (LWFA) -based light sources summarizes other viable configurations for generating energetic photons and the properties of radiation that could be obtained with a reflection to prospective applications 6 .
Multiphoton Thompson scattering was observed in a recent experiment by Yan et al 7 . The next goal is to obtain a high conversion efficiency of the electron energy into high-frequency radiation. This requires operating the source in a regime of significant radiation recoil, in classical or quantum regime of interaction. First evidence of electron slowdown in a laser field 8,9 are consistent with the classical radiation reaction predictions for scattering a LWFA electron bunch and a laser pulse 10 . However, there are still many open questions regarding the transition from the classical to the quantum radiation reaction dominated regime, that has generated keen interest in the last few years [11][12][13][14][15][16][17][18] .
The radiation reaction dominated regime (i.e. the regime where particles lose a substantial amount of energy through radiation emission) can be reached either by using higher laser intensities or more energetic electrons. The particles in this regime either emit a few very energetic photons (quantum radiation reaction) or a large number of low-energy photons successively (classical radiation reaction) 19 . Current published energy record for LWFA electrons is 4 GeV and those electrons a) Electronic mail: marija.vranic@ist.utl.pt b) Electronic mail: luis.silva@ist.utl.pt are obtained using a 16 J laser 20 . In the next generation laser facilities, electron energies on the order of 10 GeV are expected. Pairing such electron beams with intense lasers can provide access to more extreme regimes of interaction. It may even be possible to obtain multiple electron-positron pairs from electron beam-laser collision. In the two-step approximation a particle first emits a high-energy photon, which then decays into an electronpositron pair via the Breit-Wheeler (BW) process in an intense electromagnetic background. A milestone experiment was performed at SLAC 21 , where BW pairs were produced in a collision of a 46.6 GeV electron beam from a linear accelerator with a green laser at the intensity of ∼ 10 18 W/cm 2 . The next generation of laser facilities is expected to deliver a much higher pair yield 22,23 even using electrons with lower energies (on the order of a few GeV). Further laser-electron scattering experiments are planned, both relying on LWFA and conventional accelerators to provide the electrons.
Another way to generate γ-rays and pairs with intense lasers are QED cascades 19,24 . One of the most popular configurations is to use two colliding lasers to create an intense standing wave 25-30 whose nonlinear evolution in the presence of self-generated plasma can be studied theoretically resorting to QED-PIC simulations 28,[31][32][33][34][35][36][37][38] . New laser facilities will access extreme regimes of interaction [39][40][41][42] where we can expect an abundance of electron-positron pairs. This motivated a lot of theoretical effort to improve our understanding of QED cascades, considering configurations using multiple lasers [43][44][45] , as well as the challenges concerning the cascade seeding 46,47 . Among the proposed solutions for the seeding problem, it was proposed to use solid targets 48,49 . Configurations with nanowires were also proposed, to enhance the laser heating of the target electrons 50 .
In this manuscript, we focus on the gamma-ray emission. We revisit the two main configurations: laser-electron beam scattering and the two-laser QED cascade. The first aim is to connect the hard-photon emission in both of these scenarios and bring an intuitive understanding of why even in a QED cascade, the classical absorption is more important than the instantaneous quantum absorption. The ratio of the classical and quantum absorption has been found previously using QED calculations 51 to be proportional to the square of the local normalized vector potential a 0 , where a 0 = 0.86 I[10 18 W/cm 2 ]λ[µm] for linearly polarized lasers and a 0 = 0.61 I[10 18 W/cm 2 ]λ[µm] for circular polarization. In other words, this ratio is proportional to the laser intensity, and for high laser intensities, the classical absorption dominates. Here, we show how the electronphoton scattering in a QED cascade can be mapped to a simple laser-electron scattering scenario for which the quantum vs. classical absorption ratio has already been calculated. The second aim is to estimate how much energy is converted to hard photons in both configurations and identify the relevant regimes of interaction. We also give a brief summary of the scaling laws for evolution of the electron energy distribution function in the laserelectron beam scattering. This is of importance for planning of experiments, because the electron beam properties imprint on the emitted radiation. The ideas and scalings presented in this manuscript are relevant for the multi-petawatt laser projects such as ELI 39 , Apollon 52 and CoReLS 53 that aim to reach unprescedented laser intensities, as well as FACET-II 54 and LUXE 55,56 that plan to perform laser scattering experiments using 10 GeVclass high-quality electron beams. The manuscript is organized as follows. In the following section, we discuss the ratio of classical vs. quantum absorption in a scattering of a single electron with one wave. We mostly discuss the counter-propagating geometry, because it provides the highest energy photons. We then extend the ideas presented for a single wave-electron interaction, to a lepton interacting with a standing wave formed by two colliding lasers. We then review the scaling laws and the energy conversion expected for the electron-laser scattering configuration. We finally discuss the radiative absorption in a two-laser cascade. We distinguish two regimes: a regime of controlled radiation emission, when the wave is not severely affected by the presence of the plasma and a regime where the plasma density is high enough to disrupt the wave (this can happen because the target is dense in the beginning, or due to the considerable production of electron-positron pairs). We identify the parameters where the pair production is low enough to operate a controlled γ-ray source in a low-absorption regime with solid hydrogen targets already available in the laboratory. The results are supported with QED-PIC simulations for a range of parameters, both in 2D and 3D geometry. Quantum vs. classical laser absorption by a single electron was first considered by Meuren et al. 51,57 using the QED formalism. Here, we provide an intiutive picture for the corresponding scaling laws.

A. Laser-electron scattering
The first studies dealing with interaction of electrons with intense laser beams date back to a few years after the invention of the laser 58,59 . It was found that a mass change induced in the electron (considered initially at rest) by the external field of the laser shifts the wavelength of the scattered photons by an amount depending on the intensity of the incident beam. Furthermore, the absorption of multiple laser photons becomes possible, which facilitates the emission of harmonics. If the electron is initially relativistic, the scattering with the laser photons can result in the emission of γ-rays.
Effectively, only a finite space-time region is relevant for the photon emission process. It is characterized by the so-called formation length l f that depends on the wave intensity 60 where λ is the wavelength of the wave. Equation (1) is valid for χ e 1, where the definition of quantum parameter for electrons counter-propagating with an optical laser (λ = 1 µm) can be approximated as χ e 5 × 10 −6 γ 0 a 0 . According to Eq. (1) the radiation is formed over multiple laser cycles if a 0 1. In this regime it is useful to introduce the so called "dressed momentum" for the emitting electron, which is obtained by averaging the classical canonical momentum over a laser period 58,60,61 . Correspondingly, the classical mass change mentioned above appears.
For a 0 1 we have l f λ and the radiation process happens almost instantaneously in an effectively constant background field. The concept of dressed momentum in this regime becomes meaningless 51,57,62 . Instead, a semiclassical approximation becomes possible, i.e. the radiating electron follows a classical trajectory between localized emission events 63,64 .
In this regime we can distinguish two types of energy transfers: classical and quantum absorption 51,57 . The term quantum laser absorption LQ is defined as the direct contribution of the laser photons that scatter with the electrons during the quasi-instantaneous emission process. Any laser absorption that occurs elsewhere (the laser energy invested in the electron acceleration) can be described classically, and is incorporated into the instantaneous electron energy contribution e . The analysis of the energy transfer from electrons to photons allows estimating the importance of quantum laser absorp- Laser-electron head-on configuration; c), d) Two-laser standing wave, when interacting particle instantaneous momentum is perpendicular to the laser axis; e),f) Two-laser standing wave interacting with a particle at an oblique angle.
tion during the electron interaction with an intense laser beam.
The scattering process is affected by the presence of the intense laser field in two ways. One aspect is that there is a temporary momentum transfer from the laser to the electrons: the laser induces electron oscillations. A relativistic electron counter-propagating towards an intense laser with normalized vector potential a 0 , acquires an additional transverse momentum on the order of ∼ a 0 mc during the interaction, while the effect on the longitudinal momentum is small if a 0 γ 0 . For circular polarization, the additional electron transverse momentum is constant in magnitude but periodically changes the direction revolving around the propagation axis, while for linear polarization, the direction is fixed, but the magnitude changes periodically. The second effect of the very intense field to the scattering process is that the local field energy density is high, which means there can be many photons within the interaction volume. Such photon density allows to have frequent repeated scatterings, as well as to absorb more than one photon at a time in a single scattering.
The following paragraph illustrates intuitively, by comparing the typical timescales of acceleration and emission, how the factor a 2 0 between the classical and quantum photon absorption might arise. Let us assume that a 0 1 and the particle emits a hard photon not more than once during one full laser cycle. The photon formation time is τ f = l f /c, where c is the speed of light. The laser period is T L . As T L /τ f ∼ a 0 , the total work of the laser electric field during one period on one particle is a 0 times larger than during one emission event. However, at relativistic intensities, an electron can take much more energy from the wave, on the order of ξ ∼ a 2 0 mc 2 (assuming the electron is initially at rest). This is because a relativistic electron can have a longer effective interaction time, and one oscillation can last longer than a simple laser period T L (for example, if an electron copropagating with the wave). From there we get that T eff ∼ a 0 T L and T eff /τ f ∼ a 2 0 . One should note that the value of T eff depends on the initial particle energy and the scattering angle, but it is never smaller than T L /2. If we have n emissions during one laser cycle, then the average time between two emissions is T eff /n. The relevant ratio then becomes T eff /(nτ f ) ∼ a 2 0 /n. If a particle is counter-propagating with a wave, then most of the energy of the emitted photons is invested by the electron. However, the laser has to invest a few photons (even if not many) into the electron acceleration (classically), and also during the actual emission event in order to facilitate it (quantum absorption). One can show that the classical vs. quantum absorption of the laser photons scales as a 2 0 . As the quantum absorption is negligible for a 0 1, it is not considered in the QED-PIC algorithms. However, the classical absorption that happens due to the laser interaction with the plasma particles (i.e. particle acceleration) is intrinsically included in the PIC algorithm.
B. Mapping photon emission in a QED cascade produced by two colliding lasers with the photon emission in laser-electron scattering The ideas presented for a laser -electron scattering do not trivially port to the standing wave configurations with multiple lasers. Particles can be initially at rest. The particles first get accelerated, then lose energy due to radiation emission, and then get re-accelerated 65 by the electric field in the standing wave. All the energy radiated to high-frequencies comes from the laser field: either through accelerating the electrons (classical absorption) or by providing photons for scattering to occur (quantum absorption). The question is: which channel is dominant?
We can make use of what we know about particle dynamics in a standing wave to establish a connection between this setup and the simple laser-electron scattering discussed in the previous section. Electrons (and positrons) gain momentum on the order of a 0 that is perpendicular to the laser propagation axis 34,47,65 . In the case of circular polarization, leptons keep rotating, always remaining perpendicular to the laser propagation. This allows for a simplified consideration, as the laserphoton scattering always occurs at the right angle.
If the particle scatters with photons from only one laser, the situation is exactly the same as a scattering with one laser at 90 degrees. However, in principle, we could have a linear combination of m photons from one wave and n photons for the other participating in the scattering (not every combination is necessarily allowed, but here we are assuming a most general case). As the particle Lorentz factor is typically on the order of a 0 , the energy of an individual photon in the particle rest frame is ξ ph ∼ a 0 ξ ph regardless which laser the photon belongs to. Furthermore, the photons of both waves are practically co-propagating in the electron rest frame, as the Lorentz boost gave them a momentum in the same direction. In other words, the electron sees them almost as one wave, apart from the tiny difference in the momentum (∼ 1 eV for optical photons) perpendicular to the Boost direction (see Fig. 1 d)).
Particle dynamics is slightly different when we have linearly polarized laser pulses forming the standing wave. The electric field is perpendicular to the laser propagation axis, but the magnetic field then rotates the momentum, and the particles can be found counter-propagating with one of the waves. In the electron rest frame, counterpropagating photons are upshifted, while the photons of the co-propagating laser are downshifted. The most energetic photons in the electron rest frame are the counterpropagating ones. We can then see the analogy with the laser-electron scattering emerge naturally, as one wave becomes more important then the other. We note here that the normalized vector potential a 0 is Lorentzinvariant, which means that the wave with higher energy individual photons in the electron instantaneous rest frame has a higher energy density. There is a range of possible angles of incidence, with two limiting cases: particles being perpendicular to the laser propagation axis or counter-propagating with one of the lasers (see Fig. 1 e)). In general, any of these configurations can be mapped to an electron-laser scattering We have not discussed a decay of a hard photon into a pair here, because the number of emission events in a QED cascade surpasses the number of pair production events by orders of magnitude. They contribute little to the overall energy balance. Nevertheless, it was shown in Ref. 51 that the ratio between the classical and quantum laser absorption for pair production is also CL / LQ ∼ a 2 0 .

III. SCALING LAWS FOR LASER-ELECTRON SCATTERING
In a collision between an intense laser and an electron beam, the final electron energy after the interaction can be estimated as 10 for k < 1. The coefficient k depends on the laser duration τ 0 at FWHM and peak intensity I 0 in the following way where I 22 = I 0 10 22 W/cm 2 and θ represents the angle of interaction (for counter-propagation the cos θ −1). The counter-propagating configuration, therefore, corresponds to the strongest radiation reaction or losing a largest fraction of the electron energy. The energy converted to photons is then approximately Note that the electron beam energy bandwidth is bound to rise due to stochasticity in the quantum regime χ e 1 11,12 . But, the high-energy electrons, on average, radiate more than low-energy electrons. This tends to reduce the energy bandwidth, even in the quantum regime of interaction. The electron distribution function either spreads, or shrinks, depending on the local conditions. In the limit where the scattering is still Thompson in the electron rest frame, one can derive an expression for an instantaneous "turning point" 16 . If the standard deviation of the electron energy distribution function σ is larger than σ T , then the electron distribution function shrinks. For σ < σ T , stochasticity dominates and the bandwidth of the electron energy distribution rises. The value of σ T is given by a correction for the electron Gaunt factor 17,18 . However, the final expression then becomes more complex, and our aim here is to keep the scaling laws as simple as possible. The simplicity allows to estimate an asymptotic energy spread 16 as a function of the initial electron energy γ 0 and the laser intensity and duration: Equation (6) is useful for planning experiments, because it allows for a quick estimate of the expected final width of the electron energy distribution function. One can also predict the final divergence of the beam. Let us define the divergence as the average deflection angle θ F from the main axis of beam propagation. In this case, one can estimate this value as Comparisons of Eq. (7) with QED-PIC simulations are given in Fig. 2. Blue dashed line represents the θ F obtained using values of γ F and σ F Eqs.
(2) and (6). Red line represents the values obtained with Eq. (7), but using σ F and γ F measured in the simulation. The electron beam initial energy was 0.85 GeV, and it interacted with a circularly polarized laser of a 0 = 27. All other simulation parameters are given in the Appendix.

IV. ENERGY CONVERSION TO HIGH-FREQUENCY LIGHT IN A TWO-LASER STANDING WAVE (LINEAR POLARIZATION)
The number of photons of a plane wave with the normalized vector potential a 0 in a volume of λ 3 0 can be estimated as Here a S = mc 2 / ω 0 = 4.12×10 5 λ 0 [µm] represents the Schwinger critical field in units normalized to the laser frequency. In other words, this defines the dimensionless normalized vector potential that corresponds to a Schwinger field at a given wave frequency. For our calculation, the total number of available photons in a volume of λ 3 0 given by Eq. (8) should be multiplied by a factor of 2, as the standing wave is formed by two counterpropagating traveling waves.
The temporal structure of the standing wave causes a periodic emission with a period of T L . The constructive interference for magnetic and electric field is temporally spaced at T L /4. Periods of particle acceleration (when E is large) and rotation (when B dominates) therefore alternate every T L /4. The characteristic cycle of emission is T L /2 with the second half of the laser period repeating exactly the same particle dynamics in the opposite direction of motion.
Particle trajectories are chaotic 66,67 and stochastic emission does not allow for a general analytical estimate of the radiated energy using the trajectories alone. However, average energy absorbed per particle during the half-period T L /2 can be approximated using ideal simulations where pair production and current deposition are suppressed (see Figs. 3 and 4). In such simulations, the wave is not disturbed by the presence of the plasma. Photons do not decay into pairs, but the particles do experience quantum radiation reaction due to the hard photon emission.
Equation (9), obtained as a best fit to the available data, can then be used to estimate the importance of radiative laser absorption for non-ideal scenarios. What we mean by radiative absorption is the energy that was transferred from the laser to the particles through classical acceleration, and then radiated to hard photons. When the standing wave is interacting with a plasma target, this is one of the possible depletion channels. Let us assume that our target is 1 µm thick and the density n is much lower than the relativistic critical density a 0 n c . In that case, we do not expect the standing wave to be disturbed. Also, by taking a narrow plasma slab, all plasma electrons (and eventually positrons) are located within one laser wavelength. This limits the interaction of each section of the traveling waves with the plasma to one full laser period. The radiative laser depletion can therefore be estimated locally in a λ 3 0 volume and during one full laser cycle. The percentage of the wave turned into high-frequency radiation depends on the local plasma density and the local laser intensity. The density gives an estimate of the total number of particles in the λ 3 0 volume, while the emitted energy per particle depends on the laser intensity and is given by Eq. (9).
If n = n c and we consider an electron-positron plasma, we have about 10 9 electrons and 10 9 positrons in the λ 3 0 volume. The fraction of the laser energy absorbed locally during one laser period is then given by (10) Here, multiplication factor 4 comes from moving from half of T L to a full T L and considering two species, electrons and positrons. If we now follow any point of the individual traveling waves at the speed of light as it passes the interaction region, we note that any absorption that happens to that portion of the wave happens within one T L . The rest of the time, the wave propagates freely. Every T L , there is a fresh portion of the wave interacting with the plasma, while the previous has escaped the interaction region. This means that the ratio given by Eq.  For optical lasers with λ 0 = 1 µm, interacting with a pair plasma we get that in an undisturbed standing wave, the energy converted to high frequency radiation can be approximated as This estimate is not the whole picture, because we neglect the self-consistent fields in the plasma and pair production. However, Eq. (11) allows to estimate how strong is the high-frequency emission as a depletion channel when considering plasma densities much lower than the relativistic critical density a 0 n c . Figure 5 illustrates the strong radiative absorption limit as a function of intensity as predicted by Eq. (11). The characteristic density is calculated for each laser intensity assuming 10% laser absorption and solving for n. If we would like to stay in the regime of low absorption (undisturbed wave), we have to use the density below this limit. It is interesting to note that one does not need to reach the relativistic critical density for extreme absorption according to this. For a 0 = 1000 even n = 160n c is enough for strong absorption. Cryogenic targets available in experiments have a density on the order of 10 n c , which means it is possible to use them for controlled radiation sources, provided that the QED cascade does not increase the plasma density by an order of magnitude.
To evaluate a parameter range where this is the case, we should analyze also the pair production. Plasma density in a cascade increases exponentially: n = n 0 exp(Γt), where Γ is a growth rate that depends on the laser intensity. For a 0 < 500, and circularly polarized lasers, the growth rate has an approximate expression given by 47 The growth rate predicted by Eq. (12)  400. An upper estimate for linear polarization can be estimated by assuming the same total energy yields the same growth rate. This is equivalent to using a 0 / √ 2 instead of a 0 in Eq. (12) Another option for Gaussian laser pulses is to insert a 0 /2 instead of a 0 in Eq. (12) to account for the fact that not all particles experience the instantaneous maximum intensity Using a 0 /2 is also more consistent with the data obtained from cascade simulations at higher intensities 29,47 (a 0 > 1000). Equations (13) and (14) define a range of values for expected growth rate at each a 0 for linear polarization. We should note here that near-threshold pair production is very sensitive to the seeding 46,47,69 , so one may not be always able to reach the given growth rates. When using laser pulses, the actual growth rate is likely to be closer to Γ LP − than to Γ LP + . For example, the multiplicity we get for a 24.5 fs laser pulse with peak intensity of a 0 = 500 using Eq. (14) is 1.44, while the simulation data give 1.69. For the same data, the upper bound Γ LP + predicts a multiplicity of ∼ 100, which is two orders of magnitude higher. Using an infinite plane wave, the multiplicity in the simulation increases, but it stays on the same order as predicted by Γ LP − in Eq. (14). The growth rate estimates for a 0 < 500 given by Eqs. (12)- (14) are illustrated in Fig 6. As we mentioned before, for designing a radiation source, it may be favourable to keep the pair production rate low and operate in a more controlled emission regime. To avoid the uncertainty of seeding, it is also important to have a reliable estimate regarding how many particles are contained within the interaction volume. In this sense, using a gas jet would not be optimal, even though it has an advantage of offering a possibility to start at a low density. A dense solid target can be manufactured to high precision, but then an undisturbed standing wave interaction is not possible. The best option available to date could be to use cryogenic targets, composed of hydrogen ice. They have recently been used in ion acceleration experiments, and µm-level thickness was achieved 70,71 at a density ∼ 10 n c .
In the following paragraphs, we consider the output radiation generated by interaction of two pulses with one such target, initially composed of electrons and protons. For completeness, we compare the findings also with the case when the target is initially composed of electrons and positrons, as well as the case when the standing wave is undisturbed. A range of intensities between a 0 = 100 and a 0 = 1500 was considered. Figure 7 displays the angular distribution of radiation for all the cases, while Fig. 8 shows the conver- sion efficiency of optical light to high-frequency radiation. The differences between the electron-proton and electron-positron target in the typical radiation directions are almost negligible, although there is a difference in the energy conversion efficiency. This is not surprising, because in the electron-positron target, there are twice as many radiating leptons. The angular distribution in an undisturbed standing wave is quite different at some intensities compared to that of the 10 n c target. For a 0 = 100, the 10 n c target represents already 10% of the relativistic critical density for such a wave. Some level of discrepancy is therefore expected for the lower end of the explored intensities. At higher intensities, the discrepancies come from the pair production that increases the target density during the interaction, until the density is high enough for wave disruption. Mid-range intensity a 0 ∼ 500 seems to be the best choice for controlled emission, because the pair multiplicity is low, and at the same time the intensity is high enough not to be too disturbed by the presence of the target. Figure 9 shows a 3D simulation using a cryogenic target for a 0 = 1000. The conversion efficiency obtained in 3D simulations for a 0 = 100, a 0 = 500 and a 0 = 1000 are displayed together with 2D results in Fig. 8. The absorption is somewhat lower in 3D (as the average laser intensity is lower), but is of the same order of magnitude, so the conclusions regarding the different regimes of interaction are consistent with these results as well.

V. CONCLUSION
We have shown that it is possible to convert most of the energy from the laser (or interacting electrons in a laserelectron scattering) to γ-rays in the near-future laser experiments. To achieve a strong conversion efficiency, one should either use a solid-density thin target, or use intensities that can initiate a QED cascade that produces enough pairs to increase the plasma density during the interaction (a 0 ∼ 1000). It is also possible to obtain radi- ation emission in controlled conditions. Using cryogenic targets, this is possible for a 0 ∼ 500, where the number of particles is not expected to increase more than twice due to the pair production, and the target density n ∼ 10n c is small enough to keep the radiative absorption below 10%.

VI. APPENDIX
All simulations are performed with the QED module of OSIRIS 72 . The QED module is a Monte-Carlo module that accounts for the photon emission and Breit-Wheeler pair production, and is implemented as an addition to the standard PIC loop of OSIRIS. Energetic photons are initialized as an additional particle species. The emission rates are found in Refs [73][74][75] . OSIRIS QED has been used previously in Refs. 16,29,44,47,76 . A similar method for incorporating BW pair production is used in several other codes 28,31,32,[34][35][36][37][38]77 .
The simulations from Section III displayed in Fig. 7 are performed with parameters as in Vranic et al 16 . The electron beam initial energy is 0.85 GeV, and initial beam divergence is p ⊥ /p ∼ 0.2 mrad. The laser is transversely a plane wave with a temporal envelope. The total pulse duration is given by τ = τ flat + (τ rise + τ fall )/2, where τ flat is the constant amplitude section of the wave that was varied from 0 to 160 fs. The envelope function has a smooth rise and fall, the same for all the simulations: τ rise = τ fall = 26.6 fs. The simulations are performed in 2D, with a box size of 500 × 20 c 2 /ω 2 0 resolved with 5000 × 200 cells and the timestep dt = 0.04 ω −1 0 using 16 particles per cell (ppc).
Section IV has several types of simulations. Ideal simulations displayed in Figs. 3, 4 and 7 were performed with no current deposition and pair production: the standing wave is undisturbed by the presence of the plasma and the plasma density does not grow due to the new particle generation. The transverse boundary conditions are periodic. The simulation box size is 200 × 10 c 2 /ω 2 0 , resolved with 2000 × 100 cells and 9 ppc. The plasma slab is 1 µm thick, and is composed of electrons and positrons. The initial density is n = 0.001 n c . The two lasers have a 2cycle smooth rise and fall, and a 10-cycle flat section. All the measurements in ideal conditions were taken while all particles are fully immersed in the flat-top section of the