Nonradiating and radiating modes excited by quantum emitters in open epsilon-near-zero cavities

See allHide authors and affiliations

Science Advances  21 Oct 2016:
Vol. 2, no. 10, e1600987
DOI: 10.1126/sciadv.1600987


Controlling the emission and interaction properties of quantum emitters (QEs) embedded within an optical cavity is a key technique in engineering light-matter interactions at the nanoscale, as well as in the development of quantum information processing. State-of-the-art optical cavities are based on high quality factor photonic crystals and dielectric resonators. However, wealthier responses might be attainable with cavities carved in more exotic materials. We theoretically investigate the emission and interaction properties of QEs embedded in open epsilon-near-zero (ENZ) cavities. Using analytical methods and numerical simulations, we demonstrate that open ENZ cavities present the unique property of supporting nonradiating modes independently of the geometry of the external boundary of the cavity (shape, size, topology, etc.). Moreover, the possibility of switching between radiating and nonradiating modes enables a dynamic control of the emission by, and the interaction between, QEs. These phenomena provide unprecedented degrees of freedom in controlling and trapping fields within optical cavities, as well as in the design of cavity opto- and acoustomechanical systems.

  • Epsilon-near-zero materials
  • ENZ
  • flexible cavity
  • open cavity
  • quantum emitter
  • radiating mode
  • nonradiating mode


Cavity quantum electrodynamics (QED) is the field of research that investigates the interaction between quantum emitters (QEs), such as atoms and quantum dots, and a resonant cavity (1). This interaction is fundamentally interesting, and it could be the basis for quantum information processing (24). We emphasize that not only coherent but also dissipative interactions can be exploited with this purpose (57). QE-cavity interactions are also relevant for single-photon nonlinearities (8), lasing (911), quantum many-body systems (12, 13), and cavity optomechanics (14). At infrared and optical frequencies, dielectric microcavities (15), photonic crystals (1620), and Anderson-localized modes (21) are commonly used owing to their low losses and associated high-quality factors. Despite losses, plasmonic systems (2225) are also attractive because they provide subwavelength confinements, with cavity sizes well below the diffraction limit.

Aside from this spectrum of conventional cavities, more sophisticated responses in the emission and interaction properties of QEs could be obtained with cavities carved in more exotic materials. For instance, zero-index metamaterials [for example, epsilon-near-zero (ENZ) (26) or epsilon-mu-near-zero media (2729)] exhibit a decoupling between spatial and temporal field variations (28, 30), which enables numerous wave phenomena, including tunneling (26, 31) and geometry-invariant eigenfrequencies (32). In terms of tailoring the emission properties of QEs, the phase uniformity in zero-index metamaterials has been exploited to enhance the directivity (33, 34) and Purcell factor (3537) of single emitters, as well as to construct collective interference effects among multiple emitters (29, 3840). The wealth in wave phenomena related to metamaterials with near-zero parameters, as well as their potentiality in enhancing emission properties, has motivated us to investigate the emission properties of QEs embedded in open ENZ cavities, with a view toward their future application in cavity QED. Note that far from being a theoretical curiosity, there are several experimental demonstrations of zero-index metamaterials based on naturally available materials (41, 42), dispersion engineering in waveguides (31, 43, 44), photonic crystals (45, 46), and artificial electromagnetic materials (47, 48).

We demonstrate that the main signature of a QE embedded in an open ENZ cavity is the excitation of a nonradiating mode independently of the geometry of the external boundary of the cavity. Nonradiating modes have been investigated for a long time because of their connection to classical problems, such as models for stable atoms and elementary particles [see, for example, Devaney and Wolf (49, 50) and Marengo and Ziolkowski (49, 50)]. Furthermore, the extreme light confinement facilitated by nonradiating modes may also have practical applications in nonlinear optics, sensing, and heating (51, 52); in the storage of “bits” of quantized energy light (53); and in managing the reactive power surrounding an emitter (54). Recently, the excitation of nonradiating modes in spherical plasmonic cavities has been investigated (5153, 55). Here, we demonstrate theoretically that these apparently exotic nonradiating modes can be excited in open cavities with arbitrarily shaped boundaries and that they provide unique opportunities in controlling the emission properties of QEs.


Nonradiating modes in open ENZ cavities

