This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

thanks: Corresponding author: yessenov@knights.ucf.eduthanks: raddy@creol.ucf.edu

Ultra-compact synthesis of space-time wave packets

Murat Yessenov CREOL, The College of Optics & Photonics, University of Central Florida, Orlando, FL 32816, USA    Oussama Mhibik CREOL, The College of Optics & Photonics, University of Central Florida, Orlando, FL 32816, USA    Lam Mach CREOL, The College of Optics & Photonics, University of Central Florida, Orlando, FL 32816, USA    Tina M. Hayward Department of Electrical and Computer Engineering, University of Utah, Salt Lake City, Utah 84112, USA    Rajesh Menon Department of Electrical and Computer Engineering, University of Utah, Salt Lake City, Utah 84112, USA    Leonid Glebov CREOL, The College of Optics & Photonics, University of Central Florida, Orlando, FL 32816, USA    Ivan Divliansky CREOL, The College of Optics & Photonics, University of Central Florida, Orlando, FL 32816, USA    Ayman F. Abouraddy CREOL, The College of Optics & Photonics, University of Central Florida, Orlando, FL 32816, USA
Abstract

Space-time wave packets (STWPs) are pulsed fields in which a strictly prescribed association between the spatial and temporal frequencies yields surprising and useful behavior. However, STWPs to date have been synthesized using bulky free-space optical systems that require precise alignment. We describe a compact system that makes use of a novel optical component: a chirped volume Bragg grating that is rotated by 4545^{\circ} with respect to the plane-parallel device facets. By virtue of this grating’s unique structure, cascaded gratings resolve and recombine the spectrum without free-space propagation or collimation. We produce STWPs by placing a phase plate that spatially modulates the resolved spectrum between such cascaded gratings, with a device volume of 25×25×825\!\times\!25\!\times\!8 mm3, which is orders-of-magnitude smaller than previous arrangements.

preprint: APS/123-QED

Unlike conventional pulsed beams whose spectra are separable with respect to the spatial and temporal degrees of freedom, space-time wave packets (STWPs) Yessenov et al. (2022a) are structured pulsed fields in which a tight association is enforced between the spatial and temporal frequencies Turunen and Friberg (2010); Hernández-Figueroa et al. (2014). Unique features of STWPs stem from this prescribed spatio-temporal spectral structure, including propagation invariance in linear media Malaguti et al. (2008); Kondakci and Abouraddy (2017); Porras (2017), self-healing Kondakci and Abouraddy (2018), tunable group velocity Salo and Salomaa (2001); Wong and Kaminer (2017); Efremidis (2017); Kondakci and Abouraddy (2019), and anomalous refraction Bhaduri et al. (2020). Unlike previously theorized propagation-invariant focus-wave modes Brittingham (1983) and X-waves Saari and Reivelt (1997), which did not make experimental headway Turunen and Friberg (2010), STWPs are readily synthesizable without need for exorbitantly broad spectra nor large numerical apertures Yessenov et al. (2022a). Moreover, STWPs are part of a burgeoning effort investigating spatio-temporally structured fields, including flying-focus pulses Sainte-Marie et al. (2017); Froula et al. (2018); Jolly et al. (2020), wave packets with transverse angular momentum Hancock et al. (2019); Chong et al. (2020); Gui et al. (2021), toroidal pulses Wan et al. (2022); Zdagkas et al. (2022), among other examples Chen et al. (2022); Yessenov et al. (2022a). A shared feature in the synthesis of such fields is that the initial pulse spectrum is spatially resolved via a diffraction grating to spatially modulate each wavelength. Consequently, these optical arrangements are bulky and require careful alignment.

Chirped Bragg volume gratings (CBGs) offer an alternative route to spatially resolving the spectrum Glebov et al. (2014). At normal incidence, a CBG stretches the reflected pulse temporally Kaim et al. (2014), but at oblique incidence the reflected wavelengths are spatially resolved, which can be used in multiplexing and demultiplexing in optical communications systems Gerken and Miller (2003). Recently, we made use of such a component to produce linear spatial chirp without the spectral phase introduced by conventional gratings Yessenov et al. (2022b, c). Nevertheless, CBGs typically have a large footprint because the reflected field exits from the same input facet, so that increasing the bandwidth requires increasing both the length and the width of the device (e.g., CBG footprint in Yessenov et al. (2022b, c) was 37×\times25×\times6 mm3).

