Startseite Spin-controlled photonics via temporal anisotropy
Artikel Open Access

Spin-controlled photonics via temporal anisotropy

  • Carlo Rizza ORCID logo EMAIL logo , Giuseppe Castaldi und Vincenzo Galdi ORCID logo EMAIL logo
Veröffentlicht/Copyright: 21. Februar 2023
Artikel downloaden (PDF)
Veröffentlichen auch Sie bei De Gruyter Brill
Manuskript einreichen Informationen für Autor*innen
Nanophotonics
Aus der Zeitschrift Nanophotonics Band 12 Heft 14

Abstract

Temporal metamaterials, based on time-varying constitutive properties, offer new exciting possibilities for advanced field manipulations. In this study, we explore the capabilities of anisotropic temporal slabs, which rely on abrupt changes in time from isotropic to anisotropic response (and vice versa). Our findings show that these platforms can effectively manipulate the wave-spin dimension, allowing for a range of intriguing spin-controlled photonic operations. We demonstrate these capabilities through examples of spin-dependent analog computing and spin–orbit interaction effects for vortex generation. These results provide new insights into the field of temporal metamaterials, and suggest potential applications in communications, optical processing and quantum technologies.

Keywords: analog computing; anisotropy; metamaterials; spin-orbit interaction; time-varying

1 Introduction

In recent years, there has been mounting interest in “temporal” and “space-time” metamaterials [1–3]. These are artificial materials where the spatial modulation of the constitutive parameters is replaced by (or combined with) time-varying properties. This field of research, which has roots in longstanding theoretical foundations [4–7], has been spurred by the emergence of new physical concepts such as “time crystals” [8, 9] and major technological advances in rapidly reconfigurable material constituents across the electromagnetic (EM) spectrum [10–14].

Accessing the temporal dimension offers the potential for advanced spatial-spectral field manipulations and for surpassing fundamental limitations of linear, time-invariant systems [15]. Accordingly, a broad variety of concepts relying on space-time analogies have been put forward, ranging from relatively simple ideas such as temporal boundaries [16], interfaces [17], and slabs [18], to more sophisticated ones including gratings [19, 20], filters [21, 22], photonic time crystals [23–25], antireflection coatings [26–29], and absorbers [30]. Recent reviews and perspectives on this rapidly advancing research area can be found in Refs. [31, 32]. While there are inherent constraints and technological challenges in the temporal modulation of constitutive parameters [33], experimental studies are continuing to progress, and recent results have demonstrated feasibility in this area [34–36].

Of special interest for this study are the recent results on anisotropic temporal metamaterials, featuring abrupt temporal switching from isotropic to anisotropic responses. These include, for instance, temporal aiming [37] and Brewster angle [38], spatiotemporal isotropic-to-anisotropic meta-atoms [39], complete polarization conversion [40], nonreciprocity and Faraday rotation [41], and spin-temporal interactions [42].

Here, we investigate the capabilities of anisotropic temporal platforms to attain spin-controlled field manipulations. Specifically, within the recently proposed framework of short-pulsed temporal metamaterials [43, 44], we show that anisotropy can be leveraged to attain spin-dependent analog computing on an impinging wavepacket, and to overcome some limitations inherent of isotropic scenarios. Additionally, we demonstrate that typical spin–orbit interaction effects observed in spatial, anisotropic scenarios [45] can be translated to the temporal case, enabling efficient vortex generation. These findings highlight the potential of temporal anisotropic metamaterials for advanced field manipulations, with diverse and wide-ranging possible applications to communications, optical processing, and quantum technologies.

2 Results and discussion

2.1 Problem schematic and statement

As schematically illustrated in Figure 1, we consider a spatially unbounded, nonmagnetic medium described by the constitutive relationships

(1) D = ε 0 ε ̲ ̲ ( t ) ⋅ E , B = μ 0 H ,

which relate the electric and magnetic inductions (D and B, respectively) to the corresponding fields (E and H, respectively). In Equations (1), ɛ 0 and μ 0 are the vacuum dielectric permittivity and magnetic permeability, respectively, and

(2) ε ̲ ̲ ( t ) = ε 1 I ̲ ̲ , t < 0 , ε ̲ ̲ a = ε ⊥ e ̂ x ⊗ e ̂ x + e ̂ y ⊗ e ̂ y + ε ‖ e ̂ z ⊗ e ̂ z , 0 < t < τ , ε 2 I ̲ ̲ , t > τ ,

is a time-varying relative-permittivity tensor, with ɛ ν (ν = 1, 2, ⊥, ‖) denoting real-valued constants, I ̲ ̲ the identity tensor, and ε ̲ ̲ a a uniaxial tensor. Here and henceforth, e ̂ α indicates an α-directed unit vector, and ⊗ the dyadic product.

Figure 1: Schematic representation of an anisotropic temporal slab (details in the text).
Figure 1:

Schematic representation of an anisotropic temporal slab (details in the text).

Equation (2) describes a time-varying medium which, at time t = 0, undergoes an abrupt transition from a stationary, isotropic state to a uniaxially anisotropic response; subsequently, at the time instant t = τ, the response is abruptly switched back to an isotropic state, which is then maintained indefinitely. In analogy with previous studies [40], we refer to the configuration above as an “anisotropic temporal slab.” We assume that the EM fields experience ideal temporal boundaries, i.e., discontinuous changes of the dielectric permittivity at t = 0 and t = τ. However, in our full-wave simulations (see the Methods Section 4.3) we take into account finite rising/falling times (much shorter than the temporal slab duration τ). Moreover, as in previous studies on this subject [37–40], we assume to be far away from any material resonance, so as to neglect temporal dispersion. Considering the impact of dispersion is technically possible, but would necessitate a more sophisticated approach than the transfer-matrix method used here. Recent research has indicated that sudden variations in the plasma frequency of a Lorentzian-type dispersive medium can lead to the emergence of two shifted frequencies and the need for additional boundary conditions [46]. This extension will be addressed in future studies.

From Maxwell’s equations, the EM field dynamics can be described by the vector wave equation for the electric induction, viz.,

(3) ∂ 2 D ∂ t 2 + c 2 ∇ × ∇ × ε ̲ ̲ − 1 ( t ) ⋅ D = 0 ,

where −1 denotes the inverse operator, and c is the wavespeed in vacuum. In what follows, we will derive a general analytical solution, which will be subsequently particularized to two scenarios of interest.

2.2 General theory

We start by considering plane-wave solutions of the type D ( r , t ) = Re d ( k , t ) e i k ⋅ r , with i denoting the imaginary unit, r = x e ̂ x + y e ̂ y + z e ̂ z and k = k x e ̂ x + k y e ̂ y + k z e ̂ z the position and wave vectors, respectively, and d(k, t) a time-dependent plane-wave spectrum. From Equations (1) and (2), together with the Maxwell’s curl equation ∇ × E = −∂ t B, we identify the eigenwaves of the temporal anisotropic slab, i.e., the ordinary (s-polarized) and extraordinary (p-polarized) plane waves, with eigenfrequencies

(4) ω p = c k ⊥ 2 n ‖ 2 + k z 2 n ⊥ 2 , ω s = c n ⊥ k ,

and polarization unit vectors

(5a) e ̂ p = k x k z k k ⊥ e ̂ x + k y k z k k ⊥ e ̂ y − k ⊥ k e ̂ z ,

(5b) e ̂ s = − k y k ⊥ e ̂ x + k x k ⊥ e ̂ y ,

where n ⊥ = ε ⊥ , n ⊥ = ε ‖ , k = k , and k ⊥ = k ⊥ = k x 2 + k y 2 . Equations (5) clearly hold for k ⊥ ≠ 0; for the case k ⊥ = 0, we define instead

(6) e ̂ p = e ̂ x , e ̂ s = e ̂ y .

Next, we investigate the scattering of a time-harmonic plane wave, which, for t < 0, can be written as D ( i ) ( r , t ) = Re d ( i ) ( k ) e i k ⋅ r − ω 1 t . The interaction with the anisotropic temporal slab described by Equations (2) will generate forward (transmitted) and backward (reflected) waves, which, for t > τ, can be written as D ( t ) ( r , t ) = Re d ( t ) ( k ) e i k ⋅ r − ω 2 ( t − τ ) and D ( r ) ( r , t ) = Re d ( r ) ( k ) e i k ⋅ r + ω 2 ( t − τ ) , respectively, with ω 1,2 = ck/n 1,2 and n 1,2 = ε 1,2 denoting the angular frequencies and refractive indices, respectively, in the initial/final medium. By enforcing the conventional temporal boundary conditions (i.e., continuity of the electric and magnetic inductions at t = 0, τ) [16], we can obtain the temporal transmission and reflection matrices connecting the transmitted [d (t)] and reflected [d (r)] vector amplitudes of the electric induction in the polarization basis of Equations (5), viz., d ( t ) = T ̲ ̲ ⋅ d ( i ) , d ( r ) = R ̲ ̲ ⋅ d ( i ) , where

