Refraction of space-time wave packets: III. Experiments at oblique incidence
Abstract
The refraction of space-time (ST) wave packets at planar interfaces between non-dispersive, homogeneous, isotropic dielectrics exhibit fascinating phenomena, even at normal incidence. Examples of such refractive phenomena include group-velocity invariance across the interface, anomalous refraction, and group-velocity inversion. Crucial differences emerge at oblique incidence with respect to the results established at normal incidence. For example, the group velocity of the refracted ST wave packet can be tuned simply by changing the angle of incidence. In paper (III) of this sequence, we present experimental verification of the refractive phenomena exhibited by ST wave packets at oblique incidence that were predicted in paper (I). We also examine a proposal for ‘blind synchronization’ whereby identical ST wave packets arrive simultaneously at different receivers without a priori knowledge of their locations except that they are all located at the same depth beyond an interface between two media. A first proof-of-principle experimental demonstration of this effect is provided.
I Introduction
Snell’s law governs the change in the propagation direction of a monochromatic plane wave incident onto a planar interface between two optical media Saleh and Teich (2019). If and are the refractive indices of the two media, and and are the propagation angles with respect to the normal to the interface for the incident and refracted waves, respectively, then Snell’s law dictates that . Although this result, strictly speaking, applies to only monochromatic plane waves, nevertheless its utility is typically extended in practice to conventional pulsed beams where it can provide an adequate approximation, especially for narrowband paraxial fields in which the spatial and temporal degrees of freedom (DoFs) are uncoupled. Crucially, this entails that the group velocity of the transmitted wave packet depends solely on the local optical properties of the second medium, and is independent of the incident angle. In other words, no ‘memory’ of the incident wave packet is retained as far as the group velocity of the transmitted wave packet is concerned.
In contrast to conventional wave packets, space-time (ST) wave packets Kondakci and Abouraddy (2016); Parker and Alonso (2016) are a unique class of pulsed optical beams in which the spatial and temporal DoFs are non-separable, and are instead inextricably intertwined Reivelt and Saari (2003); Kiselev (2007); Turunen and Friberg (2010); Hernández-Figueroa et al. (2014). ST wave packets are endowed with a precise spatio-temporal structure in which each spatial frequency is tightly associated with a single wavelength (or temporal frequency) Donnelly and Ziolkowski (1993); Saari and Reivelt (2004); Longhi (2004); Kondakci et al. (2019); Hall et al. (2021). Uniquely, this spatio-temporal structure determines the group velocity of the ST wave packet Salo and Salomaa (2001); Recami et al. (2003); Valtna et al. (2007); Zamboni-Rached and Recami (2008); Kondakci and Abouraddy (2017); Efremidis (2017); Wong and Kaminer (2017); Kondakci and Abouraddy (2019); Bhaduri et al. (2019a). Consequently, the rearrangement of the field structure upon refraction at a planar interface leads to a change in the group velocity of the transmitted wave packet that depends on the group velocity of the incident wave packet. Because of this, ST wave packets – in contrast to conventional wave packets – retain a ‘memory’ of the incident wave packet, which is the basis for the fascinating refractive phenomena exhibited by ST wave packets at normal and oblique incidence Bhaduri et al. (2020).
In paper (I) of this series, we presented a theoretical study of the refraction of ST wave packets at normal and oblique incidence on a planar interface between two non-dispersive, homogeneous, isotropic dielectrics Yessenov et al. (2021). At normal incidence, the above-described memory effect is the basis for group-velocity invariance, anomalous refraction, and group-velocity inversion Bhaduri et al. (2020). These phenomena extend to oblique incidence, although the conditions for their realization change with incident angle. Crucially, a new refractive phenomenon emerges at oblique incidence: the group velocity of the transmitted wave packet changes with the incident angle when all else is held fixed. Indeed, the group velocity increases with incident angle upon refraction into a higher-index medium when the incident wave packet is subluminal, and it decreases when the wave packet is superluminal Bhaduri et al. (2020); Yessenov et al. (2021).
In paper (II), we provided experimental confirmation of the predicted phenomena at normal incidence Allende Motz et al. (2021a). Here we examine experimentally in detail the refraction of ST wave packets at oblique incidence at a planar interface between two non-dispersive, homogeneous, isotropic dielectrics using the interferometric measurement strategy described in paper (II). After a brief overview of the basic law of refraction governing baseband ST wave packets at oblique incidence, we confirm the predicted dependence at oblique incidence of the group velocity for the transmitted wave packet on that of the incident wave packet and on the incidence angle. We then verify the predicted changes that occur at oblique incidence in the conditions required to realize group-velocity invariance and group-velocity inversion. Finally, we examine a proposal that makes use of these unique characteristics to realize blind synchronization of multiple receivers via pulses emitted from a single transmitter with no a priori knowledge of the receiver positions except that they are situated at the same depth below an interface between two media. Building on our recent observation of isochronous ST wave packets that traverse a planar slab at a fixed group delay independently of the angle of incidence (despite the change in distance traversed) Allende Motz et al. (2021b), we report here a proof-of-principle realization of the blind synchronization scheme.


