Second-order MCSCF optimization revisited. I. Improved algorithms for fast and robust second-order CASSCF convergence

Spectroscopic Characterization of S3O Isomers: Potential Contributor to the Unknown UV Absorber in Venus's Atmosphere

A long-standing question about Venus's atmosphere concerns the origin and distribution of the mysterious ultraviolet (UV) absorber, which is responsible for the absorption band observed in the near-UV and blue regions of Venus's spectrum. In this work, we investigate the electronic structure of S3O isomers, their excited states, and their photoabsorption cross sections using quantum chemistry methods. Our results demonstrate that the cis- and trans-S3O isomers exhibit strong absorption in the 300-500 nm range, making them plausible candidates for the unknown near-UV absorber. The electronic spectrum of ring-S3O is featureless, whereas planar-S3O exhibits absorption around 300 nm. Additionally, spectroscopic data were generated using vibrational configuration interaction (VCI) for the most stable isomer and its isotopes to aid in their detection in the laboratory or the interstellar medium. This study provides valuable insights into the photochemistry of S3O isomers and the complex sulfur chemistry within Venus's atmosphere.

Rational Design and Reaction Mechanism Study of the Photochemical Rearrangement of Fulvene Derivatives

Photochemistry is considered one of the most efficient and reproducible techniques in organic synthesis. Recently, List and co-workers reported an efficient UV light triggered photochemical synthesis of spiro[2,4]heptadiene from fulvenes with different substituents ( Angew. Chem., Int. Ed. 2023, 62, e202303119); however, the mechanistic details remain unclear, and the intermediates have not been characterized. To facilitate the applications of this novel photochemical reaction, we theoretically designed a series of fulvene derivatives with different parent molecular skeletons for analyzing the substitution effects, and two of the representative fulvenes were synthesized for investigating the reaction mechanisms by employing time-resolved transient absorption spectroscopy (TA) experiments. It has been found that instead of density functional theory, the second-order n-electron valence state perturbation theory is necessary to acquire reliable theoretical characterization of the fulvenes examined. Our designed fulvenes were found to undergo the photorearrangement cyclopropanation reaction on the basis of photoproduct analysis. The intermediate species involved in the intramolecular hydrogen atom transfer and cyclization processes within the photorearrangement reaction were characterized by TA spectroscopy, and the full reaction pathways were proposed. Our work not only reveals the detailed mechanism of this photorearrangement reaction but also demonstrates the significance of appropriate theoretical methods for rational molecular design.

UV-photoelectron spectroscopy and MS-CASPT2/CASSCF study of the thermolysis of azidoethyl-methyl sulfide: Characterization and mechanism of the formation of S-methyl-N-sulfenylethanimine

The thermal decomposition of azidoethyl methyl sulfide was studied by real-time UV-photoelectron spectroscopy (UV-PES) at temperatures ranging from 773 to 1023 K. Different ionization energies were obtained using density functional theory calculations to assign UV-PES spectra. The complete active space self-consistent field and multistate second-order perturbation methods were used to predict the formation of different species present in the thermal decomposition process. N2 and S-methyl-N-sulfenylethanimine are generated at 773 K. The first step of the reaction is the dissociation of the molecule into nitrene and nitrogen. The spin state (singlet or triplet) of nitrene formed in the first step of the reaction is temperature-dependent. At low temperatures (T ≤ 650 K), both states are formed with almost the same probability; in contrast, at high temperatures (T ≥ 1000 K), singlet nitrene is the majority intermediate. From this singlet nitrene, three stable reaction products were detected in the experiments: an imine derivative, a four-member cyclic derivative, and a sulfenyl derivative.

Analytical SA-HCISCF Nuclear Gradients from Spin-Adapted Heat-Bath Configuration Interaction

This work reports an implementation of the analytical nuclear gradients and nonadiabatic couplings with state-averaged SCF wave functions from a spin-pure selected configuration interaction (SCI) method. At the core of the implementation lies the evaluation of the Lagrange multipliers required for the variational calculation of the nuclear gradient. Using the same code infrastructure, we developed a fully CI-coupled second-order orbital optimization method. Both the calculation of the nuclear gradient and the second-order orbital optimization make use of density fitting in order to accelerate the calculation of the two-electron integrals. We demonstrate the use of analytical nuclear gradients in excited-state geometry optimizations for conjugated molecules. In addition, the first triplet excited-state geometry of a transition-metal catalyst, Fe(PDI), was optimized with up to 30 orbitals in the active space. Our results outline the capabilities of the implemented methods as well as directions for future work.

Nuclear-Electronic Orbital Multireference Configuration Interaction for Ground and Excited Vibronic States and Fundamental Insights into Multicomponent Basis Sets

The nuclear-electronic orbital (NEO) approach incorporates nuclear quantum effects into quantum chemistry calculations by treating specified nuclei quantum mechanically, equivalently to the electrons. Within the NEO framework, excited states are vibronic states representing electronic and nuclear vibrational excitations. The NEO multireference configuration interaction (MRCI) method presented herein provides accurate ground and excited vibronic states. The electronic and nuclear orbitals are optimized with a NEO multiconfigurational self-consistent field (NEO-MCSCF) procedure, thereby including both static and dynamic correlation and allowing the description of double and higher excitations. The accuracy of the NEO-MRCI method is illustrated by computing the ground state protonic densities and excitation energies of the vibronic states for four molecular systems with the hydrogen nucleus treated quantum mechanically. In addition, revised conventional electronic basis sets adapted for quantized nuclei are developed and shown to be essential for achieving this level of accuracy. The NEO-MRCI approach, as well as the strategy for revising electronic basis sets, will play a critical role in multicomponent quantum chemistry.

Semiclassical Approach to Computing Vibrationally Resolved Ionization Cross Sections for Molecular Nitrogen

A semiclassical model based on the Gryzinski theory is used to compute vibrationally resolved electron-impact ionization cross sections for molecular nitrogen. This model extends the approach used by Wünderlich for molecular hydrogen and its isotopomeres. The multireference configuration interaction (MRCI) method in Molpro is used with complete active space self-consistent field reference wave functions to compute potential energy curves (PECs) and electronic wave functions for several states of interest. Nuclear wave functions and vibrational energy levels are computed from the MRCI PECs using the Fourier grid Hamiltonian method. The target orbital electron kinetic energies, Franck-Condon factors, and transition energies are parameters for the semiclassical model and are calculated directly from the computed electronic and nuclear wave functions and vibrational energy levels. The target electron kinetic energies are computed as the expectation value of the one-electron kinetic energy operator for the product of the appropriate nuclear vibrational wave function and the MRCI natural orbital for a particular state-to-state ionization process. From the fully vibrationally resolved ionization cross sections, lumped and total cross sections are calculated by summing the partial cross sections over the closure relationship of the Franck-Condon theory.

