key: cord-0976765-nmz23lrh
authors: Nagarajan, Anitharaj; Panchanathan, Aruna Priya; Chelliah, Pandian; Satoh, Hiroaki; Inokawa, Hiroshi
title: FDTD Study on Evolution of Trimer Silver@Silica Nanospheres to Dimer for SERS Characteristics
date: 2021-10-13
journal: Plasmonics
DOI: 10.1007/s11468-021-01554-y
sha: 8b43ca776b5f60313fd7597dc7e4cb51b503fc39
doc_id: 976765
cord_uid: nmz23lrh

Light enhancement occurs strongly within the plasmonic clusters by interaction with surface plasmons. Surface-enhanced Raman spectroscopic (SERS) characteristics of a series of silver@silica trimer core–shell (CS) nanosphere (NS) clusters are investigated in this paper. It is significant to understand the electric field (EF) enhancement mechanism behind the SERS technique. The effect of symmetry breaking is studied for the series starting from the highly symmetric trimer cluster and transformed to linear dimer geometry which progresses through the gradual reduction in the vertex NS. The optical activity such as the evolution of LSPR peak is discussed, the formation of hot spots is demonstrated and the strength of the local EF enhancement is calculated and correlated with the plasmon dipolar modes by using plasmon hybridization theory to understand the underlying physical concepts.

Localized surface plasmon resonance (LSPR) effects observed in metallic nanostructures due to the collective oscillations of conductive electrons under irradiation draw increasing interest in the applications such as solar cells, plasmonic waveguides, LSPR sensors and surface-enhanced Raman spectroscopy (SERS) [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] . In SERS technique, which is a more promising choice for COVID-19 detection than the time-consuming PCR, the significant challenge is to understand the interplay between the structure of the nanoparticle and its optical response [11] . SERS characteristics depend on the formation of hot spots, which are strong electric field (EF) confinement, occurring in the gap between the two individual nanospheres, direction of polarization and the strength of the EF [12, 13] . The features of a highly sensitive and rapidly analyzing SERS substrate should produce an enhancement factor greater than 10 5 , it should be stable under the effect of oxidation and should not be degraded under ambient environment conditions [14] .

Plasmonic molecules are small clusters constructed by the aggregation of individual metallic nanospheres in controlled manner [15] . Plasmonic clusters with metals encapsulated by a dielectric layer enable tuning of wavelength in the range and prevent the effect of oxidation and instability under ambient environment condition [10, 16] . Dimer with nearly touching nanospheres is the mostly studied plasmonic molecules, due to its capability to yield strong hot spots [17] . The complex cluster, trimer, is more interesting as they not only confine light in the gaps, but may open the window into numerous novel fascinating phenomena by breaking its symmetry [18] [19] [20] .

In our previous work, the study on the optimization and relation between the SERS enhancement characteristics and structure of silver@silica trimer cluster composed of nanospheres with equal sizes have been determined by symmetry breaking. The highly symmetric trimer arrangement is transformed to linear trimer by gently opening the vertex angle. The effect of polarization, evolution of LSPR peak and optimization for a highly sensitive SERS substrate have been analysed and correlated with the plasmon dipolar modes using plasmon hybridization theory for understanding the background mechanism. The role of shell such as tuning of wavelength and twofold enhancement has also been discussed [21] .

The current work focuses on symmetry breaking of the equilateral trimer cluster through a different scenario. By gently reducing the vertex nanosphere size, the symmetric trimer is transformed to linear dimer. The SERS enhancement factor is directly proportional to the fourth power of the ratio of the EF intensity around the NS to the EF intensity of the incident light (|E|/|E o |) 4 [22] and could be determined by evaluating enhancement factor of nanocluster ((|E|/|E o |) 2 ). For this current series, the local EF enhancement is calculated, the polarization selective LSPR evolution is discussed and the formation of hot spots is demonstrated and correlated with the plasmon dipolar modes.

The electromagnetic simulation for investigating the SERS characteristics is performed by 3D-FDTD method using the open-source software package called MEEP [23] . The size of the computational cell is kept as 480 nm 3 , the Gaussian source of central frequency, 0.9 PHz and width, 1.2 PHz. The boundary condition for x, y and z directions is made by perfectly matched layers of thickness 50 nm. The resolution of the computation is 15 cells/nm with the spatial step of 0.1 nm. Lorentz-Drude expression from Rakic's work is used to model Ag [24] . The refractive index of the silica and air used is 1.45 and 1, respectively. The FDTD simulation method is verified with the theoretical result by comparing the optical spectra of the monomer between the simulated and the Mie's theory.

The symmetric trimer of D 3h symmetry group consists of three nearly touching (gap = 2 nm) CS NS of radius, R CS = 30 nm (radius of silver, core is R C = 25 nm and thickness of silica is 5 nm, i.e. R CS / R C = 1.2). The symmetry breaking is done by gently reducing the vertex radius, R 1 = 30 nm, with R CS / R C maintained at 1.2 and the base radius, R 2 , and gap size are constant at 30 and 2 nm, respectively leading to the transformation of C 2v symmetry. Further reduction in size (R 1 = 0 nm) leads to vanishing of vertex NS and formation of linear dimer (D ∞h symmetry) (Fig. 1 ).

