Preprint
Article

Attosecond pulses from Ionization Injection Wake Field Accelerators

Altmetrics

Downloads

125

Views

54

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

13 September 2023

Posted:

14 September 2023

You are already at the latest version

Alerts
Abstract
High-quality ionization injection methods for Wake Field Acceleration driven by lasers or charged beams (LWFA/PWFA) can be optimised so as to generate high-brightness electron beams with tuneable duration in the attosecond range. We present a model of the minimum bunch duration obtainable with low-emittance ionization injection schemes, by spotting the roles of the ionization pulse duration, of the wake field longitudinal shape and of the delay of the ionization pulse position with respect to the node of the accelerating field. The model is tested for the resonant multi-pulse ionization injection (ReMPI) scheme, showing that bunches having length of about 300as can be obtained with a ionization pulse having duration of 30fs FWHM
Keywords: 
Subject: Physical Sciences  -   Other

1. Introduction

The generation of relativistic electron bunches with duration in the attosecond range can lead to pump/probe beams which can be fruitfully employed to unveil ultrafast dynamics [1]. In the context of plasma wake field acceleration either driven by laser pulses (LWFA) [2] or particle beams (PWFA) [3] several methods have been proposed to specifically generate electron beams with duration below the femtosecond scale, from the pioneering work about beam compression of beams externally injected ahead the driver laser pulse [4,5,6], dense attosecond beams with up-ramp density transitions [7], attosecond beams via density modulations [8], attosecond trains obtained by betatron quivering modulations [9,10], few-cycle TW pulses driven electron beams [11,12], attosecond trains via ionization injection [13] and high-brightness electron beams through ionization injection in hybrid LWFA/PWFA schemes [14,15]. As the disentanglement of the electron beam parameters including length, charge, average energy, energy spread and emittance is of paramount importance for the feasibility of the pump/probe attosecond source, a flexible injection/acceleration scheme should be preferred. The two-color ionization injection [16] and the Resonant Multi-Pulse Ionization injection (ReMPI) for LWFA [17], or their equivalent form for the PWFA a.k.a. the trojan-horse scheme result be extremely flexible yet capable of generating high-brightness electron beams [15]. All these schemes use a driver to excite a large amplitude plasma wave and a short wavelength ionization pulse to extract electrons from a dopant. The driver can be a single long-wavelength laser pulse (two-color), a train of resonantly delayed pulses (ReMPI) or a charged beam (trojan-horse). In any of these schemes, the electrons extracted by the low normalized amplitude amplitude ionization pulse ( a 0 , i = e A 0 , i / m e c 2 8 . 5 · 10 10 I λ i 2 1 , here I and λ i are the ionization pulse peak intensity and wavelength in W / c m 2 , and μ m , respectively) do quiver in the laser field until they slip back it out with a residual transverse momentum which will constitute the major contribution for the final beam emittance. As the electrons are accelerated and focused by the wakefield, they are eventually trapped in the bucket and further accelerated. During the slippage in the back of the bucket, the electron beam is compressed in both the longitudinal and radial directions and can reach longitudinal sizes of tens of n m , thus generating electron bunches that can reach attosecond scale duration. As the wakefield driver should not ionize the dopant, it’s maximum electric field should be well below the threshold for tunnel ionization [18] of the selected ionization process (e.g. K r 8 + 9 + , A r 8 + 9 + ) but it’s wakefield driving strength (which depends on the driver laser irradiance I d λ d 2 for LWFA schemes or on the driver beam electric field for PWFA) should be large enough to be able to excite wakefields with amplitudes above the trapping threshold for the extracted electron beam [19]. This contradictory requests have been solved in the two-color [16] scheme by employing a single driver pulse having a long wavelength (so as to increase the irradiance while keeping the pulse electric field below the dopant ionization threshold), thus rising the request of using two laser systems. The ReMPI scheme [17] can also use a single laser system (e.g. a Ti:Sa one), and it employs a train of resonantly delayed pulses, each one having the electric field below the ionization threshold. The particle driven based ionization injection scheme [15] is in a particularly favourable position here, as the electric field generated by electron beams driving large amplitude plasma waves are usually much lower than the ones in equally driving laser pulses. This opens to the possibility of employing ionization processes with low ionization threshold, thus paving the road for ultra-low emittance electron beams.

