Progress and Challenges in the Calculation of Electronic Excited States

The Analysis of Electron Densities: From Basics to Emergent Applications

The electron density determines all properties of a system of nuclei and electrons. It is both computable and observable. Its topology allows gaining insight into the mechanisms of bonding and other phenomena in a way that is complementary to and beyond that available from the molecular orbital picture and the formal oxidation state (FOS) formalism. The ability to derive mechanistic insight from electron density is also important with methods where orbitals are not available, such as orbital-free density functional theory (OF-DFT). While density topology-based analyses such as QTAIM (quantum theory of atoms-in-molecules) have been widely used, novel, vector-based techniques recently emerged such as next-generation (NG) QTAIM. Density-dependent quantities are also actively used in machine learning (ML)-based methods, in particular, for ML DFT functional development, including machine-learnt kinetic energy functionals. We review QTAIM and its recent extensions such as NG-QTAIM and localization-delocalization matrices (LDM) and their uses in the analysis of bonding, conformations, mechanisms of redox reactions excitations, as well as ultrafast phenomena. We review recent research showing that direct density analysis can circumvent certain pitfalls of the FOS formalism, in particular in the description of anionic redox, and of the widely used (spherically) projected density of states analysis. We discuss uses of density-based quantities for the construction of DFT functionals and prospects of applications of analyses of density topology to get mechanistic insight with OF-DFT and recently developed time-dependent OF-DFT.

Fully Polarizable Multiconfigurational Self-Consistent Field/Fluctuating Charges Approach

A multiscale model based on the coupling of the multiconfigurational self-consistent field (MCSCF) method and the classical atomistic polarizable fluctuating charges (FQ) force field is presented. The resulting MCSCF/FQ approach is validated by exploiting the CASSCF scheme through application to compute vertical excitation energies of formaldehyde and para-nitroaniline in aqueous solution. The procedure is integrated with molecular dynamics simulations to capture the solute's conformational changes and the dynamic aspects of solvation. Comparative analysis with alternative solvent models, gas-phase calculations, and experimental data provides insights into the model's accuracy in reproducing solute-solvent molecular interactions and spectral signals.

Subspace-Search Quantum Imaginary Time Evolution for Excited State Computations

Quantum systems in excited states are attracting significant interest with the advent of noisy intermediate-scale quantum (NISQ) devices. While ground states of small molecular systems are typically explored using hybrid variational algorithms like the variational quantum eigensolver (VQE), the study of excited states has received much less attention, partly due to the absence of efficient algorithms. In this work, we introduce the subspace search quantum imaginary time evolution (SSQITE) method, which calculates excited states using quantum devices by integrating key elements of the subspace search variational quantum eigensolver (SSVQE) and the variational quantum imaginary time evolution (VarQITE) method. The effectiveness of SSQITE is demonstrated through calculations of low-lying excited states of benchmark model systems including H2 and LiH molecules. A toy Hamiltonian is also employed to demonstrate that the robustness of VarQITE in avoiding local minima extends to its use in excited state algorithms. With this robustness in avoiding local minima, SSQITE shows promise for advancing quantum computations of excited states across a wide range of applications.

Rank-Reduced Equation-of-Motion Coupled Cluster Triples: an Accurate and Affordable Way of Calculating Electronic Excitation Energies

In the present work, we report an implementation of the rank-reduced equation-of-motion coupled cluster method with approximate triple excitations (RR-EOM-CC3). The proposed variant relies on tensor decomposition techniques in order to alleviate the high cost of computing and manipulating the triply excited amplitudes. In the RR-EOM-CC3 method, both ground-state and excited-state triple-excitation amplitudes are compressed according to the Tucker-3 format. This enables factorization of the working equations such that the formal scaling of the method is reduced to N6, where N is the system size. An additional advantage of our method is the fact that the accuracy can be strictly controlled by proper choice of two parameters defining sizes of triple-excitation subspaces in the Tucker decomposition for the ground and excited states. Optimal strategies of selecting these parameters are discussed. The developed method has been tested in a series of calculations of electronic excitation energies and compared to its canonical EOM-CC3 counterpart. Errors several times smaller than the inherent error of the canonical EOM-CC3 method (in comparison to FCI) are straightforward to achieve. This conclusion holds both for valence states dominated by single excitations and for states with pronounced doubly excited character. Taking advantage of the decreased scaling, we demonstrate substantial computational costs reductions (in comparison with the canonical EOM-CC3) in the case of two large molecules - l-proline and heptazine. This illustrates the usefulness of the RR-EOM-CC3 method for accurate determination of excitation energies of large molecules.

Assessment of the similarity-transformed equation of motion (STEOM) for open-shell organic and transition metal molecules

This study benchmarks the newly re-implemented single-reference excited-state methods, IP-EOM-CCSD, EA-EOM-CCSD, and STEOM-CCSD, in ORCA6.0, with a focus on open-shell systems. We compare STEOM against EOM-CCSD, CC3, and CCSDT across a range of systems, including small organic radicals, hydrated transition metal (TM) ions, and TM diatomic systems with both closed and open-shell configurations. For organic radicals, STEOM and EOM-CCSD show comparable performance, aligning closely with CC3 and CCSDT results. In the case of hydrated TM ions, IP-EOM closely matches DLPNO-CCSD results, while deviations from DLPNO-CCSD(T) are consistent. For open-shell TM systems, IP-EOM exhibits a blueshift relative to both the DLPNO-CCSD methods, while EA-EOM-CCSD shows better agreement. When comparing STEOM and CC3 to CCSDT, STEOM shows slightly larger deviations in closed-shell systems but shows excellent agreement in open-shell systems. Computational efficiency is also assessed, revealing a significant speedup in ORCA 6.0 compared to ORCA 5.0, with optimizations improving computation times. This study provides valuable insights into the performance and efficiency of STEOM in various chemical environments, highlighting its strengths and limitations.

A Constraint-Based Orbital-Optimized Excited State Method (COOX)

In this work, we present a novel method to directly calculate targeted electronic excited states within a self-consistent field calculation based on constrained density functional theory (cDFT). The constraint is constructed from the static occupied-occupied and virtual-virtual parts of the excited state difference density from (simplified) linear-response time-dependent density functional theory calculations (LR-TDDFT). Our new method shows a stable convergence behavior, provides an accurate excited state density adhering to the Aufbau principle, and can be solved within a restricted SCF for singlet excitations to avoid spin contamination. This also allows the straightforward application of post-SCF electron-correlation methods like MP2 or direct RPA methods. We present the details of our constraint-based orbital-optimized excited state method (COOX) and compare it to similar schemes. The accuracy of excitation energies will be analyzed for a benchmark of systems, while the quality of the resulting excited state densities is investigated by evaluating excited state nuclear forces and excited state structure optimizations. We also investigate the performance of the proposed COOX method for long-range charge transfer excitations and conical intersections with the ground-state.

Machine Learning-Assisted Mixed Quantum-Classical Dynamics without Explicit Nonadiabatic Coupling: Application to the Photodissociation of Peroxynitric Acid