(7a) T ̲ ̲ ( k ) = T p p ( k ) e ̂ p ⊗ e ̂ p + T s s ( k ) e ̂ s ⊗ e ̂ s ,

(7b) R ̲ ̲ ( k ) = R p p ( k ) e ̂ p ⊗ e ̂ p + R s s ( k ) e ̂ s ⊗ e ̂ s ,

and d ( j ) = d p ( j ) e ̂ p + d s ( j ) e ̂ s , with j = i, r, t. The expressions of the scattering coefficients (T pp , R pp , T ss , R ss ) are reported in the Methods Section 4.1, together with the possible generalization to an arbitrary polarization basis.

2.3 Representative examples

2.3.1 Short-pulsed regime: spin-dependent analog computing

As a leading first example, we consider the short-pulsed regime τ ≪ Δt, with Δt denoting a characteristic timescale of the wave dynamics. In our previous studies on isotropic configurations [43], we have shown that such regime may be interpreted as a nonlocal temporal boundary, whose response can be harnessed so as to perform elementary analog computing (e.g., derivatives) on an impinging wavepacket. Here, we explore to what extent anisotropy can be leveraged to attain spin-dependent operations.

To this aim, we assume circularly polarized plane waves propagating along the x-axis (i.e., k = k x e ̂ x ), and label with the subscripts “+” and “−” the associated spin, corresponding to the unit vectors e ̂ ± = ( − e ̂ z ± i e ̂ y ) / 2 , i.e., left- or right-handed circular (LHC or RHC) polarization, respectively, for the assumed incidence direction. In this case, it can be shown (see the Methods Section 4.1 for details) that the relevant transmission and reflection coefficients can be approximated as

(8a) T + + ( k x ) = T − − ( k x ) ≃ 1 2 1 + n 2 n 1 − i π 2 2 n 1 + n 2 n ⊥ 2 + n 2 n ‖ 2 k x K + O k x 2 K 2 ,

(8b) T + − ( k x ) = T − + ( k x ) ≃ i π 2 n 2 1 n ⊥ 2 − 1 n ‖ 2 k x K + O k x 2 K 2 ,

and

(8c) R + + ( k x ) = R − − ( k x ) ≃ 1 2 1 − n 2 n 1 − i π 2 2 n 1 − n 2 n ⊥ 2 − n 2 n ‖ 2 k x K + O k x 2 K 2 ,

(8d) R + − ( k x ) = R − + ( k x ) ≃ − i π 2 n 2 1 n ⊥ 2 − 1 n ‖ 2 k x K + O k x 2 K 2 ,

where K = 2π/(cτ) and O is the Landau symbol. Here and henceforth, T ±±, R ±± (T ∓±, R ∓±) denote the co-polar (cross-polar) transmission and reflection coefficients, respectively.

Equations (8) represent the first main results of our study, since they clearly illustrate the impact of local contributions (i.e., constant terms) and nonlocal terms (i.e., proportional to k x ) and their interactions with the wave spin. Remarkably, these results suggest how to exploit temporal anisotropy for designing unconventional spin-dependent analog computations. For instance, it is evident that the co-polar coefficients (T ±±, R ±±) are generally dominated by local terms, whereas the leading terms in the cross-polar ones (T ∓±, R ∓±) are nonlocal (∝ik x ); recalling the well-known property of the Fourier transform, these latter are amenable to first derivatives. Therefore, as schematically illustrated in Figure 2, for an impinging wavepacket with a given spin and profile (Figure 2(A)), we expect the co-polarized transmitted/reflected wavepackets to generally exhibit a similar profile, and the cross-polarized transmitted/reflected ones to be essentially proportional to its first derivative (Figure 2(B)).

Figure 2: Schematic illustration of spin-dependent analog computing. (A) A wavepacket with a positive spin (i.e., LHC) impinges in the initial medium. (B) After the short-pulsed anisotropic temporal modulation (t > τ), the co-polarized reflected (backward) and transmitted (forward) wavepacket exhibit the same profile as the impinging one, whereas the cross-polarized ones are proportional to its first derivative.
Figure 2:

Schematic illustration of spin-dependent analog computing. (A) A wavepacket with a positive spin (i.e., LHC) impinges in the initial medium. (B) After the short-pulsed anisotropic temporal modulation (t > τ), the co-polarized reflected (backward) and transmitted (forward) wavepacket exhibit the same profile as the impinging one, whereas the cross-polarized ones are proportional to its first derivative.

By comparison with isotropic configurations (ɛ ⊥ = ɛ ‖) [43], where the analog-computing capabilities emerge only in the reflection (backward) response for impedance-matching conditions (ɛ 1 = ɛ 2), we note that the anisotropy character enables (via polarization conversion) extended operations, also in the absence of impedance matching and in the transmission (forward) response. Most important, by controlling the wave polarization, it enables elementary spin-dependent analog computing.

For illustration and validation, we assume an incident wavepacket with a positive spin and Gaussian profile,

(9) D ( i ) ( x , t ) = D 0 e − x − v 1 t − t s v 1 σ t 2 e ̂ + , t < 0 ,

with D 0 denoting a constant amplitude, σ t a characteristic timescale, v 1 = c / ε 1 , and t s = −5σ t .

As a first illustrative example, we consider as initial and final relative permittivities ɛ 1 = 1 and ɛ 2 = 4, respectively, and an anisotropic temporal slab with ɛ ⊥ = 1 and ɛ ‖ = 4, of duration τ = 0.5σ t so as to fulfill the short-pulsed assumption. Under these conditions, the co-polar transmission and reflection responses in Equations (8a) and (8c), respectively, are dominated by local terms, whereas the leading terms in the cross-polar responses (see Equations (8b) and (8d)) are nonlocal (∝ ik x ). Figure 3(A) and (B) show the corresponding space-time maps computed from Equations (8), from which we observe the expected local reflection/transmission in the co-polar components (Figure 3(A)) and the emergence (for t > τ) of a cross-polar response with clearly nonlocal character (Figure 3(B)). For a more quantitative assessment, Figure 3(C) and (D) show the corresponding spatial cuts at a fixed time instant (t = 10σ t ). In the co-polar response (Figure 3(C)), both the reflected (backward) and transmitted (forward) waveform are essentially scaled copies of the incident wavepacket (shown in the inset), as typically observed in conventional temporal boundaries [16]. Conversely, the cross-polar responses (Figure 3(D)) contain scaled copies of the first derivatives. This is the first example of an elementary analog operation that is performed only on a selected wave spin. We also notice that, by comparison with isotropic scenarios [43], now the analog-computing capabilities are enabled in transmission too, and without the need for impedance matching. This latter condition implies the possibility to perform analog computing in conjunction with frequency conversion, which is not attainable in the isotropic case [43].

Figure 3: Example of spin-dependent analog computing. Anisotropic short-pulsed temporal slab with ɛ ⊥ = 1, ɛ ‖ = 3, ɛ 1 = 1, ɛ 2 = 4, and τ = 0.5σ t , excited by the incident Gaussian wavepacket in Equation (9) [see inset in panel (C)]. (A), (B) Space-time maps [normalized electric induction, computed from Equations (8)], for co-polar and cross-polar responses, respectively. The thick purple-dashed lines indicate the temporal boundaries. (C), (D) Corresponding spatial cuts at t = 10σ t , computed via full-wave simulations. The superposed blue-dashed curves indicate the expected first derivatives.
Figure 3:

Example of spin-dependent analog computing. Anisotropic short-pulsed temporal slab with ɛ ⊥ = 1, ɛ ‖ = 3, ɛ 1 = 1, ɛ 2 = 4, and τ = 0.5σ t , excited by the incident Gaussian wavepacket in Equation (9) [see inset in panel (C)]. (A), (B) Space-time maps [normalized electric induction, computed from Equations (8)], for co-polar and cross-polar responses, respectively. The thick purple-dashed lines indicate the temporal boundaries. (C), (D) Corresponding spatial cuts at t = 10σ t , computed via full-wave simulations. The superposed blue-dashed curves indicate the expected first derivatives.

Figure 4 illustrates another interesting example, where the parameters are selected so as to attain impedance matching (ɛ 1 = ɛ 2). This implies the vanishing of the local term in the co-polar reflection response (see Equation (8c)), which is therefore dominated by nonlocality. As a consequence, we obtain a different combination of spin-dependent analog operations, where a first derivative is performed for both wave spins in reflection (Figure 4(A) and (C)), and for one only in transmission (Figure 4(B) and (D)).

