A hybrid meta on-top functional for multiconfiguration pair-density functional theory.

Multiconfiguration pair-density functional theory (MC-PDFT) was proposed a decade ago, but it is still in the early stage of density functional development. MC-PDFT uses functionals that are called on-top functionals; they depend on the density and the on-top pair density. Most MC-PDFT calculations to date have been unoptimized translations of generalized gradient approximations (GGAs) of Kohn-Sham density functional theory (KS-DFT).

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 and electronic structure packages. Yet, the main focus of this work is using ISCs to calculate minimum-energy conical intersections (MECIs) by CMS-PDFT.

The OpenMolcas : 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.

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.

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.

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.

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.

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.

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.

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.

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.

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

The potential energy surface for the thermal decomposition of P → 2P was computed along the C 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.

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.