Nonadiabatic Excited-state Dynamics Research Articles

GPU-accelerated on-the-fly nonadiabatic semiclassical dynamics.

GPU-accelerated on-the-fly nonadiabatic dynamics is enabled by interfacing the linearized semiclassical dynamics approach with the TeraChem electronic structure program. We describe the computational workflow of the "PySCES" code interface, a Python code for semiclassical dynamics with on-the-fly electronic structure, including parallelization over multiple GPU nodes. We showcase the abilities of this code and present timings for two benchmark systems: fulvene solvated in acetonitrile and a charge transfer system in which a photoexcited zinc-phthalocyanine donor transfers charge to a fullerene acceptor through multiple electronic states on an ultrafast timescale. Our implementation paves the way for an efficient semiclassical approach to model the nonadiabatic excited state dynamics of complex molecules, materials, and condensed phase systems.

Machine-Learning Accelerated Description of Polarons in 2D Lead Halide Perovskites

MLFF-based MD trajectories to compute NAC and excited-state dynamics. At each stage, results based on machine learning are compared to traditional ab initio based electronic dissipative dynamics.The limits of this method are explored as we increase the timescale of the calculation beyond those allowed by ab initio molecular dynamics and to determine the effectiveness of MLFF-MD when used for excited-state nonadiabatic dynamics simulations using atomistic models that include charge injection or when an excited-stated potential energy surface is used to account for changes due to an excitation.

Ultrafast Hot Exciton Nonadiabatic Excited-State Dynamics in Green Fluorescent Protein Chromophore Analogue.

The ultrafast high-energy nonadiabatic excited-state dynamics of the benzylidenedimethylimidazolinone chromophore dimer has been investigated using an electronic structure method coupled with on-the-fly quantitative wave function analysis to gain insight into the photophysics of hot excitons in biological systems. The dynamical simulation provides a rationalization of the behavior of the exciton in a dimer after the photoabsorption of light to higher-energy states. The results suggest that hot exciton localization within the manifold of excited states is caused by the hindrance of torsional rotation due to imidazolinone (I) or phenolate (P) bonds i.e., ΦI- or ΦP-dihedral rotation, in the monomeric units of a dimer. This hindrance arises due to weak π-π stacking interaction in the dimer, resulting in an energetically uphill excited-state barrier for ΦI- and ΦP-twisted rotation, impeding the isomerization process in the chromophore. Thus, this study highlights the potential impact of the weak π-π interaction in regulating the photodynamics of the green fluorescent protein chromophore derivatives.

Ultrafast Excited-State Nonadiabatic Dynamics in Pt(II) Donor-Bridge-Acceptor Assemblies: A Quantum Approach for Optical Control.

The ultrafast nonadiabatic excited state dynamics of (PTZ-N-benzyl-acetylide) (trans-bis-trimethylphosphine) Pt(II) (acetylide-NDI-bis-methyl) 1, representative of a series of Pt(II) donor-bridge-acceptor assemblies experimentally studied by the Weinstein group, University of Sheffield, is investigated by means of wavepacket propagations based on the multiconfiguration time-dependent Hartree (MCTDH) method. On the basis of electronic structure data obtained at the time-dependent density functional theory (TD-DFT) level, the subpicosecond decay is simulated by solving an 11 electronic states multimode problem, up to 18 vibrational normal modes, including both spin-orbit coupling (SOC) and vibronic coupling. A careful analysis of the results, within the diabatic representation, provides the key features of the spin-vibronic mechanism at work in this complex, distinguishing between the spin-orbit and vibronically activated ultrafast processes within the excited states manifold. The knowledge of the key active normal modes that promote selectively the population of specific electronic excited states opens a route toward optical control by selectively exciting these modes in order to drive the associated nonadiabatic processes. Relevant simulations, over 2 ps, are proposed to assess the impact of these selective vibrational excitations on the branching ratio between the primary photoproducts, namely, bridge-acceptor charge-transfer (CT) and donor-acceptor charge-separated (CS) electronic states. Whereas the excitation of the localized acetylide bridge C≡C bond stretching does not modify drastically the population of the low-lying electronic states within the first two ps, vibrational excitation of the out-of-plane twisting motion of the N-benzyl group linked to the donor entity favors the population of the 1,3CS states at the expense of the lowest 1,3CT states. This quantum study opens the route to IR optical control experiments based on the specific alteration of vibrational normal modes that activate vibronic couplings between key electronic excited states. However, the presence of critical crossings along the PES channels associated with these normal modes and the role of concurrent SOC driven ultrafast transfers of population should not be underestimated.