2. Materials and Methods

In the following we will consider the dynamics of electrons extracted by field (tunnel) ionization and immersed into a plasma wake field driven by ultrarelativistic laser-pulses or dense charged beams. The time scale T e v o l of the driver evolution and of the background (eventual) plasma longitudinal variation will be much longer than the time T c h a r g i n g needed by the electron beam to be extracted and to be trapped by the wakefield, thus enabling us to employ the quasi-static approximation (QSA) [20]. When the particle dynamics is described in a moving window with the same speed v d of the driver (say along the z direction), the constant of motion [5]
γ ( t ) β d u z ϕ ( ξ ( t ) , k p x ( t ) ) h 0 = c o n s t a n t
unveils the longitudinal electron dynamics. In Eq.1  γ is the particle Lorentz factor, β d = v d / c , u = p / m e c and ϕ e Φ / m e c 2 is the normalized scalar potential, k p = 2 π / λ p is the plasma wave wavevector and ξ = k p ( z v d t ) is the longitudinal coordinate inside the window. As any particle extracted by field ionization with very low pulse amplitude is ejected from the parent ion with a negligible momentum (it will get quivering and drift momentum in the laser field right after the ejection), we can easily evaluate the h 0 constant from Eq. 1 as h 0 = 1 0 ϕ e , where ϕ ( ξ ( t = t e ) , k p x ( t = t e ) ) is the normalized potential in the particles position at the extraction time t e . Right after the ejection from the ion, the electron quivers in the laser field being also experiencing longitudinal and transverse ponderomotive forces as well as the longitudinal and transverse forces from the wake field. In the window comoving with the wakefield driver, the electron is seen slipping back towards the rear of the bucket and eventually being trapped by the wakefield. During the slippage, the electron beam is constituted by electrons extracted in different longitudinal positions in the comoving window.
As the Lorentz factor of the electrons increases during their slippage in the back of the bucket, the velocities of the particles approach the speed of light and therefore the electron bunch becomes stiff along the longitudinal direction. The beam compression therefore occurs during the early stages of the beam acceleration and reaches its stable point around the turning point of the longitudinal trajectory in the co-moving window [5], i.e. at the particle trapping point, where the longitudinal speed of the particles is the same of the wakefield, i.e. that of the driver. After straightforward manipulations of Eq. 1 at the trapping point occurring at t t we get
γ , t γ d ϕ t = 1 ϕ e ,
where ϕ t = ϕ ( ξ ( t = t t ) , k p x ( t = t t ) ) is the normalized potential at the trapping point and γ , t = 1 + u 2 ( t = t t ) is the transverse Lorentz factor at the trapping time. Equation 2 states that at the trapping point the normalized potential where the particle is placed ( ϕ t ) depends only on the potential where the particle was born ( ϕ e ) and on the transverse momentum of the particle (through γ , t ) at the trapping time. This link constitutes the basis for the evaluation of the beam length, which will depend on wakefield parameters and on the spread of the longitudinal and transverse position of the particles at their extraction time. We can accurately estimate this spread, as well as the rms residual transverse momentum, by using the theory in [21]. Here, we limit to the unsaturated ionization regime in the tunnel regime and to leading order in the parameter [22] Δ = 3 E 0 2 E a 1 / 2 U H U I 1 / 2 0 . 2 , where E 0 is the ionization pulse peak electric field, E a 0 . 51 T V / m is the atomic field, U I is the ionization energy of the ionization process and U H = 13 . 6 e V is the ionization energy of the Hydrogen atom. Theory in [21,22] shows that for a laser pulse of minimum waist w 0 , i , normalized amplitude a 0 , i and FWHM duration T i , the rms residual transverse momentum ( u , e ), radius ( r e ) and length ( δ z e ) of the beam can be evaluated as
σ ( u , e ) Δ a 0 , i , σ ( r e ) Δ w 0 , i , σ ( δ z e ) 1 2 Δ L 0 , i ,
where L 0 , i = c T i / 2 log 2 .
The wakefield structure can be obtained by analytical results in the blowout regime [23] and by PIC (or fluid simulation if no wave breaking occurs). We performed quasi-3D simulations for the ReMPI scheme in the quasi-linear regime by means of the FB-PIC code [24]. There, a 4.5J Ti:Sa pulse is split into two sub-pulses after a pick-up of 100’s mJ small size beam for the frequency tripled ionization injection. The driving train of two resonantly delayed and circularly polarized pulses had duration of 23 fs, waist of 30 μ m , normalized amplitude of a 0 , d = 0 . 9 , while the linearly polarized ionizing pulse had FWHM duration T i = 30 f s , waist w 0 , i = 4 . 2 μ m and amplitude a 0 , i = 0 . 4 . The ionization pulse phase position in the bucket was varied in the interval 0 . 7 ξ ¯ e 0 . 7 . A plasma target composed by Argon (pre-ionized up to 8th level) with electron density of n e = 7 . 5 · 10 17 c m 3 , corresponding to a plasma wave vector of k p = 0 . 164 μ m 1 is considered. The driver Lorentz factor was γ d = n c / n e 48 where n c = 1 . 1 × 10 21 / λ d 2 = 1 . 7 · 10 21 c m 3 is the critical density for the driver with wavelength λ d = 0 . 81 μ m . The extracted electron beams, having charge of Q = 5 . 2 p C and normalized emittance of about 80 n m r a d , were analyzed in a simulation time at which the trapping were completed. The fields in the simulated cylindrical region were sampled with N m = 3 rotational modes, 2 × 2 × 12 = 48 particles-per-cell deposited in the r , z , θ directions of the cell. The spatio-temporal resolution was d z = λ i / 24 and d r = λ i / 8 in the longitudinal and radial directions, being λ i = λ d / 3 270 n m the ionization pulse wavelength. A snapshot of the fields for the case T i = 30 f s is shown in Figure 1. Those simulations (see Figure 1)) show that the longitudinal field on axis exhibits a linear shape within the whole ionization pulse, with slope ξ E ^ z e = 0 . 38 . For the simulation shown there, the average field at the extraction is E ^ z e = 0 . 3 and the resulting average electric field at trapping is E ^ z e = 0 . 72 , which agrees with the one inferred by asking the particles are trapped close to the peak of the accelerating gradient [17]. We stress here that slight positive or negative delay of the ionization pulse with respect to the field node resulting in | ξ ¯ e | 1 will barely change the electric field at trapping. The radial structure of the electric field also exhibits, as expected, a linear shape close to the axis with gradient k p r e E ^ r e = 0 . 04 , which is much smaller than that obtainable in the blowout regime ( 1 / 4 [23]).

