Introduction

Recent years have seen a tremendous surge in development of laser-based, high-field terahertz (THz) sources for a wide range of applications, including THz-based accelerators1,2,3, control and metrology of electrons4, femtosecond X-ray pulse characterisation5,6, control of dynamics in solids7,8 and spectroscopy9. Most work has concentrated on generation of single-cycle, broadband THz pulses, resulting in a mature technology yielding close to millijoule energies and optical-to-THz conversion efficiencies near 1%10,11,12,13. Some emerging applications, however, such as waveguide-based electron acceleration14, resonant pumping of materials15 and narrowband spectroscopy9, require high-field, narrow-bandwidth THz pulses for which the technology is much less developed. In particular, relativistic acceleration of electrons calls for optically synchronised multi-millijoule pulses with MV/cm electric-field strengths, frequencies below 1 THz and bandwidths below 1%16,17, for which no technology currently exists.

Early on, periodically-poled lithium niobate (PPLN) was identified as promising for optical difference-frequency generation (DFG) of multicycle THz pulses18,19 due to its high nonlinear coefficient, d33 ~25 pm/V, and predicted conversion efficiencies in the range of several percent17. Such crystals would enable generation of THz pulses with tens-of-millijoule energies, sub-percent bandwidths and MV/cm fields using the Joule-level optical pulses available today. However, the insufficient optimisation of the narrowband optical-to-THz conversion, and the resultant low fraction of optical photons that contribute to this process, is a primary obstacle to achieving THz pulses with the above parameters. Recent work has therefore concentrated on addressing the systematic optimisation of laser and crystal parameters with a view towards making every photon count.

Initial experiments using compressed Ti:Sapphire pulses demonstrated that the THz yield can be significantly increased (×5) by optimising the optical bandwidth of the driver laser and by cryogenic cooling of the PPLN which reduces THz reabsorption20. Using short-pulse drivers, however, results in high intensities that limit the energy that can be applied to the crystal before the conversion efficiency saturates or the crystal damages. In addition, the simultaneous presence of photons spanning a broad bandwidth enables phase-matched generation of unwanted, higher-frequency THz20, which not only reduces efficiencies by robbing energy from the main process, but also contributes to damage and thus reduces the total amount of useful pump energy that can be applied to the crystal.

These issues, however, can be overcome by pumping a nonlinear crystal with temporally chirped pulses. Chirping the drive pulse reduces the peak intensity, which allows for large increases in optical energy, and also reduces the instantaneous spectral bandwidth by mapping frequency onto time. In order to provide the correct spectral content for DFG, two chirped pulses must then be combined with a delay.

This chirp-and-delay technique, pioneered originally for driving photoconducting antennae21, has since been extensively used in a broad range of scientific settings22,23,24 because it provides a simple means of creating optical pulse trains of tunable number and periodicity which can be effectively used to drive resonant processes. A complication is that chirped pulses from chirped-pulse amplification (CPA) based laser systems inherently contain high-order phase, such as third-order dispersion (TOD), which induces temporal variations in the periodicity of the pulse train (temporal beating) and can reduce effectiveness of the driver23,25. For narrowband THz generation, the TOD-induced pulse-train periodicity variations, which are connected to temporal variations in the instantaneous difference-frequency content, can severely limit the conversion efficiency. Unfortunately, the intuition to completely remove the TOD from the chirped drive pulses leads to approaches which are either practically infeasible or produce undesired changes in the pulse properties.

Here, we present an elegant and easy-to-implement method for precisely controlling the temporal variation of the optical modulation period. Instead of removing the TOD, we simply tune the relative spectral phase between the chirped pulses. The former linearises the chirp, but requires an impractically large change to both pulses, while the latter uses a very small difference in chirp between the pulses to monochromatise the difference-frequency spectrum and regularise the pulse train. Our method not only enables significant optimisation for efficiently driving narrowband resonant processes but also offers new possibilities in producing designer pulse trains with varying periodicities optimised for specific purposes. We apply this technique to the problem of optimising the conversion efficiency for narrowband, multicycle THz generation by regularising the driver pulse train. By tuning the relative spectral phase in a chirp-and-delay setup, we achieve a 13-fold relative enhancement in conversion efficiency resulting in sub-percent bandwidth, 361 GHz pulses with a conversion efficiency of 0.24%. Using a large-aperture PPLN, we demonstrate a record energy of 0.6 mJ which exceeds previously reported energies by over an order of magnitude26,27.