Figure 4: Example of spin-dependent analog computing. Anisotropic short-pulsed temporal slab with ɛ ⊥ = 1, ɛ ‖ = 3, ɛ 1 = ɛ 2 = 1, and τ = 0.5σ t , excited by the incident Gaussian wavepacket in Equation (9) [see upper inset in panel (C)]. (A), (B) Space-time maps [normalized electric induction, computed from Equations (8)], for co-polar and cross-polar responses, respectively. The thick purple-dashed lines indicate the temporal boundaries. (C), (D) Corresponding spatial cuts at t = 10σ t , computed via full-wave simulations. The superposed blue-dashed curves indicate the expected first derivatives. The lower inset in panel (C) shows a magnified view of the reflection (backward) response.
Figure 4:

Example of spin-dependent analog computing. Anisotropic short-pulsed temporal slab with ɛ ⊥ = 1, ɛ ‖ = 3, ɛ 1 = ɛ 2 = 1, and τ = 0.5σ t , excited by the incident Gaussian wavepacket in Equation (9) [see upper inset in panel (C)]. (A), (B) Space-time maps [normalized electric induction, computed from Equations (8)], for co-polar and cross-polar responses, respectively. The thick purple-dashed lines indicate the temporal boundaries. (C), (D) Corresponding spatial cuts at t = 10σ t , computed via full-wave simulations. The superposed blue-dashed curves indicate the expected first derivatives. The lower inset in panel (C) shows a magnified view of the reflection (backward) response.

As a further variation, in Figure 5, we select the parameters in such a way that both the local and first-order nonlocal terms in Equation (8c) vanish. It can be shown that this condition also implies the vanishing of the second-order nonlocality, thereby leaving the third-order term ∝ i k x 3 as the dominant one (see the Methods Section 4.1 for details). This enables a more sophisticated response, where a first- or third-order derivative is performed in reflection, depending on the wave spin. As also observed in previous studies on short-pulsed isotropic metamaterials [43, 44], by increasing the order of the derivatives, their amplitude may decrease rapidly. However, it is worth stressing that the above parameters were merely chosen for a basic illustration of the phenomenon, and the maximization of the amplitude was not a concern; in principle, higher efficiencies may be obtained, also in view of the inherently active character of our time-varying platform.

Figure 5: Example of spin-dependent analog computing. Anisotropic short-pulsed temporal slab with ɛ ⊥ = 1, ɛ ‖ = 3, ɛ 1 = ɛ 2 = 1.5, and τ = 0.5σ t , excited by the incident Gaussian wavepacket in Equation (9) [see upper inset in panel (C)]. (A), (B) Space-time maps [normalized electric induction, computed from Equations (8)], for co-polar and cross-polar responses, respectively. The thick purple-dashed lines indicate the temporal boundaries, and the color scale in panel (A) is suitably saturated so as to show the weakest waveform. (C), (D) Corresponding spatial cuts at t = 10σ t , computed via full-wave simulations. The superposed blue-dashed and cyan-dotted curves indicate the expected first and third derivatives, respectively. The lower inset in panel (C) shows a magnified view of the reflection (backward) response.
Figure 5:

Example of spin-dependent analog computing. Anisotropic short-pulsed temporal slab with ɛ ⊥ = 1, ɛ ‖ = 3, ɛ 1 = ɛ 2 = 1.5, and τ = 0.5σ t , excited by the incident Gaussian wavepacket in Equation (9) [see upper inset in panel (C)]. (A), (B) Space-time maps [normalized electric induction, computed from Equations (8)], for co-polar and cross-polar responses, respectively. The thick purple-dashed lines indicate the temporal boundaries, and the color scale in panel (A) is suitably saturated so as to show the weakest waveform. (C), (D) Corresponding spatial cuts at t = 10σ t , computed via full-wave simulations. The superposed blue-dashed and cyan-dotted curves indicate the expected first and third derivatives, respectively. The lower inset in panel (C) shows a magnified view of the reflection (backward) response.

It is worth highlighting that more sophisticated operations can be attained by tailoring the short-pulsed modulation waveform and/or via multiple, time-resolved short-pulsed temporal slabs, by extending to the anisotropic case of interest here the approaches developed in Refs. [43, 44] for isotropic scenarios.

For a basic illustration, as shown in Figure 6(A), we consider a scenario featuring two identical anisotropic short-pulsed temporal slabs with parameters as in Figure 5. Similar to the isotropic case [44], as an effect of the multiple interactions, we now observe two time-resolved waveforms in the reflection and transmission responses, with the presence of composed operations (second derivatives). However, as an effect of the anisotropy, we now obtain both the co-polar and cross-polar responses (see Figure 6(B) and (C)). This enables the computation of first- or second-order derivatives, depending on the wave spin.

Figure 6: Example of spin-dependent analog computing. (A) Two anisotropic short-pulsed temporal slabs with ɛ ⊥ = 1, ɛ ‖ = 3, ɛ 1 = ɛ 2 = 1.5, and τ = 0.5σ t , starting at t = 0 and t = 10σ t , excited by the incident Gaussian wavepacket in Equation (9). (B), (C) Spatial cuts of normalized electric induction at t = 40σ t , computed via full-wave simulations, for co-polar and cross-polar responses, respectively. The insets show some magnified views of the responses. The superposed blue-dashed and green-dotted curves indicate the expected first and second derivatives, respectively.
Figure 6:

Example of spin-dependent analog computing. (A) Two anisotropic short-pulsed temporal slabs with ɛ ⊥ = 1, ɛ ‖ = 3, ɛ 1 = ɛ 2 = 1.5, and τ = 0.5σ t , starting at t = 0 and t = 10σ t , excited by the incident Gaussian wavepacket in Equation (9). (B), (C) Spatial cuts of normalized electric induction at t = 40σ t , computed via full-wave simulations, for co-polar and cross-polar responses, respectively. The insets show some magnified views of the responses. The superposed blue-dashed and green-dotted curves indicate the expected first and second derivatives, respectively.

2.3.2 Spin–orbit interactions: vortex generation

As a second representative example, we consider the generation of optical vortices. This phenomenon falls under the category of spin–orbit interactions of light [47], which have been extensively studied in recent years due to their potential applications in various fields. In particular, vortices have been proposed for particle trapping [48], optical communications [49], quantum technologies [50], and microscopy imaging [51].

However, in temporal metamaterials, spin–orbit interactions remain hitherto largely unexplored. Recently, a temporal spin-Hall effect, i.e., a spin-dependent frequency shift, has been theoretically demonstrated at a temporal boundary between bianisotropic chiral and dielectric media [42]. This effect is the temporal analog of the spin-Hall effect of light, where a circularly polarized beam experiences a spin-dependent spatial shift [47].

In conventional spatial scenarios, optical vortices can be generated through the reflection from a slab [45] and/or the propagation along the optical axis of a homogeneous uniaxial crystal [52]. In both configurations, the spin–orbit interaction effect essentially stems from the difference between the dynamics of transverse electric and magnetic fields [45]. It appears therefore suggestive to investigate similar effects in time-varying scenarios, since temporal anisotropy likewise induces this type of asymmetry.

To investigate optical vortex generation, we now assume a beam propagating along the optical axis (i.e., z) of the temporal anisotropic slab, synthesized through the superposition of plane waves with wavevector k 0 = k 0 ⊥ cos ⁡ θ e ̂ x + sin ⁡ θ e ̂ y + k 0 z e ̂ z , where θ = 0,2 π , and k 0⊥, k 0z are pre-set wavenumbers. Accordingly, the incident beam can be written as

(10a) D ( i ) ( r , t ) = Re e i k 0 z z − c n 1 k 0 t D ̃ ( i ) ( r ⊥ ) ,

where

(10b) D ̃ ( i ) ( r ⊥ ) = ∫ 0 2 π d θ e i k 0 ⊥ r ⊥ ⁡ cos ( θ − ϕ ) d p ( θ ) e ̂ p + d s ( θ ) e ̂ s ,

k 0 = k 0 ⊥ 2 + k 0 z 2 , d p,s (θ) is the spectral angle of the p/s-polarized waves, and we have assumed polar coordinates r ⊥ = r ⊥ ( cos ⁡ ϕ e ̂ x + sin ⁡ ϕ e ̂ y ) in the transverse plane. Equation (10b) is the spectral representation of a monochromatic non-diffracting wavefield, which can be expressed as a superposition of different Bessel beams [53]. In essence, a Bessel beam does not experience diffraction since it is the result of the superposition of plane waves exhibiting the same angular frequency and longitudinal wavenumber, and can be written as the product between a plane-wave carrier and a transverse profile.

