Interferometric control of the absorption in optical patch antennas

Interferometric control of the absorption in optical patch antennas www.nature.com/scientificreports OPEN Interferometric control of the absorption in optical patch antennas Received: 2 February 2017 Caroline Lemaître, Emmanuel Centeno & Antoine Moreau Accepted: 21 April 2017 Optical patch nano-antennas possess unique absorption, field enhancement and concentration Published: xx xx xxxx capabilities – but their crosssection, as well as their response outside of normal incidence are not well understood. Here we explain the large cross-section by considering that each patch nanoantenna is a cavity excited from both sides. Such a simple physical picture allows to fully understand the influence of the angle of incidence – that odd resonances have a very high absorption cross-section which decreases when the incidence angle increases, while even resonances cannot be excited in normal incidence. A direct application would be to use these structures as an optical nanometric set-square. Gap-plasmon resonators are a new kind of optical metallic resonator , that has emerged as particularly prom- 2–6 7–9 ising , from a fundamental point of view as well as for all the applications that could be based on these 2, 10 structures given their extraordinary efficiency at concentrating light . It has actually been suggested that they 11, 12 13–15 could be used as cheap but very efficient selective absorbers , as sensors , as well as for designing reflection metasurfaces providing an almost complete control of light . Finally, the light confinement associated to their resonances is so high that an unprecedentedly high Purcell effect can be obtained using this kind of optical patch 17–20 antennas . These resonators are literally cavities for a peculiar plasmonic mode called a gap-plasmon which presents a very large effective index. This property alone explains how deeply sub-wavelength the resonators can be made, which in turn partly explains their efficiency: the absorption cross section of a single patch can be actually as high as 30 times its geometrical cross-section . The modes that can be excited in these cavities have already been 2, 21, 22 23 studied , and their importance largely underlined . However, the influence of the incidence angle, despite its 14, 24 25 practical importance for sensing for instance, is not well understood - and their extraordinary absorption cross-section hasn’t yet received a full explanation. Here we show that the influence of the incidence angle can be accurately taken into account if a single gap-plasmon resonator is modeled as a cavity excited from both ends, making the structure a clear example of 26 27 interferometrically controlled absorption , or of a coherent (perfect) absorber . We give analytic formula for the losses induced by the resonances inside the cavity and link it to the incidence angle. Such a description allows to understand why the absorption decreases when the incidence angle increases and explains the extraordinary 9, 28 absorption cross-section of the patch antenna compared to other kinds of gap-plasmon resonators . An analytic model for the absorption of a gap-plasmon resonator A selectively absorbing metasurface is typically constituted of periodically or randomly distributed patches separated from a metallic film by a thin dielectric material, as shown on Fig.  1. Here we will consider silver for all the metallic parts . The efficiency of the device benefits from the capacity of each patch to funnel light into the gap, so that when the resonators are close enough to each other, they are able to absorb all the incoming light and constitute a perfect absorber . Most of the structures that are fabricated consist in rectangular patches coupled to a film instead of the nanorods shown Fig.  1, but is largely admitted now that the physical behavior of both struc- tures is the same and that the die ff rence, for the fundamental modes, essentially reduces to a shift of the spectral 11 22 position of the resonances . For higher order modes, although they may be more complex , the fact that some are excited only for non-normal incidence is perfectly in agreement with what we expose here for 2D structures. We must underline that when the structure is periodic, surface plasmons are oen ex ft cited , and can be respon- sible for a very large absorption in very peculiar conditions. Such a mechanism is not present in disordered structures, where the surface plasmons are excited and absorb a non-negligible part of the energy but play a Université Clermont Auvergne, CNRS, Institut Pascal, F-63000, Clermont-Ferrand, France. Correspondence and requests for materials should be addressed to A.M. (email: antoine.moreau@uca.fr) Scientific Repo R ts | 7: 2941 | DOI:10.1038/s41598-017-03064-6 1 www.nature.com/scientificreports/ Figure 1. A 3D rendering of a selectively absorbing plasmonic metasurface based on gap-plasmon resonators: regularly spaced metallic nanorods coupled to a metallic substrate and separated from it by a nanometric dielectric layer. lesser role than when periodicity enhances their excitation. Here we will concentrate on the resonance of a single resonator, which is angularly and spectrally broader, is more efficient, does not depend on the periodicity and explains the response of disordered structures . In the thin layer sandwiched between the patch and the metallic film, light propagates under the form of a gap-plasmon along the x axis, the only propagating mode of a metal-dielectric-metal waveguide. This mode exists for p polarization only, with a magnetic field along the y axis. Solving Maxwell’s equations a looking for a guided mode with a pulsation ω and a wavevector k in a metal-dielectric-metal configuration actually yields a magnetic field that can be written ik () xt −ω Hx (, zP )( = ze ) y (1) with a vertical profile of the magnetic field  −κ z ez ∀> hd    Pz () = Cz cosh κ ∀∈ z − ,      κ z d ez ∀< − 2 (2)   −κ where to ensure the continuity of the magnetic field at all interfaces. Figure  2(a) shows the profile of C = κ d cosh the magnetic field associated to the gap-plasmon propagating in a nano-gap. e Th dispersion relation of the gap-plasmon is obtained using the continuity of the electric field at any inter - face and can be written   κ d κ 2 1 tanh κ += 0  2   ε 2 ε (3) dm 2 2 2 2 where κ =− αε k and κ =− αε k , α being the propagation constant in the x direction. Once the 1 m 0 2 d 0 pulsation ω is fixed, the dispersion relation allows to retrieve the complex propagation constant α . The main characteristics of the gap-plasmon is that its effective index, defined as n = R() α /k , diverges when the size of eff 0 the gap tends to zero, as shown Fig. 2(b). This plasmonic ee ff ct occurs when light propagates more in the metal than in the dielectric. This is what allows gap-plasmon resonators to be so small compared to the wavelength in vacuum. When the patch ends, the gap-plasmon can in fact no longer propagate. Instead, it is able to excite the surface plasmon and a waves propagating in free space , but as both present very low propagation constants, the outside space represents a very high impedance mismatch for the gap-plasmon. This is of course all the more so that the gap is small and the effective index of the gap-plasmon large. As a consequence, the reflection coefficient r of the 5, 11, 24 gap-plasmon when it encounters the edge of the patch can be very high . Computing this reflection coefficient is not an easy task. Its phase can be deduced from the location of the 5 24, 31 resonances themselves , and more direct determinations can be difficult too . Finally, it can be noticed that this can be done directly using a Fourier Modal Method , which relies on the computation of the modes that prop- agate in different vertically invariant layers to solve Maxwell’s equations and of the transmission and reflection coefficients. It directly gives access to r and α. It is even possible, relying on transformation optics techniques to introduce Perfect Matching Layers in order to get rid of the periodicity. Figure 2(c) shows the modulus of the computed reflection coefficient as a function of the waveguide width, computed using such techniques (see Methods). The phase of the reflection coefficient, which is critical for the resonance condition and thus to explain why the size of a resonator can be smaller than half of the gap-plasmon wavelength , is shown on Fig. 2(c). This phase increases when the the gap closes. Scientific Repo R ts | 7: 2941 | DOI:10.1038/s41598-017-03064-6 2 www.nature.com/scientificreports/ Figure 2. (a) Profile of the magnetic field associated with the gap-plasmon propagating along the x axis. (b) Effective index of the gap-plasmon as a function of the gap width d. (c) Modulus of the reflection R() α /k coefficient at the edge of the cavity as a function of the gap width d. Figure 3. Cavity model. (a) Sketch representing the cavity model including the reflection coefficients of the gap-plasmon on the edges of the patch. (b) Map of the magnetic field (modulus) for a single 75 nm wide square patch, separated from the metallic film by a 10 nm thick dielectric layer with an optical index of 1.54 illuminated from above in normal incidence at resonance. The field profile in the region which is enclosed in the white rectangle is accurately predicted by the model (center). The inset shows the map of the magnetic field (modulus) as predicted by the cavity model when both excitation coefficients A and A are taken equal, at resonance. l r Under the patch, gap-plasmons thus propagate from right to left and from left to right and are reflected on both sides of the patch, so that they form a stationary wave. This is the actual description of a cavity of size h , as shown on Fig. 3(a). This allows us to introduce a model for the gap-plasmon resonator, based on our knowledge of r and α, that are computed as explained above. Since a gap-plasmon propagates along x towards the left and another one propagates towards the right, the amplitude of the mode, as a function of x under the patch can simply be written ix αα −ix () xA =+ eB e . (4) This means that the magnetic field associated to the stationary wave can simply be written H (, xz )( =  xP )(z ). We choose here the edges of the cavity to be in ± . We call A and A the amplitude with which the l r gap-plasmons are excited respectively on the left and on the right of the patch. We know actually very little about these quantities, except for the fact that in normal incidence, they are obviously equal: the cavity is being excited on both sides exactly the same way. In non-normal incidence, the phase difference between the two can be esti- mated, as the following will show, but there must be a difference in the amplitude due to the shading by the cube that is difficult to take into account. These excitations should be considered as relatively large because when a plane wave is reflected by a metallic film, the magnetic field presents a maximum on the surface. Whatever the thickness of the dielectric layer, there is always a maximum of the magnetic field at the entrance of the cavity. This explains why the cavity resonances are much more easily excited in this configuration than for dimer , for instance. The resonances are, conversely, not likely to be excited through the patches – because they are twice as large as the skin depth. This simple cavity model, detailed in the Supplementary Information, thus gives the amplitude of the modes in the cavity as hh hh  −− ii ααxiαα ix ii ααxi −− αα ix  22 22 ee + re e re ee + e   () xA = + A . l r   22ih α 22ih α 11 − re − re   (5)   Scientific Repo R ts | 7: 2941 | DOI:10.1038/s41598-017-03064-6 3 www.nature.com/scientificreports/ 2 2iαh As shown by this expression, a resonance occurs each time the quantity 1 − r e presents a minimum and can −ih α −ih α be considered close enough to zero. The resonance condition, in that case, reduces to re  ± . When re  then the resonance is symmetric with respect to the x = 0 axis, which means it has an odd number of anti-nodes d d −ih α between − and . On the contrary, when re  − the magnetic field is anti-symmetric with respect to x = 0, 2 2 and the resonance presents an even number of anti-nodes. Keeping that in mind for later, it is then possible to get analytic expressions of all the fields, magnetic as well as electric, that are linked to the cavity and the stationary wave that exists in it - even out of resonance. When the amplitude () x is multiplied by the vertical profile P(Z) of the gap-plasmon, then the magnetic field is almost perfectly reproduced, as shown Fig. 3 where a comparison to a full simulation H (, xz )( =  xP )(z ) using a Fourier Modal Method of a single patch is shown. Using Maxwell’s equations, the electric field can be analytically calculated (See Supplementary Information). Finally, the volumetric losses p(r) inside the metal are classically written p() rE =. ωε I() ε E . 0 r (6) As shown in the Supplementary Information, the volumetric losses can be integrated over the metallic region, to compute analytically the absorption under one patch antenna, as a function of A and A l r h/2 +∞ Pp =. () r dd xz ∫∫ −− h/2 ∞ Discussion Two different things can be done starting from here. Either use the complex expression (given in the Supplementary Information) as it is, to compute the losses in every possible situation (out of resonance, for any angle and wavelength), or to compute a simplified analytic solution that brings a better understanding but only −ih α close to resonance. Let us first assume that we are close to a resonance, with re  ± . Then, as shown in the Supplementary Information, the absorption by the resonator can be written 22 ⁎ PA =+ β(2 AA ± R() A ) (7) lr lr In normal incidence, A = A since both ends of the cavity are excited in the same way. This explains why res- l r iαh onances with an even number of anti-nodes, for which r = −e cannot be excited in normal incidence, since we have strictly P = 0. On the opposite, for a resonance with an odd number of anti-nodes P = 4β|A | and the absorption is thus four times higher than if the cavity were excited from one side only. This shows the mechanism that is at the heart of the extraordinary large absorption cross-section that patch antenna have, when compared to their actual size. The interferometric reinforcement of the absorption, associated with the high effective index of the gap-plasmon which drives the reduction in size of the resonators, completely explains this effect. Now it is possible to take into account the incidence angle, as long as it is not too large by simply assuming that the modulus of the excitations A and A stay close to their maximum, that is reached in normal incidence, r l but that the main factor driving their variation will be a slight phase shift arising from the different optical paths followed by light to reach the entrances of the gap. With a incidence angle θ, this phase shift is Δφ = k h sin θ. We finally have an absorption at resonance that can be written either   π h   PA = 2c β os sin(θ)  λ   (8) for odd resonances (i.e. with an odd number of anti-nodes) and π h 22  PA = 2s β in sin(θ) l   λ   0  (9) for even resonances (i.e. with an even number of anti-nodes). Changing the angle of incidence thus lowers the absorption of odd resonances, while allowing the even resonances to be excited and to appear in the spectrum. Our findings can be summarized by a simple scheme (see Fig.  4). Odd (resp. even) resonances are actually symmetric (resp. anti-symmetric). Since each side of the patch is excited by the incoming field, the field under the patch can be seen as a superposition of two modes: one excited from the left and one excited from the right. In normal incidence, both excitations are perfectly in phase. When the resonance should be odd (and symmetric), then the two modes add up so that the amplitude of the field is exactly twice what it would be if the resonance had been excited from one side only (and the absorption is multiplied by four). When the resonances should be even (and anti-symmetric) the two modes interfere destructively, and the resonance is completely canceled. For oblique incidence, the two excitations are not in phase anymore, disturbing the interference between the modes. Odd resonances see their efficiency decrease, while even resonances can be excited. Now we can use the model described above, without any assumption on r, out of the resonance condition to try to predict, as the incidence angle increases, how the absorption linked to the resonance will evolve. Using the general expression for the absorption obtained thanks to our model, we plot the absorption as a function of the patch width (see Fig. 5) and compare it to the reflectance of the structure for different incidence angles . The accu- racy of the model can be clearly seen. It is almost perfect in normal incidence and for large angles the variation of Scientific Repo R ts | 7: 2941 | DOI:10.1038/s41598-017-03064-6 4 www.nature.com/scientificreports/ Figure 4. e in Th terferometric control of the absorption in patch antennas summarized in normal incidence: adding two anti-symmetric modes excited from left and right leads to a complete cancellation, while adding even modes leads to an enhancement. For oblique incidence, this control is lowered, so that the even modes don’t completely cancel each other, and the odd modes are less enhanced. Figure 5. Absorption spectrum for a silver nanopatch of variable width (in nm) from the model (left) and from the reference code (right). The layer is 5 nm thick, the wavelength 950 nm. the interferometric effect is not perfectly reproduced - probably because our assumptions regarding the optical path and the shading are not fully accurate. However, this definitely validates our physical interpretation of the nanopatch optical response. We underline that the destructive interference that literally kills odd resonances in normal incidence can be utilized to assess whether the angle of incidence is actually zero or not. As soon as the incidence angle differs from zero, the interferometric control will stop and the absorption linked to the resonance will sharply increase, as shown on the spectrum. This would be the optical equivalent of an optical set square at a nanometric level. Conclusion We have shown that a optical patch antenna can be considered as a Fabry-Perot cavity excited from both sides and that this explains why the absorption by even resonances is reinforced by an interferometric control effect while odd resonances are completely tamed by the same effect. Previous authors have already evoked the symmetry 24, 31 mismatch between the incident wave and the resonances , but the interferometric control pictures offers a way to quantify this effect. Using an analytic model whose parameters are computed numerically, we have been able to show how this effect varies when the incidence angle is changed, allowing odd resonances to be excited and slightly lowering the absorption efficiency of even resonances. This gives a better understanding of why the absorption of metasurfaces made with such resonators decreases when the angle of incidence increases - and shows that modifying the angle of incidence is an actual way to control the absorption cross section of each individual resonator and thus light concentration of the fluorescence enhancement for instance . We think that having in mind this fundamental mechanism of the absorption by optical patch antennas will help find new applications or new techniques. It may help to forecast the radiation pattern of film-coupled nanocubes , and shows that the angle of incidence is a way to control the absorption cross-section of a single patch. Finally, the fact that the excitation of odd resonances is strictly impossible at normal incidence could well be utilized to fabricate an optical set-square, giving an immediate application to the interferometric control of the absorption by optical patch antennas. Methods We used a Fourier Modal Method with coordinate stretching and Perfect Matching Layers to first solve the modal problem in a metal-insulator-metal structure. The period is of 930 nm with a 200 nm thick PML and a stretching parameter η = 0.99. On one side, we have 120 nm of metal and 600 nm on the other side. We take the mode with the smallest imaginary part for the propagation constant. This is the gap plasmon. We then compute the scattering matrix when the above structure is put next to a structure with the same parameters, except that the 600 nm metallic part is replaced by air. The rest of the parameters are kept the same. The mode with the lowest imaginary part for its propagation constant is then the surface plasmon. Using the boundary conditions allows to access the scattering matrix of the problem and then retrieve the coefficient of the scattering matrix corresponding to an incoming and outgoing gap-plasmon. The coefficient constitutes the reflection coefficient of the gap-plasmon on each edge of the structure. Scientific Repo R ts | 7: 2941 | DOI:10.1038/s41598-017-03064-6 5 www.nature.com/scientificreports/ References 1. Bozhevolnyi, S. & Sondergaard, T. General properties of slow-plasmon resonant nanostructures: nano-antennas and resonators. Optics Express 15, 10869 (2007). 2. Søndergaard, T. & Bozhevolnyi, S. Slow-plasmon resonant nanostructures: Scattering and field enhancements. Phys . Rev. B 75, 073402, doi:10.1103/PhysRevB.75.073402 (2007). 3. Jung, J., Søndergaard, T. & Bozhevolnyi, S. I. Gap plasmon-polariton nanoresonators: Scattering enhancement and launching of surface plasmon polaritons. Phys. Rev. B 79, 035401, doi:10.1103/PhysRevB.79.035401 (2009). 4. Chu, Y. & Crozier, K. B. Experimental study of the interaction between localized and propagating surface plasmons. Opt. Lett. 34, 244–246, http://ol.osa.org/abstract.cfm?URI=ol-34-3-244 (2009). 5. Nielsen, M. G., Gramotnev, D. K., Pors, A., Albrektsen, O. & Bozhevolnyi, S. I. Continuous layer gap plasmon resonators. Opt. Express 19, 19310–19322, http://www.opticsexpress.org/abstract.cfm?URI=oe-19-20-19310 (2011). 6. Yang, J. et al. Ultrasmall metal-insulator-metal nanoresonators: impact of slow-wave effects on the quality factor. Optics Express 20, 16880–16891 (2012). 7. Lassiter, J. B. et al. Third-harmonic generation enhancement by film-coupled plasmonic stripe resonators. Acs Photonics 1, 1212–1217 (2014). 8. Moreau, A., Ciraci, C. & Smith, D. R. Impact of nonlocal response on metallodielectric multilayers and optical patch antennas. Physical Review B 87, 045401 (2013). 9. Dechaux, M. et al. Influence of spatial dispersion in metals on the optical response of deeply subwavelength slit arrays. Physical Review B 93, 045413 (2016). 10. Nielsen, M. G., Pors, A., Albrektsen, O. & Bozhevolnyi, S. I. Efficient absorption of visible radiation by gap plasmon resonators. Optics express 20, 13311–13319 (2012). 11. Moreau, A. et al. Controlled-reflectance surfaces with film-coupled colloidal nanoantennas. Nature 492, 86–89 (2012). 12. Akselrod, G. M. et al. Large-area metasurface perfect absorbers from visible to near-infrared. Advanced Materials 27, 8028–8034 (2015). 13. Tittl, A. et al. Palladium-based plasmonic perfect absorber in the visible wavelength range and its application to hydrogen sensing. Nano letters 11, 4366–4369 (2011). 14. Cattoni, A. et al. λ 3/1000 plasmonic nanocavities for biosensing fabricated by soft uv nanoimprint lithography. Nano letters 11, 3557–3563 (2011). 15. Powell, A. W. et al. Plasmonic gas sensing using nanocube patch antennas. Advanced Optical Materials (2016). 16. Pors, A., Albrektsen, O., Radko, I. P. & Bozhevolnyi, S. I. Gap plasmon-based metasurfaces for total control of reflected light. Scientific reports 3 (2013). 17. Belacel, C. et al. Controlling spontaneous emission with plasmonic optical patch antennas. Nano letters 13, 1516–1521 (2013). 18. Akselrod, G. M. et al. Probing the mechanisms of large purcell enhancement in plasmonic nanoantennas. Nature Photonics 8, 835–840 (2014). 19. Hoang, T. B. et al. Ultrafast spontaneous emission source using plasmonic nanoantennas. Nature Communications 6, 7788 (2015). 20. Faggiani, R., Yang, J. & Lalanne, P. Quenching, plasmonic, and radiative decays in nanogap emitting devices. ACS Photonics 2, 1739–1744 (2015). 21. Cirac, C., Lassiter, J. B., Moreau, A. & Smith, D. R. Quasi-analytic study of scattering from optical plasmonic patch antennas. Journal of Applied Physics 114, 163108 (2013). 22. Bowen, P. T. & Smith, D. R. Coupled-mode theory for film-coupled plasmonic nanocubes. Phys. Rev. B 90, 195402 (2014). 23. Akselrod, G. M. et al. Leveraging nanocavity harmonics for control of optical processes in 2d semiconductors. Nano Letters 15, 3578–3584 (2015). 24. Alaee, R. et al. Deep-subwavelength plasmonic nanoresonators exploiting extreme coupling. Nano Letters 13, 3482–3486, doi:10.1021/nl4007694, PMID: 23805879 (2013). 25. Watts, C. M., Liu, X. & Padilla, W. J. Metamaterial electromagnetic wave absorbers. Advanced Materials 24 (2012). 26. Wan, W. et al. Time-reversed lasing and interferometric control of absorption. Science 331, 889–892 (2011). 27. Chong, Y. D., Ge, L., Cao, H. & Stone, A. D. Coherent perfect absorbers: Time-reversed lasers. Phys. Rev. Lett. 105, 053901 (2010). 28. Le Perchec, J., Quemerais, P., Barbara, A. & Lopez-Rios, T. Why metallic surfaces with grooves a few nanometers deep and wide may strongly absorb visible light. Physical review letters 100, 066408 (2008). 29. Rakić, A. D., Djurišić, A. B., Elazar, J. M. & Majewski, M. L. Optical properties of metallic films for vertical-cavity optoelectronic devices. Applied optics 37, 5271–5283 (1998). 30. Defrance, J. et al. Moosh: A numerical swiss army knife for the optics of multilayers in octave/matlab. Journal of Open Research Software 4 (2016). 31. Chandran, A., Barnard, E. S., White, J. S. & Brongersma, M. L. Metal-dielectric-metal surface plasmon-polariton resonators. Phys. Rev. B 85, 085416, doi:10.1103/PhysRevB.85.085416 (2012). 32. Moreau, A., Lafarge, C., Laurent, N., Edee, K. & Granet, G. Enhanced transmission of slit arrays in an extremely thin metallic film. Journal of Optics A: Pure and Applied Optics 9, 165 (2007). 33. Granet, G. Reformulation of the lamellar grating problem through the concept of adaptive spatial resolution. JOSA A 16, 2510–2516 (1999). 34. Shalaev, V. M. et al. Negative index of refraction in optical metamaterials. Optics letters 30, 3356–3358 (2005). Acknowledgements The authors want to thank Kofi Edee, Gérard Granet and Armel Pitelet for fruitful discussions. This work has been supported by the French National Research Agency, “Physics of Gap-Plasmons” project number ANR- 13-JS10-0003. Author Contributions C.L. made the theoretical calculations and wrote the codes. A.M. supervised the work providing help for the code and the physical analysis. All the authors provided scientific insight and contributed to the writing of the manuscript. Additional Information Supplementary information accompanies this paper at doi:10.1038/s41598-017-03064-6 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. Scientific Repo R ts | 7: 2941 | DOI:10.1038/s41598-017-03064-6 6 www.nature.com/scientificreports/ 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 Cre- ative 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 per- mitted 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 Scientific Repo R ts | 7: 2941 | DOI:10.1038/s41598-017-03064-6 7 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Scientific Reports Springer Journals

