Gamma photons and electron-positron pairs from ultra-intense laser-matter interaction: A comparative study of proposed con ﬁ gurations

High-energy γ -photongenerationvianonlinearComptonscatteringandelectron – positronpaircreationviatheBreit – Wheelerprocessdrivenby laser – plasma interaction are modeled, and a number of mechanisms are proposed. Owing to the small cross section, these processes require both an ultra-intense laser ﬁ eld and a relativistic electron bunch. The extreme conditions for such scenarios can be achieved through recent de-velopments in laser technology. Photon emission via nonlinear Thomson and Compton scattering has been observed experimentally. High-energy positron beams generated via a multiphoton process have recently been observed too. This paper reviews the principles of γ -ray emission and e + e − pair creation in the context of laser – plasma interaction. Several proposed experimental setups for γ -ray emission and e + e − pair creation by ultra-intense laser pulses are compared in terms of their ef ﬁ ciency and the quality of the γ -photonand positron beams produced forultrashort (15 fs) and longer (150 fs) multi-petawatt laser beams.


I. INTRODUCTION
The 2018 Nobel Prize in Physics was awarded to Ashkin, Mourou, and Strickland for their groundbreaking inventions in the field of laser physics. The chirped pulse amplification (CPA) technique proposed by Strickland and Mourou 1 provides significantly enhanced laser intensities in excess of 10 18 W cm −2 . The increasing intensity pushes the particle dynamics in the laser field into the relativistic region, thus providing a basis for the development of laserdriven accelerators. Electrons in the plasma with appropriate phases are accelerated to relativistic energies within millimeters by the wakefield or directly by the laser field. Several regimes have been proposed for the production by these relativistic electrons of radiation ranging from X rays to γ rays. State-of-the-art laser facilities 2,3 are able to deliver intensities of the order of 10 22 W cm −2 . The next generation of lasers are designed for the 10 PW region, with the corresponding intensities being expected to reach the order of 10 23 -10 24 W cm −2 or even higher. 4,5 As a result of interaction with such high-intensity fields, the emitted γ photons have a momentum comparable to the that of the electrons, and quantum electrodynamic (QED) effects then become significant. Processes generating secondary sources of particles, such as γ-photon emission, electron-positron pair creation, and QED cascades, come into play under these extreme intensity conditions, and the corresponding phenomena have been widely studied.  Experimental signatures of radiation reaction and QED processes have been reported recently. 29,30 The regimes for γ-ray emission in laser-plasma interaction include bremsstrahlung, synchrotron emission, nonlinear Thomson scattering, and nonlinear Compton scattering. Bremsstrahlung [31][32][33] dominates in the case of solid targets at laser intensities below 10 22 W cm −2 . Giulietti et al. 34 reported an intense γ-ray source around the giant dipole resonance for photonuclear absorption obtained via electron-bunch bremsstrahlung driven by a 10 TW laser. Vyskočil et al. 35 described γ-ray emission in laser-foil interactions using kinetic simulations and found that the corresponding γ-ray spectra show a linear dependence of the γ-ray temperature on the normalized laser potential. With highenergy electron beams, γ photons can be obtained via nonlinear Thomson scattering. Sarri et al. 36 reported experimental observations of ultrahigh-brilliance [10 20 s −1 mm −2 mrad −2 (0.1%BW) −1 ] multi-MeV γ rays from nonlinear relativistic Thomson scattering. High-energy electrons in the self-generated magnetic field emit multi-MeV photons, as described by Stark et al. 37 Yan et al. 38 presented measurements of high-order multiphoton Thomson scattering from a laser-driven electron beam. Recently, solid nanowires have been reported to be ideal materials for γ-ray generation. Martinez et al. 39 presented a numerical study of the synchrotron emission resulting from the interaction of a strong laser field with a nanowire array. Wang et al. 40 proposed the use of an electron wiggling regime to obtain collimated γ rays along a petawatt laser-irradiated wire. Nonlinear Compton scattering 41,42 is the main mechanism for γ-ray production in the case of relativistic electrons interacting with ultra-intense electromagnetic (EM) fields. Chang et al. 20 proposed the use of circularly polarized laser pulses to produce photon emission with high brightness in a plasma channel. According to Benedetti et al., 43 collimated γ-ray flashes with high brilliance [10 25 s −1 mm −2 mrad −2 (0.1%BW) −1 ] are obtained via plasma filaments as a result of synchrotron emission from an ultrarelativistic electron bunch traveling in a millimeter-scale conductor. Recently, high-order Laguerre-Gaussian mode laser and compound targets 44 have been used to induce nonlinear Compton scattering. A tunable γ-ray beam generation regime using an oscillating plasma mirror has also been proposed recently. 45 In the case of an extremely strong EM field interacting with γ photons, e + e − pairs are produced by the multiphoton Breit-Wheeler (BW) process 46,47 γ + nℏω 0 → e − + e + , where ω 0 is the laser frequency and n is the number of laser photons. The small cross section for photon-photon collisions means that the threshold field intensity for BW process can be as high as 10 24 W cm −2 . 7 QED cascades in the multiple counterpropagating regime have been proposed in recent theoretical studies. 9,19,21,48,49 Regimes involving pair creation with a single laser pulse have also been proposed on the basis of a strong electron self-injection effect 18 and circular oscillations in the plasma channel. 50 Normal injection of a multi-PW laser into a bunch of GeV-class electron beams has been reported to produce a low-divergence electron-positron pair beam. 51 A new mechanism based on a relativistic plasma mirror has recently been proposed 52 in which the pair creation efficiency and the production rate are significantly enhanced.
In this paper, we briefly introduce the laser-plasma interaction in the ultrarelativistic regime, together with the corresponding radiation and QED effects. The recent proposed regimes for γ-ray emission and pair creation are reviewed and compared. The scenarios for each setup are presented using particle-in-cell simulations. The advantages and disadvantages of this approach are summarized.

