Electron Acceleration by Beating of Two Intense Cross-Focused Hollow Gaussian Laser Beams in Plasma
Mahmoud Saleh T.1, Gauniyal Rakhi2, Ahmad Nafis1, Rawat Priyanka2, Purohit Gunjan2, †
Department of Physics, College of Science, UAE University, P.O. Box 15551 Al-Ain, United Arab Emirates
Laser Plasma Computational Laboratory, Department of Physics, DAV (PG) College, Dehradun, Uttarakhand-248001, India

 

† Corresponding author. E-mail: gunjan75@gmail.com

Abstract
Abstract

This paper presents propagation of two cross-focused intense hollow Gaussian laser beams (HGBs) in collisionless plasma and its effect on the generation of electron plasma wave (EPW) and electron acceleration process, when relativistic and ponderomotive nonlinearities are simultaneously operative. Nonlinear differential equations have been set up for beamwidth of laser beams, power of generated EPW, and energy gain by electrons using WKB and paraxial approximations. Numerical simulations have been carried out to investigate the effect of typical laser-plasma parameters on the focusing of laser beams in plasmas and further its effect on power of excited EPW and acceleration of electrons. It is observed that focusing of two laser beams in plasma increases for higher order of hollow Gaussian beams, which significantly enhanced the power of generated EPW and energy gain. The amplitude of EPW and energy gain by electrons is found to enhance with an increase in the intensity of laser beams and plasma density. This study will be useful to plasma beat wave accelerator and in other applications requiring multiple laser beams.

1 Introduction

Laser driven plasma-based accelerators are of great interest because of their ability to sustain extremely large acceleration gradients.[14] They have widespread applications in various fields such as particle physics, materials science, medicine (ranging from x-ray diagnostics to particle beam therapies), manufacturing industry and to produce extremely short electron bunches. In laser plasma interaction, various schemes are operative for particle acceleration such as laser wakefield accelerator (LWFA),[57] plasma beat wave accelerator (PBWA),[5, 8] multiple laser pulses,[910] self-modulated LWFA[1112] etc. in which very high intense laser beams are used. The key requirement for achieving higher acceleration is the large amplitude of the electron plasma wave. A large amplitude plasma wave (with phase velocity near the speed of light) is driven by intense laser pulse in which plasma electrons can be trapped and accelerated to relativistic energies. Therefore, the excitation of electron plasma wave by intense laser beams in plasmas has been an active field of research for charged particle acceleration.

There has been significant interest in the beat wave excitation of electron plasma wave and its application in plasma based particle accelerators. The plasma beat-wave accelerator scheme is one of the most mature methods for plasma acceleration in laser plasma interaction.[5] In this scheme, two intense laser beams of slightly different frequencies ω1 & ω2 and corresponding wave numbers k1 & k2 are used to resonantly excite an electron plasma wave. The resonance conditions for beat wave excitation of electron plasma wave are: and . If the beat frequency is close to the plasma frequency ωp0, a very large amplitude relativistic electron plasma wave can be generated, which can be used to accelerate electrons efficiently to high energies in short distances. In this process, self-focusing of an intense laser beam in plasmas plays an important role and arises due to increase of the on-axis index of refraction relative to the edge of the laser beam. When two intense laser beams with slightly different frequencies simultaneously propagates in plasma, self-focusing of one laser beam is affected by the self-focusing of another laser beam. Due to mutual interaction of two laser beams, cross-focusing takes place in plasma. This method is used to generate large amplitude electron plasma wave for ultrahigh gradient electron acceleration in PBWA scheme. A lot of theoretical and experimental work has been reported on the excitation of electron plasma wave and electron acceleration by beating of two laser beams in plasma.[5, 1337]

