Interferometry for full temporal reconstruction of laser-plasma accelerator-based seeded free electron lasers

TL;DR

This paper demonstrates how interferometry can be used to fully reconstruct the temporal and spectral properties of laser-plasma accelerator-based seeded free electron laser pulses in a single shot, aiding their development.

Contribution

It introduces a novel interferometric method for complete spectro-temporal characterization of LPA-based seeded FEL pulses using numerical simulations and analytical models.

Findings

Interferometric patterns reveal detailed FEL pulse properties.

Single-shot reconstruction of spectro-temporal content is possible.

Method can handle singular spatio-spectral distributions.

Abstract

The spectacular development of Laser-Plasma Accelerators (LPA) appears very promising for a free electron laser application. The handling of the inherent properties of those LPA beams already allowed controlled production of LPA-based spontaneous undulator radiation. Stepping further, we here unveil that the forthcoming LPA-based seeded FELs will present singular spatio-spectral distributions. Relying on numerical simulations and simple analytical models, we show how those interferometric patterns can be exploited to retrieve, in single-shot, the spectro-temporal content and source point properties of the FEL pulses.

Tables

Table 1. Table 1: Parameters for FEL simulation and modeling.

	Parameter	Value
Electron beam	Energy	200 MeV
	Relative energy spread	4 $\%$
	Charge $Q$	20 pC
	Divergence ${\sigma }_{x0,y0}^{\prime }$	1 mrad-rms
	Emittance ${ϵ}_{x0,y0}$	0.1 mm.mrad
	Chicane dispersion $R_{56,c}$	2.5 mm
Undulator	Period ${\lambda }_u$	18 mm
	Deflection parameter $K$	1.864
Seed laser	Central wavelength	160 nm
	Quadratic phase	1.15$\times {10}^{-4}$ fs-2
	Duration	370 fs-FWHM
	Peak power	10 kW
	Waist position	0 m
	Rayleigh length	0.7 m

Equations

$${\tilde{v}}_z\approx c\left[1-\frac{1}{2{\gamma }^2}\left(1+\frac{K^2}{2}\right)\right]$$

$$dt=z\frac{d}{d\gamma }\left(\frac{1}{{\tilde{v}}_z}\right)d\gamma \approx -\frac{z}{c}\left(1+\frac{K^2}{2}\right)\frac{d\gamma }{{\gamma }^3}$$

$$\Delta \lambda =c\Delta t=-z\left(1+\frac{K^2}{2}\right)\frac{1}{{\gamma }^2}\frac{\Delta \gamma }{\gamma }=R_{56,u}(z)\frac{\Delta \gamma }{\gamma }$$

$$\Delta s=R_{56,c}\frac{\Delta \gamma }{\gamma }$$

$$\begin{array}{rl}\Delta \lambda (z)& =\frac{R_{56,u}(z)}{R_{56,c}}{\lambda }_{seed}=-\frac{z\left(1+\frac{K^2}{2}\right)}{{\gamma }_0^2{R}_{56,c}}{\lambda }_{seed}\end{array}$$

$${\sigma }_{\gamma }\approx \frac{1}{R_{56,c}}\sqrt{{\sigma }_{z0}^2+2R_{522}^2{\sigma }_{x0}^{{}^{\prime }4}+2R_{544}^2{\sigma }_{y0}^{{}^{\prime }4}}$$

$$E_{\text{temp}}(t)=e^{-t^2/(4{\sigma }_t^2)}{e}^{i(\varphi +{\omega }_{0}t+C_1{t}^2)}$$

$$E_{fel}(x,y,z,t)=E_{seed}(x,y,z,t)+E_{coh}(x,y,z+\Delta z,t-\tau )$$

Full text

Free Electron Lasers (FELs) Madey (1971); Deacon et al. (1977) deliver ultrashort, narrow-band and ultrabright pulses down to the hard X-ray range Emma et al. (2010); Ishikawa et al. (2012); Weise and Decking (2017), enabling breakthrough experiments in chemical, physical, and biological sciences. These light sources rely on relativistic electron beams wiggling in the periodic magnetic field of an undulator as gain medium. Interacting with the spontaneous radiation of the undulator or an external seed, the electrons experience an energy modulation at the resonance wavelength which is further transformed into a density modulation by dispersive elements. After this “lethargy” Bonifacio et al. (1984) phase, the beam density modulation allows the emission of a coherent radiation which can then be exponentially amplified. A saturation is reached when the electrons energy loss is such that the resonance condition is violated, causing a red spectral shift of the FEL line Wang et al. (2009); Schneidmiller and Yurkov (2015).

