Pair beams unlock beyond-terawatt attosecond free-electron laser pulses

Abstract

Free-electron lasers (FELs) generate the brightest coherent X-ray pulses available, enabling atomic-resolution and femtosecond-timescale studies across physics, chemistry, and biology. Realising their full potential at extreme peak powers and attosecond pulse durations critically depends on sustaining coherent gain across the full bunch length. Yet, the quasi-static longitudinal space-charge field in the ultrahigh-current regime imprints a slice-dependent energy detuning that quenches gain growth, so that current schemes typically sustain efficient lasing only across a limited fraction of the bunch. Here we demonstrate that a quasi-neutral electron-positron pair beam cancels this self-field and enables full-bunch high-gain lasing in ultracompressed beams without external compensation. Three-dimensional particle-in-cell simulations in a single-pass, untapered undulator confirm the mechanism across operating regimes: in the soft X-ray regime, the pair beam reaches $1.85\,\mathrm{TW}$ at $345\,\mathrm{as}$ with enhanced odd-harmonic emission and improved spatial coherence, while the electron-only beam fails to saturate; and a high-harmonic pair-cascade configuration yields ${\sim}10\,\mathrm{TW}$ in isolated ${\sim}3.5\,\mathrm{as}$ spikes with coherent amplification extending to photon energies of ${\sim}177\,\mathrm{keV}$. These results establish a new operating regime for ultrahigh-power attosecond light sources and open a direct route to coherent gamma-ray emission ($\geq\!100\,\mathrm{keV}$) currently inaccessible to magnetic-undulator FELs, with broad implications for ultrafast structural, electronic, and nuclear sciences.

Figures (5)

