## Abstract

Can an optical pulse traverse a non-dispersive material at the speed of light in vacuum? Because traditional approaches for controlling the group velocity of light manipulate either the material or structural resonances, an absence of dispersion altogether appears to exclude such a prospect. Here we demonstrate theoretically and experimentally that “space–time” wave packets—pulsed beams in which the spatial and temporal degrees of freedom are tightly intertwined—can indeed traverse a non-dispersive transparent optical material at the speed of light in vacuum. We synthesize wave packets whose spatio-temporal spectra lie along the intersection of the material’s light-cone with a spectral hyperplane tilted to coincide with the vacuum light-line. By measuring the group delay interferometrically with respect to a generic reference pulse, we confirm that the wave packet group velocity in a variety of materials (including water, glass, and sapphire) is the speed of light in vacuum.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

## 1. INTRODUCTION

It has long been known that the group velocity of an optical pulse can vary substantially when traversing an optical material if the pulse spectrum overlaps with a resonance [1]. To date, systems that realize so-called slow light and fast light [2] exercise control over the group velocity through the effects of a carefully crafted wavelength-dependent refractive index in diverse physical embodiments, via material [3–5] or structural resonances [6,7]. Such control facilitates optical buffering, enables the implementation of true optical delays, and offers advantages with respect to interferometric precision [2]. However, in transparent materials at wavelengths far from resonance, the wavelength dependence of the refractive index may be ignored, especially for narrow bandwidths, in which case the group velocity is ${v}_{\mathrm{g}}\approx c/n$, where $c$ is the speed of light in vacuum, and $n$ is the refractive index. Under these assumptions, can a pulse travel in such a medium at $c$? In absence of dispersion, the large required change in group velocity from $c/n$ to $c$ does not seem possible.

Here, we demonstrate that the key to versatile control over the group velocity of a pulsed beam traveling in an optical material is to manipulate the spatial and temporal degrees of freedom *jointly*. Propagation-invariant space–time (ST) wave packets [8,9] endowed with tight spatio-temporal spectral correlations [10–14] can indeed achieve this apparently impossible goal: to travel diffraction-free and dispersion-free for extended distances at a group velocity of $c$ in linear isotropic non-dispersive materials (and also in dispersive media [15,16]). Because we do not rely on an optical resonance, this behavior is achieved *independently* of the refractive index or the spectral band, *without* amplification or attenuation, and potentially over large bandwidths. Indeed, spatio-temporal structuring produces arbitrary subluminal *or* superluminal shifts that can readily change the group velocity from $c/n$ to $c$—or to any other velocity desired in the material [17–20]. The question concerning the various definitions of group velocities has recently received careful consideration in Ref. [21]. Because we rely on measurements that are resolved in space and time, the group velocity we refer to henceforth is the average group velocity in the direction of the optical axis [21].

Prior efforts for spatio-temporal pulse shaping have yielded only limited control over the speed of light. In free space, typical velocities realized to date have been $1.00022c$ [22], $1.00012c$ [23], and $1.00015c$ [24] (but larger changes can be achieved by the different methodology of “tilted-pulse fronts” [25–27]). This falls short of the required variation needed to change the group velocity in a material from $c/n$ to $c$. Previous approaches to endow optical fields with spatio-temporal correlations have utilized techniques for producing Bessel beams (e.g., annular apertures [28] or axicons [22,23,29]), nonlinear processes such as second-harmonic generation [30] and laser filamentation [31,32], or spatio-temporal spectral filtering [33,34]. In contrast, we make use here of a spatio-temporal *synthesis* approach that relies on a phase-only spectral modulation scheme implemented by a spatial light modulator (SLM) [35,36]. We have exploited this energy-efficient technique in synthesizing ST wave packets that propagate in free space with controllable group velocity [37]. The precision achievable with this generalized pulse-shaping approach has been verified by the high quality of the fit between measurements and theory in a variety of experiments, including the observation of self-healing after opaque obstructions [38], propagation for extended distances in free space [39], and synthesis of broadband ST pulses [40]. However, to the best of our knowledge, there have been no demonstrations of the propagation of ST wave packets synthesized via this novel approach in optical materials, whereupon the coupling from free space to the material is expected to result in substantial variation in the spatio-temporal spectral correlations initially introduced in free space.

In general, the spatio-temporal spectra of propagation-invariant ST wave packets lie along the intersection of the material light-cone with a tilted spectral hyperplane [35,41] whose spectral tilt angle solely determines its group velocity, as described theoretically by Donnelly and Ziolkowski [10] and by Sheppard [9] (Chapter 18), and experimentally in Ref. [37]. Crossing the interface between two materials at normal incidence results in a change in this tilt angle. By judiciously selecting the spectral tilt angle in free space, we arrange for the spectral hyperplane in the material to coincide with the light-line in vacuum, leading to the group velocity in the material becoming exactly $c$. The requisite spatio-temporal spectral correlations are introduced into a pulsed plane wave via a SLM, and the group delay resulting from traversing an optical material is measured with respect to a generic short reference pulse via an interferometric arrangement. We confirm that the synthesized ST wave packets travel at $c$ in a variety of solid and liquid optical materials having refractive indices in the range from 1.3 to 1.76. Furthermore, we demonstrate smooth and continuous control of the ST wave packet group velocity in these materials below and above $c$.

