## Abstract

Quantum theory is expected to govern the electromagnetic properties of a quantum metamaterial, an artificially fabricated medium composed of many quantum objects acting as artificial atoms. Propagation of electromagnetic waves through such a medium is accompanied by excitations of intrinsic quantum transitions within individual meta-atoms and modes corresponding to the interactions between them. Here we demonstrate an experiment in which an array of double-loop type superconducting flux qubits is embedded into a microwave transmission line. We observe that in a broad frequency range the transmission coefficient through the metamaterial periodically depends on externally applied magnetic field. Field-controlled switching of the ground state of the meta-atoms induces a large suppression of the transmission. Moreover, the excitation of meta-atoms in the array leads to a large resonant enhancement of the transmission. We anticipate possible applications of the observed frequency-tunable transparency in superconducting quantum networks.

## Introduction

Quantum metamaterials, media built from quantum objects acting as meta-atoms, promise to display novel light-matter interaction phenomena^{1,2,3,4,5,6,7,8,9}. Nonclassical electrodynamics of such media is linked to many-body quantum mechanics of quantum simulators^{10,11,12}, ensemble quantum memory^{13,14}, generation of non-classical light states^{15,16}, etc. The only reported to date experiments with many superconducting meta-atoms in the quantum regime employed arrays of superconducting qubits weakly coupled to a microwave resonator^{17,18,19}. The observed variations of the transmission coefficient \(\left| {S_{21}} \right|^2\) were, in that case, rather small and limited to a narrow frequency range.

In this work, we aim at achieving tunable electromagnetic properties of the medium in a wide frequency range, not limited by the qubit–resonator interaction. Here we demonstrate intriguing electrodynamic properties of a superconducting quantum metamaterial fabricated from an array of meta-atoms, each consisting of a pair of superconducting loops coupled via a tunnel junction (twin flux qubit). Such meta-atoms provide a strong coupling between qubits and propagating electromagnetic waves. The peculiarity of the twin qubit structure is the field induced 0 → *π* transition of the Josephson junction phases that leads to an abrupt suppression of the microwave transmission in a broad frequency range. In a narrow frequency range, we observe a great enhancement of the microwave transmission. Such resonant transparency is controlled by an external magnetic field. The detailed quantitative analysis is in a good accord with measurements.

## Results

### Physical device and control

The studied system consists of an array of 15 twin flux qubits, each of them containing five Josephson junctions (see Fig. 1). The central Josephson junction (called *α*-junction) of each qubit is shared by two superconducting loops and embedded directly into the central electrode of an Al coplanar waveguide (CPW). This setup provides an extremely large coupling between qubits and electromagnetic waves propagating through the system. We measured the microwave transmission through the CPW as a function of frequency and an external dc magnetic field applied perpendicular to the plane of the substrate. Experiments were performed at low power of the applied microwave field to remain in the linear driving regime. The radiation power is estimated to be on the order of −125 dBm at the chip input port.

### Twin flux qubit model

Precluding the following demonstration of the experimental data and its detailed theoretical analysis, we present the essence of our idea to implement a metamaterial based on twin flux qubit. The twin flux qubit circuit has four nodes. By grounding one of the nodes (setting its nodal phase to zero) and assigning nodal phase variables to the rest of them we can completely describe the state of the system. The Hamilton’s function, which is sum of the charge and Josephson energies, reads

where \(\phi = \left( {\phi _0,\phi _1,\phi _2} \right)\) are the nodal phases and *q* = (*q*_{0}, *q*_{1}, *q*_{2}) are the corresponding canonically conjugate charges (see inset Fig. 1). Here, *Φ* is the external flux applied to a single loop of the twin qubit. Since the Josephson energy is periodic in *Φ*, we analyze the phase distribution across the twin qubit structure at applied fluxes *Φ* from 0 to *Φ*_{0} per single loop of the twin-qubit, where *Φ*_{0} = *h*/2*e* is the magnetic flux quantum. When *Φ* is close to 0, the phase across the central junction *ϕ*_{0} is zero. Simple calculation demonstrates that at *Φ* close to *Φ*_{0}/2 the minimum energy corresponds to *ϕ*_{0} = *π*. This result follows from phase quantization of the superconducting order parameter around the loop. Therefore, one can expect a sharp change of the transmission coefficient \(\left| {S_{21}} \right|^2\) at the critical value of *Φ*_{cr} (between 0 and *Φ*_{0}/2), and a jump back at the external flux \({{\Phi }} \simeq {{\Phi }}_0 - {{\Phi }}_{{{{\rm cr}}}}\). The value of *Φ*_{cr} depends on circuits parameters.

