Skip to main content

Quantum metamaterials in the microwave and optical ranges

Abstract

Quantum metamaterials generalize the concept of metamaterials (artificial optical media) to the case when their optical properties are determined by the interplay of quantum effects in the constituent ‘artificial atoms’ with the electromagnetic field modes in the system. The theoretical investigation of these structures demonstrated that a number of new effects (such as quantum birefringence, strongly nonclassical states of light, etc.) are to be expected, prompting the efforts on their fabrication and experimental investigation. Here we provide a summary of the principal features of quantum metamaterials and review the current state of research in this quickly developing field, which bridges quantum optics, quantum condensed matter theory and quantum information processing.

1 Introduction

The turn of the century saw two remarkable developments in physics. First, several types of scalable solid state quantum bits were developed, which demonstrated controlled quantum coherence in artificial mesoscopic structures [14] and eventually led to the development of structures, which contain hundreds of qubits and show signatures of global quantum coherence (see [5, 6] and references therein). In parallel, it was realized that the interaction of superconducting qubits with quantized electromagnetic field modes reproduces, in the microwave range, a plethora of effects known from quantum optics (in particular, cavity QED) with qubits playing the role of atoms (‘circuit QED’, [79]). Second, since Pendry [10] extended the results by Veselago [11], there was an explosion of research of classical metamaterials resulting in, e.g., cloaking devices in microwave and optical range [1214]. The logical outcome of this parallel development was to ask, what would be the optical properties of a ‘quantum metamaterial’ - an artificial optical medium, where the quantum coherence of its unit elements plays an essential role?

As could be expected, this question was arrived at from the opposite directions, and the term ‘quantum metamaterial’ was coined independently and in somewhat different contexts. In Refs. [1517] it was applied to the plasmonic properties of a stack of 2D layers, each of them thin enough for the motion of electrons in the normal direction to be completely quantized. Therefore ‘the wavelike nature of matter’ had to be taken into account at a single-electron level, but the question of quantum coherence in the system as a whole did not arise. In Refs. [18, 19] the starting point was the explicit requirement that the system of artificial atoms (qubits) maintained quantum coherence on the time scale of the electromagnetic pulse propagation across it, in the expectation that the coherent quantum dynamics of qubits interacting with the electromagnetic field governed the ‘optical’ properties of the metamaterial.

Currently the term ‘quantum metamaterial’ is being used in both senses (see, e.g., [2025]). We will follow the more restrictive usage and call quantum metamaterials (in the narrow sense) such artificial optical (in the broad sense) media that [23] (i) are comprised of quantum coherent unit elements with desired (engineered) parameters; (ii) quantum states of (at least some of) these elements can be directly controlled; and (iii) can maintain global coherence for the duration of time, exceeding the traversal time of the relevant electromagnetic signal. The totality of (i)-(iii) (in short: controlled macroscopic quantum coherence) that makes a quantum metamaterial a qualitatively different system, with a number of unusual properties and applications.

A conventional metamaterial can be described by effective macroscopic parameters, such as its refractive index. (The requirement that the size of a unit cell of the system be much less - in practice at least twice less - than the wavelength of the relevant electromagnetic signal, is implied in its definition as an optical medium, and is inherited by quantum metamaterials.) From the microscopic point of view, these parameters are functions of the appropriately averaged quantum states of individual building blocks. In a quantum metamaterial, these states can be directly controlled and maintain phase coherence on the relevant spatial and temporal scale.

The full treatment of such a system should start from quantum description of both the electromagnetic field and the ‘atoms’. In case when their role is played by qubits (that is, two-level quantum systems), a general enough Hamiltonian of a quantum metamaterial is given by [21, 26]

$$ \hat{H} = \sum_{j} \hbar\omega_{j} b_{j}^{+}b_{j} - \frac{1}{2} \sum _{k} \bigl(\epsilon_{k}\sigma_{z}^{(k)} + \Delta_{k} \sigma_{x}^{(k)} \bigr) + i \sum _{jk} \xi_{jk}\sigma_{z}^{(k)} \bigl(b_{j}-b_{j}^{+} \bigr), $$
(1)

where the first term describes unperturbed photon modes, the second the qubit degrees of freedom, and the third their interaction. The direct control over (at least some) qubits is realized through the qubit Hamiltonian parameters (the bias ϵ or, sometimes, the tunneling matrix element Δ).

Forming photon wave packets with characteristic size \(\Lambda\gg a\) (where a is the unit cell size) and averaging over the quantum states of qubits on the scale of Λ, one should eventually arrive at the effective equation of motion for the ‘Λ-smooth’ density matrix. It will describe the state of both the electromagnetic field and the quantum metamaterial, characterized by a nondiagonal, nonlocal, state- and position-dependent ‘refractive index’ matrix.

Following through with this program involves significant technical difficulties, and the task is not brought to conclusion yet. Nevertheless certain key effects in quantum metamaterials can be investigated at a more elementary level using an approximate wave function (e.g., [2629]), or treating the electromagnetic field classically [18, 19, 30, 31]. The latter is like the standard quasi-classical treatment of the atom-light interaction [32], but is more conveniently done using the lumped-elements description (see, e.g., [30], Section 2.3). The state of a system of M nodes connected by capacitors, inductors and, if necessary, Josephson junctions, is described by a Lagrangian \({\mathcal {L}}(\{\Phi,\dot{\Phi}\})\). Here the ‘node fluxes’ \(\Phi_{j}(t) = c\int^{t} dt' V_{j}(t')\) are related to the node voltages \(V_{j}(t)\) and completely describe the classical electromagnetic field degrees of freedom (the current-voltage distribution) in the system. The lumped-elements description is appropriate, since we are interested in signals with a wavelength much larger than the dimension of a unit element of the circuit (the condition of its serving as a metamaterial).

Qubits are introduced in this scheme through their own Hamiltonians and coupling terms (e.g., flux qubits can be coupled inductively - through the magnetic flux penetrating their loop, - or galvanically - through sharing a current-carrying conductor with the circuit), which are usually cast in the form

$$ \hat{H}_{q} = -\frac{1}{2} \bigl[\epsilon_{q} \sigma_{q}^{z} + \Delta_{q} \sigma _{q}^{x} \bigr] + J_{q}(\Phi,\dot{\Phi}) \sigma_{q}^{z}. $$
(2)

