Optically induced metal-to-dielectric transition in Epsilon-Near-Zero metamaterials

Abstract

Epsilon-Near-Zero materials exhibit a transition in the real part of the dielectric permittivity from positive to negative value as a function of wavelength. Here we study metal-dielectric layered metamaterials in the homogenised regime (each layer has strongly subwavelength thickness) with zero real part of the permittivity in the near-infrared region. By optically pumping the metamaterial we experimentally show that close to the Epsilon-Near-Zero (ENZ) wavelength the permittivity exhibits a marked transition from metallic (negative permittivity) to dielectric (positive permittivity) as a function of the optical power. Remarkably, this transition is linear as a function of pump power and occurs on time scales of the order of the 100 fs pump pulse that need not be tuned to a specific wavelength. The linearity of the permittivity increase allows us to express the response of the metamaterial in terms of a standard third order optical nonlinearity: this shows a clear inversion of the roles of the real and imaginary parts in crossing the ENZ wavelength, further supporting an optically induced change in the physical behaviour of the metamaterial.

Optically induced metal-to-dielectric transition in Epsilon-Near-Zero metamaterials

OPEN

R A

P

R. M. Kaipurath, M. Pietrzyk, L. Caspani, T. Roger, M. Clerici,, C. Rizza,, A. Ciattoni,A. Di Falco & D. Faccio

Epsilon-Near-Zero materials exhibit a transition in the real part of the dielectric permittivity from positive to negative value as a function of wavelength. Here we study metal-dielectric layered metamaterials in the homogenised regime (each layer has strongly subwavelength thickness) with zero real part of the permittivity in the near-infrared region. By optically pumping the metamaterial we experimentally show that close to the Epsilon-Near-Zero (ENZ) wavelength the permittivity exhibits a marked transition from metallic (negative permittivity) to dielectric (positive permittivity) as a function of the optical power. Remarkably, this transition is linear as a function of pump power and occurs on linearity of the permittivity increase allows us to express the response of the metamaterial in termsof a standard third order optical nonlinearity: this shows a clear inversion of the roles of the real and imaginary parts in crossing the ENZ wavelength, further supporting an optically induced change in the physical behaviour of the metamaterial.

Recent advances in metamaterial science have opened routes to unprecedented control over the optical properties of matter, with a wide array of applications and implications for novel light-matter interactions. Examples are the demonstration of negative index materials and, more recently, signicant attention has been devoted to the behaviour of light in a medium with zero dielectric permittivity. We will refer to these materials as Epsilon-Near-Zero (ENZ) materials with the implicit assumption that in all passive materials, the ENZ condition will only be met for the real part of the permittivity, , (as a result of absorption that will always imply that the imaginary part is greater than zero) and at one single wavelength (due to dispersion). Such ENZ materials may either occur naturally at the plasma frequency or may result from engineering the propagation medium, for example so that light propagates in a waveguide near the cuto frequency. Another option, investigated here, is to create a metamaterial (MM) made of deeply subwavelength alternating layers of dielectric and metal with thicknesses that are chosen such that is zero at a chosen wavelength. The linear properties of ENZ metamaterials have been investigated in depth with a range of applications for example in novel waveguiding regimes and for controlling the radiation pattern of electromagnetic sources118. The ENZ condition has also been predicted to have far-reaching consequences in terms of the eective optical nonlinearity of the metamaterial, but with limited experimental evidence1927. A compelling experimental evidence of the role of ENZ properties aecting the optical nonlinearity is the recent demonstration of efficient third harmonic generation due to the enhancement of the pump electric eld longitudinal component in a uniform lm of Indium-Tin-Oxide (ITO)28.

A dierent, yet related area of study, is the search for materials that can be optically controlled so as to exhibit a sharp and rapid transition from metallic to dielectric (or vice versa), thus implying a fundamental change in the material properties. Examples that have been investigated and observed in literature are optically induced phase transitions (for example in Vanadium Oxide and other compounds)29 or sudden increase in conductivity in glass when optically pumped with single cycle pulses, close to the breakdown damage threshold30,31. A

metal-to-dielectric transition has also been theoretically proposed in metallo-dielectric stacks, obtained by

School of Engineering and Physical Sciences, SUPA, Institute of Photonics and Quantum Sciences Heriot-Watt Dipartimento di

SCIENTIFIC REPORTS

1

www.nature.com/scientificreports/

