The Simplest Possible Approach for Simulating S0- S1 Conical Intersections with DFT/TDDFT: Adding One Doubly Excited Configuration

Symmetry-breaking dynamics of a photoionized carbon dioxide dimer

Photoionization can initiate structural reorganization of molecular matter and drive formation of new chemical bonds. Here, we used time-resolved extreme ultraviolet (EUV) pump - EUV probe Coulomb explosion imaging of carbon dioxide dimer ionC O 2 2 + dynamics, that combined with ab initio molecular dynamics simulations, revealed unexpected asymmetric structural rearrangement. We show that ionization by the pump pulse induces rearrangement from the slipped-parallel (C2h) geometry of the neutral C O 2 dimer towards a T-shaped (C2v) structure on the ~100 fs timescale, although the most stable slipped-parallel (C2h) structure of the ionic dimer. Moreover, we find that excited states of the ionized C O 2 dimer can exhibit formation of aCO 3 moiety in theC 2 O 4 + complex that can persist even after a suitably time-delayed second photoionization in a metastableC 2 O 4 2 + dication. Our results suggest that charge asymmetry plays an important role in the ionization-induced dynamics in such dimers that are present in C O 2 rich environments.

WFOT: A Wave Function Overlap Tool between Single- and Multi-Reference Electronic Structure Methods for Spectroscopy Simulation

We report the development of a novel diagnostic tool, named wave function overlap tool (WFOT), designed to evaluate the overlap between wave functions computed at single-reference [i.e., time-dependent density functional theory or configuration interaction singles (CIS)] and multireference (i.e., CASSCF/CASPT2) electronic structure levels of theory. It relies on truncating the single- and multireference WFs to CIS-like expansions spanning the same configurational space and maximizing the molecular orbital overlap by means of a unitary transformation. To demonstrate the functionality of the tool, we calculate the transient spectrum of acetylacetone by evaluating excited state absorption signals with multireference quality on top of single-reference on-the-fly dynamics simulations. Semiautomatic spectra generation is facilitated by interfacing the tool with the COBRAMM package, which also allows one to use WFOT with several quantum chemistry codes such as Gaussian, NWChem, and OpenMolcas. Other exciting possibilities for the utilization of the code beyond the simulation of transient absorption spectroscopy are eventually discussed.

A fast and smooth one-electron approach for investigating charge transfer states and D1-D0 crossings for systems with odd numbers of electrons

We propose an efficient version of ensemble Hartree-Fock/density functional theory to calculate a set of two charge-transfer states for systems with odd-numbers of electrons. The approach does require definitions of donor/acceptor fragments; however, the approach is not very sensitive to such definitions-even in the limit of very strong electronic coupling. The key ansatz is that, by mandating that the vector space spanned by the active orbitals projects equally onto the donor and acceptor fragments, such a constraint eliminates all intra-molecular local excitations and makes it far easier to generate potential energy surfaces that are smooth over a wide region of configuration space. The method is fast, working with only two electron configurations, and should be useful for ab initio non-adiabatic dynamics in the near future.

Diagrammatic multiplet sum method (MSM) density functional theory (DFT): Investigation of the transferability of integrals in "simple" DFT-based approaches to multideterminantal problems

Kohn-Sham density functional theory (DFT) typically works well for describing dynamic correlation. Two other types of correlation, arising in the cases of degenerate (static) or quasidegenerate (nondynamic) zero-order states, represent a difficult problem for DFT. When symmetry is present, multiplet sum method (MSM) DFT [Ziegler et al., Theor. Chim. Acta 4, 877 (1977)] provides one of the earliest and simplest ways to include static correlation in DFT. MSM-DFT assumes that DFT provides a good description of single-determinant energies and uses symmetry and simple ansätze to include the effects of static correlation. This is equivalent to determining the off-diagonal matrix elements in a small configuration interaction (CI) eigenvalue problem. Our ultimate goal, however, is nondynamic correlation in cases where symmetry is inadequate for fixing the dynamic-correlation limitation of DFT. To this end, we have developed a diagrammatic approach to MSM-DFT, which does not, by itself, solve the nondynamic correlation problem in DFT but which facilitates comparison with wave function CI and so allows educated guesses of off-diagonal CI matrix elements even in the absence of symmetry. In every case, an additional exchange-only ansatz (EXAN) allows the MSM-DFT formulas to be transformed into wave function formulas. This EXAN also works for transforming time-dependent DFT into time-dependent Hartree-Fock. Although not enough to uniquely guess DFT formulas from wave function formulas, the diagrammatic approach and the EXAN provide important constraints on any guesses that might be used. We illustrate how diagrammatic MSM-DFT may be used to guess a nondynamic correlation correction for the dissociation of H2 and how diagrammatic MSM-DFT may be used to guess a nonsymmetry-based coupling element in the O2 multiplet problem, which is reasonably close to a previous symmetry-derived result.

