Competition between the electronic and phonon–mediated scattering channels in the out–of–equilibrium carrier dynamics of semiconductors: an ab-initio approach

Ab Initio Predictions of Spin Relaxation, Dephasing, and Diffusion in Solids

Spin relaxation, dephasing, and diffusion are at the heart of spin-based information technology. Accurate theoretical approaches to simulate spin lifetimes (τs), determining how fast the spin polarization and phase information will be lost, are important to the understanding of the underlying mechanism of these spin processes, and invaluable in searching for promising candidates of spintronic materials. Recently, we develop a first-principles real-time density-matrix (FPDM) approach to simulate spin dynamics for general solid-state systems. Through the complete first-principles descriptions of light-matter interaction and scattering processes including electron-phonon, electron-impurity, and electron-electron scatterings with self-consistent spin-orbit coupling, as well as ab initio Landé g-factor, our method can predict τs at various conditions as a function of carrier density and temperature, under electric and magnetic fields. By employing this method, we successfully reproduce experimental results of disparate materials and identify the key factors affecting spin relaxation, dephasing, and diffusion in different materials. Specifically, we predict that germanene has long τs (∼100 ns at 50 K), a giant spin lifetime anisotropy, and spin-valley locking effect under electric fields, making it advantageous for spin-valleytronic applications. Based on our theoretical derivations and ab initio simulations, we propose a new useful electronic quantity, named spin-flip angle θ↑↓, for the understanding of spin relaxation through intervalley spin-flip scattering processes. Our method can be further applied to other emerging materials and extended to simulate exciton spin dynamics and steady-state photocurrents due to photogalvanic effect.

Real-Time GW-Ehrenfest-Fan-Migdal Method for Nonequilibrium 2D Materials

Quantum simulations of photoexcited low-dimensional systems are pivotal for understanding how to functionalize and integrate novel two-dimensional (2D) materials in next-generation optoelectronic devices. First-principles predictions are extremely challenging due to the simultaneous interplay of light-matter, electron-electron, and electron-nuclear interactions. We here present an advanced ab initio many-body method that accounts for quantum coherence and non-Markovian effects while treating electrons and nuclei on equal footing, thereby preserving fundamental conservation laws like the total energy. The impact of this advancement is demonstrated through real-time simulations of the complex multivalley dynamics in a molybdenum disulfide (MoS2) monolayer pumped above gap. Within a single framework, we provide a parameter-free description of the coherent-to-incoherent crossover, elucidating the role of microscopic and collective excitations in the dephasing and thermalization processes.

Real-Time GW: Toward an Ab Initio Description of the Ultrafast Carrier and Exciton Dynamics in Two-Dimensional Materials

We demonstrate the feasibility of the time-linear scaling formulation of the GW method [Phys. Rev. Lett. 124, 076601 (2020)PRLTAO0031-900710.1103/PhysRevLett.124.076601] for ab initio simulations of optically driven two-dimensional materials. The time-dependent GW equations are derived and solved numerically in the basis of Bloch states. We address carrier multiplication and relaxation in photoexcited graphene and find deviations from the typical exponential behavior predicted by the Markovian Boltzmann approach. For a resonantly pumped semiconductor we discover a self-sustained screening cascade leading to the Mott transition of coherent excitons. Our results draw attention to the importance of non-Markovian and dynamical screening effects in out-of-equilibrium phenomena.

Time-Dependent Screening Explains the Ultrafast Excitonic Signal Rise in 2D Semiconductors

We calculate the time evolution of the transient reflection signal in an MoS₂ monolayer on a SiO₂/Si substrate using first-principles out-of-equilibrium real-time methods. Our simulations provide a simple and intuitive physical picture for the delayed, yet ultrafast, evolution of the signal whose rise time depends on the excess energy of the pump laser: at laser energies above the A- and B-exciton, the pump pulse excites electrons and holes far away from the K valleys in the first Brillouin zone. Electron-phonon and hole-phonon scattering lead to a gradual relaxation of the carriers toward small Active Excitonic Regions around K, enhancing the dielectric screening. The accompanying time-dependent band gap renormalization dominates over Pauli blocking and the excitonic binding energy renormalization. This explains the delayed buildup of the transient reflection signal of the probe pulse, in excellent agreement with recent experimental data. Our results show that the observed delay is not a unique signature of an exciton formation process but rather caused by coordinated carrier dynamics and its influence on the screening.

