Space-Time Modulated Loaded-Wire Metagratings for Magnetless Nonreciprocity and Near-Complete Frequency ConversionPreprint: APS/123-QED

In recent years substantial efforts have been devoted towards the development of magnet-less non-reciprocal devices in electromagnetics as well as in acoustics Tanaka1965; Ayasli1989; Gallo2001; Kodera2011; Sounas2013; Wang2012; Popa2012; Manipatruni2009; Mazor2019_2; Fleury2014. In particular, space-time modulated structures have been explored theoretically and experimentally for this purpose Qin2014; Reiskarimian2016; Yu2009; Lira2012; Qin2016; Sounas2013_2; Estep2014. These designs are based on the creation of a synthetic sense of motion Mazor2019 that fundamentally enables the breach of time-reversal symmetry with respect to the guiding-wave sub-system. While the major part of these efforts has been dedicated to breaking reciprocity in guiding wave structures, some important focus has been also aimed at the violation of reciprocity in scattering and radiation Hadad2015; Hadad2016; Taravati2017; Taravati2017_2. For example in Hadad2016 a space-time modulated travelling wave antenna with different radiation patterns in transmit and receive mode was studied theoretically and experimentally. In that antenna, a few voltage varying capacitors were used to establish the required modulation. The underlying idea of this design is to amplify the effect of the spatiotemporal modulation by taking advantage of asymmetric interband transitions between different space-time harmonics of a mode inside and outside the light-cone. Thus, in a sense, the natural sharp filtering property of the light cone has been used to enhance the nonreciprocal effect in this low-Q leaky wave antenna system. In another work Hadad2015, a metasurface that consists of a space-time modulated impedance operator has been explored theoretically. In that proposal non-reciprocal Wood’s anomaly has been achieved through a careful design of the space-time modulation of the surface characteristics. However, unlike the antenna in Ref. Hadad2016 which required a small number of modulation actuators, for the emulation of the effective surface impedance in Hadad2015, an unrealistically, extremely complex configuration of voltage varying capacitors and dense wiring system is required to achieve the many different modulation regions.

A possible venue to overcome this issue may be based on using metagratings which are metasurfaces with electrically large unit cells. Metagratings Chalabi2017; Radi2017; Epstein2017; Rabinovich2018; Popov2018; Popov2019; Popov2019b; Epstein2018; Diaz-Rubio2017; Taravati2019, a close relative of frequency selective surfaces and diffraction gratings, have been explored recently as another means to overcome the challenge of perfect anomalous reflection and refraction by an electrically thin metasurface Yu2011; Monticone2013; Li2014; Pfeiffer2013; Selvanayagam2013; Selvanayagam2013_2; Pfeiffer2013_2; Sievenpiper2002; Hasman2003; Lin2014; Wong2014; Kim2014; Epstein2016; Epstein2014; Epstein2016_2; Dorrah2018; Sounas_Estakhri2016. As opposed to the traditional metasurface approach in which electrically deep subwavelength unit cells are involved, the unit cell of the metagrating is ${\cal O}(\lambda/2)$ or more and thus fundamentally gives rise to propagating diffraction orders. Then, a proper design of the metagrating’s unit-cell enables a significant control over the balance of the various propagating diffraction orders, and in particular makes it possible to nullify the zero order harmonic which corresponds to regular reflection and refraction. Moreover, recently, an interesting proposal for a dielectric slab metagrating with space-time modulated refractive index has been explored demonstrating a peculiar relations between the space-time harmonics in this system Taravati2019. In particular, non-reciprocal scattering has been theoretically demonstrated. Nonetheless, a realistic implementation of space-time modulated refractive index is challenging and requires either nonlinear media with significant pump power or the use of effective medium, e.g., using a dense array of unit-cells that are loaded by voltage-varying capacitors. Thus, unfortunately raising again the challenge of realizability.

Here, inspired by the recent progress in metagratings Chalabi2017; Radi2017; Epstein2017; Rabinovich2018; Popov2018; Popov2019; Popov2019b; Epstein2018; Diaz-Rubio2017; Taravati2019, we propose a solution to this problem and design efficient non-reciprocal metasurfaces by using only a minimal number of three distinct modulation regions and small modulation parameters. Specifically, we develop a rigourous theory for metagratings with time-modulated resonant elements. Using this approach we explore analytically a two-dimensional lattice of resonant, capacitively loaded, wires that are subject to space-time modulation. To that end we first develop a generalization of polarizability theory for time-modulated particles, and later exploit it together with the discrete-dipole approximation Draine1994; Tretyakov_Book; Hadad2013; Hadad2010; Mazor2012; Mazor2015 and the proper Green’s function to explore analytically the scattering from an infinite space-time modulated metagratings. We particularly show that non-reciprocal anomalous reflection can be achieved with nearly perfect coupling efficiency between the incident wave and the desired scattered wave. As opposed to anomalous reflection obtained by a stationary metagratings, here, the anomalous reflection process involves also an efficient frequency conversion process.