We have devised a hybrid quantum-classical scheme utilizing machine-learned potential energy surfaces (PES), which circumvents the need for explicit computation of nonadiabatic coupling elements. The quantities necessary to account for the nonadiabatic effects are directly obtained from the PESs. The simulation of dynamics is based on the fewest-switches surface-hopping method. We applied this scheme to model the photodissociation of both N-O and O-O bonds in a conformer of peroxynitric acid (HO2NO2). Adiabatic PES data for the six lowest states of this molecule were computed at the CASSCF level for various nuclear configurations. These served as the training data for the machine-learning models for the PESs. The dynamics simulation was initiated on the lowest optically bright singlet excited state (S4) and propagated along the two Jacobi coordinates J → 1 and J → 2 while accounting for the nonadiabatic effects through transitions between PESs. Our analysis revealed that there is a very high chance of dissociation of the N-O bond leading to the HO2 and NO2 fragments.

Restricted open-shell time-dependent density functional theory with perturbative spin-orbit coupling

When using quantum chemical methods to study electronically excited states of open-shell molecules, it is often beneficial to start with wave functions that are spin eigenfunctions. For excited states of molecules containing heavy elements, spin-orbit coupling (SOC) is important and needs to be included as well. An efficient approach is to include SOC perturbatively on top of a restricted open-shell Kohn-Sham (ROKS) time-dependent density functional theory, which can be combined with the Tamm-Dancoff approximation (TDA) to suppress numerical instabilities. We implemented and assessed the potential of such a ROKS-TDA-SOC method, also featuring the possibility of calculating transition dipole moments between states to allow for full spectrum simulation. Our study shows that the ROKS-TDA-SOC formalism yields a clear and easy-to-use method to obtain electronically excited states of open-shell molecules that are of moderate size and contain heavy elements.

An Improved Penalty-Based Excited-State Variational Monte Carlo Approach with Deep-Learning Ansatzes

We introduce several improvements to the penalty-based variational quantum Monte Carlo (VMC) algorithm for computing electronic excited states of Entwistle et al. [Nat. Commun. 14, 274 (2023)] and demonstrate that the accuracy of the updated method is competitive with other available excited-state VMC approaches. A theoretical comparison of the computational aspects of these algorithms is presented, where several benefits of the penalty-based method are identified. Our main contributions include an automatic mechanism for tuning the scale of the penalty terms, an updated form of the overlap penalty with proven convergence properties, and a new term that penalizes the spin of the wave function, enabling the selective computation of states with a given spin. With these improvements, along with the use of the latest self-attention-based ansatz, the penalty-based method achieves a mean absolute error below 1 kcal/mol for the vertical excitation energies of a set of 26 atoms and molecules, without relying on variance matching schemes. Considering excited states along the dissociation of the carbon dimer, the accuracy of the penalty-based method is on par with that of natural-excited-state (NES) VMC, while also providing results for additional sections of the potential energy surface, which were inaccessible with the NES method. Additionally, the accuracy of the penalty-based method is improved for a conical intersection of ethylene, with the predicted angle of the intersection agreeing well with both NES-VMC and multireference configuration interaction.

Perspective on Theoretical and Experimental Advances in Atmospheric Photochemistry

Research that explores the chemistry of Earth's atmosphere is central to the current understanding of global challenges such as climate change, stratospheric ozone depletion, and poor air quality in urban areas. This research is a synergistic combination of three established domains: earth observation, for example, using satellites, and in situ field measurements; computer modeling of the atmosphere and its chemistry; and laboratory measurements of the properties and reactivity of gas-phase molecules and aerosol particles. The complexity of the interconnected chemical and photochemical reactions which determine the composition of the atmosphere challenges the capacity of laboratory studies to provide the spectroscopic, photochemical, and kinetic data required for computer models. Here, we consider whether predictions from computational chemistry using modern electronic structure theory and nonadiabatic dynamics simulations are becoming sufficiently accurate to supplement quantitative laboratory data for wavelength-dependent absorption cross-sections, photochemical quantum yields, and reaction rate coefficients. Drawing on presentations and discussions from the CECAM workshop on Theoretical and Experimental Advances in Atmospheric Photochemistry held in March 2024, we describe key concepts in the theory of photochemistry, survey the state-of-the-art in computational photochemistry methods, and compare their capabilities with modern experimental laboratory techniques. From such considerations, we offer a perspective on the scope of computational (photo)chemistry methods based on rigorous electronic structure theory to become a fourth core domain of research in atmospheric chemistry.

In Silico Chemical Experiments in the Age of AI: From Quantum Chemistry to Machine Learning and Back

Computational chemistry is an indispensable tool for understanding molecules and predicting chemical properties. However, traditional computational methods face significant challenges due to the difficulty of solving the Schrödinger equations and the increasing computational cost with the size of the molecular system. In response, there has been a surge of interest in leveraging artificial intelligence (AI) and machine learning (ML) techniques to in silico experiments. Integrating AI and ML into computational chemistry increases the scalability and speed of the exploration of chemical space. However, challenges remain, particularly regarding the reproducibility and transferability of ML models. This review highlights the evolution of ML in learning from, complementing, or replacing traditional computational chemistry for energy and property predictions. Starting from models trained entirely on numerical data, a journey set forth toward the ideal model incorporating or learning the physical laws of quantum mechanics. This paper also reviews existing computational methods and ML models and their intertwining, outlines a roadmap for future research, and identifies areas for improvement and innovation. Ultimately, the goal is to develop AI architectures capable of predicting accurate and transferable solutions to the Schrödinger equation, thereby revolutionizing in silico experiments within chemistry and materials science.

Excited State-Specific CASSCF Theory for the Torsion of Ethylene

State-specific complete active space self-consistent field (SS-CASSCF) theory has emerged as a promising route to accurately predict electronically excited energy surfaces away from molecular equilibria. However, its accuracy and practicality for chemical systems of photochemical interest have yet to be fully determined. We investigate the performance of the SS-CASSCF theory for the low-lying ground and excited states in the double bond rotation of ethylene. We show that state-specific approximations with a minimal (2e,2o) active space provide comparable accuracy to state-averaged calculations with much larger active spaces, while optimizing the orbitals for each excited state significantly improves the spatial diffusivity of the wave function. However, the incorrect ordering of state-specific solutions causes excited state solutions to coalesce and disappear, creating unphysical discontinuities in the potential energy surface. Our findings highlight the theoretical challenges that must be overcome to realize practical applications of state-specific electronic structure theory for computational photochemistry.

Mechanisms of photoisomerization of the prenylated flavin mononucleotide cofactor: a theoretical study

