Fully Atomistic Real-Time Simulations of Transient Absorption Spectroscopy

Transient-absorption spectroscopy of dendrimers via nonadiabatic excited-state dynamics simulations

The efficiency of light-harvesting and energy transfer in multi-chromophore ensembles underpins natural photosynthesis. Dendrimers are highly branched synthetic multi-chromophoric conjugated supra-molecules that mimic these natural processes. After photoexcitation, their repeated units participate in a number of intramolecular electronic energy relaxation and redistribution pathways that ultimately funnel to a sink. Here, a model four-branched dendrimer with a pyrene core is theoretically studied using nonadiabatic molecular dynamics simulations. We evaluate excited-state photoinduced dynamics of the dendrimer, and demonstrate on-the-fly simulations of its transient absorption pump-probe (TA-PP) spectra. We show how the evolutions of the simulated TA-PP spectra monitor in real time photoinduced energy relaxation and redistribution, and provide a detailed microscopic picture of the relevant energy-transfer pathways. To the best of our knowledge, this is the first of this kind of on-the-fly atomistic simulation of TA-PP signals reported for a large molecular system.

Simulating ultrafast transient absorption spectra from first principles using a time-dependent configuration interaction probe

Transient absorption spectroscopy (TAS) is among the most common ultrafast photochemical experiments, but its interpretation remains challenging. In this work, we present an efficient and robust method for simulating TAS signals from first principles. Excited-state absorption and stimulated emission (SE) signals are computed using time-dependent complete active space configuration interaction (TD-CASCI) simulations, leveraging the robustness of time-domain simulation to minimize electronic structure failure. We demonstrate our approach by simulating the TAS signal of 1'-hydroxy-2'-acetonapthone (HAN) from ab initio multiple spawning nonadiabatic molecular dynamics simulations. Our results are compared to gas-phase TAS data recorded from both jet-cooled (T ∼ 40 K) and hot (∼403 K) molecules via cavity-enhanced TAS (CE-TAS). Decomposition of the computed spectrum allows us to assign a rise in the SE signal to excited-state proton transfer and the ultimate decay of the signal to relaxation through a twisted conical intersection. The total cost of computing the observable signal (∼1700 graphics processing unit hours for ∼4 ns of electron dynamics) was markedly less than that of performing the ab initio multiple spawning calculations used to compute the underlying nonadiabatic dynamics.

Prediction challenge: First principles simulation of the ultrafast electron diffraction spectrum of cyclobutanone

Computer simulation has long been an essential partner of ultrafast experiments, allowing the assignment of microscopic mechanistic detail to low-dimensional spectroscopic data. However, the ability of theory to make a priori predictions of ultrafast experimental results is relatively untested. Herein, as a part of a community challenge, we attempt to predict the signal of an upcoming ultrafast photochemical experiment using state-of-the-art theory in the context of preexisting experimental data. Specifically, we employ ab initio Ehrenfest with collapse to a block mixed quantum-classical simulations to describe the real-time evolution of the electrons and nuclei of cyclobutanone following excitation to the 3s Rydberg state. The gas-phase ultrafast electron diffraction (GUED) signal is simulated for direct comparison to an upcoming experiment at the Stanford Linear Accelerator Laboratory. Following initial ring-opening, dissociation via two distinct channels is observed: the C3 dissociation channel, producing cyclopropane and CO, and the C2 channel, producing CH2CO and C2H4. Direct calculations of the GUED signal indicate how the ring-opened intermediate, the C2 products, and the C3 products can be discriminated in the GUED signal. We also report an a priori analysis of anticipated errors in our predictions: without knowledge of the experimental result, which features of the spectrum do we feel confident we have predicted correctly, and which might we have wrong?

On-the-fly simulation of time-resolved fluorescence spectra and anisotropy

We combine on-the-fly trajectory surface hopping simulations and the doorway-window representation of nonlinear optical response functions to create an efficient protocol for the evaluation of time- and frequency-resolved fluorescence (TFRF) spectra and anisotropies of the realistic polyatomic systems. This approach gives the effective description of the proper (e.g., experimental) pulse envelopes, laser field polarizations, and the proper orientational averaging of TFRF signals directly from the well-established on-the-fly nonadiabatic dynamic simulations without extra computational cost. To discuss the implementation details of the developed protocol, we chose cis-azobenzene as a prototype to simulate the time evolution of the TFRF spectra governed by its nonadiabatic dynamics. The results show that the TFRF is determined by the interplay of several key factors, i.e., decays of excited-state populations, evolution of the transition dipole moments along with the dynamic propagation, and scaling factor of the TFRF signals associated with the cube of emission frequency. This work not only provides an efficient and effective approach to simulate the TFRF and anisotropies of realistic polyatomic systems but also discusses the important relationship between the TFRF signals and the underlining nonadiabatic dynamics.

