---
title: Helical Multifrequency Undulators
url: https://www.emergentmind.com/topics/helical-multifrequency-undulators
type: topic
---

# Helical Multifrequency Undulators

Searching arXiv for recent and foundational papers on helical multifrequency undulators, twisted-photon radiation, and helical undulator hardware.
A helical multifrequency undulator is a magnetic structure in which the transverse undulator field is formed as a superposition of several helical or elliptical components with different periods, amplitudes, chiralities, and phase offsets. In the modern theoretical formulation, the field is written as
\[
H_x=\sum_{i=1}^M H_x^i\sin \tilde{}_i,\qquad H_y=\sum_{i=1}^M H_y^i\cos \tilde{}_i,\qquad H_z=0,
\]
with
\[
\tilde_i=\pm \frac{2\pi z}{l_i}+\chi_i,
\]
so that each subundulator contributes its own frequency, polarization content, and phase to the overall device [2507.20242, 2602.07942]. Such systems are studied for two closely related reasons. First, helical geometry provides circular or elliptical polarization and strong transverse coupling in both \(x\) and \(y\). Second, multifrequency composition makes the radiation spectrum, harmonic content, and, in twisted-photon formulations, the projection of total angular momentum (TAM) and the relative phases between modes directly controllable by undulator design [2507.20242, 2602.07942]. At the same time, helical operation imposes nontrivial angular-redistribution constraints on harmonic emission, most notably the suppression of on-axis nonlinear harmonic generation in the ideal helical case [1105.1725].

## 1. Definition and field structure

In the multifrequency theory, the electron is ultrarelativistic and moves mainly along the \(z\)-axis through a stationary magnetic field built as a superposition of \(M\) elliptical undulator sections. The weak-angle and undulator assumptions are
\[
\gamma\gg 1,\qquad K_i/\gamma\ll 1,
\]
with
\[
K_i^2=\frac{\omega_i^2}{2}(a_i^2+b_i^2),\qquad a_i=\frac{eH_y^i}{m_e\omega_i^2},\qquad b_i=-\frac{eH_x^i}{m_e\omega_i^2},
\]
and
\[
\beta_3=1-\frac{1+K^2}{2\gamma^2},\qquad K^2=\sum_i K_i^2
\]
[2507.20242]. The central physical point is that the charged particle executes a superposition of oscillations at multiple frequencies, and the radiation field therefore acquires a correspondingly rich harmonic structure.

For the three-frequency helical specialization, the field is written as
\[
H_x=H_x^1\sin\tilde_1+H_x^2\sin\tilde_2+H_x^3\sin\tilde_3,\qquad
H_y=-H_x^1\cos\tilde_1-H_x^2\cos\tilde_2-H_x^3\cos\tilde_3,\qquad H_z=0,
\]
with
\[
a_i=b_i=r_i\qquad (i=1,2,3),
\]
so that the electron trajectory becomes
\[
x=x_0+\sum_{i=1}^3 r_i\cos_i,\qquad
y=y_0+\sum_{i=1}^3 r_i\sin_i,\qquad
z=z_0+\beta_z t+\sum_{i<j}d_{ij}\sin(\tilde_i-\tilde_j)
\]
[2602.07942]. This form shows that helical multifrequency operation is not merely a superposition of independent circular motions. Longitudinal mixing terms appear through the coefficients
\[
c_{ij}=\frac{\omega_i\omega_j}{2}\frac{a_ia_j-b_ib_j}{\omega_i+\omega_j},\qquad
d_{ij}=-\frac{\omega_i\omega_j}{2}\frac{a_ia_j+b_ib_j}{\omega_i-\omega_j},
\]
which encode inter-frequency coupling in the trajectory itself [2602.07942].

A common misconception is that any helical undulator is automatically a multifrequency undulator. The recent hardware papers on short-period NdFeB helices explicitly do not present experiments on multifrequency tuning, simultaneous multi-harmonic operation, or multiple radiation frequencies; rather, they provide fabrication, magnetization, assembly, and field characterization that are enabling and foundational for more advanced helical multifrequency concepts [2508.00232, 2509.06186].

## 2. Radiation spectrum and generalized harmonic conditions

The multifrequency radiation spectrum is governed by an energy-selection condition that generalizes the ordinary undulator harmonic relation. For twisted-photon radiation in an \(M\)-frequency undulator, the main harmonics satisfy
\[
k_0=\frac{\sum_{i=1}^M n_i\omega_i}{1-\tilde n_3\beta_3}
\]
[2507.20242]. In the three-frequency case, this becomes
\[
k_0=\frac{n_1\omega_1+n_2\omega_2+n_3\omega_3}{1-n_zz}
\]
[2602.07942]. The integers \(n_i\) play the role of generalized harmonic indices associated with exchange processes involving the different subundulator frequencies.