The ‘quantum Routhian’,

$$ \hat{\mathcal {R}}\bigl(\{\sigma\},\{\Phi,\dot{\Phi}\}\bigr) = \sum _{q} \hat{H}_{q} - {\mathcal {L}}\bigl(\{\Phi,\dot{\Phi} \}\bigr), $$
(3)

plays the role of the Hamiltonian for qubits (in Heisenberg representation), and its expectation value \({\mathcal {R}} =\langle\hat{\mathcal {R}} \rangle\) produces the equations of motion for the classical field variables Φ, where we have included the dissipative function \(\mathcal {Q}\) to take into account resistive losses in the circuit (see Refs. [30], Section 2.3.3, and [33]):

$$ \frac{d}{dt} \frac{\partial{ \mathcal {R}}}{\partial\dot{\Phi}_{j}} - \frac {\partial{ \mathcal {R}}}{\partial\Phi_{j}} = \frac{\partial{ \mathcal {Q}}}{\partial \dot{\Phi}_{j}}. $$
(4)

If one approximates the quantum state of the qubit subsistem by a factorized function, one can take the continuum limit and, switching to the Schrödinger representation, obtain coupled systems of equations for the field variable \(\Phi(\mathbf{r},t)\) and the qubit two-component ‘macroscopic wave function’ \(\hat{\Psi}(\mathbf{r},t)\). At any convenient stage the electromagnetic modes can be canonically quantized via \(\Phi \to\hat{Q} \sim(b+b^{\dagger})\), \(\dot{\Phi} \to\hat{P} \sim i(b-b^{\dagger })\), leading back to a certain approximation of an approach based directly on Eq. (1).

2 Superconducting quantum metamaterials

The above scheme is most often applied - though not limited - to the case of superconducting metamaterials based on various types of superconducting qubits (see, e.g., Ref. [30], Chapter 2). The simplest case is given by the experimental setup of Ref. [34] (Figure 1(a)). There a single artificial atom (a flux qubit) is placed in a transmission line, and transmission and reflection coefficients for the microwave signal are measured. Here the equations (4) for the field in the continuum limit will yield free telegraph equations for the voltage and current everywhere except the point \(x=0\), where the qubit is situated:

$$\begin{aligned} \frac{\partial V(x,t)}{\partial x} = \frac{\tilde{L}}{c^{2}} \frac {\partial I(x,t)}{\partial t}; \end{aligned}$$
(5)
$$\begin{aligned} \frac{\partial I(x,t)}{\partial x} = \tilde{C} \frac{\partial V(x,t)}{\partial t}, \end{aligned}$$
(6)

where , are the inductance and capacitance per unit length of the transmission line. Of course, for such a simple structure these equations can be written down directly. The influence of the flux qubit, which is coupled to the transmission line through the effective mutual inductance M (taking into account both magnetic and kinetic inductance) is through the matching conditions at \(x=0\),

$$\begin{aligned}& V(+0,t) = V(-0,t) - \frac{M}{c}\frac{\partial\langle\hat{I}_{q} (t) \rangle}{\partial x}; \end{aligned}$$
(7)
$$\begin{aligned}& I(+0,t) = I(-0,t). \end{aligned}$$
(8)

The qubit current operator \(\hat{I}_{q} = I_{p}\sigma^{z}\) is governed by the qubit Hamiltonian (2) with the coupling term \(J_{q} = I(0,t)M/c^{2}\). An explicit solution for the reflection/transmission amplitudes was found to be in a very good agreement with the experimental data in Ref. [34] (Figure 1(b)).

Figure 1
figure 1

Interaction of electromagnetic wave with an artificial atom in one dimension. (a) A ‘proto-metamaterial’: a single artificial atom (flux qubit) in a transmission line. (b) Reflection coefficient of a microwave signal in the line measured (left) and calculated (right) for the structure of (a). After Ref. [34].

For a structure, where the role of the ‘artificial atom’ between 1D transmission lines is played by a single qubit surrounded by an array of N coupled photonic cavities [29] the calculations in the one-excitation approximation (with electromagnetic modes treated quantum mechanically) show that in such a structure arise long-living quasi-bound states of photons and the qubit, manifested as ultra-narrow resonances in the transmission coefficient.

Going from these ‘proto-metamaterials’ to QMMs containing many artificial atoms, we return to the classical treatment of the electromagnetic field. Though 2D and 3D versions of superconducting QMMs are feasible [23, 24] (Figure 2(a)), they do not operate yet in a quantum coherent regime. Most of the research at the moment concentrates on the 1D case, which already promises interesting results. In theoretical papers [18, 31, 35] is considered a QMM formed by a set of superconducting charge qubits placed in the transmission line (Figure 2(b)). The equations of motion for the field and qubits were solved, using a factorized approximation of the quantum state vector of the qubit subsystem. For the realistic choice of qubit and transmission line parameters, the figures of merit \(\beta= \sqrt{E_{EM}/ \Delta} \sim30\) and \(\nu= \hbar\max(\Gamma_{qb},\Gamma_{TL})/\Delta \sim10^{-3}\) (Ref. [18]) ensure that the continuum approximation is justified, and that in the first approximation the effects of decoherence can be neglected. Here \(E_{EM}\) is the electromagnetic field energy per unit cell, Δ (Eq. (1)) gives the qubit energy scale, β is the dimensionless signal velocity in the transmission line (units cells per \(\hbar/E_{J}\)), and \(\Gamma_{qb}\), \(\Gamma_{TL}\) are the decoherence rates in a qubit and in the transmission line respectively. Some of the results are shown in Figure 3.

Figure 2
figure 2

Superconducting quantum metamaterials. (a) A 2D superconducting quantum metamaterial [24] (Figure 3 of [24]). (b) A realization of a 1D superconducting quantum metamaterial: charge qubits in a transmission line [31].

Figure 3
figure 3

Electromagnetic wave propagation in a one-dimensional quantum metamaterial (theory). (a) Breathing photonic crystal [18]. (b) Creating a photonic crystal without a direct control of a QMM [35]. (c) Lasing in a QMM [31].

Dimensionless equations of motion for the field vector-potential in the lowest order in field-qubits interaction yield the wave equation

