Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.cplire.ru/html/lab234/pubs/2014_02.pdf Äàòà èçìåíåíèÿ: Mon Mar 21 06:11:21 2016 Äàòà èíäåêñèðîâàíèÿ: Sun Apr 10 06:40:32 2016 Êîäèðîâêà: ISO8859-5

AR TICLE
Received 13 Dec 2013 | Accepted 26 Mar 2014 | Published 28 Apr 2014
DOI: 10.1038/ncomms4730

Multistability and switching in a superconducting metamaterial
P. Jung1, S. Butz1, M. Marthaler2, M.V. Fistul
3,4,5,

? J. Leppakangas6, V.P. Koshelets4,7 & A.V. Ustinov

1,4

The field of metamaterial research revolves around the idea of creating artificial media that interact with light in a way unknown from naturally occurring materials. This is commonly achieved using sub-wavelength lattices of electronic or plasmonic structures, so-called metaatoms. One of the ultimate goals for these tailored media is the ability to control their properties in situ. Here we show that superconducting quantum interference devices can be used as fast, switchable meta-atoms. We find that their intrinsic nonlinearity leads to simultaneously stable dynamic states, each of which is associated with a different value and sign of the magnetic susceptibility in the microwave domain. Moreover, we demonstrate that it is possible to switch between these states by applying nanosecond-long pulses in addition to the microwave-probe signal. Apart from potential applications for this all-optical metamaterial switch, the results suggest that multistability can also be utilized in other types of nonlinear meta-atoms.

1 Physikalisches Institut, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany. 2 Institut fur Theoretische Festkorperphysik, Karlsruhe Institute of ? ? Technology, 76131 Karlsruhe, Germany. 3 Theoretische Physik III, Ruhr-Universitat Bochum, 44801 Bochum, Germany. 4 Laboratory for Superconducting ? Metamaterials, National University of Science and Technology MISIS, Moscow 119049, Russia. 5 Theoretical Physics and Quantum Technologies Department, National University of Science and Technology MISIS, Moscow 119049, Russia. 6 Microtechnology and Nanoscience, MC2, Chalmers University of Technology, Goteborg 412 96, Sweden. 7 Laboratory of Superconducting Devices for Signal Detection and Processing, Kotel'nikov Institute of Radio ? Engineering and Electronics, Moscow 125009, Russia. Correspondence and requests for materials should be addressed to P.J. (email: philipp.jung@kit.edu).

NATURE COMMUNICATIONS | 5:3730 | DOI: 10.1038/ncomms4730 | www.nature.com/naturecommunications

1

& 2014 Macmillan Publishers Limited. All rights reserved.


AR TICLE
n recent years, the idea of using superconducting quantum interference devices (SQUIDs) and their quantum analogues (qubits) as meta-atoms in superconducting microwave metamaterials has become increasingly appealing as the attention of the community shifts towards tunability and parametric effects1-4. In the quantum limit, efforts to build many-qubit systems blur the border between circuit quantum electrodynamics and metamaterials5-8 while even the classical case holds a multitude of interesting and yet unexplored possibilities. Implementations of these structures, for instance in the form of parametric amplifiers9, have already attracted a great deal of attention and continue to do so as the field of superconducting metamaterials evolves. In the present work, we are focusing on an effect found in single-junction (radio frequency (rf-))SQUIDs. To some extent, they are similar to split-ring resonators, one of the most common meta-atom designs, as they consist of a (superconducting) loop interrupted by a Josephson tunnelling barrier instead of a capacitor. Although superconductivity and the Josephson effect are quantum phenomena, the way the SQUIDs are built and operated allows us to describe them using classical dynamic equations, as one would model other resonant metaatoms. Unlike other resonant meta-atoms, however, they exhibit a unique dependence of their resonance frequency on a constant, external magnetic field. It stems from the nonlinear nature of the Josephson junction. This effect has been widely discussed10-12 and demonstrated in different waveguide geometries13-15. Although this truly remarkable feature is a direct effect of the Josephson nonlinearity, the SQUIDs in the cited experiments were mostly operated at excitation powers low enough to treat the junction as a quasi-linear inductive element. With increasing power, this simple model fails as the nonlinearity starts to influence the resonant behaviour. So far, this regime has not been explored in great detail experimentally. The essence of this work is the dynamic multistability which occurs in non-hysteretic rf-SQUIDs at intermediate power levels where the linear treatment starts to fail. This has already been observed in single junctions16 and investigated numerically for rfSQUID meta-atoms17. It should be noted that here, we deal with a purely dynamic phenomenon that is not related to the multistability known from hysteretic SQUIDs18,19. For a certain choice of parameters, the dynamic multistability manifests itself as a small number of simultaneously stable states, each of which corresponds to a different value of the SQUIDs' magnetic-flux susceptibility wf at the driving frequency. More importantly, a
Arg (S21) Power (a.u.) (rad)

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4730

