Real-time exciton dynamics in two-dimensional materials under ultrashort laser pulses

This study theoretically investigates the real-time exciton dynamics in two-dimensional h-BN and GeS monolayers under ultrashort laser pulses by combining a full-electron LAPW+lo method with a time-dependent adiabatic $GW$ approximation to elucidate the critical role of many-body effects in shaping ultrafast optical responses.

The Gist

Imagine a world made of incredibly thin, flat sheets of material—so thin they are essentially two-dimensional. In this microscopic world, electrons (the tiny particles that carry electricity) and "holes" (the empty spaces they leave behind) don't just float around freely. Because the sheet is so thin, they are stuck together by a strong magnetic-like pull, forming pairs called excitons. Think of an exciton like a dance couple: the electron is the lead, the hole is the follow, and they are holding hands so tightly they move as one unit.

This paper is about watching these dance couples in real-time when you hit them with a super-fast, intense flash of light (a laser pulse).

The Problem: Why is this hard to study?

Usually, when scientists try to predict how these materials behave, they use a "lazy" method. They assume the dancers are just moving on their own, ignoring how tightly they are holding hands. This is like watching a ballroom dance and pretending the partners aren't touching. It works okay for slow, simple movements, but if you play loud, fast music (an ultrashort laser pulse), the dancers react in complex, synchronized ways that the "lazy" method misses completely.

In real life, these materials (like Hexagonal Boron Nitride or GeS) are used in future electronics and solar cells. To make better devices, we need to understand exactly how these "dance couples" react to light.

The Solution: A High-Definition Simulation

The authors of this paper built a super-accurate computer simulation to watch these dances happen in real-time.

1. The Stage: They chose two specific materials to study:

   • h-BN (Hexagonal Boron Nitride): Often called "white graphene." It's a very sturdy, wide-gap insulator (like a very high wall that electrons can't easily jump over).
   • GeS (Germanium Sulfide): A newer material that is great at moving electricity, promising for future solar panels.

2. The Camera: Instead of just guessing, they used a method called TD-aGW. Imagine this as a high-speed, slow-motion camera that doesn't just see the dancers, but also calculates the invisible tension in their arms (the electron-hole interaction) as they spin. This is much more accurate than previous methods.

3. The Music (The Laser): They hit the materials with two types of "music" (laser pulses):

   • One-Photon Beat: A single, strong flash of light.
   • Two-Photon Beat: A slightly different flash where the material has to absorb two tiny bits of light at once to get excited.

What They Discovered

1. The Quantum "Beats" (The Wobble)
When they hit the h-BN material with the right color of light, the excitons didn't just sit there. They started "beating" or wobbling back and forth.

• The Analogy: Imagine two pendulums connected by a spring. If you push one, they start swinging in a complex rhythm, sometimes in sync, sometimes out of sync. The paper shows that the electrons and holes are doing this exact thing. They are swapping energy back and forth between different "dance moves" (quantum states) about every 16 femtoseconds (a femtosecond is a quadrillionth of a second!).
• Why it matters: This "wobble" creates a specific electrical signal. If we can control this wobble, we could build computers that process information at the speed of light.

2. The Shape of the Dance
They found that the shape of the laser pulse matters.

• If they used a simple, straight-line laser pulse, the dancers stayed in a predictable pattern.
• If they changed the pulse to be more complex (using two photons), the "dance floor" got crowded with more types of dancers. The neat, rhythmic wobble got messy and distorted. This tells scientists that to control these materials, they need to be very precise with the "music" they play.

3. The "Ghost" Dancers
In the two-photon experiment, they noticed that some dancers that are usually invisible (called "dark states") suddenly started moving. It's like a magician pulling a rabbit out of a hat that you didn't know was there. This is crucial because these "ghost" dancers might hold the key to new types of light-emitting devices.

The Big Picture

This paper is like a masterclass in choreography for the microscopic world. By using a highly advanced mathematical "camera," the authors showed us that when you shine a fast laser on these 2D materials, the electrons and holes don't just jump; they perform a complex, synchronized ballet.

Why should you care?
Understanding this ballet helps engineers design:

• Faster Computers: Using light instead of electricity to switch bits.
• Better Solar Cells: Capturing more energy from the sun by understanding how light creates these electron pairs.
• New Sensors: Devices that can detect light with extreme sensitivity.

In short, the authors didn't just look at the material; they listened to the rhythm of the electrons and learned how to conduct the orchestra for the technology of tomorrow.

Technical Summary

Here is a detailed technical summary of the paper "Real-time exciton dynamics in two-dimensional materials under ultrashort laser pulses" by Tumakov and Popova-Gorelova.

1. Problem Statement

The optical response of two-dimensional (2D) materials is dominated by excitonic effects due to reduced dielectric screening, leading to strong electron-hole Coulomb interactions and stable excitons at room temperature. While ultrafast spectroscopic techniques (e.g., pump-probe, high-harmonic generation) can probe these dynamics, the fundamental mechanisms governing exciton formation and real-time evolution—particularly in the nonlinear response regime—remain poorly understood.

Existing theoretical methods face trade-offs:

• Frequency-domain methods (e.g., solving the Bethe-Salpeter Equation, BSE) are accurate for linear response but struggle with strong-field, time-dependent nonlinear dynamics.
• Real-time methods (e.g., RT-TDDFT) are computationally efficient but often lack the rigorous many-body treatment required to capture strong excitonic binding energies and correlations accurately.

There is a need for a computationally efficient, all-electron, real-time framework that incorporates many-body electron-hole interactions (beyond the independent particle approximation) to simulate ultrafast charge dynamics in 2D materials under intense laser pulses.

