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: [email protected]

Refraction of space-time wave packets: III. Experiments at oblique incidence

Murat Yessenov CREOL, The College of Optics & Photonics, University of Central Florida, Orlando, FL 32816, USA    Alyssa M. Allende Motz CREOL, The College of Optics & Photonics, University of Central Florida, Orlando, FL 32816, USA    Basanta Bhaduri 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

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.

preprint: APS/123-QED

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 n1n_{1} and n2n_{2} are the refractive indices of the two media, and ϕ1\phi_{1} and ϕ2\phi_{2} 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 n1sinϕ1=n2sinϕ2n_{1}\sin{\phi_{1}}\!=\!n_{2}\sin{\phi_{2}}. 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.

Refer to caption
Figure 1: (a) Normal and (b) oblique incidence of a ST wave packet on a planar interface.
Refer to caption
Figure 2: (a) Representation of the spectral support domain of a ST wave packet in a medium of refractive index n1n_{1} in the (x,z)(x,z) coordinate system that is aligned with the propagation direction of the ST wave packet. The spectrum lies at the intersection of the light-cone with a plane that is parallel to the kxk_{x}-axis and makes an angle θ1\theta_{1} with respect to the kzk_{z}-axis in (kx,kz,ωc)(k_{x},k_{z},\tfrac{\omega}{c}) space. The zz-axis is aligned with the plane-wave component at ω=ωo\omega\!=\!\omega_{\mathrm{o}} and kx=0k_{x}\!=\!0. (b) The spectral support domain at oblique incidence represented in the coordinate system (x,z)(x^{\prime},z^{\prime}) that is aligned with the interface but not with the propagation direction of the wave packet. The spectral support domain is rotated by ϕ1\phi_{1} (the incident angle) with respect to the ωc\tfrac{\omega}{c}-axis in (kx,kz,ωc)(k_{x}^{\prime},k_{z}^{\prime},\tfrac{\omega}{c}) space. (c) The spectral projection onto the (kx,ωc)(k_{x}^{\prime},\tfrac{\omega}{c})-plane from (b) is invariant across the interface at oblique incidence (dotted parabolic curve centered at kx0k_{x}^{\prime}\!\neq\!0). At normal incidence, the spectral projection onto the (kx,ωc)(k_{x},\tfrac{\omega}{c})-plane from (a) is invariant across the interface (solid parabolic curve centered at kx=0k_{x}\!=\!0). (d) Spectral support domain of the transmitted ST wave packet after oblique incidence. (e) Representation of the spectral support domain of a ST wave packet in a medium of refractive index n2n_{2}. The new spectral tilt angle θ2\theta_{2} is related to θ1\theta_{1} via Eq. 2.

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 kxk_{x} is associated with a single temporal frequency ω\omega; here xx and zz are the transverse and axial coordinates, respectively, and kxk_{x} and kzk_{z} are the corresponding components of the wave vector. We identify the propagation direction with the zz-axis, we assume a spatial spectrum that is symmetric around kx=0k_{x}\!=\!0, and we assign the temporal frequency ωo\omega_{\mathrm{o}} to the spatial frequency kx=0k_{x}\!=\!0. The spectral support domain for a ‘baseband’ ST wave packet Yessenov et al. (2019a) in a non-dispersive medium of index nn lies at the intersection of the light-cone kx2+kz2=n2(ωc)2k_{x}^{2}+k_{z}^{2}\!=\!n^{2}(\tfrac{\omega}{c})^{2} with a tilted plane ω=ωo+(kznko)ctanθ\omega\!=\!\omega_{\mathrm{o}}+(k_{z}-nk_{\mathrm{o}})c\tan{\theta} that is parallel to the kxk_{x}-axis and makes an angle θ\theta (the spectral tilt angle) with respect to the kzk_{z}-axis, where ko=ωo/ck_{\mathrm{o}}\!=\!\omega_{\mathrm{o}}/c and cc 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 v~=ctanθ=c/n~\widetilde{v}\!=\!c\tan{\theta}\!=\!c/\widetilde{n}, where n~=cotθ\widetilde{n}\!=\!\cot{\theta} is the group index. In other words, the group velocity is determined by an internal DoF (the spectral tilt angle θ\theta) that is related to the spatio-temporal field structure. Note that the ST wave packet reverts to a plane-wave pulse when n~=n\widetilde{n}\!=\!n. Within this framework, ωo\omega_{\mathrm{o}} is either the minimum temporal frequency when the wave packet is superluminal (v~>c/n\widetilde{v}\!>\!c/n), or the maximum temporal fewquency when the wave packet is subluminal (v~<c/n\widetilde{v}\!<\!c/n) Yessenov et al. (2019a).