The enzymatic decarboxylation of α,β-unsaturated acids using the ferulic acid decarboxylase (Fdc1) enzyme and prenylated flavin mononucleotide (prFMN) cofactor is a potential, environmentally friendly reaction for the biosynthesis of styrene and its derivatives. However, experiments showed that the enzyme activity of Fdc1 depends on the ring structure of prFMN, namely, the iminium and ketimine forms, and the loss of enzyme activity results from prFMNim → prFMNket photoisomerization. To obtain insight into this photochemical process and to improve the enzyme efficiency of Fdc1, two proposed photoisomerization mechanisms with different proton sources for the acid-base reaction were studied herein using theoretical methods. The potential energy surfaces calculated using the density functional theory method with the Becke, 3-parameter, and Lee-Yang-Parr hybrid functionals and DZP basis set (DFT/B3LYP/DZP) and TD-DFT/B3LYP/DZP methods confirmed that the light-dependent reaction occurs in the rate-determining proton transfer process and that the mechanism involving intermolecular proton transfer between prFMNim and Glu282 (external base) is energetically more favorable than that involving intramolecular proton transfer in prFMNim (internal base). The thermodynamic results obtained from the transition state theory method suggested that the exothermic relaxation energy in the photo-to-thermal process can promote the spontaneous formation of a high-energy-barrier transition state, and an effective enzymatic decarboxylation could be achieved by slowing down the formation of the undesirable thermodynamically favorable product (prFMNket). Because the rate constant for formation of the high-energy-barrier transition state varies exponentially over the temperature range of 273-298 K, and experimental results have shown that incubating Fdc1 on ice results in a complete loss of enzyme activity, it is recommended to perform the decarboxylation reaction at 285 K to strike a balance between minimizing enzyme stability loss at 273 K and mitigating the effects of UV irradiation. The computational strategy and fundamental insights obtained in this study could serve as guidelines for future theoretical and experimental investigations on the same and similar photochemical systems.

State-Specific Coupled-Cluster Methods for Excited States

We reexamine ΔCCSD, a state-specific coupled-cluster (CC) with single and double excitations (CCSD) approach that targets excited states through the utilization of non-Aufbau determinants. This methodology is particularly efficient when dealing with doubly excited states, a domain in which the standard equation-of-motion CCSD (EOM-CCSD) formalism falls short. Our goal here to evaluate the effectiveness of ΔCCSD when applied to other types of excited states, comparing its consistency and accuracy with EOM-CCSD. To this end, we report a benchmark on excitation energies computed with the ΔCCSD and EOM-CCSD methods for a set of molecular excited-state energies that encompasses not only doubly excited states but also doublet-doublet transitions and (singlet and triplet) singly excited states of closed-shell systems. In the latter case, we rely on a minimalist version of multireference CC known as the two-determinant CCSD method to compute the excited states. Our data set, consisting of 276 excited states stemming from the quest database [Véril et al., WIREs Comput. Mol. Sci. 2021, 11, e1517], provides a significant base to draw general conclusions concerning the accuracy of ΔCCSD. Except for the doubly excited states, we found that ΔCCSD underperforms EOM-CCSD. For doublet-doublet transitions, the difference between the mean absolute errors (MAEs) of the two methodologies (of 0.10 and 0.07 eV) is less pronounced than that obtained for singly excited states of closed-shell systems (MAEs of 0.15 and 0.08 eV). This discrepancy is largely attributed to a greater number of excited states in the latter set exhibiting multiconfigurational characters, which are more challenging for ΔCCSD. We also found typically small improvements by employing state-specific optimized orbitals.

Many-Body Excited States with a Contracted Quantum Eigensolver

Calculating ground and excited states is an exciting prospect for near-term quantum computing applications, and accurate and efficient algorithms are needed to assess viable directions. We develop an excited-state approach based on the contracted quantum eigensolver (ES-CQE), which iteratively attempts to find a solution to a contraction of the Schrödinger equation projected onto a subspace and does not require a priori information on the system. We focus on the anti-Hermitian portion of the equation, leading to a two-body unitary ansatz. We investigate the role of symmetries, initial states, constraints, and overall performance within the context of the model strongly correlated rectangular H4 system. We show that the ES-CQE achieves near-exact accuracy across the majority of states, covering regions of strong and weak electron correlation, while also elucidating challenging instances for two-body unitary ansatz.

Predicting the photodynamics of cyclobutanone triggered by a laser pulse at 200 nm and its MeV-UED signals-A trajectory surface hopping and XMS-CASPT2 perspective

This work is part of a prediction challenge that invited theoretical/computational chemists to predict the photochemistry of cyclobutanone in the gas phase, excited at 200 nm by a laser pulse, and the expected signal that will be recorded during a time-resolved megaelectronvolt ultrafast electron diffraction (MeV-UED). We present here our theoretical predictions based on a combination of trajectory surface hopping with XMS-CASPT2 (for the nonadiabatic molecular dynamics) and Born-Oppenheimer molecular dynamics with MP2 (for the athermal ground-state dynamics following internal conversion), coined (NA+BO)MD. The initial conditions were sampled from Born-Oppenheimer molecular dynamics coupled to a quantum thermostat. Our simulations indicate that the main photoproducts after 2 ps of dynamics are CO + cyclopropane (50%), CO + propene (10%), and ethene and ketene (34%). The photoexcited cyclobutanone in its second excited electronic state S2 can follow two pathways for its nonradiative decay: (i) a ring-opening in S2 and a subsequent rapid decay to the ground electronic state, where the photoproducts are formed, or (ii) a transfer through a closed-ring conical intersection to S1, where cyclobutanone ring opens and then funnels to the ground state. Lifetimes for the photoproduct and electronic populations were determined. We calculated a stationary MeV-UED signal [difference pair distribution function-ΔPDF(r)] for each (interpolated) pathway as well as a time-resolved signal [ΔPDF(r,t) and ΔI/I(s,t)] for the full swarm of (NA+BO)MD trajectories. Furthermore, our analysis provides time-independent basis functions that can be used to fit the time-dependent experimental UED signals [both ΔPDF(r,t) and ΔI/I(s,t)] and potentially recover the population of photoproducts. We also offer a detailed analysis of the limitations of our model and their potential impact on the predicted experimental signals.

Is this how bromine spiropyran salt is converted to merocyanine under UV irradiation? A look through the prism of quantum chemical calculations

The stability of merocyanine forms formed under UV irradiation of a solution of a spiropyran salt, in which an organic part acts as a cation and a compact bromide ion as an anion, their photophysical properties, and the formation mechanism are studied in this work using time-dependent density functional theory. Theoretical calculations show that TTC and CTT are the most stable open forms (the difference in stability energies is 10.5 and 12.0 kcal mol-1 relative to the formation energy of spiropyran, respectively). The simulated absorption bands in the UV spectrum of the merocyanine forms are observed both in the UV region at 308-366 nm and in the visible region at 544-757 nm due to n → π* and π → π* type transitions. We found that the isomerisation mechanism of spiropyran into merocyanine forms includes two key stages: the ring opening to form cisoid merocyanine forms (except unstable TCC) through conical intersection and their subsequent isomerisation to form stable transoid isomers. The length of the Cspiro-O bond is 1.97 Å and the C1'-C2'-C3'-C4' angle is 70° in the structure close to conical intersection. The stage that determines the rate of this process is the isomerisation between transoid forms, as in the case of transformation of open merocyanine forms into spiropyran.

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.

A Local Diabatisation Method for Two-State Adiabatic Conical Intersections