The nature of carotenoid S* state and its role in the nonphotochemical quenching of plants

In plants, light-harvesting complexes serve as antennas to collect and transfer the absorbed energy to reaction centers, but also regulate energy transport by dissipating the excitation energy of chlorophylls. This process, known as nonphotochemical quenching, seems to be activated by conformational changes within the light-harvesting complex, but the quenching mechanisms remain elusive. Recent spectroscopic measurements suggest the carotenoid S* dark state as the quencher of chlorophylls’ excitation. By investigating lutein embedded in different conformations of CP29 (a minor antenna in plants) via nonadiabatic excited state dynamics simulations, we reveal that different conformations of the complex differently stabilize the lutein s-trans conformer with respect to the dominant s-cis one. We show that the s-trans conformer presents the spectroscopic signatures of the S* state and rationalize its ability to accept energy from the closest excited chlorophylls, providing thus a relationship between the complex’s conformation and the nonphotochemical quenching.

Isomerization of surface functionalized SWCNTs and the critical influence on photoluminescence: static calculations and excited-state dynamics simulations.

Single-walled carbon nanotubes (SWCNTs) functionalized with sparse surface chemical groups are promising for a variety of optical applications such as quantum information and bio-imaging. However, the luminescence efficiencies and stability, two key aspects, undoubtedly govern their practical usage. Herein, we assess the surface migration of oxygen and triazine groups on as-modified SWCNT fragments by adopting transition state theory and explore the de-excitation of oxygen-functionalized SWCNT fragments by performing non-adiabatic excited-state dynamics simulations. According to the predicted moderate or even small reaction barriers, the migration of both oxygen and triazine groups is feasible from an sp3 defect configuration forming an energetically more stable sp2 configuration at moderate or even room temperatures. Such isomerization leads to drastically different light emission capabilities as indicated by the large or zero oscillator strengths. During the dynamics simulations, the lowest excited singlet (S1) state rapidly decays in energy within 20 fs and then fluctuates until the end, providing insights into the emission mechanism of SWCNTs. This study highlights the potential intrinsic limitations of surface-functionalized SWCNTs for luminescence applications.

Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics.

When molecules are coupled to an optical cavity, new light-matter hybrid states, so-called polaritons, are formed due to quantum light-matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light-matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light-matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule-cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.

Tuning UV Pump X-ray Probe Spectroscopy on the Nitrogen K Edge Reveals the Radiationless Relaxation of Pyrazine: Ab Initio Simulations Using the Quasiclassical Doorway-Window Approximation.

Transient absorption UV pump X-ray probe spectroscopy has been established as a versatile technique for the exploration of ultrafast photoinduced dynamics in valence-excited states. In this work, an ab initio theoretical framework for the simulation of time-resolved UV pump X-ray probe spectra is presented. The method is based on the description of the radiation-matter interaction in the classical doorway-window approximation and a surface-hopping algorithm for the nonadiabatic nuclear excited-state dynamics. Using the second-order algebraic-diagrammatic construction scheme for excited states, UV pump X-ray probe signals were simulated for the carbon and nitrogen K edges of pyrazine, assuming a duration of 5 fs of the UV pump and X-ray probe pulses. It is predicted that spectra measured at the nitrogen K edge carry much richer information about the ultrafast nonadiabatic dynamics in the valence-excited states of pyrazine than those measured at the carbon K edge.

Chlorine Adsorption on TiO2(110)/Water Interface: Nonadiabatic Molecular Dynamics Simulations for Photocatalytic Water Splitting

Chloride is one of the most abundant ions in sea water, which is more available than fresh water. Due to lack of H2O adsorbate states near the valence band maximum (VBM) edge, the difficulty of water dissociation incidents has been reported on the rutile TiO2 surface as the excitation energy is around the band gap energy of TiO2. It is interesting whether the extra chloride can be a benefit to the water dissociation or not. In this study, the models of chlorine adatoms placed on the rutile TiO2 (110)/water interface are constructed using ab initio methods. The time-dependent spatial charges, bond-lengths of water molecules, and Hirshfeld charges are calculated by real-time time-dependent density functional theory and the Ehrenfest dynamics theory for investigating the excited state nonadiabatic dynamics of water dissociation. This study presents two photoinduced water-splitting pathways related to chlorine and analyzes the photogenerated hole along the reactions. The first step of water dissociation relies on the localized competition of oxygen charges between the dissociated water and the bridge site of TiO2 for transforming the water into hydroxyl and hydrogen by photoinduced driving force.

