Universal route to optimal few- to single-cycle pulse generation in hollow-core fiber compressors.

Gas-filled hollow-core fiber (HCF) pulse post-compressors generating few- to single-cycle pulses are a key enabling tool for attosecond science and ultrafast spectroscopy. Achieving optimum performance in this regime can be extremely challenging due to the ultra-broad bandwidth of the pulses and the need of an adequate temporal diagnostic. These difficulties have hindered the full exploitation of HCF post-compressors, namely the generation of stable and high-quality near-Fourier-transform-limited pulses. Here we show that, independently of conditions such as the type of gas or the laser system used, there is a universal route to obtain the shortest stable output pulse down to the single-cycle regime. Numerical simulations and experimental measurements performed with the dispersion-scan technique reveal that, in quite general conditions, post-compressed pulses exhibit a residual third-order dispersion intrinsic to optimum nonlinear propagation within the fiber, in agreement with measurements independently performed in several laboratories around the world. The understanding of this effect and its adequate correction, e.g. using simple transparent optical media, enables achieving high-quality post-compressed pulses with only minor changes in existing setups. These optimized sources have impact in many fields of science and technology and should enable new and exciting applications in the few- to single-cycle pulse regime.

Introduction

The chirped pulse amplification (CPA) technique applied to Titanium Sapphire lasers has made intense near-infrared (NIR) ultrashort pulses in the 20−100 fs range widely available for scientific, biomedical and industrial applications. Special efforts have been devoted to generate even shorter pulses, in the few- and single-cycle regime, due to a number of interesting applications. In particular, such pulses have paved the way for attosecond physics and metrology[1-6], via the extreme ultraviolet (XUV) attosecond pulse trains and isolated attosecond pulses[7] that can be obtained by high-harmonic generation (HHG). The use of few-cycle optical pulses with durations close to or shorter than 10 fs in the near-infrared, visible and near-ultraviolet spectral regions has been extended in recent years to a wide range of spectroscopic techniques such as impulsive vibrational spectroscopy[8-10], time-resolved stimulated Raman spectroscopy[11-14], and ultrafast pump-probe absorption spectroscopy[15,16]. Few-cycle optical pulses have also become an interesting tool for transient absorption microscopy[17], near-field imaging techniques[18] and for generating ultrashort terahertz radiation[19].

While it is possible to obtain sub-10 fs pulses from CPA[20] or from optical parametric amplification[21] systems, the former is not easy to accomplish, and the latter is not commonplace. Therefore, post-compression techniques are usually employed for the generation of intense few- and even single-cycle pulses in the near- and mid-infrared spectral regions[22,23]. In order to post-compress ultrafast pulses down to the few-cycle regime, two steps are usually needed. First, nonlinear processes broaden the pulse spectrum, thus decreasing the Fourier-limited pulse duration. In a second step, the spectral phase resulting from the previous stage is compensated, typically using chirped mirrors, gratings, prisms, or other dispersive systems, resulting in a temporally compressed pulse. This scheme was first proposed in the context of optical fibers in the 1980s[24], and enabled achieving 6 fs pulses when compensating simultaneously the outcoming group delay dispersion (GDD) and third-order dispersion (TOD)[25]. The scheme was successfully expanded in 1996 by Nisoli and coworkers to the ultra-intense laser pulse regime (mJ-level pulses) by using hollow-core fibers (HCF) filled with gases[26]. Using the latter scheme, together with chirped mirrors as the compression system, few-cycle pulses in the hundreds of μJ energy range with 0.1 TW peak power were obtained[27]. Later, similar results were obtained using the light filamentation process in the spectral broadening stage[28].

In spite of requiring a finer input beam alignment than filamentation-based compressors, HCF compressors are today the most widely used high-energy post-compressed sources, in part due to their intrinsic spatial filtering properties which result in very high quality beam profiles. Furthermore, their long-term stability can be greatly improved by using, e.g., piezo-driven mirror mounts to ensure stable and constant spatial coupling of the input laser pulses into the HCF via a simple feedback loop. Hollow-core fiber post-compressed pulses have shown a great potential in a wide range of applications, such as pump-probe experiments in conjunction with attosecond pulses[29], ultrafast measurement of electrical and optical properties of solids[30,31], time-resolved studies of Coulomb explosion dynamics[32], ultrafast spectroscopy techniques[33-39], and very recently a new generation of compact kHz laser-plasma accelerators based on single-cycle pulses[40].