Replacing the nodal phase and charge variables by the respective operators in Hamilton’s function (see Eq. (1)) yields the Hamiltonian operator. By making use of numeric diagonalisation of the Hamiltonian the eigenvalues and the eigenstates are found (see Supplementary Note 1). The magnetic field dependence of the ground state energy of a twin flux qubit is presented in Fig. 2a. For the chosen experimental parameters of the twin flux qubit, in the vicinity of zero-applied flux, the ground-excited state transition energy *hf*_{01} is quite large (see Fig. 2b). However, if the external flux is close to the value of *Φ*_{0}/2 per each loop, the energy *hf*_{01} is about 10 GHz and has an extremum (“sweet spot”). At this point, the twin-qubit transition energy reaches a local maximum rather than minimum characteristic for the standard single-loop three-junction flux qubit^{20}. Moreover, the flux dependence of *f*_{01} is much weaker than that of the usual flux qubit. This makes the twin-qubit more immune to the well known flux qubit problem—strong sensitivity of the transition energy to the external magnetic field fluctuations leading to unavoidable dephasing.

### Measurements of the Transmission Coefficient

The measured dependence of the transmission coefficient \(\left| {S_{21}} \right|\) on both the microwave frequency *f* and the externally applied magnetic flux is shown in Fig. 3a. First, as it is expected, one notices that \(\left| {S_{21}} \right|\) is symmetric around zero field and periodic with respect to the external magnetic flux *Φ* (outside the range of Fig. 3) with period of *Φ*_{0} per single qubit loop. Second, a spectacular result of the measurement is the presence of two regions in the magnetic flux with drastically different transmission coefficients \(\left| {S_{21}} \right|\), see the orange and blue areas in Fig. 3a. The transition between these two regimes is abrupt and occurs at magnetic flux values *Φ*_{cr} ≃ ±0.3*Φ*_{0}, see Fig. 3c. In the range 0.3*Φ*_{0} ≤ *Φ* ≤ 0.7*Φ*_{0}, the transmission amplitude \(\left| {S_{21}} \right|\) changes weakly (from −10 to −5 dB) in frequency band 6–14 GHz outside the resonance (the red feature). The intriguing field-dependent resonant features are characterized by a large enhancement of the transmission amplitude (appearing in Fig. 3a in red) over a relatively narrow frequency range between 11 and 14 GHz. The two symmetric resonant features reach their maximum frequency near *Φ* ~ ±*Φ*_{0}/2. The induced resonant transmission frequency is smoothly decreasing on either side from the maximum and disappears at around 11 GHz.

In order to understand qualitatively the observed modulation of the microwave transparency, let us consider the dependence of the qubit state on the applied magnetic field. At low temperatures, the metamaterial is in its ground state. Therefore, the phases across qubit junctions can be reconstructed as functions of magnetic fluxes applied through qubit loops. Since the external magnetic flux is applied symmetrically to the twin qubits structure, induced current across central junction is zero, and its phase, consequently, is either zero or *π*. Along with increasing the applied magnetic flux, a circulating current through external junctions is induced, providing a monotonic increase of the Josephson phases of the outer Josephson junctions and, in turn, the effective qubit inductance. This causes a smooth decrease of the transparency, see Fig. 3c.

