Shape of Multireference, Equation-of-Motion Coupled-Cluster, and Density Functional Theory Potential Energy Surfaces at a Conical Intersection

The Photophysical Path to the Triplet State in Light-Oxygen-Voltage (LOV) Domains

Upon blue-light absorption, LOV domains efficiently undergo intersystem crossing (ISC) to the triplet state. Several factors potentially contribute to this efficiency. One often proposed in the literature is the heavy atom effect of the nearby (and eventually adduct-forming) cysteine. However, some LOV domain derivatives that lack the cysteine residue also undergo ISC efficiently. Using hybrid multireference quantum mechanical/molecular mechanical (QM / MM) models, we investigated the effect of the electrostatic environment in a prototypal LOV domain, Arabidopsis thaliana Phototropin 1 LOV2 (AtLOV2), compared to the effect of the dielectric of an aqueous solution. We find that the electrostatic environment of AtLOV2 is especially well tuned to stabilize a triplet( n N , π * ) ${(n_{\rm{N}}, \pi ^{\ast} )}$ state, which we posit is the state involved in the ISC step. Other low-lying triplet states that have( π , π * ) ${(\pi, \pi ^{\ast} )}$ and( n O , π * ) ${(n_{\rm{O}}, \pi ^{\ast} )}$ character are ruled out on the basis of energetics and/or their orbital character. The mechanistic picture that emerges from the calculations is one that involves the ISC of photoexcited flavin to a triplet( n N , π * ) ${(n_{\rm{N}}, \pi ^{\ast} )}$ state followed by rapid internal conversion to a triplet( π , π * ) ${(\pi, \pi ^{\ast} )}$ state, which is the state detected spectroscopically. This insight into the ISC mechanism can provide guidelines for tuning flavin's photophysics through mutations that alter the protein electrostatic environment and potentially helps to explain why ISC (and subsequent flavin photochemistry) does not occur readily in many classes of flavin-binding enzymes.

Conical Intersections Studied by the Configuration-Interaction-Corrected Tamm-Dancoff Method

Conical intersections directly mediate the internal energy conversion in photoinduced processes in a wide range of chemical and biological systems. Because of the Brillouin theorem, many conventional electronic structure methods, including configuration interaction with single excitations from a Hartree-Fock reference and time-dependent density functional theory in either the linear response approximation (TDDFT) or Tamm-Dancoff approximation (DFT-TDA), have the wrong dimensionality for conical intersections between the ground state (S0) and the first excited state (S1) of the same multiplicity. This leads to unphysical state crossings. Here, we implement and assess the configuration-interaction-corrected Tamm-Dancoff approximation (CIC-TDA) that restores the correct dimensionality of conical intersections by including the coupling between the reference state and the intersecting excited state. We apply the CIC-TDA method to the S1/S0 conical intersections in ammonia (NH3), ethylene (C2H4), bithiophene (C8H6S2), azobenzene (C12H10N2), and 11-cis retinal protonated Schiff base (PSB11) in vacuo. We show that this black-box approach can produce potential energy surfaces (PESs) of comparable accuracy to multireference wave function methods. The method validated here can allow cost-efficient explorations of photoinduced electronically nonadiabatic dynamics, especially for large molecules and complex systems.

On the Meaning of De-Excitations in Time-Dependent Density Functional Theory Computations

De-excitations play a prominent role within the mathematical formalism of time-dependent density functional theory (TDDFT) and other excited-state response methods. However, their physical meaning remains largely unexplored and poorly understood. It is the purpose of this work to shed new light on this issue. The main thesis developed here is that de-excitations are not a peculiarity of TDDFT but that they are a more fundamental property of the underlying wave functions reflecting how electrons are excited between partially occupied orbitals. The paraquinodimethane (pQDM) molecule is chosen as a convenient model system whose open-shell character can be modulated via twisting of its methylene groups. Using the one-electron transition density matrix as a rigorous basis for our analysis, we highlight qualitative and quantitative parallels in the way that de-excitations are reflected in multireference wave function and TDDFT computations. As a physically observable consequence, we highlight a lowering of the transition dipole moment that derives from destructive interference between the excitation and de-excitation contributions. In summary, we hope that this work will shed new light on formal and practical aspects regarding the application of TDDFT to excited-state computations, especially of diradicaloid systems.

Analytical nuclear gradient and derivative coupling theories for multireference perturbation methods

Electron correlations should be appropriately included in quantum chemistry calculations to accurately describe the energy and wave functions. In multiconfigurational methods, the reference functions are written as linear combinations of multiple electronic configurations to describe static correlations. Using the multiconfigurational reference functions, it is also possible to correct for dynamical correlations using various methods. Geometry optimizations and dynamics simulations are among the most prominent applications of quantum chemistry methods. Such applications become much more straightforward when analytical nuclear gradients are available. Many efficient algorithms for computing analytical nuclear gradients and derivative coupling using multireference perturbation theories (MRPTs) have recently been developed. This work aims to provide a comprehensive and easy-to-follow review of analytical gradient theories and the properties of methods for obtaining analytical gradients and derivative coupling methods using MRPTs. We also briefly review the practical applications of these methods in performing nonadiabatic dynamics simulations.

Spin-Orbit Coupling and Admixture Coefficients in SA-CASSCF and MS-CASPT2, and Triplet Excitation Yield Simulated via Trajectory Surface Hopping and Calibrated SA-CASSCF in 1,2-Dioxetane Derivatives

The energy gaps, spin-orbit coupling (SOC), and admixture coefficients over a series of the configurations are evaluated by the SA-CASSCF/6-31G, SA-CASSCF/6-31G*, SA-CASSCF/ANO-RCC-VDZP, and MS-CASPT2/ANO-RCC-VDZP to reveal the extent of the inaccuracy of the SA-CASSCF. By comparing the mean absolute errors for the energy gaps and the admixture coefficient magnitudes (ACMs) measured between the SA-CASSCF/6-31G, SA-CASSCF/6-31G*, or SA-CASSCF/ANO-RCC-VDZP and the MS-CASPT2/ANO-RCC-VDZP, the SA-CASSCF/6-31G is selected as the electronic structure method in the nonadiabatic molecular dynamics simulation. The major components of the ACMs of the SA-CASSCF/6-31G and MS-CASPT2/ANO-RCC-VDZP are identified and compared; we find that the ACMs are underestimated by the SA-CASSCF/6-31G, which is verified by the reasonable triplet quantum yield simulated by the trajectory surface hopping and the calibrated SA-CASSCF/6-31G. The magnitude of the singlet-triplet mixing positively correlates to the hopping probability between the mixed singlet and triplet states, which is confirmed by the computed S-T transition probability.

