Revealing Low-Radiative Modes of Nanoresonators with Internal Raman Scattering

Revealing Low-Radiative Modes of Nanoresonators with Internal Raman Scattering

K. V. Baryshnikova    K. Frizyuk    G. Zograf    S. Makarov    M. A.  Baranov    D. Zuev    V. A. Milichko    I. Mukhin    M. Petrov    and A. B. Evlyukhin e-mail: Department of Nanophotonics and Metamaterials, ITMO University, St. Petersburg, Russia
Laboratory of Renewable Energy Sources, St. Petersburg Academic University, St. Petersburg, Russia
Institute of Quantum Optics, Leibniz Universität Hannover, 30167 Hannover, Germany

Revealing hidden non-radiative (dark) of resonant nanostructures using optical methods such as dark-field spectroscopy often becomes a sophisticated problem due to a weak coupling of these modes with a far-field radiation, whereas methods of dark-modes spectroscopy, e.g. cathodoluminescence or elastic energy losses, are not always convenient in use. Here, we suggest an approach for experimental determining the mode structure of a nanoresonator basing on utilizing intrinsic incoherent Raman scattering. We theoretically predict the efficiency of this approach and realize it experimentally for silicon nanoparticle resonators possessing strong Raman line at 520 cm. With this method we studied a silicon nanoparticle placed on a gold substrate and reveal the spectral position of a low-radiative magnetic quadrupole mode which is hardly observable with common dark-field optical spectroscopy.


Revealing Low-Radiative Modes of Nanoresonators… \sodtitleRevealing Low-Radiative Modes of Nanoresonators with Inherent Spontaneous Raman Scattering \rauthorK. V. Baryshnikova, K. Frizyuk, G. Zograf … \sodauthorBaryshnikova, Frizyuk, Zograf, Makarov, Baranov, Zuev, Milichko, Mukhin, Petrov, Evlyukhin \dates18 April 2019 \PACS78.30.−j, 61.46.+w

1 Introduction

The resonantly enhanced optical interaction of light with metallic, dielectric, or semiconductor nanoparticles, which are the basic components of various nanooptical devices, is at the forefront of modern nanophotonics [1, 2, 3, 4]. The examples of vital applications include on-chip integrated photonic systems [5, 6, 7], lab-on-chip technologies [8], medical biosensing [9] and others. However, understanding the modal content of resonant nanostructures is non-trivial problem due to their substantially subwavelength sizes and typically complicated geometries. This problem has very limited solution within the common optical spectroscopy methods as some of the modes can be low-radiative because of their weak coupling to freely propagating waves, and in this sense are often known as dark-modes [10, 11, 12]. Thus, studying the modal spectrum of optical nanostructures with account on dark modes immediately becomes a complicated problem.

Figure 1: The elastic and inelastic scattering by bright and dark resonant modes in silicon nanoparticle on the top of gold substrate. (a) The excitation of the bright nanoparticle mode results in both elastic scattering and inelastic (Raman) scattering. (b) The excitation of the dark mode does not contribute to elastic scattering, but can be detected via Raman scattering signal.

The solution of this problem usually relies on sophisticated optical methods such as total internal reflectance microscopy (TIRM) [13] and its combination with atomic-force microscopy [14, 15]. The methods based on resonant mode excitation through accelerated electron scattering showed their efficiency in observing the dark-modes of nanostructures. In particular, the elastic electron loss spectroscopy (EELS) is very powerful method for detecting the modal structure of plasmonic resonators by measuring both radiative and non-radiative losses of electrons interacting with the electron excitations in the nanostructures [16, 17]. The cathodoluminescence (CL) method similarly to EELS uses accelerated electrons, but excites the dark mode via CL radiation generated inside a nanoresonator. In CL technique, the electronically excited luminescence signal inside the nanostructures is enhanced by the resonant modes, allowing spatial and spectral mapping of dark modes in both plasmonic [18, 19] and non-metal [20, 21, 22] structures.