Results

Monochromatising the difference-frequency spectrum

The idealised case for narrowband difference-frequency generation (DFG) calls for two narrow-line optical sources with a well-defined frequency difference, Δω. Unfortunately, these sources are not yet available with appreciable energies. Chirp-and-delay offers a way of obtaining, effectively, a two-line narrowband source from a broad-band, linearly chirped laser by combining two pulses with a delay Δt. The resulting instantaneous spectral features are narrow and have a constant frequency difference Δω. For this reason, chirp-and-delay is also known as spectral focusing28. By adjusting the delay, Δω can be tuned to match the design frequency ΩTHz of a PPLN crystal, enabling effective conversion of the provided difference-frequency content into a THz pulse (Fig. 1a).

Fig. 1
figure 1

Chirp-and-delay concept. a Two broadband but linearly chirped and delayed pulses provide narrow spectral features with constant instantaneous difference frequency Δω, which is converted in a PPLN crystal into a THz pulse of frequency ΩTHz = Δω. b Due to higher order dispersion Δω varies along the pulse and limits the range where the provided difference frequency fulfils the phase matching condition in the PPLN. c Slightly tilting one of the pulses in time-frequency space regularises Δω and thereby maximises THz generation

A complication of this concept is that chirped pulses from CPA-based laser systems include not only group delay dispersion (GDD), i.e., linear chirp, but also non-negligible third-order dispersion (TOD) which adds a quadratic component, i.e., curvature, to the chirp (Fig. 1b). The result is a linear variation of the difference frequency with time (i.e., difference-frequency chirp) which detunes a fraction of the drive pulse energy away from the resonance condition (e.g., phase-matching) of the driven process. This effect becomes increasingly severe as the bandwidth of the resonance narrows, such as for difference frequency generation of narrowband THz pulses with a large number of cycles.

An intuitive solution would be to completely remove the TOD, thereby linearising the chirp. Unfortunately, this solution is generally impractical because most dispersive systems couple TOD with large amounts of GDD and thereby also produce undesired changes to the pulse duration. Although chirped mirrors can be designed to provide pure TOD, they typically compensate only of order 104 fs3 per bounce, while CPA systems usually have at least 106 fs3 of TOD. It is, however, not necessary to completely remove the TOD. What is actually important is to monochromatise the difference frequency, Δω, which can be achieved much more elegantly by simply adding a small amount of spectral phase to one of the two pulses (see Methods). The effect is to slightly tilt one pulse in time-frequency space (Fig. 1c) so that the temporal dependence of Δω is eliminated despite the chirp curvature. In this way, Δω remains within the narrow acceptance bandwidth of the PPLN around ΩTHz over the full duration of the pulse, optimising the efficiency of the process.

Equivalently, this optimisation of the THz conversion efficiency can be described in the time domain: when the periodicity of the optical intensity modulation matches the THz carrier period, the THz generated by successive pulses in the train add coherently. Irregularities in the pulse-train periodicity induce destructive interferences which become more severe as the number of THz cycles increases. Monochromatising the difference frequency Δω, which determines the temporal beat frequency, ensures a regular pulse train and coherent addition of THz waves.

Millijoule-scale narrowband THz generation

To demonstrate this concept experimentally, we use a modified Mach-Zehnder interferometer (shown in Fig. 2a) in which uncompressed pulses from a Ti:Sapphire system are split into two identical copies and recombined after spectral phase manipulation.

Fig. 2
figure 2