I

closer analysis reveals that in many cases coexisting stable states have small-positive as well as large-positive and large-negative values of the real part of wf. This means that, depending on the current state of the meta-atom, it can either be magnetically almost transparent, strongly enhance the local-magnetic field by oscillating in phase with the drive, or even counteract the same by oscillating out of phase. This is a highly desirable feature for a meta-atom provided one can reliably switch it from one to another state in situ. Similar to the quantum switch20 and unlike most other switchable meta-atoms, the process of switching in this approach is not associated with a change of the static mechanical21 or electronic22 properties. Here, we demonstrate the switchability of the SQUIDs experimentally and provide an analytical model describing the observed behaviour. Finally, we show how a chain of these metaatoms can be implemented as a metamaterial device. Results Switching by microwave pulses. The most prominent feature associated with the effect described below is that it is possible to switch between the states reproducibly by applying a shortmicrowave pulse to the metamaterial. To verify this experimentally, we used a coplanar waveguide (CPW) loaded with a one-dimensional (1D) metamaterial consisting of 54 rf-SQUIDs. A row of 27 SQUIDs is located in each of the two gaps of the CPW. The coupling between the waveguide and the SQUIDs is designed to be primarily magnetic, which is why a change in the magnetic susceptibility seen by the waveguide can be measured as a change in its scattering characteristics. Figure 1 shows an exemplary transmission measurement recorded in such an experiment. Initially, the system is in a state with a susceptibility close to zero which results in little scattering and thus a large transmission through the loaded CPW. By applying the microwave pulse, the system is driven into a state with |wf|c1 which results in both enhanced reflection and absorption of the probe signal. Together, this causes the amount of transmitted power to drop by one order of magnitude (10 dB) for this sample. The system remains stable in the new state until it is reset. One noteworthy difference when comparing this system with other switchable metamaterials is that here, the switching process is rather fast because it only relies on the intrinsic properties of the nonlinear resonator. Thus the time required to switch is proportional to the decay time td Ì 2RC which is of the order of a

a

d
Via JJ

0.1 0.05 0 0 -5 -10

b
Capacitor

|S21| (dB)

c
100 m 0 0.1 0.2 0.3 0.4 0.5 0.6 Time t (s) 0.7 0.8 0.9 1

Figure 1 | Switching of a 54 SQUID array. (a) Sketch of the pulse sequence used to switch the 1D metamaterial. The actual measurement, indicated by two vertical dashed lines, takes place after initializing the system in a small amplitude state. (b) Transmission phase and (c) magnitude measurement of the SQUID-loaded waveguide recorded as a function of time. The steepness of the jump is limited by measurement bandwidth. The base-line has been subtracted for both curves. The pulse switches the system into another stable, dynamic state which is visible as change in transmission. (d) Microscope image of a section of the loaded CPW and a close up sketch of a single SQUID. The arrow in the round inset labelled JJ marks the position of the Josephson junction.
2
NATURE COMMUNICATIONS | 5:3730 | DOI: 10.1038/ncomms4730 | www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.


NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4730

AR TICLE
Here, g is the potential in the rotating frame23:
! pffiffiffiffiffiffiffiffiffiffiffiffiffiffi Ñ 1 o2 Ð 2 gÌ x2 ? y2 ; 1 Ð bL 2 x ? y2 Ð jea x Ð 2bL J0 op 2

