BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

C. J. Umrigar

C. J. Umrigar appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.chem-phphysics.comp-phcond-mat.str-elcond-mat.mtrl-scicond-mat.mes-hallcond-mat.otherquant-phcond-matphysics.atom-phatom-phchem-phphysics.atm-clus

Catalog footprint

What is connected

24works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

24 published item(s)

preprint2022arXiv

The chromium dimer: closing a chapter of quantum chemistry

The complex electronic structure and unusual potential energy curve of the chromium dimer have fascinated scientists for decades, with agreement between theory and experiment so far elusive. Here, we present a new ab initio simulation of the potential energy curve and vibrational spectrum that significantly improves on all earlier estimates. Our data support a shift in earlier experimental assignments of a cluster of vibrational frequencies by one quantum number. The new vibrational assignment yields an experimentally derived potential energy curve in quantitative agreement with theory across all bond lengths and across all measured frequencies. By solving this long-standing problem, our results raise the possibility of quantitative quantum chemical modeling of transition metal clusters with spectroscopic accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accurate many-body electronic structure near the basis set limit: application to the chromium dimer

We describe a method for computing near-exact energies for correlated systems with large Hilbert spaces. The method efficiently identifies the most important basis states (Slater determinants) and performs a variational calculation in the subspace spanned by these determinants. A semistochastic approach is then used to add a perturbative correction to the variational energy to compute the total energy. The size of the variational space is progressively increased until the total energy converges to within the desired tolerance. We demonstrate the power of the method by computing a near-exact potential energy curve (PEC) for a very challenging molecule -- the chromium dimer.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Chemistry of the spin-1/2 kagome Heisenberg antiferromagnet

We believe that a necessary first step in understanding the ground state properties of the spin-${\scriptstyle\frac{1}{2}}$ kagome Heisenberg antiferromagnet is a better understanding of this model's very large number of low energy singlet states. A description of the low energy states that is both accurate and amenable for numerical work may ultimately prove to have greater value than knowing only what these properties are, in particular when these turn on the delicate balance of many small energies. We demonstrate how this program would be implemented using the basis of spin-singlet dimerized states, though other bases that have been proposed may serve the same purpose. The quality of a basis is evaluated by its participation in all the low energy singlets, not just the ground state. From an experimental perspective, and again in light of the small energy scales involved, methods that can deliver all the low energy states promise more robust predictions than methods that only refine a fraction of these states.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Efficient heat-bath sampling in Fock space

We introduce an algorithm for sampling many-body quantum states in Fock space. The algorithm efficiently samples states with probability approximately proportional to an arbitrary function of the second-quantized Hamiltonian matrix element connecting the sampled state to the current state. We apply the new sampling algorithm to the recently-developed Semistochastic Full Configuration Interaction Quantum Monte Carlo method (S-FCIQMC), a semistochastic implementation of the power method for projecting out the ground state energy in a basis of Slater determinants. The heat-bath sampling requires modest additional computational time and memory compared to uniform sampling but results in newly-spawned weights that are approximately of the same magnitude, thereby greatly improving the efficiency of projection. A comparison in efficiency between uniform and approximate heat-bath sampling is performed on the all-electron nitrogen dimer at equilibrium in Dunning's cc-pVXZ basis sets with X in {D, T, Q, 5}, demonstrating a large gain in efficiency that increases with basis set size. In addition, a comparison in efficiency is performed on three all-electron first-row dimers, B_2, N_2, and F_2, in a cc-pVQZ basis, demonstrating that the gain in efficiency compared to uniform sampling also increases dramatically with the number of electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Introduction to the variational and diffusion Monte Carlo methods

We provide a pedagogical introduction to the two main variants of real-space quantum Monte Carlo methods for electronic-structure calculations: variational Monte Carlo (VMC) and diffusion Monte Carlo (DMC). Assuming no prior knowledge on the subject, we review in depth the Metropolis-Hastings algorithm used in VMC for sampling the square of an approximate wave function, discussing details important for applications to electronic systems. We also review in detail the more sophisticated DMC algorithm within the fixed-node approximation, introduced to avoid the infamous Fermionic sign problem, which allows one to sample a more accurate approximation to the ground-state wave function. Throughout this review, we discuss the statistical methods used for evaluating expectation values and statistical uncertainties. In particular, we show how to estimate nonlinear functions of expectation values and their statistical uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Comparison of polynomial approximations to speed up planewave-based quantum Monte Carlo calculations