Charge-Transfer Luminescence in a Molecular Donor-Acceptor Complex: Computational Insights.

Donor-acceptor molecular complexes are a popular class of materials utilizing charge-transfer states for practical applications. A recent class of donor-acceptor dyads based on the fluorescent BODIPY functionalized with triphenylamine (TPA) shows the peculiar property of dual fluorescence. It is hypothesized that instead of the sensitized charge-transfer state being optically dark, it provides an additional bright radiative pathway. Here we use time-dependent density functional theory to characterize the energetic alignment of excitonic and charge-transfer states in a BODIPY-TPA molecular complex. We observe that using a long-range exchange corrected functional in combination with state-specific solvation scheme gives a qualitatively correct alignment of the exciton and charge-transfer states and an enhancement in oscillator strength for the equilibrium solvated charge-transfer state, in agreement with experiment. This work provides rationalization of charge-transfer state emission and provides a foundation to explore charge-transfer using ab initio excited-state nonadiabatic dynamics.

Nonadiabatic Excited State Dynamics of Organic Chromophores: Take-Home Messages.

Nonadiabatic excited state dynamics are important in a variety of processes. Theoretical and experimental developments have allowed for a great progress in this area, while combining the two is often necessary and the best approach to obtain insight into the photophysical behavior of molecules. In this Feature Article we use examples of our recent work combining time-resolved photoelectron spectroscopy with theoretical nonadiabatic dynamics to highlight important lessons we learned. We compare the nonadiabatic excited state dynamics of three different organic molecules with the aim of elucidating connections between structure and dynamics. Calculations and measurements are compared for uracil, 1,3-cyclooctadiene, and 1,3-cyclohexadiene. The comparison highlights the role of rigidity in influencing the dynamics and the difficulty of capturing the dynamics accurately with calculations.

Infinitene: Computational Insights from Nonadiabatic Excited State Dynamics.

Progress in organic synthesis opens exploration of a rich diversity of molecules with interesting new structural topologies. This is the case of a recently synthesized helically twisted figure-eight molecule coined infinitene. The molecule belongs to a numerous family of looped polyarenes, where the degree of π-conjugation is controlled by high strain energies and steric hindrances. A particular balance of these ingredients leads to unusual optoelectronic properties potentially suitable for a range of applications in nanoelectronics and photonics. Due to its recent discovery, the photophysical properties of infinitene remain unexplored. In this Letter, atomistic nonadiabatic excited state molecular dynamics modeling unveils unique features of intramolecular electronic and vibrational energy relaxation and redistribution that take place after molecular photoexcitation. Our results detail relationships between optical and electronic properties providing useful knowledge for future molecular designs related to infinitene.

Ultrafast Excited-State Dynamics of Carotenoids and the Role of the SX State.

Carotenoids are natural pigments with multiple roles in photosynthesis. They act as accessory pigments by absorbing light where chlorophyll absorption is low, and they quench the excitation energy of neighboring chlorophylls under high-light conditions. The function of carotenoids depends on their polyene-like structure, which controls their excited-state properties. After light absorption to their bright S2 state, carotenoids rapidly decay to the optically dark S1 state. However, ultrafast spectroscopy experiments have shown the signatures of another dark state, termed SX. Here we shed light on the ultrafast photophysics of lutein, a xanthophyll carotenoid, by explicitly simulating its nonadiabatic excited-state dynamics in solution. Our simulations confirm the involvement of SX in the relaxation toward S1 and reveal that it is formed through a change in the nature of the S2 state driven by the decrease in the bond length alternation coordinate of the carotenoid conjugated chain.

Theoretical investigation on the fluorescent sensing mechanism for recognizing formaldehyde: TDDFT calculation and excited-state nonadiabatic dynamics

