# Slow waves in locally resonant metamaterials line defect waveguides

Nadège Kaina, Alexandre Causier, Yoan Bourlier, Mathias Fink, Thomas

Berthelot, Geoffroy Lerosey

► To cite this version:

Nadège Kaina, Alexandre Causier, Yoan Bourlier, Mathias Fink, Thomas Berthelot, et al. Slow waves in locally resonant metamaterials line defect waveguides. Scientific Reports, 2017, 7 (1), pp.15105. 10.1038/s41598-017-15403-8. cea-03811636

HAL Id: cea-03811636

https://cea.hal.science/cea-03811636

Submitted on 12 Oct 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Slow waves in locally resonant metamaterials line defect waveguides

Received: 24 July 2017

Accepted: 26 October 2017

Published online: 08 November 2017

Nadège Kaina<fnref n="1" />, Alexandre Causier², Yoan Bourlier², Mathias Fink<fnref n="1" />, Thomas Berthelot² & Geoffroy Lerosey<fnref n="1" />

Many efforts have been devoted to wave slowing, as it is essential, for instance, in analog signal computing and is one prerequisite for increased wave/matter interactions. Despite the interest of many communities, researches have mostly been conducted in optics, where wavelength-scaled structured composite media are promising candidates for compact slow light components. Yet their structural scale prevents them from being transposed to lower frequencies. Here, we propose to overcome this limitation using the deep sub-wavelength scale of locally resonant metamaterials. We experimentally show, in the microwave regime, that introducing coupled resonant defects in such metamaterials creates sub-wavelength waveguides in which wave propagation exhibit reduced group velocities. We qualitatively explain the mechanism underlying this slow wave propagation and demonstrate how it can be used to tune the velocity, achieving group indices as high as 227. We conclude by highlighting the three beneficial consequences of our line defect slow wave waveguides: (1) the sub-wavelength scale making it a compact platform for low frequencies (2) the large group indices that together with the extreme field confinement enables efficient wave/matter interactions and (3) the fact that, contrarily to other approaches, slow wave propagation does not occur at the expense of drastic bandwidth reductions.

Being able to temporally control the propagation of waves, and particularly achieving low group velocities, is one of the current challenges in wave physics, regarding its many related outcomes in both applied and fundamental physics. First, there are obvious benefits for the development of devices and the simplification of information processing through the design of delay lines and buffers. From a more fundamental point of view, slow waves are of primary importance to observe and exploit many fascinating, though weak, effects, of either quantum or classical physics. It indeed permits to enhance the wave-matter interaction probability by extending the interaction time in the propagation medium. Most commonly studied in the optical field, those interactions can for instance benefit low-power nonlinear physics<fnref n="1" />, efficient solar energy harvesting², miniature lasers³ or be used to control the spontaneous emission rate of atoms⁴. However, lower frequency wave-matter interactions are also relevant and hence widely studied, going from quantum QED⁵, masers⁶,⁷, electronic spin transition in NV centers⁸,⁹ or nonreciprocal devices. Moreover, the interest of slow wave propagation in applied physics is not the prerogative of short wavelengths since it finds numerous applications at lower frequencies, from the near infrared to microwaves, as in the fields of elastic waves and phononics. It is indeed required to enhance the efficiency of sensors, for the detection of hazardous products, for the development of phononic based circuits¹⁰ or for analog processing of microwaves in telecommunications¹¹, which actually often relies on modulations to optical frequencies¹² hence resulting in very high insertion losses.

It is nonetheless in optics, thanks to the extensive researches on structured based composite materials, namely photonic crystals¹³, that most of the components were developed. Those Bragg interference based media are indeed considered as natural candidates due to their ability to support line defect waveguides, named photonic crystal waveguides, or defect cavities which can be weakly coupled to create photonic crystal cavity waveguides¹⁴,¹⁵, both leading to slow wave propagation¹⁶⁻¹⁸ though implying different mechanisms. While the first ones rely on the dispersion induced in the waveguide¹⁹, the second, that are one particular implementation of

<footnote n="1">¹Institut Langevin, ESPCI ParisTech & CNRS, Paris, France. ²CEA Saclay, IRAMIS, NIMBE, LICSEN, UMR 3685, F-91191, Gif sur Yvette, France. Correspondence and requests for materials should be addressed to N.K. (email: nadegekaina@free.fr)</footnote>
coupled resonators optical waveguides (CROWs), exploit the strong dispersion brought by the resonances of its elementary constituents, namely, cavities²⁰⁻²².

Yet the use of these Bragg interference based media or CROWs sets several constraints, leading to two fundamental limitations. First, the overall dimensions of related components, as well as the wave spatial confinement, are set by the wavelength. Indeed, photonic crystal properties mainly rely on their structure, which is by definition wavelength scaled, as are the dimensions of elementary resonators in CROWs. Hence, if they remain reasonably small in optics, they obviously cannot be used as compact components at lower frequencies, as in the microwave domain for instance where they would end up on meter scaled components. Second, in both systems, slow wave propagation always comes at the expense of a bandwidth narrowing²³, coming either as a less trivial consequence of the wavelength scale of photonic crystal or the use of extremely resonant unit cells (equivalently very high quality factor resonators) for CROWs.

In this article, we investigate how both of these constraints can be released by transposing the previous slow wave concepts into the metamaterial field. Metamaterials, composite media generally composed of resonators organized at deep-subwavelength scales, indeed present effective properties that can be very high, a property that was used to image or focus waves below the diffraction limit²⁴⁻³⁰, near-zero³¹ or negative³²,³³. In a previous work³⁴, we evidenced that metamaterials negative effective properties can also be understood as bandgaps that originate from the interference between the plane waves propagating in the host matrix and the fields radiated by their resonant elements. This purely interference based interpretation of the origin of the metamaterials bandgap, that emphasizes analogies with the physics of photonic crystals, allowed us to demonstrate that local modifications of the material, analogous to photonic crystals, induce defect cavities able to confine waves to the defect region. Contrarily to structured composite media however, for which the bandgap is created by interferences on the periodic structure, and hence the defects are structural ones with sizes the order of the wavelength, the hybridization bandgap in metamaterial mainly stems from the resonant nature of its constituents³⁵,³⁶. Defect cavities then typically consist in modifying the resonance frequency of one element so that it falls within the hybridization bandgap of the surrounding medium. It results in wave confinements on much smaller dimensions and in very high Purcell factors³⁷, due to the typical deep sub-wavelength metamaterial scales. Following this idea, we further proved that inserting sequences of those resonant defects enables the creation of waveguides that can mold the flow of waves, again at deep subwavelength scales³⁴,³⁸. Through this initial proof of concept study, it was shown that the first constraint regarding the scale of components can be largely relaxed using metamaterials for which defect cavities and waveguides dimensions are independent of the wavelength.