To access the Fourier Limit of a pulse after a nonlinear propagation process one has to deal with the complex phase that the pulse acquires due to the interplay between different linear and nonlinear effects. In general, researchers optimize their HCF compressors by empirically adjusting several key parameters, such as gas type and pressure, input pulse characteristics and coupling conditions, with the final result usually involving a delicate compromise between output efficiency, amount of spectral broadening and achievable degree of compression (pulse duration and quality) for their particular system and chirped mirror set. The ability to measure and quantify the achieved degree of compression is therefore paramount to identify the main characteristics of the output pulse and to further optimize its compression. Both spectral phase oscillations and the overall spectral phase of a pulse can be visualized in a very straightforward way using the dispersion-scan (d-scan)[41] technique, which has been extensively used in the last years to characterize many state-of-the-art few-cycle pulse sources around the world[6,38-40,42].

D-scan is a recent approach for the simultaneous measurement and compression of femtosecond laser pulses. Its experimental setup is fully in-line, without any beam splitting, recombination or temporal delaying of the pulses. Experimentally it involves the measurement of the spectrum of a nonlinear signal, such as second-harmonic generation (SHG), as a function of dispersion applied to the pulse. This can be performed with pulse compression setups, such as a chirped mirror (CM) and glass wedge compressor, where the amount of glass traversed by the pulse is an independent variable controlled by insertion of one of the wedges: while the CMs impart negative dispersion, the variable positive dispersion introduced by the wedges will vary the total dispersion experienced by the pulse to be measured. In second-harmonic generation d-scan (SHG d-scan), measuring the SHG signal after the compressor provides a two-dimensional trace of the SHG spectrum vs. insertion. An optimization algorithm is then used to retrieve the spectral phase of the pulse from the measured d-scan trace and calibrated linear spectrum[41]. A recent approach to d-scan retrieval can also be used to obtain the pulse amplitude and phase from the measured trace[43], but in this case, the trace itself must be calibrated.

D-scan has been successfully demonstrated with few-cycle pulses since its inception[22,41], and over octave-spanning single-cycle pulses have been measured directly with SHG d-scan[6,40,43-45]. Apart from its robustness and performance, another important advantage of d-scan is the fact that it directly results in very intuitive traces that provide useful guiding information on the quality of the achieved pulse compression, which motivates our present use of the d-scan trace as a diagnostic tool. For instance, a flat and thin trace is indicative of excellent compression, since this means that for a particular position of the compressor, all spectral components are equally compressed and hence their SHG signal is maximized. If the trace has a tilt, this means that different parts of the spectrum are being compressed for different applied dispersions. In other words, we have a frequency-dependent chirp in the pulses, i.e., the pulses have a residual third-order phase due to uncompensated third-order dispersion. A curved, parabolic-like trace would be indicative of fourth-order dispersion, and so on. In the case of phase oscillations, these will produce spectral modulations and a wavy appearance in the resulting d-scan trace (for a more detailed description of d-scan traces and their interpretation, the reader is referred to refs[22,41]).

In this paper we look for a route to obtain the optimum stable post-compressed pulse from a HCF using the d-scan technique as the compression device and the d-scan trace to univocally identify this route. Numerical simulations and experimental measurements reveal that, in quite general conditions, the best performance results in the post-compressed pulses typically exhibiting a residual third-order phase of nonlinear origin, i.e., a signature tilt (residual TOD) intrinsic to optimum nonlinear propagation within the fiber, in agreement with experimental observations made in several laboratories around the world. Overdriving the HCF above this optimum regime invariably results in an increasingly complex nonlinear spectral phase which renders compression very hard to optimize. These results and behavior have been obtained for different gases and different setup parameters (including HCF length, gas pressure and pulse energy) showing the universality of the existence of the optimum regime. We have also probed the spatial quality of the optimum output pulse, showing a high spatio-spectral homogeneity. Moreover, we have been able to develop a simple theoretical model that explains the main nonlinear effects underlying the optimum regime, which is very useful for finding out the proper parameters to achieve the desired few- to single-cycle pulses. The understanding of the optimum propagation regime, of its intrinsic TOD and its subsequent correction using, e.g., unusual transparent optical media with adequate ratio between second- and third-order dispersion, enables achieving optimized high-quality post-compressed pulses with only minor changes in existing setups.

Results and Discussion

Regimes of nonlinear propagation in the hollow-core fiber. Identification of the optimum d-scan trace