A methodology to locally characterize conical intersections (CIs) between two adiabatic electronic states for which no nonadiabatic coupling (NAC) vectors are available is presented. Based on the Hessian and gradient at the CI, the branching space coordinates are identified. The potential energy surface around the CI in the branching space is expressed in the diabatic representation, from which the NAC vectors can be calculated in a wave-function-free, energy-based approach. To demonstrate the universality of the developed methodology, the minimum-energy CI (MECI) between the first (S1) and second (S2) singlet excited states of formamide is investigated at the state-averaged complete active space self-consistent field (SA-CASSCF) and extended multistate complete active space second-order perturbation theory (XMS-CASPT2) levels of theory. In addition, the asymmetrical MECI between the ground state (S0) and S1 of cyclopropanone is evaluated using SA-CASSCF, as well as (ME)CIs between the S1 and S2 states of benzene using SA-CASSCF and time-dependent density functional theory (TDDFT). Finally, a CI between the S1 and S2 excited states of thiophene was analyzed using TDDFT.

The energies and charge and spin distributions in the low-lying levels of singlet and triplet N2V defects in diamond from direct variational calculations of the excited states

This paper reports the energies and charge and spin distributions of the low-lying excited states in singlet and triplet N2V defects in diamond from direct Δ-SCF calculations based on Gaussian orbitals within the B3LYP, PBE0, and HSE06 functionals. They assign the observed absorption at 2.463 eV, first reported by Davies et al. [Proc. R. Soc. London 351, 245 (1976)], to the excitation of a N(sp3) lone-pair electron in the singlet and triplet states, respectively, with estimates of ∼1.1 eV for that of the unpaired electrons, C(sp3). In both cases, the excited states are predicted to be highly local and strongly excitonic with 81% of the C(sp3) and 87% of the N(sp3) excited charges localized at the three C atoms nearest neighbor (nn) to the excitation sites. Also reported are the higher excited gap states of both the N lone pair and C unpaired electron. Calculated excitation energies of the bonding sp3 hybrids of the C atoms nn to the four inner atoms are close to that of the bulk, which indicates that the N2V defect is largely a local defect. The present results are in broad agreement with those reported by Udvarhelyi et al. [Phys. Rev. B 96, 155211 (2017)] from plane wave HSE06 calculations, notably for the N lone pair excitation energy, for which both predict an energy of ∼2.7 eV but with a difference of ∼0.5 eV for the excitation of the unpaired electron.

Natural Products Dereplication: Databases and Analytical Methods

The development of efficient methods for dereplication has been critical in the re-emergence of the research in natural products as a source of drug leads. Current dereplication workflows rapidly identify already known bioactive secondary metabolites in the early stages of any drug discovery screening campaign based on natural extracts or enriched fractions. Two main factors have driven the evolution of natural products dereplication over the last decades. First, the availability of both commercial and public large databases of natural products containing the key annotations against which the biological and chemical data derived from the studied sample are searched for. Second, the considerable improvement achieved in analytical technologies (including instrumentation and software tools) employed to obtain robust and precise chemical information (particularly spectroscopic signatures) on the compounds present in the bioactive natural product samples. This chapter describes the main methods of dereplication, which rely on the combined use of large natural product databases and spectral libraries, alongside the information obtained from chromatographic, UV-Vis, MS, and NMR spectroscopic analyses of the samples of interest.

What Controls the Quality of Photodynamical Simulations? Electronic Structure Versus Nonadiabatic Algorithm

The field of nonadiabatic dynamics has matured over the last decade with a range of algorithms and electronic structure methods available at the moment. While the community currently focuses more on developing and benchmarking new nonadiabatic dynamics algorithms, the underlying electronic structure controls the outcome of nonadiabatic simulations. Yet, the electronic-structure sensitivity analysis is typically neglected. In this work, we present a sensitivity analysis of the nonadiabatic dynamics of cyclopropanone to electronic structure methods and nonadiabatic dynamics algorithms. In particular, we compare wave function-based CASSCF, FOMO-CASCI, MS- and XMS-CASPT2, density-functional REKS, and semiempirical MRCI-OM3 electronic structure methods with the Landau-Zener surface hopping, fewest switches surface hopping, and ab initio multiple spawning with informed stochastic selection algorithms. The results clearly demonstrate that the electronic structure choice significantly influences the accuracy of nonadiabatic dynamics for cyclopropanone even when the potential energy surfaces exhibit qualitative and quantitative similarities. Thus, selecting the electronic structure solely on the basis of the mapping of potential energy surfaces can be misleading. Conversely, we observe no discernible differences in the performance of the nonadiabatic dynamics algorithms across the various methods. Based on the above results, we discuss the present-day practice in computational photodynamics.

Variational Active Space Selection with Multiconfiguration Pair-Density Functional Theory

The selection of an adequate set of active orbitals for modeling strongly correlated electronic states is difficult to automate because it is highly dependent on the states and molecule of interest. Although many approaches have shown some success, no single approach has worked well in all cases. In light of this, we present the "discrete variational selection" (DVS) approach to active space selection, in which one generates multiple trial wave functions from a diverse set of systematically constructed active spaces and then selects between these wave functions variationally. We apply this DVS approach to 207 vertical excitations of small-to-medium-sized organic and inorganic molecules (with 3 to 18 atoms) in the QUESTDB database by (i) constructing various sets of active space orbitals through diagonalization of parametrized operators and (ii) choosing the result with the lowest average energy among the states of interest. This approach proves ineffective when variationally selecting between wave functions using the density matrix renormalization group (DMRG) or complete active space self-consistent field (CASSCF) energy but is able to provide good results when variationally selecting between wave functions using the energy of the translated PBE (tPBE) functional from multiconfiguration pair-density functional theory (MC-PDFT). Applying this DVS-tPBE approach to selection among state-averaged DMRG wave functions, we obtain a mean unsigned error of only 0.17 eV using hybrid MC-PDFT. This result matches that of our previous benchmark without the need to filter out poor active spaces and with no further orbital optimization following active space selection of the SA-DMRG wave functions. Furthermore, we find that DVS-tPBE is able to robustly and effectively select between the new SA-DMRG wave functions and our previous SA-CASSCF results.

Linearized Pair-Density Functional Theory for Vertical Excitation Energies

Multiconfiguration pair-density functional theory (MC-PDFT) is a computationally efficient method that computes the energies of electronic states in a state specific or state average framework via an on-top functional. However, MC-PDFT does not include state interaction among these states since the final energies do not come from the diagonalization of an effective model-space Hamiltonian. Recently, multistate extensions such as linearized PDFT (L-PDFT) have been developed to accurately model the potentials near conical intersections and avoided crossings. However, there has not been any systematic study evaluating their performance for predicting vertical excitations at the equilibrium geometry of a molecule, when the excited states are generally well separated. In this paper, we report the performance of L-PDFT on the extensive QUESTDB data set of vertical excitations using a database of automatically selected active spaces. We show that L-PDFT performs well on all these excitations and successfully reproduces the performance of MC-PDFT. These results further demonstrate the potential of L-PDFT, as its scaling is constant with the number of states included in the state-average manifold, whereas MC-PDFT scales linearly in this regard.

Quantification of the Ionic Character of Multiconfigurational Wave Functions: The Qat Diagnostic

