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

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

The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry.

The developments of the open-source OpenMolcas 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 OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Application of Multiconfiguration Pair-Density Functional Theory to the Diels-Alder Reaction.

Transition states for Diels-Alder reactions are strongly correlated, as evidenced by high-to-very-high M diagnostics, and therefore they require treatment by multireference methods. Multiconfiguration pair-density functional theory (MC-PDFT) combines a multiconfiguration wave function with a functional of the electron density and the on-top pair density to calculate the electronic energy for strongly correlated systems at a much lower cost than wave function methods that do not employ density functionals. Here we apply MC-PDFT to the Diels-Alder cycloaddition reaction of 1,3-butadiene with ethylene, where two kinds of reaction paths have been widely studied: concerted synchronous paths and diradical stepwise paths. The lowest-energy reaction path is now known to be a concerted synchronous one, and a method's ability to predict this is an important test. By comparison to the best available theoretical results in the literature, we test the accuracy of MC-PDFT with several choices of on-top functional for geometries and enthalpies of stable structures along both paths and for the transition state geometries. We also calculate the Arrhenius activation energies for both paths and compare these to experiment. We also compare to Kohn-Sham density functional theory (KS-DFT) with selected exchange-correlation functionals. CAS-PDFT gives consistently good energies and geometries for both the concerted and stepwise mechanisms, but none of the KS-DFT functionals gives accurate activation energies for both. The stepwise transition state is very strongly correlated, and MC-PDFT can treat it, but KS-DFT (which involves a single-configuration treatment) has larger errors. The results confirm that using a multiconfigurational reference function for strongly correlated transition states can significantly improve the reliability and that MC-PDFT can provide good accuracy at a much lower computational cost than competing multireference methods.

Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory.

Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn-Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin-orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals.

Zero-Field Splitting Calculations by Multiconfiguration Pair-Density Functional Theory.

Zero-field splitting (ZFS) is a fundamental molecular property that is especially relevant for single-molecule magnets (SMMs), electron paramagnetic resonance spectra, and quantum computing. Developing a method that can accurately predict ZFS parameters can be very powerful for designing new SMMs. One of the challenges is to include external correlation in an inherently multiconfigurational open-shell species for the accurate prediction of magnetic properties. Previously available methods depend on expensive multireference perturbation theory calculations to include external correlation. In this paper, we present spin-orbit-inclusive multiconfiguration and multistate pair-density functional theory (MC-PDFT) calculations of ZFSs; these calculations have a cost comparable to complete-active-space self-consistent field (CASSCF) theory, but they include correlation external to the active space. We found that combining a multistate formulation of MC-PDFT, namely, compressed-state multistate pair-density functional theory, with orbitals optimized by weighted-state-averaged CASSCF, yields reasonably accurate ZFS results.

Multiconfiguration Pair-Density Functional Theory.

Kohn-Sham density functional theory with the available exchange-correlation functionals is less accurate for strongly correlated systems, which require a multiconfigurational description as a zero-order function, than for weakly correlated systems, and available functionals of the spin densities do not accurately predict energies for many strongly correlated systems when one uses multiconfigurational wave functions with spin symmetry. Furthermore, adding a correlation functional to a multiconfigurational reference energy can lead to double counting of electron correlation. Multiconfiguration pair-density functional theory (MC-PDFT) overcomes both obstacles, the second by calculating the quantum mechanical part of the electronic energy entirely by a functional, and the first by using a functional of the total density and the on-top pair density rather than the spin densities. This allows one to calculate the energy of strongly correlated systems efficiently with a pair-density functional and a suitable multiconfigurational reference function. This article reviews MC-PDFT and related background information.

Compressed-State Multistate Pair-Density Functional Theory.