By exploiting the formalism in Equations (7), we can work out analytically the expressions of the transmitted (forward) and reflected (backward) beams (see the Methods Section 4.1 for details). In particular, by assuming the circularly polarized basis e ̂ ± = ( e ̂ x ± i e ̂ y ) / 2 , the transverse components of the incident electric induction can be written as

(11) D ̃ ⊥ ( i ) ( r ⊥ ) = ∫ 0 2 π d θ e i k 0 ⊥ r ⊥ ⁡ cos ( θ − ϕ ) × U + ( i ) e − i θ e ̂ + + U − ( i ) e i θ e ̂ − ,

where we have defined U ± ( i ) = d p k 0 z / k 0 ∓ i d s / 2 . Accordingly, we obtain for the transverse components of the transmitted and reflected electric inductions D ⊥ ( t ) ( r , t ) = Re e i k 0 z z − c n 2 k 0 ( t − τ ) D ̃ ⊥ ( t ) ( r ⊥ ) and D ⊥ ( r ) ( r , t ) = Re e i k 0 z z + c n 2 k 0 ( t − τ ) D ̃ ⊥ ( r ) ( r ⊥ ) , respectively, with

(12) D ̃ ⊥ ( j ) ( r ⊥ ) = ∫ 0 2 π d θ e i k 0 ⊥ r ⊥ ⁡ cos ( θ − ϕ ) × U + ( j ) e − i θ e ̂ + + U − ( j ) e i θ e ̂ − ,

and j = r, t. Here,

(13a) U ± ( t ) = 1 2 T p p k 0 ⊥ , k 0 z d p k 0 z k 0 ∓ i T s s k 0 ⊥ , k 0 z d s ,

(13b) U ± ( r ) = 1 2 R p p k 0 ⊥ , k 0 z d p k 0 z k 0 ∓ i R s s k 0 ⊥ , k 0 z d s ,

with the dependence on k 0⊥ and k 0z (rather than on k 0) in the transmission/reflection coefficients being a consequence of the axial symmetry of the considered beam and anisotropic temporal slab.

As an example, we assume that the incident beam has positive spin (LHC) and is devoid of topological charge, i.e.,

(14) U + ( i ) = D 0 2 π e i θ , U − ( i ) = 0 ,

with D 0 denoting a real-valued normalization constant; this implies that the spectral amplitudes fulfill the relationships d p = k 0 D 0 e i θ / ( 2 π k 0 z 2 ) , d s = i D 0 e i θ / ( 2 π 2 ) . By substituting these assumptions in Equations (11), and performing the angular integration, we obtain

(15a) D ̃ ⊥ ( i ) ( r ⊥ ) = D 0 J 0 ( k 0 ⊥ r ⊥ ) e ̂ + ,

(15b) D ̃ ⊥ ( t ) ( r ⊥ ) = D 0 T + + k 0 ⊥ , k 0 z J 0 ( k 0 ⊥ r ⊥ ) e ̂ + − T − + × k 0 ⊥ , k 0 z J 2 ( k 0 ⊥ r ⊥ ) e i 2 ϕ e ̂ − ,

(15c) D ̃ ⊥ ( r ) ( r ⊥ ) = D 0 R + + k 0 ⊥ , k 0 z J 0 ( k 0 ⊥ r ⊥ ) e ̂ + − R − + × k 0 ⊥ , k 0 z J 2 ( k 0 ⊥ r ⊥ ) e i 2 ϕ e ̂ − ,

where J m (⋅) denote the mth-order Bessel functions [54], and T ±±, R ±± (T ∓±, R ∓±) are co-polarized (cross-polarized) transmission and reflection coefficients, respectively, for the chosen circularly polarized basis e ̂ ± = ( e ̂ x ± i e ̂ y ) / 2 .

Equations (15) clearly illustrate the spin–orbit interaction effect that can occur in an anisotropic temporal slab, with the polarization conversion giving rise to a variation of the orbital angular momentum. It appears thus possible to generate a vortex beam with topological charge ℓ = 2 in the cross-polarized reflection or transmission response. For instance, as schematically illustrated in Figure 7, it is possible to tailor the anisotropic slab parameters so that R ++ = T ++ = 0 (see the Methods Section 4.1 for details), i.e., so that an impinging Bessel-type beam with a positive spin and no topological charge (ℓ = 0, Figure 7(A)) is converted into reflected and transmitted cross-polarized vortex beams with topological charge ℓ = 2 (Figure 7(B)).

Figure 7: Schematic illustration of the spin–orbit interactions and vortex generation. (A) A circularly polarized Bessel-type beam with positive spin (i.e., LHC) and no topological charge (ℓ = 0) impinges in the initial medium. (B) After the anisotropic temporal modulation (t > τ), the beam is generally converted into cross-polarized reflected (backward) and transmitted (forward) vortex beams with topological charge ℓ = 2. If the initial and final permittivities are different (ɛ 1 ≠ ɛ 2), frequency conversion (from ω 1 to ω 2) is attained too. Note the different incidence conditions by comparison with the scenario in Figure 2.
Figure 7:

Schematic illustration of the spin–orbit interactions and vortex generation. (A) A circularly polarized Bessel-type beam with positive spin (i.e., LHC) and no topological charge (ℓ = 0) impinges in the initial medium. (B) After the anisotropic temporal modulation (t > τ), the beam is generally converted into cross-polarized reflected (backward) and transmitted (forward) vortex beams with topological charge ℓ = 2. If the initial and final permittivities are different (ɛ 1 ≠ ɛ 2), frequency conversion (from ω 1 to ω 2) is attained too. Note the different incidence conditions by comparison with the scenario in Figure 2.

The above mechanism is similar to that occurring in conventional spatial scenarios involving uniaxial crystals [52], with the important difference that the temporal configuration is not necessarily bound by power conservation for EM signals, and the vortex-beam generation is also accompanied by frequency conversion if the initial and final permittivities do not coincide (ɛ 1 ≠ ɛ 2). Specifically, by calculating the time-averaged power flow of the beams before and after the temporal slab, we obtain (see the Methods Section 4.2 for further details)

(16) P + ( t ) + P − ( t ) + P + ( r ) + P − ( r ) P + ( i ) = n 1 n 2 2 ,

where, as usual, the superscripts i, r, t denote the incident, reflected and transmitted beams, respectively, whereas the subscripts +, − denote positive and negative spin, respectively. Similar to the scenario of a single temporal boundary, the power conservation holds only for the impedance-matching case (n 1 = n 2) [55]. Furthermore, it is natural to define the vortex generation efficiency η as the fraction of the incident power coupled with the generated vortex beam (i.e., η = P − ( t ) / P + ( i ) ) [45]. In the assumed conditions where R ++ = T ++ = 0, we obtain (see the Methods Section 4.2 for further details)

(17) η = n 1 n 2 n 1 n 2 + | R − + | 2 .

From Equation (17), it is evident that the vortex generation efficiency can become greater than one as a consequence of the lack of power conservation. In the impedance-matching scenario (n 1 = n 2), it is also possible to attain R −+ = 0 (see the Methods Section 4.1 for details), so that the impinging beam is perfectly converted into a cross-polarized transmitted vortex beam with unit efficiency and without frequency conversion.

This latter scenario is exemplified in Figure 8, with parameters (given in the caption) tailored so as to perfectly convert an impinging Bessel-type beam devoid of topological charge (i.e., ℓ = 0) into a cross-polarized transmitted vortex beam with topological charge ℓ = 2. Figure 8(A) and (B) show the full-wave computed transverse wavefronts (at a fixed time) of the impinging and cross-polarized transmitted beams (with all other scattering terms being below ∼ 1 0 − 7 in the normalized scale, and not shown for brevity), which confirm our theoretical predictions.

Figure 8: Example of spin–orbit interaction effects (vortex generation). Anisotropic temporal slab with ɛ ⊥ = 8, ɛ ‖ = 1.946, ɛ 1 = ɛ 2 = 1 and τ = 3 2 T $\tau =3\sqrt{2}T$ , excited by a time-harmonic Bessel-type beam as in Equations (10), with period T, wavelength λ = cT, and characteristic wavenumbers k 0 = 2π/λ and k 0⊥ = π/λ. (A), (B) Transverse wavefronts (normalized electric induction) of impinging and transmitted (forward) cross-polarized beams, respectively, computed via full-wave simulations at t = −T and t = 19.75T, respectively. All other (co- and cross-polarized) scattering terms are below ∼ 1 0 − 7 ${\sim} 1{0}^{-7}$ in the normalized scale.
Figure 8:

Example of spin–orbit interaction effects (vortex generation). Anisotropic temporal slab with ɛ ⊥ = 8, ɛ ‖ = 1.946, ɛ 1 = ɛ 2 = 1 and τ = 3 2 T , excited by a time-harmonic Bessel-type beam as in Equations (10), with period T, wavelength λ = cT, and characteristic wavenumbers k 0 = 2π/λ and k 0⊥ = π/λ. (A), (B) Transverse wavefronts (normalized electric induction) of impinging and transmitted (forward) cross-polarized beams, respectively, computed via full-wave simulations at t = −T and t = 19.75T, respectively. All other (co- and cross-polarized) scattering terms are below ∼ 1 0 − 7 in the normalized scale.

3 Conclusions

To sum up, we have shown that temporal anisotropy, in the form of abrupt transitions from isotropic to anisotropic dielectric permittivity (and vice versa), can be harnessed to attain spin-controlled photonic operations. These include, for instance, spin-dependent analog computing on an impinging wavepacket and spin–orbit interaction effects for vortex generation. Overall, our outcomes indicate new pathways for advanced control of light–matter interactions, which may find potential applications in various scenarios, ranging from telecommunications to optical and quantum computing.

From the implementation viewpoint, temporal anisotropy poses technological challenges that are comparable to those encountered in isotropic temporal metamaterials (see, e.g., the discussions in Refs. [31–33] and the recent experimental results in Refs. [34–36]); the additional anisotropic character could be implemented via a suitable (asymmetric) design of the reconfigurable meta-atoms, and should not pose insurmountable obstacles.

As possible extensions and follow-up studies, we are currently investigating more general short-pulsed variations (by generalizing our previous results in Ref. [43] for isotropic scenarios) as well as the combination of temporal anisotropy and conventional spatial modulations. Also of great interest are the modeling of temporal dispersion and the study of possible applications of temporally induced spin-controlled photonics to quantum technologies.

4 Methods

4.1 Analytical modeling

For the anisotropic temporal slab described by Equations (2), the temporal scattering problem can be solved analytically in closed form. By particularizing to our scenario the approach proposed in Ref. [40], the temporal reflection and transmission coefficients for the ordinary (s-polarized) and extraordinary (p-polarized) plane waves can be expressed as

(18a) T j j ( k ) = 1 2 1 + ω 1 ω 2 cos ω j τ − i ω 1 ω j + ω j ω 2 sin ω j τ ,

(18b) R j j ( k ) = 1 2 1 − ω 1 ω 2 cos ω j τ − i ω 1 ω j − ω j ω 2 sin ω j τ ,

with j = p, s.

In a generic polarization basis, parameterized by the angles φ and δ,

(19a) e ̂ 1 = cos ⁡ φ e ̂ p + e i δ ⁡ sin ⁡ φ e ̂ s ,

(19b) e ̂ 2 = sin ⁡ φ e ̂ p − e i δ ⁡ cos ⁡ φ e ̂ s ,

the temporal transmission and reflection matrices can be written as

(20a) M ̲ ̲ = M 11 M 12 M 21 M 22 ,

where

(20b) M 11 = M pp ⁡ cos 2 ⁡ φ + M s s ⁡ sin 2 ⁡ φ ,

(20c) M 22 = M pp ⁡ sin 2 ⁡ φ + M s s ⁡ cos 2 ⁡ φ ,

(20d) M 12 = M 21 = M pp − M s s 2 sin ( 2 φ ) ,

with M = T, R.

In Section 2.3, we have considered circular polarized waves, i.e., e ̂ ± = ( e ̂ p ± i e ̂ s ) / 2 , corresponding to φ = π/4 and δ = π/2 in Equations (19). This yields,

(21a) T + + ( k ) = T − − ( k ) = T pp ( k ) + T s s ( k ) 2 ,

(21b) T + − ( k ) = T − + ( k ) = T pp ( k ) − T s s ( k ) 2 ,

(21c) R + + ( k ) = R − − ( k ) = R pp ( k ) + R s s ( k ) 2 ,

(21d) R + − ( k ) = R − + ( k ) = R pp ( k ) − R s s ( k ) 2 .

More specifically, in the example discussed in Section 2.3.1, we have assumed pulsed plane waves propagating along the x-axis, e ̂ ± = ( − e ̂ z ± i e ̂ y ) / 2 , corresponding to k = k x e ̂ x in Equations (5). From Equations (21), the expressions in Equations (8) readily follow via Maclaurin series expansions up to the second order in τ.

Moreover, with specific reference to the example in Figure 5, the Maclaurin series expansion of the coefficient R ++ up to the third order in τ can be written as

(22) R + + ( k x ) ≃ 1 2 1 − n 2 n 1 − i π 2 2 n 1 − n 2 n ⊥ 2 − n 2 n ‖ 2 k x K − π 2 2 × 1 − n 2 n 1 1 n ⊥ 2 + 1 n ‖ 2 k x 2 K 2 + i π 3 3 1 n 1 1 n ⊥ 2 + 1 n ‖ 2 − n 2 1 n ⊥ 4 + 1 n ‖ 4 × k x 3 K 3 + O k x 4 K 4 .

From the above expression, it appears evident that the parameter choice as in Figure 5 implies the vanishing of the terms up to the second order, thereby leaving the third order ∝ i k x 3 as the dominant one.

In Section 2.3.2, we have instead considered non-diffracting beams propagating along the z-axis. Accordingly, we have expressed the transverse component of the electric induction in the basis e ̂ ± = ( e ̂ x ± i e ̂ y ) / 2 . The considered beams are a suitable superposition of plane waves with wavevector k 0 = k 0 ⊥ cos ⁡ θ e ̂ x + sin ⁡ θ e ̂ y + k 0 z e ̂ z , and θ = 0,2 π . Consequently, the temporal transmission and reflected coefficients, T + ± k 0 ⊥ , k 0 z , R + ± k 0 ⊥ , k 0 z appearing in Equations (15) correspond to the matrix entries reported in Equations (21) with k = k 0; these terms do not depend on θ in view of the axial symmetry of the considered system.

From Equations (21), we observe that the vanishing of the co-polar response (T ++ = R ++ = 0) can be attained by enforcing T pp = −T ss and R pp = −R ss , which is obtained by selecting ω j τ = m j π (j = p, s), with m p and m s being odd and even positive integers, respectively, or vice versa. This yields

(23a) k 0 ⊥ 2 n ‖ 2 + k 0 z 2 n ⊥ 2 = π m p c τ 2 ,

(23b) k 0 ⊥ 2 + k 0 z 2 = n ⊥ 2 π m s c τ 2 ,

which are consistent with the parameters chosen in Figure 8 (with m p = 3 and m s = 4). If, as in Figure 8, we also assume impedance matching, i.e., ɛ 1 = ɛ 2, it can readily be verified from Equation (18b) (with ω 1 = ω 2) that R −+ = 0, i.e., the impinging Bessel beam is perfectly converted into a cross-polarized transmitted vortex beam, without frequency conversion.

4.2 Power flow in vortex generation

To assess the (lack of) power conservation, we evaluate the time-averaged power flowing through a disk of radius ρ, lying in the z = 0 plane and centered at the origin, for the Bessel-type beams in Equations (15). For the incident (i) and transmitted (t) beams, we obtain

(24) P ± ( j ) = π ∫ 0 ρ d r ⊥ r ⊥ e ̂ z ⋅ Re E ω , ± ( j ) × H ω , ± ( j ) * ,

where E ± ( j ) = Re ( E ω , ± ( j ) e − i ω t ) , H ± ( j ) = Re ( H ω , ± ( j ) e − i ω t ) are the electric and magnetic fields with positive/negative spin, and ω = ω 1, ω 2 for j = i, t, respectively. Similarly, the reflected power flow is

(25) P ± ( r ) = π ∫ 0 ρ d r ⊥ r ⊥ e ̂ z ⋅ Re E ω 2 , ± ( r ) × H ω 2 , ± ( r ) * ,

with E ± ( r ) = Re ( E ω 2 , ± ( r ) e i ω 2 t ) , H ± ( r ) = Re ( H ω 2 , ± ( r ) e i ω 2 t ) . By using Equations (15) along with Equations (1) and Maxwell’s curl equation ∇ × E = −∂ t B, we obtain

(26a) P ± ( t ) P + ( i ) = n 1 n 2 ξ 2 ξ 1 | T ± + k 0 ⊥ , k 0 z | 2 ,

(26b) P ± ( r ) P + ( i ) = − n 1 n 2 ξ 2 ξ 1 | R ± + k 0 ⊥ , k 0 z | 2 ,

where ξ m = ∫ 0 ρ d r ⊥ r ⊥ J m 2 ( k 0 ⊥ r ⊥ ) with m = 0, 2. Furthermore, by exploiting Equations (18) and (21), it is straightforward to verify that