The Electronic Structure and Bonding in Some Small Molecules

The electronic structures of the first- and second-row homonuclear diatomics, XeF2, and the weakly bound dimers of nitric oxide and nitrogen dioxide molecules in their ground states are discussed in terms of molecular orbital (MO) theory and, where possible, valence bond theories. The current work is extended and supported by restricted and unrestricted Hartree-Fock (RHF and UHF) self-consistent field (SCF), complete active space SCF (CASSCF), multi-reference configuration interaction (MRCI), coupled cluster CCSD(T), and unrestricted Kohn-Sham (UKS) density functional calculations using a polarized triple-zeta basis. The dicarbon (C2) molecule is especially poorly described by RHF theory, and it is argued that the current MO theories taught in most undergraduate courses should be extended in recognition of the fact that the molecule requires at least a two-configuration treatment.

Ab initio electronic structures and total internal partition sums of FeH+/2

In the present work, we studied 27 FeH+ and 6 FeH2+ electronic states using multireference configuration interaction (MRCI), Davidson-corrected MRCI (MRCI+Q), and coupled cluster singles doubles and perturbative triples [CCSD(T)] wavefunction theory (WFT) calculations conjoined with large quadruple-ζ and quintuple-ζ quality correlation consistent basis sets. We report their potential energy curves (PEC), energy related properties, spectroscopic parameters, and spin-orbit couplings. Dipole moment curves (DMC) and transition dipole moment curves (TDMC) of several low-lying electronic states of FeH+ and FeH2+ are also introduced. The ground state of FeH+ is a single-reference X5Δ (6σ27σ13π21δ3) with an adiabatic D0 of ∼52 kcal mol-1, which is in agreement with the experimental value. The states with the largest adiabatic binding energies of FeH2+ (4Π and 4Δ) are multireference in nature with an approximate D0 of 22 kcal mol-1. We used CCSD(T) μ of the FeH+(X5Δ) to assess the density functional theory (DFT) errors associated with a series of functionals that span multiple rungs of Jacob's ladder of density functional approximation (DFA) and observed a general trend of improving μ when moving to more expensive functionals at the higher rungs. We expect weak spectral bands to be produced from the low-lying electronic states of FeH2+ and FeH+ due to their lower transition μ values. Lastly, we present results for the total internal partition function sums (TIPS) of FeH+ and FeH2+, which have not been presented in the literature before.

Electronic Spectroscopy of the Halogenated Criegee Intermediate, ClCHOO: Experiment and Theory

A chlorine-substituted Criegee intermediate, ClCHOO, is photolytically generated using a diiodo precursor, detected by VUV photoionization at 118 nm, and spectroscopically characterized via ultraviolet-visible (UV-vis)-induced depletion of m/z = 80 under jet cooled conditions. UV-vis excitation resonant with a π* ← π transition yields a significant ground state depletion, indicating a strong electronic transition and rapid photodissociation. The broad absorption spectrum peaks at 350 nm and is attributed to contributions from both syn (∼70%) and anti (∼30%) conformers of ClCHOO based on spectral simulations using a nuclear ensemble method. Electronic structure theory shows significant differences in the vertical excitation energies of the two conformers (330 and 371 nm, respectively) as well as their relative stabilities in the ground and excited electronic states associated with the π* ← π transition. Natural bond orbital analysis reveals significant and nonintuitive nonbonding contributions to the relative stabilities of the syn and anti conformers.

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.

On the ground and excited electronic states of LaCO and AcCO

High-level ab initio electronic structure analysis of correlated lanthanide- and actinide-based species is laborious to perform and consequently limited in the literature. In the present work, the ground and electronically excited states of LaCO and AcCO molecules were explored utilizing the multireference configuration interaction (MRCI), Davidson corrected MRCI (MRCI+Q), and coupled cluster singles doubles and perturbative triples [CCSD(T)] quantum chemical tools conjoined with correlation consistent triple-ζ and quadruple-ζ quality all-electron Douglas-Kroll (DK) basis sets. The full potential energy curves (PECs), dissociation energies (Des), excitation energies (Tes), bond lengths (res), harmonic vibrational frequencies (ωes), and chemical bonding patterns of low-lying electronic states of LaCO and AcCO are introduced. The ground electronic state of LaCO is a 4Σ- (1σ11π2) which is a product of the reaction between excited La(4F) versus CO(X1Σ+), whereas the ground state of AcCO is a 12Π (1σ21π1) deriving from ground state fragments Ac(2D) + CO(X1Σ+). The spin-orbit ground states of LaCO (14Σ-3/2) and AcCO (12Π1/2) bear ∼13 and 5 kcal mol-1D0 values, respectively. At the MRCI level, the spin-orbit curves, the spin-orbit mixing, and the Tes of spin-orbit states of LaCO and AcCO were also analyzed. Lastly, the density functional theory (DFT) calculations were performed applying 16 exchange-correlation functionals that span three rungs of "Jacob's ladder" of density functional approximations (DFAs) to assess DFT errors associated on the De and ionization energy (IE) of LaCO.

Ab initio exploration of low-lying electronic states of linear and bent MNX+ (M = Ca, Sr, Ba, Ra; X = O, S, Se, Te, Po) and their origins

High-level multireference and coupled cluster quantum calculations were employed to analyze low-lying electronic states of linear-MNX+ and side-bonded-M[NX]+ (M = Ca, Sr, Ba, Ra; X = O, S, Se, Te, Po) species. Their full potential energy curves (PECs), dissociation energies (Des), geometric parameters, excitation energies (Tes), and harmonic vibrational frequencies (ωes) are reported. The first three chemically bound electronic states of MNX+ and M[NX]+ are 3∑-, 1Δ, 1∑+ and 3A″, 1A', 1A″, respectively. The 3∑-, 1Δ, 1∑+ of MNX+ originate from the M+(2D) + NX(2Π) fragments, whereas the 3A″, 1A', 1A″ states of M[NX]+ dissociate to M+(2S) + NX(2Π) as a result of avoided crossings. The MNX+ and M[NX]+ are real minima on the potential energy surface and their interconversions are possible. The M2+NX-/M2+[NX]- ionic structure is an accurate representation for their low-lying electronic states. The Des of MNX+ species were found to depend on the dipole moment (μ) of the corresponding NX ligands and a linear relationship between these two parameters was observed.

Permutation symmetry in spin-adapted many-body wave functions