We recently introduced a new optical component that lifts these conventional restrictions: a rotated chirped Bragg volume grating (r-CBG), in which the volume grating structure is rotated by 4545^{\circ} with respect to the plane-parallel device facets Mhibik et al. (2022). Such a device has several salutary features: (1) the input beam is incident normally on one facet, and the spectrally resolved field exits a different facet also normally; (2) the emerging spectrum is collimated without requiring a lens and does not require further free-space propagation; (3) the link between the length and width of the device is severed, so that its cross-section can take the form of a narrow rectangle rather than a large square; and (4) a pair of r-CBGs can be cascaded to resolve and then recombine the spectrum, thus reconstituting the input pulse at the output.

Here, we construct an ultra-compact optical system for synthesizing STWPs comprising a cascaded pair of r-CBGs. The configuration is based conceptually on that in our previous work employing transmissive phase plates for spatially modulating the spectrally resolved wavefront of an ultrafast pulse produced by diffraction gratings Kondakci et al. (2018). Here, we sandwich such a phase plate between two appropriately cascaded r-CBGs to spatially modulate the spectrum with no need for any free-space propagation or collimation. This new spatio-temporal synthesis system has an ultra-compact footprint (25×25×825\times 25\times 8 mm3) whose volume is orders-of-magnitude smaller with respect to that in Kondakci et al. (2018), is easily aligned, and is significantly more robust, thus offering opportunities for producing STWPs in moving platforms.

The basic concept of an r-CBG is illustrated in Fig. 1. Traditionally, a pulse normally incident on a conventional CBG Glebov et al. (2014); Kaim et al. (2014) (in which the Bragg structure is parallel to the input facet) is temporally stretched [Fig. 1(a)]. At oblique incidence, each wavelength reflects from a different depth within the CBG, and spatial dispersion resolves the spectrum of a collimated pulse, which is reflected back obliquely from the same input facet [Fig. 1(b)]. It can be shown that the ratio of the CBG width WW and length LL is WL=2sinφno2sin2φ\tfrac{W}{L}\!=\!2\tfrac{\sin\varphi}{\sqrt{n_{\mathrm{o}}^{2}-\sin^{2}\varphi}}, where φ\varphi is the external incident angle with respect to the normal to the CBG entrance facet, non_{\mathrm{o}} is the average refractive index. Increasing the bandwidth at a fixed chirp rate requires increasing LL, which therefore implies simultaneously increasing its width WW Mhibik et al. (2022).

Refer to caption
Figure 1: Concept of an r-CBG. (a) A pulse normally incident on a CBG is temporally stretched but not spectrally resolved. (b) A pulse obliquely incident on a CBG is spectrally resolved. (c) Spectrally resolving a pulse normally incident on an r-CBG. (d) Photograph of an r-CBG resolving the spectrum of white light. The r-CBG is designed for operation at λo800\lambda_{\mathrm{o}}\!\approx\!800 nm, so operation in the visible requires oblique incidence.

The change in geometry associated with rotating the Bragg structure by 4545^{\circ} with respect to the plane-parallel device facets (r-CBG) produces a new configuration for spectral analysis. Here, the collimated pulse is incident normally on the input facet (width WW), and the resolved spectrum exits from the orthogonal facet (length LL) also normally [Fig. 1(c)]. This convenient configuration simplifies the alignment and the cascading of further optical elements (including other r-CBGs). Moreover, the width WW and length LL are now independent of each other: LL can be extended to increase the resolved bandwidth while holding WW fixed. The dimensions L×WL\times W for an r-CBG can thus be smaller than those for a CBG of the same chirp rate and bandwidth. In our experiments, we make use of r-CBGs of dimensions 25×12.5×825\times 12.5\times 8 mm3 designed for operation in the vicinity of λo800\lambda_{\mathrm{o}}\!\sim\!800 nm. The dimension 12.5 mm can be reduced to 4\sim\!4 mm without a change in the device performance. A photograph of the r-CBG is shown in Fig. 1(d), where it is obliquely oriented to operate in the visible.