When the ST wave packet is normally incident on the interface, the wave vector associated with ω=ωo\omega\!=\!\omega_{\mathrm{o}} is perpendicular to the interface [Fig. 1(a)], whereas all other temporal frequencies ω\omega are not. Two quantities are conserved across the interface: the transverse momentum kxk_{x} and the temporal frequency ω\omega. If the refractive indices of the media are n1n_{1} and n2n_{2}, and the group indices of the incident and transmitted wave packets are n~1\widetilde{n}_{1} and n~2\widetilde{n}_{2}, respectively, then it can be shown that:

n1(n1n~1)=n2(n2n~2).n_{1}(n_{1}-\widetilde{n}_{1})=n_{2}(n_{2}-\widetilde{n}_{2}). (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 n(nn~)n(n-\widetilde{n}) is related to the curvature of the spatio-temporal spectrum ωωoωo=kx2/ko22n(nn~)\tfrac{\omega-\omega_{\mathrm{o}}}{\omega_{\mathrm{o}}}\!=\!\tfrac{k_{x}^{2}/k_{\mathrm{o}}^{2}}{2n(n-\widetilde{n})} when projected onto the (kx,ωc)(k_{x},\tfrac{\omega}{c})-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 ωo\omega_{\mathrm{o}} is incident at an angle ϕ1\phi_{1} 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):

n1(n1n~1)cos2ϕ1=n2(n2n~2)cos2ϕ2,n_{1}(n_{1}-\widetilde{n}_{1})\cos^{2}{\phi_{1}}=n_{2}(n_{2}-\widetilde{n}_{2})\cos^{2}{\phi_{2}}, (2)

where ϕ2\phi_{2} is obtained from Snell’s law n1sinϕ1=n2sinϕ2n_{1}\sin{\phi_{1}}\!=\!n_{2}\sin{\phi_{2}}. Equation 2 indicates the invariance of a new spectral curvature n(nn~)cos2ϕn(n-\widetilde{n})\cos^{2}{\phi} at oblique incidence. The two curves corresponding to Eq. 1 and Eq. 2 intersect at the luminal point n~1=n1\widetilde{n}_{1}\!=\!n_{1} and n~2=n2\widetilde{n}_{2}\!=\!n_{2}, where the ST wave packets revert to plane-wave pulses [Fig. 3], and Eq. 2 reverts back to Eq. 1 at normal incidence ϕ1=ϕ2=0\phi_{1}\!=\!\phi_{2}\!=\!0.

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 (x,z)(x,z) – which is aligned with the propagation direction of the wave packet – lies at the intersection of the light-cone with the plane 𝒫B(θ)\mathcal{P}_{\mathrm{B}}(\theta) [Fig. 2(a)]. At normal incidence, the projection of this spectrum onto the (kx,ωc)(k_{x},\tfrac{\omega}{c})-plane is invariant across the interface Bhaduri et al. (2019a, 2020). At oblique incidence, the transverse wave number kxk_{x} is not parallel to the interface, and therefore is not conserved. However, in the coordinate system (x,z)(x^{\prime},z^{\prime}) that is aligned with the interface rather than with the propagation direction of the ST wave packet, the transverse wave number kx=kxcosϕkzsinϕk_{x}^{\prime}\!=\!k_{x}\cos{\phi}-k_{z}\sin{\phi} is conserved. In this coordinate system, the spectral support domain is rotated an angle ϕ1\phi_{1} around the ωc\tfrac{\omega}{c}-axis [Fig. 2(b)], and the projection of this new spectrum onto the (kx,ωc)(k_{x}^{\prime},\tfrac{\omega}{c})-plane is invariant across the interface [Fig. 2(c)]. It can be shown that kxcosϕk_{x}\cos{\phi} is invariant to first order (in the small-angle approximation Δkxko\Delta k_{x}\!\ll\!k_{\mathrm{o}}), and the new invariant spectral curvature at oblique incidence is n(nn~)cos2ϕn(n-\widetilde{n})\cos^{2}{\phi}. Because the opening angle of the light-cone changes from arctan(n1)\arctan{(n_{1})} in the first medium to arctan(n2)\arctan{(n_{2})} in the second, the invariance of the spectral projection onto the (kx,ωc)(k_{x}^{\prime},\tfrac{\omega}{c})-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 ϕ2\phi_{2} (rather than ϕ1\phi_{1}) around the ωc\tfrac{\omega}{c}-axis returns the spectral support domain to the (x,z)(x,z) coordinate system aligned with the propagation direction of the transmitted wave packet [Fig. 2(e)].

