Photoexcitation in Solids: First-Principles Quantum Simulations by Real-Time TDDFT

Correlated electron-nuclear dynamics of photoinduced water dissociation on rutile TiO2

Elucidating the mechanism of photoinduced water splitting on TiO2 is important for advancing the understanding of photocatalysis and the ability to control photocatalytic surface reactions. However, incomplete experimental information and complex coupled electron-nuclear motion make the microscopic understanding challenging. Here we analyse the atomic-scale pathways of photogenerated charge carrier transport and photoinduced water dissociation at the prototypical water-rutile TiO2(110) interface using first-principles dynamics simulations. Two distinct mechanisms are observed. Field-initiated electron migration leads to adsorbed water dissociation via proton transfer to a surface bridging oxygen. In the other pathway, adsorbed water dissociation occurs via proton donation to a second-layer water molecule coupled to photoexcited-hole transfer promoted by in-plane surface lattice distortions. Two stages of non-adiabatic in-plane lattice motion-expansion and recovery-are observed, which are closely associated with population changes in Ti3d orbitals. Controlling such highly correlated electron-nuclear dynamics may provide opportunities for boosting the performance of photocatalytic materials.

Efficient exact exchange using Wannier functions and other related developments in planewave-pseudopotential implementation of RT-TDDFT

The plane-wave pseudopotential (PW-PP) formalism is widely used for the first-principles electronic structure calculation of extended periodic systems. The PW-PP approach has also been adapted for real-time time-dependent density functional theory (RT-TDDFT) to investigate time-dependent electronic dynamical phenomena. In this work, we detail recent advances in the PW-PP formalism for RT-TDDFT, particularly how maximally localized Wannier functions (MLWFs) are used to accelerate simulations using the exact exchange. We also discuss several related developments, including an anti-Hermitian correction for the time-dependent MLWFs (TD-MLWFs) when a time-dependent electric field is applied, the refinement procedure for TD-MLWFs, comparison of the velocity and length gauge approaches for applying an electric field, and elimination of long-range electrostatic interaction, as well as usage of a complex absorbing potential for modeling isolated systems when using the PW-PP formalism.

Condensed Matter Systems Exposed to Radiation: Multiscale Theory, Simulations, and Experiment

This roadmap reviews the new, highly interdisciplinary research field studying the behavior of condensed matter systems exposed to radiation. The Review highlights several recent advances in the field and provides a roadmap for the development of the field over the next decade. Condensed matter systems exposed to radiation can be inorganic, organic, or biological, finite or infinite, composed of different molecular species or materials, exist in different phases, and operate under different thermodynamic conditions. Many of the key phenomena related to the behavior of irradiated systems are very similar and can be understood based on the same fundamental theoretical principles and computational approaches. The multiscale nature of such phenomena requires the quantitative description of the radiation-induced effects occurring at different spatial and temporal scales, ranging from the atomic to the macroscopic, and the interlinks between such descriptions. The multiscale nature of the effects and the similarity of their manifestation in systems of different origins necessarily bring together different disciplines, such as physics, chemistry, biology, materials science, nanoscience, and biomedical research, demonstrating the numerous interlinks and commonalities between them. This research field is highly relevant to many novel and emerging technologies and medical applications.

The impact of a dopant atom on the distribution of hot electrons and holes in Au-doped Ag nano-clusters

The generation of hot carriers (HCs) through the excitation of localized surface plasmon resonance (LSPR) in metal nanostructures is a fascinating phenomenon that fuels both fundamental and applied research. In this study, we employ first principles real-time time-dependent density-functional theory (rt-TDDFT) calculations to elucidate the creation and distribution of HCs within Au-doped Ag nanoclusters: Ag11Cl3P7H21, Ag10AucoreCl3P7H21, and Ag10AusurfCl3P7H21 nanoclusters. Our findings indicate that adjustments in HC distribution are achievable through the Au dopant atom, and precise control of HC distribution is possible by manipulating the location of the Au dopant atom. When employing a Gaussian laser pulse tailored to match the LSPR frequency, a substantial accumulation of HCs in the Ag-P bond is observed. This finding suggests a weakening of the Ag-P bonds and, consequently, the initiation of bond stretching. We propose that these findings open up possibilities for tuning HCs in Au-doped chemically functionalized Ag nanoclusters.

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.

Coherent Phonon Assisted Ultrafast Order-Parameter Reversal and Hidden Metallic State in Ta_{2}NiSe_{5}

The nonequilibrium dynamics during photoinduced insulator-to-metal transition (IMT) in the excitonic insulator (EI) candidate Ta_{2}NiSe_{5} have been investigated, which reproduce the timescale and spectral features of the ultrafast switch and reveal intricate many-body interactions involving multidegrees of freedom. The key role of lattice order parameter (OP) reversal, occurring on a timescale comparable to that of purely electronic processes (<100  fs), is identified. This reversal is enabled by the anharmonic interactions between EI-OP-coupled phonons and the conventional coherent phonons, leading to a modified potential energy landscape and a high-frequency mode up-conversion. The phonon excitation depends on the dynamics of photocarriers distributed around the Fermi level, and thus intertwines with the excitonic quenching and the complete gap collapse. These findings provide a comprehensive understanding of exciton-phonon dynamics in correlated quantum materials.