On the description of conical intersections between excited electronic states with LR-TDDFT and ADC(2)

Conical intersections constitute the conceptual bedrock of our working understanding of ultrafast, nonadiabatic processes within photochemistry (and photophysics). Accurate calculation of potential energy surfaces within the vicinity of conical intersections, however, still poses a serious challenge to many popular electronic structure methods. Multiple works have reported on the deficiency of methods like linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) or algebraic diagrammatic construction to second-order [ADC(2)]-approaches often used in excited-state molecular dynamics simulations-to describe conical intersections between the ground and excited electronic states. In the present study, we focus our attention on conical intersections between excited electronic states and probe the ability of AA LR-TDDFT and ADC(2) to describe their topology and topography, using protonated formaldimine and pyrazine as two exemplar molecules. We also take the opportunity to revisit the performance of these methods in describing conical intersections involving the ground electronic state in protonated formaldimine-highlighting in particular how the intersection ring exhibited by AA LR-TDDFT can be perceived either as a (near-to-linear) seam of intersection or two interpenetrating cones, depending on the magnitude of molecular distortions within the branching space.

Computational insights of excited state intramolecular proton transfer (ESIPT) based fluorescent detection and imaging of γ-glutamytranspeptidase activity

γ-Glutamytranspeptidase (GGT) is an important tumor biomarker that widely appears in the tumor cells. Therefore, accurate imaging and detection of GGT activity in live cells, serum and pathological cells grasp great importance for the diagnosis, management, and treatment of cancer. Herein, 2-(2-hydroxyl-phenyl)-6-chloro-4-(3H)-quinazolinone (HPQ) is considered as the fluorophore probe for the detection of GGT activity, which is known for the typical mechanism of excited-state intramolecular proton transfer (ESIPT). All the simulations adopted to evaluate the sensing mechanism were carried out via DFT and TDDFT calculations at CAM-B3LYP/TZVP level of theory. The emission properties of HPQ and HPQ-TD are thoroughly studied to understand the photoinduced electron transfer (PET) and excited state intramolecular proton transfer (ESIPT) process. The results reveal that the fluorescence quenching of HPQ (enol form) is assigned to the PET process, whereas the large Stokes shift in fluorescence emission of HPQ (keto form) is related with ESIPT mechanism. The obtained results are further cross validated by frontier molecular orbital (FMO) analysis, geometric analysis, and potential energy curve (PEC) scanning. Our calculations provide powerful evidence for the ESIPT based sensing mechanism of HPQ (keto-enol form) for GGT activity.

Analytical gradients and derivative couplings for the TDDFT-1D method

We derive and implement analytic gradients and derivative couplings for time-dependent density functional theory plus one double (TDDFT-1D) which is a semiempirical configuration interaction method whereby the Hamiltonian is diagonalized in a basis of all singly excited configurations and one doubly excited configuration as constructed from a set of reference Kohn-Sham orbitals. We validate the implementation by comparing against finite difference values. Furthermore, we show that our implementation can locate both optimized geometries and minimum-energy crossing points along conical seams of S₁/S₀ surfaces for a set of test cases.

Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method

Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. The nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.

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.

Variational Density Functional Calculations of Excited States: Conical Intersection and Avoided Crossing in Ethylene Bond Twisting

Theoretical studies of photochemical processes require a description of the energy surfaces of excited electronic states, especially near degeneracies, where transitions between states are most likely. Systems relevant to photochemical applications are typically too large for high-level multireference methods, and while time-dependent density functional theory (TDDFT) is efficient, it can fail to provide the required accuracy. A variational, time-independent density functional approach is applied to the twisting of the double bond and pyramidal distortion in ethylene, the quintessential model for photochemical studies. By allowing for symmetry breaking, the calculated energy surfaces exhibit the correct topology around the twisted-pyramidalized conical intersection even when using a semilocal functional approximation, and by including explicit self-interaction correction, the torsional energy curves are in close agreement with published multireference results. The findings of the present work point to the possibility of using a single determinant time-independent density functional approach to simulate nonadiabatic dynamics, even for large systems where multireference methods are impractical and TDDFT is often not accurate enough.