III Confirming the law of refraction at oblique incidence

Refer to caption
Figure 3: Verification of the law of refraction for ST wave packets at oblique incidence (Eq. 2) for ϕ1=20\phi_{1}\!=\!20^{\circ}, compared to normal incidence (Eq. 1) ϕ1=0\phi_{1}\!=\!0^{\circ}. The inset highlights the transition from normal refraction to anomalous refraction. We also identify the positions along the curves corresponding to group-velocity invariance v~=v~1\widetilde{v}\!=\!\widetilde{v}_{1} and group-velocity inversion v~2=v~1\widetilde{v}_{2}\!=\!-\widetilde{v}_{1}.

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 ϕ1\phi_{1} and sweep the spectral tilt angle θ1\theta_{1} of the ST wave packet incident from free space onto a 5-mm-thick layer of sapphire. For each value of θ1\theta_{1}, we measure the group delay incurred by the wave packet across the layer, from which we estimate θ2\theta_{2}. The measurements are carried out at normal incidence ϕ1=0\phi_{1}\!=\!0 and at oblique incidence ϕ1=20\phi_{1}\!=\!20^{\circ}. 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 (θ1,θ2)=(θ1L,θ2L)(\theta_{1},\theta_{2})\!=\!(\theta_{1}^{\mathrm{L}},\theta_{2}^{\mathrm{L}}). The curve for oblique incidence is shifted below that for normal incidence in the superluminal regime (θ1>θ1L\theta_{1}\!>\!\theta_{1}^{\mathrm{L}} and v~>c/n\widetilde{v}\!>\!c/n) and above it in the subluminal regime (θ1<θ1L\theta_{1}\!<\!\theta_{1}^{\mathrm{L}} and v~<c/n\widetilde{v}\!<\!c/n). We next explore the far-reaching ramifications of this change.

IV Dependence of the group velocity on the angle of incidence

Refer to caption
Figure 4: (a) Calculated group index n~2\widetilde{n}_{2} of the transmitted ST wave packet as a function of the angle of incidence ϕ1\phi_{1} and the group index n~1\widetilde{n}_{1} of the incident ST wave packet from free space (n1=1n_{1}\!=\!1) onto sapphire (n2=1.76n_{2}\!=\!1.76). The dashed line is the locus of luminal plane-wave pulses whose spatial and temporal DoFs are separable; here n~2=n2=1.76\widetilde{n}_{2}\!=\!n_{2}\!=\!1.76 is independent of ϕ1\phi_{1}. (b) Measured change in n~2\widetilde{n}_{2} at oblique incidence onto MgF2 (n2=1.38n_{2}\!=\!1.38) from free space with respect to normal incidence, Δn~2=n~2(ϕ1)n~2(0)\Delta\widetilde{n}_{2}\!=\!\widetilde{n}_{2}(\phi_{1})-\widetilde{n}_{2}(0). Measurements are carried out for a ST wave packet in the subluminal regime θ1=30\theta_{1}\!=\!30^{\circ} (v~=0.58c\widetilde{v}\!=\!0.58c) and in the superluminal regime θ1=109\theta_{1}\!=\!109^{\circ} (v~=2.9c\widetilde{v}\!=\!-2.9c). (c) Same as (a) for incidence from free space onto sapphire.