II. GENERAL DESCRIPTION OF ULTRA-INTENSE AND ULTRASHORT LASER-MATTER INTERACTION
The intensity at which the electromagnetic field exceeds the binding strength of an electron to an atom is about I a ce 2 /(8πa 4 B ) ≈ 3.5310 16 W cm −2 , where a B ℏ 2 /(m e e 2 ) is the Bohr radius and other physical constants have their usual meanings (c is the speed of light, e the electron charge, m e the electron mass, and ℏ the reduced Planck's constant). An ultra-intense laser pulse can ionize any target even without account being taken of multiphoton and tunneling ionization effects, which actually reduce the threshold value to the order of 10 14 W cm −2 . The prepulse or the pedestal of an ultra-intense laser is strong enough to ionize the target material. The strength of the laser pulse is scaled by the dimensionless parameter a 0 eE 0 /m e ωc, where E 0 is the laser electric field and ω is the angular frequency. The relation between the laser intensity and the normalized amplitude is a 0 ≈ 8.57310 −10 I(W cm −2 ) λ(μm). Laserplasma interactions in the intensity range of 10 18 -10 23 W cm −2 are discussed below.

A. Relativistic laser-plasma interaction and particle acceleration
When the dimensionless amplitude a 0 1, this means that the field provides work equal to m e c 2 within a distance λ/2π. Relativistic effects must be considered when a 0 ≫ 1, which implies that the electron can be accelerated to relativistic energy in the distance of one laser period. In such a relativistic laser-plasma interaction, electron thermal motion can be neglected, since it is much weaker than the motion driven by the laser field. The transverse and longitudinal motion of a free electron in a planar laser field satisfies u ⊥ a and u ∥ γ − 1 a 2 0 /2. Here, a eA/mc 2 is the instantaneous normalized vector potential and u ⊥ and u ∥ are respectively the transverse and longitudinal components of the electron momentum normalized by m e c. In the case of a ≫ 1, the longitudinal momentum is much larger than the transverse momentum. Therefore, the response of an electron to a relativistic laser field is dominated by its longitudinal motion driven by the ponderomotive force. It is also clear that there is no energy gain by an electron passing through a plane wave, since the vector potential is zero before and after the pulse. Energy gain is possible, for example, with a focused laser pulse in which the electron experiences different field strengths at different distances from the focal point.
Based on single-electron dynamics, we discuss the interaction of a laser and an underdense plasma with a density lower than the critical density n cr m e ω 2 /4πe 2 . A plasma with electron density higher than the critical density becomes opaque to an EM field with frequency ω. Taking the relativistically induced transparency effect into account, the cutoff density is increased to γn cr as the electron mass increases by the Lorentz factor γ. The ponderomotive force of the laser field propagating in an underdense plasma and the restoring electrostatic force in the plasma provide the possibility of creating a largeamplitude wakefield. This laser wakefield excitation is an important process in which electrons can be trapped and accelerated in an efficient way. A sinusoidal type of density perturbation δn/n 0 is induced by a laser pulse in a plasma with an initial density of n 0 . The perturbations are accompanied by an electrostatic field (wakefield) with maximum amplitude E max ∼ 2(γ − 1) . 53 When the optimal conditions are satisfied, the maximum energy gain for the electrons trapped inside the wakefield is proportional to γ } cτ P L √ , where τ and P L are the pulse length and peak power of the laser.
Electrons in the plasma can also be trapped by the laser field and heated via the direct laser acceleration (DLA) mechanism. 54 In this case, the electrons experience betatron oscillations in the transverse direction driven by the laser field. In the ultrarelativistic case, the electric and magnetic components of the Lorentz force are of the same order of magnitude, since the corresponding transverse velocity is of the order of the speed of light. The electrons are accelerated by the u ⊥ /γ 3 B force in the longitudinal direction and comoving with the laser field.