We start by considering the emission properties of a QE located at the center of a spherical vacuum bubble of radius r0, which is itself immersed in an unbounded ENZ environment (see Fig. 1A). The insulating bubble is required to avoid the singularity that arises from the direct contact of a source with an absorbing medium (56, 57). The QE is modeled as a point-like two-level system with dipole moment p = p, intrinsic nonradiative decay rate ΓNR, and transition frequency ω0 (58), which is assumed to be centered at the ENZ frequency of the background medium (ω0 ≈ ωp). Moreover, we set ω0 ≈ ωp ≈ 2π × 29.08 × 1012 rad/s (that is, λ0 ≈ λp ≈ 10.32 μm), which corresponds to the plasma frequency of silicon carbide (SiC) (41), to facilitate future experimental demonstrations. Naturally, either SiC or any other practical implementations of ENZ media must be necessarily lossy and dispersive, and the effects predicted for the ENZ frequency will only be observed for sufficiently small losses and on a bandwidth of finite extent. Therefore, although the theory introduced here is derived in the ENZ limit for the sake of simplicity, numerical simulations, including the effect of losses and frequency dispersion, are included later in the paper.

Fig. 1 Spatially electrostatic fields in unbounded ENZ media.

(A) Geometry and sketch of a QE located at the center of a vacuum spherical bubble of radius r0, embedded in an unbounded ENZ medium, ε(ω0) ≈ 0. (B) Simulated electric and magnetic field magnitude distributions (r0 = 0.25 μm, λ0 = λp = 10.31 μm). The numerical simulations ratify the excitations of fields with a spatially electrostatic distribution and zero magnetic field in the ENZ region. (C) Analytically and numerically computed effective dipole moment enhancement factor at resonance (r00 = 0.715) as a function of losses (imaginary part of the permittivity of the ENZ region).

By solving the pertinent boundary value problem, we have shown (see Supplementary Note 1) that the spatial distribution of the classical fields excited by the aforementioned insulated dipole immersed in an unbounded ENZ medium is equal to those of an electrostatic dipole [although the fields oscillate in time at frequency ω0, exp(−iω0t) time convention], with effective dipole moment peff, that isEmbedded Image(1)andEmbedded Image(2)

Figure 1B depicts the electric and magnetic field magnitude distributions obtained with a full-wave numerical solver (see Materials and Methods), for a QE embedded in a spherical vacuum bubble of radius r0 = 0.25 μm. The figure ratifies that the magnetic field in the ENZ region is zero, and it is trapped within the vacuum bubble. Therefore, we find that insulated QEs embedded in ENZ media can be effectively treated as sources of spatially electrostatic fields with a time-harmonic variation. The strength of the effective electrostatic dipole moment peff is determined by the properties of the vacuum bubble (see Supplementary Note 1), and it is given byEmbedded Image(3)with k0 = ω0/c and Ĵ1(x) ≡ √ (πx/2) J3/2(x), where Jn(x) is the cylindrical Bessel function of the first kind and order n (59, 60). For example, the effective dipole moment is three times the original dipole moment for deeply subwavelength bubbles, that is, k0r0 ≪ 1→peff ≈ 3p, in consistence with the quasi-static solution to the problem (see Supplementary Note 2). Furthermore, it is apparent from Eq. 3 that the cavity conformed by the vacuum bubble becomes resonant at the zeros of Ĵ1(k0r0). Moreover, because the spatially electrostatic nature of the fields prevents radiation losses, the magnitude of peff is only limited in practice by dissipation losses. By taking the appropriate limiting procedure (see Supplementary Note 3), we find that at resonance [that is, at Ĵ1(k0r0) = 0], peff can be asymptotically written asEmbedded Image(4)where Ĵ1′(x) ≡ ∂x{Ĵ1(x)} and ε″ is the imaginary part of the relative permittivity of the ENZ host medium. This simple analytical rule is validated in Fig. 1C, which depicts peff as computed with a full-wave numerical solver (see Materials and Methods). Note that even with an amount of loss comparable to that of naturally available ENZ materials [ε″ ≈ 0.1 (41)], the magnitude of the effective dipole is 46 times larger than that of the original dipole. Equivalently, the electric field intensity in the vicinity of the bubble is enhanced by a factor of (46)2 = 2116, which suggests applications in enhancing the coupling between the QE and its environment, as well as in triggering nonlinear phenomena.

