Pulse chirping effect on controlling the transverse cavity oscillations in nonlinear bubble regime
Vosoughian H1, 3, 4, †, Riazi Z2, Afarideh H3, Sarri G4
Laser and Optic Institute, Nuclear Science and Technology Research Institute, Tehran, Iran
Physics and Accelerator School, Nuclear Science and Technology Research Institute, Tehran, Iran
Department of Energy Engineering and Physics, Amirkabir University of Technology, P. O. Box 15875-4413, Tehran, Iran
School of Mathematics and Physics, The Queen’s University of Belfast, BT7 1NN Belfast, United Kingdom

 

† Corresponding author. E-mail: h.vosoughian@qub.ac.uk

Abstract

The propagation of an intense laser pulse in an under-dense plasma induces a plasma wake that is suitable for the acceleration of electrons to relativistic energies. For an ultra-intense laser pulse which has a longitudinal size shorter than the plasma wavelength, λ p, instead of a periodic plasma wave, a cavity free from cold plasma electrons, called a bubble, is formed behind the laser pulse. An intense charge separation electric field inside the moving bubble can capture the electrons at the base of the bubble and accelerate them with a narrow energy spread. In the nonlinear bubble regime, due to localized depletion at the front of the pulse during its propagation through the plasma, the phase shift between carrier waves and pulse envelope plays an important role in plasma response. The carrier–envelope phase (CEP) breaks down the symmetric transverse ponderomotive force of the laser pulse that makes the bubble structure unstable. Our studies using a series of two-dimensional (2D) particle-in-cell (PIC) simulations show that the frequency-chirped laser pulses are more effective in controlling the pulse depletion rate and consequently the effect of the CEP in the bubble regime. The results indicate that the utilization of a positively chirped laser pulse leads to an increase in rate of erosion of the leading edge of the pulse that rapidly results in the formation of a steep intensity gradient at the front of the pulse. A more unstable bubble structure, the self-injections in different positions, and high dark current are the results of using a positively chirped laser pulse. For a negatively chirped laser pulse, the pulse depletion process is compensated during the propagation of the pulse in plasma in such a way that results in a more stable bubble shape and therefore, a localized electron bunch is produced during the acceleration process. As a result, by the proper choice of chirping, one can tune the number of self-injected electrons, the size of accelerated bunch and its energy spectrum to the values required for practical applications.

1. Introduction

Since the generation of high power femtosecond laser pulses, numerous investigations have been devoted to the study of the interaction of high-intensity lasers with plasma. One application of the laser–plasma interaction is the generation of the high brightness and the high-energy electron beams in which electrons with energy up to 4.2 GeV, 6-pC charge, and 0.3-mrad divergence have been produced from a 9-cm-long capillary discharge waveguide.[112] The propagation of an intense laser pulse in an under-dense plasma induces a plasma wake that is suitable for the acceleration of the electrons to relativistic energies.[1316] The electrons with sufficient initial kinetic energy can ride the plasma wake, staying in phase with the longitudinal electric field inside the wake and gaining energy. Because of the small structure of the wake field, i.e., its micron-scale longitudinal and transverse size, the bunches of a duration less than one plasma period have to be injected into the laser wake field in order to get high-quality electron bunches. However, the effective acceleration mechanism occurs when the ambient plasma electrons contribute to the injection process (self-injection). The self-injection process is explained by the two physical mechanisms: wave-breaking (the longitudinal self-injection),[17,18] and the bubble regime (transverse self-injection).[1921] The bubble regime is a multidimensional effect, where the symmetric transverse ponderomotive force of the laser expels the electrons radially from the propagation axis and creates a spherical region (bubble) devoid of the plasma electrons behind the laser. This bubble moves with a velocity equal to the group velocity of the laser pulse in the plasma. An intense charge separation electric field inside the moving bubble can capture the electrons at the base of the bubble and accelerate them with a narrow energy spread.The quality of the accelerated electrons is mainly determined by the evolution of the bubble structure. Recently, there have been several theoretical[22,23] and experimental[24] investigations to study the influence of the bubble evolution on the wake field structure and hence on the quality of the electron beam.It has been shown that in the fully relativistic regime in which the normalized vector potential of laser a 0 = eE 0/m e 0 approaches 1 (where E 0 and ω0 are the peak electric field and the frequency of the laser), the interaction of the laser with the under-dense plasma can lead to several nonlinear processes such as laser self-steepening.[25] In laser self-steepening, the peak region of the pulse with high-intensity moves faster than those with low intensity at the pulse head. This leads to optical shock formation and then a nonlinear wake is excited. As a result,the plasma wavelength elongates and uninterrupted self-injection follows.[2628] The laser pulse is depleted due to the excitation of a plasma wave. For a situation where the laser is ultra-intense (a 0 ≫ 1), the pump deletion occurs at a narrow region in front of the pulse corresponding to a few cycles of the laser front.[29] A consequence of the localized pulse depletion is induction of the transverse oscillations in the bubble sheath leading to multi-injection, a large energy spread and collective transverse oscillation of the electron beams.[30] Therefore, it is important to control these instabilities to achieve high-quality electron beams. In this work the interaction of a chirped laser pulse with the under-dense plasma is studied via two-dimensional (2D) particle-in-cell (PIC) simulations. The aim of this work is to propose a new approach to suppress the induced transverse oscillations in the bubble structure at ultra-relativistic regimes (a 0 ≈ 10) using the frequency-chirped laser pulse. Based on this technique, a considerable difference in pulse depletion process is observed for positively and negatively chirped laser pulses. This difference can be used to tune the bubble structure and hence, the quality of the injected bunch in the nonlinear bubble regime.

