Search

Guided by perturbative analysis, we improve the accuracy of Aufbau suppressed coupled cluster theory in simple single excitations, multi-configurational single excitations, and charge transfer excitations while keeping the cost of its leading-order terms precisely in line with ground state coupled cluster. Combining these accuracy improvements with a more efficient implementation based on spin-adaptation, we observe high accuracy in a large test set of single excitations, and, in particular, a mean unsigned error for charge transfer states that outperforms equation-of-motion coupled cluster theory by 0.25 eV. We discuss how these results are achieved via a systematic identification of which amplitudes to prioritize for single- and multi-configurational excited states, and how this prioritization differs in important ways from the ground state theory. In particular, our data show that a partial linearization of the theory increases accuracy by mitigating unwanted side effects of Aufbau suppression., Comment: 22 pages, 5 figures, 3 tables

We present an approach for augmenting Gaussian atomic orbitals with correct nuclear cusps. Like the atomic orbital basis set itself, and unlike previous cusp corrections, this approach is independent of the many-body method used to prepare wave functions for quantum Monte Carlo. Once the basis set and molecular geometry are specified, the cusp-corrected atomic orbitals are uniquely specified, regardless of which density functionals, quantum chemistry methods, or subsequent variational Monte Carlo optimizations are employed. We analyze the statistical improvement offered by these cusps in a number of molecules and find them to offer similar advantages as molecular-orbital-based approaches while maintaining independence from the choice of many-body method., Comment: 11 pages, 5 figures, 2 tables, CGAOWS repository: https://github.com/eneuscamman/cgaows

We introduce an approach to improve single-reference coupled cluster theory in settings where the Aufbau determinant is absent from or plays only a small role in the true wave function. Using a de-excitation operator that can be efficiently hidden within a similarity transform, we create a coupled cluster wave function in which de-excitations work to suppress the Aufbau determinant and produce wave functions dominated by other determinants. Thanks to an invertible and fully exponential form, the approach is systematically improvable, size consistent, size extensive, and, interestingly, size intensive in a granular way that should make the adoption of some ground state techniques such as local correlation relatively straightforward. In this initial study, we apply the general formalism to create a state-specific method for orbital-relaxed singly excited states. We find that this approach matches the accuracy of similar-cost equation-of-motion methods in valence excitations while offering improved accuracy for charge transfer states. We also find the approach to be more accurate than excited-state-specific perturbation theory in both types of states., Comment: 17 pages, 4 tables, 1 figure

The photo-ferrioxalate system (PFS), [Fe(III)(C$_2$O$_4$)]$^{3-}$, more than an exact chemical actinometer, has been extensively applied in wastewater and environment treatment. Despite many experimental efforts to improve clarity, important aspects of the mechanism of ferrioxalate photolysis are still under debate. In this paper, we employ the recently developed W$\Gamma$-CASSCF to investigate the ligand-to-metal charge-transfer states key to the ferrioxalate photolysis. This investigation provides a qualitative picture of these states and key potential energy surface features related to the photolysis. Our theoretical results are consistent with the prompt charge transfer picture seen in recent experiments and clarify some features that are not visible in experiments. Two ligand-to-metal charge-transfer states contribute to the photolysis of ferrioxalate, and the avoided crossing barrier between them is low compared to the initial photoexcitation energy. Our data also clarify that one Fe-O bond cleaves first, followed by the C-C bond and the other Fe-O bond., Comment: 20 pages, 7 figures, 3 tables

We theoretically investigate the influence of diradical electron spin coupling on the time-resolved X-ray absorption spectra of the photochemical ring opening of furanone. We predict geometry dependent carbon K-edge signals involving transitions from core orbitals to both singly and unoccupied molecular orbitals. The most obvious features of the ring opening come from the carbon atom directly involved in the bond breaking, through its transition to both the newly formed SOMO and the available LUMO state. In addition to this primary feature, the singlet spin coupling of four unpaired electrons that arises in the core-to-LUMO states creates additional geometry dependence in some spectral features, with both oscillator strengths and relative excitation energies varying observably as a function of the ring opening. We attribute this behavior to a spin-occupancy-induced selection rule, which occurs when singlet spin coupling is enforced in the diradical state. Notably, one of these geometry-sensitive core-to-LUMO transitions excites core electrons from a backbone carbon not involved in the bond breaking, providing a novel non-local X-ray probe of chemical dynamics arising from electron spin coupling., Comment: 52 Pages, 13 Figures

