Open access peer-reviewed chapter - ONLINE FIRST

Tip-Enhanced Raman Spectroscopy of 2D Semiconductors

By Mahfujur Rahaman and Dietrich R.T. Zahn

Submitted: February 16th 2021Reviewed: August 5th 2021Published: October 6th 2021

DOI: 10.5772/intechopen.99817

Downloaded: 1


Two-dimensional (2D) semiconductors are one of the most extensively studied modern materials showing potentials in large spectrum of applications from electronics/optoelectronics to photocatalysis and CO2 reduction. These materials possess astonishing optical, electronic, and mechanical properties, which are different from their bulk counterparts. Due to strong dielectric screening, local heterogeneities such as edges, grain boundaries, defects, strain, doping, chemical bonding, and molecular orientation dictate their physical properties to a great extent. Therefore, there is a growing demand of probing such heterogeneities and their effects on the physical properties of 2D semiconductors on site in a label-free and non-destructive way. Tip-enhanced Raman spectroscopy (TERS), which combines the merits of both scanning probe microscopy and Raman spectroscopy, has experienced tremendous progress since its introduction in the early 2000s and is capable of local spectroscopic investigation with (sub-) nanometer spatial resolution. Introducing this technique to 2D semiconductors not only enables us to understand the effects of local heterogeneities, it can also provide new insights opening the door for novel quantum mechanical applications. This book chapter sheds light on the recent progress of local spectroscopic investigation and chemical imaging of 2D semiconductors using TERS. It also provides a basic discussion of Raman selection rules of 2D semiconductors important to understand TERS results. Finally, a brief outlook regarding the potential of TERS in the field of 2D semiconductors is provided.


  • TERS
  • Raman spectroscopy
  • 2D semiconductors
  • TMDC
  • MMC
  • plasmonics
  • nanoscale chemical imaging
  • local heterogeneities

1. Introduction

The ability of isolating stable, atomically thin monolayers of layered materials stimulated a new field of atomic-scale interface physics, with tremendous potential for novel quantum optoelectronic applications [1]. Graphene, the first isolated 2D material consisting of an atomically thin carbon sheet provides much stronger mechanical strength and electrical or heat conductivity compared to its bulk counterpart graphite [2, 3, 4]. Thus, graphene was predicted to have great impact on devices with ballistic charge transport or quantum anomalous Hall effect [5]. However, the lack of a bandgap coupled with challenges and consequences associated to the attempts of bandgap opening hinders its application as an active material in semiconducting devices [6]. In this context, 2D semiconductors have attracted significant attention owing to their suitable bandgap required for optoelectronic applications. Some of these semiconductors already possess carrier mobilities that can outperform existing Si CMOS technology at the scaling limit, although they are affected by some other limiting factors [7]. Therefore, extensive research is going on to gain both fundamental understanding of these materials and to explore new 2D semiconductors for potential applications.

One of the novel aspects of 2D semiconductors is that their optical properties can be tuned as a function of layer number [8]. The most famous 2D semiconductors in this respect are the transition metal dichalcogenides (TMDCs). For example, in bulk TMDCs are indirect bandgap semiconductors, which can be tuned upward ∼1 eV with decreasing layer number down to monolayer. More importantly, the nature of the bandgap also changes from indirect to direct at the monolayer thickness [9]. Density functional theory (DFT) predicts that the direct excitonic transition energy at the Brillouin zone (BZ) K point in TMDCs hardly shows any dependence with respect to the layer thickness [9]. However, the indirect transition along the ΓQ(valence band – conduction band) direction (the smallest energy gap in the bulk, i.e. the bandgap) increases monotonically as the layer number decreases. Therefore, at monolayer thickness the direct transition becomes smaller than the indirect gap and hence the TMDC becomes a direct bandgap semiconductor. Black phosphorous (BP) also exhibits a similar behavior (from 0.3 eV in bulk to 1.5 eV in monolayer) whereas the nature of the bandgap remains direct due to strong interlayer coupling [10]. Other 2D semiconductors such as, metal monochalcogenides (MMCs) are direct bandgap semiconductors in bulk and turn into indirect bandgap semiconductor at some thicknesses as the layer number decreases [11, 12].

Even though these semiconductors are few atoms thick, some of them can absorb up to 15% of light in the visible range due to strong light-matter interaction [13]. The dramatically reduced dielectric screening in the out-of-plane direction creates strongly bound excitons with binding energies in the range of few hundreds of meV [14, 15]. Therefore, their optical response is dominated by excitons or multi-particle complexes at room temperature (RT) [16, 17]. Atomically flat dangling bond free surfaces of these materials are free from career scattering caused by surface roughness, which can lead to high performance optoelectronic devices with large on–off ratio (∼108), photoresponsivity (∼ 107 mAW−1), and career mobility (103 cm2V−1 s−1) at RT [7, 18, 19]. Both few layer BP and InSe based field effect transistors show comparable career mobility (103 cm2V−1 s−1) at RT suggesting they can outperform Si based CMOS at the scaling limits [20]. However, BP suffers from poor stability in air [21] and the polar nature of InSe leads to layer dependent Fröhlich interaction [22] and thus optical phonon limited career mobility.