As soon as the external magnetic flux reaches the value of *Φ*_{cr}, it is energetically favorable to increase the energy of the central junction and to lower the energies of the outer junctions. Quantitatively this change is dictated by the generalized Bohr–Sommerfeld quantization rule ensuring a single-valued wave-function of the condensate across qubits loops. As a result at \(\left| {{\Phi }} \right| = {{\Phi }}_{{\mathrm{cr}}}\) the phase shift of the central junction *ϕ*_{0} from zero to *π* is accompanied by the phase shift of the outer junctions by *π*/2, see Supplementary Fig. 4, and the transmission amplitude sharply drops, see Fig. 3c. We call the ground state of the twin qubit in this magnetic field range the *π*-state as opposed to the zero-state at \(\left| {{\Phi }} \right| < {{\Phi }}_{{\mathrm{cr}}}\). In the *π*-state, the currents in the outer branches of the twin flux qubit also change according to the new values of the phases. In particular, at the “sweet” spot *Φ* = *Φ*_{0}/2 phases of all five Josephson junctions of the twin qubit are either zero or *π* which implies zero currents flowing through all junctions, causing the absence of current flowing in all branches of the twin flux qubit.

This nonresonant effect is observed for the entire frequency range used in our experiments and caused by the field-dependent phase distribution in the twin-qubit metamaterial in its ground state. The experimentally obtained value of *Φ*_{cr} is consistent with the calculated dependence on *α*, see Supplementary Note 2. Important to note that widely discussed Electromagnetically Induced Transparency (EIT) has a different nature than the phenomenon described here. “Conventional” EIT is caused by the quantum interference between different scattering channels. It immediately demands the presence of many (at least three) energy levels for joint photon–qubit systems^{21,22,23,24,25,26}. At contrary in the case considered here the classical electromagnetic field interacts with a regular array of qubits, and EIT occurs due to changing of the ground state of a twin flux qubits arrangement.

In the *π*-state region of external magnetic flux, the transition energy between the ground and first excited states of the twin qubit *hf*_{01} can be determined experimentally. The sharp peaks in Fig. 3b and inverted red parabolas between 11 and 14 GHz in Fig. 3a correspond to resonant enhancement of transmission amplitude resulting from the \(\left| 0 \right\rangle \to \left| 1 \right\rangle\) transition. These experimental data in Fig. 3b are fitted by the transition frequency *f*_{01} obtained from diagonalisation of the Hamiltonian (1) and excellent agreement is observed, also see Supplementary Fig. 2.

### Analysis of the transmission coefficient

In spite of the quantum nature of the twin qubit and the strong coupling between qubits and applied microwaves, here we first describe the electromagnetic field as a classical wave. Below we quantitatively analyze the wave propagation through the transmission line which is strongly inductively coupled to a periodic array of quantum oscillators. The transmission line is characterized by coordinate and time dependent charge distribution *Q*(*x*, *t*). Here, *x* is the coordinate along the transmission line. In the presence of a periodic array of twin qubits, which can be considered as superconducting SQUID-like oscillators, the electromagnetic wave propagation is determined by the equation as

where the parameter *β*{*Q*(*x*_{
n
}, *t*)} is the effective scattering strength of a single twin qubit located at *x* = *x*_{
n
}. The parameter *γ*_{tl} describes the intrinsic dissipation of the transmission line (*γ*_{tl} ≪ 1), and *c*_{0} is the velocity of electromagnetic waves in the transmission line.

In the particular case of a twin flux qubit (see Fig. 1) the parameter *β* is obtained from the quantum-mechanical average of the time-dependent Josephson phase of the *α*-junction *φ*_{0}(*t*), i.e.,

where *L* is the transmission line inductance per unit length. The quantum dynamics of the Josephson phase *ϕ*_{0} is determined by the stationary Hamiltonian (1) and the time-dependent perturbation Hamiltonian, \(H_{{\mathrm{per}}} = \left( {E_{\mathrm{J}}{\mathrm{/}}I_{\mathrm{c}}} \right)\dot Q(x,t)\phi _0\).

Let us consider an electromagnetic wave of a given frequency *ω*, *Q*(*x*,*t*) = *Q*(*x*)*e*^{iωt}. In this case, Eq. (2) can be re-written as

Following the method elaborated for the solution of the Schrödinger equation with the Kronig–Penney potential^{27}, we can present the charge distribution *Q*(*x*) in the following form:

where the wave vector \(k = \sqrt {\omega ^2 + i\gamma _{{\mathrm{tl}}}\omega } {\mathrm{/}}c_0\) and *Q*_{
n
} = *Q*(*x*_{
n
}) is the amplitude of propagating charge distribution at point *x*_{
n
}. By making use of the properties of the *δ*-function, we obtain a set of discrete equations for *Q*_{
n
}:

Here, *a* is the effective distance between adjacent point-type scatterers of the array. Taking into account that, in the microwave range, \(ka \ll 1\), we substitute the discrete values of *Q*_{
n
} with a smooth function \(\tilde Q(x)\) satisfying the equation

For low amplitude of the waves, we obtain an explicit expression for the effective scattering strength *β*:

where *Ω* = Δ*E*/*ħ* is the characteristic frequency of the oscillator that can be used to approximate the lowest-lying energy levels of the twin flux qubit. For small external flux values \(\left| {{\Phi }} \right| < {{\Phi }}_{{\mathrm{cr}}}\), where *Φ*_{cr} = *Φ*_{0}/*π* arcsin *α*, the 0-state is the unique stable state of the twin qubit. For the qubits with *α* value used in this paper, *Φ*_{cr} ≈ 0.3*Φ*_{0}. The characteristic frequency *Ω* = *Ω*_{+} in the range \(\left| {{\Phi }} \right| < {{\Phi }}_{{\mathrm{cr}}}\) is relatively large (see Fig. 2b).

Notice here, that the frequency *Ω*_{+} smoothly decreases with the applied magnetic flux *Φ* and results in a parabolic non-resonant decrease of \(\left| {S_{21}} \right|\) as observed in experiments, see the central part of the \(\left| {S_{21}} \right|\) curve in Fig. 3b.

The characteristic frequency of small-amplitude quantum oscillations in the *π*-state for \({{\Phi }}_{{\mathrm{cr}}} < \left| \Phi \right| < {{\Phi }}_0 - {{\Phi }}_{{{{\rm cr}}}}\) can be written explicitly as (see the central part of Fig. 2b).

Since the characteristic frequency \({{\Omega }}_ - \ll {{\Omega }}_ +\) one expects strong changes of the magnitude of the transmission coefficient \(\left| {S_{21}} \right|\) at magnetic flux values *Φ* = *Φ*_{cr} and *Φ* = *Φ*_{0} − *Φ*_{cr}.

Substituting (8) in (7) we obtain an equation for the propagation of electromagnetic waves in this quantum metamaterial:

where the parameter *A* ≈ *ħ*/(2*eLaI*_{c}) depends on the properties of both the transmission line and the Josephson junctions. In our case \(A = 400 \gg 1\) and in this limit we can obtain the transmission coefficient \(S_{21} = \left| {Q(l){\mathrm{/}}Q(0)} \right|\) (*l* is the length of a transmission line) in the following form:

## Discussion

Equation (12) shows a strong dependence of the transmission coefficient \(\left| {S_{21}} \right|\) on the characteristic frequency *Ω*. In particular, there is a non-resonant suppression of \(\left| {S_{21}} \right|\) due to a great enhancement of *K*(*ω*) in the 0-state (see Eq. (8)). However, as *ω* gets close to the *Ω*_{−}, the resonant condition is satisfied, i.e., \(kl\left| {\sqrt {AK(\omega )} } \right| = \pi\), and the qubit array becomes transparent. The typical magnetic field dependence of *S*_{21} for two different frequencies is shown by red lines in Fig. 3b, c. We stress here that in order to observe the both physical effects, the two essential conditions has to be satisfied, namely, a large coupling \(\left( {A \gg 1} \right)\) between qubits and the transmission line, and a large number of qubits in the metamaterial \(\left( {l{\mathrm{/}}a \gg 1} \right)\).

The resonant propagation of electromagnetic field has to be rather sensitive to the disorder in qubits parameters, and therefore, the observation of a single large resonance in the dependence of *S*_{21} on magnetic field, strongly indicates the presence of collective behavior in systems of many qubits.

