## Abstract

Space-time wave packets are a class of pulsed optical beams that are diffraction-free and dispersion-free in free space by virtue of introducing a tight correlation between the spatial and temporal degrees of freedom of the field. Such wave packets have been recently synthesized in a novel configuration that makes use of a spatial light modulator to realize the required spatio-temporal correlations. This arrangement combines pulse-modulation and beam-shaping to assign one spatial frequency to each wavelength according to a prescribed correlation function. Relying on a spatial light modulator results in several limitations by virtue of their pixelation, small area, and low energy-handling capability. Here we demonstrate the synthesis of space-time wave packets with one spatial dimension kept uniform – that is, light sheets – using transparent transmissive phase plates produced by a gray-scale lithography process. We confirm the diffraction-free behavior of wave packets having a bandwidth of 0.25 nm (filtered from a typical femtosecond Ti:sapphire laser) and 30 nm (a multi-terawatt femtosecond laser). This work paves the way for developing versatile high-energy light bullets for applications in nonlinear optics and laser machining.

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

## 1. Introduction

Freely propagating optical fields diffract: their spatial extent increases and their fine features blur with propagation [1]. Combating diffraction in free space has typically relied on exploiting beams having specific transverse profiles that are propagation-invariant solutions of the wave equation. Examples of monochromatic diffraction-free beams include those whose profiles conform to Bessel [2], Mathieu [3], Weber [4], and Airy [5] functions (see [6] for a review). In the case of pulsed beams (wave packets or optical bullets), different functional forms of the field in space and time have been found to be propagation-invariant (both diffraction-free and dispersion-free). The earliest examples are the Brittingham focus-wave mode [7] and MacKinnon’s wave packets [8], and subsequently X-waves [9–11], Bessel pulses [12], among other examples [13], in addition to more recent proposals [14, 15]. Despite the wide diversity of specific propagation-invariant wave packets identified theoretically [16,17], underlying them all is a *single fundamental characteristic*: in the spatio-temporal spectral domain, the spatial frequencies (that establish the transverse beam profile in space) are tightly correlated with the temporal frequencies (that determine the pulse profile in time) [14–18]. We thus call all such optical fields ‘space-time’ (ST) wave packets. The spatio-temporal spectra of all such free-space wave packets lie at the intersection of the light-cone with spectral hyperplanes, and the various families of previously identified ST wave packets can be generated by tilting the hyperplane with respect to the light-cone [18, 19]. In other words, instead of the spatio-temporal spectrum of the wave packet occupying a patch on the light-cone as is typical for pulsed beams (where the spatial and temporal degrees of freedom are almost independent), the spectra of ST wave packets are confined to reduced-dimensionality trajectories.

Synthesizing such ST wave packets presents considerable difficulties, especially in ensuring the realization of a tight correlation between spatial and temporal frequencies. We have recently provided an experimental solution to this challenge that combines elements from beam-shaping [20,21] and ultrafast pulse-shaping [22,23]: we use a two-dimensional (2D) spatial light modulator (SLM) to modulate the phase in the *orthogonal* direction to that of the wavelength spread (similarly to so-called 4− *f* pulse shaping [24–26]) to imprint the spatial frequency that is to be associated with each particular wavelength. Using this approach we produced diffraction-free ST wave packets of transverse width of 14 *μ*m that propagate for 100 mm, for example, in addition to hollow pulsed beams with a diffraction-free axial null [18]. More recently, we have produced *non*-accelerating Airy ST wave packets that accelerate instead in space-time; that is, this Airy wave packet travels in a straight line without the usual transverse displacement, and instead accelerates in the local frame of the propagating pulse [19]. Moreover, several very recent theoretical studies have further enriched this newly developing field of ST wave packets [27–31].

As described above, our previous experimental realizations of ST wave packets have made use of a SLM to impress a 2D phase distribution on the impinging spectrally resolved beam in order to introduce the desired spatio-temporal spectral correlations. SLMs have an obvious advantage with respect to their capability for dynamically varying the implemented phase distribution when compared to a static phase plate. Nevertheless, relying on SLMs presents some restrictions that stem from their fundamental limitations. First, because most SLMs make use of liquid-crystal-based technology, their low energy-handling capability precludes them from being utilized to synthesize high-energy ST wave packets. Second, the pixel size sets a *lower* limit on the width of the uncertainty of the spectral correlations (which dictates the diffraction-free propagation distance), while the finite SLM aperture size sets an *upper* limit on the utilizable spectral bandwidth as well as the quasi-diffraction-free propagation distance.