Narrowband THz setup. a The driver laser output is split into two identical copies, whose GDD and TOD are independently tuned using a variable amount of high-dispersive SF11 glass (adding GDD), and a set of dispersion compensating mirrors (DCM), each adding a discrete amount of TOD. The recombined pulses generate narrowband THz pulses in large-aperture cryogenically cooled magnesium-doped PPLN crystals, shown in the photograph (b). To characterise the THz pulses, a pyroelectric detector is used either directly (pulse energy) or as the detector in a longitudinal interferometer (frequency). BS: beam splitter; FF: far field camera; HWP: half-wave plate

Typical CPA systems have a positive chirp curvature, so monochromatising Δω requires adding positive GDD and negative TOD (see Methods). However, negative TOD is generally not available. While it is conceptually simpler to apply dispersion compensation to only one pulse, as illustrated in Fig. 1, an equivalent solution is to add GDD to one pulse and TOD with reversed sign to the other. We therefore split the dispersion compensation between the two interferometer arms, which allows us to use custom-designed dispersion compensating mirrors (DCMs), available with positive TOD, to tune the second pulse. This design also compactifies the setup. A tunable amount of GDD is added by inserting a variable thickness of high-dispersion glass (SF11) and TOD is added incrementally by replacing high-reflectors with DCMs.

In order to maximise the pump energy that we can apply to the DFG process, we employ large-aperture (shown in Fig. 2b), magnesium-doped PPLN (LA-PPMgLN) crystals which were custom fabricated by the Institute for Molecular Science in Japan29. Since Mach-Zehnder interferometers necessarily provide two identical outputs, we send each onto independent, cryogenically cooled crystals. This configuration provides a natural way to effectively double the crystal aperture, further increasing the allowable incident energy and making use of energy that would otherwise be wasted. Total optical pulse energies in excess of 1 Joule are thus safely applied to the crystals. The resulting THz pulse energy is then measured with a pyroelectric detector and the spectrum is characterised using a longitudinal interferometer.

Here we investigate the performance of crystals with a poling period of 330 µm and a corresponding emission frequency of 361 GHz, which is ideally suited for powering a THz-based accelerator16. Details of the performance of a PPLN crystal designed for a higher frequency of 558 GHz (i.e., with a periodicity of 212 µm), are presented in Supplementary Note 1. For our driver laser, which has a GDD of 2.05 × 106 fs2 and a TOD of −4 × 106 fs3, we calculate (see Methods) an optimum compensation of ΔGDD = 9100 fs2 and ΔTOD = 54,000 fs3 for the 330 µm crystal. Note that these phase compensation amounts are much smaller than the original spectral phase (i.e., ΔGDD GDD and ΔTOD TOD) which ensures the pulse duration is not significantly affected and is a key factor in the practicality of our technique.

Figure 3a shows the THz signal as a function of the delay between the driver pulses. In the uncompensated case (red dots), THz is generated with low signal levels over a broad range of delays. Here, each delay providing THz signal corresponds to matching of the instantaneous frequency difference, Δω, with the PPLN design frequency, ΩTHz, at a temporally and spectrally distinct part of the pulse. Adding GDD by increasing the amount of SF11 in the beam enhances the THz signal (blue triangles) by allowing more parts of the pulse to contribute to THz generation simultaneously. As a result, the range of THz-generating delays narrows until most of the pulse is engaged at a single delay. The experimentally determined optimum (green squares) of ΔGDD = 9200 fs2 (or 49 mm SF11) is in very good agreement with our calculation and results in an 11-fold increase in THz output.

Fig. 3
figure 3

Optimisation of the relative spectral phase. a Scan of pulse delay to match Δω to the PPLN frequency ΩTHz. Adding ΔGDD, by inserting SF11 into one of the interferometer arms, the THz signal increases up to a factor of 11 (green squares), and the range of THz-generating delays narrows. b Conversion efficiency enhancement vs ΔGDD for varying amounts of ΔTOD. For an optimum ΔTOD, the signal increases by another 20% relative to no ΔTOD correction, yielding a total increase in THz output by a factor of 13

