Real-time Time-dependent Density Functional Theory Research Articles

Electrical conductivity of warm dense hydrogen from Ohm’s law and time-dependent density functional theory

Abstract Understanding the electrical conductivity of warm dense hydrogen is critical for both fundamental physics and applications in planetary science and inertial confinement fusion. We demonstrate how to calculate the electrical conductivity using the continuum form of Ohm’s law, with the current density obtained from real-time time-dependent density functional theory. This approach simulates the dynamic response of hydrogen under warm dense matter conditions, with temperatures around 30,000 K and mass densities ranging from 0.02 to 0.98 g/cm³. We systematically address finite-size errors in real-time time-dependent density functional theory, demonstrating that our calculations are both numerically feasible and reliable. Our results show good agreement with other approaches, highlighting the effectiveness of this method for modeling electronic transport properties from ambient to extreme conditions.

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.

Ultrafast structural transition and electron-phonon/phonon–phonon coupling in antimony revealed by nonadiabatic molecular dynamics

Real-time time-dependent density-functional theory molecular dynamics (rt-TDDFT-MD) reveals the nonadiabatic dynamics of the ultrafast photoinduced structural transition in a typical phase-change material antimony (Sb) with Peierls distortion (PD). As the excitation intensity increases from 3.54% to 5.00%, three distinct structural transition behaviors within 1 ps are observed: no PD flipping, nonvolatile-like PD flipping, and nonstop back-and-forward PD flipping. Analyses on electron-phonon and phonon-phonon couplings indicate that the excitation-activated coherent A1gphonon mode by electron-phonon coupling drives the structural transition within several hundred femtoseconds. Then, the energy of coherent motions are transformed into that of random thermal motions via phonon-phonon coupling, which prevents the A1g-mode-like coherent structure oscillations. The electron-phonon coupling and coherent motions will be enhanced with increasing the excitation intensity. Therefore, a moderate excitation intensity that can balance the coherent and decoherent thermal movements will result in a nonvolatile-like PD flipping. These findings illustrate important roles of nonadiabatic electron-phonon/phonon-phonon couplings in the ultrafast laser-induced structural transitions in materials with PD, offering insights for manipulating their structures and properties by light.

Optical spectra of silver clusters and nanoparticles from 4 to 923 atoms from the TDDFT+U method

The localized surface-plasmon resonances of coinage-metal clusters and nanoparticles enable many applications, the conception and necessary optimization of which require precise theoretical description and understanding. However, for the size range from few-atom clusters through nanoparticles of a few nanometers, where quantum effects and atomistic structure play a significant role, none of the methods employed previously has been able to provide high-quality spectra for all sizes. The main problem is the description of the filled shells of d electrons which influence the optical response decisively. We show that the DFT+U method, employed with real-time time-dependent density-functional theory calculations (RT-TDDFT), provides spectra in good agreement with experiment for silver clusters ranging from 4 to 923 atoms, the latter representing a nanoparticle of 3 nm. Both the electron-hole-type discrete spectra of the smallest clusters and the broad plasmon resonances of the larger sizes are obtained. All calculations use the value of the effective U parameter that provides good results in bulk silver. The agreement with experiment for all sizes shows that the U parameter is surprisingly transferable. Our results open the pathway for calculations of many practically relevant systems including clusters coupled to bio-molecules or to other nano-objects.

Direct Hydrogen Production Promoted by Laser-Induced Water Plasma.

Hydrogen, as a clean energy carrier, plays an important role in addressing the current energy and environmental crisis. However, conventional hydrogen production technologies require extreme reaction conditions, such as high temperature, high pressure, and catalysts. Herein, we study the microscopic mechanism of laser-induced water plasma and subsequent H2 production with real-time time-dependent density functional theory simulations and ab initio molecular dynamics simulations. The results demonstrate that intense laser excites liquid water to generate nonequilibrium plasma in a warm-dense state, which constitutes a superior reaction environment. Subsequent annealing leads to the recombination of energetic reactive particles to generate H2, O2, and H2O2 molecules. Annealing rate and laser wavelength are shown to modulate the product ratio, and the energy conversion efficiency can reach ∼9.2% with an annealing rate of 1.0 K/fs. This work reveals the nonequilibrium atomistic mechanisms of hydrogen production from laser-induced water plasma and shows far-reaching implications for the design of optically controllable hydrogen technology.