Here we synthesize ST wave packets using transmissive phase plates that are lithographically inscribed in a transparent polymer film via a gray-scale lithography process. We thus overcome the major limitations of SLMs and implement phase distributions at a higher spatial resolution (∼ 3 *μ*m) and cover a larger area (∼ 25 × 25 mm^{2}). Consequently, we are now in a position to synthesize ST wave packets with larger spectral bandwidth than previously possible. In Refs. [18, 19] we produced ST wave packets with a bandwidth of < 1 nm. Here, we report on the synthesis of diffraction-free ST wave packets with bandwidths up to ~30 nm. In the experiments, we make use of a typical femtosecond Ti:sapphire laser for a smaller-bandwidth experiment and a Ti:sapphire-based multi-terawatt femtosecond laser for the higher bandwidth cases. Furthermore, the use of a transparent transmissive phase plate opens the path towards exploiting high-energy pulses to synthesize ST wave packets (indeed, phase plates have been used previously with high-energy pulses to produce nondiffracting double helical beams and helical filaments [32, 33]). In all cases, we synthesize the ST wave packets – for simplicity – in the form of light sheets with an approximate Gaussian beam profile at the pulse center, where the diffraction-free behavior is observed along one dimension, while the field is uniform along the second transverse dimension [18, 19].

The paper is organized as follows. We first present a brief overview of the concept of ST wave packets and delineate the role of ST correlations in producing propagation-invariant pulsed beams. Next we describe the design and fabrication of the phase plates before describing the optical setup used for wave-packet synthesis. We then present our measurement results confirming the diffraction-free behavior of the produced ST wave packets before outlining our conclusions.

## 2. Propagation-invariant ST light sheets

We consider here for simplicity one-dimensional (1D) optical beams or light sheets with electric field *E* (*x, z, t*), where *x* is the transverse coordinate, *z* is the axial coordinate, and *t* is time (the field is uniform along the other transverse coordinate *y*). The initial field at *z* = 0, *E* (*x*, 0, *t*), has a spatio-temporal spectrum *Ẽ*(*k _{x}*,

*ω*), where

*k*is the spatial frequency or the transverse component of the wave vector and

_{x}*ω*is the temporal frequency, which are subject to the dispersion relationship ${k}_{x}^{2}+{k}_{z}^{2}={\left(\omega /c\right)}^{2}$ that represents the surface of a cone known as the light-cone [Fig. 1(a)]; here

*k*is the axial component of the wave vector and

_{z}*c*is the speed of light in vacuum. Upon free-propagation [1], the field evolves axially according to

The basic concept underlying the propagation-invariance of a ST wave packet is the tight correlation between the spatial and temporal frequencies *k _{x}* and

*ω*, respectively. An optical

*beam*has a finite spatial bandwidth Δ

*k*, whereas an optical

_{x}*pulse*has a finite temporal bandwidth Δ

*ω*. A

*pulsed*optical

*beam*(or wave packet) combines finite spatial and temporal bandwidths, and thus the associated spatio-temporal spectrum occupies in general a 2D patch

*Ẽ*(

*k*) on the side of the light-cone.

_{x}, ωIf correlations are introduced between the spatial and temporal degrees of freedom such that each spatial frequency *k _{x}* is associated with a temporal frequency

*ω*[Fig. 1(b)] and thus

*Ẽ*(

*k*) →

_{x}, ω*Ẽ*(

*k*)

