Coherent Wiggler Radiation

• Coherent Wiggler Radiation (CWR) is the process where a prebunched charged-particle beam traversing a wiggler emits phase-locked radiation, with intensity scaling as the square of the particle number.
• The phenomenon is modeled through formulations such as spontaneous and stimulated superradiance and wake field impedance, detailing resonance conditions and coherence criteria.
• CWR plays a critical role in storage-ring dynamics and novel crystalline undulators, where tuning parameters like bunch length, wiggler period, and magnetic field optimize coherent emission.

Coherent Wiggler Radiation (CWR) denotes the coherent radiation emitted when a short or prebunched charged-particle beam traverses a wiggler or undulator and the emitted fields add with fixed phase rather than randomly. In the superradiant formulation, spontaneous emission from a random beam is proportional to the number of particles, whereas a prebunched beam can emit spontaneously coherent radiation proportional to the number of particles squared; in storage-ring dynamics, the same phenomenon appears as a longitudinal impedance of damping wigglers that can lower the microwave-instability threshold (Gover et al., 2018, He et al., 30 Jul 2025). The subject includes coherent spontaneous emission, stimulated superradiance in seeded systems, wake and impedance formulations for shielded and unshielded wigglers, self-interaction of a bunch with its own coherent field, and crystalline-undulator realizations (Stupakov et al., 2016, MacArthur et al., 2019, Gevorgyan et al., 2024).

1. Physical meaning and coherence criteria

In the review formulation of superradiance, the terms “wiggler” and “undulator” are used interchangeably, and coherent wiggler radiation is identified with coherent spontaneous radiation from a prebunched beam. The basic single-bunch coherence condition is

$2\sigma_{tb}<T=\frac{2\pi}{\omega},$

with $\sigma_{tb}$ the rms bunch duration and $T$ the radiation period. For a Gaussian bunch,

$$M_b(\omega)=\int_{-\infty}^{\infty} e^{i\omega t}f(t)\,dt =e^{-\omega^2\sigma_{tb}^2/2},$$

so coherence is preserved only while the bunching factor $M_b$ remains substantial. For a finite train of equally spaced bunches, the macrobunch form factor is

$$M_M(\omega)=\frac{\sin(N_M\pi\omega/\omega_b)}{N_M\sin(\pi\omega/\omega_b)},$$

which yields narrow comb lines at $\omega=n\omega_b$ (Gover et al., 2018).

A complementary formulation, used for microbunched beams in crystalline undulators, writes the total radiation as a coherent-plus-incoherent sum,

$$N_{\text{tot}}(\omega,\vartheta)=N_{\text{incoh}}(\omega,\vartheta)+N_{\text{coh}}(\omega,\vartheta),$$

with

$$N_{\text{incoh}}(\omega,\vartheta)=N_b(1-F)\,N_{ph}(\omega,\vartheta),\qquad N_{\text{coh}}(\omega,\vartheta)=N_b^2F\,N_{ph}(\omega,\vartheta).$$

Here $F=F_ZF_R$ is the bunch form factor. For a Gaussian bunch,

$\sigma_{tb}$0

At zero angle,

$\sigma_{tb}$1

so an unmodulated bunch radiates coherently only when $\sigma_{tb}$2. For a modulated bunch with density modulation $\sigma_{tb}$3,

$\sigma_{tb}$4

and the coherent term dominates when

$\sigma_{tb}$5

This is the standard prebunched-beam CWR mechanism in explicit form (Gevorgyan et al., 2024).

Short localized structures inside a longer bunch can also satisfy the coherence condition. In a six-period high-$\sigma_{tb}$6 magnetic wiggler, a current spike in the bunch tail that is shorter than the resonant wavelength emits coherently at the fundamental wavelength $\sigma_{tb}$7; the whole bunch need not be shorter than $\sigma_{tb}$8 (MacArthur et al., 2019). In insertion-device treatments of short-bunch coherent radiation, most coherent radiation occurs at wavelengths approximately equal to the bunch length and at small angles, while backward radiation occurs with wavelength about twice the undulator or wiggler period (Mikhailichenko et al., 2011).

2. Resonance, superradiance, and impedance formulations

In the mode-expansion description of coherent undulator or wiggler radiation, resonance is written through the detuning parameter

$\sigma_{tb}$9

with synchronism at

$T$0

