Discrete & Continuous Models 2019 27 (2) 133—142
. „ . http://journals.rudn.ru/miph
& Applied Computational Science
Research article
UDC 533.9
DOI: 10.22363/2658-4670-2019-27-2-133-142
On the radiation losses during motion of an electron in the field of intense laser radiation
Ekaterina V. Dobrova, Vladimir P. Milantiev
Institute of Physical Research and Technology Peoples' Friendship University of Russia (RUDN university) 6, Miklukho-Maklaya St., Moscow, 117198, Russian Federation
(received: April 20, 2019; accepted: November 14, 2019)
Motion of the relativistic electron in the field of intense laser pulse of the arbitrary shape is considered. The pulse dimension is supposed to be of the order of the Gaussian laser beam dimension in the focal plane. It is supposed that the pulse is propagating along the external constant magnetic field. In the paraxial approximation the corrections of the first order to the vectors of the field of radiation as well as the force of the radiation friction are taken into account. Averaged relativistic equations of motion of electron are obtained with the help of averaging over the fast oscillations of the laser radiation. It is shown that with taking into account corrections of the first order to the field vectors an averaged force arises. This force is defined by pulsed character of radiation and proportional to the intensity but not to gradient of intensity. It is shown that radiation losses are of little importance in the transverse plane but may considerably act on the longitudinal motion of electron.
Key words and phrases: relativistic electron, intense laser pulse, paraxial approximation, Gaussian beam, radiation friction
1. Introduction
When a relativistic charged particle moves in an electromagnetic field, radiation friction forces can play a significant role [1-3]. Various aspects of the radiation friction forces problem, including the refinement of the classical Lorentz-Abraham-Dirac expression, were discussed in [4-7]. The problem of radiation losses has become particularly relevant in connection with the creation of powerful laser radiation sources [8-10]. The charged particle motion pattern substantially depends on the properties of electromagnetic radiation. In the case of a sufficiently weak radiation intensity, the parameter
g = ——1, which is the ratio of the particle momentum oscillations amplitude
umc
in the wave field to a rest momentum, is small: g < 1. Here, |—| is the amplitude of the electric wave field, u is the wave frequency, e is the charge, m — the mass of the particle, c is the light speed in vacuum. Availability of
© Dobrova E.V., Milantiev V.P., 2019
This work is licensed under a Creative Commons Attribution 4.0 International License http://creativecommons.org/licenses/by/4.0/
a small parameter allows one to find a solution to the equations of motion for a particle in the form of expansions over this parameter and averaging over fast oscillations. When a particle moves in the high-power laser radiation field which is pulsed by its nature and of high intensity, the parameter g becomes large. In the case of an electron, this parameter is represented as
g = 0.85 • 10-9\Vl, (1)
where I[W/cm2] is the laser pulse intensity, and A [^m] is the wavelength. This shows that the parameter g > 1 at an intensity of I > 1018 W/cm2. Intensities of such an order and even higher by their value are well achieved in modern laser devices. In this case, a solution of the equations of motion for a particle in the form of expansions in the parameter (1) is impossible. Therefore, they usually use numerical methods for solving the problem. Analytical solution turns out to be possible if propagation of laser radiation is adequately described within the paraxial approximation in the form of Gaussian beams [11-17]. In this approximation, there is a small parameter
2 a ..
» = Va = 7<<1, (2)
fX/ L-h Asy
where a is the laser beam radius at focus (waist of the laser beam), zr = ka2/2 is the Rayleigh length which determines the diffraction divergence of the laser beam, k = 2^/A is the wave number, and A is the radiation wavelength. The presence of such a small parameter (2) allows finding the solution of the equations of motion for a particle in the form of expansions in this parameter and averaging over fast oscillations of the radiation field. In the case of sufficiently long pulses, the radiation field vectors in the zero approximation of the expansion of Maxwell's equations in parameter (2) are represented in the form of Gaussian beams of different modes as a solution to the parabolic equation [11,14]. Upon that, the longitudinal components of the field vectors in the form of first-order corrections, and the corrections to the transverse components of the field vectors in the form of second-order quantities arise. In the case of sufficiently short pulses, the length of which is of the order of the beam waist, the corrections to the transverse components of the field vectors are first order quantities [13,17]. These amendments were not taken into account in papers [8-10]. Note that in the case of highly-focused laser radiation, paraxial approximation is inapplicable [18].
This work is devoted to the study of a relativistic electron motion in the field of high-power pulsed laser radiation in the form of a Gaussian beam with the circularly polarized basic mode, taking into account first-order corrections to field vectors and radiation friction force.
2. Initial equations and problem statement
In case of laser radiation with the circularly polarized basic mode, as a zero approximation of expansions in a small parameter, the field vectors are described by the formulas [17]:
E°x = e^ cos(^ + 0) = A cos(^ + 0),
VI + z2 (3)
E°u = - №(0) e sin(<^ + 0) = -A sin(<^ +
In the formulas written out: d = t) is the "fast" wave phase,
P2 Z C
^(p, Z) = ^ — arctan(Z) is the slow phase, /(a) is the function which
takes into account the pulsed nature of the radiation, and where the parameter
(i — *)
a = — c , and At is pulse duration. Dimensionless quantities are used for the longitudinal coordinate Z = —, the distance from the Gaussian beam
axis in the transverse plane is p = —, where r = a/x2 + y2. In the first
a
approximation for expansions of the field vectors with respect to the small parameter not only the longitudinal field components, but also corrections to the transverse components arise in the case of sufficiently short pulses [17]:
— 1 =__^__— (0)/' /(1 — p2)2 + Z2e 1+12 x
— 2^cAi(1 + Z 2 )3/^(1 x
x cos (arctan —^^ +—^ ^ — 3 arctan(Z) + d ) = —A cos(^1 + 0), V 1 — p2 1 + Z2 J
— 1 = __— (0)/' /(1 — n2)2 + Z2er+|2 x
y 2^cAi(1 + Z 2 )3/2V(1 x
x sin (arctan —+—^ ^ — 3 arctan(Z) + 0) = A sin(^1 + 0), \ 1 — p2 1 + Z2 J
, 2/p—_ (0) ^ / Zp2 y /r7, „
—1 = J xK ei+z2 sin ( —— arctan — — 3arctan(Z) + Q z A;a(1 + Z2) + ^2 x
(4)
Here /» = .
The relationship between the transverse components of the electric and magnetic vectors is determined by the relations:
= , = . (5)
There is also a longitudinal component of the magnetic field of laser radiation:
= (0) e 1+22 sin ( Zp - arctan - - 3 arctan(Z) + tf) . (6)
2 &a(1 + Z2 ) + ^2 x J
The relativistic electron motion in an electromagnetic field is described by the equations:
d^ = -e(E+—[pH])+fr, (7)
at V mc^f J
p2
where p is a particle momentum vector, m is its mass; 7 = \l 1 +--is
V m2c2
a relativistic factor (dimensionless particle energy), fr is a radiation friction force.
The laser radiation field is given by formulas (3-6). The constant magnetic field is directed along the z axis: H0 = (0,0,H0).
The system of equations (7) must be supplemented by an equation for the wave phase 9 which is considered as a "fast" variable:
^ = ) = —-G. (8)
di V c) 7 w
Here vz is the longitudinal velocity of the particle; and the quantity G has the following meaning
G = j — —. (9)
mc
It follows from equation (8) that the phase d becomes a "slow" variable
V V
with 1--z «0. If 1--z » n, then the phase 9 can be considered as a "fast"
c c
variable over which it is possible to carry out averaging. This is assumed to perform in the future. According to the equations of motion (7), the value G is described by the equation:
" = )*" (10)
This shows that in the first approximation, the quantity G is on average preserved and is determined by the initial conditions. It is usually assumed that G = 1 [1]. From the equation (10) it follows that the quantity G can be represented as an expansion:
G = Go + + .... (11)
Here G0 & 1 is the value averaged over the fast phase, and Gi are quickly oscillating components of the quantity G.
It follows directly from the equations of motion for a particle that the relativistic factor satisfies the equation:
= P(-eE + fr). (12)
di ^f(mc)2
The main task of this work is to obtain averaged equations for electron motion using the Bogolyubov's method [19], by averaging over the fast oscillations of laser radiation, taking into account the averaged radiation friction force. In their sense, the radiation friction forces are small compared
with the effects of laser radiation. Therefore, we first consider the motion of a particle without taking radiation losses into account.
3. Averaged equations of motion for a charged particle
The initial equations of motion for a particle (7) contain members oscillating with the laser frequency, the amplitude of which is determined by parameter (1). In the conditions under consideration, this parameter is large. Therefore, it is impossible to carry out averaging according to the Bo-golyubov's method [19] using expansions in the parameter g. However, the equations of motion also contain parameter (2) which is small. In principle, this allows to seek a solution in the form of expansions in the parameter fi. However, it is necessary that the quickly oscillating members have small amplitude in order to apply the averaging method. This can be achieved if we firstly exclude in the equations of motion the quickly oscillating members with large amplitudes [20,21]. In this regard, we will make the replacement of the transverse components of the particle momentum vector, which is convenient to use in the complex form:
p — n + Ce^+9). (13)
Here, p — px + ipy is the complex transverse momentum of a particle, and K is the complex generalized momentum. Substituting (13) into system (7), we obtain the equation:
dE + ^iW) + + iC^S**» —
di di di di
- e(l- ^) Aei(*+e) + ^(tt + Cei(v+e)) (l + ^)■ (14) V c ) j V HoJ
In the equation (14) the first member on the right is quickly oscillating with large amplitude. It can be eliminated by appropriate selection of the amplitude C in the transformation (13). In its meaning, the amplitude C must be determined by the amplitude of the radiation field in zero approximation. Then, equating the members of zero approximation in the equation (14)
considering equation (8) and expansion (11), we get:
c' =^ (15)
eH
Here
ujr,r, —-is a
classic cyclotron frequency. Thus, for the found value
mc
of C, equation (14) for the variable n contains only members of the first order in the parameter ¡i. We can obtain the following expression from the definition of the relativistic factor, taking into account (9)
1 + G2 +
jL
i =-^^, (16)
where p2 = p2 Given the replacement (13), we get
= + c 2 _ sin(^ + 0-£).
Here: = + , ^ = ^ cos ff, = sin ff.
(17)
The value of G =
represents the amplitude of electron mo-
mentum oscillations in a field of laser radiation propagating along a constant magnetic field. Amplitude A is determined by formulas (3). After a series of transformations, formula (16) takes the form:
1 Gr
7 =
+ it +
2G0 ' 2 ' 2G0m2c2 ' 2
G2 Gi G^ sin(^ + в _ ff) С2G!
+ —---^-—---(18)
G0 m2 с2
2G0m2 с2'
1 1 / G \ Here the expansion —-— « — ( 1 — —1 ) is taken into account. It is
G0_+ Gx G0 \ G0/
also assumed that < C. This is a natural limitation in particle acceleration problems. From the formula (18) it follows that the relativistic factor contains "constant" and fast oscillating parts:
7 = r + 7i,
where
Г =
1 + G2 +
с2
Gi
2Gn
7i =
G2G1 G^ sin(p + 0_0)
2G2 m2c2
G0 m2 c2
(19)
(20)
Let us further consider the equation of longitudinal motion from system (7):
^ = e£i + — di г +
mcy
Pa
(0)-e 1+22 cos(^ + 0) + Я1
л/Г+Z2
' (0) .л/ГТ^2
After a series of transformations we get: dp
e i+z2 sin(^ + Q) + H
dT = eRïf§)e sin _ ^ctan(z) _ f + +
2e
m2c2 + G0m2 c2 + G2
G0 _
G0G1 m2c2 _ G2G1 _ 2G0G^ sin(^ + 0 _ ff) m2 c2 + G2 m2c2 + G2
x
x
cos(^ + в + ff) _ ^ cos(^1 + в +
z
_ AG sin + 0) + G sin (arctan --^ _ 2 arctan(Z)) +
22
+
+ A1 Csin
2p2z , Z
— 4 arctan(z) + arctan
1 + Z 2
1 — P2
+ 29
(21)
We omit the members containing , and select the members that are independent on the fast phase Q. As a result, we obtain the averaged equation of the longitudinal motion:
dPz dt
e2 X0Z
E2X (0)ff '
1
V(i — p2 )2 + z 2 e ^ x z
2mnc2 AtT (1 + Z2)2 u + uco/G0
x sin ( arctan
1 — P2
— 2 arctan(Z)) . (22)
Here pz is the average value for the longitudinal momentum of the particle. The right side of this equation, which has the first order of smallness in the parameter describes the averaged force with which high-power pulsed laser radiation acts on an electron in the direction of the magnetic field. A similar expression was obtained in [15] without taking into account the external magnetic field. However, due to the large difference between the cyclotron frequency and the laser radiation frequency, the influence of an external magnetic field is insignificant. As can be seen from (22), the averaged force is
df
due to the pulsed nature of the radiation: f ' = — + 0.
oa
4. Radiation friction force
The radiation friction force problem is one of the frequently discussed in the literature [3-10]. In [6,7], some shortcomings were noted and refinements of the classical expression for the radiation friction force were made [1,2]. In the present work, the classical Lorenz-Abraham-Dirac formula is used for the radiation friction force [1-3]:
ffr =
2e4 i ^ ^ 1 ' [ÊH] + -c
3m2 c
H[Hv]
2e4-f 3mc3
1 - \2 1 -E+-[vH]) -- (Ev)2
(23)
For the longitudinal motion, the radiation friction force when electron is moving in the laser field (3), taking into account the transformation (13), is represented as;
fz = { A2 cos2 (<p + 6) +A2 sin2 (<p + 6) + ^°-
3m2 c4 y cmj
—AC sin2 (^ + 9)+
+ AC cos2 ((f + 2e41
Pz
mcy
cos2 (if + d) + sin2 (if + 2
pz {A2 cos2 (v + d) + A2 sin2 (v + ff)--ACH0 sin2 (^ + 6)—
3m2 c3 I mcj
— ACH0 cos2(^ + 0)+ pz(^2 cos2(^ + d) + A2 sin2 (^ + 0))
c2 m27 2
A2C2 sin4 (p + 0) + A2C2 cos4 (p + 0) - 2AC#0^ sin2 (p + 0)+ + ^ cos2(^ + 0) + 2A2C2 sin2 (^ + 0) cos2 (p + 6
Averaging over the fast oscillations of the radiation field and holding only the main members, we obtain the expression for the radiation friction force in the direction of the laser pulse propagation:
f o4 /2£2(0) f^/i Q Z2(0)e4r
= 26 3m2c4(1 + Z2)61+22 {1 - mcT} - 23m2c3(1 + Z2)6^
f 2pz e2/2(0) r2p2l , ^
X <1--—-----^-e i+z2 } (24)
\ mer m2 c2 r2(1 + Z2 )(w- ^ )2 C J' ( )
Expression (24) is still quite complicated. For a rough estimate, we will assume that the particle velocity is close enough to the velocity of light. Then
the main member of this expression is--. From here, it can be
3m2 c3
concluded that the radiation reaction force can be sufficiently large in the field of high-intensity laser radiation. The calculations showed that the averaged transverse components of the radiation friction force are absent:
f =0 f =0
Javx ' Javy
This means that when a particle is moving in the transverse plane, the radiation losses on average play no role. When a particle is moving in the direction of radiation propagation, they can be significant in the field of super-intense laser radiation.
5. Conclusion
An expression for the averaged force acting on a relativistic electron in the field of an intense short laser pulse propagating along an external constant magnetic field is obtained. This force is due to the pulsed nature of the radiation, the description of which requires the inclusion of first-order members in expansions with respect to a small parameter of the paraxial approximation. It was taken into account that in the case of high-power laser radiation, the ratio of the amplitude of the particle oscillatory velocity in a wave field to the velocity of light can be large. To use the averaging over the fast oscillations of laser radiation, the transverse components of the particle momentum vector have been transformed. Thus, in the equations of motion, quickly oscillating members with large amplitude determined by the radiation field were eliminated. An expression was obtained for the averaged radiation friction force determined by the classical Lorentz-Abraham-Dirac formula. It was shown that this force in the transverse plane is zero and can be significant in the direction of laser pulse propagation.
Acknowledgments
The publication has been prepared with the support of the "RUDN University Program 5-100".
References
1. L. D. Landau, E. M. Lifshitz, Field theory [Teoriya polya], Nauka, Moscow, 1988, in Russian.
2. J. D. Jackson, Classical electrodynamics, J. Wiley, NY-L., 1962.
3. V. L. Ginzburg, Theoretical physics and astrophysics [Teoreticheskaya fizika i astrofizika], Nauka, Moscow, 1975, in Russian.
4. N. P. Klepikov, Radiation damping forces and radiation from charged particles, Soviet Physics Uspekhi 28 (1985) 506-520. doi:10.1070/PU1985v028n06ABEH005205.
5. V. S. Krivitskii, V. N. Tsytovich, Average radiation-reaction force in quantum electrodynamics, Soviet Physics Uspekhi 234 (1991) 250-258. doi:10.1070/PU1991v034n03ABEH002352.
6. G. F. Efremov, Radiative damping of a relativistic electron in classical electrodynamics, Journal of Experimental and Theoretical Physics 89 (1998) 899-904. doi:10.1134/1.558738.
7. I. V. Sokolov, Renormalization of the Lorentz-Abraham-Dirac equation for radiation reaction force in classical electrodynamics, Journal of Experimental and Theoretical Physics 109 (2009) 207-212. doi:10.1134/S1063776109080044.
8. A. L. Galkin, Dynamics of an electron in a relativistically intense laser field including radiation reaction, Journal of Experimental and Theoretical Physics 115 (2012) 201-207. doi:10.1134/S1063776112070072.
9. A. V. Bashinov, A. A. Gonoskov, A. V. Kim, M. Marklund, G. Mourou, A. M. Sergeev, Electron acceleration and emission in a field of a plane and converging dipole wave of relativistic amplitudes with the radiation reaction force taken into account, Quantum Electronics 43 (2013) 291-299. doi:10.1070/QE2013v043n04ABEH015101.
10. K. Seto, H. Nagamoto, J. Koga, K. Mima, Equations of motion with radiation reaction in ultrarelativistic laser-electron interactions, Physics of Plasmas 18 (2011) 123404. doi:10.1063/1.3663843.
11. L. W. Davis, Theory of electromagnetic beams, Physical Review A 19 (1979) 1177-1179. doi:10.1103/PhysRevA.19.1177.
12. D. Bauer, P. Mulser, W. Steeb, Relativistic ponderomotive force, uphill acceleration and transition to chaos, Physical Review Letters 75 (1995) 4622-4625. doi:10.1103/PhysRevLett.75.4622.
13. B. Quesnel, P. Mora, Theory and simulation of the interaction of ultraintense laser pulses with electrons in vacuum, Physical Review E 58 (1998) 3719-3732. doi:10.1103/PhysRevE.58.3719.
14. A. M. Goncharenko, Gaussian light beams [Gaussovy puchki sveta], Nauka i tekhnika, Minsk, 1997, in Russian.
15. D. R. Bituk, M. V. Fedorov, Relativistic ponderomotive forces, Journal of Experimental and Theoretical Physics 89 (1999) 640-646.
16. N. B. Narozhny, M. S. Fofanov, Scattering of relativistic electrons by a focused laser pulse, Journal of Experimental and Theoretical Physics 90 (2000) 753-768. doi:10.1134/1.559160.
17. V. P. Milant'ev, S. P. Karnilovich, Y. N. Shaar, Description of high-power laser radiation in the paraxial approximation, Quantum Electronics 45 (2015) 1063-1068. doi:10.1070/QE2015v045n11ABEH015800.
18. S. G. Bochkarev, V. Y. Bychenkov, Acceleration of electrons by tightly focused femtosecond laser pulses, Quantum Electronics 37 (2007) 273-284. doi:10.1070/QE2007v037n03ABEH013462.
19. N. N. Bogoljubov, Y. A. Mitropolskij, Asymptotic methods in the theory of nonlinear oscillations [Asimptoticheskiye metody v teorii nelineynykh kolebaniy], Nauka, Moscow, 1974, in Russian.
20. V. P. Milant'ev, A. J. Castillo, On the theory of the relativistic motion of a charged particle in the field of intense electromagnetic radiation, Journal of Experimental and Theoretical Physics 116 (2013) 558-566. doi:10.1134/S1063776113040067.
21. A. J. Castillo, V. P. Milant'ev, Relativistic ponderomotive forces in the field of intense laser radiation, Technical Physics 59 (2014) 1261-1266. doi:10.1134/S1063784214090138.
For citation:
E. V. Dobrova, V. P. Milantiev, On the radiation losses during motion of an electron in the field of intense laser radiation, Discrete and Continuous Models and Applied Computational Science 27 (2) (2019) 133-142. DOI: 10.22363/26584670-2019-27-2-133-142.
Information about the authors:
Ekaterina V. Dobrova (Russian Federation) — student of Institute of Physical Research and Technology of Peoples' Friendship University of Russia (RUDN University) (e-mail: dobrova03@icloud.com, phone: +7(499) 2483057)
Vladimir P. Milantiev (Russian Federation) — Professor, Doctor of Physical and Mathematical Sciences, professor of Institute of Physical Research and Technology of Peoples' Friendship University of Russia (RUDN University) (e-mail: vmilant@mail.ru, phone: +7(499) 2483057, Scopus Author ID: 6506148930)