It has been observed that the focusing of laser beams in plasma and the yield of electron acceleration depends on the spatial profile of laser beams and the nonlinearities associated with plasma. Most of the earlier investigations on excitation of electron plasma wave and electron acceleration have been carried out by taking Gaussian intensity distribution of laser beams with ponderomotive/relativistic nonlinearity. In contrast with Gaussian profile of laser beams, currently, the hollow Gaussian intensity profile of laser beams with minimum field intensity at the center has attracted much more attentions because of their unique physical properties and applications.[3840] The main feature of considered hollow Gaussian laser beams is having the same power at different beam orders with null intensity at the centre. The propagation dynamics of hollow Gaussian beams (HGBs) have been widely studied both experimentally and theoretically.[4146] Various techniques[4750] have been used to generate HGBs. These beams can be expressed as a superposition of a series of Lagurerre-Gaussian modes.[41] Moreover, when an ultra-intense laser beams propagates in plasma, both relativistic and ponderomotive nonlinearities arise simultaneously due to electron mass variation and electron density perturbations respectively, which depend on the time scale of the laser pulse.[5152] Therefore, in comparison with only relativistic/ponderomotive nonlinearity, the dynamics of the propagation of laser beams in plasma is expected to be drastically affected by cumulative effects of relativistic-ponderomotive nonlinearity. It is also important to mention that the relativistic effect and ponderomotive nonlinearity contribute to focusing on a femtosecond time scale at very high intensity. Cross-focusing of HGBs with relativistic/ponderomotive or relativistic-ponderomotive nonlinearities in plasma have been investigated in detail[5355] but no one has studied the excitation of EPW and electron acceleration. Moreover, donut wakefields generation for particle acceleration by Laguerre-Gaussian laser pulses (carrying a finite amount of orbital angular momentum) using particle-in-cell simulations have also been reported in under-dense plasma.[5657]

In the present study, we have considered the propagation of two intense hollow Gaussian beams in collisionless plasma. The intensity distribution of the beams along the axis is zero and the maximum is away from the axis. We have studied the cross-focusing of two intense hollow Gaussian laser beams in collisionless plasma with relativistic and ponderomotive nonlinearities and further its effect on the excitation of EPW and electron acceleration. The paraxial-ray approximation[5859] is used to describe the focal region of the laser beam where all the relevant parameters correspond to a narrow range around the maximum irradiance of the HGBs. Section 2 presents the relativistic-ponderomotive focusing of two HGBs in plasma under paraxial-ray approximation. The effect of the cross-focusing of the HGBs on the excitation of the electron plasma wave and electron acceleration is studied in Sec. 3. Section 4 deals with the results and discussion. The conclusions drawn from the results of present investigation have been summarized in Sec. 5.

2 Propagation of Hollow Gaussian Laser Beams in Collisional Plasma

We consider two intense linearly polarized co-axial hollow Gaussian laser beams of frequencies ω1 and ω2 propagating along the z-direction in the collisional plasma. The initial electric field distribution of the beams are given by

where is the radial coordinate of the cylindrical coordinate system, and are the initial spot size of the laser beams, and are the amplitude of HGBs maximum , is the order of the HGBs ( corresponds to the fundamental Gaussian laser beam) and a positive integer, characterizing the shape of the HGBs an position of its irradiance maximum.

For hollow Gaussian laser beams, transforming the coordinate in to () coordinate by the relation[44]

where η is a reduced radial coordinate, is the beam width of HGBs, is the position of the maximum irradiance for the propagating beams. Most of the power of HGBs are concentrated in the region around .

The ponderomotive force acting on the electrons in the presence of relativistic nonlinearity is given by[5253]

where γ is the relativistic Lorentz factor and is given by

here

This force modifies the background electron density because electrons are expelled away from the region of higher electric field. The modified electron density due to relativistic-ponderomotive force is given as[52]

Using Eq. (3) into Eq. (7), the total electron density can be written as

Here and are the intensities of first and second laser beams.

The dielectric constant of the plasma is given by

where is the plasma frequency with e is the charge of an electron, is its rest mass, and is the density of plasma electrons in the absence of laser beam. The effective dielectric constant () of the plasma (in the presence of relativistic-ponderomotive nonlinearity) at frequencies ω1 and ω2 respectively is given by

where is the linear part of dielectric constant of the plasma. The dielectric constant may be expended around the maximum of HGBs i.e. at η = 0. For HGBs, one can express the dielectric constant , in terms of . Expanding dielectric constant around η = 0 by Taylor expansion, one can write

where

with , , and are the coefficients associated with η0 and η2 in the expansion of around .

The wave equation governing the electric vectors of the two laser beams in plasma can be written as

In Eq. (10) the second term can be neglected, which is justified as long as in . One can thus write the wave equation for the two laser beams in plasma is as,

The solution of Eq. (11) can be written as[58]

