Researcher profile

Eric Neuscamman

Eric Neuscamman contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

Applying generalized variational principles to excited-state-specific complete active space self-consistent field theory

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.

preprint2020arXiv

A generalized variational principle with applications to excited state mean field theory

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.

preprint2020arXiv

A hybrid approach to excited-state-specific variational Monte Carlo and doubly excited states

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.

preprint2020arXiv

An N$^5$-scaling excited-state-specific perturbation theory

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ø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.

preprint2020arXiv

Excited State Mean-Field Theory without Automatic Differentiation

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.

preprint2020arXiv

Improving excited state potential energy surfaces via optimal orbital shapes

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.

preprint2020arXiv

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion Quantum Monte Carlo

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.

preprint2019arXiv

A Density Functional Extension to Excited State Mean-Field Theory

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.

preprint2013arXiv

A Jastrow factor coupled cluster theory for weak and strong electron correlation

We present a Jastrow-factor-inspired variant of coupled cluster theory that accurately describes both weak and strong electron correlation. Compatibility with quantum Monte Carlo allows for variational energy evaluations and an antisymmetric geminal power reference, two features not present in traditional coupled cluster that facilitate a nearly exact description of the strong electron correlations in minimal-basis N2 bond breaking. In double-zeta treatments of the HF and H2O bond dissociations, where both weak and strong correlations are important, this polynomial cost method proves more accurate than either traditional coupled cluster or complete active space perturbation theory. These preliminary successes suggest a deep connection between the ways in which cluster operators and Jastrow factors encode correlation.

preprint2013arXiv

The Jastrow antisymmetric geminal power in Hilbert space: theory, benchmarking, and application to a novel transition state

The Jastrow-modified antisymmetric geminal power (JAGP) ansatz in Hilbert space successfully overcomes two key failings of other pairing theories, namely a lack of inter-pair correlations and a lack of multiple resonance structures, while maintaining a polynomially scaling cost, variational energies, and size consistency. Here we present efficient quantum Monte Carlo algorithms that evaluate and optimize the JAGP energy for a cost that scales as the fifth power of the system size. We demonstrate the JAGP's ability to describe both static and dynamic correlation by applying it to bond stretching in H2O, C2, and N2 as well as to a novel, multi-reference transition state of ethene. JAGP's accuracy in these systems outperforms even the most sophisticated single-reference methods and approaches that of exponentially-scaling active space methods.

preprint2012arXiv

A correlator product state study of molecular magnetism in the giant Keplerate Mo72Fe30

We have studied the properties of the giant Keplerate molecular magnet Mo72Fe30, as a function of applied magnetic field, using the correlator product state (CPS) tensor network ansatz. The magnet is modeled with an S = 5/2 antiferromagnetic Heisenberg Hamiltonian on the 30-site icosidodecahedron lattice, a model for which exact diagonalization is infeasible. The CPS ansatz produces significant improvements in variational energies relative to previous studies using the density matrix renormalization group, a result of its superior ability to handle strong correlation in two dimensional spin systems. The CPS results reaffirm that the ground state energies adhere qualitatively to the parabolic progression of the rotational band model (RBM), but show important deviations near 1/3 of the saturation field. These deviations predict anomalous behavior in the differential magnetization and heat capacity that cannot be explained by the RBM alone. Finally, we show that these energetic deviations originate from a qualitative change in the ground state that resembles a finite size analogue of a phase transition.

preprint2012arXiv

An exactly size consistent geminal power via Jastrow factor networks in a local one particle basis

The accurate but expensive product of geminals ansatz may be approximated by a geminal power, but this approach sacrifices size consistency. Here we show both analytically and numerically that a size consistent form very similar to the product of geminals can be recovered using a network of location specific Jastrow factors. Upon variational energy minimization, the network creates particle number projections that remove the charge fluctuations responsible for size inconsistency. This polynomial cost approach captures strong many-electron correlations, giving a maximum error of just 1.8 kcal/mol during the double-bond dissociation of H2O in an STO-3G atomic orbital basis.

preprint2012arXiv

Striped spin liquid crystal ground state instability of kagome antiferromagnets

The Dirac spin liquid ground state of the spin 1/2 Heisenberg kagome antiferromagnet has potential instabilities[1-4]. This has been suggested as the reason why it is not strongly supported in large-scale numerical calculations[5]. However, previous attempts to observe these instabilities have failed. We report on the discovery of a projected BCS state with lower energy than the projected Dirac spin liquid state which provides new insight into the stability of the ground state of the kagome antiferromagnet. The new state has three remarkable features. First, it breaks both spatial symmetry in an unusual way that may leave spinons deconfined along one direction. Second, it breaks the U(1) gauge symmetry down to $Z_2$. Third, it has the spatial symmetry of a previously proposed "monopole" suggesting that it is an instability of the Dirac spin liquid. The state described herein also shares a remarkable similarity to the distortion of the kagome lattice observed at low Zn concentrations in Zn-Paratacamite suggesting it may already be realized in these materials.

preprint2011arXiv

Optimizing large parameter sets in variational quantum Monte Carlo

We present a technique for optimizing hundreds of thousands of variational parameters in variational quantum Monte Carlo. By introducing iterative Krylov subspace solvers and by multiplying by the Hamiltonian and overlap matrices as they are sampled, we remove the need to construct and store these matrices and thus bypass the most expensive steps of the stochastic reconfiguration and linear method optimization techniques. We demonstrate the effectiveness of this approach by using stochastic reconfiguration to optimize a correlator product state wavefunction with a pfaffian reference for four example systems. In two examples on the two dimensional Hubbard model, we study 16 and 64 site lattices, recovering energies accurate to 1% in the smaller lattice and predicting particle-hole phase separation in the larger. In two examples involving an ab initio Hamiltonian, we investigate the potential energy curve of a symmetrically dissociated 4x4 hydrogen lattice as well as the singlet-triplet gap in free base porphin. In the hydrogen system we recover 98% or more of the correlation energy at all geometries, while for porphin we compute the gap in a 24 orbital active space to within 0.02eV of the exact result. The numbers of variational parameters in these examples range from 4x10^3 to 5x10^5, demonstrating an ability to go far beyond the reach of previous formulations of stochastic reconfiguration.

preprint2010arXiv

Further developments in correlator product states: deterministic optimization and energy evaluation

Correlator product states (CPS) are a class of tensor network wavefunctions applicable to strongly correlated problems in arbitrary dimensions. Here, we present a method for optimizing and evaluating the energy of the CPS wavefunction that is non-variational but entirely deterministic. The fundamental assumption underlying our technique is that the CPS wavefunction is an exact eigenstate of the Hamiltonian, allowing the energy to be obtained approximately through a projection of the Schrödinger equation. The validity of this approximation is tested on two dimensional lattices for the spin-1/2 antiferromagnetic Heisenberg model, the spinless Hubbard model, and the full Hubbard model. In each of these models, the projected method reproduces the variational CPS energy to within 1%. For fermionic systems, we also demonstrate the incorporation of a Slater determinant reference into the ansatz, which allows CPS to act as a generalization of the Jastrow-Slater wavefunction.