In the domain of exchange-coupled polynuclear transition-metal (PNTM) clusters, local emergent symmetries exist which can be exploited to greatly increase the sparsity of the configuration interaction (CI) eigensolutions of such systems. Sparsity of the CI secular problem is revealed by exploring the site permutation space within spin-adapted many-body bases, and highly compressed wave functions may arise by finding optimal site orderings. However, the factorial cost of searching through the permutation space remains a bottleneck for clusters with a large number of metal centers. In this work, we explore ways to reduce the factorial scaling, by combining permutation and point group symmetry arguments, and using commutation relations between cumulative partial spin and the Hamiltonian operators, . Certain site orderings lead to commuting operators, from which more sparse wave functions arise. Two graphical strategies will be discussed, one to rapidly evaluate the commutators of interest, and one in the form of a tree search algorithm to predict how many and which distinct site permutations are to be analyzed, eliminating redundancies in the permutation space. Particularly interesting is the case of the singlet spin states for which an additional reversal symmetry can be utilized to further reduce the number of distinct site permutations.

Ground and Excited Electronic Structure Analysis of FeH with Correlated Wave Function Theory and Density Functional Approximations

FeH is one of the most challenging diatomic molecules to study under electronic structure theory. Here, we have successfully studied 22 electronic states of FeH using ab initio multireference configuration interaction (MRCI), Davidson-corrected MRCI (MRCI+Q), and coupled cluster singles, doubles, and perturbative triples [CCSD(T)] levels of theory. We report their potential energy curves (PECs), excitation energies, dissociation energies, equilibrium electronic configurations, and a series of spectroscopic constants with the use of augmented triple-ζ, quadruple-ζ, and quintuple-ζ quality correlation consistent basis sets. The scalar relativistic effects and active space and core electron correlation contribution on the properties of FeH are also explored. The use of a large CASSCF active space that includes 4s, 4p, 3d, and 4d orbitals of Fe and the 1s of H is critical for producing accurate full PECs with proper dissociations and predicting the exact order of the electronic states. Our findings are in harmony with the experimental results available in the literature and will serve as reference values for future studies of FeH. Furthermore, with the use of PECs, the total internal partition function sum (TIPS) of FeH was calculated across a range of temperatures. Finally, we exploited the single-reference nature of the a6Δ of FeH and its ionized product FeH+ (X5Δ) to evaluate the associated density functional theory (DFT) errors on their dissociation energies and spectroscopic parameters.

Fulminic acid: a quasibent spectacle

Fulminic acid (HCNO) played a critical role in the early development of organic chemistry, and chemists have sought to discern the structure and characteristics of this molecule and its isomers for over 200 years. The mercurial nature of the extremely flat H-C-N bending potential of fulminic acid, with a nearly vanishing harmonic vibrational frequency at linearity, remains enigmatic and refractory to electronic structure theory, as dramatic variation with both orbital basis set and electron correlation method is witnessed. To solve this problem using rigorous electronic wavefunction methods, we have employed focal point analyses (FPA) to ascertain the ab initio limit of optimized linear and bent geometries, corresponding vibrational frequencies, and the HCN + O(3P) → HCNO reaction energy. Electron correlation treatments as extensive as CCSDT(Q), CCSDTQ(P), and even CCSDTQP(H) were employed, and complete basis set (CBS) extrapolations were performed using the cc-pCVXZ (X = 4-6) basis sets. Core electron correlation, scalar relativistic effects (MVD1), and diagonal Born-Oppenheimer corrections (DBOC) were all included and found to contribute significantly in determining whether vibrationless HCNO is linear or bent. At the all-electron (AE) CCSD(T)/CBS level, HCNO is a linear molecule with ω5(H-C-N bend) = 120 cm-1. However, composite AE-CCSDT(Q)/CBS computations give an imaginary frequency (51i cm-1) at the linear optimized geometry, leading to an equilibrium structure with an H-C-N angle of 173.9°. Finally, at the AE-CCSDTQ(P)/CBS level, HCNO is once again linear with ω5 = 45 cm-1, and inclusion of both MVD1 and DBOC effects yields ω5 = 32 cm-1. A host of other topics has also been investigated for fulminic acid, including the dependence of re and ωi predictions on a variety of CBS extrapolation formulas, the question of multireference character, the N-O bond energy and enthalpy of formation, and issues that give rise to the quasibent appellation.

Wavefunction theory and density functional theory analysis of ground and excited electronic states of TaB and WB

Several low-lying electronic states of TaB and WB molecules were studied using ab initio multireference configuration interaction (MRCI), Davidson corrected MRCI (MRCI+Q), and coupled cluster singles doubles and perturbative triples [CCSD(T)] methods. Their full potential energy curves (PECs), equilibrium electron configurations, equilibrium bond distances (res), dissociation energies (Des), excitation energies (Tes), harmonic vibrational frequencies (ωes), and anharmonicities (ωexes) are reported. The MRCI dipole moment curves (DMCs) of the first 5 electronic states of both TaB and WB are also reported and the equilibrium dipole moment (μ) values are compared with the CCSD(T) μ values. The most stable 13Π (1σ22σ23σ11π3) and 15Δ (1σ22σ23σ11π21δ1) electronic states of TaB lie close in energy with ∼62 kcal mol-1De with respect to the Ta(4F) + B(2P) asymptote. However, spin-orbit coupling effects make the 15Δ0+ state the true ground state of TaB. The ground electronic state of WB (16Π) has the 1σ22σ13σ11π31δ2 electron configuration and is followed by the excited 16Σ+ and 14Δ states. Finally, the MRCI De, re, ωe, and ωexe values of the 13Π state of TaB and 16Π and 14Δ states of WB are used to assess the density functional theory (DFT) errors on a series of exchange-correlation functionals that span multiple-rungs of the Jacob's ladder of density functional approximations (DFA).

Surface hopping molecular dynamics simulation of ultrafast methyl iodide photodissociation mapped by Coulomb explosion imaging

We present a highly efficient approach to directly and reliably simulate photodissociation followed by Coulomb explosion of methyl iodide. In order to achieve statistical reliability, more than 40 000 trajectories are calculated on accurate potential energy surfaces of both the neutral molecule and the doubly charged cation. Non-adiabatic effects during photodissociation are treated using a Landau-Zener surface hopping algorithm. The simulation is performed analogous to a recent pump-probe experiment using coincident ion momentum imaging [Ziaee et al., Phys. Chem. Chem. Phys., 2023, 25, 9999-10010]. At large pump-probe delays, the simulated delay-dependent kinetic energy release signals show overall good agreement with the experiment, with two major dissociation channels leading to I(2P3/2) and I*(2P1/2) products. At short pump-probe delays, the simulated kinetic energy release differs significantly from the values obtained by a purely Coulombic approximation or a one-dimensional description of the dicationic potential energy surfaces, and shows a clear bifurcation near 12 fs, owing to non-adiabatic transitions through a conical intersection. The proposed approach is particularly suitable and efficient in simulating processes that highly rely on statistics or for identifying rare reaction channels.

Ab initio electronic structure analysis of ground and excited states of HfN0,