Multiconfiguration pair-density functional theory (MC-PDFT) is a multireference method that can be used to calculate excited states. However, MC-PDFT potential energy surfaces have the wrong topology at conical intersections because the last step of MC-PDFT is not a diagonalization of a model-space Hamiltonian matrix, as done in, for example, multistate second-order perturbation theory (MS-CASPT2). We have previously proposed methods that solve this problem by diagonalizing a model-space effective Hamiltonian matrix, where the diagonal elements are MC-PDFT energies for intermediate states, and the off-diagonal elements are evaluated by wave function theory. One previous method is called variational multistate PDFT (VMS-PDFT), whose intermediate states maximize the trace of the effective Hamiltonian, namely, the sum of the MC-PDFT energies of the model-space states; the VMS-PDFT is very robust but is more computationally expensive than another method, extended multistate PDFT (XMS-PDFT), in which the transformation to intermediate states is accomplished without needing any density functional evaluations. However, although VMS-PDFT was accurate in all cases tested, XMS-PDFT was accurate in only some of them. In the present paper, we propose a new method, called compressed-state multistate PDFT (CMS-PDFT), that is as efficient as XMS-PDFT and as accurate as VMS-PDFT. The new method maximizes the trace of the classical Coulomb energy of the intermediate states such that the electron densities of the intermediate states are compressed. We show that CMS-PDFT performs robustly even where XMS-PDFT fails.

Multi-state pair-density functional theory.

Multi-configuration pair-density functional theory (MC-PDFT) has previously been applied successfully to carry out ground-state and excited-state calculations. However, because they include no interaction between electronic states, MC-PDFT calculations in which each state's PDFT energy is calculated separately can give an unphysical double crossing of potential energy surfaces (PESs) in a region near a conical intersection. We have recently proposed state-interaction pair-density functional theory (SI-PDFT) to treat nearly degenerate states by creating a set of intermediate states with state interaction; although this method is successful, it is inconvenient because two SCF calculations and two sets of orbitals are required and because it puts the ground state on an unequal footing with the excited states. Here we propose two new methods, called extended-multi-state-PDFT (XMS-PDFT) and variational-multi-state-PDFT (VMS-PDFT), that generate the intermediate states in a balanced way with a single set of orbitals. The former uses the intermediate states proposed by Granovsky for extended multi-configuration quasi-degenerate perturbation theory (XMC-QDPT); the latter obtains the intermediate states by maximizing the sum of the MC-PDFT energies for the intermediate states. We also propose a Fourier series expansion to make the variational optimizations of the VMS-PDFT method convenient, and we implement this method (FMS-PDFT) both for conventional configuration-interaction solvers and for density-matrix-renormalization-group solvers. The new methods are tested for eight systems, exhibiting avoided crossings among two to six states. The FMS-PDFT method is successful for all cases for which it has been tested (all cases in this paper except O3 for which it was not tested), and XMS-PDFT is successful for all eight cases except the mixed-valence case. Since both XMS-PDFT and VMS-PDFT are less expensive than XMS-CASPT2, they will allow well-correlated calculations on much larger systems for which perturbation theory is unaffordable.

OpenMolcas: From Source Code to Insight.

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 postcalculation analysis and visualization, a multiscale simulation option using frozen-density embedding theory, and new electronic and muonic basis sets.

Automatic Active Space Selection for Calculating Electronic Excitation Energies Based on High-Spin Unrestricted Hartree-Fock Orbitals.

Multireference perturbation theory (MR-PT2, e.g., MS-CASPT2) and multiconfiguration pair-density functional theory (MC-PDFT) use a multiconfigurational wave function as the reference wave function, and this can be generated, for example, by the complete active space self-consistent field (CASSCF) method or restricted active space self-consistent field (RASSCF) method. The MR-PT2 and MC-PDFT methods have proved successful in many previous studies, but their performance depends on the quality of the reference wave function, and the quality of a CASSCF or RASSCF wave function depends on the active space. Even for a given number of active electrons and active orbitals, the orbitals obtained at the end of the self-consistent field iterations of the CASSCF calculation may be different for different sets of initial guess orbitals. Consequently, it is a worthwhile goal to automate active space selection, including the choice of orbitals to start the iterations, and here we examine the question of whether we can devise a broadly applicable automatic scheme for producing guess orbitals that lead to an active space that gives accurate excitation energies. Such a scheme depends on the target number of excitation energies that one is trying to calculate, and here our target is the first two spin-conserving excited electronic states of a closed-shell molecule or a doublet radical. For this target, we propose a scheme called ABC2 to automatically select an active space for MC-PDFT calculations and MS-CASPT2 calculations; the scheme uses high-spin-state unrestricted Hartree-Fock (UHF) natural orbitals as guess orbitals. A novel feature of the ABC2 scheme is that it contains a reliability test that should be passed for the results to be considered reliable. This test requires that the MC-PDFT or CASPT2 excitation energy be within 1.1 eV of the CASSCF excitation energy for the result to be considered reliable. We evaluated the performance of ABC2 for the two lowest excitation energies (or only the lowest when an accurate value of the second lowest state is not available) of 16 singlet systems and 10 doublet systems. We compared the performance of ABC2 to that of our previous ABC scheme, which finds active orbitals through a more expensive (multiconfigurational) excited-state calculation, to the use of ground-state UHF natural orbitals, and to the default scheme in OpenMolcas. We found that the new scheme is more robust than previous methods for a variety of systems.