On the inclusion of one double within CIS and TDDFT

We present an improved approach for generating a set of optimized frontier orbitals (HOMO and LUMO) that minimizes the energy of one double configuration. We further benchmark the effect of including such a double within a rigorous configuration interaction singles or a parameterized semi-empirical time-dependent density functional theory Hamiltonian for a set of test cases. Although we cannot quite achieve quantitative accuracy, the algorithm is quite robust and routinely delivers an enormous qualitative improvement to standard single-reference electronic structure calculations.

Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections

Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.

CAS without SCF-Why to use CASCI and where to get the orbitals

The complete active space self-consistent field (CASSCF) method has seen broad adoption due to its ability to describe the electronic structure of both the ground and excited states of molecules over a broader swath of the potential energy surface than is possible with the simpler Hartree-Fock approximation. However, it also has a reputation for being unwieldy, computationally costly, and un-black-box. Here, we discuss a class of alternatives, complete active space configuration interaction (CASCI) methods, paying particular attention to their application to electronic excited states. The goal of this Perspective is fourfold. First, we argue that CASCI is not merely an approximation to CASSCF, in that it can be designed to have important qualitative advantages over CASSCF. Second, we present several insights drawn from our experience experimenting with different schemes for computing orbitals to be employed in CASCI. Third, we argue that CASCI is well suited for application to nanomaterials. Finally, we reason that, with the rise in new low-scaling approaches for describing multireference systems, there is a greater need than ever to develop new methods for defining orbitals that provide an efficient and accurate description of both static correlation and electronic excitations in a limited active space.

Nonadiabatic Dynamics at Metal Surfaces: Fewest Switches Surface Hopping with Electronic Relaxation

A new scheme is proposed for modeling molecular nonadiabatic dynamics near metal surfaces. The charge-transfer character of such dynamics is exploited to construct an efficient reduced representation for the electronic structure. In this representation, the fewest switches surface hopping (FSSH) approach can be naturally modified to include electronic relaxation (ER). The resulting FSSH-ER method is valid across a wide range of coupling strengths as supported by tests applied to the Anderson-Holstein model for electron transfer. Future work will combine this scheme with ab initio electronic structure calculations.

Analytic gradients and derivative couplings for configuration interaction with all single excitations and one double excitation-En route to nonadiabatic dynamics

We present analytic gradients and derivative couplings for the simplest possible multireference configuration interaction method, CIS-1D, an electronic structure Ansatz that includes all single excitations and one lone double excitation on top of a Hartree-Fock reference state. We show that the resulting equations are numerically stable and require the evaluation of a similar number of integrals as compared to standard CIS theory; one can easily differentiate the required frontier orbitals (h and l) with minimal cost. The resulting algorithm has been implemented within the Q-Chem electronic structure package and should be immediately useful for understanding photochemistry with S₀-S₁ crossings.

Dynamical kernels for optical excitations

We discuss the physical properties and accuracy of three distinct dynamical (i.e., frequency-dependent) kernels for the computation of optical excitations within linear response theory: (i) an a priori built kernel inspired by the dressed time-dependent density-functional theory kernel proposed by Maitra et al. [J. Chem. Phys. 120, 5932 (2004)], (ii) the dynamical kernel stemming from the Bethe-Salpeter equation (BSE) formalism derived originally by Strinati [Riv. Nuovo Cimento 11, 1-86 (1988)], and (iii) the second-order BSE kernel derived by Zhang et al. [J. Chem. Phys. 139, 154109 (2013)]. The principal take-home message of the present paper is that dynamical kernels can provide, thanks to their frequency-dependent nature, additional excitations that can be associated with higher-order excitations (such as the infamous double excitations), an unappreciated feature of dynamical quantities. We also analyze, for each kernel, the appearance of spurious excitations originating from the approximate nature of the kernels, as first evidenced by Romaniello et al. [J. Chem. Phys. 130, 044108 (2009)]. Using a simple two-level model, prototypical examples of valence, charge-transfer, and Rydberg excited states are considered.