A pair of r-CBGs [Fig. 2(a)] can be cascaded to resolve and then recombine the spectrum. A collimated pulse is incident normally on the first (r-CBG1) on the facet of width WW, and the spatially resolved spectrum emerges from the facet of length LL. This spectrum then impinges on the corresponding facet of the appropriately oriented second device (r-CBG2), thereby combining the spectrum and reconstituting the pulse that emerges from the facet of width WW. Physical realizations are shown in Fig. 2(b,c). We confirm that the input beam impinging on r-CBG1 and the output beam emerging from r-CBG2 have the same spatial structure [Fig. 2(d)] and the same spectrum [Fig. 2(e)]. Between r-CBG1 and r-CBG2, the spectrum has been spread over a large spatial width [Fig. 2(d,e)]. In (d-e) the input source is a mode-locked Ti:Sapp laser having 10-nm-bandwidth (FWHM), centered at a wavelength \approx800 nm, and a beam width of 1.5\sim\!1.5 mm.

Refer to caption
Figure 2: (a) Cascaded r-CBGs resolve and then recombine the spectrum. (b) Photograph of cascaded r-CBGs with a monochromatic red input laser. (c) Schematic of the cascaded configuration, where we capture (d) the intensity profile and (e) the spectrum at 3 different locations: (1) at the input, (2) after r-CBG1, and (3) after r-CBG2.
Refer to caption
Figure 3: (a) Representation of the spatio-temporal spectrum of an STWP on the surface of the light-cone at its intersection with a tilted spectral plane. The spectral projection onto (kx,ωc)(k_{x},\tfrac{\omega}{c})-domain is a conic section. (b) Time-averaged axial intensity profile I(x,z)I(x,z) of an STWP with θ=44.97\theta\!=\!44.97^{\circ} and beam size Δx40\Delta x\!\approx\!40 μ\upmum. (c) Conventional free-space setup for producing STWPs and (d) the phase plate that introduces the requisiste spatio-temporal spectral structure into the field. (e) Schematic of the ultra-compact synthesis system for STWPs based on r-CBGs, and (f) its physical realization. Note that we also show the r-CBG-based system in (c) at the same scale as the conventional setup for comparison. G: Diffraction grating; L: cylindrical lens; PP: phase plate.

Such a device can help dramatically reduce the size of the conventional spatio-temporal synthesis setup for STWPs developed over the course of the past 5\sim\!5 years Yessenov et al. (2022a). The key feature of an STWP is that a particular form of angular dispersion – a wavelength-dependent propagation angle φ(ω)\varphi(\omega) – is introduced into a generic plane-wave pulse, such that the axial wave number kzk_{z} and temporal frequency ω\omega are related via a linear relationship Ω=(kzko)ctanθ\Omega\!=\!(k_{z}-k_{\mathrm{o}})c\tan{\theta} in free space, where Ω=ωωo\Omega\!=\!\omega-\omega_{\mathrm{o}}, ωo\omega_{\mathrm{o}} is a fixed frequency associated with on-axis propagation, cc is the speed of light in vacuum, ko=ωock_{\mathrm{o}}\!=\!\tfrac{\omega_{\mathrm{o}}}{c}, and θ\theta is the spectral tilt angle Kondakci and Abouraddy (2017). If we consider only one transverse coordinate xx (and hold the field uniform along yy), then the spatio-temporal spectrum of the STWP is the conic section at the intersection of the light-cone kx2+kz2=(ωc)2k_{x}^{2}+k_{z}^{2}\!=\!(\tfrac{\omega}{c})^{2} with a plane that is parallel to the kxk_{x}-axis and makes an angle θ\theta with the kzk_{z}-axis [Fig. 3(a)]; where kxk_{x} is the transverse wave number along xx. The STWP field E(x,z;t)=ei(kozωot)ψ(x,z;t)E(x,z;t)\!=\!e^{i(k_{\mathrm{o}}z-\omega_{\mathrm{o}}t)}\psi(x,z;t) has a propagation-invariant envelope ψ(x,z;t)=ψ(x,0;tz/v~)\psi(x,z;t)\!=\!\psi(x,0;t-z/\widetilde{v}) that travels rigidly in free space at a group velocity v~=ctanθ\widetilde{v}\!=\!c\tan{\theta} determined solely by the spectral tilt angle θ\theta Kondakci and Abouraddy (2019), so the time-averaged intensity I(x,z)=𝑑t|E(x,z;t)|2=I(x,0)I(x,z)\!=\!\int\!dt\,|E(x,z;t)|^{2}\!=\!I(x,0) is independent of zz [Fig. 3(b)].