The complete active space self-consistent field (CASSCF) method is a cornerstone in modern excited-state quantum chemistry providing the starting point for most common multireference computations. However, CASSCF, when used with a minimal active space, can produce significant errors (>2 eV) even for the excitation energies of simple hydrocarbons if the states of interest possess ionic character. After illustrating this problem in some detail, we present a diagnostic for ionic character, denoted as Q at, that is readily computed from the transition density. A set of 11 molecules is considered to study errors in vertical excitation energies. State-averaged CASSCF obtains a mean absolute error (MAE) of 0.87 eV for the 34 singlet states considered. We highlight a strong correlation between the obtained errors and the Q at diagnostic, illustrating its power to predict problematic cases. Conversely, using multireference configuration interaction with single and double excitations and Pople's size extensivity correction (MR-CISD+P), excellent results are obtained with an MAE of 0.11 eV. Furthermore, correlations with the Q at diagnostic disappear. In summary, we hope that the presented diagnostic will facilitate reliable and user-friendly multireference computations on conjugated organic molecules.

Reference Vertical Excitation Energies for Transition Metal Compounds

To enrich and enhance the diversity of the quest database of highly accurate excitation energies [Véril, M.; et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], we report vertical transition energies in transition metal compounds. Eleven diatomic molecules with a singlet or doublet ground state containing a fourth-row transition metal (CuCl, CuF, CuH, ScF, ScH, ScO, ScS, TiN, ZnH, ZnO, and ZnS) are considered, and the corresponding excitation energies are computed using high-level coupled-cluster (CC) methods, namely, CC3, CCSDT, CC4, and CCSDTQ, as well as multiconfigurational methods such as CASPT2 and NEVPT2. In many cases, to provide more comprehensive benchmark data, we also provide full configuration interaction estimates computed with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. Based on these calculations, theoretical best estimates of the transition energies are established in both the aug-cc-pVDZ and aug-cc-pVTZ basis sets. This allows us to accurately assess the performance of the CC and multiconfigurational methods for this specific set of challenging transitions. Furthermore, comparisons with experimental data and previous theoretical results are also reported.

On the Prediction of Spectroscopic Fingerprints of Co2+ Complexes Relevant for the ZIF Nucleation Process

The nucleation process of zeolitic imidazolate frameworks (ZIFs) is to date not completely understood. Recently, it has been found that, during the formation of Co-ZIF-67, after mixing imidazole-type ligands with octahedral precursors containing oxygen-coordinated ligands, a metal-organic pool with a diversity of transition metal complexes (TMCs) is formed showing fingerprints of octahedral and tetrahedral Co2+ complexes with both types of ligands [Filez, M. Cell Rep. Phys. Sci. 2021, 2, 100680]. In order to further unravel this process, we aim to characterize the d-d transitions of the TMCs and focus on their number, intensity, and position, which change during the process and can thus serve as a fingerprint for the formed species. It was previously shown that the number of ligands and symmetry has a detrimental influence on the ground state properties of Co2+ TMCs. Herein, we investigate how far excited state properties of TMCs relevant during nanoporous formation processes can be predicted by time-dependent density functional theory (TDDFT) and ligand field density functional theory (LFDFT). As TMCs are known to be challenging systems with possibly degenerate ground states and double excitations, we first investigate the performance of both techniques on first-row octahedral aqua-complexes. With this knowledge, we then focus on tetrahedral Co2+ complexes with aqua and imidazole-type ligands in order to investigate in how far we can propose a spectroscopic fingerprint that allows us to follow the Co2+ complexes during the formation of Co-ZIF-67. The results of TDDFT and LFDFT are qualitatively in agreement and provide complementary information. We found that various features can be used to distinguish between the species. However, as LFDFT is not suited for TMCs possessing the extended imidazole-type ligands and double and spin-flip states are not included in TDDFT, both techniques need to be complemented with more advanced methods to obtain complete insight into the d-d excitations of TMCs with imidazole ligands. Therefore, we particularly explored ab initio ligand field theory, which is capable of describing double excitations and is, in contrast to LFDFT, suitable for TMCs with extended ligands.

Stochastic Real-Time Second-Order Green's Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green's function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(Ne3), where Ne is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

The three kingdoms-Photoinduced electron transfer cascades controlled by electronic couplings

Excited states are the key species in photocatalysis, while the critical parameters that govern their applications are (i) excitation energy, (ii) accessibility, and (iii) lifetime. However, in molecular transition metal-based photosensitizers, there is a design tension between the creation of long-lived excited (triplet), e.g., metal-to-ligand charge transfer (3MLCT) states and the population of such states. Long-lived triplet states have low spin-orbit coupling (SOC) and hence their population is low. Thus, a long-lived triplet state can be populated but inefficiently. If the SOC is increased, the triplet state population efficiency is improved-coming at the cost of decreasing the lifetime. A promising strategy to isolate the triplet excited state away from the metal after intersystem crossing (ISC) involves the combination of transition metal complex and an organic donor/acceptor group. Here, we elucidate the excited state branching processes in a series of Ru(II)-terpyridyl push-pull triads by quantum chemical simulations. Scalar-relativistic time-dependent density theory simulations reveal that efficient ISC takes place along 1/3MLCT gateway states. Subsequently, competitive electron transfer (ET) pathways involving the organic chromophore, i.e., 10-methylphenothiazinyl and the terpyridyl ligands are available. The kinetics of the underlying ET processes were investigated within the semiclassical Marcus picture and along efficient internal reaction coordinates that connect the respective photoredox intermediates. The key parameter that governs the population transfer away from the metal toward the organic chromophore either by means of ligand-to-ligand (3LLCT; weakly coupled) or intra-ligand charge transfer (3ILCT; strongly coupled) states was determined to be the magnitude of the involved electronic coupling.

Benchmarking Aspects of Ab Initio Fragment Models for Accurate Excimer Potential Energy Surfaces

While Coupled-Cluster methods have been proven to provide an accurate description of excited electronic states, the scaling of the computational costs with the system size limits the degree for which these methods can be applied. In this work different aspects of fragment-based approaches are studied on noncovalently bound molecular complexes with interacting chromophores of the fragments, such as π-stacked nucleobases. The interaction of the fragments is considered at two distinct steps. First, the states localized on the fragments are described in the presence of the other fragment(s); for this we test two approaches. One method is founded on QM/MM principles, only including the electrostatic interaction between the fragments in the electronic structure calculation with Pauli repulsion and dispersion effects added separately. The other model, a Projection-based Embedding (PbE) using the Huzinaga equation, includes both electrostatic and Pauli repulsion and only needs to be augmented by dispersion interactions. In both schemes the extended Effective Fragment Potential (EFP2) method of Gordon et al. was found to provide an adequate correction for the missing terms. In the second step, the interaction of the localized chromophores is modeled for a proper description of the excitonic coupling. Here the inclusion of purely electrostatic contributions appears to be sufficient: it is found that the Coulomb part of the coupling provides accurate splitting of the energies of interacting chromophores that are separated by more than 4 Å.

F12+EOM Quartic Force Fields for Rovibrational Predictions of Electronically Excited States