II Law of refraction for space-time wave packets at oblique incidence
A ST wave packet is a pulsed beam having finite spatial and temporal bandwidths in which each spatial frequency is associated with a single temporal frequency ; here and are the transverse and axial coordinates, respectively, and and are the corresponding components of the wave vector. We identify the propagation direction with the -axis, we assume a spatial spectrum that is symmetric around , and we assign the temporal frequency to the spatial frequency . The spectral support domain for a ‘baseband’ ST wave packet Yessenov et al. (2019a) in a non-dispersive medium of index lies at the intersection of the light-cone with a tilted plane that is parallel to the -axis and makes an angle (the spectral tilt angle) with respect to the -axis, where and is the speed of light in vacuum Yessenov et al. (2019b, a, c). This construction results in a propagation-invariant wave packet transported rigidly at a group velocity , where is the group index. In other words, the group velocity is determined by an internal DoF (the spectral tilt angle ) that is related to the spatio-temporal field structure. Note that the ST wave packet reverts to a plane-wave pulse when . Within this framework, is either the minimum temporal frequency when the wave packet is superluminal (), or the maximum temporal fewquency when the wave packet is subluminal () Yessenov et al. (2019a).
When the ST wave packet is normally incident on the interface, the wave vector associated with is perpendicular to the interface [Fig. 1(a)], whereas all other temporal frequencies are not. Two quantities are conserved across the interface: the transverse momentum and the temporal frequency . If the refractive indices of the media are and , and the group indices of the incident and transmitted wave packets are and , respectively, then it can be shown that:
(1) |
This law of refraction for ST wave packets was derived in paper (I) Yessenov et al. (2021) and explored experimentally in paper (II) Allende Motz et al. (2021a); see also Ref. Bhaduri et al. (2020). The quantity is related to the curvature of the spatio-temporal spectrum when projected onto the -plane. We have denoted this newly identified refractive invariant the ‘spectral curvature’. The law of refraction in Eq. 1 thus expresses the invariance of the spectral curvature across the planar interface at normal incidence.
At oblique incidence, the wave vector associated with is incident at an angle with respect to the normal to the interface [Fig. 1(b)]. The group indices of the incident and transmitted ST wave packets satisfy a modified law Bhaduri et al. (2020); Yessenov et al. (2021):
(2) |
where is obtained from Snell’s law . Equation 2 indicates the invariance of a new spectral curvature at oblique incidence. The two curves corresponding to Eq. 1 and Eq. 2 intersect at the luminal point and , where the ST wave packets revert to plane-wave pulses [Fig. 3], and Eq. 2 reverts back to Eq. 1 at normal incidence .
The rationale for the modified oblique-incidence spectral curvature can be elucidated by reference to Fig. 1(b) and Fig. 2. The spectral support domain of the incident ST wave packet in the coordinate system – which is aligned with the propagation direction of the wave packet – lies at the intersection of the light-cone with the plane [Fig. 2(a)]. At normal incidence, the projection of this spectrum onto the -plane is invariant across the interface Bhaduri et al. (2019a, 2020). At oblique incidence, the transverse wave number is not parallel to the interface, and therefore is not conserved. However, in the coordinate system that is aligned with the interface rather than with the propagation direction of the ST wave packet, the transverse wave number is conserved. In this coordinate system, the spectral support domain is rotated an angle around the -axis [Fig. 2(b)], and the projection of this new spectrum onto the -plane is invariant across the interface [Fig. 2(c)]. It can be shown that is invariant to first order (in the small-angle approximation ), and the new invariant spectral curvature at oblique incidence is . Because the opening angle of the light-cone changes from in the first medium to in the second, the invariance of the spectral projection onto the -plane produces a new spectral support domain for the transmitted wave packet on the surface of the light-cone [Fig. 2(d)]. A rotation through an angle (rather than ) around the -axis returns the spectral support domain to the coordinate system aligned with the propagation direction of the transmitted wave packet [Fig. 2(e)].
III Confirming the law of refraction at oblique incidence

