## Abstract

By coupling controllable quantum systems into larger structures we introduce the concept of a quantum metamaterial. Conventional meta-materials represent one of the most important frontiers in optical design, with applications in diverse fields ranging from medicine to aerospace. Up until now however, metamaterials have themselves been classical structures and interact only with the classical properties of light. Here we describe a class of dynamic metamaterials, based on the quantum properties of coupled atom-cavity arrays, which are intrinsically lossless, reconfigurable, and operate fundamentally at the quantum level. We show how this new class of metamaterial could be used to create a reconfigurable quantum superlens possessing a negative index gradient for single photon imaging. With the inherent features of quantum superposition and entanglement of metamaterial properties, this new class of dynamic quantum metamaterial, opens a new vista for quantum science and technology.

©2011 Optical Society of America

## 1. Introduction

By offering material properties beyond that which occurs in nature, artificially engineered metamaterials are of intense interest. Typically fabricated with periodic features spaced closer than the operating wavelength, the system acts as a homogenous material. The earliest introductions were in the area of magnetic resonance imaging where conducting elements were used to produce *artificial magnetism* [1–5]. Negative index materials (NIMs) with simultaneous negative permittivity and permeability have also been engineered. A remarkable property of such negative index materials (NIMs) is the possibility of negative refraction [6]. Negative refraction arises due to the ability to design an interface between concave and convex surfaces in the material bandstructure. This interface leads to a novel dispersion relation not found in naturally occurring materials, namely all-angle negative refraction [7]. Simply put, conservation of energy forces a reversal of the wave-vector in one dimension, whilst preserving the wave-vector in the other dimension. Typically, the appropriate bandstructure engineering requires control over the electromagnetic properties on scales significantly smaller than the wavelength of the radiation, although any system that realizes the appropriate bandstructure engineering will exhibit negative refraction. NIMs with negative refraction have been demonstrated in the microwave regime in structures that consists of interlocking metal strips and conducting split-ring resonators [8–11]. By applying transformation optics, metamaterials become a rich platform for the control of electromagnetic waves. A striking consequence of such control has been the realization of an invisibility cloak [12]. Apart from electromagnetic metamaterials, acoustic [13,14] and seismic [15] metamaterials are also areas of intense research.

Here we introduce a new class of metamaterials distinguished from conventional metamaterials in that it operates in the quantum regime and is easily reconfigurable. The medium, a cavity array metamaterial (CAM), comprises of a network of coupled atom-optical cavities. Under rotating-wave approximation and in the tight-binding regime, the medium can be treated with the Jaynes-Cummings Hubbard (JCH) model [16–23]. Atom-optical cavities have already been shown to demonstrate such quantum effects as quantum collapse and revival [24], Rabi oscillations [25], photon blockade [26], and electromagnetically induced transparency [27]. The coupling of these cavities have theoretically been shown to exhibit quantum phase transitions [17, 23] and quantum chaos [28]. They have also been proposed to act as Q-swtiches [29], single photon switches [30], many-body [31] and semiconductor simulators [23]. Importantly, by electrical control of the atomic resonances, the metamaterial properties can be dynamically varied, a significant development compared to the more usual static implementations. The term quantum metamaterial was first introduced in the context of superconducting charge qubits inside a superconducting resonator [32]. Here we extend the concept of a quantum metamaterial by providing an alternative quantum platform based on cavity quantum-electrodynamics. This platform goes beyond classical electromagnetic properties and allows the exploitation of a quantum plurality of metamaterial properties.

As an illustration we present the design for a reconfigurable *quantum superlens* based on the CAM. We envisage a configuration depicted schematically in Fig. 1 where the JCH system is manipulated to produce a perfect image using single photons. The aim here is not build a better superlens, but to use the superlens as an casestudy of how CAMs can be used to exhibit metamaterial properties such as negative refraction. However as with photonic crystal (PhC) implementations, CAMs have the advantage of being ideally lossless and do not require the operating wavelength to be larger than the constitutent element spacing. Further, because the transition energy of each atom in the system can be individually controlled via a Stark shifting control voltage, CAMs have the distinct advantage of being highly tunable and reconfigurable: features not possible using conventional designs.

To form a perfect image requires the lossless convergence of the propagating and evanescent light components. Conventional lenses only focus the propagating fields, and so the resolution of the image is fundamentally limited to features greater than the optical wavelength. Subwavelength features are carried by the high spatial-frequency components encoded by the evanescent fields. The loss of the evanescent components lead to the *diffraction limit*. Near-field scanning optical microscopy overcomes this problem by scanning a probe in close proximity to the object, but this is often undesirable for applications such as optical lithography and sensing. It has been proposed that a lens built from NIMs can produce *perfect* far field imaging, exhibiting all-angle negative refraction (AANR) and evanescent wave enhancement (EWE) [33]. Because of their ability to overcome the diffraction limit, and the lack of optical axis and curved surfaces that AANR affords, such a lens is termed *superlens* [33, 34].

Negative indexing however is not a prerequisite for superlensing. A different class of metamaterial from that of NIM, is formed by photonic crystals (PhCs) which uses Bragg scattering. PhCs have also been shown to exhibit AANR and EWE [35–37], and have the advantage of low loss. The disadvantage of PhCs as a superlens is that not all the evanescent components can be uniformly amplified. In relation to the PhC as a metamaterial, it is of note that since the size and periodicity of the scattering elements in PhCs are on the order of the operating wavelength, the medium cannot be considered as homogeneous, which is a necessary condition to identify a meaningful permeability and permittivity.

The theoretical developments of superlenses have also been matched by experimental efforts. Superlensing has been demonstrated with microwaves in both NIMs [38] and PhCs [39, 40], as well as other platforms such as silver films in visible light [41]. More recently, near-field microscopy using a SiC-based superlens at mid-infrared frequency has successfully imaged features smaller than the illumination wavelength [42].