Using DCMs to tune the relative TOD between the two driver pulses, the THz signal is further enhanced by 20% yielding in total a 13-fold increase in output compared to the initial configuration (Fig. 3b). The measured optimum of ΔTOD = 20,000 fs3 is different from the predicted optimum, which may be attributable to the presence of fourth and higher orders in the pump spectral phase. In contrast to the GDD correction, the discrete nature of the TOD compensation makes fine-tuning of the TOD difficult.

The above optimisation of the THz yield can be understood as an optimisation of the temporal range over which the time-dependent frequency difference Δω(t) falls within the THz bandwidth ΩTHz ± ΔΩTHz accepted by the phase-matching process (Fig. 4a). The portion of the optical pulse contributing to THz generation is then simply estimated as the fraction satisfying | Δω(t)ΩTHz | ≤ ΔΩTHz. More precisely, this fraction is the integrated product of the phase-matching efficiency curve in the time domain and the temporal intensity profile of the optical pulse (Fig. 4a). This quantity serves as a good measure predicting the variation in THz output with delay and with GDD and TOD corrections, such as the enhancement of THz at positive delays and the reduction of THz at negative delays (Fig. 4b).

Fig. 4
figure 4

Asymmetric chirp compensation. a The 330 µm poling period PPLN has a central frequency of ΩTHz = 361 GHz with a 1% bandwidth. Initially, the difference frequency Δω varies linearly in time (red). As spectral phase compensation is applied adding ΔGDD (dashed blue) and further adding ΔTOD (solid green) the fraction of the pulse falling within the THz bandwidth ΩTHz ± ΔΩTHz, and thus the generated THz energy, increases. The top inset shows the PPLN bandwidth multiplied by the difference frequency. The integrated product of this phase-matching efficiency with the pulse envelope (grey) is a measure for the fraction of the pulse contributing to THz generation. b Adding ΔGDD (green) creates an asymmetry between the two pulses in time that reduces the THz signal at negative delays and enhances the THz at positive delays. For identical pulses (red) the THz signal is equal for positive and negative delays. A calculation (solid lines) based on the concept discussed in the main text and illustrated in panel a predicts the THz output in good agreement with the measurements (squares, dots)

To characterise the performance of the crystal, we use two definitions of conversion efficiency which highlight different quantities. The “extracted” efficiency, plotted in Fig. 5a, represents the ratio of useful THz energy extracted from the crystal to the driver incident energy. The peak extracted efficiency of 0.15% is achieved at a fluence of 180 mJ/cm2. This figure corresponds to an “internal” efficiency of 0.24%, which represents the ratio of THz and optical energies within the crystal. The internal efficiency serves as a useful metric for the optimisation of the DFG process independent from the practical complexities of input and output coupling. The fine-tuning of the TOD not only increases the conversion efficiency, but also increases the fluence at which optimum conversion efficiency occurs, leading to large increases in the total THz yield.

Fig. 5
figure 5

High THz energy. a Fine tuning the spectral phase increases the fluence at which the maximum conversion efficiency occurs and reaches 0.15% at 180 mJ/cm2 for crystal #1. Note, that with simultaneous GDD and TOD correction, the optimum determined ΔGDD is slightly different than the value reported in Fig. 3a. b Using GDD and TOD compensation, an input fluence of 350 mJ/cm2 generates individual THz pulse energies up to 458 μJ (392 µJ average) and 232 μJ (212 μJ average) in crystal #1 and #2, respectively. The difference in THz output is attributed to differences in the preparation of the crystal surfaces (see Methods). Significant red-shift of the driver spectrum is observed for high-energy THz generation indicating a high degree of optical cascading. c Interferometric correlation of the generated THz pulses and its frequency spectrum (d), confirming the central frequency of 361 GHz and a sub-1% bandwidth. The interferometer (inset) replaces the energy measurement setup, Fig. 2a. Note in c that each data point is averaged over 10 shots. The sampling of the frequency, however, accentuates the visual effect of the shot-to-shot fluctuations, while indeed the centre frequency is captured with high precision (see Supplementary Note 2). Error bars correspond to the standard deviation of 20 measurements in the case of a and b