## 2. THEORY OF SPACE–TIME WAVE PACKETS IN NON-DISPERSIVE MEDIA

We start our theoretical analysis from the decomposition of pulsed beams in a non-dispersive medium of refractive index $n$ into monochromatic plane waves ${e}^{i({k}_{x}x+{k}_{z}z-\omega t)}$, each represented by a point on the surface of the material’s light-cone ${k}_{x}^{2}+{k}_{z}^{2}={n}^{2}{(\frac{\omega}{c})}^{2}$ [Fig. 1(a)]. Here ${k}_{x}$ and ${k}_{z}$ are the components of the wave vector along the transverse $x$- and longitudinal $z$-axes, respectively, and $\omega $ is the temporal frequency. The spatio-temporal spectra of propagation-invariant ST wave packets lie along conic sections at the intersection of the light-cone with the spectral hyperplane $\mathcal{P}(\theta )$ given by the equation $\frac{\omega}{c}={k}_{\mathrm{o}}+({k}_{z}-n{k}_{\mathrm{o}})\mathrm{tan}\text{\hspace{0.17em}}\theta $ [35,41], where the spectral tilt angle $\theta $ is with respect to the ${k}_{z}$-axis, ${\omega}_{\mathrm{o}}$ is an optical carrier frequency, and ${k}_{\mathrm{o}}={\omega}_{\mathrm{o}}/c$. This implies that each temporal frequency $\omega $ is associated with a spatial frequency pair $\pm {k}_{x}$. The field $E(x,z,t)={e}^{i(n{k}_{\mathrm{o}}z-{\omega}_{\mathrm{o}}t)}\psi (x,z,t)$ has a slowly varying envelope whose plane-wave expansion is

The key to synthesizing the associated ST wave packet in *free space*, which produces the targeted wave packet once coupled into the material, is exploiting the invariance of ${k}_{x}$ and $\omega $ across planar interfaces between different media. Given a target plane $\mathcal{P}({\theta}_{\mathrm{m}})$ defining the ST wave packet in the material, the projection of the hyperbola at its intersection with the light cone [Fig. 1(a)] onto the $({k}_{x},\frac{\omega}{c})$-plane is also a hyperbola. This projected spatio-temporal spectrum is invariant across the planar interface between the material and free space [Fig. 1(b)]. Whatever form is taken by the spatio-temporal spectrum on the free-space light-cone, its projection onto the $({k}_{x},\frac{\omega}{c})$-plane *must* be identical with this hyperbola [Fig. 1(c)]. It can be shown (Appendix A) that this projection is consistent with the intersection of the free-space light-cone with a *non-planar* surface that is parallel to the ${k}_{x}$-axis. The projection onto the $({k}_{z},\frac{\omega}{c})$-plane is hence *not* a straight line, thereby indicating dispersive behavior. Nevertheless, for narrow temporal bandwidths $\mathrm{\Delta}\omega \ll {\omega}_{\mathrm{o}}$, the projection onto the $({k}_{z},\frac{\omega}{c})$-plane of the free-space ST wave packet can be approximated to first order as a straight line making a spectral tilt angle ${\theta}_{\mathrm{a}}$ with respect to the ${k}_{z}$-axis (Appendix A); [Fig. 1(c)].

This simplification amounts to approximating the exact hyperbola in the vicinity of ${k}_{x}=0$ with a parabola. Within this approximation, the transition of the ST wave packet from free space to the material changes the spectral tilt angle from ${\theta}_{\mathrm{a}}$ to ${\theta}_{\mathrm{m}}$ (Appendix A) such that

*medium*at the speed of light in vacuum (${v}_{\mathrm{g}}=c$ or $\mathrm{tan}\text{\hspace{0.17em}}{\theta}_{\mathrm{m}}=1$), one must first synthesize a ST wave packet in

*free space*for which Because $\mathrm{tan}\text{\hspace{0.17em}}{\theta}_{\mathrm{a}}>1$ for $n>1$, the ST wave packet in air must be superluminal. For example, a ST wave packet in air with ${\theta}_{\mathrm{a}}\approx 77.1\xb0$ will propagate in BK7 glass ($n\approx 1.51$) at the speed of light in vacuum instead of $c/n$. For a ST wave packet to travel at $c$ in materials with higher indices, it may require first synthesizing a ST wave packet in free space with a negative group velocity (${\theta}_{\mathrm{a}}>90\xb0$). Interestingly, the transition from positive to negative values of the required group velocity when incident from air—to propagate in the material at $c$—occurs at a refractive index equal to the golden mean:

## 3. EXPERIMENT

The experimental synthesis of ST wave packets utilizes a two-dimensional pulse shaper that simultaneously modulates the transverse spatial profile and temporal linewidth of a pulsed plane wave from a Ti:Sapphire femtosecond laser [35]; (Fig. 2). This process endows the field with the requisite spatio-temporal spectral correlations and realizes the prescribed spectral tilt angle ${\theta}_{\mathrm{a}}$. A reflective phase-only SLM assigns the appropriate spatial-frequency pair $\pm {k}_{x}$ to each wavelength $\lambda $ after spatially resolving the spectrum via a diffraction grating (Fig. 2, inset). The retro-reflected wave front returns to the grating that reconstitutes the pulse and produces the ST wave packet (Appendix A). We make use of a temporal bandwidth of $\mathrm{\Delta}\lambda \approx 0.2\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}$ to avoid dispersive effects in the materials investigated. The spatio-temporal spectrum is recorded after a double Fourier transform in space and time [Fig. 3(a)] from which we can confirm that the targeted value of ${\theta}_{\mathrm{a}}$ has been realized [Fig. 3(b)]. A charge coupled device (CCD) camera is scanned axially along $z$ to capture the diffraction-free propagation of the time-averaged intensity $I(x,z)=\int \mathrm{d}tI(x,z,t)$ of the ST wave packet [Fig. 3(c)].

Measurements of the group delay of ST wave packets are carried out interferometrically with respect to a short reference pulse [37] (Fig. 2). We first reconstruct the spatio-temporal intensity profile of the wave packet (pulse width, 9 ps) by superposing it with a shorter plane-wave reference pulse (100 fs). Observing high-visibility spatially resolved fringes indicates that the two wave packets overlap in space and time, whereupon the fringe visibility is proportional to the spatio-temporal intensity distribution [37]. The spatio-temporal intensity $I(x,z=0,\tau )={|E(x,0,\tau )|}^{2}$ is traced by varying the delay $\tau $ in the path of the reference pulse [Fig. 3(d)]. Of course, using this linear interferometric technique, only the amplitude of the spatio-temporal envelope is reconstructed and not the phase, but that suffices for our purpose here of estimating ${v}_{\mathrm{g}}$. Previous measurements of the spatio-temporal profiles of propagation invariant wave packets have made use of spatially resolved ultrafast pulse measurement techniques [23] or self-referenced interferometry [33,35,36].

We measure the group velocity of the ST wave packet in free space by introducing an extra propagation distance $L$ in the common path of the ST wave packet and reference pulse, whose group delays are ${\tau}_{\mathrm{a}}=L/(c\text{\hspace{0.17em}}\mathrm{tan}\text{\hspace{0.17em}}{\theta}_{\mathrm{a}})$ and ${\tau}_{\mathrm{o}}=L/c$, respectively. When ${\tau}_{\mathrm{o}}-{\tau}_{\mathrm{a}}$ is larger than the ST wave packet pulse width, the fringe visibility is lost, but it can be restored by introducing a delay ${\tau}_{\mathrm{d}}={\tau}_{\mathrm{a}}-{\tau}_{\mathrm{o}}$ in the reference path. Measuring ${\tau}_{\mathrm{d}}$ yields ${\tau}_{\mathrm{a}}$, and thus further confirmation of the value of ${\theta}_{\mathrm{a}}$.

## 4. MEASUREMENT RESULTS

We exploit a similar strategy to estimate the group velocity of the ST wave packet in a material. Placing a layer of length $L$ and index $n$ in the common path results in a loss of interference, which is restored by inserting a delay length $\mathrm{\Delta}\ell $ in the reference path. The group delay of the ST wave packet traversing the material ${\tau}_{\mathrm{m}}$ normalized by the group delay ${\tau}_{\mathrm{o}}$ of a traditional pulse traversing a length $L$ of free space is

*only*on the material refractive index, which provides a simple test to confirm that the ST wave packet propagates in the material at the speed of light in vacuum.