We will for the first time discuss cavity arrays as a metamaterial. The inter-cavity hopping mechanism of the CAM are markedly distinguished from NIM and Bragg scattering, which have been the hallmark of metamaterials until now, and therefore represent another class of metamaterial. In particular we will discuss the JCH system as a CAM and investigate it as a medium for a quantum-based reconfigurable superlensing device. The tunability of the system enables post fabrication dynamical control over the focal point of the lens, and by scanning over resonances, importantly allows uniform enhancement of evanescent modes, which is not possible with passive PhC-based devices. We further show a gradient negative index lens which minimizes reflection from the device.

## 2. Jaynes-Cummings-Hubbard Hamiltonian

We consider a uniform two dimensional (2D) periodic array of coupled optical cavities, embedded with single two-level atomic systems. Our treatment is implementation independent and could be realized in any coupled-cavity array system that realizes a two-dimensional square lattice topology. However it is useful to consider a concrete implementation. We specifically consider a coupled-cavity array realized in a two-dimensional photonic crystal membrane. In such a structure the thickness of the membrane effects confinement in the third dimension, whilst the photonic crystal structure defines the cavity array. Such a structure has been considered previously in the context of photonic quantum emulators [17, 20, 23], and one-dimensional coupled cavity arrays have been fabricated [43]. The cavity array is a series of quantum oscillators coupled through the overlapping of the photonic modes of adjacent cavities, such that it leads to a tight-binding Hubbard-like model. Each cavity couples to its atom via the Jaynes-Cummings (JC) interaction. In terms of the atomic (photonic) raising and lowering operators
${\sigma}_{r}^{+}$,
${\sigma}_{r}^{-1}$ (
${a}_{r}^{\u2020}$, *a _{r}*) at site

*r*, the total JCH Hamiltonian reads (

*h̄*= 1),.

_{〈}

_{r,s}_{〉}is the sum over all nearest-neighbor cavities,

*κ*> 0 is the hopping frequency (requiring that the ground eigenstate to be symmetric and the first excited eigenstate to be anti-symmetric also means that

*κ*> 0),

*ε*is the atomic transition energy,

*ω*is the cavity resonance frequency,

*β*is the single-photon Rabi frequency, and the rotating wave approximation is assumed. The onsite terms can be diagonalized in a basis of mixed photonic and atomic excitations called dressed states or polaritons, |±,

*n*〉

*= sin Θ*

_{r}*|*

_{n}*g*,

*n*〉

*+ cos Θ*

_{r}*|*

_{n}*e,n*− 1〉

*, with energy ${E}_{n}^{\pm}\hspace{0.17em}=\hspace{0.17em}n\omega \hspace{0.17em}-\hspace{0.17em}\Delta /2\hspace{0.17em}\pm \hspace{0.17em}\sqrt{n{\beta}^{2}\hspace{0.17em}+\hspace{0.17em}{(\Delta /2)}^{2}}$, and mixing angle ${\Theta}_{n}\hspace{0.17em}=\hspace{0.17em}\frac{1}{2}\hspace{0.17em}\text{arctan}[-2\sqrt{n}\beta /\Delta ]$, where Δ ≡ ω – ε.*

_{r}Using a Bloch state analysis (see appendix) the band structure in the one-excitation manifold is given by,

*K*≡ 4

*κ*cos(

*k*) cos(

_{x}d*k*) and for the unrotated lattice (inset of Fig. 2(b))

_{y}d*K*≡ 2

*κ*[cos(

*k*) + cos(

_{x}d*k*)], where $d\hspace{0.17em}\equiv \hspace{0.17em}|{\overrightarrow{d}}_{s}\hspace{0.17em}-\hspace{0.17em}{\overrightarrow{d}}_{r}|/\sqrt{2}$. Equation (2) is an exact solution for a periodic lattice. In our numerical simulations we will verify that boundary effects are negligible. Note that Eq. (1) describes the underlying connectivity of the lattice, and not explicitly its geometry. It is possible for geometrically distinct configurations to have the same Hamiltonian, and in the reciprocal lattice the same Bloch analysis would hold, resulting in equivalent lensing properties. However, the orientation of the lattice relative to the interface is of importance for superlensing, as is explained below.

_{y}dIn the following sections we will demonstrate how traditional metamaterial techniques can be applied to the JCH Hamiltonian, to show that a suitably designed CAM can exhibit the hallmarks of a superlens, namely all-angle negative refraction and evanescent wave ehancement.

## 3. All-angle negative refraction

A useful tool for analyzing light refraction at an interface is the isoenergy map plotted in *k*-space. Given a dispersion relation *E*(*k _{x}, k_{y}*), an isoenergy contour defines the curves over which the energy is constant. In this representation the gradient of the energy surface is the vector field of group velocities, which points normal to the isoenergy contour as illustrated by the colored arrows in Fig. 2(a). Using Eq. (2), the group velocity
${\overrightarrow{v}}_{g}^{\pm}\hspace{0.17em}\equiv \hspace{0.17em}{\nabla}_{\overrightarrow{k}}{E}^{\pm}$ is expressed as,

AANR of single photons can take place at the interface between free-space and a JCH lattice. Consider the band structure of a rotated lattice shown in Fig. 2(a) which is superimposed with the lattice and free-space isoenergy contours of matching energy. The contours associated with free-space are circles with radius equal to its energy (in natural units). By requiring that the surface parallel wave vector *k _{y}* is conserved at the interface, the group velocities associated with these contours determine the refraction angle. In the illustration, an incident photon with wavevector

*k⃗*

_{1}= (

*k*

_{1,}

*,*

_{x}*k*

_{1,}

*) and velocity*

_{y}*v⃗*

_{1}will couple to an allowed mode of the lattice, and propagate with

*k⃗*

_{2}= (

*k*

_{2,}

*,*

_{x}*k*

_{2,}

*) and*

_{y}*v⃗*

_{2}. The refraction angle is