2. Theory and simulation parameters

The interaction of the short-pulse and the ultra-high intense lasers (I > 1019 W/cm2) with the undersense plasma presents different features compared to moderately high intensities (I < 1018 W/cm2). At moderately high intensity, the nonlinear instabilities like relativistic self-focusing,[31] relativistic self-phase modulation,[32] Raman backward and forward scattering,[33] and envelope self-modulation[34] cause the laser pulse to be absorbed and scattered in a distance of the order of a Rayleigh length. At ultra-high intensities, due to relativistic effects,[3538] the growth rate for the above-mentioned instabilities decreases. This leads to the depletion of only the leading edge of the pulse and the main body of the pulse undergoes little modification. Therefore, the laser pulse eventually evolves to a state in which the leading edge becomes a step function. The phase and the group velocity of the laser in a dispersive medium such as plasma are not equal. Therefore, propagation of a laser pulse through plasma induces a phase shift between the carrier wave and the envelope known as the carrier–envelope phase (CEP). The time variation of the electric field for extremely short laser pulses of few-oscillations depends on the carrier–envelope phase. The localized depletion of a laser pulse at ultra-high intensities causes the varying electric field of the laser at the front of the depleted pulse to act like a few-cycle pulse[39,40] in which the asymmetric transverse ponderomotive force leads to the occurrence of the transverse oscillations in the bubble structure.[41] To stabilize the bubble shape, the depletion of the pulse leading edge must be mitigated in order to avoid the effect of the carrier phase shift on the transverse ponderomotive force. Since the depletion velocity of the pulse, i.e., ,[29] depends on the laser frequency, i.e., ω(ξdepl), it is possible to postpone (speed up) the pulse depletion rate by locating the higher frequency (lower frequency) at the front of the pulse. Here, ξdepl is the position of the localized depletion within the laser in dimensionless delayed time ξ, coordinates and ω p is the plasma frequency. By this technique, we would be able to significantly suppress/enhance the transverse cavity oscillations when the pulse propagates through the plasma. In order to demonstrate the proposed scheme, the 2D PIC simulations have been carried out by EPOCH[42] code to show the effect of pulse chirping on the evolution of the bubble structure. In these simulations, the total simulation box is 200 µm (x) × 70 µm (y) which corresponds to a window with 2000 × 900 cells and 6 × 106 particles. The under-dense hydrogen plasma (T e = 0) is initially located in the region of 30 µm < x < 175 µm, –30 µm < y < 30 µm, with the density of ne = 0.02nc, where nc = 1.1 × 1021 cm–3 is the critical density for the laser pulse with a wave length of λLaser = 1 µm. The electric field of the chirped laser pulse that propagates along the x axis can be expressed as[43]