Quartic force fields (QFFs) constructed using a sum of ground-state CCSD(T)-F12b energies with EOM-CCSD excitation energies are proposed for computation of spectroscopic properties of electronically excited states. This is dubbed the F12+EOM approach and is shown to provide similar accuracy to previous methodologies at lower computational cost. Using explicitly correlated F12 approaches instead of canonical CCSD(T), as in the corresponding (T)+EOM approach, allows for 70-fold improvement in computational time. The mean percent difference between the two methods for anharmonic vibrational frequencies is only 0.10%. A similar approach is also developed herein which accounts for core correlation and scalar relativistic effects, named F12cCR+EOM. The F12+EOM and F12cCR+EOM approaches both match to within 2.5% mean absolute error of experimental fundamental frequencies. These new methods should help in clarifying astronomical spectra by assigning features to vibronic and vibrational transitions of small astromolecules when such data are not available experimentally.

A theoretical study on the photochemical generation of phenylborylene from phenyldiazidoborane

Organic borylenes are a kind of highly reactive species, which play important roles in a lot of reactions as vigorous intermediates. In this work, we investigated the photochemical generation mechanisms of phenylborylene (PhB) together with the side product N-phenylnitrenoiminoborane (PhNBN) from phenyldiazidoborane (PhBN₆) by extrusion of dinitrogen in the two lowest electronic singlet states (S₀ and S₁) based on the complete active space self-consistent field (CASSCF) and its second-order perturbation (CASPT2) methods combined with time-dependent density functional theory (TD-DFT) calculations. Our results show that the reaction PhBN₆ → PhB + 3N₂ involves stepwise N₂ extrusion three times and the azido region rearrangement. Moreover, we found that the studied photo-induced processes are kinetically feasible because the highest energy barrier is only 0.36 eV and excitation with light of wavelength 254 nm can provide enough excess energy to overcome these energy barriers. Importantly, we revealed that several conical intersections between S₁ and S₀ states participate and facilitate the studied photochemical processes. Our results not only clarify the experimental observations, (H. F. Bettinger, J. Am. Chem. Soc. 2006, 128, 2534), but also provide valuable insights into borylene chemistry.

Target State Optimized Density Functional Theory for Electronic Excited and Diabatic States

A flexible self-consistent field method, called target state optimization (TSO), is presented for exploring electronic excited configurations and localized diabatic states. The key idea is to partition molecular orbitals into different subspaces according to the excitation or localization pattern for a target state. Because of the orbital-subspace constraint, orbitals belonging to different subspaces do not mix. Furthermore, the determinant wave function for such excited or diabatic configurations can be variationally optimized as a ground state procedure, unlike conventional ΔSCF methods, without the possibility of collapsing back to the ground state or other lower-energy configurations. The TSO method can be applied both in Hartree-Fock theory and in Kohn-Sham density functional theory (DFT). The density projection procedure and the working equations for implementing the TSO method are described along with several illustrative applications. For valence excited states of organic compounds, it was found that the computed excitation energies from TSO-DFT and time-dependent density functional theory (TD-DFT) are of similar quality with average errors of 0.5 and 0.4 eV, respectively. For core excitation, doubly excited states and charge-transfer states, the performance of TSO-DFT is clearly superior to that from conventional TD-DFT calculations. It is shown that variationally optimized charge-localized diabatic states can be defined using TSO-DFT in energy decomposition analysis to gain both qualitative and quantitative insights on intermolecular interactions. Alternatively, the variational diabatic states may be used in molecular dynamics simulation of charge transfer processes. The TSO method can also be used to define basis states in multistate density functional theory for excited states through nonorthogonal state interaction calculations. The software implementing TSO-DFT can be accessed from the authors.

Machine Learning Models for Predicting Molecular UV-Vis Spectra with Quantum Mechanical Properties

Accurate understanding of ultraviolet-visible (UV-vis) spectra is critical for the high-throughput synthesis of compounds for drug discovery. Experimentally determining UV-vis spectra can become expensive when dealing with a large quantity of novel compounds. This provides us an opportunity to drive computational advances in molecular property predictions using quantum mechanics and machine learning methods. In this work, we use both quantum mechanically (QM) predicted and experimentally measured UV-vis spectra as input to devise four different machine learning architectures, UVvis-SchNet, UVvis-DTNN, UVvis-Transformer, and UVvis-MPNN, and assess the performance of each method. We find that the UVvis-MPNN model outperforms the other models when using optimized 3D coordinates and QM predicted spectra as input features. This model has the highest performance for predicting UV-vis spectra with a training RMSE of 0.06 and validation RMSE of 0.08. Most importantly, our model can be used for the challenging task of predicting differences in the UV-vis spectral signatures of regioisomers.

Photodynamics of the Molecular Ruby [Cr(ddpd)2]3

The introduction of strong-field ligands can enable luminescence in first-row transition-metal complexes. In this way, earth-abundant near-infrared emitters can be obtained using early 3d metals. A prime example is the molecular ruby [Cr(ddpd)₂]³⁺ (ddpd = N,N'-dimethyl-N,N'-dipyridin-2-ylpyridine-2,6-diamine) that can achieve high phosphorescence quantum yields at room temperature in aqueous solution. To understand these remarkable properties, here, we simulate its photodynamics in water using trajectory surface hopping on linear vibronic coupling potentials parametrized from multiconfigurational CASSCF/CASPT2 calculations. We find that after excitation to the second absorption band, a relaxation cascade through metal-centered states occurs. After an initial back-and-forth intersystem crossing with higher-lying doublet states, the complex relaxes through a manifold of quartet metal-centered states to the low-lying doublet metal-centered states which are responsible for the experimentally observed emission. These electronic processes are driven by an elongation of the Cr-ligand bond lengths as well as the twisting motion of the trans-coordinated pyridine units in the ddpd ligands. The low-lying doublet states are reached within 1-2 ps and are close in geometry to the doublet minima, thus explaining the high phosphorescence quantum yield of the molecular ruby [Cr(ddpd)₂]³⁺.

UV-vis spectroscopy of gas-phase ions

Photodissociation action spectroscopy has made a great progress in expanding investigations of gas-phase ion structures. This review deals with aspects of gas-phase ion electronic excitations that result in wavelength-dependent dissociation and light emission via fluorescence, chiefly covering the ultraviolet and visible regions of the spectrum. The principles are briefly outlined and a few examples of instrumentation are presented. The main thrust of the review is to collect and selectively present applications of UV-vis action spectroscopy to studies of stable gas-phase ion structures and combinations of spectroscopy with ion mobility, collision-induced dissociation, and ion-ion reactions leading to the generation of reactive intermediates and electronic energy transfer.

Power of the Sine Hamiltonian Operator for Estimating the Eigenstate Energies on Quantum Computers

Quantum computers have been shown to have tremendous potential in solving difficult problems in quantum chemistry. In this paper, we propose a new classical-quantum hybrid method, named as power of sine Hamiltonian operator (PSHO), to evaluate the eigenvalues of a given Hamiltonian (Ĥ). In PSHO, for any reference state |φ0⟩, the normalized energy of the sinn(H^τ)|φ0⟩ state can be determined. With the increase of the power, the initial reference state can converge to the eigenstate with the largest |sin(Eiτ)| value in the coefficients of the expansion of |φ0⟩, and the normalized energy of the sinn(H^τ)|φ0⟩ state converges to E_i. The ground- and excited-state energies of a Hamiltonian can be determined by taking different τ values. The performance of the PSHO method is demonstrated by numerical calculations of the H₄ and LiH molecules. Compared with the current popular variational quantum eigensolver method, PSHO does not need to design the ansatz circuits and avoids the complex nonlinear optimization problems. PSHO has great application potential in near-term quantum devices.

