## Abstract

Cavity array metamaterials (CAMs), composed of optical microcavities in a lattice coupled via tight-binding interactions, represent a novel architecture for engineering metamaterials. Since the size of the CAMs’ constituent elements are commensurate with the operating wavelength of the device, it cannot directly utilise classical transformation optics in the same way as traditional metamaterials. By directly transforming the internal geometry of the system, and locally tuning the permittivity between cavities, we provide an alternative framework suitable for tight-binding implementations of metamaterials. We develop a CAM-based cloak as the case study.

© 2013 OSA

## 1. Introduction

Metamaterials, as an arrangement of artificial elements designed to achieve novel electromagnetic properties, represent a major advancement in material science. Typically constructed with conductive elements smaller than the operating wavelength, these structures have demonstrated novel electromagnetic behaviour, such as negative refraction [1] and cloaking [2]. Photonic crystals (PhCs) may also mimic metamaterial behaviour [3–5]. Tight-binding models have recently been investigated as a novel metamaterial architecture [6, 7]. A realization of this is the cavity array metamaterial (CAM) – an optical lattice composed of evanescently-coupled microcavities [6]. The quantum mechanical properties of the metamaterial are produced by coupling each cavity to single atoms that can be individually controlled. By tailoring the system’s dispersion properties, it can show effects generally found in classical metamaterials, such as negative refraction, evanescent wave amplification, and anomalous waveguiding properties, but with the distinct exception that it operates fundamentally at the quantum level. Tuning the atomic properties opens up the possibility of reconfigurable metamaterials with quantum properties.

Pivotal in the development of metamaterials was the advent of transformation optics (TO) as a methodology for designing novel electromagnetic devices, paving the way to the realization of a cloaking device [2]. Critically, because CAMs operate on a wavelength commensurate with the size of its constituent elements, they cannot directly utilise classical TO in the same way as traditional metamaterial systems. Here we develop a framework for TO using CAMs.

Traditional TO works by applying a covariant co-ordinate transformation to Maxwell’s
equations and then finding the required permittivity, *ε*, and permeability,
*μ*, for light to follow geodesics in this transformed co-ordinate system
which in the original co-ordinate system would be curved paths [8, 9]. An attempt has been made
to find a quantum equivalent that involves an invariant transformation of the Schrödinger
equation, predicting cloaking of matter waves [10]. As a fundamental technique this work may find applications in systems such as
ultracold rubidium atoms in optical lattices where the wavelength of the atom is many times greater
than the lattice spacing so the system can be treated as an effective medium, although there are
significant technical difficulties that need to be overcome first.

PhCs as a metamaterial have the advantage that electromagnetic waves can propagate through them with relatively low loss. Classically, PhCs can effect TO by either varying the dielectric composition [4] or geometrically transforming the PhC [5]. Our approach is in the spirit of the latter, in that we physically deform the medium. However our system differs from [5] because we are working with a quantum mechanical platform. Unlike [10] we work with photons and not matter waves. In particular, we will look at the propagation of a single photon.

In Sec. 2 we discuss the guiding of light in CAMs to effect TO, using the cloaking of a circular region as a case study. We then calculate in Sec. 3 the required permittivity distribution in a coupled optical array to implement the circular cloak.

## 2. Cloaking

In CAMs, light propagation can be directly controlled by the lattice geometry and intersite hopping frequency. Therefore, unlike the continuous medium of classical metamaterials, coordinate transformations can be directly represented by changing the lattice configuration. The goal then is to keep the Hamiltonian invariant whilst the lattice configuration changes according to the required propagation response. Thus TO applied here to CAMs is a literal transformation of the medium space.

To demonstrate TO in CAMs we consider an array of coupled optical cavities, the simplest of the CAMs, and show how it can be used as a cloak. This is schematically shown in Fig. 1(a). This choice of system, which has become increasingly accessible in experiments [11], suffices to demonstrate the effects of modifying geometry. The methods here are transferable to other more complex CAMs such as the Jaynes-Cummings-Hubbard (JCH) system [6].

