Introduction

Accurate knowledge of the on-target pulse intensity is essential for the correct interpretation of high-field experiments involving highly nonlinear processes. The on-target intensity can be estimated directly via different gas ionization processes, such as the intensity scaling of the above-threshold ionization1,2, the yield of single or multiple charged ions3,4, the photoelectron/photoion momentum distribution from ionization by circularly polarized laser fields5 and similar. However, these techniques are limited to intensities significantly below the relativistic limit (\(1.37 \times 10^{18}\;{\mathrm{W}}/{\mathrm{cm}}^2\) at a central wavelength of 1 µm) and require a precise theory of the ionization process as well as knowledge of the temporal and spatial pulse shape to achieve reasonable precision of the averaging of the electron/ion impact over the intensity distribution5. Therefore, in high-intensity systems, the value of the on-target intensity is commonly estimated from the spatial and temporal power/fluence distributions, which are measured separately6,7. While there are reliable approaches for the detection of pulse energy and spatial distribution8, precise characterization of the on-target temporal structure remains a challenge, especially for intense, few-cycle pulses.

Achieving reliable detection of the spatial intensity distribution is simpler than achieving that of the temporal distribution, as it can be done directly with a CCD camera in combination with a microscope objective8, and the result is not affected by dispersive optical elements in the beam path (at least when they are of good enough quality and the beam is properly attenuated to prevent nonlinear distortions). One should note that a reliable measurement requires a high dynamic range (HDR) image to detect the energy contained in a possible background around the main peak. This image can be obtained with a combination of an HDR CCD camera and the additional extension of the dynamic range by the sequential measurement of unsaturated and saturated images with different attenuations.

In contrast, the temporal characterization of ultrashort optical pulses is always indirect and suffers from dispersive optics in the diagnostic beam path. It can be performed using different techniques involving reconstruction algorithms, such as FROG (frequency-resolved optical gating)9, SPIDER (spectral phase interferometry for direct electric-field reconstruction)10, a dispersion scan11, spatially resolved Fourier transform spectrometry12, and their modifications. The HDR temporal characterization is as important as the HDR focal beam profile diagnostics. Common and reliable methods of measuring the HDR temporal structure and determining the energy content in the main peak are third-order autocorrelation13, cross-correlation14 and self-referenced spectral interferometry, which allows single-shot implementation15. Note that the temporal contrast diagnostics and the pulse duration diagnostics are complementary. The former determines the energy content in the main pulse as well as in pre- and post-pulses and in the background but has poor temporal resolution (in the best cases, ~20–50 fs). Therefore, it cannot be used for the peak power evaluation of pulses with durations below ~100 fs. Thus, pulse duration diagnostics is required for the complete measurement of the temporal structure and corresponding peak power.

The problem of most pulse duration characterization approaches is that they cannot be applied at high intensities, which makes their results not always representative for the full-energy pulses. In addition to the difficulty of attenuation itself, the problem here lies in the dispersion, which is introduced in beamlines between the attenuation and diagnostics. Usually, it is necessary to split the beam somewhere before the target, strongly attenuate it, and guide it outside a vacuum chamber through a window, as most diagnostic devices are not vacuum compatible or require fine alignment, which makes their in-vacuum realization extremely complicated. All these procedures normally introduce some dispersion that cannot always be precisely traced back or post-compensated for accurate evaluation of the on-target pulse duration. For multi-cycle pulses (>10 fs pulse duration at ~1 µm central wavelength), generated, for example in Ti:Sa laser systems, this uncertainty is often tolerable. However, it is certainly not tolerable for few-cycle pulses generated in the emerging high-peak-power OPCPA (optical parametric chirped pulse amplification) systems16,17,18. Thus, many high-field experiments, especially with these systems, will profit from a direct on-target temporal characterization technique compatible with relativistic intensities 19,20.