When the frequency ratios are rational, the spectrum reduces to an equidistant ladder. If
\[
\eta_i=\omega_i/\omega_1=\frac{h_i}{g_i},
\]
then one may introduce a common effective base frequency, and the spectrum is rewritten as
\[
k_0=\frac{\omega n}{1-\tilde n_3\beta_3},\qquad n\in\mathbb Z
\]
[2507.20242]. In the three-frequency construction the reduction is achieved by defining
\[
g=\operatorname{lcm}(g_2,g_3),\qquad \tilde g_i=g\eta_i,\qquad d=\gcd(\tilde g_1,\tilde g_2,\tilde g_3),
\]
and then
\[
\omega=\omega_1\frac{d}{g},\qquad \ell_i=\frac{\omega_i}{\omega},
\]
so that the admissible integer triples solve the Diophantine equation
\[
\ell_1 n_1+\ell_2 n_2+\ell_3 n_3=n
\]
[2602.07942].

This number-theoretic reduction is one of the distinctive features of helical multifrequency theory. It implies that several different tuples \((n_1,\ldots,n_M)\) may contribute to the same emitted energy. The radiated state at fixed energy is therefore naturally a coherent superposition rather than a single harmonic line in the one-frequency sense [2507.20242, 2602.07942].

A further notable consequence, emphasized for the rational-frequency case, is that the lowest effective frequency is smaller than the smallest subundulator frequency unless all \(\omega_i\) are integer multiples of one of them [2507.20242]. This suggests that multifrequency composition can generate spectral structure not obtainable by simply concatenating one-frequency helical sections.

## 3. Twisted photons, TAM selection rules, and composite states

The most developed contemporary application of helical multifrequency undulators is as a source of twisted photons. In the semiclassical formulation, a twisted photon is labeled by helicity \(s=\pm1\), TAM projection \(m\), longitudinal and transverse momenta \(k_3,k_\perp\), and energy \(k_0=\sqrt{k_3^2+k_\perp^2}\). The general average number of emitted twisted photons is
\[
dP(s,m,k_3,k_\perp)=e^2\left|\int d\tau\, e^{-i[k_0x^0(\tau)-k_3x_3(\tau)]} \Big\{\tfrac12[\dot x_+ a_-+\dot x_- a_+]+\dot x_3 a_3\Big\}\right|^2 n_\perp^3\,\frac{dk_3dk_\perp}{16\pi^2}
\]
with
\[
x_\pm=x\pm i y
\]
[2507.20242].

For trajectories close to the axis, the general multifrequency elliptical selection rule is
\[
m+n_1+\cdots+n_M \;\text{is even},
\]
whereas the helical specialization yields the sharper condition
\[
m=n_1+n_2+\cdots+n_M
\]
[2507.20242]. In the three-frequency helical case, the near-axis condition is stated as
\[
m=n_1+n_2+n_3
\]
[2602.07942]. The interpretation given in both works is in terms of virtual photons exchanged between the electron and the different helical components of the undulator field; the TAM transferred through those exchanges appears in the emitted twisted photon.

For rational frequencies and fixed harmonic \(n\), the TAM spectrum is constrained by the Diophantine structure. In the general helical multifrequency treatment, one obtains a generalized TAM rule of the form
\[
m=n\sum_{i=1}^M n_i^0+\sum_{a>b=1}^M(\epsilon_a-\epsilon_b)\,k_{(a,b)}
\]
[2507.20242]. In the explicitly three-frequency helical case, one representation is
\[
m=n(n_1^0+n_2^0+n_3^0)+\ell_{21}k_{21}+\ell_{31}k_{31}+\ell_{32}k_{32},
\]
and, when \(\gcd(\ell_{21},\ell_{31})=\mathrm d_{123}\),
\[
m=n(n_1^0+n_2^0+n_3^0)+\mathrm d_{123}\,r,\qquad r\in\mathbb Z
\]
[2602.07942]. Thus the allowed \(m\)-values at fixed emitted energy form an arithmetic progression.

The three-frequency theory goes beyond selection rules and treats the emitted radiation as a composite twisted state. Its explicit claim is that such undulators can generate photons in linear superpositions of modes with definite projections of the total angular momentum, amplitudes, relative phases, and polarizations, and that these quantities can be governed in a predictable way by adjusting the parameters of the multifrequency helical undulator [2602.07942]. In the two-frequency and general \(M\)-frequency formulations, the same idea appears as a coherent superposition of several TAM eigenstates at a fixed harmonic, with amplitudes controlled by \(a_i,b_i\), harmonic content controlled by \(\omega_i\), and relative phases controlled by \(\chi_i\) [2507.20242].

