# Nonlinear Magnetoinductive Transmission Lines

## Abstract

Power transmission in one-dimensional nonlinear magnetic metamaterials driven at one end is investigated numerically and analytically in a wide frequency range. The nonlinear magnetic metamaterials are composed of varactor-loaded split-ring resonators which are coupled magnetically through their mutual inductances, forming thus a magnetoiductive transmisison line. In the linear limit, significant power transmission along the array only appears for frequencies inside the linear magnetoinductive wave band. We present analytical, closed form solutions for the magnetoinductive waves transmitting the power in this regime, and their discrete frequency dispesion. When nonlinearity is important, more frequency bands with singificant power transmission along the array may appear. In the equivalent circuit picture, the nonlinear magnetoiductive transmission line driven at one end by a relatively weak electromotive force, can be modeled by coupled resistive-inductive-capacitive (RLC) circuits with voltage-dependent capacitance. Extended numerical simulations reveal that power transmisison along the array is also possible in other than the linear frequency bands, which are located close to the nonlinear resonances of a single nonlinear RLC circuit. Moreover, the effectiveness of power transmission for driving frequencies in the nonlinear bands is comparable to that in the linear band. Power transmission in the nonlinear bands occurs through the linear modes of the system, and it is closely related to the instability of a mode that is localized at the driven site.

## I Introduction

The *metamaterial* concept usually refers to a periodic arrangement of
artificially structured elements, designed to achieve advantageous and/or unusual
electromagnetic properties compared to those of naturally occuring materials
Smith et al. (2004).
For example, the first metamaterial ever realized Shelby et al. (2001),
exhibited negative refraction index in a narrow frequency band around 5 Gigahertz.
The negative refraction index property requires that both the dielectric permittivity
and the magnetic permeability are simultaneously negative. The most widely used
elements for constructing a negative refractive index metamaterial are the
electrically small resonant ’particles’ called
*split-ring resonators* (SRRs) and metallic wires.
While the array of wires gives a negative permittivity below the plasma frequency,
the array of the SRRs gives a negative magnetic permeability
above their resonance frequency.
Thus, the refraction index can be negative in a narrow frequency band.
Since their discovery, metamaterials have been the target of intensive
study. Although the developement of metamaterials at microwave frequencies has
progressed to the point where scientists and engineers are now pursuing
applications, research on metamaterials that operate at higher frequencies
is still in an early stage Linden et al. (2006); Soukoulis et al. (2007); Litchinitser & Shalaev (2008).

There is a wide subclass of metamaterials, that exhibit significant magnetic
properties and even negative magnetic permeability at Terahertz and optical
frequencies Yen et al. (2004); Katsarakis et al. (2005); Linden et al. (2006), even though they are made
of non-magnetic materials. Magnetism at such high
frequencies is particularly important for the implementation of devices such as
tunable lenses, adaptive mirrors, etc., since there are only a few natural
materials that respond magnetically above microwaves. Moreover, magnetism in those
natural materials is usually weak and within narrow frequency bands,
limiting thus their use in possible Terahertz devices.
The most common realization of such a *magnetic metamaterial* is composed
of periodically arranged SRRs in a one- or two-dimentional lattices.
In its simplest version an SRR is just a highly conducting
Linden et al. (2004); Enkrich et al. (2005); Katsarakis et al. (2005) or superconducting Gu et al. (2010)
ring with a slit. The SRR structure was originally proposed by Pendry Pendry et al. (1999),
who theoretically predicted that a periodic SRR-based structure could exhibit a negative
magnetic permeability frequency band when irradiated by an alternating electromagnetic
field of appropriate polarization. Pendry also realized that the SRR structure
has considerable potential for enchancing nonlinear phenomena.
When the SRR dimensions are much smaller than the wavelength of the incident field,
the SRR can be regarded as an electric circuit consisting of inductive, capacitive
and ohmic elements connected in series. The ring forms the inductance , while the
slit can be considered as a capacitor of capacitance .
The ohmic element may represent all possible losses of the structure.

The next step was naturally to built nonlinearity in the SRR structure. Nonlinearity implies real-time tunability and multistability that is certainly a desired property of possible future devices. The SRRs can become nonlinear either by the insertion of a strongly nonlinear dielectric Hand & Cummer (2008) or a nonlinear electronic component (i.e., a varactor) Shadrivov et al. (2006); Powell et al. (2007); Wang et al. (2008) into their slits. Both ways would lead to effectively field dependent magnetic permeability for an array of nonlinear SRRs. Recently, the dynamic tunability by an external field for a two-dimensional array of varactor-loaded SRRs was demonstrated experimentally Shadrivov et al. (2008). The SRRs in such an array are coupled through magnetic and/or electric dipole-dipole interactions, whose strength depends on the relative orientation of the SRRs and their slits in the array Sydoruk et al. (2006); Hesmer et al. (2007); Penciu et al. (2008); Feth et al. (2010).