The computational cost of quantum Monte Carlo (QMC) calculations of realistic periodic systems depends strongly on the method of storing and evaluating the many-particle wave function. Previous work [A. J. Williamson et al., Phys. Rev. Lett. 87, 246406 (2001); D. Alfè and M. J. Gillan, Phys. Rev. B 70, 161101 (2004)] has demonstrated the reduction of the O(N^3) cost of evaluating the Slater determinant with planewaves to O(N^2) using localized basis functions. We compare four polynomial approximations as basis functions -- interpolating Lagrange polynomials, interpolating piecewise-polynomial-form (pp-) splines, and basis-form (B-) splines (interpolating and smoothing). All these basis functions provide a similar speedup relative to the planewave basis. The pp-splines have eight times the memory requirement of the other methods. To test the accuracy of the basis functions, we apply them to the ground state structures of Si, Al, and MgO. The polynomial approximations differ in accuracy most strongly for MgO and smoothing B-splines most closely reproduce the planewave value for of the variational Monte Carlo energy. Using separate approximations for the Laplacian of the orbitals increases the accuracy sufficiently to justify the increased memory requirement, making smoothing B-splines, with separate approximation for the Laplacian, the preferred choice for approximating planewave-represented orbitals in QMC calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Energy Density Functionals From the Strong-Coupling Limit Applied to the Anions of the He Isoelectronic Series

Anions and radicals are important for many applications including environmental chemistry, semiconductors, and charge transfer, but are poorly described by the available approximate energy density functionals. Here we test an approximate exchange-correlation functional based on the exact strong-coupling limit of the Hohenberg-Kohn functional on the prototypical case of the He isoelectronic series with varying nuclear charge $Z<2$, which includes weakly bound negative ions and a quantum phase transition at a critical value of $Z$, representing a big challenge for density functional theory. We use accurate wavefunction calculations to validate our results, comparing energies and Kohn-Sham potentials, thus also providing useful reference data close to and at the quantum phase transition. We show that our functional is able to bind H$^-$ and to capture in general the physics of loosely bound anions, with a tendency to strongly overbind that can be proven mathematically. We also include corrections based on the uniform electron gas which improve the results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spectroscopic accuracy directly from quantum chemistry: application to ground and excited states of beryllium dimer

We combine explicit correlation via the canonical transcorrelation approach with the density matrix renormalization group and initiator full configuration interaction quantum Monte Carlo methods to compute a near-exact beryllium dimer curve, {\it without} the use of composite methods. In particular, our direct density matrix renormalization group calculations produce a well-depth of $D_e$=931.2 cm$^{-1}$ which agrees very well with recent experimentally derived estimates $D_e$=929.7$\pm 2$~cm$^{-1}$ [Science, 324, 1548 (2009)] and $D_e$=934.6~cm$^{-1}$ [Science, 326, 1382 (2009)]], as well the best composite theoretical estimates, $D_e$=938$\pm 15$~cm$^{-1}$ [J. Phys. Chem. A, 111, 12822 (2007)] and $D_e$=935.1$\pm 10$~cm$^{-1}$ [Phys. Chem. Chem. Phys., 13, 20311 (2011)]. Our results suggest possible inaccuracies in the functional form of the potential used at shorter bond lengths to fit the experimental data [Science, 324, 1548 (2009)]. With the density matrix renormalization group we also compute near-exact vertical excitation energies at the equilibrium geometry. These provide non-trivial benchmarks for quantum chemical methods for excited states, and illustrate the surprisingly large error that remains for 1$^1Σ^-_g$ state with approximate multi-reference configuration interaction and equation-of-motion coupled cluster methods. Overall, we demonstrate that explicitly correlated density matrix renormalization group and initiator full configuration interaction quantum Monte Carlo methods allow us to fully converge to the basis set and correlation limit of the non-relativistic Schrödinger equation in small molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Zigzag Phase Transition in Quantum Wires