Consider a two-dimensional (2D) array of interacting identical cavities, described by the
Hamiltonian (*h*̄ = 1),

*i*, ∑

_{〈}

_{i, j}_{〉}is the sum over all nearest-neighbor cavities,

*ω*is the cavity resonance frequency, and

*κ*is the inter-cavity coupling rate. We have used a second quantized description as it provides a convenient way to express the system, and it allows for a direct translation to more complex quantum mechanical systems, such as the JCH model. A classical electromagnetic description of the problem however, would be equally valid. Note that as the system is linear, and we are working explicitly with single photons, our results apply equally to many photons, where Hong-Ou-Mandel effects [12] are ignored.

For a uniform square lattice, Bloch’s theorem can be applied to find an analytical expression for the dispersion relation [13],

where*k⃗*=

*k*̂ +

_{x}x*k*is the wave vector associated with the

_{y}ŷ*crystal momentum*and

*d*is the distance between neighboring sites. The dispersion surface of the system is shown depicted in Fig. 2. The behavior of the propagating modes is then governed by the group velocity,

*r*′,

*ϕ*′),

*r*∈ (0,

*b*), where

*a*is the radius of the hidden region and

*b*−

*a*is the radial width of the cloak. The transformation directly relocates the relative positions of the cavities within this region, as shown in Fig. 3(b).

To illustrate the cloaking property of the system, we solve the time evolution of the system, *ψ̇*(*t*) = −*iℋψ*(*t*), in the one-excitation (single-photon) manifold, with two different initial conditions. We use isofrequency contours and group velocities given by Eqs. (2) and (3) to identify operating frequencies. In Fig. 3(a), we initialize the excitation in a superposition of all the *k*-modes at an isofrequency contour near the band edge (white curve in Fig. 2), as this is where it is most circular, to mimic a point-like source: we plot the probability distribution |*ψ*(*t*)|^{2}. In Fig. 3(c), the excitation is initialized as a Gaussian pulse in a single *k*-mode for a square contour (black bold curve in Fig. 2) to minimize dispersion, and plot the superimposed distribution ∑* _{t}* |

*ψ*(

*t*)|

^{2}. For both cases, the photon propagates around the hole and continues as if the hole was not there. Any detector outside the annulus region of deformation would not be able to distinguish between the original system and the cloaking system, thereby allowing anything inside the hole to be hidden from view.

Unlike traditional transformation optics where a change in local *ε* and *μ* are required to effect the co-ordinate transformation, in CAMs the co-ordinate transformation is a direct transformation of the physical configuration of the system, and the inter-site properties are tuned to ensure the invariance of the Hamiltonian (Eq. 1). In other words, the inter-cavity coupling for different inter-cavity spacing needs to remain the same for all sites *i* and *j* (*κ _{ij}* =

*κ*, ∀

*). In the next section we show how this can be achieved by tuning the intercavity permittivity.*

_{i,j}## 3. Permittivity detuning

Consider a 1D coupled cavity system with individual cavity width *w*, and permittivity profile described by step-functions such that *ε*(*x*) = *ε _{a}* inside cavities and

*ε*in-between cavities, as depicted in Fig. 1(b). Solving Maxwell’s equation for a single cavity ∇

_{b}^{2}

**E**

_{Ω}= −

*ε*(

*x*)

*ω*

^{2}/

*c*

^{2}

**E**

_{Ω}yields the fundamental cavity eigenmode with,

*A*satisfies ∫

*ε*

_{Ω}(

*x*)|

*E*

_{Ω}(

*x*)|

^{2}

*dx*= 1. In a tight-binding model, the nearest-neighbour coupling is given by [14],

*ℋ*invariant from lattice deformation, we solve Eq. (6) for a constant inter-cavity coupling while allowing

*ε*to vary.