Here, we concentrate on the temporal properties of such line defect waveguides in locally resonant metamaterials. While experimentally measuring the delays of pulses propagating in the waveguides, we first prove that the latter can slow down waves very effectively, with experimental group indices as high as $n_g = 227$, challenging the performances of photonic crystal based delay lines. We show that this value is furthermore largely tunable through modifications of the geometric properties of the surrounding bandgap medium. Moreover, recording both the transmission spectra and the delays enables us accessing the essential figure of merit that characterizes the efficiency of any slow wave device: the normalized delay-bandwidth product (NDBP), a dimensionless number independent of the length of the device that underlines how broadband the latter can operate given its slowing 'strength', i.e. its group index ng. It is defined as $NDBP = n_g * \Delta f/f_0$, where $\Delta f$ is the bandwidth and $f_0$ the central frequency of the device. In other composite media based slow wave devices, this NDBP is intrinsically limited: the slower the wave, the narrower the bandwidth³⁹. Though efforts have been made¹⁹,⁴⁰ to overcome this limitation, NDBP typical values remain less than unity. Strikingly, we observe in the defect line waveguides in metamaterials NDBPs at least one order of magnitude higher than proposals from the state of the art, coming mostly from the photonic crystal community⁴¹. More importantly, we show that, owing to their deep subwavelength scale, the NDBP in locally resonant metamaterials is no longer limited but can be easily increased by tailoring the density of the metamaterial line defect waveguides. This in turn overcomes the second usual intrinsic limitation of slow wave devices. We confirm this property via simulations and prove that the NDBP evolves as the inverse of the metamaterial typical spatial scale, independently of the wavelength, a property unique to the resonance driven physical mechanism of these media. Finally, we conclude by highlighting the three essential consequences of using line defect waveguides in metamaterials to slow down microwaves, namely deep sub-wavelength confinement, large group indices and extended bandwidth, while discussing how these can be of interest for both device implementations and wave-matter interactions, especially for low frequency waves.

## Results

### Measuring the spectral, spatial and temporal properties of a metamaterial line defect waveguide.
To support our demonstration, we experimentally and numerically study microwave samples based on arrays of quarter-wavelength resonant metallic wires, of given height L, sitting on a ground plane. For a defect created by shortening a wire to a length $L_d$ (or equivalently increasing its resonance frequency³⁴), this collection of resonators acts as a bandgap medium. By shortening several adjacent wires along a line, a metamaterial line defect waveguide is inserted, for which we measure both the temporal and spectral properties. To first explore the overall characteristics of the propagation in such metamaterial line defect waveguides, we manufacture a 3D printed periodic subwavelength square array (period a = 5 mm ~ $\lambda_0$/13) of L = 16 mm long copper-plated wires mounted on a ground plane (see SI for the 3D printing and copper-plating process). Such wires have a resonance frequency close to $f_0 = 4.5$ GHz, so that this medium opens a bandgap above $f_0$³⁴. The waveguide is then designed by replacing some of the wires along a tortuous path with $L_d = 14$ mm long defect wires, as depicted in Fig. 1a. To quantitatively measure the waveguide properties, two homemade input and output antennas, each one plugged to a network analyzer, feed the latter. Those antennas are designed to be impedance matched to the defect wires to maximize the coupling to the waveguide and avoid reflections (see SI). The acquired spectrum (Fig. 1c) first highlights the sudden drop of transmission above $f_0$, corresponding to the bandgap created by the long wires. It then
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Subwavelength waveguiding demonstration in a tortuous metamaterial waveguide. (a) Photograph of a 3D printed and chemically metallized wire array sitting on a ground plane comprising a tortuous defect line made out of shorter wires. To highlight the presence of the defect line, we represent in this picture a waveguide with defects separated by $a_g = 1.5$ mm, a distance three time smaller than in the sample we actually investigated the properties. (b) Real part of the experimental E-field at $f = 4.795$ GHz, mapped in superimposition with a top view of the photograph of the device. Dark circles represent wires in the surrounding medium while white ones highlight the defects waveguide path. Schematic representation of the measurement set-up for the spectral (resp. temporal) properties in red (resp. blue). (c) Spectrum of the transmission (dB). (d) (up) Temporal propagation of a short pulse centered on $f_c$ with both emitted (grey) and received pulses (blue) and (down) same for a 20 MHz wide Gaussian filtered pulse.</figcaption></figure>

evidences a transmission band of 6% bandwidth ($\Delta f = 300$ MHz) around a central frequency $f_c = 4.88$ GHz, close to the defects resonance. This band can be attributed to the presence of the line defect waveguide, as the spectral transmission of an array of all identical wires only displays a bandgap for this frequency range³⁴ (see SI). Note that for this band, the transmission spectrum is not completely flat, as we could have expected using impedance matched antennas. This is due to some small reflections that are experienced by the wave around the numerous corners of this meandering propagation path, that can however be minimized by adapting the waveguide geometry around the bends, as in photonic crystal waveguides⁴². To further assess that the defect line indeed acts as a deep-subwavelength waveguide, we record the electric field by scanning the medium over a wide frequency range with a homemade antenna mounted on a 2D moving stage, while the antenna at the output port of the waveguide is plugged in to a standard 50 Ω load (see SI). A map of one of the modes within the defects transmission band (Fig. 1b) highlights the effective deep subwavelength spatial confinement of the electric field around the defect wires (white circles). We see that the confinement in xy the propagation plane for the dimension transverse to the propagation is limited by the distance to the first neighbors in the surrounding medium (black circles). It corresponds to twice the lattice constant, i.e. $\lambda_0/6.5$. Note here that due to the defects resonance and the deep-subwavelength packing of the resonators; the field in the waveguide is also evanescent in the out of plane z-direction. This result confirms our first experimental proof of concept³⁴ while further demonstrating two things: on the one hand waves can be molded along arbitrary complex paths whereas on the other hand the coupling of adequate feeding antennas to the defect line waveguide ensures one of the essential prerequisite to turn it into an effective component. Finally, the temporal properties (the delay or equivalently the group velocity) of the microwaves propagating in the waveguide can be estimated. To do so, a short pulse is first sent in the input antenna and received at the output antenna (Fig. 1d, up), showing a strong delay of about 20 times the FWHM of the sent pulse though the waveguide is only a few wavelengths long. To quantitatively measure the group velocity, which is dispersive, we then choose to filter this pulse on a 20 MHz wide bandwidth centered on $f_c$ (Fig. 1d, down). While logically artificially broadening the signals, the filtering operation allows to evidence that the output pulse is received after a 55 ns delay that can be attributed to the propagation through the $L_g = 29.5$ cm long meandering path. This delay, equivalent to a propagation on 16.5 m in air, corresponds to a group index $n_g = 56$, and hence demonstrates that waves are strongly slowed down in the metamaterial defect waveguide, which is what we will now focus on.

