Toward a Stochastic Complete Active Space Second-Order Perturbation Theory.

In this work, an internally contracted stochastic complete active space second-order perturbation theory, stochastic-CASPT2, is reported. The method relies on stochastically sampled reduced density matrices (RDMs) up to rank four and contractions thereof with the generalized Fock matrix. A new protocol for calculating higher-order RDMs in full configuration interaction quantum Monte Carlo (FCIQMC) has been designed based on (1) restricting sampling of the corresponding excitations to a deterministic subspace, (2) averaging the RDMs from independent dynamics and (3) projecting them onto the closest positive semi-definite matrix. Our protocol avoids previously encountered numerical conditioning problems in the orthogonalization of the perturber overlap matrix stemming from numerical noise. The chromium dimer CASSCF(12,12)/CASPT2 binding curve is computed as a proof of concept.

Exploiting Locality in Full Configuration Interaction Quantum Monte Carlo for Fast Excitation Generation.

In this paper, we propose an improved excitation generation algorithm for the full configuration interaction quantum Monte Carlo method, which is particularly effective in systems described by localized orbitals. The method is an extension of the precomputed heat-bath strategy of Holmes et al., with more effective sampling of double excitations and a novel approach for nonuniform sampling of single excitations. We demonstrate the effectiveness of the algorithm for a chain of 30 hydrogen atoms with atom-localized orbitals, a stack of benzene molecules, and an Fe(II)-porphyrin model complex, whereby we show an overall efficiency gain by a factor of two to four, as measured by variance reduction per wall-clock time.

NECI: N-Electron Configuration Interaction with an emphasis on state-of-the-art stochastic methods.

We present NECI, a state-of-the-art implementation of the Full Configuration Interaction Quantum Monte Carlo (FCIQMC) algorithm, a method based on a stochastic application of the Hamiltonian matrix on a sparse sampling of the wave function. The program utilizes a very powerful parallelization and scales efficiently to more than 24 000 central processing unit cores. In this paper, we describe the core functionalities of NECI and its recent developments. This includes the capabilities to calculate ground and excited state energies, properties via the one- and two-body reduced density matrices, as well as spectral and Green's functions for ab initio and model systems. A number of enhancements of the bare FCIQMC algorithm are available within NECI, allowing us to use a partially deterministic formulation of the algorithm, working in a spin-adapted basis or supporting transcorrelated Hamiltonians. NECI supports the FCIDUMP file format for integrals, supplying a convenient interface to numerous quantum chemistry programs, and it is licensed under GPL-3.0.

Separated-pair approximation and separated-pair pair-density functional theory.

Multi-configuration pair-density functional theory (MC-PDFT) has proved to be a powerful way to combine the capabilities of multi-configuration self-consistent-field theory to represent the an electronic wave function with a highly efficient way to include dynamic correlation energy by density functional theory. All applications reported previously involved complete active space self-consistent-field (CASSCF) theory for the reference wave function. For treating large systems efficiently, it is necessary to ask whether good accuracy is retained when using less complete configuration interaction spaces. To answer this question, we present here calculations employing MC-PDFT with the separated pair (SP) approximation, which is a special case (defined in this article) of generalized active space self-consistent-field (GASSCF) theory in which no more than two orbitals are included in any GAS subspace and in which inter-subspace excitations are excluded. This special case of MC-PDFT will be called SP-PDFT. In SP-PDFT, the electronic kinetic energy and the classical Coulomb energy, the electronic density and its gradient, and the on-top pair density and its gradient are obtained from an SP approximation wave function; the electronic energy is then calculated from the first two of these quantities and an on-top density functional of the last four. The accuracy of the SP-PDFT method for predicting the structural properties and bond dissociation energies of twelve diatomic molecules and two triatomic molecules is compared to the SP approximation itself and to CASSCF, MC-PDFT based on CASSCF, CASSCF followed by second order perturbation theory (CASPT2), and Kohn-Sham density functional theory with the PBE exchange-correlation potential. We show that SP-PDFT reproduces the accuracy of MC-PDFT based on the corresponding CASSCF wave function for predicting C-H bond dissociation energies, the reaction barriers of pericyclic reactions and the properties of open-shell singlet systems, all at only a small fraction of the computational cost.

Controversial electronic structures and energies of Fe2, Fe 2 (+), and Fe 2 (-) resolved by RASPT2 calculations.

The diatomic molecule Fe2 was investigated using restricted active space second-order perturbation theory (RASPT2). This molecule is very challenging to study computationally because predictions about the ground state and excited states depend sensitively on the choice of the quantum chemical method. For Fe2 we show that one needs to go beyond a full-valence active space in order to achieve even qualitative agreement with experiment for the dissociation energy, and we also obtain a smooth ground-state potential curve. In addition we report the first multireference study of Fe 2 (+), for which we predict an (8)Σu (-) ground state, which was not predicted by previous computational studies. By using an active space large enough to remove the most serious deficiencies of previous theoretical work and by explicitly investigating the interpretations of previous experimental results, this study elucidates previous difficulties and provides - for the first time - a qualitatively correct treatment of Fe2, Fe 2 (+), and Fe 2 (-). Moreover, this study represents a record in terms of the number or active electrons and active orbitals in the active space, namely 16 electrons in 28 orbitals. Conventional CASPT2 calculations can be performed with at most 16 electrons in 16 orbitals. We were able to overcome this limit by using the RASPT2 formalism.