Bo Zhang1,2, Zhi-Meng Zhang3,4, Zhi-Gang Deng3, Wei Hong3,4, Jian Teng3,4, Shu-Kai He3, Wei-Min Zhou3,4, Yu-Qiu Gu5,6. 1. Department of High Energy Density Physics, Research Center of Laser Fusion, 621900, Mianyang, Sichuan, People's Republic of China. zhangbolfrc@caep.cn. 2. Laboratory of Science and Technology on Plasma Physics, Research Center of Laser Fusion, 621900, Mianyang, Sichuan, People's Republic of China. zhangbolfrc@caep.cn. 3. Department of High Energy Density Physics, Research Center of Laser Fusion, 621900, Mianyang, Sichuan, People's Republic of China. 4. Laboratory of Science and Technology on Plasma Physics, Research Center of Laser Fusion, 621900, Mianyang, Sichuan, People's Republic of China. 5. Department of High Energy Density Physics, Research Center of Laser Fusion, 621900, Mianyang, Sichuan, People's Republic of China. yqgu@caep.cn. 6. Laboratory of Science and Technology on Plasma Physics, Research Center of Laser Fusion, 621900, Mianyang, Sichuan, People's Republic of China. yqgu@caep.cn.
Abstract
Electron radiation and γ photon annihilation are two of the major processes in ultra intense lasers (UIL). Understanding their behavior in one coherence interval (CI) is the basis for UIL-matter interaction researches. However, most existing analytic formulae only give the average over many CIs. Present understanding of these two multi-photon processes in one CI usually assume that they emit forward and their spectra have a cutoff at the energy of the electron/γ. Such assumptions ignore the effects of involved laser photons (EILP). We deduced the formulae for these two processes in one CI with EILP included and give the conditions for the EILP to be significant. Strong EILP introduces new behaviors into these two processes in one CI, such as large angle emission and emit particles above the usually assumed cutoff. Simulations show that the EILP would be significant when laser intensity reaches 2 × 1022 W/cm2, which is within the reach of state-of-art lasers.
Electron radiation and γ photon annihilation are two of the major processes in ultra intense lasers (UIL). Understanding their behavior in one coherence interval (CI) is the basis for UIL-matter interaction researches. However, most existing analytic formulae only give the average over many CIs. Present understanding of these two multi-photon processes in one CI usually assume that they emit forward and their spectra have a cutoff at the energy of the electron/γ. Such assumptions ignore the effects of involved laser photons (EILP). We deduced the formulae for these two processes in one CI with EILP included and give the conditions for the EILP to be significant. Strong EILP introduces new behaviors into these two processes in one CI, such as large angle emission and emit particles above the usually assumed cutoff. Simulations show that the EILP would be significant when laser intensity reaches 2 × 1022 W/cm2, which is within the reach of state-of-art lasers.
Laser is presently the most intense electromagnetic field in laboratory. On next generation 10–100 petawatt laser facilities such as ELI[1], laser intensity is anticipated to reach 1023–25W/cm2, which is several orders higher than the present record of 2 × 1022 W/cm2 [2]. Such strong lasers will open a gate for many fundamental and fantastic strong field quantum electrodynamics (SFQED) phenomena[3-7]. For recent reviews, see[8-10].At such intensity, laser-matter interaction enters a new regime where two multi-photon SFQED processes, the quantum process of electron radiation (nonlinear Compton scattering, NCS)and γ photon annihilation into an electron-positron pair (nonlinear Breit-Wheeler process, NBW)become significant, i. e., can take a large fraction of laser energy (p, k, etc. are 4 vector momenta).In ultra intense lasers, the CI of NCS and NBW is very short, where the normalized field amplitudee and m are electron charge and mass, a, E and k0 are the 4-vector, electric field and photon energy of the laser. Short CI means NCS and NBW happen in different CIs are independent from each other. NCS/NBW within a CI is then governed by the local field F and instantaneous energy-momentum of the electron/γ.Secondary interactions and back reactions are usually very strong in UILs, many emitted γ photons can produce e− − e+ pairs and the produced charges can radiate, etc. NCSs can also change the electron and positron energy and momentum significantly within a CI. Due to short CI, strong secondary interactions and strong back reactions, the influence of NCS and NBW on UIL-matter interactions usually have to be investigated with simulations rather than analytic methods. In such simulations, NCS and NBW are treated as instantaneous and local processes.Behavior of NCS/NBW in one CI is then the basis for UIL-matter interaction researches. Although NCS and NBW have already been studied extensively and a lot of references have discussed EILP, what most of these studies gave are the average over CIs of a period or a pulse with secondary interactions and back reactions ignored. Corresponding formulae cannot describe the behavior of NCS and NBW in one CI.On the other hand, unfortunately, existing analytic studies of NCS and NBW within one CI only gave lightcone results. The electron emission probability[11] in one CI in UILs () iswhere , is the modified Bessel function of the νth order, α is the fine structure constant, p0 is the electron energy and . The pair production rate in one CIwhere , and . These lightcone results are functions of and , respectively. Hence many important features of NCS and NBW in one CI such as spectrum and emission directions are obscure.As an alternative, present UIL-matter interaction researches usually assume the γ photon emitted by NCS in one CI is along the instantaneous electron forward direction. The basis for this assumption on such a quantum process is that the emission angle of an ultra relativistic charge ~1/γ according to classical electrodynamics is very small[8,12,13]. The spectrum of NCS is also assumed to cut-off at the instantaneous electron energy p0. The basis for the second assumption is that, when the involved laser photon number n is small enough, considering the laser photon energy , would be the upper limit for both and .Based on these two assumptions, present understanding of NCS in one CI usually take the emitted γ photon momentum and the changed electron momentum [14-16]. Similar assumptions for NBW give and [14-16].Present understanding of NCS and NBW in one CI is a combination of the lightcone probabilities in Eqs 4 and 5 and the two assumptions. The model for NCS and NBW based on such understanding is widely accepted and applied in UIL-matter interaction researches[17-26]. As a consequence, existing simulation studies of UIL-matter interactions ignored EILP.However, the numbers of laser photons involved in these two processes, n and n′ scale as ~[8,11,27,28], therefore the energy and momentum of involved laser photons grow nonlinearly with laser intensity. When laser intensity reaches I ~ 1024 W/cm2, the energy of laser photons involved in a single NCS nk0 or in a single NBW n′k0 would reach the scale of electron vibration energy a0m or the laser wake field acceleration record of 4.2 GeV[29]. According to energy and momentum conservation, the emission angle and spectrum of these two processes must be strongly influenced at such intensity. A simple diagram for EILP in NCS is shown in Fig. 1.
Figure 1
EILP and momentum conservation of NCS in a short CI in UILs. The left depicts present model and the right includes EILP.
EILP and momentum conservation of NCS in a short CI in UILs. The left depicts present model and the right includes EILP.To investigate EILP in UIL-matter interactions, EILP on NCS and NBW in one CI is very important. The number of laser photons involved in a single NCS in UILs has a distribution peaked atIts dispersion is very small in UILs.Without loss of generality, we put the electron momentum p on the z axis and k on x–z plane as shown in the left panel of Fig. 2. The angle between p and k is θ and the angle between a and p − k plane is .
Figure 2
Geometry of NCS (left) and NBW (right) in general frames.
Geometry of NCS (left) and NBW (right) in general frames.Solve the energy-momentum conservation equation of NCS, Eq. (6), the definition of u and an additional transverse distribution approximation simultaneously (see the supplemental information for details), where , the momenta of the deflected electron and the emitted γ photon areandwhereand . The ± subscripts are the signs ofThe chances for + and − are equal because the differential probability is a even function of [11]. Emitted photon energy and the emission angle areFigure 3 shows the emission angle of a 1 GeV electron when it radiates at the peak of an UIL. Figure 3(a) is the head-on case, it is larger than 30° at I = 1024 W/cm2. Figure 3(b,c) show emission angles θ and θ when . Figure 3(d) gives the emission angle when k, p and a are on the same plane. It shows that when is small, the emission angle can be very large. This is quite different from present understanding of NCS in one CI in UILs where forward emission (vanishing emission angle) is assumed[14-26].
Figure 3
Emission angle of NCS in one CI, θ when k and p are antiparallel (a), θ when (b), θ when (c) and θ when and (d). Corresponding parameters are , (λ = 1 μm) and W/cm2.
Emission angle of NCS in one CI, θ when k and p are antiparallel (a), θ when (b), θ when (c) and θ when and (d). Corresponding parameters are , (λ = 1 μm) and W/cm2.This non-vanishing fixed emission angle with respect to instantaneous electron momentum is different from that in[27] which is with respect to the electron momentum at infinity where laser vanishes, i. e., before it enters the laser. Additionally, the emission angle in one CI is fixed in arbitrary frames while that in ref.[27] is fixed only in specific frames that k and the electron momentum at infinity are anti-parallel. This angle is also different from those induced by stochastic effects and consecutive radiations in refs[24,26].The spectrum of the emitted photon is also modified by EILP. Comparison of energy integrated spectra of a 1 GeV electron at the moment it appears at the peak of an UIL of I = 1024 W/cm2 including EILP and excluding it is shown in Fig. 4. As given in Eq. (11), EILP blue shifts half of the spectrum and red shifts the other half. Hence the spectrum gets a two-stage structure and the blue shifted part can surpass the instantaneous electron energy of p0 = 1 GeV, which is the cutoff in present model of single NCS in UILs.
Figure 4
Energy integrated spectrum of a 1 GeV electron including and excluding EILP (blue solid, red dashed). W/cm2, , for the left panel and , for the right.
Energy integrated spectrum of a 1 GeV electron including and excluding EILP (blue solid, red dashed). W/cm2, , for the left panel and , for the right.EILP has similar influences on NBW in one CI. The involved laser photon number n′ is approximatelyand its dispersion is of scale. As shown in the right panel of Fig. 2, without loss of generality, we fix the γ photon momentum k″ on the z axis and k on x–z plane. The angle between k″ and k is θ and the angle between a and k″ − k plane is .Similar to the NCS case, solve the energy-momentum conservation of NBW, Eq. (12), definition of δ and again the additional approximation simultaneously (see the supplemental information for details), the emitted electron momentum isand the emitted positron momentum iswhereandAgain, ± subscripts are the signs ofand both signs take half of the chance.EILP on a NBW in UILs is similar but different. When laser intensity is comparatively low, the e− spectrum with EILP agrees with the present model as Fig. 5(a) shows, note that the spectrum of e+ is the same. With the growth of laser intensity, as shown in (b) and (c), the EILP introduces a lower limit and this is due to the minimal number of laser photons involved in a NBW that can be easily deduced from Eq. (12). Finally, when laser intensity is very high, as (d) shows, the cutoff at disappears.
Figure 5
EILP on the spectrum of electron produced by NBW (red lines). The parameters are: , , GeV, eV (λ = 1 μm) and W/cm2 is 1022 (a), 1023 (b), 3 × 1023 (c) and 1024 W/cm2 (d), respectively. For comparison, spectra given by present NBW model (black dotted) were also shown.
EILP on the spectrum of electron produced by NBW (red lines). The parameters are: , , GeV, eV (λ = 1 μm) and W/cm2 is 1022 (a), 1023 (b), 3 × 1023 (c) and 1024 W/cm2 (d), respectively. For comparison, spectra given by present NBW model (black dotted) were also shown.Different from the NCS case, the non vanishing emission angle of NBW in one CI introduced by EILP is not fixed. Figure 6 shows the emission angle θ of e− in UILs, it strongly depends on δ. Note that the emission angle of e+ is symmetric with respect to δ = 0.5. The differential probabilities are also present to show the physical significant ranges of δ. When laser intensity is low, e− and e+ are highly possible to be emitted along directions close to k″. The probability of large angle e− and e+ emission grows with laser intensity, and becomes dominant when W/cm2.
Figure 6
Emission angle θ (red solid line for + and black dashed for −) of NBW with respect to instantaneous γ photon momentum k″ in UILs. Parameters are: , , GeV, eV (λ = 1 μm) and W/cm2 is 1022 (a), 1023 (b), 1024 (c) and 1025W/cm2 (d), respectively. The blue dotted lines are normalized differential probabilities.
Emission angle θ (red solid line for + and black dashed for −) of NBW with respect to instantaneous γ photon momentum k″ in UILs. Parameters are: , , GeV, eV (λ = 1 μm) and W/cm2 is 1022 (a), 1023 (b), 1024 (c) and 1025W/cm2 (d), respectively. The blue dotted lines are normalized differential probabilities.To explore signals of EILP on NCS and NBW on future laser facilities, Monte-Carlo simulations of interactions between an electron bunch and UIL pulses are carried out. Eqs (4), (7) and (8) are employed to describe NCSs and NBWs including EILP are described by Eqs (5), (13) and (14). Between emissions, classical equations of motion are applied to describe electron propagation. Other processes such as higher order radiations[30] are neglected for their much weaker effects under the laser conditions concerned.In the simulation, as shown in Fig. 7, a 1.022 GeV quasi mono-energetic electron bunch collides head-on with a tightly focused linearly polarized short laser pulse. The laser field is given by an approximate solution of Maxwell’s equations to the first order of and , where w0 and τ0 are the waist radius and pulse duration[31]. Electromagnetic force between electrons is ignored for it is 7–8 magnitudes weaker than the laser Lorentz force.
Figure 7
Schematics of the simulations.
Schematics of the simulations.Applied parameters are the following: peak intensity of the laser pulse is and 1024 W/cm2, the laser wave length λ = 1 μm, beam waist w0 = 1 μm and duration fs. The electron bunch includes electrons, which uniformly distribute in a R = 1 μm sphere. The mean initial electron energy is 1.022 GeV (), and both the energy dispersion and angular dispersion θ are 0.001. The initial longitudinal distance between the electron bunch and laser pulse centers is 20 μm, the simulation lasts 330 fs, which allows most particles to escape the pulse. The time step is fs, which is fine enough for the results to converge.Simulation results of the W/cm2 case are shown in Fig. 8(a,c). Radiation intensity angular distribution of emitted γ photons including and excluding EILP both concentrate around the initial electron bunch direction, EILP extends the angular divergence (FWHM) from 2° to 4°. Note it is different from that induced by consecutive stochastic emissions[24],
Figure 8
Angular distribution of radiation intensity when W/cm2 (a) and 1024 W/cm2 (b) and energy integrated spectra of emitted photons for (c) and 1024 W/cm2 (d), (d) only includes photons in the backward hemisphere of initial electron bunch. All normalized by the initial electron number N, blue solid lines include EILP and red dashed lines exclude EILP.
Angular distribution of radiation intensity when W/cm2 (a) and 1024 W/cm2 (b) and energy integrated spectra of emitted photons for (c) and 1024 W/cm2 (d), (d) only includes photons in the backward hemisphere of initial electron bunch. All normalized by the initial electron number N, blue solid lines include EILP and red dashed lines exclude EILP.When the laser intensity is increased to 1024 W/cm2, as shown in Fig. 8(b), EILP creates a strong, new peak in emission intensity angular distribution which is very close to the laser forward direction. Note that the mechanism for this peak is different from the stochastic effect peak in[26] and the parameters are quite different ( in our case while in ref.[26]). EILP on the photon spectrum is also significant. Figure 8(d) shows that, the spectrum of photons emitted in the backward hemisphere excluding EILP is bounded below 0.5 GeV while that including EILP extends 1.5 GeV and higher, and the intensity of radiation is also several times higher.In conclusion, EILP on NCS and NBW in one CI in UILs is investigated. Present understanding of these two processes is based on two assumptions, forward emission and spectrum cutoff. These assumptions that exclude EILP are good approximations when the energy scale of involved laser photon is much lower than that of the electron/γ photon.When the total energy of involved laser photons becomes comparable to that of the electron/γ photon, NCS/NBW can have large emission angles and the high energy part of spectrum can surpass the instantaneous electron/γ photon energy. Furthermore, the spectrum of NBW gets an additional lower limit.Simulation results demonstrate that EILP is very important in future SFQED experiments, it dominates single NCS/NBW when laser intensity is close to 1024 W/cm2. Its effects is significant at a much lower intensity of 2 × 1022 W/cm2, which is within the reach of state-of-art lasers.Supplemental Information
Authors: V Yanovsky; V Chvykov; G Kalinchenko; P Rousseau; T Planchon; T Matsuoka; A Maksimchuk; J Nees; G Cheriaux; G Mourou; K Krushelnick Journal: Opt Express Date: 2008-02-04 Impact factor: 3.894
Authors: W P Leemans; A J Gonsalves; H-S Mao; K Nakamura; C Benedetti; C B Schroeder; Cs Tóth; J Daniels; D E Mittelberger; S S Bulanov; J-L Vay; C G R Geddes; E Esarey Journal: Phys Rev Lett Date: 2014-12-08 Impact factor: 9.161
Authors: A Gonoskov; A Bashinov; I Gonoskov; C Harvey; A Ilderton; A Kim; M Marklund; G Mourou; A Sergeev Journal: Phys Rev Lett Date: 2014-07-02 Impact factor: 9.161