Analytical gradients of complete active space self-consistent field energies using Cholesky decomposition: Geometry optimization and spin-state energetics of a ruthenium nitrosyl complex

Analytical Evaluation of Ground State Gradients in Quantum Electrodynamics Coupled Cluster Theory

Analytical gradients of potential energy surfaces play a central role in quantum chemistry, allowing for molecular geometry optimizations and molecular dynamics simulations. In strong coupling conditions, potential energy surfaces can account for strong interactions between matter and the quantized electromagnetic field. In this paper, we derive expressions for the ground state analytical gradients in quantum electrodynamics coupled cluster theory. We also present a Cholesky-based implementation for the coupled cluster singles and doubles model. We report timings to show the performance of the implementation and present optimized geometries to highlight cavity-induced molecular orientation effects in strong coupling conditions.

Interpreting the Cu-O2 Antibonding Nature in Two Cu-O2 Complexes from Cu L-Edge X-ray Absorption Spectra

Cu-O2 structures play important roles in bioinorganic chemistry and enzyme catalysis, where the bonding between the Cu and O2 parts serves as a fundamental research concern. Here, we performed a multiconfigurational study on the copper L2,3-edge X-ray absorption spectra (XAS) of two copper enzyme model complexes to gain a better understanding of the antibonding nature from the clearly interpreted structure-spectroscopy relation. We obtained spectra in good agreement with the experiments by using the restricted active space second-order perturbation theory (RASPT2) method, which facilitated reliable chemical analysis. Spectral feature interpretations were supported by computing the spin-orbit natural transition orbitals. All major features were assigned to be mainly from Cu 2p to antibonding orbitals between Cu 3d and O2 π*, Cu 3d-πO-O* (type A), and a few also to mixed antibonding/bonding orbitals between Cu 3d and O2 π, Cu 3d ± πO-O (type M). Our calculations provided a clear illustration of the interactions between Cu 3d and O2 π*/π orbitals that are carried in the metal L-edge XAS.

Computation of NMR shieldings at the CASSCF level using gauge-including atomic orbitals and Cholesky decomposition

We present an implementation of coupled-perturbed complete active space self-consistent field (CP-CASSCF) theory for the computation of nuclear magnetic resonance chemical shifts using gauge-including atomic orbitals and Cholesky decomposed two-electron integrals. The CP-CASSCF equations are solved using a direct algorithm where the magnetic Hessian matrix-vector product is expressed in terms of one-index transformed quantities. Numerical tests on systems with up to about 1300 basis functions provide information regarding both the computational efficiency and limitations of our implementation.

Efficient implementation of molecular CCSD gradients with Cholesky-decomposed electron repulsion integrals

We present an efficient implementation of ground and excited state coupled cluster singles and doubles (CCSD) gradients based on Cholesky-decomposed electron repulsion integrals. Cholesky decomposition and density fitting are both inner projection methods, and, thus, similar implementation schemes can be applied for both methods. One well-known advantage of inner projection methods, which we exploit in our implementation, is that one can avoid storing large V3 O and V4 arrays by instead considering three-index intermediates. Furthermore, our implementation does not require the formation and storage of Cholesky vector derivatives. The new implementation is shown to perform well, with less than 10% of the time spent calculating the gradients in geometry optimizations. Furthermore, the computational time per optimization cycle is significantly lower compared to other implementations based on an inner projection method. We showcase the capabilities of the implementation by optimizing the geometry of the retinal molecule (C20H28O) at the CCSD/aug-cc-pVDZ level of theory.

Second-Order CASSCF Algorithm with the Cholesky Decomposition of the Two-Electron Integrals

In this contribution, we present the implementation of a second-order complete active space–self-consistent field (CASSCF) algorithm in conjunction with the Cholesky decomposition of the two-electron repulsion integrals. The algorithm, called norm-extended optimization, guarantees convergence of the optimization, but it involves the full Hessian and is therefore computationally expensive. Coupling the second-order procedure with the Cholesky decomposition leads to a significant reduction in the computational cost, reduced memory requirements, and an improved parallel performance. As a result, CASSCF calculations of larger molecular systems become possible as a routine task. The performance of the new implementation is illustrated by means of benchmark calculations on molecules of increasing size, with up to about 3000 basis functions and 14 active orbitals.

Modelling Photoionisation in Isocytosine: Potential Formation of Longer-Lived Excited State Cations in its Keto Form