2. Methodology

The authors developed and implemented a Real-Time Time-Dependent Adiabatic GW (RT-TD-aGW) approach within the all-electron exciting package.

• Theoretical Framework:

  • Effective Schrödinger Equation: They utilize an effective Schrödinger equation approach based on the "modern theory of polarization" (Berry phase). This avoids the explicit propagation of two-particle Green's functions or the full BSE kernel, offering a computationally efficient real-time scheme.
  • Many-Body Interactions: Electron-hole interactions are incorporated via a non-local self-energy operator derived from Many-Body Perturbation Theory (MBPT). Specifically, they use the Time-Dependent adiabatic GW (TD-aGW) approximation (equivalent to the Time-Dependent Hartree + Screened Exchange, or HSEX, approximation).
  • Hamiltonian: The system Hamiltonian includes a field-free part ($H_0$) with a scissor operator to correct the bandgap (mimicking GW quasiparticle corrections) and a time-dependent interaction term ($V_{ext}$) coupling the system to the external electric field.
  • Gauge: The length gauge is used for field coupling, allowing non-local operators (like the self-energy) to be added directly without gauge transformation complications.

• Numerical Implementation:

  • Basis Set: The Linearized Augmented Plane Wave plus Local Orbitals (LAPW+lo) basis set is employed. This allows for a highly accurate description of both delocalized and localized states, as well as core electrons (crucial for future XAS applications).
  • Propagation: The time-dependent Kohn-Sham-Bloch states are expanded in terms of ground-state orbitals. The expansion coefficients are evolved using a 4th-order Runge-Kutta (RK4) method.
  • Parallelization: Real-time communication between MPI processes is optimized by pre-calculating overlap matrices and exchanging only coefficient matrices, reducing computational overhead.
  • Materials Studied:
    1. Hexagonal Boron Nitride (h-BN): A wide-gap insulator (~6-7 eV) with strong excitonic effects ("white graphene").
    2. Germanium Sulfide (GeS): A moderate-gap semiconductor (~2.3 eV) with high carrier mobility.

3. Key Contributions

• Implementation of RT-TD-aGW in exciting: This is the first implementation of a real-time, many-body excitonic dynamics approach using the LAPW+lo all-electron basis set. Previous implementations (e.g., in YAMBO) used plane waves and pseudopotentials.
• Bridging Linear and Nonlinear Regimes: The method successfully captures excitonic effects in both linear response (where it reproduces BSE results) and strong-field nonlinear regimes (one- and two-photon excitation).
• Analysis of Quantum Beats: The study provides a detailed theoretical analysis of charge migration and quantum beats between different excitonic states (e.g., 1s, 2s, 2p) driven by ultrashort pulses.

4. Results

A. Ground State and Linear Response

• h-BN: The calculated dielectric function shows three distinct excitonic peaks (labeled a, b, c) within the bandgap, corresponding to 1s, 2p, and 2s-like states. The TD-aGW results are nearly indistinguishable from direct BSE solutions, validating the method's accuracy in the linear regime.
• GeS: Similar to h-BN, the TD-aGW approach accurately reproduces the excitonic spectrum, capturing peaks that are absent in the Independent Particle Approximation (IPA).

B. Real-Time Charge Dynamics

• h-BN (One-Photon Excitation):
  • Excitation with a 6.12 eV pulse (resonant with 2s and 2p excitons) induces quantum beats in the macroscopic polarization with a period of ~16 fs.
  • This corresponds to charge migration between 2s and 2p excitonic states in both momentum and coordinate space.
  • The charge density oscillates between Boron and Nitrogen sites.
• h-BN (Two-Photon Excitation):
  • Using a 3.06 eV pulse (two-photon resonance) excites different selection rules, including dark $A'_1$ states.
  • The quantum beat pattern is distorted due to the involvement of more than two excitonic states, and the polarization phase shifts by $\pi/2$ relative to the linear response.
• GeS Dynamics:
  • Excitation at 2.01 eV (linear) and 1 eV (nonlinear) reveals that the energy difference between excitonic states is comparable to their energy from the ground state.
  • Unlike h-BN, the GeS dynamics do not show clear beatings; instead, a "trembling" distribution is observed where multiple states are simultaneously populated.

C. Role of Electron-Hole Interactions

• Comparing TD-aGW (with interactions) vs. IPA (without interactions):
  • In the linear regime, IPA predicts negligible transitions below the bandgap, whereas TD-aGW shows significant absorption due to excitonic states.
  • In the nonlinear regime, excitonic effects significantly alter the amplitude and phase of the macroscopic polarization compared to IPA.
  • The imaginary part of the dielectric function manifests as a phase shift in the oscillation of the polarization.

5. Significance

• Fundamental Understanding: The study elucidates how many-body effects shape ultrafast excitonic processes, specifically demonstrating how symmetry breaking in the driving field can suppress or modify quantum beat scenarios.
• Methodological Advancement: By combining the high precision of the LAPW+lo basis with the efficiency of the real-time Berry-phase approach, the authors provide a robust tool for simulating ultrafast phenomena in 2D materials. This is particularly important for interpreting emerging experiments in attosecond science and high-harmonic generation.
• Application Potential: The ability to model core-electron effects and accurate bandgaps makes this approach suitable for interpreting future ultrafast X-ray absorption spectroscopy (XAS) experiments, aiding the design of optoelectronic and photonic devices based on 2D materials.

In conclusion, this work establishes a rigorous, all-electron theoretical framework for simulating non-equilibrium exciton dynamics, successfully capturing the complex interplay between light and matter in 2D materials beyond the independent-particle picture.