Figure 1. Pictures of the sample and linear characterisation. (a) Top scanning electron microscope (SEM) image of the silver layer, showing the smooth surface obtained with the Ge wetting layer. (b) SEM image of the metamaterial sample showing the alternating layers of metal (Ag, 5nm) and dielectric (SiO2, 70nm) with 35nm top and bottom layers of SiO2. The laser exciting the material nonlinear response is directed along the surface normal and is polarised parallel to the surface. (c) The measured metamaterial linear dielectric permittivity (retrieved by tting the experimental transmission and reection spectra): real (, red dashed curve) and imaginary (, blue dashed curve) parts. The solid curves show the theoretical eective-medium predictions for and , as given by the formula

= + +

l l l l

( )/( )

eff m Ag d SiO m d

2 , where lm is the thickness of the Ag layer and ld is the thickness of the SiO2 layer41. = 0 for a wavelength of =885nm. (d) Same as in (c) with slightly tuned SiO2 thickness so that the ENZ wavelength is shied to 820nm.

optically pumping close to the ENZ wavelength27 with relatively fast ~1 ps response times. Here we experimentally investigate the optical behaviour of MMs made of deeply subwavelength alternating layers of fused silica glass and silver with thicknesses that are chosen such that = 0 in the near-infrared region, as shown in Fig.1. We show that by optically pumping the MM far from the ENZ wavelength, it is possible to induce a marked and rapid (on the same time scale of the 100 fs pump pulse) transition from metallic ( < 0) to dielectric ( > 0). By changing the layer thickness of the MM we can also control the wavelength at which this transition occurs. Remarkably we nd that the metal-dielectric transition occurs linearly as a function of the pump intensity. We show that this fact allows to describe the optical behaviour in terms of a standard third order nonlinear susceptibility. Our measurement technique provides the full complex amplitude across a wide spectral range centred at the ENZ wavelength.

Results

We fabricated the metamaterial samples by electron beam deposition of alternating layers of Ag and SiO2 on a thick (1 mm) SiO2 substrate with a total of 10 layers. The Ag thickness in each layer is kept at 5 nm whilst the SiO2 thickness is the same in each layer and tuned to 8070 nm in order to provide the ENZ condition around 820890 nm, in the centre of our laser tuning region. To obtain smooth continuous layers of silver on silica below the standard percolation limit, we seeded the deposition of each metal layer with 0.7nm of Germanium32.

Figure1(a,b) show a photograph of one of the samples and an SEM image of the multilayer structure, respectively. More details on the fabrication process are provided in the Methods section.

Linear response. The linear response (real and imaginary part of , and , respectively) was measured by a standard reection/transmission measurement (see Methods for details) and is shown as the dashed lines in Fig.1(c,d) for two dierent samples with SiO2 thickness equal to 80nm (referred to as sample A in the following) and 70nm (sample B), respectively. is measured to be zero at 885nm (sample A) and 820nm (sample B).

Under the approximation of deeply subwavelength lms, light polarised parallel to the lm does not interact with each individual layer of the multilayer structure but rather with an eective homogenised medium whose complex dielectric susceptibility is given by e=(ldd+lmm)/(ld+lm) where may represent either the linear

SCIENTIFIC REPORTS

2

www.nature.com/scientificreports/

Figure 2. Experimental results: nonlinear (pumped) reectivity as a function of the probe wavelength.

Pump and probe reectivity variation for the two metamaterial samples [sample A (a) and sample B in (b)] around the ENZ wavelength (indicated with a vertical dashed line). The pump wavelength is kept xed at 785nm whilst the probe wavelength is varied (horizontal axis). The pump intensity is 20GW/cm2 and 10GW/ cm2, for sample A and B, respectively. The data show a step-like behaviour in proximity to the ENZ wavelength.

susceptibility, (1), or also the third order nonlinear susceptibility, (3)33. l indicates the thickness of the individual layer and the subscripts d and m refer to the dielectric and metal layers, respectively.

The solid lines in Fig.1(c,d) show the predictions for a homogenised material and are in excellent agreement with the measured data between 750nm and 920nm, thus indicating that in this wavelength region the metamaterial is indeed behaving as an eectively uniform and homogenised medium.