Here the electric field amplitude of the laser pulse is defined as E 0 which is a constant value. The local frequency of the linearly chirped pulse, the longitudinal laser length, and the laser beam waist are denoted by ω (ξ) = (1 – bξ), σ x and σ y respectively. Here, b > 0 (b < 0) denotes the dimensionless chirp parameter for the positively (negatively) chirped pulse and ξ = [(xct) ω0/c] is the dimensionless delayed time in which c, m e, ω0 and e are the speed of light in vacuum, the plasma electron mass, the local frequency in the center of the pulse, and the absolute charge of the plasma electron, respectively. Also, t and x represent the time and spatial coordinates, respectively. In this work the laser beam waist, the longitudinal laser length, and the normalized incident laser amplitude are set equal to σy = 8 µm, σx = 10 µm, and a0 = eE 0/m e cω0 = 10 respectively.

3. Results and discussion

The bubble structure at the three different instances during the propagation of the laser pulse through the under-dense plasma are shown in Fig. 1. Here the unchirped laser pulse with a normalized amplitude of a 0 = 10, a longitudinal size of σ x = 10 µm, and a spot size of σ y = 8 enters the simulation box from the left boundary.

Fig. 1. (color online) The evolution of THE bubble structure when the un-chirped laser pulse propagates through the under-dense plasma (n e = 0.02n c) at the times of: (a) t = 300 fs, (b) t = 500 fs, and (c) t = 600 fs. The laser pulse with the normalized amplitude of a 0 = 10, the longitudinal size of σ x = 10 µm, and the spot size of σ y = 8 µm enters the simulation box from the left boundary.

As is clear from Fig. 1(a), at the initial stage of the bubble regime, i.e., at t = 300 fs, the symmetric bubble structure is established whereas, at the later times, the instabilities appear in the surrounding sheath of the bubble as seen in Figs. 1(b) and 1(c). When the ultra-intense laser pulse interacts with the plasma, a density spike is produced by the ponderomotive force of the front of the laser pulse and subsequently a large plasma wake is excited behind the laser. The longitudinal electric field, the plasma density at the front of the pulse, and the transverse electric field of laser are plotted in Figs. 2(a), 2(b), and Fig. 3, respectively, for t = 300 fs, t = 500 fs, and t = 600 fs.

As shown in Figs. 2(a) and 2(b), the amplitude of wake and also the density spike continuously enhances during the time. As a result, the laser energy absorption by the electrons in the spike or pulse depletion will increase too. Viewed from a photon description, the photons localized within the density spike lose energy as they excite the wake and therefore their wavelength increases leading to frequency red-shifting at the front of the pulse (Fig. 3).

Fig. 2. (color online) (a) The excited wake field behind the un-chirped pulse, (b) The excited density spike at the pulse head of the un-chirped laser pulse when it propagates through the under-dense plasma (n e = 0.02n c) at the times of t = 300 fs, t = 500 fs, and t = 600 fs.
Fig. 3. (color online) The evolution of the laser’s electric field of the unchirped laser pulse during its propagation through the under-dense plasma at the times of t = 300 fs, t = 500 fs, and t = 600 fs.

In this situation the front of the pulse is approximately converted to a one-cycle profile. Once the leading edge of the pulse erodes, the bubble structure becomes susceptible to transverse instabilities. In fact, considering the plasma as a dispersive medium, due to the induced phase shift between the carrier waves and the pulse envelope, the role of the time average of the laser field at the front of the depleted pulse becomes dominant on the intensity envelope of the pulse in laser–plasma interactions. Therefore, changing the field direction at the pulse front, as it is depicted by the red and blue arrows in Fig. 3, causes the electrons to be asymmetrically expelled from the laser propagation axis leading to transverse oscillations in the bubble structure. In order to control the evolution of the bubble structure we benefit from the negatively chirped laser pulse to reduce the pulse depletion rate at the pulse’s leading edge. The interaction of a negatively chirped laser pulse with an under-dense plasma (n e = 0.02n c) is compared with a positively chirped one in Fig. 4.

