Optical beat notes assisted attosecond soft X-ray pulse generation in high-gain free electron lasers

Abstract

Attosecond soft X-ray pulses are of great importance for the study of ultrafast electronic phenomena. In this paper, a feasible method is proposed to generate isolated fully coherent attosecond soft X-ray free electron laser via optical frequency beating. Two optical lasers with the opposite frequency chirps are used to induce a gradient frequency energy modulation, which helps to generate a gradually varied spacing electron pulse train. Subsequently, the undulator sections with electron beam delay lines are used to amplify the target ultra-short radiation. Numerical start-to-end simulations have been performed and the results demonstrate that an isolated soft X-ray pulse with the peak power of 330 GW and pulse duration of 620 as can be achieved by the proposed technique.

1 Introduction

The natural timescales of electron motion in molecules and solids systems occur at the several hundred of attoseconds level, and the relevant core-level absorption edges for light elements are found at several hundred eV^{[ 1 ]}. Consequently, the generation of fully coherent soft X-ray laser pulses with the pulse duration shorter than 1 femtosecond (fs) and peak power at the gigawatt (GW) level is of great significance for the study of ultrafast electronic phenomena. High-harmonic generation (HHG)^{[ 2 – 4 ]}, which employs a strong infrared laser to drive electrons in an atomic or molecular gas and generate high-order harmonic ultra-short radiation, provides an important tool for ultrafast science research. However, the relatively low pulse energy (nJ level) and the large output fluctuation limit the application.

The free electron laser (FEL)^{[ 5 ]} is another candidate for the study of ultrafast phenomena because of the highly coherent, high peak power and ultra-short output pulses. In the past two decades, FELs have been witnessed an impressive development worldwide^{[ 6 – 15 ]}. To date, self-amplified spontaneous emission (SASE)^{[ 16 , 17 ]} has been adopted as the major operation mode by most X-ray FEL facilities^{[ 6 – 8 , 11 – 13 ]}, which can generate X-ray pulses with peak power at GW level and pulse duration of approximately 100 fs. Recently, advanced experiments have put new demands on X-ray FELs with shorter pulse durations and higher peak powers. Tens of GW-level attosecond pulses may open a new regime for the study of ultrafast X-ray sciences, such as single-molecule imaging and nonlinear electron dynamics with X-ray spectroscopy methods^{[ 18 , 19 ]}, which are very sensitive to the photon flux. The continuous increase of the pulse energy would significantly increase the number of nonlinear X-ray interactions and extend the realm of ultrafast processes that can be detected.

To meet these requirements, various new techniques either based on SASE^{[ 20 – 25 ]} or seeded FELs^{[ 26 , 27 ]} have been developed. For most of these methods, the output pulse duration and peak power will be limited by the slippage effect. To further shorten the output duration, the mode-locking technique^{[ 28 ]} is proposed to generate attosecond pulse trains. In addition, the pulse-compression scheme^{[ 29 ]} is presented, in which an external laser is employed to generate comb-like current distributions and amplify an isolated pulse via a series of chicane-undulator segments based on the super-radiance process^{[ 30 , 31 ]}. However, the X-ray delay module makes the scheme quite complex. Subsequently, the irregularly spaced current distribution scheme, achieved by the chirped laser^{[ 32 ]}, pulse-stacking or a dedicated modulator^{[ 33 ]} technique, is proposed. In these schemes, irregularly spaced current peaks are generated and the isolated attosecond X-ray pulse is amplified with multi-stage amplification. However, an ultra-short few-cycle laser with a carrier-envelope phase-locking system is required for the generation of the large frequency chirp and stabilizing the FEL output, making this method technically challenging.

In this paper, we propose a simple technique to enhance the peak power of the attosecond radiation pulse while maintaining the ultra-short duration. In the proposed scheme, the optical beating technique is utilized to compress the electron beam^{[ 34 ]}. Different from other methods, two infrared lasers with opposite chirps are adopted in the proposed technique to generate gradually varied spacing current spikes in the electron beam. The target isolated pulse will be amplified in the following chicane-undulator section (radiator) by properly setting the delay lines. In this scheme, a commercially available laser is utilized for the electron beam manipulation, and the frequency beating structure is naturally stable as two infrared lasers are achieved from one laser, which makes the proposed technique feasible and easily implemented in existing X-ray FEL facilities.