Figure 5b shows the corresponding extracted THz energies. The highest individual THz pulse energy of 458 μJ and the highest average 2-crystal output of 0.6 mJ are achieved at a fluence of 350 mJ/cm2. The uniform red-shift of the optical spectrum (inset) shows that pump depletion has occurred across the entire spectrum, which confirms that the full optical pulse engages in the DFG process. The spectrum of the generated radiation is measured using an interferometric autocorrelation of the THz pulse (Fig. 5c). The Fourier transform of this trace shows a central frequency of 361 GHz and confirms a sub-percent bandwidth (Fig. 5d). If the pulses were to be combined, the 0.6 mJ and roughly 200 ps pulses could produce a field strength of ~18 MV/m when focused to a 1/e2 waist of 2 mm (2.4 × λTHz).

Note that the pump fluence for peak yield is well above the fluence for peak conversion efficiency, which suggests benefits from even larger aperture crystals. The difference in performance of the two crystals (Fig. 5b) is partly due to anti-reflection measures, applied to only one crystal, that reduce Fresnel losses and improve the extraction efficiency by 35% (see Methods).

Our measurements show that the need for spectral-phase tuning is driven by the narrow acceptance bandwidth of the process. For a similar degree of tuning, we therefore expect a lower conversion efficiency for a crystal with a narrower phase-matching bandwidth. In fact, this effect is confirmed by the lower efficiency measured with the 212 µm periodicity crystal which has a larger number of poling cycles and therefore a narrower acceptance bandwidth. Accounting for variations in absorption coefficient, crystal length and poling period, we predict a conversion efficiency of 0.08% at 558 GHz which is in very good agreement with the measured value of 0.07% (see Supplementary Note 1).

Despite the large efficiency enhancements obtained by the relative spectral phase tuning, the maximum absolute conversion efficiency reported here is only twice larger than that reported by us previously27. An exact comparison between these results is not possible, however, since nonlinear crystals of different aperture, length and vendor were used.

In general, a proper comparison of efficiencies across different experiments requires a careful evaluation of the effects of differences in the laser, THz and crystal parameters, which can be very complex (see Supplementary Note 3). For example, the DFG efficiency intrinsically scales quadratically with THz frequency, due to higher THz photon energies, so one can expect higher efficiencies for higher frequencies. On the other hand, THz absorption also increases with frequency, mitigating, to varying degrees, the photon-energy effect. Differences in the crystal length, poling period, aperture and crystal manufacturer all also strongly affect the conversion process. The effect of crystal length on efficiency is especially complex, since length governs not only the THz bandwidth and therefore the sensitivity of the conversion process to TOD, but also the growth of cascading17 and amounts of pump depletion and THz reabsorption. Differences in the crystal aperture can lead to significant variations in efficiency via inhomogeneities in the poling structure29. The effects of crystal material and poling quality on multicycle THz generation are so far nearly unexplored. A detailed characterisation of these effects is thus an important part of the further optimisation of the multicycle THz generation process.

What we unambiguously demonstrate here is the role of TOD in detuning narrowband THz generation using the chirp-and-delay setup. The degree of relative enhancement achievable by tuning the relative spectral phase is robustly predictable from the amount of TOD and the phase-matching bandwidth, independent of other factors which may affect the efficiency. Mitigation of the effects of TOD is thus a required ingredient for optimising narrowband THz generation using chirp-and-delay.

Discussion

We demonstrate, for the first time, the sensitivity of narrowband THz generation to TOD in the chirp-and-delay setup and show that this sensitivity can be attributed to temporal variations in the difference-frequency content which cause a very strict phase-matching condition not to be satisfied for the majority of the optical pulse. The effect of TOD can be understood in the time domain as inducing pulse-train irregularities which cause a destructive superposition of THz waves generated by the various pulses within the train. To address this issue, we introduce a powerful new technique based on relative spectral-phase tuning of the chirped pulses. The technique provides, for the first time, a simple and easily-implemented method to control and fine tune the temporal dependence of the difference-frequency spectrum, and hence the pulse-train periodicity, without changes to the laser architecture or to the pulse-train intensity and duration.