An example of the delayed pulse profile after a sample (N-BK7) is plotted in Fig. 3(e) along with the initial pulse. In Fig. 4(a), we report the measured values of the normalized group delay ${\tau}_{\mathrm{m}}/{\tau}_{\mathrm{o}}$ for glass ($n\approx 1.51$ and $L=48\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$) and distilled water ($n\approx 1.33$ and $L=50\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$). Making use of Eq. (3), the required spectral tilt angle for ${v}_{\mathrm{g}}=c$ in these materials is ${\theta}_{\mathrm{a}}=77.14\xb0$ for glass and ${\theta}_{\mathrm{a}}=60.64\xb0$ for distilled water. The results confirm that we have achieved the condition ${v}_{\mathrm{g}}=c$ in both samples at the expected spectral tilt angles ${\theta}_{\mathrm{a}}$, in addition to demonstrating precise control over ${v}_{\mathrm{g}}$ above and below $c$. We also confirm the expected linear increase in the group delay with length of the optical material via measurements carried out on the BK7 glass samples of increasing thickness from $L=12$ to 60 mm, as shown in Fig. 4(b), while holding ${\theta}_{\mathrm{a}}$ at 77.1° to maintain ${v}_{\mathrm{g}}=c$ in the sample. Finally, we confirm the relation in Eq. (6) that applies to any optical material when the group velocity of the ST wave packet traversing is $c$. The measured values of $\frac{\mathrm{\Delta}\ell}{L}$ for all the materials investigated are plotted in Fig. 5, and the theoretical expectation of ${(n-1)}^{2}$ provides an excellent fit.

## 5. DISCUSSION AND CONCLUSIONS

The effect described here does *not* exploit material or structural resonances to modify the propagation of light. Only spatio-temporal correlations that can be readily encoded into the wave packet spectrum are required. Spatio-temporal structuring of light is a powerful tool enabling unprecedented control over optical propagation, effectively rendering the group velocity in a material independent of the refractive index or wavelength range. Careful spectral engineering coupled with the use of high-resolution phase plates in lieu of a SLM [40] yields ST wave packets that can propagate for extended distances [39], thus making them candidates for remote sensing and standoff detection. Moreover, our approach may potentially be extended to accelerating or decelerating wave packets [42].

The versatile control over the group velocity in optical materials enabled by the proposed wave-shaping approach can help solve many critical problems. Controlling the group velocity of optical pulses is essential for efficient frequency conversion [17,25,27,43]. Our versatile approach paves the way to perform group-velocity matching to boost the efficiency of a larger set of processes, such as non-degenerate two-photon absorption, crucial for the detection of infrared radiation and two-photon lasers [44,45].

Moreover, the spatio-temporal spectral correlations that we have studied here in the context of non-dispersive optical materials can be modified to ensure propagation invariance in dispersive materials at arbitrary group velocities. Finally, the bandwidth that can be exploited in synthesizing ST wave packets that travel in optical materials at a prescribed group velocity is—in principle—extended (enabling the utilization of ultrashort pulses). This is in contradistinction to the narrow linewidths associated with material and structural resonances exploited in traditional approaches to slow- and fast-light, which consequently enforce the utilization of only long pulses in such scenarios.

## APPENDIX A

## 1. Space–Time Wave Packets in a Non-Dispersive Material

In a non-dispersive material of index $n$, the light-cone is ${k}_{x}^{2}+{k}_{z}^{2}={n}^{2}{(\frac{\omega}{c})}^{2}$. Here we assume that the field is uniform along $y$, and thus ${k}_{y}$ is negligible. The spatio-temporal spectrum of a propagation-invariant ST wave packet in such a material lies at the intersection of this light-cone with the spectral hyperplane $\mathcal{P}(\theta )$ given by $\frac{\omega}{c}={k}_{\mathrm{o}}+({k}_{z}-n{k}_{\mathrm{o}})\mathrm{tan}\text{\hspace{0.17em}}\theta $. The projection of this intersection onto the $({k}_{x},\frac{\omega}{c})$-plane is a conic section given by

Both $\omega $ and ${k}_{x}$ are invariant across planar interfaces between different materials. Therefore, the spatio-temporal spectrum projected onto the $({k}_{x},\frac{\omega}{c})$-plane in Eq. (A1) remains invariant when going from the material to air. However, such a projected spatio-temporal spectrum *is not consistent* with the intersection of the free-space light-cone with a spectral hyperplane. Instead, it corresponds to the intersection of the free-space light-cone with a *non-planar* surface, whose projection onto the $({k}_{z},\frac{\omega}{c})$-plane is *not* a straight line, but is instead a curve given by the following equation:

By limiting the bandwidth of the ST wave packet $\mathrm{\Delta}\omega \ll {\omega}_{\mathrm{o}}$, the curved trajectory in the $({k}_{z},\frac{\omega}{c})$-plane given in Eq. (A2) becomes approximately a straight line, $\frac{\omega}{c}={k}_{\mathrm{o}}+({k}_{z}-{k}_{\mathrm{o}})\mathrm{tan}\text{\hspace{0.17em}}{\theta}_{a}$. To the first order, the projection onto the $({k}_{x},\frac{\omega}{c})$-plane given in Eq. (A1) can be approximated as a parabola in free space and in the material:

## 2. Measuring the Group Velocity of theSpace–Time Wave Packet

