Discrete & Continuous Models 2021 29 (2) 105—113
. „ . http://journals.rudn.ru/miph
& Applied Computational Science
Research article UDC 533.9
PACS 42.62.-b, 52.27.Ng, 52.35.Mw DOI: 10.22363/2658-4670-2021-29-2-105-113
On the possibility of averaging the equations of an electron motion in the intense laser radiation
Vladimir P. Milant'ev
Peoples' Friendship University of Russia (RUDN University) 6, Miklukho-Maklaya St., Moscow 117198, Russian Federation
(received: March 12, 2021; accepted: May 25, 2021)
The problem of averaging of the relativistic motion equations of electron in the intense laser radiation, caused by the decreasing of the rate of wave phase change due to the Doppler's effect, is considered. As a result the phase can go from the "fast" to "slow" variables of the motion, so averaging over the phase becomes impossible. An analysis is presented of the conditions which are necessary for averaging of the relativistic equations of motion over the "fast" phase of the intense laser radiation on the base of the general principles of the averaging method. Laser radiation is considered in the paraxial approximation, where the ratio of the laser beam waist to the Rayleigh length is accepted as a small parameter. It is supposed that the laser pulse duration is of the order if the laser beam waist. In this case first-order corrections to the vectors of the laser pulse field should be taken into account. The general criterion for the possibility of the averaging of the relativistic motion equations of electron in the intense laser radiation is obtained. It is shown that an averaged description of the relativistic motion of an electron is possible in the case of a fairly moderate (relativistic) intensity and relatively wide laser beams. The known in the literature analogical criterion has been obtained earlier on the base of the numerical results.
Key words and phrases: intense laser pulse, relativistic electron, equations of motion, averaging of equations, criterion for averaged description of motion
1. Introduction
The nature of the motion of electrons in the field of electromagnetic waves substantially depends on the wave intensity, which is characterized by the di-mensionless parameter g = eE/umec. Here E is the electric field amplitude of the wave, u is its angular frequency, e and me are the electron charge and mass, respectively, c is the velocity of light in vacuum. The first papers [1], [2] were devoted to the nonrelativistic motion of an electron in a high-frequency electromagnetic field of low intensity (parameter g < 1). It was shown by averaging over fast field oscillations, and expansions in terms of the parameter g,
© Milant'ev V.P., 2021
This work is licensed under a Creative Commons Attribution 4.0 International License http://creativecommons.org/licenses/by/4.0/
#
that the particle was subjected to the action of an averaged (ponderomo-tive) force. Later, the relativistic generalization of the ponderomotive force was considered under the condition that the parameter g was small [3], [4]. It was also noted that relativistic effects lead to various features of the averaged force [4]. In the field of high-power laser radiation, the parameter g is large (g > 1). So expansions in terms of the parameter g become impossible. In the case of electrons, the parameter g = 1, when the electric field strength Er (V/cm) = mecu/e = 3.21 ■ 1010/X (^m), where A (^m) is the wavelength. Radiation with electric field strength E > Er is called rel-ativistically strong [5]. The parameter g is commonly represented in the form:
g = 0.855 ■ 10-g\Vl, (1)
where I = (cE2/8k) [W/cm2] is the intensity of the laser pulse. The parameter g is small in the case of a relatively weak field, when I < Ir. Here Ir = m"2c3u2/8^e2 is the relativistic intensity determined by the electric field strength Er. Intensity of modern lasers can reach I > 1018 W/cm2 [6]-[8].
In the study of particle motion, an adequate description of the laser radiation field plays an important role. When describing laser radiation, the paraxial approximation and its modifications is often used [9]-[13] which are based on the expansion of field vectors in terms of a small parameter
p = a/ZR = 2/ka<1. (2)
Here a is the size of the laser beam in focus (beam waist), ZR = ka2/2 is the Rayleigh length, k = 2n/\ = u/c is the wave number. We assume that the laser field propagates in the ^-direction. From the Maxwell equations, one can find the expressions for the transverse components of the radiation field vectors E°m, B°m of the zero approximation in the form of Gaussian beams of various modes m [9]-[13]. Longitudinal components Elm, Blm also arise, which are of the first-order quantities. The parameter (2) establishes the relation between the wavelength of radiation and the size of the focal spot. Powerful laser radiation also has a characteristic scale — the length (or duration At) of the pulse. In the case of extended pulses, the corrections to the transverse components of the radiation field are second-order quantities [9]. If the pulse length cAt and the size of the focal spot a are of the same order cAt ~ a, then the first-order corrections to the transverse components of the field vectors E\m, Bj_m appear [10]-[13]. In this case, the pulsed character of the radiation is specified by a fairly smooth pulse function f(a), where the parameter a = (t — z/c)/At.
In the case of tightly focused laser radiation with the intensity I > 1022 W/cm2, the size of the focal spot can be equal to or smaller than the wavelength. In this case, the parameter (2) is not small, so that the paraxial approximation is not applicable and an exact solution of the Maxwell equations is necessary [14].
The presence of a small parameter (2) in the equations of the electron motion allows us to use the perturbation theory and perform averaging over fast oscillations of radiation. When they derive the ponderomotive force of a laser pulse, it is usually assumed that the wave amplitude varies slowly
with respect to the wave phase (for example [15], [16]). However, the specific conditions for the relative change of these parameters are not considered. Meanwhile, the absence of such an analysis can lead to the misuse of averaging of the equations of motion. The fact is that during relativistic motion, the frequency of the radiation that the particle "sees" decreases due to the Doppler shift: u' = u(1 — vz/c). Here vz is the component of the particle velocity in the direction of the laser pulse propagation. Doppler frequency shift slows down the rate of wave phase change. Therefore, at a sufficiently high longitudinal velocity of the particle, the rate of phase change may turn out to be comparable with the change in the wave amplitude. This problem was partially touched upon in the paper [4]. However, the conditions under which the averaging of the equations of motion is permissible were not discussed in detail. It was verified in the work [10] by numerical calculations that the domain of validity of the averaged description of the electron motion in the ultraintense laser pulse was given by the condition 1 — vz/c » e, where e = ^/2. However, the meaning of this condition and its validity was not discussed.
This paper is devoted to the detailed analysis of the conditions for averaging the relativistic equations of the electron motion in the field of high-power laser radiation with a sufficiently long pulse duration such that A < cAt ~ a. In this case, the existence of a small parameter (2) is assumed as in the work [10].
2. Basic relations
The motion of an electron is described by the following equations:
dpa
= —(1 — vz/c)eExm — eBlmPy/me c-í, + PR
dt ^
dt dp y
= —(1 — Vz/c)eE + eBlm px/meC^, (3)
dpz ^
e(Px ^xm + Py ^yrn )/me ,
dr p
dt me 7 d^f e
(4)
pE. (5)
dt (me c)2^
Here p = (px ,py,pz) is the electron momentum vector, 7 is the relativistic factor (dimensionless energy). The vectors of the laser field Ern, Bm of an arbitrary mode m taking into account first-order terms are determined by the formulas [13] (or [10]). So Em = E°rn + E^ and Bm = B°rn + B^. Along with the equations (3), it is necessary to use the equation for the wave phase d:
f = -u(1 - vz/c) = -uG/1. (6)
Here
G = ^-pz fme c. (7)
Let's note that the phase that the particle "senses" in the laser field differs from the phase d by additional small terms [9]-[13]. However, in the case under consideration, these terms are not significant. The equations of motion (3) in the field of high-power laser radiation are very complicated for an analytical solution. Therefore, numerical methods of solution are often used that allow one to study some features of an electron motion in a laser field [10], [17]-[19]. In this case, most often, laser radiation is specified in the form of a Gaussian beam of the fundamental mode, even at g ^ 1, which, in general, is incorrect due to the following reasons: Solution of the Maxwell equations in the form of Gaussian (or Hermite-Gaussian) laser beams is the result of expansion of the field strength over the parameter (1). In the case of the ultra-intense and ultra-short laser pulses the size of the focal spot can be comparable with the wavelength [14]. So, the relation (1) is violated and description of the laser radiation in the form of the Gaussian beams becomes invalid. In this case exact solution of the Maxwell equations should be found [14].
3. Conditions for relativistic equations of motion
averaging
A simplified description of electron interaction with a laser is achieved by averaging the equations of motion over the wave phase. Various versions of the averaged equations of motion have been considered in many works [20]-[24].
To average equations (3) over the phase d, it must be a "rapidly" changing quantity [25]. It follows from the equation (6) that this depends on the difference 1 — vz/c = A, where 0 < A < 1 (if the particle moves in the direction of wave propagation). At A ~ 1 the phase changes "quickly" and averaging over the phase is possible. The last is a general condition for averaging the equations of motion. In the ultrarelativistic limit (A < 1) the phase d becomes a "slow" (or "semi-fast") variable as well as the wave amplitude. In this case, the electron motion changes significantly and becomes more complicated [10]. In the presence of a unique small parameter (2), it is quite natural to present the above general averaging criterion in the modified form [10]:
A~1» (8)
It follows from (6), that the difference A = G/^y, where the quantity G, according to equations (3), satisfies the equation:
dG
— = e(1 — Vz/c)E1xm/me c+.... (9)
One can see that the quantity G is an integral of motion only in the case of a plane electromagnetic wave in vacuum (E^m =0). In this case, the value G is determined by the initial conditions: G = 7(0) — pz(0)/mec. If one considers the particles at rest at the initial instant of time, then G = 1. In the case of laser radiation, the longitudinal field E'^m always exists and plays
an essential role in the motion of electrons. So, in general, the quantity G contains a slowly changing part as well as quickly oscillating corrections with small amplitudes. However, it is sometimes believed that G = 1 also in the case of laser radiation [11].
Let us further consider the averaging condition (8). In this case, the value G can be represented as the expansion in terms of the parameter
G = G0 + G1 + ..., (10)
where G0 (the averaged value of the quantity G) does not depend on the wave phase 9, while Gi are periodic functions. According to equation (9), the quantity G0 remains constant up to the first-order terms. It follows from the condition (8) that averaging of the equations of motion is possible if the following inequality is fulfilled:
JV«Go -1 (11)
or
^«1/^ = ^a/X. (12)
This condition must be satisfied both during the injection of particles in the radiation field, and during their further movement. It follows from (11), (12) that an averaged description of motion is allowed when the energy of an accelerating particle is limited. Typically [10] the parameter ^ < 6.4- 10-2. So the energy of the particle is restricted by the condition ^y « 17.
Let us consider the relativistic factor 7 = ^1 + p2/(mc)2. With the definition (7), it is easy to obtain the following equation:
1 = [1 + G2 +pl/(mec)2 ]/2G. (13)
Here p\ = pX + p2. It follows from the system of equations (3) that p± — gmec. Then from (13), we obtain the following estimate: 7 — 1 + g2/2. Given the inequality (12), we conclude that the averaging of electron motion equations in the field of relativistically intense laser radiation is possible if a rather stringent condition is satisfied:
1 + g2/2 « na/\. (14)
Thus, the averaging of the equations of motion is possible in the case of fairly moderate intensity of laser radiation and a relatively wide laser beam (a/X » 1). In the case of ultra-intense radiation (g » 1), as it was already noted, the wavelength of the laser beam may be comparable with its size in the focus. Then the parameter in (2) turns out to be large, and expansion in the terms of this parameter becomes impossible. Moreover, at g » 1 the difference 1 — vz/c = [1 + p]_/(mec)2}/2^2 — g-2 « 1, and motion of an electron becomes very complicated as it was noted in the paper [10]. That means that the concept of the relativistic ponderomotive force has rather restricted domain of validity.
4. Conclusion
It is shown that the condition 1 — vz/c » ^/2, obtained by computer calculations in the paper [10], really corresponds to the general criterion (8) for averaging of the classical relativistic equations of electron motion in the intense laser beam. It leads to the conclusion that averaged description of the relativistic electron motion is possible at limited electron energy and limited intensity of the laser radiation, as it is established by the inequalities (12), (14). So averaging the equations of an electron motion over the wave phase seems to be possible in the case of a fairly moderate intensity and a relatively wide laser beam. Therefore, in general, it is impossible to consider the problem of the ponderomotive acceleration of electrons at very high intensity of the laser radiation. It should be particularly emphasized that for the averaging procedure it is necessary to take into account not only intensity but also other characteristics of the laser pulse.
Acknowledgement
The reported study was funded by RFBR, project No. 18-29-21041.
References
[1] H. Boot and R.-S. Harvie, "Charged particles in a non-uniform radiofre-quency field," Nature, vol. 180, p. 1187, 1957. DOI: 10.1038/1801187a0.
[2] A. Gaponov and M. Miller, "Use of moving high-frequency potential wells for the acceleration of charged particles," Soviet Physics JETP-USSR, vol. 7, pp. 515-516, 1958.
[3] T. Kibble, "Mutual refraction of electrons and photons," Physical Review, vol. 150, p. 1060, 1966. DOI: 10.1103/PhysRev.150.1060.
[4] D. R. Bituk and M. V. Fedorov, "Relativistic ponderomotive forces," JETP, vol. 89, pp. 640-646, 4 1999. DOI: 10.1134/1.559024.
[5] G. A. Mourou, T. Tajima, and S. V. Bulanov, "Optics in the relativistic regime," Reviews of Modern Physics, vol. 78, p. 309, 2006. DOI: 10.1103/ RevModPhys.78.309.
[6] N. M. Naumova, J. A. Nees, and G. A. Mourou, "Relativistic attosecond physics," Physics of Plasmas, vol. 12, p. 056 707, 2005. DOI: 10.1063/1. 1880032.
[7] A. V. Korzhimanov et al., "Horizons of petawatt laser technology," Physics-Uspekhi, vol. 54, p. 9, 2011. DOI: 10.3367/UFNe.0181.201101c. 0009.
[8] C. Danson, D. Hillier, N. Hopps, and D. Neely, "Petawatt class lasers worldwide," High Power Laser Science and Engineering, vol. 3, 2015. DOI: 10.1017/hpl.2019.36.
[9] L. W. Davis, "Theory of electromagnetic beams," Physical Review A, vol. 19, pp. 1177-1179, 3 1979. DOI: 10.1103/PhysRevA.19.1177.
[10] B. Quesnel and P. Mora, "Theory and simulation of the interaction of ultraintense laser pulses with electrons in vacuum," Physical Review E, vol. 58, pp. 3719-3732, 3 1998. DOI: 10.1103/PhysRevE.58.3719.
[11] G. V. Stupakov and M. Zolotorev, "Ponderomotive laser acceleration and focusing in vacuum for generation of attosecond electron bunches," Physical Review Letters, vol. 86, p. 5274, 2001. DOI: 10.1103/PhysRevLett. 86.5274.
[12] W. Wang, J. Xia L.and Xiong, H. Fang Z.and An, Z. Xie, W. Pei, and S. Fu, "Field shaping and electron acceleration by center-depressed laser beams," Physics of Plasmas, vol. 26, p. 093109, 2019. DOI: 10.1063/1. 5099508.
[13] V. P. Milant'ev, S. P. Karnilovich, and Y. N. Shaar, "Description of high-power laser radiation in the paraxial approximation," Quantum Electronics, vol. 45, pp. 1063-1068, 11 2015. DOI: 10. 1070/ QE2015V045N11ABEH015800.
[14] S. G. Bochkarev and V. Y. Bychenkov, "Acceleration of electrons by tightly focused femtosecond laser pulses," Quantum Electronics, vol. 37, pp. 273-284, 3 2007. DOI: 10.1070/QE2007v037n03ABEH013462.
[15] E. Startsev and C. McKinstrie, "Multiple scale derivation of the rela-tivistic ponderomotive force," Physical Review E, vol. 55, p. 7527, 1997. DOI: 10.1103/PhysRevE.55.7527.
[16] P. Mora and J. T. M. Antonsen, "Kinetic modeling of intense, short laser pulses propagating in tenuous plasmas," Physics of Plasmas, vol. 4, pp. 217-229, 1997. DOI: 10.1063/1.872134.
[17] Y. I. Salamin, G. R. Mocken, and C. H. Keitel, "Electron scattering and acceleration by a tightly focused laser beam," Physical Review Special Topics-Accelerators and Beams, vol. 5, p. 101301, 2002. DOI: 10.1103/PhysRevSTAB.5.101301.
[18] P. Wang et al., "Characteristics of laser-driven electron acceleration in vacuum," Journal of Applied Physics, vol. 91, pp. 856-866, 2002. DOI: 10.1063/1.1423394.
[19] A. Galkin, V. Korobkin, M. Y. Romanovsky, and O. Shiryaev, "Dynamics of an electron driven by relativistically intense laser radiation," Physics of Plasmas, vol. 15, p. 023104, 2008. DOI: 10.1063/1.2839349.
[20] D. Bauer, P. Mulser, and W. H. Steeb, "Relativistic ponderomotive force, uphill acceleration and transition to chaos," Physical Review Letters, vol. 75, pp. 4622-4625, 25 1995. DOI: 10.1103/PhysRevLett.75.4622.
[21] I. Y. Dodin, N. J. Fisch, and G. M. Fraiman, "Drift lagrangian for a relativistic particle in an intense laser field," Journal of Experimental and Theoretical Physics Letters, vol. 78, pp. 202-206, 2003. DOI: 10. 1134/1.1622032.
[22] N. B. Narozhny and M. S. Fofanov, "Scattering of relativistic electrons by a focused laser pulse," JETP, vol. 90, pp. 753-768, 5 2000. DOI: 10.1134/1.559160.
[23] A. J. Castillo and V. P. Milant'ev, "Relativistic ponderomotive forces in the field of intense laser radiation," Technical Physics, vol. 59, pp. 12611266, 9 2014. DOI: 10.1134/S1063784214090138.
[24] V. P. Milant'ev and A. J. Castillo, "On the theory of the relativistic motion of a charged particle in the field of intense electromagnetic radiation," JETP, vol. 116, pp. 558-566, 4 2013. DOI: 10. 1134/ S1063776113040067.
[25] N. N. Bogoljubov and Y. A. Mitropolskij, Asymptotic methods in the theory of nonlinear oscillations. New York: CRC Press, 1961, vol. 10, p. 537. DOI: 10.1007/BF01056172.
For citation:
V. P. Milant'ev, On the possibility of averaging the equations of an electron motion in the intense laser radiation, Discrete and Continuous Models and Applied Computational Science 29 (2) (2021) 105-113. DOI: 10.22363/26584670-2021-29-2-105-113.
Information about the authors:
Vladimir P. Milant'ev — Doctor of Physical and Mathematical Sciences, Professor of Institute of Physical Research and Technology of Peoples' Friendship University of Russia (RUDN University) (e-mail: milantiev-vp@rudn.ru, phone: +7(499) 2483057, ORCID: https://orcid.org/0000-0003-4686-4229, ResearcherlD: B-9335-2016, Scopus Author ID: 6602468700)
УДК 533.9
PACS 42.62.-b, 52.27.Ng, 52.35.Mw DOI: 10.22363/2658-4670-2021-29-2-105-113
О возможности усреднения релятивистских уравнений движения электрона в поле мощного лазерного
излучения
В. П. Милантьев
Российский университет дружбы народов ул. Миклухо-Маклая, д. 6, Москва, 117198, Россия
Рассмотрена проблема усреднения релятивистских уравнений движения электрона в поле мощного лазерного излучения, вызванная уменьшением скорости изменения фазы волны из-за эффекта Доплера. Вследствие этого фаза может перейти из числа «быстрых» в число «медленных» переменных движения, так что усреднение по фазе становится невозможным. На основе общих принципов метода усреднения проведён анализ условий, при которых допустимо усреднение уравнений движения по «быстрой» фазе излучения. Лазерное излучение рассматривается в параксиальном приближении, в котором малым параметром является отношение сужения лазерного пучка к рэлеевской длине. Предполагается, что протяжённость импульса сопоставима с порядком сужения лазерного пучка. В этом случае необходимо учитывать поправки первого порядка к векторам поля лазерного импульса. Получен общий критерий, определяющий возможность усреднения релятивистских уравнений движения частицы в поле мощного лазерного излучения. Показано, что усреднённое описание релятивистского движения электрона возможно в случае достаточно умеренной (релятивистской) интенсивности и относительно широких лазерных пучков. Известный в литературе аналогичный критерий был получен ранее на основе численных расчётов.
Ключевые слова: мощный лазерный импульс, релятивистский электрон, уравнения движения, усреднение уравнений, критерий для усреднённого описания движения