Synthesizing such an STWP requires an arrangement that (1) spatially resolves the spectrum of a generic pulse; (2) modulates the spatial distribution of each frequency ω\omega to associate with it the requisite spatial frequency kx(ω)k_{x}(\omega); before (3) recombining the spectrum to reconstitute the pulse. The first and final steps are typically realized with a diffraction grating, whereas the spatial modulation is achieved with a transmissive phase plate Kondakci et al. (2018) or reflective spatial light modulator (SLM) in a folded system Kondakci and Abouraddy (2017). A photograph of a typical synthesis system that makes use of a transmissive phase plate is shown in Fig. 3(c). The requisite phase distribution imparted by the transmissive phase plate to produce an STWP with θ=44.97\theta\!=\!44.97^{\circ} over a bandwidth Δλ15\Delta\lambda\!\approx\!15 nm is shown in Fig. 3(d). The corresponding system using a reflective SLM would cut the volume in half.

The spatio-temporal synthesis system based on r-CBGs makes use of the same strategy. However, rather than the combination of a diffraction grating and collimating lens (which occupy a large volume and require careful alignment), the resolved spectrum is formed at the output facet of r-CBG1, the phase plate abuts that facet, and r-CBG2 is placed immediately after the phase plate. In other words, the two r-CBGs sandwich the phase plate as shown in Fig. 3(e). The two r-CBGs are designed to yield a linear spectral chirp of 0.7\approx\!0.7 nm/mm over a length of L=25L\!=\!25 mm. To bring about clearly the dramatic reduction in footprint of this synthesis arrangement based on r-CBGs, we place it alongside the conventional setup in Fig. 3(c); a magnified view is provided in Fig. 3(f). The total volume of the new system is only 25×25×825\times 25\times 8 mm3.

Refer to caption
Figure 4: Layout of the setup for the ultra-compact synthesis of STWPs via a pair of r-CBGs and a phase plate, along with the characterization setup.
Refer to caption
Figure 5: Diffraction-free propagation and spatio-temporal spectrum of the STWPs synthesized with the new system. (a) Subluminal STWP with θ=44.97\theta\!=\!44.97^{\circ} and Δx40\Delta x\!\approx\!40 μ\upmum; and (b) superluminal STWP with θ=46\theta\!=\!46^{\circ} and Δx13\Delta x\!\approx\!13 μ\upmum. The columns show the following quantities: (i) the spatio-temporal spectral intensity projected onto the (kx,λ)(k_{x},\lambda)-plane; (ii) the time-averaged intensity in the transverse (x,y)(x,y)-plane at z=0z\!=\!0; (iii) the time-avaraged intensity I(x,z)I(x,z), and (iv) I(x,0)I(x,0).

A detailed layout of the setup for the ultra-compact synthesis and the characterization of STWPs is depicted in Fig. 4. The pulses from Ti:Sapphire laser mentioned above are normally incident on the input facet of r-CBG1, and the spatially resolved spectrum exits r-CBG1 normally to the output facet along the zz-axis. We place the transmissive phase plate [Fig. 3(d)] abutting this facet to introduce the requisite spatio-temporal spectral structure associated with STWPs [Fig. 3(a)]. The phase plate is fabricated using gray-scale lithography Mohammad et al. (2017) and was designed to synthesize STWPs with tilt angle θ=46\theta\!=\!46^{\circ} (superluminal STWP) over a bandwidth of Δλ17\Delta\lambda\!\approx\!17 nm when oriented such that the longest wavelength is associated with kx=0k_{x}\!=\!0. When the phase plate is flipped horizontally so that kx=0k_{x}\!=\!0 is associated with the shortest wavelength in the spectrum, an STWP with θ=44\theta\!=\!44^{\circ} (subluminal STWP) is produced. Placing r-CBG2 immediately after the phase plate then reconstitutes the spectrum, and the beam exiting r-CBG2 is an STWP with the desired spatio-temporal spectral structure. We place a telescope system in the path of the beam after r-CBG2. For the superluminal STWP of spectral tilt angle θ=46\theta\!=\!46^{\circ} in Fig. 5(b), the telescope system images the output plane of r-CBG2 with unity magnification (the focal lengths of the two lenses are f1=f2=5f_{1}\!=\!f_{2}\!=5 cm). For the subluminal STWP of spectral tilt angle θ=44\theta\!=\!44^{\circ} in Fig. 5(a), the telescope introduces 3×3\times magnification (f2=3f1=15f_{2}\!=\!3f_{1}\!=15 cm), which changes the spectral tilt angle to θ=44.97\theta\!=\!44.97^{\circ} and yields a larger beam width.