The paper is organized as follows. We describe the principle of the proposed scheme in Section 2. The 3D start-to-end simulations based on the Shanghai soft X-ray FEL facility^{[ 14 , 15 ]} with higher beam energy (SXFEL-HE) are carried out and shown in Section 3. Finally, we give the conclusion in Section 4.

2 Principle

In order to compress the electron beam, an energy chirp is usually introduced by either the radiofrequency (RF) or an external laser. After passing through a bunch compressor (BC), the high-energy particles move forward while the low-energy ones move backward, indicating the compression or stretching of the electron beam. For a given strength of the BC, the compression factor depends on the energy chirp of the electron beam, which is determined by the amplitude and frequency of the RF wave or laser. In the FEL facilities, the electron beam is usually compressed to several kA by the BC in the linear particle accelerator (linac). The RF used for bunch compression ranges from 1.3 to 12 GHz^{[ 35 ]}. For advanced electron beam manipulation techniques, the energy chirp is usually induced by an external laser with peak power of approximately 10 GW and pulse energy of approximately 10 mJ in a short undulator (modulator). The required pulse energy of the laser is proportional to the beam length and square of the beam energy^{[ 21 , 22 ]}. The frequency of the laser is about two orders of magnitude higher than that of the RF wave, resulting in a two orders of magnitude larger energy chirp in the electron beam with the same additional energy spread. As a result, the peak current of the electron beam could be enhanced to several tens of kilo-ampere, making the laser-based method much more efficient for generating energy chirp in the electron beam than the RF wave.

In the modulator, the energy modulation is introduced through the laser–electron beam interaction. The energy modulation is converted into density modulation after the beam passes through the magnetic chicane downstream of the modulator. Now we consider that the electron has the initial coordinate $\left(s,p\right)$ , with $s=z-{\overline{v}}_zt$ being the relative coordinate, $z$ the longitudinal coordinate along the electron beam, ${\overline{v}}_z$ the average velocity of the electron beam along the $z$ -axis and $t$ the time. Besides, we define a relative energy $p=\left(\gamma - {\gamma}_0\right)/\gamma$ , where $\gamma$ is the Lorentz factor of the electron and ${\gamma}_0$ denotes the average value for the electron beam. The coordinate of the electron at the end of the magnetic chicane is as follows:

(1) $$\begin{align}{s}^{\prime}=s+{R}_{56}{p}^{\prime },\end{align}$$

(2) $$\begin{align}{p}^{\prime}=p+{p}_0\sin \left(\omega t+{\phi}_0\right),\end{align}$$

where $\omega =2\pi c/\lambda$ is the frequency of the external laser at the wavelength of $\lambda$ , $c$ is the light speed, ${p}_0$ is the energy modulation amplitude in the modulator, ${\phi}_0$ is the initial phase and ${R}_{56}$ is the longitudinal dispersion of the magnetic chicane. In order to maximize the peak current, the optimum dispersion $\rho$ should be as follows:

(3) $$\begin{align}\rho =\frac{\lambda }{2\pi {p}_0},\end{align}$$

where the condition $\mathrm{d}s/ \mathrm{d}p=0$ is satisfied. The longitudinal phase space evolution of the electron beam after passing through the modulator and dispersion section is given in Figure 1 (left).

Figure 1 Longitudinal phase space of the electron beam with normal laser modulation (left), frequency beating with the unchirped laser (middle) and with the chirped laser (right).

Now we consider the frequency beating cases. The electron is transported into the modulator to interact with two unchirped lasers with the central frequencies ${\omega}_1$ , ${\omega}_2$ and energy modulation amplitudes ${p}_1$ , ${p}_2$ , assuming ${\omega}_1>{\omega}_2$ . The coordinates of the electron at the end of the magnetic chicane are as follows:

(4) $$\begin{align}{s}^{\prime}=s+{R}_{56}{p}^{\prime },\end{align}$$