Interferometric control of the absorption in optical patch antennas

Free
7 pages

Loading next page...
 
/lp/springer_journal/interferometric-control-of-the-absorption-in-optical-patch-antennas-JSPVRfWPbT
Publisher
Springer Journals
Copyright
Copyright © 2017 by The Author(s)
Subject
Science, Humanities and Social Sciences, multidisciplinary; Science, Humanities and Social Sciences, multidisciplinary; Science, multidisciplinary
eISSN
2045-2322
D.O.I.
10.1038/s41598-017-03064-6
Publisher site
See Article on Publisher Site

Abstract

www.nature.com/scientificreports OPEN Interferometric control of the absorption in optical patch antennas Received: 2 February 2017 Caroline Lemaître, Emmanuel Centeno & Antoine Moreau Accepted: 21 April 2017 Optical patch nano-antennas possess unique absorption, field enhancement and concentration Published: xx xx xxxx capabilities – but their crosssection, as well as their response outside of normal incidence are not well understood. Here we explain the large cross-section by considering that each patch nanoantenna is a cavity excited from both sides. Such a simple physical picture allows to fully understand the influence of the angle of incidence – that odd resonances have a very high absorption cross-section which decreases when the incidence angle increases, while even resonances cannot be excited in normal incidence. A direct application would be to use these structures as an optical nanometric set-square. Gap-plasmon resonators are a new kind of optical metallic resonator , that has emerged as particularly prom- 2–6 7–9 ising , from a fundamental point of view as well as for all the applications that could be based on these 2, 10 structures given their extraordinary efficiency at concentrating light . It has actually been suggested that they 11, 12 13–15 could be used as cheap but very efficient selective absorbers , as sensors , as well as for designing reflection metasurfaces providing an almost complete control of light . Finally, the light confinement associated to their resonances is so high that an unprecedentedly high Purcell effect can be obtained using this kind of optical patch 17–20 antennas . These resonators are literally cavities for a peculiar plasmonic mode called a gap-plasmon which presents a very large effective index. This property alone explains how deeply sub-wavelength the resonators can be made, which in turn partly explains their efficiency: the absorption cross section of a single patch can be actually as high as 30 times its geometrical cross-section . The modes that can be excited in these cavities have already been 2, 21, 22 23 studied , and their importance largely underlined . However, the influence of the incidence angle, despite its 14, 24 25 practical importance for sensing for instance, is not well understood - and their extraordinary absorption cross-section hasn’t yet received a full explanation. Here we show that the influence of the incidence angle can be accurately taken into account if a single gap-plasmon resonator is modeled as a cavity excited from both ends, making the structure a clear example of 26 27 interferometrically controlled absorption , or of a coherent (perfect) absorber . We give analytic formula for the losses induced by the resonances inside the cavity and link it to the incidence angle. Such a description allows to understand why the absorption decreases when the incidence angle increases and explains the extraordinary 9, 28 absorption cross-section of the patch antenna compared to other kinds of gap-plasmon resonators . An analytic model for the absorption of a gap-plasmon resonator A selectively absorbing metasurface is typically constituted of periodically or randomly distributed patches separated from a metallic film by a thin dielectric material, as shown on Fig.  1. Here we will consider silver for all the metallic parts . The efficiency of the device benefits from the capacity of each patch to funnel light into the gap, so that when the resonators are close enough to each other, they are able to absorb all the incoming light and constitute a perfect absorber . Most of the structures that are fabricated consist in rectangular patches coupled to a film instead of the nanorods shown Fig.  1, but is largely admitted now that the physical behavior of both struc- tures is the same and that the die ff rence, for the fundamental modes, essentially reduces to a shift of the spectral 11 22 position of the resonances . For higher order modes, although they may be more complex , the fact that some are excited only for non-normal incidence is perfectly in agreement with what we expose here for 2D structures. We must underline that when the structure is periodic, surface plasmons are oen ex ft cited , and can be respon- sible for a very large absorption in very peculiar conditions. Such a mechanism is not present in disordered structures, where the surface plasmons are excited and absorb a non-negligible part of the energy but play a Université Clermont Auvergne, CNRS, Institut Pascal, F-63000, Clermont-Ferrand, France. Correspondence and requests for materials should be addressed to A.M. (email: antoine.moreau@uca.fr) Scientific Repo R ts | 7: 2941 | DOI:10.1038/s41598-017-03064-6 1 www.nature.com/scientificreports/ Figure 1. A 3D rendering of a selectively absorbing plasmonic metasurface based on gap-plasmon resonators: regularly spaced metallic nanorods coupled to a metallic substrate and separated from it by a nanometric dielectric layer. lesser role than when periodicity enhances their excitation. Here we will concentrate on the resonance of a single resonator, which is angularly and spectrally broader, is more efficient, does not depend on the periodicity and explains the response of disordered structures . In the thin layer sandwiched between the patch and the metallic film, light propagates under the form of a gap-plasmon along the x axis, the only propagating mode of a metal-dielectric-metal waveguide. This mode exists for p polarization only, with a magnetic field along the y axis. Solving Maxwell’s equations a looking for a guided mode with a pulsation ω and a wavevector k in a metal-dielectric-metal configuration actually yields a magnetic field that can be written ik () xt −ω Hx (, zP )( = ze ) y (1) with a vertical profile of the magnetic field  −κ z ez ∀> hd    Pz () = Cz cosh κ ∀∈ z − ,      κ z d ez ∀< − 2 (2)   −κ where to ensure the continuity of the magnetic field at all interfaces. Figure  2(a) shows the profile of C = κ d cosh the magnetic field associated to the gap-plasmon propagating in a nano-gap. e Th dispersion relation of the gap-plasmon is obtained using the continuity of the electric field at any inter - face and can be written   κ d κ 2 1 tanh κ += 0  2   ε 2 ε (3) dm 2 2 2 2 where κ =− αε k and κ =− αε k , α being the propagation constant in the x direction. Once the 1 m 0 2 d 0 pulsation ω is fixed, the dispersion relation allows to retrieve the complex propagation constant α . The main characteristics of the gap-plasmon is that its effective index, defined as n = R() α /k , diverges when the size of eff 0 the gap tends to zero, as shown Fig. 2(b). This plasmonic ee ff ct occurs when light propagates more in the metal than in the dielectric. This is what allows gap-plasmon resonators to be so small compared to the wavelength in vacuum. When the patch ends, the gap-plasmon can in fact no longer propagate. Instead, it is able to excite the surface plasmon and a waves propagating in free space , but as both present very low propagation constants, the outside space represents a very high impedance mismatch for the gap-plasmon. This is of course all the more so that the gap is small and the effective index of the gap-plasmon large. As a consequence, the reflection coefficient r of the 5, 11, 24 gap-plasmon when it encounters the edge of the patch can be very high . Computing this reflection coefficient is not an easy task. Its phase can be deduced from the location of the 5 24, 31 resonances themselves , and more direct determinations can be difficult too . Finally, it can be noticed that this can be done directly using a Fourier Modal Method , which relies on the computation of the modes that prop- agate in different vertically invariant layers to solve Maxwell’s equations and of the transmission and reflection coefficients. It directly gives access to r and α. It is even possible, relying on transformation optics techniques to introduce Perfect Matching Layers in order to get rid of the periodicity. Figure 2(c) shows the modulus of the computed reflection coefficient as a function of the waveguide width, computed using such techniques (see Methods). The phase of the reflection coefficient, which is critical for the resonance condition and thus to explain why the size of a resonator can be smaller than half of the gap-plasmon wavelength , is shown on Fig. 2(c). This phase increases when the the gap closes. Scientific Repo R ts | 7: 2941 | DOI:10.1038/s41598-017-03064-6 2 www.nature.com/scientificreports/ Figure 2. (a) Profile of the magnetic field associated with the gap-plasmon propagating along the x axis. (b) Effective index of the gap-plasmon as a function of the gap width d. (c) Modulus of the reflection R() α /k coefficient at the edge of the cavity as a function of the gap width d. Figure 3. Cavity model. (a) Sketch representing the cavity model including the reflection coefficients of the gap-plasmon on the edges of the patch. (b) Map of the magnetic field (modulus) for a single 75 nm wide square patch, separated from the metallic film by a 10 nm thick dielectric layer with an optical index of 1.54 illuminated from above in normal incidence at resonance. The field profile in the region which is enclosed in the white rectangle is accurately predicted by the model (center). The inset shows the map of the magnetic field (modulus) as predicted by the cavity model when both excitation coefficients A and A are taken equal, at resonance. l r Under the patch, gap-plasmons thus propagate from right to left and from left to right and are reflected on both sides of the patch, so that they form a stationary wave. This is the actual description of a cavity of size h , as shown on Fig. 3(a). This allows us to introduce a model for the gap-plasmon resonator, based on our knowledge of r and α, that are computed as explained above. Since a gap-plasmon propagates along x towards the left and another one propagates towards the right, the amplitude of the mode, as a function of x under the patch can simply be written ix αα −ix () xA =+ eB e . (4) This means that the magnetic field associated to the stationary wave can simply be written H (, xz )( =  xP )(z ). We choose here the edges of the cavity to be in ± . We call A and A the amplitude with which the l r gap-plasmons are excited respectively on the left and on the right of the patch. We know actually very little about these quantities, except for the fact that in normal incidence, they are obviously equal: the cavity is being excited on both sides exactly the same way. In non-normal incidence, the phase difference between the two can be esti- mated, as the following will show, but there must be a difference in the amplitude due to the shading by the cube that is difficult to take into account. These excitations should be considered as relatively large because when a plane wave is reflected by a metallic film, the magnetic field presents a maximum on the surface. Whatever the thickness of the dielectric layer, there is always a maximum of the magnetic field at the entrance of the cavity. This explains why the cavity resonances are much more easily excited in this configuration than for dimer , for instance. The resonances are, conversely, not likely to be excited through the patches – because they are twice as large as the skin depth. This simple cavity model, detailed in the Supplementary Information, thus gives the amplitude of the modes in the cavity as hh hh  −− ii ααxiαα ix ii ααxi −− αα ix  22 22 ee + re e re ee + e   () xA = + A . l r   22ih α 22ih α 11 − re − re   (5)   Scientific Repo R ts | 7: 2941 | DOI:10.1038/s41598-017-03064-6 3 www.nature.com/scientificreports/ 2 2iαh As shown by this expression, a resonance occurs each time the quantity 1 − r e presents a minimum and can −ih α −ih α be considered close enough to zero. The resonance condition, in that case, reduces to re  ± . When re  then the resonance is symmetric with respect to the x = 0 axis, which means it has an odd number of anti-nodes d d −ih α between − and . On the contrary, when re  − the magnetic field is anti-symmetric with respect to x = 0, 2 2 and the resonance presents an even number of anti-nodes. Keeping that in mind for later, it is then possible to get analytic expressions of all the fields, magnetic as well as electric, that are linked to the cavity and the stationary wave that exists in it - even out of resonance. When the amplitude () x is multiplied by the vertical profile P(Z) of the gap-plasmon, then the magnetic field is almost perfectly reproduced, as shown Fig. 3 where a comparison to a full simulation H (, xz )( =  xP )(z ) using a Fourier Modal Method of a single patch is shown. Using Maxwell’s equations, the electric field can be analytically calculated (See Supplementary Information). Finally, the volumetric losses p(r) inside the metal are classically written p() rE =. ωε I() ε E . 0 r (6) As shown in the Supplementary Information, the volumetric losses can be integrated over the metallic region, to compute analytically the absorption under one patch antenna, as a function of A and A l r h/2 +∞ Pp =. () r dd xz ∫∫ −− h/2 ∞ Discussion Two different things can be done starting from here. Either use the complex expression (given in the Supplementary Information) as it is, to compute the losses in every possible situation (out of resonance, for any angle and wavelength), or to compute a simplified analytic solution that brings a better understanding but only −ih α close to resonance. Let us first assume that we are close to a resonance, with re  ± . Then, as shown in the Supplementary Information, the absorption by the resonator can be written 22 ⁎ PA =+ β(2 AA ± R() A ) (7) lr lr In normal incidence, A = A since both ends of the cavity are excited in the same way. This explains why res- l r iαh onances with an even number of anti-nodes, for which r = −e cannot be excited in normal incidence, since we have strictly P = 0. On the opposite, for a resonance with an odd number of anti-nodes P = 4β|A | and the absorption is thus four times higher than if the cavity were excited from one side only. This shows the mechanism that is at the heart of the extraordinary large absorption cross-section that patch antenna have, when compared to their actual size. The interferometric reinforcement of the absorption, associated with the high effective index of the gap-plasmon which drives the reduction in size of the resonators, completely explains this effect. Now it is possible to take into account the incidence angle, as long as it is not too large by simply assuming that the modulus of the excitations A and A stay close to their maximum, that is reached in normal incidence, r l but that the main factor driving their variation will be a slight phase shift arising from the different optical paths followed by light to reach the entrances of the gap. With a incidence angle θ, this phase shift is Δφ = k h sin θ. We finally have an absorption at resonance that can be written either   π h   PA = 2c β os sin(θ)  λ   (8) for odd resonances (i.e. with an odd number of anti-nodes) and π h 22  PA = 2s β in sin(θ) l   λ   0  (9) for even resonances (i.e. with an even number of anti-nodes). Changing the angle of incidence thus lowers the absorption of odd resonances, while allowing the even resonances to be excited and to appear in the spectrum. Our findings can be summarized by a simple scheme (see Fig.  4). Odd (resp. even) resonances are actually symmetric (resp. anti-symmetric). Since each side of the patch is excited by the incoming field, the field under the patch can be seen as a superposition of two modes: one excited from the left and one excited from the right. In normal incidence, both excitations are perfectly in phase. When the resonance should be odd (and symmetric), then the two modes add up so that the amplitude of the field is exactly twice what it would be if the resonance had been excited from one side only (and the absorption is multiplied by four). When the resonances should be even (and anti-symmetric) the two modes interfere destructively, and the resonance is completely canceled. For oblique incidence, the two excitations are not in phase anymore, disturbing the interference between the modes. Odd resonances see their efficiency decrease, while even resonances can be excited. Now we can use the model described above, without any assumption on r, out of the resonance condition to try to predict, as the incidence angle increases, how the absorption linked to the resonance will evolve. Using the general expression for the absorption obtained thanks to our model, we plot the absorption as a function of the patch width (see Fig. 5) and compare it to the reflectance of the structure for different incidence angles . The accu- racy of the model can be clearly seen. It is almost perfect in normal incidence and for large angles the variation of Scientific Repo R ts | 7: 2941 | DOI:10.1038/s41598-017-03064-6 4 www.nature.com/scientificreports/ Figure 4. e in Th terferometric control of the absorption in patch antennas summarized in normal incidence: adding two anti-symmetric modes excited from left and right leads to a complete cancellation, while adding even modes leads to an enhancement. For oblique incidence, this control is lowered, so that the even modes don’t completely cancel each other, and the odd modes are less enhanced. Figure 5. Absorption spectrum for a silver nanopatch of variable width (in nm) from the model (left) and from the reference code (right). The layer is 5 nm thick, the wavelength 950 nm. the interferometric effect is not perfectly reproduced - probably because our assumptions regarding the optical path and the shading are not fully accurate. However, this definitely validates our physical interpretation of the nanopatch optical response. We underline that the destructive interference that literally kills odd resonances in normal incidence can be utilized to assess whether the angle of incidence is actually zero or not. As soon as the incidence angle differs from zero, the interferometric control will stop and the absorption linked to the resonance will sharply increase, as shown on the spectrum. This would be the optical equivalent of an optical set square at a nanometric level. Conclusion We have shown that a optical patch antenna can be considered as a Fabry-Perot cavity excited from both sides and that this explains why the absorption by even resonances is reinforced by an interferometric control effect while odd resonances are completely tamed by the same effect. Previous authors have already evoked the symmetry 24, 31 mismatch between the incident wave and the resonances , but the interferometric control pictures offers a way to quantify this effect. Using an analytic model whose parameters are computed numerically, we have been able to show how this effect varies when the incidence angle is changed, allowing odd resonances to be excited and slightly lowering the absorption efficiency of even resonances. This gives a better understanding of why the absorption of metasurfaces made with such resonators decreases when the angle of incidence increases - and shows that modifying the angle of incidence is an actual way to control the absorption cross section of each individual resonator and thus light concentration of the fluorescence enhancement for instance . We think that having in mind this fundamental mechanism of the absorption by optical patch antennas will help find new applications or new techniques. It may help to forecast the radiation pattern of film-coupled nanocubes , and shows that the angle of incidence is a way to control the absorption cross-section of a single patch. Finally, the fact that the excitation of odd resonances is strictly impossible at normal incidence could well be utilized to fabricate an optical set-square, giving an immediate application to the interferometric control of the absorption by optical patch antennas. Methods We used a Fourier Modal Method with coordinate stretching and Perfect Matching Layers to first solve the modal problem in a metal-insulator-metal structure. The period is of 930 nm with a 200 nm thick PML and a stretching parameter η = 0.99. On one side, we have 120 nm of metal and 600 nm on the other side. We take the mode with the smallest imaginary part for the propagation constant. This is the gap plasmon. We then compute the scattering matrix when the above structure is put next to a structure with the same parameters, except that the 600 nm metallic part is replaced by air. The rest of the parameters are kept the same. The mode with the lowest imaginary part for its propagation constant is then the surface plasmon. Using the boundary conditions allows to access the scattering matrix of the problem and then retrieve the coefficient of the scattering matrix corresponding to an incoming and outgoing gap-plasmon. The coefficient constitutes the reflection coefficient of the gap-plasmon on each edge of the structure. Scientific Repo R ts | 7: 2941 | DOI:10.1038/s41598-017-03064-6 5 www.nature.com/scientificreports/ References 1. Bozhevolnyi, S. & Sondergaard, T. General properties of slow-plasmon resonant nanostructures: nano-antennas and resonators. Optics Express 15, 10869 (2007). 2. Søndergaard, T. & Bozhevolnyi, S. Slow-plasmon resonant nanostructures: Scattering and field enhancements. Phys . Rev. B 75, 073402, doi:10.1103/PhysRevB.75.073402 (2007). 3. Jung, J., Søndergaard, T. & Bozhevolnyi, S. I. Gap plasmon-polariton nanoresonators: Scattering enhancement and launching of surface plasmon polaritons. Phys. Rev. B 79, 035401, doi:10.1103/PhysRevB.79.035401 (2009). 4. Chu, Y. & Crozier, K. B. Experimental study of the interaction between localized and propagating surface plasmons. Opt. Lett. 34, 244–246, http://ol.osa.org/abstract.cfm?URI=ol-34-3-244 (2009). 5. Nielsen, M. G., Gramotnev, D. K., Pors, A., Albrektsen, O. & Bozhevolnyi, S. I. Continuous layer gap plasmon resonators. Opt. Express 19, 19310–19322, http://www.opticsexpress.org/abstract.cfm?URI=oe-19-20-19310 (2011). 6. Yang, J. et al. Ultrasmall metal-insulator-metal nanoresonators: impact of slow-wave effects on the quality factor. Optics Express 20, 16880–16891 (2012). 7. Lassiter, J. B. et al. Third-harmonic generation enhancement by film-coupled plasmonic stripe resonators. Acs Photonics 1, 1212–1217 (2014). 8. Moreau, A., Ciraci, C. & Smith, D. R. Impact of nonlocal response on metallodielectric multilayers and optical patch antennas. Physical Review B 87, 045401 (2013). 9. Dechaux, M. et al. Influence of spatial dispersion in metals on the optical response of deeply subwavelength slit arrays. Physical Review B 93, 045413 (2016). 10. Nielsen, M. G., Pors, A., Albrektsen, O. & Bozhevolnyi, S. I. Efficient absorption of visible radiation by gap plasmon resonators. Optics express 20, 13311–13319 (2012). 11. Moreau, A. et al. Controlled-reflectance surfaces with film-coupled colloidal nanoantennas. Nature 492, 86–89 (2012). 12. Akselrod, G. M. et al. Large-area metasurface perfect absorbers from visible to near-infrared. Advanced Materials 27, 8028–8034 (2015). 13. Tittl, A. et al. Palladium-based plasmonic perfect absorber in the visible wavelength range and its application to hydrogen sensing. Nano letters 11, 4366–4369 (2011). 14. Cattoni, A. et al. λ 3/1000 plasmonic nanocavities for biosensing fabricated by soft uv nanoimprint lithography. Nano letters 11, 3557–3563 (2011). 15. Powell, A. W. et al. Plasmonic gas sensing using nanocube patch antennas. Advanced Optical Materials (2016). 16. Pors, A., Albrektsen, O., Radko, I. P. & Bozhevolnyi, S. I. Gap plasmon-based metasurfaces for total control of reflected light. Scientific reports 3 (2013). 17. Belacel, C. et al. Controlling spontaneous emission with plasmonic optical patch antennas. Nano letters 13, 1516–1521 (2013). 18. Akselrod, G. M. et al. Probing the mechanisms of large purcell enhancement in plasmonic nanoantennas. Nature Photonics 8, 835–840 (2014). 19. Hoang, T. B. et al. Ultrafast spontaneous emission source using plasmonic nanoantennas. Nature Communications 6, 7788 (2015). 20. Faggiani, R., Yang, J. & Lalanne, P. Quenching, plasmonic, and radiative decays in nanogap emitting devices. ACS Photonics 2, 1739–1744 (2015). 21. Cirac, C., Lassiter, J. B., Moreau, A. & Smith, D. R. Quasi-analytic study of scattering from optical plasmonic patch antennas. Journal of Applied Physics 114, 163108 (2013). 22. Bowen, P. T. & Smith, D. R. Coupled-mode theory for film-coupled plasmonic nanocubes. Phys. Rev. B 90, 195402 (2014). 23. Akselrod, G. M. et al. Leveraging nanocavity harmonics for control of optical processes in 2d semiconductors. Nano Letters 15, 3578–3584 (2015). 24. Alaee, R. et al. Deep-subwavelength plasmonic nanoresonators exploiting extreme coupling. Nano Letters 13, 3482–3486, doi:10.1021/nl4007694, PMID: 23805879 (2013). 25. Watts, C. M., Liu, X. & Padilla, W. J. Metamaterial electromagnetic wave absorbers. Advanced Materials 24 (2012). 26. Wan, W. et al. Time-reversed lasing and interferometric control of absorption. Science 331, 889–892 (2011). 27. Chong, Y. D., Ge, L., Cao, H. & Stone, A. D. Coherent perfect absorbers: Time-reversed lasers. Phys. Rev. Lett. 105, 053901 (2010). 28. Le Perchec, J., Quemerais, P., Barbara, A. & Lopez-Rios, T. Why metallic surfaces with grooves a few nanometers deep and wide may strongly absorb visible light. Physical review letters 100, 066408 (2008). 29. Rakić, A. D., Djurišić, A. B., Elazar, J. M. & Majewski, M. L. Optical properties of metallic films for vertical-cavity optoelectronic devices. Applied optics 37, 5271–5283 (1998). 30. Defrance, J. et al. Moosh: A numerical swiss army knife for the optics of multilayers in octave/matlab. Journal of Open Research Software 4 (2016). 31. Chandran, A., Barnard, E. S., White, J. S. & Brongersma, M. L. Metal-dielectric-metal surface plasmon-polariton resonators. Phys. Rev. B 85, 085416, doi:10.1103/PhysRevB.85.085416 (2012). 32. Moreau, A., Lafarge, C., Laurent, N., Edee, K. & Granet, G. Enhanced transmission of slit arrays in an extremely thin metallic film. Journal of Optics A: Pure and Applied Optics 9, 165 (2007). 33. Granet, G. Reformulation of the lamellar grating problem through the concept of adaptive spatial resolution. JOSA A 16, 2510–2516 (1999). 34. Shalaev, V. M. et al. Negative index of refraction in optical metamaterials. Optics letters 30, 3356–3358 (2005). Acknowledgements The authors want to thank Kofi Edee, Gérard Granet and Armel Pitelet for fruitful discussions. This work has been supported by the French National Research Agency, “Physics of Gap-Plasmons” project number ANR- 13-JS10-0003. Author Contributions C.L. made the theoretical calculations and wrote the codes. A.M. supervised the work providing help for the code and the physical analysis. All the authors provided scientific insight and contributed to the writing of the manuscript. Additional Information Supplementary information accompanies this paper at doi:10.1038/s41598-017-03064-6 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. Scientific Repo R ts | 7: 2941 | DOI:10.1038/s41598-017-03064-6 6 www.nature.com/scientificreports/ 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 Cre- ative 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 per- mitted 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 Scientific Repo R ts | 7: 2941 | DOI:10.1038/s41598-017-03064-6 7

Journal

Scientific ReportsSpringer Journals

Published: Jun 7, 2017

References

You’re reading a free preview. Subscribe to read the entire article.


DeepDyve is your
personal research library

It’s your single place to instantly
discover and read the research
that matters to you.

Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.

All for just $49/month

Explore the DeepDyve Library

Search

Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly

Organize

Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.

Access

Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.

Your journals are on DeepDyve

Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.

All the latest content is available, no embargo periods.

See the journals in your area

DeepDyve

Freelancer

DeepDyve

Pro

Price

FREE

$49/month
$360/year

Save searches from
Google Scholar,
PubMed

Create lists to
organize your research

Export lists, citations

Read DeepDyve articles

Abstract access only

Unlimited access to over
18 million full-text articles

Print

20 pages / month

PDF Discount

20% off