3. Results

In the following we will suppose the ionization pulse size is much smaller than that of the wakefield, i.e.  c T i λ p , and w 0 , i R where T i and w 0 , i are the FWHM pulse duration and minimum waist and R is the wakefield radius. These assumptions will assure us that all the extracted electrons will lie in a similar potential at the extraction time. The constraint w 0 , i R , along with the constrain a 0 , i 1 , will also limit the extent of the transverse dynamics of the extracted electrons, thus maintaining the normalized transverse momentum in the non relativistic range | u , t | 1 . Simulations with the ReMPI scheme show that at the trapping point the bunch transverse size is a fraction of μ m (and usually smaller than the transverse size at the extraction time), while the transverse momentum is close to the one right after the ionization pulse passage. We will also consider non evolving wakefields in the temporal window of the electrons extraction and trapping process, so as to ensure that the potential function ϕ ( ξ , k p r ) is a constant of time. This will limit the amount of the charge for the trapped beam down to the pC level, as beam-loading effects might change the wakefield around the trapped beam position.
For a wakefield with rotational symmetry, the 2nd order Taylor expansion of the normalized potential close to the axis turns out be
ϕ ( ξ 0 + δ ξ , k p r ) ϕ ( ξ 0 , 0 ) E ^ z ( ξ 0 , 0 ) δ ξ 1 2 ξ E ^ z ( ξ , 0 ) | ξ 0 δ ξ 2 + k p r E ^ r ( ξ 0 , k p r ) | 0 ( k p r ) 2 ,
where E ^ z , r = ξ , k p r ϕ = E z , r / E 0 are the longitudinal and radial electric field components normalized to the Dawson field E 0 = m c 2 k p / e .
We can now consider a set of electrons emitted by field ionization at slightly different longitudinal and transverse position ( ξ e , k p r e ) , which will get the trapping point in the position ( ξ t , k p r t ) . The center of mass (phase) position of the beam at the extraction time is ( ξ ¯ e , 0 ) (where ξ ¯ e ξ e ), which lies in the potential ϕ ¯ e = ϕ ( ξ ¯ e , 0 ) . By implicitly defining the reference position ( ξ ¯ t , 0 ) of the beam at its trapping point through the relation 2
γ ¯ , t / γ d ϕ ¯ t = 1 ϕ ¯ e
where ϕ ¯ t = ϕ ( ξ ¯ t , 0 ) and γ ¯ , t = γ , t , we can access the information of the longitudinal spread of the particle position by inverting the nonlinear equation
γ , t γ d ϕ ( ξ ¯ t + δ ξ t , k p r t ) = 1 ϕ ( ξ ¯ e + δ ξ e , k p r e )
on the variable δ ξ t . As we seek at ultrashort pulses of length δ z t σ ( δ ξ t ) / k p 1 (here σ ( x ) is the rms value of the random variable x), at the trapping point the longitudinal gradient of the normalized potential should dominate over the radial one. As a result, particles should be trapped where the accelerating gradient is as high as possible. This is accomplished by setting the wakefield amplitude and the extraction position so as the strong trapping condition is reached [17]. As the radial contribution of the normalized potential ϕ spread at the trapping point is negligible over the longitudinal one, we can express Eq. 7 by using the expansion in Eqs. 4,2 we get:
δ γ , t γ d + E ^ z t δ ξ t E ^ z e δ ξ e + 1 2 ξ E ^ z e δ ξ e 2 + k p r e E ^ r e ( k p r e ) 2 .
where E ^ z t = E ^ z ( ξ ¯ t , 0 ) , E ^ z e = E ^ z ( ξ ¯ e , 0 ) , ξ E ^ z e = ξ E ^ z ( ξ , 0 ) | ξ e , k p r e E ^ r e = k p r e E ^ r ( ξ e , k p r ) | 0 and δ γ , t = γ , t γ ¯ , t 1 2 ( u 2 u 2 ) , which is valid for u 2 1 . The average longitudinal field at extraction ( E ^ z e ), trapping ( E ^ z t ) and field slopes at extraction ( ξ E ^ z e , k p r e E ^ r e ) depend on the phase of the ionization pulse within the bucket and on the wakefield’s regime. They can be analytically inferred in the blowout regime [23] or obtained by simulations. Typical values of ξ E ^ z e and k p r e E ^ r e for the quasi-linear regime at the threshold for the strong trapping or in the blowout regimes can be found in Table 1 (see also Figure 1). We refer there to quasi-3D simulations we performed for a ReMPI setup in the quasi-linear regime (at least during the beam charging/trapping time). Those simulations will be discussed in detail below.
To simplify Eq. 7, we employ the linearity of the longitudinal field in the vicinity of the node (where the ionization pulse is placed), thus substituting the average field E ^ z e at the extraction with E ^ z e ξ E ^ z e · ξ ¯ e for | ξ ¯ e | 1 :
E ^ z t δ ξ t ξ E ^ z e ξ ¯ e δ ξ e + 1 2 δ ξ e 2 + 1 2 k p r e E ^ r e ( k p r e ) 2 δ γ , t γ d .
Let’s first evaluate the average of the phase positions spread at trapping ξ t ¯ = ξ t by averaging Eq. 8 on the random variables δ ξ e and r r , which are independently distributed
E ^ z t δ ξ t 1 2 ξ E ^ z e σ 2 ( ξ e ) + k p r e E ^ r e k p 2 σ 2 ( r e ) .
We aim at evaluating the variance of δ ξ t so we need to average the square of the right hand side of Eq. 8. This involves a terms a δ γ , t 2 / γ d 2 1 4 σ 2 ( u ) / γ d 2 3 4 Δ 2 a 0 , i 2 / γ d 2 which can be neglected as γ d 1 , a 0 , i 2 1 and Δ 2 = O ( 10 3 ) . In the following we will assume γ d be sufficiently large to neglect all the terms containing δ γ 2 / γ d 2 . The evaluation of the δ ξ t variance σ 2 ( δ ξ t 2 ) = δ ξ t 2 δ ξ r 2 finally leads to:
2 ( E ^ z t ) 2 σ 2 ( δ ξ t ) ( ξ E ^ z e ) 2 2 ξ ¯ e 2 σ 2 ( δ ξ e ) + σ 4 ( δ ξ e ) + ( k p r e E ^ r e ) 2 k p 4 σ 4 ( r e ) .
We can immediately see from Eq. 10 that if all the wakefield parameters and all the ionization pulse parameters (but the ionization pulse phase ξ ¯ e ) are fixed, the rms value of the longitudinal positions at trapping σ ( z t ) get its minimum value when the ionization pulse is placed in the vicinity of the node (where ξ ¯ e = k p z ¯ e 0 ).
We evaluate the terms of Eq. 10 for the reference case given by the ReMPI simulations we run and for the case of the fully evacuated bubble. For the blowout case we will use the analytical results from [23] and we will assume that both a background plasma and an ionization pulse with the same parameters of the ReMPI PIC simulations will be employed. We also evaluated the neglected term related to δ γ t / γ d , to check the accuracy of the approximation we made.
Results are shown in Table 2, where it is apparent that the transverse contribution are much smaller than the longitudinal ones in the case of tightly focused ionization pulses ( c T i w 0 , i ). Also, the terms related to the transverse momentum spread at the trapping point is negligible for the cases considered above. The rms longitudinal beam size at trapping can be directly expressed in terms of physical quantities as
σ ( δ z f ) k p Δ 2 1 2 E z t ( ξ E ^ z e ) 2 1 4 + z ¯ e Δ L i 2 L i 4 + ( k p r e E ^ r e ) 2 w 0 , i 4
where z ¯ e is the position of the ionization pulse with respect to the node of the electric field (where we set ξ = 0 in Figure 1). Equation Eq. 11 gives us the route to obtain ultrashort electron beams. Firstly, as stated above the electric field at the trapping point should be as large as possible. If a quasi-linear regime is employed, the strong trapping condition [17] assures that the electron beam is trapped at the peak of the electric field. Secondly, as the beam length scales as k p n e 1 / 2 , tenuous plasmas might be preferred also because they led to higher values for γ d in the case of a laser driver. Finally, as the term z ¯ e / ( Δ L i ) in Eq. 11 refers to the ratio between the distance of the ionization pulse to the node ( z ¯ e ) and the longitudinal extension of the beam at the extraction time ( Δ L i / 2 ), we can obtain the shortest electron beams by placing the ionization pulse on the node of the accelerating gradient, with an acceptable jitter of the scale of a fraction of the ionization pulse length, i.e. | z ¯ e | ( Δ L i ) .
Results for dependence the rms beam duration σ ( δ t e ) = σ ( δ z e ) / c on the delay of the ionization pulse from the accelerating gradient node z ¯ e / c are shown for the case of the ReMPI setup with a 200TW Ti:Sa laser system selected here and with T = 30 f s long ionization pulse in third harmonics. Predictions from the model (Eq. 11) do overlap with the observed beam duration, which show that electron beams with duration as short as 300 as can be generated in this way. As both the ionization pulse and the driver train are generated by the same amplified pulse, the relative time jitter between them is only due to mechanical vibration of the optical elements after the pick-up and it can therefore be limited down to a few femtosecond scale.
Figure 2. Beam duration rms at the trapping point for the set of ReMPI simulations with ionization pulse duration of 30 f s FWHM vs the delay z ¯ e / c of the ionization pulse with respect to the node of the electric field. The orange line refers to the prediction by Eq. 11.
Figure 2. Beam duration rms at the trapping point for the set of ReMPI simulations with ionization pulse duration of 30 f s FWHM vs the delay z ¯ e / c of the ionization pulse with respect to the node of the electric field. The orange line refers to the prediction by Eq. 11.
Preprints 85045 g002