Analytic First-Order Derivatives of CASPT2 Combined with the Polarizable Continuum Model

The complete active space second-order perturbation theory (CASPT2) is valuable for accurately predicting electronic structures and transition energies. However, optimizing molecular geometries in the solution phase has proven challenging. In this study, we develop analytic first-order derivatives of CASPT2 using an implicit solvation model, specifically the polarizable continuum model, within the open-source package OpenMolcas. Analytic gradients and nonadiabatic coupling vectors are computed by solving a modified Z-vector equation. Comparisons with existing theoretical and experimental results demonstrate that the solvent effects can be qualitatively captured using the developed method.

OpenQP: A Quantum Chemical Platform Featuring MRSF-TDDFT with an Emphasis on Open-Source Ecosystem

The OpenQP (Open Quantum Platform) is a new open-source quantum chemistry library developed to tackle sustainability and interoperability challenges in the field of computational chemistry. OpenQP provides various popular quantum chemical theories as autonomous modules such as energy and gradient calculations of HF, DFT, TDDFT, SF-TDDFT, and MRSF-TDDFT, thereby allowing easy interconnection with third-party software. A scientifically notable feature is the innovative mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) and its customized exchange-correlation functionals such as the DTCAM series of VAEE, XI, XIV, AEE, and VEE, which significantly expand the applicability scope of DFT and TDDFT. OpenQP also supports parallel execution and is optimized with BLAS and LAPACK for high performance. Future enhancements such as extended Koopman's theorem (EKT)-MRSF-TDDFT and spin-orbit coupling (SOC)-MRSF-TDDFT will further expand OpenQP's capabilities. Additionally, a Python wrapper PyOQP is provided that performs tasks such as geometry optimization, conical intersection searches, and nonadiabatic coupling calculations, among others, by prototyping the modules of the OpenQP library in combination with third-party libraries. Overall, OpenQP aligns with modern trends in high-performance scientific software development by offering flexible prototyping and operation while retaining the performance benefits of compiled languages like Fortran and C. They enhance the sustainability and interoperability of quantum chemical software, making OpenQP a crucial platform for accelerating the development of advanced quantum theories like MRSF-TDDFT.

UMRSF-TDDFT: Unrestricted Mixed-Reference Spin-Flip-TDDFT

An unrestricted version of Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory (UMRSF-TDDFT) was developed based on unrestricted Kohn-Sham orbitals (UKS) with a new molecular orbital (MO) reordering scheme. Additionally, a simple yet accurate method for estimating ⟨S2⟩ expectation values was devised. UMRSF-TDDFT was benchmarked against cases where DFT, TDDFT, and SF-TDDFT traditionally fail to provide accurate descriptions. In an application to the ground and excited states of a Be atom, UMRSF-TDDFT successfully recovers the degenerate states, with its energies slightly reduced compared to its RO counterpart, due to the additional variational flexibility of UKS. A clear difference between UMRSF and U-SF-TDDFT is evident in the bond breaking of the hydrogen fluoride system, as the latter misses an important configuration. In the case of the Jahn-Teller distortion of trimethylenemethane (TMM), the relative singlet energy compared to the triplet is lower by 0.1 and 0.2 eV for UMRSF and U-SF-TDDFT, respectively, than that of MRSF-TDDFT. The reduction in UMRSF energy is attributed to spatial orbital relaxations, whereas the reduction in U-SF-TDDFT energy results from spin contamination. Overall, the additional orbital relaxations afforded by unrestricted Kohn-Sham (UKS) orbitals in UMRSF-TDDFT lead to lower total system energies compared to their restricted open-shell counterparts. This enhancement adds a practical and accurate quantum chemical theory to the existing RO variant for addressing challenging systems where traditional quantum theories suffer.

Excited state electronic structure of dimethyl disulfide involved in photodissociation at ∼200 nm

Dimethyl disulfide (DMDS), one of the smallest organic molecules with an S-S bond, serves as a model system for understanding photofragmentation in polypeptides and proteins. Prior studies of DMDS photodissociation excited at ∼266 nm and ∼248 nm have elucidated the mechanisms of S-S and C-S bond cleavage, which involve the lowest excited electronic states S1 and S2. Far less is known about the dissociation mechanisms and electronic structure of relevant excited states of DMDS excited at ∼200 nm. Herein we present calculations of the electronic structure and properties of electronic states S1-S6 accessed when DMDS is excited at ∼200 nm. Our analysis includes a comparison of theoretical and experimental UV spectra, as well as theoretically predicted one-dimensional cuts through the singlet and triplet potential energy surfaces along the S-S and C-S bond dissociation coordinates. Finally, we present calculations of spin-orbit coupling constants at the Franck-Condon geometry to assess the likelihood of ultrafast intersystem crossing. We show that choosing an accurate yet computationally efficient electronic structure method for calculating the S0-S6 potential energy surfaces along relevant dissociation coordinates is challenging due to excited states with doubly excited character and/or mixed Rydberg-valence character. Our findings demonstrate that the extended multi-state complete active space second-order perturbation theory (XMS-CASPT2) balances this computational efficiency and accuracy, as it captures both the Rydberg character of states in the Franck-Condon region and multiconfigurational character toward the bond-dissociation limits. We compare the performance of XMS-CASPT2 to a new variant of equation of motion coupled cluster theory with single, double, and perturbative triple corrections, EOM-CCSD(T)(a)*, finding that EOM-CCSD(T)(a)* significantly improves the treatment of doubly excited states compared to EOM-CCSD, but struggles to quantitatively capture asymptotic energies along bond dissociation coordinates for these states.