We have presented the results of the study of the microwave transmission through the quantum metamaterial composed of twin flux qubits. The measured transmission spectra demonstrate that the transparency of the metamaterial is controlled by the external magnetic field in a broad frequency range. Observed large resonant enhancement of the transmission coefficient of the metamaterial in a relatively narrow frequency range corresponds to the qubit \(\left| 0 \right\rangle\) → \(\left| 1 \right\rangle\) transition. Additionally, an abrupt drop of the transmission induced by the external magnetic field and observed in a wide frequency range is caused by the 0 → *π* transition of the *α*-junction, corresponding to a change of the system ground state.

## Methods

### Sample fabrication

The Josephson junctions were fabricated by the shadow-evaporation technique using thin-film Al with a AlO_{
x
} tunnel barrier. The four junctions placed on the outer sides of the superconducting loops have the size of 200 × 355 nm^{2} each, and the central Josephson junction is by a factor *α* = 0.72 smaller. The resulting outer Josephson junctions are estimated to have a Josephon energy of the standard junctions *E*_{J} = 39 × *h* GHz. The ratio of the Josephson and charging energy of *α*-junction is *E*_{J}/*E*_{C} ≈ 3.8.

### Experimental setup

All measurements were carried out at a temperature of about 20 mK in a dilution refrigerator. The frequency dependent microwave transmission coefficient \(\left| {S_{21}} \right|\) was measured with a vector network analyzer.

A superconducting coil around the sample holder was employed to apply a homogeneous dc magnetic field perpendicular to the sample plane. We calibrated the applied dc magnetic field in units of the magnetic flux quantum *Φ*_{0} penetrating the individual qubit loops by making use of similar experiments performed with a single-loop reference SQUID. Ground planes and contact pads on the chip were made out of gold in order to minimize flux focusing and reduce flux noise produced by creeping Abrikosov vortices.

### Data availability

The data that support the findings of this study are available from the corresponding author on request.

## References

- 1.
Rakhmanov, A. L., Zagoskin, A. M., Savel’ev, S. & Nori, F. Quantum metamaterials: electromagnetic waves in a Josephson qubit line.

*Phys. Rev. B***77**, 144507 (2008). - 2.
Zagoskin, A. M., Felbacq, D. & Rousseau, E. Quantum metamaterials in the microwave and optical ranges.

*EPJ Quantum Technol.***3**, 2 (2016). - 3.
Zagoskin, A. M.

*Quantum Engineering: Theory and Design of Quantum Coherent Structures*. pp 272–311 (Cambridge University Press, Cambridge, 2011). - 4.
Jung, P., Ustinov, A. V. & Anlage St., M. Progress in superconducting metamaterials.

*Supercond. Sci. Technol.***27**, 073001 (2014). - 5.
Fink, J. M. et al. Dressed collective qubit states and the tavis-cummings model in circuit QED.

*Phys. Rev. Lett.***103**, 083601 (2009). - 6.
Koppenhofer, M., Marthaler, M. & Schon, G. Superconducting quantum metamaterials as an active lasing medium: effects of disorder.

*Phys. Rev. A***93**, 063808 (2016). - 7.
Ivic, Z., Lazarides, N. & Tsironis, G. P. Qubit lattice coherence induced by electromagnetic pulses in superconducting metamaterials.

*Sci. Rep.***6**, 29374 (2016). - 8.
Blais, A., Huang, R.-S. h., Wallraff, A., Girvin, S. M. & Schoelkopf, R. J. Cavity quantum electrodynamics for superconducting electrical circuits: an architecture for quantum computation.

*Phys. Rev. A***69**, 062320 (2004). - 9.
Monz, T. et al. 14-qubit entanglement: creation and coherence.

*Phys. Rev. Lett.***106**, 130506 (2011). - 10.
Houck, A. A., Tureci, H. E. & Koch, J. On-chip quantum simulation with superconducting circuits.

*Nat. Phys.***8**, 292299 (2012). - 11.
Salathe, Y. et al. Digital quantum simulation of spin models with circuit quantum electrodynamics.

*Phys. Rev. X***5**, 021027 (2015). - 12.
Zakosarenko, V. et al. Realization of a classical counterpart of a scalable design for adiabatic quantum computation.

*Appl. Phys. Lett.***90**, 022501 (2007). - 13.
Saito, S. et al. Towards realizing a quantum memory for a superconducting qubit: storage and retrieval of quantum states.

*Phys. Rev. Lett.***111**, 107008 (2013). - 14.
Wu, H. et al. Storage of multiple coherent microwave excitations in an electron spin ensemble.

*Phys. Rev. Lett.***105**, 140503 (2010). - 15.
Mukhin, S. I. & Fistul, M. V. Generation of non-classical photon states in superconducting quantum metamaterials.

*Supercond. Sci. Technol.***26**, 084003 (2013). - 16.
Hofheinz, M. et al. Generation of Fock states in a superconducting quantum circuit.

*Nature***454**, 310–314 (2008). - 17.
Macha, P. et al. Implementation of a quantum metamaterial using superconducting qubits.

*Nat. Commun.***5**, 5146 (2014). - 18.
Kakuyanagi, K. et al. Observation of collective coupling between an engineered ensemble of macroscopic artificial atoms and a superconducting resonator.

*Phys. Rev. Lett.***117**, 210503 (2016). - 19.
Shulga, K. V. et al. Observation of a collective mode of an array of transmon qubits.

*JETP Lett.***105**, 38 (2017). - 20.
Orlando, T. P. et al. Superconducting persistent-current qubit.

*Phys. Rev. B***60**, 15398 (1999). - 21.
Abdumalikov, A. A. et al. Electromagnetically induced transparency on a single artificial atom.

*Phys. Rev. Lett.***104**, 193601 (2010). - 22.
Fleischhauer, M., Imamoglu, A. & Marangos, J. P. Electromagnetically induced transparency: optics in coherent media.

*Rev. Mod. Phys.***77**, 633 (2005). - 23.
Hoi, I.-C. et al. Giant crosskerr effect for propagating microwaves induced by an artificial atom.

*Phys. Rev. Lett.***111**, 053601 (2013). - 24.
Sun, H.-C. et al. Electromagnetically induced transparency and autler-townes splitting in superconducting flux quantum circuits.

*Phys. Rev. A***89**, 063822 (2014). - 25.
Inomata, K. et al. Single microwavephoton detector using an artificial

*λ*-type threelevel system.*Nat. Commun.***7**, 12303 (2016). - 26.
Long J. et al. Electromagnetically induced transparency in circuit QED with nested polariton states. Preprint at https://arxiv.org/abs/1704.08777 (2017).

- 27.
Lifshits, I. M., Gredeskul, S. A. & Pastur, L. A.

*Introduction to the Theory of Disordered Systems*. (Wiley-VCH, New York, 1988).

## Acknowledgements

This work was partially supported by the Ministry of Education and Science of Russian Federation in the framework of Increase Competitiveness Program of the NUST MISIS, contracts no. K2-2016-063 (experiment) and K2-2016-067 (theory). Experiments and analytical tools were developed with the financial support from the Russian Science Foundation (contract no. 16-12-00095). M.V.F. acknowledges hospitality of the International Institute of Physics, Natal, Brazil where this work has been finished.

## Author information

### Affiliations

### Contributions

K.V.S. designed and performed the experiments, analyzed the data and wrote the paper; E.I. created the analytic model with contributions from O.V.A. and A.V.U., supervised its analysis and edited the manuscript; I.S.B. developed analytical tools and wrote the Supplementary Information; S.B. designed the experimental samples; M.V.F. designed the theoretical model; U.H. provided the experimental samples.

### Corresponding author

## Ethics declarations

### Competing interests

The authors declare no competing financial interests.

## Additional information

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Electronic Supplementary Material

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Shulga, K.V., Il’ichev, E., Fistul, M.V. *et al.* Magnetically induced transparency of a quantum metamaterial composed of twin flux qubits.
*Nat Commun* **9, **150 (2018). https://doi.org/10.1038/s41467-017-02608-8

Received:

Accepted:

Published:

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.