Ultrafast Relaxation of Symmetry-Breaking Photo-Induced Atomic Forces

We present a first-principles method for the calculation of the temperature-dependent relaxation of symmetry-breaking atomic driving forces in photoexcited systems. We calculate the phonon-assisted decay of the photoexcited force on the low-symmetry E_{g} mode following absorption of an ultrafast pulse in Bi, Sb, and As. The force decay lifetimes for Bi and Sb are of the order of 10 fs and in agreement with recent experiments, demonstrating that electron-phonon scattering is the primary mechanism relaxing the symmetry-breaking forces. Calculations for a range of absorbed photon energies suggest that larger amplitude, symmetry-breaking atomic motion may be induced by choosing a pump photon energy which maximizes the product of the initial E_{g} force and its lifetime. The high-symmetry A_{1g} force undergoes a partial decay to a nonzero constant on similar timescales, which has not yet been measured in experiments. The average imaginary part of the electron self-energy over the photoexcited carrier distribution provides a crude indication of the decay rate of symmetry-breaking forces.

Many-body perturbation theory calculations using the yambo code

yambo is an open source project aimed at studying excited state properties of condensed matter systems from first principles using many-body methods. As input, yambo requires ground state electronic structure data as computed by density functional theory codes such as Quantum ESPRESSO and Abinit. yambo's capabilities include the calculation of linear response quantities (both independent-particle and including electron-hole interactions), quasi-particle corrections based on the GW formalism, optical absorption, and other spectroscopic quantities. Here we describe recent developments ranging from the inclusion of important but oft-neglected physical effects such as electron-phonon interactions to the implementation of a real-time propagation scheme for simulating linear and non-linear optical properties. Improvements to numerical algorithms and the user interface are outlined. Particular emphasis is given to the new and efficient parallel structure that makes it possible to exploit modern high performance computing architectures. Finally, we demonstrate the possibility to automate workflows by interfacing with the yambopy and AiiDA software tools.

Spin States Protected from Intrinsic Electron-Phonon Coupling Reaching 100 ns Lifetime at Room Temperature in MoSe2

We present time-resolved Kerr rotation measurements, showing spin lifetimes of over 100 ns at room temperature in monolayer MoSe₂. These long lifetimes are accompanied by an intriguing temperature-dependence of the Kerr amplitude, which increases with temperature up to 50 K and then abruptly switches sign. Using ab initio simulations, we explain the latter behavior in terms of the intrinsic electron-phonon coupling and the activation of transitions to secondary valleys. The phonon-assisted scattering of the photoexcited electron-hole pairs prepares a valley spin polarization within the first few ps after laser excitation. The sign of the total valley magnetization, and thus the Kerr amplitude, switches as a function of temperature, as conduction and valence band states exhibit different phonon-mediated intervalley scattering rates. However, the electron-phonon scattering on the ps time scale does not provide an explanation for the long spin lifetimes. Hence, we deduce that the initial spin polarization must be transferred into spin states, which are protected from the intrinsic electron-phonon coupling, and are most likely resident charge carriers, which are not part of the itinerant valence or conduction band states.

Intravalley Spin-Flip Relaxation Dynamics in Single-Layer WS2

In monolayer (1L) transition metal dichalcogenides (TMDs) the valence and conduction bands are spin-split because of the strong spin-orbit interaction. In tungsten-based TMDs the spin-ordering of the conduction band is such that the so-called dark excitons, consisting of electrons and holes with opposite spin orientation, have lower energy than A excitons. The transition from bright to dark excitons involves the scattering of electrons from the upper to the lower conduction band at the K point of the Brillouin zone, with detrimental effects for the optoelectronic response of 1L-TMDs, since this reduces their light emission efficiency. Here, we exploit the valley selective optical selection rules and use two-color helicity-resolved pump-probe spectroscopy to directly measure the intravalley spin-flip relaxation dynamics in 1L-WS₂. This occurs on a sub-ps time scale, and it is significantly dependent on temperature, indicative of phonon-assisted relaxation. Time-dependent ab initio calculations show that intravalley spin-flip scattering occurs on significantly longer time scales only at the K point, while the occupation of states away from the minimum of the conduction band significantly reduces the scattering time. Our results shed light on the scattering processes determining the light emission efficiency in optoelectronic and photonic devices based on 1L-TMDs.