4. Discussion

We discussed here a simple model to infer the beam duration of electron beams obtained by high-quality ionization injection schemes (e.g. two-color, ReMPI or trojan-horse) with ultrarelativistic wakefield driver ( γ d 1 ), low amplitude ( a 0 , i 1 ), tightly focused ( w 0 , i R ) and short ( k p c T 1 / Δ ) ionization pulses. Those conditions are naturally satisfied with the trojan-horse scheme, which employs electron beams with γ = γ d 1 as drivers. As the ReMPI scheme can work with Ti:Sa pulses, at plasma densities in the range of ( 1 · 10 17 1 · 10 18 ) c m 3 the driver relativistic factor γ d = n e / n c is in the range 40 130 . The two-color scheme needs a long wavelength driver ( λ d = 5 μ m in [16]), so the achievement of the condition γ d 1 requires much lower plasma densities (a factor of five less to get the same γ d with a Ti:Sa pulse). The model was tested against a set of ReMPI simulations in the quasi-linear regime, with γ d 50 . Both the minimum value of the beam duration (about 300 as rms) and the dependence of duration on the delay of the ionization pulse from the accelerating gradient node, agree with a percent error. The model doesn’t take into account beam loading effects, curvature effects of the ξ u z trajectories at the inversion point (which decrease by employing high values of γ d ) and the eventual evolution of the wakefield during the beam charging which can be caused by driver evolution or non flat longitudinal background plasma profiles. The first effect can be mitigated by reducing the amount of the extracted charge, the second is virtually negligible in beam driven schemes and it can be mitigated in laser driven schemes by reducing the plasma density. The third one is negligible in trojan-horse like schemes with flat plasma profiles, while it can be reduced in LWFA by choosing very small waist sizes for the ionization injection pulse and/or driver pulses close to their foci. The effect of the ionization pulse-to-driver temporal jitter was also discussed. For a ReMPI scheme driven by a single Ti:Sa pulse the time jitter of the ionization pulse to the node of the accelerating gradient can be limited down to a few femtoseconds, resulting in a slight increase of the trapped beam length for the case of ionization pulses of duration T = 30 f s FWHM. A similar scenario can be envisaged in trojan-horse schemes for which the photocathode laser pulse and the ionization pulse can be obtained by the same laser system.

