State-averaged Monte Carlo configuration interaction applied to electronically excited states

Using a multistate mapping approach to surface hopping to predict the ultrafast electron diffraction signal of gas-phase cyclobutanone

Using the recently developed multistate mapping approach to surface hopping (multistate MASH) method combined with SA(3)-CASSCF(12,12)/aug-cc-pVDZ electronic structure calculations, the gas-phase isotropic ultrafast electron diffraction (UED) of cyclobutanone is predicted and analyzed. After excitation into the n-3s Rydberg state (S2), cyclobutanone can relax through two S2/S1 conical intersections, one characterized by compression of the CO bond and the other by dissociation of the α-CC bond. Subsequent transfer into the ground state (S0) is then achieved via two additional S1/S0 conical intersections that lead to three reaction pathways: α ring-opening, ethene/ketene production, and CO liberation. The isotropic gas-phase UED signal is predicted from the multistate MASH simulations, allowing for a direct comparison to the experimental data. This work, which is a contribution to the cyclobutanone prediction challenge, facilitates the identification of the main photoproducts in the UED signal and thereby emphasizes the importance of dynamics simulations for the interpretation of ultrafast experiments.

Modular Approach to Selected Configuration Interaction in an Arbitrary Spin Basis: Implementation and Comparison of Approaches

A modular selected configuration interaction (SCI) code has been developed that is based on the existing Monte-Carlo configuration interaction code (MCCI). The modularity allows various selection protocols to be implemented with ease and allows for fair comparison between wave functions built via different criteria. We have initially implemented adaptations of existing SCI theories, which are based on either energy- or coefficient-driven selection schemes. These codes have been implemented not only in the basis of Slater determinants (SDs) but also in the basis of configuration state functions (CSFs) and extended to state-averaged regimes. This allows one to take advantage of the reduced dimensionality of the wave function in the CSF basis and also the guarantee of pure spin states. All SCI methods were found to be able to predict potential energy surfaces to high accuracy, producing compact wave functions, when compared to full configuration interaction (FCI) for a variety of bond-breaking potential energy surfaces. The compactness of the error-controlled adaptive configuration interaction approach, particularly in the CSF basis, was apparent with nonparallelity errors within chemical accuracy while containing as little as 0.02% of the FCI CSF space. The size-to-accuracy was also extended to FCI spaces approaching one billion configurations.

Analytic Non-adiabatic Couplings for Selected Configuration Interaction via Approximate Degenerate Coupled Perturbed Hartree-Fock

We use degenerate perturbation theory and assume that for degenerate pairs of orbitals, the coupled perturbed Hartree-Fock coefficients are symmetric in the degenerate basis to show [Formula: see text] is the only modification needed in the original molecular orbital basis. This enables us to develop efficient and accurate analytic nonadiabatic couplings between electronic states for selected configuration interactions (CIs). Even when the states belong to different irreducible representations, degenerate orbital pairs cannot be excluded by symmetry. For various excited states of carbon monoxide and trigonal planar ammonia, we benchmark the method against the full CI and find it to be accurate. We create a semi-numerical approach and use it to show that the analytic approach is correct even when a high-symmetry structure is distorted to break symmetry so that near degeneracies in orbitals occur. For a range of geometries of trigonal planar ammonia, we find that the analytic non-adiabatic couplings for selected CI can achieve sufficient accuracy using a small fraction of the full CI space.

Efficient Computation of Two-Electron Reduced Density Matrices via Selected Configuration Interaction

We create an approach to efficiently calculate two-electron reduced density matrices (2-RDMs) using selected configuration interaction wavefunctions. This is demonstrated using the specific example of Monte Carlo configuration interaction (MCCI). The computation of the 2-RDMs is accelerated by using ideas from fast implementations of full configuration interaction (FCI) and recent advances in implementing the Slater–Condon rules using hardware bitwise operations. This method enables a comparison of MCCI and truncated CI 2-RDMs with FCI values for a range of molecules, which includes stretched bonds and excited states. The accuracy in energies, wavefunctions, and 2-RDMs is seen to exhibit a similar behavior. We find that MCCI can reach sufficient accuracy of the 2-RDM using significantly fewer configurations than truncated CI, particularly for systems with strong multireference character.