The local electric field enhancement of the trimer with R 1 = 30, 25, 15 and 0 nm is calculated under both transverse and longitudinal polarizations and the plasmonic spectra are shown in Fig. 2 . The D 3h symmetry with R 1 = 30 nm shows EF enhancement with degenerative nature at two peaks. This degeneracy of the D 3h symmetry could be discussed by exploring the structure of its plasmon mode. According to plasmon hybridization theory, D 3h group consists of six in-plane modes in which three are bonding and three are anti-bonding modes. The correlation table shown in Fig. 3 indicates three bonding dipolar modes of D 3h symmetry, two 

The resizing of vertex particle leads to the lift in the degeneracy of the plasmon modes and the corresponding transformation of the three bonding modes is clearly represented in correlation table for each trimer configuration. The polarization selective excitation of the mode and the associated electric field enhancement could be explained by using plasmon modes shown in Fig. 3 . The A 1 mode of the C 2v configuration is originated from the lift in the dipole orientation of transverse E′ modes. For the case of R 1 = 25 and 15 nm under transverse polarization, the A 1 bonding mode gradually blue shifts and the electric field intensity slightly diminishes corresponding to the decrease in its dipole moment. When R 1 is reduced to 0 nm, it corresponds to linear dimer of D ∞h symmetry, the A 1 mode is transformed to the null mode of dimer cluster, the LSPR further blue shifts and electric field is decreased well. The e-field of null mode is similar to the intensity of the monomer, because all the localized charges are on the vertex CS NS with a zero size, while the charge amplitude of the base CS NSs is essentially zero.

On the other hand, under longitudinal polarization, for C 2v configuration, the longitudinal B 2 bonding modes originated from A 2 ′ modes of D 3h symmetry are responsible for the strong peak in clusters with R 1 = 25 and 15 nm, with interestingly no change in the position of the spectral peak. For D ∞h symmetry, transformation to bright Σ + u mode occurs which results in the stronger peak. The antibonding peak is gradually diminishes and vanishes in both the cases of transverse and longitudinal polarization, when the configuration is altered from D 3h symmetry to D ∞h symmetry. The surface charge excitation is not enough to enhance the e-field and only a small bonding peak appears for D ∞h symmetry, enabling polarization selective switching of e-field enhancement in dimers. Figure 4 shows the plot of spatial distribution for CS NS trimers to demonstrate the generation of e-field hot spots for understanding the e-field localization in different trimer configurations. For longitudinal and transverse polarization in D 3h symmetry, hot spots rises at the gap between three NS corresponding to the longitudinal E′ mode and gap between the base NS corresponding to the dipole orientation of transverse E′ mode, respectively. As the degeneracy is lifted, e-field selectively diminishes with respect to the direction of polarization. In C 2v symmetry, for longitudinal polarization, hot spot arises between the base NSs corresponding to the B 2 mode with decreases e-field intensity and for transverse polarization exhibits apparent diminishing of hot spots. In D ∞h symmetry, hot spots are switched on and off corresponding to the bright Σ + u mode in longitudinal polarization and null mode in transverse polarization.

The electric field enhancement and LSPR mode of the trimer series are analysed at the Z centre , i.e. same z-position of centre of the three CS NSs and discussed in the previous sections. The effect of the trimer series on the electric field enhancement and the LSPR wavelength with respect to the Fig. 5a and 5b, the electric field enhancement and LSPR wavelength of the trimer series are calculated at Z centre (solid lines) and Z top (dotted lines) for longitudinal, E x (square), and transverse, E y (circle), polarizations. In the case of D 3h symmetry configurations, there is no change in the geometry of the cluster with respect to the change in the z-position at the top surface; therefore, no change is observed for the LSPR mode and field enhancement. However, the geometry of the C 2v symmetry (R 1 = 15 and 25 nm) is broken, i.e. the position of vertex NS is shifted when the top surface of the vertex NS is aligned with the same z-position of the two base NSs. For Ex polarization, the effect of change in z-position does not affect the LSPR wavelength as the plasmonic effect is caused by the bonding mode, B 2 mode, that occurs between the base NSs which lies along x direction; however, a small increase in the field enhancement is observed for the Z top . For E y polarization, predominant change is observed for the C 2v symmetry due to coupling of vertex NS with base NSs along y direction. The LSPR peak is blue shifted and the field enhancement shows a high difference. In the case of Z top , C 2v symmetry is broken and does not obey the plasmon hybridization theory. In the case of D ∞h symmetry, there is no change in the geometry with respect to the Z top as the vertex NS is absent in this linear dimer structure; therefore, no change is observed for the LSPR mode and field enhancement similar to the case of D 3h symmetry.

In this paper, the SERS characteristics for the transformation of trimer arrangement of silver@silica CS NS cluster to dimer arrangement are numerically investigated for the first time and correlated with the plasmon modes for understanding the background physical concepts. The e-field localization corresponding to the dipolar orientation of the plasmon modes is clearly evident by the hot spot generation as represented in the spatial map. The D 3h symmetry exhibits high enhancement with degenerative nature. A polarization selective decrement is observed as the degeneracy in the plasmon dipolar mode is lifted for the C 2v and D ∞h symmetry. The optical behaviour of different series is shaped by the structure of the plasmon modes. The origin of the plasmonic behaviour of the dimer is observed from plasmon mode evolution from trimer to dimer. The effect of z-position on centre of vertex NS shows the significance of arrangement of NSs on symmetry breaking. This study paves way for designing SERS substrate with complex clusters such as tetramers and pentamers to study their SERS behaviour by breaking their symmetry into different geometrical configurations for polarization selective light coupling. 

Design of Polarization Independent SERS Substrate with Raman Gain Evaluated Using Purcell Factor