A suitable approach can be a known temporal characterization technique involving second or higher-order harmonic generation, with the harmonic source being based on the relativistic surface high harmonic generation instead of a nonlinear crystal. Relativistic surface harmonics21,22 are compatible with very high pulse intensity and thus fulfill the main goal. There is, for example, a theoretical proposal of implementing an autocorrelator based on the detection of the third harmonic at relativistic intensities20. The temporal characterization approach to be used should ideally be reliable, robust and easy to implement but at the same time should be accurate and support ultrashort pulses. The dispersion scan (D-scan)11 method fulfills all these requirements and thus seems to be one of the best candidates. D-scan is based on the measurement of the dependence of the second harmonic spectrum on the varying dispersion in the fundamental beam path. Unlike autocorrelation approaches19,20, it requires no beam splitting and delaying of pulse replicas relative to each other and has been demonstrated to be suitable for the accurate characterization of few-cycle pulses11,23. Since our approach uses the D-scan scheme combined with relativistic surface second harmonic generation (SHG), we call it the relativistic surface second harmonic generation dispersion scan (RSSHG-D-scan).

This paper introduces the RSSHG-D-scan (Fig. 1). We demonstrate both experimental and PIC (particle in cell) simulation results, proving the feasibility of the technique, and discuss optimum experimental parameters and limitations of the method.

Fig. 1: Scheme of the experimental setup.
figure 1

The p-polarized pulses are focused onto a BK7 target with a gold-coated ø1ʺ f/1.6 off-axis parabolic mirror. The reflected fundamental radiation and the generated harmonics are re-collimated with an aluminum-coated ø2ʺ f/1.3 off-axis parabolic mirror and guided with a set of mirrors and wedges outside the vacuum chamber for diagnostics. The inset shows a typical measured spectrum

Results

SHG saturation intensity

Saturation must be avoided in any measurement setup, including the RSSHG-D-scan. In a standard D-scan, SHG saturation needs to be avoided to ensure the quadratic scaling of the second harmonic signal with the intensity of the fundamental radiation, which allows for precise pulse reconstruction. This is normally automatically ensured by using a very thin nonlinear crystal, which is already required for broadband phase matching. In a RSSHG-D-scan, saturation of the second harmonic signal needs to be avoided as well. To determine the intensity range with quadratic SHG scaling, we measured the dependence of the generated second harmonic signal on the intensity of the incoming pulse (Fig. 2a). The experimental results demonstrate the quadratic scaling of the SHG efficiency for the normalized vector potential of an incoming field of \(a_{0} \,{\lesssim}\, 1\) (see Methods for determination and more information); however, for higher intensities, the SHG saturates, and the ratio becomes linear. This experimental finding is in good agreement with the results of 1D PIC simulations24 (Fig. 2b), where both quadratic and linear scaling regions as well as the turning point between them are well reproduced.

Fig. 2: Nonlinearity of the RSSHG process.
figure 2

Experimental (a) and PIC simulation (b) results show very good agreement with each other. The simulation is performed for the plasma scale length of 0.2λ. The axes are normalized to \(I_{\omega \_{\mathrm{in}}}(a_0 = 1)\); thus, the x-axis is \(I_{\omega \_{\mathrm{in}}}/I_{\omega \_{\mathrm{in}}}(a_0 = 1)\), and the y-axis is \(I_{2\omega \_{\mathrm{out}}}/I_{\omega \_{\mathrm{in}}}(a_0 = 1)\). The corresponding vector potential is \(a_0 = \sqrt {I_{\omega \_{\mathrm{in}}}/I_{\omega \_{\mathrm{in}}}(a_0 = 1)}\), and \(I_{\omega \_{\mathrm{in}}}(a_0 = 1)\) is the intensity corresponding to a0(ω) = 1. Both the experimental and simulation results demonstrate quadratic scaling (orange lines) for the normalized vector potential of the incoming field of \(a_0 \,{\lesssim}\, 1\); however, for higher incoming intensities (\(a_0 \,{\gtrsim}\, 2\)), the SHG saturates, and the ratio becomes linear (green lines)