$$ \ddot{\alpha}(\xi,\tau) - \beta^{2} \frac{\partial^{2}\alpha(\xi,\tau )}{\partial\xi^{2}} + V(\xi,\tau) \alpha(\xi,\tau) = 0, $$
(9)

where \(V(\xi,\tau) \propto\langle\hat{\Psi}(\xi,\tau)|\cos\phi|\hat {\Psi}(\xi,\tau)\rangle\). The tunneling matrix element of the charge qubit Hamiltonian (ϕ being the qubit superconducting phase) is determined by the qubit state \(|\hat{\Psi}(\xi,\tau)\rangle\) at the given point. The dispersion law is thus directly dependent on the QMM quantum state, as expected. For example, if the qubit state is a periodic function of the coordinate, \(\hat{\Psi}(x+\Lambda)=\hat{\Psi }(x)\), the QMM behaves as a photonic crystal, with gaps opening in the electromagnetic spectrum [18, 35], which can be manipulated by controlling the quantum state of qubits. For example, if qubits are placed in a spatially periodic superposition of their eigenstates,

$$ \bigl|\hat{\Psi}(x,t)\bigr\rangle = A(x)|g\rangle+ B(x)|e\rangle e^{-i\omega t}, $$
(10)

where \(A(x)\), \(B(x)\) are periodic functions and ħω is the energy splitting between the ground (\(|g\rangle\)) and excited (\(|e\rangle\)) state of the qubit. Quantum beats between these states will produce a ‘breathing’ photonic band structure (Figure 3(a)); external control of qubit states allows to trap a portion of radiation in a pocket of such a structure and move it across the QMM at a desired speed [18]. In the absence of direct control over individual qubit states (apart from their initialization in the ground state) it is still possible to create a photonic crystal structure [35] by sending into the QMM specially shaped ‘priming’ electromagnetic pulses from the opposite directions (Figure 3(b)). Exiting the QMM, they leave behind a spatially periodic pattern of the probability of finding a qubit in the ground (excited) state.

Solving the coupled equations for the classical field and qubits numerically (still in the approximation of factorized qubit state) allows to investigate lasing in a QMM [31]. If the qubits are initialized in the excited state (e.g. by sending a priming pulse through the QMM), an initial pulse triggers a coherent transition of energy from qubits to the electromagnetic field (Figure 3(c)). Remarkably, not only the process has a precipitous character, but its onset starts the sooner the greater the amplitude of the triggering pulse: \(\sqrt{\mathrm{field}\ \mathrm{amplitude}} \times\tau_{\mathrm{onset}} \approx\mathrm{const}\).

In equilibrium a fully quantum treatment of a superconducting QMM becomes possible, which allows the investigation of phase transitions in the photon system [26]. The chosen model of a QMM (a series of RF SQUIDS coupled to the transmission line and considered in two-level approximation) lead to a generic Hamiltonian (1), with only parameters being model dependent. Using the instanton approach, the effective action of the photon subsystem was obtained as a function of the photon field momentum P (in imaginary time).

In case when \(P_{0}\) is independent on the imaginary time τ, the photon system may undergo a second order classical phase transition: above a critical temperature \(T^{\star}\) the momentum \(P=0\), while below it the system can choose between two values \(\pm P_{0}\). The transition details depend on the level of disorder in the qubit chain, which is modeled here by the random distribution of tunneling matrix elements Δ. In the case of a low disorder, \(\Delta_{k} \approx \Delta_{0}\), the transition occurs at the critical temperature

$$ T^{\star}_{n} = \frac{m\eta^{2}N}{k_{B}}, $$
(11)

where m is the transmission line inductance per unit cell, η parametrizes the qubit-field interaction, and N is the total number of qubits in the QMM. This phase transition occurs only if \(\Delta_{0} < k_{B}T^{\star}_{n}\). In the case of strong disorder with \(\Delta_{k}\) distributed from zero to some \(\Delta_{0}\), the transition temperature

$$ T^{\star}_{w} = \frac{\Delta_{0} \exp [-\Delta_{0}/m\eta^{2}N ]}{k_{B}}, $$
(12)

and the transition occurs only if \(\Delta_{0} > m\eta^{2}N\). As the authors of Ref. [26] note, the first case is similar to the metal-ferromagnet phase transition, while the second one is reminiscent of the normal metal-superconductor or Peierls metal-insulator transition. In either case a coherent state of the photon field emerges, with nonzero value of the order parameter \(P_{0}\), which in the case of low disorder is proportional to the number of qubits in the system, N.

It turns out that there is the possibility of a quantum transition as well, into a state representing a superposition of semiclassical states \(\pm P_{0}\). Then the order parameter becomes a periodic function of τ. In the case of low disorder and strong field-qubit coupling it occurs below the temperature (Figure 4)

$$ T^{QOP}_{n} \propto \biggl[\frac{\Delta_{0}^{2}m\eta^{2}N}{\pi^{2}} \biggr]^{1/3} \frac{1}{k_{B}}, $$
(13)

while in the opposite case it is given by

$$ T^{QOP}_{n} = \frac{\Delta_{0}}{k_{B}\pi} \exp \bigl[- \Delta_{0}/m\eta^{2}N \bigr]. $$
(14)
Figure 4
figure 4

Classical and quantum phase transition in the state of photons in a 1D superconducting quantum metamaterial [ 26 ]. The photon order parameter \(P_{0}(\tau)\) as a function of imaginary time τ is either constant (classical phase transition, red line, \(P_{0} \propto N\)) or proportional to \(\sin[P_{0}(T) \tau]\) (quantum phase transition, blue line, \(P \propto N^{1/3}\)). The low disorder case is shown: \(\Delta_{k} \approx\Delta_{0} = 4\)K; \(T^{\star}_{n} = 20\)K.

The estimates of Ref. [26] for the transition temperatures \(T^{\star} \sim0.1, \ldots,50\mbox{ K}\) let us hope for direct observation of photon phase transitions in a superconducting QMM, when the number of qubits, the homogeneity of their parameters and their coherence times are improved as compared to the existing QMM prototype [36, 37]. The prototype (Figure 5(a)) consists of 20 flux qubits placed in a coplanar waveguide resonator. The inductive coupling of qubits to the resonator and each other was comparable, but strong decoherence (decoherence time of the order of a few nanoseconds) effectively suppressed qubit-qubit coupling. Nevertheless the transmission measurements in the resonant regime showed the formation of three ensembles of interacting qubits (two of four qubits each and one of eight qubits). This interaction through the electromagnetic modes is the key element of the operation of a QMM. In order to improve the operation of a QMM it is suggested to increase the qubit-qubit coupling, in order to counteract the effects of qubit parameter dispersion.