For all the extraordinary optoelectronic properties of 2D semiconductors, the influence of local heterogeneities such as edges, grain boundaries, defects, strain, doping, chemical bonding etc.are very crucial due to the strong Coulomb screening in these systems [23, 24]. Therefore, a powerful analytical technique, which is capable of studying chemical structures, phase, stress, impurities, molecular interaction the nanoscale, is essential to understand the underlying physics important for device applications. Raman spectroscopy meets all criteria when it comes to measure these heterogeneities at the micro-scale. It has been frequently applied as a non-destructive and straightforward to measure technique to study 2D semiconductors over the last decade [25]. However, the major disadvantages of this technique are its low sensitivity and optical diffraction limited spatial resolution. While typically one in every 105–107 scattered photons is counted for Raman scattering, the typical spatial resolution of conventional micro-Raman scattering is in the range of 0.5 to 1 μm. Hence, both of these characteristics limit Raman spectroscopy from acquiring the local heterogeneous information stated above. Therefore, tip-enhanced Raman spectroscopy (TERS) has become a popular analytical technique for studying 2D semiconductors in recent years [26]. It combines the benefits of both scanning probe microscopy (SPM) and Raman spectroscopy, thus permits a very good spatial resolution and high sensitivity for local spectroscopic investigation and chemical imaging. In this chapter, we start with discussing the principle of TERS by explaining the technical aspects of it. Then, recent developments and applications of TERS of 2D semiconductors are discussed. Finally, conclusions and future directions of TERS of 2D semiconductors are addressed.


2. Principle of TERS

Raman scattering is at the core of TERS, which is light inelastically scattered by elementary excitations such as vibrational modes (phonons) in the material probed. Consequently photons are emitted that are shifted in energy with respect to the energy of the exciting light. The shifts are directly correlated to the phonon frequencies of the material. This effect was first observed by C. V. Raman in 1928 and thus named after him [27]. The fundamentals and the various aspects of Raman scattering are comprehensively covered the book series “Light Scattering in Solids I-IX” [28, 29]. Here, we would like to point out a few facts, which are relevant for TERS. First, the Raman process in matter involves the instantaneous absorption of a photon, excitation of the electronic system, scattering of electron or hole by a vibrational mode or phonon, and relaxation of the excited state accompanied by the emission of a photon. This higher order process has an extremely weak Raman scattering cross section (typically 10−26–10−31 cm2) [30, 31]. The intensity of the inelastically scattered Raman light is thus orders of magnitude weaker than the elastically scattered Rayleigh light, which imposes the need of strong reduction of the elastically scattered light by the Raman spectrometers. The weakness of Raman scattering is especially critical if the amount of matter and/or the available scattering volume are very small as is the case for low-dimensional structures such as quantum dots or 2D semiconductors. For such samples, enhancement mechanisms are employed. An important one is resonance enhancement, for which the exciting light is tuned to a real electronic transition of the sample leading to strong optical absorption and consequently stronger Raman scattering (typically two orders of magnitude). However, this may not be sufficient for low-dimensional systems and thus further enhancement using surface-enhanced Raman scattering (SERS) is required. SERS has seen tremendous increase in applications in recent years [32]. Compared to conventional Raman cross sections SERS cross section can be up to ten orders of magnitude larger allowing e.g.single molecule detection [33, 34, 35]. SERS uses metallic nanostructures to create a very strong electric field via localized surface plasmon resonances (LSPR) or lightening rod effects also known as plasmonic effects [36]. Besides a local electromagnetic (EM) field enhancement, photo-driven charge transfer (CT) mechanisms can also lead to an additional enhancement factor of up to 103 [34]. Despite the remarkable improvement compared to conventional Raman scattering, SERS still suffers from the optically diffraction limit and is restricted to resolve spatial chemical information beyond 200 nm [26]. TERS can also be esteemed as the ultimate case of SERS with the SPM tip acting as the single metallic nanostructure enhancing the Raman cross section. The photon emitted in the Raman process is shifted in energy with respect to the photon energy of the exciting light towards higher (anti-Stokes) or lower energy (Stokes) by the energy of the elementary excitation (phonon) depending on whether a phonon is absorbed or emitted. The Stokes-shifted part of the Raman spectrum is predominantly probed in Raman experiments. As the Raman shift directly reflects the energy of the elementary excitation (phonon), it can directly be used to identify materials. Moreover, the intensity (under non-resonant conditions) is related to the number of scatterers (scattering volume), while the energy position and the lineshape can be influenced by parameters such as temperature, strain, doping, and crystallinity. Finally, the Raman scattering process obeys symmetry selection rules, which are represented by the so-called Raman tensors (derivative of the polarization tensors with respect to the phonon elongation). By choosing the polarization of the incident and scattered light it is possible to probe distinct components of the Raman tensor. The Raman selection rules can be altered by e.g.electric fields, strain, and by reduced dimensionality as well as by the presence of the SPM tip inducing strong electric field gradients.

TERS works using the same principal as SERS, i.e. it utilizes the plasmonic enhancement around a sharp metallic tip. However, unlike SERS it employs a single SPM tip to collect local phonon information by taking advantage of the lateral resolution of SPM [37, 38]. When excited by a suitable photon energy, the collective oscillation of the conduction band electrons in the sharp metallic tip amplifies and confines the optical field in the vicinity of the tip apex. The resulting enhancement of the EM field then produces Raman scattering from a nanoscopic volume of the sample under the apex. The general schematic of the TERS mechanism is shown in Figure 1a. For simplification the tip apex can be assumed as shown in the schematic to be a metal sphere acting as the metal nanostructure prerequisite for the plasmonic enhancement. When an EM wave with Efield parallel to the tip long axis excites the tip, a strong plasmonic field (also known as local electric field or simply local field) is created around the apex since the metal nanostructure can facilitate both LSPR and lightening rod effect. A finite element method (FEM) simulation of the plasmonic field enhancement around the Au tip apex with a radius of 10 nm at 638 nm excitation is presented in Figure 1b. The local field amplitude is enhanced by a factor of 5 at the tip apex. Since, the Raman scattering is approximately proportional to the fourth power of the excitation field, TERS can thus amplify the local Raman information by a factor of 625 under these simulated conditions.