For high-energy narrowband THz generation, where the priority is optimising the optical-to-THz conversion efficiency, we adjust the spectral phase to monochromatise the difference-frequency spectrum and to create a regularised pulse-train (see Supplementary Notes 4 and 5), and thus obtain a 13-fold enhancement in THz yield (at 361 GHz). By using a Joule-class laser as the chirped-pulse source to pump two custom, large-aperture PPLN crystals (LA-PPMgLN29), we achieve a record energy of 0.6 mJ which exceeds previous results in optical generation of multicycle THz by over an order of magnitude. The internal conversion efficiency we measure of 0.24% is also nearly twice larger than previously reported for narrowband sources20,26,27. This result is especially significant in light of the phase-matching difficulties associated with sub-percent THz bandwidths and the additional complexities associated with scaling to centimetre apertures and Joule-level pump energies29.

Millijoule-scale narrowband THz pulses represent a break-through which demonstrates the feasibility of the emerging technology of THz-based acceleration that holds high promise for revolutionising electron accelerators and the light sources based on them16. It is expected that 0.6 mJ of THz in a 200-ps-long pulse as produced in this work could create a field strength of 50 MV/m when coupled in to a waveguide designed for matched acceleration. With acceleration over the roughly 4 cm allowed by such a pulse, sub-relativistic electrons could be accelerated beyond 2 MeV, an energy already relevant for ultrafast electron diffraction. This can be compared to past acceleration in a waveguide with broadband THz, where the acceleration length was limited to 3 mm and the energy gain was limited to 7 keV1.

However, our results are equally applicable for enhancing the yield of narrowband THz pulses at any energy level. For the average researcher with a mJ- or sub mJ-scale laser system, making every photon count is an imperative. Order-of-magnitude-scale efficiency enhancements are thus a game changer for enabling experiments requiring high-field, narrowband THz pulses, such as pumping of Josephson plasmon resonances in layered superconductors30,31,32. Spectroscopic applications may be accessible with our source parameters9, and we have demonstrated the technique at multiple frequencies. However, the lack of demonstrated tunability of the PPLN technology may limit the practical application to spectroscopic scenarios that do not require a large number of frequencies.

As chirp-and-delay is ubiquitously used for generating pulse trains, tuning the pulse-train is also relevant for a wide range of applications that require resonant excitation, including excitation of plasma waves33, high-order harmonic generation22, laser modulation of electron beams23, coherent anti-Stokes Raman scattering microscopy of biological structures24, excitation of atoms34, molecules and solids and driving photocathode injectors for conventional accelerators35. These applications have so far been limited by the presence of high-order phase in the chirped-pulse drivers and have suffered from the lack of practical possibilities to remove it. Our technique addresses this need. In addition, the technique opens up the possibility of generating pulse trains with customised temporal periodicity profiles with novel applications such as generation of pre-chirped multicycle THz pulses that compress under the action of waveguide dispersion, and thus address a difficult problem in the delivery of high-field THz to an interaction point36.

Methods

Laser system

We use a 260-ps FWHM uncompressed, high-power CPA-based Ti:Sapphire laser providing 1.2 J at 5 Hz repetition rate. Inherent to the Öffner stretcher design, the laser pulses have a GDD of ϕ2 = 2.05 × 106 fs2 and a TOD of ϕ3 = −4 × 106 fs3 with additional uncharacterised fourth and higher order spectral phase. The beam has a super-Gaussian profile of 13 mm 1/e2 intensity diameter, and a spectral width of more than 30 nm FWHM centred at 800 nm.

Dispersion management