(27) | T + + | 2 + | T − + | 2 − | R + + | 2 − | R − + | 2 = n 1 n 2 .

Recalling that ξ 2/ξ 1 → 1 in the limit ρ → ∞, and by combining Equations (26a) with Equations (27), we obtain Equation (16) reported in Section 2.3.2. From Equation (26a), we obtain for the vortex efficiency

(28) η = lim ρ → + ∞ P − ( t ) P + ( i ) = n 1 n 2 | T − + k 0 ⊥ , k 0 z | 2 ,

which, recalling Equation (27) with T ++ = R ++ = 0, becomes Equation (17) reported in Section 2.3.2.

4.3 Full-wave solution

Our full-wave solution is obtained via a rigorous numerical approach that generalizes the method introduced in Refs. [43, 56] for isotropic scenarios. In essence, the method operates on the spatially algebrized Maxwell’s curl equations

(29a) i k ε 0 × ε ̲ ̲ − 1 t ⋅ d k , t = − d b k , t d t ,

(29b) i k × b k , t = μ 0 d d k , t d t ,

i.e., a system of six ordinary differential equations in time, where the wavevector k is considered as a fixed parameter. For given initial conditions and wavevector, the above system is solved numerically by means of the NDSolve routine available in Mathematica TM [57], which applies adaptively several numerical methods (e.g., Runge–Kutta, predictor-corrector, implicit backward differentiation). In our implementation, we rely on default settings and parameters. Moreover, to favor numerical convergence and to assess the effects of non-ideal temporal boundaries, we model the abrupt changes in the permittivity by means of an analytical, smooth unit-step function U s (t) = [tanh(t/T s ) + 1]/2, with T s = 10−4 T.

For a given impinging wavefield [Equation (9) for the example in Section 2.3.1, and Equations (10) for the example in Section 2.3.2], we compute a time-dependent plane-wave spectrum, from which we derive a set of discretized wavevectors and corresponding initial conditions. For each wavevector, we then solve the system in Equations (29). From these numerical solutions, we finally synthesize the physical observable of interest (electric induction) in terms of a plane-wave spectral integral, numerically implemented via fast-Fourier-transform by means of the Fourier routine available in MathematicaTM [57].

Corresponding authors: Carlo Rizza, Department of Physical and Chemical Sciences, University of L’Aquila, L’Aquila, I-67100, Italy, E-mail: ; and Vincenzo Galdi, Department of Engineering , Fields & Waves Lab, University of Sannio, Benevento, I-82100, Italy, E-mail: .

Carlo Rizza and Giuseppe Castaldi contributed equally to this work.

Funding source: University of Sannio

Award Identifier / Grant number: FRA 2021

Funding source: Italian Ministry of University and Research (MUR)

Funding source: National Innovation Ecosystem

Funding source: Recovery and Resilience National Plan (PNRR) Italy

Funding source: Vitality

Award Identifier / Grant number: CUP D73C2200084000

  1. Author contributions: C. R., G. C., and V. G. conceived the idea. C. R. and G. C. developed the analytical framework. G. C. carried out the numerical simulations. V. G. supervised the study. All authors participated in the data analysis and interpretation. C. R. and V. G. wrote the manuscript, with inputs and feedback from all authors.

  2. Research funding: This work was supported by the Italian Ministry of University and Research (MUR), National Innovation Ecosystem, Recovery and Resilience National Plan (PNRR) Italy, Vitality, CUP D73C2200084000. G. C. and V. G. acknowledge partial support from the University of Sannio via the FRA 2021 Program.

  3. Conflict of interest statement: The authors declare no conflicts of interest regarding this article.

References

[1] N. Engheta, “Metamaterials with high degrees of freedom: space, time, and more,” Nanophotonics, vol. 10, no. 1, pp. 639–642, 2021, https://doi.org/10.1515/nanoph-2020-0414.Suche in Google Scholar

[2] C. Caloz and Z. Deck-Léger, “Spacetime metamaterials—part I: general concepts,” IEEE Trans. Antenn. Propag., vol. 68, no. 3, pp. 1569–1582, 2020, https://doi.org/10.1109/tap.2019.2944225.Suche in Google Scholar

[3] C. Caloz and Z. Deck-Léger, “Spacetime metamaterials—part II: theory and applications,” IEEE Trans. Antenn. Propag., vol. 68, no. 3, pp. 1583–1598, 2020, https://doi.org/10.1109/tap.2019.2944216.Suche in Google Scholar

[4] F. R. Morgenthaler, “Velocity modulation of electromagnetic waves,” IRE Trans. Microw. Theor. Tech., vol. 6, no. 2, pp. 167–172, 1958, https://doi.org/10.1109/tmtt.1958.1124533.Suche in Google Scholar

[5] A. A. Oliner and A. Hessel, “Wave propagation in a medium with a progressive sinusoidal disturbance,” IRE Trans. Microw. Theor. Tech., vol. 9, no. 4, pp. 337–343, 1961, https://doi.org/10.1109/tmtt.1961.1125340.Suche in Google Scholar

[6] L. Felsen and G. Whitman, “Wave propagation in time-varying media,” IEEE Trans. Antenn. Propag., vol. 18, no. 2, pp. 242–253, 1970, https://doi.org/10.1109/tap.1970.1139657.Suche in Google Scholar

[7] R. Fante, “Transmission of electromagnetic waves into time-varying media,” IEEE Trans. Antenn. Propag., vol. 19, no. 3, pp. 417–424, 1971, https://doi.org/10.1109/tap.1971.1139931.Suche in Google Scholar

[8] F. Wilczek, “Quantum time crystals,” Phys. Rev. Lett., vol. 109, no. 16, p. 160401, 2012, https://doi.org/10.1103/physrevlett.109.160401.Suche in Google Scholar PubMed

[9] K. Sacha and J. Zakrzewski, “Time crystals: a review,” Rep. Prog. Phys., vol. 81, no. 1, p. 016401, 2017, https://doi.org/10.1088/1361-6633/aa8b38.Suche in Google Scholar PubMed

[10] N. Kamaraju, A. Rubano, L. Jian, et al.., “Subcycle control of terahertz waveform polarization using all-optically induced transient metamaterials,” Light Sci. Appl., vol. 3, no. 2, p. e155, 2014, https://doi.org/10.1038/lsa.2014.36.Suche in Google Scholar

[11] M. Z. Alam, I. De Leon, and R. W. Boyd, “Large optical nonlinearity of indium tin oxide in its epsilon-near-zero region,” Science, vol. 352, no. 6287, pp. 795–797, 2016, https://doi.org/10.1126/science.aae0330.Suche in Google Scholar PubMed

[12] A. Kord, M. Tymchenko, D. L. Sounas, H. Krishnaswamy, and A. Alù, “CMOS integrated magnetless circulators based on spatiotemporal modulation angular-momentum biasing,” IEEE Trans. Microw. Theor. Tech., vol. 67, no. 7, pp. 2649–2662, 2019, https://doi.org/10.1109/tmtt.2019.2915074.Suche in Google Scholar

[13] J. B. Khurgin, M. Clerici, V. Bruno, et al.., “Adiabatic frequency shifting in epsilon-near-zero materials: the role of group velocity,” Optica, vol. 7, no. 3, pp. 226–231, 2020, https://doi.org/10.1364/optica.374788.Suche in Google Scholar

[14] Y. Zhou, M. Z. Alam, M. Karimi, et al.., “Broadband frequency translation through time refraction in an epsilon-near-zero material,” Nat. Commun., vol. 11, no. 1, p. 2180, 2020, https://doi.org/10.1038/s41467-020-15682-2.Suche in Google Scholar PubMed PubMed Central

[15] Z. Hayran and F. Monticone, “Challenging fundamental limitations in electromagnetics with time-varying systems,” arXiv:2205.07142, 2022.Suche in Google Scholar

[16] Y. Xiao, D. N. Maywar, and G. P. Agrawal, “Reflection and transmission of electromagnetic waves at a temporal boundary,” Opt. Lett., vol. 39, no. 3, pp. 574–577, 2014, https://doi.org/10.1364/ol.39.000574.Suche in Google Scholar PubMed

[17] L. Stefanini, S. Yin, D. Ramaccia, A. Alù, A. Toscano, and F. Bilotti, “Temporal interfaces by instantaneously varying boundary conditions,” Phys. Rev. B, vol. 106, no. 9, p. 094312, 2022, https://doi.org/10.1103/physrevb.106.094312.Suche in Google Scholar

