Phonon-induced Markovian and non-Markovian effects on absorption spectra of moiré excitons in twisted transition metal dichalcogenide bilayers

This theoretical study reveals that phonon-induced effects on the absorption spectra of moiré excitons in twisted transition metal dichalcogenide bilayers transition from non-Markovian dynamics and phonon sidebands at small twist angles to Markovian broadening at larger angles, while intraband scattering significantly suppresses higher-energy absorption bands when their bandwidth exceeds the optical phonon energy.

The Gist

Imagine you have a piece of fabric made of two ultra-thin, transparent sheets of a special material (like a high-tech silk). If you stack them perfectly on top of each other, they look like one sheet. But, if you twist one sheet slightly relative to the other, a new pattern emerges on the surface: a giant, honeycomb-like grid called a Moiré pattern.

In the world of physics, this twisting creates a "trap" for tiny particles of light and electricity called excitons. These excitons are like little couples dancing on the fabric. The way they dance depends entirely on how much you twist the sheets.

This paper is a deep dive into how vibrations in the fabric (called phonons) affect these dancing couples, and how the "twist angle" changes the rules of the dance.

Here is the breakdown using simple analogies:

1. The Two Dance Floors: Tight vs. Loose

The researchers found that the twist angle creates two very different environments for the excitons:

• The "Tiny Room" (Small Twist Angle): When you twist the sheets just a tiny bit, the Moiré pattern creates huge, deep traps. The excitons get stuck in one spot, like a dancer trapped in a tiny, cozy room. They can't move around much.

  • The Vibe: Because they are stuck, they interact with the fabric's vibrations in a very complex, "remembering" way. If the floor vibrates, the dancer feels it, reacts, and the floor remembers that reaction for a while. This is called Non-Markovian behavior. It's like a conversation where everyone keeps interrupting and remembering what was said five minutes ago.
  • The Result: The light they absorb looks like a single, sharp peak with a fuzzy "halo" around it (sidebands).

• The "Open Ballroom" (Large Twist Angle): When you twist the sheets more, the traps get shallower and wider. The excitons are no longer stuck; they can run freely across the fabric.

  • The Vibe: Now, the interaction with the floor vibrations is quick and forgetful. The dancer bumps into a vibration, gets a little push, and moves on immediately without looking back. This is Markovian behavior. It's like a fast-paced conversation where you only care about what is being said right now.
  • The Result: The light absorption peak gets wider and looks lopsided (asymmetric), similar to how light behaves in a single sheet of this material.

2. The Two Types of Floor Vibrations

The fabric vibrates in two main ways, and the paper explains how each affects the dancers:

• The "Rumble" (Acoustic Phonons): These are low-frequency, long-wavelength vibrations, like a slow rumble of the floor.

  • In the Tiny Room, this rumble creates a complex, symmetrical blur around the dancer's light.
  • In the Open Ballroom, the rumble causes the dancer to stumble and lose energy quickly, making the light peak stretch out to one side.

• The "Clap" (Optical Phonons): These are high-frequency, short-wavelength vibrations, like a sudden, sharp clap of the floor.

  • These create distinct "echoes" or side-peaks in the light spectrum.
  • The Big Discovery: The paper found a "magic trick" with these claps. If the dancer is in a high-energy state (a higher band) and the "room" is wide enough (large twist angle), the dancer can use a single "clap" to jump down to a lower energy state instantly.
  • The Consequence: This jump is so efficient that it completely erases the light signal from that higher-energy dancer. It's as if the dancer suddenly becomes invisible because they are constantly jumping down to the floor before anyone can see them. This explains why some expected light peaks disappear in experiments.

3. The "Magic Angle"

Just like in the famous "Magic Angle" graphene discovery, this paper identifies specific twist angles where things get weird.

• There is a specific angle where the "rumble" of the floor perfectly matches the energy needed for the dancer to jump from one spot to another. At this angle, the interaction becomes super strong, and the dancer's behavior changes drastically.

Why Does This Matter?

Think of these materials as the future building blocks for quantum computers and super-efficient solar cells.

• To build a quantum computer, you need particles that stay in one place (the "Tiny Room") and don't get confused by noise.
• To build a solar cell, you want particles that can move freely (the "Open Ballroom") to carry energy.

This paper gives engineers a "recipe book." It tells them: "If you want sharp, stable light for quantum tech, twist the sheets by 1 degree. If you want to absorb light efficiently for solar power, twist them by 5 degrees. And watch out for the 'claps' (optical phonons), because they might make your higher-energy signals disappear!"

In a nutshell: By simply twisting two sheets of material, we can switch the physics from a complex, memory-filled dance to a simple, fast-paced one, and we can even make certain light signals vanish entirely. This gives us precise control over how these materials interact with light.

Technical Summary

1. Problem Statement

Transition metal dichalcogenide (TMDC) bilayers, when twisted relative to each other, form a moiré superlattice. This superlattice creates a periodic potential that localizes excitons, with the degree of localization depending on the twist angle ($\theta$).

• Small twist angles: Excitons are quasi-localized with flat energy dispersions (0D-like behavior).
• Large twist angles: Excitons become delocalized with curved dispersions (2D-like behavior).

While the influence of phonons on 0D emitters (like quantum dots) and 2D materials (monolayers) is well-studied separately, the twist-angle-dependent crossover between these regimes in moiré systems remains poorly understood. Specifically, the interplay between Markovian (memoryless, leading to line broadening) and non-Markovian (memory-dependent, leading to sidebands and dephasing) phonon effects on the optical absorption spectra of these tunable systems was not theoretically quantified.

2. Methodology

The authors developed a comprehensive theoretical framework to model the polarization dynamics and absorption spectra of intralayer moiré excitons in a twisted MoSe$_2$/WSe$_2$ heterobilayer.