The identification of the optimum parameters of a nonlinear process to obtain a desired output pulse is not an easy task. Fortunately, in the context of nonlinear propagation of ultrashort laser pulses, we have helpful theoretical models to guide us. We used a nonlinear spatio-temporal model (see the Methods section) to simulate the nonlinear propagation of laser pulses in a static HCF filled with different gases (argon, neon and air) and in different conditions (gas pressure; input pulse energy and duration), and their subsequent compression with a d-scan system based on chirped mirrors and glass wedges. All cases studied are representative of real experimental situations and present qualitatively similar dynamics.

In Fig. 1 we present the theoretical (top row) and experimental (bottom row) d-scan traces obtained in a pressure-scan experiment. The left, middle and right column represent the low pressure cases (weak interaction regime), the optimum pressure cases (optimum interaction regime) and the high pressure cases (strong interaction regime), respectively (see the figure caption for the detailed set of parameters). The theoretical d-scan traces are calculated from the on-axis field obtained at the end of the HCF. The d-scan compressor is composed of chirped mirrors followed by a pair of BK7 glass wedges, as in the d-scan setup commonly used for few-cycle pulse compression[41], but taking only into account their group delay dispersion effect and neglecting higher order terms. We assume ultra-broadband chirped mirrors, which introduce −120 fs2 of pure GDD, with the BK7 wedges introducing 46.6 fs2/mm at the central wavelength 780 nm. As we will see in the comparison between simulations and experimental results, neglecting higher order terms in the numerical d-scan stage is a very good approximation because the high-order contributions of CMs and wedges cancel out. This simplification helps us identify the origin of the higher order dispersion terms that may appear in the pulse phase, which must come from the propagation in the HCF because the compression/measurement stage does not introduce them. Therefore, and based on the explanation given in the Introduction section, we will be able to identify any high order phase terms of the output pulse directly from the d-scan trace structure.

Figure 1

The theoretical d-scan traces (top row) correspond to the propagation of a Fourier-limited 25 fs full-width at half-maximum (FWHM) Gaussian pulse, with 0.2 mJ of energy and with the spectrum centered at 780 nm, in a 25 cm long HCF with a core diameter of 250 micron and filled with Ar at 0.2 bar (left), 0.6 bar (middle) and 1.0 bar (right). All the simulations are done assuming that the beam couples perfectly into the fundamental mode of the HCF. The experimental d-scan traces (bottom row) correspond to the propagation of ∼25 fs FWHM pulses with 1 mJ energy centered at 780 nm in a 1.0 meter long HCF with a core diameter of 250 micron filled with argon at 0.15 bar (left), 0.265 bar (middle) and 0.9 bar (right).

The observed similarity between the theoretical and experimental pressure scans indicates that the theoretical model, including the calculation of d-scan traces, takes into account the main effects occurring in the experiment. The results presented in Fig. 1 also help us identify the optimum d-scan traces (those in the middle column). It is clear that the weak interaction regime traces (left column) represent situations where the pulses did not experience enough nonlinear interaction to broaden their spectra, showing narrow d-scan traces in the spectral coordinate. In contrast, the intricate modulations of the strong interaction regime traces (right column) indicate that they correspond to situations where the pulses had such intense nonlinear interaction that, although the spectral broadening was large, the complex spectral phase makes the output pulse useless for most applications. Moreover, this regime is easily identified in the experiments due to a poor output stability. We can therefore define the optimum interaction regime as the one that generates d-scan traces similar to those of the middle column, which represent situations where the pulses are spectrally broadened but in a regime in which the acquired spectral phase can be adequately compensated for. The retrieved pulses obtained from the optimal d-scan traces presented in Fig. 1 have durations of ∼4.1 fs FWHM for the theoretical case (top row) and ∼4.0 fs FWHM for the experimental case (bottom row). These optimum d-scan traces already correspond to sub-two-cycle pulses, and additional compensation of their intrinsic residual TOD[44] allows obtaining high-quality pulses in the single-cycle regime[6,43,44].

The exact limits of the three identified propagation regimes are somewhat flexible but, as shown in Fig. 1, they can be easily identified using the d-scan as the guiding tool. There is a large number of parameters that one can use to reach the optimum regime in the laboratory: input pulse energy, input pulse temporal duration (chirp), coupling conditions (focusing, mode and numerical aperture matching), gas type and pressure. Using some of them one can get into the optimum regime to achieve the optimal post-compressed output pulse, which will not be the one with the broadest spectrum, but the shortest among the less temporally-structured obtained pulses, which we have been able to univocally identify through the resulting d-scan trace.