This perspective is closely related to earlier work showing that higher harmonics of a one-frequency helical undulator carry orbital angular momentum. That work emphasized that all harmonics with \(n>1\) carry OAM, whereas the first harmonic is the dominant ordinary helical-undulator line concentrated on axis and does not carry OAM in the same sense [1109.1603]. The multifrequency theory extends that single-frequency geometric argument into an explicit mode-engineering framework.

## 4. Angular distribution, harmonic suppression, and practical flux limits

A central constraint on helical multifrequency operation comes from the angular structure of harmonic radiation. In a helical undulator, unlike the planar case, nonlinear harmonic generation is strongly suppressed on axis and is redistributed off axis [1105.1725]. The paraxial Maxwell analysis of harmonic emission from a helical undulator gives, for the \(h\)-th harmonic,
\[
\frac{dP_h}{d\Omega} = \frac{\pi}{c} \left[\frac{b_h I(t)\, h N_w K \bar\gamma_z}{((h-1)!(1+K^2))}\right]^2 \left[ \frac{h^2 L_w K^2 \omega_r}{4\pi N_w c(1+K^2)} (\theta_x^2+\theta_y^2) \right]^{h-1}
\]
\[
\times \mathrm{sinc}^2\!\left[ \frac{hL_w\omega_r}{4c}(\theta_x^2+\theta_y^2) \right] \exp\!\left[ -\frac{hL_w\omega_r}{c}(N_x\theta_x^2+N_y\theta_y^2) \right],
\]
with
\[
N_{x,y}=\frac{h\omega_r\sigma_{x,y}^2}{cL_w}
\]
[1105.1725]. For \(h>1\), the factor \((\theta_x^2+\theta_y^2)^{h-1}\) forces the emission to vanish on axis in the ideal helical case.

The experimental benchmark at Elettra confirmed this structure. In coherent harmonic generation mode at \(260\) nm, the measured intensity peaks on axis. In nonlinear harmonic generation mode, with the radiator tuned to \(520\) nm and coherent emission observed at \(260\) nm as the second harmonic, the on-axis signal is much weaker, and the off-axis distribution is consistent with the theoretical prediction that the harmonic emission is concentrated away from the axis [1105.1725]. The comparison of the CHG peak to the residual on-axis NHG signal gave about \(60\) experimentally and about \(140\) in the acceptance-smeared theory, which the authors regarded as satisfactory agreement given alignment uncertainty, beam fluctuations, and model simplifications [1105.1725].

This has direct significance for multifrequency helical devices. A plausible implication is that multifrequency design cannot be assessed solely from harmonic resonance conditions or TAM selection rules; it must also account for angular flux redistribution, finite acceptance, and the rapid reduction of usable off-axis harmonic power. The same study explicitly noted a FERMI@Elettra example in which \(3\times10^{12}\) on-axis photons/pulse at \(5\) nm would correspond to only about \(10^8\) photons/pulse at the second harmonic and far fewer at the third harmonic [1105.1725]. That estimate does not describe a multifrequency device directly, but it quantifies the practical difficulty of using pure helical harmonic radiation as a high-flux operational output.

Related exact work on an extended helical undulator reinforces the importance of angular structure. The exact cylindrical-multipole solution gives the Doppler relation
\[
\omega_m=\frac{m\omega_0}{1-B_z\cos\theta},
\qquad
\cos\theta_m=\frac{1}{B_z}\left(1-\frac{m\omega_0}{\omega}\right),
\]
and identifies the fundamental mode \(m=1\) as the one responsible for radiation exactly on axis, while higher modes contribute near the axis with smaller relative weight [2410.05730]. This exact framework is not a multifrequency theory, but it is directly relevant to controlling harmonic overlap, angular acceptance, and quasi-monochromatic operation in helical devices [2410.05730].

## 5. Hardware realizations and enabling technologies

Recent prototype work has established practical fabrication routes for compact helical undulators built from single-piece rare-earth helices. One approach uses wire electrical discharge machining in combination with a flat tool and rotary workpiece motion to cut helical NdFeB magnets directly from a cylinder [2508.00232]. In the reported \(20\) mm-period prototype, an assembly of two oppositely longitudinally magnetized helices with a half-period shift created a measured on-axis field of \(0.53\) T, which, for a \(20\) mm period, gives \(K\approx0.99\) [2508.00232]. The same paper calculated that Halbach-type helical micro-undulators with periods of \(3\)–\(6\) mm built from four helices could provide a field of \(1\) T and \(K=0.28\)–\(0.6\) on axis [2508.00232].