The ST wave packet and the reference pulse have group velocities $c\text{\hspace{0.17em}}\mathrm{tan}\text{\hspace{0.17em}}{\theta}_{\mathrm{a}}$ and $c$ in air, respectively. Maximal interference occurs when the delays encountered by both wave packets, when they reach the detector, are equal. If the ST wave packet propagates a distance ${L}_{1}$, and the reference pulse propagates a distance ${L}_{2}$, then ${L}_{1}/(c\text{\hspace{0.17em}}\mathrm{tan}\text{\hspace{0.17em}}{\theta}_{\mathrm{a}})={L}_{2}/c$. Inserting a distance $L$ of free space in the common path results in a loss of interference. The spatially resolved high-visibility fringes are regained by introducing a distance $\mathrm{\Delta}\ell $ in the reference path, to reach a new temporal balance $({L}_{1}+L)/(c\text{\hspace{0.17em}}\mathrm{tan}\text{\hspace{0.17em}}{\theta}_{\mathrm{a}})=({L}_{2}+L+\mathrm{\Delta}\ell )/c$, such that $\mathrm{\Delta}\ell /L=\mathrm{cot}\text{\hspace{0.17em}}{\theta}_{\mathrm{a}}-1$, from which ${\theta}_{\mathrm{a}}$ and hence ${v}_{\mathrm{g}}=c\text{\hspace{0.17em}}\mathrm{tan}\text{\hspace{0.17em}}{\theta}_{\mathrm{a}}$ in air can be obtained, and which can be rewritten as ${\tau}_{\mathrm{a}}={\tau}_{\mathrm{o}}+{\tau}_{\mathrm{d}}$, where ${\tau}_{\mathrm{a}}=L/(c\text{\hspace{0.17em}}\mathrm{tan}\text{\hspace{0.17em}}{\theta}_{\mathrm{a}})$, ${\tau}_{\mathrm{o}}=L/c$, and ${\tau}_{\mathrm{d}}=\mathrm{\Delta}\ell /c$.

Introducing a material of length $L$ and index $n$ in the common path diminishes the observed interference fringes. The group velocities of the ST wave packet and the reference pulse in the material are $c\text{\hspace{0.17em}}\mathrm{tan}\text{\hspace{0.17em}}{\theta}_{\mathrm{m}}$ and $\frac{c}{n}$, respectively. However, adding an extra free-space length $\mathrm{\Delta}\ell $ in the path of the reference pulse can restore the temporal balance such that $\mathrm{\Delta}\ell /L=\mathrm{cot}\text{\hspace{0.17em}}{\theta}_{\mathrm{m}}-\mathrm{cot}\text{\hspace{0.17em}}{\theta}_{\mathrm{a}}-(n-1)$. Substituting for $\mathrm{tan}\text{\hspace{0.17em}}{\theta}_{\mathrm{m}}$ in terms of $\mathrm{tan}\text{\hspace{0.17em}}{\theta}_{\mathrm{a}}$ from Eq. (A4), we obtain—within the parabolic approximation valid for narrow bandwidths—the following relationship:

## 3. Synthesis of Space–Time Wave Packets

The experimental setup combines spatial beam modulation with ultrafast pulse shaping. We start with a horizontally polarized pulsed laser beam from a Ti:Sapphire femtosecond laser (Tsunami, Spectra Physics) having a bandwidth of 8.5 nm and a central wavelength of 800 nm (pulse width, 100 fs). The beam is expanded spatially and collimated to produce a 25-mm-diameter plane-wave front. A reflective diffraction grating G (1200 lines/mm, area $25\times 25\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{2}$; Newport 10HG1200-800-1) disperses the spectrum in space, and the second diffraction order is directed to a reflective spatial light modulator (Hamamatsu X10468-02) through a cylindrical lens ${\mathrm{L}}_{1-y}$ in a $2f$ arrangement ($f=50\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{cm}$), which results in an estimated spectral uncertainty of $\delta \lambda \approx 30\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{pm}$ [Fig. 3(a)].

The SLM introduces controllable spatio-temporal spectral correlations. We implement on the SLM at the position associated with each wavelength $\lambda $ the linear phase distribution that corresponds to the required spatial frequency ${k}_{x}$ after introducing a reduction in scale by a factor of $4\times $ (i.e., we replace ${k}_{x}$ with ${k}_{x}/4$). We then subsequently compensate for this reduction at the output via a $4f$ telescope system that introduces a demagnification by a factor of $4\times $ along the $x$-direction (cylindrical lenses ${\mathrm{L}}_{2-x}$ and ${\mathrm{L}}_{3-x}$ of focal lengths 40 cm and 10 cm, respectively). This system implements the transformation ${k}_{x}\to 4{k}_{x}$, which compensates for the transformation ${k}_{x}\to {k}_{x}/4$ implemented at the SLM. An aperture A is introduced to reduce the temporal bandwidth to 0.2 nm. The modulated wave front is retro-reflected back to the grating through ${\mathrm{L}}_{1-y}$, where the pulse is reconstituted.

## 4. Measurement Procedure for the Space–Time Wave Packets