**Origin of the slow wave propagation in metamaterial line defect waveguides.** To understand the origin of this very high group index and determine if it can be further enhanced, it is necessary to understand the physics underlying the propagation of waves in these line defect waveguides. The shorter wires forming the waveguide are resonant defects within the bandgap created by the resonance of the longer wires constituting the surrounding medium. The field is then necessarily tightly confined around each defect. The subwavelength transverse size of the waveguide furthermore proscribes the presence of propagating waves. Each of these defects can hence only be, in first approximation, evanescently coupled to its neighbors. Waves then tunnel from defect to defect through a tight binding like interaction, characterized by a coupling strength $\kappa$. This mechanism is very similar to that of CROWs in optics⁴³, which have been experimentally realized using ultra-high Q cavities such as ring resonators⁴⁴,⁴⁵ or defect cavities in photonic crystals¹⁵,⁴⁶. In microwave, we can name magneto-inductive
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Effect of the bandgap medium density on the transmission properties of the slow wave line defect waveguides. (a,b) Schematic representation of waveguides embedded in a bandgap medium with different densities and photographs of the four investigated devices. (c–g) Experimental results for all samples with densities ranging from a = 2 mm to a = 5 mm comprising: (c) the envelope of the filtered sent (black) and propagated (color) pulses at central frequency $f_c$, (d) the group index at central frequency $f_c$ as a function of $1/a^2$, (e) the transmission spectra (dB) (f), the dispersion relations displayed with normalized wavenumbers and with the frequency centered on the central frequency and normalized by the bandwidth for each waveguide and (g), The attenuation of the transmission plateau as a function of $a^{-2}$. Note that for the a = 2 mm medium, the attenuation is higher than expected due to a lower quality of the copper plating, inducing extra losses.</figcaption></figure>

resonator chains or capacitively coupled resonant microstrip lines as possible implementation of slow wave components based on this propagation mechanism. Those are however quite different from the proposed metamaterial line defect waveguides since our resonant wires do not naturally exhibit any capacitive or inductive coupling. In this case, the evanescent coupling occurs not between wire's resonances but rather between the fields that are confined due to the defects insertion. The second major difference is that metamaterials are structured at deep subwavelength scales, much smaller than that of CROWs whose resonant elements can be several wavelengths large, or even that of photonic crystals waveguides, which are typically wavelength scaled. We will see later that this feature considerably enhances the slowing properties of these waveguides and especially their corresponding NDBP. Knowing the propagation mechanism in the line defect waveguide, we now have a clear insight on how to further increase the group index $n_g$, because the group velocity is partly driven by the coupling strength $\kappa$. Indeed, the lower the coupling, i.e. the more difficult the tunneling, the slower the propagation and the higher the group index $n_g$. The idea is then simply to modulate $\kappa$, which can be done while acting on one of the two characteristics affecting the tunneling efficiency: the distance between defects (i.e. the period in the waveguide) or the field confinement around each defect. The latter is in practice set by the bandgap attenuation efficiency that can easily be modified by adjusting the density of resonators in the surrounding medium. The denser this medium is, the more efficient the bandgap attenuation becomes (Fig. 2a). The field around the defects in the waveguide is then more tightly confined, decreasing the coupling strength $\kappa$ and resulting in smaller group velocities. Note that such a modification of the surrounding medium's structure is a unique prerogative offered by metamaterials, since they are subwavelength scaled and governed by a local resonance rather than by Bragg interferences as in structured materials. The latter indeed have a determined periodicity which is constrained by the wavelength and that cannot be freely scaled without scaling the operating frequency as well.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Comparison of the delays and bandwidths of tight-binding dispersion relations chains in the case of coupled cavities with wavelength-scaled (photonic crystal/CROWs) in blue and deeply subwavelength (metamaterial) in red periodicities. (a) Comparison of the expected group index in case both materials present the same bandwidth. (b) Comparison of the expected bandwidth in case both materials present the group index around the central frequency.</figcaption></figure>

**Tuning n<sub>g</sub> through modifications of the defect-defect tunneling strength.** For the sake of simplicity, the group velocity tailoring is now demonstrated in straight waveguides. They are composed of 60 of the previous defect wires separated by a period a<sub>g</sub> = 1.5 mm ~ λ<sub>0</sub>/40. The surrounding medium's wires are packed in a periodic array of lattice constant a, varying from 5 mm to 2 mm in the four investigated samples (Fig. 2a), corresponding to periodicities ranging from λ<sub>0</sub>/13 to λ<sub>0</sub>/30. We first implement temporal measurements in order to evaluate the group velocity in our waveguides. For each sample, we send a 20 MHz wide pulse centered on the transmission band central frequency and measure the temporal delay of the received signal (see Fig. S7 in SI for short pulse measurements). The results, presented in Fig. 2c, evidence that shrinking the medium periodicity from a = 5 mm to a = 2 mm indeed strongly modulates the delay, that extends from 20 ns to 68 ns corresponding to group indices from n<sub>g</sub> = 65 to n<sub>g</sub> = 227. This strong delay is equivalent to free space propagation of 5.8 m to 20.4 m, although the waveguide's physical length is only L<sub>g</sub> = 9 cm. Such high n<sub>g</sub> is already comparable to the best performances obtained in optics with structured composite materials based slow light devices. Interestingly, the group index presents a linear dependence on the inverse of the medium density a<sup>2</sup> (Fig. 2d) and is hence highly tunable and improvable, since the lattice constant is theoretically solely limited by the minimum spacing between the wires, namely their radius. After successfully proving that the group velocity can be tailored, it is important to determine how the transmission bandwidth is affected by the velocity reduction, which is a natural concern regarding slow wave propagation. The measured transmission spectra for each waveguide (Fig. 2e) show that while increasing the group index, the transmission bandwidth accordingly shrinks from Δf = 610 MHz to Δf = 250 MHz, that is from 10% to 2% of the central frequency. This is not surprising since in a tight-binding model, the bandwidth Δf is directly proportional to the coupling strength κ. As we decreased κ to slow down microwaves, through the surrounding medium density, the bandwidth is naturally highly impacted. Note that for this linear geometry of waveguides, the transmission bandwidths are flat, confirming the impedance matching of the antennas.

