Capillary Discharge Waveguides

• Capillary discharge waveguides are plasma-filled channels in dielectric capillaries, formed by a pulsed electrical discharge that produces a parabolic electron density profile for low-loss optical guiding.
• They enable precise matching of optical modes, with typical spot sizes of 30–200 µm, facilitating the guidance of high-power lasers and the focusing of relativistic particle beams over multi-centimeter lengths.
• Advanced designs employ laser-heater techniques and engineered geometries to enhance plasma channel control, supporting high-repetition operation (up to 400 Hz) and stable beam transport.

A capillary discharge waveguide is a transient, plasma-filled waveguide formed within a dielectric or ceramic capillary by a pulsed electrical discharge in a background gas. These structures produce long, low-density plasma channels with tailored radial electron density profiles that enable optical guiding of high-power, ultrashort laser pulses or focusing of relativistic charged-particle beams. Capillary discharge waveguides play a central role in laser-wakefield acceleration (LWFA), plasma wakefield acceleration (PWFA), active plasma lenses, and advanced beam transport applications. The discharge-driven plasma channel is typically parabolic near the axis, facilitating low-loss guiding of fundamental electromagnetic or particle modes over multi-centimeter to decimeter lengths, and can be manipulated with pulsed laser heaters or via capillary geometry for advanced functionality.

1. Principles of Plasma Channel Formation in Capillary Discharges

A capillary discharge waveguide is established by filling a dielectric or ceramic capillary (e.g., sapphire, alumina, Shapal, or Macor) with a low-pressure gas (typically H₂, Ar), followed by a high-voltage, high-current electrical discharge pulse. The current pulse (hundreds of amps to several kA, with rise times ≈100–500 ns) Ohmically heats the gas, leading to full or partial ionization and rapid formation of a plasma column along the capillary axis (Bagdasarov et al., 2017, Turner et al., 2020, Crincoli et al., 24 Jul 2025).

After breakdown, the plasma reaches a quasi-steady-state in which Ohmic heating is balanced by thermal conduction losses to the capillary wall. Under uniform pressure and full ionization, the radial power balance is

$\nabla\cdot(\kappa \nabla T_e) + j^2/\sigma = 0,$

where $\kappa$ is the electron thermal conductivity, $j$ the current density, and $\sigma$ the Spitzer (or modified) conductivity. The resulting temperature profile is quadratic:

$T_e(r) = T_e(0) - \Delta T (r/R)^2,$

with $R$ the capillary radius. For $p = n_e k_B T_e(r) = \text{const}$, the electron density is

$$n_e(r) = n_{e,\rm axis} + \Delta n_0 \left(\frac{r}{R}\right)^2,$$

with

$$\Delta n_0 \simeq \frac{p}{k_B T_e(0)} \frac{j^2 R^2}{4 \kappa T_e(0)}.$$

Typical on-axis densities are $n_{e,\rm axis} \approx 10^{17}\ \mathrm{cm}^{-3}$ for $\kappa$0 and $\kappa$1 A (Bobrova et al., 2013, Mewes et al., 19 Jun 2025).

A parabolic density minimum on-axis yields a refractive-index well for laser guiding, described to leading order by $\kappa$2. In the vicinity of the axis, this index profile closely matches the requirements for low-loss optical guiding of high-intensity laser pulses.

2. Optical Guiding and Matched Mode Properties

The plasma channel acts as a graded-index medium for electromagnetic waves:

$\kappa$3

where $\kappa$4 is the critical density for laser frequency $\kappa$5.

For a parabolic $\kappa$6 profile, the matched spot size $\kappa$7 for a lowest-order Gaussian mode is

$\kappa$8

with $\kappa$9 the on-axis density and $j$0 the channel depth (Turner et al., 2020). For typical capillaries ($j$1m, $j$2–$j$3, $j$4), matched mode radii are $j$5 30–200 $j$6m, enabling guiding of relativistic laser pulses over tens of centimeters (Bagdasarov et al., 2017, Crincoli et al., 24 Jul 2025).