In this Letter, we suggest detecting an intrinsic Raman signal, which corresponds to inelastic photon scattering via optical phonons excitation in a crystalline lattice of a nanostructured solid to determine its electromagnetic modal structure. The Raman signal generated by lattice vibrations is distributed over the volume of the nanoresonator and effectively couples to all resonant modes including dark ones. By exciting the optical phonons and their consequent spontaneous decay, the incoherent Stokes signal is generated inside the nanoresonator under study. The intensity of the Raman signal is governed by the optical density of states inside the nanoresonator at the frequency of the Raman emission, and additionally enhanced at the excitation frequency by resonant excitation of the eigenmodes. We show that the approach allows to reveal magnetic quadrupole mode in a rather complicated system of a silicon nanoparticle arranged on a gold substrate as well as an anapole state, demonstrating the applicability of our approach even when dark-field optical spectroscopy methods are not efficient. Silicon as one of the main materials used in all-dielectric nanophotonics [1, 23, 3] has a well pronounced peak in Raman spectrum at  cm, which is defined by the light coupling with optical phonons at the band central point [24, 25, 26]. The system under study consists of a subwavelength silicon nanoparticle (SiNP) placed over a gold substrate as show in Figure 1. We note that the effect of a metallic substrate on Mie resonances of SiNP has recently attracted much attention due to the effective optical bi-anisotropy of such system [27, 28, 29, 30], surface plasmon polariton launching[31, 32] and biosensing[33]. It should be stressed that approach can be applied to detect any resonant mode excited in any structures that have inherent Raman scattering at visible range.

2 Results and discussion

Figure 2: (a) Experimentally measured and numerically calculated dark field spectra of nanoparticle with radius 125 nm placed on glass substrate. The angle of incidence is 65. The collection aperture and the SEM image of silicon nanoparticle placed on the glass substrate are shown in figure by insets. (b) Numerically calculated multipoles contributions to scattering of nanoparticle with radius 125 nm placed on the glass substrate. (c) The same as for (a) but for case of gold substrate. (d) The same as for (b) but for case of gold substrate.

We started our consideration with the spectral properties of elastic scattering on the spherical nanoparticles. The nanoparticles of different radii  nm were fabricated by using laser ablation technique [34, 35] and then deposited on glass and gold substrates. The distances between neighboring particles was much larger than the excitation wavelength in order to avoid their interaction. Than their geometrical parameters were investigted with scanning electron microscopy (SEM). Thus, SEM picture of a nanoparticle with radius about 125 nm is shown at inset in Figure 2(a). Next, the dark-field spectral measurements were performed in the geometry of oblique incidence, when the first objective (Mitutoyo M Plan Apo NIR x10 0.26 NA) is used to illuminate the sample by polarized light at the angle of 67 to the surface normal with broadband light source (HL-2000 halogen lamp), and the second one (Mitutoyo M Plan Apo NIR x50 0.42 NA) is used to collect the scattered light. Inelastic Raman scattering on the SiNPs was studied on the same setup but with the He-Ne laser as a light source with a wavelength of 632.8 nm and intensity which was cut by a notch filter in collected light. The scattered signal subsequently was analyzed by a confocal setup with a spectrometer (HORIBA LabRam HR) and a cooled charge-couled device (CCD) Camera (Andor DU 420A-OE 325) equipped with a 150 g/mm (dark-field experiments) and a 600 g/mm (Raman scattering) diffraction gratings. The position of the nanoparticles was controlled optically via camera (Canon 400D). The measured (solid curves) typical dark-field spectra along with numerical modeling (dashed curves) are shown in Figure 2(a) for nanoparticle of radii 125 nm placed on glass substrate and in Figure 2(b) for nanoparticle placed on gold substrate. Multipole decompositions for both types of substrate are shown also in Figure 2(c,d).

We note that the optical response of nanoparticle on gold substrate significantly differs from the case of nanoparticle in vacuum or on glass substrate due to a strong coupling with the metal substrate [27]. This is clearly seen after calculating multipoles excited in nanoparticles [30]. One can see that the electric dipole mode gives a predominant contribution to the dark field signal, and its resonances are coincide with the peaks of dark field spectra. At the same time, in this geometry the magnetic mode gives small contribution into the scattering. Magnetic quadrupole (MQ) moment, which was calculated numerically using Comsol Multyphysics software, is of very high-quality in the vicinity of 700 nm, however its presence on dark field spectra is almost neglected. The minimum in the scattering in the vicinity of MQ resonance can be explained with the multipole decomposition spectra; as one can see, other multipoles are also suppressed in this region, especially electric dipole, which corresponds to the case of so-called anapole state[29, 36, 37]. Weak scattering of MQ stems from high quality factor and absorbance in silicon, which is reflected in the spectrum of the electromagnetic energy stored inside the nanoparticle but not in the scattering spectra. Recently, the strong field enhancement at the anapole frequency was utilized for observing a number of effect such as third harmonic generation [36, 38] or nanolasing [39, 40]. The enhancement of electric energy inside the nanoparticle at the MQ mode with anapole state has a key importance for Raman signal enhancement. Indeed, the Raman scattering can be described semi-classically by introducing the Raman tensor  [25]. The subscript index corresponds to three independent polarizations and of phonon modes corresondingly, which are mutually incoherent. Then polarization vector at the point , which generates Raman signal emission will be as follows [25]:


where is exciting field inside the nanoparticle, and the Raman tensor has the following non-zero components:


where is the Raman tensor amplitude.

The Raman emission is a spontaneous process, and, thus, it is defined not only by the amplitude of Raman polarization, but also by local density of electromagnetic states, which is known as Purcell effect [41, 42]. Then the Raman intensity [43]


where is the Purcell factor calculated at the frequency of Stocks line , which is lower than the exciting frequency for the frequency of optical phonon , the integral is taken over the whole volume of nanoparticle. The considered SiNP over gold substrate manifest itself as an open resonator, and the Purcell effect in such systems has been intensively discussed lately [44, 42, 45]. We have used numerical simulations in Comsol Multyphysics package to model the Raman emission spectrum. The results for SiNP of 115 nm radius are plotted in Figure 3 (a) for different excitation wavelengths along with the electromagnetic energy spectrum. Nanoparticle radius was chosen for realization of MQ resonance at 633 nm, which corresponds to the excitation wavelength in the experiment. It should to be stressed that wavelength of MQ resonance stays almost the same for both cases of glass and gold substrate. One can see that the Raman emission intensity generally follows the electromagnetic energy and is enhanced by the resonances of the SiNP. In particular, it is strongly enhanced at the magnetic quadrupole resonance. However, the peak of the Raman emission resonance at the MQ resonance is split as the Purcell factor enhancement is shifted for optical phonon frequency relatively to the peak of the exciting field enhancement at . For this spectral range the optical phonon frequency corresponds to  nm spectral shift for HeNe 633 nm laser excitation, (520cm), which is comparable to MQ width, that provides the splitting of the Raman intensity peak. This effect is smeared for electric dipole mode which is spectrally much wider than the Raman shift. Herewith total electromagnetic power inside the nanoparticle placed on glass substrate is much smaller than this one for case of gold substrate. The same behaviour is observed for simulated Raman signal. It origins from the high enhancement of fields inside the resonator in the vicinity of metal substrate and the following effect of bianisotropy.

Figure 3: (a) The modeled spectra of electromagnetic energy inside the nanoparticle with radius  nm placed on gold (dashed red line) or glass (solid blue line) substrate. The modelled dependencies of Raman scattering are shown as function of the excitation wavelength with scattered marks. Inset: typical Raman scattering spectra for a single SiNP placed on glass (blue line) or gold (red line) substrate. (b) The maximal Raman intensity measured from SiNP of different diameters placed on gold substrate at 633 nm excitation wavelength (red dots). The red solid line corresponds to numerical simulations.

To verify these assumptions experimentally, we measured the cold Raman emission intensity from SiNP of different sized and arranged on a top of gold substrate. The typical Raman spectrum from silicon nanoparticle placed on gold and glass substrate with sharp maximum around 520 cm is shown in Figure 3(b) by inset which fully prove our assumptions. The excitation wavelength used in the equipment was fixed at 633 nm not allow making spectral sweep of the maximal Raman intensity as shown in Figure 3(a). However, in order to study the effect of Raman enhancement at different resonances we considered a set of nanoparticles of different radii  nm. Due to the properties of Mie resonances their resonant wavelengths shift linearly to the long wavelength region with increasing nanoparticle size and the multipole resonances coincide with the excitation wavelength one after the another. Finally, for different size of SiNP we obtained dependence shown in Figure 3(b), which is obtained after normalizing the Raman intensity over nanoparticle’s volume. The measured experimental data (red circles) show a good correspondence to the results of numerical modeling (red solid line). The important result here is that the specific contribution of MQ mode in Raman intensity is much higher than for dipole resonances, contrary to the case of elastic scattering, when the MQR is almost undetectable in total scattering spectra.

3 Conclusions

In conclusion, we have suggested an efficient approach for experimental determining the spectral content of nanoresonator basing on utilizing an incoherent Raman scattering, inherent to the resonator material. The approach was confirmed theoretically and experimentally for silicon nanoparticle resonators possessing a strong intrinsic Raman response. We have shown that the suggested method allows us to reveal the magnetic quadrupole resonance combining with the anapole state in the system of a silicon nanoparticle arranged on a gold substrate, demonstrating the applicability of our approach for broad range of non-plasmonic resonant nanostructures.