• Figure 1: SASE FEL simulation scheme and diagnostic layout for the $\mathrm{e^-/e^+}$ pair beam.a Schematic of the single-pass undulator stage. b Transverse profile of the Poynting flux amplitude, sampled in the laboratory frame at three stations along the undulator axis corresponding to initialisation, linear gain, and saturation. c Longitudinal profile of the electric field amplitude recorded in the bunch rest frame at the same three stations; the one-to-one correspondence with panel b follows directly from the inverse Lorentz transformation. d Far-field radiation signal amplitude at the downstream detector. Together, panels b--d form a complementary diagnostic set (near-field intensity, rest-frame field structure, and far-field emission) used consistently across all FEL configurations in this work.
• Figure 2: Ellipsoidal versus pancake bunch geometry in the UV FEL: PIC simulation results. All panels compare the quasineutral $\mathrm{e^-/e^+}$ pair beam (orange) with the $\mathrm{e^-}$-only beam (blue); the top row (a--c) shows the ellipsoidal bunch ($L_z = 100\,\mu\mathrm{m}$, $\sigma_\perp = 60\,\mu\mathrm{m}$) and the bottom row (d--f) the pancake bunch ($L_z = 10\,\mu\mathrm{m}$, $\sigma_\perp = 200\,\mu\mathrm{m}$). a, d Longitudinal macro-particle density profiles at saturation, used as a proxy for microbunching extent. b, e FEL power growth curves; insets show the kinetic energy change distribution at saturation, with dashed lines marking the mean energy shift per species, confirming synchronous energy transfer to the radiation field. c, f Far-field power spectra up to the 7th harmonic with odd-harmonic intensity ratios below each spectrum. The legend in c and f applies to all panels in the respective row. g Cross-geometry comparison of the far-field spectrum between the $\mathrm{e^-}$-only ellipsoidal and $\mathrm{e^-/e^+}$ pancake beams, highlighting the harmonic enhancement from DC-field neutralisation.
• Figure 3: Pair-beam SASE in the soft X-ray regime. Fixed parameters: $E_b = 2\,\mathrm{GeV}$, $R_b \approx 90\,\mu\mathrm{m}$, $\varepsilon_n = 0.3\,\mu\mathrm{m\,rad}$, $\sigma_\delta = 0.1\%$, $\lambda_u = 3\,\mathrm{cm}$, $K = 3$ (single, untapered undulator; $L_{\mathrm{sat}} \lesssim 6\,\mathrm{m}$), $\lambda_{\mathrm{r}} \approx 5.4\,\mathrm{nm}$. Two charge--length settings are shown: $(Q_b,\,L_b) = (165\,\mathrm{pC},\,2\,\mu\mathrm{m})$ in panel a and $(330\,\mathrm{pC},\,1\,\mu\mathrm{m})$ in panel b, corresponding to peak currents of ${\sim}\,25\,\mathrm{kA}$ and ${\sim}\,100\,\mathrm{kA}$; all other parameters are held fixed. a, b Temporal power envelope at the undulator exit for the $\mathrm{e^-/e^+}$ pair beam (orange) and the $\mathrm{e^-}$-only beam (blue). c Initial longitudinal momentum phase space in the laboratory frame, confirming beam quality within experimentally demonstrated ranges. d Longitudinal momentum phase space at saturation in the bunch rest frame; the pair beam retains a compact resonant core while the $\mathrm{e^-}$-only beam shows phase-space heating consistent with LSC-driven phase mixing. e, f Transverse particle distribution in the central slice of the pair beam (e) and the $\mathrm{e^-}$-only beam (f) at saturation, with fundamental bunching factors $|b_1| \approx 0.67$ and $|b_1| \approx 0.18$, respectively; the pair beam yields terawatt-class attosecond SASE spikes.
• Figure 4: Coherent $\gamma$-ray lasing with the HHPC-FEL pair-beam concept. Fixed parameters: $E_{\mathrm{slice}} = 48\,\mathrm{GeV}$, $Q_{\mathrm{slice}} \approx 5\,\mathrm{pC}$, $L_{\mathrm{slice}} \approx 4\,\mathrm{nm}$, $R_{\mathrm{slice}} \approx 15\,\mu\mathrm{m}$, $\varepsilon_n = 0.1\,\mu\mathrm{m\,rad}$, $\sigma_\delta = 0.03\%$, $\lambda_u = 3\,\mathrm{cm}$, $K = 2.5$ (single, untapered undulator; $L_{\mathrm{sat}} \approx 18\,\mathrm{m}$), $\lambda_{\mathrm{r}} \approx 7\,\mathrm{pm}$. a, b Rest-frame radiation intensity snapshots at $t_{\mathrm{sim}} = 1.6\,\mathrm{ps}$ and $3.1\,\mathrm{ps}$, capturing sequential gain states of the $\mathrm{e^-/e^+}$ slice comb. At $t = 1.6\,\mathrm{ps}$ the leading slice has saturated while the trailing slice is still in exponential growth, confirming quasi-independent lasing; slice envelopes are traced by the cyan macro-particle sample. The apparent longitudinal compression of the slice train arises because the beam outside the undulator outruns the beam inside, an effect not observable in the laboratory frame. c Zoomed intensity snapshot of the mid-slice at $t_{\mathrm{sim}} = 3.1\,\mathrm{ps}$. The white envelope marks the slippage-broadened FEL pulse ($\mathrm{FWHM} \sim 3.5\,\mathrm{as}$ in the laboratory frame); the $\mathrm{e^-}$-only signal (green) fluctuates incoherently about zero, confirming the absence of gain without DC-field neutralisation. d Far-field power spectrum showing coherent amplification to the 5th harmonic, with the fundamental at ${\sim}\,177\,\mathrm{keV}$ and odd-harmonic intensity ratios. e, f Transverse radiation intensity profiles downstream of the 2nd slice for the $\mathrm{e^-/e^+}$ pair beam (e) and the $\mathrm{e^-}$-only beam (f), illustrating the spatial-coherence contrast at saturation.
• Figure 5: Saturation-power landscape as a function of total charge and bunch length. Colour map: estimated saturation power $P_{\mathrm{sat}}$ (Eq. \ref{['eq:Psat']}) across the $(Q_{\mathrm{beam}},\,L_{\mathrm{beam}})$ plane for fixed parameters $K = 1.5$, $R_b \simeq 250\,\mu\mathrm{m}$, $E_b \simeq 10\,\mathrm{GeV}$. Dashed contours show $\Lambda_{\mathrm{LSC}}$ at three thresholds: negligible suppression ($\Lambda_{\mathrm{LSC}} = 0.2$, black), onset of gain reduction ($\Lambda_{\mathrm{LSC}} = 0.7$, red), and breakdown of high-gain lasing ($\Lambda_{\mathrm{LSC}} = 2.0$, blue). Although $P_{\mathrm{sat}}$ increases with peak current toward the upper-left corner, the contours impose hard boundaries for an electron-only beam: beyond the blue contour, LSC detuning ultimately destroys coherent gain. A quasi-neutral pair beam cancels the DC self-field, driving $\Lambda_{\mathrm{LSC}} \to 0$ regardless of bunch charge or length, and unlocks the full high-power domain across the ultrashort, high-charge operating space.