A D-scan measurement at intensities in the saturation region results in a trace that cannot be accurately retrieved by a standard algorithm since it does not take saturation into account and thus fails. The simulated dependence of the reconstruction accuracy on the incoming pulse intensity and thus on the saturation level is presented in Table 1. The simulations were performed for Lp = 0.05λ. According to these data, increasing the incoming field strength over a0 = 1 causes an increase in the peak intensity reconstruction error. There are a few approaches to overcome this problem and to implement the suggested technique in an ultra-relativistic laser system. One option is a modification of the retrieval approach by taking SHG saturation into account, but this requires substantial theoretical support and will have problems with unambiguity of the retrieval result, especially at very high intensities (a01). Therefore, a better option is a well-controlled and characterized scaling down of the on-target intensity in the characterization setup. This can be realized, for example, by measuring the RSSHG-D-scan signal in the focus of another parabolic mirror with a longer focal length, or the incoming beam can be attenuated with a reflective wedge (or a set of wedges) to fulfill the condition of \(a_0 \lesssim 1\). Unlike other contemporary temporal characterization approaches, the diagnostics can be implemented directly in a target chamber without an additional beamline, which would introduce dispersion.

Table 1 Simulated dependence of the pulse reconstruction accuracy on the peak intensity of the incoming pulse and thus on the SHG saturation level (Fig. 2)

Another outcome of the presented results (Fig. 2) is the SHG efficiency optimum being at a0 ≈ 2–3. A further increase in the intensity will not improve the SHG efficiency, but it will result in energy transfer into higher harmonics, which saturate later25.

Simulation results

The theoretical feasibility of the RSSHG-D-scan was numerically investigated by comparing the RSSHG-D-scan with the standard SHG-D-scan using an input pulse with a nontrivial spectral phase (Fig. 3a). The standard SHG-D-scan was performed assuming perfect phase-matching conditions and using the following equation:

$$I_{2\omega }(\omega ,L) = \left| {{\int} {E_{2\omega }(t,L)e^{ - i\omega t}dt } } \right|^2$$
(1)

with

$$E_{2\omega }(t,L) = \left( {{\int} {E_\omega (\omega )e^{i\omega t + iLk(\omega )}d\omega } } \right)^2$$

where E2ω and I2ω are the field and intensity of the second harmonic signal, respectively, Eω is the spectral field distribution of the incoming fundamental signal, L is the thickness of the material that is being scanned, and k(ω) = n(ω)ω/c is the wave vector, with n(ω) being the material dispersion. The RSSHG-D-scan was simulated with a 1D PIC code24. Both results, i.e., the second harmonic spectral intensity as a function of the inserted material length, are shown in Fig. 3b, c and demonstrate a very good agreement with the RMS difference of just G = 0.03 (see Methods section for definition), which represents a good level of consistency26. This consistency infers that all retrieval algorithms developed for the standard D-scan are directly applicable to the RSSHG-D-scan. Therefore, the RSSHG-D-scan method can directly profit from the advances in the retrieval approaches of the D-scan technique that have been developed over the last decade and are very well established.

Fig. 3: Theoretical proof of the feasibility of the RSSHG-D-scan by comparing it with the conventional SHG-D-scan method in simulation.
figure 3

a Spectral intensity and phase (in terms of group delay) used in the simulations. b Simulation of the SHG-D-scan (see main text for details). c PIC simulation of the RSSHG-D-scan (simulation parameters: a0 = 1 for the best compressed pulse, Lp = 0.05λ). L is the thickness of the fused silica, determining the dispersion introduced in the fundamental pulse

Experimental results

The experimental test of the RSSHG-D-scan was performed using the setup schematically shown in Fig. 1 and described in more detail in the Methods section. To demonstrate the consistency of the RSSHG-D-scan with the standard SHG-D-scan, we characterized the same pulse using both techniques in parallel. As mentioned in the introduction, the measurement of the direct on-target pulse duration with a SHG-D-scan setup typically requires more effort. In particular, any dispersive medium needs to be avoided, which means that the setup has to be completely in vacuum and that no dispersive optical elements should be used. For this reason, only uncoated or metal-coated reflective optics were used in the implemented SHG-D-scan diagnostic setup. Second harmonic generation was performed in the thinnest commercially available BBO (β barium borate) crystal, with a 5 µm thickness, to ensure very good broadband phase matching. The RSSHG-D-scan, on the other hand, requires no special precautions, apart from the use of appropriate optics in the diagnostic beamline, which have no strong absorption in the spectral range of interest. The obtained experimental results are presented in Figs. 4 and 5. The relatively rough step size in the scan was mainly determined by technical issues caused by the necessity of re-optimization of the OPCPA stages for each material thickness (see Methods section for more details).