Chicanes as dispersive elements are used to speed up the electron beam energy to density modulation conversion Csonka (1978); Vinokurov and Skrinskii (1978). Associated to electron beam energy chirping in the Radio-Frequency Accelerators (RFAs) structures producing those beams and then to frequency chirped seeds, they opened the door to FEL wavelength tunability Shaftan and Yu (2005), two-color operation De Ninno et al. (2013), amplitude and phase control Wu et al. (2007a); Lutman et al. (2009), spectro-temporal shaping Gauthier et al. (2015), chirped-pulse compression Wu et al. (2007b); Gauthier et al. (2016) and temporal reconstruction De Ninno et al. (2015).

In high-gain FELs Bonifacio et al. (1984), the initial radiation can be amplified by several orders of magnitudes up to saturation in one single pass. The gain of those systems is mainly driven by the electron beam properties. This is the reason why first FELs were developed using classical RFAs which presently deliver the requisite high-quality beams. However since those accelerators are typically several hundreds of meters long for a lasing in the X-ray range, the new Laser Plasma Accelerators (LPAs) Tajima and Dawson (1979) have been proposed as an alternative towards miniaturised FELs Grüner et al. (2007).

However from proposal to demonstration, the path is tortuous. In LPAs, the intrinsic large divergence ($\approx$ few milliradians) and energy spread ($\approx$ few percent) of the electron beam Couperus et al. (2017); Tsai et al. (2018) dramatically limit the FEL gain. Several techniques can be implemented to reach the mandatory Pierce parameter $\rho_{fel}$ Bonifacio et al. (1984) qualifying the FEL gain: use plasma lenses Kuschel et al. (2016); Thaury et al. (2015); van Tilborg et al. (2015); Chiadroni et al. (2018) and/or high gradient quadrupoles in the immediate vicinity of the LPA source to handle the large divergence Eichner et al. (2007); Mihara et al. (2018); Ghaith et al. (2018), use a chicane Couprie et al. (2014); Maier et al. (2012) or a transverse gradient undulator Smith et al. (1979); Huang et al. (2012) to lower the slice energy spread, and/or use a chromatic matching Loulergue et al. (2015) to synchronize the beam maximum density with the radiation pulse propagation inside the undulator. Few years after the first observation of LPA-based undulator radiation Schlenvoigt et al. (2008); Fuchs et al. (2009); Lambert et al. (2012); Anania et al. (2014), significative achievements in terms of electron beam transport and manipulation André et al. (2018) turned a corner, allowing for control of the synchrotron radiation spectral properties Ghaith et al. (2019). LPA based FEL amplification is now within reach at least in the visible or ultra-violet range. To ease this next step, a seeded configuration should be considered, since it enables to speed up the lethargy phase, i.e. requires shorter undulator lengths and relaxes the transport issues. Given the present LPAs performance, the first demonstration is more likely to be obtained in the low-gain Elias et al. (1976) regime as proposed in Liu et al. (2019).

In this letter, we reveal that those forthcoming LPA based seeded FELs can present original spatio-spectral distributions that enable the full reconstruction of the temporal FEL pulse. Indeed, under simple conditions, the emitted radiation exhibits a clear red shift with respect to the input seed wavelength, leading to spectral interference fringes between this additional emitted radiation and the input seed. In the case of LPAs, the electron beam energy chirp, i.e. longitudinal sorting in energy, is not induced by the accelerating structure but by a dispersive transport of the high energy spread beam through a chicane. The resulting particles velocity gradient is so steep that the undulator dispersion drives the energy to density modulation conversion while simultaneously stretching the density modulation period (ahead travelling faster than tail particles). This effect is usually negligible on RFA based FELs, even though observed on a self-seeding configuration Inoue et al. (2019). In the case of LPA beams however, the stretching is such that the final FEL wavelength can be red shifted by more than one seed linewidth. Due to the low-gain configuration, the red-shifted coherent emission is observed together with the seed, leading to strong intereferences in the spectral domain. As in De Ninno et al. (2015), this interferometric pattern can then be used for temporal reconstruction of the FEL pulses.

