Electron Dynamics and Thomson Scattering for Ultra-Intense Lasers: Elliptically Polarized and OAM Beams

We investigated the classical nonlinear Thomson scattering (TS), from a single relativistic electron, generated by either: (a) an incoming plane wave monochromatic laser radiation and general elliptical polarization or (b) incoming radiations with intrinsic orbital angular momentum (OAM). Both (a) and (b) propagate along the z direction, with wave vector k0, frequency ω0, and initial phase φ0≠0 and have any intensity. Item (a) enables obtaining general electron TS Doppler frequencies and other quantities, for fusion plasmas. We explored the possibility of approximating nonlinear TS with OAM beams (Item (b)) by means of nonlinear TS with plane wave beams (Item (a)). For Item (a), a general explicit solution of the Lorentz relativistic equation and the subsequent TS are given in terms of ζ=ω0t−k0z (t denoting time). In particular, it includes the cases for linear and circular polarizations and φ0≠0 for fusion plasmas, thereby extending previous studies for φ0=0. The explicit solutions give rise to very efficient computations of electron TS Doppler frequencies, periods of trajectories, and drift velocities, and the comparisons with ab initio numerical solutions (for Item (a)) yield an excellent match. The approximate approach, using explicit solutions for Item (a), towards TS OAM (employing ab initio numerical computations for Item (b)), extending previously reported ones) yields a quite satisfactory agreement over time spans including several optical cycles, for a wide range of laser intensities, polarizations, and electron energies. The role of φ0≠0 was analyzed. A simple quantitative criterion to predict whether the agreement between the two approaches (a) and (b) would be observed over a given time span is discussed.