High-level ab initio electronic structure analysis of third-row transition metal (TM)-based diatomic species is challenging and has been perpetually lagging. In this work, fourteen and eighteen electronic states of HfN and HfN+ respectively are studied, employing multireference configuration interaction (MRCI) and coupled cluster singles doubles and perturbative triples [CCSD(T)] theories under larger correlation-consistent basis sets. Their potential energy curves (PECs), energetics, and spectroscopic parameters are reported. Core electron correlation effects on their properties are also investigated. Chemical bonding patterns of several low-lying electronic states are introduced based on the equilibrium electron configurations. The ground state of HfN (X2Σ+) has the 1σ22σ23σ11π4 electronic configuration, and the ionization of the 3σ1 electron produces the ground state of HfN+ (X1Σ+). Ground states of both HfN and HfN+ are triple bonded in nature and bear 124.86 and 109.10 kcal mol-1 binding energies with respect to their ground state fragments. The findings of this work agree well with the limited experimental literature available and provide useful reference values for future experimental analysis of HfN and HfN+.

Geometric Optimization of Restricted-Open and Complete Active Space Self-Consistent Field Wave Functions

We explore Riemannian optimization methods for Restricted-Open-shell Hartree-Fock (ROHF) and Complete Active Space Self-Consistent Field (CASSCF) methods. After showing that ROHF and CASSCF can be reformulated as optimization problems on so-called "flag manifolds", we review Riemannian optimization basics and their application to these specific problems. We compare these methods to traditional ones and find robust convergence properties without fine-tuning of numerical parameters. Our study suggests that Riemannian optimization is a valuable addition to orbital optimization for ROHF and CASSCF, warranting further investigation.

Towards reliable and efficient modeling of [Cu2O2]2+-based compound electronic structures with the partially fixed reference space protocols

This work reports a computationally efficient approach for reliable modeling of complex electronic structures based on [Cu2O2]2+ moieties. Specifically, we explore the recently developed partially fixed reference space (PFRS) protocol to minimize the active space size, taking into account the double d-shell effects. We show that the ground-state electronic structure of the core [Cu2O2]2+ model system is dominated by the d9/d10 occupations. The PFRS-crafted active spaces are further used to generate the reference wave functions for the multi-reference coupled cluster, configuration interaction, and multi-reference perturbation theory calculations. Specifically, we demonstrate that the bare [Cu2O2]2+ core can be modeled qualitatively using active spaces as small as CAS(2,2)PFRS. To obtain quantitative agreement with the reference DMRG(32,62)CI calculations, the CAS(4,4) has to be used in conjunction with the MRCCSD correction on top of it. This reliable and computationally efficient protocol is further used to model the electronic structure and properties of ammonia coordinated [Cu2O2]2+ complexes. Finally, based on the large amount of available experimental data regarding the oxo-peroxo equilibrium of [Cu2O2]2+-based systems, it is possible to formulate educated guesses regarding the effect of each experimental variable over each d-occupancy-specific state. With a large sample size of d-occupancy-specific state dependence with ligands and solvents, it should be possible to propose new ligands with specific d-occupancy and, therefore, oxidative properties based on the d-occupancy energy gaps of relatively low-cost calculations.

Anharmonic Vibrational Spectroscopy of Germanium-Containing Clusters, GexC4-x and GexSi4-x (x = 0-4), for Interstellar Detection

An extensive, high-level theoretical study on tetra-atomic germanium carbide/silicide clusters is presented. Accurate harmonic and anharmonic vibrational frequencies and rotational constants are calculated at the CCSD(T)-F12a(b)/cc-pVT(Q)Z-F12 levels of theory. With growing capabilities to discern more of the chemical composition of the interstellar medium (ISM), an accurate database of reference material is required. The presence of carbon is ubiquitous in the ISM, and silicon is known to be present in interstellar dust grains; however, germanium-containing molecules remain elusive. To begin understanding the presence and role of germanium in the ISM, we present this study of the vibrational and rotational spectroscopic properties of various germanium-containing molecules to aid in their potential identification in the ISM with modern observational tools such as the James Webb Space Telescope. Structures studied herein include rhomboidal (r-), diamond (d-), and trapezoidal (t-) tetra-atomic molecules of the form GexC4-x and GexSi4-x, where x = 0-4. The most promising structure for detection is r-Ge2C2 via the ν4 mode with a frequency of 802.7 cm-1 (12.5 μm) and an intensity of 307.2 km mol-1. Other molecules that are potentially detectable, i.e., through vibrational modes or rotational transitions, include r-Ge3C, r-GeSi3, d-GeC3, r-GeC3, and t-Ge2C2.

MBE-CASSCF Approach for the Accurate Treatment of Large Active Spaces

We present a novel implementation of the complete active space self-consistent field (CASSCF) method that makes use of the many-body expanded full configuration interaction (MBE-FCI) method to incrementally approximate electronic structures within large active spaces. On the basis of a hybrid first-order algorithm employing both Super-CI and quasi-Newton strategies for the optimization of molecular orbitals, we demonstrate both computational efficacy and high accuracy of the resulting MBE-CASSCF method. We assess the performance of our implementation on a set of established numerical tests before applying MBE-CASSCF in the investigation of the triplet-quintet spin gap of iron(II) porphyrin with active spaces as large as 50 electrons in 50 orbitals.

Exploring the photochemistry of OAlOH: Photodissociation pathways and electronic spectra

This study was focused on the photochemistry of OAlOH and three possible pathways, which were studied with high-level multireference configuration interaction ab initio calculations. We computed cuts of the six-dimensional potential energy surfaces for the ground, the lowest singlet and triplet excited states, and probed the photodissociation mechanisms and the stabilities. The OAlOH electronic spectrum, with an energy reaching 7.15 eV, contained four prominent peaks. Photodissociation to AlO, OH, and AlOH constituted a plausible mechanism within the deep-UV range (λ = 250.4 nm). Our data indicated the photostability of OAlOH in the near-UV‒Vis region, so detection with laser-induced fluorescence is possible. Fluorescence and phosphorescence may occur upon excitation at 363.5 nm. The roles of OAlOH in the photochemical reactions of Al-bearing molecules in the upper atmosphere and VY Canis Majoris are discussed.

Self-Consistent Field Approach for the Variational Quantum Eigensolver: Orbital Optimization Goes Adaptive