_{b}Optical microcavities such as microdisks, pillars and toroids in general, and 2D PhC resonators in particular, represent a versatile platform for realizing this scheme. PhC cavities are created by introducing a local inhomogeneity in a periodic dielectric lattice. By placing such defects in close proximity, they become evanescently coupled to form a coupled-cavity array. Alternatively, this coupling can be achieved through a line-defect, which has demonstrated resonant transmission across over 100 cavities in a thin silicon PhC membrane [11]. In both cases, the effective inter-site permittivity (*ε _{b}*) is related to the refractive indices of the dielectric medium and hole infiltrations [15], and also the holes’ size, geometry and position [16] in the crystal regions between the cavities.

That the index of liquid crystal infiltrated into air pores of PhC can be altered substantially using an applied electric field [15] [Fig. 1(c)], represents a means to dynamically tune the local inter-cavity coupling. Mechanisms to change the index of the background material have also been developed including local heating [17], chemical etching [18], electrical carrier injections [19], near-field tip perturbation [20], photosensitive material illumination [21], carbon dot deposition [22], fiber-taper probing [23], ion-beam irradiation [24], and nanoparticle intrusion [25]. Controlling the atomic-cavity detuning in coupled atom-optical cavities has also been shown to play a similar role to local index modification [6].

As a specific case study, we consider an infiltrated silicon-based PhC [Fig. 1(d)] implementation operating at *λ* = 1.5 *μ*m with typical setup characteristics of *ε _{a}* = 11.7,

*ε*≈ 2.3 (prior to transformation),

_{b}*κ*= 10

^{14}rad/s, and

*w*=

*λ*/2. To achieve the above cloaking prescription with uniform

*ω*and

*ε*, we require the distribution of the inter-site permittivity to take up values shown in Fig. 4 for the example in Fig. 3. The permittivity in-between neighboring sites is assumed to be constant. Each plot point in Fig. 4b) represents an inter-site position, and is characterized by a radial distance and a permittivity value. Near the center of the cloak the deformation stretches the physical distance between nearest neighbors. To maintain the same hopping frequency, the intersite permittivity is decreased. Near the outer radius of the cloaking annulus, the sites are squashed closer together, and correspondingly the intersite permittivity here is greater than the baseline permittivity outside the cloaking region.

_{a}## 4. Conclusion and outlook

CAMs offers the potential of metamaterials with quantum properties. Here we have shown how transformation optics is implemented in the context of CAMs: a spatial transformation of the physical sites and careful tuning of the intersite coupling. As a case study, we showed how a single photon device can cloak a circular region. This work provides the framework for future investigation into quantum mechanical behaviour such as superposition and entanglement, in this system. These types of quantum mechancal behaviour could, for example, be achieved by coherently manipulating the states of atoms or atom-like systems coupled to the cavity sites.

## Acknowledgments

We acknowledge useful discussions with Andrew. M. Martin, Andrew C. Hayward, and Lloyd C.L. Hollenberg. J.Q.Q. would like to thank Susan M. Heywood for support and general discussion. This work was supported by the Australian Research Council (ARC) under the Centre of Excellence scheme. A.D.G. acknowledges the financial support of the ARC under Project DP0880466.

## References and links

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

**2. **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]

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

**4. **Y. A. Urzhumov and D. R. Smith, “Transformation optics with photonic band gap media,” Phys. Rev. Lett. **105**, 163901–163905 (2010). [CrossRef]

**5. **Z. Liang and J. Li, “Scaling two-dimensional photonic crystals for transformation optics,” Opt. Express **19**16821–16829 (2011). [CrossRef] [PubMed]

**6. **J. Q. Quach, C.-H. Su, A. M. Martin, A. D. Greentree, and L. C. L. Hollenberg, “Reconfigurable quantum metamaterials”, Opt. Express **19**, 11018–11033 (2011). [CrossRef] [PubMed]

**7. **C.-H. Su, “Novel quantum technology based on atom-cavity physics”, Ph.D. thesis, The University of Melbourne, Victoria (2010).

**8. **U. Leonhardt, “Optical conformal mapping”, Science **312**, 1777–1780 (2006). [CrossRef] [PubMed]