few nanoseconds for our samples. Here, R and C are the intrinsic resistance and capacitance of the SQUID meta-atom, respectively. Experimentally we were able to verify that the switching takes place for pulses as short as 2 ns. Analytical modelling. To better understand and predict the existence of the multiple stable, dynamic states, we developed an analytical model that is able to quantitatively reproduce most of the experimental observations. We start from the basic differential equation that defines the evolution of the total gaugeinvariant phase difference j over the Josephson junction of the SQUID when driven externally by jext. " # 1 1 _ ?1î j ? bL sinj ? j ? 2 j Ì jext ; oc op
2c h where bL Ì 2pLIc , oc Ì 2pRIc , o2 Ì CpI0 , F0 Ì 2e, Ic is the maxp F0 F0 F imum supercurrent that can be carried by the Josephson junction, C is the capacitance shunting the junction, L is the inductance of the loop, and R is the damping. Assuming a harmonic drive jext Ì jea cosot, we can make an ansatz for the solution of the equation of motion: j Ì x cos(ot) Ð y sin(ot). This corresponds pffiffiffiffiffiffiffiffiffiffiffiffiffiffi to all harmonic solutions that have an amplitude A Ì x2 ? y2 and a phase d Ì arctan(y/x) with respect to the drive. The magnetic-flux susceptibility can then be defined as the time average of the ratio between the flux in the loop and the driving flux which can be expressed as wj Ì /j/jextSt Ð 1. It should be noted at this point, that it is not obvious that this approach will describe the dynamics of the system sufficiently well, because it does not contain contributions from higher harmonics. A detailed treatment of why this approximation is applicable in our case is provided in the methods section. Using the mentioned ansatz, one can write down the equations of motion in the rotating-wave approximation (RWA) together with the stability condition:

? 4î

bL

o @g o! _ Ð bL y Ì 0 yÌ Ð o2 @x oc p o @g o! _ xÌ Ð bL x Ì 0 o2 @y oc p

?2î

where J0 is the zeroth order Bessel function of the first kind. These equations yield both stable and unstable solutions which can be distinguished from each other through a stability analysis. Figure 2a,b shows the real and imaginary part of the calculated magnetic-flux susceptibility for a fixed frequency as a function of the amplitude of the drive. For a certain range of driving power, multiple stable (black) and unstable (grey) solutions coexist. Moreover, the simultaneously stable states can have a smallpositive as well as a large-negative and a large-positive values of the real part of wf. These results can then be used in a transmission-line approach to calculate the transmitted power through the sample as seen in Fig. 2c (black line). Qualitatively, this effect exists for a wide range of frequencies and magnetic field values, however with varying magnitude. Here, we focus on a region slightly above the LC-resonance frequency fLC Ì Ð pffiffiffiffiffiffiÑ Ð 1 where the effect is most pronounced. 2p LC The fact that multistability only occurs in a limited range of driving power24 can be understood from the particular form of the Josephson nonlinearity sin j in equation (1). For very small-driving amplitudes, this term can be linearized around a static phase value jconst which is in turn related to the direct current (DC) magnetic flux seen by the SQUID loop. Therefore, the SQUID is usually treated as a DC-flux tunable, quasi-linear resonator in this driving regime. For very large-driving amplitudes the nonlinear term in equation (1) is strongly suppressed. As a result, the system can be described as a DCflux independent LC circuit resonating at fLC in this regime. Consequently, the only case in which the nonlinearity contributes to the dynamic behaviour of the system is in the intermediate range of driving amplitudes. Our analysis also shows that both the number of observable states and their amplitudes strongly depend on the difference between driving frequency and fLC. Multistability in a single SQUID. To test the analytical model and verify the existence of the dynamic states, we demonstrate two ways of measuring them: by sweeping power up and down thus going through a hysteresis loop and by exciting the system

b

L

?3î

a
Real (
)

100 50 0 -50

c
-1.4 -1.5 |S21| (dB) -1.6 -1.7 -1.8 -1.9 -65

-100 10-3

10-2

10

-1

10

0

b
-Imag (
)

200 150 100 50 0 10-3 10-2 10-1 Driving-flux amplitude/0 10
0

-60

-55

-50 -45 -40 Input power (dBm)

-35

-30

-25

Figure 2 | Multistability of a single SQUID. (a) Real and (b) imaginary part of the magnetic-flux susceptibility of a single SQUID meta-atom calculated from equations (2)-(4). Stable solutions are shown in black, while unstable solutions are in grey. (c) Calculated (black) and measured (blue and red) transmission through the sample arrangement containing only a single SQUID. The red data show a hysteresis loop from low to high power and back. The red arrows indicate the direction of the sweep. The blue data are a collection of different power sweeps of varying length and initial conditions. The overall transmission of less than Ð 1.4 dB is caused by the experimental setup and calibration and not related to the performance of the meta-atom.
NATURE COMMUNICATIONS | 5:3730 | DOI: 10.1038/ncomms4730 | www.nature.com/naturecommunications