We explore the performance of a recently-introduced $N^5$-scaling excited-state-specific second order perturbation theory (ESMP2) on the singlet excitations of the Thiel benchmarking set. We find that, without regularization, ESMP2 is quite sensitive to $\pi$ system size, performing well in molecules with small $\pi$ systems but poorly in those with larger $\pi$ systems. With regularization, ESMP2 is far less sensitive to $\pi$ system size and shows a higher overall accuracy on the Thiel set than CC2, EOM-CCSD, CC3, and a wide variety of time-dependent density functional approaches. Unsurprisingly, even regularized ESMP2 is less accurate than multi-reference perturbation theory on this test set, which can in part be explained by the set's inclusion of some doubly excited states but none of the strong charge transfer states that often pose challenges for state-averaging. Beyond energetics, we find that the ESMP2 doubles norm offers a relatively low-cost way to test for doubly excited character without the need to define an active space.

We present an excited-state-specific coupled-cluster approach in which both the molecular orbitals and cluster amplitudes are optimized for an individual excited state. The theory is formulated via a pseudoprojection of the traditional coupled-cluster wavefunction that allows correlation effects to be introduced atop an excited state mean field starting point. The approach shares much in common with ground state CCSD, including size extensivity and an $N^6$ cost scaling. Preliminary numerical tests show that, when augmented with $N^5$-cost perturbative corrections for key terms, the method can improve over excited-state-specific second order perturbation theory in valence, charge transfer, and Rydberg states., Comment: 17 pages, 2 figures, 9 tables Updated to reflect changes made during peer review process, conclusions remain the same

Modern quantum Monte Carlo (QMC) methods often capture electron correlation through both explicitly correlating Jastrow factors and small to mid-sized configuration interaction (CI) expansions. Here, we study the additional optimization difficulty created by including increasing numbers of CI parameters. We find evidence that the quality of Variational Mone Carlo (VMC) optimization can be limited by the ability to statistically resolve the CI parameters in the presence of a Jastrow factor. Although using larger statistical samples can mitigate this issue and bring an optimization closer to its true minimum, this approach is computationally intensive. We present evidence that similar gains to optimization quality can be had without increasing the sample size by avoiding CI parameters that are statistically the most difficult to resolve. Our findings suggest that, in addition to the value of using traditional selected configuration interaction (sCI) methods to prepare VMC wave functions, sCI-like methods can play an important role within VMC by improving the effectiveness of stochastic energy minimization., Comment: 15 pages, 6 figures, 3 tables

We investigate the issue of optimization stability in variance-based state-specific variational Monte Carlo, discussing the roles of the objective function, the complexity of wave function ansatz, the amount of sampling effort, and the choice of minimization algorithm. Using a small cyanine dye molecule as a test case, we systematically perform minimizations using variants of the linear method as both a standalone algorithm and in a hybrid combination with accelerated descent. We demonstrate that adaptive step control is crucial for maintaining the linear method's stability when optimizing complicated wave functions and that the hybrid method enjoys both greater stability and minimization performance. As a verification of variance minimization's practical utility, we report an excitation energy in the cyanine dye that is in good agreement with both benchmark quantum chemistry values and results obtained from state-averaged energy-based variational Monte Carlo., Comment: 30 pages, 15 figures, 2 tables

We present a discussion of recent progress in excited-state-specific quantum chemistry and quantum Monte Carlo alongside a demonstration of how a combination of methods from these two fields can offer reliably accurate excited state predictions across singly excited, doubly excited, and charge transfer states. Both of these fields have seen important advances supporting excited state simulation in recent years, including the introduction of more effective excited-state-specific optimization methods, improved handling of complicated wave function forms, and ways of explicitly balancing the quality of wave functions for ground and excited states. To emphasize the promise that exists at this intersection, we provide demonstrations using a combination of excited-state-specific complete active space self-consistent field theory, selected configuration interaction, and state-specific variance minimization. These demonstrations show that combining excited-state-specific quantum chemistry and variational Monte Carlo can be more reliably accurate than either equation of motion coupled cluster theory or multi-reference perturbation theory, and that it can offer new clarity in cases where existing high-level methods do not agree., Comment: 21 pages, 7 figures, 16 tables