**Delay-bandwidth products in locally resonant metamaterial-based line defects waveguides.** Despite this bandwidth reduction, one substantial observation is that Δf remains very large compared to what can be obtained in conventional composite structure slow wave devices. This is made particularly clear while calculating the corresponding NDBP for each waveguide, ranging from 7.5 (a = 5 mm) to 10.6 (a = 2 mm), values that are more than 10 times higher than any prior proposition based on CROWs or PC. This feature can be explained as a direct consequence of the deep subwavelength nature of metamaterials, as evidenced in Fig. 3. We represent here the schematic sinusoidal dispersion relation typical of tight binding interactions that take place in both defect cavities photonic crystal waveguides/CROWs and our line defect waveguides in metamaterials. These dispersion relations inform simultaneously about the group indices (through the slope of the curves, for instance around the central frequency) and about the total bandwidth that can be achieved in both kind of defect chains. To compare the two kind of materials, that are structured media and metamaterials, we consider two basic cases. In the first one (Fig. 3a), both linear chains of defects are supposed to exhibit the same bandwidth Δf, corresponding to an equivalent coupling strength κ. As it is pictured, the group index in the case of metamaterial line defect waveguides is then largely enhanced compared to the structured material case. This can be attributed to the difference of the first Brillouin zone extension, i.e. the maximum values of wavenumbers that can be reached for waves propagating in the waveguide (for the lower frequency branch we are interested in), in each kind of material. These wavenumbers are indeed only governed by the periodicity in the waveguide, through π/a<sub>g</sub>. In photonic crystal or CROWs structures (blue curve), a<sub>g</sub> between the defects is typically of the order of the wavelength while in metamaterials (red curve) this periodicity is deeply subwavelength. Our line defect waveguides can then support much higher wavenumber modes than other structured media based slow wave devices. This implies that for a given bandwidth, the group index can be much higher in metamaterials, which is geometrically observable through the difference of slopes in the dispersion relations close to the central frequency f<sub>c</sub> (Fig. 3a). In the second case, we assume that the group indices of both type of chains are now identical (the slopes are the same around the central frequency). Again, the schematic view of Fig. 3b demonstrates that the difference of first Brillouin zone extension implies that the bandwidth achieved in metamaterials is significantly larger than the one in the structured medium chain.
From the previous measurements, we however observe that the NDBP, though slightly increasing with the surrounding medium density, is far from being on a par with the very large group index modulation. This can easily be understood considering that the NDBP = n<sub>g</sub>*Δf/f<sub>0</sub> and the n<sub>g</sub> modulation is more or less compensated by the bandwidth shrinking. For a tight-binding based propagation indeed, assuming that the total bandwidth Δf is simply proportional to the linear effective bandwidth δf, the definition of the group index n<sub>g</sub> = c*δk/2πδf leads to NDBP ~ c*δk/f<sub>0</sub>. In this case, the NDBP is then independent of the bandwidth of the device and is solely proportional to the wavenumber of the modes propagating in the line defect waveguide. This explains why this quantity is restricted to small and relatively constant values in more conventional composite structures, because the wave vector is constrained by the wavelength.

In our case, though the wavenumbers are the same for all waveguides, we still observe a small variation in the NDBP. This can be attributed to a slight deviation from the tight-binding model that is highlighted by the experimental determination of the dispersion relations of the wave propagating in the waveguides (see SI for protocol). Indeed, having a closer glimpse at those dispersion curves (Fig. 2f), we see that they are not purely sinusoidal, as expected in CROWs for instance, especially for the largest densities. This highlights that the propagation mechanism, while remaining principally tight binding, also takes into account a small propagative part whose contribution depends on the strength of the field confinement around the defects. It gets more important when the bandgap attenuation efficiency is decreased, leading to dispersion relations presenting some features of a polariton, which is the expected dispersion relation when resonant wires are solely coupled through a propagating wave. Beyond these propagation mechanism considerations, we observe that these curves further confirm the group index modulation through the change in the slope near the central frequency.

Despite the small increment of the NDBP while increasing the surrounding medium density, it remains relatively limited, due to the drastic bandwidth reduction. This leads to the following question: can the NDBP be further enhanced or, in other words, can the group index be largely increased while limiting the bandwidth shrinking? In the previous samples, we only tuned the group index via the coupling strength κ, i.e. via the bandwidth δf but a<sub>g</sub>, and hence δk, remained constant, and so did the NDBP, as in structured media. Contrarily to the latter however, the wavenumber δk in metamaterials can be as well adjusted, providing a second lever that directly influences the NDBP, which is what we demonstrate in the following section.