_{x}*δ*(

*ω*−

*f*(

*k*)), the spatio-temporal spectrum is then confined to a 1D trajectory instead of a 2D patch; here (

_{x}*ω*=

*f*(

*k*) is the functional form of the correlation between

_{x}*ω*and

*k*. This relationship is taken here to be one-to-one between |

_{x}*k*| and

_{x}*ω*. If this spatio-temporal spectral trajectory defined by

*ω*=

*f*(

*k*) is the result of the intersection of the light-cone with a plane that is parallel to the

_{x}*k*-axis and orthogonal to the (

_{x}*k*)-plane as shown in Fig. 1(a), then

_{z}, ω*any*wave packet whose spatio-temporal spectrum lies along this trajectory can be shown to be diffraction-free

*and*dispersion-free [18]. This can be understood by noting that the projection of this trajectory onto the (

*k*)-plane takes the form of a straight line

_{z}, ω*ω*/

*c*=

*k*

_{o}+ (

*k*−

_{z}*k*

_{o})

*v*

_{g}/

*c*where

*k*

_{o}is a constant wave number and

*v*

_{g}is the group velocity of the wave packet, which is related to the tilt angle

*φ*of the spectral plane with the horizontal (

*k*−

_{x}*k*

_{z})-plane,

*v*

_{g}=

*c*tan

*φ*. The resulting wave packet thus takes the form

*k*and

_{z}*ω*described above leads to rigid translation of the wave packet along

*z*at a group velocity

*v*

_{g}without diffraction or dispersion [14, 34–36]. A typical spatio-temporal profile |

*E*(

*x, z, t*)|

^{2}of such a propagation-invariant ST wave packet is shown in Fig. 1(c) at a fixed value of

*z*.

## 3. Synthesis of space-time wave packets

#### 3.1. Optical setup

We depict schematically in Fig. 2 the optical setup to synthesize the ST wave packets using a transmissive phase plate that modulates the wave front such as to introduce the desired spatio-temporal correlations. The system can be divided into two major sections that carry out the tasks of field *synthesis* and *analysis*, with the latter performed in physical space (*x, z, t*) and in the spatio-temporal spectral domain (*k _{x}, k_{z}, ω*). As described in the previous section, the synthesis of a ST wave packet requires associating each wavelength

*λ*with a particular spatial frequency

*k*. To achieve this, we start from a pulsed plane wave whose spatio-temporal spectrum is thus separable, having approximately the form

_{x}*Ẽ*(

*k*) ∼

_{x}, λ*δ*(

*k*(

_{x}Ẽ*λ*). The pulse spectrum

*Ẽ*(

*λ*) having a bandwidth of Δ

*λ*is spread spatially along one dimension via a diffraction grating G

_{1}and directed to the transmissive phase plate (that replaces the SLM in [18, 19]) by a cylindrical lens L

_{1}in a 2

*f*configuration. The phase plate then imparts a 2D phase distribution to the spectrally resolved field before we reconstitute the pulse via a second cylindrical lens L

_{2}(identical to L

_{1}) and a second diffraction grating G

_{2}(identical to G

_{1}). At this stage, the spatio-temporal spectrum has the form

*Ẽ*(

*k*) ∼

_{x}, λ*Ẽ*(

*k*)

_{x}*δ*(

*λ*−

*f*(

*k*)). Note that the finite sizes of the apertures and the minimum feature size in the phase plate prevent the realization of an idealized delta-function correlation. Instead, the delta function is replaced by a narrow function that introduces a finite spectral uncertainty

_{x}*δλ*, which determines the maximum diffraction-free propagation length [14].

The resulting ST wave packet is split into two paths by a beam splitter for analysis. First, the diffraction-free behavior is confirmed in physical space by scanning a CCD camera (CCD_{1}) along the propagation axis after a 4 *f* telescope. The measured intensity is thus the time-averaged value *I*(*x, z*) = ∫*dt* |*E* (*x, z, t*)|^{2}. Second, the spatio-temporal spectrum is measured by sampling the field *before* G_{2} (i.e., while the field is still spectrally resolved) and implementing a spatial Fourier transform via a lens in a 2 *f* configuration to resolve the spatial spectrum along *k _{x}*. The spatio-temporal spectrum |

*Ẽ*(

*k*)|

_{x}, λ^{2}is then recorded by a camera (CCD

_{2}).

Note that this system is similar conceptually to that implemented in [18, 19] where the arrangement was folded back on itself by employing a *reflective* SLM, such that the gratings G_{1} and G_{2} coincide and the lenses L_{1} and L_{2} coincide.

#### 3.2. Phase plate design and fabrication

The phase plates were fabricated using gray-scale lithography [37]. First, a positive-tone photoresist (Shipley 1813) [38] was spin-coated on an RCA-cleaned 2-inch-diameter D263 glass wafer at 1000 rpm for 60 seconds, which gives a film thickness of ≈ 2.6 *μ*m. The sample was then baked in an oven at 110 °C for 30 minutes. A calibration was conducted to determine the depth of the photoresist after development for a given exposure dose as described previously [39, 40]. A laser pattern generator (Heidelberg Instruments) [41] was used to write the design. The exposed sample was then developed in AZ 1:1 solution [42] for 90 seconds followed by DI water rinse.

We produced two phase plates with this process. The first phase plate was designed to accommodate a bandwidth of 2 nm, that of a typical femtosecond Ti:Sa laser, with dimensions 15×12 mm^{2}. By spatially filtering the spectrum incident on the phase plate, different bandwidths can be utilized. In our experiment, we restrict its use to a bandwidth of ~0.25 nm. The second plate had dimensions of 25 × 25 mm^{2} and was designed to accommodate a larger bandwidth of ∼40 nm. The multi-terawatt femtosecond laser used here has a bandwidth of ∼30 nm. Optical micrographs of the fabricated devices are given in Fig. 4(b) and 6(b).

## 4. Space-time wave-packet synthesis and characterization

#### 4.1. Space-time wave packet with 0.25-nm bandwidth

The experimental arrangement for synthesizing and characterizing ST wave packets, which applies for bandwidths of Δ*λ* ∼0.25 nm or 30 nm, is shown in Fig. 2; however, different values for the component parameters are selected for the two bandwidths targeted. The laser used with the first phase plate for a bandwidth of Δ*λ* ∼0.25 nm is filtered from a generic pulsed Ti:Sa laser with a pulse width of ∼100 fs (FWHM) at a wavelength of ~800 nm (Tsunami, Spectra Physics). We first characterize the initial laser pulse in the spectral and temporal domains using the FROG technique (employing a GRENOUILLE 8-50, Swamp Optics); see Fig. 3 for measurement results. Crucially, minimal spectral phase and temporal chirping are observed over the spectral bandwidth exploited in our experiment prior to the phase modulation applied by the phase plate. The beam is directed to a reflective diffraction grating (G_{1}) with 1200 lines/mm and dimensions 25 × 25 mm^{2} (Newport 10HG1200-800-1) at an incidence angle of 68° to spread the pulse spectrum in space. We select the second diffraction-order at 81° that provides a spatial dispersion of 0.064 nm/mrad. The diffracted beam is normally incident on the transmissive phase plate after a cylindrical lens L_{1−}* _{y}* of focal length

*f*= 50 cm in a 2

*f*configuration (we denote the direction of the spatially dispersed spectrum

*y*).

The phase plate is designed to synthesize ST wave packets whose spatio-temporal spectral loci lie along the intersection of the light-cone with a plane that makes an angle of *φ* = 90° with the *k _{z}*-axis; that is, the izo-

*k*plane [Fig. 1(a)]. The phase pattern thus implements the dispersion relationship ${k}_{\mathrm{o}}^{2}={\left(2\pi /\lambda \right)}^{2}-{k}_{\mathrm{x}}^{2}$, which produces the required correlation between

_{z}*k*and

_{x}*λ*. Here

*k*

_{o}is a fixed axial wave number that corresponds to the maximum wavelength used and is associated with the spatial frequency

*k*=0. For experimental convenience, we implement the spatial frequencies

_{x}*k*on the phase plate with a reduction in scale by a factor 4× and then subsequently compensate for that reduction at the output via a 4

_{x}*f*telescope system that introduces a demagnification by a factor of 4× along the

*x*direction [18].

After traversing the phase plate, the field is directed through a cylindrical lens L_{2−}* _{y}* (that is identical to L

_{1−}

*) in a 2*

_{y}*f*configuration to a beam splitter BS

_{1}to split the beam in two paths for characterization. In one path, we confirm that the desired spatio-temporal correlation between

*k*and

_{x}*λ*has been realized by measuring the spatio-temporal spectrum |

*Ẽ*(

*k*)|

_{x}, λ^{2}. In this path, we arrange for a spherical lens L

_{5−s}with a focal length of

*f*=10 cm in a 2

*f*configuration to collimate the spectrum along

*y*and carry out a Fourier transform spatially along

*x*. A CCD

_{2}camera (The ImagingSource, DMK 33UX178) placed in the focal plane of L

_{5−s}thus directly images the spatio-temporal spectrum |

*Ẽ*(

*k*)|

_{x}, λ^{2}, which is presented in Fig. 4(a). The data confirms that we obtain the targeted hyperbolic correlation function between

*k*and

_{x}*λ*that is imposed by the phase distribution [Fig. 4(b)] provided by the phase plate.

In the other path after BS_{1}, we characterize the ST wave packet in physical space to confirm its diffraction-free behavior by reconstituting the pulse via a diffraction grating G_{2} (that is identical to G_{1}). After G_{2}, an imaging 4 *f* system translates the beam and implements the above-described demagnification factor 4× (the values of all the focal lengths of the lenses are provided in the inset of Fig. 2). Starting from the output plane of this imaging system, we capture the transverse intensity distribution with a camera CCD_{1} (The ImagingSource, DMK 33UX178) while scanning its axial position along the *z* direction. The measurement corresponds to the time-averaged intensity *I*(*x, z* )= ∫*dt* |*E* (*x, z, t*)|^{2} (the intensity distribution is uniform along *y*). The measured intensity distribution is plotted in Fig. 4(c) for a scanned axial distance scanned of 25 mm, confirming the expected diffraction-free behavior. The detailed structure of such a ST wave packet in space and time is presented in [18].

#### 4.2. Space-time wave packet with 30-nm bandwidth

The experimental setup for synthesizing and characterizing the ST wave packets for the case of the larger bandwidth Δ*λ* ∼ 30 nm is also that shown in Fig. 2, except for the change in the values of the focal lengths of the lenses used (provided in the inset of Fig. 2). The laser source of the initial pulses here is the Multi-Terawatt Femtosecond Laser (MTFL) [43, 44] housed at the University of Central Florida. This Ti:Sa-based chirped pulse amplification system consists of 3 stages of amplification and is capable of delivering 500-mJ-energy pulses at a wavelength of 800 nm having a pulse width of ∼50 fs (FWHM) at a repetition rate of 10 Hz. The beam used in this experiment was extracted directly from the regenerative amplifier (first stage of amplification) at an energy of only 1 mJ, which was further attenuated to avoid saturating the CCD camera. We have characterized these pulses in the temporal and spectral domains using the FROG technique as described in the previous section, and the measurements are presented in Fig. 5. We note in Fig. 5(b) the existence of a spectral phase extending across the bandwidth we utilize. Once the spatio-temporal spectral correlations are introduced between *k _{x}* and

*λ*by the phase plate, this spectral phase will be shared between the spatial and temporal degrees of freedom. We thus expect a slight tilt in the propagating wave front.

Because of the larger spectral bandwidth involved here, we modify the experimental parameters for spreading the spectrum in space and also use a phase plate of larger dimensions. The input pulsed beam is directed at an incidence angle of 51° to the grating G_{1} that has in this case 1800 lines/mm and dimensions 25×25 mm^{2} (Thorlabs GR25-1850), and select the first diffraction-order at an angle of 41°, which is directed to a phase plate having dimensions 25 × 25 mm^{2}. The field at the phase plate extends over a width of ≈24 mm. The phase plate is designed to produce a ST wave packet whose spatio-temporal spectrum that lies along the hyperbolic trajectory resulting from the intersection of the light cone with a spectral plane tilted an angle *φ*=45.3° with respect to the *k _{z}*-axis [Fig. 1(a)].

The synthesized wave packet is characterized in physical space *I*(*x, z*) and in the spatio-temporal domain |*Ẽ*(*k _{x}, λ*)|

^{2}using the same previously described methodology. Measurements of the spatio-temporal spectrum |

*Ẽ*(

*k*)|

_{x}, λ^{2}are shown in Fig. 6(a). The data confirms that the targeted hyperbolic correlation function between

*k*and

_{x}*λ*has been realized over the full bandwidth after making use of the phase distribution depicted in Fig. 6(b). The measured time-averaged intensity

*I*(

*x, z*) while scanning the camera CCD

_{1}axially is plotted in Fig. 6(c). The diffraction-free behavior is observed over an axial distance of ∼45 mm. Note that a slight tilt in the wave front is seen in the data as expected from the spectral phase of the initial femtosecond pulses [Fig. 5(b)].

## 5. Conclusion

In conclusion, we have synthesized broadband one-dimensional diffraction-free ST wave packets (or light sheets) using phase plates produced by a gray-scale lithography process. The phase plates are larger in overall size (25×25 mm^{2}) than a typical SLM, have smaller minimal feature size (∼ 3 *μ*m), and have higher diffraction efficiency. These characteristics have allowed us to synthesize ST wave packets of two types. First, we exploited the fine spatial features of the phase plate to produce a ST wave packet with narrow bandwidth Δ*λ ~* 0.25. Second, we made use of the additional advantage of the large area of the phase plate to produce a ST wave packet of considerably large bandwidth Δ*λ* 30 nm compared to previous demonstrations that exploited SLMs (Δ*λ* ∼ 1 nm in Refs. [18, 19]). In both cases, we maintain a low spectral uncertainty *δλ* necessary for the diffraction-free propagation. The possibility of exploiting transparent transmissive phase plates paves the way to implementing high-energy ST wave packets that can be useful in nonlinear optical applications, including laser filamentation studies and laser machining or additive manufacturing. Future work will be directed at increasing the diffraction-free propagation distance, realizing ST wave packets with control over both transverse dimensions (ST needles in lieu of light sheets), and synthesizing high-energy ST wave packets. Finally, we note that the correlations introduced between the spatial and temporal degrees of freedom here are the continuous analog of the correlations between discretized optical degrees of freedom known as ‘classical entanglement’ [45, 46] that has shown significant recent utility in studying optical coherence theory.

## Funding

U.S. Office of Naval Research (ONR) (N00014-17-1-2458); U.S. Office of Naval Research (ONR) (N66001-10-1-4065); U.S. Army Research Office (ARO) MURI program on Light Filamentation Science (W911NF1110297); HEL JTO (FA9550-11-1-0001).

## Acknowledgments

We thank the staff at the University of Utah nanofabrication facility, especially Brian Baker and Steve Pritchett, for useful suggestions to improve the phase-plate fabrication process.

## References and links

**1. **B. E. A. Saleh and M. C. Teich, *Fundamentals of Photonics* (Wiley, 2007).

**2. **J. Durnin, J. J. Miceli, and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. **58**, 1499–1501 (1987). [CrossRef] [PubMed]

**3. **M. A. Bandres, J. C. Gutiérrez-Vega, and S. Chávez-Cerda, “Parabolic nondiffracting optical wave fields,” Opt. Lett. **29**, 44–46 (2004). [CrossRef] [PubMed]

**4. **B. M. Rodríguez-Lara, “Normalization of optical Weber waves and Weber-Gauss beams,” J. Opt. Soc. Am. A **27**, 327–332 (2010). [CrossRef]

**5. **G. A. Siviloglou and D. N. Christodoulides, “Accelerating finite energy Airy beams,” Opt. Lett. **32**, 979–981 (2007). [CrossRef] [PubMed]

**6. **U. Levy, S. Derevyanko, and Y. Silberberg, “Light modes of free space,” Prog. Opt. **61**, 237–281 (2016). [CrossRef]

**7. **J. N. Brittingham, “Focus wave modes in homogeneous Maxwell’s equations: Transverse electric mode,” J. Appl. Phys. **54**, 1179–1189 (1983). [CrossRef]

**8. **L. Mackinnon, “A nondispersive de Broglie wave packet,” Found. Phys. **8**, 157–176 (1978). [CrossRef]

**9. **J.-Y. Lu and J. F. Greenleaf, “Nondiffracting X waves–exact solutions to free-space scalar wave equation and their finite aperture realizations,” IEEE Trans. Ultrason. Ferroelec. Freq. Control **39**, 19–31 (1992). [CrossRef]

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

**11. **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] [PubMed]

**12. **C. J. R. Sheppard, “Generalized Bessel pulse beams,” J. Opt. Soc. Am. A **19**, 2218–2222 (2002). [CrossRef]

**13. **D. N. Christodoulides, N. K. Efremidis, P. Di Trapani, and B. A. Malomed, “Bessel X waves in two- and three-dimensional bidispersive optical systems,” Opt. Lett. **29**, 1446–1448 (2004). [CrossRef] [PubMed]

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

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

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

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

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

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

**20. **C. Maurer, A. Jesacher, S. Bernet, and M. Ritsch-Marte, “What spatial light modulators can do for optical microscopy,” Laser Photon. Rev. **5**, 81–101 (2011). [CrossRef]

**21. **A. Forbes, A. Dudley, and M. McLaren, “Creation and detection of optical modes with spatial light modulators,” Adv. Opt. Photon. **8**, 200–227 (2016). [CrossRef]

**22. **A. M. Weiner, “Femtosecond pulse shaping using spatial light modulators,” Rev. Sci. Instrum. **71**, 1929–1960 (2000). [CrossRef]

**23. **A. M. Weiner, *Ultrafast Optics* (Wiley, Weinheim, 2009). [CrossRef]

**24. **R. M. Koehl, T. Hattori, and K. A. Nelson, “Automated spatial and temporal shaping of femtosecond pulses,” Opt. Commun. **157**, 57–61 (1998). [CrossRef]

**25. **T. Feurer, J. C. Vaughan, R. M. Koehl, and K. A. Nelson, “Multidimensional control of femtosecond pulses by use of a programmable liquid-crystal matrix,” Opt. Lett. **27**, 652–654 (2002). [CrossRef]

**26. **B. J. Sussman, R. Lausten, and A. Stolow, “Focusing of light following a 4–f pulse shaper: Considerations for quantum control,” Phys. Rev. A **77**, 043416 (2008). [CrossRef]

**27. **L. J. Wong and I. Kaminer, “Abruptly focusing and defocusing needles of light and closed-form electromagnetic wavepackets,” ACS Photon. **4**, 1131–1137 (2017). [CrossRef]

**28. **L. J. Wong and I. Kaminer, “Ultrashort tilted-pulse-front pulses and nonparaxial tilted-phase-front beams,” ACS Photon. **4**, 2257–2264 (2017). [CrossRef]

**29. **A. Sainte-Marie, O. Gobert, and F. Quere, “Controlling the velocity of ultrashort light pulses in vacuum through spatio-temporal couplings,” Optica **4**, 1298–1304 (2017). [CrossRef]

**30. **M. A. Porras, “Gaussian beams diffracting in time,” Opt. Lett. **42**, 4679–4682 (2017). [CrossRef] [PubMed]

**31. **N. K. Efremidis, “Spatiotemporal diffraction-free pulsed beams in free-space of the Airy and Bessel type,” Opt. Lett. **23**, 5038–5041 (2017). [CrossRef]

**32. **N. Barbieri, M. Weidman, G. Katona, M. Baudelet, Z. Roth, E. Johnson, G. Siviloglou, D. Christodoulides, and M. Richardson, “Double helical laser beams based on interfering first-order Bessel beams,” J. Opt. Soc. Am. A **28**, 1462–1469 (2011). [CrossRef]

**33. **N. Barbieri, Z. Hosseinimakarem, K. Lim, M. Durand, M. Baudelet, E. Johnson, and M. Richardson, “Helical filaments,” Appl. Phys. Lett. **104**, 261109 (2014). [CrossRef]

**34. **P. A. Bélanger, “Packetlike solutions of the homogeneous-wave equation,” J. Opt. Soc. Am. A **1**, 723–724 (1984). [CrossRef]

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

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

**37. **C. McKenna, K. Walsh, M. Crain, and J. Lake, “Maskless direct write grayscale lithography for MEMS applications,” in Micro/Nano Symposium (UGIM), 2010 18th Biennial University/Government/Industry, pp. 1–4. IEEE (2010).

**38. **Data sheet of Shipley 1813 photoresist: http://www.microchem.com/PDFs_Dow/S1800.pdf

**39. **N. Mohammad, M. Meem, X. Wan, and R. Menon, “Full-color, large area, transmissive holograms enabled by multi-level diffractive optics,” Sci. Rep. **7**, 5789 (2017). [CrossRef] [PubMed]

**40. **P. Wang, J. A. Dominguez-Caballero, D. J. Friedman, and R. Menon, “A new class of multi-bandgap high-efficiency photovoltaics enabled by broadband diffractive optics,” Prog. Photovolt. **23**, 1073–1079 (2015). [CrossRef]

**41. **Data sheet of Heidelberg μ PG 101: http://www.himt.de/index.php/upg-101.html

**42. **Data sheet of AZ developer: http://www.microchemicals.com/micro/az_developer.pdf

**43. **B. Webb, J. Bradford, K. Lim, N. Bodnar, A. Vaupel, E. McKee, M. Baudelet, M. M. Durand, L. Shah, and M. Richardson, “Compact 10 TW laser to generate multi-filament arrays,” CLEO: 2014, OSA Technical Digest (Optical Society of America, 2014), paper SM1F.6.

**44. **B. Webb, “Design and Engineering of Ultrafast Amplification Systems,” Ph.D. thesis, University of Central Florida (2016).

**45. **K. H. Kagalwala, G. Di Giuseppe, A. F. Abouraddy, and B. E. A. Saleh, “Bell’s measure in classical optical coherence,” Nat. Photon. **7**, 72–78 (2013). [CrossRef]

**46. **C. Okoro, H. E. Kondakci, A. F. Abouraddy, and K. C. Toussaint, “Demonstration of an optical-coherence converter,” Optica **4**, 1052–1058 (2017). [CrossRef]