To verify the law of refraction in Eq. 2 for baseband ST wave packets at oblique incidence, we make use of the interferometric procedure outlined in paper (II); see also Kondakci and Abouraddy (2019); Bhaduri et al. (2019a). The planar slabs used in our measurements are placed in the common path of the co-aligned ST wave packet and the reference pulse, and we tilt the sample with respect to the incident wave packets. We plot the measurement results in Fig. 3. Here we fix the angle of incidence and sweep the spectral tilt angle of the ST wave packet incident from free space onto a 5-mm-thick layer of sapphire. For each value of , we measure the group delay incurred by the wave packet across the layer, from which we estimate . The measurements are carried out at normal incidence and at oblique incidence . The measurements verify for the first time the change in the law of refraction for ST wave packets with respect to that at normal incidence, which was demonstrated in Bhaduri et al. (2020). It is clear that the two curves for normal and oblique incidence in Fig. 3 intersect at the point corresponding to the luminal condition . The curve for oblique incidence is shifted below that for normal incidence in the superluminal regime ( and ) and above it in the subluminal regime ( and ). We next explore the far-reaching ramifications of this change.
IV Dependence of the group velocity on the angle of incidence

Figure 3 yields the relationship between the group indices and for the incident and refracted ST wave packets, respectively, at a fixed incident angle . We next examine the change in the group index for the transmitted wave packet as we vary the incidence angle while holding fixed the group index for the incident wave packet. In Fig. 4(a), we plot the calculated while varying both and for incidence from free space onto sapphire . This calculation highlights the above-mentioned ‘memory’ effect. The group index of the transmitted wave packet depends not only the refractive index of the second medium , but also on the characteristics of the incident wave packet: its group index and incident angle . The luminal condition for the incident ST wave packet is , in which case this wave packet is simply a plane-wave pulse, and the refracted wave packet is also a luminal plane-wave pulse . For subluminal incidence ST wave packets , decreases with ; i.e., the group velocity of the transmitted subluminal wave packet increases with in the higher-index medium. Conversely, for superluminal incident ST wave packets , increases with ; i.e., the group velocity of the transmitted superluminal wave packet increases with in the higher-index medium. In the case of a conventional wave packet whose spatial and temporal DoFs are separable, is independent of .
This behavior can be understood on the basis of Eq. 2. Substituting for in terms of from Snell’s law yields:
(3) |
where the factor depends on the indices and through
(4) |
The behavior of depends on the ratio . When , drops monotonically from its initial value to at . Because the term is positive in the subluminal regime , the decrease of with results in a decrease of with . In the superluminal regimes, the term is negative, and increase with . On the other hand, when , increases monotonically with and reaches a singularity when , corresponding to the critical angle (total internal reflection) at the interface. Consequently, the opposite trends for with ensue.
Our measurements confirm this predicted behavior. In Fig. 4(b,c) we plot the change in the refracted group index with with respect to that at normal incidence, , for incidence from free space onto two materials: MgF2 in Fig. 4(b) and sapphire in Fig. 4(c). For each material, we carry out measurements with ST wave packets synthesized in free space in the subluminal () and superluminal () regimes. We note several general observations about the results that are expected from our analysis. First, deviates monotonically from the normal-incidence value . Second, the group index drops at oblique incidence in the subluminal regime, and increases in the superluminal regime. That is, the group velocity of a refracted subluminal ST wave packet increases with incident angle in the higher-index medium. Conversely, the group velocity of a refracted superluminal ST wave packet decreases with . These features are critical for the blind synchronization scheme we investigate below.
V Impact of the angle of incidence on group-velocity invariance and inversion
The law of refraction at normal incidence (Eq. 1) predicts three phenomena that occur for any pair of media: (1) group velocity invariance , which occurs when ; (2) anomalous refraction whereby when , which occurs when ; and (3) group-velocity inversion , which occurs when ; see Fig. 3. These three phenomena were confirmed experimentally in Ref. Bhaduri et al. (2020) and in more detail in paper (II) Allende Motz et al. (2021a). Despite the difference between the law of refraction at oblique incidence (Eq. 2) from that at normal incidence (Eq. 1), these three phenomena are still realizable, albeit with modifications to the conditions to be satisfied.
First, group-velocity invariance occurs at oblique incidence when Yessenov et al. (2021):
(5) |
where . We plot in Fig. 5 measurements of the group delay incurred upon traversing equal lengths ( mm) of free space () and MgF2 (). At normal incidence corresponding to [Fig. 5(a)], and at oblique incidence we have , corresponding to [Fig. 5(b)]. In both cases the ST wave packets accrue approximately the same delay in equal lengths of free space and MgF2. Whereas anomalous refraction occurs at normal incidence when , this regime is expanded for oblique incidence to at .