B. Radiation reaction and QED
With increasing intensity and particle energy, the charged particle dynamics becomes dissipative owing to the emission of highenergy photons. Here, we introduce a dimensionless parameter ϵ rad 4πr e /3λ, where r e e 2 /m e c 2 is the classical electron radius. Nonlinear Thomson scattering has emission power proportional to the fourth power of the electron energy: P γ ≈ ϵ rad m e c 2 ωγ 4 . Radiation effects become important when the energy gain per unit time (P e a 0 ωm e c 2 ) is balanced by the radiation loss, i.e., a 0 > ε −1/3 rad . The corresponding laser intensity is I R ≈ 2.65 3 10 23 (1 μm/λ) 4/3 W cm −2 . 55 The electron dynamics is significantly affected by the radiation reaction due to the recoil momentum from the emitted photons, which is equivalent to the radiation damping force in the Landau-Lifshitz form 56 In the case of the electron DLA mechanism, the radiation damping force prevents the electrons from escaping from the pulse center.
When the emitted photon energy is comparable to the electron kinetic energy, i.e., ℏω γ ≈ γ e m e c 2 , QED effects must be taken into account. The threshold for these effects is determined by the dimensionless relativistic and gauge-invariant parameter the EM field, p ] is the electron 4-momentum, and E s m 2 e c 3 /eZ is the critical QED electric field, i.e., the Schwinger limit field. In the case of a plane EM wave, the invariant parameter can be expressed as χ e (E/E s )(γ − p ∥ /m e c). Considering γ ≫ 1 under the ultrarelativistic condition, the corresponding parameter is maximized in a counterpropagating geometry of EM field and electrons, χ ↑↓ e ≈ 2γ(E/E s ). However in the copropagating case, it dramatically decreases to χ ↑↑ e ≈ (2γ) −1 (E/E s ). High energy photons produce e + e − pairs in the extremely strong field via the multiphoton BW process γ + nℏω → e − + e + . An important parameter in calculating the probability of this QED process is where ℏk ] (ℏω γ , ℏk γ ) is the 4-momentum of the high-energy photon. It can be expressed in the threedimensional geometry as Similar to the case of γ-photon emission discussed above, χ γ can be calculated in the case of photons interacting with the plane EM field: The counterpropagating and copropagating photons have χ ↑↓ γ ≈ 2(Zω/m e c 2 )(E/E s ) and χ ↑↑ γ ≈ 0, respectively.