(5) $$\begin{align}{p}^{\prime}=p+{p}_1\sin {\omega}_1t+{p}_2\sin \left({\omega}_2t+\Delta \phi \right),\end{align}$$

where $\Delta \phi$ is the phase difference between two lasers. We take $\Delta \phi =0$ and assume that the two lasers have the same envelope with ${p}_1={p}_2={p}_0$ ; then the energy modulation could be written as follows:

(6) $$\begin{align}{p}^{\prime}=p+2{p}_0\sin \left(\frac{\omega_1+{\omega}_2}{2}t\right)\cos \left(\frac{\omega_1-{\omega}_2}{2}t\right).\end{align}$$

It is obvious that the right-hand side of Equation (6) can be separated as the fast-oscillating term and the slow-oscillating term. The frequency beating node, indicating the slow-oscillating term and equal to the wavelength of the current spike, is dominated by the wave number difference of two lasers. The longitudinal phase space evolution of the electron beam along the modulator and dispersion section is shown in Figure 1 (middle), from which one can see that the current spikes are equispaced. The wavelength of the seed laser is much shorter than the distance between the current spikes.

Furthermore, we consider the case that two linearly-chirped lasers are adopted to interact with the electron beam in the modulator with the central frequencies ${\omega}_1$ and ${\omega}_2$ . An unchirped laser passes through the dispersion component and the phase modulation $\phi \left(\omega \right)$ can be expanded in the Taylor series at ${\omega}_0$ as follows^{[ 36 ]}:

(7) $$\begin{align}\phi \left(\omega \right)=\sum \limits_{m=0}^{\infty}\frac{{\left(\omega -{\omega}_0\right)}^m}{m!}{\phi}_{{m}},\end{align}$$

where the dispersion coefficient ${\phi}_{{m}}$ evaluates the mth derivative of $\phi \left(\omega \right)$ with respect to $\omega$ . We consider the linear chirp and ignore the third and higher order terms of Equation (7). The coordinate of the electron at the end of the magnetic chicane is as follows:

(8) $$\begin{align}{p}^{\prime}=p+{p}_1\sin \left({\omega}_1t+{\alpha}_1{t}^2+{\phi}_1\right)+{p}_2\sin \left({\omega}_2t+{\alpha}_2{t}^2+{\phi}_2\right),\end{align}$$

where the chirp parameter $\left|{\alpha}_{{i}}\right|\cong 1/{\sigma}_0{\sigma}_{{i}}$ , and pulse duration ${\sigma}_{{i}}={\sigma}_0{\left(1+4{\phi}_{{i}}^2/{\sigma}_0^4\right)}^{1/2}$ . Here, ${\sigma}_0$ is the input transform-limited pulse duration. We assume that the two chirped lasers have the same central wavelength ${\omega}_1={\omega}_2={\omega}_0$ , the same envelope ${p}_1={p}_2={p}_0$ and opposite frequency chirps ${\alpha}_1=-{\alpha}_2=\alpha$ . We take ${\phi}_1={\phi}_2=0$ . The energy modulation can be written as follows:

(9) $$\begin{align}{p}^{\prime}&=p+{p}_0\sin \left({\omega}_0t+\alpha {t}^2\right)+{p}_0\sin \left({\omega}_0t-\alpha {t}^2\right)\nonumber\\ &=p+2{p}_0\sin ({\omega}_0t)\cos (\alpha {t}^2).\end{align}$$

The evolution of the longitudinal phase space of the electrons is shown in Figure 1 (right), from which one can see that the gradually varying spacing current enhancement along the electron beam is achieved, which could be used to generate the isolated attosecond X-ray pulses in the following chicane-undulator segments.

The schematic layout of the proposed scheme is shown in Figure 2. In the radiator, the radiation pulse train, which reflects the longitudinal current distribution of the electron beam, is generated in the first undulator segment. The undulator length should be optimized to make the peak power high enough to suppress the shot noise but far from saturation to prevent the degradation of the electron quality. After the first delay line, the target radiation pulse is shifted forward to the upstream current peak and the microbunching (nm level) in the electron beam is washed out. Due to the unequal interval of the current spikes (several hundred nm), other radiation pulses will miss the current peak and the radiation gain processes are suppressed. The ultra-short target pulse is continuously amplified in the downstream undulator segments by repeating the above process.