3

& 2014 Macmillan Publishers Limited. All rights reserved.


AR TICLE
using microwave pulses. We first demonstrate the hysteresis (red curve in Fig. 2c) for a single SQUID meta-atom to obtain results that can be compared directly with the single SQUID analytic solution. At a low initial-power level, the systems starts out in a transparent state (small absolute value of wf) resulting in a large value of the transmission |S21| through the SQUID-loaded CPW. When increasing the power, it remains in this state until it is no longer stable (E Ð 37 dBm in this case). The system then settles in one of the excited states (large-absolute value of wf), which results in a drop in the transmitted power. Following a series of small jumps, the system eventually reaches a stable level. Upon sweeping the power back, the observed transmission takes a different path by following several of the highly excited states until they become unstable then jumping to the closest stable one. Although being consistent with the model, this method only reveals a small fraction of the predicted-stable states. For a more complete picture one can prepare the system in each of the states using the hysteresis and then follow each path in both directions until it becomes unstable. The resulting curves are shown as blue lines in Fig. 2c. One noteworthy aspect in this figure is that all but one of the highly excited states correspond to solutions with a negative-real part of wf as seen in Fig. 2a. The only solution with a large-positive real part of wf in Fig. 2a corresponds to the line that crosses the base-line at about Ð 50 dBm driving power in Fig. 2c. Although the data agree with the theoretical prediction to a quantitative level in many aspects, there is still a slight discrepancy between the two. The most noticeable difference is that even by varying the initial conditions for each sweep, we are not able to prepare the system in all states predicted by the model. Multistability in SQUID arrays. As mentioned above, this approach can be extended from a single meta-atom to a metamaterial. In this case, however, even without any disorder, the multistability will not have such a clear signature in transmission measurements any more. This is owing to the fact that the probe signal gets absorbed and reflected as it passes through the 1D metamaterial so that every meta-atom sees a different effective input power. In case of the 1D 54 SQUID array, one can still observe the characteristic hysteresis loop (blue curve in Fig. 3) albeit without the distinct response of the different states. Figure 1 shows that in this case, it is possible to decrease the transmission by up to 10 dB by applying a short-microwave pulse. With the knowledge of the
1 Transmission magnitude (dB) 0 -2 -4 -6 -8 -10 -60 -55 -50 Input power (dBm) -45 -40

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4730

theory and the hysteresis curve, the pulse sequence shown in Fig. 1a becomes clear: the first part is necessary to reset the system to the transparent state. The intermediate power level (usually just referred to as input or driving power) moves the system to the multi-stable power region, while keeping it in the transparent branch. Finally, the pulse excites the system into a different branch of the hysteresis. The decrease in transmission thus depends on the input-power level and is subject to a statistical spread. A more detailed overview of the switching as a function of the input power is given in Fig. 3. The black points show the mean magnitude of the difference in |S21| before and after applying the microwave pulse and its standard deviation. Each of the points contain information from 150 measurements using the sequence outlined in Fig. 1. Apart from the driving-power level, the switching also depends on duration and amplitude of the pulse, all of which influence the number of switched SQUIDs in the array and thus the resulting change in transmission. Ideally, the magnitude of the switching should follow the difference between the up- and down-sweep branches of the hysteresis. For the combination of parameters chosen in Fig. 3 (pulse duration: 200 ns, pulse power at the sample: Ð 38 dBm), however, the switching only works reliably in part of the input-power range (between Ð 58 dBm and Ð 52 dBm). Discussion The presented results demonstrate the general possibility of creating an on-chip metamaterial switch from switchable metaatoms. As such, its performance (on/off ratio of about 10:1 and limited dynamic range) can be questioned. With the current design of identical SQUIDs there is indeed a limit to those quantities. The performance of the device is related to the hysteresis curve of the single SQUID. This relation can qualitatively be summarized as follows: when all SQUIDs are in the transparent state, they all see the same amount of driving power since the scattering at every meta-atom is small. However, as soon as some of the SQUIDs are switched to a different state, they absorb and reflect part of the microwave drive, so the ones located behind them (in propagation direction) feel a weaker drive. In a long chain of SQUIDs a situation will eventually be reached in which the elements further back in the chain are biased outside of their hysteresis range. Thus the limit of the switching magnitude of the array is given by the width of the hysteresis curve of a single SQUID as seen in Fig. 2. Some options to improve upon the basic design include optimizing the individual SQUIDs (in terms of width and onset of the hysteresis) and making them non-identical. Despite the aforementioned restrictions we believe that this approach presents a viable alternative to other switching techniques for the following reasons: first, the device can be switched by comparatively short, low-power pulses and remains stable with a high fidelity. Second, it can easily be integrated with other technology based on Josephson junctions such as qubits or Josephson parametric amplifiers. Third, the concept can directly be extended into two or three dimensions by arranging SQUIDs accordingly thus making a pulse-switchable surface or 3D body. We would also like to point out that, although the nature of the Josephson junction is responsible for this particular shape of the multistability, the way it is implemented here is of a more general nature. Switching the oscillatory state of a meta-atom by means of pulses can be extended to other types of nonlinearites and different technologies that do not require superconductivity25. In conclusion, we have implemented and demonstrated a 1D metamaterial consisting of switchable meta-atoms. By operating non-hysteretic, single-junction SQUIDs in an intermediate power range, we showed that they can be switched between different