We present a self-consistent field (SCF) approach within the adaptive derivative-assembled problem-tailored ansatz variational quantum eigensolver (ADAPT-VQE) framework for efficient quantum simulations of chemical systems on near-term quantum computers. To this end, our ADAPT-VQE-SCF approach combines the idea of generating an ansatz with a small number of parameters, resulting in shallow-depth quantum circuits with a direct minimization of an energy expression that is correct to second order with respect to changes in the molecular orbital basis. Our numerical analysis, including calculations for the transition-metal complex ferrocene [Fe (C5H5)2], indicates that convergence in the self-consistent orbital optimization loop can be reached without a considerable increase in the number of two-qubit gates in the quantum circuit by comparison to a VQE optimization in the initial molecular orbital basis. Moreover, the orbital optimization can be carried out simultaneously within each iteration of the ADAPT-VQE cycle. ADAPT-VQE-SCF thus allows us to implement a routine analogous to the complete active space SCF, a cornerstone of state-of-the-art computational chemistry, in a hardware-efficient manner on near-term quantum computers. Hence, ADAPT-VQE-SCF paves the way toward a paradigm shift for quantitative quantum-chemistry simulations on quantum computers by requiring fewer qubits and opening up for the use of large and flexible atomic orbital basis sets in contrast to earlier methods that are predominantly based on the idea of full active spaces with minimal basis sets.

Convergent ab initio analysis of the multi-channel HOBr + H reaction

High-level potential energy surfaces for three reactions of hypobromous acid with atomic hydrogen were computed at the CCSDTQ/CBS//CCSDT(Q)/complete basis set level of theory. Focal point analysis was utilized to extrapolate energies and gradients for energetics and optimizations, respectively. The H attack at Br and subsequent Br-O cleavage were found to proceed barrierlessly. The slightly submerged transition state lies -0.2 kcal mol-1 lower in energy than the reactants and produces OH and HBr. The two other studied reaction paths are the radical substitution to produce H2O and Br with a 4.0 kcal mol-1 barrier and the abstraction at hydrogen to produce BrO and H2 with an 11.2 kcal mol-1 barrier. The final product energies lie -37.2, -67.9, and -7.3 kcal mol-1 lower in energy than reactants, HOBr + H, for the sets of products OH + HBr, H2O + Br, and H2 + BrO, respectively. Additive corrections computed for the final energetics, particularly the zero-point vibrational energies and spin-orbit corrections, significantly impacted the final stationary point energies, with corrections up to 6.2 kcal mol-1.

Sterically Hindered Tetradentate [Pt(O^N^C^N)] Emitters with Radiative Decay Rates up to 5.3 × 105 s-1 for Phosphorescent Organic Light-Emitting Diodes with LT95 Lifetime over 9200 h at 1000 cd m-2

Described here are sterically hindered tetradentate [Pt(O^N^C^N)] emitters (Pt-1, Pt-2, and Pt-3) developed for stable and high-performance green phosphorescent organic light-emitting diodes (OLEDs). These Pt(II) emitters exhibit strong saturated green phosphorescence (λmax = 517-531 nm) in toluene and mCP thin films with emission quantum yields as high as 0.97, radiative rate constants (kr) as high as 4.4-5.3 × 105 s-1 and reduced excimer emission, and with a preferential horizontally oriented transition dipole ratio of up to 84%. Theoretical calculations show that p-(hetero)arene substituents at the periphery of the ligand scaffolds in Pt-1, Pt-2, and Pt-3 can i) enhance the spin-orbit coupling (SOC) between the lower singlet excited states and the T1 state, and S0→Sn (n = 1 or 2) transition dipole moment, and ii) introducing additional SOC activity and the bright 1ILCT[π(carbazole)→π*(N^C^N)] excited state (Pt-2 and Pt-3), which are the main contributors to the increased kr values. Utilizing these tetradentate Pt(II) emitters, green phosphorescent OLEDs are fabricated with narrow-band electroluminescence (FWHM down to 36 nm), high external quantum efficiency, current efficiency up to 27.6% and 98.7 cd A-1, and an unprecedented device lifetime (LT95) of up to 9270 h at 1000 cd m-2 under laboratory conditions.

Economical quasi-Newton unitary optimization of electronic orbitals

We present an efficient quasi-Newton orbital solver optimized to reduce the number of gradient evaluations and other computational steps of comparable cost. The solver optimizes orthogonal orbitals by sequences of unitary rotations generated by the (preconditioned) limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm equipped with trust-region step restriction. The low-rank structure of the L-BFGS inverse Hessian is exploited when solving the trust-region problem. The efficiency of the proposed "Quasi-Newton Unitary Optimization with Trust-Region" (QUOTR) solver is compared to that of the standard Roothaan-Hall approach accelerated by the Direct Inversion of Iterative Subspace (DIIS), and other exact and approximate Newton solvers for mean-field (Hartree-Fock and Kohn-Sham) problems.

Atmospheric Gas-Phase Formation of Methanesulfonic Acid

Despite its impact on the climate, the mechanism of methanesulfonic acid (MSA) formation in the oxidation of dimethyl sulfide (DMS) remains unclear. The DMS + OH reaction is known to form methanesulfinic acid (MSIA), methane sulfenic acid (MSEA), the methylthio radical (CH3S), and hydroperoxymethyl thioformate (HPMTF). Among them, HPMTF reacts further to form SO2 and OCS, while the other three form the CH3SO2 radical. Based on theoretical calculations, we find that the CH3SO2 radical can add O2 to form CH3S(O)2OO, which can react further to form MSA. The branching ratio is highly temperature sensitive, and the MSA yield increases with decreasing temperature. In warmer regions, SO2 is the dominant product of DMS oxidation, while in colder regions, large amounts of MSA can form. Global modeling indicates that the proposed temperature-sensitive MSA formation mechanism leads to a substantial increase in the simulated global atmospheric MSA formation and burden.

Quantum Information-Assisted Complete Active Space Optimization (QICAS)

We propose an effective quantum information-assisted complete active space optimization scheme (QICAS). What sets QICAS apart from other correlation-based selection schemes is (i) the use of unique measures from quantum information that assess the correlation in electronic structures in an unambiguous and predictive manner and (ii) an orbital optimization step that minimizes the correlation discarded by the active space approximation. Equipped with these features, QICAS yields, for smaller correlated molecule, sets of optimized orbitals with respect to which the complete active space configuration interaction energy reaches the corresponding complete active space self-consistent field (CASSCF) energy within chemical accuracy. For more challenging systems such as the chromium dimer, QICAS offers an excellent starting point for CASSCF by greatly reducing the number of iterations required for numerical convergence. Accordingly, our study validates a profound empirical conjecture: the energetically optimal nonactive spaces are predominantly those that contain the least entanglement.

Excited-state van der Waals potential energy surfaces for the NO A2Σ+ + CO2X1Σg+ collision complex