Electronic Structure Methods for Simulating Flavin's Spectroscopy and Photophysics: Comparison of Multi-reference, TD-DFT, and Single-Reference Wave Function Methods

The use of flavins and flavoproteins in photocatalytic, sensing, and biotechnological applications has led to a growing interest in computationally modeling the excited-state electronic structure and photophysics of flavin. However, there is limited consensus regarding which computational methods are appropriate for modeling flavin's photophysics. We compare the energies of low-lying excited states of flavin computed with time-dependent density functional theory (TD-DFT), equation-of-motion coupled cluster (EOM-EE-CCSD), scaled opposite-spin configuration interaction [SOS-CIS(D)], multiconfiguration pair-density functional theory (MC-PDFT), and several multireference perturbation theory (MR-PT2) methods. In the first part, we focus on excitation energies of the first singlet excited state (S1) of five different redox and protonation states of flavin, with the goal of finding a suitable active space for MR-PT2 calculations. In the second part, we construct two sets of one-dimensional potential energy surfaces connecting the S0 and S1 equilibrium geometries (S0-S1 path) and the S1 (π,π*) and S2 (n,π*) equilibrium geometries (S1-S2 path). The first path therefore follows a Franck-Condon active mode of flavin while the second path maps crossings points between low-lying singlet and triplet states in flavin. We discuss the similarities and differences in the TD-DFT, EOM-EE-CCSD, SOS-CIS(D), MC-PDFT and MR-PT2 energy profiles along these paths. We find that (TD-)DFT methods are suitable for applications such as simulating the spectra of flavins but are inconsistent with several other methods when used for some geometry optimizations and when describing the energetics of dark (n,π*) states. MR-PT2 methods show promise for the simulation of flavin's low-lying excited states, but the selection of orbitals for the active space and the number of roots used for state averaging must be done carefully to avoid artifacts. Some properties, such as the intersystem crossing geometry and energy between the S1 (π,π*) and T2 (n,π*) states, may require additional benchmarking before they can be determined quantitatively.

On the Topological Phase around Conical Intersections with Tamm-Dancoff Linear-Response Time-Dependent Density Functional Theory

Regions of nuclear-configuration space away from the Franck-Condon geometry can prove problematic for some electronic structure methods, given the propensity of such regions to possess conical intersections, i.e., (highly connected) points of degeneracy between potential energy surfaces. With the likelihood (perhaps even inevitability) for nonadiabatic dynamics simulations to explore molecular geometries in close proximity to conical intersections, it is vital that the performance of electronic structure methods is routinely examined in this context. In a recent paper [Taylor, J. T. J. Chem. Phys. 2023, 159, 214115.], the ability of linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) to provide a proper description of conical intersections, in terms of their topology and topography, was investigated, with particular attention paid to conical intersections between two excited electronic states. For the same prototypical molecules, protonated formaldimine and pyrazine, we herein consider whether AA LR-TDDFT can correctly reproduce the topological phase accumulated by the adiabatic electronic wave function upon traversing a closed path around an excited-to-excited state conical intersection despite not using the appropriate quadratic-response nonadiabatic coupling vectors. Equally, we probe the ability of the ground-to-excited state intersection ring exhibited by AA LR-TDDFT in protonated formaldimine to give rise to a similar topological phase in spite of its incorrect dimensionality.

Multiconfigurational Correlation at DFT + U Cost: On-Site Electron-Electron Interactions Yield a Block-Localized Configuration Interaction Hamiltonian

This work presents a first-principles wavefunction-in-DFT approach based on the Hubbard density functional theory (DFT) + U method. This approach begins with the standard DFT reference system of noninteracting electrons and introduces an electron-electron interaction projected onto DFT+U-type atomic states. The reference system's configuration interaction Hamiltonian is block-localized to these states and can be expressed in terms of state occupation numbers, state self-energies (which correspond to unscreened Hubbard U values), and the promotion energies of doubly excited Slater determinants. Simple approximations for the promotion energies provide multiconfigurational correlation energies without requiring explicit orbital localization/transform. Numerical results for fractionally occupied chromium atom, bonded chromium dimer, dissociating covalent bonds, and large active spaces show that the approach provides beyond-zero-sum accuracy at computational cost comparable to standard DFT+U.

Quantum simulation of conical intersections

We explore the simulation of conical intersections (CIs) on quantum devices, setting the groundwork for potential applications in nonadiabatic quantum dynamics within molecular systems. The intersecting potential energy surfaces of H3+ are computed from a variance-based contracted quantum eigensolver. We show how the CIs can be correctly described on quantum devices using wavefunctions generated by the anti-Hermitian contracted Schrödinger equation ansatz, which is a unitary transformation of wavefunctions that preserves the topography of CIs. A hybrid quantum-classical procedure is used to locate the seam of CIs. Additionally, we discuss the quantum implementation of the adiabatic to diabatic transformation and its relation to the geometric phase effect. Results on noisy intermediate-scale quantum devices showcase the potential of quantum computers in dealing with problems in nonadiabatic chemistry.

Minimum-Energy Conical Intersections by Compressed Multistate Pair-Density Functional Theory

Compressed multistate pair-density functional theory (CMS-PDFT) is a multistate version of multiconfiguration pair-density functional theory that can capture the correct topology of coupled potential energy surfaces (PESs) around conical intersections. In this work, we develop interstate coupling vectors (ISCs) for CMS-PDFT in the OpenMolcas and PySCF/mrh electronic structure packages. Yet, the main focus of this work is using ISCs to calculate minimum-energy conical intersections (MECIs) by CMS-PDFT. This is performed using the projected constrained optimization method in OpenMolcas, which uses ISCs to restrain the iterations to the conical intersection seam. We optimize the S1/S0 MECIs for ethylene, butadiene, and benzene and show that CMS-PDFT gives smooth PESs in the vicinities of the MECIs. Furthermore, the CMS-PDFT MECIs are in good agreement with the MECI calculated by the more expensive XMS-CASPT2 method.