Next, as illustrated in Fig. 2A (see also fig. S1), the QE is placed within an ENZ environment of finite size, defining an open cavity. In the most general case, the cavity can be of any geometry, including the presence of other dielectric bodies within it. In this generalized scenario, the fields excited in the ENZ host will be those of the unbounded case, plus the fields “scattered” at the interface of the cavity with the unbounded external space and the dielectric bodies. Because the sources of the problem initially excite spatially electrostatic fields within the ENZ region, a valid solution to the scattering problem is the spatial distribution determined by the solution of the spatially electrostatic problem. Owing to the uniqueness of the solution, this is the spatial distribution of the fields excited in the time-harmonic case. In this manner, the electric field in the ENZ region can be written as the gradient of a scalar potential E(r) = −∇V(r), and consequently, it corresponds to the solution of the Laplace equation ∇2V(r) = 0 subject to the boundary conditions imposed by the continuity of both V and ε ∂V/∂n.

Fig. 2 Nonradiating modes in open ENZ cavities.

(A) Sketch of a bubble-insulated QE (shown as red arrow) embedded within an open ENZ cavity of arbitrary shape (shown as gray background) with several vacuum bubbles (shown in green). Simulation results for the (B) electric and (C) magnetic field magnitude distributions inside and outside the cavity, when the QE is at the center of its vacuum bubble. The electric field is trapped within the open cavity, whereas the magnetic field is confined to the vacuum bubble.

Crucially, and as it is noted by Landau and Lifshitz (61), the Laplace equation is independent of the background permittivity, which only appears in the solution through the boundary conditions. Only the ratio between the permittivities at each side of the interface matters. That is, the field distribution of a body of permittivity ε2 immersed in a background of permittivity ε1 is identical to that of a body of permittivity ε21 immersed in vacuum (61). Therefore, when the background medium is ENZ, ε1 ≈ 0, all material bodies behave as effectively perfect electric conductors, ε21→∞, for spatially electrostatic fields. Consequently, the fields excited by the QE will be unable to either escape the cavity or penetrate any dielectric body within this ENZ region. This effect is illustrated in Fig. 2B, which depicts the electric and magnetic field magnitude distributions in the XZ-plane cut. The magnetic field is again trapped within the vacuum bubble containing the QE at its center. By contrast, the electric field penetrates within the ENZ cavity in the form of a time-varying electrostatic field, but it is nonetheless unable to either escape the cavity or enter the other dielectric bodies.

This geometry-invariant confinement can also be understood by noting that the bound charges excited at the interface of ENZ and vacuum have distributions identical to those that would be excited at the interface of vacuum and a perfect electrical conductor (PEC) of analogous geometry. The former are given by σb = −∇∙P = ε0(ε − 1)En1 = −ε0En1, whereas the latter are given by σf = ∇∙D = −ε0En1, where En1 is the electric field normal to the interface in the first medium. Because identical charge distributions give rise to the same electric field, it is clear that, in terms of finding the solution to the scattered field, all boundaries in contact to the ENZ medium effectively behave as PEC boundaries. Thus, QEs embedded in open ENZ cavities excite nonradiating modes that are confined within the extent of the cavity independently of its geometry and even when the cavity itself is open to an unbounded vacuum space. We emphasize that although nonradiating modes have already been predicted in open spherical plasmonic cavities (5153, 55), these modes exist in open ENZ cavities independently of the geometry of its external boundary and that this effect is empowered by the spatially electrostatic nature of the fields, as demonstrated here.

Switching between nonradiating and radiating modes

The above analysis is valid as long as the fields excited by the QE-vacuum bubble system are dominated by the electric dipole mode in the bubble. This is exactly the case when the QE is at the center of a spherical bubble, and it is a very accurate estimation for subwavelength bubbles. At the same time, bubbles with larger sizes can efficiently excite different, possibly radiating, modes. Far from being a limitation, this opens up the possibility of switching between nonradiating and radiating modes, and hence activating/deactivating the interaction of QEs with a system external to the cavity. Take for example, as schematically depicted in Fig. 3A (see also fig. S2), an open ENZ cavity with arbitrarily shaped external boundary containing a spherical vacuum bubble, whose radius has been tuned to be resonant with a radiating mode (for example, a magnetic dipolar mode); because of the symmetry of the fields, only the electric dipolar mode is excited when the QE is at the center of the sphere, resulting in a nonradiating cavity mode. However, this symmetry is broken as the QE is shifted from the center of the bubble, and the QE excites the resonant radiating mode. In practice, the position of the QE can be controlled with different mechanisms, for example, sound waves, microelectromechanical systems (MEMSs) (62), and optomechanical techniques. One of these options is schematically illustrated in Fig. 3A, in which the QE is assumed to be attached to a membrane in the cavity. Thus, if the cavity were excited by an external optical or acoustic wave, then the membrane would vibrate and the position of the QE and, hence, its emission properties would oscillate in synchrony with the vibrational mode of the cavity.