Fig. 4. (color online) The evolution of bubble structure when the chirped laser pulse propagates through the under-dense plasma (n e = 0.02n c). The evolution of the bubble in the case of the negatively chirped pulse with the chirped parameter of b = –0.08 at the times of: (a1) t = 300 fs, (a2) t = 500 fs, and (a3) t = 600 fs. The evolution of the bubble in the case of the positively chirped pulse with the chirped parameter of b = 0.08 at the times of: (b1) t = 300 fs, (b2) t = 500 fs, and (b3) t = 600 fs. The other laser parameters are the same as the un-chirped case.

As can be seen for the negatively chirped pulse the symmetric structure of the moving bubble remains stable with time (Figs. 4(a1)4(a3)), whereas as Figures 4(b1)4(b3) show the cavity oscillations in the surrounding sheath of the bubble for the positively chirped pulse occur even sooner than the bubble for the un-chirped pulse. The density spike at the pulse head and the longitudinal electric field for both the negatively and the positively chirped pulse are depicted in Fig. 5 and Fig. 6, respectively for t = 300 fs, t = 500 fs, and t = 600 fs.

Fig. 5. (color online) The electron density at the pulse head of: (a) the negatively chirped pulse, (b) the positively chirped pulse, when propagating through the under-dense plasma (n e = 0.02n c) at the times of t = 300 fs, t = 500 fs, and t = 600 fs.

In fact, with the higher frequencies located in front of the pulse, the plasma transparency increases at the leading part of the pulse that results in the excitation of a small density spike (Fig. 5(a)) and also a weaker wake field compared with the un-chirped pulse (Fig. 6(a)), whereas the low frequency cycles at the front of the positively chirped pulse initiate the largest density spike (Fig. 5(b)) and the strongest wake field behind the laser (Fig. 6(b)).

Fig. 6. (color online) The excited wake field behind the chirped laser pulse propagating through the under-dense plasma: (a) the negatively chirped laser pulse, (b) the positively chirped laser pulse at the times of t = 300 fs, t = 500 fs, and t = 600 fs.
Fig. 7. (color online) The evolution of the laser’s electric field: (a) The negatively chirped pulse, (b) The positively chirped pulse, during its propagation through the under-dense plasma at the times of t = 300 fs, t = 500 fs, and t = 600 fs.

As a result, the pulse depletion rate at the leading edge of the negatively chirped (positively chirped) pulse decreases (increases) and thus, the carrier phase shift becomes negligible (significant) in the leading edge of the pulse as depicted in Fig. 7(a) (Fig. 7(b)). Therefore, the ponderomotive force of laser keeps (breaks down) its symmetric properties in transverse direction during the pulse propagation.

The longitudinal electric field together with the plasma density are shown in Figs. 8(a)8(c) for negatively chirped pulse at times t = 410, 510, and 570 fs. The electrostatic fields are denoted by the blue curves and the density distributions of injected electrons are specified by the black curves. Also, the energy spectrum of the injected electrons at the corresponding time is plotted in Fig. 8(d). Figures 9(a)9(d) are the same as Figs. 8(a)8(d), but for positively chirped pulse. In these figures the electron-accelerating phase buckets (where E x < 0) and the electron-deceleration buckets (where E x > 0) of the wake wave are crosshatched by the green lines and the red lines, respectively. It is seen from Fig. 8(a) that the electrons are trapped in plasma wake field, as pointed out by the blue arrow. These trapped electrons stay in the accelerating phase of the field as can be seen in Figs. 8(b) and 8(c). This leads to continuous acceleration of electrons (Fig. 8(d)) to higher energies.

