|Home | About | Journals | Submit | Contact Us | Français|
Signatures of stochastic effects in the radiation of a relativistic electron beam interacting with a counterpropagating superstrong short focused laser pulse are investigated in a quantum regime when the electron’s radiation dominates its dynamics. We consider the electron-laser interaction at near-reflection conditions when pronounced high-energy gamma-ray bursts arise in the backward-emission direction with respect to the initial motion of the electrons. The quantum stochastic nature of the gamma-photon emission is exhibited in the angular distributions of the radiation and explained in an intuitive picture. Although, the visibility of the stochasticity signatures depends on the laser and electron beam parameters, the signatures are of a qualitative nature and robust. The stochasticity, a fundamental quantum property of photon emission, should thus be measurable rather straightforwardly with laser technology available in near future.
The next generation petawatt laser systems1, 2 will open a door not only to novel regimes of laser-matter interaction3, 4, but also to new perspectives for the investigation of fundamental problems5–9. In ultrastrong laser fields the quantum properties of electron radiation, the discrete and probabilistic character of photon emission, can be conspicuous. While the discreteness of the radiation photon energy is known to be observed straightforwardly, e.g., in the Compton scattering as a shift of the emission frequency10–12, the signatures of the stochastic character of photon emission are more subtle and elaborate for observation. The latter has impact on the radiation back-action to the electron dynamics and should be more apparent in the so-called radiation dominated regime (RDR) of interaction8, 13, 14, when multiple emission of photons by an electron becomes probable. One of the conceptual consequences of stochasticity effects (SE) in photon emissions, i.e., the probabilistic nature of photon emission, is the broadening of the energy spread of an electron beam in a plane laser field15, 16, while similar effect can cause electron stochastic heating in a standing laser field17. However, in an experiment in a focused laser beam, competing effects may arise, e.g., an additional energy spreading of the electron beam due to the difference of radiative losses of electrons in the electron-beam cross section. Another SE signature is the so-called electron straggling effect during radiation in strong fields18–21, when the electrons propagate a long distance without radiation due to SE, resulting in the increase of the yield of high-energy photons. However, we will show below that the straggling effect is rather weak in the considered regime. Different signatures of quantum radiation reaction have been discussed22–35 which include all quantum effects, such as photon recoil, stochasticity, and interferences35. However, they do not isolate information on the specific role of stochasticity. Unfortunately, so far unequivocal SE have not been observed in an experiment, which would also be important as a clear test for the underlying theories. We underline that in this paper SE refers to the quantum probabilistic nature of photon emission, but not to the chaotic classical dynamics as, e.g., in ref. 36.
The quantum effects in strong laser fields are determined by the invariant parameter 37, 38, where F μν is the field tensor, p ν=(ε, p) the incoming electron 4-momentum, and e and m are the electron charge and mass, respectively (Planck units ħ=c=1 are used throughout). The RDR, when the radiation losses during a laser period are comparable with the electron initial energy, is characterized by the parameter 8, where α is the fine structure constant, ξ|e|E 0/(mω 0) the invariant laser field parameter, while E 0 and ω 0 are the laser field amplitude and frequency, respectively. In the quantum RDR the stochastic nature of photon emission is a fundamental quantum property and has to be taken into account during the high-energy photon radiation39–42, that significantly affects the electron dynamics and the high-energy photon emission.
In this paper, we investigate signatures of the stochastic nature of photon emission in the nonlinear Compton scattering in the quantum RDR during the interaction of a superstrong short focused laser pulse with a counterpropagating relativistic electron beam. We consider the interaction in the electron near-reflection regime14, where the front electrons are reflected and emit gamma-rays in the near-backward direction29: due to the combined effects of the laser focusing and radiation reaction, the front electrons of the electron beam are reflected and emit ultrashort gamma-rays in the near-backward direction with respect to the initial electron motion29. In the considered case with , the straggling effect is rather weak, however, the electron near-reflection regime offers a possibility to observe the SE in the angular distribution of the radiation. In fact, the photon emissions in different laser cycles are essentially modified due to the SE in the RDR, and the latter is mapped into the broad backward-emission angles when the electron is at the near-reflection condition. We calculate angle-resolved radiation intensity and photon numbers and show that, due to the stochastic nature of photon emission, the radiation’s angular distribution (RAD) has a single prominent peak in the backward direction, which is broad and easily observable in an experiment. In contrast, when the SE are ignored, the backward radiation yields an angular distribution with several peaks. Furthermore, we investigate the influences of the laser and electron-beam parameters (the laser focal radius, the laser pulse duration and the electron initial energy) on the visibility of those SE signatures, to optimize parameters for a future experimental setup.
The considered quantum RDR requires the invariant parameters and , while the electron near-reflection regime does γ ∼ ξ/2, where γ is the Lorenz factor of the electron, and β the electron velocity scaled by the light speed in vacuum. The two conditions above demand γ ∼ ξ ∼ 103, i.e., an electron beam of GeV energies and laser intensities of 1023–1024W/cm2 anticipated in next generation facilities1, 2.
The calculation of the radiation is based on Monte-Carlo simulations employing QED theory for the electron radiation and classical equations of motion for the propagation of electrons between photon emissions40–42. In superstrong laser fields ξ ≫ 1, the coherence length of the photon emission is much smaller than the laser wavelength and the typical size of the electron trajectory38, 43. Then, the photon emission probability is determined by the local electron trajectory, consequently, by the local value of the parameter χ 44. The radiation in the quantum regime ignoring the SE are calculated by employing the Sokolov equation22–24, when the radiation is emitted continuously along the electron trajectory and modifies accordingly the classical equations of motion. In this work we consider the regime when the quantum invariant field parameter χ is not small, . Therefore, the quantum recoil in photon emission is significant even in the case when the probabilistic character (stochasticity) of photon emission is neglected. Therefore radiation intensity and the radiation reaction should be calculated using quantum corrections. The Sokolov equation is a phenomenologically derived equation of motion for an electron in the ξ ≫ 1 limit, which is based on the energy-momentum conservation within the system of the electron and emitted photons at each formation length of radiation [however, it is argued in ref. 45 that this equation cannot be derived from QED by taking the classical limit]. The Sokolov equation adapts the electron energy-momentum at each step according to the radiated energy-momentum including quantum corrections. It provides a natural way to treat the electron dynamics in the external field classically but to take into account in the radiation-reaction quantum-recoil corrections (but not stochasticity). In contrast, the conceptually more systematic Lorentz-Abraham-Dirac or Landau-Lifshitz equations account for the radiation reaction only classically (i.e. without recoil). The simulation methods including and excluding the SE are described in detail in the section “Methods”.
We employ a linearly polarized tightly focused laser pulse with a Gaussian temporal profile, which propagates along +z-direction and polarizes in x-direction. The spatial distribution of the electromagnetic fields takes into account up to the ε 3-order of the nonparaxial solution46, 47, where ε=w 0/z r, w 0 is the laser focal radius, the Rayleigh length with the laser wave vector k 0=2π/λ 0, and λ 0 the laser wavelength. The expressions of the laser fields are represented in the section “Methods”.
A typical angular distribution of radiation which carries the signature of the stochastic nature of photon emission, is illustrated in Fig. 1; ϕ=0° and ±180° correspond to the positive and negative directions of the laser polarization, respectively. The peak intensity of the 6-cycle (FWHM) laser pulse is I≈4.9×1023W/cm2 (ξ=600), λ 0=1μm, and w 0=2μm. The electron beam, with radius w e=λ 0, length L e=6λ 0, and density n e≈1015cm−3, initially counterpropagates with the laser pulse, i.e. . The initial mean kinetic energy of the electron beam is ε 0=180MeV (γ 0≈353, the maximum value of χ during interaction ), and the energy and angular spread are Δε/ε 0=Δθ=0.02. The electron-beam parameters are typical for current laser-plasma acceleration setups5–7.
Figures 1(a,b) demonstrate RAD including and excluding SE, respectively. The radiation is most significant along the strongest field component E x of the linearly polarized laser field. However, in a tightly focused laser beam, other components of the electric field, E y and E z, are not negligible and play a significant role in the electron dynamics and the photon emission. Consequently, the electrons radiate continuously over the azimuthal angle ϕ, with a Gaussian radiation distribution with respect to ϕ corresponding to the laser transverse profile. Moreover, the radiation sweeps from the polar angle down to θ≈11°. The electrons initially counterpropagate with the laser pulse, emit forwards, and create a high-intensity spectral region around θ≈180°. Due to radiative losses, the electron energy decreases, facilitating the reflection condition γ≈ξ/2 in the strong-field region and inducing electron reflection, during which the radiation flash sweeps down to θ≈11°. The emission angle after the reflection with respect to the electron average motion δθ ∼ ξ/γ is determined by the values of ξ and γ in situ after the reflection. In the region of θ≈0°, the radiation intensity is vanishing, because of the very low χ value for the co-propagating electron: the radiation energy ε R χ γξ(1−βcosθ)≈0 at θ≈0°.
The radiation distributions with respect to the polar angle θ with and without SE are essentially different in Fig. 1. While with SE the RAD is smoothly peaking at θ=180° and at small angles, without SE it shows a band structure corresponding to the radiation emerging from different laser cycles. A similar behaviour is seen in RAD for photon numbers (See the Supplemental Materials for the details, which includes the impacts of different parameters on the radiation angular distributions).
A quantitative comparison between RADs including and excluding SE is represented in Fig. 2. We focus on the strongest radiation domain along the polarization plane in the region of −15°≤ϕ≤+15°, analysing the radiation energy d/[dθsin(θ)]=dϕ dε R/dΩ and the photon number d/[dθsin(θ)]=dϕ dN R/dΩ in this domain. The stochastic nature of photon emission is clearly discernible in RAD: a single broad high-intensity gamma-photon peak is formed in the near-reflection direction when SE is included, while in the case without SE multiple radiation peaks emerge corresponding to the emission from different laser cycles. Therefore, the main detectable difference of RADs with and without SE is the peak number of the radiation: the previous has only one radiation peak, and the latter has several peaks corresponding to the laser-cycle structure.
We proceed discussing the role of SE in shaping RAD. The dynamics and radiation of a sample electron is analysed in Fig. 3 (an electron at the beam center, z=L e/2 and x=y=0, is considered). We follow the mean longitudinal momentum and the mean motion direction of the electron during interaction with the laser pulse in dependence on the laser phase , see Fig. 3(a,b), when SE are included or excluded, respectively. To reproduce SE we repeat the simulation 200 times for the same initial conditions of the electron. The single-electron radiation angle is (). The point corresponds to the electron reflection, when the emission angle sharply changes from the forward into the backward direction.
In the case without SE, the single-electron radiation angle is well defined at each moment during interaction. In contrast, in the case with SE one observes spreading of the average longitudinal momenta and the emission angles at a certain laser phase , see Fig. 3(a,b), which stems from the probabilistic character of the photon-emission process. In each possible trajectory, the electron emits photons at different moments and photons of different energies. Consequently, the electron has different energy and emission angle at the same laser phase (same ξ).
The dynamics of the photon emission in the SE case described above is illustrated in Fig. 3(c,d) for two sample trajectories, having a small and a large energy after the reflection in Fig. 3(a). Although both of the trajectories begin with the same initial condition, in the first case the number of photon emissions before the reflection () is larger. In the second case the electron “straggles” (does not emit large number of photons) most of the time before the reflection, only emitting a large photon near the reflection (See the Supplemental Materials for the details, which includes the impacts of different parameters on the radiation angular distributions). Because of the latter, the parameter χ after the reflection is larger for the first trajectory with respect to the second, yielding significantly more radiative loss for the first trajectory than for the second one, and correspondingly, smaller energy (larger emissison angle) for the first trajectory than for the second one after the reflection (See the Supplemental Materials for the details, which includes the impacts of different parameters on the radiation angular distributions). Thus, different probabilistic dynamics of the photon emission, i.e. SE, induces spreading of the photon emission angle at each laser phase.
The angle-resolved radiation intensity and radiation energy are shown in Fig. 3(e,f), respectively, for the case without SE. In each cycle strongest radiation arises near the peaks of the cycles at a certain angle. Between adjacent radiation peaks, there is a gap in emission angle corresponding to the weak-field part of the laser cycle. Therefore, the RAD reveals the laser-cycle structure when SE are neglected.
The RAD with SE for a single electron is shown in Fig. 3(g,h). Since with SE the radiation angle in each laser cycle has a very broad spread shown in Fig. 3(b,g), the gaps in emission angle between adjacent radiation peaks of each cycle are filled out. Consequently, the radiation intensity in this case shows a single gamma-radiation peak corresponding to the peak of the laser pulse. Note that the discussed qualitative features of the RAD for a sample electron do not depend on the initial position of the electron in the electron beam. The variation of the initial position introduces only a slight modification of a quantitative character (See the Supplemental Materials for the details, which includes the impacts of different parameters on the radiation angular distributions).
The electron straggling effect in our setup is estimated in Fig. 2(c,d) for the emitted photon number integrated over the whole 2π-region of the azimuthal angle. The straggling effect is observed as an enhancement of the emission of high-energy photons when the SE are included, see Fig. 2(d) for the emitted photon number above 50MeV. However, the electron straggling effect in the considered regime is insignificant and not relevant for the observation of SE.
We investigate the optimization of the laser (focal radius, pulse duration) and electron (energy) parameters for the best identification of SE. The influence of the laser focusing effect on RAD is studied in Fig. 4. The case of w 0=2λ 0 is optimal, being the same as in Figs 1 and and2,2, when the SE are exhibited by a broad smooth peak of the gamma radiation in backward direction. For w 0=λ 0=w e, the near-reflection radiation is rather weak, and there is no laser-cycle structure in the “without SE” case. For w 0=3λ 0, the main (highest) gamma-photon peaks in the region of are much stronger apparently than those for w 0=2λ 0, since all electrons are much closer to the laser-intensity center, which also weakens the shift of the main peaks in polar angle when the SE are excluded. Consequently, higher angular resolution is required to distinguish the “without SE” case from the one with SE. In the extreme case of a plane-wave laser pulse the structure of multiple gamma-photon peaks in the “without SE” case cannot be observed any more. Therefore, a tightly focused laser pulse with a focal radius of roughly 2 times of the electron-beam radius applies well for the observation of SE.
The dependence of RAD on the laser-pulse length is shown in Fig. 5. As the laser-pulse length increases from 4λ 0 to 10λ 0, the main (highest) radiation peak in the region of gradually decreases, and the plateaus between the peaks gradually increase. In principle, for a longer laser pulse there are more peaks corresponding to the laser cycles in the “without SE” case, which is helpful to distinguish SE. However, for a longer laser pulse, e.g., L L=10λ 0, the laser-field gradient is smaller than for shorter laser pulses, consequently, the ratio between the peaks and the plateaus becomes smaller, and the visibility of the signatures becomes worse. Then, a short laser pulse of L L≈6λ 0 is optimal: it has a more-cycle structure than for ultrashort pulses, and the visibility of the laser-cycle structure is better than longer pulses.
The role of the initial kinetic energy of the electron beam on the radiation angular distribution is analysed in Fig. 6. As the electron kinetic energy ε 0 increases from 120MeV to 240MeV, the radiation energy ε R χ ξγ increases as well since γ ∼ ε0 and ξ being constant is assumed. The peaks in RAD shift right to a larger polar angle since the emission angle after reflection satisfies . Therefore, larger initial kinetic energies are favourable for generation of more clear SE signatures, but one has to keep in mind that for the discussed SE features the reflection condition should be held in advance.
Finally, the role of the colliding angle between the laser pulse and the electron beam is investigated, as shown in Fig. 7. The figure shows that the variation of the colliding angle of 5° does not change the qualitative SE signatures. This is the strength of the qualitative SE signatures that they are robust with respect to the variation of the process parameters.
We have revealed signatures of the stochastic nature of gamma-photon emission during nonlinear Compton scattering of a superstrong short focused laser pulse by a counterpropagating electron beam in the quantum radiation-dominated regime. The signatures are manifested in the qualitative features of the angular distribution of the radiation in the near-backward direction with respect to the initial electron motion, which arises when the electron energy is at the reflection condition. When the stochasticity effects are included, a single broad gamma-photon peak is formed in the near-reflection radiation angular distribution, while several gamma-photon peaks would arise if there were no stochasticity in the photon emission. The signatures are enhanced with tightly focused and short laser pulses. They are robust with respect to variation of the laser and electron parameters and can be observed in near-future laser facilities, such as ELI and XCELS.
In this work, we employ a linearly polarized focused Gaussian laser pulse propagating along +z direction and polarized in x direction. The expressions of the electromagnetic fields are presented in the following refs 46 and 47:
E 0 is the amplitude of the laser fields with normalization factor F n=i to keep at the focus, yielding the scaled coordinates
where is the Rayleigh length, , η=ω 0 t−k 0 z, and ϕ CEP the carrier-envelop phase. Note that ϕ CEP=0 is employed throughout.
where is the normalized emitted photon energy, k 0, k and p are the four-vectors of the driving laser photon, the emitted photon and the electron, respectively, and ω r=2ω/3(1−χω). The applicability of the constant crossed-field approximation for the emission probability in Eq. (9) related to the laser intensity has been estimated and justified (See the Supplemental Materials for the details, which includes the impacts of different parameters on the radiation angular distributions). The photon emission of electrons is considered to be a Monte-Carlo stochastic process40–42. During the electron-laser interaction, for each propagation coherent length Δη, the photon emission will take place if the condition (dW/dη)Δη≥N r is fulfilled, where N r is a uniformly distributed random number in . Herein, the coherent length Δη is inversely proportional to the invariant laser field parameter ξ, i.e., . However, to keep the total photon emission energy consistent, i.e., to exclude numerical error of the simulation of photon emission, we choose . The photon emission probability
where ω min is the driving laser photon energy, and ω max equals the electron instantaneously kinetic energy. In addition, the emitted photon energy ω R is determined by the relation:
where, is an another independent uniformly distributed random number in . Between the photon emissions, the electron dynamics in the laser field is governed by classical equations of motion:
Given the smallness of the emission angle ~1/γ for an ultrarelativistic electron, the photon emission is assumed to be along the electron velocity8, 38. The photon emission induces the electron momentum change p f≈(1−ω R/cp)p i, where p i,f are the electron momentum before and after the emission, respectively. In this simulation, the interference effects between emissions in adjacent coherent lengths are negligible since the laser fields employed are superstrong, i.e., ξ ≫ 1. Therefore, the photon emissions happening in each coherent length are independent of each other. However, as the interference effects are significant, i.e. , the photon emission probabilities in adjacent coherent lengths are summed up to compare with the random number N r until the photon emission happens.
As the stochastic effects are not taken into account, the rate of the electron radiation loss ℐ can be considered as
where is the proper time,
the four-vector of the coordinate, τ c2e 2/(3m), and ℐc = 2αω2ξ2 the classical radiation loss rate. The radiation energy and photon number at each step from the initial phase η i to the final phase η f read:
Then, the photon number of the radiation energy ε R above an arbitrary energy ε arb at each step yields:
where ω arb=ε arb.
To prove explicitly the convergence of our Monte-Carlo simulations, we re-simulate the radiation of Fig. 2 via employing a smaller step size, as shown in Fig. 8. The blue-solid curves present the same radiation as in Fig. 2 with the simulation step number N step=2.7×104, and the red-dash curves show the radiation with a smaller simulation step size, i.e., larger step number of N step=5.4×104. The simulations clearly show that the behaviour of the photon emission angular distribution does not change with decreasing the simulation step size, i.e., increasing the simulation step number, which indicates the convergence of our results. Besides, the simulations based on the “Sokolov” equation have been compared with those based on a modified “Landau-Lifshitz” equation (rescaled radiation but no stochasticity effects)48, and they demonstrate good coincidences (See the Supplemental Materials for the details, which includes the impacts of different parameters on the radiation angular distributions).
J.-X.L., K.Z.H. and C.H.K. conceived the project, J.-X.L. and Y.-Y.C. developed the simulation model and carried out the calculations. All contributed to the development of ideas, discussion of the results, and preparation of the manuscript.
The authors declare that they have no competing interests.
Electronic supplementary material
Supplementary information accompanies this paper at 10.1038/s41598-017-11871-0.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Jian-Xing Li, Email: ed.gpm.dh-ipm@iL.gniX-naiJ.
Karen Z. Hatsagortsyan, Email: email@example.com.