Excited state van der Waals (vdW) potential energy surfaces (PESs) of the NO A2Σ+ + CO2X1Σg+ system are thoroughly investigated using coupled cluster theory and complete active space perturbation theory to second order (CASPT2). First, it is shown that pair natural orbital coupled cluster singles and doubles with perturbative triples yields comparable accuracy compared to CCSD(T) for molecular properties and vdW-minima at a fraction of computational cost of the latter. Using this method in conjunction with highly diffuse basis sets and counterpoise correction for basis set superposition error, the PESs for different intermolecular orientations are investigated. These show numerous vdW-wells, interconnected for all geometries except one, with a maximum depth of up to 830 cm-1; considerably deeper than those on the ground state surface. Multi-reference effects are investigated with CASPT2 calculations. The long-range vdW-surfaces support recent experimental observations relating to rotational energy transfer due the anisotropy in the potentials.

Excited Electronic States of Sr2: Ab Initio Predictions and Experimental Observation of the 21Σu+ State

Despite its apparently simple nature with four valence electrons, the strontium dimer constitutes a challenge for modern electronic structure theory. Here we focus on excited electronic states of Sr₂, which we investigate theoretically up to 25000 cm^-1 above the ground state, to guide and explain new spectroscopic measurements. In particular, we focus on potential energy curves for the 1¹Σ_u⁺, 2¹Σ_u⁺, 1¹Π_u, 2¹Π_u, and 1¹Δ_u states computed using several variants of ab initio coupled-cluster and configuration-interaction methods to benchmark them. In addition, a new experimental study of the excited 2¹Σ_u⁺ state using polarization labeling spectroscopy is presented, which extends knowledge of this state to high vibrational levels, where perturbation by higher electronic states is observed. The available experimental observations are compared with the theoretical predictions and help to assess the accuracy and limitations of employed theoretical models. The present results pave the way for future more accurate theoretical and experimental spectroscopic studies.

Excited States, Symmetry Breaking, and Unphysical Solutions in State-Specific CASSCF Theory

State-specific electronic structure theory provides a route toward balanced excited-state wave functions by exploiting higher-energy stationary points of the electronic energy. Multiconfigurational wave function approximations can describe both closed- and open-shell excited states and avoid the issues associated with state-averaged approaches. We investigate the existence of higher-energy solutions in complete active space self-consistent field (CASSCF) theory and characterize their topological properties. We demonstrate that state-specific approximations can provide accurate higher-energy excited states in H₂ (6-31G) with more compact active spaces than would be required in a state-averaged formalism. We then elucidate the unphysical stationary points, demonstrating that they arise from redundant orbitals when the active space is too large or symmetry breaking when the active space is too small. Furthermore, we investigate the singlet-triplet crossing in CH₂ (6-31G) and the avoided crossing in LiF (6-31G), revealing the severity of root flipping and demonstrating that state-specific solutions can behave quasi-diabatically or adiabatically. These results elucidate the complexity of the CASSCF energy landscape, highlighting the advantages and challenges of practical state-specific calculations.

Kylin 1.0: An ab-initio density matrix renormalization group quantum chemistry program

The accurate evaluation of electron correlations is highly necessary for the proper descriptions of the electronic structures in strongly correlated molecules, ranging from bond-dissociating molecules, polyradicals, to large conjugated molecules and transition metal complexes. For this purpose, in this paper, a new ab-initio quantum chemistry program Kylin 1.0 for electron correlation calculations at various quantum many-body levels, including configuration interaction (CI), perturbation theory (PT), and density matrix renormalization group (DMRG), is presented. Furthermore, fundamental quantum chemistry methods such as Hartree-Fock self-consistent field (HF-SCF) and the complete active space SCF (CASSCF) are also implemented. The Kylin 1.0 program possesses the following features: (1) a matrix product operator (MPO) formulation-based efficient DMRG implementation for describing static electron correlation within a large active space composed of more than 100 orbitals, supporting both U 1 n × U 1 S z $$ \mathrm{U}{(1)}_{\mathrm{n}}\times \mathrm{U}{(1)}_{{\mathrm{S}}_{\mathrm{z}}} $$ and U 1 n × SU 2 S $$ \mathrm{U}{(1)}_{\mathrm{n}}\times \mathrm{SU}{(2)}_{\mathrm{S}} $$ symmetries; (2) an efficient second-order DMRG-self-consistent field (SCF) implementation; (3) an externally contracted multi-reference CI (MRCI) and Epstein-Nesbet PT with DMRG reference wave functions for including the remaining dynamic electron correlation outside the large active spaces. In this paper, we introduce the capabilities and numerical benchmark examples of the Kylin 1.0 program.

Second-Order Self-Consistent Field Algorithms: From Classical to Quantum Nuclei

This work presents a general framework for deriving exact and approximate Newton self-consistent field (SCF) orbital optimization algorithms by leveraging concepts borrowed from differential geometry. Within this framework, we extend the augmented Roothaan-Hall (ARH) algorithm to unrestricted electronic and nuclear-electronic calculations. We demonstrate that ARH yields an excellent compromise between stability and computational cost for SCF problems that are hard to converge with conventional first-order optimization strategies. In the electronic case, we show that ARH overcomes the slow convergence of orbitals in strongly correlated molecules with the example of several iron-sulfur clusters. For nuclear-electronic calculations, ARH significantly enhances the convergence already for small molecules, as demonstrated for a series of protonated water clusters.

Photochemistry of 1-Phenyl-1-diazopropane and Its Diazirine Isomer: A CASSCF and MS-CASPT2 Study

In this work, we studied the wavelength (520 or 350 nm) dependence of the photochemical decomposition of 1-phenyl-1-diazopropane (PDP) and 1-phenyl-1-propyl diazirine (PED) by means of high-level ab initio quantum chemical calculations (CASSCF and MS-CASPT2) to obtain qualitative and quantitative results. It is found that the photochemistry of PDP is governed by nonradiative deactivation processes that can involve one or two S1/S0 conical intersections (CI1 and CI2) depending on the wavelength of the radiation; CI2 is only accessible at the shortest wavelength. It is demonstrated that the main intermediate of the photochemistry of the titled compounds is 1-ethyl-1-phenyl carbene (EPC). Upon irradiation of PDP with the 520 nm light, the carbene is always generated in its ground state as closed-shell singlet carbene. In contrast, the 350 nm radiation can directly decompose PDP into S1 carbene (open shell) and N2 when the conical intersection CI2 is avoided. Once the carbene is formed in the S1 state, it can experience excited state intramolecular proton transfer along a seam of crossing (ESIPT-SC) of the S1 and S0 states to yield the alkene derivative; that is, the proton transfer reaction takes places on a degenerate potential energy surface where the two electronic states have equal energy. In addition, it is found that EPC absorbs at 350 nm (double excitations); therefore, there is another possible route that can induce as well a slightly different photochemistry in changing the wavelength of the radiation because the shortest wavelength (when it is intense enough) decreases the amount of available EPC or generates a highly vibrationally excited state of the carbene; that is, after 350 nm excitation, the carbene intermediate can deactivate via radiation emission or can decay through a cascade of conical intersections to its first excited state (S1), where ESIPT-SC is operative again.

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

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