Numerical simulations and ivestigation of anapole states was financially supported by Russian Science Foundation (project 17-72-10230). The experimental part was supported by Russian Science Foundation (project 18-79-00338).


  • [1] Kuznetsov, A. I.; Miroshnichenko, A. E.; Brongersma, M. L.; Kivshar, Y. S.; Lukyanchuk, B. Science 2016, 354, aag2472–aag2472
  • [2] Schuller, J. A.; Barnard, E. S.; Cai, W.; Jun, Y. C.; White, J. S.; Brongersma, M. L. Nature Materials 2010, 9, 193–204
  • [3] I. Staude; J. Schilling, Nature Photonics 2017, 11, 274
  • [4] Evlyukhin, A. B.; Novikov, S. M.; Zywietz, U.; Eriksen, R. L.; Chichkov, B. N. Nano letters 2012, 12, 3749–3755
  • [5] Rios, C.; Stegmaier, M.; Hosseini, P.; Wang, D.; Scherer, T.; Wright, C. D.; Bhaskaran, H.; Pernice, W. H. P. Nature Photonics 2015, 9, 725–732
  • [6] Khasminskaya, S.; Pyatkov, F.; Slowik, K.; Ferrari, S.; Kahl, O.; Kovalyuk, V.; Rath, P.; Vetter, A.; Hennrich, F.; Kappes, M. M.; Gol’tsman, G.; Korneev, A.; Rockstuhl, C.; Krupke, R.; Pernice, W. H. P. Nature Photonics 2016, 10, 727–732
  • [7] Liu, S.; Saravi, S.; Keeler, G. A.; Sinclair, M. B.; Yang, Y.; Reno, J.; Pertsch, T.; Brener, I. Nano Letters 2016, acs.nanolett.6b01816
  • [8] Schwarz, B.; Reininger, P.; Ristanić, D.; Detz, H.; Andrews, A. M.; Schrenk, W.; Strasser, G. Nature Communications 2014, 5, 4085
  • [9] Brolo, A. Nature Photonics 2012, 6, 709–713
  • [10] Schmidt, F.-P.; Ditlbacher, H.; Hohenester, U.; Hohenau, A.; Hofer, F.; Krenn, J. R. Nano Letters 2012, 12, 5780–5783
  • [11] Nordlander, P. Nature Nanotechnology 2013, 8, 76–77
  • [12] Yan, J.; Liu, P.; Lin, Z.; Wang, H.; Chen, H.; Wang, C.; Yang, G. ACS Nano 2015, 9, 2968–2980
  • [13] Yang, S. C.; Kobori, H.; He, C. L.; Lin, M. H.; Hung-Ying, C.; Li, C.; Kanehara, M.; Teranishi, T.; Gwo, S. Nano Letters 2010, 10, 632–637
  • [14] Lahiri, B.; Holland, G.; Aksyuk, V.; Centrone, A. Nano Letters 2013, 13, 3218–3224
  • [15] Chae, J.; Lahiri, B.; Kohoutek, J.; Holland, G.; Lezec, H.; Centrone, A. Optics Express 2015, 23, 25912
  • [16] Chu, M.-W.; Myroshnychenko, V.; Chen, C. H.; Deng, J.-P.; Mou, C.-Y.; García de Abajo, F. J. Nano letters 2009, 9, 399–404
  • [17] Koh, A. L.; Fernandez-Dominguez, A.; Maier, S.; Yang, J.; McComb, D. Nano letters 2011, 17, 764–765
  • [18] Coenen, T.; Vesseur, E. J. R.; Polman, A.; Koenderink, A. F. Nano Letters 2011, 11, 3779–3784
  • [19] Barnard, E. S.; Coenen, T.; Vesseur, E. J. R.; Polman, A.; Brongersma, M. L. Nano Letters 2011, 11, 4265–4269
  • [20] Coenen, T.; Van De Groep, J.; Polman, A. ACS Nano 2013, 7, 1689–1698
  • [21] van de Groep, J.; Coenen, T.; Mann, S. a.; Polman, A. Optica 2016, 3, 93
  • [22] Brenny, B. J. M.; Beggs, D. M.; van der Wel, R. E. C.; Kuipers, L.; Polman, A. ACS Photonics 2016, DOI: 10.1021/acsphotonics.6b00557
  • [23] Krasnok, A.; Makarov, S.; Petrov, M.; Savelev, R.; Belov, P.; Kivshar, Y. Proceedings of SPIE 2015, 9502, 950203
  • [24] Dmitriev, P. A.; Baranov, D. G.; Milichko, V. A.; Makarov, S. V.; Mukhin, I. S.; Samusev, A. K.; Krasnok, A. E.; Belov, P. A.; Kivshar, Y. S. Nanoscale 2016, 8, 9721–9726
  • [25] Yu, P. Y.; Cardona, M. Fundamentals of Semiconductors; Springer-Verlag Berlin Heidelberg, 2010; Vol. 1; p 778
  • [26] Cao, L.; Nabet, B.; Spanier, J. E. Physical Review Letters 2006, 96, 157402
  • [27] Sinev, I.; Iorsh, I.; Bogdanov, A.; Permyakov, D.; Komissarenko, F.; Mukhin, I.; Samusev, A.; Valuckas, V.; Kuznetsov, A. I.; Luk’yanchuk, B. S.; Miroshnichenko, A. E.; Kivshar, Y. S. Laser and Photonics Reviews 2016, 10, 799–806
  • [28] Evlyukhin, A. B.; Bozhevolnyi, S. I. Physical Review B 2015, 92, 245419
  • [29] Miroshnichenko, A. E.; Evlyukhin, A. B.; Kivshar, Y. S.; Chichkov, B. N. ACS Photonics 2015, 2, 1423–1428
  • [30] Markovich, D. L.; Ginzburg, P.; Samusev,A. K.; Belov, P. a.; Zayats, a. V. Optics Express 2014, 22, 10693–10702
  • [31] Krasnok, A.; Li, S.; Lepeshov, S.; Savelev, R.; Baranov, D. G.; Alú, A. Phys. Rev. Applied 2018, 9, 014015
  • [32] Sinev, I. S.; Komissarenko, F. E.; Mukhin, I. S.; Petrov, M. I.; Iorsh, I. V., Belov, P. A., Samusev, A. K. Nanosystems: Physics, Chemistry, Mathematics 2018, 9, 609–613
  • [33] Krasilin, A. A.; Volodina, K.; Sukhova, A. A.; Petrov, M. I.; Zuev, D. A.; Dyachuk, V. A.; Milichko, V. A. Journal of biophotonics 2018, 11, e201700322
  • [34] Zywietz, U.; Evlyukhin, A. B.; Reinhardt, C.; Chichkov, B. N. Nature Communications 2014, 5, 3402
  • [35] Dmitriev, P. A.; Makarov, S. V.; Milichko, V. A.; Mukhin, I. S.; Gudovskikh, A. S.; Sitnikova, A. A.; Samusev, A. K.; Krasnok, A. E.; Belov, P. A. Nanoscale 2016, 8, 5043–5048
  • [36] Grinblat, G.; Li, Y.; Nielsen, M. P.; Oulton, R. F.; Maier, S. A. Nano Letters 2016, 16, 4635–4640
  • [37] Papasimakis, N.; Fedotov, V. A.; Savinov, V.; Raybould, T. A.; Zheludev, N. I. Nature Materials 2016, 15, 263–271
  • [38] Grinblat, G.; Li, Y.; Nielsen, M. P.; Oulton, R. F.; Maier, S. A. ACS Nano 2016, acsnano.6b07568
  • [39] Gongora, J. S. T.; Miroshnichenko, A. E.; Kivshar, Y. S.; Fratalocchi, A. Nat. Commun. 2017, 8, 15535
  • [40] Baryshnikova, K. V.; Smirnova, D. A.; Luk’yanchuk, B. S.; Kivshar, Y. S. Advanced Optical Materials 2019, 1801350
  • [41] Purcell, E. Physical Review 1946, 69, 674–674
  • [42] Krasnok, A. E.; Slobozhanyuk, A. P.; Simovski, C. R.; Tretyakov, S. A.; Poddubny, A. N.; Miroshnichenko, A. E.; Kivshar, Y. S.; Belov, P. A. Scientific reports 2015, 5, 12956
  • [43] Frizyuk, K.; Hasan, M.; Krasnok, A.; Alú, A.; Petrov, M. Phys. Rev. B 2018, 97, 085414
  • [44] Zambrana-puyalto, X.; Bonod, N. Physical Review B 2015, 91, 195422
  • [45] Sauvan, C.; Hugonin, J. P.; Maksymov, I. S.; Lalanne, P. Physical Review Letters 2013, 110, 237401
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description