Fig. 4: Experimental results of the SHG-D-scan using a 5-µm-thick BBO crystal.
figure 4

Measured (a) and retrieved (b) trace. The dispersion was changed by varying the thickness (L) of the fused silica substrates, which is measured relative to an arbitrarily defined 0 position, with smaller values corresponding to less material. c Retrieved spectral phase. The light-blue area around the main curve represents the uncertainty of the reconstructed spectral phase determined by the standard deviation of the retrieved function. d Reconstructed temporal pulse shape

Fig. 5: Experimental results of the RSSHG-D-scan.
figure 5

Measured (a) and retrieved (b) trace. The intensity of the best compressed pulse in the experiment corresponds to a0 ≈ 1. The dispersion was changed by varying the thickness (L) of the fused silica substrates, which is measured relative to an arbitrarily defined 0 position, with smaller values corresponding to less material. c Retrieved spectral phase. The light-blue area around the main curve represents the uncertainty of the reconstructed spectral phase determined by the standard deviation of the retrieved function. The shape of the phase resembles the optical parametric phase16; thus, it is caused by a slight mismatch between the final experimental parameters and the estimated parametric phase included in the design dispersion of the chirped mirrors. d Reconstructed temporal pulse shape. The retrieved full width at half maximum (FWHM) pulse duration is just 2% larger than the transform-limited one

The retrieval approach based on the differential evolution algorithm26 was used. The RMS reconstruction error in both cases was 0.04. The uncertainty of the reconstructed spectral phase (Figs. 4c and 5c) determined by the standard deviation of the retrieved phase corresponds to the pulse duration uncertainty of 1%. The retrieved phases are in good agreement with each other as well as with the SHG-FROG results of our previous work on the laser system itself16. Note that this latter result has an opposite sign due to the sign ambiguity of the SHG-FROG method. The shape of the residual spectral phase resembles the optical parametric phase16 and is most likely caused by the matching of the dispersion of the chirped mirror compressor to the entire system, which is quite well optimized, but it is not perfect.

Although, as shown above, it is possible (in principle) to perform a reliable measurement with the SHG-D-scan, it has many drawbacks compared with the RSSHG-D-scan. First, the RSSHG-D-scan is a more direct technique, ensuring that the pulse duration is measured directly on-target. Second, the efficiency of the second harmonic generation in a thin crystal is small, causing detection problems, especially in the presence of a parasitic SHG background. For example, nonlinear crystals used in OPA systems quite often support some collinear second harmonic generation that contains not more than a fraction of a percent of the total energy of the amplified pulse and thus causes no problems for experiments, but this parasitic SHG signal from a few mm-thick OPA crystals can be higher than the SHG-D-scan signal from the few µm-thick diagnostic crystal needed for good phase matching. As will be discussed in more detail below, RSSHG shows a very high conversion efficiency, which makes it very tolerant to any parasitic background. Note that the problem of the parasitic SHG background can be avoided in a standard D-scan using the cross-polarized wave XPW-D-scan27, but this approach suffers from other problems, such as dispersion of the polarizers and polarization purity of the incoming beam. Third, a potential error in the estimation of the introduced dispersion will result in a pronounced increase in the reconstruction error, which will indicate that something is wrong. In contrast, in any out-of-vacuum diagnostic approach, an incorrect estimate of the beamline dispersion directly affects the measured pulse duration without any indication of the inconsistency. In addition, in our setup, the RSSHG-D-scan was simpler to design, build and align than the classical SHG-D-scan setup.

RSSHG conversion efficiency

RSSHG can support high conversion efficiency for ultrashort high-energy pulses28, as one can directly see from the typical detected spectrum presented in Fig. 1. At the optimum pre-plasma scale length, the integrated detected energy in the reflected beam in the second harmonic range was of the same order as in the fundamental range. In particular, the measured energy was 5 mJ for both the fundamental and second harmonic beams after taking into account the beamline transmission. An FEL700 filter was used to measure the fundamental energy, and a BG40 filter, to measure the second harmonic. The incoming pulse energy was 22 mJ, which is equivalent to a0 ≈ 3. Thus, the measured energies correspond to 23% of the second harmonic conversion efficiency. The rest of the incoming energy is partially converted into higher harmonics that are outside the detection range of the used spectrometers and might be partially absorbed or transmitted by the target.