**Enhanced delay-bandwidth products through wavenumber modulations.** To modulate the wavenumber of the modes within the waveguide, we again take advantage of metamaterials unique properties to modify the density of defects in the waveguide, namely a<sub>g</sub>. Decreasing this period indeed leads to higher wavenumbers and consequently increases the group index. However, putting the defects closer also has one opposite consequence on the velocity: it facilitates the propagation through the augmentation of the coupling strength, hence decreasing the group index. As we previously demonstrated though, this coupling can as well be controlled by the surrounding medium wires' density independently of the waveguide geometrical parameters. In order to limit the impact of the waveguide density on the coupling strength, we propose to implement a scaling of the device's dimensions that amounts to simultaneously modify the medium density (as in the previous paragraph) and the waveguide lattice constant. We then manufacture two components (Fig. 4a); the reference one (Sc = 1) is a straight waveguide with the lattice constants a = 10 mm (λ<sub>0</sub>/6.5) for the medium and a<sub>g</sub> = 6.97 mm (λ<sub>0</sub>/8.5) for a 18 cm long waveguide. The second one (Sc = 0.5) has a and a<sub>g</sub> reduced by half and so is the waveguide length L<sub>g</sub>. The spectral and temporal measurements demonstrate a striking result. While this scaling operation leads to a group delay multiplied by a factor 2.4, the transmission bandwidth reduces only by 25% (Fig. 4b,c,e,f). This value is far beyond the 50% bandwidth reduction that would be expected for an equivalent slowing strength when the coupling κ is alone modulated, keeping the wavenumber constant (Fig. 2e). As a direct consequence of the increase of the wavenumber δk of the slow wave modes, resulting from the geometrical scaling of the metamaterial, the NDBP is then multiplied by a factor 1.7, going from 4 to 7.5. This result completely overcomes the previous preconception of a systematic dramatic bandwidth reduction for large delay propagation in slow wave components. To complete the study, we again experimentally measure the dispersion relations and compare it with the theoretical tight-binding ones (Fig. 4d,g). As for the previous waveguides, we observe a slight deviation that can be attributed to a partly polaritonic-like behavior. This happens particularly for the Sc = 0.5 sample for which the smaller period in the waveguide may not completely prevent from direct propagative coupling. This could be avoided by decreasing the coupling between neighbors that is with higher medium density samples or larger periodicities in the waveguide.

**Delay-bandwidth products as a function of the metamaterial spatial scale.** To confirm these experimental results, we finally prove that the scaling operation can be further implemented and that the denser we scale the metamaterial samples, that is the larger the wavenumber of the modes are, the higher the delay-bandwidth products are. Because our fabrication process has a limited resolution imposed by our 3D printer, we simulate, using the software CST microwave studio, twelve of the previous devices with scaling factors ranging from Sc = 2 to Sc = 0.125. We observe the same trend as in experiment: when Sc decreases, while the normalized transmission bandwidth shrinks by a factor of 6, from 18% to 3% of the central frequency (Fig. 5a,b), the group index increases by more than two orders of magnitude (Fig. 5b). It results in NDBPs with very large values, up to 16, that seem to follow a linear dependence to a power of Sc close to -1 (Fig. 5c). It is consequently far from being limited to constant values. Again, we explain this by the fact that the NDBP value principally depends on the wavenumber of the waves propagating in the waveguide. The slight difference from the expected power -1 could be the consequence of the deviation of our dispersion relation from the pure tight-binding one that introduces a small dependence of the NDBP on the bandwidth. Note again that such a scaling operation is simply impossible to realize with other structured media since their spatial scale imposes their operating frequency, hence setting a fixed maximum value for δk.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Effect of the device scaling on the transmission properties of the slow wave line defect waveguides. (a) Schematic representations and pictures of the scaled samples Sc = 1 (left) and Sc = 0.5 (right). (b–d) Transmission spectrum (dB), envelope of the filtered sent (black) and propagated (color) pulses at central frequency $f_c$ and profiles of the dispersion relations along with the theoretical tight-binding model (dotted black) for the Sc = 1 sample. (e–g) Same for the Sc = 0.5 sample. The wavenumbers in **d** and **g** are normalized by the waveguide periodicity $a_g$ in the sample Sc = 1.</figcaption></figure>

Discussion

Studying these line defect metamaterial waveguides have enabled evidencing three key features that can together be beneficial to both spatial and temporal control of wave propagation, especially because they exhibit certain advantages over other typical implementations of slow-wave devices. First, it allows a wave guiding with transverse dimensions that are deeply subwavelength. This stems from the fact that the hybridization bandgap in locally resonant metamaterial, that is here exploited to create defects, is rather due to the resonance of its constituents than to their spatial organization, as it is the case in structured media. The main consequence is that it leads to an extreme compactness, with dimensions actually independent of the wavelength and whose sole limitation is the resonators size, providing though that there is no direct coupling between the latter. Second, the slowing strength of these line defect waveguides is very effective, with experimental group indices as high as 227 that can be increased by further densifying the bandgap medium. In other words, the more compact it is, the slower the wave propagates, without affecting too much the operating frequency, which again is a unique prerogative of metamaterial line defect waveguides. These first two characteristics combined, deep subwavelength wave confinement and long interaction time, perfectly fulfill the two requirements for increased wave-matter interactions. Finally, using metamaterials to create very dense line defect waveguides permits to obtain unprecedented and unrestricted delay-bandwidth products. As explained in Fig. 3, this comes as a less trivial consequence of the deep subwavelength scale of metamaterials allowing very large wavenumbers for the slow wave modes.

At the light of these findings, we believe that the concept we propose, line defect waveguides in metamaterials to route and slow down waves, demonstrated here in microwaves, can be of substantial interest, especially for low frequency waves. These domains indeed suffer a lack of practical solutions for both compact broadband components and wave-matter interaction enhancement platforms, despite the many applications they could benefit to. We can for instance mention microwave/atom interaction based phenomena such as NMR, widely used in medical imaging, MASERS<sup>6,7</sup>, quantum electrodynamics implemented with microwave cavities<sup>5</sup>, or even NV centers<sup>8,9</sup> and spin manipulation that are serious candidate for quantum computing. At larger scales, microwave/bulk material interaction are fundamental to observe phenomena stemming from magneto-electric or nonlinear effects for instance. Amongst potential applications are non-reciprocal devices, magnetic field sensors or data storage. From a more device-oriented point of view, in microwave to THz, compact delay lines as well as flexible wave guiding are required for all analog signals routing, radars, sensors, detection of hazardous products etc.
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Simulation of the effect of the sample scaling on the delay-bandwidth products of the slow wave line defect waveguides. (a) Simulated transmissions (dB) for six out of twelve devices with scaling factors from $Sc = 2$ to $Sc = 0.125$. (b) Evolution of the group index $n_g$ at the central frequency (y-axis left in blue, dotted blue line and black points) and the transmission bandwidth in % of the central frequency (y-axis right in red, dotted red line and black points) as a function of the inverse of the scaling factor $1/Sc$. The colored points on the two curves correspond to the values of scaling factors represented in (a,c). Normalized delay-bandwidth product as a function of $Sc^{-\alpha}$, where $\alpha = 0.84$ gives a linear dependence. Black dotted line is the linear fit and colored points correspond to the values of scaling factors represented in a.</figcaption></figure>