where is the complex amplitude of the electric field , and is the dielectric function, equivalent to the highest electric field on the wave front of the HGBs. Substituting from Eq. (12) into Eq. (11), the wave equation becomes

The complex amplitude may be expressed as

where is a complex function of space and is the eikonal associated with the HGBs. Substituting Eq. (14) into Eq. (13) and separating the real and imaginary parts of the resulting equation, the following set of equations is obtained.

In the paraxial approximation i.e. when , Eqs. (15) and (16) lead to the maintenance of the HGBs character during propagation. Using Eq. (3), Eqs. (15) and (16) in terms of η can be expressed as

The solution of Eq. (17) can be written as[5859]

where Here is an arbitrary function of z and is the beam width parameter for first and second HGBs respectively. Substituting from Eqs. (19) and (20) into Eq. (17) and equating the coefficients of η0 and η2 on both sides of the resulting equation, one obtains

where is the dimensionless distance of propagation is the dimensionless initial beam width.

3 Beat Wave Excitation of Plasma Wave and Electron Acceleration

Due to beating of two hollow Gaussian laser beams in plasma and modification in plasma density via relativistic-ponderomotive force, electron plasma wave (EPW) is generated. In this process, the contribution of ions is negligible because they only provide a positive background, i.e., only plasma electrons are responsible for excitation of EPW. The amplitude of electron plasma wave (which depends on the background electron density) gets strongly coupled to the laser beams. To study the effect of relativistic and ponderomotive nonlinearities on the generation of the plasma wave by the beat wave process in paraxial region, we start with the following fluid equations:

(i) The continuity equation:

(ii) The momentum equation:

(iii) Poisson’s equation:

where Ne is the total electron density, V is the sum of the drift velocities of the electron in the EM field and self-consistent field, is the Landau damping factor,[60] E and B are the self-consistent electric and magnetic fields of laser beams corresponding to the electron plasma wave; and other symbol have their usual meaning.

Using Eq. (22) to Eq. (24), we obtain the following equation governing the electron plasma wave

where is the electron thermal speed. At the difference frequency , Eq. (25) reduces to

where is the component of electron density oscillating at frequency . The drift velocities of electron in the pump field at the frequency are . Therefore, one can write

Equation (26) contains two plasma waves (both have different frequencies). The first one is supported by hot plasma and the second by the source term at the different frequency. If , the phase velocity is almost equal to the thermal velocity of the electron and Landau damping occur. But as , the phase velocity is very large as compared to the electron thermal velocity and Landau damping is negligible. The solution of Eq. (9) with in the WKB approximation can be expressed as

where , N10 and N20 are the slowly varying real function of the space coordinates. Substituting for N1 from Eq. (28) into Eq. (26) the equation for N10 and N20 can be obtained by equating the coefficient of and in the resulting equation.

After separating the real and imaginary parts of Eq. (29), we obtain

The equation for N20 is given by

The solution of Eqs. (31) and (32) can be written as[58, 59]

where , a0 and B are the constants to be determined by the boundary condition that the magnitude of the generated plasma wave at z = 0 is zero. This gives

and . In Eq. (34), fp is the dimensionless beam width parameter of the plasma wave given by

Using Poisson’s equation one can obtain the electric vector of the plasma wave excited at the difference frequency:

where

The power associated with electron wave is given by[29]

where is the group velocity. Equation (38) has been solved numerically with the help of Eqs. (21), (36), and (37).

The excited electron plasma wave transfers its energy to accelerate the electrons at the difference frequency of laser beams. The equation governing electron momentum and energy are

where is the energy gain by the electron or electron energy. This equation can be written as

where is given by Eqs. (37) and (39) gives the electron energy. Equation (39) has been solved numerically, where we have used f1 and f2 by Eq. (21).

4 Numerical Results and Discussion

In order to study the cross-focusing of two intense hollow Gaussian laser beams in collisonal plasma with dominant relativistic-ponderomotive nonlinearity and its effect on the generation of electron plasma wave and electron acceleration in paraxial-ray approximation, numerical computation of Eqs. (19), (21), (38), and (39) has been performed for the typical laser plasma parameters. The following set of the parameters has been used in the numerical calculation: , , , , , , , 2.4, 2.8, and , 2, 3. For initial plane wave fronts of the beams, the initial conditions for are and at .

