Getting full control over light radiation requires materials with a vast range of characteristics. However, nature only provides us with a limited set of materials that can be utilized in order to manipulate electromagnetic waves. Artificial materials fill this gap by providing a rich platform to realize effective material properties not readily available in nature. Metasurfaces, as an important class of artificial materials, are two-dimensional arrays of artificially engineered subwavelength scatterers that provide a versatile platform to engineer the flow of electromagnetic waves, with promising applications from microwaves to optics [
1,
2,
3,
4,
5,
6]. Space-gradient metasurfaces are a special subclass of these electrically thin layers, tailored to exhibit spatially varying optical responses, thus providing a powerful tool to engineer the direction and wavefront of an incident wave in reflection and/or transmission. These surfaces have found many applications, including bending the flow of electromagnetic waves, cloaking, and sensing [
7,
8,
9,
10,
11,
12,
13].
Although these devices have been driving significant excitement in the optics community, most of the relevant proposals to date provide responses that are strictly bound by Lorentz reciprocity. Breaking reciprocity opens a whole new set of functionalities that are impossible with reciprocal structures [
14,
15]. Different approaches can be employed to break reciprocity in these structures [
16]. Magnetic gyrotropy (magnetized ferrite) or electric gyrotropy (magnetized plasma) have been among the first proposals to design nonreciprocal gradient metasurfaces (see, e.g., [
17]). However, nonreciprocal metasurfaces based on these materials usually lead to bulky, expensive, and difficult-to-integrate devices. As an alternative approach, nonlinear gradient metasurfaces have been proposed for nonreciprocal anomalous refraction or reflection (see, e.g., [
18]). However, these structures are intensity-dependent and may, as a result, distort the signal and not appeal several application areas.
Recently, the emerging field of time-varying metasurfaces has provided several interesting opportunities to expand the realm of exotic light-matter interactions with metasurfaces [
19,
20,
21], including nonreciprocal responses [
14,
15]. In particular, it was shown that time-modulated gradient metasurfaces open new degrees of freedom for engineered nonreciprocal responses. This concept has recently been employed to design nonreciprocal surfaces for optical isolation and circulation [
22,
23,
24]. Although surfaces with spatio-temporal modulation address many of the challenges related to the other approaches mentioned above, their functionalities are fundamentally limited by several issues: because of temporal modulation, we can generally expect frequency mixing with associated signal distortion [
22,
23,
24]. Furthermore, dynamic modulation is an active approach that requires external energy applied to impart the required modulation signals, and correspondingly bulky feeding networks in the case of electro-optical modulations.
In this letter, we introduce nonreciprocal metagratings based on media in synthetic motion, which overcome the abovementioned challenges. We first derive the minimum physical requirements to realize such a platform, and then propose a design to implement unitary-efficiency free-space circulation. In 1908, three years after Einstein presented the special theory of relativity, Minkowski formulated the electrodynamics of moving media [
25,
26], which, at the time, was a mysterious problem within the framework of classical physics. The first part of Minkowski’s theory deals with the invariant nature of Maxwell’s equations, stating that Maxwell’s equations for a stationary medium are also valid for a medium that moves uniformly at velocity
with respect to a stationary (fixed) reference frame. Then, constitutive relations for moving media are determined utilizing the special theory of relativity, assuming that the constitutive relations for the medium at rest are known [
25,
26].
As expected, these constitutive relations are nonreciprocal in nature, showing opportunities to be employed as building blocks of nonreciprocal surfaces. However, in order to get considerable nonreciprocal response, the medium should be moving with a speed comparable to the speed of light which is not practical. However, analogous nonreciprocal constitutive relations in stationary artificial composite media, i.e., metamaterials at rest, may be obtained [
27] by artificially engineering the inclusion bianisotropy to synthesize the electromagnetic behavior of a moving medium. The required coupling between electric and magnetic phenomena in this class of stationary engineered materials was termed bianisotropic “moving” coupling. Based on this concept, we can make stationary inclusions, which, as a function of the direction from which you look at them, respond electromagnetically as if the particles move towards you or against you (or in another direction) with a controllable speed. In the following, we present a platform based on bianisotropic moving particles for rerouting electromagnetic waves in a nonreciprocal fashion. Furthermore, we show how such nonreciprocal surfaces can be realized using stationary magnet-free inclusions.
Let us consider a two-dimensional array of bianisotropic moving particles located at a distance
from a ground plane (
Figure 1). We assume that the particles considered here can only provide electric moment in the
direction and magnetic moment along
. The structure is illuminated with a transverse magnetic (TM) plane wave with electric field component
. Here,
is the free-space wave impedance,
is the free-space wavenumber,
is the incidence angle, and
is the amplitude of the incident electric field. Because of this illumination, surface electric and magnetic currents will be excited on the structure, which sustain the scattering from the array:
where
;
is the Dirac’s delta function; and
,
are the periodicities in the
and
directions, respectively, and
is the area of a unit cell. Here,
and
are the surface electric and magnetic current densities on each unit cell. The radiated fields can be written in terms of Floquet modal expansion with electric field component
, where
. Here,
is the wavenumber in the
direction for the
mode (the first and second indices refer to the Floquet orders with respect to the
and
axes),
and
are the surface currents, which can be written in terms of effective electric
, magnetic
, and bianisotropic moving polarizabilities
as
and
, where
(
) is the electric (magnetic) field in the absence of the array. Assuming no Ohmic loss, passivity requires that the total radiated power carried away from the metasurface through all Floquet channels equals the extinction power
, where the radiated and extinction powers read
and
, respectively [
28]. Applying these conditions will impose the following constraints on the effective polarizabilities of the inclusions:
where
.
Now, we can tailor the polarizability response to define the desired nonreciprocal light-matter interactions. Depending on the defined functionality, bianisotropic particles with different polarizabilities, satisfying (3), will be required. As an example, here we design a metagrating that can operate as an ideal (i.e., with unitary efficiency) free-space circulator, realizing a surface that directs all the incident wave from direction A to B, all the incident wave from B to C, and all the incident wave from C back to A (
Figure 1). First, we choose the periodicity so as to align three diffraction orders with the directions of interest (i.e.,
, where
is the wavelength and
is the angle of Floquet channels from the normal to the surface). In order to reroute all normally incident waves from A towards B, first we need to nullify the reflection into the specular direction, which can be achieved if
. Next, we need to make sure that no scattering trickles towards diffraction order C, which can be guaranteed if
. It can be shown that the required polarizabilities for such functionality read.
Substituting these polarizabilities into the third condition of (3), we can derive the design equation . As a second design step, we need to make sure that, if the layer is illuminated from direction B, all the energy is transferred into the Floquet channel C and no scattering happens towards A. It can be shown that in this case we need (when the excitation is from port B, 00, 0(-1), and 0(-2) Floquet channels are aligned with ports C, A, and B, respectively). This corresponds to requiring that, when the metagrating is excited from port B, the bianisotropic moving particles do not get excited at all. In this case, the incident wave only sees the ground plane and, as a result, all the illuminated wave gets reflected into the Floquet channel C as initially planned. It is quite interesting to note that, although the bianisotropic moving particles are stationary, from an electromagnetic standpoint, they seem to move with the speed of the incident wave away from B in the same direction, so that the incident wave never gets a chance to interact with them. Finally, we need to make sure that, if we excite the metagrating from port C, the illuminated wave will be fully directed into port A and no scattering will happen towards B or back to C. Using (4), it can be shown that for this case and (when the excitation is from port C, 00, 01, and 02 Floquet channels are aligned with B, A, and C, respectively). Under these conditions, from an electromagnetic point of view, the layer seems to move in the direction antiparallel to the incident wave and, as it interacts with it, all the illuminated wave gets redirected into port A.
To understand how the metagrating behaves versus frequency, here we model the polarizabilities using the conventional Lorentz dispersion model which, near its resonance, can adequately describe the dipolar response of small inclusions, i.e.,
,
, and
, where
,
, and
are the amplitude coefficients and
,
, and
are the loss factors. Here, we assume that the resonance frequencies are the same for all polarizabilities (i.e.,
). As an example, here we design the metagrating to circulate the free space wave between three Floquet channels located 45 degrees apart (see
Figure 1). Amplitude coefficients
,
, and
can be calculated by substituting these polarizability expressions in (4) (note that, here, as an example, we assume
).
Figure 1b–d show the amplitude of the reflected wave into different Floquet channels under excitation from different Floquet channels. It can be seen that at resonance (i.e.,
), the surface is capable of circulating the space wave with unitary efficiency, which is impossible to achieve using gradient metasurfaces.
This concept is quite general, and different approaches can be taken to realize stationary particles with the required bianisotropic moving coupling. A simple approach can be based on ferrite inclusions, in order to get the required magnetic gyrotropy to break reciprocity, and combined with a proper set of electric dipoles in order to emulate the electromagnetic properties of a real moving inclusion [
27,
28,
29,
30]. However, as mentioned above, such inclusions require a magnetic bias, resulting in bulky and poorly integrable structures. Another recently proposed approach has been based on mimicking Faraday rotation using dc current bias [
31]. Following this approach, recently magnet-free nonreciprocal bianisotropic inclusions were designed and utilized as the building blocks of nonreciprocal metasurfaces [
32]. Here, inspired by these works, we design an artificial stationary inclusion with bianisotropic moving coupling. The inclusion is composed of a ring loaded with an ideal isolator, which allows rotational travelling wave only in one direction beneath the ring. An electric dipole in the form of a wire strip has been located underneath the ring that, in conjunction with the directional rotating magnetic dipole, provides the bianisotropic moving coupling. Based on this inclusion design, we simulate a nonreciprocal metagrating that is capable of ideally circulating space waves between three Floquet channels located 45 degrees apart. The designed metagrating is shown in
Figure 2a. The results in
Figure 2b show that the structure nicely circulates the space wave between three ports. As it can be seen from the simulated scattering parameters, despite the presence of realistic absorption loss due to the isolators used in the loops, the performance is close to ideal. As it was shown in the theoretical analysis, the nonreciprocal metagrating, in principle, can be designed to circulate space wave with unitary efficiency. Here, for the proof of concept, we considered the case in which the ports are located 45 degrees apart, however it should be noted that the involved ports can be located with almost 90-degree angular separation with respect to each other while providing unitary efficiency. More efficient inclusions can be utilized in order to provide the required nonreciprocal functionalities in a more efficient way.
To conclude, in this paper we have introduced a new conceptual platform for nonreciprocal wave manipulation with unitary efficiency. The proposed concept is based on nonreciprocal metagratings composed of stationary inclusions synthesizing bianisotropic moving coupling, while being stationary. As an example of nonreciprocal wave manipulation, we have shown that these metagratings can perform nonreciprocal circulation for space waves with unitary efficiency. There are a few important points that distinguish the proposed concept from other platforms presented for nonreciprocal wave manipulation. The proposed platform does not require a continuous profile of spatio-temporal modulation, but it simply consists of a periodic array of bianisotropic inclusions, where all the building blocks are identical. They are stationary, and simply biased by a dc current to break reciprocity. It was recently revealed how gradient metasurfaces suffer from fundamental bounds on conversion efficiency [
12,
33,
34], yet, in the proposed design there is no fundamental limit on coupling to one diffraction order to the other [
35], enabling a versatile platform for nonreciprocal wave manipulation with large efficiency. In contrast to structures based on spatio-temporal modulation, in the proposed design, no up- or down-conversion of frequency arises. The period in the array is comparable to the wavelength, avoiding severe challenges to translate these concepts to higher frequencies. We also envision similar concepts for nonreciprocal manipulation of acoustic waves based on similar principles.