Finally, since our concept is general, it could be easily transferred in the field of acoustic and phononics, using Helmholtz resonators⁴⁷, pillars⁴⁸ or vibrating rods⁴⁹. Wave/matter interaction in the field of phononics are indeed equivalently important, with the development of acoustic analog signal processing using surface acoustic waves¹⁰. Hence, transposing to the metamaterial field concepts usually exploited with photonic crystal/CROWs, for which research has been very active in optics, we showed two things. First, the opening of new perspectives for low frequency wave manipulation (from microwave to THz, acoustics and elastic waves) and second, the evidence that the common preconception of severe bandwidth restrictions when slowing down waves is actually no fatality and is no more relevant while using metamaterials.

Naturally, the line defect waveguides suffer from some limitations. For instance, as many slow wave devices, they are dispersive (Figs 2f and 4d,g) so that the effective operating bandwidth is in practice restrained. We can however stress that the dispersion is linear (i.e. the propagation dispersionless) over quite a wide range of frequencies near the central frequency, in agreement with the principal tight-binding contribution in the propagation mechanism. This in turn differs drastically from photonic crystal waveguides slow wave propagation occurring only on the edge of the operating bandwidth¹⁹ or “rainbow trapping”⁵⁰,⁵¹ and slow wave approaches based on polaritonic dispersion systems as spoof plasmons, for which very large group indices only occur on the restrained flat band edge of the dispersion relation⁵²⁻⁵⁵.

Moreover, if any concrete application was to be designed, attention should be taken to minimize the losses that are particularly problematic when dealing with slow waves, given the long propagation time spent in the material. In our specific microwave design, the losses that can be mainly attributed to the metal intrinsic properties can be evaluated for the different components. To do so, we display (Fig. 2g) the transmission amplitude variation with the inverse of the density of the medium for the waveguides of Fig. 2. This amounts, given the linear variation of $n_g$ with $a^{-2}$, to evaluate the attenuation per time unit spent in the $Lg = 1.6 \lambda$ long media. We observe that the wave attenuation increases rather linearly with the group index, attesting of a 0.1 dB/ns loss. Though not negligible, the performance of the metamaterial line defect waveguides in terms of attenuation is then comparable to the best
coaxial cables while largely outperforming by one order of magnitude microstrip line based delay lines (see SI). It is moreover important to stress that, on the contrary to these common microwave devices, the fair attenuation performances of our system are coupled with an extreme compactness of the waveguide. The required propagation length to reach delays of tens of nanoseconds is indeed the order of the wavelength (few centimeters) while coaxial cables would have for instance been hundreds of wavelengths long (few meters). Usually, this bulk issue is overcome by implementing microwave photonics<sup>56</sup> components that require modulating microwave signals forth and back to optical frequencies before delaying them in compact photonic delay lines. Though quite functional, this method undergoes several sources of substantial losses (insertion loss, conversion...) while providing typical optical delays of only few picoseconds that barely fit microwave requirements. Using metamaterial line defect waveguides then seems to be conceptually of interest to slow down microwave signals, in terms of large delays, compactness and reasonable losses.

Of course, the metal intrinsic properties prevent from transposing this specific design (metallic resonant wires) directly in optics, where material losses would severely hinder propagation. However, we can imagine that the concept itself might as well be adaptable to higher frequencies, provided the use of more suitable materials or resonators (as cold atoms for instance) taking into account the inherent specificities of optics.

## Conclusion

To conclude, we demonstrated a new approach to slow down waves based on metamaterial subwavelength line defect waveguides. Owing to the subwavelength structure of metamaterials, we were able to demonstrate both very large group indices, here experimentally up to 227 and numerically up to 600, and relatively large operating bandwidths. Consequently, we measured unprecedented normalized delay-bandwidth products, more than one order a magnitude higher than prior propositions based on structured materials and that can be further increased by a simple scaling of the physical dimensions of the devices. We believe that our concept that achieves both extremely high confinement and overcomes the current limitations related to the delay-bandwidth product of slow wave devices can be of great interest in many fields of research relying on wave/matter interactions, especially for low frequency waves. It also paves the way to the design of ultra-compact components that spatially and temporally control the propagation of waves at scales independent of the wavelength. Finally, we would like to emphasize that though demonstrated here in microwave with a very specific design, our approach is very general. If it seems to suit best the requirements of low frequency regimes for now, the concept might in the future be transposed to other frequency ranges (up to THz and eventually optics) provided the use of appropriate materials. Furthermore, the idea can largely be adapted to other kind of resonators, as long as near-field interactions between unit cells can be safely neglected, or even for different type of waves such as acoustic or elastic waves.

## Methods

**Sample fabrication.** The samples were 3D printed in an ABS-like polymer resin using an Objet30 Pro 3D printer. The surface of the polymer was first annealed then submitted to an acidic treatment to produce chemical functions for the copper metallic ions chelation. The latter were chemically reduced to create copper nanoparticles or clusters onto the surface that acts as a seed layer catalyzing the copper electroless plating. To reach the copper bulk electrical properties and hence minimize the dissipation of waves in the polymer matrix, the thickness of the copper electroless plated layer was further increased over 20 µm (see SI for thickness control) by performing copper electroplating of the resulting conductive samples.

**Microwave measurements.** Measuring the spectral properties of the waveguide required homemade antennas consisting in cut stain wires soldered on SMA connectors, whose length was manually adapted in order to get the best possible coupling to the waveguide (see SI, Fig. S4). Spectral measurements were conducted while connecting both antennas to a beforehand calibrated network analyzer (Agilent Technologies N5230C) for frequencies ranging from 3 GHz to 7 GHz. Each measure was averaged 10 times.