Unraveling Hidden Charge Density Wave Phases in 1T-TiSe_{2}

The unexpected chiral order observed in 1T-TiSe_{2} represents an exciting area to explore chirality in condensed matter, while its microscopic mechanism remains elusive. Here, we have identified three metastable collective modes-the so-called single-q modes-in single layer TiSe_{2}, which originate from the unstable phonon eigenvectors at the zone boundary and break the threefold rotational symmetry. We show that polarized laser pulse is a unique and efficient tool to reconstruct the transient potential energy surface, so as to drive phase transitions between these states. By designing sequent layers with chiral stacking order, we propose a practical means to realize chiral charge density waves in 1T-TiSe_{2}. Further, the constructed chiral structure is predicted to exhibit circular dichroism as observed in recent experiments. These facts strongly indicate the chirality transfer from photons to the electron subsystem, meanwhile being strongly coupled to the lattice degree of freedom. Our work provides new insights into understanding and modulating chirality in quantum materials that we hope will spark further experimental investigation.

Circularly polarized attosecond light generation from OCS molecules irradiated by the combination of linear polarized infrared and orthogonal terahertz fields

Researching ultrafast dynamics and creating coherent light sources will both benefit significantly from the establishment of polarization control in high-order harmonic generation (HHG). By employing the time-dependent density functional theory method, we investigate HHG of carbonyl sulfide molecules using a combination of a linear polarized infrared (IR) laser and a weaker orthogonal Terahertz (THz) field. Our findings show that by adjusting the amplitude of the THz field, the movement scale of electrons in the THz direction can be tuned, thereby one can control the harmonic intensity in the IR laser direction. This method allows for the creation of near-circularly polarized attosecond pulses. Furthermore, the ellipticity of the attosecond pulse may be changed by modifying the carrier-envelope phase of the IR laser pulse.

Coherent ultrafast photoemission from a single quantized state of a one-dimensional emitter

Femtosecond laser-driven photoemission source provides an unprecedented femtosecond-resolved electron probe not only for atomic-scale ultrafast characterization but also for free-electron radiation sources. However, for conventional metallic electron source, intense lasers may induce a considerable broadening of emitting energy level, which results in large energy spread (>600 milli-electron volts) and thus limits the spatiotemporal resolution of electron probe. Here, we demonstrate the coherent ultrafast photoemission from a single quantized energy level of a carbon nanotube. Its one-dimensional body can provide a sharp quantized electronic excited state, while its zero-dimensional tip can provide a quantized energy level act as a narrow photoemission channel. Coherent resonant tunneling electron emission is evidenced by a negative differential resistance effect and a field-driven Stark splitting effect. The estimated energy spread is ~57 milli-electron volts, which suggests that the proposed carbon nanotube electron source may promote electron probe simultaneously with subangstrom spatial resolution and femtosecond temporal resolution.

Giant acceleration of polaron transport by ultrafast laser-induced coherent phonons

Polaron formation is ubiquitous in polarized materials, but severely hampers carrier transport for which effective controlling methods are urgently needed. Here, we show that laser-controlled coherent phonon excitation enables orders of magnitude enhancement of carrier mobility via accelerating polaron transport in a prototypical material, lithium peroxide (Li2O2). The selective excitation of specific phonon modes, whose vibrational pattern directly overlap with the polaronic lattice deformation, can remarkably reduce the energy barrier for polaron hopping. The strong nonadiabatic couplings between the electronic and ionic subsystem play a key role in triggering the migration of polaron, via promoting phonon-phonon scattering in q space within sub-picoseconds. These results extend our understanding of polaron transport dynamics to the nonequilibrium regime and allow for optoelectronic devices with ultrahigh on-off ratio and ultrafast responsibility, competitive with those of state-of-the-art devices fabricated based on free electron transport.

Coherent-Phonon-Driven Intervalley Scattering and Rabi Oscillation in Multivalley 2D Materials

Resolving the complete electron scattering dynamics mediated by coherent phonons is crucial for understanding electron-phonon couplings beyond equilibrium. Here we present a time-resolved theoretical investigation on strongly coupled ultrafast electron and phonon dynamics in monolayer WSe_{2}, with a focus on the intervalley scattering from the optically "bright" K state to "dark" Q state. We find that the strong coherent lattice vibration along the longitudinal acoustic phonon mode [LA(M)] can drastically promote K-to-Q transition on a timescale of ∼400  fs, comparable with previous experimental observation on thermal-phonon-mediated electron dynamics. Further, this coherent-phonon-driven intervalley scattering occurs in an unconventional steplike manner and further induces an electronic Rabi oscillation. By constructing a two-level model and quantitatively comparing with ab initio dynamic simulations, we uncover the critical role of nonadiabatic coupling effects. Finally, a new strategy is proposed to effectively tune the intervalley scattering rates by varying the coherent phonon amplitude, which could be realized via light-induced nonlinear phononics that we hope will spark experimental investigation.