In the equivalent circuit picture, the dynamics of a periodic nonlinear SRR array
irradiated by an electromagnetic field, can be described by a set of coupled and driven
ordinary differential equations with on-site nonlinearity. The discreteness, which is
inherent in those SRR-based magnetic metamaterials, along with the nonlinearity and the
weak coupling between their elements,
allow for the generation of nonlinear excitations in the form of intrinsic localized
modes or *discrete breathers* Sievers & Takeno (1988).
These nonlinear modes appear generically in discrete and nonlinear extended systems.
Recent theoretical work in one- and two-dimensional nonlinear SRR lattices have
demonstrated the
existence and the stability of discrete breathers both in energy-conserving and
dissipative
systems Lazarides et al. (2006); Eleftheriou et al. (2008); Lazarides et al. (2008); Eleftheriou et al. (2009).
It has also demonstrated that breathers may be formed spontaneously through
modulational instability in binary nonlinear magnetic metamaterials
Molina et al. (2009); Lazarides et al. (2009, 2010).
Moreover, domain-wall excitations Shadrivov et al. (2006) and envelope solitons
Kourakis et al. (2007); Cui et al. (2009); Tsironis et al. (2010); Cui et al. (2010)
may as well be excited in those systems, which
seem to be stable even in the presence of noise Tsurumi (2008).

The SRR-based magnetic metamaterials support a new kind of electromagnetic waves, the magnetoinductive waves, which exhibit phonon-like dispersion curves, and they can transfer energy along the array Shamonina & Solymar (2004); Syms et al. (2006, 2010). It is thus possible to fabricate a contact-free data and power transfer device, a magnetoinductive transmission line, which make use of the unique properties of the magnetic metamaterial structure, and may function as a frequency-selective communication channel for devices via their magneto-inductive wave modes Stevens et al. (2010). In the present work we investigate the power transmission along a one-dimensional varactor-loaded SRR-based magnetic metamaterial, in the two possible geometries, which is driven at one end by a sinusoidal power sourse. In the linear limit, energy tranfer along the array occurs only for driving frequencies in the linear magnetoinductive wave band. However, the nonlinearity could generate more frequency bands where efficient power transmission along the array is possible. Power transmission in chains of coupled anharmonic oscillators for driving frequencies in the band gap of the linear spectrum has been recently investigated, and that effect is referred to as self-induced transparency Maniadis et al. (2006) or supratransmission Geniet & Leon (2002). For frequencies inside the linear magnetoinductive wave band the power is transmitted along the array from the magnetoinductive wave modes, slightly modified by the presence of nonlinearity. However, for frequencies in the pass-bands resulting from nonlinearity, the power is transmitted along the array from nonlinear, breather-like excitations. In the next Section we derive the model equations for a varactor-loaded SRR-based magnetic metamaterial, where the nonlinearity has the form a polynomial expansion on the main variable, like that investigated in Wang et al. (2008). We also present bifurcation diagrams for the shortest possible end-driven array, i.e., an array with only two varactor-loaded SRRs, from which only one is directly driven. In Section III we obtain analytical solutions for the linear magnetoinductive waves of the end-driven, finite SRR array when losses are neglected, along with their discrete dispersion relation. In Section IV we present numerical results for power transmission and discuss their dependence on the model parameters for a short, end-driven varactor-loaded SRR array with and without an absorbing boundary at the non-driven end. We finish with concluding remarks in Section V.

## Ii Varactor-Loaded SRR Array Model

Consider a ring-shaped or a squared split-ring resonator with a hyperabrupt tuning varactor mounted onto its slit. The varactor type could be selected from a large variety of available varactors, according to the needs of the experiment. In recent experimental works on varactor-loaded SRR-based metamaterials Wang et al. (2008); Poutrina et al. (2010); Larouche et al. (2010); Huang et al. (2010), the varactor selected was a Skyworks SMV1231-079, whose voltage-dependent capacitance ( is the voltage across the diode) is given by

(1) |

where is the DC rest capacitance, is the intrinsic potential, and is a parameter. The values of those parameters are provided by the manufacturer SPICE model for that varactor to be

(2) |

From Eq. (1) we can determine the voltage-dependence of the normalized charge

(3) |

with being the charge in the diode, as a function of the normalized diode voltage as

(4) |

Assuming that , the earlier equation can be solved for the diode voltage as a function of the charge, i.e.,

(5) |

When , the voltage in the earlier equation can be expanded in a Taylor series inpowers of around zero. Thus, neglecting term of order or higher, we get