Signatures of Conical Intersections in Extreme Ultraviolet Photoelectron Spectra of Furan Measured with 15 fs Time Resolution

Ultrafast internal conversion of furan upon deep UV excitation at 200 nm is studied by using extreme ultraviolet time-resolved photoelectron spectroscopy with a time resolution of 15 fs. Ballistic nuclear wavepacket motion from the 1B2(ππ*) state to the ground state is fully observed using 21.7 eV probe pulses. Through the performance of a comparison with the results of electronic structure calculations at the MS(3)-CASPT2(10,10)/cc-pVTZ level of theory, the photoelectron signals from the conical intersection regions are identified.

Which Electronic Structure Method to Choose in Trajectory Surface Hopping Dynamics Simulations? Azomethane as a Case Study

Nonadiabatic dynamics simulations have become a standard approach to explore photochemical reactions. Such simulations require underlying potential energy surfaces and couplings between them, calculated at a chosen level of theory, yet this aspect is rarely assessed. Here, in combination with the popular trajectory surface hopping dynamics method, we use a high-accuracy XMS-CASPT2 electronic structure level as a benchmark for assessing the performances of various post-Hartree-Fock methods (namely, CIS, ADC(2), CC2, and CASSCF) and exchange-correlation functionals (PBE, PBE0, and CAM-B3LYP) in a TD-DFT/TDA context, using the isomerization around a double bond as test case. Different relaxation pathways are identified, and the ability of the different methods to reproduce their relative importance and time scale is discussed. The results show that multireference electronic structure methods should be preferred, when studying nonadiabatic decay between excited and ground states. If not affordable, TD-DFT with TDA and hybrid functionals and ADC(2) are efficient alternatives but overestimate the nonradiative decay yield and thus may miss deexcitation pathways.

Visible light-induced palladium-carbon bond weakening in catalytically relevant T-shaped complexes

Triggering one-electron redox processes during palladium catalysis holds the potential to unlock new reaction mechanisms and synthetic methods not previously accessible in the typical two-electron reaction manifolds that dominate palladium catalysis. We report that T-shaped organopalladium(ii) complexes coordinated by a bulky monophosphine, a class of organometallic intermediate featured in a range of contemporary catalytic reactions, undergo blue light-promoted bond weakening leading to mild and efficient homolytic cleavage of strong Pd(ii)-C(sp3) bonds under ambient conditions. The origin of light-triggered radical formation in these systems, which lack an obvious ligand-based chromophore (i.e., π-systems), was investigated using a combination of DFT calculations, photoactinometry, and transient absorption spectroscopy. The available data suggest T-shaped organopalladium(ii) complexes manifest unusual blue light-accessible Pd-to-C(sp3) transition. The quantum efficiency and excited state lifetime of this process were unexpectedly superior compared to a prototypical (α-diimine)Pd(ii) complex featuring a low-lying, ligand-centered LUMO (π*). These results suggest coordinatively-unsaturated organopalladium(ii) compounds, catalysts in myriad catalytic processes, have untapped potential for one-electron reactivity under visible light excitation.

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.

Assessment of the Electron Correlation Treatment on the Quantum-Classical Dynamics of Retinal Protonated Schiff Base Models: XMS-CASPT2, RMS-CASPT2, and REKS Methods