Figure 1.

Schematic representation of TERS (a) and FEM simulation of the plasmonic field enhancement at the tip apex for an Au tip (b).

The lightning rod effect is introduced by the anisotropy of the tip geometry (the larger the anisotropy, the stronger the field enhancement) and thus independent of the excitation wavelength. However, the LSPR is created due to the collective oscillation of conduction band electrons at the metal surface. Therefore, several factors influence the LSPR energy and oscillator strength, such as material, size, shape, and dielectric interface [36, 39, 40]. Gold and silver are the two most widely used noble metals for TERS tips. Both metals reveal small dielectric loss, stability in air, and tunability of the LSPR in the visible spectrum. Importantly, they both have some advantages and disadvantages, which become critical deciding factors for TERS applications. For example, silver shows a more pronounced plasmonic effect but gold provides a better environmental stability. Thus, the latter makes gold the more popular material for TERS tips. Additionally, for the same size and shape of the apex, the LSPR of a tip made of Ag lies at a higher energy than that of an Au tip [36]. Therefore, Ag tips are more commonly in use for green excitation and Au tips are more suitable for red and near-IR TERS applications.

The TERS schematic shown in Figure 1a is also known as conventional TERS, in which the local electric field is created by the tip dipole oscillation. The TERS sensitivity (i.e.enhancement factor, EF, and spatial resolution, SR) is moderate in a conventional TERS geometry. Hence, a usual practice of increasing TERS sensitivity multifold is to introduce a metal substrate, which then creates an image dipole inside the substrate [41]. Thus, the combined tip-substrate system acts as a dimer and enhances the sensitivity to a great extent such that the EF can reach more than three order of magnitude higher values than that in the conventional TERS [42]. This way of enhancing the TERS sensitivity is also known as gap-mode TERS configuration.

Figure 2a displays the schematic of the gap-mode TERS configuration. The strength of the dimer coupling and consequently the gap-mode TERS sensitivity depend on the tip-substrate distance d. A comparison of generated local electric field distribution around an Au tip apex in conventional TERS and gap-mode TERS is presented in Figure 2b and c. As can be seen, the enhancement of the local electric field amplitude in the case of conventional TERS is almost 20 times smaller than that in the gap-mode TERS configuration. The spatial resolution of gap-mode TERS is also superior to that of the conventional TERS. The full width at half maximum (FWHM) of the local electric field Elocdistribution around the tip apex in conventional TERS is proportional to (∼ 1.345(R + d)), where Ris the tip radius and dis the distance between the sample surface and the tip apex [41]. However, in the case of gap-mode TERS the FWHM of the Elocdistribution is given by 2Rdas shown in Figure 2d [37]. Since the TERS intensity is proportional to the fourth power of Eloc, the TERS SR also turns out to be much smaller than that in the conventional geometry. The analysis of SR can be performed using the schematic shown in Figure 2a. The potential drop between the tip sphere and the image sphere can be written as ΔV=Elocd; whereas the potential difference between two sites in the absence of the these two spheres can be expressed as ΔV=E02R+d, E0is the incident field. Since, the spheres are equipotential, we can thus write

Figure 2.

Schematic of the gap-mode TERS configuration (a). Comparison of local electric field distribution between conventional (b) and gap-mode TERS (c). Local electric field amplitude (d) and TERS enhancement (e) in a gap-mode configuration shown in (c) calculated for a tip-substrate distance of 1 nm. TERS EF (f) and SR (g) in a gap-mode TERS configuration as a function of tip diameter. The tip-substrate distance was kept at 1 nm. Figure (a) to (e) reproduced from ref. [37] with permission from the Royal Society of Chemistry and (f) and (g) from ref. [38].


For a certain tip-substrate geometry both Rand dare fixed at the smallest distance from the tip apex to the sample surface being d, which is along the center of the electric field lines. Therefore, the lateral offset of the electric field from the center can be written as [38].


For the TERS intensity Eq. (2) becomes


At x = 0, Eloc4has a maximum of ΔVd4. Hence, at half maximum we can write


For a very small dand R> > d, we can neglect the terms containing d2. Thus, the SR of gap-mode TERS becomes


A more detailed theoretical study of the TERS enhancement and spatial resolution in gap-mode TERS as a function of tip diameter is shown in Figure 2fg. The tip-substrate distance dwas kept fixed at 1 nm in the simulation. Both the calculated TERS EF and SR under 638 nm excitation show two distinctive regimes, one below and the other above a tip size of 28 nm. The sharp increase of the TERS EF and SR below 28 nm of the tip diameter is due to the lightning rod effect. On contrary, the slow increment of both EF and SR above 28 nm is the result of the LSPR of the tip apex. The red asterisks in the SR graph (Figure 2g) represent the results of Eq. (5) in very good agreement to the SR calculated from the LSPR contribution of the tip apex. Eq. (5) clearly indicates that for a given tip radius both EF and SR can be amplified by reducing the distance d. However, at very small d, the quantum mechanical effects such as non-local screening effects and electron tunneling weaken the TERS sensitivity and impose an unavoidable limiting factor avoiding a singularity at d = 0 [43, 44].


3. Surface selection rules of TERS in 2D semiconductors