Second, group-velocity inversion occurs at oblique incidence when Yessenov et al. (2021):
(6) |
where when . We plot in Fig. 6 the measured group delay incurred upon traversing equal lengths ( mm) of sapphire and MgF2 separately, and after traversing a bilayer of the two media. The required group index for a ST wave packet in free space that experiences group-velocity invariance in sapphire and MgF2 is (); see Fig. 6(a). At oblique incidence (), the required free-space group index to realize group-velocity inversion is (); see Fig. 6(b). In both cases, the group delays measured in the individual layers are approximately equal in magnitude but opposite in sign. Consequently, the group delay in the bilayer almost vanishes as shown.
VI Blind synchronization
VI.1 Concept of blind synchronization using ST wave packets

Consider the configuration shown in Fig. 7. A transmitter (Tx) in the first medium at a height above the interface directs an ST wave packet at different angles to receivers (Rx) at different locations but at the same depth beneath the interface in the second medium. The optical delay in the first medium from Tx to the interface is , where is the group delay at normal incidence. The group velocity of the ST wave packet is the same in all directions independently of because the medium is homogeneous and isotropic. Therefore, the group delay increases with . The group delay in the second medium is , where is the group index of the transmitted wave packet, determined from Eq. 2. In contrast to the incident wave packet, the group velocity of the transmitted wave packet depends on .
The relative group delay between obliquely and normally incident ST wave packets in the first medium is:
(7) |
and in the second medium is:
(8) |
We define the total relative group delay , where is the total group delay incurred from Tx to Rx across the interface. Blind synchronization requires independently of , such that the ST wave packet arrives simultaneously at all receivers at the depth despite the different physical distances covered along the different trajectories.
We first ascertain whether realizing is physically feasible. Note that in the first medium is always positive, because the group delay along the longer oblique path is always larger than that along the shorter normal path. Blind synchronization therefore requires that in the second medium be negative. That is, the group velocity of the transmitted wave packet must increase at oblique incidence to compensate for the longer propagation distances in both media. In other words, a necessary (but not sufficient) condition is that . From the analysis and measurements reported above, this scenario occurs for subluminal ST wave packets when . We next consider whether the blind synchronization condition can be met quantitatively.

VI.2 Calculated group delay of ST wave packets between transmitter and receivers
The total relative group delay is given by:
(9) | |||||
where we made use of the law of refraction in Eq. 2, and we normalized with respect to . We plot in Fig. 8 calculations for as a function of the incident angle and the group index of the incident wave packet in two different configurations. In the first configuration we have , , and ; and in the second configuration we have , , and . In both cases we find the contour occurs at a specific group index over a range of incident angles. Such ST wave packets arrive simultaneously at receivers located at the same depth within this angular range despite following paths of very different lengths, thereby realizing blind synchronization. We find in general that meeting the requirement for blind synchronization is more favorable when and/or increase.
For the sake of comparison, we plot in Fig. 9(a) the relative group delay for ST wave packets having different group indices alongside that for a conventional wave packet whose spatial and temporal DoFs are separable. The conventional wave packet propagates in the first medium at a group velocity and in the second at (assuming non-dispersive media). The relative group delay is given by:
(10) |
The relative group delay is always positive for conventional pulses because the pulses traveling longer distances inevitably accrues larger group delays. In contrast, varying for ST wave packets changes from positive values (similarly to a conventional pulse where longer paths incur longer delays) to anomalously negative values (longer paths incur shorter delays) passing through the desired target of blind synchronization . In Fig. 9(b) we highlight this transition from to as is changed in the vicinity of .