Figure 2 Schematic layout of the proposed scheme.

3 Simulation

To show the performance of the proposed technique, 3D start-to-end simulations have been carried out based on the parameters of the SXFEL-HE (under design), which consists of a 3-GeV linac and two undulator lines. A high-quality electron beam with the charge of 500 pC is generated in the photoinjector and then compressed to about 1 kA in the linac with two-stage BC chicanes. The longitudinal phase space of the electron beam at the exit of the linac, simulated by the tracking code ASTRA^{[ 37 ]} in the injector and ELEGANT^{[ 38 ]} in the linac, is shown in Figure 3 (left). The bunch length of the electron beam is about 400 fs by the end of the linac. The normalized slice emittance is about 1.0 μm⋅rad and the slice energy spread is about 50 keV. These simulation results fit quite well with the measurements at the SXFEL facility^{[ 39 , 40 ]}.

Figure 3 Longitudinal phase space of the electron beam at the end of the linac (left) and at the exit of the ESASE section (right). The bunch head is to the right.

In the modulation section, a transform-limited laser pulse (longitudinal Gaussian distribution) with the central wavelength of 800 nm and pulse duration of 100 fs (full width at half maximum (FWHM)) is split and sent into two branches. The separated two optical pulses are stretched by a dispersive medium with opposite dispersions, resulting in two 1.1 ps (FWHM) long laser pulses with opposite frequency chirps. Finally, these two chirped laser pulses are sent through a splitter to recombine to a single laser pulse, as shown in the left-hand part of Figure 2, and sent into the modulator to interact with the electron beam to generate the energy modulation. The peak power of each chirped laser pulse is about 15 GW with laser waist of about 300 μm and pulse energy of about 16 mJ. The energy modulation is converted into density modulation in the downstream dispersion section, generating the gradually varied spacing pulse train. The longitudinal phase space of the electron beam at the exit of the modulation section, simulated by the 3D algorithm^{[ 41 ]}, is shown in Figure 3 (right). One can find that the peak current of about 5 kA is achieved with the current spacing of around 40 fs. The head and the tail parts of the electron beam, containing 5.1% particles in total, are cut during the simulation, which has little contribution to the radiation process but will result in the distortion of the electron beam in the longitudinal direction. Subsequently, the electron pulse train is transported to the downstream radiator.

It is worth pointing out that the longitudinal space charge effect is not negligible due to the high peak current^{[ 42 , 43 ]}. The longitudinal space charge field, according to Ref. [43], could be calculated by the following:

(10) $$\begin{align}Ez\approx -\frac{Z_0{I}^{\prime }(s)}{4\pi {\overline{\gamma}}_z^2}\left(2\ln \frac{{\overline{\gamma}}_z{\sigma}_z}{r_{\mathrm{b}}}+1-\frac{r^2}{r_{\mathrm{b}}^2}\right),\end{align}$$

where ${Z}_0=377\ \Omega$ is the free space impedance, ${I}^{\prime }(s)= \mathrm{d}I/ \mathrm{d}s$ is the derivative of the electron current profile with respect to the longitudinal bunch coordinate $s$ , ${\overline{\gamma}}_z=\gamma /\sqrt{1+{K}^2/2}$ , $K$ is the undulator parameter, ${\sigma}_z$ is the root mean square (rms) bunch length of the current spike, ${r}_{\mathrm{b}}$ is the beam radius of a uniform transverse distribution and $r=\sqrt{x^2+{y}^2}$ . Here, we take $K=2.67$ , $\gamma =5871$ , ${r}_{\mathrm{b}}=100$ μm and ${\sigma}_z=293$ nm. With these parameters, the accumulated energy modulation $\Delta E$ after the 39 m undulator is as given in Figure 4. One can see that the longitudinal space charge produces additional energy chirps in the current spikes with a peak-to-peak energy variation of about 20 MeV. The longitudinal space charge effect can significantly degrade the FEL performance and is taken into account in the simulations based on GENESIS^{[ 44 ]}.

