Direct Laser-Driven Electron Acceleration and Energy Gain in Helical Beams

Etele Molnár; Dan Stutman

doi:10.1155/2021/6645668

Direct Laser-Driven Electron Acceleration and Energy Gain in Helical Beams

Published online by Cambridge University Press: 01 January 2024

Etele Molnár

and

Dan Stutman

Show author details

Etele Molnár*: Affiliation:
Extreme Light Infrastructure–Nuclear Physics ELI-NP, Horia Hulubei National Institute for Physics and Nuclear Engineering, 30 Reactorului Street, RO-077125, Bucharest-Magurele, Romania
Dan Stutman: Affiliation:
Extreme Light Infrastructure–Nuclear Physics ELI-NP, Horia Hulubei National Institute for Physics and Nuclear Engineering, 30 Reactorului Street, RO-077125, Bucharest-Magurele, Romania
*: Correspondence should be addressed to Etele Molnár; etele.molnar@eli-np.ro

Article contents

Abstract
Introduction
Direct Laser-Driven Electron Acceleration in Helical Beams
Results
Conclusions
Data Availability
Conflicts of Interest
References

Abstract

A detailed study of direct laser-driven electron acceleration in paraxial Laguerre–Gaussian modes corresponding to helical beams LG0m with azimuthal modes m = {1,2,3,4,5} is presented. Due to the difference between the ponderomotive force of the fundamental Gaussian beam LG00 and helical beams LG0m, we found that the optimal beam waist leading to the most energetic electrons at full width at half maximum is more than twice smaller for the latter and corresponds to a few wavelengths Δw0 = {6,11,19}λ0 for laser powers of P0 = {0.1, 1,10} PW. We also found that, for azimuthal modes m ≥ 3, the optimal waist should be smaller than Δw0 < 19λ0. Using these optimal values, we have observed that the average kinetic energy gain of electrons is about an order of magnitude larger in helical beams compared to the fundamental Gaussian beam. This average energy gain increases with the azimuthal index m leading to collimated electrons of a few 100 MeV energy in the direction of the laser propagation.

Information

Type: Research Article
Information: Laser and Particle Beams , Volume 2021 , 2021 , e11

DOI: https://doi.org/10.1155/2021/6645668 [Opens in a new window]
Creative Commons: This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright: Copyright © 2021 Etele Molnár and Dan Stutman.

1. Introduction

The well-known Laguerre–Gaussian (LG_pm) modes of paraxial light with integer radial and azimuthal, p and m, indices have a phase distribution of $e^{- i m θ}$ , where θ is the azimuthal angle. Here, m corresponds to the azimuthal or rotational order of the mode with a well-defined orbital angular momentum (OAM) of mℏ per photon [Reference Allen, Beijersbergen, Spreeuw and Woerdman1]. Such higher-order modes lead to twisted light by forming an optical vortex in the direction of propagation associated with the optical phase singularity. The spiral shape of the wavefront is formed by |m| interconnected helices; hence, ${LG}_{p, m \neq 0}$ beams are called helical beams.

These distinctive properties of higher-order Laguerre–Gaussian modes set forth to an ever-increasing amount of applications in different areas of optical communications, imaging techniques, quantum information technologies, and other topics; see [Reference Grier2–Reference Krenn, Handsteiner and Fink6] and references therein.

In recent years with the development of ultraintense multi-PW laser technology, Laguerre–Gaussian laser modes have other potential applications in the fields of plasma accelerators and inertial confinement fusion and in the generation of X-rays and γ-rays with OAM [Reference Zhang, Shen and Yin7–Reference Zhang, Chen and Luo14]. In particular, it has been shown that, in laser wakefields driven by LG₀₁ helical pulses, the wakefield shows a donut-like structure with a ring-shaped hollow electron beam [Reference Zhang, Chen and Luo14]. Furthermore, for a lower-density plasma or a smaller laser spot size, besides the donut-like wakefield, a central bell-like wakefield forms in the center of the donut-like wake. On the contrary, further reducing the plasma density or laser spot size leads to an on-axis electron beam acceleration only.

It was also demonstrated that Laguerre–Gaussian beams transfer a part of their OAM to electrons through the dephasing process similar to the direct electron acceleration (DLA) in Gaussian beams [Reference Nuter, Korneev, Dmitriev, Thiele and Tikhonchuk15]. Furthermore, the propagation of optical beams with OAM leads to plasma waves that may also carry OAM which couple to the plasma electrons and involve Landau damping and particle acceleration accompanied with the generation of quasi-static axial and azimuthal magnetic fields [Reference Blackman, Nuter, Korneev and Tikhonchuk16]. When Laguerre–Gaussian plasma waves are subjected to Landau damping, a higher azimuthal mode number leads to a larger OAM transfer to particles traveling close to the phase velocity of the plasma wave [Reference Blackman, Nuter, Korneev and Tikhonchuk17]. Vacuum-based charge acceleration with Laguerre–Gaussian beams has also been studied very recently [Reference Vaziri, Golshani, Sohaily and Bahrampour18, Reference Akou and Firouzjaei19] showing that it is possible to generate GeV high-quality electron bunch with low spread in energy and radial deflection.

Motivated by these interesting results, in this paper, we will study direct electron acceleration in vacuum in various helical Laguerre–Gaussian LG_0m laser pulses corresponding to helical modes m = {0,1,2,3,4,5}. Laser beams with p = 0 and m = 0 define the fundamental Gaussian mode, while beams with p = 0 and m ≥ 1 have a hollow ring-like transverse intensity profile with zero intensity at the center. The central hollowness corresponds to a potential well which confines/accelerates the electric charges through the transverse/longitudinal ponderomotive forces.

This type of ponderomotive trap is also realized using higher-order transverse electromagnetic modes (TEMs), i.e., the Hermite–Gaussian modes, such as TEM_1,0 or the combination of TEM_1,0 with TEM_0,1. This ponderomotive potential is similarly axisymmetric and has a minimum on the axis [Reference Chaloupka, Fisher, Kessler and Meyerhofer20–Reference Li, Gu and Zhu24]. This is certainly expected since orthogonal Hermite–Gaussian modes may be decomposed into Hermite–Gaussian modes with a phase difference, and vice versa; see [Reference Allen, Beijersbergen, Spreeuw and Woerdman1, Reference Abramochkin and Volostnikov25] for the general formulas.

Here, we specifically focus our study on the energy gain in femtosecond lasers of low to very high intensity, thereby gaining valuable estimates about relevant parameters for lasers operating at ELI-NP [Reference Tanaka, Spohr and Balabanski26]. Our results are based on the 3-dimensional numerical solution of the relativistic equations of motion for free electrons in paraxial laser fields. Henceforth, similarly as in very low-density plasmas, in DLA, we also observe collimated and intertwining electron beams in the direction of the laser propagation, while in addition, we also show that the average energy gain increases with the azimuthal mode index m.

