Three Dimensional Theory of the Ion Channel Laser

TL;DR

The paper develops a complete 3D theory of the Ion Channel Laser for the planar off-axis configuration, deriving a paraxial Hamiltonian, pendulum dynamics, Maxwell-Klimontovich equations, and a dispersion relation to capture 3D growth and diffraction effects. It then solves the coupled system with a Van Kampen normal-mode expansion and a Crank-Nicolson numerical scheme to obtain radiation growth rates and transverse profiles, explicitly incorporating diffraction, detuning, and nonzero energy and undulator parameter spreads. The results quantify how diffraction via the Fresnel parameter and beam quality constraints limit the 3D gain, and they provide concrete emittance bounds for several phase-space distributions, informing feasibility and beam-conditioning requirements. These findings guide future work toward 3D PIC simulations and proof-of-principle experiments at facilities such as FACET-II, with implications for compact, high-gain, short-wavelength radiation sources.

Abstract

The ion channel laser (ICL) is a plasma-based alternative to the free electron laser (FEL) that uses the electric field of a uniform-density ion channel rather than the magnetic field of an undulator to induce transverse oscillations of electrons in an ultrarelativistic bunch and thereby produce coherent radiation via a collective electromagnetic instability. The powerful focusing of the ion channel generally yields significantly higher gain parameters in the ICL as compared to the FEL. This permits lasing in extremely short distances using electron bunches with an energy spread as large as a few percent; a value readily achievable with current plasma-based accelerators. ICLs, however, impose stringent transverse phase space requirements on the electron bunch beyond what is required in FELs. In this work, we present a novel 3D theory of the planar off-axis configuration of the ICL that accounts for a number of effects including diffraction, transverse radiation profile, frequency and betatron phase detuning, and nonzero spread in energy and undulator parameter. We derive the ICL pendulum and field equations, which we use to write down the 3D Maxwell-Klimontovich equations. After linearizing, we obtain an integro-differential equation describing the $z$-evolution of the radiation field. The 3D ICL dispersion relation is obtained using a Van Kampen normal mode expansion. We numerically solve the $z$-evolution equation to compute radiation power growth rates and transverse radiation profiles over a range of different ICL parameters. We examine the gain reduction due to 3D effects, energy spread, and emittance. Electron bunch phase space and emittance requirements for lasing are derived. Finally, we make general observations about the performance and feasibility of the ICL and discuss future prospects.

Figures (13)

• Figure 1: Dependence of the Fresnel parameter $\mathcal{F}_D$ defined in Eq. (\ref{['eq:fresnel']}) on the undulator parameter $K_r$, assuming $|\Delta \nu| \ll 1$. $\mathcal{F}_D$ quantifies the level of diffraction in an ICL with smaller values corresponding to a greater degree of diffraction. Note that $\mathcal{F}_D \times (-I / \gamma_r I_A)^{-1/3}$ is plotted instead of $\mathcal{F}_D$ directly to pull out the dependence on $I$ and $\gamma_r$
• Figure 2: Simulated radiation power gain of an ICL with the "base" parameters described in Section \ref{['sec:numericalresults']}.
• Figure 3: Simulated transverse radiation profiles at different $\hat{z}$ positions for an ICL with the "base" parameters described in Section \ref{['sec:numericalresults']}. (a) Normalized two-dimensional field profile (b) Lineouts at $\hat{y}=0$ (solid) and $\hat{x}=0$ (dashed). Figures show the evolution of the radiation profile from the initial seed, through the startup process, to the final steady state.
• Figure 4: Dependence of the relative 3D gain parameter (a) and the transverse radiation sizes in $\hat{x}$ (solid lines) and $\hat{y}$ (dashed lines) (b) on the Cold 1D gain parameter $\rho_0$ for various values of the undulator parameter: $K_r=1$ (blue), $K_r=2$ (green), and $K_r\rightarrow\infty$ (yellow). Results are shown for a cold electron bunch with no frequency detuning ($\Delta \nu = 0$) and are computed using the algorithm described in Section \ref{['sec:numericalalgorithm']}, except for the dashed lines in figure (a) which are computed using the method from Davoine et al.davoine2018.
• Figure 5: Effect of finite Gaussian energy/undulator spread spread on the gain parameter $\rho$ for $K_r=1$ (blue), $K_r=2$ (green), and $K_r\rightarrow\infty$ (yellow), and for $\rho_0=0.01$ (a), $\rho_0=0.02$ (b), and $\rho_0 =0.03$ (c). $\Sigma$ is defined in terms of $\sigma_K$, $\sigma_{\gamma}$, and $\rho_0$ in Eq. (\ref{['eq:bigsigmadef']}).