VI.3 Blind synchronization
We plot in Fig. 10 results for a proof-of-principle test of the blind synchronization effect. The overall system configuration is illustrated in Fig. 10(a). We consider the configuration in Fig. 8(a) for propagation from free space () to sapphire () with the propagation distances at normal incidence having the ratio ; we use here mm and mm. We measure the group delay of the ST wave packet in free space [Fig. 10(b)] and in the sapphire layer [Fig. 10(c)] separately. The separate delays for the oblique-incidence paths are measured by tilting the sapphire layer and in free space by extending the distance over which the delay is measured.
We synthesize a subluminal ST wave packet in free space with (). According to Fig. 8(a), blind synchronization is realized by such a wave packet over a range of incidence angles . We plot in Fig. 10(b) the measured group delay in free space over distances extending from mm (corresponding to ) to mm (corresponding to ), with . We next plot in Fig. 10(c) the measured group delay in the sapphire layer as we increase the angle of incidence. Here, the group delay decreases with the increases distance traveled in the layer at larger . As shown in Fig. 4, this is a general feature at oblique incidence for subluminal ST wave packets when . For this particular ST wave packet with , the decrease in delay in the layer (increase in group velocity) with is sufficient to counterbalance the increase in delay in free space and the increase in propagation distance in the sapphire layer, such that the total delays in free space and the layer are approximately equal [Fig. 10(d)] for all trajectories corresponding to incident angles from to .
VII Discussion and Conclusions
The proof-of-principle experimental demonstration of blind synchronization described here extends our recent report on isochronous ST wave packets Allende Motz et al. (2021b). These are wave packets that traverse a material layer after accruing a fixed group delay independently of the angle of incidence, despite traveling larger distances in the layer (, where is the layer thickness). This can be seen as a special case of blind synchronization where we set . An increase in group velocity with incident angle for isochronous ST wave packets is required to compensate for the increase in propagation distance within the layer. This is achieved, as we have done here for blind synchronization, by employing subluminal incident ST wave packets. However, the blind synchronization scenario is more stringent because the necessary increase in group velocity with incident angle is larger: this change must counterbalance the increase in propagation distance within both media. As such, an incident ST wave packet with larger (deeper within the subluminal regime) is required for blind synchronization compared to isochronous ST wave packets. Pursuing the concept of blind synchronization further requires increasing the propagation distance over which this effect is realized by reducing the so-called ‘spectral uncertainty’ Yessenov et al. (2019b); Bhaduri et al. (2019b); Kondakci et al. (2019), which is the unavoidable ‘fuzziness’ in the association between spatial and temporal frequencies in the wave packet spectrum. In general, reducing the spectral uncertainty to increase the propagation distance requires increasing the system’s numerical aperture.
The work reported in this paper sequence, in addition to Bhaduri et al. (2020), has examined the refraction of ST wave packets in the simplest scenario: at a planar interface between two non-dispersive, homogeneous, isotropic dielectrics. Future work will be directed to an exploration of other important configurations. Paramount amongst these are the cases of dispersive and anisotropic media, which will be critical as a prelude to studying nonlinear interactions involving ST wave packets. Additionally, it will be useful to study the refraction of ST wave packets in which both transverse dimensions are included (ST needles rather than ST sheets), and ST wave packets in which a parameter is controlled along the propagation axis, such as wave packets that accelerate or decelerate Yessenov and Abouraddy (2020), or are endowed with axial spectral encoding Allende Motz et al. (2021c).
In conclusion, we have verified experimentally the law of refraction for ST wave packets obliquely incident on a planar interface between two non-dispersive, homogeneous, isotropic dielectrics (Eq. 2). This law was verified first by changing the group velocity of the incident ST wave packet at a fixed angle of incidence, and, second, by changing the angle of incidence while holding the group velocity of the incident wave packet fixed. Our measurements reveal that the group velocity of the transmitted wave packet increases with the incident angle in the subluminal regime (when the index of the second medium is higher than that of the first); conversely, the group velocity of the transmitted wave packet decreases with incident angle in the superluminal regime. Furthermore, we observed the shift that occurs in the conditions for group-velocity invariance and group-velocity inversion at oblique incidence with respect to those at normal incidence. These findings made possible the first observation of blind synchronization, whereby ST wave packets produced from a fixed source (Tx) are arranged to arrive simultaneously at different receivers (Rx) at a priori unknown positions, except that they are located at a fixed depth beyond an interface.
Funding
U.S. Office of Naval Research (ONR) contract N00014-17-1-2458.
Disclosures. The authors declare no conflicts of interest.
References
- Saleh and Teich (2019) B. E. A. Saleh and M. C. Teich, Fundamentals of Photonics (Wiley, 2019).
- Kondakci and Abouraddy (2016) H. E. Kondakci and A. F. Abouraddy, Opt. Express 24, 28659 (2016).
- Parker and Alonso (2016) K. J. Parker and M. A. Alonso, Opt. Express 24, 28669 (2016).
- Reivelt and Saari (2003) K. Reivelt and P. Saari, arxiv:physics/0309079 (2003).
- Kiselev (2007) A. P. Kiselev, Opt. Spectrosc. 102, 603 (2007).
- 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).
- Donnelly and Ziolkowski (1993) R. Donnelly and R. W. Ziolkowski, Proc. R. Soc. Lond. A 440, 541 (1993).
- Saari and Reivelt (2004) P. Saari and K. Reivelt, Phys. Rev. E 69, 036612 (2004).
- Longhi (2004) S. Longhi, Opt. Express 12, 935 (2004).
- Kondakci et al. (2019) H. E. Kondakci, M. A. Alonso, and A. F. Abouraddy, Opt. Lett. 44, 2645 (2019).
- Hall et al. (2021) L. A. Hall, M. Yessenov, and A. F. Abouraddy, Opt. Lett. 46, 1672 (2021).
- Salo and Salomaa (2001) J. Salo and M. M. Salomaa, J. Opt. A 3, 366 (2001).
- Recami et al. (2003) E. Recami, M. Zamboni-Rached, K. Z. Nóbrega, and C. A. Dartora, IEEE J. Sel. Top. Quantum Electron. 9, 59 (2003).
- Valtna et al. (2007) H. Valtna, K. Reivelt, and P. Saari, Opt. Commun. 278, 1 (2007).
- Zamboni-Rached and Recami (2008) M. Zamboni-Rached and E. Recami, Phys. Rev. A 77, 033824 (2008).
- Kondakci and Abouraddy (2017) H. E. Kondakci and A. F. Abouraddy, Nat. Photon. 11, 733 (2017).
- Efremidis (2017) N. K. Efremidis, Opt. Lett. 42, 5038 (2017).
- Wong and Kaminer (2017) L. J. Wong and I. Kaminer, ACS Photon. 4, 2257 (2017).
- Kondakci and Abouraddy (2019) H. E. Kondakci and A. F. Abouraddy, Nat. Commun. 10, 929 (2019).
- Bhaduri et al. (2019a) B. Bhaduri, M. Yessenov, and A. F. Abouraddy, Optica 6, 139 (2019a).
- Bhaduri et al. (2020) B. Bhaduri, M. Yessenov, and A. F. Abouraddy, Nat. Photon. 14, 416 (2020).
- Yessenov et al. (2021) M. Yessenov, B. Bhaduri, and A. F. Abouraddy, accompanying paper (2021).
- Allende Motz et al. (2021a) A. M. Allende Motz, M. Yessenov, B. Bhaduri, and A. F. Abouraddy, accompanying paper (2021a).
- Allende Motz et al. (2021b) A. M. Allende Motz, M. Yessenov, and A. F. Abouraddy, arXiv:2102.10505 (2021b).
- Yessenov et al. (2019a) M. Yessenov, B. Bhaduri, H. E. Kondakci, and A. F. Abouraddy, Phys. Rev. A 99, 023856 (2019a).
- Yessenov et al. (2019b) M. Yessenov, B. Bhaduri, L. Mach, D. Mardani, H. E. Kondakci, M. A. Alonso, G. A. Atia, and A. F. Abouraddy, Opt. Express 27, 12443 (2019b).
- Yessenov et al. (2019c) M. Yessenov, B. Bhaduri, H. E. Kondakci, and A. F. Abouraddy, Opt. Photon. News 30, 34 (2019c).
- Bhaduri et al. (2019b) B. Bhaduri, M. Yessenov, D. Reyes, J. Pena, M. Meem, S. R. Fairchild, R. Menon, M. C. Richardson, and A. F. Abouraddy, Opt. Lett. 44, 2073 (2019b).
- Yessenov and Abouraddy (2020) M. Yessenov and A. F. Abouraddy, Phys. Rev. Lett. 125, 233901 (2020).
- Allende Motz et al. (2021c) A. M. Allende Motz, M. Yessenov, and A. F. Abouraddy, Phys. Rev. Appl. 15, 024067 (2021c).