To measure the spatio-temporal spectrum ${|\tilde{E}({k}_{x},\lambda )|}^{2}$, we sample a portion of the retro-reflected field from the SLM after passing through the lens ${\mathrm{L}}_{1-y}$ via a beam splitter ${\mathrm{BS}}_{3}$. The field is directed through a spherical lens ${\mathrm{L}}_{4-\mathrm{s}}$ of focal length $f=7.5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{cm}$ to a CCD camera (${\mathrm{CCD}}_{2}$, the Imaging Source, DMK 72AUC02) in a $2f$ configuration to collimate the spectrum along $y$ and carry out a Fourier transform spatially along $x$.

The reference pulse is obtained from the initial laser pulses via beam splitter ${\mathrm{BS}}_{1}$. A neutral density filter adjusts the power, and the beam is expanded and spatially filtered using two spherical lenses (focal lengths 50 cm and 10 cm) and a pinhole (diameter 30 μm). The reference pulses traverse an optical delay line and are then combined with the ST wave packet via beam splitter ${\mathrm{BS}}_{4}$, and the resulting interferogram is recorded with a CCD camera (${\mathrm{CCD}}_{1}$, the Imaging Source, DMK 33UX178). The optical-material samples are placed in the common path after ${\mathrm{BS}}_{4}$. The solid samples (UV fused silica, BK7, and sapphire) are in the form of flat windows, whereas the liquid samples (distilled water, sea water, tetrahydrofuran, DMSO, benzene, and ${\mathrm{CS}}_{2}$) are placed in a cylindrical cell (Starna Cells 34-Q-50, interior length 50 mm, and diameter 19 mm). To accommodate the samples, ${\mathrm{CCD}}_{1}$ was kept at a fixed distance of 60 mm from ${\mathrm{BS}}_{4}$.

## 5. Measurement of Refractive Index of the Samples

For the liquid samples, interference is first recorded with an empty cylindrical cell placed in the reference path, and a ST wave packet of the spectral tilt angle $\theta =60\xb0$ is synthesized (any other value of $\theta $ can be used). The cylindrical cell is next filled with a liquid, resulting in a loss of the interference. The delay line in the reference arm is *reduced* to regain the interference. Using the measured delay and the inner length of the cylindrical cell (50 mm), we estimate the refractive index. The sea water was initially filtered using an Acrodisc PTFE CR syringe filter (25 mm diameter, 0.45 μm pore size) to remove any impurities. A similar approach was implemented for the solid samples. In Table 1, we list all the materials used and their relevant parameters.

## Funding

Office of Naval Research (ONR) (N00014-17-1-2458).

## Acknowledgment

We thank Scott Webster and H. Esat Kondakci for assistance.

## REFERENCES

**1. **L. Brillouin, *Wave Propagation and Group Velocity* (Academic, 1960).

**2. **R. W. Boyd and D. J. Gauthier, “Controlling the velocity of light pulses,” Science **326**, 1074–1077 (2009). [CrossRef]

**3. **L. V. Hau, S. E. Harris, Z. Dutton, and C. Behroozi, “Light speed reduction to 17 m per second in an ultracold atomic gas,” Nature **397**, 594–598 (1999). [CrossRef]

**4. **L. J. Wang, A. Kuzmich, and A. Dogariu, “Gain-assisted superluminal light propagation,” Nature **406**, 277–279 (2000). [CrossRef]

**5. **K. Y. Song, M. G. Herráez, and L. Thévenaz, “Gain-assisted pulse advancement using single and double Brillouin gain peaks in optical fibers,” Opt. Express **13**, 9758–9765 (2005). [CrossRef]

**6. **T. Baba, “Slow light in photonic crystals,” Nat. Photonics **2**, 465–473 (2008). [CrossRef]

**7. **K. L. Tsakmakidis, O. Hess, R. W. Boyd, and X. Zhang, “Ultraslow waves on the nanoscale,” Science **358**, eaan5196 (2017). [CrossRef]

**8. **J. Turunen and A. T. Friberg, “Propagation-invariant optical fields,” Prog. Opt. **54**, 1–88 (2010). [CrossRef]

**9. **H. E. Hernández-Figueroa, E. Recami, and M. Zamboni-Rached, eds., *Non-diffracting Waves* (Wiley-VCH, 2014).

**10. **R. Donnelly and R. Ziolkowski, “Designing localized waves,” Proc. R. Soc. Lond. A **440**, 541–565 (1993). [CrossRef]

**11. **S. Longhi, “Gaussian pulsed beams with arbitrary speed,” Opt. Express **12**, 935–940 (2004). [CrossRef]

**12. **P. Saari and K. Reivelt, “Generation and classification of localized waves by Lorentz transformations in Fourier space,” Phys. Rev. E **69**, 036612 (2004). [CrossRef]

**13. **H. E. Kondakci and A. F. Abouraddy, “Diffraction-free pulsed optical beams via space-time correlations,” Opt. Express **24**, 28659–28668 (2016). [CrossRef]