The experimentally demonstrated high SHG conversion efficiency is quite attractive for many applications requiring broadband high-energy second or even higher-order harmonic radiation29,30,31,32,33, especially in combination with their good phase-locking properties34,35. For example, RSSHG can be used for the generation of very short optical wave forms by synthesizing the fundamental radiation with the generated second harmonic and potentially higher harmonics in a way similar to the generation of light transients from a supercontinuum after hollow-core fibers 30.

The SHG conversion efficiency depends on the pre-plasma scale length, as was demonstrated for higher harmonics in other experiments on relativistic surface high harmonic generation36,37,38,39,40, and the 23% efficiency presented above was achieved only under optimum conditions. The performed pre-pulse scan (Fig. 6) shows a clear optimum around the plasma scale length of \(L_{\mathrm{p}} \approx 0.1\lambda - 0.15\lambda\), which agrees well with the predictions of the PIC simulations (Fig. 6), having a maximum at \(L_{\mathrm{p}} \approx 0.1\lambda\). Therefore, the experiments, aiming at efficient frequency doubling, should be performed at a plasma scale length of ~0.1λ.

Fig. 6: Measured (blue round dots) and simulated (green squares) dependence of the SHG efficiency on the pre-pulse delay and corresponding plasma scale length.
figure 6

SHG efficiency is determined as \(U_{2\omega \_{\mathrm{out}}}/U_{\omega \_{\mathrm{in}}}\), with \(U_{\omega \_{\mathrm{in}}}\) being the incoming fundamental energy and \(U_{2\omega \_{\mathrm{out}}}\) being the reflected energy in the SHG spectral region. Error bars originate from the standard deviation of multiple measurements. The results demonstrate the excellent agreement between the experiment and simulation and have a clear optimum at the plasma scale length of Lp ≈ 0.1λ

Discussion

We have presented RSSHG-D-scan—a technique for the direct on-target temporal characterization of high-energy, few-cycle optical pulses at relativistic intensity. The feasibility of the approach for the temporal characterization of few-cycle pulses was proved both experimentally and in PIC simulations. The direct on-target temporal characterization of 6 fs pulses with few-TW peak power was experimentally demonstrated. RSSHG-D-scan will be especially suitable for pulse characterization in emerging high-energy, few-cycle OPCPA-based systems and can play an important role in the optimization of pulse compression in applications such as high-order harmonic generation and laser-based particle acceleration.

In addition, RSSHG can be used for efficient second harmonic generation. A conversion efficiency of 23% was experimentally demonstrated, which is in good agreement with the results of PIC simulations. This high SHG conversion efficiency offers possibilities for many applications, including the synthesis of single-cycle optical pulses at an energy level that has not been accessible thus far.

Materials and methods

Experimental setup

The scheme of the experimental setup is shown in Fig. 1. The relativistic surface harmonics are driven by the pulses generated in the OPCPA system, providing a 6–7 fs pulse duration at a 900 nm central wavelength with up to 25 mJ of energy on a target and a contrast of better than 1011 at a rise time of only 1 ps16. The typical (RMS) pulse energy instability was 3%. As demonstrated in previous works36,37,38,39,40, the performance of relativistic surface harmonics is very sensitive to the pre-plasma gradient, which needs to be optimized. To do so, a pre-pulse with an intensity of ~1015 W/cm2 and an adjustable delay (τ) is introduced before the main pulse36 to pre-ionize the target and initiate the plasma expansion with a speed of Cs ≈ 80 nm/ps41, leading to an approximately exponential42 plasma density gradient at the solid–vacuum interface (\(N_{\mathrm{e}} = N_0e^{ - x/L_{\mathrm{p}}}\)) with a scale length of Lp = Cs × τ. The radiation reflected from the target (including generated harmonics and the rest of the fundamental harmonic) is re-collimated with an f/1.3 90o off-axis Al-coated parabolic mirror and guided out of the vacuum chamber, where the spectrum is recorded with a set of two spectrometers (NIRQuest-512 and HR2000 + from Ocean Optics) covering the entire fundamental and second harmonic spectral range.