As for the temporal measurement, one of the antennas was connected to a waveform generator (Tektronix AWG7102, 20 GS/s when interleaving channels) thus acting as a source, while the second was connected to an oscilloscope (Tektronix TDS6604B, 20GS/s). In order to probe the properties over the whole bandwidth of the waveguide, we sent for each device a wide pulse centered on the latter central frequency. The sent and transmitted pulses were then filtered on a 20 MHz bandwidth around the central frequency of each waveguide. The group index was retrieved from the delay between the sent and received filtered pulses that was measured from the maxima of the envelope of those pulses. The reference velocity was set to $3e^8$ m/s.

To measure the dispersion relations, we scanned the E-field on a line along the waveguide from the first wire to the last one, around 1 mm away from the top of the wires, with a step of $a_g/4$. The waveguide was fed by one of the previously described homemade antennas connected to one port of the network analyzer. The field probing was implemented with a homemade antenna consisting of a very short, hence inefficiently radiative that mainly probes the evanescent field, stain wire mounted on a SMA connector. The latter was plugged to the second port of the network analyzer. This antenna was mounted on a 2D translational stage (Newport M-IMS400PP). In order to prevent from reflections at the end of the waveguide that would create stationary waves, the second antenna in the sample was maintained (Fig. S5) but loaded with a 50 load. A spatial Fourier transform of the spectrum was performed and for each frequency, the wave vector with the maximum contribution in the FFT is kept to plot the dispersion curves as they are presented in Figs 2 and 4.

## References

1. Soljacic, M. & Joannopoulos, J. D. Enhancement of nonlinear effects using photonic crystals. *Nat. Mater.* **3**, 211–219 (2004).

2. Duché, D. et al. Slow Bloch modes for enhancing the absorption of light in thin films for photovoltaic cells. *Appl. Phys. Lett.* **92**, 1–4 (2008).

3. Painter, O. et al. Two-Dimensional Photonic Band-Gap Defect Mode Laser. *Science* **284**, 1819–1821 (1999).
4. Lodahl, P., Driel, A. F., Van, Nikolaev, I. S. & Irman, A. Controlling the dynamics of spontaneous emission from quantum dots by photonic crystals. *Nature* **430**, 8–11 (2004).

5. Wallraff, A. et al. Strong coupling of a single photon to a superconducting qubit using circuit quantum electrodynamics. *Nature* **431**, 162–167 (2004).

6. Oxborrow, M., Breeze, J. D. & Alford, N. M. Room-temperature solid-state maser. *Nature* **488**, 353–6 (2012).

7. Gordon, J., Zeiger, H. J. & Townes, C. H. The Maser-New Type of Microwave Amplifier, Frequency Standard, and Spectrometer. *Phys. Rev. Lett.* **99** (1955).

8. Jelezko, F., Gaebel, T., Popa, I., Gruber, A. & Wrachtrup, J. Observation of Coherent Oscillations in a Single Electron Spin. *Phys. Rev. Lett.* **92**, 76401 (2004).

9. Gaebel, T. et al. Room-temperature coherent coupling of single spins in diamond. *Nat. Phys.* **2**, 408–413 (2006).

10. Campbell, C. *Surface Acoustic Wave Devices and Their Signal Processing Applications*. (Elsevier, 2012).

11. Caloz, C., Gupta, S., Zhang, Q. & Nikfal, B. Analog Signal Processing: A Possible Alternative or Complement to Dominantly Digital Radio Schemes. *IEEE Microw. Mag.* **14**, 87–103 (2013).

12. Capmany, J. & Novak, D. Microwave photonics combines two worlds. *Nat. Photonics* **1**, 319–330 (2007).

13. Joannopoulos, J. D., Johnson, S. G., Winn, J. N. & Meade, R. D. *Photonic Crystals: Molding the Flow of Light (Second Edition)*. (Princeton University Press, 2008).

14. Notomi, M., Kuramochi, E. & Tanabe, T. Large-scale arrays of ultrahigh-Q coupled nanocavities. *Nat. Photonics* **2**, 741–747 (2008).

15. Matsuda, N., Kuramochi, E., Munro, W. J., Notomi, M. & Takesue, H. An on-chip coupled resonator optical waveguide single-photon buffer. *Nat. Commun.* **4**, 1–7 (2013).

16. Krauss, T. F. Why do we need slow light? *Nat. Photonics* **2**, 448–450 (2008).

17. Baba, T. & Kondo, K. Dynamic Control of Slow Light Pulses in Photonic Crystal Waveguides. *Transparent Opt. Networks (ICTON), 2014 16th Int. Conf. on, IEEE* **1–4** (2014).

18. Ek, S. et al. Slow-light-enhanced gain in active photonic crystal waveguides. *Nat. Commun.* **5**, 5039 (2014).

19. Baba, T. Slow light in photonic crystals. *Nat. Photonics* **2**, 465–473 (2008).

20. Yariv, A., Xu, Y., Lee, R. K. & Scherer, A. Coupled-resonator optical waveguide: a proposal and analysis. *Opt. Lett.* **24**, 711–713 (1999).

21. Morichetti, F., Ferrari, C., Canciamilla, A. & Mellononi, A. The first decade of coupled resonator optical waveguides: bringing slow light to applications. *Laser Photon. Rev.* **6**, 74–96 (2012).

22. Scheuer, J., Paloczi, G. T., Poon, J. K. S. & Yariv, A. Coupled Resonator Optical Waveguides Toward the Slowing & Storage of Light. *Opt. Photonics News* **16**, 36–40 (2005).

23. Miller, D. A. B. Fundamental limit to linear one-dimensional slow light structures. *Phys. Rev. Lett.* **99**, 203903 (2007).

24. Lemoult, F., Fink, M. & Lerosey, G. Revisiting the wire medium: an ideal resonant metalens. *Waves in Random and Complex Media* **21**(4), 591–613 (2011).

25. Lemoult, F., Fink, M. & Lerosey, G. Far-field sub-wavelength imaging and focusing using a wire medium based resonant metalens. *Waves in Random and Complex Media* **21**(4), 614–627 (2011).

26. Lemoult, F., Lerosey, G., De~Rosny, J. & Fink, M. Resonant Metalenses for Breaking the Diffraction Barrier. *Phys. Rev. Lett.* **104**, 203901 (2010).

27. Belov, P. A., Simovski, C. R. & Ikonen, P. Canalization of subwavelength images by electromagnetic crystals. *Phys. Rev. B* **71**, 193105 (2005).

28. Belov, P. A. & Silveirinha, G. Resolution of sub-wavelength transmission devices formed by a wire medium. 1–10 (2006).