**14. **K. J. Parker and M. A. Alonso, “The longitudinal iso-phase condition and needle pulses,” Opt. Express **24**, 28669–28677 (2016). [CrossRef]

**15. **S. Szatmári, P. Simon, and M. Feuerhake, “Group-velocity-dispersion-compensated propagation of short pulses in dispersive media,” Opt. Lett. **21**, 1156–1158 (1996). [CrossRef]

**16. **H. Sõnajalg and P. Saari, “Suppression of temporal spread of ultrashort pulses in dispersive media by bessel beam generators,” Opt. Lett. **21**, 1162–1164 (1996). [CrossRef]

**17. **R. Danielius, A. Piskarskas, P. Di Trapani, A. Andreoni, C. Solcia, and P. Foggi, “Matching of group velocities by spatial walk-off in collinear three-wave interaction with tilted pulses,” Opt. Lett. **21**, 973–975 (1996). [CrossRef]

**18. **K. Reivelt and P. Saari, “Optical generation of focus wave modes,” J. Opt. Soc. Am. A **17**, 1785–1790 (2000). [CrossRef]

**19. **M. A. Porras, G. Valiulis, and P. Di Trapani, “Unified description of Bessel X waves with cone dispersion and tilted pulses,” Phys. Rev. E **68**, 016613 (2003). [CrossRef]

**20. **C. J. Zapata-Rodríguez and M. A. Porras, “X-wave bullets with negative group velocity in vacuum,” Opt. Lett. **31**, 3532–3534 (2006). [CrossRef]

**21. **P. Saari, “Reexamination of group velocities of structured light pulses,” Phys. Rev. A **97**, 063824 (2018). [CrossRef]

**22. **F. Bonaretti, D. Faccio, M. Clerici, J. Biegert, and P. Di Trapani, “Spatiotemporal amplitude and phase retrieval of Bessel-X pulses using a Hartmann-Shack sensor,” Opt. Express **17**, 9804–9809 (2009). [CrossRef]

**23. **P. Bowlan, H. Valtna-Lukner, M. Lõhmus, P. Piksarv, P. Saari, and R. Trebino, “Measuring the spatiotemporal field of ultrashort Bessel-X pulses,” Opt. Lett. **34**, 2276–2278 (2009). [CrossRef]

**24. **K. B. Kuntz, B. Braverman, S. H. Youn, M. Lobino, E. M. Pessina, and A. I. Lvovsky, “Spatial and temporal characterization of a Bessel beam produced using a conical mirror,” Phys. Rev. A **79**, 043802 (2009). [CrossRef]

**25. **J. Hebling, G. Almási, I. Z. Kozma, and J. Kuhl, “Velocity matching by pulse front tilting for large-area THz-pulse generation,” Opt. Express **10**, 1161–1166 (2002). [CrossRef]

**26. **J. Hebling, K.-L. Yeh, M. C. Hoffmann, B. Bartal, and K. A. Nelson, “Generation of high-power terahertz pulses by tilted-pulse-front excitation and their application possibilities,” J. Opt. Soc. Am. B **25**, B6–B19 (2008). [CrossRef]

**27. **J. A. Fülöp and J. Hebling, “Applications of tilted-pulse-front excitation,” in *Recent Optical and Photonic Technologies*, K. Y. Kim, ed. (InTech, 2010).

**28. **P. Saari and K. Reivelt, “Evidence of X-shaped propagation-invariant localized light waves,” Phys. Rev. Lett. **79**, 4135–4138 (1997). [CrossRef]

**29. **I. Alexeev, K. Y. Kim, and H. M. Milchberg, “Measurement of the superluminal group velocity of an ultrashort Bessel beam pulse,” Phys. Rev. Lett. **88**, 073901 (2002). [CrossRef]

**30. **P. Di Trapani, G. Valiulis, A. Piskarskas, O. Jedrkiewicz, J. Trull, C. Conti, and S. Trillo, “Spontaneously generated X-shaped light bullets,” Phys. Rev. Lett. **91**, 093904 (2003). [CrossRef]

**31. **D. Faccio, M. A. Porras, A. Dubietis, F. Bragheri, A. Couairon, and P. Di Trapani, “Conical emission, pulse splitting, and X-wave parametric amplification in nonlinear dynamics of ultrashort light pulses,” Phys. Rev. Lett. **96**, 193901 (2006). [CrossRef]

**32. **D. Faccio, A. Averchi, A. Couairon, M. Kolesik, J. Moloney, A. Dubietis, G. Tamosauskas, P. Polesana, A. Piskarskas, and P. Di Trapani, “Spatio-temporal reshaping and X wave dynamics in optical filaments,” Opt. Express **15**, 13077–13095 (2007). [CrossRef]