Optimum d-scan trace properties: Universality, Spatial and Spectral structure, and Physical origin

In the previous results, we have been able to identify the optimum d-scan trace when using argon. We now apply the theoretical model to verify that similar optimum d-scan traces appear using different gases and experimental parameters. Figure 2 shows the theoretical optimum d-scan traces obtained at the end of a HCF using argon (left), neon (middle) and air (right). The retrieved output pulse duration for the three cases is ∼4.1 fs, ∼4.0 fs and ∼4.7 fs FWHM, respectively. All these traces were identified after simulating pressure scans as the one presented in the top row of Fig. 1. When using argon or neon we obtained almost identical optimum traces, while when using air there are some variations due to the presence of the Raman effect, which is not present when using atomic gases. Surprisingly, regardless of the different conditions (fiber length, pulse energy, gas pressure and gas type), we have always been able to identify the optimum d-scan trace at the end of the HCF. This universal behavior, together with the particular structure of the optimum d-scan trace as a univocal fingerprint, is a very useful tool both for optimizing existing systems and to promote the spreading of few-cycle pulse systems for many more applications that can benefit from a stable and reproducible source.

Figure 2

Simulated d-scan traces for a pulse with 0.2 mJ energy propagating in a 25 cm long HCF filled with argon at 0.6 bar (left), a pulse with 0.4 mJ propagating in a 50 cm long HCF filled with neon at 2.0 bar (middle), and a pulse with 0.15 mJ propagating in a 75 cm long HCF filled with air at 0.5 bar (right). All the simulations assume a perfect coupling into the fundamental mode of a HCF with a core diameter of 250 micron, and an input pulse with 25 fs FWHM, spectrally centered at 780 nm.

For completeness we present a movie (see Supplementary Material) showing the evolution of the d-scan trace at the end of the HCF during a pressure scan with argon, together with the evolution of the wedge-insertion-dependent post-compressed pulse. The left plot in the movie shows the particular features of the d-scan trace along all three stages. The right plot shows the final compressed pulse that can be obtained, displaying the continuous evolution of the structure of the output pulse.

The structure of the optimum d-scan trace gives us important information about the output pulse obtained from the post-compression system. The first distinctive feature of the optimum d-scan trace is a slightly negative slope which corresponds to a negative TOD. This is the remaining residual TOD accumulated during the nonlinear propagation of the pulse inside the HCF. According to the simulations, this TOD mostly comes from the self-steepening effect, similarly to what occurs in post-compression setups for pulses in the mid-infrared spectral region[46]. This tilt of the d-scan trace is very helpful to find experimentally the optimum regime to obtain single-cycle post-compressed pulses. To our knowledge, this remaining TOD is clearly present in experimental HCF post-compressed pulses obtained in different laboratories around the world when optimizing for the output pulses (e.g., United Kingdom[45], Portugal[44,45,47], France[48], Sweden[42,49], Germany[50], Canada[51], China[52], Austria[46] and United States of America[6,39], as some examples). Those experiments were done using HCFs with different characteristics (length and/or core diameter), filled with different gases (Ar, Ne or He) and with lasers in the NIR and mid-IR regions. All these observations support that the route to obtain an optimum post-compressed pulse described here is universal.