Hot electron relaxation dynamics in semiconductors: assessing the strength of the electron-phonon coupling from the theoretical and experimental viewpoints

The rapid development of the computational methods based on density functional theory, on the one hand, and of time-, energy-, and momentum-resolved spectroscopy, on the other hand, allows today an unprecedently detailed insight into the processes governing hot-electron relaxation dynamics, and, in particular, into the role of the electron-phonon coupling. Instead of focusing on the development of a particular method, theoretical or experimental, this review aims to treat the progress in the understanding of the electron-phonon coupling which can be gained from both, on the basis of recently obtained results. We start by defining several regimes of hot electron relaxation via electron-phonon coupling, with respect to the electron excitation energy. We distinguish between energy and momentum relaxation of hot electrons, and summarize, for several semiconductors of the IV and III-V groups, the orders of magnitude of different relaxation times in different regimes, on the basis of known experimental and numerical data. Momentum relaxation times of hot electrons become very short around 1 eV above the bottom of the conduction band, and such ultrafast relaxation mechanisms are measurable only in the most recent pump-probe experiments. Then, we give an overview of the recent progress in the experimental techniques allowing to obtain detailed information on the hot-electron relaxation dynamics, with the main focus on time-, energy-, and momentum-resolved photoemission experiments. The particularities of the experimental approach developed by one of us, which allows to capture time-, energy-, and momentum-resolved hot-electron distributions, as well as to measure momentum relaxation times of the order of 10 fs, are discussed. We further discuss the main advances in the calculation of the electron-phonon scattering times from first principles over the past ten years, in semiconducting materials. Ab initio techniques and efficient interpolation methods provide the possibility to calculate electron-phonon scattering times with high precision at reasonable numerical cost. We highlight the methods of analysis of the obtained numerical results, which allow to give insight into the details of the electron-phonon scattering mechanisms. Finally, we discuss the concept of hot electron ensemble which has been proposed recently to describe the hot-electron relaxation dynamics in GaAs, the applicability of this concept to other materials, and its limitations. We also mention some open problems.

Ab Initio Calculations of Ultrashort Carrier Dynamics in Two-Dimensional Materials: Valley Depolarization in Single-Layer WSe2

In single-layer WSe₂, a paradigmatic semiconducting transition metal dichalcogenide, a circularly polarized laser field can selectively excite electronic transitions in one of the inequivalent K^± valleys. Such selective valley population corresponds to a pseudospin polarization. This can be used as a degree of freedom in a "valleytronic" device provided that the time scale for its depolarization is sufficiently large. Yet, the mechanism behind the valley depolarization still remains heavily debated. Recent time-dependent Kerr experiments have provided an accurate way to visualize the valley dynamics by measuring the rotation of a linearly polarized probe pulse applied after a circularly polarized pump pulse. We present here a clear, accurate and parameter-free description of the valley dynamics. By using an atomistic, ab initio approach, we fully disclose the elemental mechanisms that dictate the depolarization effects. Our results are in excellent agreement with recent time-dependent Kerr experiments. We explain the Kerr dynamics and its temperature dependence in terms of electron-phonon-mediated processes that induce spin-flip intervalley transitions.

Photo-Induced Bandgap Renormalization Governs the Ultrafast Response of Single-Layer MoS2

Transition metal dichalcogenides (TMDs) are emerging as promising two-dimensional (2D) semiconductors for optoelectronic and flexible devices. However, a microscopic explanation of their photophysics, of pivotal importance for the understanding and optimization of device operation, is still lacking. Here, we use femtosecond transient absorption spectroscopy, with pump pulse tunability and broadband probing, to monitor the relaxation dynamics of single-layer MoS2 over the entire visible range, upon photoexcitation of different excitonic transitions. We find that, irrespective of excitation photon energy, the transient absorption spectrum shows the simultaneous bleaching of all excitonic transitions and corresponding red-shifted photoinduced absorption bands. First-principle modeling of the ultrafast optical response reveals that a transient bandgap renormalization, caused by the presence of photoexcited carriers, is primarily responsible for the observed features. Our results demonstrate the strong impact of many-body effects in the transient optical response of TMDs even in the low-excitation-density regime.