29. Shvets, G., Trendafilov, S., Pendry, J. B. & Sarychev, A. Guiding, focusing, and sensing on the subwavelength scale using metallic wire arrays. *Phys. Rev. Lett.* **99**, 53903 (2007).

30. Kaina, N., Lemoult, F., Fink, M. & Lerosey, G. Negative refractive index and acoustic superlens from multiple scattering in single negative metamaterials. *Nature* **525**, 77–81 (2015).

31. Alü, A., Silveirinha, M., Salandrino, A. & Engheta, N. Epsilon-near-zero metamaterials and electromagnetic sources: Tailoring the radiation phase pattern. *Phys. Rev. B* **75**, 155410 (2007).

32. Fang, N. et al. Ultrasonic metamaterials with negative modulus. *Nat. Mater.* **5**, 452–6 (2006).

33. Alü, A., Salandrino, A. & Engheta, N. Negative effective permeability and left-handed materials at optical frequencies. *Opt. Express* **14**, 1557–67 (2006).

34. Lemoult, F., Kaina, N., Fink, M. & Lerosey, G. Wave propagation control at the deep subwavelength scale in metamaterials. *Nat. Phys.* **9**, 55–60 (2013).

35. Liu, Z. et al. Locally Resonant Sonic. *Materials. Science* **289**, 1734–1736 (2000).

36. Cowan, M. L., Page, J. H. & Sheng, P. Ultrasonic wave transport in a system of disordered resonant scatterers: Propagating resonant modes and hybridization gaps. *Phys. Rev. B* **84**, 94305 (2011).

37. Kaina, N., Lemoult, F., Fink, M. & Lerosey, G. Ultra small mode volume defect cavities in spatially ordered and disordered metamaterials. *Appl. Phys. Lett.* **102**, 144104 (2013).

38. Gao, Z., Gao, F. & Zhang, B. Guiding, bending, and splitting of coupled defect surface modes in a surface-wave photonic crystal. *Appl. Phys. Lett.* **108**, 41105 (2016).

39. Vlasov, Y. A., Hamann, H. F. & Mcnab, S. J. Active control of slow light on a chip with photonic crystal waveguides. *Nature* **438**, 65–69 (2005).

40. Yanik, M. F. & Fan, S. Stopping light all optically. *Phys. Rev. Lett.* **92**, 83901 (2004).

41. Hao, R. et al. Improvement of delay-bandwidth product in photonic crystal slow-light waveguides. *Opt. Express* **18**, 16309–16319 (2010).

42. Mekis, A. et al. High Transmission through Sharp Bends in Photonic Crystal Waveguides. *Phys. Rev. Lett.* **77**, 3787–3790 (1996).

43. Bayindir, M., Temelkuran, B. & Ozbay, E. Tight-binding description of the coupled defect modes in three-dimensional photonic crystals. *Phys. Rev. Lett.* **84**, 2140–3 (2000).

44. Xia, F., Sekaric, L. & Vlasov, Y. Ultracompact optical buffers on a silicon chip. *Nat. Photonics* **1**, 65–71 (2006).

45. Poon, J. K. S., Zhu, L., Derose, G. A. & Yariv, A. Transmission and group delay of microring coupled-resonator optical waveguides. *Opt. Lett.* **31**, 456–458 (2006).

46. Kuramochi, E. et al. Ultrahigh-Q photonic crystal nanocavities realized by the local width modulation of a line defect. *Appl. Phys. Lett.* **88**, 41112 (2006).

47. Theocharis, G., Richoux, O., Romero Garcia, V., Merkel, A. & Tournat, V. Limits of slow sound propagation and transparency in lossy, locally resonant periodic structures. *New J. Phys.* **16** (2014).

48. Achaoui, Y., Khelif, A., Benchabane, S., Robert, L. & Laude, V. Experimental observation of locally-resonant and Bragg band gaps for surface guided waves in a phononic crystal of pillars. *Phys. Rev. B* **83**, 104201 (2011).

49. Rupin, M., Lemoult, F., Lerosey, G. & Roux, P. Experimental Demonstration of Ordered and Disordered Multiresonant Metamaterials for Lamb Waves. *Phys. Rev. Lett.* **112**, 234301 (2014).

50. Tsakmakidis, K. L., Boardman, A. D. & Hess, O. 'Trapped rainbow' storage of light in metamaterials. *Nature* **450**, 397–401 (2007).

51. Gan, Q., Fu, Z., Ding, Y. J. & Bartoli, F. J. Ultrawide-bandwidth slow-light system based on thz plasmonic graded metallic grating structures. *Phys. Rev. Lett.* **100**, 1–4 (2008).
52. Williams, C. R. et al. Highly confined guiding of terahertz surface plasmon polaritons on structured metal surfaces. *Nat. Photonics* **2**, 175–179 (2008).

53. Hibbins, A. P., Evans, B. R. & Sambles, J. R. Experimental Verification of Designer Surface Plasmons. *Science* **308**, 670–672 (2005).

54. Fernández-Dominguez, A. I., Moreno, E., Martín-Moreno, L. & García-Vidal, F. J. Terahertz wedge plasmon polaritons. *Opt. Lett.* **34**, 2063–2065 (2009).

55. Gao, Z. & Zhang, B. Broadband wave manipulation in surface-wave photonic crystal. *ArXiv e-prints* 1–11 (2016).

56. Seeds, A. J. Microwave Photonics. *IEEE Trans. Microw. Theory Tech.* **50**, 877–887 (2002).

## Acknowledgements

N.K. acknowledges funding from French "Direction Generale de l'Armement". This work is supported by LABEX WIFI (Laboratory of Excellence within the French Program "Investments for the Future") under references ANR-10-LABX-24 and ANR-10-IDEX-0001-02 PSL*.

## Author Contributions

G.L. initiated and supervised the project, N.K. performed the simulations and experiments, G.L. and N.K. analyzed the results and developed the theory. A.C., Y.B. and T.B. developed the sample fabrication process, T.B. supervised the sample fabrication. All authors discussed the results and wrote the manuscript.

## Additional Information

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-017-15403-8.

**Competing Interests:** The authors declare that they have no competing interests.

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

<figure><img src="image_8.png" /></figure>

**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/.

© The Author(s) 2017