**9. **P. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields”, Science **312**, 1780–1782 (2006). [CrossRef] [PubMed]

**10. **S. Zhang, D. A. Genov, C. Sun, and X. Zhang, “Cloaking of matter waves”, Phys. Rev. Lett. **100**, 123002 (2008). [CrossRef] [PubMed]

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

**12. **C. K. Hong, Z. Y. Ou, and L. Mandel, “Measurement of subpicosecond time intervals between two photons by interference”, Phys. Rev. Lett. **59**, 2044–2046, (1987). [CrossRef] [PubMed]

**13. **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]

**14. **A. Yariv, Y. Xu, R. K. Lee, and A. Scherer, “Coupled-resonator optical waveguide: a proposal and analysis,” Opt. Lett. **24**, 771–713 (1999). [CrossRef]

**15. **K. Busch and S. John, “Liquid-crystal photonic-band-gap materials: the tunable electromagnetic vacuum” Phys. Rev. Lett. **83**, 967 (1999). [CrossRef]

**16. **E. Kuramochi, M. Notomi, S. Mitsugi, A. Shinya, T. Tanabe, and T. Watanabe, “Ultrahigh-Q photonic crystal nanocavities realized by the local width modulation of a line defect,” Appl. Phys. Lett. **88**, 041112 (2006). [CrossRef]

**17. **A. Faraon and J. Vučković, “Local temperature control of photonic crystal devices via micron-scale electrical heaters,” Appl. Phys. Lett. **95**, 043102 (2009). [CrossRef]

**18. **D. Dalacu, S. Frédérick, P. J. Poole, G. C. Aers, and R. L. Williams, “Postfabrication fine-tuning of photonic crystal microcavities in InAs/InP quantum dot membranes,” Appl. Phys. Lett. **87**, 151107 (2005). [CrossRef]

**19. **D. Englund, B. Ellis, E. Edwards, T. Sarmiento, J. S. Harris, D. A. B. Miller, and J. Vuckovic, “Electrically controlled modulation in a photonic crystal nanocavity,” Opt. Express **17**, 15409–15419 (2009). [CrossRef] [PubMed]

**20. **G. Le Gac, A. Rahmani, C. Seassal, E. Picard, E. Hadji, and S. Callard, “Tuning of an active photonic crystal cavity by an hybrid silica/silicon near-field probe,” Opt. Express **17**, 21672–21679 (2009). [CrossRef] [PubMed]

**21. **A. Faraon, D. Englund, D. Bulla, B. Luther-Davies, B. J. Eggleton, N. Stoltz, P. Petroff, and J. Vuckovic, “Local tuning of photonic crystal cavities using chalcogenide glasses,” Appl. Phys. Lett. **92**, 043123 (2008). [CrossRef]

**22. **M.-K. Seo, H.-G. Park, J.-K. Yang, J.-Y. Kim, S.-H. Kim, and Y.-H. Lee, “Controlled sub-nanometer tuning of photonic crystal resonator by carbonaceous nano-dots,” Opt. Express **16**, 9829–9837 (2008). [CrossRef] [PubMed]

**23. **G. Shambat, K. Rivoire, J. Lu, F. Hatami, and J. Vučković, “Tunable-wavelength second harmonic generation from GaP photonic crystal cavities coupled to fiber tapers,” Opt. Express **18**, 12176–12184 (2010). [CrossRef] [PubMed]

**24. **S. Tomljenovic-Hanic, A. D. Greentree, C. M. de Sterke, and S. Prawer, “Flexible design of ultrahigh-Q microcavities in diamond-based photonic crystal slabslexible design of ultrahigh-Q microcavities in diamond-based photonic crystal slabs,” Opt. Express **18**, 6465–6475 (2009). [CrossRef]

**25. **S. Tomljenovic-Hanic, A. D. Greentree, B. C. Gibson, T. J. Karle, and S. Prawer, “Nanodiamond induced high-Q resonances in defect-free photonic crystal slabs,” Opt. Express **19**, 22219–22226 (2011). [CrossRef] [PubMed]