Nonlinear response. We measured the nonlinear response of the metamaterial by monitoring the changes in reectivity, R = R Rlin, and transmissivity, T = T Tlin, in a pump and probe experiment (where Rlin and Tlin are the linear without the pump reectivity and transmissivity, respectively). The pump (with a xed wavelength of 785 nm, pulse duration 100 fs, horizontally polarised) is at normal incidence on the sample and the probe (wavelength tuned in the 7001000nm region, 100 fs pulse duration, 50 Hz repetition rate, vertically polarised) is incident at a small ~1 deg. angle with respect to the pump. The pump-probe delay was adjusted to maximise the nonlinear eect i.e. was zero within the precision of the pulse duration. The probe power is always kept extremely low so that alone it does not induce any nonlinear eects whilst the pump power is varied between zero (pump blocked) and ~30GW/cm2.

Two sets of measurements for pump intensities of 20 GW/cm2 and 10 GW/cm2 are shown in Fig.2(a,b) for samples A and B, respectively. We note that there is a clear step-like increase in the normalised R/Rlin as the probe is tuned across the ENZ wavelength. We did not observe any measurable dierence in transmissivity between the pumped and not pumped case (within the noise limit of our detectors, a few percent). This can be explained by a simultaneous change in absorption that eventually balance a variation in the light transmitted into the sample. We thus consider T=0.

While the sudden change in reectivity is clearer for sample A, we still observe a signicant variation also for sample B, yet with a broader response. As we will show in the following, dierent eects contribute to the change in reectivity R/Rlin, including the dispersion of the (3) and the transition through the ENZ wavelength. A more thorough analysis of the variation of the permittivity with the pump intensity and probe wavelength is therefore required to unveil the underlying processes. We thus extend the same method followed to extract the linear permittivity from Rlin and Tlin, to also extract in the presence of the pump: the values of the reectivity and transmissivity in the nonlinear (pumped) case allow to retrieve the nonlinear (pumped) value of the permittivity. In particular, here we are interested in the behaviour just above the ENZ wavelength: as can be seen in Fig.1, here the unpumped is negative. In the presence of a positive and sufficiently large increase in due to the optical pump we may predict that the permittivity will transition from below to above zero. Figure3 shows as a function of pump intensity for sample A and sample B, measured at 890nm and 825nm, respectively (in both cases, 5nm above the ENZ wavelength). As can be seen, in both cases the permittivity transitions from negative to positive, thus indicating a transition of the medium from metallic to dielectric. The total variation ~0.05 is of the same order of the absolute value of the permittivity itself, implying a relatively large bandwidth of ~10nm over which the optically-induced metal-dielectric transition occurs for the maximum pump power (limited by material damage). In the inset (a) to Fig.3 we also show the corresponding imaginary parts of the permittivity that also increase with pump intensity. Inset (b) shows as a function of the relative pump-probe delay measured on sample A for a pump power of 17GW/cm2: the rise time is of the order of the 100fs pump pulse duration (followed by a decay time of a few ps that is typical for Ag). We see that by tuning the delay, it is possibly to tune the precise value of , crossing from metal to dielectric and back again.

A notable feature of this data is the clear linear dependence of both and with pump power (in disagreement for example with the theoretical predictions of Husakou et al.27): the dashed lines in Fig.3 represent linear ts to the data, which are seen to pass through the value measured in the absence of the pump (and reported in Fig.1, as expected). We also note that this linear behaviour was observed over a wide range of wavelengths (700nm to 1000nm, data not shown). This feature is remarkable as it allows us to relate the behaviour of the MM and the transition from metal to dielectric in terms of a standard third order nonlinear susceptibility,

SCIENTIFIC REPORTS

3

www.nature.com/scientificreports/

Figure 3. Nonlinear (pumped) behaviour of the permittivity. (a) Real part of () versus input pump intensity for sample A (blue circles) and B (red squares). The dashed lines show linear ts to the experimental data. Panel (b) shows instead the imaginary party of (). (c) Variation in versus pump-probe delay for sample A measured at 17GW/cm2 pump intensity. pump=785nm, probe = 890 and 825nm, for sample A and B, respectively. The error bars have been evaluated as standard deviation over 10 samples.

(3)34. Indeed, because of this linear behaviour we can extract the nonlinear susceptibility tensor element as (see Methods for details):

=

n c I

I

p p p

0 pr

where pr and p are the probe and pump frequencies, respectively, and np is the real part of the medium refractive index at the pump frequency.