When two intense hollow Gaussian laser beams of slightly different frequencies simultaneously propagate through the plasma, the background electron density of plasma is modified due to ponderomotive force. Equation (19) describes the intensity profile of HGBs in plasma along the radial direction in presence of relativistic and ponderomotive nonlinearities, while Eq. (21) determines the focusing/defocusing of the hollow Gaussian beams along the distance of propagation in the plasma, where the first term on the right hand side shows the diffraction behaviour of laser beams and remaining terms on the right hand side show the converging behaviour that arises due to the relativistic and ponderomotive nonlinearity. The focusing/defocusing behaviour of laser beams depends on the magnitudes of the nonlinear coupling terms. It is clear from Eq. (19), the intensity profile of both laser beams depends on the focusing/defocusing of the laser beams i.e. beam width parameters .

The results of Eqs. (19) and (21) are presented in Figs. 1 and 2 respectively. Figures 1(a)1(c) illustrate the focusing behavior of two hollow Gaussian laser beams in plasma with normalized distance of propagation. It is clear that both of the beams show oscillatory self focusing. One can see from Fig. 1(a) as the order of the hollow Gaussian beam (n) increases focusing of both beams also increases, which implies that with the increase in n the hollow space across the beam also increases. It is observed from Figs. 1(b) and 1(c) that with an increase in the intensity of incident laser beam and plasma density there is an increase in the extent of self-focusing of both the laser beams. This is due to the fact that with an increase in plasma density the number of electrons contributing to the pondermotive-relativistic nonlinearity also increases. These results indicate that focusing of one beam is significantly affected by the presence of another beam. Figures 2(a) and 2(b) represent the normalized intensity of first and second laser beam in plasma respectively along with the normalized distance of propagation for different n. It is observed from the figures that with an increase in n, normalized intensity of both laser beams also increases. It has also been observed that the intensity of laser beams is higher at higher plasma density and incident laser beam intensity (results not shown here). This is due to the strong focusing of both HGBs at higher values of n and plasma density respectively. Such kind of highly self-focused intense laser beams interact with each other and generate large amplitude electron plasma wave.

Fig. 1 (Color online) Variation of the beam width parameters and with the normalized distance of propagation (ξ) of first and second HGBs: (a) For different orders of n (red, black and blue colors are for n = 1, 2 and 3) with , and . (b) For different values of (red, black, and blue colors are for , 2.4, 2.8) with , and . (c) For different values of ωp0 (red, black, and blue colors are for , , and ) with , . and , respectively, when relativistic and ponderomotive nonlinearities are operative. The solid line represents and the dotted line represents .
Fig. 2 (Color online) Variation of the normalized intensity of HGBs with the normalized distance of propagation (ξ) for different orders of n (red, black and blue colors are for , 2 and 3) keeping , , , when relativistic and ponderomotive nonlinearities are operative. (a) and (b) are for first and second laser beam.

To see the effect of cross-focusing of two intense hollow Gaussian laser beams on the generation of electron plasma wave, we have studied the excitation of the electron plasma wave by the beat wave process in the presence of relativistic and ponderomotive nonlinearities. It is important to mention here that beating between two co-propagating intense laser beams in plasma can generate a longitudinal electron plasma wave with a high electric field and a relativistic phase velocity. The beat plasma wave is very efficient for electron acceleration up to ultrarelativistic energies. Equation (38) describes the normalized power of electron plasma wave generated as a result of beating of the two laser beams. It is evident from Eq. (38), the power of electron plasma wave depends upon the focusing of laser beams in plasma , focusing of electron plasma wave and the electric field associated with electron plasma wave . The same set of parameters has been used for numerical calculations as in Figs. 12. Figures 3(a) and 3(b) show the effect of changing the order of HGBs and the initial intensity of second laser beam on the power of generated electron plasma wave, at the position of optimum irradiance of the beams (η = 0). By increasing laser beam orders and the intensity of second laser beam, the power of generated electron plasma wave increases. This is because the power of generated electron plasma wave depends on the focusing behaviour of laser beams , and intensity of laser beams.

Fig. 3 (Color online) Variation of the normalized power (P) of electron plasma wave with normalized distance (ξ): (a) For different orders of n (red, black and blue colors are for n = 1, 2 and 3) with , , and . (b) For different values of (red, black, and blue colors are for , 2.4, 2.8) with , , and , respectively, when relativistic and ponderomotive nonlinearities are operative.