Figure 3 | SQUID-array results. Hysteresis curve for a CPW loaded with 54 SQUIDs (blue line) normalized to the value at Ð 65 dBm. The black dots show the mean magnitude and standard deviation (150 samples) of the difference in |S21| before and after applying a microwave pulse. Up to about Ð 58 dBm this quantity is well-described by the difference between the lower and upper branch of the hysteresis. Above that threshold, the switching magnitude drops to a lower value.
4

NATURE COMMUNICATIONS | 5:3730 | DOI: 10.1038/ncomms4730 | www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.


NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4730

AR TICLE
equation of motion: 0Ì Ð @g 3o Ð bL y3 : @ x3 op ?11î

dynamic states corresponding to different values of their magnetic susceptibility. For individual meta-atoms, we were able to map these states directly in a transmission experiment. We developed an analytical model based on the rotating-wave approximation which shows excellent agreement with the measured data. Using an array containing 54 SQUID metaatoms, we showed that one can switch the transmission through this 1D metamaterial by applying a nanosecond-long microwave pulse. Unlike in other switchable meta-atom implementations, we employ an intrinsic nonlinearity which makes these SQUIDs very fast-switching meta-atoms. Methods
The samples. The SQUIDs are thin-film Nb structures with a Nb-AlOx-Nb tunnel barrier. They are B45 Ò 49 mm2 in size with a 2 mm2 junction area which is shunted by a large Nb2O5 capacitor. They are positioned in the gaps of a CPW on a 4 Ò 4 mm Si chip. The values used to model the SQUIDs are op Ì 5.276 Ñ 1010s Ð 1, oc Ì 5.438 Ñ 1012s Ð 1, bL Ì 0.45. Experimental setup. All measurements were performed in liquid 4He at a temperature of 4.2 K using cryogenic amplification. The sweeps were recorded with a vector network analyzer. Pulses were generated using a microwave mixer and a pulse generator. Data calibration. The measurements were thru-calibrated at the connector closest to the sample. Additional scatterers that could not be calibrated out this way (that is, bond-connections) were accounted for whenever theory and experiment were compared. Validity of the RWA. In general, the solution for the phase j is of the form, j Ì x0 ?
1 XÒ Ó xn cos?not î? yn sin?not î : nÌ1

0Ì

@g 3o Ð bL x3 : @ y3 op

?12î

Solving this, we get for the amplitude of the oscillations of the third mode pffiffiffiffiffiffiffiffiffiffiffiffiffiffi qffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2 2 b j J3 ? x1 ? y1 î j 2 2 ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi : A3 Ì x3 ? y3 Ì rL 1 Ð 9b
o2 L o2 p 2

?13î

? 9b

2o Lo

2 2 c

On resonance we have b

Ì 1, thus we can simplify the expression to qffiffiffiffiffiffiffiffiffiffiffiffiffiffi b b 2 2 A3 % L J3 x1 ? y1 Ì L jJ3 ?A1 îj: 8 8

o2 L o2 p

?14î

