Theory of interlayer exciton dynamics in 2D TMDCs Heterolayers under the influence of strain reconstruction and disorder

Marten Richter

DOI: https://doi.org/10.1103/PhysRevB.109.125308

2024-02-27

Abstract:Monolayers of transition metal dichalcogenides (TMDC) became one of the most studied nanostructures in the last decade. Combining two different TMDC monolayers results in a heterostructure whose properties can be individually tuned by the twist angle between the lattices of the two van-der-Waals layers and the relative placement of the layers, leading to Moiré cells. For small twist angles, lattice reconstruction leads to strong strain fields in the Moiré cells. In this paper, we combine an existing theory for lattice reconstruction with a quantum dynamic theory for interlayer excitons and their dynamics due to exciton phonon scattering using a polaron transformation. The exciton theory is formulated in real space instead of the commonly used quasi-momentum space to account for imperfections in the heterolayer breaking lattice translational symmetry. We can analyze the structure of the localized and delocalized exciton states and their exciton-phonon scattering rates for single phonon processes using Born-Markov approximation and multi-phonon processes using a polaron transformation. Furthermore, linear optical spectra and exciton relaxation Green functions are calculated and discussed. A P-stacked MoSe$_2$/WSe$_2$ heterolayer is used as an illustrative example. It shows excitons localized in the potential generated through the Moiré-pattern and strain and a delocalized continuum. The exciton-phonon relaxation times vary depending on the strain and range from sub-pico seconds up to nanoseconds.

Mesoscale and Nanoscale Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to understand the interlayer exciton dynamics in two - dimensional transition - metal dichalcogenide (TMDC) heterostructures, which are affected by strain reconstruction and disorder. Specifically, the authors combined the existing strain - reconstruction theory and quantum - dynamics theory to study interlayer excitons and their dynamic behaviors under exciton - phonon scattering. By introducing the polaron transformation, the authors were able to describe these phenomena in real space, thus considering the imperfections in the heterostructures (such as the breaking of lattice translational symmetry) more accurately. In addition, the authors also analyzed the structures of localized and non - localized exciton states and their exciton - phonon scattering rates, and calculated and discussed the linear optical spectra and exciton - relaxation Green's functions.

### Main problems:
1. **Interlayer exciton dynamics**: How to describe the dynamic behaviors of interlayer excitons in the presence of strain reconstruction and disorder?
2. **Exciton - phonon interaction**: How to quantitatively describe the interaction between excitons and phonons, especially the influence of multi - phonon processes?
3. **Optical properties**: How to predict and explain the linear optical spectra and exciton - relaxation characteristics of these systems through the theoretical framework?

### Research methods:
- **Strain - reconstruction theory**: Calculate the strain field through energy minimization, considering supercells of different scales and directions.
- **Quantum - dynamics theory**: Use the Schrödinger equation to describe the interlayer exciton wave function, and separate the relative motion and the center - of - mass motion through the Wannier equation.
- **Polaron transformation**: Introduce the polaron transformation to deal with multi - phonon processes and modify the Fermi golden rule to include multi - phonon scattering.

### Key formulas:
1. **Total energy expression**:
\[
E=\int_{S}[U(\mathbf{r}) + W_{s}(\mathbf{r})]d^{2}\mathbf{r}
\]
where \(U(\mathbf{r})\) is the intralayer contribution and \(W_{s}(\mathbf{r})\) is the interlayer contribution.

2. **Intralayer strain energy**:
\[
U(\mathbf{r})=\sum_{l = t,b}\left[\frac{\lambda_{l}}{2}(u_{ii}^{(l)})^{2}+\mu_{l}(u_{ij}^{(l)})^{2}\right]
\]

3. **Interlayer strain energy**:
\[
W_{s}(\mathbf{r})=V_{s}(\mathbf{r})-\epsilon Z_{s}^{2}(\mathbf{r})
\]

4. **Schrödinger equation for the interlayer exciton wave function**:
\[
E\Psi(\mathbf{r}_{e},\mathbf{r}_{h})=\left[-\frac{\hbar^{2}\nabla_{e}^{2}}{2m_{e}}-\frac{\hbar^{2}\nabla_{h}^{2}}{2m_{h}}+V_{eh}(\mathbf{r}_{e}-\mathbf{r}_{h})+V_{h}(\mathbf{r}_{h})+V_{e}(\mathbf{r}_{e})\right]\Psi(\mathbf{r}_{e},\mathbf{r}_{h})
\]

5. **Coulomb interaction potential**:
\[
V_{eh}(r)=\frac{e^{2}}{4\pi\epsilon_{0}\epsilon_{r}r_{\text{eff}}}\left[\log\left(\frac{r}{r + r_{\text{eff}}}\right)+(\gamma-\log(2))e^{-r/r_{\text{eff}}}\right]
\]
where \(r_{\text{eff}}=r^{*}+r'^{*}+d\), \(r^{*}\) and \(r'^{*}\) are material parameters, and \(\gamma\) is the Euler constant.

6. **Single - particle potential**