(6) |

where, by using Eq. (2),

(7) |

The expression Eq. (6) can then be used in the voltage equation, obtained from Kirchhoff’s voltage law, which represents the varactor-loaded SRR as an effective resistive-inductive-capacitive (RLC) circuit with voltage-dependent capacitance

(8) |

where

(9) |

, is the induced electromotive force of amplitude and frequency resulting from an applied electromagnetic field with the same frequency, of the form

(10) |

where is the inductance of the ring, is the sum of the Ohmic resistance of the SRR ring and the series resistance of the diode, and the temporal variable. The eigenfrequency of the circuit described by Eq. (8) is given, in the linear limit without losses and driving, by

(11) |

That approximate model of the varactor-loaded SRR is found to be adequate for varactor voltages not exceeding ( in normalized units), when the dissipative current of the varactor can be neglected Wang et al. (2008).

Using Eq. (3) and the relations

(12) |

Eq. (8) can be written in the normalized form

(13) |

where

(14) |

is the normalized loss coefficient. The varactor-loaded SRR is a nonlinear oscillator, which exhibits multistability, hysteretic effects, and secondary resonances, when driven by a sufficiently strong external field. A typical current amplitude - frequency curve for the equivalent circuit model Eq. (II), like that shown in Fig. 1, exhibits all those characteristics. Notice the hysteresis loop at normalized frequency around , close to the resonance frequency of the linear system, and the very strong resonance at its second harmonic, i.e., at , where a smaller hysteresis loop also appears. Furthermore, we observe a subharmonic resonance at , as well as weaker resonances at , and . For frequency intervals lying between the boundaries of the hysteresis loops there are two different states of the oscillator, with low and high current amplitude, which are simultaneously stable.

Consider a one-dimensional periodic array of *identical* varactor-loaded SRRs,
which are coupled
magnetically to their first neighbors through their mutual inductances.
The array can be formed in two possible geometries, according to the relative
orientation of its elementary units (i.e., the SRRs) in the array.
As can be seen in Fig. 2, the SRRs can be arranged either in the planar geometry,
where all of them are lying on the same plane, or in the axial geometry,
where the axes of all the SRRs are lying on the same line.
The orientation of the incident electromagnetic field,
which is also shown in the figure, is such that the magnetic component can excite
an electromotive force in each SRR.
Below we consider arrays that are driven only at one end (let us say the left end).
Thus only one SRR receives energy directly from the driver. However, the energy can
be transmitted along the array due to the interaction between its elements.
Including the coupling of each SRR to its nearest neighbors, we get for the
charge of the th varactor the dynamic equations

(15) |

where , with being the total number of SRRs in the array, is the current flowing in the th SRR, is the approximate voltage across the th varactor, and (not to be confused with the varactor parameter) is the mutual inductance between two neighboring SRRs. Using the relations from Eq. (12), the coupled nonlinear Eqs. (15) can be written in normalized form as

(16) |

where is the coupling coefficient between neighboring SRRs. The function is unity for , while it is zero for any other . The coupling parameter may assume both positive and negative values, corresponding to planar and axial geometry of the varactor-loaded SRR array, respectively.

For a particular array, the values of the coupling and loss coefficients and , respectively, can be estimated from its geometrical and material parameters. Specifically, the mutual inductance , the ring inductance , and the Ohmic resistance of the ring can be calculated from the following expressions

(17) | |||||

(18) | |||||

(19) |

where is the magnetic permeability in vacuum, is the average radius of each SRR, the diameter of the wire of each SRR, is the center to center distance between neighboring SRRs, and and is the conductivity and skin depth of the SRR wire. Eqs. (17)–(19) refer to circular SRRs with circular cross-section. However, they can also be used for circular SRRs with square cross-section by defining an equivalent diameter of the SRR wire , with and being the depth and width, respectively, of the SRR wire. For , , , we get from Eqs. (17) and (18), respectively, that and , and consequently . For and (for an SRR made of copper wire at at room temperature) we get from Eq. (19) that . The diode series resistance of the varactor used in Ref. Wang et al. (2008) is around ; that is, most of the Ohmic resistance of the varactor-loaded SRR, comes from the varactor. By substituting in the second of Eq. (14), along with and , we get that . The resonance frequency of a single varactor-loaded SRR is . In the numerical calculations below we used the estimated value of that indicates weakly coupled SRRs. However, we use a smaller value of , which could be achieved in practice by using a metalic wire with higher conductivity and an ultra-low resistance varactor with similar capacitance-voltage relation.