Fig. 8. (color online) The superposition of one-dimensional (1D) diagrams of the wake field structure and plasma density during the propagation of the negatively chirped pulse in plasma at the times of: (a) t = 410 fs, (b) t = 510 fs, and (c) t = 570 fs. The blue curves show the electrostatic field of wake field and the black curves specify the density distribution of injected electrons. Also the green region determines the acceleration phase of the wake field. (d) The energy spectrum of injected electrons at the corresponding times of: t = 410 fs, t = 510 fs, and t = 570 fs.
Fig. 9. (color online) The superposition of 1D diagrams of the wake field structure and the plasma density during the propagation of the positively chirped pulse in plasma at the times of: (a) t = 400 fs, (b) t = 510 fs, and (c) t = 570 fs. The blue curves show the electrostatic field of the wake field and the black curves specify the density distribution of the injected electrons. Also the green and red regions determine the acceleration phase and deceleration phase of the wake field, respectively. (d) The energy spectrum of the injected electrons at the corresponding times of: t = 410 fs, t = 510 fs, and t = 570 fs.

However, as can be clearly seen in Figs. 9(a)9(d), for the positively chirped pulse, due to the evolution of the bubble base and the loading of a large number of electrons in the moving bubble, the accelerating region of the wake field turns into a deceleration region frequently and the injected electrons sometimes enter the deceleration region and lose energy.

By comparing the energy spectrum of the electrons for the negatively and the positively chirped pulse (Figs. 8(d) and 9(d)), we can see that higher-energy electrons are obtained from the former and a higher number of electrons is obtained from the latter due to the multiple self-injection.

4. Conclusion

In the nonlinear bubble regime, the transverse oscillations of the cavity structure as a result of the depletion of the pulse leading edge can be controlled by applying the chirped laser pulses. The simulation results show that by locating the higher frequencies in front of the pulse, the pulse depletion process is postponed, thus the cavity structure stays stable for a longer time and a single electron bunch with a higher energy can be generated. In the case in which the lower frequencies are located at the front of the pulse, the pulse depletion speeds up, thus the transverse instability of the bubble appears in an earlier stage of laser plasma interaction. The unstable bubble results in the accelerating phase of the wake field being turned into the decelerating phase, frequently. As a consequence the energy loss of the injected electrons causes the generation of the electron bunches with lower mean energy. Therefore, the utilizations of the chirped laser pulse provides a method to control the quality of the generated bunch in the nonlinear bubble regime.