III. MECHANISMS OF γ-RAY AND ELECTRON-POSITRON PAIR PRODUCTION
In this section, we review several recently proposed mechanisms for generating γ rays and e + e − pairs. As mentioned above, the probabilities for γ-ray emission and e + e − pair creation are both inversely proportional to the Schwinger limit field E s ≈ 1.32 3 10 18 V m −1 .
Assuming that the electric field strength reaches ∼10 15 V m −1 , corresponding to a laser intensity of 10 24 W cm −2 , the ratio E/E s ∼ 10 −3 . This implies that the cross sections for the QED processes are still exponentially small. The way to increase the parameters χ e and χ γ to χ e,γ ∼ 1 is to use an optimal geometry and enhance the energy of the electrons and the photons. To satisfy the optimal conditions, various setups have been proposed and studied intensively in recent years. We classify these mechanisms according to the setup and summarize the achievements. For each category, we present our simulations using the 2D PIC code EPOCH 57,58 to illustrate the main results. The laser wavelength λ 1 μm and the power chosen in the simulations is 70 PW for the mechanism employing a single pulse or 35 PW in the case of two laser pulses. The focal spot size is 3 μm, and thus the peak laser intensities are 5 3 10 23 W cm −2 and 2.5 3 10 23 W cm −2 , respectively. Two different pulse lengths are compared for each mechanism: a short pulse of 15 fs and a long pulse of 150 fs.

A. Multicolliding laser pulses
Pair creation by multicolliding pulses in vacuum was proposed in Ref. 9, where it was theoretically predicted that the pair creation threshold could be reduced by concentrating the EM field energy into a small volume. Since then, there have been several presentations of the use of multiple counterpropagating pulses interacting with a plasma target. 19,21,48,49,59 Electrons accelerated in a plasma generate a large number of γ photons via nonlinear Compton scattering, and these photons interact with the strong field of the standing wave formed by the laser pulses from both directions. BW processes occur, and e + e − pairs are created.
In our case, two counterpropagating linearly p-polarized laser pulses interact with a circular target lying at their common focal spot [ Fig. 1(a)]. The target is composed of electrons and deuterium ions, and its diameter is 3 μm. The simulation box of dimensions 40 μm 3 100 μm is resolved with 1200 3 3000 cells, each containing 100 simulation particles. As the two counterpropagating laser pulses overlap, a transient standing wave is created. The electrons oscillating in the electromagnetic standing wave serve as a source of γ photons. These emitted photons then interact with photons of the laser field, which results in e + e − pair creation via the BW process. For a given laser pulse, there exists an optimal initial target density such that the maximum number of seed electrons remain in the focal spot until the highest intensity of the standing wave is established. In our case, the highest number of positrons created per seed particle is achieved when a target of initial density 50n cr (200n cr ) is used for the interaction with two colliding laser pulses having duration 15 fs (150 fs). Since both pulses have the same peak intensity, the longer one delivers more energy, which leads to more efficient positron production. In this case, 0.01 positrons per seed electron are created, while in the case of a short pulse, the number of positrons is ten times lower. The absolute number of created photons and positrons is much higher in the case of a longer laser pulse [ Table I and Fig. 1(b)]. This is due to the four times higher number of seed particles and also to the much longer duration of the established standing wave. Therefore, a greater fraction of the laser energy is converted into photons and positrons in the long-pulse case. The probability of pair production depends on both the magnitude of the electric field and the photon momentum, as well as on their mutual orientation. Since the intensity and the mean photon energy are comparable in both cases, the mean positron energy is expected to be similar for both laser pulse durations (Table I) As the electric field of the standing wave oscillates in the direction of the laser pulse polarization, the created positrons are accelerated predominantly in this direction, as shown in Figs. 1(d) and 1(f). As the number of created positrons is much lower in the case of a shorter laser pulse, the actual phase of the oscillating standing wave is more pronounced in Fig. 1(d), resulting in an asymmetric distribution. After one-half of a laser period later, the electric field of the standing wave will acquire the opposite phase, and consequently will push positrons in the upward direction.