Fig. 3 Switching between nonradiating and radiating modes.

(A) Geometry and sketch of an arbitrarily shaped open ENZ cavity, with a vacuum spherical bubble (r0 = 5.27 μm), containing a QE attached to a membrane, so that its position may be displaced along the x axis because of an external stimulus, for example, the vibrational modes excited by an external optical/acoustic wave. (B) Simulated (at λ0 = λp = 10.31 μm) electric field magnitude distributions for displacements: Δx = 0 (nonradiating mode) and Δx = 3.5 μm (radiating mode) in the XZ-plane cut. Excitation rate Γexc (normalized to the free-space excitation rate) (C) and quantum yield ηrad (D), as a function of emitter displacement, for different amounts of loss in the ENZ medium (ε″) and assuming an intrinsic quantum yield ηint = 0.5.

This effect can be appreciated in Fig. 3B, which depicts the simulated electric field distribution in the XZ-plane cut for the QE emitter at two different emitter positions (see Materials and Methods). As anticipated, a nonradiating mode is excited for Δx = 0, and the field is confined within the cavity. By contrast, the field is strongly radiated outside the cavity for Δx = 3.5 μm. In our numerical simulation, the cavity has been considered with an arbitrary (not particularly designed) external boundary, although the radius of the internal spherical vacuum bubble has been numerically optimized to r0 = 5.27 μm to trigger the magnetic dipolar resonance when the QE is displaced away from the bubble’s center (see Supplementary Note 4 and fig. S3).

We make use of the nanoantenna formalism developed by Novotny and van Hulst (63) and Bharadwaj et al. (64) to quantitatively assess the emission properties of a QE embedded in this specific cavity. In particular, we calculate the quantum yield ηrad (or radiation efficiency, that is, the ratio of radiative to total decay rates, where the latter includes both the dissipation decay in the cavity and the intrinsic nonradiative decay rate), as well as the normalized excitation rate Γexcexcfree (that is, the rate of excitation via spontaneous emission of a receiver located outside the cavity, at the position schematically depicted in Fig. 3A, and normalized to the free-space excitation rate) with Γexc = Γrad Drad, where Embedded Image dS and Embedded Image are the radiative decay rate and radiation directivity, respectively, and S = 1/2 Re{E × H*} (63, 64). Both quantities are depicted in Fig. 3 (C and D) as a function of the emitter position displacement, Δx, and for different amounts of loss of the ENZ host medium ε″. For illustrative purposes, we assume that the intrinsic quantum yield is 0.5. It is evident from Fig. 3 (C and D) that both the excitation rate and the radiation efficiency are suppressed when the dipole is at the center of the vacuum bubble, consistent with the excitation of nonradiating modes. Moreover, both figures of merit increase as the QE is shifted from the center, and they reach a maximum at the specific displacement Δx = 3.5 μm. This optimal value corresponds to the maxima of the magnetic dipolar coefficient as computed from the addition theorem (see Supplementary Note 5 and fig. S4) (65). Naturally, the excitation rate and radiation efficiency are limited by the losses of the ENZ medium. However, even with losses ε″ = 0.1, the excitation rate is enhanced by a factor of 5 with respect to that of free space (see Fig. 3C). This result indicates that the effect of switching between radiating and nonradiating modes could potentially be experimentally validated at λ ~ 10 μm using SiC (41) or at telecommunication wavelengths (λ ~ 1550 nm) using transparent conducting oxides (42, 66). Better performances and even higher frequencies of operation might be obtained with synthetic ENZ media [such as waveguides at cutoff frequencies (31, 43, 44), photonic crystals (45, 46), and/or metamaterials (47, 48)]. Furthermore, we emphasize that a similar effect takes place for near-field interactions because, as shown in Fig. 3B, the near field outside the cavity is also suppressed when a nonradiating mode is excited.

Coupling/decoupling of QEs

