Search

We modify the reweighting factor of the projector used in diffusion Monte Carlo to reduce the time-step error of the total energy. Further, we present a reweighting scheme that has the desirable feature that it is exactly size-consistent, i.e, the energy of a system containing widely separated fragments is the same as the sum of the energies of the individual fragments. The practical utility of the latter improvement is that it reduces the time-step error of the binding energies of some weakly interacting systems., Comment: 15 pages, 3 figures; expanded discussion of results, discussion of time-step error for all-electron systems added to supplementary material

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.

The semistochastic heat-bath configuration interaction (SHCI) method is a selected configuration interaction plus perturbation theory method that has provided near-full configuration interaction (FCI) levels of accuracy for many systems with both single- and multi-reference character. However, obtaining accurate energies in the complete basis set limit is hindered by the slow convergence of the FCI energy with respect to basis size. Here we show that the recently developed basis-set correction method based on range-separated density-functional theory can be used to significantly speed up basis-set convergence in SHCI calculations. In particular, we study two such schemes that differ in the functional used, and apply them to transition metal atoms and monoxides to obtain total, ionization, and dissociation energies well converged to the complete-basis-set limit within chemical accuracy., Comment: 5 figures, supplemental material

We study several approaches to orbital optimization in selected configuration interaction plus perturbation theory (SCI+PT) methods, and test them on the ground and excited states of three molecules using the semistochastic heatbath configuration interaction (SHCI) method. We discuss the ways in which the orbital optimization problem in SCI resembles and differs from that in complete active space self-consistent field (CASSCF). Starting from natural orbitals, these approaches divide into three classes of optimization methods according to how they treat coupling between configuration interaction (CI) coefficients and orbital parameters, namely uncoupled, fully coupled, and quasi-fully coupled methods. We demonstrate that taking the coupling into account is crucial for fast convergence and recommend two quasi-fully coupled methods for such applications: accelerated diagonal Newton and Broyden-Fletcher-Goldfarb-Shanno (BFGS)., Comment: 5 figures

We present a version of the T-moves approach for treating nonlocal pseudopotentials in diffusion Monte Carlo which has much smaller time-step errors than the existing T-moves approaches, while at the same time preserving desirable features such as the upper-bound property for the energy. In addition, we modify the reweighting factor of the projector used in diffusion Monte Carlo to reduce the time-step error. The latter is applicable not only to pseudopotential calculations but to all-electron calculations as well., Comment: 25 pages, 4 figures; additional figures added, expanded discussion of results, added VMC energies to suppl.pdf; We discovered a programming error which had led us to recommend a reweighting factor which gives a suboptimal time-step error. The reweighting factor in Eq. 9 has been replaced with a better recommendation, and the results have been updated. The conclusions of the paper are unchanged

We report on the findings of a blind challenge devoted to determining the frozen-core, full configuration interaction (FCI) ground state energy of the benzene molecule in a standard correlation-consistent basis set of double-$\zeta$ quality. As a broad international endeavour, our suite of wave function-based correlation methods collectively represents a diverse view of the high-accuracy repertoire offered by modern electronic structure theory. In our assessment, the evaluated high-level methods are all found to qualitatively agree on a final correlation energy, with most methods yielding an estimate of the FCI value around $-863$ m$E_{\text{H}}$. However, we find the root-mean-square deviation of the energies from the studied methods to be considerable (1.3 m$E_{\text{H}}$), which in light of the acclaimed performance of each of the methods for smaller molecular systems clearly displays the challenges faced in extending reliable, near-exact correlation methods to larger systems. While the discrepancies exposed by our study thus emphasize the fact that the current state-of-the-art approaches leave room for improvement, we still expect the present assessment to provide a valuable community resource for benchmark and calibration purposes going forward., Comment: 29 pages, 1 figure, 2 tables. SI as an ancillary file