Figure 3 yields the relationship between the group indices n~1=cotθ1\widetilde{n}_{1}\!=\!\cot{\theta_{1}} and n~2=cotθ2\widetilde{n}_{2}\!=\!\cot{\theta_{2}} for the incident and refracted ST wave packets, respectively, at a fixed incident angle ϕ1\phi_{1}. We next examine the change in the group index n~2\widetilde{n}_{2} for the transmitted wave packet as we vary the incidence angle ϕ1\phi_{1} while holding fixed the group index n~1\widetilde{n}_{1} for the incident wave packet. In Fig. 4(a), we plot the calculated n~2\widetilde{n}_{2} while varying both n~1\widetilde{n}_{1} and ϕ1\phi_{1} for incidence from free space n1=1n_{1}\!=\!1 onto sapphire n2=1.76n_{2}\!=\!1.76. This calculation highlights the above-mentioned ‘memory’ effect. The group index n~2\widetilde{n}_{2} of the transmitted wave packet depends not only the refractive index of the second medium n2n_{2}, but also on the characteristics of the incident wave packet: its group index n~1\widetilde{n}_{1} and incident angle ϕ1\phi_{1}. The luminal condition for the incident ST wave packet is n~1=n1=1\widetilde{n}_{1}\!=\!n_{1}\!=\!1, in which case this wave packet is simply a plane-wave pulse, and the refracted wave packet is also a luminal plane-wave pulse n~2=n2=1.76\widetilde{n}_{2}\!=\!n_{2}\!=\!1.76. For subluminal incidence ST wave packets n~1>n1\widetilde{n}_{1}\!>\!n_{1}, n~2\widetilde{n}_{2} decreases with ϕ2\phi_{2}; i.e., the group velocity of the transmitted subluminal wave packet increases with ϕ1\phi_{1} in the higher-index medium. Conversely, for superluminal incident ST wave packets n~1<n1\widetilde{n}_{1}\!<\!n_{1}, n~2\widetilde{n}_{2} increases with ϕ1\phi_{1}; i.e., the group velocity of the transmitted superluminal wave packet increases with ϕ1\phi_{1} in the higher-index medium. In the case of a conventional wave packet whose spatial and temporal DoFs are separable, n~2\widetilde{n}_{2} is independent of ϕ1\phi_{1}.

This behavior can be understood on the basis of Eq. 2. Substituting for ϕ2\phi_{2} in terms of ϕ1\phi_{1} from Snell’s law yields:

n~2(ϕ1)n2=1+(n~1n11)η(ϕ1),\frac{\widetilde{n}_{2}(\phi_{1})}{n_{2}}=1+\left(\frac{\widetilde{n}_{1}}{n_{1}}-1\right)\eta(\phi_{1}), (3)

where the factor η(ϕ1)\eta(\phi_{1}) depends on the indices n1n_{1} and n2n_{2} through

η(ϕ1)=cos2ϕ1(n2n1)2sin2ϕ1.\eta(\phi_{1})=\frac{\cos^{2}{\phi_{1}}}{(\frac{n_{2}}{n_{1}})^{2}-\sin^{2}{\phi_{1}}}. (4)

The behavior of η(ϕ1)\eta(\phi_{1}) depends on the ratio n1/n2n_{1}/n_{2}. When n1<n2n_{1}\!<\!n_{2}, η(ϕ1)\eta(\phi_{1}) drops monotonically from its initial value to 0 at ϕ1=90\phi_{1}\!=\!90^{\circ}. Because the term (n~1n11)(\tfrac{\widetilde{n}_{1}}{n_{1}}-1) is positive in the subluminal regime n~1>n1\widetilde{n}_{1}\!>\!n_{1}, the decrease of η(ϕ1)\eta(\phi_{1}) with ϕ1\phi_{1} results in a decrease of n~2(ϕ1)\widetilde{n}_{2}(\phi_{1}) with ϕ1\phi_{1}. In the superluminal regimes, the term (n~1n11)(\tfrac{\widetilde{n}_{1}}{n_{1}}-1) is negative, and n~2(ϕ1)\widetilde{n}_{2}(\phi_{1}) increase with ϕ1\phi_{1}. On the other hand, when n1>n2n_{1}\!>\!n_{2}, η(ϕ1)\eta(\phi_{1}) increases monotonically with ϕ1\phi_{1} and reaches a singularity when sinϕ1=n2/n1\sin{\phi_{1}}\!=\!n_{2}/n_{1}, corresponding to the critical angle (total internal reflection) at the interface. Consequently, the opposite trends for n~2\widetilde{n}_{2} with ϕ1\phi_{1} ensue.