In free space, $T$1, giving

$T$2

Finite interaction length produces the usual resonance envelope through $T$3, and in free space the spectral width is $T$4 (Gover et al., 2018).

The same review gives the canonical $T$5 versus $T$6 scaling. For a random beam,

$T$7

whereas for a perfectly short bunch,

$T$8

With finite bunching,

$T$9

and in the presence of a seed field the stimulated-superradiant interference term is

$$M_b(\omega)=\int_{-\infty}^{\infty} e^{i\omega t}f(t)\,dt =e^{-\omega^2\sigma_{tb}^2/2},$$0

This formalism places CWR, spontaneous superradiance, and stimulated superradiance on the same footing (Gover et al., 2018).

A complementary theory treats CWR as coherent synchrotron radiation from the oscillatory reference orbit of a wiggler. For an ultrarelativistic short line bunch on a plane orbit between perfectly conducting parallel plates, the longitudinal impedance is

$$M_b(\omega)=\int_{-\infty}^{\infty} e^{i\omega t}f(t)\,dt =e^{-\omega^2\sigma_{tb}^2/2},$$1

Under the high-frequency and paraxial approximations this becomes

$$M_b(\omega)=\int_{-\infty}^{\infty} e^{i\omega t}f(t)\,dt =e^{-\omega^2\sigma_{tb}^2/2},$$2

For a finite wiggler the impedance is split into $$M_b(\omega)=\int_{-\infty}^{\infty} e^{i\omega t}f(t)\,dt =e^{-\omega^2\sigma_{tb}^2/2},$$3, representing source and observer both inside the wiggler, entrance transient, and exit transient. For an infinitely long free-space wiggler, the low-frequency asymptotics are

$$M_b(\omega)=\int_{-\infty}^{\infty} e^{i\omega t}f(t)\,dt =e^{-\omega^2\sigma_{tb}^2/2},$$4

with $$M_b(\omega)=\int_{-\infty}^{\infty} e^{i\omega t}f(t)\,dt =e^{-\omega^2\sigma_{tb}^2/2},$$5. In this framework, CWR is precisely the CSR generated when the reference orbit is the oscillatory trajectory of a wiggler rather than a single bend (Stupakov et al., 2016).

3. Shielding, chamber geometry, and finite-length structure

Free-space CWR estimates are not generally sufficient in realistic insertion devices. In short-bunch insertion-device radiation, the emitted wavelength for harmonic $$M_b(\omega)=\int_{-\infty}^{\infty} e^{i\omega t}f(t)\,dt =e^{-\omega^2\sigma_{tb}^2/2},$$6 is

$$M_b(\omega)=\int_{-\infty}^{\infty} e^{i\omega t}f(t)\,dt =e^{-\omega^2\sigma_{tb}^2/2},$$7

so backward observation $$M_b(\omega)=\int_{-\infty}^{\infty} e^{i\omega t}f(t)\,dt =e^{-\omega^2\sigma_{tb}^2/2},$$8 gives

$$M_b(\omega)=\int_{-\infty}^{\infty} e^{i\omega t}f(t)\,dt =e^{-\omega^2\sigma_{tb}^2/2},$$9

and for the fundamental $M_b$0. The same treatment states that most coherent radiation occurs at wavelengths approximately equal to the bunch length and at small forward angles, whereas backward radiation can become fully coherent because its wavelength is much longer (Mikhailichenko et al., 2011).

Vacuum chambers modify this picture through cutoff and waveguide dispersion. For a chamber with characteristic transverse size $M_b$1, the guided wavelength satisfies

$M_b$2

The consequence is that the chamber can suppress, redirect, or resonantly enhance coherent wiggler radiation. In the arbitrary-cross-section simulations of insertion-device coherent radiation, long wavelengths may be below cutoff in the central chamber but propagate in side slits; the coherent field is therefore sensitive to exact slit dimensions and to the resonant properties of those slits (Mikhailichenko et al., 2011).

The parallel-plate impedance formulation makes the same point in a different language. Shielding selects odd vertical waveguide modes,

$M_b$3

and the shielded infinite-wiggler impedance shows resonant-like spikes produced by synchronism between the wiggling beam and chamber waveguide modes. In a rectangular chamber the resonance condition is

$M_b$4