The large amplitude electron plasma wave can be used to accelerate the electrons in plasma beat wave accelerator scheme. Equation (39) gives the expression for energy gain by the electrons. This equation has been solved numerically with the help of Eq. (37) i.e. energy gain depends on the electric field associated with excited electron plasma wave. Figures 4(a) and 4(b) respectively illustrate the effect of changing the value of n and the intensity of second laser beam on energy gain by the electrons. These figures clearly indicate that the maximum energy gain by electrons is significantly increased by increasing the value of b and . Although, the nature of these results are similar to Figs. 3(a) and 3(b) due to the same reasons as discussed above. Thus, we see that cross-focusing of laser beams plays a crucial role for efficient generation of electron plasma wave and electron acceleration.

Fig. 4 (Color online) Variation of the (γ) of the electron with normalized distance (ξ): (a) For different orders of n (red, black and blue colors are for , 2 and 3) with , , and . (b) For different values of (red, black, and blue colors are for , 2.4, 2.8) with , , and , respectively, when relativistic and ponderomotive nonlinearities are operative.
5 Conclusions

In conclusion, plasma wave generation and electron acceleration by beating of two intense HGBs in collisionless plasma with dominant relativistic-ponderomotive nonlinearity has been discussed. Paraxial-ray approximation has been used to establish the given formulation. Effects of laser and plasma parameters on the focusing of two HGBs in plasma at difference frequency, generation of electron plasma wave and particle acceleration process is examined. Following conclusions are made from the results:

(i) The order of the hollow Gaussian laser beams plays an important role in plasma beat wave accelerator scheme.

(ii) Focusing of both hollow Gaussian laser beams in plasma is enhanced by increasing the laser beam orders, incident laser intensity, and plasma density.

(iii) The intensity of both laser beams in plasma is higher for higher orders of the HGBs.

(iv) Maximum power of generated electron plasma wave depends on the extent of focusing of laser beams in plasma and the electric field associated with electron plasma wave respectively.

(v) The power of electron plasma wave increases by increasing the laser beam orders and the initial intensity of second/first laser beams.

(vi) The maximum energy gain also depends on the electric field associated with electron plasma wave and increases by increasing the laser beam orders and the initial intensity of second/first laser beams.

The results of the present investigation are relevant to laser beat wave based particle accelerators, terahertz generation and in other applications requiring multiple laser beams.