Electronic energies from coupled fermionic "Zombie" states' imaginary time evolution

Zombie states are a recently introduced formalism to describe coupled coherent fermionic states that address the fermionic sign problem in a computationally tractable manner. Previously, it has been shown that Zombie states with fractional occupations of spin orbitals obeyed the correct fermionic creation and annihilation algebra and presented results for real-time evolution [D. V. Shalashilin, J. Chem. Phys. 148, 194109 (2018)]. In this work, we extend and build on this formalism by developing efficient algorithms for evaluating the Hamiltonian and other operators between Zombie states and address their normalization. We also show how imaginary time propagation can be used to find the ground state of a system. We also present a biasing method, for setting up a basis set of random Zombie states, that allows much smaller basis sizes to be used while still accurately describing the electronic structure Hamiltonian and its ground state and describe a technique of wave function "cleaning" that removes the contributions of configurations with the wrong number of electrons, improving the accuracy further. We also show how low-lying excited states can be calculated efficiently using a Gram-Schmidt orthogonalization procedure. The proposed algorithm of imaginary time propagation on biased random grids of Zombie states may present an alternative to the existing quantum Monte Carlo methods.

Reinforcement Learning Configuration Interaction

Selected configuration interaction (sCI) methods exploit the sparsity of the full configuration interaction (FCI) wave function, yielding significant computational savings and wave function compression without sacrificing the accuracy. Despite recent advances in sCI methods, the selection of important determinants remains an open problem. We explore the possibility of utilizing reinforcement learning approaches to solve the sCI problem. By mapping the configuration interaction problem onto a sequential decision-making process, the agent learns on-the-fly which determinants to include and which to ignore, yielding a compressed wave function at near-FCI accuracy. This method, which we call reinforcement-learned configuration interaction, adds another weapon to the sCI arsenal and highlights how reinforcement learning approaches can potentially help solve challenging problems in electronic structure theory.

Investigation of challenging spin systems using Monte Carlo configuration interaction and the density matrix renormalization group

We investigate if a range of challenging spin systems can be described sufficiently well using Monte Carlo configuration interaction (MCCI) and the density matrix renormalization group (DMRG) in a way that heads toward a more "black box" approach. Experimental results and other computational methods are used for comparison. The gap between the lowest doublet and quartet state of methylidyne (CH) is first considered. We then look at a range of first-row transition metal monocarbonyls: MCO when M is titanium, vanadium, chromium, or manganese. For these MCO systems we also employ partially spin restricted open-shell coupled-cluster (RCCSD). We finally investigate the high-spin low-lying states of the iron dimer, its cation and its anion. The multireference character of these molecules is also considered. We find that these systems can be computationally challenging with close low-lying states and often multireference character. For this more straightforward application and for the basis sets considered, we generally find qualitative agreement between DMRG and MCCI. © 2017 Wiley Periodicals, Inc.

A Deterministic Projector Configuration Interaction Approach for the Ground State of Quantum Many-Body Systems