Consider a periodically loaded perfect electrically conducting (PEC) wire that is co-aligned with the $z$-axis as illustrated in Fig. 1(a). The wire radius is $r_{0}$, the loading impedance is $Z_{L}$, and the loading periodicity is small on the wavelength, i.e., $\Delta\ll\lambda$. Now let ${\bf{E}}(\omega)=\hat{z}E(\omega)$ be the electric field at the wire location but in the absence of the wire itself. Then, the induced current on the wire is given by Tretyakov_Book

with the effective susceptibilityComm $\alpha$,

where

is the susceptibility of an unloaded PEC wire. In Eq. (3), $\eta=120\pi\Omega$ and $k=\omega/c$ are the free space impedance and wavenumber, respectively. Also, $\omega$ is the radial frequency, $c$ is the speed of light in vacuum, and $H_{0}^{(2)}$ is zero order Hankel function of the second kind. To simplify subsequent notations, in the following we use

Then, for capacitive periodic loading we can write

where $C$ is the per-unit-length loading capacitance.

Up to this point the formulation has been strictly carried out in the frequency domain. However, in the following we shall introduce temporal modulation of the loading capacitor. We assume

where $\omega_{m}$ is the modulation frequency, $m$ is the modulation index and $\phi$ is the modulation phase shift between the wires in the meta-grating structure. To model the scattering process in this case it is required to repeat the previous derivation, but in the time-domain. To that end we use the convolution property of the Fourier transform on Eq. (1) with Eq. (4). We have

where ${\cal E}(t),\gamma(t)$, and ${\cal I}(t)$ are the time domain counterparts of $E(\omega),\gamma(\omega)$, and $I(\omega)$, respectively, and $*$ denotes convolution as defined in Eq. (7). By explicitly expanding Eq. (7) using Eq. (5), we get

Here ${\cal F}^{-1}$ stand for the inverse Fourier transform. In the last equality we used the identity ${\cal F}^{-1}\{I(\omega)/j\omega\}=\int_{-\infty}^{t}{\cal I}(\tau)d\tau-0.5I(\omega)|_{\omega=0}$, together with the fact $I(\omega)|_{\omega=0}=0$ (implying no DC currents in the problem).

Assuming that the exciting field is harmonic at frequency $\omega$, then, due to modulation at $\omega_{m}$, the induced current and the electric field in the problem may be expressed by

where ${\cal X}$ stands for either ${\cal I}$ or ${\cal E}$, and ${\cal X}_{n}$ denotes the complex coefficient of n’th harmonic at $\omega_{n}=\omega+n\omega_{m}$. Next, in order to obtain a difference equation that relates between the different harmonics of the exciting electric field and the induced current, we apply Eq. (10) in Eq. (8). Beginning with the last term in Eq. (8), we find after a simple derivation,

Next, we evaluate the first term in Eq. (8). To that end we use the convolution identity

with the Fourier transform of Eq. (10) (here, ${\cal X}={\cal I}$)

Therefore leading to

Finally, by combining Eqs. (14) and (II) in Eq. (8) and with Eq. (10), we get

Using Eq. (15) it is possible to define a “generalized” effective susceptibility $\alpha$ for time-modulated wires. In this case, $\alpha$ becomes a matrix that relates between the electric field and the induced currents in different harmonics. For example, if we choose to take into account only three harmonics, $n=0,\pm 1$, that is at $\omega_{0}=\omega,\omega_{\pm 1}=\omega\pm\omega_{m}$ we get

where $\tilde{{\bf{E}}}=[{\cal E}_{-1},{\cal E}_{0},{\cal E}_{1}]^{T}$ and $\tilde{{\bf{I}}}=[{\cal I}_{-1},{\cal I}_{0},{\cal I}_{1}]^{T}$, and

Consider a linear array of capacitively loaded time-modulated wires, located at $y=h$, above an infinite PEC plate that is placed on the $y=0$ plane. The inter-wire spacing is $a$ and the wires are loaded with capacitors $C(t)$, as in Eq. (6), with

where $l\in\mathbb{Z}$ is the wire index and $N$ is the spatial periodicity. For simplicity, in the following we shall assume that $N=3$, which is the minimal spatial periodicity that is required to achieve the effect of a synthetic motion, and thus, non-reciprocity. An illustration of the structure in this case is shown in Fig. 1(b).