Large-Scale Benchmarking of Multireference Vertical-Excitation Calculations via Automated Active-Space Selection

We have calculated state-averaged complete-active-space self-consistent-field (SA-CASSCF), multiconfiguration pair-density functional theory (MC-PDFT), hybrid MC-PDFT (HMC-PDFT), and n-electron valence state second-order perturbation theory (NEVPT2) excitation energies with the approximate pair coefficient (APC) automated active-space selection scheme for the QUESTDB benchmark database of 542 vertical excitation energies. We eliminated poor active spaces (20-40% of calculations) by applying a threshold to the SA-CASSCF absolute error. With the remaining calculations, we find that NEVPT2 performance is significantly impacted by the size of the basis set the wave functions are converged in, regardless of the quality of their description, which is a problem absent in MC-PDFT. Additionally, we find that HMC-PDFT is a significant improvement over MC-PDFT with the translated PBE (tPBE) density functional and that it performs about as well as NEVPT2 and second-order coupled cluster on a set of 373 excitations in the QUESTDB database. We optimized the percentage of SA-CASSCF energy to include in HMC-PDFT when using the tPBE on-top functional, and we find the 25% value used in tPBE0 to be optimal. This work is by far the largest benchmarking of MC-PDFT and HMC-PDFT to date, and the data produced in this work are useful as a validation of HMC-PDFT and of the APC active-space selection scheme. We have made all the wave functions produced in this work (orbitals and CI vectors) available to the public and encourage the community to utilize this data as a tool in the development of further multireference model chemistries.

Rationalizing the Unexpected Sensitivity in Excited State Lifetimes of Adenine to Tautomerization by Nonadiabatic Molecular Dynamics

The remarkable photostability of canonical nucleobases makes them ideal building blocks for DNA and RNA. Even minor structural changes are expected to lead to drastic alteration of their subpicosecond excited state lifetimes. However, it is interesting to note that while the 9H- and 7H-amino tautomers of adenine possess drastically different lifetimes, 9H- and 7H-keto guanine possess similar excited state lifetimes. With an aim to explain this unexpected difference in sensitivity of lifetimes to tautomerization, we have investigated the excited state relaxation mechanism of UV-excited adenine and guanine tautomers using surface hopping based nonadiabatic molecular dynamics. We find that internal conversion in both guanine tautomers is almost barrierless while both adenine tautomers encounter significant barriers before they can deactivate. Moreover, the major deactivation channel (C2-puckering) in 9H-amino adenine is overall more efficient than the one (C6-puckering) in the 7H-amino form. We trace this difference to the frequent rotation of the amino group which disrupts its conjugation with the heterocyclic ring thereby reducing the strength of nonadiabatic coupling and, hence, delaying internal conversion.

Accurate Calculation of Excited-State Absorption for Small-to-Medium-Sized Conjugated Oligomers: Multiconfigurational Treatment vs Quadratic Response TD-DFT

Excited-state absorption (ESA) spectra of π-conjugated compounds are frequently calculated by (quadratic response) time-dependent density functional theory, (QR) TD-DFT, often giving a reasonable representation of the experimental results despite the (known) incomplete electronic description. To investigate whether this is inherent to the method, we calculate here the ESA spectra of small-to-medium-sized oligophenylenevinylenes (nPV) and oligothiophenes (nT) using QR TD-DFT as well as CASPT2 based on CASSCF geometries. CASPT2 gives indeed a reliable, theoretically correct description of the ESA features for all compounds; the computational effort can be reduced without significant loss of accuracy using TD-DFT geometries. QR TD-DFT, based on BHandHLYP and CAM-/B3LYP functionals, fails on short nTs but provides a reasonable description for spectral positions of nPVs and long nTs. The failure on short nTs is, however, only partly due to the incomplete configuration description but, in particular, related to an improper MO description, resulting in an asymmetric energy spacing of the occupied vs unoccupied MOs in the DFT scheme. Longer nTs, on the other side, adapt approximately the MO scheme for alternant hydrocarbons just like in nPVs, while contributions by two triplet excitations combined to a singlet (which inhibits an accurate treatment of polyenes with standard TD-DFT) do not play a relevant role in the current case. For such "well-behaved" systems, a reasonable representation of ESA spectra is found at the QR TD-DFT level due to the rather small energy shifts when including higher-order excitations.

Nucleic Acids as a Playground for the Computational Study of the Photophysics and Photochemistry of Multichromophore Assemblies

ConspectusThe interaction between light and multichromophoric assemblies (MCAs) is the primary event of many fundamental processes, from photosynthesis to organic photovoltaics, and it triggers dynamical processes that share remarkable similarities at the molecular scale: light absorption, energy and charge transfer, internal conversions, emission, and so on. Those events often involve many chromophores and different excited electronic states that are coupled on an ultrafast time scale. This Account aims to discuss some of the chemical physical effects ruling these processes, a fundamental step toward their control, based on our experience on nucleic acids.In the last 15 years, we have, indeed, studied the photophysics and photochemistry of DNA and its components. By combining different quantum mechanical methods, we investigated the molecular processes responsible for the damage of the genetic code or, on the contrary, those preventing it by dissipating the excess energy deposited in the system by UV absorption. Independently of its fundamental biological role, DNA, with its fluctuating closely stacked bases stabilized by weak nonbonding interactions, can be considered a prototypical MCA. Therefore, it allows one to tackle within a single system many of the conceptual and methodological challenges involved in the study of photoinduced processes in MCA.In this Account, by using the outcome of our studies on oligonucleotides as a guideline, we thus highlight the most critical modellistic issues to be faced when studying, either experimentally or computationally, the interaction between UV light and DNA and, at the same time, bring out their general relevance for the study of MCAs.We first discuss the rich photoactivated dynamics of nucleobases (the chromophores), highlighting the main effects modulating the interplay between their excited states and how the latter can affect the photoactivated dynamics of the polynucleotides, either providing effective monomer-like nonradiative decay routes or triggering reactive processes (e.g., triplet generation).We then tackle the reaction paths involving multiple bases, showing that in the DNA duplex the most important ones involve two stacked bases, forming a neutral excimer or a charge transfer (CT) state, which exhibit different spectral signatures and photochemical reactivity. In particular, we analyze the factors affecting the dynamic equilibrium between the excimer and CT, such as the fluctuations of the backbone or the rearrangement of the solvent.Next, we highlight the importance of the effects not directly connected to the chromophores, such as the flexibility of the backbone or the solvent effect. The former, affecting the stacking geometry of the bases, can determine the preference between different deactivation paths. The latter is particularly influential for CT states, making very important an accurate treatment of dynamical solvation effects, involving both the solvent bulk and specific solute-solvent interactions.In the last section, we describe the main methodological challenges related to the study of polynucleotide excited states and stress the benefits derived by the integration of complementary approaches, both computational and experimental. Only exploiting different point of views, in our opinion, it is possible to shed light on the complex phenomena triggered by light absorption in DNA, as in every MCA.