Introduction
The scattering of electromagnetic radiation by free charges is a standard diagnostic of their distribution functions (see for example [1,2] for applications in fusion plasmas) and also for the diagnosis of intense or ultra-intense laser beams. The diagnoses in both cases are based on the properties (spectrum, polarization, etc.) of the scattered radiation.
Ultra-intense pulsed lasers have been developed dramatically with the invention of CPA technology [3]. Intensities in the 10 22 W/cm 2 range are now feasible with nearinfrared petawatt or multipetawatt lasers [4]. Such laser intensities ionize and accelerate the released electrons to relativistic speeds, the threshold for that being about 10 18 W/cm 2 . Hence, at 10 22 W/cm 2 , electron motion is highly relativistic. At those limits, the radiation reaction starts playing a role, and the electron dynamics has to be calculated according to that [5]. New realms for light-matter interactions have been opened by these developments that were simply not available at lower laser peak powers.
Laser transverse patterns and/or polarization states can have a very rich structure, the simplest one being the well-known Gaussian (or TEM 00 ) mode. Laser polarization is often associated with the spin angular momentum of the photon. Laser beams may also have orbital angular momentum (OAM); see for example [6], the chapter by Allen and Padgett [7], and also [8][9][10][11][12][13] for some experiments with light in states having OAM. Those beams have remarkable properties because of their structure so different from the ordinary plane wave ones. For example, it is now clear that the photoelectric effect selection rules are valid for OAM photons provided that the photon orbital angular momentum is added in the selection rules [14], so the typical dipole transition rules are violated, and an exchange of more that one unit angular momentum per photon in that process is allowed.
In modern petawatt laser sources, it seems possible to include a variety of transverse structures, having a strong influence on the laser-atom or laser-electron interaction. It seems worthwhile to study a general Laguerre-Gauss mode (with a step, or other, function to describe its time-pulsed structure). Although it is difficult to obtain such OAM pulses at the multi-joule level, it seems very reasonable to consider various possible physical processes generated by the former, one of them being the motion of a fully relativistic electron in such OAM pulses. In this paper, we considered two cases that are of growing interest: amplifying mechanisms leading to petawatt peak powers generate linearly polarized laser fields in most cases. That explains why so far, most of the research has been on such a polarization state or in the other situation of circular polarization; although, at ultrahigh intensities, circular polarization is difficult to achieve. However, in many cases, the mirrors of the beam transport from compressor to target can induce some degree of elliptical polarization (intentional or not). Therefore, the present model describing an arbitrary degree of polarization can be of broad interest to model future experiments. On the other hand, there is a growing interest in OAM beams because of the new possibilities they bring. Until recently, ultra-intense OAM beams were considered very difficult to achieve; however, in a recent paper [15], it was shown that off axis spiral phase mirrors can produce vortex beams of extremely high intensity, so the interest in those beams in the ultrahigh frontier is growing quickly, and the present model can be of relevance for the understanding of some features of those experiments as well.
The pulse diagnostic is in itself a challenge: at extreme powers, atoms are ionized during the turn on of the pulse; therefore, the only remaining particles are electrons and positive ions (protons or other ions). Ions move much slower than electrons due to the mass/charge difference. Thus, the acceleration of the electrons is the dominant process for light scattering: laser-driven electrons move relativistically and radiate due to Thomson scattering (TS), and in turn, this scattering can be used to obtain information on the laser field itself.
Recently, the authors extended the computation of the electron dynamics and its TS radiation to the case of ultra-intense beams with net orbital angular momentum [16]. In that work, the study of TS in OAM beams (heavily) relied on numerical computations. Analytical tools that contribute to an understanding of these numerical results are worth investigating. Accordingly, we considered TS with incoming monochromatic plane waves in the most general state of polarization, i.e., elliptical polarization, and with OAM, both with a nontrivial initial phase dependence. In the particular case of initial phase dependence, our results extend previous ones (see for example [17] and the references therein) where it was clearly shown that in the ultra-intense domain, such an influence has to be taken into account to fully understand electron dynamics. As shown in this paper, the results for elliptically polarized plane-wave radiation, besides their own interest, provide an excellent approximation to the OAM dynamics for short time spans (typically covering a few tens of optical cycles). We also discuss the conditions for the approximations to hold.
Besides generalizing previous results, another reason to study such plane-wave solutions is to compute TS quantities (Doppler frequencies, spectral or spatial intensity distribution, etc.), with potential interest for fusion plasmas/laser beam diagnostics. In plasma physics, classical incoherent Thomson scattering (TS) is the standard diagnostic technique for electron temperature and density measurement; see for instance [2,[18][19][20][21][22][23][24].
Our interest here focused exclusively on incoherent TS by relativistic electrons (which behave independently of each other). For collective scattering of laser radiation by electrondensity fluctuations with relativistic phase velocities, see [25] for experimental observations and [26] for a fully relativistic theoretical analysis.
There is a difference in emphasis when TS is used for diagnostic purposes (for example in plasma physics) and TS as issued from the interaction with ultra-intense lasers; in plasma fusion diagnosis, the change of the scattering frequency at the detector is mainly (or nearly completely) due to the initial distribution of velocities in the electron population and the scattering angle via the Doppler effect (the laser being important in setting the expected position of scattered frequencies). In TS enacted by an ultra-intense laser, however, the properties of the scattered radiation should include its effect as well, because it substantially and nonlinearly perturbs the electron motion. Providing theoretical and/or numerical tools for a unified understanding of Thomson scattering, either as a diagnostic of the electron distribution or as a diagnostic of the laser state, could contribute to exploiting to the limits the application of a technique with a long and fruitful history.
The present paper is organized as follows. Section 2 is devoted to the dynamical equations of motion, the asymptotic Liénard-Wiechert retarded radiation fields, and the spectral representation of the latter, in general. In Section 3, incoming elliptically polarized monochromatic plane-wave radiation is considered: the full solution of the dynamical equations of motion and several physically interesting quantities. Section 4 deals with the TS spectral representation of the asymptotic Liénard-Wiechert retarded radiation fields. For the same incoming radiation as in Section 3, the TS spectral representation of the Liénard-Wiechert fields and Doppler frequencies are given in Section 4. Section 5 deals with TS for incoming radiation with OAM and detailed comparisons of it with TS for incoming plane-wave radiations. Section 6 contains the conclusions and discussions. Appendix A collects certain useful coefficients that are defined and employed in Section 3. Appendix B discusses in outline form a new generalization of the Bessel functions, appearing in the spectral representation of the asymptotic Liénard-Wiechert retarded radiation fields for elliptic polarization. Appendix C summarizes the essentials of the spectral representation of the asymptotic Liénard-Wiechert retarded radiation fields with an initial phase dependence, for linear and circular polarizations. Appendix D summarizes the Gauss-Laguerre polynomials and electric and magnetic fields for incoming radiation with OAM.

Classical Equations of Motion for the Electron
We used the MKS system of units (see, for instance, [34]) throughout the paper. We considered a classical relativistic electron in R 3 with position x(t) and momentum p(t) at time t. The electron interacts with a classical electromagnetic field in vacuum. The incoming electric and magnetic fields are E i and B i , which by assumption, correspond to a (in general, non-monochromatic) radiation. The total electric and magnetic fields are, in general, the sums of the incoming ones and the dynamical fields radiated by the electron itself (after it interacts with the input beam field).
In order to solve the dynamical problem, an approximation method was used throughout this work, based on the assumption that the total electric and magnetic fields be replaced, respectively, by E i and B i . Then, the equations of motion of the relativistic electron, subject to the Lorentz force of the incoming electric and magnetic fields E i = E i (y = x(t), t) and B i = B i (y = x(t), t), read [27] e is the charge of the electron. c is the velocity of light in vacuum. Notice that dx/dt = c 2 p/[m 2 c 4 + c 2 p 2 ] 1/2 . With the above understanding of E i (y = x(t), t) and B i (y = x(t), t), the dynamical problem boils down to solving the nonlinear Equation (1) for the electron position x(t) and momentum p(t) at time t. Six suitable variables related to the initial position and momentum of the electron give the initial conditions to solve Equation (1). Let β 0 = v 0 /c be the normalized initial velocity of the electron before being affected by the laser. The detailed solution (analytical or numerical) is used later in the radiated Liénard-Wiechert fields. In the above approximation scheme, the effect of the radiation reaction (becoming increasingly appreciable, for a near-infrared intense laser pulse, beyond 10 23 W/cm 2 ) and the controversial possibility of runaway solutions (a feature associated with the Abraham-Lorentz equation) [27,34,37] were disregarded from the outset.

Asymptotic Liénard-Wiechert Retarded Radiation Fields
Let E(y, t) rr,∞ and B(y, t) rr,∞ be the asymptotic Liénard-Wiechert retarded radiation fields (subscripts rr and ∞) at the position y at ("detection") time t [34,35] for any incoming (in general, non-monochromatic) radiation. Both E(y, t) rr,∞ and B(y, t) rr,∞ depend nonlinearly on x(t ), the solution of Equation (1) at the "radiation" time t , and hence, on the input beam fields. For the TS situations treated here, it is a good approximation to consider that the observation point (the point where the detector lies) is located far away from the radiating electron. In other words, the electron is not far from the origin of the coordinates initially and does not separate too much from it in the course of its interaction with the laser. The asymptotic fields read: with R = y , n = R −1 y and β ≡ c −1 (dx/dt ) = (β 1 , β 2 , β 3 ). 0 is the dielectric permittivity of the vacuum. n (the scattering unit vector) indicates the direction at which a detector is located at y. Both β and the acceleration of the relativistic electron d 2 x/dt 2 in Equations (2) and (3) are taken at time t . For given detector position y and "detection" time t, for the asymptotic fields, the time t is obtained from t − c −1 n · x(t ) = t − c −1 R and R x(t ) . One has: n = (n 1 , n 2 , n 3 ) with n 2 1 + n 2 2 + n 2 3 = 1. With standard unit vectors i, j, k, n is typically chosen to be orthogonal to j: n = n 0 = sin θ 0 i + cos θ 0 k (hence, the scattering vector lies in the (y 1 , y 3 ) plane). A general scattering vector n will also be considered later.
Once the solution of Equation (1) for an incoming (in general, non-monochromatic) radiation has been obtained (analytically or numerically), the computation of B(y, t) rr,∞ in Equation (3) x(t ) numerically, as was done in [35]. On the other hand, Equations (5)-(8) enable computing the spectral fieldB(y, ω) just by evaluating those integrals over ξ (numerically or analytically), without needing to solve for the implicit equation t − c −1 n · x(t ) = t − c −1 R. We carried out exact analytic computations with Equations (5)-(8) for incoming elliptically, linearly, and circularly polarized plane-wave radiation in Sections 3, 4 and Appendix C, respectively. Some numerical computations for them will also be carried out there, in order to check the feasibility of the latter. Numerical computations for Equations (5)-(8) with incoming radiation with OAM, for which no analytical solution of Equation (1) is available, will be considered in Section 5.

Representation of an Arbitrary Non-Monochromatic Plane Wave through the Vector Potential
The radiation (or Coulomb) gauge was employed. Let A i be the transverse vector potential of the incoming electric and magnetic fields E i and B i . E i , B i , and A i depend on a three-dimensional position y = (y 1 , y 2 , y 3 ) and on time (t). One has in the radiation gauge: Let the input beam field, corresponding to an arbitrary non-monochromatic plane wave and fulfilling the charge-free Maxwell equations in vacuum, propagate along the with A i ≡ dA i /dξ and k being the standard vector along the y 3 axis. Variables with overbars always denote two-dimensional vectors in the (y 1 , y 2 )-plane. A general case was considered in which β 0 is not necessarily orthogonal to A i , that is β 0 · A i = 0 was not imposed.
The incoming electromagnetic field is described by unique electric and magnetic fields (E i and B i ), and the dynamics and trajectory of the relativistic electron are also uniquely determined. On the other hand, the vector potential A i orĀ i is not unique: the gives rise to the same E i and B i (gauge invariance of the electromagnetic field). The electron trajectory is independent of δĀ i . See [36].

The Vector Potential and the Solution of the Electron Dynamics for Monochromatic Elliptically Polarized Radiation
The (real) electric field for a plane wave with elliptic polarization is given by: where Re denotes the real part; the amplitudes E 0 , E 1 , the unit vectors 1 and 2 , and the phase ϕ 0 are all real. 1 and 2 are: θ d is a real angle that gives the orientation of the principal axes of the polarization ellipse with respect to the (x, y) axes. ω 0 (real and > 0) is the frequency. The vector potential, with ζ = ω 0 ξ, reads: whereā co = E 0 1 ,ā si = E 1 2 . Let: As an initial condition for A i , we always imposed that A i = 0 at ζ = ζ 0 (typically and without loss of generality, ζ 0 = 0). Then: Notice that this is a particular choice of gauge; see [36]. The phase ϕ 0 has a physical meaning, and it is gauge-independent. Recall that the physical (gauge-independent) timedependent momentum isp(ξ) =f − eĀ i . Then, with the present gauge choice, the constant of motionf (gauge-dependent) coincides with the physicalp(ξ = 0).x(ξ), We solved the equations of motion for the vector potential given by A i (ξ) by following [27]. The initial conditions for the electron coordinates are x 10 , x 20 , x 30 . The initial conditions for the electron momentum p 10 , p 20 , p 30 follow from the three integration constantsf = ( f 1 , f 2 ) and γ 1 ; see [36]. Introducing (x 1 = k 0 x 1 , and so on) and (p 1 = p 1 /mc, and so on), we obtained: and:p The parameters appearing in the formulas above are given in Appendix A. The Lorentz factor γ, the normalized velocity β = (β 1 , β 2 , β 3 ) = v/c, and the acceleration a = (a 1 , a 2 , a 3 ) = dv/dt can also be computed in closed form in terms of the parameter ζ, namely: Coordinate time (t ) and also the time at the detector (t) admit closed (although involved) expressions in terms of ζ: from: and performing the integration of the left-hand side of the last equation (which involves only elementary functions), t is obtained as an explicit function of ζ, namely: The differential relationship that connects trajectory time t and the time at the detector t ("detection" time) is given by dt = (1 − n) · β)dt ; see [32,36]. Taking into account that in radiation problems, it is (in general) a very good approximation to put n = constant, n being a constant scattering vector, and making use of the differential relationship between dt and dζ, one obtains after elementary transformations that: In this case, we considered a unit scattering vector n(= (n 1 , n 2 , n 3 ) = (n, n 3 )) that is not lying necessarily in the (y 1 , y 3 )-plane, that being a generalization of the treatment presented in our previous paper [36]. The integration of the left-hand side of Equation (34) can also be carried out using elementary functions, and t is also an explicit function of ζ: The coefficients τ l , τ c , τ s , τ 2s , τ l , τ c , τ s , τ 2s are given in Appendix A. t 0 and t 0 are integration constants that satisfy t = t 0 and t = t 0 at ζ = ζ 0 .
All this algebraic machinery allows also obtaining explicit formulas for the Liénard-Wiechert radiated fields at the detector, using their definition from Equations (2) and (3), with all the variables as functions of the parameter ζ. The parametric formulas, together with the relationship between ζ and t already mentioned, allow for a complete description of Liénard-Wiechert fields in the time domain in a computationally very efficient way. This representation is the counterpart to the frequency domain representation that will be dealt with in the next section and contains, of course, the same information.