Accurate computed singlet-triplet energy differences for cobalt systems: implication for two-state reactivity

Accurate singlet-triplet energy differences for cobalt and rhodium complexes were calculated by using several wave function methods, such as MRCISD, CASPT2, CCSD(T) and BCCD(T). Relaxed energy differences were obtained by considering the singlet and triplet complexes, each at the minimum of their potential energy surfaces. Active spaces for multireference calculations were carefully checked to provide accurate results. The considered systems are built by increasing progressively the first coordination sphere around the metal. We included in our set two CpCoX complexes (Cp = cyclopentadienyl, X = alkenyl ligand), which have been suggested as intermediates in cycloaddition reactions. Indeed, cobalt systems have been used for more than a decade as active species in this kind of transformations, for which a two-state reactivity has been proposed. Most of the considered systems display a triplet ground state. However, in the case of a reaction intermediate, while a triplet ground state was predicted on the basis of Density Functional Theory results, our calculations suggest a singlet ground state. This stems from the competition between the exchange term (stabilising the triplet) and the accessibility of an intramolecular coordination (stabilising the singlet). This finding has an impact on the general mechanism of the cycloaddition reaction. Analogous rhodium systems were also studied and, as expected, they have a larger tendency to electron pairing than cobalt species.

Near-Exact Nuclear Gradients of Complete Active Space Self-Consistent Field Wave Functions

In this paper, we study the nuclear gradients of heat bath configuration interaction self-consistent field (HCISCF) wave functions and use them to optimize molecular geometries for various molecules.We show that the HCISCF nuclear gradients are fairly insensitive to the size of the "selected" variational space, which allows us to reduce the computational cost without introducing significant error.The ability of HCISCF to treat larger active spaces combined with the flexibility for users to control the computational cost makes the method very attractive for studying strongly correlated systems which require a larger active space than possible with complete active space self-consistent field (CASSCF).Finally, we study the realistic catalyst, Fe(PDI), and highlight some of the challenges this system poses for density functional theory (DFT).We demonstrate how HCISCF can clarify the energetic stability of geometries obtained from DFT when the results are strongly dependent on the functional.We also use the HCISCF gradients to optimize geometries for this species and study the adiabatic singlet-triplet gap. During geometry optimization, we find that multiple near-degenerate local minima exist on the triplet potential energy surface.

A combined first- and second-order optimization method for improving convergence of Hartree–Fock and Kohn–Sham calculations

In this work, we investigate the optimization of Hartree–Fock (HF) orbitals with our recently proposed combined first- and second-order (SO-SCI) method, which was originally developed for multi-configuration self-consistent field (MCSCF) and complete active space SCF (CASSCF) calculations. In MCSCF/CASSCF, it unites a second-order optimization of the active orbitals with a Fock-based first-order treatment of the remaining closed-virtual orbital rotations. In the case of the single-determinant wavefunctions, the active space is replaced by a preselected “second-order domain,” and all rotations involving orbitals in this subspace are treated at second-order. The method has been implemented for spin-restricted and spin-unrestricted Hartree–Fock (RHF, UHF), configuration-averaged Hartree–Fock (CAHF), as well as Kohn–Sham (KS) density functional theory (RKS, UKS). For each of these cases, various choices of the second-order domain have been tested, and appropriate defaults are proposed. The performance of the method is demonstrated for several transition metal complexes. It is shown that the SO-SCI optimization provides faster and more robust convergence than the standard SCF procedure but requires, in many cases, even less computation time. In difficult cases, the SO-SCI method not only speeds up convergence but also avoids convergence to saddle-points. Furthermore, it helps to find spin-symmetry broken solutions in the cases of UHF or UKS. In the case of CAHF, convergence can also be significantly improved as compared to a previous SCF implementation. This is particularly important for multi-center cases with two or more equal heavy atoms. The performance is demonstrated for various two-center complexes with different lanthanide atoms.

Dipolar 1,3-cycloaddition of thioformaldehyde S-methylide (CH2 SCH2 ) to ethylene and acetylene. A comparison with (valence) isoelectronic O3 , SO2 , CH2 OO and CH2 SO

Methods rooted in the density functional theory and in the coupled cluster ansatz were employed to investigate the cycloaddition reactions to ethylene and acetylene of 1,3-dipolar species including ozone and the derivatives issued from replacement of the central oxygen atom by the valence-isoelectronic sulfur atom, and/or of one or both terminal oxygen atoms by the isoelectronic CH₂ group. This gives rise to five different 1,3-dipolar compounds, namely ozone itself (O₃ ), sulfur dioxide (SO₂ ), the simplest Criegee intermediate (CH₂ OO), sulfine (CH₂ SO), and thioformaldehyde S-methylide (CH₂ SCH₂ , TSM). The experimental and accurate theoretical data available for some of those molecules were employed to assess the accuracy of two last-generation composite methods employing conventional or explicitly correlated post-Hartree-Fock contributions (jun-Cheap and SVECV-f12, respectively), which were then applied to investigate the reactivity of TSM. The energy barriers provided by both composite methods are very close (the average values for the two composite methods are 7.1 and 8.3 kcal mol^-1 for the addition to ethylene and acetylene, respectively) and comparable to those ruling the corresponding additions of ozone (4.0 and 7.7 kcal mol^-1 , respectively). These and other evidences strongly suggest that, at least in the case of cycloadditions, the reactivity of TSM is similar to that of O₃ and very different from that of SO₂ .

A trust-region augmented Hessian implementation for state-specific and state-averaged CASSCF wave functions

In this work, we present a one-step second-order converger for state-specific (SS) and state-averaged (SA) complete active space self-consistent field (CASSCF) wave functions. Robust convergence is achieved through step restrictions using a trust-region augmented Hessian (TRAH) algorithm. To avoid numerical instabilities, an exponential parameterization of variational configuration parameters is employed, which works with a nonredundant orthogonal complement basis. This is a common approach for SS-CASSCF and is extended to SA-CASSCF wave functions in this work. Our implementation is integral direct and based on intermediates that are formulated in either the sparse atomic-orbital or small active molecular-orbital basis. Thus, it benefits from a combination with efficient integral decomposition techniques, such as the resolution-of-the-identity or the chain-of-spheres for exchange approximations. This facilitates calculations on large molecules, such as a Ni(II) complex with 231 atoms and 5154 basis functions. The runtime performance of TRAH-CASSCF is competitive with the other state-of-the-art implementations of approximate and full second-order algorithms. In comparison with a sophisticated first-order converger, TRAH-CASSCF calculations usually take more iterations to reach convergence and, thus, have longer runtimes. However, TRAH-CASSCF calculations still converge reliably to a true minimum even if the first-order algorithm fails.