The recently developed semistochastic heat-bath configuration interaction (SHCI) method is a systematically improvable selected configuration interaction plus perturbation theory method capable of giving essentially exact energies for larger systems than is possible with other such methods. We compute SHCI atomization energies for 55 molecules which have been used as a test set in prior studies because their atomization energies are known from experiment. Basis sets from cc-pVDZ to cc-pV5Z are used, totaling up to 500 orbitals and a Hilbert space of $10^{32}$ Slater determinants for the largest molecules. For each basis, an extrapolated energy well within chemical accuracy (1 kcal/mol or 1.6 mHa/mol) of the exact energy for that basis is computed using only a tiny fraction of the entire Hilbert space. We also use our almost exact energies to benchmark coupled-cluster [CCSD(T)] energies. The energies are extrapolated to the complete basis set limit and compared to the experimental atomization energies. The extrapolations are done both without and with a basis-set correction based on density-functional theory. The mean absolute deviations from experiment for these extrapolations are 0.46 kcal/mol and 0.51 kcal/mol, respectively. Orbital optimization methods used to obtain improved convergence of the SHCI energies are also discussed., Comment: 6 figures, 4 supplemental files

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., Comment: 9 pages, 6 figures, 1 table

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., Comment: 4 figures

We report that a recent active space model of the nitrogenase FeMo cofactor, proposed in the context of quantum simulations, is not representative of the electronic structure of the FeMo cofactor ground-state. Although quantum resource estimates, outside of the cost of adiabatic state preparation, will not be much affected, conclusions should not be drawn from the complexity of classical simulations of the electronic structure of this system in this active space. We provide a different model active space for the FeMo cofactor that contains the basic open-shell qualitative character, which may be useful as a benchmark system for making classical and quantum resource estimates., Comment: 3 pages, 3 figures

We extend our recently-developed heat-bath configuration interaction (HCI) algorithm, and our semistochastic algorithm for performing multireference perturbation theory, to the calculation of excited-state wavefunctions and energies. We employ time-reversal symmetry, which reduces the memory requirements by more than a factor of two. An extrapolation technique is introduced to reliably extrapolate HCI energies to the Full CI limit. The resulting algorithm is used to compute the twelve lowest-lying potential energy surfaces of the carbon dimer using the cc-pV5Z basis set, with an estimated error in energy of 30-50 {\mu}Ha compared to Full CI. The excitation energies obtained using our algorithm have a mean absolute deviation of 0.02 eV compared to experimental values. We also calculate the complete active-space (CAS) energies of the S0, S1, and T0 states of tetracene, which are of relevance to singlet fission, by fully correlating active spaces as large as 18 electrons in 36 orbitals.

We present the extension of variational Monte Carlo (VMC) to the calculation of electronic excitation energies and oscillator strengths using time-dependent linear-response theory. By exploiting the analogy existing between the linear method for wave-function optimisation and the generalised eigenvalue equation of linear-response theory, we formulate the equations of linear-response VMC (LR-VMC). This LR-VMC approach involves the first-and second-order derivatives of the wave function with respect to the parameters. We perform first tests of the LR-VMC method within the Tamm-Dancoff approximation using single-determinant Jastrow-Slater wave functions with different Slater basis sets on some singlet and triplet excitations of the beryllium atom. Comparison with reference experimental data and with configuration-interaction-singles (CIS) results shows that LR-VMC generally outperforms CIS for excitation energies and is thus a promising approach for calculating electronic excited-state properties of atoms and molecules., Comment: Advances in Quantum Chemistry, 2017, Novel Electronic Structure Theory: General Innovations and Strongly Correlated Systems