According to the FEM simulations, the overall electric field intensity beneath the tip apex is zero for an incident beam polarized along the XY plane as shown in Figure 3a. Therefore, for an incident beam polarized along the z-direction (along the tip long axis) the intensity enhancement matrix can be written as [26].

Figure 3.

Incident light polarization dependentElocdistribution around the tip apex (a). Schematic of the gap-mode TERS in the side-illumination geometry (b). The dashed line projects the normal to the basal plane of 2D semiconductor at any given geometrical situation, which is inclined at an angleθto thez-axis.


The above enhancement condition is also applicable for the gap-mode TERS configuration. Hence, for a side-illumination geometry at an incident angle, which satisfies the tip polarization condition, a tip in gap-mode TERS configuration confines and enhances the Elocalong z-direction as shown in Figure 3b. Thus, the TERS intensity can be expressed as


Gωinand Gωscare the intensity enhancement matrices for incident and scattered light. αzzis the Raman tensor involved in a vibrational mode along the z-direction, which can be derived from the 3 x 3 matrix of the anisotropic polarizability tensor associated to the phonon vibrations as follows


The idea of surface selection rules was first conceptualized by Moskovits on flat metallic surfaces [45]. Due to the adsorption of molecules on a flat metal surface, they will experience a different electric field, since incident light is reflected off the metal substrate. Assuming both s-and p-polarized incident light causing Raman scattering, the scattered radiation at any given angle of incidence, ϕconsists of superposition of two waves as given by [45]


The primes are associated to the scattered radiation, rsand rpare the respective Fresnel reflection coefficient, and piis the polarization inside the molecule expressed by p=αE, where αcan be derived from Eq. (8). The four components associated to the Raman experiment are ss, sp, ps, and pp, the first and second letters refer to the incident and scattered light polarization, respectively. In the limit of a very good reflector, all the components except ppreduce to zero, since rs1and rp1. Therefore, the effective Raman scattering intensity becomes


According to Eq. (10) the effective Raman scattering depends on the z-component of the Raman tensor, normally not measured in ordinary Raman spectroscopy. However, in the gap-mode TERS configuration the coupling between the phonon with tensor elements having z-component and the Elocchanges the TERS intensity dramatically leading to so-called TERS selection rules.

In a real situation, the normal to the 2D material basal plane is parallel to the direction of Eloc(i.e. θ = 0 in Figure 3b) in a Raman backscattering geometry. Therefore, one can expect that the TERS signal is only associated to the αzzcomponent. However, if the normal to the 2D basal plane is tilted by an angle θdue to local structural deformation, other components of the Raman tensor also contribute to the overall TERS response.

It is also worth noting that above mentioned assumption is a simplified model, which does not consider interface chemistry. For example, in the case of chemisorption the molecular geometry or orientation may alter significantly leading to the breaking of the Raman selection rules. In such cases, some Raman inactive modes can become Raman active [46]. Moreover, light-plasmon coupling in a nanocavity between the tip and the metal substrate can lead to the alteration of classical Raman selection due to photon tunneling through perturbation of the evanescent field [47].

In polarization dependent Raman measurements we observe the scattered light intensity as a function of polarization directions of both the incident and scattered light. The Raman intensity can be written as


Here, esand eiare the unit vectors of the electric field of the scattered and incident light and Ris a second rank tensor also called Raman tensor derived from Eq. (8). There is plenty of literature including text books available discussing polarization and Raman tensors in conventional Raman measurements [25, 50, 51, 52]. In this chapter, we revisit the Raman tensors of first order Raman modes in 2D semiconductors, which are important to understand the TERS selection rules. Group theory analysis tells us that the Raman active phonon modes belong to the irreducible representations, the basis function of which contain quadratic terms of x, y, z. Thus, to have non-zero Raman intensity at a given esand eigeometry, Rshould have a non-zero matrix element. The Raman tensors of three first order Raman modes in 2D semiconductors (D3hpoint group) are presented below [25, 50].


There is one more doubly degenerate Raman mode with E(E1g) symmetry available in 2D semiconductors, which is forbidden in backscattering geometry. Therefore, the discussion of this mode is omitted here. Interested readers can learn more about this phonon feature in this review [50]. From Eq. (12) it is clear that the out-of-plane, A1(A1g) phonon mode has a Raman tensor component along z-direction (αzz). Therefore, the A1Raman mode is strongly sensitive to the TERS geometry. Since other modes are active in cross polarization geometry, they can still be observed in side-illumination TERS configuration, even though the TERS intensity is not as enhanced as for the out-of-plane Raman mode. An experimental demonstration of TERS selection rules on MoS2 is shown in Figure 4 [48, 49]. In the case of few layer MoS2 on an Au surface and resonantly excited by a 638 nm laser the far-field contribution on the substrate close to the flake shows comparable intensity of the in-plane E2gand out-of-plane A1gmode (Figure 4a). However, as the tip approaches the flake, the out-of-plane A1gmode enhances significantly [48]. Similarly, a marked TERS enhancement of the A1gmode is observed from the monolayer MoS2 deposited on Au nanodisc arrays on Si substrate and excited by 785 nm excitation (Figure 4b). Due to the below bandgap excitation (bandgap of monolayer MoS2 is ∼660 nm) only the in-plane E2gmode is observed in micro-Raman spectra. In contrast, TERS spectra of monolayer MoS2 on Au nanodisc show a significant modification with a giant enhancement of the A1gmode. Together with other factors (as discussed later in this chapter), the TERS selection rules amplify multifold the out-of-plane phonon modes [49].

Figure 4.