Multiconfigurational Calculations and Photodynamics Describe Norbornadiene Photochemistry

Storing solar energy is a vital component of using renewable energy sources to meet the growing demands of the global energy economy. Molecular solar thermal (MOST) energy storage is a promising means to store solar energy with on-demand energy release. The light-induced isomerization reaction of norbornadiene (NBD) to quadricyclane (QC) is of great interest because of the generally high energy storage density (0.97 MJ kg-1) and long thermal reversion lifetime (t1/2,300K = 8346 years). However, the mechanistic details of the ultrafast excited-state [2 + 2]-cycloaddition are largely unknown due to the limitations of experimental techniques in resolving accurate excited-state molecular structures. We now present a full computational study on the excited-state deactivation mechanism of NBD and its dimethyl dicyano derivative (DMDCNBD) in the gas phase. Our multiconfigurational calculations and nonadiabatic molecular dynamics simulations have enumerated the possible pathways with 557 S2 trajectories of NBD for 500 fs and 492 S1 trajectories of DMDCNBD for 800 fs. The simulations predicted the S2 and S1 lifetimes of NBD (62 and 221 fs, respectively) and the S1 lifetime of DMDCNBD (190 fs). The predicted quantum yields of QC and DCQC are 10 and 43%, respectively. Our simulations also show the mechanisms of forming other possible reaction products and their quantum yields.

Photoinduced Dynamics at the Water/TiO_{2}(101) Interface

We present a femtosecond time-resolved optical pump-soft x-ray probe photoemission study in which we follow the dynamics of charge transfer at the interface of water and anatase TiO_{2}(101). By combining our observation of transient oxygen O 1s core level peak shifts at submonolayer water coverages with Ehrenfest molecular dynamics simulations we find that ultrafast interfacial hole transfer from TiO_{2} to molecularly adsorbed water is completed within the 285 fs time resolution of the experiment. This is facilitated by the formation of a new hydrogen bond between an O_{2c} site at the surface and a physisorbed water molecule. The calculations fully corroborate our experimental observations and further suggest that this process is preceded by the efficient trapping of the hole at the surface of TiO_{2} by hydroxyl species (-OH), that form following the dissociative adsorption of water. At a water coverage exceeding a monolayer, interfacial charge transfer is suppressed. Our findings are directly applicable to a wide range of photocatalytic systems in which water plays a critical role.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Ultrafast Internal Conversion Dynamics through the on-the-Fly Simulation of Transient Absorption Pump-Probe Spectra with Different Electronic Structure Methods

An on-the-fly surface-hopping simulation protocol is developed for the evaluation of transient absorption (TA) pump-probe (PP) signals of molecular systems exhibiting internal conversion to the electronic ground state. We study the nonadiabatic dynamics of azomethane and the associating TA PP spectra at three levels of the electronic-structure theory, OM2/MRCI, SA-CASSCF, and XMS-CASPT2. The impact of these methods on the population dynamics and time-resolved TA PP signals is substantially different. This difference is attributed to the strong non-Condon effects that must be taken into account for the proper understanding and interpretation of time-resolved TA PP signals of nonadiabatic polyatomic systems. This shows that the combination of the dynamical and spectral simulations definitely provides more accurate and detailed information on the microscopic mechanisms of photophysical and photochemical processes. Hence the simulation of time-resolved spectroscopic signals provides another important dimension to examine the accuracy of quantum chemistry methods.

Probing interfacial charge transfer in the heterojunctions for photocatalysis

Photocatalytic reactions can sustainably employ inexhaustible solar energy for environmental remediation and conversion of photon energy into chemical energy, and thereby show great potential in alleviating the environmental stress and...

Spectral Fingerprint of Excited-State Energy Transfer in Dendrimers through Polarization-Sensitive Transient-Absorption Pump-Probe Signals: On-the-Fly Nonadiabatic Dynamics Simulations

The time-resolved polarization-sensitive transient-absorption (TA) pump-probe (PP) spectra are simulated using on-the-fly surface-hopping nonadiabatic dynamics and the doorway-window representation of nonlinear spectroscopy. A dendrimer model system composed of two linear phenylene ethynylene units (2-ring and 3-ring) is taken as an example. The ground-state bleach (GSB), stimulated emission (SE), and excited-state absorption (ESA) contributions as well as the total TA PP signals are obtained and carefully analyzed. It is shown that intramolecular excited-state energy transfer from the 2-ring unit to the 3-ring unit can be conveniently identified by employing pump and probe pulses with different polarizations. Our results demonstrate that time-resolved polarization-sensitive TA PP signals provide a powerful tool for the elucidation of excited-state energy-transfer pathways, notably in molecular systems possessing several optically bright nonadiabatically coupled electronic states with different orientations of transition dipole moments.