Inspired by the activity-based sensing method, the hydrazine-modified naphthalene derivative (Naph1) was synthesized and used as a fluorescent probe to detect formaldehyde (FA) in living cells. Through the condensation reaction between the probe Naph1 and analyte FA, researchers observed a ∼14 folds enhancement of fluorescent signal around 510 nm in an experiment, realizing the high selectivity and sensitivity detection of FA. However, a theoretical understanding of the sensing mechanism was not provided in the experimental work. Given this, the light-up fluorescent detecting mechanism was in-depth unveiled by performing the time-dependent density functional theory (TDDFT) and the complete active space self-consistent field (CASSCF) theoretical calculations on excited-state intramolecular proton transfer (ESIPT) and non-adiabatic excited-state dynamics simulation. The deactivation channel of S1/T2 intersystem crossing (ISC) was turned off to successfully recognize FA. Insight into the ESIPT-based fluorescent detecting mechanism indicated that ESIPT was essential to light-up fluorescent probes. This work would provide a new viewpoint to develop ESIPT-based fluorescent probes for detecting reactive carbon species in vivo or vitio.

Excited state non-adiabatic dynamics of large photoswitchable molecules using a chemically transferable machine learning potential

Light-induced chemical processes are ubiquitous in nature and have widespread technological applications. For example, photoisomerization can allow a drug with a photo-switchable scaffold such as azobenzene to be activated with light. In principle, photoswitches with desired photophysical properties like high isomerization quantum yields can be identified through virtual screening with reactive simulations. In practice, these simulations are rarely used for screening, since they require hundreds of trajectories and expensive quantum chemical methods to account for non-adiabatic excited state effects. Here we introduce a diabatic artificial neural network (DANN), based on diabatic states, to accelerate such simulations for azobenzene derivatives. The network is six orders of magnitude faster than the quantum chemistry method used for training. DANN is transferable to azobenzene molecules outside the training set, predicting quantum yields for unseen species that are correlated with experiment. We use the model to virtually screen 3100 hypothetical molecules, and identify novel species with high predicted quantum yields. The model predictions are confirmed using high-accuracy non-adiabatic dynamics. Our results pave the way for fast and accurate virtual screening of photoactive compounds.

Resolving competing conical intersection pathways: time-resolved X-ray absorption spectroscopy of trans-1,3-butadiene.

Time-resolved X-ray absorption spectroscopy is emerging as a uniquely powerful tool to probe coupled electronic-nuclear dynamics in photo-excited molecules. Theoretical studies to date have established that time-resolved X-ray absorption spectroscopy is an atom-specific probe of excited-state wave packet passage through a seam of conical intersections (CIs). However, in many molecular systems, there are competing dynamical pathways involving CIs of different electronic and nuclear character. Discerning these pathways remains an important challenge. Here, we demonstrate that time-resolved X-ray absorption spectroscopy (TRXAS) has the potential to resolve competing channels in excited-state non-adiabatic dynamics. Using the example of 1,3-butadiene, we show how TRXAS discerns the different electronic structures associated with passage through multiple conical intersections. trans-1,3-Butadiene exhibits a branching between polarized and radicaloid pathways associated with ethylenic "twisted-pyramidalized" and excited-state cis-trans isomerization dynamics, respectively. The differing electronic structures along these pathways give rise to different XAS signals, indicating the possibility of resolving them. Furthermore, this indicates that XAS, and other core-level spectroscopic techniques, offer the appealing prospect of directly probing the effects of selective chemical substitution and its ability to affect chemical control over excited-state molecular dynamics.

Nonadiabatic excited-state dynamics of ReCl(CO)3(bpy) in two different solvents.

We present a study of excited-states relaxation of the complex ReCl(CO)3(bpy) (bpy = 2,2-bipyridine) using a nonadiabatic TD-DFT dynamics on spin-mixed potential energy surfaces in explicit acetonitrile (ACN) and dimethylsulfoxide (DMSO) solutions up to 800 fs. ReCl(CO)3(bpy) belongs to a group of important photosensitizers which show ultrafast biexponential subpicosecond fluorescence decay kinetics. The choice of solvents was motivated by the different excited-state relaxation dynamics observed in subpicosecond time-resolved IR (TRIR) experiments. Simulations of intersystem crossing (ISC) showed the development of spin-mixed states in both solvents. Transformation of time-dependent populations of spin-mixed states enabled to monitor the temporal evolution of individual singlet and triplet states, fitting of bi-exponential decay kinetics, and simulating the time-resolved fluorescence spectra that show only minor differences between the two solvents. Analysis of structural relaxation and solvent reorganization employing time-resolved proximal distribution functions pointed to the factors influencing the fluorescence decay time constants. Nonadiabatic dynamics simulations of time-evolution of electronic, molecular, and solvent structures emerge as a powerful technique to interpret time-resolved spectroscopic data and ultrafast photochemical reactivity.