Photodissociation dynamics of N3+

The photodissociation dynamics of [Formula: see text] excited from its (linear) 3[Formula: see text]/(bent) 3A″ ground to the first excited singlet and triplet states is investigated. Three-dimensional potential energy surfaces for the 1A′, 1A″, and 3A′ electronic states, correlating with the 1Δg and 3Πu states in linear geometry, for [Formula: see text] are constructed using high-level electronic structure calculations and represented as reproducing kernels. The reference ab initio energies are calculated at the MRCI+Q/aug-cc-pVTZ level of theory. For following the photodissociation dynamics in the excited states, rotational and vibrational distributions P( v′) and P( j′) for the N2 product are determined from vertically excited ground state distributions. Due to the different shapes of the ground state 3A″ potential energy surface and the excited states, appreciable angular momentum j′ ∼ 60 is generated in diatomic fragments. The lifetimes in the excited states extend to at least 50 ps. Notably, results from sampling initial conditions from a thermal ensemble and from the Wigner distribution of the ground state wavefunction are comparable.

Importance of dynamical electron correlation in diabatic couplings of electron-exchange processes

We demonstrate the importance of the dynamical electron correlation effect in diabatic couplings of electron-exchange processes in molecular aggregates. To perform a multireference perturbation theory with large active space of molecular aggregates, an efficient low-rank approximation is applied to the complete active space self-consistent field reference functions. It is known that kinetic rates of electron-exchange processes, such as singlet fission, triplet-triplet annihilation, and triplet exciton transfer, are not sufficiently explained by the direct term of the diabatic couplings but efficiently mediated by the low-lying charge transfer states if the two molecules are in close proximity. It is presented in this paper, however, that regardless of the distance of the molecules, the direct term is considerably underestimated by up to three orders of magnitude without the dynamical electron correlation, i.e., the diabatic states expressed in the active space are not adequate to quantitatively reproduce the electron-exchange processes.

State-averaged CASSCF with polarizable continuum model for studying photoreactions in solvents: Energies, analytical nuclear gradients, and non-adiabatic couplings

This paper presents state-averaged complete active space self-consistent field in polarizable continuum model (PCM) for studies of photoreactions in solvents. The wavefunctions of the solute and the PCM surface charges of the solvent are optimized simultaneously such that the state-averaged free energy is variationally minimized. The method supports both fixed weights and dynamic weights where the weights are automatically adjusted based on the energy gaps. The corresponding analytical nuclear gradients and non-adiabatic couplings are also derived. Furthermore, we show how the new method can be entirely formulated in terms of seven basic operations, which allows the implementation to benefit from existing high-performance libraries on graphical processing units. Results demonstrating the accuracy and performance of the implementation are presented and discussed. We also apply the new method to the study of minimal conical intersection search and photoreaction energy pathways in solvents. Effects from the polarity of the solvents and different formulas of dynamic weights are compared and discussed.

Iterative subspace algorithms for finite-temperature solution of Dyson equation

One-particle Green’s functions obtained from the self-consistent solution of the Dyson equation can be employed in the evaluation of spectroscopic and thermodynamic properties for both molecules and solids. However, typical acceleration techniques used in the traditional quantum chemistry self-consistent algorithms cannot be easily deployed for the Green’s function methods because of a non-convex grand potential functional and a non-idempotent density matrix. Moreover, the optimization problem can become more challenging due to the inclusion of correlation effects, changing chemical potential, and fluctuations of the number of particles. In this paper, we study acceleration techniques to target the self-consistent solution of the Dyson equation directly. We use the direct inversion in the iterative subspace (DIIS), the least-squared commutator in the iterative subspace (LCIIS), and the Krylov space accelerated inexact Newton method (KAIN). We observe that the definition of the residual has a significant impact on the convergence of the iterative procedure. Based on the Dyson equation, we generalize the concept of the commutator residual used in DIIS and LCIIS and compare it with the difference residual used in DIIS and KAIN. The commutator residuals outperform the difference residuals for all considered molecular and solid systems within both GW and GF2. For a number of bond-breaking problems, we found that an easily obtained high-temperature solution with effectively suppressed correlations is a very effective starting point for reaching convergence of the problematic low-temperature solutions through a sequential reduction of temperature during calculations.

Nitrene formation is the first step of the thermal and photochemical decomposition reactions of organic azides

In this work, the decomposition of a prototypical azide, isopropyl azide, both in the ground and excited states, has been investigated through the use of multiconfigurational CASSCF and MS-CASPT2 electronic structure approaches. Particular emphasis has been placed on the thermal reaction starting at the S₀ ground state surface. It has been found that the azide thermally decomposes via a stepwise mechanism, whose rate-determining step is the formation of isopropyl nitrene, which is, in turn, the first step of the global mechanism. After that, the nitrene isomerizes to the corresponding imine derivative. Two routes are possible for such a decomposition: (i) a spin-allowed path involving a transition state; and (ii) a spin-forbidden one via a S₀/T₀ intersystem crossing. Both intermediates have been determined and characterised. Their associated relative energies have been found to be quite similar, 45.75 and 45.52 kcal mol^-1, respectively. To complete this study, the kinetics of the singlet and triplet channels are modeled with the MESMER (Master Equation Solver for Multi-Energy Well Reactions) code by applying the RRKM and Landau-Zener (with WKB tunnelling correction) theories, respectively. It is found that the canonical rate-coefficients of the singlet path are 2-orders of magnitude higher than the rate-coefficients of the forbidden reaction. In addition, the concerted mechanism has been investigated that would lead to the formation of the imine derivative and nitrogen extrusion in the first step of the decomposition. After a careful analysis of CASSCF calculations with different active spaces and their comparison with single electronic configuration methods (MP2 and B3LYP), the concerted mechanism is discarded.

iCISCF: An Iterative Configuration Interaction-Based Multiconfigurational Self-Consistent Field Theory for Large Active Spaces

An iterative configuration interaction (iCI)-based multiconfigurational self-consistent field (SCF) theory, iCISCF, is proposed to handle systems that require large active spaces. The success of iCISCF stems from three ingredients: (1) efficient selection of individual configuration state functions spanning the active space while maintaining full spin symmetry; (2) the use of Jacobi rotation for optimization of the active orbitals in conjunction with a quasi-Newton algorithm for the core/active-virtual and core-active orbital rotations; (3) a second-order perturbative treatment of the residual space left over by the selection procedure (i.e., iCISCF(2)). Several examples that go beyond the capability of CASSCF are taken as showcases to reveal the efficacy of iCISCF and iCISCF(2), facilitated by iCAS for imposed automatic selection and localization of active orbitals.

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.