HTML
-
To quantitatively probe plasmonic enhancement using SERS, structural parameters that strongly affect SERS EF calculations should be defined as clearly as possible26, 27. First, the geometry of plasmonic antennas should be precisely characterized, particularly the exact shape of the nanogaps. The gap region of a typical 1L MoS2-NPOM is displayed in transmission electron microscopy (TEM) cross-sectional images (Fig. 1c). The structure consists of an ultrasmooth gold film with 0.32 nm root-mean-square roughness (Supplementary Fig. S1), a 1 L MoS2 with a precisely defined thickness of 0.62 nm and a 50 nm AuNP (see Materials and Methods for sample fabrication). The atomic migration effect induced by the interaction between the Au and the MoS2 flattens the bottom region of the AuNP34, leading to an average facet size of 19.4 nm (Supplementary Fig. S2). The citrate molecules loosely covering the AuNP surface are nearly invisible to TEM, as the measured average gap distance of the nanocavity equals the thickness of the 1L MoS2 (see detailed discussions in Supplementary Fig. S2). Therefore, the MoS2-NPOM system provides a robust nanocavity that allows for the quantitative determination of the area of the "hotspot" and the gap distance.
Second, the SERS probes should lie exactly in the "hotspot", with known orientations with respect to the nanogap axis. The NPOM geometry20, 30, 34-36 combined with a two-dimensional atomic crystal probe of MoS2 can perfectly address these issues. On one hand, the MoS2 probe uniformly fills the entire gap of the nanocavity, where the local field is maximized. The exact area of the MoS2 probe (corresponding to the number of conventional probe molecules) inside the "hotspot" is naturally determined by the size of the nanocavity. On the other hand, the single atomic layer of MoS2 has a definite lattice orientation such that its out-of-plane (A1g) and in-plane (${{E}}_{2{{g}}}^1$) lattice vibrations are aligned with the vertical and horizontal local fields of the MoS2-NPOM, respectively (Fig. 1b-d). The MoS2-NPOM enhanced Raman scattering intensity ISERS can expressed as $I_{{{SERS}}} \propto | {\mathop {\alpha }\limits^ \leftrightarrow (\omega _{{R}}, \omega){\bf{E}}(\omega)} |^2$, where E(ω) is the local electric field induced by the incident field, and $\mathop {\alpha }\limits^ \leftrightarrow (\omega _{{R}}, \omega)$ is the Raman polarizability tensor. For 1L MoS2, the expression reads as follows:
$$ \mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle \leftrightarrow$}} \\ \over \alpha } ^{{E}}\left( {\omega _{{R}},\omega } \right) = \left( {\begin{array}{*{20}{c}} {\alpha _{xx}^{{E}}} & {\alpha _{xy}^{{E}}} & 0 \\ {\alpha _{yx}^{{E}}} & {\alpha _{yy}^{{E}}} & 0 \\ 0 & 0 & 0 \end{array}} \right)\,{{and}}\,\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle \leftrightarrow$}} \\ \over \alpha } ^{{A}}\left( {\omega _{{R}},\omega } \right) \\ = \left( {\begin{array}{*{20}{c}} {\alpha _{xx}^{{A}}} & 0 & 0 \\ 0 & {\alpha _{yy}^{{A}}} & 0 \\ 0 & 0 & {\alpha _{zz}^{{A}}} \end{array}} \right) $$ with $\alpha _{xx}^{{E}} = \alpha _{xy}^{{E}} = \alpha _{yx}^{{E}} = - \alpha _{yy}^{{E}}$ for the ${{E}}_{2{{g}}}^1$ mode and $\alpha _{xx}^{{A}} = \alpha _{yy}^{{A}}$ for the A1g mode37. The equations immediately imply that the Raman intensity of the in-plane phonon ${{E}}_{2{{g}}}^1$ is solely determined by the horizontal local fields (Ex and Ey), whereas the out-of-plane phonon A1g is contributed mainly by Ez because the vertical local field is dominant in the MoS2-NPOM. Based on our further analysis (see Supplementary Note S1), the horizontal (vertical) plasmonic enhancement of the MoS2-NPOM can be quantitatively obtained by the SERS EF from the ${{E}}_{2{{g}}}^1$ (A1g) phonon. Figure 1e shows the Raman spectra obtained from a 1L MoS2 on quartz, a 1L MoS2 on ultrasmooth gold film and a 1L MoS2-NPOM with a 32-nm-thick Al2O3 surface coating. An enlarged view of the spectral region around the distinct phonons is shown in Fig. 1f to better compare the peak positions and shapes. When few-layer MoS2 contacts well with the gold film, a doping effect and a local mechanical strain effect will occur38, 39. As a result, both the ${{E}}_{2{{g}}}^1$ and A1g are softened, with two new peaks appearing near the red side of the ${{E}}_{2{{g}}}^1$ and the A1g, labeled the ${{E}}_{2{{g}}}^1$′ and A1g′ modes (Fig. 1f and Supplementary Fig. S3a). Based on handedness-resolved Raman measurements (Supplementary Fig. S3), the ${{E}}_{2{{g}}}^1$′ (A1g′) is considered split from the ${{E}}_{2{{g}}}^1$ (A1g). Therefore, we sum the ${{E}}_{2{{g}}}^1$′ (A1g′) and ${{E}}_{2{{g}}}^1$ (A1g) and label them ${{E}}_{2{{g}}}^1$+${{E}}_{2{{g}}}^1$′ (A1g+A1g′) modes in the following analysis. Comparing the Raman spectra of the 1L MoS2-NPOM with that of the 1L MoS2 on an ultrasmooth gold film, we find that both the A1g+A1g′ and ${{E}}_{2{{g}}}^1$+${{E}}_{2{{g}}}^1$′ are largely enhanced in the MoS2-NPOM, whereas the spectral intensity of the former is approximately two orders of magnitude higher than that of the latter (Fig. 1e). These results guarantee the validity of using the ${{E}}_{2{{g}}}^1$+${{E}}_{2{{g}}}^1$′ and A1g+A1g′ phonons in probing the horizontal and vertical field enhancements in the MoS2-NPOM nanocavity (Supplementary Note S1). Raman imaging of the A1g+A1g′ (396 ± 10 cm−1) shows a bright spot at the position of the AuNP, confirming that the enhancement originates from the MoS2-NPOM nanocavity (the inset of Fig. 1e).
-
To optimize the SERS efficiency, the polarization of incident light should be aligned with the nanogap axis. In the NPOM geometry20, 30, 34-36, the nanoparticle interacts with its electromagnetic image induced under the metallic mirror, in analogy with a nanoparticle dimer antenna with a fixed nanogap axis along the surface normal. In our SERS measurements, the MoS2-NPOM is excited by a slightly focused 785 nm laser beam illuminating from the side with an angle of 80° from the normal to the sample (see Materials and methods and Supplementary Fig. S4a). Here, the incident light with polarization parallel to the incident plane, labeled p-polarized light, is expected to excite the nanogap plasmons efficiently. To clarify the contributions from the plasmonic gap modes excited in the SERS measurements, we performed dark-field scattering measurements on the MoS2-NPOMs using the same configuration, except with the laser replaced with white light (Materials and methods, the inset of Fig. 2a and Supplementary Fig. S4b).
Fig. 2 Plasmonic modes analysis.
a, b Measured dark-field scattering (a) and simulated absorption (b) spectra of a 1 L MoS2-NPOM with 32-nm-thick Al2O3 surface coating excited by oblique incident white light with different polarizations. The inset in a shows the excitation configuration both applied in experiments and simulations. The inset in b with red (blue) square shows the surface charge distribution at the D (M) peak excited by p- (s-) polarized light, where the nanocavity is plotted with opened geometry in order to see the charges clearly (black arrows)Figure 2a shows the dark-field scattering spectra of the same 1L MoS2-NPOM shown in Fig. 1e with un-polarized, p-polarized, and s-polarized incident light. The spectra manifest two peaks at ~785 nm, which are labeled M mode (743 nm) and D mode (798 nm). The M mode can be activated by both p- and s-polarized light, whereas the D mode can only be excited under p-polarized excitation. The abovementioned far-field behaviors of the D and M modes are also confirmed by the corresponding simulated polarization-dependent absorption spectra with a similar excitation configuration (Materials and Methods and Fig. 2b). The insets of Fig. 2b show the surface charge distributions of the D (M) mode excited by p- (s-) polarized light. For the D mode, one charge is localized on the entire bottom facet, and the opposite charge is distributed on the remaining part of the AuNP surface, as well as on the gold film surface beneath the AuNP. This charge pattern, which can only be activated by p-polarized light, suggests that the D mode is a strongly radiating antenna mode2, 40. In the case of the M mode, whereas the top surface of the AuNP is left nearly blank, both positive and negative charges are mainly concentrated on the two halves of the AuNP's bottom facet with a blank node at the centre. A similar pattern with opposite charges is induced on the surface of the subjacent gold mirror. The M mode can be understood as the lowest frequency cavity mode15, 40, 41, which appears in parallel flat terminals formed in the MoS2-NPOM geometry. This tightly confined charge distribution in the gap region makes the M mode insensitive to the refractive index change on top of the AuNP but ultrasensitive to that inside the nanocavity15.
The electric field distributions of the 1L MoS2-NPOM at 785 nm excited by p-polarized and s-polarized light are shown in Fig. 3a, b. For p-polarization excitation, both the D and M modes are excited, but the latter shows greater detuning with respect to the laser wavelength. The modes provide intense electric field enhancement, with a maximum enhancement of 240-fold, where the contribution of the M mode is negligible in the SERS process. For s-polarization excitation, the D mode cannot be excited, while the M mode can still be effectively excited, providing a 20-fold maximum electric field enhancement. The results agree well with the measured polarization-dependent SERS from the MoS2-NPOMs shown in Supplementary Fig. S4d. We determined the average vertical (horizontal) plasmonic enhancement $\bar g_z({\bar g_{xy}})$ based on the approximation that the SERS EF is proportional to the fourth order of the local electric field enhancement in the MoS2-NPOM28. Specifically, The enhancement is associated with the surface-averaged vertical or horizontal SERS EF as follows (see details in Supplementary Note S1):
Fig. 3 Local field distributions of MoS2-NPOM and effective areas of MoS2 probe.
a, b Electric field distributions at 785 nm of a 32-nm-thick Al2O3 coated 1 L MoS2-NPOM excited by p-polarized (a) and s-polarized (b) light. The MoS2-NPOM cross-sections are taken from a plane 5 nm away from the xz plane. c, d Vertical (c) and horizontal (d) components of the electric field distribution in the 1 L MoS2-NPOM with p-polarization excitation. e Fraction fz (fxy) as a function of integral radius ρz (ρxy). The inset in e shows the definition of integral radius ρ. The effective vertical (horizontal) local field area (${{\pi }}R_z^2$ ) is determined by setting ρz= Rz (ρxy=Rxy) for fz (Rz) = fxy (Rxy) = 1-e−5 (~99.3%)${{\pi }}R_{xy}^2$ $$ \overline {{{EF}}} _z = \frac{{I_{{{SERS}}}^{{A}}{{/}}S_{{{SERS}}}^z}}{{I_{{{Ref}}}^{{A}}{{/}}S_{{{Ref}}}}} \approx \left| {\bar g_z} \right|^4 $$ (1) $$ \overline {{{EF}}} _{xy} = \frac{{I_{{{SERS}}}^{{E}}{{/}}S_{{{SERS}}}^{{{xy}}}}}{{I_{{{Ref}}}^{{E}}{{/}}S_{{{Ref}}}}} \approx \left| {\bar g_{xy}} \right|^4 $$ (2) where $I_{{{SERS}}}^{{A}}$ ($I_{{{SERS}}}^{{E}}$) and $S_{{{SERS}}}^z$ ($S_{{{SERS}}}^{xy}$) are the Raman scattering intensity and the effective "hotspot" area of the A1g+A1g′ (${{E}}_{2{{g}}}^1$+${{E}}_{2{{g}}}^1$′) in the MoS2-NPOM, and $I_{{{Ref}}}^{{A}}$ ($I_{{{Ref}}}^{{E}}$) and SRef are the Raman scattering intensity and the excitation area of the A1g (${{E}}_{2{{g}}}^1$) for the same layer MoS2 on quartz. SRef is established as the collection area (~3.2 μm2) because the excitation beam is larger than the collection area in our experiments. Because the "hotspot" dominates the SERS signal, the effective local field areas $S_{{{SERS}}}^z$ and $S_{{{SERS}}}^{xy}$ should be much smaller than the collection area and can be determined based on the local field distribution of the MoS2-NPOM. Figure 3c, d shows the vertical and horizontal electric field distributions of a MoS2-NPOM with structural parameters corresponding to those in Fig. 1. The results suggest that the vertical local field is mostly localized in the area below the AuNP's circular bottom facet and decays rapidly outside the nanocavity region (Fig. 3c). In contrast, the horizontal local field is distributed almost entirely outside the facet region, with only 9.7% SERS EF contributing from the area directly below the bottom facet of the AuNP and the rest from the outer region with a longer decay length (Fig. 3d). Here, the $S_{{{SERS}}}^z$ ($S_{{{SERS}}}^{xy}$) is treated as a circular area with radius Rz (Rxy), which can be calculated by defining a fraction fz (fxy) as the ratio of the vertical (horizontal) SERS EF contributed from a circular area ${{\pi }}\rho _z^2$ (${{\pi }}\rho _{xy}^2$) centered at the nanogap region (the inset of Fig. 3e) to that from the collection area (considered as infinity):
$$ f_z \left(\rho _z\right) = {\int}_{ 0}^{{{\pi }}\rho _z^2} {\left| {E_z} \right|^4/\left| {E_{0z}} \right|^4{{d}}s/{{\int}_{ 0}^\infty}\left| {E_z} \right|^4/{\left| {E_{0z}} \right|^4{{d}}s} } $$ (3) $$ f_{xy}{{(}}\rho _{xy}{{) = }}{\int}_{ 0}^{{{\pi }}\rho _{xy}^2} {\left| {E_{xy}} \right|^4/\left| {E_{0xy}} \right|^4{{d}}s} {{/}}{\int}_{ 0}^\infty {\left| {E_{xy}} \right|^4/\left| {E_{0xy}} \right|^4{{d}}s}\hskip-14pt $$ (4) where Ez (Exy) and E0z (E0xy) are the vertical (horizontal) electric components of the local and incident field, respectively. fz (fxy) is depicted as a function of ρz (ρxy) in Fig. 3e, demonstrating that fz (fxy) saturates quickly with increasing ρz (ρxy). $S_{{{SERS}}}^z$ can be obtained by setting ρz = Rz for fz(Rz) = 1−e−5 (~99.3%, see Fig. 3e), meaning that this probe area contributes the total SERS signal. In this case, the effective radius of the vertical local field area Rz is 11 nm, only slightly larger than the bottom facet radius of the AuNP. A similar procedure can be applied to the horizontal direction, which yields an effective radius Rxy of 42 nm. After inserting these values into equations (1) and (2), $\overline {{{EF}}} _z$ ($\bar g_z$) and $\overline {{{EF}}} _{xy}$ ($\bar g_{xy}$) can be determined in a straightforward manner.
-
To obtain the maximum SERS signal, the plasmonic resonance wavelength of the nanoantenna should overlap with the incident wavelength as well as the outgoing Raman wavelength. In experiments, satisfying this condition is nontrivial, e.g., by wavelength-scanned SERS spectroscopy7, 25, 42. However, this technique requires expensive narrow-line lasers with multiple/tunable wavelengths and careful calibration schemes, preventing its widespread application. Here, we performed a plasmon-scanned SERS measurement to achieve a similar goal: plasmon resonance is gradually redshifted by successively depositing Al2O3 layers onto the sample surface to match the excitation laser (and the outgoing Raman light) at a fixed wavelength43. The effect of the Al2O3 coating on the near-field distributions of the MoS2-NPOM system are shown in Supplementary Fig. S5.
Dark-field scattering spectra of a 1 L MoS2-NPOM with different Al2O3 coating thicknesses are shown in Fig. 4a. The D and M modes are overlapping at ~700 nm without the Al2O3 coating. The modes show continuous redshift and broadening as the thickness of the Al2O3 layer increases, resulting from the dielectric screening effect induced by the high refractive index layer. As the Al2O3 thickness exceeds 20 nm, the D and M modes start splitting into two peaks. The separation originates from the larger spectral shift of the D mode relative to that of the M mode in response to the equal-thickness dielectric coating. This feature can be understood as the result of the distribution of a higher proportion of surface charges on the outer AuNP surface for the D mode than for the M mode (the inset of Fig. 2b). This coating dependence of the resonance wavelength, together with the polarization-dependent dark-field spectra shown in Fig. 2, guarantees an unambiguous identification of the plasmonic resonances in experiments. For comparison, we calculated the corresponding far-field absorption spectra of the MoS2-NPOM with p-polarization excitation (Fig. 4b, see the corresponding scattering spectra in Supplementary Fig. S6). The results show similar redshift and broadening behaviors of the D peak, followed by the appearance of the M peak on the blue side of the D peak. The measured and simulated scattering spectra with an Al2O3 coating exceeds 50 nm; the results of further analysis are shown in Supplementary Fig. S6.
Fig. 4 Plasmon-scanned SERS measurement.
a, b Measured dark-field scattering (a) and simulated absorption (b) spectra of a 1 L MoS2-NPOM with the Al2O3 coating thickness varied from 0 nm to 50 nm. The gray dotted lines in a and b show the laser line at 785 nm. c Intensity map of the SERS spectra of the 1 L MoS2-NPOM as the Al2O3 coating thickness increases from 8 nm to 102 nmSERS measurements were performed after each Al2O3 coating under p-polarization excitation below 4 μW/μm2 to avoid damage to the sample or possible nonlinear effects. SERS spectra from the same 1L MoS2-NPOM are mapped in Fig. 4c, showing a pronounced enhancement of the spectral intensity when the wavelength of the plasmon resonance overlaps that of the excitation laser. For each Al2O3 coating, we can determine the corresponding peak position of the D mode, $\overline {{{EF}}} _z$ and $\overline {{{EF}}} _{xy}$ via the dark-field scattering spectrum, equations (1) and (2), respectively. Then, $\overline {{{EF}}} _z$ and $\overline {{{EF}}} _{xy}$ can be plotted as a function of the resonance wavelength of the D mode, as shown in Fig. 5a. $\overline {{{EF}}} _z$ increases rapidly as the D peak wavelength approaches the excitation wavelength and then reaches the maximum value of 1.5 × 108 when the D peak wavelength is between the incoming laser wavelength (785 nm) and the outgoing Raman wavelength (~810 nm). $\overline {{{EF}}} _z$ drops quickly when the D peak is red-detuned to the Raman wavelength. This D-peak-scanned $\overline {{{EF}}} _z$ profile almost follows the far-field shape of the D mode (Fig. 5a). The corresponding results in the horizontal direction show similar behavior, but the maximum $\overline {{{EF}}} _{xy}$ is 6.3 × 104, which is ~2400-fold weaker than that of $\overline {{{EF}}} _z$. Similar performances are also observed in the statistical values of $\overline {{{EF}}} _z$ and $\overline {{{EF}}} _{xy}$ in response to the D peaks (Fig. 5b) obtained from 17 different 1L MoS2-NPOMs (Supplementary Fig. S7). As the gap distance in our MoS2-NPOM system is precisely determined by the thickness of the MoS2, the error bars of $\overline {{{EF}}} _z$ mainly originate from the derivation of the diameter and the bottom facet size from different measured AuNPs (Supplementary Fig. S2 and Supplementary Fig. S8). It should be noted that the maximum SERS EFs of several 1L MoS2-NPOMs are not included in Fig. 5b because their D peaks become invisible in the dark-field spectra as the Al2O3 coating thickness exceeds ~70 nm (see details in Supplementary Fig. S6).
Fig. 5 Vertical and horizontal SERS EFs in response to the D peak position.
a and$\overline {{{EF}}} _z$ of a 1 L MoS2-NPOM in response to the resonance wavelength of its D mode, along with its dark-field (DF) scattering spectrum at 36-nm-thick Al2O3 coating (where the$\overline {{{EF}}} _{xy}$ reaches maximum). The light orange (blue) shadow area represents the far-field shape of the D mode (M mode) by fitting the DF scattering spectrum with Lorentz curves. b Statistical average of the$\overline {{{EF}}} _z$ and$\overline {{{EF}}} _z$ as a function of the D peak position, obtained from 17 individual 1 L MoS2-NPOMs. The errors bars represent the standard deviations of the$\overline {{{EF}}} _{xy}$ and$\overline {{{EF}}} _z$ calculated by averaging the data points of the D peak dependent$\overline {{{EF}}} _{xy}$ and$\overline {{{EF}}} _z$ with 5 nm wavelength a step$\overline {{{EF}}} _{xy}$ Classical electromagnetic simulations predict a smooth curve for the Al2O3 thickness-dependent $\overline {{{EF}}} _z$ (Supplementary Fig. S9), whereas curves in Fig. 5a contain one or several dips near the peak region. A similar feature is observed in other MoS2-NPOMs (Supplementary Fig. S7) and the previous wavelength-scanned SERS measurement results7, 25, 42. In conventional single-molecule SERS experiments, large-intensity fluctuation is a typical observation, due to the possible random movements or photodamage of the probe molecules around the 'hotspot'. In contrast, benefiting from the higher photodamage threshold of the two-dimensional atomic crystal and its strict lattice arrangement in the 'hotspot', the SERS spectra of our MoS2-NPOM system exhibit high stability and repeatability over time, demonstrating 4.3% intensity fluctuation over 150 minutes (Supplementary Fig. S10a). Therefore, we could exclude the possibility of large-intensity fluctuation during the SERS measurements. Power-dependent SERS measurements show a linear relationship between the incident power and the SERS intensity (Supplementary Fig. S10b), indicating that the process is unlikely to be a stimulated Raman scattering one. Another possible explanation is associated with the phonon nature of the Raman scattering, which may also involve an optomechanical mechanism32, 33, 44. Further studies are required to clarify this feature, which is beyond the scope of the current work.
Plasmon-scanned SERS measurement can obtain the maximum SERS enhancement of an individual MoS2-NPOM, thus enabling quantitative probing of the limits of plasmonic enhancement in nanogaps. We repeated these measurements on six individual bilayer (2 L) MoS2-NPOMs and four individual trilayer (3 L) MoS2-NPOMs to obtain their maximum SERS EFs. The largest $\overline {{{EF}}} _z$ values for 1 L, 2 L, and 3 L MoS2 probes reach up to 5.1 × 108, 1.7 × 108, and 8.5 × 106, respectively (Supplementary Fig. S7). The statistical average of the maximum $\overline {{{EF}}} _z$ as a function of the gap distance is plotted in Fig. 6, showing that $\overline {{{EF}}} _z$ rapidly increases as the gap distance narrows. For comparison, we first calculated the total SERS EFs using the E4 model (see Materials and method, and Supplementary Fig. S9), which can approximate the $\overline {{{EF}}} _z$ due to the dominated field enhancement in the vertical direction. The measured and simulated $\overline {{{EF}}} _z$ are in good agreement at gap distances of 1.24 nm and 1.86 nm. Thus, our measured $\overline {{{EF}}} _z$ can be well described by the electromagnetic enhancement without introducing the concept of chemical enhancement applied in most SERS experiments27 and thereby demonstrating the reliability of our designs in probing the plasmonic enhancement. However, this precise prediction starts to overestimate the measured plasmonic enhancement at narrower gap distances. As the gap distance decreases to 0.62 nm, the measured $\bar g_z$ ≈ 114 ($\overline {{{EF}}} _z$ = 1.7 × 108) becomes ~38.4% (~6.9-fold) lower than the $\bar g_z$ ≈ 185 ($\overline {{{EF}}} _z$ = 1.17 × 109) predicted by the E4 model. Because the gap distance and the orientation of the probe in our antenna system are well-defined and robust, the reduction of the measured plasmonic enhancement in the subnanometer gap is most likely caused by the emergence of previously predicted quantum mechanical effects18, 19, 22. Note that these values are the average values over the total "hotspot" area. The maximum out-of-plane plasmonic enhancement $g_z^{{{max}}}$ (at the "hottest" position) can be obtained if the field distribution is precisely known. Taking the geometry variations in the experiments into account, the MoS2-NPOM provides a satisfactory linear relationship between the maximum and the surface-averaged SERS EFs: ${{EF}}_z^{{max}} = 3.24 \times {{10}}^8 + 2.5\overline {{EF}} _z$ (Supplementary Fig. S8). Therefore, the measured ${{EF}}_z^{{max}}$ ($g_z^{{{max}}}$) is evaluated to be as large as 4.93 × 108 (148).
Fig. 6 Measured and simulated maximum $\overline {{{EF}}} _{{z}}$ in response to the gap distance.
The error bars represent the standard deviations of the measured maximum averaged from seventeen 1 L, six 2 L and four 3 L MoS2-NPOMs. The simulated SERS EFs are calculated by the E4 model and two-study model (TSM), respectively (see Materials and Methods, Supplementary Note S2). The simulated SERS EFs based on the pure classical electromagnetic theories match well with the measured SERS EFs as the gap distance is no less than 1.24 nm, while significantly overestimate the field enhancement at shorter gap distance$\overline {{{EF}}} _{{z}}$ To further clarify the issue regarding the quenching of the field enhancement, we also introduced a two-study model (TSM) to calculate the SERS EFs (Fig. 6) of the A1g phonon in the MoS2-NPOMs and the same layer MoS2 on quartz (see Materials and Methods, Supplementary Note S2). Using a polarization current as the source at the Raman frequency, the emission enhancement of the MoS2-NPOM can be fully captured45. A key difference between molecule vibrations and lattice phonons is the coherence of the Raman signal from different locations, which is characterized by the correlation length. This effect can be qualitatively considered in the TSM by setting the diameter of the MoS2 probe equal to the correlation length, as the polarization current is coherent within a MoS2 sheet of finite size. Upon assuming the correlation length of the A1g modes in 1 L, 2 L, and 3 L MoS2 to be 24 nm, 28 nm, and 34 nm, the SERS EFs obtained by the TSM match those predicted by the E4 model very well (Fig. 6). These assumed correlation lengths are comparable with the measured value of ~30 nm for optical phonons in graphene46. Additionally, the decrease in the correlation length with the decrease in the number of layers is reasonable because the 1 L MoS2 may have local wrinkles that disturb the coherence of phonons. Therefore, we can conclude that as the gap distance exceeds 1.24 nm, the two models based purely on classical Maxwell's descriptions can predict the behaviors of the measured plasmonic enhancement quite well, but they start to overestimate the field enhancement at narrower gap distances. These performances regarding the narrow gap match well with the prediction of the quantum-corrected model, suggesting the probable emergence of electron tunneling across the 1 L MoS2 layer.