Optical force and torque in near-field excitation of C3H6: A first-principles study using RT-TDDFT.

Optical trapping is an effective tool for manipulating micrometer-sized particles, although its application to nanometer-sized particles remains difficult. The field of optical trapping has advanced significantly, incorporating more advanced techniques such as plasmonic structures. However, single-molecule trapping remains a challenge. To achieve a deeper understanding of optical forces acting on molecular systems, a first-principles approach to analyze the optical force on molecules interacting with a plasmonic field is crucial. In our study, the optical force and torque induced by the near-field excitation of C3H6 were investigated using real-time time-dependent density functional theory calculations on real-space grids. The near field from the scanning tunneling probe was adopted as the excitation source for the molecule. The optical force was calculated using the polarization charges induced in the molecule based on Lorentz force. While the optical force and torque calculated as functions of the light energy were in moderate agreement with the oscillator strengths obtained from the far-field excitation of C3H6, a closer correspondence was achieved with the power spectrum of the induced dipole moment using near-field excitation. Time-domain analysis of the optical force suggests that the simultaneous excitation of multiple excited states generally weakens the force because of mismatches between the directions of the induced polarization and the electric field. This study revealed a subtle damping mechanism for the optical force arising from intrinsic electronic states and the influence of beating.

Origins of Electromagnetic Radiation from Spintronic Terahertz Emitters: A Time-Dependent Density Functional Theory plus Jefimenko Equations Approach.

Microscopic origins of charge currents and electromagnetic (EM) radiation generated by them in spintronic THz emitters-such as, femtosecond laser pulse-driven single magnetic layer or its heterostructures with a nonmagnetic layer hosting strong spin-orbit coupling (SOC)-remain poorly understood despite nearly three decades since the discovery of ultrafast demagnetization. We introduce a first-principles method to compute these quantities, where the dynamics of charge and current densities is obtained from real-time time-dependent density functional theory, which are then fed into the Jefimenko equations for properly retarded electric and magnetic field solutions of the Maxwell equations. By Fourier transforming different time-dependent terms in the Jefimenko equations, we unravel that in the 0.1-30THz range the electric field of far-field EM radiation by the Ni layer, chosen as an example, is dominated by charge current pumped by demagnetization, while often invoked magnetic dipole radiation from the time-dependent magnetization of a single magnetic layer is a negligible effect. Such an effect of charge current pumping by a time-dependent quantum system, whose magnetization is shrinking while its vector does not rotate, does not require any spin-to-charge conversion via SOC effects. In the Ni/Pt bilayer, EM radiation remains dominated by the charge current within the Ni layer, whose magnitude is larger than in the case of a single Ni layer due to faster demagnetization, while often invoked spin-to-charge conversion within the Pt layer provides an additional but smaller contribution. By using the Poynting vector and its flux, we also quantify the efficiency of conversion of light into emitted EM radiation, and its angular distribution.

Controlling Ultrafast Magnetization Dynamics via Coherent Phonon Excitation in a Ferromagnet Monolayer.

Exploring ultrafast magnetization control in 2D magnets via laser pulses is established, yet the interplay between spin dynamics and the lattice remains underexplored. Utilizing real-time time-dependent density functional theory (rt-TDDFT) coupled with Ehrenfest dynamics and nonadiabatic molecular dynamics (NAMD) simulations, we systematically investigate the laser-induced spin-nuclei dynamics with pre-excited A1g and E2g coherent phonons in the 2D ferromagnet Fe3GeTe2 (FGT) monolayer. Selective pre-excitation of coherent phonons under ultrafast laser irradiation significantly alters the local spin moment of FGT, consequently inducing additional spin loss attributed to the nuclear motion-induced asymmetric interatomic charge transfer. Excited spin-resolved charge undergoes a bidirectional spin-flip between spin-down and spin-up states, characterized by a subtle change in the spin moment within approximately 100 fs, followed by unidirectional spin-flip, which will further contribute to the spin moment loss of FGT within tens of picoseconds. Our results shed light on the coupling of coherent phonons with magnetization dynamics in 2D limit.