TERS spectra of few layer MoS2 deposited on an Au coated glass slide acquired using 638 nm excitation (a). Micro-Raman and TERS spectra of monolayer MoS2 deposited on an Au nanodisc array on a Si substrate acquired using 785 nm excitation (b). The micro-Raman spectrum was obtained over an area of ∼1 μm2 (spot size) and the TERS spectrum was taken locally on top of Au nanodisc with a spatial resolution of 2.3 nm. In both cases, together with theA1gRaman mode, the infra-red activeA2uis enhanced significantly. Figure (a) is adapted from ref. [48] (Copyright © 2017, IEEE) and (b) is adapted from ref. [49] with permission from the Royal Society of Chemistry.

Another important aspect of TERS selection rules is demonstrated by the enhancement of the infra-red (IR) active out-of-plane A2umode in both works. As shown in Figures 2 and 3, the tip-substrate dimer system creates a strong dipole moment in the nanocavity (or even sub-nanocavity) resulting in a strong out-of-plane Eloc. Therefore, IR active phonon modes, which induce a transition dipole moment parallel to the Eloc, can also be strongly enhanced. Since the IR active A2uphonon vibration creates a transition dipole moment perpendicular to the MoS2 basal plane, it couples to the Elocat the tip-substrate nanocavity polarized along z-direction. This is remarkable since the A2uphonon mode is Raman forbidden and only observed in resonant Raman conditions. Therefore, Voronine et al. [48] could still see this mode on the substrate as a far field contribution. However, the strong plasmonic field created at the tip-nanodisc nanocavity is able to not only excite this Raman forbidden IR mode, but also amplify it significantly under non-resonant 785 nm excitation as shown in Figure 4b [49].


4. Local phonon properties in 2D semiconductors

4.1 Strain induced local bandgap modulation

Strain plays an important role in manipulating optoelectronic properties of 2D semiconductors. The remarkable strength of 2D semiconductors such as TMDCs compared to other conventional semiconductors makes them well suited for straintronic applications. For example, a MoS2 monolayer can sustain a large biaxial strain >10%; whereas bulk silicon usually breaks at a strain of ∼1.5% [53]. This extraordinary property of 2D semiconductors has led to burgeoning research of strain engineered optoelectronic properties in recent times [53, 54, 55]. Among them, most of the studies provide macro/microscale information of strain induced optoelectronic manipulation using conventional optical spectroscopy. In contrast, the local modification of band structures due to spatially inhomogeneous strain distribution on the nanoscale is much less explored. Like the band structure, phonons in 2D semiconductors are directly affected by a wide range of parameters including strain [56, 57]. Therefore, Raman spectroscopy, one of the most frequently used non-destructive spectroscopies, becomes a powerful technique to investigate the different components of strain and their effect on optoelectronic properties of these semiconductors. Due to energy and momentum conservation, phonons participating in first order Raman scattering processes need to satisfy the wavevector condition q0(i.e.near the Γpoint in the crystal). Though, phonons with non-zero qcan also contribute to the Raman signals. However, they require two or more scattering processes assisted by a real electronic transition, a process called resonant-Raman scattering [58]. Hence, monitoring both the first order and the higher order Raman modes can provide quantitative knowledge of strain and its impact on the electronic structure of 2D semiconductors. Additionally, both tip enhanced photoluminescence (TEPL) and TERS can be performed simultaneously and can be used as complementary techniques to investigate local strain effects on the excitonic response. Though, special care needs to be taken to interpret the results since PL intensity, lineshape, and spectral position also depend on other local heterogeneities [24].

The two well-studied first order Raman modes in TMDCs are in-plane E2gand out-of-plane A1g. Using the combination of both experiments and ab initiocalculations shows that the in-plane E2gmode is more sensitive to uni- or bi-axial strain than the out-of-plane A1gmode [56, 60]. TERS reveals similar behavior for multilayer and monolayer TMDCs when probing highly localized strain at the nanoscale. Rahaman et al. studied local strain heterogeneities in 3 L MoS2 caused by a nanostructured gold substrate using TERS with a spatial resolution of less than 25 nm [59]. TERS measurements were performed in a side-illumination geometry under 638 nm excitation. When deposited on hexagonal periodic arrays of gold nanotriangles, 3 L MoS2 underwent local deformations on top and in between the nanotriangles as shown in Figure 5a. Hence, this is an ideal system to study the local flexibility at the nanoscale. Looking at the AFM topography and the 2nd order derivative of it, which represents the local curvature image (Figure 5b), it is clear that a variety of deformations caused strain inhomogeneities in MoS2. The center of the valley is unstrained and the twisted areas between two gold nanotriangles is most strained. The corresponding TERS map together with selective TERS spectra of the MoS2/Au heterostructure are presented in Figure 5c and d. As explained in the previous section the out-of-plane A1gmode is found to be more enhanced than the in-plane E2gmode due to TERS selection rules. After deconvoluting all the spectra using voigt functions, the frequency shift of the in-plane E2gmode agrees well with the local curvature image and a maximum strain of 1.4% is determined in the twisted regions from the frequency shift of 4.2 cm−1. However, at the corner of the nanotriangles the strain is calculated to be ∼0.9% from a frequency shift of 2.6 cm−1 of the E2gmode. The same authors then studied monolayer MoS2 on top of gold nanodiscs and observed a similar amount of strain at the rim of the discs [49].

Figure 5.