The complex values of the permittivity, (pr, Ip), were retrieved from the reectivity and transmissivity measurements at dierent probe wavelengths and pump intensities. We used the transfer matrix approach to determine the value of permittivity that results in the measured reectivity and transmissivity35 (see Methods for further details). The nonlinear susceptibility is then calculated from these values using Eq.(1).

We note that the linearity observed in the variation of with pump intensity (as shown in Fig.3) implies that Eq.(1) is consistent, as it provides us with (3) values that are constants (do not depend on Ip). We remark that this method diers from the simple retrieval of the permittivity in the pumped case, as it exploits the found linear behaviour of versus intensity to interpret the nonlinear mechanism in term of a third-order nonlinearity, also allowing to extract the complex value of the (3) tensor at dierent pump and probe wavelengths.

This derivation neglects the variation of the pump intensity inside the sample along the propagation direction (due to absorption). Averaging the intensity over the sample thickness would lead to a small correction factor ~2 for the values of (3), thus here we chose the simplied formulation.

We use Eq.(1), applied separately to the real and imaginary part of to plot the real and imaginary parts of (3). Figure4(a,c) show the real (solid blue line) and imaginary (dashed red line) third order nonlinear coefficients r(3) and i(3) for samples A and B, respectively. Figure4(b,d) show the same data but plotted as the absolute value |(3)| (solid blue line) and phase, (dashed red line). The notable feature of these results is that as the probe wavelength crosses from the dielectric-like region ( > 0) to the metallic-like region ( < 0), the |(3)| changes in nature from predominantly real to predominantly imaginary.

Discussion

The maximum measured |(3)| is of the same order or slightly larger than that measured by other means in bulk silver |(3)|=2.81019m2/V234,36. However, it is interesting to note that whereas silver cannot be used in transmission in samples thicker than ~40nm (all light is reected from the surface at these thicknesses), the metamaterial is signicantly thicker and more than two orders of magnitude thicker than the total thickness of silver present in the sample. This therefore provides a much longer eective interaction length for light with the high nonlinearity of the metal, as compared to the bare material. Note that the increased transmission is a known result for alternating layers of dielectric and metals, and it is for example exploited in one-dimensional metal-dielectric photonic crystal37,38.

We also underline that the large imaginary component of the (3) tensor measured in the metallic wavelength region does not imply larger nonlinear losses. Indeed, third order nonlinearities lead to nonlinear phase shifts and nonlinear absorption that are determined by the real and imaginary part of nonlinear refractive index, respectively (and not of the susceptibility). These are given by (in the degenerate case)39:

= +

n c n n D

3/(4 )[ ]/

2 r r r i i

0

SCIENTIFIC REPORTS

(3) ( , ) 3

( , ) ,

pr p

(1)

2 2 .

The imaginary part of n2 is usually associated to what is known as the nonlinear absorption coefficient, 2=4n2i/, where is the vacuum wavelength. We underline that n2r and 2 depend respectively, on the sums and dierences of the real and imaginary parts of of the (3) tensor. This implies for example that a large i(3) will

enhance Kerr eects such as phase modulation that are associated to the n2r coefficient whilst simultaneously minimising nonlinear absorption (the data in Fig.4 shows that actually 2~0 close to the ENZ wavelength).

In conclusion, ENZ metamaterials allow to tailor and access novel optical propagation regimes. The interplay between the linear and nonlinear, real and imaginary propagation constants leads to metamaterials that exhibit nonlinearities with amplitudes similar to those of metals yet spread over material thicknesses two orders of

(3) (3) where = +

D n n n

( )

r r i

(3) (3) and

0

=

n c n n D

3/(4 )[ ]/

2 i r r i i

4

www.nature.com/scientificreports/

Figure 4. Experimental values of the third-order nonlinear susceptibility. Real and imaginary parts of the third order nonlinearity r(3) (solid blue line) and i(3) (dashed red line), respectively for sample A (a) and

sample B (c). The same data is also plotted as the absolute value |(3)| (solid blue line) and phase (dashed red line) for sample A in (b) and for sample B in (d). The vertical dashed black line represents the ENZ wavelength. The maximum pump intensity used for extracting the (3) values was 21.5GW/cm2 for sample A and 14.6GW/ cm2 for sample B.

magnitude larger and with substantially reduced losses. We have also shown that it is possible to optically induce a metal-dielectric transition over a wide bandwidth that can support ~100fs laser pulses. Further optimisation by e.g. reducing losses (and hence the linear refractive index) by introducing a gain medium25 may enable new forms of efficient switches for light and even a platform for non-perturbative nonlinear optics at low light intensities.