Studying the effects of UV and VUV radiation on non-canonical DNA/RNA nucleobases allows us to compare how they release excess energy following absorption with respect to their canonical counterparts. This has attracted much research attention in recent years because of its likely influence on the origin of our genetic lexicon in prebiotic times. Here we present a CASSCF and XMS-CASPT2 theoretical study of the photoionisation of non-canonical pyrimidine nucleobase isocytosine in both its keto and enol tautomeric forms. We analyse their lowest energy cationic excited states including 2 π + , 2 n O + and 2 n N + and compare these to the corresponding electronic states in cytosine. Investigating lower-energy decay pathways we find - unexpectedly - that keto-isocytosine+ presents a sizeable energy barrier potentially inhibiting decay to its cationic ground state, whereas enol-isocytosine+ features a barrierless and consequently ultrafast pathway analogous to the one previously found for the canonical (keto) form of cytosine+ . Dynamic electron correlation reduces the energy barrier in the keto form substantially (by ∼1 eV) but it is nevertheless still present. We additionally compute the UV/Vis absorption signals of the structures encountered along these decay channels to provide spectroscopic fingerprints to assist future experiments in monitoring these intricate photo-processes.

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

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

Analytic gradients for multiconfiguration pair-density functional theory with density fitting: Development and application to geometry optimization in the ground and excited states

Density fitting reduces the computational cost of both energy and gradient calculations by avoiding the computation and manipulation of four-index electron repulsion integrals. With this algorithm, one can efficiently optimize the geometries of large systems with an accurate multireference treatment. Here, we present the derivation of multiconfiguration pair-density functional theory for energies and analytic gradients with density fitting. Six systems are studied, and the results are compared to those obtained with no approximation to the electron repulsion integrals and to the results obtained by complete active space second-order perturbation theory. With the new approach, there is an increase in the speed of computation with a negligible loss in accuracy. Smaller grid sizes have also been used to reduce the computational cost of multiconfiguration pair-density functional theory with little effect on the optimized geometries and gradient values.

The bimetallic and the anchoring group effects on both optical and charge transport properties of hexaphyrin amethyrin

Bimetallic Cu(ii)-hexaphyrin amethyrin proposed as a molecular switch operated by the application of an external magnetic field.

Modern quantum chemistry with [Open]Molcas

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

Structural and Optical Properties of Metal-Nitrosyl Complexes

The electronic, structural and optical properties (including Spin–Orbit Coupling) of metal nitrosyl complexes [M(CN)₅(NO)]²⁻ (M = Fe, Ru or Os) are investigated by means of Density Functional Theory, TD-DFT and MS-CASPT2 based on an RASSCF wavefunction. The energy profiles connecting the N-bound (η¹-N), O-bound (η¹-O) and side-on (η²-NO) conformations have been computed at DFT level for the closed shell singlet electronic state. For each structure, the lowest singlet and triplet states have been optimized in order to gain insight into the energy profiles describing the conformational isomerism in excited states. The energetics of the three complexes are similar—with the N-bound structure being the most stable—with one exception, namely the triplet ground state of the O-bound isomer for the iron complex. The conformation isomerism is highly unfavorable in the S₀ electronic state with the occurrence of two energy barriers higher than 2 eV. The lowest bands of the spectra are assigned to MLCT_NO/LLCT_NO transitions, with an increasing MLCT character going from iron to osmium. Two low-lying triplet states, T1 (MLCT_NO/LLCT_NO) and T2 (MLCT_NO/IL_NO), seem to control the lowest energy profile of the excited-state conformational isomerism.

Evaluation of full valence correlation energies and gradients

Complete-active-space self-consistent field (CASSCF) wave functions are central to understanding strongly correlated molecules as they capture the entirety of electronic interactions within a subset of the orbital space. The most interesting case for CASSCF is the full valence limit, where all bonding and an equal number of virtual orbitals are included in the active space, and no approximation is made in selecting the important valence orbitals or electrons. While conventional algorithms require exponential computational time to evaluate full valence CASSCF, this article shows that the method of increments can do the same with polynomial effort, in a new method denoted iCASSCF. The method of increments can also provide density matrices and other necessary ingredients for the construction of the nuclear gradient. These goals are met through a many-body expansion that breaks the problem into smaller pieces that are subsequently reassembled to form close approximations of conventional CAS results. Practical demonstrations on a number of medium-sized molecules, with up to 116 valence electrons correlated in 116 orbitals, show the power of this methodology.

Excited-State Dynamics of Thienoguanosine, an Isomorphic Highly Fluorescent Analogue of Guanosine