Acknowledgments

Extreme Light Infrastructure Nuclear Physics (ELI-NP) Phase II, is a project co-financed by the Romanian Government and the European Union through the European Regional Development Fund and the Competitiveness Operational Programme (1/07.07.2016, COP, ID 1334). This work was supported by the contract sponsored by the Romanian Ministry of Research and Innovation: PN 23 21 01 05, and the IOSIN funds for research infrastructures of national interest. Access to HPC through PRACE EHPC-BEN-2023B05-023 is also acknowledged.

References

  1. Krausz, F.; Ivanov, M. Attosecond physics. Reviews of modern physics 2009, 81, 163.
  2. Tajima, T.; Dawson, J.M. Laser Electron Accelerator. Phys. Rev. Lett. 1979, 43, 267–270. [CrossRef]
  3. Chen, P.; Dawson, J.M.; Huff, R.W.; Katsouleas, T. Acceleration of electrons by the interaction of a bunched electron beam with a plasma. Physical review letters 1985, 54, 693.
  4. Khachatryan, A.G. Trapping, compression, and acceleration of an electron beam by the laser wake wave. Journal of Experimental and Theoretical Physics Letters 2001, 74, 371–374.
  5. Khachatryan, A.G. Trapping, compression, and acceleration of an electron bunch in the nonlinear laser wakefield. Physical Review E 2002, 65, 046504.
  6. Khachatryan, A.; Van Goor, F.; Boller, K.J.; Reitsma, A.; Jaroszynski, D. Extremely short relativistic-electron-bunch generation in the laser wakefield via novel bunch injection scheme. Physical review special topics-Accelerators and Beams 2004, 7, 121301.
  7. Li, F.; Sheng, Z.; Liu, Y.; Meyer-ter Vehn, J.; Mori, W.; Lu, W.; Zhang, J. Dense attosecond electron sheets from laser wakefields using an up-ramp density transition. Physical review letters 2013, 110, 135002.
  8. Tooley, M.; Ersfeld, B.; Yoffe, S.; Noble, A.; Brunetti, E.; Sheng, Z.; Islam, M.; Jaroszynski, D. Towards attosecond high-energy electron bunches: Controlling self-injection in laser-wakefield accelerators through plasma-density modulation. Physical review letters 2017, 119, 044801.
  9. Luttikhof, M.; Khachatryan, A.; Van Goor, F.; Boller, K.J. Generating ultrarelativistic attosecond electron bunches with laser wakefield accelerators. Physical review letters 2010, 105, 124801.
  10. Horny, V.; Krus, M.; Yan, W.; Fulop, T. Attosecond betatron radiation pulse train. Scientific Reports 2020, 10, 15074.
  11. Zhu, X.L.; Liu, W.Y.; Chen, M.; Weng, S.M.; He, F.; Assmann, R.; Sheng, Z.M.; Zhang, J. Generation of 100-MeV attosecond electron bunches with terawatt few-cycle laser pulses. Physical Review Applied 2021, 15, 044039.
  12. Ferri, J.; Hornỳ, V.; Fülöp, T. Generation of attosecond electron bunches and x-ray pulses from few-cycle femtosecond laser pulses. Plasma Physics and Controlled Fusion 2021, 63, 045019.
  13. Deng, A.; Li, X.; Luo, Z.; Li, Y.; Zeng, J. Generation of attosecond micro bunched beam using ionization injection in laser wakefield acceleration. Optics Express 2023, 31, 19958–19967.
  14. Hidding, B.; Rosenzweig, J.; Xi, Y.; O’Shea, B.; Andonian, G.; Schiller, D.; Barber, S.; Williams, O.; Pretzler, G.; Königstein, T.; others. Beyond injection: Trojan horse underdense photocathode plasma wakefield acceleration. AIP Conference Proceedings. American Institute of Physics, 2012, Vol. 1507, pp. 570–575.
  15. Hidding, B.; Assmann, R.; Bussmann, M.; Campbell, D.; Chang, Y.Y.; Corde, S.; Cabadağ, J.C.; Debus, A.; Döpp, A.; Gilljohann, M.; others. Progress in hybrid plasma wakefield acceleration. Photonics. MDPI, 2023, Vol. 10, p. 99.
  16. Yu, L.L.; Esarey, E.; Schroeder, C.B.; Vay, J.L.; Benedetti, C.; Geddes, C.G.R.; Chen, M.; Leemans, W.P. Two-Color Laser-Ionization Injection. Phys. Rev. Lett. 2014, 112, 125001. [CrossRef]
  17. Tomassini, P.; De Nicola, S.; Labate, L.; Londrillo, P.; Fedele, R.; Terzani, D.; Gizzi, L.A. The resonant multi-pulse ionization injection. Physics of Plasmas 2017, 24, 103120. [CrossRef]
  18. Ammosov, M.V.; Delone, N.B.; Krainov, V.P. Tunnel Ionization Of Complex Atoms And Atomic Ions In Electromagnetic Field. High Intensity Laser Processes; Alcock, J.A., Ed. International Society for Optics and Photonics, SPIE, 1986, Vol. 0664, pp. 138 – 141. [CrossRef]
  19. Esarey, E.; Pilloff, M. Trapping and acceleration in nonlinear plasma waves. Physics of Plasmas 1995, 2, 1432–1436.
  20. Sprangle, P.; Esarey, E.; Ting, A. Nonlinear interaction of intense laser pulses in plasmas. Phys. Rev. A 1990, 41, 4463–4469. [CrossRef]
  21. Tomassini, P.; Massimo, F.; Labate, L.; Gizzi, L.A. Accurate electron beam phase-space theory for ionization-injection schemes driven by laser pulses. High Power Laser Science and Engineering 2022, 10, e15.
  22. Schroeder, C.; Vay, J.L.; Esarey, E.; Bulanov, S.; Benedetti, C.; Yu, L.L.; Chen, M.; Geddes, C.; Leemans, W. Thermal emittance from ionization-induced trapping in plasma accelerators. Physical Review Special Topics-Accelerators and Beams 2014, 17, 101301.
  23. Pukhov, A.; Gordienko, S.; Kiselev, S.; Kostyukov, I. The bubble regime of laser–plasma acceleration: monoenergetic electrons and the scalability. Plasma physics and controlled fusion 2004, 46, B179.
  24. Lehe, R.; Kirchen, M.; Andriyash, I.A.; Godfrey, B.B.; Vay, J.L. A spectral, quasi-cylindrical and dispersion-free Particle-In-Cell algorithm. Computer Physics Communications 2016, 203, 66–82. [CrossRef]