Time-dependent equation-of-motion coupled-cluster simulations with a defective Hamiltonian

Simulations of laser-induced electron dynamics in a molecular system are performed using time-dependent (TD) equation-of-motion (EOM) coupled-cluster (CC) theory. The target system has been chosen to highlight potential shortcomings of truncated TD-EOM-CC methods [represented in this work by TD-EOM-CC with single and double excitations (TD-EOM-CCSD)], where unphysical spectroscopic features can emerge. Specifically, we explore driven resonant electronic excitations in magnesium fluoride in the proximity of an avoided crossing. Near the avoided crossing, the CCSD similarity-transformed Hamiltonian is defective, meaning that it has complex eigenvalues, and oscillator strengths may take on negative values. When an external field is applied to drive transitions to states exhibiting these traits, unphysical dynamics are observed. For example, the stationary states that make up the time-dependent state acquire populations that can be negative, exceed one, or even complex-valued.

Polarization Dependence of Laser-Induced Dynamics on Non-Flat Metal Surfaces: A Time-Dependent Density Functional Theory Approach

Real-time time-dependent density functional theory was used to study the laser-pulse-induced ion dynamics on metal surfaces featuring rows of atomic ridges. In contrast to atomically flat surfaces, the rows of atomic ridges induce anisotropy on the surface even in surface-parallel directions. This anisotropy causes the laser-induced ion dynamics to depend on the orientation of the laser polarization vector in the surface-parallel directions. This polarization dependence occurs for both copper (111) and aluminum (111) surfaces, indicating that the existence of localized d orbitals in the electronic system does not play a crucial role. The difference in kinetic energies between ions on the ridges and those on the planar surface reached a maximum when the laser polarization vector was perpendicular to the rows of ridges but parallel to the surface. A simple mechanism for the polarization dependence and some potential applications in laser processing are discussed.

A comprehensive investigation of the plasmonic-photocatalytic properties of gold nanoparticles for CO2 conversion to chemicals

Understanding the interactions between plasmonic gold (Au) nanoparticles and the adsorbate is essential for photocatalytic and plasmonic applications. However, it is often challenging to identify a specific reaction mechanism in the ground state and to explore the optical properties in the excited states because of the complicated pathways of carriers. In this study, photocatalytic reduction of carbon dioxide (CO₂) to C1 products (for example, CO and CH₄) on the Au(111) nanoparticle (NP) surface was studied based on reaction pathway analysis, adsorbate reactivity, and its ability to stabilize or deactivate the surface. The calculated reaction Gibbs free energies and activation barriers revealed that the first step in CO reduction via a direct hydrogen transfer mechanism on Au(111) is the formation of formyl (*CHO) instead of hydroxymethylidyne (*COH). Furthermore, the size enhanced and symmetry sensitive optical responses of cuboctahedral Au(111) NPs on localized surface plasmon resonance (LSPR) were investigated by using time-dependent DFT (TDDFT) calculations. Although near field enhancement around cuboctahedral Au(111) NPs is only weakly dependent on the morphology of NPs, it was observed that corner sites stabilize *C-species to drive the CO₂ reduction to CO. The density of active surface states interacting with the adsorbate states near the Fermi level gradually decreases from the (111) on-top site toward the corner site of the Au(111) NP-CO system, which strongly affects the molecule's binding on catalytic sites and, in particular, electronic excitation. Finally, the spatial distribution of the charge oscillations was determined as a guide for the fabrication of Au NPs with an optimal LSPR response.

Accurate Relativistic Real-Time Time-Dependent Density Functional Theory for Valence and Core Attosecond Transient Absorption Spectroscopy

First principles theoretical modeling of out-of-equilibrium processes observed in attosecond pump-probe transient absorption spectroscopy (TAS) triggering pure electron dynamics remains a challenging task, especially for heavy elements and/or core excitations containing fingerprints of scalar and spin-orbit relativistic effects. To address this, we formulate a methodology for simulating TAS within the relativistic real-time, time-dependent density functional theory (RT-TDDFT) framework, for both the valence and core energy regimes. Especially for TAS, full four-component (4c) RT simulations are feasible but computationally demanding. Therefore, in addition to the 4c approach, we also introduce the atomic mean-field exact two-component (amfX2C) Hamiltonian accounting for one- and two-electron picture-change corrections within RT-TDDFT. amfX2C preserves the accuracy of the parent 4c method at a fraction of its computational cost. Finally, we apply the methodology to study valence and near-L_2,3-edge TAS processes of experimentally relevant systems and provide additional physical insights using relativistic nonequilibrium response theory.

Efficient All-Electron Time-Dependent Density Functional Theory Calculations Using an Enriched Finite Element Basis