In this work we propose a novel approach to solve the Schrödinger equation which combines projection onto the ground state with a path-filtering truncation scheme. The resulting projector configuration interaction (PCI) approach realizes a deterministic version of the full configuration interaction quantum Monte Carlo (FCIQMC) method [Booth, G. H.; Thom, A. J. W.; Alavi, A. J. Chem. Phys. 2009, 131, 054106]. To improve upon the linearized imaginary-time propagator, we develop an optimal projector scheme based on an exponential Chebyshev expansion in the limit of an infinite imaginary time step. After writing the exact projector as a path integral in determinant space, we introduce a path filtering procedure that truncates the size of the determinantal basis and approximates the Hamiltonian. The path filtering procedure is controlled by one real threshold that determines the accuracy of the PCI energy and is not biased toward any determinant. Therefore, the PCI approach can equally well describe static and dynamic electron correlation effects. This point is illustrated in benchmark computations on N2 at both equilibrium and stretched geometries. In both cases, the PCI achieves chemical accuracy with wave functions that contain less than 0.5% determinants of full CI space. We also report computations on the ground state of C2 with up to quaduple-ζ basis sets and wave functions as large as 200 million determinants, which allow a direct comparison of the PCI, FCIQMC, and density matrix renormalization group (DMRG) methods. The size of the PCI wave function grows modestly with the number of unoccupied orbitals, and its accuracy may be tuned to match that of FCIQMC and DMRG.

Investigating Multireference Character and Correlation in Quantum Chemistry

We review a range of multireference diagnostics for quantum chemistry and discuss them in terms of choices of the molecular orbitals. We show how an approach1 of P.-O. Löwdin can also be viewed as quantifying the electron correlation via the spatial entanglement relative to a single determinant. We consider three example systems from quantum chemistry that exhibit three different combinations of multireference character and correlation: not strongly multireference and not strongly correlated, strongly multireference but not strongly correlated, and strongly multireference together with strong correlation. We find that a multireference measure (MR) does not change substantially with the cutoff used for a Monte Carlo configuration interaction calculation and investigate the effect of using natural orbitals. We see that a coupled-cluster singles and doubles diagnostic and a density-functional theory diagnostic give a correct general prediction of the multireference character for these systems. We also look at the issue of multireference character for a collection of noninteracting hydrogen molecules and the effect of basis size on the multireference character of a stretched hydrogen molecule.

Approaching exact hyperpolarizabilities via sum-over-states Monte Carlo configuration interaction

We propose using sum-over-states calculations with the compact wavefunctions of Monte Carlo configuration interaction to approach accurate values for higher-order dipole properties up to second hyperpolarizabilities in a controlled way. We apply the approach to small systems that can generally be compared with full configuration interaction (FCI) results. We consider hydrogen fluoride with a 6-31g basis and then look at results, including frequency dependent properties, in an aug-cc-pVDZ basis. We extend one calculation beyond FCI by using an aug-cc-pVTZ basis. The properties of an H4 molecule with multireference character are calculated in an aug-cc-pVDZ basis. We then investigate this method on a strongly multireference system with a larger FCI space by modelling the properties of carbon monoxide with a stretched geometry. The behavior of the approach with increasing basis size is considered by calculating results for the neon atom using aug-cc-pVDZ to aug-cc-pVQZ. We finally test if the unusual change in polarizability between the first two states of molecular oxygen can be reproduced by this method in a 6-31g basis.

Adaptive multiconfigurational wave functions

A method is suggested to build simple multiconfigurational wave functions specified uniquely by an energy cutoff Λ. These are constructed from a model space containing determinants with energy relative to that of the most stable determinant no greater than Λ. The resulting Λ-CI wave function is adaptive, being able to represent both single-reference and multireference electronic states. We also consider a more compact wave function parameterization (Λ+SD-CI), which is based on a small Λ-CI reference and adds a selection of all the singly and doubly excited determinants generated from it. We report two heuristic algorithms to build Λ-CI wave functions. The first is based on an approximate prescreening of the full configuration interaction space, while the second performs a breadth-first search coupled with pruning. The Λ-CI and Λ+SD-CI approaches are used to compute the dissociation curve of N2 and the potential energy curves for the first three singlet states of C2. Special attention is paid to the issue of energy discontinuities caused by changes in the size of the Λ-CI wave function along the potential energy curve. This problem is shown to be solvable by smoothing the matrix elements of the Hamiltonian. Our last example, involving the Cu2O2(2+) core, illustrates an alternative use of the Λ-CI method: as a tool to both estimate the multireference character of a wave function and to create a compact model space to be used in subsequent high-level multireference coupled cluster computations.