In long ($j$740 cm) capillaries, the electron density and matched spot size can be reproduced shot-to-shot to better than 1% and 0.2% variations, respectively (Turner et al., 2020). Laser beam propagation in channels with higher-order ($j$8) deviations from parabolicity maintains diffraction-limited performance when the injected beam is approximately matched ($j$9), but non-parabolic terms produce measurable beam distortions and ring-like modes in "plasma telescope" configurations ($\sigma$0).

The scaling of waveguide properties with capillary and discharge parameters approximately follows (Bagdasarov et al., 2017):

• $\sigma$1
• $\sigma$2

3. Extension and Manipulation of Discharge Channels

a) Advanced Thermal Control: Laser-Heater Techniques

A secondary, ns-scale laser pulse ("laser heater") can precisely deepen and narrow the on-axis plasma channel by selective electron heating, transiently doubling the index curvature and reducing $\sigma$3 by a factor $\sigma$4 for several ns (Bobrova et al., 2013). The governing equations for laser-induced heating are:

$\sigma$5

with $\sigma$6 the inverse bremsstrahlung absorption. The resulting additional curvature ($\sigma$7) scales as $\sigma$8 (heater pulse energy/spot squared) and supports smaller matched spots ($\sigma$9–90 $T_e(r) = T_e(0) - \Delta T (r/R)^2,$0m for $T_e(r) = T_e(0) - \Delta T (r/R)^2,$1 nm laser pulses).

b) Shaping and Focusing of Particle Beams

The same azimuthal discharge current $T_e(r) = T_e(0) - \Delta T (r/R)^2,$2 that sustains the plasma channel establishes a linear $T_e(r) = T_e(0) - \Delta T (r/R)^2,$3 profile:

$T_e(r) = T_e(0) - \Delta T (r/R)^2,$4

This magnetic field can focus relativistic electron beams (active plasma lens). Engineering rectangular or square capillary cross-sections allows tuning of focusing asymmetry, favoring the production and simultaneous focusing of flat beams with high emittance ratios ($T_e(r) = T_e(0) - \Delta T (r/R)^2,$5), as demonstrated in 2D MHD simulations and analysis (Bagdasarov et al., 2017, Bagdasarov et al., 2017).

Curved capillary discharges ("active bending plasma") can confine and guide relativistic bunches along curved trajectories, with minimal distortion of the longitudinal phase space compared to conventional dipoles (Pompili et al., 2017). The bending angle is

$T_e(r) = T_e(0) - \Delta T (r/R)^2,$6

and preservation of bunch length is substantially improved in the capillary channel.

4. Capillary Materials, High Repetition Rate Operation, and Longevity

Material choice governs the waveguide's thermal resilience, machinability, and lifetime under high-fluence, high repetition rate operation. Shapal Hi M Soft™ (AlN-BN hybrid) and Macor® (glass-ceramic) have been implemented in 1 mm $T_e(r) = T_e(0) - \Delta T (r/R)^2,$7 2 mm, $T_e(r) = T_e(0) - \Delta T (r/R)^2,$8–10 cm capillaries for 10–400 Hz operation (Crincoli et al., 24 Jul 2025). Shapal exhibits high thermal conductivity ($T_e(r) = T_e(0) - \Delta T (r/R)^2,$9 92 W/m·K), high melting point ($R$01900°C), and machinability compatible with complex geometries.

Quantitative studies show that with $R$1–$R$2 J and $R$3–$R$4 A, the wall temperature remains within safe limits up to $R$5 Hz, provided heat is efficiently transported to water-cooled holders. Plasma properties ($R$6, $R$7) are maintained to within 5% over $R$8 pulses, with negligible capillary erosion and stable guiding characteristics (Crincoli et al., 24 Jul 2025).