We study the quantum phase transition of interacting electrons in quantum wires from a one-dimensional (1D) linear configuration to a quasi-1D zigzag arrangement using quantum Monte Carlo methods. As the density increases from its lowest values, first, the electrons form a linear Wigner crystal; then, the symmetry about the axis of the wire is broken as the electrons order in a quasi-1D zigzag phase; and, finally, the electrons form a disordered liquid-like phase. We show that the linear to zigzag phase transition is not destroyed by the strong quantum fluctuations present in narrow wires; it has characteristics which are qualitatively different from the classical transition.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Approaching Chemical Accuracy with Quantum Monte Carlo

A quantum Monte Carlo study of the atomization energies for the G2 set of molecules is presented. Basis size dependence of diffusion Monte Carlo atomization energies is studied with a single determinant Slater-Jastrow trial wavefunction formed from Hartree-Fock orbitals. With the largest basis set, the mean absolute deviation from experimental atomization energies for the G2 set is 3.0 kcal/mol. Optimizing the orbitals within variational Monte Carlo improves the agreement between diffusion Monte Carlo and experiment, reducing the mean absolute deviation to 2.1 kcal/mol. Moving beyond a single determinant Slater-Jastrow trial wavefunction, diffusion Monte Carlo with a small complete active space Slater-Jastrow trial wavefunction results in near chemical accuracy. In this case, the mean absolute deviation from experimental atomization energies is 1.2 kcal/mol. It is shown from calculations on systems containing phosphorus that the accuracy can be further improved by employing a larger active space.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Semistochastic Projector Monte Carlo Method

We introduce a semistochastic implementation of the power method to compute, for very large matrices, the dominant eigenvalue and expectation values involving the corresponding eigenvector. The method is semistochastic in that the matrix multiplication is partially implemented numerically exactly and partially with respect to expectation values only. Compared to a fully stochastic method, the semistochastic approach significantly reduces the computational time required to obtain the eigenvalue to a specified statistical uncertainty. This is demonstrated by the application of the semistochastic quantum Monte Carlo method to systems with a sign problem: the fermion Hubbard model and the carbon dimer.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Basis set construction for molecular electronic structure theory: Natural orbital and Gauss-Slater basis for smooth pseudpotentials

A simple yet general method for constructing basis sets for molecular electronic structure calculations is presented. These basis sets consist of atomic natural orbitals from a multi-configurational self-consistent field calculation supplemented with primitive functions, chosen such that the asymptotics are appropriate for the potential of the system. Primitives are optimized for the homonuclear diatomic molecule to produce a balanced basis set. Two general features that facilitate this basis construction are demonstrated. First, weak coupling exists between the optimal exponents of primitives with different angular momenta. Second, the optimal primitive exponents for a chosen system depend weakly on the particular level of theory employed for optimization. The explicit case considered here is a basis set appropriate for the Burkatzki-Filippi-Dolg pseudopotentials. Since these pseudopotentials are finite at nuclei and have a Coulomb tail, the recently proposed Gauss-Slater functions are the appropriate primitives. Double- and triple-zeta bases are developed for elements hydrogen through argon. These new bases offer significant gains over the corresponding Burkatzki-Filippi-Dolg bases at various levels of theory. Using a Gaussian expansion of the basis functions, these bases can be employed in any electronic structure method. Quantum Monte Carlo provides an added benefit: expansions are unnecessary since the integrals are evaluated numerically.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantum Monte Carlo facing the Hartree-Fock symmetry dilemma: The case of hydrogen rings

When using Hartree-Fock (HF) trial wave functions in quantum Monte Carlo calculations, one faces, in case of HF instabilities, the HF symmetry dilemma in choosing between the symmetry-adapted solution of higher HF energy and symmetry-broken solutions of lower HF energies. In this work, we have examined the HF symmetry dilemma in hydrogen rings which present singlet instabilities for sufficiently large rings. We have found that the symmetry-adapted HF wave function gives a lower energy both in variational Monte Carlo and in fixed-node diffusion Monte Carlo. This indicates that the symmetry-adapted wave function has more accurate nodes than the symmetry-broken wave functions, and thus suggests that spatial symmetry is an important criterion for selecting good trial wave functions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantum Monte Carlo with Jastrow-valence-bond wave functions