We employ a generalized variational principle to improve the stability, reliability, and precision of fully excited-state-specific complete active space self-consistent field theory. Compared to previous approaches that similarly seek to tailor this ansatz's orbitals and configuration interaction expansion for an individual excited state, we find the present approach to be more resistant to root flipping and better at achieving tight convergence to an energy stationary point. Unlike state-averaging, this approach allows orbital shapes to be optimal for individual excited states, which is especially important for charge transfer states and some doubly excited states. We demonstrate the convergence and state-targeting abilities of this method in LiH, ozone, and MgO, showing in the latter that it is capable of finding three excited state energy stationary points that no previous method has been able to locate., Comment: 17 pages, 7 figures, 6 tables

We extend our hybrid linear-method/accelerated-descent variational Monte Carlo optimization approach to excited states and investigate its efficacy in double excitations. In addition to showing a superior statistical efficiency when compared to the linear method, our tests on the carbon dimer and cyclopentadiene show good energetic agreement with benchmark methods and experiment, respectively. We also demonstrate the ability to treat double excitations in systems that are too large for a full treatment by selective configuration interaction methods via an application to 4-aminobenzonitrile. Finally, we investigate the stability of state-specific variance optimization against collapse to other states' variance minima and find that symmetry, ansatz quality, and sample size all have roles to play in achieving stability., Comment: 44 pages, 9 figures, 7 tables plus supplementary material, accepted by The Journal of Chemical Physics

We test the efficacy of excited state mean field theory and its excited-state-specific perturbation theory on the prediction of K-edge positions and X-ray peak separations. We find that the mean field theory is surprisingly accurate, even though it contains no accounting of differential electron correlation effects. In the perturbation theory, we test multiple core-valence separation schemes and find that, with the mean field theory already so accurate, electron-counting biases in one popular separation scheme become a dominant error when predicting K-edges. Happily, these appear to be relatively easy to correct for, leading to a perturbation theory for K-edge positions that is lower scaling and more accurate than coupled cluster theory and competitive in accuracy with recent high-accuracy results from restricted open-shell Kohn Sham theory. For peak separations, our preliminary data show excited state mean field theory to be exceptionally accurate, but more extensive testing will be needed to see how it and its perturbation theory compare to coupled cluster peak separations more broadly., Comment: 7 pages, 2 figures, 3 tables

We present a systematically-improvable approach to core excitations in variational Monte Carlo. Building on recent work in excited-state-specific Monte Carlo, we show how a straightforward protocol, starting from a quantum chemistry guess, is able to capture core state's strong orbital relaxations, maintain accuracy in the near-nuclear region during these relaxations, and explicitly balance accuracy between ground and core excited states. In water, ammonia, and methane, which serve as prototypical representatives for oxygen, nitrogen, and carbon core states, respectively, this approach delivers accuracies on par with the best available theoretical methods even when using relatively small wave function expansions., Comment: 10 pages, 4 figures, 1 table

We demonstrate that, rather than resorting to high-cost dynamic correlation methods, qualitative failures in excited-state potential energy surface predictions can often be remedied at no additional cost by ensuring that optimal molecular orbitals are used for each individual excited state. This approach also avoids the weighting choices required by state-averaging and dynamic weighting and obviates their need for expensive wave function response calculations when relaxing excited state geometries. Although multi-state approaches are of course preferred near conical intersections, other features of excited-state potential energy surfaces can benefit significantly from our single state approach. In three different systems, including a double bond dissociation, a biologically relevant amino hydrogen dissociation, and an amino-to-ring intramolecular charge transfer, we show that state-specific orbitals offer qualitative improvements over the state-averaged status quo., Comment: 6 pages, 6 figures, 1 table

We show that, as in Hartree Fock theory, the orbitals for excited state mean field theory can be optimized via a self-consistent one-electron equation in which electron-electron repulsion is accounted for through mean field operators. In addition to showing that this excited state ansatz is sufficiently close to a mean field product state to admit a one-electron formulation, this approach brings the orbital optimization speed to within roughly a factor of two of ground state mean field theory. The approach parallels Hartree Fock theory in multiple ways, including the presence of a commutator condition, a one-electron mean-field working equation, and acceleration via direct inversion in the iterative subspace. When combined with a configuration interaction singles Davidson solver for the excitation coefficients, the self consistent field formulation dramatically reduces the cost of the theory compared to previous approaches based on quasi-Newton descent., Comment: 6 pages, 3 tables, 1 figure, plus supplementary material