Designs targeting next-generation facilities (e.g., EuPRAXIA) specify operation at 100–400 Hz, longevity $R$9 shots, and stable $p = n_e k_B T_e(r) = \text{const}$0 cm$p = n_e k_B T_e(r) = \text{const}$1 with $p = n_e k_B T_e(r) = \text{const}$2–50 µm.

5. Numerical Modeling, Diagnostics, and Cross-Sectional Effects

Comprehensive modeling of capillary discharges employs single- and multi-fluid, two- or three-temperature MHD, or hydrodynamic codes (e.g., MARPLE, FLASH, NPINCH, INF&RNO, COMSOL Multiphysics) (Bagdasarov et al., 2017, Diaw et al., 2022, Mewes et al., 19 Jun 2025). The choice of electron transport model (Spitzer, Epperlein-Haines, Davies–Wen–Ji–Held) introduces 10–40% variations in predicted electron temperature and channel profile, directly impacting guided laser mode size and focusing strength (Diaw et al., 2022). The Davies–Wen–Ji–Held closure is preferred for accuracy and smooth parameter dependence, with benchmarking against Stark-broadened $p = n_e k_B T_e(r) = \text{const}$3 and (ideally) interferometric data for $p = n_e k_B T_e(r) = \text{const}$4 (Mewes et al., 19 Jun 2025).

Square- and rectangular-section capillaries yield plasma and $p = n_e k_B T_e(r) = \text{const}$5-field profiles on axis that match circular capillaries to within a few percent for $p = n_e k_B T_e(r) = \text{const}$6, making them suitable for laser guiding and enabling asymmetric plasma lensing (Bagdasarov et al., 2017, Bagdasarov et al., 2017). For guiding and focusing, azimuthal inhomogeneity is negligible for laser spots $p = n_e k_B T_e(r) = \text{const}$7, allowing the use of square capillaries to enhance optical diagnostic access and shape electron beams.

6. Dielectric Wakefields and Coupling to Relativistic Bunches

Dielectric-walled capillaries (fused silica, UV-grade polymers) not only support plasma formation but also act as slow-wave structures in the absence of plasma, sustaining longitudinal and transverse (dipole) wakefields when traversed by an off-axis relativistic bunch (Verra et al., 2024). The transverse wake potential induced by misalignment scales as $p = n_e k_B T_e(r) = \text{const}$8 (with $p = n_e k_B T_e(r) = \text{const}$9 misalignment and $$n_e(r) = n_{e,\rm axis} + \Delta n_0 \left(\frac{r}{R}\right)^2,$$0 the aperture), leading to slice-dependent deflection ($$n_e(r) = n_{e,\rm axis} + \Delta n_0 \left(\frac{r}{R}\right)^2,$$1).

Plasma formation screens dielectric wakefields; once the plasma density is sufficiently high ($$n_e(r) = n_{e,\rm axis} + \Delta n_0 \left(\frac{r}{R}\right)^2,$$2 beam-wall distance), the net effect vanishes, and only the plasma's focusing/wakefield properties dominate (Verra et al., 2024). For high-quality beam transport in PWFA, specifications on alignment ($$n_e(r) = n_{e,\rm axis} + \Delta n_0 \left(\frac{r}{R}\right)^2,$$3), aperture, and plasma density are explicit.

7. Mode Structure, Coupling Efficiency, and Waveguide Loss

In plasma channels and metallic capillaries with $$n_e(r) = n_{e,\rm axis} + \Delta n_0 \left(\frac{r}{R}\right)^2,$$4, the fundamental guided mode is an $$n_e(r) = n_{e,\rm axis} + \Delta n_0 \left(\frac{r}{R}\right)^2,$$5 circularly polarized "R" mode, with an attenuation length $$n_e(r) = n_{e,\rm axis} + \Delta n_0 \left(\frac{r}{R}\right)^2,$$6 (Tuev et al., 2020). Optimal coupling of a linearly polarized Gaussian beam into the capillary is achieved for $$n_e(r) = n_{e,\rm axis} + \Delta n_0 \left(\frac{r}{R}\right)^2,$$7 (spot size), resulting in $$n_e(r) = n_{e,\rm axis} + \Delta n_0 \left(\frac{r}{R}\right)^2,$$8 power transfer into the fundamental. Modal attenuation over tens of centimeters enables practical low-loss guiding for LWFA and related applications.