A later implementation study reported two \(6\) mm-period helical microundulators with a \(1\) mm bore [2509.06186]. The first, composed of two oppositely magnetized longitudinal helices, produced a transverse on-axis field exceeding \(0.93\) T. The second, a hybrid design using two oppositely longitudinally pre-magnetized rare-earth helices alternating with two unmagnetized steel helices, reached about \(1.5\) T on axis [2509.06186]. In that work, the steel helices function as flux shapers approximating an ideal helical Halbach-type arrangement while avoiding the need for radially magnetized rare-earth helices.

The standard undulator parameter relation used in these hardware studies is
\[
K = \frac{e B_u \lambda_u}{2\pi m_e c} \approx 0.934\, B[\mathrm{T}]\, \lambda_u[\mathrm{cm}]
\]
[2508.00232], and in the microundulator study
\[
K = \frac{eBd}{2\pi mc},\qquad
\lambda \approx \frac{d}{2\gamma^2}(1+K^2)
\]
[2509.06186]. These relations connect field strength and period directly to the compact-FEL parameter regime relevant for multifrequency or harmonic applications.

The current limitation is conceptual rather than mechanical. Both hardware papers explicitly emphasize that they do not directly demonstrate multifrequency emission. Their role is enabling and foundational: they show that precisely fabricated NdFeB helices can achieve strong on-axis fields and useful \(K\) values in short-period geometries, which is the kind of platform needed for more advanced helical undulator architectures, including multifrequency variants [2508.00232, 2509.06186].

A distinct enabling direction is the laser-based helical undulator. A Bessel light beam carrying orbital angular momentum has been proposed as an ultra-short-period helical undulator whose effective period can reach the sub-millimeter range through electron dephasing in a superluminal-phase structure [1701.00910]. That proposal is not a multifrequency theory, but it broadens the hardware concept of what may serve as a helical undulator.

## 6. Modeling frameworks, related regimes, and outlook

Three-dimensional and exact modeling frameworks are important because helical multifrequency devices combine polarization control, harmonic structure, transverse mode content, and, in many cases, off-axis emission. A general time-dependent FEL framework capable of modeling planar, helical, and elliptical undulators represents the optical field with Gaussian modes and variable polarization, uses full Newton–Lorentz particle tracking, and includes an analytic APPLE-II model to treat arbitrary elliptical polarizations [1611.01649]. In that treatment, the APPLE-II phase shift \(\phi\) interpolates continuously between planar and helical limits, with
\[
u_e = \frac{1-\cos\phi}{1+\cos\phi}
\]
for \(0\le \phi \le \pi/2\), and
\[
u_e = \frac{1+\cos\phi}{1-\cos\phi}
\]
for \(\pi/2 \le \phi \le \pi\), so that \(\phi=0\) is planar and \(\phi=\pi/2\) is helical [1611.01649]. This framework does not demonstrate explicit multifrequency helical output, but it provides the numerical machinery needed to study polarization-resolved harmonic and sideband dynamics.

At long wavelengths, a THz study based on Liénard–Wiechert fields found that helical undulators typically exhibit reduced side lobes and a more focused main radiation peak, whereas planar undulators are more harmonic-rich but show broader angular structure [2404.16503]. This does not establish a multifrequency helical source in the strict sense, but it clarifies a recurring tradeoff: cleaner directionality and circular polarization in helical geometry versus richer harmonic content in planar geometry.

Twisted-photon work in dispersive media adds another related regime. In undulators filled with a homogeneous dielectric dispersive medium, the helical selection rule
\[
m=n
\]
is preserved, plasma media can suppress lower harmonics, and inverse-polarization domains can produce larger orbital angular momentum at a given harmonic [2004.01443]. These results concern one-frequency helical wigglers in a medium rather than vacuum multifrequency devices, but they show that OAM selection structures remain robust under non-vacuum modifications of the spectrum.

Taken together, these results define the present state of the subject. The mature theoretical core lies in the multifrequency twisted-photon treatments, where helical superpositions give explicit control over emitted energies, TAM spectra, amplitudes, and relative phases through \(\omega_i\), \(K_i\), and \(\chi_i\) [2507.20242, 2602.07942]. The principal practical constraint is that helical harmonic radiation is strongly structured angularly and may be severely flux-limited on axis [1105.1725]. The principal engineering advance is the emergence of compact, high-field helical prototypes fabricated from single-piece rare-earth helices, which provide realistic physical platforms for future multifrequency implementations [2508.00232, 2509.06186]. This suggests that the near-term development of helical multifrequency undulators will depend on combining precise field synthesis, exact or full-3D modeling, and application-specific choices about whether the priority is structured twisted states, polarization control, spectral extension, or high-flux operation.

Source: https://www.emergentmind.com/topics/helical-multifrequency-undulators