Our measurements confirm this predicted behavior. In Fig. 4(b,c) we plot the change in the refracted group index n~2\widetilde{n}_{2} with ϕ1\phi_{1} with respect to that at normal incidence, Δn~2=n~2(ϕ1)n~2(0)\Delta\widetilde{n}_{2}\!=\!\widetilde{n}_{2}(\phi_{1})-\widetilde{n}_{2}(0), 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 (θ1<45\theta_{1}\!<\!45^{\circ}) and superluminal (θ1>45\theta_{1}\!>\!45^{\circ}) regimes. We note several general observations about the results that are expected from our analysis. First, n~2(ϕ1)\widetilde{n}_{2}(\phi_{1}) deviates monotonically from the normal-incidence value n~2(0)\widetilde{n}_{2}(0). Second, the group index drops at oblique incidence n~2(ϕ1)<n~2(0)\widetilde{n}_{2}(\phi_{1})\!<\!\widetilde{n}_{2}(0) in the subluminal regime, and increases n~2(ϕ1)>n~2(0)\widetilde{n}_{2}(\phi_{1})\!>\!\widetilde{n}_{2}(0) 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 ϕ1\phi_{1}. 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 v~2=v~1\widetilde{v}_{2}\!=\!\widetilde{v}_{1}, which occurs when n~1=n~th=n1+n2\widetilde{n}_{1}\!=\!\widetilde{n}_{\mathrm{th}}\!=\!n_{1}+n_{2}; (2) anomalous refraction whereby v~2>v~1\widetilde{v}_{2}\!>\!\widetilde{v}_{1} when n2>n1n_{2}\!>\!n_{1}, which occurs when n~1>n~th\widetilde{n}_{1}\!>\!\widetilde{n}_{\mathrm{th}}; and (3) group-velocity inversion v~2=v~1\widetilde{v}_{2}\!=\!-\widetilde{v}_{1}, which occurs when n~1=n1n2\widetilde{n}_{1}\!=\!n_{1}-n_{2}; 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):

n~1=n~th(ϕ1)=n1+n21+n1n2sin2ϕ1,\widetilde{n}_{1}=\widetilde{n}_{\mathrm{th}}(\phi_{1})=\frac{n_{1}+n_{2}}{1+\frac{n_{1}}{n_{2}}\sin^{2}{\phi_{1}}}, (5)