Magnitudes Derived from the Parametric Representation of Trajectories
The parametric representation of dynamics obtained in the previous subsection allows for an analytical computation of several quantities of interest associated with the trajectory in a rather straightforward manner, namely the fundamental period/frequency of the trajectory, the fundamental period/frequency at detector position, the drift velocity (defined as the velocity at which the particle is at rest on average), the extremal values for the momentum components and the γ factor. From (33) and (35) and taking two ζ values 2π apart from each other, one obtains for the fundamental period of the trajectory, the fundamental period, and the fundamental frequency (ω f = 2π/T f ) at the detector the expressions: Due to the relevance of Ω f in nonlinear TS problems, let us write it down explicitly as: In this form, the dependencies of Ω f on the initial momentum (throughf , s/c, γ 1 ), the scattering vector (through n), and the parameters characterizing the laser (c 1 ,Ā 0 ) are apparent. SinceĀ 0 depends on the initial phase of the laser field ϕ 0 , the fundamental Thomson scattering frequency also depends on it in the high intensity regime. It is clear that this dependency on the phase is absent ifĀ 0 vanishes; moreover, if c 1 andĀ 0 vanish or are very small (the case of a low-intensity laser), one recovers the standard Doppler frequency.
The dependency of Ω f on several kinematical (i.e., associated with the initial position and the initial momentum of electron) and laser parameters (intensity, polarization, helicity) makes it possible (to some extent) to select the frequency at an specific detection point, that is an engineering of the radiated spectrum seems feasible by a judicious choice of the control parameters. Let us illustrate this in Figure 1 with a simple example with the following setup. In the upper panel, the electron is initially at the origin of the coordinates with E kin0 = 0. The reference electric field of the incoming laser is chosen E r = 5.0 × 10 12 V/m. E 0 , E 1 are defined from E r and the angle θ s as E 0 = E r cos θ s , E 1 = E r sin θ s . In the lower panel, E kin0 = 200.0 keV, and the angles giving the orientation of the initial momentum are θ = π/2, ϕ = π/4; the laser parameters are the same. The point at which the scattered radiation is measured was assumed to lie 0.5 m along the positive direction of the y axis. The opportunities for scattered frequency engineering are clearly displayed in the two panels, by playing with either kinematical or laser parameters, including, in a non-trivial way, the dependency on the initial phase (ϕ 0 ) and the polarization/helicity of the laser beam. Please, note that in the figures of the paper, "adim." stands for "adimensional", meaning a magnitude with no dimensions, i.e., a pure number.   The fundamental period (T f und ) of the trajectory (and hence, its fundamental frequency) also has a non-trivial dependency on ϕ 0 and the laser polarization state, inherited from the parameters τ l , τ c , τ s and τ 2s . This is shown (Figure 2) in a case where the laser has a reference peak electric field of E r = 4.257 × 10 12 V/m, E 0 = E r cos θ s , E 1 = E r sin θ s , and the electron has an initial kinetic energy of E kin0 = 10.0 keV, its initial velocity being oriented with respect to the laser propagation direction (z− axis) at angles θ = 0.807433, ϕ = 4.32799. Three different polarizations corresponding to θ s = 0, θ s = ±π/4, θ s = ±π/8 are displayed. The vertical axis of this figure, labeled as "Fundamental Period (adim.)", is defined as ω 0 × T f und .  The important question of the field amplitudes (and hence, the power) associated with the different frequency components at the detector, i.e., the Thomson scattering spectrum, is perhaps more conveniently addressed in the frequency domain and shall be dealt with in the next section.
The equations for p 1 , p 2 extremal values are easily solved; the one for the extremal values of p 3 or γ (which are the same due to the relationship mcγ = γ 1 + p 3 , γ 1 constant) is slightly more involved, but it is nevertheless amenable to efficient numerical computation. As an example, the non-trivial dependency of the maximum of γ factor on ϕ 0 is shown in Figure 4, where, for the same conditions as mentioned for Figure 2, the initial phase of the laser beam is systematically varied in the interval [−π, π].
We omitted a lengthy analytical computation. We found for Λ(ξ) E : All g n,E 's are gauge-independent. The parameter s was defined in [36]. With the definitions given above, the coefficients g n,E , n = 0, 1, 2, 3 are dimensionless. The treatment given here includes the cases in which β 0 · A i may be = 0 or = 0.