We show that by working in a basis similar to that of the natural transition orbitals and using a modified zeroth order Hamiltonian, the cost of a recently-introduced perturbative correction to excited state mean field theory can be reduced from seventh to fifth order in the system size. The (occupied)$^2$(virtual)$^3$ asymptotic scaling matches that of ground state second order M{\o}ller-Plesset theory, but with a significantly higher prefactor because the bottleneck is iterative: it appears in the Krylov-subspace-based solution of the linear equation that yields the first order wave function. Here we discuss the details of the modified zeroth order Hamiltonian we use to reduce the cost as well as the automatic code generation process we used to derive and verify the cost scaling of the different terms. Overall, we find that our modifications have little impact on the method's accuracy, which remains competitive with singles and doubles equation-of-motion coupled cluster., Comment: 11 pages, 2 figures, 7 tables, accepted by the Journal of Chemical Theory and Computation

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing accuracy. Advances in real space methods include techniques for accurate computation of band gaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods including GW and density functional based techniques. To provide an improved foundation for these calculations we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

We present a formulation of excited state mean-field theory in which the derivatives with respect to the wave function parameters needed for wave function optimization (not to be confused with nuclear derivatives) are expressed analytically in terms of a collection of Fock-like matrices. By avoiding the use of automatic differentiation and grouping Fock builds together, we find that the number of times we must access the memory-intensive two-electron integrals can be greatly reduced. Furthermore, the new formulation allows the theory to exploit existing strategies for efficient Fock matrix construction. We demonstrate this advantage explicitly via the shell-pair screening strategy, with which we achieve a cubic overall cost scaling. Using this more efficient implementation, we also examine the theory's ability to predict charge redistribution during charge transfer excitations. Using coupled cluster as a benchmark, we find that by capturing orbital relaxation effects and avoiding self-interaction errors, excited state mean field theory out-performs other low-cost methods when predicting the charge density changes of charge transfer excitations.

We present a generalization of the variational principle that is compatible with any Hamiltonian eigenstate that can be specified uniquely by a list of properties. This variational principle appears to be compatible with a wide range of electronic structure methods, including mean-field theory, density functional theory, multi-reference theory, and quantum Monte Carlo. Like the standard variational principle, this generalized variational principle amounts to the optimization of a nonlinear function that, in the limit of an arbitrarily flexible wave function, has the desired Hamiltonian eigenstate as its global minimum. Unlike the standard variational principle, it can target excited states and select individual states in cases of degeneracy or near-degeneracy. As an initial demonstration of how this approach can be useful in practice, we employ it to improve the optimization efficiency of excited state mean field theory by an order of magnitude. With this improved optimization, we are able to demonstrate that the accuracy of the corresponding second-order perturbation theory rivals that of singles-and-doubles equation-of-motion coupled cluster in a substantially broader set of molecules than could be explored by our previous optimization methodology.

We investigate an extension of excited state mean-field theory in which the energy expression is augmented with density functional components in an effort to include the effects of weak electron correlations. The approach remains variational and entirely time-independent, allowing it to avoid some of the difficulties associated with linear response and the adiabatic approximation. In particular, all of the electrons' orbitals are relaxed state specifically and there is no reliance on Kohn-Sham orbital energy differences, both of which are important features in the context of charge transfer. Preliminary testing shows clear advantages for single-component charge transfer states, but the method, at least in its current form, is less reliable for states in which multiple particle-hole transitions contribute significantly.

We present a comparison between a number of recently introduced low-memory wave function optimization methods for variational Monte Carlo in which we find that first and second derivative methods possess strongly complementary relative advantages. While we find that low-memory variants of the linear method are vastly more efficient at bringing wave functions with disparate types of nonlinear parameters to the vicinity of the energy minimum, accelerated descent approaches are then able to locate the precise minimum with less bias and lower statistical uncertainty. By constructing a simple hybrid approach that combines these methodologies, we show that all of these advantages can be had at once when simultaneously optimizing large determinant expansions, molecular orbital shapes, traditional Jastrow correlation factors, and more nonlinear many-electron Jastrow factors.