with $M_b$5 odd and $M_b$6 even. This is why shielding does not merely attenuate CWR; it structures it spectrally and changes the short-range wake even when the long-range impedance spectrum differs strongly from the free-space case (Stupakov et al., 2016).

4. Self-interaction, superradiant emission, and beam manipulation

CWR is not only a radiative output channel. In a short, high-$M_b$7 planar wiggler, the coherent field emitted by one region of a bunch can slip forward and act back on another region of the same bunch. In the experimentally demonstrated phase-stable self-modulation regime, the fundamental resonant wavelength is

$M_b$8

and the relevant coherence condition is that the tail current spike be shorter than $M_b$9. Radiation generated by that tail slips over the bunch and modulates the core. The measured outcome was a six-period carrier-envelope-phase stable infrared field with gigawatt power and a few MeV, phase-stable modulation in the beam core (MacArthur et al., 2019).

The same work gives a compact paraxial field model,

$$M_M(\omega)=\frac{\sin(N_M\pi\omega/\omega_b)}{N_M\sin(\pi\omega/\omega_b)},$$0

and, after transverse averaging, the mean relative energy modulation in the beam core becomes

$$M_M(\omega)=\frac{\sin(N_M\pi\omega/\omega_b)}{N_M\sin(\pi\omega/\omega_b)},$$1

This explicitly attributes the quasi-single-cycle energy modulation to diffraction and finite wiggler length rather than to a single-cycle emitted pulse (MacArthur et al., 2019).

Benchmarking of OPAL-FEL against LCLS and AWA experiments shows two limiting collective regimes in short wigglers. In the radiation-dominated LCLS case, a short strong wiggler produced a chirp in the bunch center and a single-cycle energy modulation. In the lower-energy AWA case, the dominant effect was space-charge-enhanced longitudinal interaction and energy-spread growth. A practical regime discriminator introduced there is the radiation diffraction limit

$$M_M(\omega)=\frac{\sin(N_M\pi\omega/\omega_b)}{N_M\sin(\pi\omega/\omega_b)},$$2

with radiative effects becoming negligible when $$M_M(\omega)=\frac{\sin(N_M\pi\omega/\omega_b)}{N_M\sin(\pi\omega/\omega_b)},$$3 (Albà et al., 2021). This suggests that CWR in short wigglers should be treated as a continuum between propagating coherent radiation and near-field longitudinal self-interaction, rather than as a purely far-field source term.

5. Storage-ring impedance and microwave instability

In low-emittance storage rings, CWR is treated as a high-frequency longitudinal impedance source associated with damping wigglers. In the SuperKEKB LER study of CSR-driven microwave instability, the impedance model explicitly included the CSR impedance of bending magnets and damping wigglers, the latter identified as coherent wiggler radiation. The CSRZ code was used to calculate the impedances, and Vlasov-Fokker-Planck simulations showed that above the microwave-instability threshold the high-frequency CSR and CWR impedances significantly drive microbunching, leading to an additional increase in bunch length and energy spread. In that case, CSR from regular bends was the most significant source in the high-frequency region, but CWR still contributed notably (Dastan et al., 2024).

A dedicated collider study gives explicit CWR models and threshold scalings. In the low-frequency free-space steady-state regime,

$$M_M(\omega)=\frac{\sin(N_M\pi\omega/\omega_b)}{N_M\sin(\pi\omega/\omega_b)},$$4

with

$$M_M(\omega)=\frac{\sin(N_M\pi\omega/\omega_b)}{N_M\sin(\pi\omega/\omega_b)},$$5

valid for

$$M_M(\omega)=\frac{\sin(N_M\pi\omega/\omega_b)}{N_M\sin(\pi\omega/\omega_b)},$$6

For threshold analysis the paper uses

$$M_M(\omega)=\frac{\sin(N_M\pi\omega/\omega_b)}{N_M\sin(\pi\omega/\omega_b)},$$7

and for low-frequency CWR obtains

$$M_M(\omega)=\frac{\sin(N_M\pi\omega/\omega_b)}{N_M\sin(\pi\omega/\omega_b)},$$8

The same work also develops transient parallel-plate and free-space models, denoted PP-TR and FS-TR, and shows that in the high-frequency regime the free-space impedance recovers the $$M_M(\omega)=\frac{\sin(N_M\pi\omega/\omega_b)}{N_M\sin(\pi\omega/\omega_b)},$$9 scaling familiar from CSR (He et al., 30 Jul 2025).