In this section we use the theory developed above in order to demonstrate one possible functionality of the time-modulated metagratings surface. We consider a space-time modulated metagratings structure with the following parameters: reference frequency $f_{r}$, reference wavelength $\lambda_{r}=c/f_{r}$, inter-wire distance $a=0.6\lambda_{r}$, wires distance from the PEC plate $h$, wire loading capacitance $C_{0}=1$pF, loading periodicity $\Delta=0.1\lambda_{r}$, and wire radius $r_{0}=0.5$mm. We introduce space-time modulation of the wires in the lattice such that $\omega_{m}=0.1\omega_{r}$ ($\omega_{r}=2\pi f_{r}$) with modulation index $m=0.1$, and super-cell size $N=3$. A stationary metagratings structure, namely, in the absence of the space-time modulation, has been shown to be able to completely nullify the specular reflection Chalabi2017; Radi2017; Epstein2017; Rabinovich2018; Popov2018; Popov2019; Popov2019b; Epstein2018; Diaz-Rubio2017. Instead, the energy is directed to one of the grating lobes of the structure. This functionality has been recently termed - anomalous reflection. Nonetheless, up to now this effect has been demonstrated only in reciprocal structures, with the exception in Taravati2019 which explored a metagrating dielectric slab with a continuous space-time modulation of the refractive index. Due to reciprocity, if an impinging wave at $\theta^{i}$ is fully reflected to $\theta^{r}$, then a wave arriving at $-\theta^{r}$ will be reflected with the same efficiency to $-\theta^{i}$, where $\theta^{i}$ and $\theta^{r}$ are measured with respect to the $y$ axis, and considered positive when corresponding to waves with positive $k_{x}$ component of the wave number. Here, we show that upon space-time modulation this effect can be turned to be strongly non-reciprocal, and moreover, we demonstrate that even with a small modulation index it is possible to achieve a near complete frequency conversion by the temporal modulation.

With the parameters listed above, in Fig. 2 we show the suppression of the specular reflection as a function of the incidence angle $\theta^{i}$ and the distance between the lattice and the PEC plate $h$.

The figure shows the specular reflection intensity, color-coded in logarithmic scale. For large range of parameters nearly full reflection takes place. However, at certain regions in the reflection intensity map inter-harmonic resonances take place leading to full suppression on the specular reflection. From the asymmetry in the figure, it is clear that this effect is non-reciprocal, leading to a very different power transmission for impinging waves that are coming from complementary direction, $\theta^{i}$ or $-\theta^{i}$. Let us consider a specific example, with $h=0.7775\lambda_{r}$ and incidence angle $\theta^{i}=-30^{\circ}$. With these parameters, using Eq. (III.4) it is clear that besides the specular reflection (fundamental harmonic $(m,n)=(0,0)$), also the higher order $(m,n)=(1,-1)$ and $(m,n)=(0,1)$ harmonics are propagating. This means, that the incident wave power that is not specularly reflected has to efficiently couple to one or two of these additional propagating diffraction orders. This is demonstrated in Fig. 3 below. In the figure, the power flux density along the $y$ direction, $S_{mn}$, is shown for the propagating harmonics. The specular reflection wave due to an incident wave at $\theta^{i}=-30^{\circ}$ is shown in blue. At the desired reference frequency, specular reflection is completely suppressed. The incident wave energy experiences practically complete conversion to the $(m,n)=(1,-1)$ harmonic that propagates at frequency $\omega_{-1}=\omega_{r}-\omega_{m}$ which is reflected towards $\theta_{(1,-1)}\approx 42^{\circ}$. As opposed to that, the $(m,n)=(0,1)$ harmonic which is also propagating is very weakly excited. On the contrary, in the absence of inter-harmonic resonances as shown in Fig. 2, the wave that impinges at the complementary direction, i.e., at $\theta^{i}=30^{\circ}$ experiences nearly complete specular reflection, and therefore, practically, no additional propagating harmonics will be excited.

A few words on practical realization of the proposed design. In the example considered above only $N=3$ modulation regions are required per unit cell. The distance between the wires is large on the wavelength thus enabling a convenient wiring system. Furthermore, we note that all the capacitors that are distributed on each wire can be modulated simultaneously by modulation voltage that is applied at the two ends of each wire. Moreover, an additional simplification of the practical implementation may be achieved if each family of wires that are subject to the same modulation (namely, all the wires of the same color in Fig. 1) will be placed at a different height with respect to the ground plane. In that case, the modulation scheme may be designed using even a smaller number of modulation drivers. Note though that the analysis carried out in this paper should be slightly adjusted for this particular case.

To conclude, in this paper we have developed a theoretical model based on the discrete-dipole-approximation and polarizability theory for the scattering by a two-dimensional lattice of space-time modulated resonant capacitively loaded wires. Using this simple methodology we demonstrate the design of a significantly non-reciprocal anomalous reflection using metagrating structure with electrically large unit cells that contain a single wire each. Thus, significantly simplifying the practical requirements for the modulation system and hence opening a realistic possibility for actual fabrication of such devices. The non-reciprocal scattering process described here is also associated with efficient frequency conversion of the anomalous refracted beam. While the model explored in this paper involves capacitively loaded wires, the same approach may be also augmented to the optical frequency regime using fast modulation techniques Phare2015; Liu2011. And with possible implications in radio-frequency devices as well as for more efficient photovoltaic processes Green2012; Zhu2014. Moreover, we note that the model developed in this paper may be readily augmented to solve the excitation response due to a localized source by taking an approach akin to Hadad2011; Hadad2013_2.