Nonequilibrium Electron Dynamics Research Articles

Hot carrier transfer from plasmon decay in Ag20at H-Si(111) surface: real-time TDDFT simulation in Wannier gauge.

Plasmon decay is believed to play an essential role in inducing hot carrier transfer at the interfaces between plasmonic nanoparticles and semiconductor surfaces. In this work, we employ real-time time-dependent density functional theory (RT-TDDFT) simulation in the Wannier gauge to gain quantum-mechanical insights into the nonlinear dynamics of the plasmon decay in the Ag20nanoparticle at a semiconductor surface. The first-principles simulations show that the plasmon decay is more than two times faster when the Ag20nanoparticle is adsorbed on a hydrogen-terminated Si(111) surface, taking place within 100 femtoseconds of the plasmon excitation. Hot carrier transfer across the interface is observed as the plasmon decay takes place, and nearly 30% of holes are generated deep in the valence band of the semiconductor surface. The use of Wannier gauge in RT-TDDFT simulation is particularly convenient for gaining quantum-mechanical insights into non-equilibrium electron dynamics in complex heterogeneous systems.

Hot Electron Cooling in n-Doped Colloidal Nanoplatelets Following Near-Infrared Intersubband Excitation.

Intersubband transition was recently discovered in colloidal nanoplatelets, but the associated intersubband carrier relaxation dynamics remains poorly understood. In particular, it is crucial to selectively excite the intersubband transition and to follow the hot electron dynamics in the absence of valence-band holes. This is achieved herein by exciting the predoped electrons in CdSe/ZnS nanoplatelets using near-infrared femtosecond pulses and monitoring nonequilibrium electron dynamics using broad-band visible pulses. We find that the n = 2 electrons relax to the n = 1 subband and establish a Fermi-Dirac distribution within 200 fs, and finally reach an equilibrium with the lattice within a few ps. The cooling dynamics depend mainly on the excitation fluence but weakly on the doping density and the lattice temperature. These characteristics are well captured by our numerical simulation that explicitly accounts for the state occupation effect and optical phonon scattering.

Toward Direct Exploration of the Few-Femtosecond Dynamics of Electronic Coherence and Correlation in Quantum Materials Using Time- and Angle-Resolved Photoemission Spectroscopy

Over the last two decades, time- and angle-resolved photoemission spectroscopy (trARPES) has become a mature and established experimental technique for the study of ultrafast electronic and structural dynamics in materials. To date, most trARPES investigations have focused on the investigation of processes occurring on time scales of ≳30 fs, in particular, relaxation and thermalization, and have therefore been blind to the initial sub-10 fs dynamics related to electronic coherence and correlation effects. In this article, we illustrate how current trARPES setups reach their limits when it comes to addressing such extraordinarily short time scales and present an experimental configuration that provides the time, energy, and momentum resolutions required to monitor few-femtosecond dynamics on the relevant energy and momentum scales. We discuss the potential capabilities of such an experiment to study the electronic response of materials in the strong-field interaction regime at PHz frequencies and finally review a theoretical concept that may in the future even overcome the competing resolution limitations of trARPES experiments, as imposed by the time–bandwidth product of the probing laser pulse. Our roadmap for ultrafast trARPES indicates a path to break new experimental ground in quantum nonequilibrium electronic dynamics, from which new possibilities for ultrafast control of optical and electronic signals in quantum materials can be explored.

Real-Time Time-Dependent Density Functional Theory for Simulating Nonequilibrium Electron Dynamics.

The explicit real-time propagation approach for time-dependent density functional theory (RT-TDDFT) has increasingly become a popular first-principles computational method for modeling various time-dependent electronic properties of complex chemical systems. In this Perspective, we provide a nontechnical discussion of how this first-principles simulation approach has been used to gain novel physical insights into nonequilibrium electron dynamics phenomena in recent years. Following a concise overview of the RT-TDDFT methodology from a practical standpoint, we discuss our recent studies on the electronic stopping of DNA in water and the Floquet topological phase as examples. Our discussion focuses on how RT-TDDFT simulations played a unique role in deriving new scientific understandings. We then discuss existing challenges and some new advances at the frontier of RT-TDDFT method development for studying increasingly complex dynamic phenomena and systems.

Real-time observation of non-equilibrium phonon-electron energy and angular momentum flow in laser-heated nickel.

Identifying the microscopic nature of non-equilibrium energy transfer mechanisms among electronic, spin, and lattice degrees of freedom is central to understanding ultrafast phenomena such as manipulating magnetism on the femtosecond timescale. Here, we use time- and angle-resolved photoemission spectroscopy to go beyond the often-used ensemble-averaged view of non-equilibrium dynamics in terms of quasiparticle temperature evolutions. We show for ferromagnetic Ni that the non-equilibrium electron and spin dynamics display pronounced variations with electron momentum, whereas the magnetic exchange interaction remains isotropic. This highlights the influence of lattice-mediated scattering processes and opens a pathway toward unraveling the still elusive microscopic mechanism of spin-lattice angular momentum transfer.

First-Principles Approach for Coupled Quantum Dynamics of Electrons and Protons in Heterogeneous Systems.

The coupled quantum dynamics of electrons and protons is ubiquitous in many dynamical processes involving light-matter interaction, such as solar energy conversion in chemical systems and photosynthesis. A first-principles description of such nuclear-electronic quantum dynamics requires not only the time-dependent treatment of nonequilibrium electron dynamics but also that of quantum protons. Quantum mechanical correlation between electrons and protons adds further complexity to such coupled dynamics. Here we extend real-time nuclear-electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) to periodic systems and perform first-principles simulations of coupled quantum dynamics of electrons and protons in complex heterogeneous systems. The process studied is an electronically excited-state intramolecular proton transfer of o-hydroxybenzaldehyde in water and at a silicon (111) semiconductor-molecule interface. These simulations illustrate how environments such as hydrogen-bonding water molecules and an extended material surface impact the dynamical process on the atomistic level. Depending on how the molecule is chemisorbed on the surface, excited-state electron transfer from the molecule to the semiconductor surface can inhibit ultrafast proton transfer within the molecule. This Letter elucidates how heterogeneous environments influence the balance between the quantum mechanical proton transfer and excited electron dynamics. The periodic RT-NEO-TDDFT approach is applicable to a wide range of other photoinduced heterogeneous processes.

A multi-fragment real-time extension of projected density matrix embedding theory: Non-equilibrium electron dynamics in extended systems.

In this work, we derive a multi-fragment real-time extension of the projected density matrix embedding theory (pDMET) designed to treat non-equilibrium electron dynamics in strongly correlated systems. As in the previously developed static pDMET, the real time pDMET partitions the total system into many fragments; the coupling between each fragment and the rest of the system is treated through a compact representation of the environment in terms of a quantum bath. The real-time pDMET involves simultaneously propagating the wavefunctions for each separate fragment-bath embedding system along with an auxiliary mean-field wavefunction of the total system. The equations of motion are derived by (i) projecting the time-dependent Schrödinger equation in the fragment and bath space associated with each separate fragment and by (ii) enforcing the pDMET matching conditions between the global 1-particle reduced density matrix (1-RDM) obtained from the fragment calculations and the mean-field 1-RDM at all points in time. The accuracy of the method is benchmarked through comparisons to time-dependent density-matrix renormalization group and time-dependent Hartree-Fock (TDHF) theory; the methods were applied to a one- and two-dimensional single-impurity Anderson model and multi-impurity Anderson models with ordered and disordered distributions of the impurities. The results demonstrate a large improvement over TDHF and rapid convergence to the exact dynamics with an increase in fragment size. Our results demonstrate that the real-time pDMET is a promising and flexible method that balances accuracy and efficiency to simulate the non-equilibrium electron dynamics in heterogeneous systems of large size.

Excited electron and spin dynamics in topological insulator: A perspective from ab initio non-adiabatic molecular dynamics

We perform an ab initio non-adiabatic molecular dynamics simulation to investigate the non-equilibrium spin and electron dynamics in a prototypical topological insulator (TI) Bi2Se3. Different from the ground state, we reveal that backscattering can happen in an oscillating manner between time-reversal pair topological surface states (TSSs) in the non-equilibrium dynamics. Analysis shows the phonon excitation induces orbital composition change by electron-phonon interaction, which further stimulates spin canting through spin-orbit coupling. The spin canting of time-reversal pair TSSs leads to the non-zero non-adiabatic coupling between them and then issues in backscattering. Both the spin canting and backscattering result in ultrafast spin relaxation with a timescale around 100 fs. This study provides critical insights into the non-equilibrium electron and spin dynamics in TI at the ab initio level and paves a way for the design of ultrafast spintronic materials.

Extracting quantitative dielectric properties from pump-probe spectroscopy

Optical pump-probe spectroscopy is a powerful tool for the study of non-equilibrium electronic dynamics and finds wide applications across a range of fields, from physics and chemistry to material science and biology. However, a shortcoming of conventional pump-probe spectroscopy is that photoinduced changes in transmission, reflection and scattering can simultaneously contribute to the measured differential spectra, leading to ambiguities in assigning the origin of spectral signatures and ruling out quantitative interpretation of the spectra. Ideally, these methods would measure the underlying dielectric function (or the complex refractive index) which would then directly provide quantitative information on the transient excited state dynamics free of these ambiguities. Here we present and test a model independent route to transform differential transmission or reflection spectra, measured via conventional optical pump-probe spectroscopy, to changes in the quantitative transient dielectric function. We benchmark this method against changes in the real refractive index measured using time-resolved Frequency Domain Interferometry in prototypical inorganic and organic semiconductor films. Our methodology can be applied to existing and future pump-probe data sets, allowing for an unambiguous and quantitative characterisation of the transient photoexcited spectra of materials. This in turn will accelerate the adoption of pump-probe spectroscopy as a facile and robust materials characterisation and screening tool.

Nonlocal and nonadiabatic Pauli potential for time-dependent orbital-free density functional theory

Time-dependent orbital-free density functional theory (TD-OFDFT) is an efficient ab-initio method for calculating the electronic dynamics of large systems. In comparison to standard TD-DFT, it computes only a single electronic state regardless of system size, but it requires an additional time-dependent Pauli potential term. We propose a nonadiabatic and nonlocal Pauli potential whose main ingredients are the time-dependent particle and current densities. Our calculations of the optical spectra of metallic and semiconductor clusters indicate that nonlocal and nonadiabatic TD-OFDFT performs accurately for metallic systems and semiquantitatively for semiconductors. This work opens the door to wide applicability of TD-OFDFT for nonequilibrium electron and electron-nuclear dynamics of materials.

Subpicosecond metamagnetic phase transition in FeRh driven by non-equilibrium electron dynamics

Femtosecond light-induced phase transitions between different macroscopic orders provide the possibility to tune the functional properties of condensed matter on ultrafast timescales. In first-order phase transitions, transient non-equilibrium phases and inherent phase coexistence often preclude non-ambiguous detection of transition precursors and their temporal onset. Here, we present a study combining time-resolved photoelectron spectroscopy and ab-initio electron dynamics calculations elucidating the transient subpicosecond processes governing the photoinduced generation of ferromagnetic order in antiferromagnetic FeRh. The transient photoemission spectra are accounted for by assuming that not only the occupation of electronic states is modified during the photoexcitation process. Instead, the photo-generated non-thermal distribution of electrons modifies the electronic band structure. The ferromagnetic phase of FeRh, characterized by a minority band near the Fermi energy, is established 350 ± 30 fs after the laser excitation. Ab-initio calculations indicate that the phase transition is initiated by a photoinduced Rh-to-Fe charge transfer.

Understanding the mechanism of tunable-barrier single-electron pumping: Mechanism crossover and optimal accuracy

Understanding nonequilibrium electron dynamics in a tunable-barrier dynamic quantum dot is vital to achieve high-accuracy single-electron pumping, which has applications to metrological standards and quantum information devices. However, the dynamic mechanism determining the pumping process has not been fully understood. In this paper, we study and clarify the physical mechanisms of single-electron pumping by analyzing and calculating master equations with a realistic model in detail. Our focus is mainly on the regime which lies at the crossover between the two extreme physical models of single-electron pumping. The mechanism crossover strongly depends on a capacitive coupling between a gate for tuning the barrier and the quantum dot. We find that there is an optimal value of the coupling with the best pumping accuracy and the mechanism crossover depends on gate voltage in some range of the coupling values. Our results offer a guideline for evaluation of pumping characteristics and for optimization of the device structure, which is an important step toward understanding high-speed nonequilibrium electron dynamics and realizing reproducible high-accuracy single-electron pumps.

First-principles calculations for attosecond electron dynamics in solids

Nonequilibrium electron dynamics in solids is an important subject from both fundamental and technological points of view. The recent development of laser technology has enabled us to study ultrafast electron dynamics in the time domain. First-principles calculation is a powerful tool for analyzing such complex electron dynamics and clarifying the physics behind the experimental observation. In this article, we review the recent development of the first-principles calculation for light-induced electron dynamics in solids by revising its application to recent attosecond experiments. The electron dynamics calculations offer an accurate description of static and transient optical properties of solids and provide physics insight into light-induced electron dynamics. Furthermore, the microscopic decomposition of transient properties of nonequilibrium systems has been developed to extract microscopic information from the simulation results. The first-principles analysis opened a novel path to analyze the nonequilibrium electron dynamics in matter and to provide the fundamental understanding complementarily with the sophisticated experimental technique.

Perturbo: A software package for ab initio electron–phonon interactions, charge transport and ultrafast dynamics

Perturbo is a software package for first-principles calculations of charge transport and ultrafast carrier dynamics in materials. The current version focuses on electron–phonon interactions and can compute phonon-limited transport properties such as the conductivity, carrier mobility and Seebeck coefficient. It can also simulate the ultrafast nonequilibrium electron dynamics in the presence of electron–phonon scattering. Perturbo uses results from density functional theory and density functional perturbation theory calculations as input, and employs Wannier interpolation to reduce the computational cost. It supports norm-conserving and ultrasoft pseudopotentials, spin–orbit coupling, and polar electron–phonon interactions for bulk and 2D materials. Hybrid MPI plus OpenMP parallelization is implemented to enable efficient calculations on large systems (up to at least 50 atoms) using high-performance computing. Taken together, Perturbo provides efficient and broadly applicable ab initio tools to investigate electron–phonon interactions and carrier dynamics quantitatively in metals, semiconductors, insulators, and 2D materials. Program summaryProgram Title:PerturboCPC Library link to program files:https://doi.org/10.17632/34m2p6v79t.1Developer’s repository link:https://perturbo-code.github.ioLicensing provisions: GNU General Public Licence 3.0Programming language:Fortran, PythonExternal routines/libraries:LAPACK, HDF5, MPI, OpenMP, FFTW, Quantum-ESPRESSO, Wannier90Nature of problem: Computing transport properties from first-principles in materials, including the electrical conductivity, carrier mobility and Seebeck coefficient; Simulating ultrafast nonequilibrium electron dynamics, such as the relaxation of excited carriers via interactions with phonons.Solution method: We implement the first-principles Boltzmann transport equation, which employs materials properties such as the electronic structure, lattice dynamics, and electron–phonon collision terms computed with density functional theory and density functional perturbation theory. The Boltzmann transport equation is solved numerically to compute charge transport and simulate ultrafast carrier dynamics. Wannier interpolation is employed to reduce the computational cost.Additional comments: Hybrid MPI plus OpenMP parallelization is implemented to run large calculations and take advantage of high-performance computing. Most results are output to HDF5 file format, which is portable and convenient for post-processing using high-level languages such as Python and Julia.

Dynamical transition orbitals: A particle-hole description in real-time TDDFT dynamics.

We expand the concept of natural transition orbitals in the context of real-time time-dependent density functional theory (RT-TDDFT) and show its application in practical calculations. Kohn-Sham single-particle wavefunctions are propagated in RT-TDDFT simulation, and physical properties remain invariant under their unitary transformation. In this work, we exploit this gauge freedom and expand the concept of natural transition orbitals, which is widely used in linear-response TDDFT, for obtaining a particle-hole description in RT-TDDFT simulation. While linear-response TDDFT is widely used to study electronic excitation, RT-TDDFT can be employed more generally to simulate non-equilibrium electron dynamics. Studying electron dynamics in terms of dynamic transitions of particle-hole pairs is, however, not straightforward in the RT-TDDFT simulation. By constructing natural transition orbitals through projecting time-dependent Kohn-Sham wave functions onto occupied/unoccupied eigenstate subspaces, we show that linear combinations of a pair of the resulting hole/particle orbitals form a new gauge, which we refer to as dynamical transition orbitals. We demonstrate the utility of this framework to analyze RT-TDDFT simulations of optical excitation and electronic stopping dynamics in the particle-hole description.