We present an efficient and systematically convergent approach to all-electron real-time time-dependent density functional theory (TDDFT) calculations using a mixed basis, termed as enriched finite element (EFE) basis. The EFE basis augments the classical finite element basis (CFE) with a compactly supported numerical atom-centered basis, obtained from atomic ground-state DFT calculations. Particularly, we orthogonalize the enrichment functions with respect to the classical finite element basis to ensure good conditioning of the resultant basis. We employ the second-order Magnus propagator in conjunction with an adaptive Krylov subspace method for efficient time evolution of the Kohn-Sham orbitals. We rely on a priori error estimates to guide our choice of an adaptive finite element mesh as well as the time step to be used in the TDDFT calculations. We observe close to optimal rates of convergence of the dipole moment with respect to spatial and temporal discretizations. Notably, we attain a 50-100 times speedup for the EFE basis over the CFE basis. We also demonstrate the efficacy of the EFE basis for both linear and nonlinear responses by studying the absorption spectra in sodium clusters, the linear to nonlinear response transition in the green fluorescence protein chromophore, and the higher harmonic generation in the magnesium dimer. Lastly, we attain good parallel scalability of our numerical implementation of the EFE basis for up to ∼1000 processors, using a benchmark system of a 50-atom sodium nanocluster.

Decoupled ultrafast electronic and structural phase transitions in photoexcited monoclinic VO2

Photoexcitation has emerged as an efficient way to trigger phase transitions in strongly correlated materials. There are great controversies about the atomistic mechanisms of structural phase transitions (SPTs) from monoclinic (M₁-) to rutile (R-) VO₂ and its association with electronic insulator-metal transitions (IMTs). Here, we illustrate the underlying atomistic processes and decoupling nature of photoinduced SPT and IMT in nonequilibrium states. The photoinduced SPT proceeds in the order of dilation of V-V pairs and increase of twisting angles after a small delay of ~40 fs. Dynamic simulations with hybrid functionals confirm the existence of isostructural IMT. The photoinduced SPT and IMT exhibit the same thresholds of electronic excitations, indicating similar fluence thresholds in experiments. The IMT is quasi-instantaneously (<10 fs) generated, while the SPT takes place with time a constant of 100 to 300 fs. These findings clarify some key controversies in the literature and provide insights into nonequilibrium phase transitions in correlated materials.

Laser-induced enhancement of vertical polarization in ferroelectric bilayer WTe2

Light-matter interaction is one of the key means to manipulate the structural and electronic properties of materials, especially in two-dimensional (2D) layered materials, which are optically accessible due to their atomic thickness. We propose that an ultrashort laser pulse could drastically enhance the ferroelectric polarization of bilayer WTe₂by our real-time time-dependent density functional theory simulations. It is noted that bilayer WTe₂is a 2D sliding ferroelectric material recently discovered whose vertical polarization can be controlled by a slight horizontal displacement. We demonstrate that interlayer sliding and compression are simultaneously achieved upon illumination of linearly polarized near-infrared laser pulse, leading to an ultrafast electric polarization enhancement by ∼230% within hundreds of femtosecond. Two major contributions have been identified: (a) the piezoelectric effect due to laser-induced interlayer compression, caused by interlayer charge transfer and dipole-dipole interaction; (b) the interlayer sliding along the opposite direction of ferroelectric switching, induced by inhomogeneous excited carrier distribution and specific electron-phonon couplings. This work provides new insights on controlling ferroelectricity of layered materials, which may extend to other van der Waals bilayers and even bulk materials.

Orbital-Free Methods for Plasmonics: Linear Response

Plasmonic systems, such as metal nanoparticles, are widely used in different application areas, going from biology to photovoltaics.The modeling of the optical response of such systems is of fundamental importance to analyze their behavior and to design new systems with required properties.When the characteristic sizes/distances reach a few nanometers, non-local and spill-out effects become relevant and conventional classical electrodynamics models are no more appropriate. Methods based on the Time-Dependent Density-Functional Theory (TD-DFT) represent the current reference for the description of quantum effects. However, TD-DFT is based on knowledge of all occupied orbitals whose calculation is computationally prohibitive to model large plasmonic systems of interest for applications.On the other hand, methods based on the Orbital-Free (OF) formulation of TD-DFT, can scale linearly with the system size.In this Review, OF methods ranging from semiclassical models to the quantum hydrodynamic theory, will be derived from the linear response TD-DFT, so that the key approximations and properties of each method can be clearly highlighted. The accuracy of the various approximations will be then validated for the linear optical properties of jellium nanoparticles, the most relevant model system in plasmonics. OF methods can describe the collective excitations in plasmonic systems with great accuracy andwithout system-tuned parameters. The accuracy on these methods depends only on the accuracy on the (universal) kinetic energy functional of the ground-state electronic density. Current approximations and future development directions will be indicated.

Traversing Double-Well Potential Energy Surfaces: Photoinduced Concurrent Intralayer and Interlayer Structural Transitions in XTe2 (X = Mo, W)