We compare the performance of three different multiconfigurational wave function-based electronic structure methods and two implementations of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method. The study is characterized by three features: (i) it uses a small set of quantum-classical trajectories rather than potential energy surface mapping, (ii) it focuses, exclusively, on the photoisomerization of retinal protonated Schiff base models, and (iii) it probes the effect of both methyl substitution and the increase in length of the conjugate π-system. For each tested method, the corresponding analytical gradients are used to drive the quantum-classical (Tully's FSSH method) trajectory propagation, including the recent multistate XMS-CASPT2 and RMS-CASPT2 gradients. It is shown that while CASSCF, XMS-CASPT2, and RMS-CASPT2 yield consistent photoisomerization dynamics descriptions, REKS produces, in some of these systems, qualitatively different behavior that is attributed to a flatter and topographically different excited state potential energy surface. The origin of this behavior can be traced back to the effect of the employed density functional approximation. The above studies are further expanded by benchmarking, at the CASSCF and REKS levels, the electronic structure methods using a QM/MM model of the visual pigment rhodopsin.

Doubly Tuned Exchange-Correlation Functionals for Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory

It is demonstrated that significant accuracy improvements in MRSF-TDDFT can be achieved by introducing two different exchange-correlation (XC) functionals for the reference Kohn-Sham DFT and the response part of the calculations, respectively. Accordingly, two new XC functionals of doubly tuned Coulomb attenuated method-vertical excitation energy (DTCAM-VEE) and DTCAM-AEE were developed on the basis of the "adaptive exact exchange (AEE)" concept in the framework of the Coulomb-attenuating XC functionals. The values by DTCAM-VEE are in excellent agreement with those of Thiel's set [mean absolute errors (MAEs) and the interquartile range (IQR) values of 0.218 and 0.327 eV, respectively]. On the other hand, DTCAM-AEE faithfully reproduced the qualitative aspects of conical intersections (CIs) of trans-butadiene and thymine and the nonadiabatic molecular dynamics (NAMD) simulations on thymine. The latter functional also remarkably exhibited the exact 1/R asymptotic behavior of the charge-transfer state of an ethylene-tetrafluoroethylene dimer and the accurate potential energy surfaces (PESs) along the two torsional angles of retinal protonated Schiff base model with six double bonds (rPSB6). Overall, DTCAM-AEE generally performs well, as its MAE (0.237) and IQR (0.41 eV) are much improved as compared to BH&HLYP. The current idea can also be applied to other XC functionals as well as other variants of linear response theories, opening a new way of developing XC functionals.

Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory: Multireference Advantages with the Practicality of Linear Response Theory

The density functional theory (DFT) and linear response (LR) time-dependent (TD)-DFT are of the utmost importance for routine computations. However, the single reference formulation of DFT suffers in the description of open-shell singlet systems such as diradicals and bond-breaking. LR-TDDFT, on the other hand, finds difficulties in the modeling of conical intersections, doubly excited states, and core-level excitations. In this Perspective, we demonstrate that many of these limitations can be overcome by recently developed mixed-reference (MR) spin-flip (SF)-TDDFT, providing an alternative yet accurate route for such challenging situations. Empowered by the practicality of the LR formalism, it is anticipated that MRSF-TDDFT can become one of the major workhorses for general routine tasks.

Organic compounds for solid state luminescence enhancement/aggregation induced emission: a theoretical perspective

Organic luminophores displaying one or more forms of luminescence enhancement in solid state are extremely promising for the development and performance optimization of functional materials essential to many modern key technologies. Yet, the effort to harness their huge potential is riddled with hurdles that ultimately come down to a limited understanding of the interactions that result in the diverse molecular environments responsible for the macroscopic response. In this context, the benefits of a theoretical framework able to provide mechanistic explanations to observations, supported by quantitative predictions of the phenomenon, are rather apparent. In this perspective, we review some of the established facts and recent developments about the current theoretical understanding of solid-state luminescence enhancement (SLE) with an accent on aggregation-induced emission (AIE). A description of the macroscopic phenomenon and the questions it raises is accompanied by a discussion of the approaches and quantum chemistry methods that are more apt to model these molecular systems with the inclusion of an accurate yet efficient simulation of the local environment. A sketch of a general framework, building from the current available knowledge, is then attempted via the analysis of a few varied SLE/AIE molecular systems from literature. A number of fundamental elements are identified offering the basis for outlining design rules for molecular architectures exhibiting SLE that involve specific structural features with the double role of modulating the optical response of the luminophores and defining the environment they experience in solid state.

High-performance strategies for the recent MRSF-TDDFT in GAMESS

Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. This paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios.

Analytic first-order derivatives of CASPT2 with IPEA shift

Complete active space second-order perturbation theory (CASPT2) is useful for accurately predicting properties of complex electronic structures, but it is well known that it systematically underestimates excitation energies. The underestimation can be corrected using the ionization potential-electron affinity (IPEA) shift. In this study, analytic first-order derivatives of CASPT2 with the IPEA shift are developed. CASPT2-IPEA is not invariant with respect to rotations among active molecular orbitals, and two additional constraint conditions are necessary in the CASPT2 Lagrangian to formulate analytic derivatives. The method developed here is applied to methylpyrimidine derivatives and cytosine, and minimum energy structures and conical intersections are located. By comparing energies relative to the closed-shell ground state, we find that the agreement with experiments and high-level calculations is indeed improved by the inclusion of the IPEA shift. The agreement of geometrical parameters with high-level calculations may also be improved in some cases.

Exact non-adiabatic coupling vectors for the time-dependent density functional based tight-binding method

We report on non-adiabatic coupling vectors between electronic excited states for the time-dependent-density functional theory based tight-binding (TD-DFTB) method. The implementation includes orbital relaxation effects that have been previously neglected and covers also the case of range-separated exchange-correlation functionals. Benchmark calculations with respect to first principles TD-DFT highlight the large dependence of non-adiabatic couplings on the functional. Closer investigations of the topology around a conical intersection between excited states show that TD-DFTB delivers near-exact values of the Berry phase, which paves the way for consistent non-adiabatic molecular dynamics simulations for large systems.

Assessment of State-Averaged Driven Similarity Renormalization Group on Vertical Excitation Energies: Optimal Flow Parameters and Applications to Nucleobases

We present a comprehensive excited-state benchmark for the state-averaged (SA) driven similarity renormalization group (DSRG) [Li, C.; Evangelista, F. A. J. Chem. Phys. 2018, 148, 124106]. Following the QUEST database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; Loos, P.-F. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2021, 11, e1517], 280 vertical transition energies of 35 medium-sized molecules are computed using the SA-DSRG derived second- and third-order perturbation theories (PT2/PT3) along with a nonperturbative approach [sq-LDSRG(2)]. Comparing to the theoretical best estimates, the optimal flow parameter is found to be 0.35 and 2.0 Eh-2 for SA-DSRG-PT2 and SA-DSRG-PT3, respectively. For SA-sq-LDSRG(2), a flow parameter of 1.5 Eh-2 provides converged equations without compromising the accuracy. We then assess the accuracy of the SA-DSRG hierarchy using these parameters. The SA-DSRG-PT2 scheme outperforms the level-shifted CASPT2 by 0.10 eV in mean absolute error (MAE), yet this accuracy is slightly inferior than that of CASPT2 with the ionization-potential-electron-affinity shift. Both SA-DSRG-PT3 and SA-sq-LDSRG(2) yield a MAE of 0.10 eV, which is comparable to that of CASPT3 (0.09 eV). Finally, we compute vertical excitation energies of several low-lying singlet states of nucleobases. The SA-sq-LDSRG(2) approach provides highly accurate results for π → π* excitations, while n → π* transitions are better described by SA-DSRG-PT3.

Good Vibrations: Calculating Excited-State Frequencies Using Ground-State Self-Consistent Field Models

The use of Δ-self-consistent field (SCF) approaches for studying excited electronic states has received a renewed interest in recent years. In this work, the use of this scheme for calculating excited-state vibrational frequencies is examined. Results from Δ-SCF calculations for a set of representative molecules are compared with those obtained using configuration interaction with single substitutions (CIS) and time-dependent density functional theory (TD-DFT) methods. The use of an approximate spin purification model is also considered for cases where the excited-state SCF solution is spin-contaminated. The results of this work demonstrate that an SCF-based description of an excited-state potential energy surface can be an accurate and cost-effective alternative to CIS and TD-DFT methods.

Conical Intersection in Chemiluminescence of Cyclic Peroxides