Figure 1. Field structure for one of the ReMPI simulations in the quasi-linear regime and for a case in which the ionization pulse is delayed from its optimal position in ξ = ξ ¯ e = 0 . a) Line-out on the axis: the accelerating field (orange) is in units of E 0 , the laser pulses transverse field (blue) is in arbitrary units. The ionization pulse is placed in the phase ξ ¯ e = 0 . 7 , with an average longitudinal field at extraction E ^ z e = 0 . 3 . The longitudinal field at trapping ( ξ ¯ t = 2 . 3 ) is E ^ z e = 0 . 72 . The dashed (green) line shows a linear fit of the accelerating gradient at ξ = 0 . b) Transverse lineout of the transverse field normalized to E 0 in the node of the longitudinal field (blue) and on the trapping point (orange). The dashed (green) line shows a linear fit of the transverse gradient at ξ = 0 . c) 2D map of n ( ξ , k p r ) / n e 1 showing that a quasi-linear regime is obtained. The dashed line represents the accelerating gradient in arb. units.
Figure 1. Field structure for one of the ReMPI simulations in the quasi-linear regime and for a case in which the ionization pulse is delayed from its optimal position in ξ = ξ ¯ e = 0 . a) Line-out on the axis: the accelerating field (orange) is in units of E 0 , the laser pulses transverse field (blue) is in arbitrary units. The ionization pulse is placed in the phase ξ ¯ e = 0 . 7 , with an average longitudinal field at extraction E ^ z e = 0 . 3 . The longitudinal field at trapping ( ξ ¯ t = 2 . 3 ) is E ^ z e = 0 . 72 . The dashed (green) line shows a linear fit of the accelerating gradient at ξ = 0 . b) Transverse lineout of the transverse field normalized to E 0 in the node of the longitudinal field (blue) and on the trapping point (orange). The dashed (green) line shows a linear fit of the transverse gradient at ξ = 0 . c) 2D map of n ( ξ , k p r ) / n e 1 showing that a quasi-linear regime is obtained. The dashed line represents the accelerating gradient in arb. units.
Preprints 85045 g001
Table 1. Longitudinal and radial gradients of the electric field at the node of the longitudinal electric field (see Figure 1).
Table 1. Longitudinal and radial gradients of the electric field at the node of the longitudinal electric field (see Figure 1).
Parameter Quasi-linear Spherical bubble (theory)
ξ E ^ z e 0.38 0.50
k p r e E ^ r e 0.04 0.25
Table 2. Estimation of the terms in Eq. 10.
Table 2. Estimation of the terms in Eq. 10.
Parameter Quasi-linear Spherical bubble (theory)
[ ξ E ^ z e σ ( δ ξ e ) 2 ] 2 3 · 10 4 4 · 10 4
[ k p r e E ^ r e k p 2 σ ( r e ) 2 ] 2 1 · 10 6 5 · 10 5
[ σ ( δ γ t ) / γ d ] 2 4 · 10 8 4 · 10 8
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated