This work presents a theoretical analysis of the generation of twisted terahertz (THz) radiation using laser-bunched relativistic electron beams in a magnetic wiggler. By employing a laser-bunched relativistic electron beam, which introduces a transverse modulation to the electron beam density, and a magnetic wiggler, which induces a transverse deflection to the electron trajectories, the generation of twisted THz radiation is achieved. The interaction between the modulated electron beam and the magnetic field leads to the emission of THz photons with a twisted phase structure. The findings of this study provide valuable insights into the generation and manipulation of twisted THz radiation contributing to the advancement of THz technology and its diverse applications.

THz radiation has garnered a great deal of attention recently due to its potential applications in spectroscopy, imaging, and communications, among other fields. Its frequency range extends from infrared to microwaves.1–3 The ability to manipulate the properties of THz radiation opens up new avenues for technological advancement. One such route is the generation of twisted THz radiation also known as vortex beams, which has an orbital angular momentum (OAM) that is precisely specified. The field of orbital angular momentum (OAM) has gained significant attention since 1992, when Allen et al.4 demonstrated that Laguerre–Gaussian (LG) beams have an orbital angular momentum of l per photon, where l is the topological charge of the OAM mode. The characteristics of such laser beams are an annular spatial profile with a central dark spot and a wavefront that spirals in a corkscrew-like manner along the beam's propagation direction.5,6 Twisted THz radiation has enormous potential for applications that rely on the special features of OAM-carrying beams.7–9 For example, in communication systems, utilizing the OAM of THz radiation can dramatically boost data capacity and improve security measures.10,11 Furthermore, the ability to adjust the OAM of THz radiation can provide increased structural and chemical research capabilities in spectroscopy and imaging.12,13 Recently, many different strategies have been presented to generate twisted THz waves.14–16 Sobhani et al. studied the cross-focusing of two twisted coaxial laser beams to create a resonant vortex THz beam. Vortex THz radiation can be produced more effectively by adjusting the laser and plasma parameters. It is also shown that as the beating frequency approaches the plasma frequency, the THz radiation amplitude increases.17 Laser-plasma vortex interactions provide a new path toward THz vortex beams, where emission qualities are tuned and angular momentum transfer is guided by the helical density.18,19 Shobani et al. further investigated the generation of twisted THz radiation by two nonlinear currents that are triggered by a static electric field. By structuring the radiation into a Laguerre–Gaussian envelope, one can create new potential for THz applications by modifying the plasma density.20 

On the other hand, it has also been shown that Bunched beams present a compelling solution to the growing need for compact and controllable THz sources.21–23 Kumar et al. provided a formalism for the controllable generation of coherent THz radiation through the use of a magnetic wiggler and a relativistic electron beam regulated by two laser beams. At around 1.9 THz, they were able to attain a maximum value of conversion efficiency of the order 10 7 .24 Hasanbeigi studied the production of coherent THz radiation from the interaction of a helical wiggler pumped with a bunched relativistic electron beam. They found that the injected beam's radius, axial velocity, and ion-channel density all increased with the highest THz power.25 Li et al. reported the simulation of a single-pass free electron laser (FEL) powered by a THz-pulse-train photoinjector, which produces a strong narrow-band THz source with a wide tuning range of frequency of radiation. According to simulation results, this FEL can provide power outputs in the few tens of megawatts and pulse energy in the several tens of micro-joules at frequencies between 0.5 and 5 THz.26 Zhen et al. suggested a technique for producing strong subpicosecond density bunching in high-intensity relativistic electron beams that is based on slice energy spread modulation. The adjustable beam's radiation can cover the frequency range of 1–10 THz with high power, as demonstrated by both theory and simulations.27 

In this research paper, we explore the generation of twisted THz radiation using a laser-bunched relativistic electron beam in a magnetic wiggler. Incorporating a magnetic wiggler further enhances the efficiency and control over the generation process enabling the precise manipulation of the radiation's OAM. Two laser pulses interact with a modulator within the electron beam path. These pulses transfer energy to the electrons through the modulator, modulating their energy levels. This energy modulation results in a corresponding modulation of their velocities, with higher energies translating to faster speeds. As a consequence, electrons with similar velocities tend to cluster together, forming bunches with higher densities compared to surrounding regions. This prebunched electron beam then enters a magnetic wiggler. The wiggler's alternating magnetic field exerts a transverse force on the electrons, causing them to oscillate perpendicular to their direction of travel. This oscillatory motion generates electromagnetic radiation in the THz range.

Section II establishes expressions for velocity modulation and density modulation of an electron beam under the influence of twisted laser beams. Section III delves into the derivation of an expression for the nonlinear current generated within this modulated beam. Section IV outlines the calculation of radiated power resulting from this interaction. Finally, Secs. V and VI, respectively, present the key results and draw conclusions from the analysis.

Let us consider two high-power vortex laser beams having frequency ω 1 and ω 2 and wave vectors k 1 and k 2, which propagates through a modulator. The electric field of the laser is assumed to be polarized along the x-axis. The radial index p and the phase singularity on axis with strength l, often known as the optical vortex charge number, can be used to characterize laser beams with LG distribution. Here, the laser beams are expected to exhibit a doughnut shape p = 0 , and the evolution of the beam width is disregarded. The electric field of such laser beam can be written as
(1)
where
where j = 1 , 2. E 0 is the initial amplitude of electric field, l j + 1 arctan Z Z R is the Gouy phase, r is the radial coordinate, r 0 represents the beam width, β z = z 1 + Z R Z 2 is the radius of curvature of the wavefront, and Z R is the Rayleigh length. Figure 1 depicts the wavefront, spatial, and phase characteristics of vortex laser beams and Gaussian beams. As the value of vortex charge number increases, a greater twist in the wavefront of the beams can be seen.
An electron beam with velocity v 0 b and density n 0 b passes through a modulator of width d. The electrons in the beam acquire oscillatory velocity due to laser when they stream axially,
(2)
where γ 0 = 1 v 0 b 2 c 2 1 / 2 .
They are subjected to ponderomotive force F P = m / 2 v 1 . v 2 * exerted by lasers at beat frequency ω = ω 1 ω 2 and wave number k = k 1 k 2,
(3)
where
where l = l 1 l 2. The beam electrons' response to the ponderomotive force in the modulator is determined by the relativistic equation of motion given as
(4)
The velocity of electron can be obtained by writing v = v 0 b z ̂ + v 1 z ̂ and γ = γ 0 + v 0 b γ 0 3 v 1 c 2 in Eq. (4). After linearizing, we get
(5)
By putting z = v 0 b ( t t 0 ) and integrating within limits t 0 to t 1 = t 0 + d / v 0 b, one can get velocity modulation,
(6)
where
(6a)
(6b)
t 0 is the time at which electron arrives at the modulator, and d is the length of the modulator. The electron beam's velocity is modulated in accordance with the energy modulation. Higher energy electrons move more quickly, while lower energy electrons move more slowly. The electron density varies periodically as a result, resulting in high and low density zones along the beam. The initial energy and velocity changes become a well-defined density modulation as the modulated electron beam propagates. This results in density bunching of electrons. The beam density modulation can be reduced to a simpler form by using the technique described in Ref. 24,
(7)
where δ = ω ( d + L ) v 0 b , t 2 = t 0 + d v 0 b + L v 0 b 1 + ξ ψ r , φ , z sin ω t 0 + θ b g is the electron's time of escape from the region of drift space and L is the length of drift space.

A prebunched electron beam is introduced into a magnetic wiggler, having magnetic field, B w = B 0 ( x ̂ + i y ̂ ) e i k w z as shown in Fig. 2.

The oscillating magnetic field of the wiggler applies a transverse force that causes the electrons to oscillate perpendicular to their direction of passage. The oscillatory velocity acquired by the wiggler is
(8)
This will give nonlinear current density at ω , k + k w,
(9)
On substituting t 2 = t z / v 0 b in Eq. (9), we get
(10)
The mathematical representation of the vector potential, as observed at a distant location, generated by the electron bunch takes the form of
(11)
For a narrow beam, one may write r r r r . z r 2 r z cos θ,
(12)
On further solving, we get
(13)
where L b is the bunch length. The magnetic field corresponding to this B r , t = × A i ω / c r ̂ × A is
(14)
The average Poynting vector is
(15)
where θ = ω L b 2 c v 0 b c cos θ + c k w ω.
The normalized THz power is
(16)
where

The governing equation of twisted THz radiation using laser-bunched relativistic electron beams in a magnetic wiggler is modeled analytically. For a given set of parameters, d = 400 μ m , L = 1 cm , v 0 b = 0.90 c , γ 0 = 2.29 , L b = 0.40 cm , L w = 0.40 cm , λ w = 1.8 mm , e E 0 / m ω 1 c = e E 0 / m ω 2 c = 0.05 , λ 1 = 1.06 μ m , λ 2 = 1.06 ( 1 + λ 1 / λ THz ) μ m , r b = 18 μ m , r 0 = 2 μ m , ω p 0 / ω 1 = 10 4 , and a range of topological charge number of both lasers equation (16) has been calculated to demonstrate the THz radiation in a helical wiggler magnetic field. Figure 3 illustrates the normalized THz power behavior as a function of normalized frequency ω / ω 1 for l 1 = l 2 = 1 , l = 0. The power initially exhibits a positive correlation with frequency, reaching a maximum value ( 2 × 10 3) at frequency 1.9 THz, but as frequency gets high, a limit sets in due to the Bessel function. As frequency increases beyond a certain point, the system becomes less efficient in converting power and the transmitted power starts to decline. By increasing the beam velocity from 0.90 c to 0.95 c, the THz power output can be enhanced by a factor of 1.9 with the peak shifting toward a higher frequency 4.2 THz. Figure 4 plots the relationship between THz radiation power and frequency ω / ω 1 for l 1 = l 2 = 2 , l = 0. The frequency and power are both normalized, making the plot more general. As the vortex charge number increases (from l 1 = l 2 = 1 to l 1 = l 2 = 2), the plot reveals that the peak power of the THz radiation gets significantly stronger. Not only does the power increase, but it also shifts toward higher frequencies. This is because at higher OAM, the intensity of the LG beam can be concentrated in a smaller region within the material. This intense region can trigger nonlinear optical effects that further enhance the generation of THz radiation. Figures 5 and 6 show the variation of normalized THz power as a function of normalized frequency ω / ω 1 for l 0. The curve in this instance is increasing exponentially because the term ψ ( r , φ , z ) in Eq. (16) has a dominant exponential component for l 0. THz power drops as beam velocity rises. THz power decreases as the number of lasers with increasing topological charge increases. This is because the relative phase between the lasers has a significant impact on the production of THz radiation for l 0. As Figs. 7 and 8 make clear, the relative Gouy phase causes a decrease in the amplitude of the THz radiation emission. The Gouy phase is the extra phase shift that a concentrated laser beam's finite size causes in various regions of the beam. A distinctive phase curvature is produced throughout the beam profile as a result. Because the Gouy phase introduces undesired phase deviations throughout the beam profile, it might interfere with phase matching. Because of this mismatch, there is less overlap between effective interaction zones, which lowers conversion efficiency and produces fewer THz. In simple terms, matching the topological numbers of the lasers maximizes their interaction and eliminates a major source of inefficiency, which results in a large increase in the amplitude of the generated THz radiation.

The generation of twisted THz radiation using a laser-bunched relativistic electron beam in a magnetic wiggler is explored in this study. Two laser pulses and a modulator are in contact within the electron beam path. These pulses modify the energy levels of the electrons by imparting energy through the modulator. This energy modulation, which also influences their velocities, converts higher energies into quicker speeds. This leads to the formation of bunches with higher densities than the surrounding areas when electrons traveling at similar speeds gather together. After being prebunched, the electron beam travels via a magnetic wiggler. The electrons oscillate perpendicular to their initial path of travel due to the transverse force of the wiggler's alternating magnetic field. This repetitive motion generates electromagnetic radiation in the THz range. In a laser-bunched relativistic electron beam system, our study uncovers the essential parameters controlling the production of twisted THz radiation, including the vital function of vortex charge number and beam velocity. As a result of this study, twisted THz radiation creation will be better understood and may be used practically in a number of scientific and technological domains. This research lays the groundwork for future developments and breakthroughs in the field of twisted THz radiation, which has a wide range of possible uses in biological spectroscopy, security screening, and THz imaging.

The authors have no conflicts to disclose.

Himani Juneja: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Software (equal); Supervision (equal); Writing – original draft (equal). Anuraj Panwar: Supervision (equal). Prashant Chauhan: Supervision (equal).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
B.
Ferguson
and
X. C.
Zhang
, “
Materials for terahertz science and technology
,”
Nat. Mater.
1
,
26
33
(
2002
).
2.
E.
Pickwell
and
V. P.
Wallace
, “
Biomedical applications of terahertz technology
,”
J. Phys. D
39
,
R301
R310
(
2006
).
3.
Z.
Sheng
,
K.
Mima
,
J.
Zhang
, and
H.
Sanuki
, “
Emission of electromagnetic pulses from laser wakefields through linear mode conversion
,”
Phys. Rev. Lett.
94
,
095003
(
2005
).
4.
L.
Allen
,
M. W.
Beijersbergen
,
R. J. C.
Spreeuw
, and
J. P.
Woerdman
, “
Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes
,”
Phys. Rev. A
45
,
8185
(
1992
).
5.
A.
Allen
,
S. M.
Barnett
, and
M.
Padgett
,
Optical Angular Momentum
(
Institute of Physics
,
Bristol, UK
,
2003
).
6.
S.
Franke-Arnold
,
L.
Allen
, and
M.
Padgett
, “
Advances in optical angular momentum
,”
Laser Photonics Rev.
2
,
299
313
(
2008
).
7.
K.
Liu
,
Y.
Cheng
,
Z.
Yang
,
H.
Wang
,
Y.
Qin
, and
X.
Li
, “
Orbital-angular-momentum-based electromagnetic vortex imaging
,”
IEEE Antennas Wireless Propag. Lett.
14
,
711
714
(
2014
).
8.
M. T.
Lin
,
Y.
Gao
,
P. G.
Liu
, and
J. B.
Liu
, “
Improved OAM-based radar targets detection using uniform concentric circular arrays
,”
Int. J. Antennas Propag.
2016
,
1852659
.
9.
M. T.
Lin
,
Y.
Gao
,
P. G.
Liu
, and
J. B.
Liu
, “
Super-resolution orbital angular momentum based radar targets detection
,”
Electron. Lett.
52
,
1168
(
2016
).
10.
W.
Cheng
,
W.
Zhang
,
H.
Jing
,
S.
Gao
, and
H.
Zhang
, “
Orbital angular momentum for wireless communications
,”
IEEE Wireless Commun.
26
,
100
107
(
2019
).
11.
X. N.
Hui
,
S. L.
Zheng
,
Y. L.
Chen
,
Y. P.
Hu
,
X. F.
Jin
,
H.
Chi
, and
X. M.
Zhang
, “
Multiplexed millimeter wave communication with dual orbital angular momentum (OAM) mode antennas
,”
Sci. Rep.
5
,
10148
(
2015
).
12.
L. A.
Downes
,
D. J.
Whiting
,
C. S.
Adams
, and
K. J.
Weatherill
, “
Rapid readout of terahertz orbital angular momentum beams using atom-based imaging
,”
Opt. Lett.
47
(
22
),
6001
6004
(
2022
).
13.
S.
Fu
,
Y.
Zhai
,
J.
Zhang
et al, “
Universal orbital angular momentum spectrum analyzer for beams
,”
PhotoniX
1
,
19
(
2020
).
14.
H.
He
,
M. E.
Friese
,
N. R.
Heckenberg
, and
H.
Rubinsztein-Dunlop
, “
Direct observation of transfer of angular momentum to absorptive particles from a laser beam with a phase singularity
,”
Phys. Rev. Lett.
75
(
5
),
826
829
(
1995
).
15.
J.
He
,
X.
Wang
,
D.
Hu
,
J.
Ye
,
S.
Feng
,
Q.
Kan
, and
Y.
Zhang
, “
Generation and evolution of the terahertz vortex beam
,”
Opt. Express
21
,
20230
20239
(
2013
).
16.
R.
Imai
,
N.
Kanda
,
T.
Higuchi
,
K.
Konishi
, and
M.
Kuwata-Gonokami
, “
Generation of broadband terahertz vortex beams
,”
Opt. Lett.
39
,
3714
3717
(
2014
).
17.
H.
Sobhani
,
M.
Vaziri
,
H.
Rooholamininejad
et al, “
Terahertz twisted beams generation in plasma
,”
Eur. Phys. J D
70
,
168
(
2016
).
18.
H.
Sobhani
,
H.
Rooholamininejad
, and
A.
Bahrampour
, “
Creation of twisted terahertz waves carrying orbital angular momentum via a plasma vortex
,”
J. Phys. D
49
,
295107
(
2016
).
19.
H.
Sobhani
, “
Creation of twisted terahertz carrying orbital angular momentum via stimulated Raman scattering in a plasma vortex
,”
Laser Phys.
27
,
096001
(
2017
).
20.
H.
Shobani
,
S.
Feili
, and
M.
Taraz
, “
Generation of twisted terahertz radiation through plasma in the presence of a static electric field
,”
Laser Phys.
28
,
106001
(
2018
).
21.
G. L.
Carr
,
M. C.
Martin
,
W. R.
McKinney
et al, “
High-power terahertz radiation from relativistic electrons
,”
Nature
420
,
153
156
(
2002
).
22.
S.
Khan
,
K.
Holldack
,
T.
Kachel
,
R.
Mitzner
, and
T.
Quast
, “
Femtosecond undulator radiation from sliced electron bunches
,”
Phys. Rev. Lett.
97
,
074801
(
2006
).
23.
H.
Saberi
,
F.
Jahangiri
, and
A. R.
Niknam
, “
Terahertz wave amplification by a laser-modulated relativistic electron beam
,”
Phys. Rev. Accel. Beams
24
,
120703
(
2021
).
24.
M.
Kumar
and
V. K.
Tripathi
, “
Terahertz radiation from a laser bunched relativistic electron beam in a magnetic wiggler
,”
Phys. Plasmas
19
,
073109
(
2012
).
25.
A.
Hasanbeigi
,
H.
Mehdian
, and
P.
Gomar
, “
Enhancement of terahertz radiation power from a prebunched electron beam using helical wiggler and ion-channel guiding
,”
Phys. Plasmas
22
,
123116
(
2015
).
26.
H.
Li
,
Y.
Lu
,
Z.
He
et al, “
Generation of intense narrow-band tunable terahertz radiation from highly bunched electron pulse train
,”
J. Infrared, Millimeter, Terahertz Waves
37
,
649
657
(
2016
).
27.
Z.
Zhang
,
L.
Yan
,
Y.
Du
,
W.
Huang
,
C.
Tang
, and
Z.
Huang
, “
Generation of high-power, tunable terahertz radiation from laser interaction with a relativistic electron beam
,”
Phys. Rev. Accel. Beams
20
,
050701
(
2017
).
28.
E.
Karimi
, “
Generation and manipulation of laser beams carrying orbital angular momentum for classical and quantum information applications
,” Ph.D. thesis (
Università degli Studi di Napoli Federico II
,
2009
).