Incoming Laser Beam with OAM: Vector Potential and Electromagnetic Fields
The input beam field was chosen to be an electromagnetic wave in vacuum that propagates along the z axis, from −∞ towards +∞, and corresponds to a Gauss-Laguerre mode in the (x, y)-plane and, so, having some given orbital angular momentum (OAM). The fields of the OAM mode are not strictly transverse: they have small components in the direction of propagation (namely, the z axis). See [6,16,40]. The radiation (or Coulomb) gauge will be employed. Let A i be the vector potential of the incoming electric (E i ) and magnetic (B i ) fields. E i , B i , and A i depend on a three-dimensional position y = (x, y, z) and on time (t). Equation (9) continues to hold in the radiation gauge. By assumption, the (real) vector potential A i lies in the (x, y)-plane (namely, the (y 1 , y 2 )-plane): A i (y; t) = A x i + A y j, with: (62) As usual, c.c. denotes the complex conjugate of the preceding term. A 0 is a real amplitude, ρ = (x 2 + y 2 ) 1/2 . η 0 and η 1 are complex numbers, with | η 0 | 2 + | η 1 | 2 = 1. We write: π 0 , π 1 , σ 0 , and σ 1 are real constants fulfilling the condition π 2 0 + π 2 1 + σ 2 0 + σ 2 1 = 1 and describing the polarization state of the field. Although we concentrated on the case of an incoming laser beam linearly polarized along the x axis, other polarization states can be treated as well by a suitable selection of π 0 , π 1 , σ 0 , and σ 1 . The Laguerre-Gauss (LG) functions LG l,p (ρ, z) are given in Appendix D. On the other hand: with φ = arctan(y/x). χ 0 is a constant phase. w 0 is the beam-waist parameter. z 0 is the Rayleigh range, defined by z 0 = k 0 w 2 0 /2 with k 0 = 2π/λ 0 . ω 0 (real and > 0) is the frequency, and ω 0 = k 0 c.
For an incoming laser beam linearly polarized along the x axis, one has π 1 = σ 0 = σ 1 = 0 and E y = 0 for any t, (x, y, z). Then, A y = 0, B x = 0, while E x , B y , and B z are non-vanishing.