Excited-State Properties for Extended Systems: Efficient Hybrid Density Functional Methods

Time-dependent density functional theory has become state-of-the-art for describing photophysical and photochemical processes in extended materials because of its affordable cost. The inclusion of exact exchange was shown to be essential for the correct description of the long-range asymptotics of electronic interactions and thus a well-balanced description of valence, Rydberg, and charge-transfer excitations. Several approaches for an efficient treatment of exact exchange have been established for the ground state, while implementations for excited-state properties are rare. Furthermore, the high computational costs required for excited-state properties in comparison to ground-state computations often hinder large-scale applications on periodic systems with hybrid functional accuracy. We therefore propose two approximate schemes for improving computational efficiency for the treatment of exact exchange. Within the auxiliary density matrix method (ADMM), exact exchange is estimated using a relatively small auxiliary basis and the introduced basis set incompleteness error is compensated by an exchange density functional correction term. Benchmark results for a test set of 35 molecules demonstrate that the mean absolute error introduced by ADMM is smaller than 0.3 pm for excited-state bond lengths and in the range of 0.02-0.04 eV for vertical excitation, adiabatic excitation, and fluorescence energies. Computational timings for a series of covalent-organic frameworks demonstrate that a speed-up of at least 1 order of magnitude can be achieved for excited-state geometry optimizations in comparison to conventional hybrid functionals. The second method is to use a semiempirical tight binding approximation for both Coulomb and exchange contributions to the excited-state kernel. This simplified Tamm-Dancoff approximation (sTDA) achieves an accuracy comparable to approximated hybrid density functional theory when referring to highly accurate coupled-cluster reference data. We find that excited-state bond lengths deviate by 1.1 pm on average and mean absolute errors in vertical excitation, adiabatic excitation, and fluorescence energies are in the range of 0.2-0.5 eV. In comparison to ADMM-approximated hybrid functional theory, sTDA accelerates the computation of broad-band excitation spectra by 1 order of magnitude, suggesting its potential use for large-scale screening purposes.

Reference Energies for Cyclobutadiene: Automerization and Excited States

Cyclobutadiene is a well-known playground for theoretical chemists and is particularly suitable to test ground- and excited-state methods. Indeed, due to its high spatial symmetry, especially at the D_4h square geometry but also in the D_2h rectangular arrangement, the ground and excited states of cyclobutadiene exhibit multiconfigurational characters and single-reference methods, such as standard adiabatic time-dependent density-functional theory (TD-DFT) or standard equation-of-motion coupled cluster (EOM-CC), are notoriously known to struggle in such situations. In this work, using a large panel of methods and basis sets, we provide an extensive computational study of the automerization barrier (defined as the difference between the square and rectangular ground-state energies) and the vertical excitation energies at D_2h and D_4h equilibrium structures. In particular, selected configuration interaction (SCI), multireference perturbation theory (CASSCF, CASPT2, and NEVPT2), and coupled-cluster (CCSD, CC3, CCSDT, CC4, and CCSDTQ) calculations are performed. The spin-flip formalism, which is known to provide a qualitatively correct description of these diradical states, is also tested within TD-DFT (combined with numerous exchange-correlation functionals) and the algebraic diagrammatic construction [ADC(2)-s, ADC(2)-x, and ADC(3)]. A theoretical best estimate is defined for the automerization barrier and for each vertical transition energy.

Adiabatic models for the quantum dynamics of surface scattering with lattice effects

In this contribution, we review models for the lattice effects in quantum dynamics calculations on surface scattering, which is important to modeling heterogeneous catalysis for achieving an interpretation of experimental measurements. Unlike dynamics models for reactions in the gas phase, those for heterogeneous reactions have to include the effects of the surface. For manageable computational costs in calculations, the effects of static surface (SS) are firstly modeled as this is simply and easily implemented. Then, the SS model has to be improved to include the effects of the flexible surface, that is the lattice effects. To do this, various surface models have been designed where the coordinates of the surface atoms are introduced in the Hamiltonian operator, especially those of the top surface atom. Based on this model Hamiltonian operator, extensive multi-dimension quantum dynamics calculations can be performed to recover the lattice effects. Here, we first review an overview of the techniques in constructing the Hamiltonian operator, which is a sum of the kinetic energy operator (KEO) and potential energy surface (PES). Since the PES containing the coordinates of the surface atoms in a cell is still expensive, the SS model is often accepted. We consider a mathematical model, called the coupled harmonic oscillator (CHO) model, to introduce the concepts of adiabatic and diabatic representations for separating the molecule and surface. Under the adiabatic model, we further introduce the expansion model where the potential function is Taylor expanded around the optimized geometry of the surface. By an expansion model truncated at the first and second order, various coupling surface models between the molecule and surface are derived. Moreover, by further and deeply understanding the adiabatic representation, an effective Hamiltonian operator is obtained by optimizing the total wave function in factorized form. By this factorized form of wave function and effective Hamiltonian operator, the geometry phase of the surface wave function is theoretically found. This theoretical prediction may be measured by carefully designing experiments. Finally, discussions on the adiabatic representation, the PES construction, and possibility of the classical-dynamics solutions are given. Based on these discussions, a simple outlook on the dynamics of photocatalytics is finally given.

A Mountaineering Strategy to Excited States: Revising Reference Values with EOM-CC4

In the framework of the computational determination of highly accurate vertical excitation energies in small organic compounds, we explore the possibilities offered by the equation-of-motion formalism relying on the approximate fourth-order coupled-cluster (CC) method, CC4. We demonstrate, using an extended set of more than 200 reference values based on CC including up to quadruples excitations (CCSDTQ), that CC4 is an excellent approximation to CCSDTQ for excited states with a dominant contribution from single excitations with an average deviation as small as 0.003 eV. We next assess the accuracy of several additive basis set correction schemes, in which vertical excitation energies obtained with a compact basis set and a high-order CC method are corrected with lower-order CC calculations performed in a larger basis set. Such strategies are found to be overall very beneficial, though their accuracy depends significantly on the actual scheme. Finally, CC4 is employed to improve several theoretical best estimates of the QUEST database for molecules containing between four and six (nonhydrogen) atoms, for which previous estimates were computed at the CCSDT level.

The Role of Local DFT+U Minima in the First-Principles Modeling of the Metal-Insulator Transition in Vanadium Dioxide

The DFT+U method is frequently employed to improve the first-principles description of strongly correlated materials. However, it is prone to deliver metastable electronic minima. While these local minima of the DFT+U method are often considered to be computational artifacts, their physical meaning and relationship to true excited states remains unclear. In this work, the possibility of theoretically modeling transformations in the solid state that require thermal or optical excitations of electrons is explored, taking into account the metastable states of the computationally undemanding DFT+U formalism. For this purpose, we choose to examine the example of the VO₂ metal-insulator transition. Metastable states that are located on different electronic potential energy surfaces are found to correspond to experimentally observed VO₂ phases. The identified metastable electronic states can be used to model the collapse of the VO₂ band gap at elevated temperatures and upon photoexcitation as well as other monoclinic-monoclinic phase transformations. The results suggest that local DFT+U minima can indeed carry physical meaning, while they remain under-reported in theoretical literature on transition metal oxides like VO₂.