[18] D. Ramaccia, A. Toscano, and F. Bilotti, “Light propagation through metamaterial temporal slabs: reflection, refraction, and special cases,” Opt. Lett., vol. 45, no. 20, pp. 5836–5839, 2020, https://doi.org/10.1364/ol.402856.Suche in Google Scholar PubMed

[19] E. Galiffi, Y. T. Wang, Z. Lim, J. B. Pendry, A. Alù, and P. A. Huidobro, “Wood anomalies and surface-wave excitation with a time grating,” Phys. Rev. Lett., vol. 125, no. 12, p. 127403, 2020, https://doi.org/10.1103/physrevlett.125.127403.Suche in Google Scholar PubMed

[20] S. Taravati and G. V. Eleftheriades, “Generalized space-time-periodic diffraction gratings: theory and applications,” Phys. Rev. Appl., vol. 12, no. 2, p. 024026, 2019, https://doi.org/10.1103/physrevapplied.12.024026.Suche in Google Scholar

[21] D. Ramaccia, A. Alù, A. Toscano, and F. Bilotti, “Temporal multilayer structures for designing higher-order transfer functions using time-varying metamaterials,” Appl. Phys. Lett., vol. 118, no. 10, p. 101901, 2021, https://doi.org/10.1063/5.0042567.Suche in Google Scholar

[22] G. Castaldi, M. Moccia, N. Engheta, and V. Galdi, “Herpin equivalence in temporal metamaterials,” Nanophotonics, vol. 11, no. 20, pp. 4479–4488, 2022, https://doi.org/10.1515/nanoph-2022-0338.Suche in Google Scholar

[23] J. S. Martínez-Romero, O. M. Becerra-Fuentes, and P. Halevi, “Temporal photonic crystals with modulations of both permittivity and permeability,” Phys. Rev. A, vol. 93, no. 6, p. 063813, 2016, https://doi.org/10.1103/physreva.93.063813.Suche in Google Scholar

[24] E. Lustig, Y. Sharabi, and M. Segev, “Topological aspects of photonic time crystals,” Optica, vol. 5, no. 11, pp. 1390–1395, 2018, https://doi.org/10.1364/optica.5.001390.Suche in Google Scholar

[25] M. Lyubarov, Y. Lumer, A. Dikopoltsev, E. Lustig, Y. Sharabi, and M. Segev, “Amplified emission and lasing in photonic time crystals,” Science, vol. 377, no. 6604, pp. 425–428, 2022, https://doi.org/10.1126/science.abo3324.Suche in Google Scholar PubMed

[26] A. Shlivinski and Y. Hadad, “Beyond the Bode–Fano bound: wideband impedance matching for short pulses using temporal switching of transmission-line parameters,” Phys. Rev. Lett., vol. 121, no. 20, p. 204301, 2018, https://doi.org/10.1103/physrevlett.121.204301.Suche in Google Scholar

[27] V. Pacheco-Peña and N. Engheta, “Antireflection temporal coatings,” Optica, vol. 7, no. 4, pp. 323–331, 2020, https://doi.org/10.1364/optica.381175.Suche in Google Scholar

[28] G. Castaldi, V. Pacheco-Peña, M. Moccia, N. Engheta, and V. Galdi, “Exploiting space-time duality in the synthesis of impedance transformers via temporal metamaterials,” Nanophotonics, vol. 10, no. 14, pp. 3687–3699, 2021, https://doi.org/10.1515/nanoph-2021-0231.Suche in Google Scholar

[29] E. Galiffi, S. Yin, and A. Alù, “Tapered photonic switching,” Nanophotonics, vol. 11, no. 16, pp. 3575–3581, 2022, https://doi.org/10.1515/nanoph-2022-0200.Suche in Google Scholar

[30] H. Li and A. Alù, “Temporal switching to extend the bandwidth of thin absorbers,” Optica, vol. 8, no. 1, pp. 24–29, 2021, https://doi.org/10.1364/optica.408399.Suche in Google Scholar

[31] E. Galiffi, R. Tirole, S. Yin, et al.., “Photonics of time-varying media,” Adv. Photon., vol. 4, no. 1, p. 014002, 2022, https://doi.org/10.1117/1.ap.4.1.014002.Suche in Google Scholar

[32] V. Pacheco-Peña, D. M. Solís, and N. Engheta, “Time-varying electromagnetic media: opinion,” Opt. Mater. Express, vol. 12, no. 10, p. 3829, 2022, https://doi.org/10.1364/ome.471007.Suche in Google Scholar

[33] Z. Hayran, J. B. Khurgin, and F. Monticone, “ℏω versus ℏk: dispersion and energy constraints on time-varying photonic materials and time crystals,” Opt. Mater. Express, vol. 12, no. 10, pp. 3904–3917, 2022, https://doi.org/10.1364/ome.471672.Suche in Google Scholar

[34] H. Moussa, G. Xu, S. Yin, E. Galiffi, Y. Radi, and A. Alù, “Observation of temporal reflections and broadband frequency translations at photonic time-interfaces,” arXiv:2208.07236, 2022.10.1038/s41567-023-01975-ySuche in Google Scholar

[35] X. Wang, M. S. Mirmoosa, V. S. Asadchy, C. Rockstuhl, S. Fan, and S. A. Tretyakov, “Metasurface-based realization of photonic time crystals,” arXiv:2208.07231, 2022.10.1126/sciadv.adg7541Suche in Google Scholar PubMed PubMed Central

[36] T. Liu, J. Y. Ou, K. F. MacDonald, and N. I. Zheludev, “Photonic analogue of a continuous time crystal,” arXiv:2209.00324, 2022.Suche in Google Scholar

[37] V. Pacheco-Peña and N. Engheta, “Temporal aiming,” Light Sci. Appl., vol. 9, no. 1, p. 129, 2020, https://doi.org/10.1038/s41377-020-00360-1.Suche in Google Scholar PubMed PubMed Central

[38] V. Pacheco-Peña and N. Engheta, “Temporal equivalent of the Brewster angle,” Phys. Rev. B, vol. 104, no. 21, p. 214308, 2021, https://doi.org/10.1103/physrevb.104.214308.Suche in Google Scholar

[39] V. Pacheco-Peña and N. Engheta, “Spatiotemporal isotropic-to-anisotropic meta-atoms,” New J. Phys., vol. 23, no. 9, p. 095006, 2021, https://doi.org/10.1088/1367-2630/ac21df.Suche in Google Scholar

[40] J. Xu, W. Mai, and D. H. Werner, “Complete polarization conversion using anisotropic temporal slabs,” Opt. Lett., vol. 46, no. 6, pp. 1373–1376, 2021, https://doi.org/10.1364/ol.415757.Suche in Google Scholar

[41] H. Li, S. Yin, and A. Alù, “Nonreciprocity and Faraday rotation at time interfaces,” Phys. Rev. Lett., vol. 128, no. 17, p. 173901, 2022, https://doi.org/10.1103/physrevlett.128.173901.Suche in Google Scholar

[42] M. H. Mostafa, M. S. Mirmoosa, and S. A. Tretyakov, “Spin-temporal interactions of light,” arXiv:2204.01574, 2022.Suche in Google Scholar

[43] C. Rizza, G. Castaldi, and V. Galdi, “Short-pulsed metamaterials,” Phys. Rev. Lett., vol. 128, no. 25, p. 257402, 2022, https://doi.org/10.1103/physrevlett.128.257402.Suche in Google Scholar PubMed

[44] G. Castaldi, C. Rizza, N. Engheta, and V. Galdi, “Multiple actions of time-resolved short-pulsed metamaterials,” Appl. Phys. Lett., vol. 122, no. 2, p. 021701, 2023, https://doi.org/10.1063/5.0132554.Suche in Google Scholar

[45] A. Ciattoni, A. Marini, and C. Rizza, “Efficient vortex generation in subwavelength epsilon-near-zero slabs,” Phys. Rev. Lett., vol. 118, no. 10, p. 104301, 2017, https://doi.org/10.1103/physrevlett.118.104301.Suche in Google Scholar

[46] D. M. Solís, R. Kastner, and N. Engheta, “Time-varying materials in the presence of dispersion: plane-wave propagation in a Lorentzian medium with temporal discontinuity,” Photon. Res., vol. 9, no. 9, pp. 1842–1853, 2021, https://doi.org/10.1364/prj.427368.Suche in Google Scholar

[47] K. Y. Bliokh, F. J. Rodríguez-Fortuño, F. Nori, and A. V. Zayats, “Spin-orbit interactions of light,” Nat. Photon., vol. 9, no. 12, pp. 796–808, 2015, https://doi.org/10.1038/nphoton.2015.201.Suche in Google Scholar