We consider the use in quantum Monte Carlo calculations of two types of valence bond wave functions based on strictly localized active orbitals, namely valence bond self-consistent-field (VBSCF) and breathing-orbital valence bond (BOVB) wave functions. Complemented by a Jastrow factor, these Jastrow-Valence-Bond wave functions are tested by computing the equilibrium well depths of the four diatomic molecules C$_2$, N$_2$, O$_2$, and F$_2$ in both variational Monte Carlo (VMC) and diffusion Monte Carlo (DMC). We show that it is possible to design compact wave functions based on chemical grounds that are capable of describing both static and dynamic electron correlation. These wave functions can be systematically improved by inclusion of valence bond structures corresponding to additional bonding patterns.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Compact and Flexible Basis Functions for Quantum Monte Carlo Calculations

Molecular calculations in quantum Monte Carlo frequently employ a mixed basis consisting of contracted and primitive Gaussian functions. While standard basis sets of varying size and accuracy are available in the literature, we demonstrate that reoptimizing the primitive function exponents within quantum Monte Carlo yields more compact basis sets for a given accuracy. Particularly large gains are achieved for highly excited states. For calculations requiring non-diverging pseudopotentials, we introduce Gauss-Slater basis functions that behave as Gaussians at short distances and Slaters at long distances. These basis functions further improve the energy and fluctuations of the local energy for a given basis size. Gains achieved by exponent optimization and Gauss-Slater basis use are exemplified by calculations for the ground state of carbon, the lowest lying excited states of carbon with $^5S^o$, $^3P^o$, $^1D^o$, $^3F^o$ symmetries, carbon dimer, and naphthalene. Basis size reduction enables quantum Monte Carlo treatment of larger molecules at high accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase transformation in Si from semiconducting diamond to metallic beta-Sn phase in QMC and DFT under hydrostatic and anisotropic stress

Silicon undergoes a phase transition from the semiconducting diamond phase to the metallic beta-Sn phase under pressure. We use quantum Monte Carlo calculations to predict the transformation pressure and compare the results to density functional calculations employing the LDA, PBE, PW91, WC, AM05, PBEsol and HSE06 exchange-correlation functionals. Diffusion Monte Carlo predicts a transition pressure of 14.0 +- 1.0 GPa slightly above the experimentally observed transition pressure range of 11.3 to 12.6 GPa. The HSE06 hybrid functional predicts a transition pressure of 12.4 GPa in excellent agreement with experiments. Exchange-correlation functionals using the local-density approximation and generalized-gradient approximations result in transition pressures ranging from 3.5 to 10.0 GPa, well below the experimental values. The transition pressure is sensitive to stress anisotropy. Anisotropy in the stress along any of the cubic axes of the diamond phase of silicon lowers the equilibrium transition pressure and may explain the discrepancy between the various experimental values as well as the small overestimate of the quantum Monte Carlo transition pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantum Monte Carlo calculations of electronic excitation energies: the case of the singlet $n \to π^*$ (CO) transition in acrolein

We report state-of-the-art quantum Monte Carlo calculations of the singlet $n \to π^*$ (CO) vertical excitation energy in the acrolein molecule, extending the recent study of Bouabça {\it et al.} [J. Chem. Phys. {\bf 130}, 114107 (2009)]. We investigate the effect of using a Slater basis set instead of a Gaussian basis set, and of using state-average versus state-specific complete-active-space (CAS) wave functions, with or without reoptimization of the coefficients of the configuration state functions (CSFs) and of the orbitals in variational Monte Carlo (VMC). It is found that, with the Slater basis set used here, both state-average and state-specific CAS(6,5) wave functions give an accurate excitation energy in diffusion Monte Carlo (DMC), with or without reoptimization of the CSF and orbital coefficients in the presence of the Jastrow factor. In contrast, the CAS(2,2) wave functions require reoptimization of the CSF and orbital coefficients to give a good DMC excitation energy. Our best estimates of the vertical excitation energy are between 3.86 and 3.89 eV.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Localization in an Inhomogeneous Quantum Wire