Thienoguanosine (^th G) is an isomorphic analogue of guanosine with promising potentialities as fluorescent DNA label. As a free probe in protic solvents, ^th G exists in two tautomeric forms, identified as the H1, being the only one observed in nonprotic solvents, and H3 keto-amino tautomers. We herein investigate the photophysics of ^th G in solvents of different polarity, from water to dioxane, by combining time-resolved fluorescence with PCM/TD-DFT and CASSCF calculations. Fluorescence lifetimes of 14.5-20.5 and 7-13 ns were observed for the H1 and H3 tautomers, respectively, in the tested solvents. In methanol and ethanol, an additional fluorescent decay lifetime (≈3 ns) at the blue emission side (λ≈430 nm) as well as a 0.5 ns component with negative amplitude at the red edge of the spectrum, typical of an excited-state reaction, were observed. Our computational analysis explains the solvent effects observed on the tautomeric equilibrium. The main radiative and nonradiative deactivation routes have been mapped by PCM/TD-DFT calculations in solution and CASSCF in the gas phase. The most easily accessible conical intersection, involving an out-of plane motion of the sulfur atom in the five-membered ring of ^th G, is separated by a sizeable energy barrier (≥0.4 eV) from the minimum of the spectroscopic state, which explains the large experimental fluorescence quantum yield.

Ultrafast and radiationless electronic excited state decay of uracil and thymine cations: computing the effects of dynamic electron correlation

In this article we characterise the radiationless decay of the first few electronic excited states of the cations of DNA/RNA nucleobases uracil and thymine, including the effects of dynamic electron correlation on energies and geometries (optimised with XMS-CASPT2). In both systems, we find that one state of ²n and another two of ²π⁺ character can be populated following photoionisation, and their different minima and interstate crossings are located. We find strong similarities between uracil and thymine cations: with accessible conical intersections suggesting that depopulation of their electronic excited states takes place on ultrafast timescales in both systems, suggesting that they are photostable in agreement with previous theoretical (uracil⁺) evidence. We find that dynamic electron correlation separates the energy levels of the "3-state" conical intersection (D₂/D₁/D₀)_CI previously located with CASSCF for uracil⁺, which will therefore have a different geometry and higher energy. Simulating the electronic and vibrational absorptions allows us to characterise spectral fingerprints that could be used to monitor these cation photo-processes experimentally.

Valence ππ* Excitations in Benzene Studied by Multiconfiguration Pair-Density Functional Theory

We explore the valence singlet and triplet ππ* excitations of benzene with complete active pace self-consistent field (CASSCF) theory, complete active space perturbation theory (CASPT2), and multiconfiguration pair-density functional theory (MC-PDFT) for four different choices of active space. We propose a new way to quantify the covalent and ionic character of the electronic states in terms of the components of the total electronic energy. We also explore the effect of scaling the exchange and correlation components of the on-top density functional used in MC-PDFT; we observe that increasing the exchange contribution improves the MC-PDFT excitation energies for benzene.

Photophysical properties and the NO photorelease mechanism of a ruthenium nitrosyl model complex investigated using the CASSCF-in-DFT embedding approach

A complete state-averaged active space self-consistent field (SA-CASSCF) calculation by means of the SA-CASSCF(18,14)-in-BP86 Miller-Manby embedding approach was performed to explore the ground and excited electronic states of the trans-[RuCl(NO)(NH₃)₄]²⁺ complex. Insights into the NO photodissociation mechanism and Ru-NO bonding properties are provided. In addition, spin-orbit (SO) interactions were taken into account to describe and characterize the spin-forbidden transitions observed at the low-energy regions of the trans-[RuCl(NO)(NH₃)₄]²⁺ UV-Vis spectrum. The SA-CASSCF(18,14)-in-BP86 electronic spectrum is in great agreement with the experimental data of Schreiner [Schreiner et al., Inorg. Chem., 1972, 11, 880].

Multireference Perturbation Theory with Cholesky Decomposition for the Density Matrix Renormalization Group

We present a second-order N-electron valence state perturbation theory (NEVPT2) based on a density matrix renormalization group (DMRG) reference wave function that exploits a Cholesky decomposition of the two-electron repulsion integrals (CD-DMRG-NEVPT2). With a parameter-free multireference perturbation theory approach at hand, the latter allows us to efficiently describe static and dynamic correlation in large molecular systems. We demonstrate the applicability of CD-DMRG-NEVPT2 for spin-state energetics of spin-crossover complexes involving calculations with more than 1000 atomic basis functions. We first assess, in a study of a heme model, the accuracy of the strongly and partially contracted variant of CD-DMRG-NEVPT2 before embarking on resolving a controversy about the spin ground state of a cobalt tropocoronand complex.

Multiconfigurational Second-Order Perturbation Theory with Frozen Natural Orbitals Extended to the Treatment of Photochemical Problems

