Real-Time Exciton Dynamics with Time-Dependent Density-Functional Theory

Stability of the Long-Range Corrected Exchange-Correlation Functional and the Proca Procedural Functional in Time-Dependent Density-Functional Theory

Excitonic effects in the optical absorption spectra of solids can be described with time-dependent density-functional theory (TDDFT) in the linear-response regime, using a simple class of approximate, long-range corrected (LRC) exchange-correlation functionals. It was recently demonstrated that the LRC approximation can also be employed in real-time TDDFT to describe exciton dynamics. Here, we investigate the numerical stability of the time-dependent LRC approach using a two-dimensional model solid. It is found that the time-dependent Kohn-Sham equation with an LRC vector potential becomes more and more prone to instabilities for increasing exciton binding energies. The origin of these instabilities is traced back to time-averaged violations of the zero-force theorem, which leads to a simple and robust numerical stabilization scheme. This explains and justifies a recently proposed method by Dewhurst, J. K. [ Phys. Rev. B 2025, 111, L060302] to stabilize the LRC vector potential, known as the Proca procedural functional.

Machine-Learning Electron Dynamics with Moment Propagation Theory: Application to Optical Absorption Spectrum Computation Using Real-Time TDDFT

We present an application of our new theoretical formulation of quantum dynamics, moment propagation theory (MPT) (Boyer et al., J. Chem. Phys. 160, 064113 (2024)), for employing machine-learning techniques to simulate the quantum dynamics of electrons. In particular, we use real-time time-dependent density functional theory (RT-TDDFT) simulation in the gauge of the maximally localized Wannier functions (MLWFs) for training the MPT equation of motion. Spatially localized time-dependent MLWFs provide a concise representation that is particularly convenient for the MPT expressed in terms of increasing orders of moments. The equation of motion for these moments can be integrated in time, while the analytical expressions are quite involved. In this work, machine-learning techniques were used to train the second-order time derivatives of the moments using first-principles data from the RT-TDDFT simulation, and this MPT enabled us to perform electron dynamics efficiently. The application to computing optical absorption spectrum for various systems was demonstrated as a proof-of-principles example of this approach. In addition to isolated molecules (water, benzene, and ethene), condensed matter systems (liquid water and crystalline silicon) were studied, and we also explored how the principle of the nearsightedness of electrons can be employed in this context.

Real-Time Dyson-Expansion Scheme: Efficient Inclusion of Dynamical Correlations in Nonequilibrium Spectral Properties

Time-resolved photoemission spectroscopy is the key technique to probe the real-time nonequilibrium dynamics of electronic states. Theoretical predictions of the time dependent spectral function for realistic systems is however, a challenge. Employing the Kadanoff-Baym equations to find this quantity results in a cubic scaling in the total number of time steps, quickly becoming prohibitive and often fail quantitatively and even qualitatively. In comparison, mean-field methods have more favorable numerical scaling both in the number of time steps and in the complexity associated with the cost of evolving for a single time step, however they miss key spectral properties such as emergent spectral features. Here we present a scheme that allows for the inclusion of dynamical correlations to the spectral function while maintaining the same scaling in the number of time steps as for mean-field approaches, while capturing the emergent physics. Further, the scheme can be efficiently implemented on top of equilibrium real-time many-body perturbation theory schemes and codes. We see excellent agreement with exact results for test systems. Furthermore, we exemplify the method on a periodic system and demonstrate clear evidence that our proposed scheme produces complex spectral features including excitonic band replicas, features that are not observed using static mean-field approaches.

High-Temperature Excitonic Condensation in 2D Lattice

Exploration of high-temperature bosonic condensation is of significant importance for the fundamental many-body physics and applications in nanodevices, which, however, remains a huge challenge. Here, in combination of many-body perturbation theory and first-principles calculations, a new-type spatially indirect exciton can be optically generated in two-dimensional (2D) Bi2S2Te because of its unique structure feature. In particular, the spin-singlet spatially indirect excitons in Bi2S2Te monolayer are dipole/parity allowed and reveal befitting characteristics for excitonic condensation, such as small effective mass and satisfied dilute limitation. Based on the layered Bi2S2Te, the possibility of the high-temperature excitonic Bose-Einstein condensation (BEC) and superfluid state in two dimensions, which goes beyond the current paradigms in both experiment and theory, are proved. It should be highlighted that record-high phase transition temperatures of 289.7 and 72.4 K can be theoretically predicted for the excitonic BEC and superfluidity in the atomic thin Bi2S2Te, respectively. It therefore can be confirmed that Bi2S2Te featuring bound bosonic states is a fascinating 2D platform for exploring the high-temperature excitonic condensation and applications in such as quantum computing and dissipationless nanodevices.

Reformulation of Time-Dependent Density Functional Theory for Nonperturbative Dynamics: The Rabi Oscillation Problem Resolved

Rabi oscillations have long been thought to be out of reach in simulations using time-dependent density functional theory (TDDFT), a prominent symptom of the failure of the adiabatic approximation for nonperturbative dynamics. We present a reformulation of TDDFT which requires response quantities only, thus enabling an adiabatic approximation to predict such dynamics accurately because the functional is evaluated on a density close to the ground state, instead of on the fully nonperturbative density. Our reformulation applies to any real-time dynamics, redeeming TDDFT far from equilibrium. Examples of a resonantly-driven local excitation in a model He atom, and charge-transfer in the LiCN molecule are given.

Subcycle dynamics of excitons under strong laser fields

