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.

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][34][35][36][37][38][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][39][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][44][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 secondand 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 fs 2 of pure GDD, with the BK7 wedges introducing 46.6 fs 2 /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.
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.
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  . 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.
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 Δω max = 0.86Δω 0 B max , where Δω 0 and B max 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 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 in , temporal duration t p , and spectrum centered at λ 0 , propagating in a HCF of length L F and with a core of radius r F . 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 R )n 2 , where x R indicates the amount of Raman effect (x R = 0 for atomic gases and x 0 6 R ∼ . for molecular gases like air). The maximum accumulated B-integral then reads: 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 ( L 1 F α  ), we have B I x n L 2 ( 1 ) . 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 in → η 2 E in , r F → ηr F , n 2 → η −2 n 2 and L F → η 2 L F 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,  B 4 sim (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  B 7 max , when using ΔT i = 25 fs and ΔT f 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 sim , B 4 opt sim = , 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 F and core radius r F , filled with some gas with nonlinear refractive index ⁎ − x n (1 ) R 2 (being ⁎ n 2 the nonlinear parameter at 1 bar pressure), with an input pulse of energy E in and a temporal duration t p : 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 B 9 exp  (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 exp > B sim . 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 fs 3 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 fs 3 (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.

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 SCIeNtIFIC REPORTS | (2018) 8:2256 | DOI:10.1038/s41598-018-20580-1 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.

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 g , being v g 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 .Th e first part of the propagation equation, L , represents the linear propagation effects: diffraction, dispersion and linear losses. The other part, N , 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 EH 1m modes of the hollow-core fiber 54 , 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, k 0 = n 0 2π/λ 0 , with λ 0 the central wavelength of the pulse. n 2 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/τ K exp(−(T − t)/τ K ), where x is the ratio between the SPM and the Raman effect, and τ K fs the characteristic time for the Raman response. For the ionization, we have ioniz C 0 0 1 where σ is the cross section for the inverse Bremsstrahlung that depends on the collision time (τ C ), 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 )(ρ at − ρ), where W(|A| 2 ) is the ionization rate calculated using the PPT model 60 and ρ at is the atomic density of the medium. Finally, the absorption term is where U i 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.