*k*

_{2}

*is given by Eq. (13). Since the isoenergy contours of the lattice are*

_{,x}*convex*, we have negative refraction (

*θ*< 0), and since the lattice contour is larger than the air contour, this occurs for all incident angles.

_{R}To converge the light the isoenergy contour needs to be as circular as possible. This occurs at the energy band extrema, and so for a sharp focus it is preferable to work as close to these frequencies as possible.

Refraction is not invariant under lattice rotation. This is because there is a change in the air-lattice interface under rotation. From the conservation of the *k _{y}* condition, Fig. 2(b) shows that the unrotated lattice does not exhibit AANR. When AANR does occur, the propagating modes in free space can be brought into a focus to form an image on the other side of the lattice even with a planar lattice slab. Such a device therefore satisfies the first criterion of a superlens. In contrast, conventional lenses that rely on positive refraction must have curvatures to converge light.

An exact numerical simulation to test the predicted negative refraction (Eqs. (3) and (4)) for an air-lattice interface requires knowledge of, and is dependent on, the specific light-cavity coupling mechanism at the interface of a physical implementation. To demonstrate the underlying principles discussed above without recourse to a specific coupling mechanism, we conduct our numerical simulation at an interface between two JCH lattices, with no loss in generality.

We use a segmented lattice as our platform (depicted in Fig. 3(a)), where our single photon source is initialized in the ‘source’ region, and there is an identical ‘image’ region which acts as the image plane. The sandwiched ‘lens’ region negatively refracts the excitation and brings it to focus in the image plane. The band structures for the source and lens region are shown in Figs. 3(b) and 3(c) respetively. For the chosen operating frequency, indicated by the bold isoenergy contours, the system exhibits AANR.

Due to the dielectric mismatch of the interface, there must be a finite probability of reflection. At the interface, the discrete scattering eigenequation can be used to derive the reflection coefficient (see appendix),

*T*≡ 1 −

*R*. Comparing Eq. (4) and Eq. (5), there is a trade-off between refraction and reflection, i.e, large negative refraction is accompanied by large reflection. This trade-off is illustrated in a comparison plot (Fig. 4) of the refraction angle and reflection coefficient for different incident angles and varying detuning.

The propagation of the field in the lattice is governed by the Schrödinger equation |*ψ*(*t*)〉 = *e ^{iℋt}*|

*ψ*(0)〉. We consider the case when the source is initialized in an equal superposition of atomic and photonic modes. It is instructive to use a directional pulse by specifying an initial state with a normalized Gaussian momentum distribution around

*k⃗*

_{1}= (

*π*/4,±

*π*/4) so that it is incident on the lens at ±45°, as shown in Fig. 5(a). The superposition of the two

*k⃗*modes manifests in a coherent interference pattern in the

*y*-direction.

The lens atomic detuning is set to −5.27*κ*, which by Eq. (4), predicts a refraction angle of *θ _{R}* = −25°. Superimposing different time instances, the incident, reflected and refracted pulses in Fig. 5(a) follow the predicted refraction angle and the trajectory predicted by Eq. (3), to converge at a location on the image plane. The reflection and transmission coefficients are also found to be in good agreement with Eq. (20). The incident and reflected polariton (±

*π*/4,

*π*/4) coherently interfere near the interface to give an interference pattern along the

*x*-direction. Note that there is considerable reflection, so that the population density has been multiplied by a factor

*M*in Fig. 5(a) and Fig. 5(b) for clearer representation.

An important property of our system, distinct from the existing PhC superlens implementations, is the ability to tune the atomic transition energy, *ε*, after fabrication. Such manipulations can be achieved dynamically by, for example, a controlled external electric field via Stark shift. This control allows one to tailor the dispersion relation, and hence the light guiding properties and focal point of the lens. The effect of changing *ε* is demonstrated in Fig. 5(b), where decreasing Δ_{2} by 0.73*κ* shifts the focus 56 sites to the right.