Chemiluminescence (CL) utilizing chemiexcitation for energy transformation is one of the most highly sensitive and useful analytical techniques. The chemiexcitation is a chemical process of a ground-state reactant producing an excited-state product, in which a nonadiabatic event is facilitated by conical intersections (CIs), the specific molecular geometries where electronic states are degenerated. Cyclic peroxides, especially 1,2-dioxetane/dioxetanone derivatives, are the iconic chemiluminescent substances. In this Perspective, we concentrated on the CIs in the CL of cyclic peroxides. We first present a computational overview on the role of CIs between the ground (S₀) state and the lowest singlet excited (S₁) state in the thermolysis of cyclic peroxides. Subsequently, we discuss the role of the S₀/S₁ CI in the CL efficiency and point out misunderstandings in some theoretical studies on the singlet chemiexcitations of cyclic peroxides. Finally, we address the challenges and future prospects in theoretically calculating S₀/S₁ CIs and simulating the dynamics and chemiexcitation efficiency in the CL of cyclic peroxides.

Applying Generalized Variational Principles to Excited-State-Specific Complete Active Space Self-consistent Field Theory

We employ a generalized variational principle to improve the stability, reliability, and precision of fully excited-state-specific complete active space self-consistent field theory. Compared to previous approaches that similarly seek to tailor this ansatz's orbitals and configuration interaction expansion for an individual excited state, we find the present approach to be more resistant to root flipping and better at achieving tight convergence to an energy stationary point. Unlike state-averaging, this approach allows orbital shapes to be optimal for individual excited states, which is especially important for charge-transfer states and some doubly excited states. We demonstrate the convergence and state-targeting abilities of this method in LiH, ozone, and MgO, showing in the latter that it is capable of finding three excited-state energy stationary points that no previous method has been able to locate.

Ab initio trajectory surface-hopping dynamics studies of excited-state proton-coupled electron transfer reactions in trianisoleheptazine-phenol complexes

The excited-state proton-coupled electron-transfer (PCET) reaction in hydrogen-bonded complexes of trianisoleheptazine (TAHz), a chromophore related to polymeric carbon nitrides widely used in hydrogen-evolution photocatalysis, with several phenol derivatives were recently studied by Schlenker and coworkers with time-resolved photoluminescence quenching and pump-probe experiments. A pronounced dependence of the PCET reactivity on the electron-donating/electron-withdrawing character of the substituents on phenol was found, with indications of a barrierless or nearly barrierless PCET reaction for the most strongly electron-donating substituent, methoxy. In the present work, the excited-state PCET dynamics was explored with first-principles nonadiabatic dynamics simulations using the TDDFT/ωB97X-D electronic-structure model for two selected complexes, TAHz-phenol and TAHz-methoxyphenol. The qualitative reliability of the TDDFT/ωB97X-D electronic-structure model was assessed by extensive benchmarking of excitation energies and potential-energy profiles against a wave-function-based ab initio method, the algebraic-diagrammatic construction of second order (ADC(2)). The nonadiabatic dynamics simulations provide temporally and structurally resolved insights into paradigmatic PCET reactions in TAHz-phenol complexes. The radiationless relaxation of the photoexcited bright ¹ππ* state to the long-lived dark S₁ state of TAHz occurs in less than 100 fs. The ensuing PCET reaction on the adiabatic S₁ surface is faster in TAHz-methoxyphenol complexes than in TAHz-phenol complexes due to a lower H-atom-transfer barrier, as observed in the experiments. The relaxation of the complexes to the electronic ground state is found to occur exclusively via PCET within the 250 fs time window covered by the present simulations, confirming the essential role of the hydrogen bond for the fluorescence quenching process. The absolute values of the computed PCET time constants are significantly shorter than those extracted from time-resolved photoluminescence measurements for mixtures of TAHz with phenolic substrates in toluene. The possible origins of this discrepancy are discussed.

Analytic First-Order Derivatives of (X)MS, XDW, and RMS Variants of the CASPT2 and RASPT2 Methods

Crossings between states involve complex electronic structures, making the accurate characterization of the crossing point difficult. In this study, the analytic derivatives of three complete active space second-order perturbation theory (CASPT2) variants as well as an extension of the restricted active space (RASPT2) are developed. These variants are applied to locating minimum energy conical intersections. Our results demonstrate that the three CASPT2 variants predict qualitatively similar results, but a recently developed variant, the rotated multistate CASPT2 (RMS-CASPT2), is least sensitive to the number of states considered in the calculation. We demonstrate that CASPT2 and the reference self-consistent field calculations predict qualitatively different energetics and bond lengths.

InteraChem: Exploring Excited States in Virtual Reality with Ab Initio Interactive Molecular Dynamics

InteraChem is an ab initio interactive molecular dynamics (AI-IMD) visualizer that leverages recent advances in virtual reality hardware and software, as well as the graphical processing unit (GPU)-accelerated TeraChem electronic structure package, in order to render quantum chemistry in real time. We introduce the exploration of electronically excited states via AI-IMD using the floating occupation molecular orbital-complete active space configuration interaction method. The optimization tools in InteraChem enable identification of excited state minima as well as minimum energy conical intersections for further characterization of excited state chemistry in small- to medium-sized systems. We demonstrate that finite-temperature Hartree-Fock theory is an efficient method to perform ground state AI-IMD. InteraChem allows users to track electronic properties such as molecular orbitals and bond order in real time, resulting in an interactive visualization tool that aids in the interpretation of excited state chemistry data and makes quantum chemistry more accessible for both research and educational purposes.

Excited-state double proton transfer of 1,8-dihydroxy-2-naphthaldehyde: A MS-CASPT2//CASSCF study