A typical setup for a LPA based seeded FEL is shown in Fig. 1. In state-of-the-art LPAs, nonlinear plasma waves are driven in a gas by femtosecond pulses of a multi-TW laser system, allowing for the trapping and further acceleration of electron bunches. The properties of the generated electron beams strongly depend on the LPA implementation. We choose a set of parameters (see Table 1) compatible with the most recent achievements of the groups targeting both quality and reliability, mandatory for an FEL application Couperus et al. (2017); Weingartner et al. (2012); Buck et al. (2013); Wang et al. (2016); Liu et al. (2019). With a short period undulator Benabderrahmane et al. (2017), these typical 200 MeV beams can be used to target the vacuum ultra-violet range. We therefore consider a seed at 160 nm (see Table 1) generated in gas Salières et al. (1995); Lambert et al. (2008) using a branch of the main laser Lambert et al. (2008). The electron beam is transported along a dedicated line allowing for a three steps manipulation Loulergue et al. (2015): refocusing at the LPA exit using high gradient quadrupoles Marteau et al. (2017), longitudinal stretching and energy chirping using a chicane and chromatic matching using standard quadrupoles. Locally driving the beam off-axis, the chicane is also used to insert a mirror for the seed injection on-axis. Experimentally, an imaging spectrometer can provide the spatio-spectral distribution of the FEL radiation. Combined with adequate focusing optics, one can observe either the near- or far-field spectral distribution.

Using parameters of Table 1, a 6D electron beam is generated with a flat top profile in energy and Gaussian profiles in other dimensions. The distribution is then transported to the undulator entrance using second order transport formalism Brown (1972). The FEL amplification is simulated with GENESIS Reiche (1999) using this 6D electron beam distribution together with an amplitude and phase temporal description of the seed while the undulator magnetic field is defined by its period $\lambda_{u}$ and deflection parameter $K$ (see Table 1). The output radiation is then analyzed using Fourier optics (see supplemental materials) to predict the observable near- and far-field spatio-spectral distributions. Typical simulation outputs are presented in Fig. 2. To identify the seed contribution, the charge $Q$ is first forced to zero. As expected, Gaussian distributions are obtained in the near- and far-fields (Fig. 2(a,b)). The seed power $P_{seed}$ is then forced to zero to distinguish the spontaneous emission contribution (Fig. 2(c,d)) expected at a resonance wavelength $\lambda_{R}=\lambda_{u}/(2\gamma^{2})(1+K^{2}/2+\theta^{2}\gamma^{2})$ for a mono-energetic beam with a Lorentz factor $\gamma$, $\theta$ being the observation angle. Hence the observed broad band emission results from the large energy spread (different $\gamma$) while the “moon shapes” in the far-field result from the off-axis contribution in $\theta^{2}$ Giannessi et al. (2013). When the seed overlaps the electron beam, i.e. in the FEL case, a significant red shift ($\approx$1.5 nm) appears (Fig. 2(e,f)), together with high contrast fringes.

The standard FEL red spectral shift due to beam energy loss can only be responsible here for sub-nanometer spectral shifts (0.002 nm in the case of Table 1). In order to understand the process underlying here, we reconsider the analytical treatment of the electron beam longitudinal motion along the undulator. For a relativistic electron of energy $\gamma mc^{2}$ travelling in a purely vertical undulator magnetic field, the Lorentz equation in the horizontal direction can be integrated to derive the average longitudinal velocity Huang and Kim (2007):

[TABLE]

The time for an electron to travel the distance $z$ in the undulator is simply $t=z/\tilde{v}_{z}$ so that the difference of travel time $dt$ for electrons of different energies $d\gamma$ close to $\gamma$ is:

[TABLE]

Thus, electrons with different energies induce a stretching $\Delta\lambda$ of the beam due to the longitudinal dispersion $R_{56,u}(z)$ of the undulator:

[TABLE]

In a first approximation, the longitudinal energy sorting of the electrons introduced in the chicane is given by:

[TABLE]

with $R_{56,c}$ the chicane strength or dispersion. Considering now a beam portion of one seed laser period, i.e. $\Delta s=\lambda_{seed}$, and inserting Eq. (4) into Eq. (3), we find that this beam portion is increased by:

[TABLE]

The initial $\lambda_{seed}$ density modulation is stretched proportionnally to the undulator and chicane dispersion ratio up to $\lambda_{seed}+\Delta\lambda(z)$. Since the coherent radiation is emitted at the final density modulation period, the FEL central wavelength is found shifted by $\Delta\lambda$. Because the density modulation is not present immediately at the undulator entrance, i.e. at $z=0$, this model can be further refined replacing $z$ in Eq. (5) by $z-z_{i}$, where $z_{i}$ is the “effective” starting point of the red shifting. Systematic analysis of the simulations reveals that $z_{i}$ could be defined as the location for which the coherent emission intensity becomes $\approx 0.5$ times the initial seed intensity.