To show how a point-like source would converge with AANR onto the image plane, we specify a source with an initial superposition of Gaussian momentum distribution summed over *k⃗* in Fig. 5(c). Taking a snapshot of the propagation at time, *t* = 300/*κ*, it shows that as the point source propagates into the planar lens, all components are negatively refracted, so that an image of the point source is successfully formed on the image plane. As with Fig. 5(a) and Fig. 5(b) there is considerable reflection at the lens’ interfaces.

In our examples so far, the lenses are homogenous and there is an abrupt change at the interface. As a result, there is considerable reflection such that the total transmission through the lens is less than 25%. As expressed by Eq. (5), reflection increases with the greater the change in *k _{x}*. This can be minimized if we provide an adiabatic spatial change of the atomic transition energy within the lens, in effect producing a gradient-index (GRIN) structure. In Fig. 5(d) the detuning distribution follows the form,

*x*′ is the number of sites from the interface,

*w*is the width of the GRIN region and

*W*is the total width of the lens. By fine tuning the GRIN region, the level of reflection can be made arbitrarily small, although the physical trade-off is a larger lens. The removal of reflection losses is an important development, demonstrating fine control of propagation possible in our system.

## 4. Evanescent wave enhancement

The ability of lenses to resolve images is limited by the wavelength of the light source because the high-spatial-frequency modes that describe the subwavelength features are non-propagating and do not reach the image plane. To see this, the dispersion relation in free space,
${k}_{x}\hspace{0.17em}=\hspace{0.17em}\sqrt{{\omega}^{2}\hspace{0.17em}-\hspace{0.17em}{k}_{y}^{2}}$, implies that the modes with *k _{y}* >

*ω*exponentially decay away from the source along the

*x*-axis. Existing superlens proposals overcome this diffraction limit by amplifying the evanescent wave (EW) components.

PhC-based evanescent wave enhancement (EWE) devices can be regarded as a type of resonator. At resonance, the transmission of the evanescent components is divergent. The total transmission across a lens of width *W*, derived from taking the summation of the multiple scattering events at the left and right interfaces is,

*T*(

_{ij}*R*) is the transmission (reflection) amplitude at the interface between region

_{ij}*i*and

*j*. At the resonance condition $\text{exp}(-2i{k}_{2,x}W)\hspace{0.17em}-\hspace{0.17em}{R}_{23}^{2}\hspace{0.17em}=\hspace{0.17em}0$, transmission is divergent. The resonant condition is just the condition for total internal reflection where the accumulated phase shift in a round trip is a multiple of 2

*π*.

The resonant *bound modes* [36] allow a build-up of these bound states to produce an amplified evanescent tail on the image side of the lens. These resonant modes are discrete, therefore they will not amplify all evanescent modes. NIM-based superlenses do not have this limitation, but because their fabrication is based on conducting elements they suffer the problem of loss.

Since JCH-based superlenses can dynamically shift the resonant points, they can overcome the limitation of PhC-based superlenses and amplify a contiguous range of evanescent modes (but not simultaneously).

Resonant bound modes can either exist along the interface (surface bound modes) or in the bulk (bulk bound modes). As both mechanisms follow the same underlying resonance principle, we demonstrate only the latter.

A quantum equivalent of an EW is the evanescent tail of a stationary state. We prepare our system such that the detuning in regions 1 and 3, Δ_{1} and Δ_{3} respectively, are sufficiently different from the detuning in region 2, Δ_{2} (see Fig. 6(a)), so that we can setup an eigenstate where evanescent tails exists in region 1 and 3. This setup is analogous to that of a square well. Region 4 is the lens, and region 5 our image plane.

The resonant condition for EWE can be achieved by tuning the atomic transition energy of the atoms in the lens region or the width of the lens. We will only consider the former here, as it is more appropriate to the idea of a reconfigurable device. We solve for the time evolution of the system and observe in Fig. 6(b) that significant coupling to the lens occurs at distinct values of the atomic detuning, Δ_{4}, where the resonance condition is met.

When the lens is exactly on resonance, as demonstrated in Fig. 6(c) where Δ_{2} = Δ_{4} = 0, the population exchange occurs between the source and the lens is that of coupled homogenous resonators, i.e., *P*
_{4}(*t*) = sin^{2}(Ω*t*), where Ω is the characteristic mutual coupling. Following the increase in *P*
_{4}, the EW incident on the lens is transmitted amplified by 6 orders of magnitude after a time *t* = 10^{4}
*κ*
^{−1}.

When the lens is only quasi-resonant (Fig. 6(d)), the population exchange between the coupled resonators (source and lens) can be approximated by,

*η*is the difference in the eigenenergy of the source and lens. By fitting

*P*

_{4}(

*t*) to numerical results, we find that for Δ

_{4}= 0.305

*κ*,

*η*∼10

^{−3}

*κ*. This is in good agreement with the minimum energy difference between source and lens obtained by solving the lens Hamiltonian

*ℋ*directly.

Figure 6(c) shows that at *t* = 10^{3}
*κ*
^{−1} the incident EW is amplified by a third of the exact resonance case. Thus, although exact resonance is an optimal condition for EWE it is not a necessary condition for enhancement. When the lens is tuned away from resonance, the degree of EWE can quickly diminish as seen in Fig 6(e).

The diffraction limit restricts the resolution of conventional lenses to the operating wavelength, Λ_{0}. Our lens’ resolution, *δ* = 2*π*/*k*
_{max}, is determined by the maximum *k* that still satisfies the resonant condition. However we would also like to resolve all the *k*-components leading up to *k*
_{max}. This implies minimizing the bulk energy band spectrum so that the deviation from the resonant energy is always small. The drawback of this is a reduction in sharpness of focus. A better solution is to introduce surface mode resonance. This can be achieved by having a different *ε* at the lens surfaces from that of the bulk. As shown in Fig. 7, the *flatter* surface mode band (see appendix) provides the necessary minimal deviation from resonance to maximize *k*
_{max}, leaving the bulk mode to provide the AANR and focal sharpness.

Figure 7 shows that *k*
_{max} = 2*π*/*d* − *ω*
_{0}, because beyond this the evanescent modes fold back into the light cone and the associated bound modes become *leaky* states [36]. Thus the maximum resolution of our lens is,

*d*< Λ

_{0}/2, the resolution exceeds that of conventional lenses.

Inter-cavity spacing is however limited by the size of the cavity. This means that to beat the diffraction limit, one at the very least needs the cavity size to be less than Λ_{0}. Typically cavity resonant wavelength is twice the cavity size, so the subwavelength resolution condition becomes, *ω* > *ω*
_{0}. The non-linear interaction introduced by the cavity atom allows, beyond that which is available through just inter-cavity hopping, the cavity resonance frequency to be greater than the operating frequency.

Using Eq. (2) where *E*
^{±} = *ω*
_{0}, the relative resolution of our lens can be approximated by,

*δ*

_{0}

*/δ*

^{+}< 1, only the resolution from the negative energy branch,

*δ*

^{−}, can better the diffraction limit. Conventionally

*β*/

*ω*, restricted by the so-called fine structure constant limit, is of the order 0.01 (although larger values are possible for unconventional coupling mechanisms [44]), so the improvement over the diffraction limit is typically small.

## 5. Experimental feasibility and outlook

The lead time from theoretical inception to experimental realization can be long. NIMs were originally proposed in 1968 [6], but it was not physically realised until 2000 [8]. The case for the experimental feasibility of CAMs is fortunately alot clearer. Optical microcavities have been created with whispering gallery modes (WGM) [45], fabry-perot microcavities [46, 47], photonic bandgap (PBG) defects [48], and in slot waveguides [49]. The coupling of PBG nanocavities into a one-dimensional array has also been achieved [43]. Recently the coupling of Nitrogen-vacancy (NV) centers to PBG cavities [50] and WGM microdisks [51] have been experimentally verified. With the astonishing advancement in microcavities and arrays of microcavities, and their coupling to multi-level quantum systems, it is feasible that the next advancement in experimental development would be that of arrays of microcavities coupled to multi-level quantum systems, in other words, CAMs. The dynamic control of Stark tunable solid-state emitters have been been demonstrated with diamond colour centres [52], quantum dots [53] single molecules [54], and single-ion dopants [55].