A systematic model study quantifying how conical intersection topography modulates photochemical reactions

Despite their important role in photochemistry and expected presence in most polyatomic molecules, conical intersections have been thoroughly characterized in a comparatively small number of systems. Conical intersections can confer molecular photoreactivity or photostability, often with remarkable efficacy, due to their unique structure: at a conical intersection, the adiabatic potential energy surfaces of two or more electronic states are degenerate, enabling ultrafast decay from an excited state without radiative emission, known as nonadiabatic transfer. Furthermore, the precise conical intersection topography determines fundamental properties of photochemical processes, including excited-state decay rate, efficacy, and molecular products that are formed. However, these relationships have yet to be defined comprehensively. In this article, we use an adaptable computational model to investigate a variety of conical intersection topographies, simulate resulting nonadiabatic dynamics, and calculate key photochemical observables. We varied the vibrational mode frequencies to modify conical intersection topography systematically in four primary classes of conical intersections and quantified the resulting rate, total yield, and product yield of nonadiabatic decay. The results reveal that higher vibrational mode frequencies reduce nonadiabatic transfer, but increase the transfer rate and resulting photoproduct formation. These trends can inform progress toward experimental control of photochemical reactions or tuning of molecules' photochemical properties based on conical intersections and their topography.

An Excited State Intramolecular Proton Transfer-Based Fluorescent Probe with a Large Stokes Shift for the Turn-on Detection of Cysteine: A Detailed Theoretical Exploration

DFT and TDDFT calculations are adopted to study the sensing mechanism of a turn-on-type cysteine fluorescent probe (2-(1-phenyl-imidazo[1,5-α]pyridine-3-yl)phenyl acrylate, denoted as MZC-AC). The photoinduced electron transfer (PET) process of MZC-AC and the excited state intramolecular proton transfer (ESIPT) process of MZC have been investigated in detail. We demonstrate that the fluorescence quenching of MZC-AC is ascribed to the PET mechanism and the large Stokes shift fluorescence emission of MZC is the result of the ESIPT mechanism. The results have been cross-validated by geometries, frontier molecular orbital analysis, and potential energy curve scanning. As a result, our calculations completely reproduce the experimental results and give powerful evidence for the sensing mechanism of MZC-AC for cysteine.

Hole-hole Tamm-Dancoff-approximated density functional theory: A highly efficient electronic structure method incorporating dynamic and static correlation

The study of photochemical reaction dynamics requires accurate as well as computationally efficient electronic structure methods for the ground and excited states. While time-dependent density functional theory (TDDFT) is not able to capture static correlation, complete active space self-consistent field methods neglect much of the dynamic correlation. Hence, inexpensive methods that encompass both static and dynamic electron correlation effects are of high interest. Here, we revisit hole-hole Tamm-Dancoff approximated (hh-TDA) density functional theory for this purpose. The hh-TDA method is the hole-hole counterpart to the more established particle-particle TDA (pp-TDA) method, both of which are derived from the particle-particle random phase approximation (pp-RPA). In hh-TDA, the N-electron electronic states are obtained through double annihilations starting from a doubly anionic (N+2 electron) reference state. In this way, hh-TDA treats ground and excited states on equal footing, thus allowing for conical intersections to be correctly described. The treatment of dynamic correlation is introduced through the use of commonly employed density functional approximations to the exchange-correlation potential. We show that hh-TDA is a promising candidate to efficiently treat the photochemistry of organic and biochemical systems that involve several low-lying excited states-particularly those with both low-lying ππ^* and nπ^* states where inclusion of dynamic correlation is essential to describe the relative energetics. In contrast to the existing literature on pp-TDA and pp-RPA, we employ a functional-dependent choice for the response kernel in pp- and hh-TDA, which closely resembles the response kernels occurring in linear response and collinear spin-flip TDDFT.

Configuration interaction approaches for solving quantum impurity models

We develop several configuration interaction approaches for characterizing the electronic structure of an adsorbate on a metal surface (at least in model form). When one can separate the adsorbate from the substrate, these methods can achieve a reasonable description of adsorbate on-site electron-electron correlation in the presence of a continuum of states. While the present paper is restricted to the Anderson impurity model, there is hope that these methods can be extended to ab initio Hamiltonians and provide insight into the structure and dynamics of molecule-metal surface interactions.