[48] K. T. Gahagan and G. A. Swartzlander, “Optical vortex trapping of particles,” Opt. Lett., vol. 21, no. 11, pp. 827–829, 1996, https://doi.org/10.1364/ol.21.000827.Suche in Google Scholar PubMed

[49] Z. Wang, N. Zhang, and X.-C. Yuan, “High-volume optical vortex multiplexing and de-multiplexing for free-space optical communication,” Opt. Express, vol. 19, no. 2, pp. 482–492, 2011, https://doi.org/10.1364/oe.19.000482.Suche in Google Scholar

[50] G. Molina-Terriza, J. P. Torres, and L. Torner, “Management of the angular momentum of light: preparation of photons in multidimensional vector states of angular momentum,” Phys. Rev. Lett., vol. 88, no. 1, p. 013601, 2001, https://doi.org/10.1103/physrevlett.88.013601.Suche in Google Scholar PubMed

[51] L. Li and F. Li, “Beating the Rayleigh limit: orbital-angular-momentum-based super-resolution diffraction tomography,” Phys. Rev. E, vol. 88, no. 3, p. 033205, 2013, https://doi.org/10.1103/physreve.88.033205.Suche in Google Scholar PubMed

[52] X. Ling, H. Luo, F. Guan, X. Zhou, H. Luo, and L. Zhou, “Vortex generation in the spin-orbit interaction of a light beam propagating inside a uniaxial medium: origin and efficiency,” Opt. Express, vol. 28, no. 19, pp. 27258–27267, 2020, https://doi.org/10.1364/oe.403650.Suche in Google Scholar

[53] J. Durnin, J. J. Miceli, and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett., vol. 58, no. 15, pp. 1499–1501, 1987, https://doi.org/10.1103/physrevlett.58.1499.Suche in Google Scholar PubMed

[54] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions: With Formulas Graphs, and Mathematical Tables, New York, NY, Dover, 1965.10.1063/1.3047921Suche in Google Scholar

[55] W. Mai, J. Xu, and D. H. Werner, “Fundamental asymmetries between spatial and temporal boundaries in electromagnetics,” arXiv:2207.04286, 2022.10.3390/sym15040858Suche in Google Scholar

[56] C. Rizza, G. Castaldi, and V. Galdi, “Nonlocal effects in temporal metamaterials,” Nanophotonics, vol. 11, no. 7, pp. 1285–1295, 2022, https://doi.org/10.1515/nanoph-2021-0605.Suche in Google Scholar

[57] Wolfram Research, Inc., Mathematica, Version 12.3.1, Champaign, IL, 2021.Suche in Google Scholar
Received: 2022-12-28
Accepted: 2023-02-10
Published Online: 2023-02-21

© 2023 the author(s), published by De Gruyter, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 International License.

Artikel in diesem Heft

  1. Frontmatter
  2. Editorial
  3. Nanophotonics in support of Ukrainian Scientists
  4. Reviews
  5. Asymmetric transmission in nanophotonics
  6. Integrated circuits based on broadband pixel-array metasurfaces for generating data-carrying optical and THz orbital angular momentum beams
  7. Singular optics empowered by engineered optical materials
  8. Electrochemical photonics: a pathway towards electrovariable optical metamaterials
  9. Sustainable chemistry with plasmonic photocatalysts
  10. Perspectives
  11. Ukraine and singular optics
  12. Machine learning to optimize additive manufacturing for visible photonics
  13. Through thick and thin: how optical cavities control spin
  14. Research Articles
  15. Spin–orbit coupling induced by ascorbic acid crystals
  16. Broadband transfer of binary images via optically long wire media
  17. Counting and mapping of subwavelength nanoparticles from a single shot scattering pattern
  18. Controlling surface waves with temporal discontinuities of metasurfaces
  19. On the relation between electrical and electro-optical properties of tunnelling injection quantum dot lasers
  20. On-chip multivariant COVID 19 photonic sensor based on silicon nitride double-microring resonators
  21. Nano-infrared imaging of metal insulator transition in few-layer 1T-TaS2
  22. Electrical generation of surface phonon polaritons
  23. Dynamic beam control based on electrically switchable nanogratings from conducting polymers
  24. Tilting light’s polarization plane to spatially separate the ultrafast nonlinear response of chiral molecules
  25. Spin-dependent phenomena at chiral temporal interfaces
  26. Spin-controlled photonics via temporal anisotropy
  27. Coherent control of symmetry breaking in transverse-field Ising chains using few-cycle pulses
  28. Field enhancement of epsilon-near-zero modes in realistic ultrathin absorbing films
  29. Controlled compression, amplification and frequency up-conversion of optical pulses by media with time-dependent refractive index
  30. Tailored thermal emission in bulk calcite through optic axis reorientation
  31. Tip-enhanced photoluminescence of monolayer MoS2 increased and spectrally shifted by injection of electrons
  32. Quantum-enhanced interferometer using Kerr squeezing
  33. Nonlocal electro-optic metasurfaces for free-space light modulation
  34. Dispersion braiding and band knots in plasmonic arrays with broken symmetries
  35. Dual-mode hyperbolicity, supercanalization, and leakage in self-complementary metasurfaces
  36. Monocular depth sensing using metalens
  37. Multimode hybrid gold-silicon nanoantennas for tailored nanoscale optical confinement
  38. Replicating physical motion with Minkowskian isorefractive spacetime crystals
  39. Reconfigurable nonlinear optical element using tunable couplers and inverse-designed structure
https://doi.org/10.1515/nanoph-2022-0809
Schlagwörter für diesen Artikel
analog computing; anisotropy; metamaterials; spin-orbit interaction; time-varying
Creative Commons
BY 4.0

Artikel in diesem Heft

  1. Frontmatter
  2. Editorial
  3. Nanophotonics in support of Ukrainian Scientists
  4. Reviews
  5. Asymmetric transmission in nanophotonics
  6. Integrated circuits based on broadband pixel-array metasurfaces for generating data-carrying optical and THz orbital angular momentum beams
  7. Singular optics empowered by engineered optical materials
  8. Electrochemical photonics: a pathway towards electrovariable optical metamaterials
  9. Sustainable chemistry with plasmonic photocatalysts
  10. Perspectives
  11. Ukraine and singular optics
  12. Machine learning to optimize additive manufacturing for visible photonics
  13. Through thick and thin: how optical cavities control spin
  14. Research Articles
  15. Spin–orbit coupling induced by ascorbic acid crystals
  16. Broadband transfer of binary images via optically long wire media
  17. Counting and mapping of subwavelength nanoparticles from a single shot scattering pattern
  18. Controlling surface waves with temporal discontinuities of metasurfaces
  19. On the relation between electrical and electro-optical properties of tunnelling injection quantum dot lasers
  20. On-chip multivariant COVID 19 photonic sensor based on silicon nitride double-microring resonators
  21. Nano-infrared imaging of metal insulator transition in few-layer 1T-TaS2
  22. Electrical generation of surface phonon polaritons
  23. Dynamic beam control based on electrically switchable nanogratings from conducting polymers
  24. Tilting light’s polarization plane to spatially separate the ultrafast nonlinear response of chiral molecules
  25. Spin-dependent phenomena at chiral temporal interfaces
  26. Spin-controlled photonics via temporal anisotropy
  27. Coherent control of symmetry breaking in transverse-field Ising chains using few-cycle pulses
  28. Field enhancement of epsilon-near-zero modes in realistic ultrathin absorbing films
  29. Controlled compression, amplification and frequency up-conversion of optical pulses by media with time-dependent refractive index
  30. Tailored thermal emission in bulk calcite through optic axis reorientation
  31. Tip-enhanced photoluminescence of monolayer MoS2 increased and spectrally shifted by injection of electrons
  32. Quantum-enhanced interferometer using Kerr squeezing
  33. Nonlocal electro-optic metasurfaces for free-space light modulation
  34. Dispersion braiding and band knots in plasmonic arrays with broken symmetries
  35. Dual-mode hyperbolicity, supercanalization, and leakage in self-complementary metasurfaces
  36. Monocular depth sensing using metalens
  37. Multimode hybrid gold-silicon nanoantennas for tailored nanoscale optical confinement
  38. Replicating physical motion with Minkowskian isorefractive spacetime crystals
  39. Reconfigurable nonlinear optical element using tunable couplers and inverse-designed structure
Jetzt anmelden und 20% sparen
Institutioneller Zugang
Wie funktioniert Authentifizierung?
Haben Sie eine Idee, wie wir unsere Website verbessern können?
Bitte schreiben Sie uns.
© 2025 De Gruyter Brill
Heruntergeladen am 16.9.2025 von https://www.degruyterbrill.com/document/doi/10.1515/nanoph-2022-0809/html
Button zum nach oben scrollen