• Hamiltonian Formulation:

  • Started with a two-band effective mass model for electrons and holes.
  • Introduced a moiré potential ($V_m$) derived from the twist angle, which couples homogeneous exciton states via reciprocal lattice vectors.
  • Derived moiré exciton operators ($X_{n,k}$) as Bloch-like states, solving an eigenvalue equation to obtain twist-angle-dependent band structures ($\hbar\omega_{n,k}$).
  • Modeled exciton-phonon coupling (both acoustic and optical) using a linear deformation potential coupling, transforming the interaction into the moiré exciton basis.

• Dynamics Simulation (TCL Master Equation):

  • Employed a Time-Convolutionless (TCL) master equation approach to describe the microscopic polarization dynamics ($p_n(t) = \langle X_{n,0} \rangle$).
  • This approach captures non-Markovian memory effects explicitly through time-dependent dissipation coefficients $\Gamma(t)$, which are derived from the Generalized Phonon Spectral Density (gPSD).
  • The gPSD, $\rho(\Omega)$, encodes the density of states of the moiré excitons and the phonon dispersion, determining the strength of energy-conserving (Markovian) and energy-mismatched (non-Markovian) transitions.

• Absorption Calculation:

  • Calculated the linear absorption spectrum $\alpha(\omega)$ by Fourier transforming the time-evolved polarization after an ultrashort optical pulse.
  • Separated the spectral features into a Zero-Phonon Line (ZPL) and Phonon Sidebands (PSBs) using a power series expansion of the dephasing function.

• Parameters:

  • Simulated MoSe$_2$ intralayer excitons in a MoSe$_2$/WSe$_2$ bilayer.
  • Considered twist angles $\theta = 1^\circ, 3^\circ, 5^\circ$ and temperatures $T = 4, 70, 200$ K.
  • Included both acoustic (linear dispersion) and optical (dispersionless) phonon branches.

3. Key Contributions and Results

A. Twist-Angle Dependent Regimes (Acoustic Phonons)

The study identifies two distinct physical regimes governed by the twist angle:

1. Small Twist Angles (Flat Bands, e.g., $1^\circ$):

   • The system behaves like a localized 0D emitter (similar to quantum dots).
   • Non-Markovian dynamics dominate: The polarization decay is slow, and the Markovian approximation significantly underestimates the decay.
   • Spectral Features: The absorption spectrum is dominated by a broad, symmetric phonon sideband (PSB) with a suppressed ZPL at elevated temperatures. This mirrors the "independent boson model" behavior.
   • Magic Angle: A specific "magic angle" ($\theta_c \approx 2.9^\circ$ for acoustic phonons) exists where the energy mismatch for transitions between the $\gamma$-point and the van Hove singularity ($m$-point) vanishes, maximizing the Markovian decay rate.

2. Large Twist Angles (Curved Bands, e.g., $5^\circ$):

   • The system behaves like a delocalized 2D exciton (similar to monolayer TMDCs).
   • Markovian dynamics dominate: The polarization decay is well-described by a constant decay rate, and memory effects vanish quickly ($\sim 3$ ps).
   • Spectral Features: The spectrum exhibits a characteristic asymmetric peak with a tail towards lower energies, a signature of delocalized excitons coupled to acoustic phonons.

B. Optical Phonon Effects

• Sidebands: Optical phonons introduce distinct sidebands at energy offsets $\pm n\hbar\Omega_{opt}$.
• Bandwidth Connection: The width of these sidebands is directly linked to the moiré exciton bandwidth.
  • At small twist angles (flat bands), sidebands are narrow peaks.
  • At large twist angles (wide bands), sidebands broaden and eventually merge with the ZPL at high temperatures.
• Markovian Limit: Unlike acoustic phonons, optical phonons generally do not contribute to the Markovian decay rate at $\Omega=0$ unless a "magic angle" condition is met where the bandwidth allows for energy-conserving transitions at the optical phonon energy.

C. Multi-Band System and Peak Suppression

When considering multiple bright moiré exciton bands:

• Dipole Moments: The oscillator strength (dipole moment) of higher bands decreases as the twist angle increases, due to the folding of the Brillouin zone and symmetry constraints ($C_{3v}$ point group).
• Phonon-Induced Suppression: A critical finding is that intraband scattering via optical phonons can completely suppress absorption peaks of higher-lying bands.
  • If the bandwidth of a higher band exceeds the optical phonon energy ($\hbar\Omega_{opt}$), efficient energy-conserving scattering occurs.
  • This leads to a large Markovian decay rate for that specific band, effectively washing out its absorption peak in the spectrum, even if it has a non-zero dipole moment.

4. Significance

• Unified Framework: The paper provides a unified theoretical description that bridges the gap between 0D localized excitons and 2D delocalized excitons, showing how the transition is controlled by the twist angle.
• Predictive Power: It predicts specific spectral signatures (symmetric vs. asymmetric peaks, sideband widths) that can be used to experimentally determine the twist angle and the nature of exciton localization in moiré materials.
• Design Implications: The discovery of phonon-induced suppression of higher bands is crucial for designing optoelectronic devices. It suggests that simply increasing the twist angle to access higher energy bands may not yield the expected absorption features if phonon scattering becomes efficient.
• Methodological Advance: The application of the TCL master equation with gPSD analysis offers a robust tool for simulating non-Markovian dynamics in complex, tunable solid-state systems, moving beyond simple phenomenological broadening models.

In summary, this work demonstrates that the optical properties of moiré excitons are not static but are dynamically tuned by the twist angle, which dictates whether the system is governed by non-Markovian localization physics or Markovian delocalization physics, with profound consequences for absorption line shapes and the visibility of higher excitonic bands.