Solvent Effects on Ultrafast Photochemical Pathways.

Photochemical reactions are increasingly being used for chemical and materials synthesis, for example, in photoredox catalysis, and generally involve photoexcitation of molecular chromophores dissolved in a liquid solvent. The choice of solvent influences the outcomes of the photochemistry because solute-solvent interactions modify the energies of and crossings between electronic states of the chromophores, and they affect the evolving structures of the photoexcited molecules. Ultrafast laser spectroscopy methods with femtosecond to picosecond time resolution can resolve the dynamics of these photoexcited molecules as they undergo structural and electronic changes, relax back to the ground state, dissipate their excess internal energy to the surrounding solvent, or undergo photochemical reactions. In this Account, we illustrate how experimental studies using ultrafast lasers can reveal the influences that different solvents or cosolutes exert on the photoinduced nonadiabatic dynamics of internal conversion and intersystem crossing in nonradiative relaxation pathways. Although the environment surrounding a solute molecule is rapidly changing, with fluctuations in the coordination to neighboring solvent molecules occurring on femtosecond or picosecond time scales, we show that it is possible to photoexcite selectively only those molecular chromophores transiently experiencing specific solute-solvent interactions such as intermolecular hydrogen bonding.The effects of different solvation environments on the photodynamics are illustrated using four selected examples of photochemical processes in which the solvent has a marked effect on the outcomes. We first consider two aromatic carbonyl compounds, benzophenone and acetophenone, which are known to undergo fast intersystem crossing to populate the first excited triplet state on time scales of a few picoseconds. We show that the nonadiabatic excited-state dynamics are modified by transient hydrogen bonding of the carbonyl group to a protic solvent or by coordination to a metal cation cosolute. We then examine how different solvents modify the competition between two alternative relaxation pathways in a photoexcited UVA-sunscreen molecule, diethylamino hydroxybenzoyl hexyl benzoate (DHHB). This relaxation back to the ground electronic state is an essential part of the effective operation of the sunscreen compound, but the dynamics are sensitive to the surrounding environment. Finally, we consider how solvents of different polarity affect the energies and lifetimes of excited states with locally excited or charge-transfer character in heterocyclic organic compounds used as excited-state electron donors for photoredox catalysis. With these and other examples, we seek to develop a molecular level understanding of how the choice of solution environment might be used to control the outcomes of photochemical reactions.

Nonadiabatic dynamics of cobalt tricarbonyl nitrosyl for ligand dissociation induced by electronic excitation

We utilize real-time time-dependent density functional theory and Ehrenfest dynamics scheme to investigate excited-state nonadiabatic dynamics of ligand dissociation of cobalt tricarbonyl nitrosyl, Co(CO)3NO, which is a precursor used for cobalt growth in advanced technologies, where the precursor’s reaction is enhanced by electronic excitation. Based on the first-principles calculations, we demonstrate two dissociation pathways of the NO ligand on the precursor. Detailed electronic structures are further analyzed to provide an insight into dynamics following the electronic excitations. This study sheds light on computational demonstration and underlying mechanism of the electronic-excitation-induced dissociation, especially in molecules with complex chemical bonds such as the Co(CO)3NO.

Exploring the effects of quantum decoherence on the excited-state dynamics of molecular systems

Fewest-switches surface hopping (FSSH) is employed in order to investigate the nonadiabatic excited-state dynamics of thiophene and related compounds and hence to establish a connection between the electronic system, the critical points in configuration space and the deactivation dynamics. The potential-energy surfaces of the studied molecules were calculated with complete active space self-consistent field and time-dependent density-functional theory. They are analyzed thoroughly to locate and optimize minimum-energy conical intersections, which are essential to the dynamics of the system. The influence of decoherence on the dynamics is examined by employing different decoherence schemes. We find that irrespective of the employed decoherence algorithm, the population dynamics of thiophene give results which are sound with the expectations grounded on the analysis of the potential-energy surface. A more detailed look at single trajectories as well as on the excited-state lifetimes, however, reveals a substantial dependence on how decoherence is accounted for. In order to connect these findings, we describe how ensemble averaging cures some of the overcoherence problems of uncorrected FSSH. Eventually, we identify carbon–sulfur bond cleavage as a common feature accompanying electronic transitions between different states in the simulations of all thiophene-related compounds studied in this work, which is of interest due to their relevance in organic photovoltaics.