Reference
[1] Geddes C G R Toth Cs van Tilborg J Esarey E Schroeder C B Bruhwiler D Nieter C Cary J Leemans W P 2004 Nature 431 538
[2] Faure J Rechatin C Norlin A Lifschitz A Glinec Y Malka V 2006 Nature 444 737
[3] Blumenfeld I Clayton C E Decker F J Hogan M J Huang C Ischebeck R Iverson R Joshi C Katsouleas T Kirby N Lu W Marsh K A Mori W B Muggli P Oz E Siemann R H Walz D Zhou M 2007 Nature 445 741
[4] Lu W Tzoufras M Joshi C Tsung F S Mori W B Vieira J Fonseca R A Silva L O 2007 Phys. Rev. STAB 10 061301
[5] Zhao X Y Xie B S Wu H C Zhang S Hong X R Aimidula A 2012 Phys. Plasmas 19 033108
[6] Wang J W Yu W Yu M Y Xu H Ju J J Luan S X Murakami M Zepf M Rykovanov S 2016 Phys. Rev. Accel. Beams. 19 021301
[7] Esarey E Schroeder C B Leemans W P 2009 Rev. Mod. Phys. 81 1229
[8] Su H Y Huang Y S Wang N Y Tang X Z Lu W 2014 Chin. Phys. Lett. 31 075202
[9] Yin C L Wang W M Liao G Q Li M C Li Y T Zhang J 2015 Acta Phys. Sin. 64 144102 (in Chinese)
[10] Wu F J Zhou W M Shan L Q Li F Liu D X Zhang Z M Li B Y Bi B Wu B Wang W W Zhang F Gu Y Q Zhang B H 2014 Acta Phys. Sin. 63 094101 (in Chinese)
[11] Li W T Wang W T Liu J S Wang C Zhang Z J Qi R Yu C H Li R X Xu Z Z 2015 Chin. Phys. B 24 015205
[12] Leemans W P Gonsalves A J Mao H S Nakamura K Benedetti C Schroeder C B Mittelberger D E Bulanov S S Vay J L Geddes C G R Esarey E 2014 Phys. Rev. Lett. 113 245002
[13] Sprangle P Esarey E Krall J Joyce G 1992 Phys. Rev. Lett. 69 2200
[14] Esarey E Sprangle P Krall J Ting A Joyce G 1993 Phys. Fluids B 5 2690
[15] GorbunovL M Kirsanov V I 1987 Sov. Phys. JETP 66 290
[16] Esarey E Sprangle P Krall J Ting A 1996 IEEE Trans. Plasma Sci. 24 252
[17] Akhiezer A I Polovin R V 1956 Sov. Phys. JETP 3 696
[18] Dawson J M 1959 Phys. Rev. 113 383
[19] Pukhov A Meyer-terVehn J 2002 Appl. Phys. B: Lasers Opt. 74 355
[20] Zhidkov A Koga J Hosokai T Kinoshita K Uesaka M 2004 Phys. Plasmas 11 5379
[21] Lu W Huang C Zhou M Mori W B Katsouleas T 2006 Phys. Rev. Lett. 96 165002
[22] Ma Y Chen M Hafz N A M Li D Z Huang K Yan W C Dunn J Sheng Z M Zhang J 2014 Appl. Phys. Lett. 105 161110
[23] Lee S Lee T H Gupta D N Uhm H S Suk H 2015 Plasma Phys. Control. Fusion 57 075002
[24] Hafz N A M Lee S K Jeong T M Lee J 2011 Nucl. Instru. Methods Phys. Res. A 637 S51
[25] DeMartini F Gustafson T K Kelley P 1967 Phys. Rev. 164 312
[26] Kalmykov S Y Beck A Yi S A Khudik V N Downer M C Lefebvre E Shadwick B A Umstadter D P 2011 Phys. Plasmas 18 056704
[27] Faure J Glinec Y Santos J J Ewald F Rousseau J P Kiselev S Pukhov A Hosokai T Malka V 2005 Phys. Rev. Lett. 95 205003
[28] Kalmykov S Y Beck A Davoine X Lefebvre E Shadwick B A 2012 New J. Phys. 14 033025
[29] Decker C D Mori W B Tzeng K C Katsouleas T 1996 Phys. Plasmas 3 2047
[30] Popp A Vieira J Osterhoff J Major Zs Horlein R Fuchs M Weingartner R Rowlands-Rees T R Marti M Fonseca R A Martins S F Silva L O Hooker S M Krausz F Gruner F Karsch S 2010 Phys. Rev. Lett. 105 215001
[31] Sprangle P Tang C M Esarey E 1987 IEEE Trans. Plasma Sci. PS-15 145
[32] Max C E Arons J Langdon A B 1974 Phys. Rev. Lett. 33 209
[33] Mori W B Decker C D Katsouleas T Hinkel D E 1994 Phys. Rev. Lett. 72 1482
[34] Esarey E Krall J Sprangle P 1994 Phys. Rev. Lett. 72 2887
[35] Bulanov S V Inovenkov I N KirsanovV I Naumova N M Sakharov A S 1992 Phys. Fluids B 4 1935
[36] Sprangle P Esarey E 1991 Phys. Rev. Lett. 67 2021
[37] Esarey E Sprangle P 1992 Phys. Rev. A 45 5872
[38] Sakharaov A S Kirsanov V I 1994 Phys. Rev. E 49 3274
[39] Paulus G G Grasbon F Walther H Villoresi P Nisoli M Stagira S Priori E De Silvestri S 2001 Nature 414 182
[40] Cormier E Lambropoulos P 1998 Eur. Phys. J. D 2 15
[41] Nerush E N Kostyukov I Yu 2009 Phys. Rev. Lett. 103 035001
[42] See http://ccpforge.cse.rl.ac.uk/gf/project/epoch/ for “EPOCH: Extendable PIC Open Collaboration” (2011)
[43] Afhami S Eslami E 2014 Phys. Plasmas 21 063108