Excited-state double proton transfer (ESDPT) is a controversial issue which has long been plagued with theoretical and experimental communities. Herein, we took 1,8-dihydroxy-2-naphthaldehyde (DHNA) as a prototype and used combined complete active space self-consistent field (CASSCF) and multi-state complete active-space second-order perturbation (MS-CASPT2) methods to investigate ES-DPT and excited-state deactivation pathways of DHNA. Three different tautomer minima of S1-ENOL, S1-KETO-1, and S1-KETO-2 and two crucial conical intersections of S1S0-KETO-1 and S1S0-KETO-2 in.and between the S0 and S1 states were obtained. S1-KETO-1 and S1-KETO-2 should take responsibility for experimentally observing dual-emission bands. In addition, two-dimensional potential energy surfaces (2D-PESs) and linear interpolated internal coordinate paths connecting relevant structures were calculated at the MS-CASPT2//CASSCF level and confirmed a stepwise ESDPT mechanism. Specifically, the first proton transfer from S1-ENOL to S1-KETO-1 is barrierless, whereas the second one from S1-KETO-1 to S1-KETO-2 demands a barrier of ca. 6.0 kcal/mol. The linear interpolated internal coordinate path connecting S1-KETO-1 (S1-KETO-2) and S1S0-KETO-1 (S1S0-KETO-2) is uphill with a barrier of ca. 12.0 kcal/mol, which will trap DHNA in the S1 state while therefore enabling dual-emission bands. On the other hand, the S1/S0 conical intersections would also prompt the S1 system to decay to the S0 state, which could be to certain extent suppressed by locking the rotation of the C5−C8−C9−O10 dihedral angle. These mechanistic insights are not only helpful for understanding ESDPT but also useful for designing novel molecular materials with excellent photoluminescent performances.

Generalized Formulation of the Density Functional Tight Binding-Based Restricted Ensemble Kohn-Sham Method with Onsite Correction to Long-Range Correction

We present a generalized formulation for the combination of the density functional tight binding (DFTB) approach and the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS or SSR) method by considering onsite correction (OC) as well as the long-range corrected (LC) functional. The OC contribution provides more accurate energies and analytic gradients for individual microstates, while the multireference character of the SSR provides the correct description for conical intersections. We benchmark the LC-OC-DFTB/SSR method against various DFTB calculation methods for excitation energies and conical intersection structures with π/π* or n/π* characters. Furthermore, we perform excited-state molecular dynamics simulations with a molecular rotary motor with variations of LC-OC-DFTB/SSR approaches. We show that the OC contribution to the LC functional is crucial to obtain the correct geometry of conical intersections.

Fewest switches surface hopping with Baeck-An couplings

In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.

Fewest switches surface hopping with Baeck-An couplings

In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.

A Conical Intersection Influences the Ground State Rearrangement of Fulvene to Benzene

The rearrangement of fulvene to benzene is believed to play an important role in the formation of soot during hydrocarbon combustion. Previous work has identified two possible mechanisms for the rearrangement─a unimolecular path and a hydrogen-atom-assisted, bimolecular path. Computational results to date have suggested that the unimolecular mechanism faces a barrier of about 74 kcal/mol, which makes it unable to compete with the bimolecular mechanism under typical combustion conditions. This computed barrier is about 10 kcal/mol higher than the experimental value, which is an unusually large discrepancy for modern electronic structure theory. In the present work, we have reinvestigated the unimolecular mechanism computationally, and we have found a second transition state that is approximately 10 kcal/mol lower in energy than the previously identified one and, therefore, in excellent agreement with the experimental value. The existence of two transition states for the same rearrangement arises because there is a conical intersection between the two lowest singlet states which occurs in the vicinity of the reaction coordinates. The two possible paths around the cone on the lower adiabatic surface give rise to the two distinct saddle points. The lower barrier for the unimolecular mechanism now makes it competitive with the bimolecular one, according to our calculations. In support of this conclusion, we have reanalyzed some previous experimental results on anisole pyrolysis, which leads to benzene as a significant product and have shown that the unimolecular and bimolecular mechanisms for fulvene → benzene must be occurring competitively in that system. Finally, we have identified that similar conical intersections arise during the isomerizations of benzofulvene and isobenzofulvene to naphthalene.

Analytical Gradient Theory for Spin-Free State-Averaged Second-Order Driven Similarity Renormalization Group Perturbation Theory (SA-DSRG-MRPT2) and Its Applications for Conical Intersection Optimizations

Second-order multireference-driven similarity renormalization group perturbation theory (DSRG-MRPT2) provides an efficient means of correcting the dynamical correlation with the multiconfiguration reference function. The state-averaged DSRG-MRPT2 (SA-DSRG-MRPT2) method is the simplest means of treating the excited states with DSRG-MRPT2. In this method, the Hamiltonian dressed with dynamical correlation is diagonalized in the CASCI state subspace (SA-DSRG-MRPT2c) or the configuration subspace (SA-DSRG-MRPT2). This work develops analytical gradient theory for spin-free SA-DSRG-MRPT2(c) with the density-fitting approximation. We check the accuracy of the analytical gradients against the numerical gradients. We present applications for optimizing minimum energy conical intersections (MECI) of ethylene and retinal model chromophores (PSB3 and RPSB6). We investigate the dependence of the optimized geometries and energies on the flow parameters and reference relaxations. The smoothness of the SA-DSRG-MRPT2(c) potential energy surfaces near the reference (complete active space self-consistent field) MECI is comparable to the XMCQDPT2 one. These results render SA-DSRG-MRPT2(c) theory a promising approach for studies of conical intersections.

Analytical Nonadiabatic Couplings and Gradients within the State-Averaged Orbital-Optimized Variational Quantum Eigensolver

We introduce several technical and analytical extensions to our recent state-averaged orbital-optimized variational quantum eigensolver (SA-OO-VQE) algorithm (see Yalouz et al. Quantum Sci. Technol. 2021, 6, 024004). Motivated by the limitations of current quantum computers, the first extension consists of an efficient state-resolution procedure to find the SA-OO-VQE eigenstates, and not just the subspace spanned by them, while remaining in the equi-ensemble framework. This approach avoids expensive intermediate resolutions of the eigenstates by postponing this problem to the very end of the full algorithm. The second extension allows for the estimation of analytical gradients and nonadiabatic couplings, which are crucial in many practical situations ranging from the search of conical intersections to the simulation of quantum dynamics, in, for example, photoisomerization reactions. The accuracy of our new implementations is demonstrated on the formaldimine molecule CH₂NH (a minimal Schiff base model relevant for the study of photoisomerization in larger biomolecules), for which we also perform a geometry optimization to locate a conical intersection between the ground and first-excited electronic states of the molecule.

