Accurately modeling the electronic structure of systems with many unpaired electrons remains a major challenge in quantum chemistry. Qualitatively correct electronic structures generally require large active space multireference wave functions, while dynamic correlation effects beyond the active space are crucial for quantitatively accurate descriptions of magnetic, catalytic and optical properties of such systems. Here, we present an uncontracted multireference perturbation theory based on the FCIQMC imaginary-time evolution of effective Hamiltonians, built upon the generalized active space concept and Löwdin’s partitioning technique. The configurational interaction space is split into a reference space, consisting of the most important configurations, and a perturber space, containing the more numerous configurations responsible for dynamic correlation effects. The generalized active space algorithm allows the flexible partitioning of the configurational space. Löwdin’s partitioning technique is then used to construct an effective Hamiltonian which is stochastically solved. This strategy allows us to apply perturbative corrections on large active space reference wave functions, without requiring high-order reduced density matrices, which have been found the bottleneck in other perturbation theory strategies. The capabilities of the resulting method, called Stochastic-SplitGAS, are demonstrated on the triplet-quintet spin gap of an Fe(II)-porphyrin model system and the spin ladder of a [Fe(III)2S2]2– complex.
Accurately modeling the electronic structure of systems with many unpaired electrons remains a major challenge in quantum chemistry. Qualitatively correct electronic structures generally require large active space multireference wave functions, while dynamic correlation effects beyond the active space are crucial for quantitatively accurate descriptions of magnetic, catalytic and optical properties of such systems. Here, we present an uncontracted multireference perturbation theory based on the FCIQMC imaginary-time evolution of effective Hamiltonians, built upon the generalized active space concept and Löwdin’s partitioning technique. The configurational interaction space is split into a reference space, consisting of the most important configurations, and a perturber space, containing the more numerous configurations responsible for dynamic correlation effects. The generalized active space algorithm allows the flexible partitioning of the configurational space. Löwdin’s partitioning technique is then used to construct an effective Hamiltonian which is stochastically solved. This strategy allows us to apply perturbative corrections on large active space reference wave functions, without requiring high-order reduced density matrices, which have been found the bottleneck in other perturbation theory strategies. The capabilities of the resulting method, called Stochastic-SplitGAS, are demonstrated on the triplet-quintet spin gap of an Fe(II)-porphyrin model system and the spin ladder of a [Fe(III)2S2]2– complex.
An algorithm to perform stochastic generalized active space calculations, Stochastic-GAS, is presented, that uses the Slater determinant based FCIQMC algorithm as configuration interaction eigensolver. Stochastic-GAS allows the construction and stochastic optimization of preselected truncated configuration interaction wave functions, either to reduce the computational costs of large active space wave function optimizations, or to probe the role of specific electron correlation pathways. As for the conventional GAS procedure, the preselection of the truncated wave function is based on the selection of multiple active subspaces while imposing restrictions on the interspace excitations. Both local and cumulative minimum and maximum occupation number constraints are supported by Stochastic-GAS. The occupation number constraints are efficiently encoded in precomputed probability distributions, using the precomputed heat bath algorithm, which removes nearly all runtime overheads of GAS. This strategy effectively allows the FCIQMC dynamics to a priori exclude electronic configurations that are not allowed by GAS restrictions. Stochastic-GAS reduced density matrices are stochastically sampled, allowing orbital relaxations via Stochastic-GASSCF, and direct evaluation of properties that can be extracted from density matrices, such as the spin expectation value. Three test case applications have been chosen to demonstrate the flexibility of Stochastic-GAS: (a) the Stochastic-GASSCF optimization of a stack of five benzene molecules, that shows the applicability of Stochastic-GAS towards fragment-based chemical systems; (b) an uncontracted stochastic MRCISD calculation that correlates 96 electrons and 159 molecular orbitals, and uses a large (32, 34) active space reference wave function for an Fe(II)-porphyrin model system, showing how GAS can be applied to systematically recover dynamic electron correlation, and how in the specific case of the Fe(II)-porphyrin dynamic correlation further differentially stabilizes the triplet over the quintet spin state; (c) the study of an Fe4S4 cluster's spin-ladder energetics via highly truncated stochastic-GAS wave functions, where we show how GAS can be applied to understand the competing spin-exchange and charge-transfer correlating mechanisms in stabilizing different spin-states.
An algorithm to perform stochastic generalized active space calculations, Stochastic-GAS, is presented, that uses the Slater determinant based FCIQMC algorithm as configuration interaction eigensolver. Stochastic-GAS allows the construction and stochastic optimization of preselected truncated configuration interaction wave functions, either to reduce the computational costs of large active space wave function optimizations, or to probe the role of specific electron correlation pathways. As for the conventional GAS procedure, the preselection of the truncated wave function is based on the selection of multiple active subspaces while imposing restrictions on the interspace excitations. Both local and cumulative minimum and maximum occupation number constraints are supported by Stochastic-GAS. The occupation number constraints are efficiently encoded in precomputed probability distributions, using the precomputed heat bath algorithm, which removes nearly all runtime overheads of GAS. This strategy effectively allows the FCIQMC dynamics to a priori exclude electronic configurations that are not allowed by GAS restrictions. Stochastic-GAS reduced density matrices are stochastically sampled, allowing orbital relaxations via Stochastic-GASSCF, and direct evaluation of properties that can be extracted from density matrices, such as the spin expectation value. Three test case applications have been chosen to demonstrate the flexibility of Stochastic-GAS: (a) the Stochastic-GASSCF optimization of a stack of five benzene molecules, that shows the applicability of Stochastic-GAS towards fragment-based chemical systems; (b) an uncontracted stochastic MRCISD calculation that correlates 96 electrons and 159 molecular orbitals, and uses a large (32, 34) active space reference wave function for an Fe(II)-porphyrin model system, showing how GAS can be applied to systematically recover dynamic electron correlation, and how in the specific case of the Fe(II)-porphyrin dynamic correlation further differentially stabilizes the triplet over the quintet spin state; (c) the study of an Fe4S4 cluster's spin-ladder energetics via highly truncated stochastic-GAS wave functions, where we show how GAS can be applied to understand the competing spin-exchange and charge-transfer correlating mechanisms in stabilizing different spin-states.
An algorithm to perform stochastic generalized active space calculations, Stochastic-GAS, is presented, that uses the Slater determinant based FCIQMC algorithm as configuration interaction eigensolver. Stochastic-GAS allows the construction and stochastic optimization of preselected truncated configuration interaction wave functions, either to reduce the computational costs of large active space wave function optimizations, or to probe the role of specific electron correlation pathways. As for the conventional GAS procedure, the preselection of the truncated wave function is based on the selection of multiple active subspaces while imposing restrictions on the interspace excitations. Both local and cumulative minimum and maximum occupation number constraints are supported by Stochastic-GAS. The occupation number constraints are efficiently encoded in precomputed probability distributions, using the precomputed heat bath algorithm, which removes nearly all runtime overheads of GAS. This strategy effectively allows the FCIQMC dynamics to a priori exclude electronic configurations that are not allowed by GAS restrictions. Stochastic-GAS reduced density matrices are stochastically sampled, allowing orbital relaxations via Stochastic-GASSCF, and direct evaluation of properties that can be extracted from density matrices, such as the spin expectation value. Three test case applications have been chosen to demonstrate the flexibility of Stochastic-GAS: (a) the Stochastic-GASSCF optimization of a stack of five benzene molecules, that shows the applicability of Stochastic-GAS towards fragment-based chemical systems; (b) an uncontracted stochastic MRCISD calculation that correlates 96 electrons and 159 molecular orbitals, and uses a large (32, 34) active space reference wave function for an Fe(II)-porphyrin model system, showing how GAS can be applied to systematically recover dynamic electron correlation, and how in the specific case of the Fe(II)-porphyrin dynamic correlation further differentially stabilizes the triplet over the quintet spin state; (c) the study of an Fe4S4 cluster's spin-ladder energetics via highly truncated stochastic-GAS wave functions, where we show how GAS can be applied to understand the competing spin-exchange and charge-transfer correlating mechanisms in stabilizing different spin-states.
Full Configuration Interaction Quantum Monte Carlo (FCIQMC) has been effectively applied to very large configuration interaction (CI) problems and was recently adapted for use as an active space solver and combined with orbital optimization. In this work, we detail an approach within FCIQMC to allow for efficient sampling of fully internally contracted multireference perturbation theories within the same stochastic framework. Schemes are described to allow for the close control over the resolution of stochastic sampling of the effective higher-body intermediates within the active space. It is found that while complete active space second-order perturbation theory seems less amenable to a stochastic reformulation, strongly contracted N-Electron Valence second-order Perturbation Theory (NEVPT2) is far more stable, requiring a similar number of walkers to converge the sc-NEVPT2 expectation values as to converge the underlying CI problem. We demonstrate the application of the stochastic approach to the computation of sc-NEVPT2 within a (24, 24) active space in a biologically relevant system and show that small numbers of walkers are sufficient for a faithful sampling of the sc-NEVPT2 energy to chemical accuracy, despite the active space already exceeding the limits of practicality for traditional approaches. This raises prospects of an efficient stochastic solver for multireference chemical problems requiring large active spaces, with an accurate treatment of external orbitals.
The computationally expensive evaluation and storage of high-rank reduced density matrices (RDMs) has been the bottleneck in the calculation of dynamic correlation for multireference wave functions in large active spaces. We present a stochastic formulation of multireference configuration interaction and perturbation theory that avoids the need for these expensive RDMs. The algorithm presented here is flexible enough to incorporate a wide variety of active space reference wave functions, including selected configuration interaction, matrix product states, and symmetry-projected Jastrow mean field wave functions. It enjoys the usual attractive features of Monte Carlo methods, such as embarrassing parallelizability and low memory costs. We find that the stochastic algorithm is already competitive with the deterministic algorithm for small active spaces, containing as few as 14 orbitals. We illustrate the utility of our stochastic formulation using benchmark applications.
In this article we report a stochastic evaluation of the recently proposed multireference linearized coupled cluster theory [S. Sharma and A. Alavi, J. Chem. Phys. 143, 102815 (2015)]. In this method, both the zeroth-order and first-order wavefunctions are sampled stochastically by propagating simultaneously two populations of signed walkers. The sampling of the zeroth-order wavefunction follows a set of stochastic processes identical to the one used in the full configuration interaction quantum Monte Carlo (FCIQMC) method. To sample the first-order wavefunction, the usual FCIQMC algorithm is augmented with a source term that spawns walkers in the sampled first-order wavefunction from the zeroth-order wavefunction. The second-order energy is also computed stochastically but requires no additional overhead outside of the added cost of sampling the first-order wavefunction. This fully stochastic method opens up the possibility of simultaneously treating large active spaces to account for static correlation and recovering the dynamical correlation using perturbation theory. The method is used to study a few benchmark systems including the carbon dimer and aromatic molecules. We have computed the singlet-triplet gaps of benzene and m-xylylene. For m-xylylene, which has proved difficult for standard complete active space self consistent field theory with perturbative correction, we find the singlet-triplet gap to be in good agreement with the experimental values.