Furthermore, we found that the optimal beam waist leading to the most energetic electrons is more than twice smaller in case of higher-order Laguerre–Gaussian beams than in case of the fundamental Gaussian beam that was presented in [Reference Molnar, Stutman and Ticos27]. For lasers of P ₀ = {0.1, 1,10} PW power, the optimal beam waists are Δw ₀ = {6,11,19}λ₀ and correspond to I ₀ = I _0m × {0.76, 2.2, 7.6}10²¹ W/cm² peak intensities, while taking into account the intensity profiles of the LG_0m beams, these values are reduced by $I^{0 m} \approx (0.37, 0.27, 0.23, 0.19, 0.17)$ for m = {1,2,3,4,5}. Using these optimal values, we have observed that, for the same laser power, the average kinetic energy gain of electrons is about an order of magnitude larger in helical beams compared to the fundamental Gaussian beam.

This paper is organized as follows. In Section 2, we present some of the characteristic properties of helical beams, the equations of motion for electrons, and the initial conditions corresponding to our study. In Section 3.1, we have estimated the optimal values of beam waist leading to the most energetic electrons for given laser power. Applying these optimal values, we present and discuss the electron dynamics and energy gains in linearly polarized (LP) and circularly polarized (CP) helical beams in Section 3.2. The conclusions are summarized in Section 4.

2. Direct Laser-Driven Electron Acceleration in Helical Beams

2.1. Laguerre–Gaussian Beams

A well-known solution [Reference Siegman28, Reference Goldsmith29] to the paraxial wave equation is obtained in cylindrical coordinates (r, θ, z), with the help of generalized Laguerre polynomials. This solution is cylindrically symmetric around the axis of propagation z, with radius $r = \sqrt{x^{2} + y^{2}}$ and azimuth $θ = arctan (x / y)$ , expressed in Cartesian coordinates. These are the Laguerre–Gaussian (LG_pm) beams [Reference Allen, Beijersbergen, Spreeuw and Woerdman1, Reference Baumann and Pukhov10, Reference Chen, Hatsagortsyan and Keitel12], with radial index p and azimuthal index m. The general expression for the electric field distribution of a monochromatic LG_pm pulse is

(1)

\begin{matrix} E_{T, p m} (r, θ, z) = C_{p m} E_{0} \exp (- i k_{0} z + i ϕ_{0}) \\ \times \frac{w_{0}}{w (z)} \exp (- \frac{r^{2}}{w^{2} (z)} - i \frac{z r^{2}}{Z_{R} w^{2} (z)} + i \arctan (\frac{z}{Z_{R}})) \\ \times {(\frac{\sqrt{2} r}{w (z)})}^{(m)} \exp (i (2 p + (m)) \arctan (\frac{z}{Z_{R}})) \\ \times L_{p}^{(m)} (\frac{2 r^{2}}{w^{2} (z)}) \exp (- i m θ), \end{matrix}

where k ₀ = ω ₀/c is the wavenumber, $c = 1 / \sqrt{ɛ_{0} μ_{0}}$ is the speed of light in vacuum, ω ₀ is the angular frequency, E ₀ is the amplitude of the electric field, and ϕ ₀ is the initial phase. Furthermore, $Z_{R} = w_{0}^{2} k_{0} / 2$ is the Rayleigh range, $R_{c} (z) = (z^{2} + Z_{R}^{2}) / z$ is the radius of curvature, and $ϕ_{G} (z) = arctan (z / Z_{R})$ is the Gouy phase. The beam waist is defined as

(2)

\begin{matrix} w (z) = w_{0} \sqrt{1 + {(\frac{z}{Z_{R}})}^{2}}, \end{matrix}