The microscopic arrangement of atoms and molecules is the determining factor in how materials behave and perform; i.e., the structure determines the property, a traditional paradigm in materials science. Photoexcitation-driven manipulation of the crystal structure and associated electronic properties in quantum materials provides opportunities for the exploration of exotic physics and practical applications; however, a generalized mechanism for such symmetry engineering is absent. Here, by ultrafast electron diffraction, structure factor calculation, and TDDFT-MD simulations, we report the photoinduced concurrent intralayer and interlayer structural transitions in the Td and 1T' phases of XTe₂ (X = Mo, W). We discuss the modification of multiple quantum electronic states associated with the intralayer and interlayer structural transitions, such as the topological band inversion and the higher-order topological state. The twin structures and the stacking faults in XTe₂ are also identified by ultrafast structural responses. The comprehensive study of the ultrafast structural response in XTe₂ suggests the traversal of all double-well potential energy surfaces (DWPES) by laser excitation, which is expected to be an intrinsic mechanism in the field of photoexcitation-driven global/local symmetry engineering and also a critical ingredient inducing the exotic properties in the non-equilibrium state in a large number of material systems.

Calibrating Out-of-Equilibrium Electron-Phonon Couplings in Photoexcited MoS2

Nonequilibrium electron-phonon coupling (EPC) serves as a dominant interaction in a multitude of transient processes, including photoinduced phase transitions, coherent phonon generation, and possible light-induced superconductivity. Here we use monolayer MoS₂ as a prototype to investigate the variation in electron-phonon couplings under laser excitation, on the basis of real-time time-dependent density functional theory simulations. Phonon softening, anisotropic modification of the deformation potential, and enhancement of EPC are observed, which are attributed to the reduced electronic screening and modulated potential energy surfaces by photoexcitation. Furthermore, by tracking the transient deformation potential and nonthermal electronic population, we can monitor the ultrafast time evolution of the energy exchange rate between electrons and phonons upon laser excitation. This work provides an effective strategy to investigate the nonequilibrium EPC and constructs a scaffold for understanding nonequilibrium states beyond the multitemperature models.

Fano Resonance and Incoherent Interlayer Excitons in Molecular van der Waals Heterostructures

Complex van der Waals heterostructures from layered molecular stacks are promising optoelectronic materials offering the means to efficient, modular charge separation and collection layers. The effect of stacking in the electrodynamics of such hybrid organic-inorganic two-dimensional materials remains largely unexplored, whereby molecular scale engineering could lead to advanced optical phenomena. For instance, tunable Fano engineering could make possible on-demand transparent conducting layers or photoactive elements, and passive cooling. We employ an adapted Gersten-Nitzan model and real time time-dependent density functional tight-binding to study the optoelectronics of self-assembled monolayers on graphene nanoribbons. We find Fano resonances that cause electromagnetic induced opacity and transparency and reveal an additional incoherent process leading to interlayer exciton formation with a characteristic charge transfer rate. These results showcase hybrid van der Waals heterostructures as paradigmatic 2D optoelectronic stacks, featuring tunable Fano optics and unconventional charge transfer channels.

Real-Time, Time-Dependent Density Functional Theory Study on Photoinduced Isomerizations of Azobenzene Under a Light Field

The trans to cis photoisomerization of azobenzene and its reverse (i.e., the cis to trans) processes are studied using real-time propagation time-dependent density functional theory combined with molecular dynamics for ions. We show that the wavelength of the applied laser may significantly affect the transition process. The simulations also show that the photon-excited electrons play essential roles in the isomerization processes, in which the hot electrons couple to phonon modes that drive the transitions.

Optical Control of Multistage Phase Transition via Phonon Coupling in MoTe_{2}

