In this article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already includes a large number of new developments realized during the transition from the commercial MOLCAS product to the open-source platform. The paper initially describes the technical details of the new software development platform. This is followed by brief presentations of many new methods, implementations, and features of the OpenMolcas program suite. These developments include novel wave function methods such as stochastic complete active space self-consistent field, density matrix renormalization group (DMRG) methods, and hybrid multiconfigurational

Predicting the singlet-triplet splittings of divalent radicals is a challenging task for electronic structure theory. In the present work, we investigate the performance of multiconfiguration pair-density functional theory (MC-PDFT) for computing the singlet-triplet splitting for small main-group divalent radicals for which accurate experimental data are available. In order to define theoretical model chemistries that can be assessed consistently, we define three correlated participating orbitals (CPO) schemes (nominal, moderate, and extended, abbreviated as nom, mod, and ext) to define the constitution of complete active spaces, and we test them systematically. Broken-symmetry Kohn-Sham DFT calculations have also been carried out for comparison. We found that the extended CPO-PDFT scheme with translated on-top pair-density functionals have smaller mean unsigned errors than weighted-average broken-symmetry Kohn-Sham DFT with the corresponding exchange-correlation functional. The accuracy of the translated Perdew-Burke-Ernzerhof (tPBE) on-top pair-density functionals with ext-CPO active space is even better than some of the more accurately parametrized exchange-correlation density functionals that we tested; this is very encouraging for MC-PDFT theory.

Analytic gradient routines are a desirable feature for quantum mechanical methods, allowing for efficient determination of equilibrium and transition state structures and several other molecular properties. In this work, we present analytical gradients for multiconfiguration pair-density functional theory (MC-PDFT) when used with a state-specific complete active space self-consistent field reference wave function. Our approach constructs a Lagrangian that is variational in all wave function parameters.We find that MC-PDFT locates equilibrium geometries for several small-to mediumsized organic molecules that are similar to those located by complete active space second-order perturbation theory but that are obtained with decreased computational cost.

The developments of the open-source chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes can address, while showing that is an attractive platform for state-of-the-art atomistic computer simulations.

In this article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already includes a large number of new developments realized during the transition from the commercial MOLCAS product to the open-source platform. The paper initially describes the technical details of the new software development platform. This is followed by brief presentations of many new methods, implementations, and features of the OpenMolcas program suite. These developments include novel wave function methods such as stochastic complete active space self-consistent field, density matrix renormalization group (DMRG) methods, and hybrid multiconfigurational wave function and density functional theory models. Some of these implementations include an array of additional options and functionalities. The paper proceeds and describes developments related to explorations of potential energy surfaces. Here we present methods for the optimization of conical intersections, the simulation of adiabatic and nonadiabatic molecular dynamics and interfaces to tools for semiclassical and quantum mechanical nuclear dynamics. Furthermore, the article describes features unique to simulations of spectroscopic and magnetic phenomena such as the exact semiclassical description of the interaction between light and matter, various X-ray processes, magnetic circular dichroism and properties. Finally, the paper describes a number of built-in and add-on features to support the OpenMolcas platform with post calculation analysis and visualization, a multiscale simulation option using frozen-density embedding theory and new electronic and muonic basis sets.

Analytic gradients are important for efficient calculations of stationary points on potential energy surfaces, for interpreting spectroscopic observations, and for efficient direct dynamics simulations. For excited electronic states, as are involved in UV–Vis spectroscopy and photochemistry, analytic gradients are readily available and often affordable for calculations using a state-averaged complete active space self-consistent-field (SA-CASSCF) wave function. However, in most cases, a post-SA-CASSCF step is necessary for quantitative accuracy, and such calculations are often too expensive if carried out by perturbation theory or configuration interaction. In this work, we present the analytic gradients for multiconfiguration pair-density functional theory based on SA-CASSCF wave functions, which is a more affordable alternative. A test set of molecules has been studied with this method, and the stationary geometries and energetics are compared to values in the literature as obtained by other methods. Excited-state geometries computed with state-averaged pair-density functional theory have similar accuracy to those from complete active space perturbation theory at the second-order.

A novel approach to strongly contracted N-electron valence perturbation theory (SC-NEVPT2) as a means of describing dynamic electron correlation for quantum chemical density matrix renormaliza-tion group (DMRG) calculations is presented. In this approach the strongly contracted perturber functions are projected onto a renormalized Hilbert space. Compared to a straightforward implementation of SC-NEVPT2 with DMRG wavefunctions, the computational scaling and storage requirements are reduced. This favorable scaling opens up the possibility of calculations with larger active spaces. A specially designed renormalization scheme ensures that both the electronic ground state and the perturber functions are well represented in the renormalized Hilbert space. Test calculations on the N 2 and [Cu 2 O 2 (en) 2 ] 2+ demonstrate some key properties of the method and indicate its capabilities. Published by AIP Publishing. [http://dx.

The accurate description of ground- and excited-state potential energy surfaces poses a challenge for many electronic structure methods, especially in regions where strong electronic state interaction occurs. Here we introduce a new methodology, state-interaction pair-density functional theory (SI-PDFT), to target molecular systems exhibiting strong interaction of electronic states. SI-PDFT is an extension of multiconfiguration pair-density functional theory in which a set of N electronic states is generated through the diagonalization of an N × N effective Hamiltonian. We demonstrate the accuracy of the method by performing calculations on the ionic-neutral avoided crossing in lithium fluoride and the ππ-πσ avoided crossing in the H-O bond photodissociation in phenol. We show that SI-PDFT can be a useful tool in the study of photochemistry and nonadiabatic dynamics.