AFM topography (a) and 2nd order derivative of the topography (b). Both topography and 2nd derivative images show heterogeneous biaxial strain distribution in 3 L MoS2. The corresponding TERS map created for the spectral range of 350–490 cm−1 (c) and two representative TERS spectra (d) taken from two respective areas shown in (c). Site dependent peak position of two first order Raman modes for the quantification of strain (e). Adapted from ref. [59] (Copyright © 2017 American Chemical Society).

It is important to note that the above described TERS experiments were performed in resonant Raman condition with 638 nm excitation (excitation close to the Aoptical transition in MoS2) [61]. Therefore, there are as many as seven Raman features observed in the TERS spectra together with the two first order Raman modes (E2gand A1g). The Raman band around 378 cm−1 is the in-plane E2gM[62]. However, some reports argued that Davydov splitting of the in-plane E2gΓshould appear around this frequency position as well [13]. The feature around 390 cm−1 is a Raman inactive E1uexcited due to the resonant Raman condition in the TERS configuration [63]. The Raman mode around 410 cm−1 is an out-of-plane A1gphonon in the vicinity of the Mpoint of the BZ and analogous to the mode E2gM. The frequency around 420 cm−1 is attributed to a two phonon process involving a successive emission of a dispersive longitudinal quasi-acoustic (QA) phonon and a dispersionless transverse optical (TO) phonon along the c-axis [64]. The broad feature around 450 cm−1 is a convolution of minimum three modes. The first one around 440 cm−1 is E2gΓ+A1gΓ[62]. The frequency centered around 450 cm−1 involves two phonons of the longitudinal acoustic branch with opposite momenta at the Mpoint of BZ (2LAM) [65], and the third feature around 460 cm−1 is the infrared active A2umode as discussed in the previous section.

The rich variety of second order Raman features in MoS2 in resonant Raman condition also yields important information about its physical properties and the electronic band structure. Zhang et al. performed similar TERS experiments on monolayer MoS2 deposited on gold nanopyramids and by combining TERS and TEPL they probed the nanoscale variation of the electronic band structure induced by strain [60]. Under resonance excitation, a few of the second order Raman features in MoS2 were assigned to longitudinal acoustic LAand/or transverse acoustic TAphonons. These Raman bands scatter light through intervalley resonant scattering of the excited electrons by phonons when excited near the Bexciton energy. The process is known as double resonance Raman scattering (DRRS) and is similar to the process that leads to the creation of the Dand 2Dband in graphene/graphite [66, 67]. The DRRS process is explained by the schematic shown in Figure 6a. After the excitation of an e-h pair near the Kvalley of BZ, the excited electron is inelastically scattered from Kto the Kvalley by the emission of a phonon. In the relaxation process the excited electron is then inelastically scattered back to the Kvalley by the emission of a second phonon, where the e-h pair recombines and emits a photon. In MoS2, the Raman mode around 420 cm−1 and the 2 LAmode around 450 cm−1 are the results of DRRS processes. Zhang et al. observed that together with the shift in peak position of these second order features, the peak intensities are also inhomogeneously modified due to local strain [60]. Their DFT calculations of the phonon dispersion for acoustic branches show that DRRS processes are sensitive to the changes in the momentum and energy conservation constraints that govern which transitions and states in the BZ participate in the resonant interactions. Thus any distortion in the band structure caused by strain is reflected in the DRRS processes. With positive hydrostatic strain, the calculated phonon dispersion indicates a negative shift in peak position for both LAand TAmodes at the Kpoint in the BZ due to increasing lattice constants. On the other hand, the increase in peak areas associated to the DRRS processes depends on the relative position of the conduction band minima at the Kand Qpoint in the BZ. In an unstrained monolayer MoS2, the Kpoint conduction band minima is lower in energy than the Qpoint conduction band minima. Therefore, the KKtransition mediated by the qKacoustic phonons are much more prominent than the KQtransition mediated by qMphonons. However, when tensile strain is applied, the difference in energy between the Kand Qpoint conduction band minima decreases. Hence, more transitions mediated by the qMphonons satisfy the constraints on a DRRS process, which in turn increases the respective peak area.

Figure 6.

Illustration of the DRRS process in 2D semiconductors (a). AFM topography of monolayer MoS2 deposited on gold NTs (b). A line profile is created from the topography along the dashed line and presented in (c). Integrated Raman peak area of the DRRS triggered 2LAmode in MoS2 (d) taken at the arrow positions along the line profile shown in (c). Adapted from ref. [60] (© 2018 American Physical Society).

Like externally induced strain via transferring TMDCs monolayers on nanostructured substrates as discussed above, direct evaporation of a metal, for example gold on monolayer MoS2, can cause large strain (∼ 5%) at the interface due to lattice mismatch. Combining various SPM techniques and TERS Jo et al. demonstrated that large tensile strain induced at the interface in this process reduces the bandgap of MoS2 and thus influences the contact resistance as a consequence [68].

TERS has become a powerful technique to probe local strain variation in 2D semiconductors, which in turn helps us to understand the inhomogeneities in optoelectronic response of these materials.

4.2 Effect of variance in atomic registry on electronic properties