Approximating the Electron Dynamics in OAM Beams by the Plane-Wave Solution: Numerical Computations
One of the reasons to study in depth the plane-wave solution, as performed in previous sections, is of course to generalize previous results, including the most general state of polarization for this type of wave (i.e., elliptical polarization) and the non-trivial initial phase dependency of dynamical variables. Another important reason is to provide a tool for the approximate analytical (or semi-analytical) treatment of more complex laser beams, like the ones with OAM, where one has to (heavily) rely on numerical computations to obtain an understanding of their dynamics. See [16] for a description and implementation of the numerical methods used to characterize the electron dynamics and the ensuing Thomson scattering emission in laser beams with OAM.
The rule that has been found to apply to obtain a good fit between OAM and planewave solution is as follows. From the characteristics of the OAM beam and the initial position of the electron, compute the local values of the laser parameter and phase, namely α local = 2A 0 LG l,p (x 10 , x 20 , x 30 )/(mc) and ψ local = ψ(x 10 , x 20 , x 30 , t 0 ). Introduce these values in the plane-wave solution with the same polarization state as the OAM wave and use the same initial electron kinetic energy and the same values for the angles θ, ϕ.
By means of the examples to follow, we claimed here that the plane-wave solutions for the dynamics of the electron can be, in some instances to be described, a very good fit to the actual dynamics in beams with OAM. The conditions to be met for the fitting are basically (i) that local values for the adimensional laser parameter and phase are chosen in OAM, and put into the plane-wave solution, (ii) that the transverse excursion of the electron is not too large, and (iii) that the dynamics is studied for a limited number of optical cycles. If these conditions are met, the plane-wave solution can fit OAM solutions to a high degree of precision for a wide range of parameters in electron energy, laser intensity, or polarization.
The electron dynamics was computed in an OAM laser beam having the following parameters: λ 0 = 800 nm, w 0 = 100 µm, χ 0 = 0, (l, p) = (1, 0), A 0 ranging from 2.5 × 10 −3 V·s/m to 7.5 × 10 −3 V·s/m. Polarization, if not otherwise stated, was linear along the x axis. Electrons with initial kinetic energy ranging from zero to 100 keV were studied; their initial positions were randomly chosen in the intervals . Polar (θ) and azimuthal (ϕ) angles giving the initial momentum direction (with respect to the z axis or laser propagation direction) were randomly chosen in the intervals [0, π], [0, 2π]. In Figures 5-9, the color code is as follows: red, green, and blue colors are assigned to the x-, y-, or z-component of the corresponding magnitude (in the case of vector quantities). For example, the x-component of velocity will be assigned red color in the figures or the z-component of the radiated electric field will be assigned blue color, etc. Figure 5 shows the comparison of normalized velocity components, the γ factor, and the radiated E field at the detector (two orthogonal quadratures along the x and z axes). The parameters were A 0 = 2.5 × 10 −3 V·s/m (OAM beam), E kin0 = 10.0 keV, x 10 = 98.4242 µm, x 20 = −52.9688 µm, x 30 = −3.57517 µm; θ = 0.807433, ϕ = 4.32799. It can be seen that the fitting was quite good in this example, and a simple measure of it can be obtained if maximum values of the γ factor are compared: the maximum of γ OAM = 3.213 vs. the maximum of γ plane = 3.216, i.e., a relative difference of less than 10 −3 . It is clear that radiated spectra will be very similar due to the similarity of radiated fields at the detector in the time domain, although no systematic comparison is this connection was attempted in this paper.
A second example is shown ( Figure 6) in which A 0 = 2.5 × 10 −3 V·s/m (OAM beam), E kin0 = 0.0 keV, x 10 = 28.0601 µm, x 20 = −1.38741 µm, x 30 = 0.108658 µm. In this case, some drift was apparent in the velocity of the OAM solution; hence, the fitting was not perfect, since the plane-wave solution simply cannot capture the changes induced by large variations in the transverse profile of the OAM mode. It was nevertheless an overall good fit for a limited section of the trajectory, some of its main characteristics being captured by the more simple plane-wave solution. Figure 7 shows a case in which A 0 = 5.0 × 10 −3 V·s/m (OAM beam), E kin0 = 10.0 keV, x 10 = 71.0055 µm, x 20 = 30.9694 µm, x 30 = −0.262421 µm, θ = 0.0433, ϕ = 0.6058. The plane-wave solution closely fits the OAM solution in this case, which showed a large excursion of the relativistic γ factor, up to nearly seven. The maximum value of γ OAM = 6.70 should be compared with the maximum of γ plane = 6.66, i.e., a relative difference of less than 6.0 × 10 −3 .
To end this relation of comparative studies between the dynamics and radiated fields in laser beams with OAM and the corresponding analytical solution, let us present a case in which the OAM beam has a genuine elliptical polarization state, characterized by (π 0 , π 1 , σ 0 , σ 1 ) = (cos π/8, 0, 0, sin π/8), and compare it with the corresponding elliptically polarized analytical solution (Figure 8). The rest of the parameters in the comparison were the same as in Figure 7. The match of OAM vs. the analytical solution was also very good here, showing the ability of these extended solutions to also cover the case of complex OAM beams with arbitrary (in particular, genuine elliptical) polarization.    A simple quantitative estimation of the conditions in which a good match between OAM and plane-wave solutions is expected is as follows: if the transverse excursion of the electron during the integration time for a particular trajectory is much less than the size of the OAM laser beam w 0 , a good match is expected. The reason would be that in that case, the transverse variation of the laser intensity is less important than the fact that the local OAM laser field can always be faithfully fitted by a plane wave of adequate intensity, phase, and polarization. In formulas: c β 2 1,av + β 2 2,av 2πN cycles /ω 0 << w 0 , where T int = 2πN cycles /ω 0 means the integration time of the trajectory, β 1,av and β 2,av are the transverse components of the drift velocity (as obtained from the plane-wave solution), w 0 sets the transverse spatial scale of the OAM beam (w 0 = 100.0 µm is used throughout this section), and λ 0 is the laser wavelength, equal to 800.0 nm in our case.
In the course of the investigations reported here, a new feature was noticed in very long integrations (involving several hundreds of optical cycles). In such a case and for OAM beams with (l, p) = (1, 0), we observed trajectories that perform a finite, i.e., not drifting, movement in the transverse X-Y plane, while they acquired some net drift velocity in the Z-direction. We provisionally named them as trapped trajectories or solutions. They typically appeared if the initial conditions in the transverse direction were close to the beam center and were found not only for E kin0 = 0, but also for a substantial range of initial kinetic energy. The character of the dynamics was not easy to ascertain, but evidence from Fourier analysis was compatible with a quasi-periodic motion involving combinations of close, but incommensurate frequencies. A specific challenge to study these solutions is that the integration must cope with both fast and slow time scales, the finite transverse motion being typically 100 times slower than the optical frequency. Figure 9 shows an example in which 512 optical cycles were integrated with the following initial conditions for the electron and the OAM beam: x 10 = 1.03196 µm, x 20 = −4.82347 µm, x 30 = −0.107291 µm, E kin0 = 2.5 keV, initial velocity along the positive x axis direction; (l, p) = (1, 0), λ 0 = 800.0 nm, w 0 = 100.0 µm, χ 0 = 0, A 0 = 5.0 × 10 −3 V·s/m, the laser beam linearly polarized along the x axis. Please note that Figure 9 represents just a particular case and that the relative amplitude and phase of the X and Y components and the z axis drift shown in it can substantially change with the initial condition, the initial velocity, and the laser parameters. There is ongoing research to more fully characterize this kind of solution, which seems to go beyond the theoretical framework developed in this paper and which poses some numerical and theoretical challenges to their understanding.