Using the common description of spectral phase, the time-dependent frequency of a chirped pulse can be written as ω(t, ϕ2, ϕ3) = ω0 + 1/ϕ2 × t − ϕ3/2ϕ23 × t2, with ω0 the central frequency, and ϕ2 and ϕ3 the pulse GDD and TOD. Here, we neglect higher orders and assume the GDD is large. The parameters ϕ2 and ϕ3 are properties of the pulse and are typically inherent to the laser design. When there is no TOD (i.e., ϕ3 = 0), the chirp is linear and the instantaneous difference in frequency Δω is constant, determined only by the chirp rate (i.e., GDD, ϕ2) and the delay:

$${\mathrm{\Delta }}\omega \left( t \right) = \omega \left( {t,\phi _2} \right) - \omega \left( {t - {\mathrm{\Delta }}t,\phi _2} \right) = \frac{{{\mathrm{\Delta }}t}}{{\phi _2}}.$$
(1)

In the presence of TOD, however, the chirp becomes nonlinear (i.e., curved), and Δω becomes a linear function of time with a slope that is proportional to the TOD:

$${\mathrm{\Delta }}{\it{\omega }}\left( t \right) = \omega \left( {t,\phi _2,\phi _3} \right) - \omega \left( {t - {\mathrm{\Delta }}t,\phi _2,\phi _3} \right) = \left[ {\frac{{{\mathrm{\Delta }}t}}{{\phi _2}} + \frac{{\phi _3}}{{2\phi _2^3}}\Delta t^2} \right] - \left[ {\frac{{\phi _3}}{{\phi _2^3}}{\mathrm{\Delta t}}} \right] \times t.$$
(2)

Applying the condition Δω = ΩTHz shows that phase matching is only achieved at a particular time, t*, which is dependent on the delay Δt and given by:

$$t^ \ast \left( {{\mathrm{\Delta }}t} \right) = \frac{{\phi _2^2}}{{\phi _3}} + \frac{1}{2}{\mathrm{\Delta }}t - \frac{{\phi _2^3}}{{\phi _3{\mathrm{\Delta }}t}}{\mathrm{\Omega }}_{THz}.$$
(3)

In the case of two identical chirped pulses, perfect phase matching requires elimination of the time-dependence of Δω, and thus can only be achieved by complete elimination of the TOD. Adjusting the relative dispersion of the two pulses, however, provides additional options. Adding a small amount of GDD, Δϕ2, and TOD, Δϕ3, reduces the curvature of the time-frequency mapping, which is equivalent to slightly tilting the pulse in time-frequency space, shown in Fig. 1c. The instantaneous difference in frequency between two chirped pulses, where one has been delayed by Δt and tuned in spectral phase by adding Δϕ2 and Δϕ3, then reads

$${\mathrm{\Delta }}\omega \left( t \right) = \omega \left( {t,\phi _2,\phi _3} \right) - \omega (t - {\mathrm{\Delta }}t,\phi _2 + {\mathrm{\Delta }}\phi _2,\phi _3 + {\mathrm{\Delta }}\phi _3).$$
(4)

It is straightforward to show that Δω(t) is now a second order equation in t, expressible as Δω(t) = c0 + c1×t + c2×t2 with ci = ci2, ϕ3, Δϕ2, Δϕ3, Δt). More explicitly stated,

$$\begin{array}{*{20}{l}} {{\mathrm{\Delta }}\omega \left( t \right)} \hfill & = \hfill & {\left[ {\frac{1}{{\phi _2 + {\mathrm{\Delta }}\phi _2}}{\mathrm{\Delta }}t + \frac{{\phi _3 + {\mathrm{\Delta }}\phi _3}}{{2\left( {\phi _2 + {\mathrm{\Delta }}\phi _2} \right)}}{\mathrm{\Delta }}t^2} \right]} \hfill \\ {} \hfill & {} \hfill & { + \left[ {\frac{1}{{\phi _2}} - \frac{1}{{\phi _2 + {\mathrm{\Delta }}\phi _2}} - \frac{{\phi _3 + {\mathrm{\Delta }}\phi _3}}{{\left( {\phi _2 + {\mathrm{\Delta }}\phi _2} \right)^3}}{\mathrm{\Delta }}t} \right]} \hfill \\ {} \hfill & {} \hfill & { - \left[ {\frac{{\phi _3}}{{2\phi _2^3}} - \frac{{\phi _3 + {\mathrm{\Delta }}\phi _3}}{{2\left( {\phi _2 + {\mathrm{\Delta }}\phi _2} \right)^3}}} \right] \times t^2}. \hfill \end{array} \times t$$
(5)