We study interaction-induced localization of electrons in an inhomogeneous quasi-one-dimensional system--a wire with two regions, one at low density and the other high. Quantum Monte Carlo techniques are used to treat the strong Coulomb interactions in the low density region, where localization of electrons occurs. The nature of the transition from high to low density depends on the density gradient--if it is steep, a barrier develops between the two regions, causing Coulomb blockade effects. Ferromagnetic spin polarization does not appear for any parameters studied. The picture emerging here is in good agreement with measurements of tunneling between two wires.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Optimization of quantum Monte Carlo wave functions by energy minimization

We study three wave function optimization methods based on energy minimization in a variational Monte Carlo framework: the Newton, linear and perturbative methods. In the Newton method, the parameter variations are calculated from the energy gradient and Hessian, using a reduced variance statistical estimator for the latter. In the linear method, the parameter variations are found by diagonalizing a non-symmetric estimator of the Hamiltonian matrix in the space spanned by the wave function and its derivatives with respect to the parameters, making use of a strong zero-variance principle. In the less computationally expensive perturbative method, the parameter variations are calculated by approximately solving the generalized eigenvalue equation of the linear method by a nonorthogonal perturbation theory. These general methods are illustrated here by the optimization of wave functions consisting of a Jastrow factor multiplied by an expansion in configuration state functions (CSFs) for the C$_2$ molecule, including both valence and core electrons in the calculation. The Newton and linear methods are very efficient for the optimization of the Jastrow, CSF and orbital parameters. The perturbative method is a good alternative for the optimization of just the CSF and orbital parameters. Although the optimization is performed at the variational Monte Carlo level, we observe for the C$_2$ molecule studied here, and for other systems we have studied, that as more parameters in the trial wave functions are optimized, the diffusion Monte Carlo total energy improves monotonically, implying that the nodal hypersurface also improves monotonically.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Energy and variance optimization of many body wave functions

We present a simple, robust and efficient method for varying the parameters in a many-body wave function to optimize the expectation value of the energy. The effectiveness of the method is demonstrated by optimizing the parameters in flexible Jastrow factors, that include 3-body electron-electron-nucleus correlation terms, for the NO$_2$ and decapentaene (C$_{10}$H$_{12}$) molecules. The basic idea is to add terms to the straightforward expression for the Hessian that are zero when the integrals are performed exactly, but that cancel much of the statistical fluctuations for a finite Monte Carlo sample. The method is compared to what is currently the most popular method for optimizing many-body wave functions, namely minimization of the variance of the local energy. The most efficient wave function is obtained by optimizing a linear combination of the energy and the variance.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

The maximum density droplet to lower density droplet transition in quantum dots

We show that, Landau level mixing in two-dimensional quantum dot wave functions can be taken into account very effectively by multiplying the exact lowest Landau level wave functions by a Jastrow factor which is optimized by variance minimization. The comparison between exact diagonalization and fixed phase diffusion Monte Carlo results suggests that the phase of the many-body wave functions are not affected much by Landau level mixing. We apply these wave functions to study the transition from the maximum density droplet state (incipient integer quantum Hall state with angular momentum L=N(N-1)/2) to lower density droplet states (L>N(N-1)/2).

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Monte Carlo Eigenvalue Methods in Quantum Mechanics and Statistical Mechanics

In this review we discuss, from a unified point of view, a variety of Monte Carlo methods used to solve eigenvalue problems in statistical mechanics and quantum mechanics. Although the applications of these methods differ widely, the underlying mathematics is quite similar in that they are stochastic implementations of the power method. In all cases, optimized trial states can be used to reduce the errors of Monte Carlo estimates.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Monte Carlo Optimization of Trial Wave Functions in Quantum Mechanics and Statistical Mechanics

This review covers applications of quantum Monte Carlo methods to quantum mechanical problems in the study of electronic and atomic structure, as well as applications to statistical mechanical problems both of static and dynamic nature. The common thread in all these applications is optimization of many-parameter trial states, which is done by minimization of the variance of the local or, more generally for arbitrary eigenvalue problems, minimization of the variance of the configurational eigenvalue.

0 citations0 reviewsNo code linkedOpen access