Weak Interactions in Alkaline Earth Metal Dimers by Pair-Density Functional Theory.

Alkaline earth dimers have small bond energies (less than 5 kcal/mol) that provide a difficult challenge for electronic structure calculations. They are especially challenging for Kohn-Sham density functional theory (KS-DFT) using generalized gradient approximations (GGAs) as the exchange-correlation density functional because GGAs often do not provide accurate results for weak interactions. Here we treat alkaline earth dimers from six different rows of the periodic table. We show that the dominant correlating configurations are the same in all six dimers. We also show that multiconfiguration pair-density functional theory (MC-PDFT) using a fully translated GGA as the on-top density functional not only performs much better than KS-DFT with GGAs in predicting equilibrium distances and dissociation energies but also performs better than the more computationally demanding complete active space second-order perturbation theory (CASPT2) with large basis sets and performs even better than CASPT2 with smaller basis sets.

Multireference Methods for Calculating the Dissociation Enthalpy of Tetrahedral P4 to Two P2.

The potential energy surface for the thermal decomposition of P4 → 2P2 was computed along the C2 v reaction trajectory. Single-reference methods were not suitable for describing this complex bond-breaking process, so two multiconfigurational methods, namely, multistate complete active space second-order perturbation theory (MS-CASPT2) and multiconfiguration pair-density functional theory (MC-PDFT), were used with the aim of determining the accuracy and efficiency of these methods for this process. Several active spaces and basis sets were explored. It was found that the multiconfiguration pair-density functional theory method was up to 900 times faster than multistate complete active space second-order perturbation theory while providing similar accuracy.

Automatic Selection of an Active Space for Calculating Electronic Excitation Spectra by MS-CASPT2 or MC-PDFT.

Multireference methods such as multistate complete active space second-order perturbation theory (MS-CASPT2) and multiconfiguration pair-density functional theory (MC-PDFT) can be very accurate, but they have the disadvantage that they are not black-box methods, and finding a good active space for the reference wave function often requires nonsystematic procedures based on intimate knowledge of the system under study. In this paper, we propose a scheme called the ABC scheme, which is a three-step calculation using three parameters, for automatic selection (without looking at the orbitals and without using any knowledge of the specific system at hand) of the active space (including both the size of the active space and the orbitals in the active space) for MS-CASPT2 or MC-PDFT calculations, and we report tests of the scheme on the first five excitation energies for a set of ten doublet radicals. The results show that the ABC scheme is very successful for both MS-CASPT2 and MC-PDFT for all ten systems considered here.

Multiconfiguration pair-density functional theory for doublet excitation energies and excited state geometries: the excited states of CN.

Multiconfiguration pair-density functional theory (MC-PDFT) is a post multiconfiguration self-consistent field (MCSCF) method with similar performance to complete active space second-order perturbation theory (CASPT2) but with greater computational efficiency. Cyano radical (CN) is a molecule whose spectrum is well established from experiments and whose excitation energies have been used as a testing ground for theoretical methods to treat excited states of open-shell systems, which are harder and much less studied than excitation energies of closed-shell singlets. In the present work, we studied the adiabatic excitation energies of CN with MC-PDFT. Then we compared this multireference (MR) method to some single-reference (SR) methods, including time-dependent density functional theory (TDDFT) and completely renormalized equation-of-motion coupled-cluster theory with singles, doubles and noniterative triples [CR-EOM-CCSD(T)]; we also compared to some other MR methods, including configuration interaction singles and doubles (MR-CISD) and multistate CASPT2 (MS-CASPT2). Through a comparison between SR and MR methods, we achieved a better appreciation of the need to use MR methods to accurately describe higher excited states, and we found that among the MR methods, MC-PDFT stands out for its accuracy for the first four states out of the five doublet states studied this paper; this shows its efficiency for calculating doublet excited states.