Electronic structure, absorption spectra and oxidation dynamics in polyynes and dicyanopolyynes.

The advent of femtosecond to attosecond experimental tools has made now possible to study such ultrafast carrier dynamics, e.g., the spatial and temporal charge density evolution, after an initial oxidation or reduction in molecules, candidates for atomic wires like polyynes and dicyanopolyynes. Here, we study the electronic structure and hole transfer in symmetric molecules containing carbon, nitrogen and hydrogen, the first members in the series of polyynic carbynes and dicyanopolyynes, using methods based on density functional theory (DFT): constrained DFT (CDFT), time-dependent DFT (TDDFT) and real-time TDDFT (RT-TDDFT), with Löwdin population analysis, comparing many levels of theory and obtaining convergence of the results. For the same purposes, we develop a tight binding (TB) variant using all valence orbitals of all atoms. This TB variant is applied here in linear molecules, but it is also adequate for electronic structure, charge transfer and charge transport of non-linear molecules and clusters of molecules. We calculate the electronic structure, the time-dependent dipole moment and the probabilities of finding the hole at each site, their mean over time values, the mean transfer rates from the oxidation site to other sites and the frequency content (using charge as well as dipole moment oscillations). We take into account zero-point motion. The initial conditions for RT-TDDFT are obtained by CDFT. For TB, we explore different initial conditions: we place the hole at a particular orbital or distribute it among a number of orbitals; it is also possible to include phase differences between orbitals. Finally, we compare with available experimental data.

Improving carrier separation in ZnIn2S4 to boost photocatalytic degradation of metronidazole based on machine learning prediction, experimental verification and theory calculation

Carrier separation efficiency significantly impacts the photocatalyst performance for wastewater purification. However, there is no established theory to accurately guide carrier separation efficiency regulation. This study used machine learning to screen dopants to change the band gap of ZnIn2S4 (ZIS), thus precisely regulating the ZIS carrier separation efficiency. Nitrogen with 0.5 % of doping amount showed the best regulation effect. Metronidazole degradation kinetic rate with nitrogen doped ZIS (N/ZIS) was 1.47 times higher than that with ZIS. Real-time time-dependent density functional theory (rt-TDDFT) was employed to assess carrier separation at the orbital level. Experiments and rt-TDDFT confirmed that N/ZIS possessed better carrier separation efficiency than that of ZIS. An “in-situ” potential research strategy was established to determine carrier acceptors, which can compare the orbital potential of molecules and catalysts in the adsorption state. This approach revealed that O2 served as the photogenerated electron acceptor, while OH– functioned as the hole acceptor. H2O was not actively involved in carrier separation. O2 anti-bonding orbital activation effect by the catalysts was studied to unveil the influence mechanism of carrier separation efficiency on the photocatalyst activity. This study provides a new insight into carrier dynamics to design efficient photocatalysts.

QRCODE: Massively parallelized real-time time-dependent density functional theory for periodic systems

We present a new software module, QRCODE (Quantum Research for Calculating Optically Driven Excitations), for massively parallelized real-time time-dependent density functional theory (RT-TDDFT) calculations of periodic systems in the open-source Qbox software package. Our approach utilizes a custom implementation of a fast Fourier transformation scheme that significantly reduces inter-node message passing interface (MPI) communication of the major computational kernel and shows impressive scaling up to 16,344 CPU cores. In addition to improving computational performance, QRCODE contains a suite of various time propagators for accurate RT-TDDFT calculations. As benchmark applications of QRCODE, we calculate the current density and optical absorption spectra of hexagonal boron nitride (h-BN) and photo-driven reaction dynamics of the ozone-oxygen reaction. We also calculate the second and higher harmonic generation of monolayer and multi-layer boron nitride structures as examples of large material systems. Our optimized implementation of RT-TDDFT in QRCODE enables large-scale calculations of real-time electron dynamics of chemical and material systems with enhanced computational performance and impressive scaling across several thousand CPU cores.

Deeper-band electron contributions to stopping power of silicon for low-energy ions.