The ratio between the amplitude of the third harmonic to the amplitude at the driving frequency A3/A1 is a useful figure to quantify the validity of the RWA. We can find an upper limit for this quantity since the Bessel function divided by its argument can be estimated by |J3(A1)|/A1o1/9 and bLo1/2. Thus A3/A1o1/144 which means that the rotating-wave approximation is well justified. The frequency range of interest for the experiment at hand is determined by the low-power qffiffiffiffiffiffiffiffiffi range of the resonance frequency ffi tuning qffiffiffiffiffiffiffiffiwhich is bounded by o
max

Ì op

frequency range the ratio is A3/A1o1/63, which is still well within the regime of the rotating wave approximation.

1?b bL

L

from above and omin Ì o

p

1Ðb bL

L

from below. In this

References
1. Zheludev, N. I. & Kivshar, Y. S. From metamaterials to metadevices. Nat. Mater. 11, 917-924 (2012). 2. Anlage, S. M. The physics and applications of superconducting metamaterials. J. Opt. 13, 024001 (2010). 3. Zheludev, N. I. A roadmap for metamaterials. Optics Photon. News 22, 30-35 (2011). 4. Boardman, A. D. et al. Active and tunable metamaterials. Laser Photon. Rev. 5, 287-307 (2011). ` 5. Lalumiere, K. et al. Input-output theory for waveguide QED with an ensemble of inhomogeneous atoms. Phys. Rev. A 88, 043806 (2013). 6. van Loo, A. F. et al. Photon-mediated interactions between distant artificial atoms. Science 342, 1494-1496 (2013). 7. Viehmann, O., von Delft, J. & Marquardt, F. Observing the nonequilibrium dynamics of the quantum transverse-field ising chain in circuit QED. Phys. Rev. Lett. 110, 030601 (2013). 8. Macha, P. et al. Implementation of a Quantum Metamaterial. Preprint at http:// arxiv.org/abs/1309.5268 (2013). 9. Castellanos-Beltran, M. A., Irwin, K. D., Hilton, G. C., Vale, L. R. & Lehnert, K. W. Amplification and squeezing of quantum noise with a tunable Josephson metamaterial. Nat. Phys. 4, 928-931 (2008). 10. Du, C., Chen, H. & Li, S. Quantum left-handed metamaterial from superconducting quantum-interference devices. Phys. Rev. B 74, 113105 (2006). 11. Lazarides, N. & Tsironis, G. P. RF superconducting quantum interference device metamaterials. Appl. Phys. Lett. 90, 163501 (2007). 12. Maimistov, A. I. & Gabitov, I. R. Nonlinear response of a thin metamaterial film containing Josephson junctions. Opt. Commun. 238, 1633-1639 (2010). 13. Jung, P., Butz, S., Shitov, S. V. & Ustinov, A. V. Low-loss tunable metamaterials using superconducting circuits with Josephson junctions. Appl. Phys. Lett. 102, 062601 (2013). 14. Butz, S., Jung, P., Filippenko, L. V., Koshelets, V. P. & Ustinov, A. V. A onedimensional tunable magnetic metamaterial. Opt. Express 21, 22540-22548 (2013). 15. Trepanier, M., Zhang, D., Mukhanov, O. & Anlage, S. M. Realization and modelling of metamaterials made of rf superconducting quantum-interference devices. Phys. Rev. X 3, 041029 (2013). 16. Siddiqi, I. et al. Direct observation of dynamical bifurcation between two driven oscillation states of a Josephson junction. Phys. Rev. Lett. 94, 027005 (2005). 17. Lazarides, N. & Tsironis, G. P. Multistability and self-organization in disordered SQUID metamaterials. Supercond. Sci. Technol. 26, 084006 (2013). 18. Du, C., Chen, H. & Li, S. Stable and bistable SQUID metamaterials. J. Phys. Condens. Matter 20, 345220 (2008). 19. Caputo, J.-G., Gabitov, I. & Maimistov, A. I. Electrodynamics of a split-ring Josephson resonator in a microwave line. Phys. Rev. B 85, 205446 (2012). 20. Marlantoni, M. et al. Two-resonator circuit quantum electrodynamics: a superconducting quantum switch. Phys. Rev. B 78, 104508 (2008). 21. Liu, A. Q., Zhu, W. M., Tsai, D. P. & Zheludev, N. I. Micromachined tunable metamaterials: a review. J. Opt. 14, 114009 (2012).
5

?5î

In equation (1), we only keep the resonant components y1 Ì y and x1 Ì x according to the rotating-wave approximation. Since the anharmonicity is relatively large, we also want to explicitly control the validity of this approximation. If we go beyond the rotating-wave approximation not only the first mode, oscillating with frequency o, will be excited but also higher modes. The closer the higher the modes are to the first mode, the stronger they should be excited. Therefore, we keep the second and third mode oscillating at 2o and 3o, respectively. In this case, the quasi-energy Hamiltonian takes the form
gÌ ! qffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 o2 2 2 2 2 1 Ð bL 2 ?x1 ? y1 îÐ jea x1 Ð 2bL J0 x1 ? y1 cos?x0 î 2 op ! " X1 n2 o2 2 2 1 Ð b L 2 ? x n ? yn î ? 2 op nÌ2;3 qffiffiffiffiffiffiffiffiffiffiffiffiffiffi ! 12 2 2 x1 ? y1 sin?x0 ? pn=2î?xn cos?ndîÐ yn sin?ndîî ? y0 ? x Ð bL Jn bL

?6î
2 0

where y0 is the conjugate variable of x0, Jn the nth order Bessel functions of the first kind and we assumed that xn and yn remain small for n Ì 2, 3. The resulting equations of motion for nZ1, are given by b
L

no @g no _ Ð bL yn : yn Ì Ð o2 @ xn oc p no @g no _ Ð bL xn : xn Ì o2 @ yn oc p

?7î

b And for n Ì 0 we get

L

?8î

1 @g op _ Ð y0 : y0 Ì Ð oc op @ x0 1 @g _ : x0 Ì op @ y0

?9î

?10î

_ _ Equation (9) has a stationary state (y0 Ì x0 Ì 0) for x0 Ì 0 independent of x1 and y1. In this case we only have to consider the coupling between the resonant mode and x3, y3 because the first mode will not couple to the second mode. Using the lowest order of perturbation theory we calculate the solution for x1 and y1 as discussed above and find an order of magnitude estimate from the stationary

NATURE COMMUNICATIONS | 5:3730 | DOI: 10.1038/ncomms4730 | www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.


AR TICLE
22. Kafesaki, M., Shen, N. H., Tzortzakis, S. & Soukoulis, C. M. Optically switchable and tunable terahertz metamaterials through photoconductivity. J. Opt. 14, 114008 (2012). 23. Dykman, M. I., Maloney, C. M., Smelyanskiy, V. N. & Silverstein, M. Fluctuational phase-flip transitions in parametrically driven oscillators. Phys. Rev. E 57, 5202-5232 (1998). 24. Bishop, L. S., Ginossar, E. & Girvin, S. M. Response of the strongly driven Jaynes-Cummings oscillator. Phys. Rev. Lett. 105, 100505 (2010). 25. Unterreithmeier, Q. P., Faust, T. & Kotthaus, J. P. Nonlinear switching dynamics in a nanomechanical resonator. Phys. Rev. B 81, 241405(R) (2010).

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4730

Federation under contract no. 11.G34.31.0062. P.J. acknowledges financial support by the Helmholtz International Research School for Teratronics (HIRST).

Author contributions
P.J. proposed and set up the experiment and wrote the paper assisted by S.B., M.M. and A.V.U. The sample fabrication was managed by V.P.K. The experiment itself was carried out by S.B. and P.J., the theory developed by M.M., J.L. and M.V.F.

Additional information Acknowledgements
The authors would like to acknowledge interesting and productive discussions with S.M. Anlage and N. Lazarides. This work was supported by the EU project SOLID, the Deutsche Forschungsgemeinschaft (DFG), and the State of Baden-Wurttemberg ? through the DFG-Center for Functional Nanostructures (CFN) within subproject B3.5. This work was supported by the Ministry for Education and Science of Russian Competing financial interests: The authors declare no competing financial interests. Reprints and permission information is available online at http://npg.nature.com/ reprintsandpermissions/ How to cite this article: Jung, P. et al. Multistability and switching in a superconducting metamaterial. Nat. Commun. 5:3730 doi: 10.1038/ncomms4730 (2014).

6

NATURE COMMUNICATIONS | 5:3730 | DOI: 10.1038/ncomms4730 | www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.