According to this expression, a constant instantaneous frequency difference Δω = ΩTHz, as illustrated in Fig. 1c, can be recovered by setting the coefficients c1 and c2 to zero, and solving for Δϕ2, Δϕ3, and Δt

$$\begin{array}{*{20}{l}} {{\mathrm{\Delta }}\phi _2} \hfill & = \hfill & {\phi _2\left( {\frac{{\phi _2}}{{\phi _2 - 2\phi _3{\mathrm{\Omega }}_{THz}}}} \right)^{\frac{1}{2}} - \phi _2} \hfill \\ {{\mathrm{\Delta }}\phi _3} \hfill & = \hfill & {\phi _3\left( {\frac{{\phi _2}}{{\phi _2 - 2\phi _3{\mathrm{\Omega }}_{THz}}}} \right)^{\frac{3}{2}} - \phi _3} \hfill \\ {{\mathrm{\Delta }}t} \hfill & = \hfill & {\frac{{2\phi _2{\mathrm{\Omega }}_{THz}}}{{1 + \left( {1 - \frac{{2\phi _3{\mathrm{\Omega }}_{THz}}}{{\phi _2}}} \right)^{\frac{1}{2}}}}}. \hfill \end{array}$$
(6)

THz generation

The chirp-and- delay concept is implemented using a Mach-Zehnder interferometer, with the GDD and TOD compensation distributed between the two interferometer arms. SF11 adds GDD of Δϕ2/L = +187.5 fs2/mm and additional small TOD of Δϕ3/L = +126 fs3/mm. Custom dielectric mirrors each add Δϕ3 = +20,000 fs3 and no additional GDD, independently confirmed with a white-light interferometer. For experiments at high fluence, two telescopes locally increase the beam size on the DCMs. The dispersion through the lenses is about 800 fs2 and compensated with additional 4.4 mm SF11, which is not explicitly mentioned in the main text. The setup is very sensitive to the angular overlap of the two beams, which we ensure to sub-100 µrad level with a far field camera after the second beam splitter. The PPLN crystals are cooled below 100 K in a cryostat operating at 8 × 10−4 mbar. The entrance window is AR coated for the NIR driver pulses. A ceramic aperture in front of the crystals (8.5 × 14.5 mm crystal #1, 9.5 × 14.5 mm crystal #2) prevents the beam from clipping on the crystal edges. We use Teflon plates with calibrated transmission, to separate the driver IR beam from the generated THz pulses. The THz energy is detected with a pyroelectric sensor (Gentec-EO: THZ9B-BL-BNC) that had been calibrated to a device calibrated by the Physikalisch-Technische Bundesanstalt (PTB). Each data point reported in the figures is averaged over 5 pulses (long delay scan, Fig. 4b), 100 pulses (GDD and TOD tuning, Fig. 3), and 20 pulses (fluence scan, Fig. 5). When error bars are not shown the error is below 5% rms and not visible within the given figure.

PPLN crystals

We use custom 10 × 15 mm2 large aperture crystals29 with a length of 36 mm and poling periods of 330 µm and 212 µm, corresponding to THz pulses of 361 GHz and 558 GHz, respectively. For the generation of high-energy THz pulses reported in Fig. 5, the 361 GHz crystal #1 is AR coated for the driver laser and features a 100-µm fused silica wafer at the end facet of the crystal to minimise Fresnel losses and maximise THz extraction. For comparison, the 361 GHz crystal #2 is left uncoated and without a wafer installed. The AR coating and fused silica wafer together enhance the THz output of crystal #1 by about 35% compared to crystal #2. The remaining 25% difference in observed THz output is attributed to variations in the crystal manufacturing.