This study provides accurate results for the electronic stopping cross sections of H, He, N, and Ne in silicon in low to intermediate energy ranges using various non-perturbative theoretical methods, including real-time time-dependent density functional theory, transport cross section, and induced-density approach. Recent experimental findings [Ntemou et al., Phys. Rev. B 107, 155145 (2023)] revealed discrepancies between the estimates of density functional theory and the observed values. We show that these discrepancies vanish by considering the nonuniform electron density of the deeper silicon bands for ion velocities approaching zero (v → 0). This indicates that mechanisms such as "elevator" and "promotion," which can dynamically excite deeper-band electrons, are active, enabling a localized free-electron gas to emulate ion energy loss, as pointed out by Lim et al. [Phys. Rev. Lett. 116, 043201 (2016)]. The observation and the description of a velocity-proportionality breakdown in electronic stopping cross sections at very low velocities are considered to be a signature of the contributions of deeper-band electrons.

Near-field induced local excitation dynamics of Na10 and Na10-N2 from real-time TDDFT.

Electron dynamics of the Na10 chain and the Na10-N2 complex locally excited by an atomistic optical near-field are investigated using real-time time-dependent density functional theory calculations on real-space grids. Ultrafast laser pulses were used to simulate the near-field excitation under on- and off-resonance conditions. Off-resonance excitation did not lead to the propagation of the excitation through the Na10 chain. In contrast, under the resonance conditions, the excited state is delocalized over the entire Na chain. Analysis of the local dipole moment of each atom in Na10 indicates that this behavior is consistent with the transition density. Adding an N2 molecule to the opposite end of the local excitation region results in energy transfer via the Na10 chain. The energy transfer efficiency of the N2 molecule is well correlated with the absorption spectrum of Na10. The present study paves the way for realizing remote excitation and photonic devices at the atomic scale.

Modeling Intermolecular Coulombic Decay with Non-Hermitian Real-Time Time-Dependent Density Functional Theory.

In this work, we investigate the capability of using real-time time-dependent density functional theory (RT-TDDFT) in conjunction with a complex absorbing potential (CAP) to simulate the intermolecular Coulombic decay (ICD) processes following the ionization of an inner-valence electron. We examine the ICD dynamics in a series of noncovalent bonded dimer systems, including hydrogen-bonded and purely van der Waals (VdW)-bonded systems. In comparison to previous work, we show that RT-TDDFT simulations with a CAP correctly capture the ICD phenomenon in systems exhibiting a stronger binding energy. The calculated time scales for ICD of the studied systems are in the range of 5-50 fs, in agreement with previous studies. However, there is a breakdown in the accuracy of the methodology for the pure VdW-bonded systems. Overall, the presented RT-TDDFT/CAP methodology provides a powerful tool for differentiating between competing electronic relaxation pathways following inner-valence or core ionization without necessitating any a priori assumptions.

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.

Reduction reactions at the interface between CdS quantum dot and Z-type ligands driven by electron injection in the electroluminescent processes.

The efficient and stable electroluminescence of quantum dots (QDs) is of great importance in their applications in new display technologies. The short service life of blue QDs, however, hinders their development and commercialization. Different mechanisms have been proposed for the destabilization of QDs in electroluminescent processes. Based on real-time time-dependent density functional theory studies on the QD models covered by Z-type ligands (XAc2, X = Cd, Zn, Mg), the structural evolution is simulated to reveal the mechanism of the reduction reactions induced by electron injection. Our simulations reproduce the experimental observations that the reduction reactions occur at the QD-ligand interface, and the reduced Cd atom is almost in a zero valence state. However, different sites are predicted for the reactions in which the surface metal atom of the QD instead of the metal atom in the ligands is reduced. As a result, one of the arms of the chelate ligand leaves the QD, which tends to cause damage to its electroluminescent performance. Our findings contribute to a mechanistic understanding of the reduction reactions that occurred at the QD-ligand interface.

Ultrafast charge carrier dynamics of methylammonium lead iodide from first principles