B. Single laser with plasma mirror reflection
In one of the pioneering works on hard X-ray emission via nonlinear Thomson scattering, Phuoc et al. 11 proposed a regime  employing a solid target serving as the reflection mirror at the end of a gas target. Here, we demonstrate the possibility of using the reflection mirror for enhancing the nonlinear Compton scattering and the BW process, which are discussed in detail in Refs. 45 and 52. The simulation box has size 60λ 3 40λ, with a resolution of 100 cells per λ. As can be seen in Fig. 2(a), the target is composed of a near-criticaldensity (NCD) plasma and an ultrathin solid layer. The NCD target has a density of 4n cr and is 40λ in size. The main mechanism of electron acceleration depends on the target density. At NCD, direct laser acceleration dominates. The electrons are accelerated by the laser field, reaching the GeV level in a short distance. The solid target reflects the laser field, serving as a relativistic oscillating mirror. The solid target consists of a gold ion layer with a charge state of 11+ and free electrons with a corresponding density of 600n cr . The thickness of the target is 4λ. With the Doppler shift and self-focusing effects, the reflected intensity can be enhanced significantly. There is a head-on collision geometry between the reflected laser field and the highenergy accelerated electrons. The corresponding probabilities for nonlinear Compton scattering and the BW process become higher than in the case of electron self-injection. 18 The spectrum in Fig. 2(b) indicates that the number of emitted γ photons in the τ 150 fs case (black) is larger than that in the τ 15 fs case (blue). However, the short-pulse case provides more highfrequency photons, and the cutoff energy reaches more than 1.5 GeV. The long pulse accelerates more electrons, since it occupies a much larger spatial area compared with the short pulse. Therefore, more electrons have the chance to emit γ photons. On the other hand, the laser ponderomotive force is inversely proportional to pulse length. A stronger ponderomotive force is applied to the electrons in the τ 15 fs case and accelerates them to higher energy. This can be seen from the tail of the γ-photon distribution function. The total numbers of photons in the long-and short-pulse cases are 10 20 μm −1 and 1. the photons and the positrons are well collimated, with a narrow opening angle. However, the long-pulse results present a broadened and isotropic distribution of both γ photons and positrons, as can be seen in Figs. 2(e) and 2(f). The difference in beam quality comes from the field reflected by the plasma mirror. The solid target pushed by the incident pulse reflects the EM field via the relativistic oscillations of the electrons around the ion layer. The velocity of the ion layer (also called the piston velocity in the piston model) can be estimated by balancing the radiation momentum flux and the charged-particle momentum flux. Under the density conditions considered here, the velocity β f is of the order of 0.1c. The corresponding reflection process is regular in the short-pulse case owing to the minor displacement of the ion layer, δs ∼ τβ f ∼ 0.5λ. However, the long pulse significantly distorts the solid target, with a displacement as large as 5λ. Therefore, the long pulse is not favorable to being reflected by a plasma mirror in the high-intensity regime. Despite the fact that it provides relatively high numbers of γ photons and positrons, the energy and the beam quality are not comparable to those obtained in the short-pulse case.

IV. COMPARISON OF THE SCHEMES
In this section, we compare the simulation results obtained in the different regimes described above. The mean energy, the number of secondary particles and the laser energy conversion efficiency are summarized and compared in Table I. The mean energy of the γ photons and e + e − pairs are similar in all of the setups, which indicates that the mean energy is related to the amplitude of the laser field instead of the total energy. It should be mentioned that the energy of the positrons in the case with the plasma mirror and short laser pulse can be greatly enhanced in the later stage, since the positrons are captured by the reflected field, experiencing further acceleration during their propagation. The longer pulses provide a higher number of secondary particles. In the case of a standing wave formed by two colliding laser pulses, the use of long pulses is preferable, since it increases the spatial area occupied by the strong field of the standing wave. The production rates for both γ photons and positrons are much higher with the long pulse. However, in the setup using a single pulse and a plasma mirror, the simulation results show the opposite behavior. Although the absolute number is lower, the efficiencies are much higher with τ 15 fs. The main reason for this is the destruction of the plasma mirror by the long pulse, which reduces the reflection efficiency and collimation. For a given incident laser energy, the short pulse with a plasma mirror provides the highest production rate for γ photons and positrons. The beam qualities presented above are the best results for this regime too.
In conclusion, for the generation of γ photons and e + e − pairs with longer petawatt laser pulses, the standing wave configuration formed by two counterpropagating laser pulses provides better results in terms of the laser energy transformation efficiency. This results in the highest number of e + e − pairs being created in this configuration. For shorter petawatt laser pulses, the plasma mirror setup is beneficial, since it significantly enhances the production rate by more efficient laser pulse reflection and focusing. Furthermore, it provides the ideal beam qualities for further utilization of the γ-ray flash and the positron beam. The laser parameters discussed in this paper may be attained at upcoming facilities like ELI-Beamlines, 4 and thus these results are useful for designing future experiments.