Another very important feature of the output pulse related to the optimum d-scan trace is spatial homogeneity. All the theoretical d-scan traces presented until now have been calculated from the on-axis field at the end of the HCF. One would expect that, as the optimum regime corresponds to a moderate nonlinear interaction regime, the resulting output beam should have a quite good spatial homogeneity, as it has in fact been proved experimentally in similar conditions[47]. Figure 3 shows the theoretical far-field d-scan calculated from a spatial selection of the field obtained at the end of the HCF. The parameters here are the same used to obtain the optimum d-scan trace in argon (leftmost case of Fig. 2). The far-field distribution is calculated by doing the Hankel transform integrating spatially over the selected part of the beam. Then the d-scan trace is finally calculated using the far-field corresponding to k⊥ = 0, i.e. the zero divergence far-field. Figure 3a shows the d-scan trace of the far-field when taking only into account the on-axis near field, which is basically the same presented using the field on-axis (left picture of Fig. 2). Figure 3b shows the d-scan trace of the far-field when integrating the near-field up to 62.5 μm (half of the HCF core radius, which contains 94% of the output energy). As it can be observed, the d-scan trace essentially retains the same shape as before, being slightly narrower in the spectral direction and more stretched along the dispersion axis. Figure 3c shows the d-scan of the far-field obtained by integrating up to 100 μm (which accounts for 99% of the output energy). In this case the spectral narrowing of the d-scan is more pronounced, as we are taking into account the most external part of the beam that sees less nonlinearity. Note that even this last d-scan trace preserves the TOD signature related to the optimum regime, although the optimal pulse duration changes from ∼4.1 fs in the first case to ∼4.7 fs and ∼5.3 fs in the second and third cases, respectively. Finally Fig. 3d shows the beam fluence at the end of the HCF with circles indicating the integrated area used to calculate the far-field d-scan for the three cases shown, which helps to visualize the good homogeneity of the beam obtained under the optimum parameters. The smooth output spatial profile also shows that the post-compressed pulse is suitable for applications which require good beam quality.

Figure 3

Simulated far-field d-scan traces obtained by integrating the beam with an on-axis aperture with a radius of 1 μm (taking only into account the on-axis field) (a) 62.5 μm (b) and 100 μm (c). All the d-scan traces correspond to the case of a 25 fs FWHM input pulse, spectrally centered at 780 nm, with 0.2 mJ of energy, propagating in a 25 cm long HCF filled with argon at 0.6 bar. Like before, a perfect coupling into the fundamental mode of the HCF is assumed. In (d) the fluence of the beam at the end of the HCF is shown, with circles indicating the different apertures used to calculate the three far-field d-scans.

The identification of the optimum regime and the spatial homogeneity of the resulting output beam indicate that this regime is achieved for moderate nonlinear interaction conditions (not presenting, for instance, self-focusing dynamics, strong gas ionization or similar phenomena related to intense nonlinear interactions). To unveil the main nonlinear effects underlying this regime, we can use a theoretical model based on propagation equations (this model is discussed in Chapter 4 of ref.[53]). As we will show below, this model is useful for estimating the parameters of the optimum propagation regime and, consequently, for obtaining the desired few- or single-cycle post-compressed pulses. The model assumes that self-phase modulation (SPM) and linear losses are the only phenomena affecting the pulse during propagation. With these assumptions we are able to anticipate the maximum induced spectral broadening after propagating a distance L in the fiber, which in the case of a Gaussian pulse turns out to be Δω = 0.86Δω0B, where Δω0 and B represent the initial spectral width and the maximum nonlinear phase shift, respectively. Assuming now that the input and output pulses are Fourier-limited, which is valid if we are able to compensate properly the spectral phase of the output pulse, their spectral widths are directly related to their respective temporal widths by and . Thus we can obtain an estimate of the required maximum nonlinear phase shift, i.e. the maximum B-integral, for any given initial and final (compressed) pulse durations:

This basic model indicates that to broaden the spectrum of Fourier-limited input pulses of similar duration towards the single-cycle regime, we will always need to accumulate roughly the same maximum amount of the B-integral. This way, the B-integral becomes a very useful tool for predicting adequate values of some of the parameters.

We now proceed to finding the accumulated B-integral of a pulse with energy E, temporal duration t, and spectrum centered at λ0, propagating in a HCF of length L and with a core of radius r. We assume that the beam couples perfectly to the fundamental mode of the HCF, with well-known losses denoted by α[54]. Inside the HCF we have a gas at constant pressure, with nonlinear index (1 − x)n2, where x indicates the amount of Raman effect (x = 0 for atomic gases and for molecular gases like air). The maximum accumulated B-integral then reads:where is the maximum input intensity. We should point out another argument that further confirms the suitability of the B-integral as a good parameter for determining the optimal compression regime, based on the scaling properties of nonlinear wave equations for pulse propagation in gases. From Eq. 2, in the case of weak linear absorption (), we have . Following Heyl et al.[49,55], the propagation equations remain invariant if we transform the input parameters with a scaling dimensionless parameter η according to certain rules. In our case, the scaling relations E → η2E, r → ηr, n2 → η−2n2 and L → η2L make B invariant under this transformation.