This effect also empowers a dynamic control of the dipole-dipole interactions between different emitters and their associated collective effects (67, 68). As depicted in Fig. 4A, one could embed two emitters, for example, p1 and p2, in different spherical bubbles placed within an open ENZ cavity. When both emitters are at the center of their bubbles, they excite nonradiating modes and are effectively decoupled. However, they become resonantly coupled as their position separates from such a symmetric position. In this case, we can assume that the position of the QEs may be controlled by means of, for example, MEMSs placed outside the cavity (62). To this end, the cavity may again be assumed to be formed by an ENZ host with a noncanonical geometry (see Fig. 4A and fig. S5), whereas spherical bubbles, separated by a distance d = 6 μm, can be assumed to be made of silicon (Si), characterized by relative permittivity εSi = 11.7, and their radius r0 = 1.505 μm has been tuned to trigger the magnetic dipolar resonance (see Supplementary Note 4 and fig. S4). The cavity is pierced by a Si rod of cross-section 0.25 μm × 0.25 μm to illustrate that the proposed cavities are robust against the geometrical modifications that could be required to implement position control of mechanisms (for example, MEMSs).

Fig. 4 Coupling/decoupling QEs in open ENZ cavities.

(A) Geometry and sketch of an open ENZ cavity of arbitrary shape with two Si spherical bubbles (r0 = 1.505 μm, εSi = 11.7), separated by a distance of d = 6 μm, and containing QEs. The cavity is assumed to be pierced by a Si rod of cross-section 0.25 μm × 0.25 μm, whose position may be externally controlled by MEMSs. (B) Simulated (at λ0 = λp = 10.31 μm) electric field magnitude distribution excited by the QE in the left bubble for displacements: Δx = 0 (decoupled) and Δx = 1 μm (coupled) in the XZ-plane cut. (C) Individual decay rate, Γ11, decay rate related to coupling, Γ21, and cooperative Lamb shift, Δω21, normalized to the free-space individual decay rate, Γ0, as a function of frequency in the coupled (Δx = 1 μm, first row) and decoupled (Δx = 0, second row) states. The ENZ host has been modeled with a dispersive Drude model ε(ω) = 1 − ωp2/ω(ω + iωc), where ωp = 2π × 29.08 × 1012 rad/s, and for three different amounts of loss, characterized by ωc = 0.001ωp (first column), ωc = 0.01ωp (second column), and ωc = 0.1ωp (third column).

The coupling/decoupling mechanism mediated by the nonradiating and radiating modes is illustrated in Fig. 4B, which represents the simulated electric field magnitude distributions in the XZ-plane cut excited by the QE in the left bubble for displacements Δx = 0 μm (decoupled) and Δx = 1 μm (optimal coupling case according to the addition theorem; see Supplementary Note 5). Note that in our simulation, the presence of the Si rod has no appreciable impact on the radiating and nonradiating nature of the fields. A quantitative estimation of the coupling/decoupling performance is gathered in Fig. 4C. To this end, we compute the individual decay rate, Embedded Image, n = 1, 2; the decay rate due to coupling, Embedded Image; and the cooperative Lamb shift, Embedded Image (69, 70). On the one hand, in the decoupled case (Δx = 0), the decay rate and coupling parameters are orders of magnitude smaller than the free-space decay rate, and the QEs are effectively decoupled. On the other hand, the cooperative processes excited in the coupled configuration, Δx = 1 μm, are defined by the geometry of the cavity and the insulating bubbles. In this case, the dynamics of the coupling between the emitters are equivalent to those of two coupled single-mode resonators, to the splitting of the individual bubble resonances into subradiant and superradiant states. This effect can be qualitatively understood by using a simple circuit model [see Supplementary Note 6, fig. S6, and Liberal et al. (71) for the circuit model analysis of the light scattering of two coupled nanoparticles].

These simulations also serve to illustrate that, although an ENZ response must be necessarily dispersive, the bandwidth provided by one particular implementation can be sufficient for certain applications. Specifically, it is clear from Fig. 4C that the coupling between both emitters in the Δx = 0 μm (no displacement) configuration is suppressed on a reasonable bandwidth. In addition, if the linewidth of the resonators is smaller than this bandwidth, then resonant coupling (Δx = 1 μm) can still be observed without a detrimental effect because of the ENZ dispersion. Again, we emphasize that the validity of our analysis and the proposed configuration will be limited to those scenarios in which the interaction within the QE and the environment is dominated by this frequency band, in which the environment exhibits a near-zero permittivity response. This will be the case, for example, for spontaneous emission in the weak coupling regime.