The reduced dimensionality of 2D semiconductors makes them susceptible to structural variations near the internal and external boundaries. Parameters, such as impurities, defects, and grain boundaries interfere strongly with intrinsic properties resulting in highly inhomogeneous optoelectronic response [23, 69, 70]. While conventional optical spectroscopy leaves a confusing picture regarding the effect of these heterogeneities [71], near-field optical studies can help us to understand the correlation of local structural heterogeneities with nano-optical response. In recent years, the combination of TEPL and TERS has been successfully introduced to investigate these heterogeneities. In TMDCs, the out-of-plane A1gRaman mode is particularly sensitive to changes in carrier concentration (or doping) via strong electron–phonon coupling [72]. Therefore, by comparing the peak positions of both in-plane E2gand out-of-plane A1gmodes it is possible to separate the local doping effect from other heterogeneities [49]. Additionally, the Raman intensity of the out-of-plane A1gmode is also affected by the physical and chemical irregularities such as doping, atomic vacancies, bond strain, grain boundaries (GB) at the atomic level [73]. Understanding the variations of atomic registries are especially important for 2D semiconductors as Park et al. showed that the excitonic behavior differs at the edges, nucleation sites (NS), and twin boundaries (TB) [24]. Using the combination of TERS, TEPL, and atomic force local strain microscopy they studied the PL response of PVD grown WSe2 and the results suggest that the defect concentration and stoichiometry of W and Se atoms in different sites can result in different excitonic mechanisms.

Due to their extraordinary physical properties, many groups are now synthesizing 2D semiconductors by CVD and other deposition techniques with the vision of industry scale production. However, in terms of crystal quality these deposition processes still require further optimization to be able to use them in electronic applications. CVD grown 2D semiconductors often reveal several structural heterogeneities as discussed in this section. Therefore, proper understanding is the key for the optimization of the deposition process and hence improving the crystal quality. While confocal micro-Raman spectroscopy is unable to obtain locally heterogeneous information, TERS can uncover them with nanometer precision. Smithe et al.studied the growth mechanism, GB, and defects in CVD grown monolayer MoSe2 using TERS [74]. Their investigation revealed that CVD MoSe2 monolayers have nano-domains of MoO3, which are invisible in confocal Raman spectroscopy. They also observed that the TERS intensity was strongly suppressed at GBs, which was attributed to the fast degradation of the GBs in ambient conditions due to the presence of higher defect concentration.

For optoelectronic devices the physics of semiconductor–metal interfaces is critical, since they are at the core of charge separation and collection. Jariwala et al.directly probed interfaces between WSe2 and noble metals such as gold and silver to investigate the influence of intrinsic doping on electronic properties using the cross-correlated spatial mapping of contact potential difference (CPD), photocurrent, TERS, and TEPL [75]. Their TERS map acquired with 638 nm excitation showed two distinct regions with domains of 10–100 nm lateral size, which correlated with CPD images as well. The regions of lower CPD exhibited non-resonant TERS spectra with a single Raman mode around 250 cm−1. This single Raman band corresponds to the combination of E2gand A1gmodes of WSe2. In contrast, TERS spectra representing higher CPD regions exhibited resonant Raman scattering with a number of additional features. It is important to note that the laser line 638 nm overlaps with the broad shoulder on the high-energy side of the Aexciton in WSe2. Therefore, Raman spectra acquired with this laser line induces resonant Raman response in this material with a number of additional peaks including DRR features, which are not observed with other excitation wavelengths. Later, by performing photocurrent measurement they observed that both higher and lower CPD areas generated higher photocurrents, but with opposite signs. Therefore, they concluded that the areas showing resonant Raman behavior are intrinsically n-type, while non-resonant areas are p-type.

4.3 Semiconducting to metallic phase transition

The stable crystal structure of TMDCs is the 2Hsemiconducting phase. However, these materials are known to form also metallic 1Tor 1Tphases (though not stable) due to high density of charge donation or electron doping [76, 77]. Due to strong light matter interaction TMDC monolayers, when integrated with plasmonic antennas, show great potential via modulation of their optoelectronic properties [78]. Generated hot electrons from plasmonic nanostructures can be injected into the conduction band of these 2D semiconductors due to favorable Schottky barrier heights at the interface and a large oscillator strength of these hot electrons [22, 79]. The concept of the doping induced structural phase change is illustrated in Figure 7a. According to the crystal field theory, the semiconducting 2Hphase consists of M (M = Mo/W) outermost dorbitals of a D3h- MX6 (X = S, Se) unit into three groups creating the stable 2Hlattice structure [80]. The doping via hot electrons first destabilizes the 2Hphase and induces a splitting of the metal outermost dorbitals to form metallic Oh-MX6 units into two groups as shown in Figure 7a [81, 82, 83]. The unoccupied orbitals then accommodate the extra electrons crossing the Schottky barrier (Figure 8).

Figure 7.

Single layer flakes of WSe2 exfoliated onto template stripped silver. Combined TERS map (intensity of 295 cm−1 peak in blue and of the 250 cm−1 (A1g+E2g) peak in red) (a). Photocurrent map at zero bias voltage collected simultaneously with TERS map (b), CPD image of the same area (c) and averaged TERS spectra showing resonant and non-resonant behavior (d). Section analysis of the photocurrent map again showing that the value of the photocurrent generated over differently doped domains was similar in absolute value, but of the opposite direction (e). The image is reproduced from ref. [75] (© 2018 IOP Publishing Ltd).

Figure 8.

High resolution TERS map of monolayer MoS2 on gold nanodiscs acquired with a gold tip and 785 nm excitation (a). Corresponding TERS spectra along the circles shown in the map (b). Thanks to the very high spatial resolution the structural phase change from 2Hto 1Tdue to highly localized plasmonic hot electron doping can be observed. The image is reproduced from ref. [49] with permission from the Royal Society of Chemistry.