The temporal characters of laser-driven phase transition from 2H to 1T^{'} has been investigated in the prototype MoTe_{2} monolayer. This process is found to be induced by fundamental electron-phonon interactions, with an unexpected phonon excitation and coupling pathway closely related to the nonequilibrium relaxation of photoexcited electrons. The order-to-order phase transformation is dissected into three substages, involving energy and momentum scattering processes from optical (A_{1}^{'} and E^{'}) to acoustic phonon modes [LA(M)] in subpicosecond timescale. An intermediate metallic state along the nonadiabatic transition pathway is also identified. These results have profound implications on nonequilibrium phase engineering strategies.

Creation of a novel inverted charge density wave state

Charge density wave (CDW) order is an emergent quantum phase that is characterized by periodic lattice distortion and charge density modulation, often present near superconducting transitions. Here, we uncover a novel inverted CDW state by using a femtosecond laser to coherently reverse the star-of-David lattice distortion in 1T-TaSe₂. We track the signature of this novel CDW state using time- and angle-resolved photoemission spectroscopy and the time-dependent density functional theory to validate that it is associated with a unique lattice and charge arrangement never before realized. The dynamic electronic structure further reveals its novel properties that are characterized by an increased density of states near the Fermi level, high metallicity, and altered electron-phonon couplings. Our results demonstrate how ultrafast lasers can be used to create unique states in materials by manipulating charge-lattice orders and couplings.

Inq, a Modern GPU-Accelerated Computational Framework for (Time-Dependent) Density Functional Theory

We present inq, a new implementation of density functional theory (DFT) and time-dependent DFT (TDDFT) written from scratch to work on graphic processing units (GPUs). Besides GPU support, inq makes use of modern code design features and takes advantage of newly available hardware. By designing the code around algorithms, rather than against specific implementations and numerical libraries, we aim to provide a concise and modular code. The result is a fairly complete DFT/TDDFT implementation in roughly 12 000 lines of open-source C++ code representing a modular platform for community-driven application development on emerging high-performance computing architectures.

Plasmon-Induced Water Splitting on Ag-Alloyed Pt Single-Atom Catalysts

A promising route to realize solar-to-chemical energy conversion resorts to water splitting using plasmon photocatalysis. However, the ultrafast carrier dynamics and underlying mechanism in such processes has seldom been investigated, especially when the single-atom catalyst is introduced. Here, from the perspective of quantum dynamics at the atomic length scale and femtosecond time scale, we probe the carrier and structural dynamics of plasmon-assisted water splitting on an Ag-alloyed Pt single-atom catalyst, represented by the Ag₁₉Pt nanocluster. The substitution of an Ag atom by the Pt atom at the tip of the tetrahedron Ag₂₀ enhances the interaction between water and the nanoparticle. The excitation of localized surface plasmons in the Ag₁₉Pt cluster strengthens the charge separation and electron transfer upon illumination. These facts cooperatively turn on more than one charge transfer channels and give rise to enhanced charge transfer from the metal nanoparticle to the water molecule, resulting in rapid plasmon-induced water splitting. These results provide atomistic insights and guidelines for the design of efficient single-atom photocatalysts for plasmon-assisted water splitting.

All-electron real-time and imaginary-time time-dependent density functional theory within a numeric atom-centered basis function framework

Real-time time-dependent density functional theory (RT-TDDFT) is an attractive tool to model quantum dynamics by real-time propagation without the linear response approximation. Sharing the same technical framework of RT-TDDFT, imaginary-time time-dependent density functional theory (it-TDDFT) is a recently developed robust-convergence ground state method. Presented here are high-precision all-electron RT-TDDFT and it-TDDFT implementations within a numerical atom-centered orbital (NAO) basis function framework in the FHI-aims code. We discuss the theoretical background and technical choices in our implementation. First, RT-TDDFT results are validated against linear-response TDDFT results. Specifically, we analyze the NAO basis sets' convergence for Thiel's test set of small molecules and confirm the importance of the augmentation basis functions for adequate convergence. Adopting a velocity-gauge formalism, we next demonstrate applications for systems with periodic boundary conditions. Taking advantage of the all-electron full-potential implementation, we present applications for core level spectra. For it-TDDFT, we confirm that within the all-electron NAO formalism, it-TDDFT can successfully converge systems that are difficult to converge in the standard self-consistent field method. We finally benchmark our implementation for systems up to ∼500 atoms. The implementation exhibits almost linear weak and strong scaling behavior.

Attosecond charge migration following oxygen K-shell ionization in DNA bases and base pairs

Core ionization of DNA begins a cascade of events which could lead to cellular inactivation or death. The created core-hole following an impulse inner-shell ionization of molecules naturally decays in the auger timescale. We simulated charge migration (CM) phenomena following an impulsive core ionization of individual DNA bases at the oxygen K-edge which occurs before Auger decay of the oxygen. Our approach is based on real-time time dependent density functional theory (RT-TDDFT). It is shown that the pronounced hole fluctuation observed around bonds of the initial core-hole results in various valence orbital migrations. Also, the same photo-core-ionized dynamics is studied for the related base pairs. We investigate the role of base pairing and H-bonding interactions in the attosecond CM dynamics. In particular, the creation of a core-hole in the oxygen involved in H-bonding leads to an enhancement of charge migration relative to the respective single bases. Importantly, the hole oscillation of the adenine-thymine base pair upon creation of a core-hole at the oxygen, which does not contribute to the donor-acceptor interactions (not H-bonded), decreases compared to the single thymine base. Understanding the detailed dynamics of the localized core-hole initiating CM process would open the way for chemically controlling DNA damage/repair in the future.

Simulating electronic excitation and dynamics with real-time propagation approach to TDDFT within plane-wave pseudopotential formulation

We give a perspective on simulating electronic excitation and dynamics using the real-time propagation approach to time-dependent density functional theory (RT-TDDFT) in the plane-wave pseudopotential formulation. RT-TDDFT is implemented in various numerical formalisms in recent years, and its practical application often dictates the most appropriate implementation of the theory. We discuss recent developments and challenges, emphasizing numerical aspects of studying real systems. Several applications of RT-TDDFT simulation are discussed to highlight how the approach is used to study interesting electronic excitation and dynamics phenomena in recent years.

Direct treatment of interaction between laser-field and electrons for simulating laser processing of metals

Laser ablation is often simulated by the two-temperature model in which electrons are assumed to be thermalized by laser irradiation, while an explicit representation of interaction between laser-field and electrons is challenging but beneficial as being free from any adjustable parameters. Here, an ab initio method based on the time-dependent density functional theory (TDDFT) in which electron-ion dynamics under a laser field are numerically simulated is examined as a tool for simulating femtosecond laser processing of metals. Laser-induced volume expansion in surface normal directions of Cu(111) and Ni(111) surfaces are simulated by using repeating slab models. The amount of simulated volume expansion is compared between Cu(111) and Ni(111) slabs for the same laser pulse conditions, and the Ni slab is found to expand more than the Cu slab despite the smaller thermal expansion coefficient of Ni compared with Cu. The analyzed electronic excitation and lattice motion were compared to those in the two-temperature model. The threshold fluence to release surface Cu atom deduced from current TDDFT approach is found to be comparable to those of Cu ablation reported experimentally.

The 2021 ultrafast spectroscopic probes of condensed matter roadmap

In the 60 years since the invention of the laser, the scientific community has developed numerous fields of research based on these bright, coherent light sources, including the areas of imaging, spectroscopy, materials processing and communications. Ultrafast spectroscopy and imaging techniques are at the forefront of research into the light-matter interaction at the shortest times accessible to experiments, ranging from a few attoseconds to nanoseconds. Light pulses provide a crucial probe of the dynamical motion of charges, spins, and atoms on picosecond, femtosecond, and down to attosecond timescales, none of which are accessible even with the fastest electronic devices. Furthermore, strong light pulses can drive materials into unusual phases, with exotic properties. In this roadmap we describe the current state-of-the-art in experimental and theoretical studies of condensed matter using ultrafast probes. In each contribution, the authors also use their extensive knowledge to highlight challenges and predict future trends.

Ab initio nonadiabatic molecular dynamics of charge carriers in metal halide perovskites

Photoinduced nonequilibrium processes in nanoscale materials play key roles in photovoltaic and photocatalytic applications. This review summarizes recent theoretical investigations of excited state dynamics in metal halide perovskites (MHPs), carried out using a state-of-the-art methodology combining nonadiabatic molecular dynamics with real-time time-dependent density functional theory. The simulations allow one to study evolution of charge carriers at the ab initio level and in the time-domain, in direct connection with time-resolved spectroscopy experiments. Eliminating the need for the common approximations, such as harmonic phonons, a choice of the reaction coordinate, weak electron-phonon coupling, a particular kinetic mechanism, and perturbative calculation of rate constants, we model full-dimensional quantum dynamics of electrons coupled to semiclassical vibrations. We study realistic aspects of material composition and structure and their influence on various nonequilibrium processes, including nonradiative trapping and relaxation of charge carriers, hot carrier cooling and luminescence, Auger-type charge-charge scattering, multiple excitons generation and recombination, charge and energy transfer between donor and acceptor materials, and charge recombination inside individual materials and across donor/acceptor interfaces. These phenomena are illustrated with representative materials and interfaces. Focus is placed on response to external perturbations, formation of point defects and their passivation, mixed stoichiometries, dopants, grain boundaries, and interfaces of MHPs with charge transport layers, and quantum confinement. In addition to bulk materials, perovskite quantum dots and 2D perovskites with different layer and spacer cation structures, edge passivation, and dielectric screening are discussed. The atomistic insights into excited state dynamics under realistic conditions provide the fundamental understanding needed for design of advanced solar energy and optoelectronic devices.

Manipulating Weyl quasiparticles by orbital-selective photoexcitation in WTe2

Optical control of structural and electronic properties of Weyl semimetals allows development of switchable and dissipationless topological devices at the ultrafast scale. An unexpected orbital-selective photoexcitation in type-II Weyl material WTe2 is reported under linearly polarized light (LPL), inducing striking transitions among several topologically-distinct phases mediated by effective electron-phonon couplings. The symmetry features of atomic orbitals comprising the Weyl bands result in asymmetric electronic transitions near the Weyl points, and in turn a switchable interlayer shear motion with respect to linear light polarization, when a near-infrared laser pulse is applied. Consequently, not only annihilation of Weyl quasiparticle pairs, but also increasing separation of Weyl points can be achieved, complementing existing experimental observations. In this work, we provide a new perspective on manipulating the Weyl node singularity and coherent control of electron and lattice quantum dynamics simultaneously.

Dual Resonant Sum Frequency Generations from Two-Dimensional Materials

We propose dual resonant optical sum frequency generation (SFG), where the two most singular resonances could be selected, and report for the monolayer (1L-) WSe₂ when one (ω₁) of two excitation pulses is resonant to A exciton and their sum frequency (ω₁ + ω₂) to D exciton. The dual resonant SFG confirms that, under an irradiation of ω₁ and ω₂ pulses with the same fluence of ∼1.4 × 10¹⁰ W/m², its signal intensity could be enhanced about 20 times higher than the resonant SHG (i.e., 2ω₁ to the D excitonic absorption). Further, the dual resonant SFG intensity of 1L-WSe₂ is found to be 1 order of magnitude higher than the single resonant SFG intensity of 1L-WS₂ under the same condition of two-pulse irradiation. Finally, observations of the dual resonant SFG are thoroughly examined using real-time time-dependent density functional theory (rt-TDDFT), and the relevant nonlinear optical characteristics are scrutinized using the Greenwood-Kubo formalism.

How does the porphyrin-like vacancy affect the spectral properties of graphene quantum dots? A theoretical study

By processing graphene quantum dot, an ideal semiconductor material with suitable band gap and higher electron mobility can be obtained. Thus, it has a broad prospect in the application of photoelectric response materials. Here, a graphene defect with porphyrin-like structure is selected to achieve the controlable light absorption. The double five-membered-ring parallel vacancy are based on self-healing properties of popular graphene defects. Aimed to separate exciton and hole more effectively and achieve higher photoelectric conversion efficiency, the occupied orbital and unoccupied orbital of the quantum dot with objected defect structure is taken as orderly dispersion to form an obvious charge separation state under the demonstration with first principles calculation. Most importantly, a real time real space charge separation is calculated by time-dependent ab-initio quantum dynamics based on numerical atomic basis sets. The result shows the specific graphene defects can form an efficient pure graphene photoelectric response medium like porphyrin skeleton, and the vacancy will induce to adjust and control the specific wavelength of the response light and charge separated state manipulably with odevity of number of peripheral carbon rings by the calculation of ultrafast process.

Ultrafast charge ordering by self-amplified exciton-phonon dynamics in TiSe2

The origin of charge density waves (CDWs) in TiSe\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{2}$$\end{document}2 has long been debated, mainly due to the difficulties in identifying the timescales of the excitonic pairing and electron–phonon coupling (EPC). Without a time-resolved and microscopic mechanism, one has to assume simultaneous appearance of CDW and periodic lattice distortions (PLD). Here, we accomplish a complete separation of ultrafast exciton and PLD dynamics and unravel their interplay in our real-time time-dependent density functional theory simulations. We find that laser pulses knock off the exciton order and induce a homogeneous bonding–antibonding transition in the initial 20 fs, then the weakened electronic order triggers ionic movements antiparallel to the original PLD. The EPC comes into play after the initial 20 fs, and the two processes mutually amplify each other leading to a complete inversion of CDW ordering. The self-amplified dynamics reproduces the evolution of band structures in agreement with photoemission experiments. Hence we resolve the key processes in the initial dynamics of CDWs that help elucidate the underlying mechanism.

The physical origins of charge density waves in 1T-TiSe2 and their response to ultrafast excitation have long been a topic of theoretical and experimental debate. Here the authors present an ab initio theory that successfully captures the observed dynamics of charge density wave formation.

Extreme nonlinear strong-field photoemission from carbon nanotubes

Strong-field photoemission produces attosecond (10-18 s) electron pulses that are synchronized to the waveform of the incident light. This nonlinear photoemission lies at the heart of current attosecond technologies. Here we report a new nonlinear photoemission behaviour-the nonlinearity in strong-field regime sharply increases (approaching 40th power-law scaling), making use of sub-nanometric carbon nanotubes and 800 nm pulses. As a result, the carrier-envelope phase sensitive photoemission current shows a greatly improved modulation depth of up to 100% (with a total modulation current up to 2 nA). The calculations reveal that the behaviour is an interplay of valence band optical-field emission with charge interaction, and the nonlinear dynamics can be tunable by changing the bandgap of carbon nanotubes. The extreme nonlinear photoemission offers a new means of producing extreme temporal-spatial resolved electron pulses, and provides a new design philosophy for attosecond electronics and photonics.

Indirect but Efficient: Laser-Excited Electrons Can Drive Ultrafast Polarization Switching in Ferroelectric Materials

To enhance the efficiency of next-generation ferroelectric (FE) electronic devices, new techniques for controlling ferroelectric polarization switching are required. While most prior studies have attempted to induce polarization switching via the excitation of phonons, these experimental techniques required intricate and expensive terahertz sources and have not been completely successful. Here, we propose a new mechanism for rapidly and efficiently switching the FE polarization via laser-tuning of the underlying dynamical potential energy surface. Using time-dependent density functional calculations, we observe an ultrafast switching of the FE polarization in BaTiO₃ within 200 fs. A laser pulse can induce a charge density redistribution that reduces the original FE charge order. This excitation results in both desirable and highly directional ionic forces that are always opposite to the original FE displacements. Our new mechanism enables the reversible switching of the FE polarization with optical pulses that can be produced from existing 800 nm experimental laser sources.