Photoinduced charge-transfer in chromophore-labeled gold nanoclusters: quantum evidence of the critical role of ligands and vibronic couplings

The electron flow between a metallic aggregate and an organic molecule after excitation with light is a crucial step on which hybrid photovoltaic nanomaterials are based. So far, designing such devices with the help of theoretical approaches has been heavily limited by the computational cost of quantum dynamics models able to track the evolution of the excited states over time. In this article, we present the first application of the time-dependent density functional tight-binding (TD-DFTB) method for an experimental nanometer-sized gold-organic system consisting of a hexyl-protected Au₂₅ cluster labelled with a pyrene fluorophore, in which the fluorescence quenching of the pyrene is attributed to the electron transfer from the metallic cluster to the dye. The full quantum rationalization of the electron transfer is attained through quantum dynamics simulations, highlighting the crucial role of the protecting ligand shell in electron transfer, as well as the coupling with nuclear movement. This work paves the way towards the fast and accurate theoretical design of optoelectronic nanodevices.

Experimental determination of charge carrier dynamics in carbon nitride heterojunctions

Carbon nitride (CN_x) is an emerging photocatalyst with the potential to efficiently produce solar fuels. CN_x heterojunctions often show significant photocatalytic activity improvements. We review the charge carrier dynamics in a range of CN_x heterojunctions including carbon-based material, black phosphorus, Ru complexes, molybdenum sulphide and metal phosphides. Time resolved photoluminescence (TRPL) and transient absorption spectroscopy (TAS) were the most common techniques employed for experimental charge carrier dynamics measurements. The low photoluminescence quantum yield of CN_x appeared to limit the depth of conclusions from TRPL, with both lengthening and shortening of the TRPL lifetimes observed and attributed to enhanced charge separation. Overall, the charge carrier dynamics studies often showed a relative lifetime change of ∼2-fold and an activity improvement of >10-fold. We highlight the need for the use of a wider range of techniques to monitor the charge carrier dynamics for conclusive determination of photophysics-activity relationships and elucidation of improvement mechanisms.

Nature of Optical Excitations in Porphyrin Crystals: A Joint Experimental and Theoretical Study

The nature of optical excitations and the spatial extent of excitons in organic semiconductors, both of which determine exciton diffusion and carrier mobilities, are key factors for the proper understanding and tuning of material performances. Using a combined experimental and theoretical approach, we investigate the excitonic properties of meso-tetraphenyl porphyrin-Zn(II) crystals. We find that several bands contribute to the optical absorption spectra, beyond the four main ones considered here as the analogue to the four frontier molecular orbitals of the Gouterman model commonly adopted for the isolated molecule. By using many-body perturbation theory in the GW and Bethe-Salpeter equation approach, we interpret the experimental large optical anisotropy as being due to the interplay between long- and short-range intermolecular interactions. In addition, both localized and delocalized excitons in the π-stacking direction are demonstrated to determine the optical response, in agreement with recent experimental observations reported for organic crystals with similar molecular packing.

Ultrafast Dynamics of Nonrigid Zinc-Porphyrin Arrays Mimicking the Photosynthetic "Special Pair"

Conjugated porphyrin arrays are heavily investigated as efficient molecular systems for photosynthesis and photocatalysis. Recently, a series of one-, two-, and six-zinc-porphyrin arrays, noncovalently linked through benzene-based hubs, have been synthesized with the aim of mimicking the structure and function of the bacteriochlorophyll "special pair" in photosynthetic reaction centers. The excitonically coupled porphyrin subunits are expected to activate additional excited state relaxation channels with respect to the monomer. Here, we unveil the appearance of such supramolecular electronic interactions using ultrafast transient absorption spectroscopy with sub-25 fs time resolution. Upon photoexcitation of the Soret band, we resolve energy trapping within ∼150 fs in a delocalized dark excitonic manifold. Moreover, excitonic interactions promote an additional fast internal conversion from the Q-band to the ground state with an efficiency of up to 60% in the hexamer. These relaxation pathways appear to be common loss channels that limit the lifetime of the exciton states in noncovalently bound molecular aggregates.

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green's functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

Machine learning enables long time scale molecular photodynamics simulations

Photo-induced processes are fundamental in nature but accurate simulations of their dynamics are seriously limited by the cost of the underlying quantum chemical calculations, hampering their application for long time scales. Here we introduce a method based on machine learning to overcome this bottleneck and enable accurate photodynamics on nanosecond time scales, which are otherwise out of reach with contemporary approaches. Instead of expensive quantum chemistry during molecular dynamics simulations, we use deep neural networks to learn the relationship between a molecular geometry and its high-dimensional electronic properties. As an example, the time evolution of the methylenimmonium cation for one nanosecond is used to demonstrate that machine learning algorithms can outperform standard excited-state molecular dynamics approaches in their computational efficiency while delivering the same accuracy.