# Mimicking Electromagnetic Wave Coupling in Tokamak Plasma with Fishnet Metamaterials

K. Rustomji, R. Abdeddaim, J. Achard, M. Chmiaa, E. Georget, M. Goniche, W. Helou, J. Hillairet, S. Enoch, G. Tayeb

► To cite this version:

K. Rustomji, R. Abdeddaim, J. Achard, M. Chmiaa, E. Georget, et al., Mimicking Electromagnetic Wave Coupling in Tokamak Plasma with Fishnet Metamaterials. Scientific Reports, 2018, 8 (1), 10.1038/s41598-018-24250-0. [hal-01786013](https://hal.science/hal-01786013)

HAL Id: hal-01786013

https://hal.science/hal-01786013

Submitted on 19 Jun 2023

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

# Mimicking Electromagnetic Wave Coupling in Tokamak Plasma with Fishnet Metamaterials

Received: 4 October 2017

Accepted: 26 March 2018

Published online: 11 April 2018

K. Rustomji<fnref n="1" />,², R. Abdeddaim<fnref n="1" />, J. Achard³, M. Chmiaa<fnref n="1" />, E. Georget⁴, M. Goniche³, W. Helou³, J. Hillairet³, S. Enoch<fnref n="1" /> & G. Tayeb<fnref n="1" />

This paper reports a fishnet hyperbolic metamaterial that mimics the electromagnetic properties of magnetically confined plasma. These electromagnetic properties are strongly anisotropic and different from any conventional material, therefore cannot be mimicked by bulk materials. The structure is made of a stack of thin copper grids spaced by Rohacell foam. We numerically and experimentally show that this kind of structuration matches well the properties of a homogeneous plasma. This solution breaks a long-lasting bottleneck and will accelerate the development of high-frequency heating systems to be used in nuclear fusion.

Tokamaks are expected to become the standard choice for magnetic confinement devices to produce nuclear fusion, an environmentally friendly solution for future large-scale energy supply. International projects such as ITER are leading the way in this field. To transfer energy into the plasma, dedicated radio frequency antennas are required. Their development and optimization require the availability of a loading material that mimics the plasma. Until now, there exists no realistic load able to simulate a tokamak plasma, and consequently no way to test the antennas that will be used to create the plasma.

The problem is to transfer power from the antenna to the plasma, and test these antennas under relevant conditions. The power transfer is mainly related to the characteristics of the plasma at the tokamak edge (location of the antenna's mouth).

In this paper, the idea is to design and build a metamaterial to mimic the plasma. It appears legitimate to simplify the characteristics of the edge plasma by tending towards hyperbolic materials. Our design allowed us to build a metamaterial that overcomes the problem of spatial dispersion encountered in previous studies. We carried out numerical simulations and measurements, which confirm the validity of this concept. As a result, we now have a way to characterize and test the antennas. One of the other advantages is that these tests can be done at low power, without the need for elaborate sources or cooling. It should be noted that we do not intend to model with our metamaterials actual tokamak plasma, but only simplified plasma from the tokamak edge that retains the essential properties related to the coupling of the waves launched by the antenna to the plasma.

## Design and Fabrication

The coupling of Lower Hybrid Range of Frequencies (LHRF, 1–8 GHz) waves to strongly magnetized plasmas is a critical issue for tokamaks as it often limits the RF power, which can be transferred from the antenna to the plasma. Development of new types of antennas to improve the ability of the antenna to handle large power in stationary conditions, as required on a fusion reactor, is hampered by the long and costly delay between the design and the feedback from experiments on large facilities such as tokamaks. Numerical codes that model both the plasma and the antenna are now available and provide an accurate characterization of the RF coupling. Before a full-scale test of the antenna on the plasma, the test of a mock-up with a load mimicking the plasma is believed to be a step forward by reducing the risks and accelerating the development process.

The magnetic configuration of a tokamak requires a strong DC current flowing along the plasma ring<fnref n="1" />. Generally, this current is inductively driven by the solenoid located on the tokamak axis. However, in order to achieve continuous operation of the reactor, an external source of current is needed. RF waves emitted from the

<footnote n="1">¹Aix Marseille Univ, CNRS, Centrale Marseille, Institut Fresnel, Marseille, France. ²Centre for Ultrahigh Bandwidth Devices for Optical Systems and Institute of Photonics and Optical Science, School of Physics, University of Sydney, Sydney, NSW 2006, Australia. ³CEA, IRFM, F-13108, St-Paul-Lez-Durance, France. ⁴CEA-Saclay, DRF/12BM/ Neurospin/UNIRS, 91191, Gif-sur-Yvette Cedex, France. Correspondence and requests for materials should be addressed to G.T. (email: gerard.tayeb@fresnel.fr)</footnote>
<figure><img src="image_3.png" /><figcaption>Figure 1. Top: sketch of an antenna made of 6 waveguides, heating tokamak plasma (not to scale). The yellow lines depict the plasma current. For a picture of an actual set of waveguides feeding a tokamak, see the supplementary information. Bottom: schematic view of the metamaterial fishnet load.</figcaption></figure>

plasma periphery are well fitted for this task when the wavevector spectrum is properly chosen in order to transfer the wave energy to the electrons of the plasma current (Fig. 1). Among these waves, LHRF waves are a good candidate, thanks to their high current drive efficiency².

The dielectric tensor of magnetized plasma in the cold plasma approximation can be written in the form:³

$$ \epsilon = \epsilon_0 \begin{bmatrix} \epsilon_{\perp} & \epsilon_A & 0 \\ \epsilon_A & \epsilon_{\perp} & 0 \\ 0 & 0 & \epsilon_{\parallel} \end{bmatrix} $$

where the subscript/denotes the direction of the confinement magnetic field, i.e., the z-axis on Fig. 1.

In the case of LHRF waves of frequency $\omega_0/2\pi$ in the 1–8 GHz range, the electron density in the vicinity of the antenna has to be optimized to minimize the RF power reflected towards the generator and the electric field in the waveguides. This is achieved when the density $n_e$ exceeds the cut-off density ($n_{cut-off} = (\epsilon_0 m_e/e^2) \omega_0^2$) typically by a factor of 3. Under such conditions, the non-diagonal term $\epsilon_A$ is small ($\epsilon_A \sim 0.1$) and $\epsilon_{\perp}$ differs from 1 by less than 1%. Consequently, $\epsilon$ can be simplified using $\epsilon_{\perp} \approx 1$ and $\epsilon_A \approx 0$. The third parameter $\epsilon_{\parallel}$ is approximatively reduced to a Drude model $\epsilon_{\parallel} = 1 - \omega_{pe}^2/\omega_0^2 = 1 - n_e/n_{cut-off}$ where $\omega_{pe}/2\pi$ is the plasma frequency ($\omega_{pe}^2 = n_e e^2/(\epsilon_0 m_e)$).

Good coupling conditions of the wave to the plasma are achieved when the relative permittivity in the direction of the magnetic field is negative ($\epsilon_{\parallel} < 0$). The density of the plasma layer facing the antenna is far from being spatially homogeneous and constant with time. There is usually a strong radial density gradient in the vicinity of the antenna with a typical gradient length of a few centimeters, shorter than the LHRF wavelength. However RF modeling shows that the coupling properties depend weakly on the gradient and, moreover, the simplification of the dielectric tensor ($\epsilon_{\perp} \approx 1$ and $\epsilon_A \approx 0$) still holds for the plasma layer close to the antenna, which matters for the coupling of the wave. In high performance plasma, instabilities occur in the plasma; the density at the plasma edge varies periodically with time from a low value ($\epsilon_{\parallel} \sim -1$) to a high value ($\epsilon_{\parallel} \sim -10)^3$ and can depart from the optimal coupling conditions achieved for $-3 < \epsilon_{\parallel} < -2$. If a medium mimicking a homogeneous plasma is well-suited to fully determining the properties of an RF antenna, a complete characterization would require several loads to cover the range of useful densities. We present here a metamaterial load that is close to optimal coupling conditions ($\epsilon_{\parallel} = -3$) for LHRF waves.

We are led to materials whose permittivity is diagonal, in the form $\epsilon = \epsilon_0$ diag (1, 1, $\epsilon_{\parallel}$), where $\epsilon_{\parallel} < 0$. These materials are called hyperbolic materials, since waves polarized in the xy-plane (TE polarization) behave as in a vacuum, whereas TM waves obey a hyperbolic dispersion relation (their dispersion curve is a hyperboloid)⁵,⁶. In our case the field entering into the plasma from the antenna propagates in the xz-plane. See supplementary information for more details.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Equi-frequency dispersion curves at 3.7 GHz. (a) The gray hyperbola is the dispersion curve of a hyperbolic material with $ε_⊥ = 1$ and $ε_{||} = -3$. The colored curves are the dispersion curves of the grid-made metamaterial with periods $p_x = 4$ mm, $p_z = 5.1$ mm, and $p_y$ varying from 10 mm, 20 mm, 25 mm, 28 mm, 35 mm for blue, yellow, red, green and magenta curves respectively. The dashed black vertical lines correspond to the wired media. (b) Enlargement of the curve for $p_y = 28$ mm. The arrows show the average direction of energy propagation for $k_z = 2k_0$.</figcaption></figure>

One way to obtain a material with a diagonal permittivity $ε = ε_0$ diag(1, 1, $ε_{||}$) is to use an array of thin metallic wires that are parallel to the z-direction. Previous works$^{7,8}$ have shown that in such metamaterials $ε_{||}$ expresses using the Drude formula $ε_{||} = 1 - ω_p^2/ω_0^2$ with a plasma frequency $ω_p$ depending on the wire's dimensions and spacing, at least when the wavelength is much greater than the wire's spacing, when the wavevector has no component along the z-direction (i.e. $k_z = 0$), and when the electric field is parallel to the z-direction. It is shown in ref.$^9$ that these wire metamaterials have a strong spatial dispersion, i.e. $ω_p$ (and consequently $ε_{||}$ and the response of the medium) depends on the wavevector component along the wire direction.

Considering our case, it is clear from Fig. 1 that the antenna must emit waves with $k_z \approx 0$. In our configuration, we use $k_z = k_{||} \approx 2k_0 = 2ω_0/c$, which is a convenient value for current drive applications in a tokamak. To this end, we use a phased array of waveguides antennas. It means that the previous design with parallel wires cannot be used. Then, in order to mimic the plasma we need to reduce the spatial dispersion of the metamaterial. Several solutions have been proposed for this$^{10}$. We use the method of connecting the conducting wires of our metamaterial, which has been shown to limit the spatial dispersion$^9$.

We designed and tested a fishnet load composed of stacked thin metallic grid layers. The fishnet structure has the important practical advantage that it can be built layer by layer. This kind of construction was previously used by some of the authors for very low $ε_{||}$ medium to propagate a very low divergence beam in a vacuum$^{11,12}$, and more recently to control the density of states in hyperbolic metamaterials$^{13}$.

The fishnet metamaterial load is depicted in Fig. 1. It is made of a stack of thin copper grids. Each grid is parallel to the yz-plane with the following dimensions: 0.2 mm width of the conductors along y or z, 18 µm thickness along x, periods $p_y$ and $p_z$. The grids are printed on a 50 µm-thick polyimide film. They are spaced using Rohacell foam with low density (permittivity $ε = 1.04$) that gives a period $p_x$ along the x-axis.

The electromagnetic properties of the structure strongly depend on these dimensions. In order to obtain a structure with properties as close as possible to the expected hyperbolic material, we adjusted these parameters and tried to get dispersion curves as close as possible to those of a hyperbolic material with $ε_{||} = -3$. Due to the symmetry of the excitation versus y, the set of phased waveguides generates waves in the plasma that have a $k_y$ component centered in $k_y = 0$; the $k_x$ bandwidth can be reduced by increasing the dimension of the antenna along y, or by piling up several rows of waveguides. For this reason, we assume in the following that the spectrum of waves launched into the plasma has a negligible $k_y$ component. Since the excitation is at a given frequency, all the electromagnetic characteristics of the material related to wave propagation are contained in the equi-frequency dispersion curve for the components $k_x$ and $k_z$ of the wavevector.

Figure 2a shows the equi-frequency dispersion curves ($k_x$, $k_z$) for various dimensions of the fishnet cell. These curves are obtained by the eigenmode solver of CST Microwave Studio. To get these curves, we enforced the electric field to have no component along the y-axis. This is consistent with the way the material is excited (TE₁₀ mode inside the waveguides). The gray hyperbola is the dispersion curve of a homogeneous plasma with $ε_⊥ = ε_{xx} = ε_{yy} = 1$ and $ε_{||} = ε_{zz} = -3$. Note that its equation is $ε_{xx} k_x^2 + ε_{zz} k_z^2 = ε_{xx} ε_{zz} k_0^2$, i.e., $k_x^2 - 3 k_z^2 = -3 k_0^2$. The colored curves show the dispersion curves of the fishnet metamaterials for several values of the period $p_y$. The two black vertical dashed lines are obtained when the wires along the y-direction are suppressed (the only wires are along the z-direction). As mentioned before, even if the effective parameters for a null z-component of the wavevector are correct, wire media fail to reproduce the dispersion relation of the homogeneous plasma for larger $k_z$ values. For our application, we are interested in propagative waves with $k_x > k_0$, but no propagative waves exist inside the wire media for these wavevector values.

As shown in Fig. 2, for a fishnet metamaterial, the dispersion relation is not a perfect hyperbola. Thus, we had to make some choices: either choose parameters that give a curve as close as possible to the desired hyperbola in the range of $k_z$, or choose parameters that fulfill some other physical properties. We chose to build a metamaterial
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Numerical modeling of the modulus of the electric field at 3.7 GHz. The multi-junction antenna mouth is surrounded by a ground plane mimicking the tokamak vessel wall (horizontal line on these graphs). The field maps are shown in the middle of the antenna (y=0 according to Fig. 1). The multi-junction antenna is loaded with vacuum (a), an ideal plasma with $ε_⊥ = 1$ and $ε_{||} = -3$. (b), and the fishnet metamaterial (c). A high resolution version of this graph showing the entire multi-junction antenna is given in supplementary information.</figcaption></figure>

that will give a direction of energy propagation close to that of a homogeneous plasma. Note that the direction of energy propagation is normal to the dispersion curve¹⁴. Since the phased waveguide antenna launches waves with $n_{||} = k_z/k_0 \approx 2$, we retained the parameters associated with the green dispersion curve (Fig. 2b). This is a good compromise between a curve close to the hyperbola and one that gives for $k_z/k_0 \approx 2$ the same energy propagation direction as the homogeneous plasma (magenta arrows).

The built structure is shown in Fig. 1. The dimensions of the grid layers are 396 mm (along the z-direction) x 168 mm (along the y-direction). To reduce the parasitic reflections of the electromagnetic field at the boundaries of the fishnet load, we surrounded the load by a 30 mm thick RF absorber foam on five of its faces, the sixth being the one where the waveguides feed the load.

The antenna exciting the LHRF wave is composed of one input waveguide feeding 6 narrow phased waveguides stacked along the z-direction (see Fig. 3). This type of antenna, called a multi-junction antenna, allows reduction of the reflection coefficient in the input principal waveguide at the expense of an increase of the electric field in the narrow waveguides compared to a conventional antenna. The field emitted by this antenna has a spectrum centered on $n_{||} = k_z/k_0 \approx 2$. Consequently, there is no propagation inside the vacuum, a high reflection coefficient, and a high standing-wave ratio (SWR) inside the antenna, as shown in Fig. 3 left. Propagation inside the plasma is allowed according to Fig. 2(b). In this case, the SWR inside the antenna is low (Fig. 3, middle). The modeling depicted in Fig. 3 (right) shows that the fishnet metamaterial load that we designed preserves these properties: similar emitted field and low SWR. Moreover, the antenna emitted spectrum contains a secondary lobe centered on $n_{||} \approx -6$ (more details are given in supplementary). It can be seen in Fig. 2 that these waves can propagate in the ideal plasma (propagation direction on the left in Fig. 3b). In the case of the metamaterial, we can see from Fig. 2 that for $n_{||} = -6$ there is no existing propagative mode, as confirmed in Fig. 3c.

To quantify the ability of our structure to fit to the plasma behavior, we measured the reflection coefficient $S_{11}$ of the multi-junction antenna, and compared it (Fig. 4) to simulations using both CST Microwave Studio and the open-source code ALOHA¹⁵. CST is a full-wave frequency solver while ALOHA is a mode-matching code where the electromagnetic waves are fully absorbed in idealized semi-infinite plasma, which explains the small discrepancies between their results. We call $d_{vac}$ the gap distance between the antenna and the load. As expected $S_{11}$ is low when the metamaterial is close to the antenna and increases with $d_{vac}$. Considering that in this case the results given by ALOHA have an uncertainty of approximately $\pm 2$ dB, one can note that there is a good agreement between the measurements and the modeling.

The map of the x-component of the magnetic field propagating from the antenna to the medium was constructed from the measurements and compared to the simulations (Fig. 5). The main lobe of the launched spectrum ($n_{||} = 2$) is propagating along an angle $\alpha$ with the parallel direction (z-axis). The angle is given by $\alpha = \tan^{-1}\left[\frac{ε_{xx}(ε_{xx} - n_z^2)}{ε_{zz} n_z^2}\right]^{1/2} \approx 27$ deg, whereas the measurement gives $\alpha = 29$ deg. The radiation pattern of the antenna obtained by Fourier transform of the field at its mouth shows secondary $n_{||}$ lobes. The details are given in supplementary information. The highest amplitude secondary mode is for $n_{||} = -6$. This mode can be seen on Fig. 5c, and it propagates towards increasing z.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Reflection coefficient as a function of the gap between the antenna and the metamaterial load $d_{vac}$. We also give the value of $S_{11}$ when the antenna radiates in a vacuum. In the figure caption, “$\epsilon_{//} = -3$” stands for the homogeneous load with $\epsilon_{//} = -3$ and $\epsilon_{\perp} = 1$.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Magnetic field map with the metamaterial load (a) experiment, (b) modeling. The modeling of the equivalent plasma ($\epsilon_{//} = -3$) is shown in (c).</figcaption></figure>
## Discussion and Conclusion

We numerically and experimentally show that the metamaterial load we designed matches well the properties of a homogeneous plasma with equivalent hyperbolic dielectric tensor, resulting in low reflection or low insertion losses for a typical LHRF antenna with $n_{\parallel}=2$. Experimental measurements and numerical modeling of magnetic field maps show that waves propagate in the anisotropic medium ($\epsilon_{\parallel} = -3$, $\epsilon_{\perp} = 1$) along the expected directions. Thus, the metamaterial load will provide a drastic simplification for development of tokamak antennas and will help to make tokamak technology more efficient.

We envision many different versions and improvements. In a next step, two loads at 5 GHz with $\epsilon_{\parallel} = -10$ and $-1$ will be manufactured using the same scheme. These loads will be tested with a simple 4-waveguide antenna.

New concepts of antenna will also be tested with these loads. From a longer term perspective, several refinements can be proposed. For instance, a load mimicking a typical density gradient along the x-axis (radial) of a tokamak plasma near the antenna ($\Delta\epsilon_{\parallel}/\Delta x \sim -10/5$ cm) can be obtained by continuously varying the dimensions of the metallic grid across the 25 layers. Gyrotropic properties that we neglected, but could be important for some frequency ranges could also be simulated by specific metamaterials structures¹⁶,¹⁷. Moreover, active metamaterials could reproduce the dynamics and instabilities of the plasma¹⁸. Finally, such anisotropic structures can be used in optical domain to enhance spontaneous emission, negative refraction and superlensing effects⁵.

## Methods

The load is tested with a 6-waveguide multi-junction-type antenna module. The individual waveguides (dimensions 8 mm × 72 mm) are phased with built-in phase shifters ($\Delta\Phi = \pi/2$) in order to launch a wave with a parallel wave index centered on $n_{\parallel} = 2.0$. The height of the waveguides is such that all modes, except TE$_{10}$, are evanescent. Due to the limited number of waveguides the $n_{\parallel}$ spectrum is rather broad ($\Delta n_{\parallel}$)$_{mid-height} = \lambda_0/L_{\parallel,ant} \approx 1$. A detailed analysis of the emitted spectrum is presented in supplementary information and shows that in addition a secondary lobe exists for $n_{\parallel} = -6$. The distance between the antenna and the load, which is terminated by one half of a period $p_x/2$ of foam layer, needs to be accurately tuned as the wave ($n_{\parallel} > 1$) is evanescent in vacuum/air with an evanescent length in the mm range. However, because of the rather poor flatness and thickness accuracy of the layers of spacers and absorbers, the mean gap cannot be controlled with accuracy better than 0.5–1 mm.

The RF field in the load was measured in a plane at mid-height of the waveguides ($y=0$) where the amplitude of the TE$_{10}$ mode is the highest. The field was measured with a magnetic loop ($\phi=3$ mm) inserted inside the load through empty channels hollowed in the foam.

The experimental results are compared to those obtained by full wave modeling, using either an open-source code ALOHA including a plasma module¹⁵ or the commercial RF code CST Microwave Studio, which can also describe the wave propagation in homogeneous plasma.

## References

1. Ongena, J., Koch, R., Wolf, R. & Zohm, H. Magnetic-confinement fusion. *Nature Physics* **12**, 398–410 (2016).

2. Stevens, J., Ono, M., Horton, R. & Wilson, J. R. Edge density profile effects for lower hybrid waveguide coupling. *Nuclear Fusion* **21**(10), 1259 (1981).

3. Stix, T.H. *Waves In Plasmas*, AIP-Press, 1992.

4. Goniche, M. et al. SOL characterization and LH coupling measurements on JET in ITER-relevant conditions. *Plasma Phys. Control. Fusion* **51**, 044002 (2009).

5. Poddubny, A., Iorsh, I., Belov, P. & Kivshar, Y. Hyperbolic metamaterials. *Nature Photonics* **7**, 948–957 (2013).

6. Ferrari, L. et al. Hyperbolic metamaterials and their applications. *Progress in Quantum Electronics* **40**, 1–40 (2015).

7. Guida, G., Maystre, D., Tayeb, G. & Vincent, P. Mean-field theory of two-dimensional metallic photonic crystals. *J. Opt. Soc. Am. B* **15**, 2308–2315 (1998).

8. Pendry, J. B., Holden, A. J., Stewart, W. J. & Youngs, I. Extremely Low Frequency Plasmons in Metallic Mesostructures. *Phys. Rev. Lett.* **76**, 4773 (1996).

9. Belov, P. A. et al. Strong spatial dispersion in wire media in the very large wavelength limit. *Phys. Rev. B* **67**, 113103 (2003).

10. Demetriadou, A. & Pendry, J. B. Taming spatial dispersion in wire metamaterial. *J. Phys. Condens. Matter* **20**, 295222 (2008).

11. Enoch, S., Tayeb, G., Sabouroux, P., Guérin, N. & Vincent, P. A metamaterial for directive emission. *Phys. Rev. Lett.* **89**, 213902 (2002).

12. Ourir, A., Abdeddaim, R. & de Rosny, J. Planar metamaterial based on hybridization for directive emission. *Optics Express* **20**(16), 17545–17551 (2012).

13. Rustomji, K., Abdeddaim, R., de Sterke, C. M., Kuhlmey, B. & Enoch, S. Measurement and simulation of the polarization-dependent Purcell factor in a microwave fishnet metamaterial. *Phys. Rev. B* **95**, 035156 (2017).

14. Gralak, B., Enoch, S. & Tayeb, G. Anomalous refractive properties of photonic crystals. *J. Opt. Soc. Am. A* **17**, 1012–1020 (2000).

15. Hillairet, J. et al. ALOHA: an Advanced LOwer Hybrid Antenna coupling code. *Nucl. Fusion* **50**, 125010, https://github.com/jhillairet/ALOHA (2010).

16. Wang, Z. et al. Gyrotropic response in the absence of a bias field. *Proceedings of the National Academy of Sciences* **109**(33), 13194–13197 (2012).

17. Rogacheva, A. V., Fedotov, V. A., Schwanecke, A. S. & Zheludev, N. I. Giant gyrotropy due to electromagnetic-field coupling in a bilayered chiral structure. *Phys. Rev. Lett.* **97**(17), 177401 (2006).

18. Zheludev, N. I. & Kivshar, Y. S. From metamaterials to metadevices. *Nature Materials* **11**(11), 917–924 (2012).

## Acknowledgements

This work is partially supported by the Région Provence-Alpes-Côte d'Azur (Project DEB 13-1464). The project leading to this publication has received funding from Excellence Initiative of Aix-Marseille University - A*MIDEX, a French "Investissements d'Avenir" program. Research was conducted within the context of the International Associated Laboratory "ALPhA: Associated Laboratory for Photonics between France and Australia".

## Author Contributions

All authors contributed to the concept of the metamaterial load, to the analysis of numerical and experimental results and to the writing of the paper. K.R, R.A and J.H. performed numerical simulations. R.A., J.A., E.G., M.C. and M.G. conducted experiments and performed the data processing.
**Additional Information**

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-018-24250-0.

**Competing Interests:** The authors declare 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) 2018