Applying Eq. 2 to all the simulations presented in Fig. 1 we obtain similar values for the B-integral, (4.2, 4.4 and 5.6 for the simulations corresponding to argon, neon and air, respectively), as predicted by the basic model described above. In fact, Eq. 1 gives a constant value , when using ΔT = 25 fs and ΔT corresponding to a 1.5-cycle, 3.9 fs pulse, proving that this simple model includes the main effects underlying this regime: SPM and linear losses. The optimum B-integral, defined as the lower limit value of B, , provides a useful hint for finding optimum propagation parameters for input pulses of this duration. Taking into account that the nonlinearity depends linearly on the gas pressure, we have used Eq. 2 to calculate, for example, the pressure needed to obtain the optimum d-scan trace at the end of a HCF of fixed length L and core radius r, filled with some gas with nonlinear refractive index (being the nonlinear parameter at 1 bar pressure), with an input pulse of energy E and a temporal duration t:

We have checked that the parameters obtained from Eq. 3 bring us to a situation very close to the optimal d-scan trace for argon, neon and air, over a wide range of parameters, and whenever the propagation occurs in the moderate nonlinear regime. These expressions demonstrate the tunability of the optimal trace with respect to key parameters such as pressure, type of gas and input energy or pulse duration, setting this procedure as a universal route for obtaining few-cycle pulses in the visible-NIR spectral region.

We have also estimated the maximum B-integral for some experimental cases[6,39,44,47], besides the experiments presented in this work, obtaining a value (range of values from 9 to 17 with an average of 14). To calculate these values we have estimated that the energy coupled into the HCF is 40% of the available input energy and we have assumed that the input pulse is Fourier-limited. Even with all these assumptions, which are not always fulfilled in experiments, we have obtained similar values of the experimental maximum B-integral, which demonstrates that this parameter is also useful for bringing an experiment close to the optimum regime. We can see that the estimated values of the B-integral in experiments are consistently larger than in numerical simulations, i.e., B > B. This is a consequence of the less effective nonlinear interaction present in real experiments compared to the numerical models, due to factors such as the imperfect spatial intensity distribution of the laser beam, as reported in the literature since the early studies of nonlinear ultrashort pulse propagation (see ref.[56]).

TOD compensation of the optimal d-scan trace

The residual TOD which appears in the optimal propagation regime can be further compensated for, not only to reduce the output pulse duration but also to increase the amount of energy in the pulse. This can be done, for example, by using an adequate transparent medium such as water[42,44,45], z-cut KDP[43] or z-cut ADP[6], which directly enables achieving high-quality pulses in a more strict single-cycle regime. In spite of the relatively small magnitude of this TOD (e.g., −40 fs3 reported in ref.[44]), its effect on a few- and single-cycle pulse can be dramatic: its compensation enabled going from a 3.8 fs pulse (with satellite pulses due to TOD) to a near-transform-limited 3.2 fs 1.3-cycle pulse, which was accompanied by a significant improvement in pulse contrast and peak intensity[44].

In the center plot of Fig. 4 we show the d-scan trace obtained after compensating the residual TOD by adding +20 fs3 (at the central wavelength of 780 nm) in the measurement/compression stage, compared with the non-compensated case (left plot). We see that TOD compensation effectively tilts the d-scan trace so that its average negative slope practically disappears. To illustrate the improvement in intensity throughput after TOD compensation we show on the right plot of Fig. 4 the direct output pulse and the TOD-compensated pulse, together with the Fourier-limited pulse. The used glass insertion in this case is the one which results in maximum intensity in the d-scan trace. Clearly, the pulse corresponding to the TOD-corrected trace has better features: it presents an increase in peak intensity of the order of 33%, a reduction of the temporal FWHM pulse duration from ∼4.1 fs to ∼2.9 fs, which is in very good agreement with experimental results, and smaller secondary temporal structures than the non-corrected case, as also confirmed by experiments.

Figure 4

D-scan traces of (left) the optimum pulse for the simulation parameters of the top middle graph in Fig. 1 (Ar-filled HCF with 0.6 bar pumped with Fourier-limited ∼25 fs pulses) and (center) the same pulse with additional TOD compensation (+20 fs3). The dashed white lines are a visual guide to see the change of slope in the trace associated to TOD compensation. (right) Output pulse in the time domain, before (red line) and after (blue line) TOD compensation. The Fourier-limited pulse (green line) is also shown for reference.

Conclusion