A new flavor of the frozen natural orbital complete active space second-order perturbation theory method (FNO-CASPT2, Aquilante et al., J. Chem. Phys. 131, 034113) is proposed herein. In this new implementation, the virtual space in Cholesky decomposition-based CASPT2 computations (CD-CASPT2) is truncated by excluding those orbitals that contribute the least toward preserving a predefined value of the trace of an approximate density matrix, as that represents a measure of the amount of dynamic correlation retained in the model. In this way, the amount of correlation included is practically constant at all nuclear arrangements, thus allowing for the computation of smooth electronic states surfaces and energy gradients-essential requirements for theoretical studies in photochemistry. The method has been benchmarked for a series of relevant biochromophores for which large speed-ups have been recorded while retaining the accuracy achieved in the corresponding CD-CASPT2 calculations. Both vertical excitation energies and gradient calculations have been carried out to establish general guidelines as to how much correlation needs to be retained in the calculation for the results to be consistent with the CD-CASPT2 findings. Our results feature errors within a tenth of an eV for the most difficult cases and have been validated to be used for gradient computations where an up to 3-fold speed-up is observed depending on the size of the system and the basis set employed.

Biological properties of novel ruthenium- and osmium-nitrosyl complexes with azole heterocycles

Since the discovery that nitric oxide (NO) is a physiologically relevant molecule, there has been great interest in the use of metal nitrosyl compounds as antitumor pharmaceuticals. Particularly interesting are those complexes which can deliver NO to biological targets. Ruthenium- and osmium-based compounds offer lower toxicity compared to other metals and show different mechanisms of action as well as different spectra of activity compared to platinum-based drugs. Novel ruthenium- and osmium-nitrosyl complexes with azole heterocycles were studied to elucidate their cytotoxicity and possible interactions with DNA. Apoptosis induction, changes of mitochondrial transmembrane potential and possible formation of reactive oxygen species were investigated as indicators of NO-mediated damage by flow cytometry. Results suggest that ruthenium- and osmium-nitrosyl complexes with the general formula (indazolium)[cis/trans-MCl₄(NO)(1H-indazole)] have pronounced cytotoxic potency in cancer cell lines. Especially the more potent ruthenium complexes strongly induce apoptosis associated with depolarization of mitochondrial membranes, and elevated reactive oxygen species levels. Furthermore, a slight yet not unequivocal trend to accumulation of intracellular cyclic guanosine monophosphate attributable to NO-mediated effects was observed.

Benzophenone Ultrafast Triplet Population: Revisiting the Kinetic Model by Surface-Hopping Dynamics

The photochemistry of benzophenone, a paradigmatic organic molecule for photosensitization, was investigated by means of surface-hopping ab initio molecular dynamics. Different mechanisms were found to be relevant within the first 600 fs after excitation; the long-debated direct (S1 → T1) and indirect (S1 → T2 → T1) mechanisms for population of the low-lying triplet state are both possible, with the latter being prevalent. Moreover, we established the existence of a kinetic equilibrium between the two triplet states, never observed before. This fact implies that a significant fraction of the overall population resides in T2, eventually allowing one to revisit the usual spectroscopic assignment proposed by transient absorption spectroscopy. This finding is of particular interest for photocatalysis as well as for DNA damages studies because both T1 and T2 channels are, in principle, available for benzophenone-mediated photoinduced energy transfer toward DNA.

Molcas 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas-Kroll-Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC-PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large-scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

Orbital entanglement and CASSCF analysis of the Ru-NO bond in a Ruthenium nitrosyl complex

Multiconfigurational wavefunction analysis and entanglement measures based on von Neumann entropy shed light on the electronic structure of a Ru nitrosyl complex, in particular on the Ru–NO bond.

Complete active space self-consistent field (CASSCF) wavefunctions and an orbital entanglement analysis obtained from a density-matrix renormalisation group (DMRG) calculation are used to understand the electronic structure, and, in particular, the Ru–NO bond of a Ru nitrosyl complex. Based on the configurations and orbital occupation numbers obtained for the CASSCF wavefunction and on the orbital entropy measurements evaluated for the DMRG wavefunction, we unravel electron correlation effects in the Ru coordination sphere of the complex. It is shown that Ru–NO π bonds show static and dynamic correlation, while other Ru–ligand bonds feature predominantly dynamic correlation. The presence of static correlation requires the use of multiconfigurational methods to describe the Ru–NO bond. Subsequently, the CASSCF wavefunction is analysed in terms of configuration state functions based on localised orbitals. The analysis of the wavefunctions in the electronic singlet ground state and the first triplet state provides a picture of the Ru–NO moiety beyond the standard representation based on formal oxidation states. A distinct description of the Ru and NO fragments is advocated. The electron configuration of Ru is an equally weighted superposition of Ru^II and Ru^III configurations, with the Ru^III configuration originating from charge donation mostly from Cl ligands. However, and contrary to what is typically assumed, the electronic configuration of the NO ligand is best described as electroneutral.