The development of microcavities can be characterized by their size and quality factor, *Q*. WGM microcavities have experimentally produced *Q ∼* 10^{9} in microspheres [56] and 10^{8} in microtoroids [57, 58]. However their sizes are relatively large: 10^{3}
*μ*m^{3} and 180 *μ*m^{3} respectively. PBG microcavities have achieved *Q ∼* 10^{7} with cavity mode volume, *V ∼* (*λ*/2)^{3}, where *λ* is the operating wavelength [59]. NV couplings in PBG microcavities have been calculated as *β* ∼ 10^{10} Hz [60]. Assuming that photon hopping limits *Q*, we can approximate the inter-cavity tunneling frequency as *κ* = *ω*/*Q*. For the superlensing properties presented in this work, *β* = 100*κ*, requiring in the visible light regime *Q ∼* 10^{7}, which is at current experimental limits. In PBG arrays with over 100 microcavities however, only *Q ∼* 10^{6} as yet been experimentally verified [43]. Coupling subwavelength-sized slot-waveguide cavities in 2D lattice has also been discussed [61].

Atom-optical microcavities have demonstrated many quantum effects ranging from quantum collapse and revival to eletromagnetically induced transparency. It is upon these quantum effects that CAMs’ true potential lies. In particular, a transmission line formed by coupling superconducting charge qubits prepared in a coherent superposition of quantum states has been studied to show an oscillating band gap and transmittance [32]. First experimental steps to realize such materials have been reported [62–64]. Along similar lines for future investigation, coupling of a three level Λ-system (e.g. using quantum dots [65] or diamond defects [66]), to the cavity modes could yield a superposition of dispersion relations and thus create a metama-terials with a quantum superposition of material properties. One of the ground states (|*g*〉) in the individual cavity Λ-system would be strongly coupled to the cavity modes while the other (|*f*〉) is not (Fig. 8). The atom-photon coupling offers intriguing potential for entirely new quantum devices. For example, if some of the atoms are prepared in a Greenberger-Horne-Zeilinger–like (GHZ-like) state |*gg...g*〉 + *|ff...f*〉, then the system will exhibit a superposition of dispersion relations. This in turn implies a superlens with two focal points in quantum superposition.

## 6. Conclusion

In summary, we have combined the previously unrelated fields of quantum mechanics and metamaterials, by proposing cavity arrays as a new class of dynamic metamaterial. Operating at the quantum level, it opens up new possibilities for quantum optical devices. By applying traditional metamaterial techniques we showed that the CAM can exhibit the features of a superlens. In a more general sense, this work lays down the framework for local manipulation of photons, the quantum superposition of metamaterial properties, the preservation and interaction with entangled fields, and other non-local effects, in cavity array metamaterials, creating a new area of investigation in quantum transformation optical phenomena. This invites quantum technology into the realm of metamaterials.

## Appendix

1. The method described here for the derivation of the band structure solution in Eq. (2) follows Ref. [23]. In the presence of intercavity coupling, the onsite energies *E _{r}* are no longer the polaritonic energies
${E}_{n}^{\pm}$ and in general, satisfy the relation,

*H*is the Hamiltonian that relates site

_{rs}*r*to site

*s*. Employing Bloch’s theorem for periodic structures,

*d⃗*is the displacement to site

_{r}*r*, and

*k⃗*≡ (

*k*,

_{x}*k*) is the wavevector associated with the crystal momentum, Eq. (11) becomes an energy eigenequation whose eigenvalues are the energy band structure or the dispersion relation of the medium.

_{y}2. For an incident field of energy *E*
_{1} and wavevector *k⃗*
_{1} = (*k*
_{1,}
* _{x}*,

*k*

_{1,}

*) with transmitted field of*

_{y}*E*

_{2}and

*k⃗*

_{2}= (

*k*

_{2,}

*,*

_{x}*k*

_{2,}

*) at the interface, energy conservation (*

_{y}*E*

_{1}=

*E*

_{2}) and phase matching (

*k*

_{1,}

*=*

_{y}*k*

_{2,}

*=*

_{y}*k*) requires that

_{y}*k*

_{2,}

*satisfies the condition,*

_{x}*K*≡ 4

*κ*cos(

*k*

_{2,}

*) cos(*

_{x}d*k*

_{1,}

*) for the rotated lattice.*

_{y}d3. For the derivation of the reflection coefficient (Eq. (5)), the state vector can be expanded in the bare atom-photon basis,

*g,n*〉

*and*

_{r}*|e*,

*n*〉

*denote the ground and excited state respectively, with*

_{r}*n*photonic excitations at site

*r*. Given the symmetry, we consider a 5-site unit cell in X-configuration that is translational-invariant along the

*y*-direction. Using the standard eigenenergy equation

*ℋ*|

*ψ*〉 =

*E*|

*ψ*〉, one arrives at a discrete scattering equation for each region of the lattice and the interfaces. In particular in region

*j*with associated parameters (

*ω*,

_{j}*ε*,

_{j}*β*,

_{j}*κ*), we have

_{j}*p,q*) is surrounded by four nearest-neighboring sites at coordinates (

*p*± 1,

*q*± 1) and the conservation of energy requires

*E*=

_{j}*E*. At the interface centered at the origin (0,0),

*k*=

_{j,p}*k*–

_{j,x}*k*and

_{y}*k*=

_{j,q}*k*+

_{j,x}*k*. In region 2, the transmitted wave has the form,

_{y}*r*and

*t*are used here to denote reflection and transmission amplitudes respectively. Substituting these solutions in to the interface equation (Eq. (16)) and applying continuity condition

*t*= 1 +

*r*, we arrive at the reflection coefficient

*R*≡ |

*r*|

^{2},

*κ*=

_{j}*κ*, we retrieve the required expression. Finally since a polaritonic pulse has a momentum distribution

*G*(

*k⃗*), we define an effective reflection coefficient,

4. To calculate the surface mode energy band, we need to take two adjacent sites along the surface as the primitive cell. After applying Bloch’s theorem the Hamiltonian is,

*ε*and

_{s}*ε*are the surface and bulk atomic transition energies respectively. We solve for the eigenvalues to get the surface mode energy bands.

_{b}## Acknowledgments

This work was supported by the Australian Research Council (ARC) under the Centre of Excellence scheme. L.C.L.H. and A.D.G. acknowledge the financial support of the ARC under Projects DP0770715 and DP0880466, respectively.

## References and links

**1. **J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart, “Magnetism from conductors and enhanced nonlinear phenomena,” IEEE Trans. Microwave Theory Tech. **47**, 2075–2084 (1999). [CrossRef]

**2. **M. C. K. Wiltshire, J. B. Pendry, I. R. Young, D. J. Larkman, D. J. Gilderdale, and J. V. Hajnal, “Microstructured magnetic materials for RF flux guides in magnetic resonance imaging,” Science **291**, 849–851 (2001). [CrossRef] [PubMed]

**3. **M. C. K. Wiltshire, J. V. Hajnal, J. B. Pendry, D. J. Edwards, and C. J. Stevens, “Metamaterial endoscope for magnetic field transfer: near field imaging with magnetic wires,” Opt. Express **11**, 709–715 (2003). [CrossRef] [PubMed]

**4. **M. C. Wiltshire, J. B. Pendry, D. J. Larkman, D. J. Gilderdale, D. Herlihy, I. R. Young, and J. V. Hajnal, “Geometry preserving flux ducting by magnetic metamaterials,” Proc. Int. Soc. Mag. Reson. Med. **11**, 713–713 (2003).

**5. **T. J. Yen, W. J. Padilla, N. Fang1, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basovm, and X. Zhang, “Terahertz magnetic response from artificial materials,” Science **303**, 1494–1496 (2004). [CrossRef] [PubMed]

**6. **V. G. Veselago, “The electrodynamics of substances with simultaneously negative values of *ε* and *μ*,” Sov. Phys. Uspekhi-USSR **10**, 509–514 (1968). [CrossRef]

**7. **T. Paul, C. Rockstuhl, C. Menzel, and F. Lederer, “Anomalous refraction, diffraction, and imaging in metamaterials,” Phys. Rev. B **79**, 115430 (2009). [CrossRef]

**8. **D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat-Nasser, and S. Schultz, “Composite medium with simultaneously negative permeability and permittivity,” Phys. Rev. Lett. **84**, 4184–4187 (2000). [CrossRef] [PubMed]

**9. **R. Shelby, D. R. Smith, and S. Schultz, “Experimental verification of a negative index of refraction,” Science **292**, 77–79 (2001). [CrossRef] [PubMed]

**10. **A. A. Houck, J. B. Brock, and I. L. Chuang, “Experimental observations of a left-handed material that obeys Snell’s law,” Phys. Rev. Lett. **90**, 137401 (2003). [CrossRef] [PubMed]

**11. **C. G. Parazzoli, R. B. Greegor, K. Li, B. E. C. Koltenbah, and M. Tanielian, “Experimental verification and simulation of negative index of refraction using Snell’s law,” Phys. Rev. Lett. **90**, 107401 (2003). [CrossRef] [PubMed]

**12. **D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science **314**, 977–980 (2006). [CrossRef] [PubMed]

**13. **Z. L. Liu, X. Zhang, Y. Mao, Y. Y. Zhu, Z. Yang, C. T. Chan, and P. Sheng, “Locally resonant sonic materials,” Science **289**, 1734–1736 (2000). [CrossRef] [PubMed]

**14. **S. Guenneau, A. Movchan, G. Ptursson, and S. A. Ramakrishna, “Acoustic metamaterials for sound focusing and confinement,” New J. Phys. **9**, 399 (2007). [CrossRef]

**15. **M. Brun, S. Guenneau, and A. B. Movchan, “Achieving control of in-plane elastic waves,” Appl. Phys. Lett. **94**, 061903 (2009). [CrossRef]

**16. **M. J. Hartmann, F. G. S. L. Brandão, and M. B. Plenio, “Strong interacting polaritons in coupled arrays of cavities,” Nat. Phys. **2**, 849–855 (2006). [CrossRef]

**17. **A. D. Greentree, C. Tahan, J. H. Cole, and L. C. L. Hollenberg, “Quantum phase transitions of light,” Nat. Phys. **2**, 856–861 (2006). [CrossRef]

**18. **D. G. Angelakis, M. F. Santos, and S. Bose, “Photon-blockade-induced Mott transitions and *XY* spin models in coupled cavity arrays,” Phys. Rev. A **76**, 031805(R) (2007). [CrossRef]

**19. **L. Zhou, H. Dong, C. P. Sun, and F. Nori, “Quantum supercavity and atomic mirrors,” Phys. Rev. A **78**, 063827 (2008). [CrossRef]

**20. **M. I. Makin, C. Cole, C. D. Tahan, L. C. L. Hollenberg, and A. D. Greentree, “Quantum phase transitions in photonic cavities with two-level systems,” Phys. Rev. A **80**, 043842 (2009). [CrossRef]

**21. **S. Schmidt and G. Blatter, “Strong coupling theory for the Jaynes-Cummings-Hubbard model,” Phys. Rev. Lett. **103**, 086403 (2009). [CrossRef] [PubMed]

**22. **P. Pippan, H. G. Evertz, and M. Hohenadler, “Excitation spectra of strongly correlated lattice bosons and polaritons,” Phys. Rev. A **80**, 033612 (2009). [CrossRef]

**23. **J. Quach, M. I. Makin, C.-H. Su, A. D. Greentree, and L. C. L. Hollenberg, “Band structure, phase transitions and semiconductor analogs in one-dimensional solid light systems,” Phys. Rev. A **80**, 063838 (2009). [CrossRef]

**24. **G. Rempe, H. Walther, and N. Klein, “Observation of quantum collapse and revival in a one-atom maser,” Phys. Rev. Lett. **58**, 353–356 (1987). [CrossRef] [PubMed]

**25. **M. Brune, F. Schmidt-Kaler, A. Maali, J. Dreyer, E. Hagley, J. M. Raimond, and S. Haroche, “Quantum Rabi oscillation: a direct test of field quantization in a cavity,” Phys. Rev Lett. **76**, 1800–1803 (1996). [CrossRef] [PubMed]

**26. **K. M. Birnbaum, A. Boca, R. Miller, A. D. Boozer, T. E. Northup, and H. J. Kimble, “Photon blockade in an optical cavity with one trapped atom,” Nature **436**, 87–90 (2005). [CrossRef] [PubMed]

**27. **M. Mucke, E. Figueroa, J. Bochmann, C. Hahn, K. Murr, S. Ritter, C. J. Villas-Boasz, and G. Rempe, “Electromagnetically induced transparency with single atoms in a cavity,” Nature **465**, 755–758 (2010). [CrossRef] [PubMed]

**28. **A. Hayward and A. D. Greentree, “Quantum and classical chaos in kicked coupled Jaynes-Cummings cavities,” Phys. Rev. A **81**, 063831 (2010). [CrossRef]

**29. **C.-H. Su, A. D. Greentree, W. J. Munro, K. Nemoto, and L. C. L. Hollenberg, “Pulse shaping by coupled cavities: single photons and qudits,” Phys. Rev. A **80**, 033811 (2009). [CrossRef]

**30. **L. Zhou, Z. R. Gong, Y.-X. Liu, C. P. Sun, and F. Nori, “Controllable scattering of a single photon inside a one-dimensional resonator waveguide,” Phys. Rev. Lett. **101**, 100501 (2008). [CrossRef] [PubMed]

**31. **A. Tomadin and R. Fazio, “Many-body phenomena in QED-cavity arrays,” J. Opt. Soc. Am. B **27**, A130–A136 (2010). [CrossRef]

**32. **A. L. Rakhmanov, A. M. Zagoskin, S. Savelev, and F. Nori, “Quantum metamaterials: electromagnetic waves in a Josephson qubit line,” Phys. Rev. B. **77**, 144507 (2008). [CrossRef]

**33. **J. B. Pendry, “Negative refraction makes a perfect lens,” Phys. Rev. Lett. **85**, 3966–3969 (2000). [CrossRef] [PubMed]

**34. **X. Zhang and Z. Liu, “Superlenses to overcome the diffraction limit,” Nature Mater. **7**, 435–441 (2008). [CrossRef]

**35. **C. Luo, S. G. Johnson, J. D. Joannopoulos, and J. B. Pendry, “All-angle negative refraction without negative effective index,” Phys. Rev. B **65**, 201104(R) (2002). [CrossRef]

**36. **C. Luo, S. G. Johnson, J. D. Joannopoulos, and J. B. Pendry, “Subwavelength imaging in photonic crystals,” Phys. Rev. B **68**, 045115 (2003). [CrossRef]

**37. **C. Luo, S. G. Johnson, J. D. Joannopoulos, and J. B. Pendry, “Negative refraction without negative index in metallic photonic crystals,” Opt. Express **11**, 746–754 (2003). [CrossRef] [PubMed]

**38. **A. Grbic and G. V. Eleftheriades, “Overcoming the diffraction limit with a planar left-handed transmission-line lens,” Phys. Rev. Lett. **92**, 117403 (2004). [CrossRef] [PubMed]

**39. **E. Cubukcu, K. Aydin, E. Ozbay, S. Foteinopoulou, and C. M. Soukoulis, “Subwavelength resolution in a two-dimensional photonic-crystal-based superlens,” Phys. Rev. Lett. **91**, 207401 (2003). [CrossRef] [PubMed]

**40. **E. Cubukcu, K. Aydin, E. Ozbay, S. Foteinopoulou, and C. M. Soukoulis, “Electromagnetic waves: negative refraction by photonic crystals,” Nature **423**, 604–605 (2003). [CrossRef] [PubMed]

**41. **Z. Liu, N. Fang, T.-J. Yen, and X. Zhang, “Rapid growth of evanescent wave with a silver superlens,” Appl. Phys. Lett. **83**, 5184–5186 (2003). [CrossRef]

**42. **T. Taubner, D. Korobkin, Y. Urzhumov, G. Shvets, and R. Hillenbrand, “Near-field microscopy through a SiC superlens,” Science **313**, 1595–1595 (2006). [CrossRef] [PubMed]

**43. **M. Notomi, E. Kuramochi, and T. Tanabe, “Large-scale arrays of ultrahigh-Q coupled nanocavities,” Nat. Photonics **2**, 741–747 (2008). [CrossRef]

**44. **M. Devoret, S. Girvin, and R. Schoelkopf, “Circuit-QED: how strong can the coupling between a Josephson junction atom and a transmission line resonator be?,” Ann. Phys. (Leipzig) **16**, 767–779 (2007). [CrossRef]

**45. **V. Lefevre-Seguin and S. Haroche, “Towards cavity-QED experiments with silica microspheres,” Mat. Sci. Eng. B **48**, 53–58 (1997). [CrossRef]

**46. **C. J. Hood, H. J. Kimble, and J. Ye, “Characterization of high-finesse mirrors: loss, phase shifts, and mode structure in an optical cavity,” Phys. Rev. A **64**, 033804 (2001). [CrossRef]

**47. **M. Kohnen, M. Succo, P. G. Petrov, R. A. Nyman, M. Trupke, and E. A. Hinds, “An array of integrated atom-photon junctions,” Nat. Photonics **5**, 35–38 (2011). [CrossRef]

**48. **C. Reese, B. Gayral, B. D. Gerardot, A. Imamoglu, P. M. Petroff, and E. Hu, “High-Q photonic crystal microcavities fabricated in a thin GaAs membrane,” J. Vac. Sci. Technol. B **19**, 2749–2752 (2001). [CrossRef]

**49. **R. Sun, P. Dong, N.-N. Feng, C.-Y. Hong, J. Michel1, M. Lipson, and L. Kimerling, “Horizontal single and multiple slot waveguides: optical transmission at *λ* = 1550 nm,” Opt. Express **15**, 17967–17972 (2007). [CrossRef] [PubMed]

**50. **M. Barth, N. Nüsse, B. Löchel, and O. Benson, “Controlled coupling of a single-diamond nanocrystal to a photonic crystal cavity,” Opt. Lett. **34**, 1108–1110 (2009). [CrossRef] [PubMed]

**51. **P. E. Barclay, C. Santori, K.-M. Fu, R. G. Beausoleil, and O. Painter, “Coherent interference effects in a nano-assembled diamond NV center cavity-QED system,” Opt. Express **17**, 8081–8097 (2009). [CrossRef] [PubMed]

**52. **Ph. Tamarat, T. Gaebel, J. R. Rabeau, M. Khan, A. D. Greentree, H. Wilson, L. C. L. Hollenberg, S. Prawer, P. Hemmer, F. Jelezko, and J. Wrachtrup, “Stark shift control of single optical centers in diamond,” Phys. Rev. Lett. **97**, 083002 (2006). [CrossRef] [PubMed]

**53. **S. A. Empedocles and M. G. Bawendi, “Quantum-confined Stark effect in single CdSe nanocrystallite quantum dots,” Science **278**, 2114–2117 (1997). [CrossRef]

**54. **Ch. Brunel, Ph. Tamarat, B. Lounis, J. C. Woehl, and M. Orrit, “Stark effect on single molecules of Dibenzan-thanthrene in a Naphthalene crystal and in a n-Hexadecane Shpol’skii matrix,” J. Phys. Chem. A **103**, 2429–2434 (1999). [CrossRef]

**55. **A. L. Alexander, J. J. Longdell, M. J. Sellars, and N. B. Manson, “Photon echoes produced by switching electric fields,” Phys. Rev. Lett. **96**, 043602 (1006). [PubMed]

**56. **M. L. Gorodetsky, A. A. Savchenkov, and V. S. Ilchenko, “Ultimate Q of optical microsphere resonators,” Opt. Lett. **21**, 453–455 (1996). [CrossRef] [PubMed]

**57. **D. K. Armani, T. J. Kippenberg, S. M. Spillane, and K. J. Vahala, “Ultrahigh-Q toroid microcavity on a chip,” Nature **421**, 925–928 (2003). [CrossRef] [PubMed]

**58. **S. M. Spillane, T. J. Kippenberg, K. J. Vahala, K. W. Goh, E. Wilcut, and H. J. Kimble, “Ultrahigh-Q toroidal microresonators for cavity quantum electrodynamics,” Phys. Rev. A **71**, 013817 (2005). [CrossRef]

**59. **S. Noda, M. Fujita, and T. Asano, “Spontaneous-emission control by photonic crystals and nanocavities,” Nat. Photonics **1**, 449–458 (2007). [CrossRef]

**60. **M. P. Hiscocks, C.-H. Su, B. C. Gibson, A. D. Greentree, L. C. L. Hollenberg, and F. Ladouceur, “Slot-waveguide cavities for optical quantum information applications,” Opt. Express **17**, 7295–7303 (2009). [CrossRef] [PubMed]

**61. **C.-H. Su, M. P. Hiscocks, B. C. Gibson, A. D. Greentree, L. C. L. Hollenberg, and F. Ladouceur, “Coupling slot-waveguide cavities for large-scale quantum optical devices,” Opt. Express **19**, 6362–6373 (2011). [CrossRef]

**62. **A. A. Abdumalikov Jr., O. Astafiev, A. M. Zagoskin, Yu. A. Pashkin, Y. Nakamura, and J. S. Tsai, “Electromagnetically induced transparency on a single artificial atom,” Phys. Rev. Lett. **104**, 193601 (2010).

**63. **O. V. Astafiev, A. A. Abdumalikov Jr., A. M. Zagoskin, Yu. A. Pashkin, Y. Nakamura, and J. S. Tsai, “Ultimate on-chip quantum amplifier,” Phys. Rev. Lett. **104**, 183603 (2010).

**64. **O. Astafiev, A. M. Zagoskin, A. A. Abdumalikov Jr., Yu. A. Pashkin, T. Yamamoto, K. Inomata, Y. Nakamura, and J. S. Tsai, “Resonance fluorescence of a single artificial atom,” Science **327**, 840–843 (2010). [CrossRef] [PubMed]

**65. **X. Xu, B. Sun, P. R. Berman, D. G. Steel, A. S. Bracker, D. Gammon, and L. J. Sham, “Coherent population trapping of an electron spin in a single negatively charged quantum dot,” Nat. Phys. **4**, 692–695 (2008). [CrossRef]

**66. **C. Santori, P. Tamarat, P. Neumann, J. Wrachtrup, D. Fattal, R. G. Beausoleil, J. Rabeau, P. Olivero, A. D. Greentree, S. Prawer, F. Jelezko, and P. Hemmer, “Coherent population trapping of single spins in diamond under optical excitation,” Phys. Rev. Lett. **97**, 247401 (2006). [CrossRef]