For the Super Tau-Charm Facility, the 1 GeV case is the most restrictive. The design bunch current is $\omega=n\omega_b$0 mA, the quoted CSR threshold is $\omega=n\omega_b$1 mA, and the quoted CWR threshold is $\omega=n\omega_b$2 mA. At higher energies the CWR threshold rises: at $\omega=n\omega_b$3 GeV it is $\omega=n\omega_b$4 mA, at $\omega=n\omega_b$5 GeV it is $\omega=n\omega_b$6 mA, and at $\omega=n\omega_b$7 GeV it is $\omega=n\omega_b$8 mA (He et al., 30 Jul 2025).

Beam energy	Design bunch current	$\omega=n\omega_b$9
$$N_{\text{tot}}(\omega,\vartheta)=N_{\text{incoh}}(\omega,\vartheta)+N_{\text{coh}}(\omega,\vartheta),$$0 GeV	$$N_{\text{tot}}(\omega,\vartheta)=N_{\text{incoh}}(\omega,\vartheta)+N_{\text{coh}}(\omega,\vartheta),$$1 mA	$$N_{\text{tot}}(\omega,\vartheta)=N_{\text{incoh}}(\omega,\vartheta)+N_{\text{coh}}(\omega,\vartheta),$$2 mA
$$N_{\text{tot}}(\omega,\vartheta)=N_{\text{incoh}}(\omega,\vartheta)+N_{\text{coh}}(\omega,\vartheta),$$3 GeV	$$N_{\text{tot}}(\omega,\vartheta)=N_{\text{incoh}}(\omega,\vartheta)+N_{\text{coh}}(\omega,\vartheta),$$4 mA	$$N_{\text{tot}}(\omega,\vartheta)=N_{\text{incoh}}(\omega,\vartheta)+N_{\text{coh}}(\omega,\vartheta),$$5 mA
$$N_{\text{tot}}(\omega,\vartheta)=N_{\text{incoh}}(\omega,\vartheta)+N_{\text{coh}}(\omega,\vartheta),$$6 GeV	$$N_{\text{tot}}(\omega,\vartheta)=N_{\text{incoh}}(\omega,\vartheta)+N_{\text{coh}}(\omega,\vartheta),$$7 mA	$$N_{\text{tot}}(\omega,\vartheta)=N_{\text{incoh}}(\omega,\vartheta)+N_{\text{coh}}(\omega,\vartheta),$$8 mA
$$N_{\text{tot}}(\omega,\vartheta)=N_{\text{incoh}}(\omega,\vartheta)+N_{\text{coh}}(\omega,\vartheta),$$9 GeV	$$N_{\text{incoh}}(\omega,\vartheta)=N_b(1-F)\,N_{ph}(\omega,\vartheta),\qquad N_{\text{coh}}(\omega,\vartheta)=N_b^2F\,N_{ph}(\omega,\vartheta).$$0 mA	$$N_{\text{incoh}}(\omega,\vartheta)=N_b(1-F)\,N_{ph}(\omega,\vartheta),\qquad N_{\text{coh}}(\omega,\vartheta)=N_b^2F\,N_{ph}(\omega,\vartheta).$$1 mA

The same study identifies shortening the wiggler period $$N_{\text{incoh}}(\omega,\vartheta)=N_b(1-F)\,N_{ph}(\omega,\vartheta),\qquad N_{\text{coh}}(\omega,\vartheta)=N_b^2F\,N_{ph}(\omega,\vartheta).$$2 as a particularly effective mitigation. For the 1 GeV STCF design, tracking with the PP-TR model gives a threshold of $$N_{\text{incoh}}(\omega,\vartheta)=N_b(1-F)\,N_{ph}(\omega,\vartheta),\qquad N_{\text{coh}}(\omega,\vartheta)=N_b^2F\,N_{ph}(\omega,\vartheta).$$3 mA for $$N_{\text{incoh}}(\omega,\vartheta)=N_b(1-F)\,N_{ph}(\omega,\vartheta),\qquad N_{\text{coh}}(\omega,\vartheta)=N_b^2F\,N_{ph}(\omega,\vartheta).$$4 m, $$N_{\text{incoh}}(\omega,\vartheta)=N_b(1-F)\,N_{ph}(\omega,\vartheta),\qquad N_{\text{coh}}(\omega,\vartheta)=N_b^2F\,N_{ph}(\omega,\vartheta).$$5 mA for $$N_{\text{incoh}}(\omega,\vartheta)=N_b(1-F)\,N_{ph}(\omega,\vartheta),\qquad N_{\text{coh}}(\omega,\vartheta)=N_b^2F\,N_{ph}(\omega,\vartheta).$$6 m, and $$N_{\text{incoh}}(\omega,\vartheta)=N_b(1-F)\,N_{ph}(\omega,\vartheta),\qquad N_{\text{coh}}(\omega,\vartheta)=N_b^2F\,N_{ph}(\omega,\vartheta).$$7 mA for $$N_{\text{incoh}}(\omega,\vartheta)=N_b(1-F)\,N_{ph}(\omega,\vartheta),\qquad N_{\text{coh}}(\omega,\vartheta)=N_b^2F\,N_{ph}(\omega,\vartheta).$$8 m (He et al., 30 Jul 2025).