Figure 5
figure 5

Quantum metamaterial prototype [ 36 ]. (a) Twenty superconducting flux qubits placed in a resonator. The signal wavelength is of the order of the resonator length (23 mm) and greatly exceeds the qubit size (approx. \(2 \times6\) μm). The ensemble average level spacing and persistent current in the qubits are 5.6 GHz and 74 nA respectively. (b) The QMM in the resonant regime. Three distinct qubit ensembles are seen.

3 Optical quantum metamaterials

The domain of quantum metamaterials in the optical, or near IR, region of the spectrum is still in its infancy. As it has already been stated, some authors use the term quantum metamaterial to denote a structure in which quantum degrees of freedom are inserted [15]. In some other cases, it is the expression ‘quantum dots metamaterials’ that is used: this is to stress that, although quantum dots are inserted in a metamaterial, one is not interested in the quantum coherence of the dots, but rather on the gain that they provide, to counteract the losses due to the presence of metallic inclusions [38]. In other proposals, it is quantum wells that are inserted in a photonic structures. The quantum well are described electromagnetically by a permittivity allowing some control over the behavior of the structure. In Ref. [17], a layered metamaterial is investigated, in which the period comprises two GaAs quantum wells. This structure results in an effective permittivity tensor allowing to obtain a negative refraction. The effective properties strongly depend upon the 2D electron density in the quantum well. In Ref. [15] the same kind of structure is investigated in order to control plasmon propagation, allowing to obtain ultra-long propagation distances.

An original proposal was made in [39] to extend the concept of metamaterial to quantum magnetism. The idea is to use molecular engineering or organic synthesis to fabricate magnetic quantum metamaterials. It is shown theoretically, by ab initio calculations, that CuCoPc2 (a chain of copper-phtalocyanine (CuPc) and cobalt phtalocyanine (CoPc)) possesses a relatively strong ferromagnetic interaction.

In the specific meaning used in this review, a proposal was made in ([21, 28]) to study the full quantum processes that occurs between the quantized electromagnetic field and two-level atoms. The system studied there is a 2D network of coupled atom-optical cavities, called a cavity array metamaterial (CAM) (see Figure 6). The authors propose to realize the model by using a two-dimensional photonic crystal membrane. The quantum oscillators could be quantum dots or substitution centers. Under reasonable assumptions, this system can be described by a Jaynes-Cummings-Hubbard Hamiltonian [40]. The system exhibits a quantum phase transitions [21, 41] and it was proposed that it could be used as a quantum simulator. The effects of cloaking and negative refraction were also demonstrated. Quite naturally, the excitations of the system are hybrids of photonic and atomic states, namely polaritons. This kind of result is very well-known since the pioneering work of Hopfield [42]. Continuing along this way, an effective permittivity that is non-local in time and space could be derived, in exactly the same way as for natural material. This is the line followed in a series of papers by Weick where a collection of metallic nanoparticles is shown to exhibit collective plasmonic modes [4345]. However, a genuine quantum metamaterial requires more than that, namely the active coherent control of the quantum state of the ‘atoms’ inserted in the photonic structure, so as to induce a control over the collective properties of the medium.

Figure 6
figure 6

A cavity array metamaterial, taken from [ 28 ].

Such a system could be implemented by considering, as above, a photonic crystal in which quantum oscillators are inserted, under the guise of quantum dots for instance. The quantum dots can be described semi-classically by a dielectric function \(\varepsilon_{QD}\) that reads as [46]:

$$ \varepsilon_{QD}(\omega)=\varepsilon_{b}+ \bigl(f_{c}(E_{e})-f_{v}(E_{h})\bigr) \frac {a}{\omega^{2}-\omega_{0}^{2}+2i \omega\gamma}. $$
(15)

The prefactor \((f_{c}(E_{e})-f_{v}(E_{h}))\) representing the difference between the populations of the levels can be either positive or negative. In the first case, it is in the absorption regime while in the latter it is in the emission (amplifying) regime. The graph of \(\varepsilon_{QD}(\omega)\) is given in Figure 7. The quantum dots can be grown inside dielectric nanopillars, and the nanopillars can be organized into a 2D periodic array, resulting into a photonic crystal with quantum dots. The bare photonic crystal is then tuned in such a way as to present a photonic band gap at the emission frequency of the quantum dots. The idea is then to realize a pump/probe experiment, where the pump controls the state of the quantum dots (absorption or emission). When the quantum dots are in the emission regime, a transmission peak appears in the transmission spectrum of the probe (Figure 8). This somewhat simplified model shows a macroscopic property, a conduction band, results from the quantum states of the microscopic quantum components. A full quantum treatment should be performed in order to address properly the quantum coherence of the system of the field coupled to the ‘atoms’ of the metamaterial. As compared to the situations encountered in cavity quantum electrodynamics, the quantization of the electromagnetic field in open space comes with severe technical difficulties, if one follows the usual mode-decomposition path, abundantly described in all the literature devoted to field quantization. In fact, for open systems, Maxwell equations do not lead to hermitian eigenvalue problems and thus the eigenfunctions are not normalizable. A recent interesting approach is to use the so-called ‘quasi-normal modes’ of the system [47].

Figure 7
figure 7

The dielectric function of a core-shell quantum dot [ 46 ].

Figure 8
figure 8

Transmission spectrum the quantum metamaterial in the absorption regime (in blue) and in the emission regime (in red). The peak corresponds to the transition of the quantum dots.

4 A review of the theoretical tools for quantum metamaterials in optics

Here we review the tools of quantum optics for the description of the quantum dynamics of a collection of emitters (e.g. quantum dots) in a complex electromagnetic environment that can be constituted by dielectric and/or metallic elements as depicted in Figure 9.

Figure 9
figure 9

A schematic view of a possible quantum metamaterial in optics, made by several quantum dots embedded in semiconducting nanowires surrounded by metallic spheres.

For a single emitter in free space, the interaction with light is given by the minimal-coupling Hamiltonian that reads [48]:

$$H = \frac{1}{2m}(\mathbf{p}-e\mathbf{A})^{2} + eV(r) + \int d^{3}r' \mathbf {E}^{\bot2}+ \mathbf{B}^{2}, $$

where p is the electron momentum, A the vector-potential, \(V(r)\) the binding potential for the electron, B the magnetic field and \(\mathbf{E}^{\bot}\) the transverse part of the electric field, satisfying \(\mathbf{\nabla. E^{\bot}} = 0\) (cf. [49], p.254). The emitters are coupled to each other through the electromagnetic field that comprises both radiating and evanescent terms.

For neutral emitters in a complex electromagnetic environment, one may need a mesoscopic description, where the emitters are described through their multipole moments instead of the dynamics of the electric charges. Moreover, it can be convenient to describe the electromagnetic environment by polarization and magnetization fields instead of charge and current densities [50]. Last but not least, a Hamiltonian involving the physical fields E and B rather than the potentials \(\varphi(\vec{r})\) and \(\mathbf{A}(\vec{r})\) can also be more convenient. Applying the Power-Zienau-Wooley transform to the previous Hamiltonian [48, 50, 51] (see also [52], p.282) fulfills these requirements. The Power-Zienau-Wooley transform leads to the following Hamiltonian (see Milonni [53], p.121):

$$\begin{aligned} H = H_{\mathrm{mat}} + H_{\mathrm{field}} + H_{\mathrm{int}} , \end{aligned}$$
(16)

where

  • \(H_{\mathrm{mat}} = \frac{\mathbf{p}}{2m} + V(r) + \int d \vec{r} \frac {\mathbf{P}^{\bot2}}{2\varepsilon_{0}} \) is the Hamiltonian describing the dynamics of the atom variables, P is the polarization field. Even if the term \(\int d \vec{r} \mathbf{P}^{\bot2}/2\varepsilon _{0}\) is important to reproduce the correct dynamics of the emitter [48], it is usually neglected when studying the interaction between the emitter and light. It is usually argued that this term merely shifts the energy levels, an effect that can be accounted for by a correct renormalization of the emitter energy levels. Nevertheless in the ultra-strong coupling regime, this term has to be taken into account [54]. It leads to a decoupling of matter and light states because of a screening of the incident light by the polarization field P, resulting for example in a reduction of the Purcell factor, while increasing the coupling between the field and the emitter [54].

  • \(H_{\mathrm{field}} = \int \frac{ \mathbf{D}^{2}(\vec{r})}{2\varepsilon _{0}} + \frac{\mathbf{B}^{2}(\vec{r})}{2\mu_{0}}\) is the Hamiltonian describing the electromagnetic field dynamics.

  • \(H_{\mathrm{int}} = -\int \mathbf{P}(\vec{r})\cdot\mathbf{D}(\vec{r}) \,d\vec {r} \) is the Hamiltonian describing the interaction between light and matter.

\(\mathbf{D}(\vec{r})\) is the displacement vector. Note that the Hamiltonian given by the equation Eq. (16) is exact. It is completely equivalent to the minimal-coupling Hamiltonian (Cohen [52], p.298). From Hamiltonian Eq. (16), with the help of the Heisenberg equation, one can find the dynamical equations satisfied by each operators (field and matter operators).