**33. **M. Dallaire, N. McCarthy, and M. Piché, “Spatiotemporal Bessel beams: theory and experiments,” Opt. Express **17**, 18148–18164 (2009). [CrossRef]

**34. **O. Jedrkiewicz, Y.-D. Wang, G. Valiulis, and P. Di Trapani, “One dimensional spatial localization of polychromatic stationary wave-packets in normally dispersive media,” Opt. Express **21**, 25000–25009 (2013). [CrossRef]

**35. **H. E. Kondakci and A. F. Abouraddy, “Diffraction-free space-time beams,” Nat. Photonics **11**, 733–740 (2017). [CrossRef]

**36. **H. E. Kondakci and A. F. Abouraddy, “Airy wavepackets accelerating in space-time,” Phys. Rev. Lett. **120**, 163901 (2018). [CrossRef]

**37. **H. E. Kondakci and A. F. Abouraddy, “Optical space-time wave packets of arbitrary group velocity in free space,” arXiv:1810.08893 (2018).

**38. **H. E. Kondakci and A. F. Abouraddy, “Self-healing of space-time light sheets,” Opt. Lett. **43**, 3830–3833 (2018). [CrossRef]

**39. **B. Bhaduri, M. Yessenov, and A. F. Abouraddy, “Meters-long propagation of diffraction-free space-time light sheets,” Opt. Express **26**, 20111–20121 (2018). [CrossRef]

**40. **H. E. Kondakci, M. Yessenov, M. Meem, D. Reyes, D. Thul, S. R. Fairchild, M. Richardson, R. Menon, and A. F. Abouraddy, “Synthesizing broadband propagation-invariant space-time wave packets using transmissive phase plates,” Opt. Express **26**, 13628–13638 (2018). [CrossRef]

**41. **M. Yessenov, B. Bhaduri, H. E. Kondakci, and A. F. Abouraddy, “Classification of propagation-invariant space-time light-sheets in free space: Theory and experiments,” arXiv:1809.08375 (2018).

**42. **M. Clerici, D. Faccio, A. Lotti, E. Rubino, O. Jedrkiewicz, J. Biegert, and P. D. Trapani, “Finite-energy, accelerating Bessel pulses,” Opt. Express **16**, 19807–19811 (2008). [CrossRef]

**43. **J. P. Torres, M. Hendrych, and A. Valencia, “Angular dispersion: an enabling tool in nonlinear and quantum optics,” Adv. Opt. Photon. **2**, 319–369 (2010). [CrossRef]

**44. **D. A. Fishman, C. M. Cirloganu, S. Webster, L. A. Padilha, M. Monroe, D. J. Hagan, and E. W. Van Stryland, “Sensitive mid-infrared detection in wide-bandgap semiconductors using extreme non-degenerate two-photon absorption,” Nat. Photonics **5**, 561–565 (2011). [CrossRef]

**45. **M. Reichert, A. Smirl, G. Salamo, D. Hagan, and E. Van Stryland, “Observation of nondegenerate two-photon gain in GaAs,” Phys. Rev. Lett. **117**, 073602 (2016). [CrossRef]

**46. **G. M. Hale and M. R. Querry, “Optical constants of water in the 200-nm to 200-μm wavelength region,” Appl. Opt. **12**, 555–563 (1973). [CrossRef]

**47. **M. J. Dodge, “Refractive properties of magnesium fluoride,” Appl. Opt. **23**, 1980–1985 (1984). [CrossRef]

**48. **I. H. Malitson, “Interspecimen comparison of the refractive index of fused silica,” J. Opt. Soc. Am. **55**, 1205–1208 (1965). [CrossRef]

**49. **I. Z. Kozma, P. Krok, and E. Riedle, “Direct measurement of the group-velocity mismatch and derivation of the refractive-index dispersion for a variety of solvents in the ultraviolet,” J. Opt. Soc. Am. B **22**, 1479–1485 (2005). [CrossRef]

**50. **K. Moutzouris, M. Papamichael, S. C. Betsis, I. Stavrakas, G. Hloupis, and D. Triantis, “Refractive, dispersive and thermo-optic properties of twelve organic solvents in the visible and near-infrared,” Appl. Phys. B **116**, 617–622 (2014). [CrossRef]

**51. **R. Jedamzik, S. Reichel, and P. Hartmann, “Optical glass with tightest refractive index and dispersion tolerances for high-end optical designs,” Proc. SPIE **8982**, 89821F (2014). [CrossRef]

**52. **S. Kedenburg, M. Vieweg, T. Gissibl, and H. Giessen, “Linear refractive index and absorption measurements of nonlinear optical liquids in the visible and near-infrared spectral region,” Opt. Mater. Express **2**, 1588–1611 (2012). [CrossRef]

**53. **I. H. Malitson, “Refraction and dispersion of synthetic sapphire,” J. Opt. Soc. Am. **52**, 1377–1379 (1962). [CrossRef]