Summary Table: Key Properties and Typical Performance

Parameter	Typical Range	Comment
Capillary radius ($$n_e(r) = n_{e,\rm axis} + \Delta n_0 \left(\frac{r}{R}\right)^2,$$9)	100–2000 $$\Delta n_0 \simeq \frac{p}{k_B T_e(0)} \frac{j^2 R^2}{4 \kappa T_e(0)}.$$0m	Sapphire, AlN, Macor, etc.
Length ($$\Delta n_0 \simeq \frac{p}{k_B T_e(0)} \frac{j^2 R^2}{4 \kappa T_e(0)}.$$1)	1–40 cm	$$\Delta n_0 \simeq \frac{p}{k_B T_e(0)} \frac{j^2 R^2}{4 \kappa T_e(0)}.$$290 mm at BELLA
Discharge current ($$\Delta n_0 \simeq \frac{p}{k_B T_e(0)} \frac{j^2 R^2}{4 \kappa T_e(0)}.$$3)	200–2000 A	340–500 ns rise, $$\Delta n_0 \simeq \frac{p}{k_B T_e(0)} \frac{j^2 R^2}{4 \kappa T_e(0)}.$$4s FWHM
On-axis density ($$\Delta n_0 \simeq \frac{p}{k_B T_e(0)} \frac{j^2 R^2}{4 \kappa T_e(0)}.$$5)	$$\Delta n_0 \simeq \frac{p}{k_B T_e(0)} \frac{j^2 R^2}{4 \kappa T_e(0)}.$$6–$$\Delta n_0 \simeq \frac{p}{k_B T_e(0)} \frac{j^2 R^2}{4 \kappa T_e(0)}.$$7 cm$$\Delta n_0 \simeq \frac{p}{k_B T_e(0)} \frac{j^2 R^2}{4 \kappa T_e(0)}.$$8	Tunable with gas fill
Channel depth ($$\Delta n_0 \simeq \frac{p}{k_B T_e(0)} \frac{j^2 R^2}{4 \kappa T_e(0)}.$$9)	0.05–0.2	Parabolic/weakly-Gaussian
Matched spot size ($n_{e,\rm axis} \approx 10^{17}\ \mathrm{cm}^{-3}$0)	30–200 $n_{e,\rm axis} \approx 10^{17}\ \mathrm{cm}^{-3}$1m	$n_{e,\rm axis} \approx 10^{17}\ \mathrm{cm}^{-3}$2 nm
High-rep operation ($n_{e,\rm axis} \approx 10^{17}\ \mathrm{cm}^{-3}$3)	10–400 Hz	20M shots, $n_{e,\rm axis} \approx 10^{17}\ \mathrm{cm}^{-3}$4C
Guiding length	>10 × Rayleigh range	$n_{e,\rm axis} \approx 10^{17}\ \mathrm{cm}^{-3}$5

Capillary discharge waveguides are thus established as robust, highly configurable plasma channels for advanced laser–plasma acceleration and high-gradient beam manipulation, offering stability, longevity, and fine-tuned optical and focusing properties over the multi-centimeter to decimeter scale (Bobrova et al., 2013, Bagdasarov et al., 2017, Turner et al., 2020, Diaw et al., 2022, Crincoli et al., 24 Jul 2025, Bagdasarov et al., 2017, Pompili et al., 2017, Verra et al., 2024, Mewes et al., 19 Jun 2025, Tuev et al., 2020, Bagdasarov et al., 2017).