Concerning the field operators, it is assumed that they are related to each other through Maxwell equations. This assumption implies that there exists some commutators between the field operators (p.18 in [55]): \([\mathbf{D}(\vec{r},t),\mathbf{A}(\vec{r'},t)] = \frac{i \hbar}{\varepsilon_{0}} \delta^{T}(\vec{r}-\vec{r'})\) where \(\delta ^{T}(\vec{r})\) is the transverse delta distribution (see also [52], p.233). Finally, one gets the following set of equations between the field operators:

$$\begin{aligned}& \nabla\times \mathbf{E} = - \partial_{t} \mathbf{B}, \\& \nabla\times\mathbf{B} = \mu_{0} \partial_{t} \mathbf{P}+ \frac{1}{c^{2}} \partial^{2}_{t^{2}} \mathbf{E}. \end{aligned}$$

One arrives at the well-known wave equations satisfied by the electric-field operator:

$$ \nabla\times\nabla\times\mathbf{E} - \frac{1}{c^{2}} \frac{\partial^{2} \mathbf{E}}{\partial t^{2}}= - \mu_{0} \frac{\partial^{2} \mathbf {P}}{\partial t^{2}}. $$
(17)

Concerning the dynamics of the matter degrees of freedom, some usual approximations are done. We write the polarization field as a sum of a polarization field due to the atoms \(\mathbf{P}_{a}\) and a polarization field due to the electromagnetic environment \(\mathbf{P}_{\mathrm{ind}}\): \(\mathbf{P} = \mathbf{P}_{a} + \mathbf{P}_{\mathrm{ind}}\). We assume that the polarization field due to the electromagnetic environment \(\mathbf {P}_{\mathrm{ind}}\) responds linearly and locally to the electric field. It reads as [56]

$$\mathbf{P}_{\mathrm{ind}}(\vec{r},t) = \int_{-\infty}^{+\infty} \varepsilon_{0} \chi\bigl( \vec{r},t'\bigr)\mathbf{E}\bigl(\vec{r},t'\bigr) \,dt', $$

where \(\chi(\vec{r},t')\) is the susceptibility at position r⃗ and time t.

Concerning the polarization field due to emitters, we work in the usual dipole approximation and approximate the polarization field by keeping only the first term in the multipolar expansion even if this approximation can be crude for quantum dots [57]. If there are \(N_{\alpha}\) emitters, the polarization field due to the emitters is then written as \(\mathbf{P}_{a}(\vec{r}) = \sum_{\alpha=1}^{N_{\alpha}} \mathbf{d}^{\alpha} \delta(\vec{r} - \vec{r_{\alpha}})\) where \(\mathbf{d}^{\alpha}\) and \(\vec{r_{\alpha}}\) are respectively the dipole-moment operator and the position of the emitter labeled by α, and \(\delta(\vec{r})\) is the usual Dirac distribution. It is convenient to express all matter-operators with the help of the basis defined by the eigenstates of the matter Hamiltonian \(H_{\mathrm{mat}}\). The matter Hamiltonian \(H_{\mathrm{mat}}\) is written as \(H_{\mathrm{mat}} = \sum_{\alpha=1}^{N_{\alpha}} H_{\mathrm{mat}}^{\alpha}\) where \(H_{\mathrm{mat}}^{\alpha}\) is the Hamiltonian of the \(\alpha^{th}\) emitter. We note \(\{|\alpha,i\rangle \}\) the eigenbasis constructed from the eigenstates of \(H_{\mathrm{mat}}^{\alpha}\) that satisfied \(H_{\mathrm{mat}}|\alpha,i\rangle = E^{\alpha}_{i} |i\rangle \) and the completeness condition \(\sum_{i} |\alpha, i\rangle \langle \alpha,i| = I_{d}^{\alpha}\), where \(I_{d}^{\alpha}\) is the identity matrix acting on the subspace of the \(\alpha^{th}\) emitter.

We now assume that emitters are two-level systems. The ground state is labelled by \(i=-\) whereas the excited level is labelled by \(i=+\). Applying twice the completeness condition on the Hamiltonian \(H_{\mathrm{mat}}^{\alpha}\), one finds \(H_{\mathrm{mat}}^{\alpha} = \frac {E_{-}+E_{+}}{2}I_{d}^{\alpha} + \frac{\hbar\omega^{\alpha}}{2} \sigma ^{\alpha}_{z} \) where \(\hbar\omega^{\alpha} = E_{+}^{\alpha} - E_{-}^{\alpha}\) and \(\sigma_{z}^{\alpha} = |\alpha,+\rangle \langle \alpha,+| - |\alpha,-\rangle \langle \alpha,-|\). \(\sigma_{z}^{\alpha}\) acts on the subspace defined by the eigenvectors of the \(\alpha^{th}\) emitter and measured its population difference between the excited and the ground state. The first term in \(H_{\mathrm{mat}}^{\alpha}\) is a constant that can be omitted by choosing correctly the reference of the energy. We then write the matter Hamiltonian as (see [58], p.23 and [53], p.128):

$$H_{\mathrm{mat}} = \sum_{\alpha= 1} ^{N_{\alpha}} \frac{\hbar\omega^{\alpha }}{2} \sigma^{\alpha}_{z}. $$

The polarization field due to the emitters can also be written with the help of the eigenbasis of each emitter by writing the dipole-moment operator in this basis:

$$\mathbf{d}^{\alpha} = \vec{d}^{\alpha}_{+-} |\alpha,+\rangle \langle \alpha ,-| + \vec{d}^{\alpha}_{-+} |\alpha,-\rangle \langle \alpha,+|, $$

where \(\vec{d}^{\alpha}_{+-} = \langle \alpha,+| \mathbf{d}^{\alpha} |\alpha,+\rangle \) is the projection of the dipole-moment operator in the basis \(\{ |\alpha,+\rangle , |\alpha,-\rangle \}\). The diagonal elements are null because we assume that the emitters have no permanent dipole. We introduce the raising \(\sigma^{\alpha}_{+} = |\alpha,+\rangle \langle \alpha ,-|\) and lowering operators \(\sigma^{\alpha}_{-} = |\alpha,-\rangle \langle \alpha,+|\). Finally, the polarization field due to the atoms reads:

$$\mathbf{P}_{a}(\vec{r}) = \sum_{\alpha=1}^{N_{\alpha}} \bigl(\vec{d}^{\alpha }_{+-} \sigma^{\alpha}_{+} + \vec{d}^{\alpha}_{-+} \sigma^{\alpha}_{-} \bigr) \delta(\vec{r} - \vec{r_{\alpha}}). $$

If we assume that the dipole moment projection \(\vec{d}^{\alpha}_{+-} = (\vec{d}^{\alpha}_{-+})^{\star}= \vec{d}_{0}\alpha\) is real, the polarization field simplifies with the help of the Pauli matrix \(\sigma_{x}\):

$$\mathbf{P}_{a}(\vec{r}) = \sum_{\alpha=1}^{N_{\alpha}} \vec{d}^{\alpha }_{0} \sigma^{\alpha}_{x} \delta(\vec{r} - \vec{r_{\alpha}}) $$

since \(\sigma_{x} = \frac{1}{2}(\sigma^{\alpha}_{+} + \sigma^{\alpha }_{-})\) (see [58], p.24, [59], p.35).

With the help of a Fourier transform, the Hamiltonian for the free electromagnetic field can be written \(H_{\mathrm{field}} = \sum_{j,s} \hbar\omega _{j} b_{j,s}^{+}b_{j,s}\) where \(b_{j}\), s (resp. \(b^{+}_{j,s}\)) is the annihilation (resp. creation) operator of the electromagnetic mode j with polarization s. Within this decomposition the electric field on its own reads \(\mathbf{E} = \sum_{j,s} e_{j,s}(b_{j,s} - b^{+}_{j,s})\). Finally, following these successive approximations one finds the Hamiltonian given by the equation Eq. (1)

The behavior of the quantum metamaterial can be computed by solving simultaneously the equation Eq. (17) and the equations of motions for the collection of two-level ‘atoms’ given by [48], p.37 and [59]:

$$\begin{aligned}& \dot{\sigma_{x}^{\alpha}} = -\omega^{\alpha}\sigma_{y}^{\alpha}, \\& \dot{\sigma_{y}^{\alpha}} = \omega^{\alpha}\sigma_{x}^{\alpha} + \frac {2}{\hbar} \vec{d}_{0}^{\alpha}\cdot \mathbf{E}(\vec{r}_{\alpha}) \sigma _{z}^{\alpha}, \\& \dot{\sigma_{z}^{\alpha}} = - \frac{2}{\hbar} \vec{d}_{0}^{\alpha}\cdot\mathbf {E}(\vec{r}_{\alpha}) \sigma_{y}^{\alpha}, \end{aligned}$$

where all terms proportional to \(\mathbf{P}^{2}(\vec{r})\) have been neglected. These equations are non-linear coupled differential equations since \(\sigma_{x}^{\alpha}\) are sources of the electric field Eq. (17). The electromagnetic environment contributes also as a source term in Eq. (17). This source term is important since it is responsible for all the unusual effects demonstrated theoretically or experimentally with classical metamaterials. Its effect on quantum metamaterials has been barely studied since solving this system of equations is challenging, even for a very simple geometry [60]. New theoretical tools should be developed to accurately describe the behavior of a quantum metamaterial in a complex electromagnetic environment. Nevertheless more approximations can be done to solve these equations. One can use the single electromagnetic-mode approximation [21, 28, 61, 62] or perform a semi-classical approximation [18, 19, 30, 31, 63]. In the latter, it is assumed that there are no correlations [59] between the electromagnetic field and the matter degrees of freedom.

5 Conclusions

The field of quantum metamaterials research arose at the intersection of quantum optics, microwave and Josephson physics, and quantum information processing. One of its rather paradoxical feature is that, while the theoretical progress in this area still significantly outweighs the experiment, the theoretical challenges seem more significant. Indeed, the existing experimental techniques, especially in case of superconducting structures, already allow creating massive arrays. The 20-qubits prototype [36] is much smaller than a recently fabricated 1000+-qubits superconducting quantum annealer D-Wave 2X. Given a simpler structure, and less strict demands to a quantum metamaterial than to a quantum computer, making and testing quantum metamaterials on this scale is a question of time and funding. On the other hand, the theoretical analysis of quantum metamaterials produces promising results, already using simple approximations. Nevertheless the understanding of the full scale of effects which can be expected in these systems requires a more detailed analysis of large scale quantum coherences and entanglement. Because of the well-known impossibility to effectively simulate a large quantum system by classical means, a direct approach to this is currently limited to structures containing (optimistically) less than a hundred qubits. New theoretical tools need to be developed, generalizing the methods of quantum theory of solid state [5].

These challenges also present alluring opportunities. Developing and testing new theoretical methods applicable to large quantum coherent systems would be valuable for the whole field of quantum technologies, including quantum computing. Optical elements based on quantum metamaterials would provide new methods for image acquisition and processing. Last but not least, a quantum metamaterial would be a natural test bed for the investigation of quantum-classical transition, which makes this class of structures interesting also from the fundamental point of view.

References

  1. Nakamura Y, Pashkin YA, Tsai JS. Nature. 1999;398:786-8.

    Article  ADS  Google Scholar 

  2. Mooij JE, Orlando TP, Levitov L, Tian L, van der Wal CH, Lloyd S. Science. 1999;285:1036-9.

    Article  Google Scholar 

  3. Friedman JR, Patel V, Chen W, Tolpygo SK, Lukens JE. Nature. 2000;406:43-6.

    Article  ADS  Google Scholar 

  4. Hayashi T, Fujisawa T, Cheong HD, Jeong YH, Hirayama Y. Phys Rev Lett. 2003;91:226804.

    Article  ADS  Google Scholar 

  5. Zagoskin AM, Ilichev E, Grajcar M, Betouras JJ, Nori F. Front Phys. 2014; 2 33. doi:10.3389/fphy.2014.00033.

    Article  Google Scholar 

  6. Albash T, Vinci W, Mishra A, Warburton PA, Lidar DA. Phys Rev A. 2015;91:042314.

    Article  ADS  Google Scholar 

  7. Blais A, Huang R-S, Wallraff A, Girvin SM, Schoelkopf RJ. Phys Rev B. 2004;69:062320.

    Article  Google Scholar 

  8. Blais A, Huang R-S, Wallraff A, Girvin SM, Schoelkopf RJ. Phys Rev A. 2004;69:062320.

    Article  ADS  Google Scholar 

  9. Blais A, Gambetta J, Wallraff A, Schuster DI, Girvin SM, Devoret MH, Schoelkopf RJ. Phys Rev A. 2007;75(3):032329.

    Article  ADS  Google Scholar 

  10. Pendry JB. Phys Rev Lett. 2000;85:3966-9.

    Article  ADS  Google Scholar 

  11. Veselago VG. Sov Phys Usp. 1968;10:509-14.

    Article  ADS  Google Scholar 

  12. Greenleaf A, Kurylev Y, Lassas M, Uhlmann G. SIAM Rev. 2009;51:3-33.

    Article  ADS  MathSciNet  MATH  Google Scholar 

  13. Alitalo P, Tretyakov S. Mater Today. 2013;12:22-9.

    Article  Google Scholar 

  14. Fleury R, Monticone F, Alu A. Rev Phys Appl. 2015;4:037001.

    Article  Google Scholar 

  15. Plumridge J, Phillips C. Phys Rev B. 2007;76:075326.

    Article  ADS  Google Scholar 

  16. Plumridge J, Clarke E, Murray R, Phillips C. Solid State Commun. 2008;146:406-8.

    Article  ADS  Google Scholar 

  17. Plumridge JR, Steed RJ, Phillips CC. Phys Rev B. 2008;77:205428.

    Article  ADS  Google Scholar 

  18. Rakhmanov AL, Zagoskin AM, Savel’ev S, Nori F. Phys Rev B. 2008;77:144507.

    Article  ADS  Google Scholar 

  19. Zagoskin AM, Rakhmanov AL, Savel’ev S, Nori F. Phys Status Solidi B. 2009;246:955-60.

    Article  ADS  Google Scholar 

  20. Zheludev NI. Science. 2010;328:582-3.

    Article  ADS  Google Scholar 

  21. Quach JQ, Su C-H, Martin AM, Greentree AD, Hollenberg LCL. Opt Express. 2011;19:11018-33.

    Article  ADS  Google Scholar 

  22. Felbacq D, Antezza M. SPIE Newsroom. 2012. doi:10.1117/2.1201206.004296.

  23. Zagoskin AM. J Opt. 2012;14:114011.

    Article  ADS  Google Scholar 

  24. Savinov V, Tsiatmas A, Buckingham AR, Fedotov VA, de Groot PAJ, Zheludev NI. Sci Rep. 2012;2:450.

    Article  ADS  Google Scholar 

  25. Zheludev NI, Kivshar YS. Nat Mater. 2012;11:917-24.

    Article  ADS  Google Scholar 

  26. Mukhin SI, Fistul MV. Supercond Sci Technol. 2013;26:084003.

    Article  ADS  Google Scholar 

  27. Zhou L, Gong R, Liu Y-X, Sun CP, Nori F. Phys Rev Lett. 2008;101:100501.

    Article  ADS  Google Scholar 

  28. Quach J, Makin MI, Su C-H, Greentree AD, Hollenberg LCL. Phys Rev A. 2009;80:063838.

    Article  ADS  Google Scholar 

  29. Biondi M, Schmidt S, Blatter G, Tuereci HE. Phys Rev A. 2014;89:025801.

    Article  ADS  Google Scholar 

  30. Zagoskin AM. Quantum engineering: theory and design of quantum coherent structures. Cambridge: Cambridge University Press; 2011.

    Book  Google Scholar 

  31. Asai H, Savel’ev S, Kawabata S, Zagoskin AM. Phys Rev B. 2015;91:134513.

    Article  ADS  Google Scholar 

  32. Blokhintsev DI. Quantum mechanics. Berlin: Springer; 1964.

    Book  MATH  Google Scholar 

  33. Zagoskin AM. Superconducting quantum metamaterials. In: Shadrivov IV, Lapine M, Kivshar YS, editors. Nonlinear, tunable and active metamaterials. Springer series in materials science. vol. 200. Berlin: Springer; 2014. p. 255-79.

    Google Scholar 

  34. Astafiev O, Zagoskin AM, Abdumalikov AA Jr., Pashkin YA, Yamamoto T, Inomata K, Nakamura Y, Tsai JS. Science. 2010;327:840-3. doi:10.1126/science.1181918.

    Article  ADS  Google Scholar 

  35. Shvetsov A, Satanin AM, Nori F, Savel’ev S, Zagoskin AM. Phys Rev B. 2013;87:235410.

    Article  ADS  Google Scholar 

  36. Macha P, Oelsner G, Reiner J-M, Marthaler M, Andre S, Schoen G, Huebner U, Meyer H-G, Il’ichev E, Ustinov AV, Schon G. Nat Commun. 2014;5:5146.

    Article  ADS  Google Scholar 

  37. Ustinov AV. IEEE Trans Terahertz Sci Tech. 2015;5:22-6.

    Google Scholar 

  38. Decker M, Staude I, Shishkin II, Samusev KB, Parkinson P, Sreenivasan VKA, Minovich A, Miroshnichenko AE, Zvyagin A, Jagadish C, Neshev DN, Kivshar YS. Nat Commun. 2013;4:2949.

    Article  ADS  Google Scholar 

  39. Wu W. J Phys Condens Matter. 2014;26:296002.

    Article  Google Scholar 

  40. Henry RA, Quach JQ, Su CH, Greentree AD, Martin AM. Phys Rev A. 2014;90:043639.

    Article  ADS  Google Scholar 

  41. Greentree AD, Tahan C, Cole JH, Hollenberg LCL. Nat Phys. 2006;2:856-61.

    Article  Google Scholar 

  42. Hopfield JJ. Phys Rev. 1958;112:1555-67.

    Article  ADS  MATH  Google Scholar 

  43. Weick G, Woollacott C, Barnes WL, Hess O, Mariani E. Phys Rev Lett. 2013;110:106801.

    Article  ADS  Google Scholar 

  44. Weick G, Mariani E. Eur Phys J B. 2015;88:7.

    Article  ADS  MathSciNet  Google Scholar 

  45. Sturges TJ, Woollacott C, Weick G, Mariani E. 2D Mater. 2015;2:014008.

    Article  Google Scholar 

  46. Holmström P, Thylén L, Bratkovsky A. J Appl Phys. 2010;107:064307.

    Article  ADS  Google Scholar 

  47. Sauvan C, Hugonin JP, Maksymov IS, Lalanne P. Phys Rev Lett. 2013;110:237401.

    Article  ADS  Google Scholar 

  48. Ackerhalt JR, Milonni PW. J Opt Soc Am B. 1984;1(1):116-20.

    Article  ADS  Google Scholar 

  49. Jackson JD. Classical electrodynamics, 3rd ed. New York: Wiley; 1999.

    MATH  Google Scholar 

  50. Jenkins SD, Ruostekoski J. Phys Rev B. 2012;86:085116.

    Article  ADS  Google Scholar 

  51. Power EA, Thirunamachandran T. J Opt Soc Am B. 1985;2:1100-5.

    Article  ADS  Google Scholar 

  52. Cohen-Tannoudji C, Dupont-Roc J, Grynberg G. Atom-photon interactions. New York: Wiley; 2008.

    Google Scholar 

  53. Milonni PW. The quantum vacuum. San Diego: Academic Press; 1994.

    Google Scholar 

  54. De Liberato S. Phys Rev Lett. 2014;112:016401.

    Article  ADS  Google Scholar 

  55. Leonhardt U. Essential quantum optics: from quantum measurements to black holes. Cambridge: Cambridge University Press; 2012.

    Google Scholar 

  56. Vogel W, Welsch D-G. Quantum optics. New York: Wiley; 2006.

    Book  Google Scholar 

  57. Yan J-Y, Zhang W, Duan S, Zhao X-G, Govorov AO. Phys Rev B. 2008;77:165301.

    Article  ADS  Google Scholar 

  58. Carmichael H. An open systems approach to quantum optics. Berlin: Springer; 1993.

    MATH  Google Scholar 

  59. Allen L, Eberly JH. Optical resonance and two-level atoms. New York: Dover; 1975.

    Google Scholar 

  60. Braak D. Phys Rev Lett. 2011;107:100401.

    Article  ADS  Google Scholar 

  61. Quach JQ, Su C-H, Greentree AD. Opt Express. 2013;21:5575-81.

    Article  ADS  Google Scholar 

  62. Everitt MJ, Samson JH, Savel’ev SE, Wilson R, Zagoskin AM, Spiller TP. Phys Rev A. 2014;90:023837.

    Article  ADS  Google Scholar 

  63. McEnery KR, Tame MS, Maier SA, Kim MS. Phys Rev A. 2014;89:013822.

    Article  ADS  Google Scholar 

Download references

Acknowledgements

AZ was supported in part by the EPSRC grant EP/M006581/1 and by the Ministry of Education and Science of the Russian Federation in the framework of Increase Competitiveness Program of NUSTMISiS (No. K2-2014-015).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Alexandre M Zagoskin.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zagoskin, A.M., Felbacq, D. & Rousseau, E. Quantum metamaterials in the microwave and optical ranges. EPJ Quantum Technol. 3, 2 (2016). https://doi.org/10.1140/epjqt/s40507-016-0040-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1140/epjqt/s40507-016-0040-x

Keywords