Finally we would like to highlight that our results are not specic to the sample used in our experiments, rather they can be generally applied to all ENZ media. Indeed, for these materials, close to the ENZ wavelength we can expect to observe a region where r is slightly negative, and provided that the nonlinear response is fast enough and the (3) is positive, we might expect similar results (alternatively for materials with negative (3), similar results can be obtained at a wavelength displaying a slightly positive real permittivity). For example, non-structured materials exhibiting an accessible ENZ region in the near infrared, such as transparent conductive oxides (TCOs), are expected to display a similar behaviour. In aluminium-doped zinc oxide (AZO) thin lms ultrafast switching close to the ENZ wavelength at 1300 nm has been recently observed40, and this could be a suitable candidate to observe a similar dynamics.

Methods

Fabrication. The ENZ metamaterial was fabricated on 1 mm thick slides of SiO2. The substrate was thoroughly cleaned in ultrasonic assisted baths of Acetone and Isopropanol (5mins each) and blown dried with N2 ow. The metallo-dielectric stack was deposited using a EDWARD auto 306 electron beam evaporator, with a base pressure below 3106bar. The deposition rate was kept below 0.1nm/s for the metals and below 0.3nm/s for the SiO2 to grant uniformity. The thickness of the seeding Ge layer (0.7 nm) was chosen aer a thorough optimisation process, to minimise the losses of the thinnest achievable layer of Silver (5 nm). An excess of 50 independent evaporations was performed to determine the nal values. The development process was complemented by a combination of SEM, STEM and AFM measurements on sacricial test samples.

Linear characterisation. To characterise the linear properties of the MM we built a simple setup for the measurement of the reection and the transmission. An Ocean Optics HL2000 Halogen Lamp covering the visible-NIR spectrum was collimated with a telescope to a beam radius of mm

5 , with controlled polarisation.

SCIENTIFIC REPORTS

5

www.nature.com/scientificreports/

A thin lm beam splitter was used to separate the reected light from the incident one. The two beams were then focused on the tip of multimode bres and analysed with two Ocean Optics spectrum analysers. Light was impinging on the sample on the side of the multilayer. The reference for the reection was a silver mirror, for the transmission we recorded the collected light without any sample. The raw data were then used to retrieve the eective complex permittivity of the multilayer, using a standard least mean square tting procedure: the experimental reected and transmitted value vs where compared with those calculated with a transmission matrix simulation, using a test value for the permittivity. Using the same method we also characterised the linear permit-tivity of SiO2 and that of a single Ge/Ag bilayer, which were then used to calculate the eective index permittivity of the homogenised stack.

Nonlinear characterisation. The material polarisation at the probe wavelength is given by

=

P I E

( ) [ ( , ) 1] ( ), (2)

pr 0 pr pr

where pr and p are the probe and pump frequencies, respectively, E(pr) is the probe electric eld, and the relative nonlinear permittivity is dened as

p

= + + .

I In c

( , ) 1 ( ) 6 ( , ) 2 (3)

p

being np the real part of the medium refractive index at the pump frequency. It is thus clear that a linear scaling of the permittivity with the pump intensity allows an interpretation in terms of a third-order nonlinearity, where the the (3) tensor can be retrieved by deriving the relative permittivity (pr, Ip) with respect to Ip:

=

n c I

I

p p p

0 pr

The nonlinear (pumped) permittivity (pr, Ip) has been evaluated from the reectivity and transmissivity measurements in presence of the pump, for dierent pump intensities and probe wavelength. Starting from the linear (pr, 0) and exploiting a transfer matrix approach similar to the one used for the linear characterisation, we found the value (pr, Ip) resulting in a R/Rlin and T/Tlin that best matched the experimental data at the given Ip35. The dual condition on reectivity and transmissivity allows one to retrieve both the real and imaginary part of (pr, Ip). We then used Eq.(4), identical to Eq.(1) in the main text, to calculate (3)(pr, p).

Data availability. All data relevant to this work may be obtained at doi: 10.17861/37615dc3-8a2d-4242-bdf8-3d06580e3571.

References

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

2. Edwards, B., Al, A., Young, M. E., Silveirinha, M. & Engheta, N. Experimental verication of epsilon-near-zero metamaterial coupling and energy squeezing using a microwave waveguide. Phys. Rev. Lett. 100, 033903 (2008).