Our results demonstrate that bubble-insulated QEs embedded in open ENZ cavities present the unique signature of exciting nonradiating modes independently of the geometry of the external boundary of the cavity. We believe that this effect provides unprecedented degrees of freedom in controlling and trapping electromagnetic fields within an open optical cavity. From a practical point of view, the possibility of having an arbitrarily shaped external boundary offers several opportunities in engineering this class of open cavities, also suitable for foldable, moldable, and flexible photonic platforms.

In addition, our study reveals that it is possible to switch between nonradiating and radiating modes. This could provide new venues in controlling the emission properties of QEs, such as enhancing/suppressing the spontaneous emission exiting the cavity, as well as dynamically activating/deactivating the coupling between them. The fact that these effects can take place in cavities with arbitrarily shaped boundaries could be exploited to resonantly couple with other physical processes, such as sound waves, enabling the coupling with specifically designed cavity-induced vibrational modes.

The geometry of the external boundary could also be tailored to boost the emission of radiating modes and/or to facilitate the excitation and manipulation of QEs with electromagnetic waves operating at frequencies where the cavity could be transparent or resonant. We emphasize that the geometries in the examples studied in Figs. 2 to 4 were selected to illustrate the possibility of using an arbitrarily shaped boundary. Thus, these configurations represent specific examples that are by no means optimized. For example, better performances in terms of stronger coupling parameters are obtained with a simple cubic cavity (see fig. S7). The design and optimization efforts required to meet the specifications of one or another application are out of the scope of the present work. Our work has also been restricted to spherical bubbles for the sake of simplicity and to establish a clear comparison with the theory. However, a variety of enhanced dissipative and coherent coupling phenomena could be achieved (and dynamically activated/deactivated) by using (not necessarily identical) multifrequency or multimode resonant bubbles.


The commercially available full-wave electromagnetic simulator software COMSOL Multiphysics version 5.2 (72) was used to compute the dyadic Green’s functions and the field distributions displayed in all figures of the main text and the Supplementary Materials. Specifically, we carried out analysis in the frequency domain solver using a tetrahedral mesh. The QE was modeled as a point dipole source with dipole moment p. The Green’s functions and field distributions provided by the numerical solver were used to compute the following related quantities: effective electrostatic dipole moment, excitation rate, quantum efficiency, decay rate associated to coupling, and photonic Lamb shift, which were evaluated according to their definition introduced in the main text. The effective dipole moment peff depicted in Fig. 1C was numerically computed by evaluating the field at the position (0, 0, 1.05r0), where r0 is the radius of the vacuum bubble containing the QE, and normalizing it relative to that of an electrostatic dipole with the same dipole moment.


Supplementary material for this article is available at

Supplementary Note 1. Derivation of Eqs. 1 to 3.

Supplementary Note 2. Quasi-static solution to the problem.

Supplementary Note 3. Derivation of Eq. 4.

Supplementary Note 4. Magnetic dipole resonances.

Supplementary Note 5. QEs shifted from the origin of the coordinates.

Supplementary Note 6. Equivalent circuit model.

fig. S1. Sketch and dimensions of the system studied in Fig. 2.

fig. S2. Sketch and dimensions of the system studied in Fig. 3.

fig. S3. Optimal bubble radius for resonant magnetic dipole excitation.

fig. S4. Electric and magnetic dipole excitations as a function of emitter displacement.

fig. S5. Sketch and dimensions of the system studied in Fig. 4.

fig. S6. Equivalent circuit model.

fig. S7. Coupling parameters in a cubic cavity.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Funding: This work was supported in part by the U.S. Air Force Office of Scientific Research Multidisciplinary University Research Initiatives on Quantum Metaphotonics & Metamaterials (award no. FA9550-12-1-0488) and the National Security Science and Engineering Faculty Fellowship award under the grant from the Office of Naval Research (award no. N00014-16-1-2029). Author contributions: I.L. and N.E. conceived the idea. I.L. carried out the analytical modeling and numerical simulations. N.E. supervised the project. Both authors contributed to the design of the problems, interpretation and analysis of the presented results, discussion and understanding of findings, and writing of the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from N.E. (engheta{at}

Stay Connected to Science Advances

Navigate This Article