Fig. 3(a) shows a typical evolution of the FEL spectrum along the undulator. The FEL wavelength drifts away from the initial seed wavelength, following the analytical prediction of Eq. (5) (solid line). The horizontal dashed line corresponds to the starting point $z_{i}$ in the model. Figure 3(b) summarises the FEL central wavelengths obtained in simulation as a function of the chicane strength for various electron beam charges and divergences. The final wavelengths are found well within the theoretical predictions.

The analysis of the FEL spectral content can then be pushed further. The typical on-axis radiation spectra at the undulator entrance, i.e. of the seed (blue dashed line), and undulator exit, i.e. of the FEL (continuous black line), are shown in Fig. 4(a). In addition to the red shift amplitude, the FEL linewidth $\sigma_{\lambda}$ can be analysed as a function of the chicane dispersion (Fig. 4(b)). It clearly comes out that, for a given charge, the FEL linewidth can be controlled by the chicane dispersion scaling as $1/R_{56,c}$. The FEL linewidth also appears strongly correlated to the initial beam divergence (color trends) while hardly sensitive to the beam emittance (full versus empty dots). In the low gain regime, according to the resonance wavelength expression, a relative energy spread of $\sigma_{\gamma}$ can lead to an inhomogeneous broadening of the relative linewidth by 2$\sigma_{\gamma}$. In the case of LPA beams, using simple analytical formula Labat et al. (2018), the electron beam slice energy spread at the undulator entrance can be approximated by:

[TABLE]

with $R_{522,544}$ the transport matrix coefficicents, $\sigma_{z0}$ the length and $\sigma_{x0,y0}^{\prime}$ the divergence of the beam at the LPA source. Hence, after transport, the initial divergence turns out to be the dominant contribution to the electron beam slice energy spread. The slice energy spread inferred from the electron beam distribution at the undulator exit follows in satisfying agreement with this prediction (Fig. 4(c)). Such a correlation of the FEL linewidth with the initial beam divergence (Fig. 4(d)) can be seen as a special feature of the LPA based FELs.

The on-axis FEL spectra (Figs. 3(a,b)- 4(a)) present fringes that are also visible in the 2D spatio-spectral distribution (Fig. 2(e,f)). Most features of this interferometric pattern can be simply understood modeling the exit radiation as a superposition of two independent Gaussian light pulses: the seed laser ($E_{seed}$) and a coherent contribution ($E_{coh}$). Each light pulse electric field is described in 3D as the product of a transverse field $E_{\text{trans}}(x,y,z)$ and a longitudinal field $E_{\text{temp}}(t)$, where ($x$,$y$) stand for the coordinates in the transverse plane perpendicular to the propagation axis $z$ and $t$ refers to the time (see supplemental material). The longitudinal electric field can be expressed as:

[TABLE]

where $\sigma_{t}$ is the RMS pulse duration, $\phi$ a phase offset, $\omega_{0}$ the optical carrier frequency and $C_{1}$ the quadratic phase. Assuming that the FEL results from the superposition of those two fields, its electric field is finally obtained according to:

[TABLE]

with $\Delta z$ the difference between seed and coherent component waist positions and $\tau$ their relative delay. The ratio $E_{coh}/E_{seed}$ corresponds to the FEL “gain”. Using Fourier optics, calculation of both near- and far-field spectral distributions is then straightforward (see supplemental material). Assuming that the seed properties are known and using the double pulse model (Eq. (8)), the simulated near-field pattern (Fig. 5(a,b)) is fitted using a standard least mean square method. The result of the fit is shown in Fig. 5(c,d). Both near- and far-field distributions present a good agreement with the initial simulations. The few percent fit accuracy is illustrated in Fig. 5(e), comparing the near-field on-axis spectra (i.e. at $y=0$ mm). Among the fitted parameters are found key features such as the FEL source point and Rayleigh length (respectively 0.23 m before the undulator exit and 0.25 m in the case of Fig. 5(c-d)). Finally, the reconstructed FEL intensity and phase in the time-domain are displayed in Fig. 5(f), again in very good agreement with the initial simulation.

In conclusion, we presented simulation results of a feasible LPA based seeded FEL. We showed that it presents original spatio-spectral distributions, with red shifted interference fringes. Those features were explained using simple considerations, allowing a deep understanding of the physical effects underlying. Moreover, we demonstrated that our double pulse model can provide with key features on the FEL source point as well as the full temporal reconstruction of the FEL pulse on a single-shot basis. Such tools should ease the path towards the very first LPA based FEL demonstration. We verified that this work can be extended to FELs in the visible or XUV range. Fringes can be indeed obtained as long as the red shift amplitude is of the order of the seed linewidth, and as the gain remains below two orders of magnitude.