3. Liu, R. et al. Experimental demonstration of electromagnetic tunneling through an epsilon-near-zero metamaterial at microwave frequencies. Phys. Rev. Lett. 100, 023903 (2008).

4. Ziolkowski, R. W. Propagation in and scattering from a matched metamaterial having a zero index of refraction. Phys. Rev. E 70, 046608 (2004).

5. Zhou, H. et al. A high-directive patch antenna based on all-dielectric near-zero-index metamaterial superstrates. J. Electromagn. Waves Appl. 24, 13871396 (2010).

6. Yang, J., Huang, M. & Peng, J. Directive emission obtained by Mu and epsilon-near-zero metamaterials. Radioengineering 18, 124128 (2009).

7. Soric, J. C., Engheta, N., Maci, S. & Al, A. Omnidirectional metamaterial antennas based on e-near-zero channel matching. IEEE Trans. Antennas Propag. 61, 3344 (2013).

8. Jin, Y. & He, S. Enhancing and suppressing radiation with some permeability-near-zero structures. Opt. Express 18, 1658716593 (2010).

9. Silveirinha, M. & Engheta, N. Tunneling of electromagnetic energy through subwavelength channels and bends using -near-zero materials. Phys. Rev. Lett. 97, 157403 (2006).

10. Silveirinha, M. & Engheta, N. Design of matched zero-index metamaterials using nonmagnetic inclusions in epsilon-near-zero media. Phys. Rev. B 75, 075119 (2007).

11. Pan, Y. & Xu, S. Energy tunnelling through an ultrasmall epsilon-near-zero channel in circular waveguide. IET Microwaves Antennas Propag. 3, 821825 (2009).

12. Cheng, Q., Liu, R., Huang, D., Cui, T. J. & Smith, D. R. Circuit verication of tunneling eect in zero permittivity medium. Applied Physics Letters 91 (2007).

13. Al, A. & Engheta, N. Coaxial-to-waveguide matching with -near-zero ultranarrow channels and bends. IEEE Trans. Antennas Propag. 58, 328339 (2010).

14. Nguyen, V. C., Chen, L. & Halterman, K. Total transmission and total reection by zero index metamaterials with defects. Phys. Rev. Lett. 105, 233908 (2010).

15. Feng, S. Loss-induced omnidirectional bending to the normal in -near-zero metamaterials. Phys. Rev. Lett. 108, 193904 (2012).16. Fleury, R. & Al, A. Enhanced superradiance in epsilon-near-zero plasmonic channels. Phys. Rev. B 87, 201101 (2013).17. Maas, R., Parsons, J., Engheta, N. & Polman, A. Experimental realization of an epsilon-near-zero metamaterial at visible wavelengths. Nat. Photonics 7, 907912 (2013).

18. Kyoung, J. et al. Epsilon-near-zero meta-lens for high resolution wide-eld imaging. Opt. Express 22, 3187531883 (2014).19. Argyropoulos, C., Chen, P., DAguanno, N., Engheta, G. & Al, A. Boosting optical nonlinearities in -near-zero plasmonic channels. Phys. Rev. B 185, 045129 (2012).

(1)

pr

pr

(3) pr p

p p

0

(3) ( , ) 3

( , ) ,

pr p

(4)

SCIENTIFIC REPORTS

6

www.nature.com/scientificreports/

20. Ciattoni, A., Rizza, C. & Palange, E. Extreme nonlinear electrodynamics in metamaterials with very small linear dielectric permittivity. Phys. Rev. A 81, 043839 (2010).

21. Vincenti, M. A., de Ceglia, D., Ciattoni, A. & Scalora, M. Singularity-driven second- and third-harmonic generation at -near-zero crossing points. Phys. Rev. A 84, 063826 (2011).

22. Rizza, C., Ciattoni, A. & Palange, E. Two-peaked and at-top perfect bright solitons in nonlinear metamaterials with epsilon near zero. Phys. Rev. A 83, 053805 (2011).

23. Ciattoni, A., Rizza, C. & Palange, E. Transverse power ow reversing of guided waves in extreme nonlinear metamaterials. Opt. Express 18, 1191111916 (2010).

24. Kauranen, M. & Zayats, A. Nonlinear plasmonics. Nat. Photon. 6, 737 (2012).25. Rizza, C., Di Falco, A. & Ciattoni, A. Gain assisted nanocomposite multilayers with near zero permittivity modulus at visible frequencies. Applied Physics Letters 99, 221107 (2011).