where n~th(ϕ1)<n~th(0)=n1+n2\widetilde{n}_{\mathrm{th}}(\phi_{1})\!<\!\widetilde{n}_{\mathrm{th}}(0)\!=\!n_{1}+n_{2}. We plot in Fig. 5 measurements of the group delay incurred upon traversing equal lengths (L=5L\!=\!5 mm) of free space (n1=1n_{1}\!=\!1) and MgF2 (n2=1.38n_{2}\!=\!1.38). At normal incidence n~th(0)=2.38\widetilde{n}_{\mathrm{th}}(0)\!=\!2.38 corresponding to θ122.8\theta_{1}\!\approx\!22.8^{\circ} [Fig. 5(a)], and at oblique incidence ϕ1=30\phi_{1}\!=\!30^{\circ} we have n~th(30)=1.86\widetilde{n}_{\mathrm{th}}(30^{\circ})\!=\!1.86, corresponding to θ28.2\theta\!\approx\!28.2^{\circ} [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 θ1<22.8\theta_{1}\!<\!22.8^{\circ}, this regime is expanded for oblique incidence to θ1<28.2\theta_{1}\!<\!28.2^{\circ} at ϕ1=30\phi_{1}\!=\!30^{\circ}.

Refer to caption
Figure 5: Impact of oblique incidence on the condition for realizing group-velocity invariance. (a) Normal incidence, ϕ1=0\phi_{1}\!=\!0. We plot the on-axis x=0x\!=\!0 input-pulse profile I(z=0;τ)I(z\!=\!0;\tau) on the left, and the output pulse profiles I(L;τ)I(L;\tau) after traversing L=5L\!=\!5 mm of free space (dotted curve) and MgF2 (solid curve) at normal incidence for a ST wave packet with θ1=22.8\theta_{1}\!=\!22.8^{\circ} (v~=0.42c\widetilde{v}\!=\!0.42c) synthesized in free space. (b) Same as (a) for oblique incidence at ϕ1=30\phi_{1}\!=\!30^{\circ}. Here the group delays traversing equal lengths of free space and MgF2 are equal but have a smaller overall delay than in (a), despite the longer distance traveled.
Refer to caption
Figure 6: Impact of oblique incidence on group-velocity inversion. (a) At normal incidence (ϕ1=0\phi_{1}\!=\!0), a ST wave packet with θ0=145\theta_{0}\!=\!145^{\circ} traverses a 5-mm-thick layer of MgF2 with a negative group delay (dashed curve; v~1=2.6c\widetilde{v}_{1}\!=\!-2.6c and θ1=111\theta_{1}\!=\!111^{\circ} in MgF2) and a 5-mm-thick layer of sapphire with a positive group delay (dotted curve; v~2=2.6c\widetilde{v}_{2}\!=\!2.6c and θ1=69\theta_{1}\!=\!69^{\circ} in sapphire), corresponding to group delay inversion. When the ST wave packet is incident from free space on a bilayer of MgF2 and sapphire, the total group delay is zero (solid curve), corresponding to group-delay cancellation. (b) Same as (a) for oblique incidence at ϕ1=30\phi_{1}\!=\!30^{\circ}. The incident ST wave packet from free space has θ0=156\theta_{0}\!=\!156^{\circ}. The delays incurred by the two pulses in the individual layers in larger than at normal incidence because the distance traveled is larger at oblique incidence and the wave packet is in the normal refraction regime.

Second, group-velocity inversion occurs at oblique incidence when Yessenov et al. (2021):

n~1(ϕ1)=n1n21n1n2sin2ϕ1,\widetilde{n}_{1}(\phi_{1})=\frac{n_{1}-n_{2}}{1-\frac{n_{1}}{n_{2}}\sin^{2}{\phi_{1}}}, (6)

where n~1(ϕ1)<n~1(0)=n1n2\widetilde{n}_{1}(\phi_{1})\!<\!\widetilde{n}_{1}(0)\!=\!n_{1}-n_{2} when n2>n1n_{2}\!>\!n_{1}. We plot in Fig. 6 the measured group delay incurred upon traversing equal lengths (L=5L\!=\!5 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 n~01.43\widetilde{n}_{0}\!\approx\!-1.43 (θ0145\theta_{0}\!\approx\!145^{\circ}); see Fig. 6(a). At oblique incidence (ϕ1=30\phi_{1}\!=\!30^{\circ}), the required free-space group index to realize group-velocity inversion is n~02.25\widetilde{n}_{0}\approx\!-2.25 (θ0156\theta_{0}\!\approx\!156^{\circ}); 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

Refer to caption
Figure 7: Concept of blind synchronization. We depict the configuration for the blind synchronization of multiple receivers (Rx) receiving pulses from a transmitter (Tx). All the receivers are at the same depth d2d_{2} beneath the interface, but their locations are otherwise not known a priori. When using conventional wave packets, the group delay increases as the physical path length from Tx to Rx increases, which occurs when the incident angle ϕ1\phi_{1} at the interface increases. By employing ST wave packets instead, the group delay accrued along the different pathways from Tx to Rx can be held constant because the group velocity of the transmitted wave packet increases with ϕ1\phi_{1}, thereby compensating for the longer path length covered.

Consider the configuration shown in Fig. 7. A transmitter (Tx) in the first medium at a height d1d_{1} above the interface directs an ST wave packet at different angles to receivers (Rx) at different locations but at the same depth d2d_{2} beneath the interface in the second medium. The optical delay in the first medium from Tx to the interface is τ1(ϕ1)=τ1(0)/cosϕ1\tau_{1}(\phi_{1})\!=\!\tau_{1}(0)/\cos{\phi_{1}}, where τ1(0)=d1n~1/c\tau_{1}(0)\!=\!d_{1}\widetilde{n}_{1}/c is the group delay at normal incidence. The group velocity of the ST wave packet is the same in all directions independently of ϕ1\phi_{1} because the medium is homogeneous and isotropic. Therefore, the group delay τ1(ϕ1)\tau_{1}(\phi_{1}) increases with ϕ1\phi_{1}. The group delay in the second medium is τ2(ϕ1)=d2n~2(ϕ1)/(ccosϕ2)\tau_{2}(\phi_{1})\!=\!d_{2}\widetilde{n}_{2}(\phi_{1})/(c\cos{\phi_{2}}), where n~2(ϕ1)\widetilde{n}_{2}(\phi_{1}) 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 ϕ1\phi_{1}.

The relative group delay between obliquely and normally incident ST wave packets in the first medium is:

Δτ1=τ1(ϕ1)τ1(0)=d1cn~1(1cosϕ11),\Delta\tau_{1}=\tau_{1}(\phi_{1})-\tau_{1}(0)=\frac{d_{1}}{c}\widetilde{n}_{1}\left(\frac{1}{\cos{\phi_{1}}}-1\right), (7)

and in the second medium is:

Δτ2=τ2(ϕ1)τ2(0)=d2c(n~2(ϕ1)cosϕ2n~2(0)).\Delta\tau_{2}=\tau_{2}(\phi_{1})-\tau_{2}(0)=\frac{d_{2}}{c}\left(\frac{\widetilde{n}_{2}(\phi_{1})}{\cos{\phi_{2}}}-\widetilde{n}_{2}(0)\right). (8)

We define the total relative group delay Δτ=Δτ1+Δτ2=τ(ϕ1)τ(0)\Delta\tau\!=\!\Delta\tau_{1}+\Delta\tau_{2}\!=\!\tau(\phi_{1})-\tau(0), where τ(ϕ1)=τ1(ϕ1)+τ2(ϕ1)\tau(\phi_{1})\!=\!\tau_{1}(\phi_{1})+\tau_{2}(\phi_{1}) is the total group delay incurred from Tx to Rx across the interface. Blind synchronization requires Δτ=0\Delta\tau\!=\!0 independently of ϕ1\phi_{1}, such that the ST wave packet arrives simultaneously at all receivers at the depth d2d_{2} despite the different physical distances covered along the different trajectories.

We first ascertain whether realizing Δτ=0\Delta\tau\!=\!0 is physically feasible. Note that Δτ1\Delta\tau_{1} 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 Δτ2\Delta\tau_{2} 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 v~2(ϕ1)>v~2(0)\widetilde{v}_{2}(\phi_{1})\!>\!\widetilde{v}_{2}(0). From the analysis and measurements reported above, this scenario occurs for subluminal ST wave packets when n2>n1n_{2}\!>\!n_{1}. We next consider whether the blind synchronization condition can be met quantitatively.

Refer to caption
Figure 8: Plots of the relative group delay Δτ\Delta\tau (normalized with respect to d1/cd_{1}/c) as a function of the incident angle ϕ1\phi_{1} and the group index n~1\widetilde{n}_{1} of the incident ST wave packet: (a) n1=1n_{1}\!=\!1, n2=1.76n_{2}\!=\!1.76, and d2/d1=5d_{2}/d_{1}\!=\!5; and (b) n1=1n_{1}\!=\!1, n2=2n_{2}\!=\!2, and d2/d1=10d_{2}/d_{1}\!=\!10.

VI.2 Calculated group delay of ST wave packets between transmitter and receivers

The total relative group delay Δτ=Δτ1+Δτ2\Delta\tau\!=\!\Delta\tau_{1}+\Delta\tau_{2} is given by:

Δτ(ϕ1)d1/c\displaystyle\frac{\Delta\tau(\phi_{1})}{d_{1}/c}\!\!\!\! =\displaystyle= n2d2d1{(1cosϕ11)+(n~1n11)(η(ϕ1)cosϕ2η(0))}\displaystyle\!\!\!\!n_{2}\frac{d_{2}}{d_{1}}\left\{\!\left(\!\frac{1}{\cos{\phi_{1}}}-1\right)+\left(\!\frac{\widetilde{n}_{1}}{n_{1}}\!-\!1\!\right)\left(\!\frac{\eta(\phi_{1})}{\cos{\phi_{2}}}-\eta(0)\!\right)\!\right\} (9)
+\displaystyle+ n~1(1cosϕ11).\displaystyle\widetilde{n}_{1}\left(\frac{1}{\cos{\phi_{1}}}\!-\!1\!\right).

where we made use of the law of refraction in Eq. 2, and we normalized Δτ\Delta\tau with respect to d1/cd_{1}/c. We plot in Fig. 8 calculations for Δτ\Delta\tau as a function of the incident angle ϕ1\phi_{1} and the group index of the incident wave packet n~1\widetilde{n}_{1} in two different configurations. In the first configuration we have n1=1n_{1}\!=\!1, n2=1.76n_{2}\!=\!1.76, and d2/d1=5d_{2}/d_{1}\!=\!5; and in the second configuration we have n1=1n_{1}\!=\!1, n2=2n_{2}\!=\!2, and d2/d1=10d_{2}/d_{1}\!=\!10. In both cases we find the contour Δτ=0\Delta\tau\!=\!0 occurs at a specific group index n~1\widetilde{n}_{1} 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 d2/d1d_{2}/d_{1} and/or n2/n1n_{2}/n_{1} increase.

For the sake of comparison, we plot in Fig. 9(a) the relative group delay Δτ\Delta\tau 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 c/n1c/n_{1} and in the second at c/n2c/n_{2} (assuming non-dispersive media). The relative group delay is given by:

Δτconv(ϕ1)d1/c=n1(1cosϕ11)+n2d2d1(1cosϕ21).\frac{\Delta\tau_{\mathrm{conv}}(\phi_{1})}{d_{1}/c}=n_{1}\left(\frac{1}{\cos{\phi_{1}}}-1\right)+n_{2}\frac{d_{2}}{d_{1}}\left(\frac{1}{\cos{\phi_{2}}}-1\right). (10)

The relative group delay Δτ\Delta\tau is always positive for conventional pulses because the pulses traveling longer distances inevitably accrues larger group delays. In contrast, varying n~1\widetilde{n}_{1} for ST wave packets changes Δτ\Delta\tau 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 Δτ=0\Delta\tau\!=\!0. In Fig. 9(b) we highlight this transition from Δτ>0\Delta\tau\!>\!0 to Δτ<0\Delta\tau\!<\!0 as n~1\widetilde{n}_{1} is changed in the vicinity of n~1=3\widetilde{n}_{1}\!=\!3.

Refer to caption
Figure 9: (a) Comparison of the relative group delay Δτ\Delta\tau (normalized with respect to d1/cd_{1}/c) of ST wave packets (solid curves; Eq. 9) from Fig. 8(a) at n~1=1.46\widetilde{n}_{1}\!=\!1.46, 3 and 4 with Δτ\Delta\tau of a traditional pulse (dashed curve; Eq. 10). (b) Same as (a), but with a reduced range for Δτ\Delta\tau along the vertical axis.
Refer to caption
Figure 10: Proof-of-principle realization of blind synchronization. (a) Schematic of the experimental configuration. (b) Measured group delay accrued over 1 mm of free space by a ST wave packet with n~1=3\widetilde{n}_{1}\!=\!3. The delay increases with ϕ1\phi_{1}. (c) Measured group delay over 5 mm of sapphire for the ST wave packet from (b). Here, the group delay decreases with ϕ1\phi_{1}. (d) Sum of the group delays from (b) and (c), showing an approximately constant value despite the different path lengths associated with different incident angles ϕ1\phi_{1}.

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 (n1=1n_{1}\!=\!1) to sapphire (n2=1.76n_{2}\!=\!1.76) with the propagation distances at normal incidence having the ratio d2/d1=5d_{2}/d_{1}\!=\!5; we use here d1=1d_{1}\!=\!1 mm and d2=5d_{2}\!=\!5 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 d1d1/cosϕ1d_{1}\!\rightarrow\!d_{1}/\cos{\phi_{1}} over which the delay is measured.

We synthesize a subluminal ST wave packet in free space with n~1=3\widetilde{n}_{1}\!=\!3 (θ1=18\theta_{1}\!=\!18^{\circ}). According to Fig. 8(a), blind synchronization Δτ=0\Delta\tau\!=\!0 is realized by such a wave packet over a range of incidence angles ϕ1=0±35\phi_{1}\!=\!0\rightarrow\pm 35^{\circ}. We plot in Fig. 10(b) the measured group delay in free space over distances extending from d1=1d_{1}\!=\!1 mm (corresponding to ϕ1=0\phi_{1}\!=\!0^{\circ}) to d1=1.2d_{1}\!=\!1.2 mm (corresponding to ϕ1=35\phi_{1}\!=\!35^{\circ}), with d1(ϕ1)=d1/cosϕ1d_{1}(\phi_{1})\!=\!d_{1}/\cos{\phi_{1}}. 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 ϕ1\phi_{1}. As shown in Fig. 4, this is a general feature at oblique incidence for subluminal ST wave packets when n2>n1n_{2}\!>\!n_{1}. For this particular ST wave packet with n~1=3\widetilde{n}_{1}\!=\!3, the decrease in delay in the layer (increase in group velocity) with ϕ1\phi_{1} 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 ϕ1=0\phi_{1}\!=\!0^{\circ} to ϕ1=35\phi_{1}\!=\!35^{\circ}.

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 (L/cosϕ2L/\cos{\phi_{2}}, where LL is the layer thickness). This can be seen as a special case of blind synchronization where we set d1=0d_{1}\!=\!0. 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 n~1\widetilde{n}_{1} (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).