Mechanistic photophysics and photochemistry of unnatural bases and sunscreen molecules: insights from electronic structure calculations

Photophysics and photochemistry are basic subjects in the study of light-matter interactions and are ubiquitous in diverse fields such as biology, energy, materials, and environment. A full understanding of mechanistic photophysics and photochemistry underpins many recent advances and applications. This contribution first provides a short discussion on the theoretical calculation methods we have used in relevant studies, then we introduce our latest progress on the mechanistic photophysics and photochemistry of two classes of molecular systems, namely unnatural bases and sunscreens. For unnatural bases, we disclose the intrinsic driving forces for the ultrafast population to reactive triplet states, impacts of the position and degree of chalcogen substitutions, and the effects of complex environments. For sunscreen molecules, we reveal the photoprotection mechanisms that dissipate excess photon energy to the surroundings by ultrafast internal conversion to the ground state. Finally, relevant theoretical challenges and outlooks are discussed.

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.

Analytic gradients for restricted active space second-order perturbation theory (RASPT2)

The computational cost of analytic derivatives in multireference perturbation theory is strongly affected by the size of the active space employed in the reference self-consistent field calculation. To overcome previous limits on the active space size, the analytic gradients of single-state restricted active space second-order perturbation theory (RASPT2) and its complete active space second-order perturbation theory (CASPT2) have been developed and implemented in a local version of OpenMolcas. Similar to previous implementations of CASPT2, the RASPT2 implementation employs the Lagrangian or Z-vector method. The numerical results show that restricted active spaces with up to 20 electrons in 20 orbitals can now be employed for geometry optimizations.

Fast and Accurate Computation of Nonadiabatic Coupling Matrix Elements Using the Truncated Leibniz Formula and Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory

We present a fast and accurate numerical algorithm for computing the first-order nonadiabatic coupling matrix element (NACME). The algorithm employs the truncated Leibniz formula (TLF) approximation within the finite-difference method, which makes it easily applicable in connection with any wave function-based methodology. In this work, we used the algorithm in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT, MRSF for brevity). The accuracy is assessed for NACME between the singlet electronic states of a dissociating hydrogen molecule. It is demonstrated that an intermediate approximation, TLF(1), affords a negligible numeric error on the order of ∼10^-10 a.u. while enabling a fast computation of NACME. As the MRSF method yields the correct description of the dissociation curves of H₂ for all the electronic states involved, the numeric TLF(1)/MRSF NACME values are in excellent agreement with the reference analytical values obtained by the full configuration interaction. For polyatomic molecules, the MRSF NAC vectors agree very closely with the MRCISD NAC vectors. Hence, the proposed protocol is a promising tool for the evaluation of NACMEs.

Role of Triplet States in the Photodynamics of Aniline

The dynamics of excited heteroaromatic molecules is a key to understanding the photoprotective properties of many biologically relevant chromophores that dissipate their excitation energy nonreactively and thereby prevent the detrimental effects of ultraviolet radiation. Despite their structural variability, most substituted aromatic compounds share a common feature of a repulsive ¹πσ* potential energy surface. This surface can lead to photoproducts, and it can also facilitate the population transfer back to the ground electronic state by means of a ¹πσ*/S₀ conical intersection. Here, we explore a hidden relaxation route involving the triplet electronic state of aniline, which has recently been discovered by means of time-selected photofragment translational spectroscopy [J. Chem. Phys. 2019, 151, 141101]. By using the recently available analytical gradients for multiconfiguration pair-density functional theory, it is now possible to locate the minimum-energy crossing points between states of different spin and therefore compute the intersystem crossing rates with a multireference method, rather than with the less reliable single-reference methods. Using such calculations, we demonstrate that the population loss of aniline in the T₁(³ππ*) state is dominated by C₆H₅NH₂ → C₆H₅NH· + H· dissociation, and we explain the long nonradiative lifetimes of the T₁(³ππ*) state at the excitation wavelengths of 294-264 nm.

Signatures of Conical Intersection Dynamics in the Time-Resolved Photoelectron Spectrum of Furan: Theoretical Modeling with an Ensemble Density Functional Theory Method

The non-adiabatic dynamics of furan excited in the ππ* state (S2 in the Franck-Condon geometry) was studied using non-adiabatic molecular dynamics simulations in connection with an ensemble density functional method. The time-resolved photoelectron spectra were theoretically simulated in a wide range of electron binding energies that covered the valence as well as the core electrons. The dynamics of the decay (rise) of the photoelectron signal were compared with the excited-state population dynamics. It was observed that the photoelectron signal decay parameters at certain electron binding energies displayed a good correlation with the events occurring during the excited-state dynamics. Thus, the time profile of the photoelectron intensity of the K-shell electrons of oxygen (decay constant of 34 ± 3 fs) showed a reasonable correlation with the time of passage through conical intersections with the ground state (47 ± 2 fs). The ground-state recovery constant of the photoelectron signal (121 ± 30 fs) was in good agreement with the theoretically obtained excited-state lifetime (93 ± 9 fs), as well as with the experimentally estimated recovery time constant (ca. 110 fs). Hence, it is proposed to complement the traditional TRPES observations with the trXPS (or trNEXAFS) measurements to obtain more reliable estimates of the most mechanistically important events during the excited-state dynamics.

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.

Equation-of-motion coupled-cluster method with double electron-attaching operators: Theory, implementation, and benchmarks

We report a production-level implementation of the equation-of-motion (EOM) coupled-cluster (CC) method with double electron-attaching (DEA) EOM operators of 2p and 3p1h types, EOM-DEA-CCSD. This ansatz, suitable for treating electronic structure patterns that can be described as two-electrons-in-many orbitals, represents a useful addition to the EOM-CC family of methods. We analyze the performance of EOM-DEA-CCSD for energy differences and molecular properties. By considering reduced quantities, such as state and transition one-particle density matrices, we compare EOM-DEA-CCSD wave functions with wave functions computed by other EOM-CCSD methods. The benchmarks illustrate that EOM-DEA-CCSD is capable of treating diradicals, bond-breaking, and some types of conical intersections.