We present a method for finding individual excited states' energy stationary points in complete active space self-consistent field theory that is compatible with standard optimization methods and highly effective at overcoming difficulties due to root flipping and near-degeneracies. Inspired by both the maximum overlap method and recent progress in excited state variational principles, our approach combines these ideas in order to track individual excited states throughout the orbital optimization process. In a series of tests involving root flipping, near-degeneracies, charge transfers, and double excitations, we show that this approach is more effective for state-specific optimization than either the naive selection of roots based on energy ordering or a more direct generalization of the maximum overlap method. Furthermore, we provide evidence that this state-specific approach improves the performance of complete active space perturbation theory. With a simple implementation, a low cost, and compatibility with large active space methods, the approach is designed to be useful in a wide range of excited state investigations., Comment: 13 pages, submitted to JCTC

We introduce a basis of counting functions that, by cleanly tessellating three dimensional space, allows real space number counting Jastrow factors to be straightforwardly applied to general molecular situations. By exerting direct control over electron populations in local regions of space and encoding pairwise correlations between these populations, these Jastrow factors allow even very simple reference wave functions to adopt nodal surfaces well suited to many strongly correlated settings. Being trivially compatible with traditional Jastrow factors and diffusion Monte Carlo and having the same cubic per-sample cost scaling as a single determinant trial function, these Jastrow factors thus offer a powerful new route to the simultaneous capture of weak and strong electron correlation effects in a wide variety of molecular and materials settings. In multiple strongly correlated molecular examples, we show that even when paired with the simplest possible single determinant reference, these Jastrow factors allow quantum Monte Carlo to out-perform coupled cluster theory and approach the accuracy of traditional multi-reference methods., Comment: 21 pages, 22 figures, 7 tables

We combine recent advances in excited state variational principles, fast multi-Slater Jastrow methods, and selective configuration interaction to create multi-Slater Jastrow wave function approximations that are optimized for individual excited states. In addition to the Jastrow variables and linear expansion coefficients, this optimization includes state-specific orbital relaxations in order to avoid the compromises necessary in state-averaged approaches. We demonstrate that, when combined with variance matching to help balance the quality of the approximation across different states, this approach delivers accurate excitation energies even when using very modest multi-Slater expansions. Intriguingly, this accuracy is maintained even when studying a difficult chlorine-anion-to-$\pi^{*}$ charge transfer in which traditional state-averaged multi-reference methods must contend with different states that require drastically different orbital relaxations., Comment: 16 pages, 6 figures, 2 tables

We perform excited-state variational Monte Carlo and diffusion Monte Carlo calculations using a simple and efficient wave function ansatz. This ansatz follows the recent variation-after-response formalism, accurately approximating a configuration interaction singles wave function as a sum of only two non-orthogonal Slater determinants, and further including important orbital relaxation. The ansatz is used to perform diffusion Monte Carlo calculations with large augmented basis sets, comparing to benchmarks from near-exact quantum chemical methods. The significance of orbital optimization in excited-state diffusion Monte Carlo is demonstrated, and the excited-state optimization procedure is discussed in detail. Diffuse excited states in water and formaldehyde are studied, in addition to the formaldimine and benzonitrile molecules., Comment: 12 pages, 3 figures

We show that for both single-Slater-Jastrow and Jastrow geminal power wave functions, the formal cost scaling of Hilbert space variational Monte Carlo can be reduced from fifth to fourth order in the system size, thus bringing it in line with the long-standing scaling of its real space counterpart. While traditional quantum chemistry methods can reduce costs related to the two-electron integral tensor through resolution of the identity and Cholesky decomposition approaches, we show that such approaches are ineffective in the presence of Hilbert space Jastrow factors. Instead, we develop a simple semi-stochastic approach that can take similar advantage of the near-sparsity of this four-index tensor. Through demonstrations on alkanes of increasing length, we show that accuracy and overall statistical uncertainty are not meaningfully affected and that a total cost crossover is reached as early as 50 electrons., Comment: 8 pages, 7 figures