Milekhin et al.studied the extremely localized site dependent transient reversible 2Hto 1Tphase change in monolayer MoS2 using TERS. Monolayer MoS2 was deposited on periodic arrays of gold nanodiscs, thus sandwiched between the nanodisc and the tip to conduct this experiment [49]. Hence, the tip-nanodisc system creates a strongly coupled dimer, which then pumps a high concentration of electrons into the conduction band of MoS2. As demonstrated by the high resolution TERS images and later by FEM simulations by Rahaman et al.[37] the coupling is the strongest between the tip and the rim of the nanodiscs. Therefore, the rim of the nanodiscs became the strongest source of hot electron doping in MoS2, which consequently showed the transition from the 2Hto 1Tphase via the appearance of an additional Raman mode in TERS spectra. Metallic MoS2 has three additional Raman modes around 156, 226, and 333 cm−1 known as J1, J2, and J3, respectively [84]. Due to strong increasing Rayleigh background, Milekhin et al.[49] were able to observe only the J3 Raman mode around 333 cm−1. As can be seen in Figure 7b, the transition from 2Hto 1Tis extremely site dependent with an area limited by the spatial resolution, which was 2 nm.

4.4 Probing edge related properties

TMDCs have two distinct edges in the form of armchair and zigzag, which possess different electronic properties. For example, the zigzag edge in MoS2 is metallic, whereas the armchair is semiconducting [85]. Therefore, quantitative measurements or characterization of these edges are important for effective edge engineering, especially for edge contacts for electronic device applications. Huang et al. studied edge related properties of atomically thin MoS2 using TERS [69]. Their local Raman measurements aided by strong plasmonic enhancement of either gold or silver tips probed different Raman features of the edge defects in mono- and bilayer MoS2 corresponding to the unique lattice vibration and electronic properties of the respective edge.

A broad Raman feature around 220 cm−1 was observed along the edges (full of dangling bonds) and absent on the wrinkles and main body (Figure 9a). The peak was assigned to the defect (dangling bonds) induced acoustic LAwith a DRRS feature as explained in Section 4.1. More interestingly, the TERS spectra taken along the edges of the bilayer and the step between mono- and bilayer exhibits an additional feature around 396 cm−1. This particular feature is explicitly observed in bulk at cryogenic temperature and assigned to the LAM+TAMacoustic mode [86]. This Raman mode involves a DRRS process between the Kand Qvalleys with wavevector qMand thus the Raman cross section of this mode depends on the relative energy position of the conduction band minima at the Kand Qpoints in the BZ. Since in bulk the conduction band minima at the Qpoint is lower than the Kpoint, the Raman scattering efficiency of this mode becomes larger. Therefore, this mode is absent in the basal plane of mono- and bilayer MoS2. However, along the edges, especially in the case of bilayer MoS2 the modification of the electronic band structure leads to the considerable lowering of the conduction band energy at the Qpoint resulting in the activation of the DRRS process of the LAM+TAMmode.

Figure 9.

TERS of mono- and bilayer MoS2 acquired along the edges and the step between mono- and bilayer (a). AFM topography of a monolayer MoS2 showing different terminated edges (b). Fitted peak position ofA1gmode acquired along a line crossing both zigzag and armchair edges (c). DFT calculated Raman spectra of at different sites in MoS2 (d). Reproduced from ref. [69].

Another important observation of this study is the edge dependent local Raman sensitivity. The frequency position of the out-of-plane A1gmode showed either upward or downward shift compared to the basal plane when probed along the armchair or zigzag edges, respectively. The AFM topography of monolayer MoS2 presented in Figure 9b shows edges terminated with different angles. For an angle of 60° both edges have the same structure (either armchair or zigzag) and for an angle of 90° the edges are different. The TERS spectra acquired along a line crossing both zigzag and armchair edges show a clear pattern, namely that the A1gmode downshifts at one of the edges and upshifts at other one (Figure 9c). According to the DFT calculations the local strain introduced by the two edge structures cause this opposite direction of the frequency shift as shown in Figure 9d.


5. Conclusion and outlook

TERS has developed into a versatile characterization technique for a wide range of materials. In recent years we have witnessed chemical mapping of single molecules with a spatial resolution of the bond length using TERS. For 2D semiconductors, an excellent SR of ∼2 nm in TERS helped us to understand local heterogeneous behavior of monolayer MoS2 and demonstrate its tremendous capabilities. In this chapter, we explained the basic principal of TERS and described the surface selection rules in regards to the 2D semiconductors important for understanding the near-field Raman spectra. We also reviewed the recent progress of this technique in the 2D semiconductor field. The potential of TERS certainly guarantees new breakthroughs in the 2D semiconductor field in the near future. One of the recent direction of 2D semiconductors is to create heterostacks like Lego, which promises exotic physics due to the creation of moiré superlattices. The size of a moiré unit cell varies from nm to few 10s nm depending on the lattice mismatch and rotation angle. The moiré superlattice forces the constituent monolayers into phonon renormalization, also known as moiré phonons. TERS can be used to probe these moiré phonons, thus resolving the moiré supercell critical for understanding the moiré physics. Besides, there is a lot of information yet to be resolved on how local heterogeneities control the optical response such as strain induced exciton funneling, defects induced excitonic response, or single photon emitters etc. in 2D semiconductors, for which TERS can be a powerful tool to achieve a thorough understanding.



This work was supported by the German research council funded project DFG-ZA 146/44-1.


Conflict of interest

The authors declare no conflict of interest.


chapter PDF

© 2021 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Mahfujur Rahaman and Dietrich R.T. Zahn (October 6th 2021). Tip-Enhanced Raman Spectroscopy of 2D Semiconductors [Online First], IntechOpen, DOI: 10.5772/intechopen.99817. Available from:

chapter statistics

1total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us