Excitons play a key role in the linear optical response of two-dimensional (2D) materials. However, their role in the nonlinear response to intense, nonresonant, low-frequency light is often overlooked as strong fields are expected to tear the electron-hole pair apart. Using high-harmonic generation as a spectroscopic tool, we theoretically study their formation and role in the nonlinear optical response. We show that the excitonic contribution is prominent and that excitons remain stable even when the driving laser field surpasses the strength of the Coulomb field binding the electron-hole pair. We demonstrate a parallel between the behavior of strongly laser-driven excitons in 2D solids and strongly driven Rydberg states in atoms, including the mechanisms of their formation and stability. Last, we show how the excitonic contribution can be singled out by encapsulating the 2D material in a dielectric, tuning the excitonic energy and its contribution to the high-harmonic spectrum.

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.

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.

Laser-induced ultrafast spin-transfer processes in non-linear zigzag carbon chain systems

We combine the high-level quantum chemistry theory CCSD and EOM-CCSD together with local and global Λ processes to investigate the details of the laser-induced ultrafast spin manipulation scenarios in non-linear zigzag carbon chain systems Ni2@C32H32 and Ni2@C36H36. The spin density distribution, which is calculated on each many-body state using a Mulliken population analysis, fulfills the requirements to accomplish the spin dynamics processes. Various spin-flip and spin-transfer scenarios are accomplished. All the spin-dynamics processes can be achieved within subpicosecond times. Under the influence of a magnetic field, we find that the spin-transfer scenarios are preserved, while the local spin-flip scenario on a Ni atom can be significantly inhibited depending on the strength of the magnetic field. The impact of the propagation direction of the laser pulse on the spin dynamics processes by varying their polar and azimuthal angles in spherical coordinates is investigated. Additionally, we find that double laser pulses successfully induce the spin-transfer processes. Our outcomes underline the significant potential of carbon chain systems as building blocks for developing future all-optical integrated logic processing units.

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.

Computational Approaches for Organic Semiconductors: From Chemical and Physical Understanding to Predicting New Materials

While a complete understanding of organic semiconductor (OSC) design principles remains elusive, computational methods─ranging from techniques based in classical and quantum mechanics to more recent data-enabled models─can complement experimental observations and provide deep physicochemical insights into OSC structure-processing-property relationships, offering new capabilities for in silico OSC discovery and design. In this Review, we trace the evolution of these computational methods and their application to OSCs, beginning with early quantum-chemical methods to investigate resonance in benzene and building to recent machine-learning (ML) techniques and their application to ever more sophisticated OSC scientific and engineering challenges. Along the way, we highlight the limitations of the methods and how sophisticated physical and mathematical frameworks have been created to overcome those limitations. We illustrate applications of these methods to a range of specific challenges in OSCs derived from π-conjugated polymers and molecules, including predicting charge-carrier transport, modeling chain conformations and bulk morphology, estimating thermomechanical properties, and describing phonons and thermal transport, to name a few. Through these examples, we demonstrate how advances in computational methods accelerate the deployment of OSCsin wide-ranging technologies, such as organic photovoltaics (OPVs), organic light-emitting diodes (OLEDs), organic thermoelectrics, organic batteries, and organic (bio)sensors. We conclude by providing an outlook for the future development of computational techniques to discover and assess the properties of high-performing OSCs with greater accuracy.

Implementation of real-time TDDFT for periodic systems in the open-source PySCF software package

We present a new implementation of real-time time-dependent density functional theory (RT-TDDFT) for calculating excited-state dynamics of periodic systems in the open-source Python-based PySCF software package. Our implementation uses Gaussian basis functions in a velocity gauge formalism and can be applied to periodic surfaces, condensed-phase, and molecular systems. As representative benchmark applications, we present optical absorption calculations of various molecular and bulk systems and a real-time simulation of field-induced dynamics of a (ZnO)₄ molecular cluster on a periodic graphene sheet. We present representative calculations on optical response of solids to infinitesimal external fields as well as real-time charge-transfer dynamics induced by strong pulsed laser fields. Due to the widespread use of the Python language, our RT-TDDFT implementation can be easily modified and provides a new capability in the PySCF code for real-time excited-state calculations of chemical and material systems.

Theoretical Approach for Electron Dynamics and Ultrafast Spectroscopy (EDUS)

In this manuscript, we present a theoretical framework and its numerical implementation to simulate the out-of-equilibrium electron dynamics induced by the interaction of ultrashort laser pulses in condensed-matter systems. Our approach is based on evolving in real time the density matrix of the system in reciprocal space. It considers excitonic and nonperturbative light-matter interactions. We show some relevant examples that illustrate the efficiency and flexibility of the approach to describe realistic ultrafast spectroscopy experiments. Our approach is suitable for modeling the promising and emerging ultrafast studies at the attosecond time scale that aim at capturing the electron dynamics and the dynamical electron-electron correlations via X-ray absorption spectroscopy.

Analyzing Excitation-Energy Transfer Based on the Time-Dependent Density Functional Theory in Real Time

Excitation-energy transfer is a key step in processes such as photosynthesis that convert light into other forms of energy. Time-dependent density functional theory (DFT) in real time is ideal for the first-principles simulation of such processes due to its computational efficiency. We here demonstrate how real-time DFT can be used for analyzing excitation-energy transfer from first-principles. We discuss several measures of energy transfer that are based solely on the time-dependent density, are well founded in the DFT framework, allow for intuitive understanding and visualization, and reproduce important limiting cases of an analytical model. We demonstrate their usefulness in calculations for model systems, both with static nuclei and in the context of DFT-based Ehrenfest dynamics.

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.