All parts of the experimental setup, starting from the compressor of the pump pulses and the OPCPA stages, were in vacuum to prevent the degradation of the beam quality due to nonlinear distortions. The material thickness during the dispersion scan was changed before the seed pulse was coupled into the vacuum setup and therefore before the amplification chain. This resulted in a change in the delay between pump and seed pulses in the OPCPA stages, which needed to be re-optimized for each material thickness to ensure identical parameters of the amplified pulses. This was the main factor limiting the number of steps in the dispersion scan during one experimental campaign. Note that the above-mentioned problem is specific to our current beamline design and not a general one. Therefore, it can be solved in our setup after the next upgrade, and it can be avoided in other systems. In particular, wedges for the variation of the material thickness can be installed after the OPCPA (but preferably before the compressor).

During the dispersion scan, each spectrum was measured by averaging over 10 laser shots to diminish the influence of the pulse energy instability.

Effective vector potential evaluation

The laser parameters listed above correspond to a peak intensity of \(I = 4 \times 10^{19}\;{\mathrm{W}}/{\mathrm{cm}}^2\) resulting in the maximum normalized vector potential of \(a_0\,=\,4.8\), where \(a_0 = eE_0/m_{\mathrm{e}}\omega c = \sqrt {I\lambda ^2/(1.37 \times 10^{18})}\), \(\lambda\) is the central wavelength in µm, \(I\) is the intensity in W/cm2, e and \(m_{\mathrm{e}}\) are the electron charge and mass, respectively, and E0 and ω are the laser field amplitude and angular frequency, respectively. However, to compare the experimental results with the 1D PIC simulations, an effective normalized value, in our case, averaged over the full width at half maximum diameter of the focal intensity distribution, should be used. The maximum effective vector potential in the performed experiments was thus \(\approx 4.8/\sqrt 2 = 3.4\), the second power of which limits, for example, the maximum value of \(I_{\omega \_{\mathrm{in}}}\) in Fig. 2a.

PIC simulations

1D PIC simulations were performed with the LPIC + +24 code. In the simulations, the plasma density profile has an exponential layer in front of a constant slab density layer. The density of the interface layer rises from 0.1nc up to a maximum of 400nc, corresponding to the density of a glass target when fully ionized, where nc is the critical electron density. The p-polarized laser pulse is incident onto the plasma at an angle of 45°. This oblique incidence geometry is transformed into a 1D case using the Bourdier technique43.

D-scan retrieval

The retrieval of the D-scan results presented in this paper was carried out with a home-written program based on the differential evolution algorithm and a reconstruction procedure described in detail elsewhere26.

The error of the D-scan retrieval was determined using the standard equation26:

$$G = \sqrt {\frac{1}{{NM}}\mathop {\sum }\limits_{i = 1}^N \;\mathop {\sum }\limits_{j = 1}^M \;I_{2\omega }^{({\mathrm{meas}})}\;\left( {\omega _i,L_j} \right) - R_iI_{2\omega }^{({\mathrm{retr}})}(\omega _i,L_j)}$$

where \(I_{2\omega }^{({\mathrm{meas}})}\) and \(I_{2\omega }^{({\mathrm{retr}})}\) are the measured/simulated and retrieved SHG signals, respectively, which are normalized to the unity maximum value, N and M are the resolution (length of the data list) in the spectral and dispersion domains, respectively, and R accounts for (in the case of measured data) the spectrometer response and SHG conversion efficiency:

$$R_i = \frac{{\mathop {\sum }\nolimits_{j = 1}^M \;I_{2\omega }^{({\mathrm{meas}})}(\omega _i,L_j)I_{2\omega }^{({\mathrm{retr}})}(\omega _i,L_j)}}{{\mathop {\sum }\nolimits_{j = 1}^M \;I_{2\omega }^{({\mathrm{meas}})}\;\left( {\omega _i,L_j} \right)^2}}$$