Figure 4 Electron beam energy modulation from the longitudinal space charge field from the 39 m undulator.

Thirteen 3-m-long undulators in total are used with the period length of 30 mm. The numbers of undulators used in the first to ninth undulator segments are chosen as 3, 2, 2, 1, 1, 1, 1, 1 and 1, respectively. The simulation results are summarized in Figure 5. In the first undulator segment, an attosecond pulse train, which reflects the comb-like current distribution, is generated with the peak power of about several tens of megawatt. Then the attosecond pulse train is shifted forward in the downstream delay line, so that the target radiation pulse (the last one) can overlap the adjacent electron pulse while the others meet the low-current part of the electron beam because of the unequal interval of the current distribution. Meanwhile, the microbunching accumulated in the first undulator segment is smeared out when the beam passes through the delay line. An X-band deflecting cavity is located at the end of the undulator line with temporal resolution better than 10 fs, which helps one to choose the strength of the delay line during future experimental research. The relative temporal jitter between the laser pulses and the electron beam is only tens of femtoseconds at FEL facilities, for example, approximately 50 fs rms at SXFEL^{[ 39 ]}. A stable comb-like structure without obvious shot-to-shot variations had been experimentally demonstrated at FERMI^{[ 45 ]}, indicating that convenient optimization methods will be feasible for the tuning of the proposed scheme. Besides, one could increase the pulse duration of the seed laser and the electron beam to further increase the tolerance of temporal jitters.

Figure 5 The output radiation evolution and the spectrum along the undulator line at the end of the first, fourth and ninth undulator segments.

In the second undulator segment, the target radiation pulse is shifted and interacts with a ‘fresh’ part of the electron beam, leading to continuous amplification, similar to the direct seeding scheme^{[ 46 ]}. The amplification of the satellite pulses is suppressed, leading to continuous growth of the contrast, which is defined as the percentage of the target radiation pulse energy. By repeating the process of the undulator amplification and radiation delay, an isolated radiation pulse is achieved with a peak power of 330 GW and a pulse duration of about 620 as after nine undulator segments. The contrast of the target radiation pulse is over 99%. The bandwidth of the final output is about 0.5%, which is about 1.1 times that of the Fourier transform limit. The final peak power of the radiation can be further improved by optimizing the parameters of the frequency beating to generate more current spikes and employ more undulator segments. The output pulse duration can be further shortened by using shorter undulators, which helps to decrease the slippage effect.

3D simulations for SASE have also been performed, sharing the same electron beam at the end of the linac and the same undulator parameters. The comparison results are given in Figure 6. It is easy to find that the normal SASE saturates at around 35 m with the peak power of about 20 GW, which is one order of magnitude lower than that of the proposed scheme.

Figure 6 FEL gain curves for the proposed scheme and the normal SASE.

4 Conclusion

In this paper, we proposed a novel and feasible method to generate isolated attosecond radiation pulses based on the electron beam compression with optical beating technique. In the proposed method, two chirped optical lasers with the same central wavelength and opposite frequency chirps are employed to modulate the electron beam and generate the gradually varied spacing current spikes. The isolated attosecond pulse could be generated in the downstream radiator with a multi-stage amplification setup. 3D start-to-end simulations have been performed and the results show that an isolated radiation pulse at 2 nm can be achieved with the peak power of about 330 GW and pulse duration of about 620 as, with the longitudinal space charge effect and longitudinal resistive wall wakefield being taken into account. The output peak power can be further improved and the pulse duration can be shortened by optimizing the parameters of the optical laser, the electron beam and the undulators. For example, by increasing the beam energy (4 GeV) to generate shorter wavelength (1 nm) radiation and increasing durations of the laser pulses and the electron beam to get more current spikes, it is possible to further suppress the slippage effect to get shorter pulses with higher peak power. The proposed scheme is feasible and can be implemented in the existing FEL facilities. This kind of attosecond coherent X-ray light source has great potential applications in the study of ultrafast electronic phenomena.