Reference
[1] Hooker S. M. Nat. Photon. 7 2013 775
[2] Esarey E. Schroeder C. B. Leemans W. P. Rev. Mod. Phys. 81 2009 1229
[3] Malka V. Faure J. Gauduel Y. A. Nat. Phys. 4 2008 447
[4] Bingham R. Phil. Trans. R. Soc. A. 364 2006 559
[5] Tajima T. Dawson J. M. Phys. Rev. Lett. 43 1979 267
[6] Pukhov A. Meyer-ter-vehn J. Appl. Phys. B 74 2002 355
[7] Albert F. Thomas A. G. R. Mangles S. P. D. Plasma Phys. Control. Fusion 56 2014 084015
[8] Filip C. V. Narang R. Tochitsky S. Ya. Phys. Rev. E 69 2004 026404
[9] Nakajima K. Phys. Rev. A 45 1992 1149
[10] Umstadter D. Esarey E. Kim J. Phys. Rev. Lett. 72 1994 1224
[11] Leemans W. P. Catravas P. Esarey E. Phys. Rev. Lett. 89 2002 174802
[12] Chen W. Chien T. Lee C. Phys. Rev. Lett. 92 2004 075003
[13] Rosenbluth M. N. Liu C. S. Phys. Rev. Lett. 29 1972 701
[14] Joshi C. Mori W. B. Katsouleas T. Nature (London) 311 1984 525
[15] Kitagawa Y. Matsumoto T. Minamihata T. Phys. Rev. Lett. 68 1984 48
[16] Tang C. M. Sprangle P. Sudan R. N. Appl. Phys. Lett. 45 1984 375
[17] Modena A. Najmudin Z. Dangor A. E. Nature (London) 377 2002 606
[18] Clayton C. E. Joshi C. Darrow C. Umstadter D. Phys. Rev. Lett. 54 1985 2343
[19] Ebrahim N. A. Lavigne P. Aithal S. IEEE Trans. Nucl. Sci. 32 1985 3539
[20] Clayton C. E. Marsh K. A. Dyson A. Phys. Rev. Lett. 70 1993 37
[21] Everett M. Lal A. Gordon D. Nature (London) 368 1994 527
[22] Ebrahim N. A. J. Appl. Phys. 76f50 1994 7645
[23] Amiranoff F. Bernard D. Cros B. Phys. Rev. Lett. 74 1995 5220
[24] Ghizzo A. Bertrand P. Lebas J. Phys. Plasmas 5 1998 4044
[25] Tochitsky S. Y. Narang R. Filip C. V. Phys. Rev. Lett. 92 2004 095004
[26] Filip C. V. Narang R. Tochitsky S. Ya. Phys. Rev. E 69 2004 026404
[27] Walton B. Najmudin Z. Wei M. S. Phys. Plasmas 13 2006 013103
[28] Kalmykov S. Yi S. A. Shvets G. Plasma Phys. Control. Fusion 51 2009 024011
[29] Gupta D. N. Singh M. Suk H. J. Plasma Phys. 81 2015 905810324
[30] Singh A. Gupta N. Phys. Plasmas 22 2015 062115
[31] Sodha M. S. Govind Tewari D. P. J. Appl. Phys. 50 1979 158
[32] Mahmoud S. T. Pandey H. D. Sharma R. P. J. Plasma Phys. 69 2003 45
[33] Chauhan P. Mahmoud S. T. Sharma R. P. Pandey H. D. J. Plasma Phys. 73 2007 117
[34] Purohit G. Chauhan P. K. Sharma R. P. Phys. Scr. 77 2008 065503
[35] Gupta M. K. Sharma R. P. Gupta V. L. Phys. Plasmas 12 2005 123101
[36] Sharma R. P. Chauhan P. K. Phys. Plasmas 15 2008 063103
[37] Rawat P. Singh R. K. Sharma R. P. Purohit G. Eur. Phys. J. D 68 2014 57
[38] Yin J. Gao W. Zhu Y. Progress in Optics 44 2003 119
[39] Yuan Y. Cai Y. Qu J. Opt. Express 17 2009 17344
[40] Xu X. Wang Y. Jhe W. J. Opt. Soc. Am. B 17 2002 1039
[41] Cai Y. Lu X. Lin Q. Opt. Lett. 28 2003 1084
[42] Cai Y. Lin Q. J. Opt. Soc. Am. A 21 2004 1058
[43] Mei Z. Zhao D. J. Opt. Soc. Am. A 22 2005 1898
[44] Sodha M. S. Mishra S. K. Misra S. J. Plasma Phys. 75 2009 731
[45] Chaitanya N. A. Jabir M. V. Banerji J. Samanta G. K. Sci. Report 6 2016 32464
[46] Misra S. Mishra S. K. Prog. Electromagnetics Res. B 16 2009 291
[47] Herman R. M. Wiggins T. A. J. Opt. Soc. Am. A 8 1991 932
[48] Wang X. Littman M. G. Opt. Lett. 1 1993 767
[49] Paterson C. Smith R. Opt. Commun. 124 1996 121
[50] Zhao C. Cai Y. Wang F. Opt. Lett. 33 2008 1389
[51] Borisov A. B. Borovskiy A. V. Shiryaev O. B. Phys. Rev. A 45 1992 5830
[52] Brandi H. S. Manus C. Mainfray G. Lehner T. Phys. Rev. E 47 1993 3780
[53] Gupta R. Sharma P. Rafat M. Sharma R. P. Laser part. Beams 29 2011 227
[54] Sharma P. Bhardwaj A. K. Sharma R. P. J. Phys.: Conf. Ser. 534 2014 012049
[55] Sharma P. Laser Part. Beams 33 2015 755
[56] Mendonca J. T. Vieira J. Phys. Plasmas 21 2014 033107
[57] Vieira J. Mendonca J. T. Phys. Rev. Lett. 112 2014 215001
[58] Akhmanov A. S. Sukhorukov A. P. Khokhlov R. V. Sov. Phys. Usp. 10 1968 609
[59] Sodha M. S. Tripathi V. K. Ghatak A. K. Prog. Optics 13 1976 169
[60] Krall J. Ting A. Esarey E. Sprangle P. Phys. Rev. E 48 1993 2157