Ab Initio Simulation of Attosecond Transient Absorption Spectroscopy in Two-Dimensional Materials

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.

Laser-Induced Electronic and Vibronic Dynamics in the Pyrene Molecule and Its Cation

Among polycyclic aromatic hydrocarbons, pyrene is widely used as an optical probe thanks to its peculiar ultraviolet absorption and infrared emission features. Interestingly, this molecule is also an abundant component of the interstellar medium, where it is detected via its unique spectral fingerprints. In this work, we present a comprehensive first-principles study on the electronic and vibrational response of pyrene and its cation to ultrafast, coherent pulses in resonance with their optically active excitations in the ultraviolet region. The analysis of molecular symmetries, electronic structure, and linear optical spectra is used to interpret transient absorption spectra and kinetic energy spectral densities computed for the systems excited by ultrashort laser fields. By disentangling the effects of the electronic and vibrational dynamics via ad hoc simulations with stationary and moving ions, and, in specific cases, with the aid of auxiliary model systems, we rationalize that the nuclear motion is mainly harmonic in the neutral species, while strong anharmonic oscillations emerge in the cation, driven by electronic coherence. Our results provide additional insights into the ultrafast vibronic dynamics of pyrene and related compounds and set the stage for future investigations on more complex carbon-conjugated molecules.

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

Linear-response time-dependent density-functional theory (TDDFT) can describe excitonic features in the optical spectra of insulators and semiconductors, using exchange-correlation (xc) kernels behaving as -1/k^{2} to leading order. We show how excitons can be modeled in real-time TDDFT, using an xc vector potential constructed from approximate, long-range corrected xc kernels. We demonstrate, for various materials, that this real-time approach is consistent with frequency-dependent linear response, gives access to femtosecond exciton dynamics following short-pulse excitations, and can be extended with some caution into the nonlinear regime.

Simulated field-modulated x-ray absorption in titania

In this paper, we present a method to compute the x-ray absorption near-edge structure (XANES) spectra of solid-state transition metal oxides using real-time time-dependent density functional theory, including spin-orbit coupling effects. This was performed on bulk-mimicking anatase titania (TiO₂) clusters, which allows for the use of hybrid functionals and atom-centered all electron basis sets. Furthermore, this method was employed to calculate the shifts in the XANES spectra of the Ti L-edge in the presence of applied electric fields to understand how external fields can modify the electronic structure, and how this can be probed using x-ray absorption spectroscopy. Specifically, the onset of t_2g peaks in the Ti L-edge was observed to red shift and the e_g peaks were observed to blue shift with increasing fields, attributed to changes in the hybridization of the conduction band (3d) orbitals.

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light-matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).

Excitons in two-dimensional atomic layer materials from time-dependent density functional theory: mono-layer and bi-layer hexagonal boron nitride and transition-metal dichalcogenides

Time-dependent density functional theory has been applied to the calculation of absorption spectra for two dimensional atomic layer materials: mono-layer and bi-layer hexagonal boron nitride and mono-layer transition metal dichalcogenides.