We characterize the synthesized STWP in the physical domain by scanning a CCD camera (C1; TheImagingSource, DMK 27BUP031) along the propagation axis zz and monitor the axial evolution of the transverse spatial profile I(x,z)I(x,z) [Fig. 5, columns (ii-iv)]. To confirm that the desired spatio-temporal spectral correlation between kxk_{x} and λ\lambda has been realized by the r-CBG-based setup, we characterize the STWPs in the spectral domain by measuring spatio-temporal spectrum |ψ~(kx,λ)|2|\widetilde{\psi}(k_{x},\lambda)|^{2} projected onto the (kx,λ)(k_{x},\lambda)-plane [Fig. 5, column (i)]. For that, we exploit the finite diffraction efficiency of r-CBG2 and use the transmitted spatially resolved spectrum on which the phase plate has imprinted a spatio-temporal spectral structure. A CCD camera (C2; TheImagingSource, DMK 33UX178) placed in the Fourier plane of a cylindrical lens Lx5L_{x-5} captures the spatio-temporal spectrum. A telescope (focal lengths f4=f3/2.5=10f_{4}\!=\!f_{3}/2.5\!=10 cm) placed in the path of the spectrum to de-magnifies the field by 2.5×2.5\times to fit the spatio-temporal spectrum into the active area of C2C_{2}.

The diffraction-free characteristics of the produced STWPs are confirmed in Fig. 5. We first plot the spatio-temporal spectra projected onto the (kx,λ)(k_{x},\lambda)-plane for the subluminal [Fig. 5(a), column (i)] and superluminal [Fig. 5(a), column (i)] STWPs. Next we plot the time-averaged intensity profiles in the transverse plane at z=0z\!=\!0 [Fig. 5, column (ii)] and along the optical axis I(x,z)I(x,z) [Fig. 5, column (iii)]. The propagation length in Fig. 5(a) exceeds the Rayleigh range of a Gaussian beam with the same width Δx\Delta x by 50×50\times, and in Fig. 5(b) by 90×90\times.

We have drastically reduced the size of the setup needed for synthesizing STWPs by eliminating the need for diffraction gratings and collimating lenses. By utilizing instead r-CBGs, the STWP-synthesis system is compact (25×25×8\sim\!25\times 25\times 8 mm3), easily aligned, and robust. These newly developed r-CBGs can be similarly used to reduce the size of the synthesis arrangement for other spatio-temporally structured pulses. It remains an open question whether the system size can be further reduced by relying on a single metasurface or a pair of metasurfaces.

In conclusion, we have constructed an ultra-compact setup for synthesizing STWPs based on an alternative strategy for spectrally resolving the spectrum of an optical pulse. Rather than using conventional diffraction gratings, we employ chirped Bragg volume grating in which the Bragg structure is rotated by 4545^{\circ} with respect to the plane-parallel facets of the device. Because the spectrum is resolved without further need of free-space propagation or collimation by a lens, and because the spatially resolved spectrum exits normally from a facet orthogonal to that of the input, we can sandwich a phase plate that spatially modulates the spectrally resolved wavefront between a cascade of two r-CBGs. The STWPs produced using this ultra-compact configuration have the same characteristics as those produced by setups based on diffraction gratings, whose volume exceeds that of the r-CBG-based setup by several orders of magnitude.

Funding: U.S. Office of Naval Research (ONR) N00014-17-1-2458 and N00014-20-1-2789.

Disclosures: The authors declare no conflicts of interest.

