1. Introduction
In the past few years, there has been a growing interest in Transition Metal Dichalcogenides (TMDCs), specifically in the form of monolayer (ML) van der Waals materials. This has sparked significant research interest for various applications in electronics and optoelectronics [
1,
2,
3].
In comparison to conventional semiconductors, such as GaAs, the Coulomb interaction between conduction electrons and valence holes, as well as the oscillator strengths of excitons in MLs TMDCs, are significantly higher due to the two-dimensional confinement of charge carriers, heavy effective masses, and weak screening in 2D systems [
4,
5,
6,
7,
8,
9]. That is why scientists studying the physics of semiconductor nanosystems have been very interested in 2D materials during the past few years, such as graphene, hexagonal boron nitride MLs, TMDC MLs, and the heterostructures they generate [
10,
11]. The ML TDMCs that have been investigated the most are MoS
2, MoSe
2, WS
2, and WSe
2 [
12,
13].
The transport of carriers in ML 2D materials under low fields is influenced by multiple scattering mechanisms, including interactions with acoustic and optical phonons. Additionally, scattering can occur due to polar coupling with the substrate beneath or with dielectrics, introducing another factor involving remote optical phonons. Thus, Polar optical phonons, situated at the interface, play a significant role in scattering carriers in TMDCs through the Fröhlich coupling. [
14,
15,
16,
17]. Therefore, it is crucial to comprehend these scattering events by examining the coupling between surface optical phonons (SOP) and electrons in TMDCs. Developing models that can elucidate experimental results becomes essential. This coupling is typically characterized by interactions between electronic excitations and phonons, giving rise to various intriguing effects in a crystal, including the formation of polarons [
14,
15,
16,
17,
18,
19,
20].
Several studies in the literature proved the importance of the role played by the SOP coupling in the optical proprieties in ML TDMCs deposited on polar substrates. As an example, we referenced Suvodeep Paul et al. [
21], who showed that the creation of a
heterostructure leads to the coupling between electrons in
and polar phonons in
. This coupling governs the enhancement of defect-bound excitons and biexcitons. Additionally, they have performed an extensive resonant Raman analysis, varying both polarization and magnetic field, which provided further confirmation of the electron-phonon coupling in the
heterostructure.
Likewise, Colin M. Chow et al. [
22], observed a resonant Raman scattering phenomenon through cross-material exciton-phonon coupling at van der Waals interfaces. They noted that the sensitivity of excitons in monolayer materials to their phononic environments, such as those provided by
,
, and sapphire, can be exploited to deepen our understanding of atomically thin devices. Else, Bastian Miller et al. [
23], investigated exciton-phonon coupling in charge-tunable single-layer MoS
2 devices using polarization-resolved Raman spectroscopy. They found a strong defect-mediated coupling between the long-range oscillating electric field of the longitudinal optical phonon in the dipolar medium and the exciton.
Sanjay Gopalan et al. [
24] also explored the impact of the dielectric environment on electronic transport in monolayers TMDCs. By employing ab initio methods, they calculated the low-field carrier mobility in free-standing layers, considering the effects of dielectric screening on electron-phonon interactions induced by the bottom oxide and gate insulator, as well as scattering from hybrid interface optical-phonon/plasmon excitations. Their findings revealed that using insulators with a high dielectric constant can greatly improve carrier mobility.
This paper is structured as follows: first, we theoretically investigate the interaction between electrons and surface optical phonons in ML TMDCs on polar substrates, such as silicon carbide () and hexagonal boron nitride (). Furthermore, we present a theoretical examination of the polaronic oscillator strength in ML TMDCs on polar substrates. Finally, we investigate the temperature dependence of the polaronic scattering rate in ML TDMCs on and h polar substrates.
2. Electron-Surface Optical Phonon Interaction in ML TMDCs on and h Dielectric Substrates
A monolayer transition metal dichalcogenide (TMDC) consists of a central layer of metal M atoms arranged in a triangular lattice, flanked by two layers of chalcogen X atoms positioned on the same triangular lattice. The triangular Bravais lattice is defined by the basis vectors:
and
(see
Figure 1a).
Figure 1b illustrates the reciprocal lattice, defined in relation to the triangular Bravais lattice and characterized by the vectors:
and where is the lattice constant.
The two-dimensional Brillouin zone of the TMDCs exhibits a hexagonal shape, featuring high-symmetry points denoted as each defined as follows:
, ,
The effective
2×2 Hamiltonian characterizing the states of the conduction and valence bands with the parallel spins
in the vicinity of the point
is represented by the following expression [
25,
26,
27]:
(1)
Where, denotes the two-dimensional wave vector of the electron measured from the point , the parameter is directly proportional to the interband matrix element of the momentum operator, where is the electron effective mass and represents the width of the band gap.
The Hamiltonian describing a pair of spin sublevels with in the same valley has the form of Equation (1) through the substitution , where represents the sum of the spin–orbit splittings of the conduction and valence bands. The effective Hamiltonian in the , valley is derived from Equation (1) through the substitution .
The energy spectrum of the electron derived from the Hamiltonian in Equation (1) has the Dirac form:
(2)
Here, and correspond to the conduction and valence bands, respectively.
In our study, we use the assumption of homogeneous and defect-free interfaces between transition metal dichalcogenide monolayers and dielectric substrates as a simplification often used in theoretical models and simulations to make the problem tractable.
In this work we have investigated the electron surface optical phonon (SOP) interaction in ML TDMCs on
and
dielectric substrates across the long-range Fröhlich coupling. Indeed, the long-range Fröhlich coupling model provides a robust framework for understanding electron-SOP interactions in TMDCs on polar substrates. However, it relies on several approximations; like for example Born-Oppenheimer Approximation [
28]. Short-range interactions [
29] (e.g., electron-phonon interactions in non-polar materials) are not considered. It assumes a constant effective mass for the electron, non-linear interactions and multi-phonon processes are typically neglected [
30]. Impurities, defects, and other forms of disorder that can affect the electron-phonon interaction in real materials are usually not included in the idealized Fröhlich model [
28,
31]. The phonon dispersion is typically assumed to be linear, which is an approximation that might not hold for all phonon modes or substrates [
32]. Often, a single dominant phonon mode is considered, neglecting the possible contribution of multiple phonon modes [
29].
To simplify, our analysis, we consider the phonon spectrum as isotropic, implying that phonons exhibit either longitudinal or transverse polarization. The Fröhlich Hamiltonian introduces an interaction term wherein an electron scatters from to involving the emission or absorption of a phonon. In both cases the total momentum is conserved and is expressed as follows:
(3)
The term denotes the phonon energies, incorporating both the Longitudinal Optical (LO) and Surface Optical (SO) modes, and can be expressed as:
(4)
In this context, , represent the creation and annihilation operators, respectively, for the phonon characterized by the wave vector , while refers to the frequency of the phonon.
The second term
is the Hamiltonian of interaction between electron and phonon [
33]:
(5)
The Fröhlich Hamiltonian is given as follow:
(6)
The interaction between carriers in monolayer Transition Metal Dichalcogenides (TMDCs) and surface optical phonons is described by the second term in Equation (6).
The coupling element in the Fröhlich Hamiltonian
represents the interaction between the electron in TMDCs and surface optical phonon of the polar substrates. This matrix element is expressed as [
34,
35,
36]:
(7)
In the given context,
represents the magnitude of the polarization field, which is determined by the Fröhlich coupling [
37]:
(8)
While
and
are the low- and high-frequency dielectric constants of the polar substrate, (refer to
Table 1),
represents the internal distance between the TMDCs and polar substrate (refer to
Table 2). The term
denotes the energy of SO phonon of the polar substrates with two branches
The SOP energies are extracted from the bulk longitudinal optical (LO) phonons as follows [
38]:
(9
The screening of the Coulomb interaction by the polar dielectric environment is considered through
. Given the weak screening of the electric field perpendicular to the plane of the ML TMDCs,
is set to1 [
44].
On polar substrates, surface optical phonons (SOP) induce an electric field that interacts with the electrons in the neighboring ML TMDCs. Using Equations (7) and (8), the SOP coupling is expressed as:
(10)
The summation is performed over one spin and one valley, where is the area of the two-atom unit cell.
In our analysis in the present case, we have followed the same theoretical method presented in our previous calculations [
14,
15,
16]. So, to study the interactions between electron and surface optical phonons in the ML TMDCs, we have specifically considered the electronic states
and
, with electron energies
and
, respectively. We have also considered the effective 2×2 Hamiltonian characterizing the states of the conduction and valence bands with the parallel spins
in the vicinity of the point
of the hexagonal Brillouin zone.
The space of polaronic states results from a tensor product between the two subspaces of electronic and phononic states. Thus, we consider new states called polaronic states given by:
(11)
The Polaron electron energies
for the states
in ML TMDCs on polar substrates are given bellow [
14,
15,
16]:
(12)
Figure 2 depicts the SO coupling strength between the electronic states
and
versus the wave vector
in ML WS
2 on
and
h polar substrate. As shown in
Figure 2, it is evident that the coupling with SOP is significantly influenced by the type of polar substrate.
We show in
Figure 3, the Polaron electron energies versus the wave vector
in ML
on
and
h polar substrate. For comparison purposes, we have included in the same figures the energies of the non-interacting states
and
. For example, in the case of WS
2 these noninteracting levels cross near
in the case of
and near
in the case of
indicating resonant coupling (see
Table 3). These crossings indicate that the separation between electronic levels is equal to
for both
and
cases, where
respectively. In fact, the electronic level crossings are clearly replaced by significant anticrossings with energy levels approximately at ~94 meV and ~70 meV for
and
polar substrates, respectively. It can be also observed in
Figure 3, the enhancement of the Rabi splitting of the electron levels when shifting from
to
(refer
Table 4).
In these anticrossings, the wave functions of the levels become mixed, allowing for multiple transitions such as , and . This demonstrates that the interaction between electrons and surface polar phonons cannot be considered a weak coupling. The coupling between electrons and SOP leads to the Rabi splitting of the electron levels. Hence, the calculations indicate the possibility of energetically resonant coupling between the electronic sub-levels and surface vibration modes in ML TMDCs on the studied polar substrates. Furthermore, the two resulting polaron states can be expressed as:
(13)
The weight of electronic component
and the weight of the one-phonon component
of the polaron states ± vary with the polaron energies
. The expressions detailing these dependencies are as follows [
14,
15,
16]:
(14)
(15)
where
is the SO coupling strength between the electronic states
[
14,
15,
16,
17].
Figure 4, depicts the weight of the electronic components and the one-phonon components of the lower polaron state
in ML TMDCs for
and
h polar substrates as a function of the wave vector
in ML TMDCs. It can be seen for example in the case of
on
polar substrate (see
Table 3,
Figure 4) that; when the wave vector
approaches
the value of the weight of the one-phonon component for the lower polaron state
is much larger compared to that of the electronic component
. This result demonstrates that the SOP situated at the interface of the ML
on
polar substrate plays a crucial role in the resonant coupling between the noninteracting states
and
allowing the formation of the polaron states. It can be noted that the same result has been proved for the other cases of ML
on
and
h polar substrates (see
Table 3 and
Figure 4).
3. Polaronic Oscillator Strength of ML TMDCs on SiC and hBN Polar Substrates
In the following section, we theoretically investigate the polaronic oscillator strength (OS) which is another crucial quantity. Drawing an analogy to the oscillator strengths of interband transitions in quantum dots, we have computed the OS for interband transitions in ML TMDCs on polar substrates. In the strong confinement limit, the OS is linked to the overlap integral of the polaronic states,
by [
45,
46]:
(16)
where is the Kane energy and is the emission energy for one-phonon of the ML TMDCs on polar substrates, which is given by:
(17)
where is the lower polaron energy of the exciton, is the emitted photon energy and is the energy gaps of the ML TMDCs.
We have calculated the OS for the lower polaron state which is a linear combination of the two states and :
(18)
Figure 5 shows the polaronic OS of
on
and
h polar substrates versus the wave vector
in ML
.
As a result, we proved theoretically; that the polaronic OS is especially sensitive to the phonon mode of the surrounding dielectrics. In fact this result is due to the emission energy for one-phonon of the ML TMDCs for
and
h polar substrates, which is given by:
. Hence, the highest polaronic oscillator strength corresponds to the highest optical phonon energy of the polar substrates
. This is by analogy to polaronic oscillator strength in ML TMDCs which is also much higher as compared with conventional semiconductors [
6,
7,
8,
9].
Using the same method, we can easily prove that for the other ML TMDCs such as: , , and :
. (19)
This result is due to the
dielectric constant as well as the
phonon energy compared to that of
, (see
Table 1), so the polar field created near the interface should be the highest in the case of
substrate, compared to
substrate. This result leads to the highest polaronic OS in ML TDMCs on
polar substrate. Similarly, it can be concluded from the
Figure 6 that:
. (20)
Hence, this result can be explained by the difference in electron effective masses in ML TDMCs. This result actually confirms the decrement in the polaronic OS for heavy electron in the nearby ML TDMCs / dielectric substrate interface. Otherwise, in the case of light electron, the polaronic OS increases considerably at the strong confinement regime.