In this work we have presented and analyzed a route to obtain optimum hollow-core fiber post-compressed pulses using the d-scan setup as the compression and diagnostic device. We have demonstrated that by changing the gas pressure and/or input pulse energy, one is able to identify the optimal post-compression parameters to obtain the shortest, less structured and stable output pulses. This optimal setup can be univocally identified by the d-scan trace of the output pulses, which shows a marked TOD feature, whereas overdriving the HCF above this optimum regime invariably results in an increasingly complex nonlinear spectral phase, which renders compression very hard to optimize. We have verified that the optimum d-scan trace can be found for different gases and setup conditions, highlighting the universality of this phenomenon. Moreover, we have also proved the good spatio-spectral characteristics of the output pulse obtained under these conditions. For the optimum propagation regime, the remaining TOD can be further corrected, for example by simple linear propagation in an adequate transparent medium, which improves the temporal shape and the peak power of the output pulse. The identified propagation regime and approach enable the generation of stable and high-quality few- to single-cycle pulses, which has direct implications in the performance of current and new HCF pulse post-compression systems and will help improve and extend the applications of these extreme light sources in many fields of science and technology.

Methods

Simulations

To study the nonlinear propagation of a laser pulse in a HCF we have implemented the standard nonlinear envelope propagation equation[57]. The model includes the spatial and temporal pulse dynamics. We use a local frame moving with the pulse T = t − z/v, being v the group velocity of the pulse, and assume cylindrical symmetry (r radial coordinate, z axial coordinate). The propagation equation for the temporal envelope of the pulse, A(r, z, T), is

The first part of the propagation equation, , represents the linear propagation effects: diffraction, dispersion and linear losses. The other part, , represents the nonlinear propagation effects, which include self-phase modulation, Raman scattering, ionization, losses due to the ionization process and plasma absorption, and self-steepening.

To solve Eq. 4 we use a split-step method dividing each propagation step into two sub-steps[53]. The first sub-step consists in applying only the linear effects by decomposing the input pulse into the EH1 modes of the hollow-core fiber[54,58]where represents the Fourier Transform of A(r, z, T). The coefficients of the decomposition can be calculated by doing the inverse Hankel Transform of the spatial beam distribution in the core of the HCFbeing a the core radius of the HCF and α the m-zero of J0(x), where J is the Bessel function of the first kind of order ν. We solve Eq. 6 by using the discrete Hankel Transform scheme proposed in[59]. The complete linear propagation in the HCF is simulated by using the complex propagation coefficient of each mode, β(ω)[54], as shown in Eq. 7. The real and imaginary parts of β(ω) take into account all the dispersion and losses of the m-mode of the HCF, respectively

The second sub-step of the method consists in applying the nonlinear effects. Separating , a mathematical expression for each term can be obtained[57], as given below. For the self-phase modulation, Raman scattering and self-steepening, we have

In Eq. 8, k0 = n02π/λ0, with λ0 the central wavelength of the pulse. n2 the nonlinear refractive index, ω0 = 2πc/λ0 and K(T − t) representing the SPM together with the Raman scattering, which has the following form (1 − x)δ(T − t) + x/τexp(−(T − t)/τ), where x is the ratio between the SPM and the Raman effect, and τ fs the characteristic time for the Raman response. For the ionization, we havewhere σ is the cross section for the inverse Bremsstrahlung that depends on the collision time (τ), the critical density of the medium and the central frequency of the laser pulse[57]; ρ represents the ionized electron density, whose evolution is governed by ∂ρ/∂t = W(|A|2)(ρ − ρ), where W(|A|2) is the ionization rate calculated using the PPT model[60] and ρ is the atomic density of the medium. Finally, the absorption term iswhere U represents the ionization potential of the gas.

Experiments

The experiments were performed by employing 23–25 fs Fourier-transform-limited laser pulses centered at 780 nm. These pulses were generated with a 1 kHz Ti:Sapphire CPA laser system (Femtolasers FemtoPower Compact PRO HE CEP) which is part of the CLPU facility. The maximum pulse energy available was 2.5 mJ and we adjusted it to 1 mJ in our experiments. The laser pulses were focused by a spherical mirror (1 m focal length) into a hollow-core fiber (HCF) with an inner diameter of 250 micron and 1 meter length. The HCF was filled with argon gas at different pressures. The output pulses were compressed and measured using a d-scan system from Sphere Ultrafast Photonics. The scanning SHG d-scan setup includes a double-angle chirped mirror set (Ultrafast Innovations GmbH) and motorized BK7 glass wedges to induce the dispersion scan.

Data availability

The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.

Video 1