We report variational and fixed-node diffusion quantum Monte Carlo (QMC) calculations of anti-ferromagnetic iron oxide (FeO) in the ground state B1 crystal structure. The goal of this study was a systematic investigation of the sensitivity of several ground state properties to a variety of QMC wave function generation techniques including advanced wave functions such as multi-determinant expansions and backflow transformations. We found that the predicted lattice distortion was largely controlled by the choice of single particle orbitals used to construct the wave function, rather than by subsequent wave function optimization techniques within QMC. However, the absolute magnetic moment was remarkably insensitive to the method of wave function construction. QMC estimates of total spin density indicate that in addition to strong electronic correlation of the Fe $3d$ states, charge transfer may be an important but challenging piece of physics to accurately capture within existing QMC methods. Finally, we highlight the need for advanced and systematically improvable many-body wave functions suitable for accurately describing challenging real systems., Comment: Further calculations are in progress. This version does not accurately represent the results of our investigation

We present a mean field theory for excited states that is broadly analogous to ground state Hartree-Fock theory. Like Hartree-Fock, our approach is deterministic, state-specific, applies a variational principle to a minimally correlated ansatz, produces energy stationary points, relaxes the orbital basis, has a Fock-build cost-scaling, and can serve as the foundation for correlation methods such as perturbation theory and coupled cluster theory. To emphasize this last point, we pair our mean field approach with an excited state analogue of second order Moller-Plesset theory and demonstrate that in water, formaldehyde, neon, and stretched lithium fluoride, the resulting accuracy far exceeds that of configuration interaction singles and rivals that of equation of motion coupled cluster., Comment: 6 pages

We present an approach to studying optical band gaps in real solids in which quantum Monte Carlo methods allow for the application of a rigorous variational principle to both ground and excited state wave functions. In tests that include small, medium, and large band gap materials, optical gaps are predicted with a mean-absolute-deviation of 3.5% against experiment, less than half the equivalent errors for typical many-body perturbation theories. The approach is designed to be insensitive to the choice of density functional, a property we exploit in order to provide insight into how far different functionals are from satisfying the assumptions of many body perturbation theory. We explore this question most deeply in the challenging case of ZnO, where we show that although many commonly used functionals have shortcomings, there does exist a one particle basis in which perturbation theory's zeroth order picture is sound. Insights of this nature should be useful in guiding the future application and improvement of these widely used techniques., Comment: 8 pages, 5 figures, 2 tables

QMCPACK is an open source quantum Monte Carlo package for ab-initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wave functions in conjunction with a sophisticated optimizer capable of optimizing tens of thousands of parameters. The orbital space auxiliary field quantum Monte Carlo method is also implemented, enabling cross validation between different highly accurate methods. The code is specifically optimized for calculations with large numbers of electrons on the latest high performance computing architectures, including multicore central processing unit (CPU) and graphical processing unit (GPU) systems. We detail the program's capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://www.qmcpack.org .

We demonstrate that a broad class of excited state variational principles is not size consistent. In light of this difficulty, we develop and test an approach to excited state optimization that transforms between variational principles in order to achieve state selectivity, size consistency, and compatibility with quantum Monte Carlo. To complement our formal analysis, we provide numerical examples that confirm these properties and demonstrate how they contribute to a more black box approach to excited states in quantum Monte Carlo., Comment: 27 pages, 10 figures, 2 tables

We present a simple and efficient wave function ansatz for the treatment of excited charge-transfer states in real-space quantum Monte Carlo methods. Using the recently-introduced variation-after-response method [J. Chem. Phys. 145, 081103 (2016)], this ansatz allows a crucial orbital optimization step to be performed beyond a configuration interaction singles expansion, while only requiring calculation of two Slater determinant objects. We demonstrate this ansatz for the illustrative example of the stretched LiF molecule, for a range of excited states of formaldehyde, and finally for the more challenging ethylene-tetrafluoroethylene molecule., Comment: 10 pages, 5 figures

In the regime where traditional approaches to electronic structure cannot afford to achieve accurate energy differences via exhaustive wave function flexibility, rigorous approaches to balancing different states' accuracies become desirable. As a direct measure of a wave function's accuracy, the energy variance offers one route to achieving such a balance. Here, we develop and test a variance matching approach for predicting excitation energies within the context of variational Monte Carlo and selective configuration interaction. In a series of tests on small but difficult molecules, we demonstrate that the approach it is effective at delivering accurate excitation energies when the wave function is far from the exhaustive flexibility limit. Results in C$_3$, where we combine this approach with variational Monte Carlo orbital optimization, are especially encouraging., Comment: 11 pages, 7 figures, 2 tables