References

  • Yessenov et al. (2022a) M. Yessenov, A. A. Hall, K. L. Schepler,  and A. F. Abouraddy, Adv. Opt. Photon. 14, 455 (2022a).
  • Turunen and Friberg (2010) J. Turunen and A. T. Friberg, Prog. Opt. 54, 1 (2010).
  • Hernández-Figueroa et al. (2014) H. E. Hernández-Figueroa, E. Recami,  and M. Zamboni-Rached, eds., Non-diffracting Waves (Wiley-VCH, 2014).
  • Malaguti et al. (2008) S. Malaguti, G. Bellanca,  and S. Trillo, Opt. Lett. 33, 1117 (2008).
  • Kondakci and Abouraddy (2017) H. E. Kondakci and A. F. Abouraddy, Nat. Photon. 11, 733 (2017).
  • Porras (2017) M. A. Porras, Opt. Lett. 42, 4679 (2017).
  • Kondakci and Abouraddy (2018) H. E. Kondakci and A. F. Abouraddy, Opt. Lett. 43, 3830 (2018).
  • Salo and Salomaa (2001) J. Salo and M. M. Salomaa, J. Opt. A 3, 366 (2001).
  • Wong and Kaminer (2017) L. J. Wong and I. Kaminer, ACS Photon. 4, 2257 (2017).
  • Efremidis (2017) N. K. Efremidis, Opt. Lett. 42, 5038 (2017).
  • Kondakci and Abouraddy (2019) H. E. Kondakci and A. F. Abouraddy, Nat. Commun. 10, 929 (2019).
  • Bhaduri et al. (2020) B. Bhaduri, M. Yessenov,  and A. F. Abouraddy, Nat. Photon. 14, 416 (2020).
  • Brittingham (1983) J. N. Brittingham, J. Appl. Phys. 54, 1179 (1983).
  • Saari and Reivelt (1997) P. Saari and K. Reivelt, Phys. Rev. Lett. 79, 4135 (1997).
  • Sainte-Marie et al. (2017) A. Sainte-Marie, O. Gobert,  and F. Quéré, Optica 4, 1298 (2017).
  • Froula et al. (2018) D. H. Froula, D. Turnbull, A. S. Davies, T. J. Kessler, D. Haberberger, J. P. Palastro, S.-W. Bahk, I. A. Begishev, R. Boni, S. Bucht, J. Katz,  and J. L. Shaw, Nat. Photon. 12, 262 (2018).
  • Jolly et al. (2020) S. W. Jolly, O. Gobert, A. Jeandet,  and F. Quéré, Opt. Express 28, 4888 (2020).
  • Hancock et al. (2019) S. W. Hancock, S. Zahedpour, A. Goffin,  and H. M. Milchberg, Optica 6, 1547 (2019).
  • Chong et al. (2020) A. Chong, C. Wan, J. Chen,  and Q. Zhan, Nat. Photon. 14, 350 (2020).
  • Gui et al. (2021) G. Gui, N. J. Brooks, H. C. Kapteyn, M. M. Murnane,  and C.-T. Liao, Nat. Photon. 15, 608 (2021).
  • Wan et al. (2022) C. Wan, Q. Cao, J. Chen, A. Chong,  and Q. Zhan, Nat. Photon. 16, 519 (2022).
  • Zdagkas et al. (2022) A. Zdagkas, C. McDonnell, J. Deng, Y. Shen, G. Li, T. Ellenbogen, N. Papasimakis,  and N. I. Zheludev, Nat. Photon. 16, 523 (2022).
  • Chen et al. (2022) L. Chen, W. Zhu, P. Huo, J. Song, H. J. Lezec, T. Xu,  and A. Agrawal, Sci. Adv. 8, eabq8314 (2022).
  • Glebov et al. (2014) L. Glebov, V. Smirnov, E. Rotari, I. Cohanoschi, L. Glebova, O. Smolski, J. Lumeau, C. Lantigua,  and A. Glebov, Opt. Eng. 53, 051514 (2014).
  • Kaim et al. (2014) S. Kaim, S. Mokhov, B. Y. Zeldovich,  and L. B. Glebov, Opt. Eng. 53, 051509 (2014).
  • Gerken and Miller (2003) M. Gerken and D. A. B. Miller, IEEE Photon. Technol. Lett. 15, 1097 (2003).
  • Yessenov et al. (2022b) M. Yessenov, J. Free, Z. Chen, E. G. Johnson, M. P. J. Lavery, M. A. Alonso,  and A. F. Abouraddy, Nat. Commun. 13, 4573 (2022b).
  • Yessenov et al. (2022c) M. Yessenov, , Z. Chen, M. P. J. Lavery,  and A. F. Abouraddy, Opt. Lett. 47, 4131 (2022c).
  • Mhibik et al. (2022) O. Mhibik, M. Yessenov, L. Mach, L. Glebov, A. F. Abouraddy,  and I. Divliansky, arXiv 2212.09871 (2022).
  • Kondakci et al. (2018) H. E. Kondakci, M. Yessenov, M. Meem, D. Reyes, D. Thul, S. R. Fairchild, M. Richardson, R. Menon,  and A. F. Abouraddy, Opt. Express 26, 13628 (2018).
  • Mohammad et al. (2017) N. Mohammad, M. Meem, X. Wan,  and R. Menon, Scientific Reports 7, 1 (2017).