The shortest end-driven varactor-loaded SRR array is that comprised of two elements. Since only one of them is driven by the external field, we refer to this system as the asymmetrically driven nonlinear dimer. The increase of degrees of freedom leads to interesting and complicated behavior, that is marked by the appearance of chaos. In Fig. 3, the bifurcation diagram of the total energy of the dimer with varying normalized frequency is shown for positive coupling . For negative coupling with equal strength we get very similar results that share a number of characteristic features. First of all, multistability is again present, for normalized frequencies in an interval around . There, two simultaneously stable energy states, with high and low energy, co-exist, as it is shown in Fig. 3a. The energy difference between those states increases condiderably with increasing frequency. Moreover, low energy state becomes chaotic for frequencies around . It is also interesting to see how the total energy is divided between the two elements. The energy fractions and , with and being the energies of the first (driven) and second SRR of the dimer, are shown both for the high and the low energy states in Fig. 3b and 3c as a function of the normalized frequency . When the dimer is in the low energy state, most of the energy is concentrated in driven SRR whose energy fraction is close to unity (correspondingly, the energy fraction of the other SRR is close to zero). However, for frequencies around the resonances, energy can be transfered easily from the one SRR to the other, and then the energy fractions of the SRRs in the dimer attain comparable values. In the chaotic regions, the energy is tranfered irregularly from one SRR to the other and the energy fractions span the whole interval between zero and one. In the high energy state (Fig. 3c), the energy fractions are almost equal for frequencies above the linear resonance frequency, . Below that frequency, the energy is either concentrated in the driven SRR (for frequencies far from resonances), or the values of the energy fractions are comparable (close to resonances.

## Iii Driven Linear Solutions

For a weakly driven array, the nonlinear terms which are proportional to and are not so important and they can be neglected. Then, the system of Eqs. (II) reduces to a linear one

(20) |

The dispersion relation of the linear system can be obtained by the substitution in the absence of losses and applied field, i.e., for and , which gives the linear dispersion of magnetoiductive waves

(21) |

where is the (normalized) wavevector (). Eq. (21) defines a frequency band of width that is bounded by a minimum and a maximum frequency and , respectively.

However, it is also possible to calculate the dispersion, expressed as a series of resonant frequencies, and the exact form of the linear magnetoinductive modes for the finite system in the absence of losses (). For finite systems, Eq. (III) should be implemented with free-end boundary conditions, to account for the termination of the structure, i.e.,

(22) |

By substituting into Eq. (III), the stationary equations can be written as

(23) |

where

(24) |

or, in matrix form

(25) |

where and are dimensional vectors with componets and , respectively, and is the inverse of the coupling matrix . The latter is a real, symmetric tridiagonal matrix that has diagonal elements equal to unity, while all the other non-zero elements are equal to . The need to find the inverse of tridiagonal matrices like arises in many scientific and engineering applications. Recently, Huang and McColl Huang & McColl (1997) related the inversion of a general trigiagonal matrix to second order linear recurrences, and they provided a set of very simple analytical formulae for the elements of the inverse matrix. Those formulae lead immediately to closed forms for certain trigiagonal matrices, like Lazarides & Tsironis (2010). The components of the vector can be written as

(26) |

where is the element of the inverse coupling matrix , whose explicit form is given in Lazarides & Tsironis (2010). Then, the solution of the linear system Eq. (III) and (22) for any driving frequency and finite is given by

(27) |

for and , and

(28) |

for , where

(29) |

Note that the s are uniquely determined by the parameters of the system.

From the analytical solution in the linear magentoinductive wave band, Eq. (III), which corresponds to either or , we can infer the resonant frequencies from zeroing the denominator, i.e., by setting . Then we see that the solution has a resonance for

(30) |

where is an integer (). The resonant freqquencies are obtained by solving the first of Eq. (24) with respect to , and substituting the resonant values of from Eq. (30). Then we get the discrete dynamic dispersion relation of linear magnetoinductive waves in the varactor-loaded SRR array, the discrete analogue of Eq. (21), as

(31) |

where is the mode number (). In Fig. 4, a density plot of the ’s on the plane, calculated analytically from Eqs. (III) and (III), is shown for positive coupling coefficient . In the frequency interval shown, the ’s have rather large values because, due to the narrow band, all frequencies are very close to a resonance. Moreover, the analytical solutions presented above have been obtained in the absence of losses, so that there is nothing to prevent the ’s from increasing indefinitely close to a resonance (and to go to infinity exactly at resonance). The white vertical regions in these figures correspond to frequency intervals around a resonance, where the ’s have very large values. Thus, those regions indicate the frequencies of the resonances. We can also see in Fig. 4 that the solutions close to the lower bound () of the linear magnetoinductive wave band exhibit very high amplitudes. For negative coupling coefficient with equal strength the high amplitude solutions appear close to the upper bound () of the linear magnetoinductive wave band.

In Fig. 5 the profiles of the analytical solutions for five different frequencies, both inside and outside the linear magnetoinductive wave band are shown for positive . The frequencies in Figs. 5a and 5e lie outside the linear band, so that the ’s decrease exponentially with increasing site number .

## Iv Power transmission along the array.

In the following we focus on the power transmission along a varactor-loaded SRR array with elements, which is driven at the left end (). The average power dissipated in the th SRR is defined as

(32) |

where denotes time-average over time-units. The power density can be normalized to , where and the value of resistance that results from the value of the dissipation coefficient used in the simulations. For we have that and . It is convenient to express in dBm, according to the relation

(33) |

where is the normalized current in the th SRR.

The system of Eqs. (II) is integrated with a 4th order Runge-Kutta algorithm with fixed time-stepping, and boundary conditions to account for the termination of the structure. The initial conditions were set to zero. The time allowed for the elimination of transients in each value of the frequency is ( ), while the time interval over which the average power is calculated is typically . In order to convert time intervals to physical units one has to multiply by the inverse of the resonance frequency , which gives . In the linear regime, i.e., and , the average power density plot on the site number - frequency plane shown in Fig. 6 exhibits the expected behavior; that appreciable power transmission only occurs for frequencies in the linear band, i.e., around .

However, the activation of the nonlinear terms reveals that even a relatively low driving amplitude is sufficient for nonlinear power transmission bands to bew formed. Thus, power can be transmitted not only for frequencies in the linear band, but also in the otherwise forbidden frequency regions. As it can be observed in Fig. 7, a number of transmission bands have appeared due to the nonlinearity. The array thus becomes transparent in power transmission at frequency intervals around the nonlinear resonances of a single varactor-loaded SRR oscillator. That self-induced transparency is a robust effect and survives with considerably increased losses, as it is shown in Fig. 8, where the loss coefficient has been increases by an order of magnitude (from to . In the density plots shown above, significant power transmission is indicated by red color. In both Figs. 7 and 8, nonlinear power transmission occurs at , , and . The latter nonlinear band, around the second harmonic of the resonance frequency of the linear system, seems to be the wider of all.

When losses are small, the appearance of more nonlinear bands is possible, like the one seen just above the linear band in Fig. 7, at . That transmission band is due to an interplay between nonlinearity and geometrical (lattice) resonances, and it is very sensitive to changes in the parameters of the model. Note that in the density plots presented above the frequency is given in natural units, using the relation .

It is very illuminating to see the time evolution of the current flowing in each SRR of the array for some frequencies in the transmission bands. We have chosen two different frequencies, one in the linear band and one the high frequency band which is around the second harmonic of the linear resonance frequency. The time evolution of the currents in the first, fourth, seventh, and tenth (last) SRR are shown for and in Figs. 9a and 9b, respectively. The currents oscillate periodically in all SRRs, although with different amplitudes and relative phases with respect to the driver. This figure suggests that power transmission for frequencies in the forbidden band gap occurs through the linear modes of the system. For frequencies outside the linear and nonlinear bands, the currents exhibit similar oscillations, although their amplitude becomes vanishingly small with increasing site number .

The relative efficiency of power transmission in the two frequencies used in the previous figure, are shown in Fig. 10, along with the power transmission for a specific frequency () outside all linear and nonlinear bands. It is remarkable that the transmission efficiency is almost the same for frequencies in the two pass-bands, while for the other one the transmitted power practically vanish at the fourth site.

In the density plots above, the right boundary was actually a reflecting one, allowing the formation of stationary states in the array. We have also used a totally absorbing boundary by adding five more sites at the right end of the array, where the damping coefficient increases exponentially with increasing site number. Thus, the waves coming from the left are completely damped at the end of the extended array, and no reflection occurs. The total number of SRRs is now where , and the dissipation coefficient is now a function of the site number , , given by

(34) |

where , with as before. For the array with dependent dissipation coefficient, a typical average power density plot is shown in Fig. 11, where the power in each site is the average over a time interval of . We observe, at least for the first sites of the array, a pattern similar with those seen in Figs. 7 and 8. That is, the appearance of the linear and the nonlinear transmission bands at about the same frequency intervals. We also observe that almost all of the transmitted power is absorbed from the five last sites of the array. Importantly, the narrow band coming from geometrical resonances in Figs. 7 and 8, which is around , appears in Fig. 11 as well, and moreover it seems that it is capable of transmitting power to the end of the array.

## V Discussion and Concluding Remarks

We have used a simple electric circuit model of inductively coupled resistive-inductive-capacitive circuits with voltage-dependent capacitance to investigate the transmission of power in a one-dimensional, end-driven, nonlinear magnetic metamaterial comprised of varactor-loaded SRRs. We have used both a reflecting and an absorbing boundary for the end that is not driven by the external field. The importance of such discrete models for describing the dynamic behaviour of SRR-based magnetic metamaterials has been recently appreciated Liu et al. (2009), although similar discrete models have been employed earlier for the study of localized excitations in those systems Shadrivov et al. (2006); Lazarides et al. (2006)

In the present work, we have focused on the transmission of power along the array, and especially that in the nonlinear regime, which results in the formation of nonlinear pass-bands for frequencies close to the resonances of a single varactor-loaded SRR that is described by the approximate model equation (II). The approximate model holds fairly well for relatively low driving amplitudes, for which the dissipative diode current of the varactor can be neglected. It seems, thus that the properties of the individual elements of the array predominantly determine its transmission properties, at least in the case of weak magnetoinductive coupling investigated in the present work. The numerical results demonstrate that the nonlinear pass-bands transmit power with efficiently comparable to that of the linear band. For the parameter sets used in this work, which are close to those in Wang et al. (2008) for a single varactor-loaded SRR, there are three nonlinear bands, both above and below the linear band. Their width apparently depends on the frequency, i.e., the bands at higher frequencies are significanlty wider than those at lower frequencies.

Moreover, we have also observed the formation, either complete or partial, of another nonlinear pass-band at frequencies slightly above the linear band (at ). That band, which is rather narrow and can be observed both in Figs. 7 and 8 as well as in Fig. 11, seem to depend sensitively on the model parameters , , (for and fixed), and mainly on the loss coefficient . Comparing Figs. 7 and 8, which have been obtained for different loss coefficient and , respectively, we see that the power transmission in that band considerably decreases when the losses increase by an order of magnitude. That type of nonlinear band appears due to geometrical resonances of the standing waves that are formed in the array due to driving at one end. Geometrical resonances seem to offer an important mechanism for the generation of nonlinear pass-bands, which can be employed for constructing transmission lines with banded power transmission spectra with bands at desired frequency ranges. Similar conduct-free transmission lines could be formed by superconducting magnetic metamaterials, where the basic structural element (the SRR) is replaced by its direct superconducting analogue, the rf SQUID Lazarides & Tsironis (2007); Lazarides et al. (2008). The rf SQUID, where the acronym stands for radio-frequency superconductive quantum interference device, is just a superconducting ring interrupted by a Josephson junction, which makes it an intrinsically nonlinear element. A one-dimensional array of rf SQUIDs coupled through their mutual inductances forms a nonlinear magnetoinductive transmission line Lazarides et al. (2008) whose power transmission properties is a matter of future work.

## Acknowledgments

This work was supported in part by the European Office of Aerospace Research and Development, AFSOR award FA8655-10-1-3039, and by the EURYI and MEXT-CT-2006-039047 .

### References

- Cui, W., Zhu, Y.-Y., Li, H.-x. & Liu, S. [2009] “Self-induced gap solitons in nonlinear magnetic metamaterials,” Phys. Rev. E 80,036608-1–036608-5.
- Cui, W., Zhu, Y.-Y., Li, H.-X. & Liu, S. [2010] “Solitons excitations in a one-dimensional nonlinear diatomic chain of split-ring resonators,” Phys. Rev. E 81, 016604-1–016604-9.
- Eleftheriou, M., Lazarides, N. & Tsironis, G. P. [2008] “Magnetoinductive breathers in metamaterials,” Phys. Rev. E 77, 036608-1–036608-13.
- Eleftheriou, M., Lazarides, N., Tsironis, G. P. & Kivshar, Yu. S. [2009] “Surface magnetoinductive breathers in two-dimensional magnetic metamaterials,” Phys. Rev. E 80, 017601-1–017601-4.
- Enkrich, C., [2005] “Magnetic metamaterials at telecommunication and visible frequencies,” Phys. Rev. Lett. 95, 203901-1–203901-4.
- Feth, N., König, M., Husnik, M., Stannigel, K., Niegemann, J., Busch, K., Wegener, M. & Linden, S. [2010] “Electromagnetic interaction of spit-ring resonators: The role of separation and relative orientation,” Opt. Express 18, 6545–6554.
- Geniet, F. & Leon, J. [2002] “Energy transmission in the forbidden band gap of a nonlinear chain,” Phys. Rev. Lett 89, 134102-1–134102-4.
- Gu, J., Singh, R., Tian, Z., Cao, W., Xing, Q., He, M.-X., Zhang, J. W., Han, J., Chen, H. & Zhang, W. [2010] “Terahertz superconductor metamaterial,” Appl. Phys. Lett. 97, 071102-1–071102-3.
- Hand, T. H. & Cummer, S. A. [2008] “Frequency tunable electromagnetic metamaterial using ferroelectric loaded split rings,” J. Appl. Phys. 103, 066105-1–066105-3.
- Hesmer, F., Tatartschuk, E., Zhuromskyy, O., Radkovskaya, A. A., Shamonin, M., Hao, T., Stevens, C. J., Faulkner, G., Edwardds, D. J. & Shamonina, E. [2007] “Coupling mechanisms for split-ring resonators:Theory and experiment,” Phys. Stat. Sol. (B) 244, 1170–1175.
- Huang, D., Poutrina, E. & Smith, D. R. [2010] “Analysis of the power dependent tuning of a varactor-loaded metamaterial at microwave frequencies,” Appl. Phys. Lett. 96, 104104-1–104104-3.
- Huang, Y. & McColl, W. F. [1997] “Analytical inversion of general tridiagonal matrices,” J. Phys. A: Math. Gen. 30, 7919–7933.
- Katsarakis, N., [2005] “Magnetic response of split-ring resonators in the far infrared frequency regime,” Opt. Lett. 30, 1348–1350.
- Kourakis, I., Lazarides, N. & Tsironis, G. P. [2007] “Self-focusing and envelope pulse generation in nonlinear magnetic metamaterials,” Phys. Rev. E 75, 067601-1–067601-4.
- Larouche, S., Rose, A., Poutrina, E., Huang, D. & Smith, D. R. [2010] “Experimental determination of the quadratic nonlinear magnetic susceptibility of a varactor-loaded split ring resonator metamaterial,” Appl. Phys. Lett. 97, 011109-1–011109-3.
- Lazarides, N., Eleftheriou, M. & Tsironis, G. P. [2006] “Discrete breathers in nonlinear magnetic metamaterials,” Phys. Rev. Lett. 97, 157406-1–157406-4.
- Lazarides, N. & Tsironis, G. P. [2007] “rf superconducting quantum interference device metamaterials,” Appl. Phys. Lett. 90, 163501-1–163501-3.
- Lazarides, N., Tsironis, G. P. & Kivshar, Yu. S. [2008] “Surface breathers in discrete magnetic metamaterials,” Phys. Rev. E 77, 065601(R)-1–065601(R)-4.
- Lazarides, N., Tsironis, G. P. & Eleftheriou, M. [2008] “Dissipative discrete breathers in rf SQUID metamaterials,” Nonlin. Phen. Compl. Syst. 11, 250–258.
- Lazarides, N., Molina, M. I. & Tsironis, G. P. [2009] “Breather induction by modulational instability in binary metamaterials,” Acta Phys. Pol. A 116, 635-637.
- Lazarides, N. & Tsironis, G. P. [2010] “Driven linear modes: Analytical solutions for finite discrete systems,” Phys. Lett. A 374, 2179–2181.
- Lazarides, N., Molina, M. I. & Tsironis, G. P. [2010] “Breathers in one-dimensional binary metamaterial models,” Physica B 405, 3007-3011.
- Lazarides, N., Molina, M. I., Tsironis, G. P. & Kivshar, Yu. S. [2010] “Multistability and localization in coupled nonlinear split-ring resonators,” Phys. Lett. A 374, 2095-2097.
- Linden, S., Enkrich, C., Dolling, G., Klein, M. W., Zhou, J., Koschny, T., Soukoulis, C. M., Burger, S., Schmidt, F. & Wegener, M. [2004] “Magnetic response of metamaterials at 100 Terahertz,” Science 306, 1351–1353.
- Linden, S., Enkrich, C., Dolling, G., Klein, M. W., Zhou, J., Koschny, T., Soukoulis, C. M., Burger, S., Schmidt, F. & Wegener, M. [2006] “Photonic metamaterials: magnetism at optical frequencies,” IEEE J. Selec. Top. Quant. Electron. 12, 1097–1105.
- Litchinitser, N. M. & Shalaev, V. M. [2008] “Photonic metamaterials,” Laser Phys. Lett. 5, 411–420.
- Liu, N., Liu, H., Zhu, S. N. & Giessen, H. [2009] “Stereometamaterials,” Nature Photon. 3, 157– .
- Maniadis, P., Kopidakis, G. & Aubry, S. [2006] “Energy dissipation threshold and self-induced transparency in systems with discrete breathers,” Physica D 216, 121–135.
- Molina, M. I., Lazarides, N. & Tsironis, G. P. [2009] “Bulk and surface magnetoinductive breathers in binary metamaterials,” Phys. Rev. E 80, 046605-1–046605-9.
- Pendry, J. B., Holden, A. J., Robbins, D. J. & Stewart, W. J. [1999] “Magnetism from conductors and enchanced nonlinear phenomena,” IEEE Trans. Microwave Theory Tech. 47, 2075–2084.
- Penciu, R. S., Aydin, K., Kafesaki, M., Koschny, Th., Ozbay, E., Economou, E. N. & soukoulis, C. M. [2008] “Multi-gap individual and coupled split-ring resonator structures,” Opt. Express 16, 18131–18144.
- Poutrina, E., Huang, D. & Smith, D. R. [2010] “Analysis of nonlinear electromagnetic metamaterials,” New J. Phys. 12, 093010-1–093010-27.
- Powell, D. A., Shadrivov, I. V., Kivshar, Yu. S. & Gorkunov, M. V. [2007] “Self-tuning mechanisms of nonlinear split-ring resonators,” Appl. Phys. Lett. 91, 144107-1–144107-3.
- Shadrivov, I. V., Morrison, S. K. & Kivshar, Yu. S. [2006] “Tunable split-ring resonators for nonlinear negative-index metamaterials,” Opt. Express 14, 9344–9349.
- Shadrivov, I. V., Zharov, A. A., Zharova, N. A. & Kivshar Yu. S. [2006] “Nonlinear magnetoinductive waves and domain walls in composite metamaterials,” Photonics Nanostruct. Fundam. Appl. 4, 69–74.
- Shadrivov, I. V., Reznik, A. N. & Kivshar Yu. S. [2007] “Magnetoinductive waves in arrays of split-ring resonators,” Physica B 394, 180–183.
- Shadrivov, I. V., Kozyrev, A. B., van der Weide, D. W. & Kivshar, Yu. S. [2008] “Tunable transmission and harmonic generation in nonlinear metamaterials,” Appl. Phys. Lett. 93, 161903-1–161903-3.
- Shalaev, V.M. [2007] , “Optical negative-index metamaterials,” Nature Photonics 1, 41–48.
- Shamonina, E. & Solymar, L. [2004] “Magneto-inductive waves supported by metamaterials elements: components for a one-dimensinal waveguide,” J. Phys. D: Appl. Phys. 37, 362– .
- Shelby, R. A., Smith, D. R. & Schultz, S. [2001] “Experimental verification of a negative index of refraction,” Science 292, 77-79.
- Sievers, C. & Takeno, C. [1988] “Intrinsic localized modes in anharmonic crystals,” Phys. Rev. Lett. 61, 970–973.
- Smith, D. R., Pendry, J. B. & Wiltshire, M. C. K. [2004] “Metamaterials and negative refractive index,” Science 305, 788–792.
- Soukoulis, C. M., Linden, S. & Wegener, M. [2007] “Negative refractive index at optical wavelengths,” Science 315, 47–49.
- Stevens, C. J., Chan, C. W. T. Stamatis, K. & Edwards, D. J. [2010] “Magnetic metamaterials as 1-D data tranfer channels: An application for magneto-inductive waves,” IEEE Trans. Microw. Theory Techniques 58, 1248–1256.
- Sydoruk, O., Radkovskaya, A., Zhuromskyy, O., Shamonina, E., Shamonin, M., Stevens, C. J., Faulnker, G., Edwards, D. J. & Solymar, L. [2006] “Tailoring the near-field guiding properties of magnetic metamaterials with two resonant elements per unit cell,” Phys. Rev. B 73, 224406-1–224406-12.
- Syms, R. R. A., Young, I. R. & Solymar, L. [2006] “Low-loss magneto-inductive waveguides,” J. Phys. D: Appl. Phys. 39, 3945– .
- Syms, R. R. A., Solymar, L., Young, I. R. & Floume, T. [2010] “Thin-film magneto-inductive cables,” J. Phys. D: Appl. Phys. 43, 055102-1–055102-7.
- Tsironis, G. P., Lazarides, N. & Eleftheriou, M. [2010] “Discrete breathers and solitons in metamaterials,” Springer Series in Optical Sciences 150, 273–285.
- Tsurumi, T. [2008] “Soliton Propagation in Nonlinear Magnetic Metamaterials with Microscopic Disorder,” J. Phys. Soc. Jpn. 77, 074006-1–074006-8.
- Wang, B., Zhou, J., Koschny, T. & Soukoulis C. M. [2008] “Nonlinear properties of split-ring resonators,” Opt. Express 16, 16058–16064.
- Yen, T. J., Padilla, W. J., Fang, N., Vier, D. C., Smith, D. R., Pendry, J. B., Basov, D. N. & Zhang, X. [2004] “Terahertz magnetic response from artificial materials,” Science 303, 1494–1496.