26. Suchowski, H. et al. Phase Mismatch-Free Nonlinear Propagation in Optical Zero-Index Materials. Science 342, 12231226 (2013).27. Husakou, A. & Hermann, J. Steplike transmission of light through a metal-dielectric multilayer structure due to an intensity-dependent sign of the eective dielectric constant. Phys. Rev. Lett. 99, 127402 (2007).

28. Capretti, A., Wang, Y., Engheta, N. & Dal Negro, L. Enhanced third-harmonic generation in Si-compatible epsilon-near-zero indium tin oxide nanolayers. Opt. Lett. 40, 15001503 (2015).

29. Lysenko, S., Rua, A., Fernandez, F. & Liu, H. Optical nonlinearity and structural dynamics of vo2 lms. J. Appl. Phys. 105, 043502 (2009).

30. Schultze, M. et al. Controlling dielectrics with the electric eld of light. Nature 493, 75 (2013).31. Schirin, A. et al. Optical-eld-induced current in dielectrics. Nature 493, 70 (2013).32. Chen, W., Thoreson, M. D., Ishii, S., Kildishev, A. V. & Shalaev, V. M. Ultra-thin ultra-smooth and low-loss silver films on a germanium wetting layer. Optics Express 18, 51245134 (2010).

33. Sipe, J. & Boyd, R. Nonlinear susceptibility of composite optical materials in the maxwell-garnett model. Phys. Rev. A 46, 1614 (1992).

34. Boyd, R. Nonlinear Optics (Academic Press, New York, 2008).35. Owens, D. T., Fuentes-Hernandez, C., Hales, J. M., Perry, J. W. & Kippelen, B. A comprehensive analysis of the contributions to the nonlinear optical properties of thin ag lm. J. Appl. Phys. 107, 123114 (2010).

36. Bloembergen, N., Burns, W. & Matsuoka, M. Reected third harmonic generated by picosecond laser pulses. Opt. Commun. 1, 195198 (1969).

37. Bennink, R. S., Yoon, Y.-K., Boyd, R. W. & Sipe, J. E. Accessing the optical nonlinearity of metals with metaldielectric photonic bandgap structures. Opt. Lett. 24, 14161418 (1999).

38. Lepeshkin, N., Schweinsberg, A., Piredda, G., Bennink, R. & Boyd, R. Enhanced nonlinear optical response of one-dimensional metal-dielectric photonic crystals. Phys. Rev. Lett. 93, 123902 (2004).

39. Smith, D. D. et al. Z-scan measurement of the nonlinear absorption of a thin gold lm. J. Appl. Phys. 86, 6200 (1999).40. Kinsey, N. et al. Epsilon-near-zero Al-doped ZnO for ultrafast switching at telecom wavelengths. Optica 2, 616622 (2015).41. Garnett, J. M. Enhanced nonlinear optical response of composite materials. Philos. Trans. R. Soc. London 203, 385 (1904).

Acknowledgements

This work was supported by the EPSRC grant EP/J004200/1. D.F. acknowledges nancial support from the European Research Council under the European Union Seventh Framework Programme (FP/20072013)/ ERC GA 306559 and EPSRC (UK, Grant No. EP/J00443X/1). L.C. and M.C. acknowledge the support from the People Programme (Marie Curie Actions) of the European Unions FP7 Programme THREEPLE (GA 627478) and KOHERENT (GA 299522). A.C. and C.R. acknowledge support from U.S. Army International Technology Center Atlantic for nancial support (Grant No. W911NF-14-1-0315). The authors acknowledge discussions with M. Scalora.

Author Contributions

R.M.K. performed experiments with assistance from T.R., M.C. and M.P. fabricated and characterised linearly the samples, L.C. and M.P. performed data and theoretical analysis, C.R. and A.C. provided theoretical support, D.F. and A.D.F. directed the project and wrote the manuscript. All authors contributed to scientic discussions and to the manuscript.

Additional Information

Competing nancial interests: The authors declare no competing nancial interests.

How to cite this article: Kaipurath, R. M. et al. Optically induced metal-to-dielectric transition in Epsilon-Near-Zero metamaterials. Sci. Rep. 6, 27700; doi: 10.1038/srep27700 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

SCIENTIFIC REPORTS

7