where the beam waist radius at focus is w ₀ ≡ w(z = 0.

The Laguerre polynomials are denoted by $L_{p}^{m}$ , and the normalization constant, $C_{p m} = \sqrt{p! / (p + (m))!}$ , follows from the orthonormality of the Laguerre polynomials [Reference González de Alaiza Martínez, Duchateau and Chimier30]. The fundamental Gaussian beam LG₀₀ is obtained for p = m ≡ 0, where $C_{00} = L_{0}^{0} \equiv 1$ .

The components of the electric and magnetic fields are

(3)

\begin{matrix} E_{x, p m} = α_{x} E_{T, p m}, \\ E_{y, p m} = - i α_{y} E_{T, p m}, \\ B_{x, p m} = - \frac{1}{c} E_{y, p m}, \\ B_{y, p m} = + \frac{1}{c} E_{x, p m}, \end{matrix}

where $α_{x} = \sqrt{(1 + α_{P}) / 2}$ and $α_{y} = \sqrt{(1 - α_{P}) / 2}$ such that α_P = 1 or −1 in case of linear polarization along the x-axis or y-axis, respectively, and α_P = 0 for circular polarization, while elliptic polarization, otherwise.

The longitudinal components of the electric and magnetic fields are calculated from Maxwell’s equations, $\nabla \cdot \overset{⟶}{E} = \nabla \cdot \overset{⟶}{B} \equiv 0$ ; thus, in the paraxial approximation,

(4)

\begin{matrix} E_{z, p m} = - \frac{i}{k_{0}} (\frac{\partial E_{x, p m}}{\partial x} + \frac{\partial E_{y, p m}}{\partial y}) \\ = \frac{i}{k_{0} w^{2} (z)} (x E_{x, p m} + y E_{y, p m}) \\ \times (2 (1 + i \frac{z}{Z_{R}}) - (m) \frac{w^{2} (z)}{r^{2}} + \frac{4 L_{p - 1}^{(m) + 1} (2 r^{2} / w^{2} (z))}{L_{p}^{(m)} (2 r^{2} / w^{2} (z))}) \\ - \frac{m}{k_{0} r^{2}} (y E_{x, p m} - x E_{y, p m}), \\ B_{z, p m} = - \frac{i}{c k_{0}} (\frac{\partial E_{x, p m}}{\partial y} - \frac{\partial E_{y, p m}}{\partial x}) \\ = \frac{i}{c k_{0} w^{2} (z)} (y E_{x, p m} - x E_{y, p m}) \\ \times (2 (1 + i \frac{z}{Z_{R}}) - (m) \frac{w^{2} (z)}{r^{2}} + \frac{4 L_{p - 1}^{(m) + 1} (2 r^{2} / w^{2} (z))}{L_{p}^{(m)} (2 r^{2} / w^{2} (z))}) \\ + \frac{m}{c k_{0} r^{2}} (x E_{x, p m} + y E_{y, p m}) . \end{matrix}

Therefore, the electromagnetic field of Laguerre–Gaussian pulses is given by

(5)

\begin{matrix} {\overset{⟶}{E}}_{p m} (t, r, θ, z) = Re ({\overset{⟶}{E}}_{p m} (r, θ, z) g (t, z)), \end{matrix}

(6)

\begin{matrix} {\overset{⟶}{B}}_{p m} (t, r, θ, z) = Re ({\overset{⟶}{B}}_{p m} (r, θ, z) g (t, z)), \end{matrix}

where the Gaussian temporal envelope with τ₀ duration and peak intensity position at z_F reads [Reference Robinson, Arefiev and Neely31]

(7)

\begin{matrix} g (t, z) = \exp (i ω_{0} t - {(\frac{t - (z - z_{F}) / c}{τ_{0}})}^{2}) . \end{matrix}

Note that there are more appropriate choices for the temporal profile such as the hyperbolic secant $g (t, z) \sim 1 / cosh ((t - (z - z_{F}) / c) / τ_{0})$ ; see [Reference Kirk32–Reference Ong, Moritaka and Takabe34] for more details. Using this profile, we found that the energy gain may be reduced by as much as 30% compared to Gaussian.

2.2. Helical Beams

The Laguerre–Gaussian beams, LG_0m, with nonzero azimuthal modes |m| ≠ 0 contain a phase change given by $e^{- i m θ}$ . Note, however, that, for all modes, the generalized Laguerre polynomials have a contribution equal to one, i.e., $L_{0}^{m} (x) = L_{0} (x) = 1$ . These type of Laguerre–Gaussian beams, LG_0m, are the helical beams associated with the nonzero OAM of light [Reference Allen, Beijersbergen, Spreeuw and Woerdman1, Reference Longman and Fedosejevs35].

In Figure 1, the normalized intensity profiles of linearly polarized fundamental Gaussian and different helical beams are shown as a function of radial distance, in units of beam waist radius w ₀. These plots represent $I^{0 m} = {(E_{x}^{0 m} / E_{0})}^{2}$ using equation (1), for t = z ≡ 0 and k ₀ = ω ₀ = w ₀ ≡ 1, while the intensity is $I_{0}^{0 m} \sim E_{0}^{2} I^{0 m}$ .

Figure 1:

The relative intensity profiles at constant power and beam waist as a function of the radius of an LP laser: the fundamental Gaussian beam LG₀₀ with dotted black line and the helical beams LG₀₁, LG₀₂, LG₀₃, LG₀₄, and LG₀₅ with green, magenta, red, blue, and black lines correspondingly.

The intensity of LG_0m-modes is largest for m = 0 corresponding to the fundamental Gaussian beam, and it is decreasing with an increasing number of azimuthal modes. Due to the ${(\sqrt{2} r / w (z))}^{(m)}$ factor, all modes with $(m) \geq 1$ have zero intensity at the center, r = 0, that is known as an optical vortex or phase singularity on the axis.

The intensity profiles of Gaussian beams are concave functions on the whole interval. The intensity of helical beams is independent of θ and has a convex part with maxima at $r_{m} = w_{0} \sqrt{(m) / 2}$ after which the functions change from convex to concave. The width of the convex part widens with increasing azimuthal index m, while using the positions of maxima, the intensity peaks of LG_0m-modes relative to fundamental Gaussian lead, $I^{0 m} (r_{m}) = C_{0 m} {(m)}^{(m)} e^{- (m)}$ , as shown in Figure 1.

2.3. Electron Acceleration in the Electromagnetic Field

The motion of electrically charged particles in an external electromagnetic field is governed by the Lorentz force, ${\overset{⟶}{F}}_{L} \equiv d \overset{⟶}{p} / d t = q (\overset{⟶}{E} + \overset{⟶}{v} \times \overset{⟶}{B})$ , and leads to the following set of nonlinear differential equations [Reference Jackson36]:

(8)

\begin{matrix} \frac{d \overset{⟶}{x}}{d t} = c \overset{⟶}{β}, \end{matrix}

(9)

\begin{matrix} \frac{d \overset{⟶}{β}}{d t} = \frac{- e}{γ m_{e} c} (- \overset{⟶}{β} (\overset{⟶}{β} \cdot \overset{⟶}{E}) + \overset{⟶}{E} + c \overset{⟶}{β} \times \overset{⟶}{B}) . \end{matrix}

Here, q = −e is the electron’s charge, the Cartesian coordinates and normalized velocity are denoted by $\overset{⟶}{x}$ and $\overset{⟶}{β} \equiv \overset{⟶}{v} / c$ , while the Lorentz factor is $γ = 1 / \sqrt{1 - {(\overset{⟶}{β})}^{2}}$ .

The four-momentum of the electron is $p^{μ} \equiv (p^{0}, \overset{⟶}{p}) = m_{e} γ c (1, \overset{⟶}{β})$ , where $m_{e} \equiv \sqrt{p^{μ} p_{μ} / c^{2}} = 0.511 MeV / c^{2}$ is its invariant rest mass. The relativistic energy $E$ and momentum $\overset{⟶}{p}$ of electrons are expressed as

(10)

\begin{matrix} p^{0} \equiv \frac{E}{c} = γ m_{e} c, \\ \overset{⟶}{p} = γ m_{e} \overset{⟶}{v} . \end{matrix}

Laser-driven electron acceleration in vacuum is the consequence of the direct interaction of the laser pulse with electrons [Reference Hartemann, Fochs and Le Sage37–Reference He, Yu and Lu46]. At any given time, equations (8) and (9) are input for the electromagnetic field of the laser pulse, i.e., equations (5) and (6); hence, the trajectory of the propagating electric charge dynamically maps the laser pulse.

The 3-dimensional solutions to the electron trajectories and velocities are obtained by solving these coupled differential equations numerically by an adaptive time-step Runge–Kutta method with an accuracy and numerical precision up to 12 digits.

2.4. Initial Conditions

Unless stated otherwise, initially, all electrons are at rest, i.e., ${\overset{⟶}{β}}_{i} (t_{0}) = {\overset{⟶}{β}}_{0, i} \equiv \overset{⟶}{0}$ and $γ_{0, i} = 1$ . These electrons at z _0,i = 0 coordinates are uniformly distributed in the orthogonal plane, $(x_{0, i}, y_{0, i})$ , on a disk with a radius that is three times the beam waist radius, i.e., r ₀ = 3w ₀. Thus, initially, over 99% of the laser’s energy is contained within this disk. The initial position of the peak of the laser pulse is located on the longitudinal axis at $z_{F} = - 5 τ_{0} c$ behind the electrons, i.e., full pulse interaction, while all electrons are independent from each other and only interact with the laser pulse [Reference Molnar, Stutman and Ticos27].

For current purposes, we have fixed the laser wavelength to λ₀ = 800 nm. The laser pulse duration at full width at half maximum (FWHM), Δτ ₀ = 25 fs, corresponds to $τ_{0} = Δ τ_{0} / 2 \sqrt{ln 2}$ . Similarly, the beam waist at FWHM Δw ₀ leads to $w_{0} = Δ w_{0} / 2 \sqrt{ln 2}$ waist radius.

The peak intensity and peak power for a monochromatic LP Gaussian laser are $I_{0} \equiv a_{0}^{2} {(m_{e} c ω_{0} / e)}^{2} (c ε_{0} / 2)$ and $P_{0} \equiv I_{0} (π w_{0}^{2} / 2)$ , where the normalized electric field amplitude is $a_{0} = (e E_{0} / m_{e} c ω_{0})$ . For P ₀ = {0.1, 1,10} PW power and beam waists of Δw ₀ = {6,11,19}λ₀, these values are listed in Table 1. For helical beams, the field intensities at the local maxima are $a_{0}^{0 m} (r_{m}) = a_{0} \sqrt{C_{0 m} {(m)}^{(m)} e^{- (m)}}$ since the total power of the beam is constant, $P_{0} \equiv (E_{0}^{2} / m!) \int_{0}^{\infty} 2 π r {(\sqrt{2} r / w_{0})}^{2 (m)} \exp (- 2 r^{2} / w_{0}^{2}) d r = E_{0}^{2} w_{0}^{2} π / 2$ .

Table 1:The normalized field amplitudes corresponding to LP Gaussian pulses for different waist radii and laser powers.

3. Results

3.1. Optimal Beam Waist

First of all, we are interested in the values of the beam waist for different laser powers that lead to maximal energy gains, in case of the Gaussian beam LG₀₀ and the helical beams LG₀₁, LG₀₂, LG₀₃, LG₀₄, and LG₀₅.

The energy gain of electrons interacting with an LG_0m laser pulse is a function of the initial location of electrons, the laser spot size, and azimuthal mode m. To estimate the value of the beam waist that correspondingly leads to maximum energy gains for a given laser power, we have varied the initial position of electrons uniformly $x_{0, i} = y_{0, i} \equiv (0.05, 0.1, 0.15, \dots, 2.5) w_{0}$ at z _0,i = 0 and calculated the respective energy gain of a single electron for each position:

(11)

\begin{matrix} Δ E_{i} (t) \equiv E_{i} (t) - E_{i} (t_{0}) = (γ_{i} - 1) m_{e} c^{2}, \end{matrix}

where $E_{i} (t_{0}) = m_{e} c^{2}$ . Furthermore, for all discrete values of beam waist Δw ₀ = {1,2, …, 200}λ₀, we have also calculated the average of these energy gains, $Δ E = \sum_{i}^{N} Δ E_{i} / N$ , where N = 50 corresponds to the number of the initial positions of electrons in the transverse plane.

These weighted averages for different laser powers are shown in Figure 2 for both LP and CP pulses. Here, the average energy gains corresponding to laser powers of P ₀ = {0.1, 1,10} PW are plotted with red, blue, and black, as a function of initial spot size. The full and dashed lines correspond to CP and LP lasers, respectively.

Figure 2:

The average energy gain of an electron initially located at z ₀ = 0 and x ₀ = y ₀ ≡ {0.05, 0.1, …,2.5} Δw ₀ as a function of the beam waist. The red, blue, and black lines correspond to the weighted average energy of P ₀ = {0.1, 1,10} PW laser powers, respectively. The full and dotted lines correspond to CP and LP lasers. (a–c) The Gaussian beam, LG₀₀, and the helical beams LG₀₁ and LG₀₂. (d–f) The helical beams LG₀₃, LG₀₄, and LG₀₅.

First, we discuss these results qualitatively. For a given laser power, the amplitude of the electron oscillations along the polarization direction increases with intensity. Increasing the beam waist decreases the intensity, and the scattering of electrons decreases; hence, they remain confined in the pulse being able to gain more energy from the laser, until the oscillations become larger than the waist and the electron scatters out from the pulse.

On the contrary, increasing the waist of the beam also decreases the longitudinal components of the electric field and the Lorentz force, i.e., E_z and ${(\overset{⟶}{v} \times \overset{⟶}{B})}_{z}$ , and therefore reduces the net kinetic energy gain of electrons. The electrons are accelerated to larger and larger velocities in the front part of the pulse, and thus, the electron trajectories are elongated in the direction of the laser propagation, while the deceleration in the back part of the pulse becomes less efficient.

This means that the beam waist corresponding to the highest average energy gain, i.e., the highest peaks of the averages in Figure 2, represents the optimal waist for the given laser power and polarization. For larger power lasers, a wider initial beam waist is more optimal to ensure that the electrons remain confined inside the pulse to gain more energy.

Therefore, using Figure 2, we can approximate the beam waists corresponding to the peaks in net energy gain. In case of the fundamental Gaussian beam, the optimal beam waist at FWHM corresponds to few tens of wavelengths, Δw ₀ = {13,23,41}λ₀, for laser powers of P ₀ = {0.1, 1,10} PW, such that the net energy gain increases about $Δ E \approx \sqrt{10}$ for every order of magnitude increase in laser power. Even with those optimal values, the average energy gain of only a few MeV was observed in LG₀₀ beams; see [Reference Molnar, Stutman and Ticos27] for more details. Note also that, in case of the Gaussian pulse, the outcome is independent on the polarization, and the averaged results overlap. This behavior is similar for the LG₀₁ helical beam, but for higher modes, the distinction between CP and LP pulses becomes more apparent.

Now, averaging the energy gains $Δ E$ of the helical beams LG₀₁ and LG₀₂ for a circularly polarized laser from Figure 2, we have approximated the optimal beam waists at FWHM. These optimal beam waists are Δw ₀ = {6,11,19}λ₀ corresponding to increasing laser power of P ₀ = {0.1, 1,10} PW. These approximated values represent the chosen optimal values for direct electron acceleration in all helical beams of interest, LG₀₁,…, LG₀₅. Note, however, that the optimal beam waist also reduces slightly as the mode number increases (see Figure 2). Therefore, the previously chosen beam waists are suboptimal for LG₀₄ and LG₀₅ modes where even a smaller beam waist would be more favorable for larger energy gains.

The most important observation in case of CP beams is that the optimal spot sizes of helical beams are more than twice smaller than for the fundamental Gaussian beam. In case of LP beams, this difference further increases with the mode index |m| ≥ 1. Furthermore, for these relatively tight initial waists, of a few wavelengths, the helical beams might lead to almost an order of magnitude larger net energy gains compared to the fundamental Gaussian beam at high laser power.

A straightforward explanation can be formulated in terms of the ponderomotive force, $F_{P}^{0 m} = - \nabla Φ^{0 m}$ , where the ponderomotive potential, the cycle-averaged oscillation energy, is directly proportional to the intensity, $Φ^{0 m} ≃ e^{2} {({\overset{⟶}{E}}_{0 m})}^{2} / (4 m_{e} ω_{0}^{2}) \sim I_{0}^{0 m}$ , shown in Figure 1. The ponderomotive force causes the charges oscillating in an inhomogeneous electric field to drift from where the electric field is larger to where it is smaller. Therefore, the immediate consequence of the transverse ponderomotive force is the scattering of charges from “regions” of higher to lower electric field intensity. Due to the fact that the intensity of the helical beams has a wide convex region, the charges found in this region are naturally driven to the beam center with zero intensity.

In other words, the transverse ponderomotive force in Gaussian beams is always positive; hence, the electrons are scattered outwards from the pulse. In helical beams, the transverse ponderomotive force of the convex region is negative, and hence, the electrons are effectively trapped inside the “hollow” pulse (see Figure 1). This leads to less spread, better focusing, and collimated electron trajectories confined near the axis of propagation. These captured electrons are accelerated further by the longitudinal ponderomotive force while continue to gain more energy through the phase synchronization process leading to larger energy gains; see [Reference Akou and Firouzjaei19] for more details.

3.2. Energy Gain for Optimal Beam Waists

Using the previously given initial conditions together with the optimal values listed in Table 1 for both CP and LP lasers of P ₀ = {0.1, 1,10} PW power and initial beam waist of Δw ₀ = {6,11,19}λ₀, we have numerically calculated the direct laser-driven electron acceleration corresponding to the fundamental Gaussian LG₀₀ and LG_0m helical beams with azimuthal modes m = {1,2,3,4,5}.

In Figure 3, we have plotted all Cartesian components of the electric field as a function of time in units of the normalized field intensity a ₀ as mapped by N_e = 6000 accelerated electrons with lower index i. In LG₀₀ beams, the electrons are accelerated at the front part of the pulse, but without reaching the available peak intensity of the pulse, the electrons are scattered out by the transverse ponderomotive force. In helical beams, e.g., LG₀₁ and LG₀₅, the available “peak” intensities are about 2.5 and 5.5 times smaller than in LG₀₀. However, these intensity maxima are located further away from the midpoint, i.e., r = 0, the location of maximum intensity of the Gaussian beam and zero intensity for the helical beams (see Figure 1). Therefore, the electrons found in the convex part of the pulse are captured and accelerated for a much longer time, being clearly visible in Figure 3, where the time axes of the electric field components are an order of magnitude longer for the helical beams. This also means that the trajectories of some electrons captured by helical beams are up to an order of magnitude longer in the direction of the laser propagation. Furthermore, as apparent in Figure 3, the electrons also have a twisted circular motion with intertwining trajectories, similar as presented in [Reference Ju, Zhou and Jiang9].

Figure 3:

CP lasers of P ₀ = 1 PW power, with Δτ ₀ = 25 fs and Δw ₀ = 11 λ₀. (a0, b0, c0) The electric fields seen by the electrons E_x,i, E_y,i, and E_z,i in units of a ₀ as a function of time, corresponding to the LG₀₀ beam. (a1, b1, c1) They are for the helical beam LG₀₁.

To further elucidate these issues, in Figure 4, we have plotted the net energy gain of electrons as a function of the initial radial distance, $r_{0, i} = \sqrt{x_{0, i}^{2} + y_{0, i}^{2}}$ , at the origin z _i,0 = 0, for varying laser power. Here, all figures correspond to an initial beam waist of Δw ₀ = 11λ0, while CP and LP pulses are plotted with “o” and “x,” respectively. Similar to the previously presented figures, red, blue, and black correspond to P ₀ = {0.1, 1,10} PW laser power. Figures 4(a0) and 4(b0) correspond to Gaussian beams, while Figures 4(a1) and 4(b1) are for helical beams LG₀₁. Similarly, the helical beams, LG₀₂, LG₀₃ and LG₀₄, LG₀₅, with different polarizations are shown.

Figure 4:

The energy gained from laser pulses with Δτ₀ = 25 fs and Δw ₀ = 11λ₀ of P ₀ = {0.1, 1,10} PW power with red, blue, and black correspondingly as a function of the same initial radial position of electrons. (a0, b0) The energy gained from a Gaussian laser pulse corresponding to CP and LP pulses, which are plotted with “o” and “x,” respectively. (a1, b1) The energy gained as a function of the initial radial position in LG₀₁ helical CP and LP beams, (a2, b2) LG₀₂ and (a3, b3) LG₀₃ helical CP and LP beams, and (a4, b4) LG₀₄ and (a5, b5) LG₀₅ helical CP and LP beams.

Here, we observe that the net energy gain as a function of the initial radial distance also reflects the initial intensity profiles shown in Figure 1. The intensity of the Gaussian beam falls off exponentially as a function of the radius; hence, the electric charges found further away from the center gain less and less energy. In the case of the helical beams, the largest acceleration occurs within the convex part of the intensity curves. Furthermore, as the distance between the center and the intensity peak widens with the azimuthal index, the peaks in intensity of the LG_0m helical beams are also decreasing (see Figure 1). This explains why there is less and less net energy gained around the middle of the helical beams with increasing mode index m (see Figure 4).

In Figure 5, the histograms of the energy $Δ E_{i}$ and the polar angle ϕ_i of electrons are shown after the interaction with the pulse. Here, the 3-dimensional polar angle is $ϕ_{i} = arccos (z_{i} / R_{i})$ , where $R_{i} = \sqrt{x_{i}^{2} + y_{i}^{2} + z_{i}^{2}}$ is the radial distance from the origin. All histograms correspond to a CP laser of P ₀ = 10 PW power and Δw ₀ = 11λ₀ spot. In Figure 5, (a0), (b0) and (a1, b1) show the histograms of energy and polar angle corresponding to LG₀₀ and LG₀₁ beams. Similarly, Figures (a2), (b2), (a3), (b3) and (a4), (b4), (a5), (b5) show the outcome from the helical beams LG₀₂, LG₀₃ and LG₀₄, LG₀₅.

Figure 5:

The logarithmic scale histogram of net energy gained $Δ E_{i}$ and the polar angle histogram ϕ_i of electrons after interaction with a circularly polarized laser pulse with Δτ ₀ = 25 fs, Δw ₀ = 11λ₀ beam waist, and P ₀ = 10 PW power. (a0, b0) The energy gained from a Gaussian pulse and the corresponding angular distribution of electrons. (a1, b1) The energy gain and angular distribution in case of the LG₀₁ helical beam, for helical beams LG₀₂ (a2, b2) and LG₀₃ (a3, b3), and for helical beams LG₀₄ (a4, b4) and LG₀₅ (a5, b5).

These histograms once again reflect the difference in the ponderomotive force and its influence on the energy gain as well as the angular distribution of electrons. Due to the interaction with the pulse, the electrons have scattered out with ϕ_i polar angle. In case of the fundamental Gaussian beam, this polar angle is predominantly in the direction orthogonal to the direction of laser propagation, i.e., ϕ_i ≥ 80 degrees. For helical beams with increasing azimuthal mode, we observe an increasing number of electrons that are scattered parallel to the longitudinal axis, at polar angles of $0 \leq ϕ_{i} < 15$ degree, and, at the same time, fewer electrons in the orthogonal directions. This obviously means that helical beams lead to a larger number of collimated electrons than fundamental Gaussian beams, while their number is also increasing with increasing mode index.

Furthermore, the distribution of energy gain is also different in helical beams. The electrons are distributed following a multimodal distribution, leading to at least a second peak at finite energy at about 150 MeV energy for the helical beams, LG₀₃, LG₀₄, and LG₀₅. Although, here, we have only shown the results for the P ₀ = 10 PW laser, a very similar behavior is observed for the P ₀ = 1 PW laser, with a peak at about 50 MeV energy. For a P ₀ = 0.1 PW laser, the intensities of the helical beams are so weak such that these striking differences disappear.

In Figures 6–8 corresponding to initial beam waists of Δw ₀ = 6λ₀, Δw ₀ = 11λ₀, and Δw ₀ = 19λ₀, we have plotted the average net energy and the largest energy gained, $Δ E_{m}$ , both from circularly and linearly polarized LG_0m pulses. Here, the average energy is calculated from $Δ E = (1 / N_{e}) \sum_{i = 1}^{N_{e}} Δ E_{i}$ , where N_e = 6000 is the number of electrons. These averages are shown with full lines for circularly polarized and with dotted lines for linearly polarized lasers. Similar as before, red, blue, and black, respectively, represent the power of the laser P ₀ = {0.1, 1,10 PW. In Figures 6–8, (a0) and (b0) show the average energy gained within the fundamental Gaussian beams, while Figures (a1) and (b1) show the same for the helical beams LG₀₁ with given waist. Figures (a2) and (b2) correspond to LG₀₂, while Figures (a3) and (b3) correspond to LG₀₃. Similarly, Figures (a4), (b4) and (a5), (b5) present the outcome for LG₀₄ and LG₀₅ helical beams.

Figure 6:

Lasers with Δτ ₀ = 25 fs pulse duration and Δw ₀ = 6λ₀ beam waist radius. The full and dotted lines correspond to a CP and LP laser of P ₀ = {0.1, 1,10} PW power, with red, blue, and black, respectively. (a0) The time evolution of the mean net energy gain of electrons; (b0) the time evolution of the highest energy electron, $Δ E_{m}$ ; both correspond to Gaussian pulses of different powers and polarizations. (a1, b1) The energy gains corresponding to the LG₀₁ helical beam, (a2, b2) LG₀₂ and (a3, b3) LG₀₃ helical beams, and (a4, b4) LG₀₄ and (a5, b5) LG₀₅ helical beams.

Figure 7:

Similar to Figure 6. All figures correspond to lasers with Δτ ₀ = 25 fs and Δw ₀ = 11λ₀ beam waist radius. The full and dotted lines correspond to a CP and LP laser of P ₀ = {0.1, 1,10} PW power, with red, blue, and black, respectively.

Figure 8:

Similar to Figures 6 and 7. All figures correspond to lasers with Δτ ₀ = 25 fs and Δw ₀ = 19λ₀ beam waist radius. The full and dotted lines correspond to a CP and LP laser of P ₀ = {0.1, 1,10} PW power, with red, blue, and black, respectively.

We observe that the mean energy gain of electrons for any given laser power is largest for an optimal waist, and this optimal waist is increasing with increasing laser power. Therefore, for a P ₀ = 0.1 PW laser, the optimal waist size is Δw ₀ = 6λ₀, compare the mean energies, i.e., the red lines, in Figures 6–8. Similarly, by comparing the blue lines, for a P ₀ = 1 PW laser, the waist size that leads to the largest energy gains is Δw ₀ = 11λ₀, while for a P ₀ = 10 PW laser, the optimal waist size is Δw ₀ = 19λ₀, in accordance with the results of Section 3.1.

On the contrary, we also recognize that the average energy gain in helical beams is increasing with increasing azimuthal index $(m) \geq 1$ . This is evident by comparing the blue and black lines in Figures 6–8. For a larger initial waist, there is larger energy gain for higher laser power. However, in Figure 8, we also observe that, for the largest waist and largest laser power, the average energy is increasing up to $(m) \leq 3$ , while it is decreasing for higher modes. This was somewhat expected based on Figure 2, meaning that, for the highest power lasers with $(m) \geq 3$ modes, the optimal waist should be smaller than $Δ w_{0} \leq 19 λ_{0}$ .

The effect of laser polarization on energy gain is a delicate matter as already suggested by Figure 2. The average net energy gain in fundamental Gaussian beams is largely independent on the polarization. However, for helical beams, the outcome is different in each case depending on the beam waist and mode index, and therefore, the difference is from a few percent to up to 50%. We have also observed that a larger average energy gain slightly favors the circularly polarized helical beams with mode indices $(m) \geq 3$ for any beam waist and laser power. For linearly polarized helical beams, the optimal waist size should be slightly smaller than in the circularly polarized case to obtain similar energy gains, as can be concluded from the results.

The average net energy gain in Gaussian pulses, for the highest laser power and largest waist, is less than 2 MeV, while the highest energy electrons are $E_{\max} \approx 150$ MeV (see Figure 8). As previously discussed, the energy gains in fundamental Gaussian pulses may be a little higher for optimal beam waists [Reference Molnar, Stutman and Ticos27], but even so, these values are easily surpassed in helical beams leading to an average of about a few and up to 45 MeV, and, in some cases, with highest energy electrons of $E_{\max} \approx 650$ MeV or more. Such energy gains are still below than the so-called ponderomotive limit $Δ E \approx m_{e} c^{2} a_{0}^{2} / 2$ [Reference Stupakov and Zolotorev21, Reference Dodin and Fisch47], where the intensity maxima of higher-order helical beams can be more than 5 times lower than those of Gaussian beams.

4. Conclusions

In this paper, we have studied and compared the direct laser-driven electron acceleration in Laguerre–Gaussian beams with azimuthal mode indices, LG₀₀ corresponding to the fundamental Gaussian beam, and helical beams LG₀₁, LG₀₂, LG₀₃, LG₀₄, and LG₀₅.

We have found that the acceleration of electrons from rest is vastly different in helical beams compared to the fundamental Gaussian beam, mainly due to the difference in ponderomotive forces. Most importantly, for relatively tight initial waists, the helical beams lead to at least an order of magnitude larger energy gain compared to the fundamental Gaussian beam, and this energy gain also increases with the azimuthal mode $(m) \geq 1$ .

For laser powers of P ₀ = {0.1, 1,10} PW, the optimal waist of helical beams leading to the most energetic electrons is more than twice smaller, Δw ₀ = {6,11,19}λ₀, than in the case of the fundamental Gaussian beam [Reference Molnar, Stutman and Ticos27]. This also means that the beam waist that is optimal for helical beams is suboptimal for the Gaussian beam, and vice versa.

Finally, compared to Gaussian beams, the electron trajectories in helical beams are confined in the direction of the laser propagation leading to collimated electrons with trajectories that are at least an order of magnitude longer in the direction of the laser propagation. These conclusions hold for both linearly and circularly polarized lasers.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors are thankful to S. Ataman, D. Doria, J. F. Ong, K. Tanaka, and S. Tzenov for corrections and suggestions. E. Molnár thanks H. S. Ghotra for the comparison of early results and valuable discussions. D. Stutman acknowledges the support by a grant of the Ministry of Education and Research, CNCS-UEFISCDI (project no. PN-IIIP4-ID-PCCF-2016–0164), within PNCDI III. The authors are thankful for the financial support from the Nuclei Project (PN 19060105).

References

Allen, L., Beijersbergen, M. W., Spreeuw, R. J. C., and Woerdman, J. P., “Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes,” Physical Review A, vol. 45, no. 11, pp. 8185–8189, 1992.CrossRef Google Scholar PubMed

Grier, D. G., “A revolution in optical manipulation,” Nature, vol. 424, no. 6950, pp. 810–816, 2003.CrossRef Google Scholar PubMed

Gibson, G., Courtial, J., Padgett, M. J. et al., “Free-space information transfer using light beams carrying orbital angular momentum,” Optics Express, vol. 12, no. 22, pp. 5448–5456, 2004.CrossRef Google Scholar PubMed

Torres, J. P. and Torner, L., Twisted Photons: Applications of Light with Orbital Angular Momentum, Wiley VCH, Weinheim, Germany, 2011.CrossRef Google Scholar

Wang, J., Yang, J.-Y., Fazal, I. M. et al., “Terabit free-space data transmission employing orbital angular momentum multiplexing,” NaturePhoton, vol. 6, no. 7, pp. 488–496, 2012.Google Scholar

Krenn, M., Handsteiner, J., Fink, M. et al., “Twisted light transmission over 143 km,” Proceedings of the National Academy of Sciences, vol. 113, no. 48, pp. 13648–13653, 2016.CrossRef Google Scholar PubMed

Zhang, X., Shen, B., Yin, S. et al., “Generation of intense high-order vortex harmonics,” Physical Review Letters, vol. 114, Article ID 173901, 2015.CrossRef Google Scholar PubMed

Ju, L. B., Huang, T. W., Xiao, K. D. et al., “Controlling multiple filaments by relativistic optical vortex beams in plasmas,” Physical Review E, vol. 94, Article ID 033202, 2016.CrossRef Google Scholar PubMed

Ju, L. B., Zhou, C. T., Jiang, K. et al., “Manipulating the topological structure of ultrarelativistic electron beams using Laguerre-Gaussian laser pulse,” New Journal of Physics, vol. 20, Article ID 063004, 2018.CrossRef Google Scholar

Baumann, C. and Pukhov, A., “Electron dynamics in twisted light modes of relativistic intensity,” Physics of Plasmas, vol. 25, no. 8, Article ID 083114, 2018.CrossRef Google Scholar

Hu, L.-X., Yu, T.-P., Lu, Y. et al., “Dynamics of the interaction of relativistic Laguerre-Gaussian laser pulses with a wire target,” Plasma Physics and Controlled Fusion, vol. 61, Article ID 025009, 2019.CrossRef Google Scholar

Chen, Y.-Y., Hatsagortsyan, K. Z., and Keitel, C. H., “Generation of twisted γ-ray radiation by nonlinear Thomson scattering of twisted light,” Matter and Radiation at Extremes, vol. 4, no. 2, Article ID 024401, 2019.CrossRef Google Scholar

Luís Martins, J., Vieira, J., Ferri, J., and Fülöp, T., “Radiation emission in laser-wakefields driven by structured laser pulses with orbital angular momentum,” Scientific Reports, vol. 9, no. 1, p. 9840, 2019.CrossRef Google Scholar PubMed

Zhang, G.-B., Chen, M., Luo, J. et al., “Acceleration of on-axis and ring-shaped electron beams in wakefields driven by Laguerre-Gaussian pulses,” Journal of Applied Physics, vol. 119, no. 10, Article ID 103101, 2016.CrossRef Google Scholar

Nuter, R., Korneev, P., Dmitriev, E., Thiele, I., and Tikhonchuk, V. T., “Gain of electron orbital angular momentum in a direct laser acceleration process,” Physical Review E, vol. 101, Article ID 053202, 2020.CrossRef Google Scholar

Blackman, D. R., Nuter, R., Korneev, P., and Tikhonchuk, V. T., “Kinetic plasma waves carrying orbital angular momentum,” Physical Review E, vol. 100, Article ID 013204, 2019.CrossRef Google Scholar PubMed

Blackman, D. R., Nuter, R., Korneev, P., and Tikhonchuk, V. T., “Nonlinear Landau damping of plasma waves with orbital angular momentum,” Physical Review E, vol. 102, Article ID 033208, 2020.CrossRef Google Scholar PubMed

Vaziri, M., Golshani, M., Sohaily, S., and Bahrampour, A., “Electron acceleration by linearly polarized twisted laser pulse with narrow divergence,” Physics of Plasmas, vol. 22, Article ID 033118, 2015.CrossRef Google Scholar

Akou, H. and Firouzjaei, A. S., “Direct electron bunch acceleration by Laguerre-Gauss laser pulse direct electron bunch acceleration by Laguerre-Gauss laser pulse,” Physics of Plasmas, vol. 27, no. 9, Article ID 093102, 2020.CrossRef Google Scholar

Chaloupka, J. L., Fisher, Y., Kessler, T. J., and Meyerhofer, D. D., “Single-beam, ponderomotive-optical trap for free electrons and neutral atoms,” Optics Letters, vol. 22, no. 13, pp. 1021–1023, 1997.CrossRef Google Scholar PubMed

Stupakov, G. V. and Zolotorev, M. S., “Ponderomotive laser acceleration and focusing in vacuum for generation of attosecond electron bunches,” Physical Review Letters, vol. 86, no. 23, pp. 5274–5277, 2001.CrossRef Google Scholar PubMed

Kong, Q., Miyazaki, S., Kawata, S. et al., “Electron bunch trapping and compression by an intense focused pulse laser,” Physical Review E, vol. 69, Article ID 056502, 2004.CrossRef Google Scholar

Miyazaki, S., Kawata, S., Kong, Q. et al., “Generation of a microelectron beam by an intense short pulse laser in the TEM(1, 0) + TEM(0, 1) mode in vacuum,” Journal of Physics D: Applied Physics, vol. 38, no. 11, pp. 1665–1673, 2005.CrossRef Google Scholar

Li, Y. Y., Gu, Y. J., Zhu, Z. et al., “Direct laser acceleration of electron by an ultra intense and short-pulsed laser in under-dense plasma,” Physics of Plasmas, vol. 18, Article ID 053104, 2011.CrossRef Google Scholar

Abramochkin, E. and Volostnikov, V., “Beam transformations and nontransformed beams,” Optics Communications, vol. 83, no. 1-2, pp. 123–135, 1991.CrossRef Google Scholar

Tanaka, K. A., Spohr, K. M., Balabanski, D. L. et al., “Current status and highlights of the ELI-NP research program,” Matter and Radiation at Extremes, vol. 5, Article ID 024402, 2020.CrossRef Google Scholar

Molnar, E., Stutman, D., and Ticos, C., “Optimizing direct laser-driven electron acceleration and energy gain at ELI-NP,” European Physical Journal D, vol. 74, p. 229, 2020.CrossRef Google Scholar

Siegman, A. E., Lasers, University Science Books, Revised edition, 1986.Google Scholar

Goldsmith, P. F., Quasioptical Systems: Gaussian Beam Quasioptical Propogation and Applications, Wiley-IEEE Press, Hoboken, NJ, USA, 1998.CrossRef Google Scholar

González de Alaiza Martínez, P., Duchateau, G., Chimier, B. et al., “Maxwell-consistent, symmetry- and energy-preserving solutions for ultrashort-laser-pulse propagation beyond the paraxial approximation,” Physical Review A, vol. 98, Article ID 043849, 2018.CrossRef Google Scholar

Robinson, A. P. L., Arefiev, A. V., and Neely, D., “Generating “superponderomotive” electrons due to a non-wake-field interaction between a laser pulse and a longitudinal electric field,” Physical Review Letters, vol. 111, Article ID 065002, 2013.CrossRef Google Scholar

Kirk, T., “McDonald,” https://www.physics.princeton.edu//mcdonald/examples/axicon.pdf.Google Scholar

Kirk, T., “McDonald,” https://www.physics.princeton.edu//mcdonald/examples/gaussian2.pdf.Google Scholar

Ong, J. F., Moritaka, T., and Takabe, H., “The suppression of radiation reaction and laser field depletion in laser-electron beam interaction,” Physics of Plasmas, vol. 25, Article ID 033113, 2018.CrossRef Google Scholar

Longman, A. and Fedosejevs, R., “Optimal Laguerre-Gaussian modes for high-intensity optical vortices,” Journal of the Optical Society of America A, vol. 37, no. 5, p. 841, 2020.CrossRef Google Scholar PubMed

Jackson, J. D., Classical Electrodynamics, John Wiley & Sons, New York, NY, USA, 3rd edition, 1999.Google Scholar

Hartemann, F. V., Fochs, S. N., Le Sage, G. P. et al., “Nonlinear ponderomotive scattering of relativistic electrons by an intense laser field at focus,” Physical Review E, vol. 51, no. 5, pp. 4833–4843, 1995.CrossRef Google Scholar

Hartemann, F. V., Van Meter, J. R., Troha, A. L. et al., “Three-dimensional relativistic electron scattering in an ultrahigh-intensity laser focus,” Physical Review E, vol. 58, no. 4, pp. 5001–5012, 1998.CrossRef Google Scholar

Esarey, E., Sprangle, P., and Krall, J., “Laser acceleration of electrons in vacuum,” Physical Review E, vol. 52, no. 5, pp. 5443–5453, 1995.CrossRef Google Scholar PubMed

Quesnel, B. and Mora, P., “Theory and simulation of the interaction of ultraintense laser pulses with electrons in vacuum,” Physical Review E, vol. 58, p. 3791, 1998.CrossRef Google Scholar

Salamin, Y. I. and Keitel, C. H., “Electron acceleration by a tightly focused laser beam,” Physical Review Letters, vol. 88, Article ID 095005,Google Scholar

Salamin, Y. I., “Electron acceleration from rest in vacuum by an axicon Gaussian laser beam,” Physical Review A, vol. 73, Article ID 043402, 2006.CrossRef Google Scholar

Gupta, D. N., Kant, N., Kim, D. E., and Suk, H., “Electron acceleration to GeV energy by a radially polarized laser,” Physics Letters A, vol. 368, no. 5, pp. 402–407, 2007.CrossRef Google Scholar

Zhang, X. P., Wang, W., Xie, Y. J., Wang, P. X., Kong, Q., and Ho, Y. K., “Field properties and vacuum electron acceleration in a laser beam of high-order Laguerre-Gaussian mode,” Optics Communications, vol. 281, no. 15-16, pp. 4103–4108, 2008.CrossRef Google Scholar

Fortin, P.-L., Pieche, M., and Varin, C., “Direct-field electron acceleration with ultrafast radially polarized laser beams: scaling laws and optimization,” Journal of Physics. B, Atomic, Molecular and Optical Physics, vol. 43, Article ID 025401, 2010.CrossRef Google Scholar

He, F., Yu, W., Lu, P. et al., “Ponderomotive acceleration of electrons by a tightly focused intense laser beam,” Physical Review E, vol. 68, Article ID 046407, 2003.CrossRef Google Scholar PubMed

Dodin, I. Y. and Fisch, N. J., “Relativistic electron acceleration in focused laser fields after above-threshold ionization,” Physical Review E, vol. 68, Article ID 056402, 2003.CrossRef Google Scholar PubMed

Figure 1: The relative intensity profiles at constant power and beam waist as a function of the radius of an LP laser: the fundamental Gaussian beam LG00 with dotted black line and the helical beams LG01, LG02, LG03, LG04, and LG05 with green, magenta, red, blue, and black lines correspondingly.

Table 1: The normalized field amplitudes corresponding to LP Gaussian pulses for different waist radii and laser powers.

Figure 2: The average energy gain of an electron initially located at z0 = 0 and x0 = y0 ≡ {0.05, 0.1, …,2.5} Δw0 as a function of the beam waist. The red, blue, and black lines correspond to the weighted average energy of P0 = {0.1, 1,10} PW laser powers, respectively. The full and dotted lines correspond to CP and LP lasers. (a–c) The Gaussian beam, LG00, and the helical beams LG01 and LG02. (d–f) The helical beams LG03, LG04, and LG05.

Figure 3: CP lasers of P0 = 1 PW power, with Δτ0 = 25 fs and Δw0 = 11 λ0. (a0, b0, c0) The electric fields seen by the electrons Ex,i, Ey,i, and Ez,i in units of a0 as a function of time, corresponding to the LG00 beam. (a1, b1, c1) They are for the helical beam LG01.

Figure 4: The energy gained from laser pulses with Δτ0 = 25 fs and Δw0 = 11λ0 of P0 = {0.1, 1,10} PW power with red, blue, and black correspondingly as a function of the same initial radial position of electrons. (a0, b0) The energy gained from a Gaussian laser pulse corresponding to CP and LP pulses, which are plotted with “o” and “x,” respectively. (a1, b1) The energy gained as a function of the initial radial position in LG01 helical CP and LP beams, (a2, b2) LG02 and (a3, b3) LG03 helical CP and LP beams, and (a4, b4) LG04 and (a5, b5) LG05 helical CP and LP beams.

Figure 5: The logarithmic scale histogram of net energy gained ΔEi and the polar angle histogram ϕi of electrons after interaction with a circularly polarized laser pulse with Δτ0 = 25 fs, Δw0 = 11λ0 beam waist, and P0 = 10 PW power. (a0, b0) The energy gained from a Gaussian pulse and the corresponding angular distribution of electrons. (a1, b1) The energy gain and angular distribution in case of the LG01 helical beam, for helical beams LG02 (a2, b2) and LG03 (a3, b3), and for helical beams LG04 (a4, b4) and LG05 (a5, b5).

Figure 6: Lasers with Δτ0 = 25 fs pulse duration and Δw0 = 6λ0 beam waist radius. The full and dotted lines correspond to a CP and LP laser of P0 = {0.1, 1,10} PW power, with red, blue, and black, respectively. (a0) The time evolution of the mean net energy gain of electrons; (b0) the time evolution of the highest energy electron, ΔEm; both correspond to Gaussian pulses of different powers and polarizations. (a1, b1) The energy gains corresponding to the LG01 helical beam, (a2, b2) LG02 and (a3, b3) LG03 helical beams, and (a4, b4) LG04 and (a5, b5) LG05 helical beams.

Figure 7: Similar to Figure 6. All figures correspond to lasers with Δτ0 = 25 fs and Δw0 = 11λ0 beam waist radius. The full and dotted lines correspond to a CP and LP laser of P0 = {0.1, 1,10} PW power, with red, blue, and black, respectively.

Figure 8: Similar to Figures 6 and 7. All figures correspond to lasers with Δτ0 = 25 fs and Δw0 = 19λ0 beam waist radius. The full and dotted lines correspond to a CP and LP laser of P0 = {0.1, 1,10} PW power, with red, blue, and black, respectively.

Article contents

Direct Laser-Driven Electron Acceleration and Energy Gain in Helical Beams

Abstract

Information

1. Introduction

2. Direct Laser-Driven Electron Acceleration in Helical Beams

2.1. Laguerre–Gaussian Beams

2.2. Helical Beams

2.3. Electron Acceleration in the Electromagnetic Field

2.4. Initial Conditions

3. Results

3.1. Optimal Beam Waist

3.2. Energy Gain for Optimal Beam Waists

4. Conclusions

Data Availability

Conflicts of Interest

Acknowledgments

References

Save article to Kindle

Save article to Dropbox

Save article to Google Drive

Reply to: Submit a response

Your details

You have entered the maximum number of contributors

Conflicting interests