6. Crystalline and structured-mode extensions, and scope of the term

A direct crystalline analogue of CWR is developed for a modulated positron bunch in a crystalline undulator. The average motion is

$$N_{\text{incoh}}(\omega,\vartheta)=N_b(1-F)\,N_{ph}(\omega,\vartheta),\qquad N_{\text{coh}}(\omega,\vartheta)=N_b^2F\,N_{ph}(\omega,\vartheta).$$9

with $F=F_ZF_R$0 and $F=F_ZF_R$1. In this system the coherent spectrum near resonance is

$F=F_ZF_R$2

and the total coherent photon number is

$F=F_ZF_R$3

Because the dielectric permittivity is

$F=F_ZF_R$4

the crystalline medium generates not only the usual hard branch $F=F_ZF_R$5 but also a soft zero-angle branch $F=F_ZF_R$6. For the proposed diamond-undulator example with a modulated positron bunch, the predicted coherent output is $F=F_ZF_R$7 photons at $F=F_ZF_R$8 keV, and coherent radiation dominates if $F=F_ZF_R$9 (Gevorgyan et al., 2024).

Not every wiggler-related radiation paper, however, is a CWR paper in the collective sense. The planar-wiggler twisted-photon theory derives single-particle twisted-mode emission amplitudes, including recoil, and establishes the forward-radiation selection rule

$\sigma_{tb}$00

with $\sigma_{tb}$01 the projection of total angular momentum and $\sigma_{tb}$02 the harmonic number. The paper is explicit that it is not a paper on coherent wiggler radiation in the collective FEL or bunched-beam sense; its relevance is that the single-particle twisted-mode amplitudes provide an input for any future many-electron coherent twisted-photon theory (Bogdanov et al., 2019).

Related generalizations extend CWR concepts beyond magnetic wigglers. In a counterpropagating intense laser pulse, the laser acts as a traveling electromagnetic wiggler, and the bunch spectrum factorizes as

$\sigma_{tb}$03

with the coherent component emitted backward with a narrow cone and an angle-integrated spectrum scaling approximately as $\sigma_{tb}$04 below the coherence cutoff $\sigma_{tb}$05 (Gelfer et al., 2023). By contrast, the laser-driven wire-wiggler $\sigma_{tb}$06-ray source is explicitly described as a self-generated plasma wiggler producing directional synchrotron $\sigma_{tb}$07-rays, but the radiation is not shown to be coherent and should not be treated as evidence for CWR (Wang et al., 2017). Likewise, coherent radiation in photonic crystals involves a periodic electromagnetic structure, bunch form factors, and self-consistent beam-field coupling, but it is not a standard magnetic-wiggler mechanism (Baryshevsky, 2022).

A recurring misconception is therefore to use “CWR” as a generic label for every intense wiggler-like source. The literature distinguishes at least three cases: coherent spontaneous or stimulated emission from a prebunched beam in a magnetic or crystalline wiggler; single-particle structured-mode emission in a wiggler, which can be a building block for coherent theories but is not itself collective CWR; and adjacent plasma, laser-wiggler, or photonic-crystal sources, which may share periodic-radiator physics without satisfying the usual magnetic-wiggler collective-emission definition (Gover et al., 2018, Gevorgyan et al., 2024, Bogdanov et al., 2019).