We extend the recently proposed heat-bath configuration interaction (HCI) method [Holmes, Tubman, Umrigar, J. Chem. Theory Comput. 12, 3674 (2016)], by introducing a semistochastic algorithm for performing multireference Epstein-Nesbet perturbation theory, in order to completely eliminate the severe memory bottleneck of the original method. The proposed algorithm has several attractive features. First, there is no sign problem that plagues several quantum Monte Carlo methods. Second, instead of using Metropolis-Hastings sampling, we use the Alias method to directly sample determinants from the reference wavefunction, thus avoiding correlations between consecutive samples. Third, in addition to removing the memory bottleneck, semistochastic HCI (SHCI) is faster than the deterministic variant for many systems if a stochastic error of 0.1 mHa is acceptable. Fourth, within the SHCI algorithm one can trade memory for a modest increase in computer time. Fifth, the perturbative calculation is embarrassingly parallel. The SHCI algorithm extends the range of applicability of the original algorithm, allowing us to calculate the correlation energy of very large active spaces. We demonstrate this by performing calculations on several first row dimers including F2 with an active space of (14e, 108o), Mn-Salen cluster with an active space of (28e, 22o), and Cr2 dimer with up to a quadruple-zeta basis set with an active space of (12e, 190o). For these systems we were able to obtain better than 1 mHa accuracy with a wall time of merely 55 seconds, 37 seconds, and 56 minutes on 1, 1, and 4 nodes, respectively.

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.

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., Comment: Advances in Quantum Chemistry, 2015, Electron Correlation in Molecules -- ab initio Beyond Gaussian Quantum Chemistry, pp.0000

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\Sigma^-_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\"odinger equation in small molecules.

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., Comment: Accepted for the JCP Special Topic Issue "Advances in DFT Methodology"

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\`e 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.

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., Comment: 5 pages, 5 figures. Published version: several clarifications, such as comment on finite size effects

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., Comment: 5 pages, 5 figures

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., Comment: 6 pages, 5 figures

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., Comment: 5 pages, 4 figures, suggested PACS numbers 02.70.Ss, 71.10.Fd, 31.15.-p

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., Comment: 6 pages, 3 figures, 2 tables, to appear in "Advances in Quantum Monte Carlo", ACS

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.

We report state-of-the-art quantum Monte Carlo calculations of the singlet $n \to \pi^*$ (CO) vertical excitation energy in the acrolein molecule, extending the recent study of Bouab\c{c}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., Comment: 6 pages, 1 figure, 2 tables, to appear in Progress in Theoretical Chemistry and Physics

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.

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., Comment: 9 pages, 7 figures

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., Comment: 8 Pages, 2 Figures, 9 Tables

The accuracy of the energy landscape of silicon systems obtained from various density functional methods, a tight binding scheme and force fields is studied. Quantum Monte Carlo results serve as quasi exact reference values. In addition to the well known accuracy of DFT methods for geometric ground states and metastable configurations we find that DFT methods give a similar accuracy for transition states and thus a good overall description of the energy landscape. On the other hand, force fields give a very poor description of the landscape that are in most cases too rugged and contain many fake local minima and saddle points or ones that have the wrong height.

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., Comment: 4 pages; 2 new figures, substantial revisions and clarifications

We pursue the development and application of the recently-introduced linear optimization method for determining the optimal linear and nonlinear parameters of Jastrow-Slater wave functions in a variational Monte Carlo framework. In this approach, the optimal parameters are found iteratively by diagonalizing the Hamiltonian matrix in the space spanned by the wave function and its first-order derivatives, making use of a strong zero-variance principle. We extend the method to optimize the exponents of the basis functions, simultaneously with all the other parameters, namely the Jastrow, configuration state function and orbital parameters. We show that the linear optimization method can be thought of as a so-called augmented Hessian approach, which helps explain the robustness of the method and permits us to extend it to minimize a linear combination of the energy and the energy variance. We apply the linear optimization method to obtain the complete ground-state potential energy curve of the C_2 molecule up to the dissociation limit, and discuss size consistency and broken spin-symmetry issues in quantum Monte Carlo calculations. We perform calculations of the first-row atoms and homonuclear diatomic molecules with fully optimized Jastrow-Slater wave functions, and we demonstrate that molecular well depths can be obtained with near chemical accuracy quite systematically at the diffusion Monte Carlo level for these systems., Comment: 15 pages, 3 figures, to appear in Journal of Chemical Physics

Diffusion Monte Carlo (DMC) calculations are performed on the monocyclic and bicyclic forms of m-benzyne, which are the equilibrium structures at the CCSD(T) and CCSD levels of coupled cluster theory. We employed multi-configuration self-consistent field trial wave functions which are constructed from a carefully selected 8-electrons-in-8-orbitals complete active space [CAS(8,8)], with CSF coefficients that are reoptimized in the presence of a Jastrow factor. The DMC calculations show that the monocyclic structure is lower in energy than the bicyclic structure by 1.9(2) kcal/mole, in excellent agreement with the best coupled cluster results., Comment: 5 pages, 2 figures. to be published in JCP

We argue that Coulomb blockade phenomena are a useful probe of the cross-over to strong correlation in quantum dots. Through calculations at low density using variational and diffusion quantum Monte Carlo (up to r_s ~ 55), we find that the addition energy shows a clear progression from features associated with shell structure to those caused by commensurability of a Wigner crystal. This cross-over (which occurs near r_s ~ 20 for spin-polarized electrons) is, then, a signature of interaction-driven localization. As the addition energy is directly measurable in Coulomb blockade conductance experiments, this provides a direct probe of localization in the low density electron gas., Comment: 4 pages, published version, revised discussion

We construct improved quantum Monte Carlo estimators for the spherically- and system-averaged electron pair density (i.e. the probability density of finding two electrons separated by a relative distance u), also known as the spherically-averaged electron position intracule density I(u), using the general zero-variance zero-bias principle for observables, introduced by Assaraf and Caffarel. The calculation of I(u) is made vastly more efficient by replacing the average of the local delta-function operator by the average of a smooth non-local operator that has several orders of magnitude smaller variance. These new estimators also reduce the systematic error (or bias) of the intracule density due to the approximate trial wave function. Used in combination with the optimization of an increasing number of parameters in trial Jastrow-Slater wave functions, they allow one to obtain well converged correlated intracule densities for atoms and molecules. These ideas can be applied to calculating any pair-correlation function in classical or quantum Monte Carlo calculations., Comment: 13 pages, 9 figures, published version

We study the development of electron-electron correlations in circular quantum dots as the density is decreased. We consider a wide range of both electron number, N<=20, and electron gas parameter, r_s<18, using the diffusion quantum Monte Carlo technique. Features associated with correlation appear to develop very differently in quantum dots than in bulk. The main reason is that translational symmetry is necessarily broken in a dot, leading to density modulation and inhomogeneity. Electron-electron interactions act to enhance this modulation ultimately leading to localization. This process appears to be completely smooth and occurs over a wide range of density. Thus there is a broad regime of ``incipient'' Wigner crystallization in these quantum dots. Our specific conclusions are: (i) The density develops sharp rings while the pair density shows both radial and angular inhomogeneity. (ii) The spin of the ground state is consistent with Hund's (first) rule throughout our entire range of r_s for all 4

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., Comment: 18 pages, 8 figures, final version

We present a simple, robust and highly efficient method for optimizing all parameters of many-body wave functions in quantum Monte Carlo calculations, applicable to continuum systems and lattice models. Based on a strong zero-variance principle, diagonalization of the Hamiltonian matrix in the space spanned by the wav e function and its derivatives determines the optimal parameters. It systematically reduces the fixed-node error, as demonstrated by the calculation of the binding energy of the small but challenging C$_2$ molecule to the experimental accuracy of 0.02 eV.

Properties of the "electron gas" - in which conduction electrons interact by means of Coulomb forces but ionic potentials are neglected - change dramatically depending on the balance between kinetic energy and Coulomb repulsion. The limits are well understood. For very weak interactions (high density), the system behaves as a Fermi liquid, with delocalized electrons. In contrast, in the strongly interacting limit (low density), the electrons localize and order into a Wigner crystal phase. The physics at intermediate densities, however, remains a subject of fundamental research. Here, we study the intermediate-density electron gas confined to a circular disc, where the degree of confinement can be tuned to control the density. Using accurate quantum Monte Carlo techniques, we show that the electron-electron correlation induced by an increase of the interaction first smoothly causes rings, and then angular modulation, without any signature of a sharp transition in this density range. This suggests that inhomogeneities in a confined system, which exist even without interactions, are significantly enhanced by correlations., Comment: final version, modified introduction and clarifications, 4 pages

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)., Comment: 8 pages, 5 figures, accepted for publication in Phys. Rev. B

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., Comment: 4 pages, 4 figures, minor corrections of inexact statements, missing 2

We address the issue of accurately treating interaction effects in the mesoscopic regime by investigating the ground state properties of isolated irregular quantum dots. Quantum Monte Carlo techniques are used to calculate the distributions of ground state spin and addition energy. We find a reduced probability of high spin and a somewhat larger even/odd alternation in the addition energy from quantum Monte Carlo than in local spin density functional theory. In both approaches, the even/odd effect gets smaller with increasing number of electrons, contrary to the theoretical understanding of large dots. We argue that the local spin density approximation over predicts the effects of interactions in quantum dots., Comment: Final Version, to appear in PRB as a Rapid Comm

Several of the energies in Tables I and II of our previous paper are incorrect. Most of the errors are due to incorrectly inputting the symmetry of some of the states, others are due to incorrectly converting or transcribing the energies. In particular, for N=4 we had a near degeneracy and a violation of Hund's first rule. The corrected result has $|L=0,S=1>$ as the ground state, as predicted by Hund's rule. Hund's first rule is satisfied for all $N$ for the confining potential used. We present corrected versions of Tables I and II., Comment: 2 pages. Erratum to cond-mat/9912166 (PRB62,8120 (2000)), submitted to PRB

We present ground and excited state energies obtained from Diffusion Monte Carlo (DMC) calculations, using accurate multiconfiguration wave functions, for $N$ electrons ($N\le13$) confined to a circular quantum dot. We analyze the electron-electron pair correlation functions and compare the density and correlation energies to the predictions of local spin density approximation theory (LSDA). The DMC estimated change in electrochemical potential as function of the number of electrons in the dot is compared to that from LSDA and Hartree-Fock (HF) calculations., Comment: 7 pages, 4 eps figures. To be published in Phys. Rev. B, September 15th 2000. See erratum cond-mat/0305710

In order to find the equilibrium geometries of molecules and solids and to perform ab initio molecular dynamics, it is necessary to calculate the forces on the nuclei. We present a correlated sampling method to efficiently calculate numerical forces and potential energy surfaces in diffusion Monte Carlo. It employs a novel coordinate transformation, earlier used in variational Monte Carlo, to greatly reduce the statistical error. Results are presented for first-row diatomic molecules., Comment: 5 pages, 2 postscript figures

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.

We consider two perturbative schemes to calculate excitation energies, each employing the Kohn-Sham Hamiltonian as the unperturbed system. Using accurate exchange-correlation potentials generated from essentially exact densities and their exchange components determined by a recently proposed method, we evaluate energy differences between the ground state and excited states in first-order perturbation theory for the Helium, ionized Lithium and Beryllium atoms. It was recently observed that the zeroth-order excitations energies, simply given by the difference of the Kohn-Sham eigenvalues, almost always lie between the singlet and triplet experimental excitations energies, corrected for relativistic and finite nuclear mass effects. The first-order corrections provide about a factor of two improvement in one of the perturbative schemes but not in the other. The excitation energies within perturbation theory are compared to the excitations obtained within $\Delta$SCF and time-dependent density functional theory. We also calculate the excitation energies in perturbation theory using approximate functionals such as the local density approximation and the optimized effective potential method with and without the Colle-Salvetti correlation contribution.

Compact and accurate wave functions can be constructed by quantum Monte Carlo methods. Typically, these wave functions consist of a sum of a small number of Slater determinants multiplied by a Jastrow factor. In this paper we study the importance of including high-order, nucleus-three-electron correlations in the Jastrow factor. An efficient algorithm based on the theory of invariants is used to compute the high-body correlations. We observe significant improvements in the variational Monte Carlo energy and in the fluctuations of the local energies but not in the fixed-node diffusion Monte Carlo energies. Improvements for the ground states of physical, fermionic atoms are found to be smaller than those for the ground states of fictitious, bosonic atoms, indicating that errors in the nodal surfaces of the fermionic wave functions are a limiting factor., Comment: 9 pages, no figures, Latex

We calculate the electronic structure of several atoms and small molecules by direct minimization of the Self-Interaction Corrected Local Density Approximation (SIC-LDA) functional. To do this we first derive an expression for the gradient of this functional under the constraint that the orbitals be orthogonal and show that previously given expressions do not correctly incorporate this constraint. In our atomic calculations the SIC-LDA yields total energies, ionization energies and charge densities that are superior to results obtained with the Local Density Approximation (LDA). However, for molecules SIC-LDA gives bond lengths and reaction energies that are inferior to those obtained from LDA. The nonlocal BLYP functional, which we include as a representative GGA functional, outperforms both LDA and SIC-LDA for all ground state properties we considered., Comment: 14 pages, 5 figures

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., Comment: 27 pages to appear in " Recent Advances in Quantum Monte Carlo Methods" edited by W.A. Lester