Bypassing the Structural Bottleneck in the Ultrafast Melting of Electronic Order.

Impulsive optical excitation generally results in a complex nonequilibrium electron and lattice dynamics that involves multiple processes on distinct timescales, and a common conception is that for times shorter than about 100fs the gap in the electronic spectrum is not seriously affected by lattice vibrations. Here, however, by directly monitoring the photoinduced collapse of the spectral gap in a canonical charge-density-wave material, the blue bronze Rb_{0.3}MoO_{3}, we find that ultrafast (∼60 fs) vibrational disordering due to efficient hot-electron energy dissipation quenches the gap significantly faster than the typical structural bottleneck time corresponding to one half-cycle oscillation (∼315 fs) of the coherent charge-density-wave amplitude mode. This result not only demonstrates the importance of incoherent lattice motion in the photoinduced quenching of electronic order, but also resolves the perennial debate about the nature of the spectral gap in a coupled electron-lattice system.

Nonequilibrium Electron and Lattice Dynamics of Strongly Correlated Quantum Materials

Nonequilibrium Electron and Lattice Dynamics of Strongly Correlated Quantum Materials

Tunable Onset of Hydrogen Evolution in Graphene with Hot Electrons.

Here, we show that the turn-on voltage for the hydrogen evolution reaction on a graphene surface can be tuned in a semiconductor-insulator-graphene (SIG) device immersed in a solution. Specifically, it is shown that the hydrogen evolution reaction (HER) onset for the graphene can shift by >0.8 V by application of a voltage across a graphene-Al2O3-silicon junction. We show that this shift occurs due to the creation of a hot electron population in graphene due to tunneling from the Si to graphene. Through control experiments, we show that the presence of the graphene is necessary for this behavior. By analyzing the silicon, graphene, and solution current components individually, we find an increase in the silicon current despite a fixed graphene-silicon voltage, corresponding to an increase in the HER current. This additional silicon current appears to directly drive the electrochemical reaction, without modifying the graphene current. We term this current "direct injection current" and hypothesize that this current occurs due to electrons injected from the silicon into graphene that drives the HER before any electron-electron scattering occurs in the graphene. To further determine whether hot electrons injected at different energies could explain the observed total solution current, the nonequilibrium electron dynamics was studied using a 2D ensemble Monte Carlo Boltzmann transport equation (MCBTE) solver. By rigorously considering the key scattering mechanisms, we show that the injected hot electrons can significantly increase the available electron flux at high energies. These results show that semiconductor-insulator-graphene devices are a platform which can tune the electrochemical reaction rate via multiple mechanisms.

Vibrational wavepacket dynamics in Fe carbene photosensitizer determined with femtosecond X-ray emission and scattering

The non-equilibrium dynamics of electrons and nuclei govern the function of photoactive materials. Disentangling these dynamics remains a critical goal for understanding photoactive materials. Here we investigate the photoinduced dynamics of the [Fe(bmip)2]2+ photosensitizer, where bmip = 2,6-bis(3-methyl-imidazole-1-ylidine)-pyridine, with simultaneous femtosecond-resolution Fe Kα and Kβ X-ray emission spectroscopy (XES) and X-ray solution scattering (XSS). This measurement shows temporal oscillations in the XES and XSS difference signals with the same 278 fs period oscillation. These oscillations originate from an Fe-ligand stretching vibrational wavepacket on a triplet metal-centered (3MC) excited state surface. This 3MC state is populated with a 110 fs time constant by 40% of the excited molecules while the rest relax to a 3MLCT excited state. The sensitivity of the Kα XES to molecular structure results from a 0.7% average Fe-ligand bond length shift between the 1 s and 2p core-ionized states surfaces.

Picosecond detection of electron motion.

The motion of a single electron can now be sampled with picosecond resolution, which helps to characterize and understand non-equilibrium electron dynamics in nanoscale conductors.