Conclusions and Final Comments
Regarding the interrelated research Items (a) and (b) stated in the Abstract, the following results were obtained.
(a) The classical nonlinear incoherent TS from a single relativistic electron with an incoming plane wave monochromatic laser radiation with frequency ω 0 , wave vector k 0 , initial phase ϕ 0 , elliptical polarization, and any intensity was investigated analytically and numerically. A complete explicit solution in terms of the parameter ζ = ω 0 t − k 0 z and with full account of the non-trivial influence of ϕ 0 = 0 in dynamical variables, and hence in the scattered spectrum, was given. In particular, the cases of TS with linear and circular polarizations and ϕ 0 = 0 were treated, thereby extending previous studies [36] for ϕ 0 = 0.
In so doing, we also carried out a detailed comparison in the TS framework between the analytical solution and ab initio numerical solutions of the Lorentz relativistic equation. Specifically, in Equations (5)-(7), we used successively the analytical solution and the numerical solution: an excellent match was found between both approaches. This contributes to confirming the reliability of Equations (5)-(7) when ab initio numerical solutions are fed in.
Detailed formulas were provided for the fundamental TS (Doppler) frequency for a general scattering vector, with full account of the influence of the laser intensity, polarization, and initial phase ϕ 0 . The dependences of TS frequency on the kinematical and laser parameters (in particular, on ϕ 0 ) were displayed as a means to engineer the scattered spectrum. The explicit solutions for general elliptic polarization gave rise to very efficient computations of the electron TS Doppler frequencies, the period of trajectories, the drift velocities and relativistic factors for electrons, etc., with potential interest for TS in fusion plasmas. The dependences on ϕ 0 are displayed in Figure 1 (electron TS Doppler frequencies for linear, circular, and elliptical polarizations), Figure 2 (periods of electron trajectories for linear, circular, and elliptical polarizations), Figure 3 (electron drift velocities), and Figure 4 (electron relativistic factor).
(b) The analytical results for incoming plane-wave radiations obtained in this paper were successfully employed to approximate (for the sake of a tentative understanding) TS for the more complex laser beams with intrinsic orbital angular momentum (OAM). Previous numerical TS results for such beams [16] and new and extended numerical computations for them were used. The results were compared with the ones for the analytical solutions with plane waves presented here. They were found to be in very good agreement with each other for time spans of several optical cycles, provided that the local laser parameter and the local phase read from the OAM beam were used in the analytical solution. This agreement was shown to hold for a wide range of laser intensities, electron energies, and for several polarization states (Figures 5-7), in particular and specifically for genuine elliptical polarization (Figure 8). We also analyzed and stressed the non-trivial role played by the initial phase ϕ 0 on several dynamical variables and in the scattered spectrum. According to the results presented in the paper, the local phase of the OAM laser beams and the local (transverse) variation of the laser intensity explained to a large extent the scattering spectra computed and discussed in [16].
The two cases considered in the paper are of growing current interest: amplifying mechanisms leading to petawatt peak powers generate linearly polarized laser fields in most cases. That explains why, so far, most of the research has been on such a polarization state or in the case of circular polarization, although at ultrahigh intensities, circular polarization is difficult to achieve. However, in many cases, the mirrors of the beam transport from compressor to target can induce some degree of elliptical polarization (intentional or not). Therefore, the present model describing an arbitrary degree of polarization can be of broad interest to model future experiments. On the other hand, there is a growing interest in OAM beams because of the new possibilities they bring. Until recently, ultra-intense OAM beams were considered very difficult to achieve; however, it has been shown that off-axis spiral phase mirrors can produce vortex beams of extremely high intensity, so the interest in those beams in the ultrahigh frontier is growing quickly, and the present model can be of relevance for the understanding of some features of those experiments as well.
A new feature was observed, for OAM beams with (l, p) = (1, 0), in very long integrations (involving several hundreds of optical cycles). Specifically, we found trajectories that perform a finite, i.e., not drifting, movement in the transverse X-Y plane, while they acquire some net drift velocity in the Z-direction. See Figure 9. Such trajectories typically appeared if the initial conditions in the transverse direction were close to the beam center and were found not only for E kin0 = 0, but for a substantial range of initial kinetic energy.
For TS with an incoming plane-wave radiation with elliptical polarization, the analytical expressions obtained for the scattered spectrum were given in terms of infinite series of new generalizations of the standard Bessel function. Such new generalizations (so far unknown and, hence, not previously reported, to the best of our knowledge) constitute an extension of the already known generalized Bessel functions [39].
Finally, it is perhaps worth saying that petawatt lasers are a reality (both for shot pulses, femtoseconds, and for long pulses, picoseconds) and that multi-petawatt lasers are expected to be fully operative in a few years. Therefore, the results of the paper dealt with a realistic experimental situation. OAM beams are easily obtained from those lasers, using for example off-axis spiral phase mirrors; see for example [15]. Off-axis spiral mirrors of more than 10 cm in diameter and with high precision in the phase exist now. The technique is quite promising, and larger spiral phase mirrors are in construction; therefore, the study of the dynamics of electrons driven by such intense OAM beams can be of high interest. J G,l (x, y; ϕ 2s,E − ϕ 1s,E ) is a natural generalization, arising in our present study of TS for an incoming elliptically polarized plane wave, of the generalized Bessel function J l (x, y) [39], which, in turn, has appeared upon treating TS for an incoming linearly polarized plane wave. The above properties of J G,l (x, y; ϕ 2s,E − ϕ 1s,E ) generalize those of the generalized Bessel function J l (x, y) [39].
The generalized Doppler-like formula is: ω L,l = l ω 0 g 1,L (A56) for the given g 1,L (see (A51)) and for all the integers l (= 0, ±1, ±2, ±3, ...). If ω varies in a finite interval, then only a discrete set of values for l contributes. The inclusion of the initial phase ϕ 0 makes the coefficients g 1,L and g 2,L depend on it, and since they account, respectively, for the scattered frequency and the amplitude of the scattered spectrum, the latter two quantities acquire a non-trivial dependence on ϕ 0 as well. Such an effect was not accounted for in our previous work [36].
describing some electromagnetic wave with OAM associated with the prescribed l, p.