Methylammonium lead iodide (MAPbI3) has been a major focus of photovoltaic research for the last decade. The unique interplay between the structural and electronic properties of this material contributes to its exciting optical properties especially under the action of an ultrafast laser pulse. First-principles methods like real-time time-dependent density functional theory (RT-TDDFT) enable performing corresponding simulations without the aid of empirical parameters: the gained knowledge can be applied to future studies of other complex materials. In this work, we investigate the ultrafast charge-carrier dynamics and the nonlinear optical response of MAPbI3 excited by a resonant pulse above the gap. First, we examine the electronic and optical properties in the static regime. Next, we impinge the system with a femtosecond field of varying intensity and follow the evolution of the photoexcited carrier density. A pronounced intensity-dependent response is observed, manifested by high-harmonic generation and nonlinear trends in the number of excited electrons and excitation energy. Our results provide relevant indications about the behavior of MAPbI3 under strong and coherent radiation and confirm that RT-TDDFT is a viable tool to simulate the photo-induced dynamics of complex materials from first principles.

A theoretical study on the electronic and optical properties of M-TiO2/ZnO (M=Li, Na, K, Fe) toward application in photocatalysis

This study employs the tight-binding density functional method GFN1-xTB to investigate the structural and electronic properties of TiO2-ZnO and TiO2-ZnO composite materials modified by Li, Na, K, and Fe metals. Computational analyses reveal the formation of weak covalent bonds between the TiO2 and ZnO components within the composite system. Doping TiO2-ZnO with metals induces significant alterations in the electronic structure, particularly in terms of ionization energy and global electrophilic index. The UV-VIS spectra are calculated using real-time time-dependent density functional theory with xTB Hamiltonians. The obtained results demonstrate minimal impact of metal presence on the absorption spectrum of TiO2-ZnO, but significant influence on the recombination potential of photogenerated charge carriers. It is suggested that Fe/TiO2-ZnO and Li/TiO2-ZnO will demonstrate higher photocatalytic activity than the pristine TiO2-ZnO material.

Real-Time Simulation of Ultrafast Electronic Dynamics of Nanoscale Systems Involving an Organic Molecule and a Nanoparticle Dimer.

Understanding and predicting the behavior of nanomaterials composed of plasmons interacting with quantum emitters at ultrafast timescales is crucial for the better manipulation of light at the nanoscale and advancing technologies like ultrafast communication and computing. Here we perform a simulation of the "real-time" electronic dynamics of a coupled molecule-metal nanoparticle dimer interacting with an ultrashort resonant laser pulse by combining the real-time time-dependent density functional theory (RT-TDDFT) approach with the time-domain frequency-dependent fluctuating charge (TD-ωFQ) model, an atomistic electromagnetic (AEM) model for the dynamic plasmonic response of nanoparticles. It is shown that the induced dipoles evolve from an exponential decay pattern to a beat pattern with an increase in coupling strength, which is altered by changing the molecular orientation relative to the dimer axis. It is further shown that in the strong coupling regime, both the excited molecule and the plasmon relax rapidly due to the molecule-plasmon interaction, and the efficient coherent energy exchange between the interacting molecule and plasmon modes occurs on a femtosecond (fs) timescale. This work provides guidance on manipulating light-matter interaction and studying molecular plasmonics at extremely fast timescales.

Propagated (fragment) Pipek-Mezey Wannier functions in real-time time-dependent density functional theory.

Localization procedures are an important tool for analysis of complex systems in quantum chemistry, since canonical molecular orbitals are delocalized and can, therefore, be difficult to align with chemical intuition and obscure information at the local level of the system. This especially applies to calculations obeying periodic boundary conditions. The most commonly used approach to localization is Foster-Boys Wannier functions, which use a unitary transformation to jointly minimize the second moment of the orbitals. This procedure has proven to be robust and fast but has a side effect of often mixing σ- and π-type orbitals. σ/π-separation is achieved by the Pipek-Mezey Wannier function (PMWF) approach [Lehtola and Jónsson, J. Chem. Theory Comput. 10, 642 (2014) and Jónsson et al., J. Chem. Theory Comput. 13, 460 (2017)], which defines the spread functional in terms of partial charges instead. We have implemented a PMWF algorithm in the CP2K software package using the Cardoso-Souloumiac algorithm to enable their application to real-time time-dependent density functional theory. The method is demonstrated on stacked CO2 molecules, linear acetylenic carbon, boron and nitrogen co-doped graphene, and nitrogen-vacancy doped diamond. Finally, we discuss its computational scaling and recent efforts to improve it with fragment approaches.