Rydberg states with quantum Monte Carlo

Calculations on Rydberg states are performed using quantum Monte Carlo methods. Excitation energies and singlet-triplet splittings are calculated for two model systems, the carbon atom (3P and 1P) and carbon monoxide ((1Sigma and 3Sigma). Kohn-Sham wave functions constructed from open-shell localized Hartree-Fock orbitals are used as trial and guide functions. The fixed-node diffusion quantum Monte Carlo (FN-DMC) method depends strongly on the wave function's nodal hypersurface. Nodal artefacts are investigated for the ground state of the carbon atom. Their effect on the FN-DMC results can be analyzed quantitatively. FN-DMC leads to accurate excitation energies but to less accurate singlet-triplet splittings. Variational Monte Carlo calculations are able to reproduce the experimental results for both the excitation energies and the singlet-triplet splittings.     

An optimized initialization algorithm to ensure accuracy in quantum Monte Carlo calculations

Quantum Monte Carlo (QMC) calculations require the generation of random electronic configurations with respect to a desired probability density, usually the square of the magnitude of the wavefunction. In most cases, the Metropolis algorithm is used to generate a sequence of configurations in a Markov chain. This method has an inherent equilibration phase, during which the configurations are not representative of the desired density and must be discarded. If statistics are gathered before the walkers have equilibrated, contamination by nonequilibrated configurations can greatly reduce the accuracy of the results. Because separate Markov chains must be equilibrated for the walkers on each processor, the use of a long equilibration phase has a profoundly detrimental effect on the efficiency of large parallel calculations. The stratified atomic walker initialization (STRAW) shortens the equilibration phase of QMC calculations by generating statistically independent electronic configurations in regions of high probability density. This ensures the accuracy of calculations by avoiding contamination by nonequilibrated configurations. Shortening the length of the equilibration phase also results in significant improvements in the efficiency of parallel calculations, which reduces the total computational run time. For example, using STRAW rather than a standard initialization method in 512 processor calculations reduces the amount of time needed to calculate the energy expectation value of a trial function for a molecule of the energetic material RDX to within 0.01 au by 33%.     

Dynamic correlations with time-dependent quantum Monte Carlo

In this paper, we solve quantum many-body problem by propagating ensembles of trajectories and guiding waves in physical space. We introduce the "effective potential" correction within the recently proposed time-dependent quantum Monte Carlo methodology to incorporate the nonlocal quantum correlation effects between the electrons. The associated correlation length is calculated by adaptive kernel density estimation over the walker distribution. The general formalism is developed and tested on one-dimensional helium atom in laser field of different intensities and carrier frequencies. Good agreement with exact results for the atomic ionization is obtained.     

Vanadium oxide compounds with quantum Monte Carlo

Calculations with the diffusion quantum Monte Carlo method are presented for vanadium oxide molecules VO0/+0(n) with n = 1-4 and for V2O5. Atomization and ionization energies are calculated as well as oxygen abstraction energies. The fixed-node approximation is compared for guide functions with orbitals from B3LYP and BP86 calculations and higher accuracy was obtained with the latter orbitals. Additionally, all-electron and pseudopotential calculations are compared for the oxygen atom. The overall accuracy is found to be comparable to CCSD(T) calculations where experimental data is available.     

Bond breaking with auxiliary-field quantum Monte Carlo

Bond stretching mimics different levels of electron correlation and provides a challenging test bed for approximate many-body computational methods. Using the recently developed phaseless auxiliary-field quantum Monte Carlo (AF QMC) method, we examine bond stretching in the well-studied molecules BH and N(2) and in the H(50) chain. To control the sign/phase problem, the phaseless AF QMC method constrains the paths in the auxiliary-field path integrals with an approximate phase condition that depends on a trial wave function. With single Slater determinants from unrestricted Hartree-Fock as trial wave function, the phaseless AF QMC method generally gives better overall accuracy and a more uniform behavior than the coupled cluster CCSD(T) method in mapping the potential-energy curve. In both BH and N(2), we also study the use of multiple-determinant trial wave functions from multiconfiguration self-consistent-field calculations. The increase in computational cost versus the gain in statistical and systematic accuracy are examined. With such trial wave functions, excellent results are obtained across the entire region between equilibrium and the dissociation limit.     

DFT-LDA pseudopotentials in quantum Monte Carlo

We investigate the portability of standard norm-conserving pseudopotentials outside the density functional theory-local density approximation (DFT-LDA) framework, i.e., their use and interpretation as electron-ion effective potentials in valence-only diffusion Monte Carlo simulations. While first-principles many-body pseudopotentials are not available in the literature yet, the use of approximate pseudopotentials in quantum Monte Carlo simulations is becoming widespread. Here we attempt a systematic analysis of the portability of norm-conserving pseudopotentials generated within DFT-LDA, focusing on a model many-body system, the two-electron valence-only ion. Our results indicate that the portability is good in most cases, hence the use of pseudopotentials in quantum Monte Carlo simulations is in general a reasonable approximation but suggest that in some cases this approximation may be relevant. © 1997 John Wiley & Sons, Inc.     

Prospects for release-node quantum Monte Carlo

We perform release-node quantum Monte Carlo simulations on the first row diatomic molecules in order to assess how accurately their ground-state energies can be obtained. An analysis of the fermion-boson energy difference is shown to be strongly dependent on the nuclear charge, Z, which in turn determines the growth of variance of the release-node energy. It is possible to use maximum entropy analysis to extrapolate to ground-state energies only for the low Z elements. For the higher Z dimers beyond boron, the error growth is too large to allow accurate data for long enough imaginary times. Within the limit of our statistics we were able to estimate, in atomic units, the ground-state energy of Li(2) (-14.9947(1)), Be(2) (-29.3367(7)), and B(2)(-49.410(2)).     

Calculating expectations with time-dependent perturbations in quantum Monte Carlo

We show that a small perturbation periodic in imaginary time can be used to compute expectation values of nondifferential operators that do not commute with the Hamiltonian within the framework of quantum diffusion Monte Carlo. Some results for the harmonic oscillator and the helium atom are presented showing the validity of the proposed method.     

The accuracy of diffusion quantum Monte Carlo simulations in the determination of molecular equilibrium structures

For a test set of 17 first-row small molecules, the equilibrium structures are calculated with Ornstein-Uhlenbeck diffusion quantum Monte Carlo simulations guiding by trial wave functions constructed from floating spherical Gaussian orbitals and spherical Gaussian geminals. To measure performance of the Monte Carlo calculations, the mean deviation, the mean absolute deviation, the maximum absolute deviation, and the standard deviation of Monte Carlo calculated equilibrium structures with respect to empirical equilibrium structures are given. This approach is found to yield results having a uniformly high quality, being consistent with empirical equilibrium structures and surpassing calculated values from the coupled cluster model with single, double, and noniterative triple excitations [CCSD(T)] with the basis sets of cc-pCVQZ and cc-pVQZ. The nonrelativistic equilibrium atomization energies are also presented to assess performance of the calculated methods. The mean absolute deviations regarding experimental atomization energy are 0.16 and 0.21 kcal/mol for the Monte Carlo and CCSD(T)/cc-pCV(56)Z calculations, respectively.     

Energy-consistent pseudopotentials for quantum Monte Carlo calculations

The authors present scalar-relativistic energy-consistent Hartree-Fock pseudopotentials for the main-group elements. The pseudopotentials do not exhibit a singularity at the nucleus and are therefore suitable for quantum Monte Carlo (QMC) calculations. They demonstrate their transferability through extensive benchmark calculations of atomic excitation spectra as well as molecular properties. In particular, they compute the vibrational frequencies and binding energies of 26 first- and second-row diatomic molecules using post-Hartree-Fock methods, finding excellent agreement with the corresponding all-electron values. They also show their pseudopotentials give superior accuracy than other existing pseudopotentials constructed specifically for QMC. Finally, valence basis sets of different sizes (VnZ with n=D,T,Q,5 for first and second rows, and n=D,T for third to fifth rows) optimized for our pseudopotentials are also presented.     

Time-dependent quantum Monte Carlo and the stochastic quantization

We examine the relation between the recently proposed time-dependent quantum Monte Carlo (TDQMC) method and the principles of stochastic quantization. In both TDQMC and stochastic quantization, particle motion obeys stochastic guidance equations to preserve quantum equilibrium. In this way the probability density of the Monte Carlo particles corresponds to the modulus square of the many-body wave function at all times. However, in TDQMC, the motion of particles and guide waves occurs in physical space unlike in stochastic quantization where it occurs in configuration space. Hence, the practical calculation of time evolution of many-body fully correlated quantum systems becomes feasible within the TDQMC methodology. We illustrate the TDQMC technique by calculating the symmetric and antisymmetric ground state of a model one-dimensional helium atom, and the time evolution of the dipole moment when the atom is irradiated by a strong ultrashort laser pulse.     

Correlated electron dynamics with time-dependent quantum Monte Carlo: three-dimensional helium

Here the recently proposed time-dependent quantum Monte Carlo method is applied to three dimensional para- and ortho-helium atoms subjected to an external electromagnetic field with amplitude sufficient to cause significant ionization. By solving concurrently sets of up to 20,000 coupled 3D time-dependent Schro?dinger equations for the guide waves and corresponding sets of first order equations of motion for the Monte Carlo walkers we obtain ground state energies in close agreement with the exact values. The combined use of spherical coordinates and B-splines along the radial coordinate proves to be especially accurate and efficient for such calculations. Our results for the dipole response and the ionization of an atom with un-correlated electrons are in good agreement with the predictions of the conventional time-dependent Hartree-Fock method while the calculations with correlated electrons show enhanced ionization that is due to the electron-electron repulsion.     

Assessing the accuracy of quantum Monte Carlo and density functional theory for energetics of small water clusters

We present a detailed study of the energetics of water clusters (H(2)O)(n) with n ≤ 6, comparing diffusion Monte Carlo (DMC) and approximate density functional theory (DFT) with well converged coupled-cluster benchmarks. We use the many-body decomposition of the total energy to classify the errors of DMC and DFT into 1-body, 2-body and beyond-2-body components. Using both equilibrium cluster configurations and thermal ensembles of configurations, we find DMC to be uniformly much more accurate than DFT, partly because some of the approximate functionals give poor 1-body distortion energies. Even when these are corrected, DFT remains considerably less accurate than DMC. When both 1- and 2-body errors of DFT are corrected, some functionals compete in accuracy with DMC; however, other functionals remain worse, showing that they suffer from significant beyond-2-body errors. Combining the evidence presented here with the recently demonstrated high accuracy of DMC for ice structures, we suggest how DMC can now be used to provide benchmarks for larger clusters and for bulk liquid water.     

Electron-nucleus cusp correction and forces in quantum Monte Carlo

A simple method is presented which ensures the electron-nucleus cusp condition is satisfied by the Slater-Jastrow wavefunctions commonly employed in quantum Monte Carlo simulations. The method is applied in variational energy calculations of the neon atom and a selection of molecules using both Gaussian and Slater basis sets. In addition, we discuss the relationship between the electron-nucleus cusps and the variance of forces, and investigate the sensitivity of forces to the quality of the cusps for various diatomic molecules.     

Optimum and efficient sampling for variational quantum Monte Carlo

Quantum mechanics for many-body systems may be reduced to the evaluation of integrals in 3N dimensions using Monte Carlo, providing the Quantum Monte Carlo ab initio methods. Here we limit ourselves to expectation values for trial wave functions, that is to variational quantum Monte Carlo. Almost all previous implementations employ samples distributed as the physical probability density of the trial wave function, and assume the central limit theorem to be valid. In this paper we provide an analysis of random error in estimation and optimization that leads naturally to new sampling strategies with improved computational and statistical properties. A rigorous lower limit to the random error is derived, and an efficient sampling strategy presented that significantly increases computational efficiency. In addition the infinite variance heavy tailed random errors of optimum parameters in conventional methods are replaced with a Normal random error, strengthening the theoretical basis of optimization. The method is applied to a number of first row systems and compared with previously published results.     

Random number generators tested on quantum Monte Carlo simulations

We have tested and compared several (pseudo) random number generators (RNGs) applied to a practical application, ground state energy calculations of molecules using variational and diffusion Monte Carlo metheds. A new multiple recursive generator with 8th‐order recursion (MRG8) and the Mersenne twister generator (MT19937) are tested and compared with the RANLUX generator with five luxury levels (RANLUX‐[0–4]). Both MRG8 and MT19937 are proven to give the same total energy as that evaluated with RANLUX‐4 (highest luxury level) within the statistical error bars with less computational cost to generate the sequence. We also tested the notorious implementation of linear congruential generator (LCG), RANDU, for comparison. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Kinetic Monte Carlo modeling of chemical reactions coupled with heat transfer

In this paper, we describe two types of effective events for describing heat transfer in a kinetic Monte Carlo (KMC) simulation that may involve stochastic chemical reactions. Simulations employing these events are referred to as KMC-TBT and KMC-PHE. In KMC-TBT, heat transfer is modeled as the stochastic transfer of "thermal bits" between adjacent grid points. In KMC-PHE, heat transfer is modeled by integrating the Poisson heat equation for a short time. Either approach is capable of capturing the time dependent system behavior exactly. Both KMC-PHE and KMC-TBT are validated by simulating pure heat transfer in a rod and a square and modeling a heated desorption problem where exact numerical results are available. KMC-PHE is much faster than KMC-TBT and is used to study the endothermic desorption of a lattice gas. Interesting findings from this study are reported.     

Auxiliary-field quantum Monte Carlo calculations of molecular systems with a Gaussian basis

We extend the recently introduced phaseless auxiliary-field quantum Monte Carlo (QMC) approach to any single-particle basis and apply it to molecular systems with Gaussian basis sets. QMC methods in general scale favorably with the system size as a low power. A QMC approach with auxiliary fields, in principle, allows an exact solution of the Schrodinger equation in the chosen basis. However, the well-known sign/phase problem causes the statistical noise to increase exponentially. The phaseless method controls this problem by constraining the paths in the auxiliary-field path integrals with an approximate phase condition that depends on a trial wave function. In the present calculations, the trial wave function is a single Slater determinant from a Hartree-Fock calculation. The calculated all-electron total energies show typical systematic errors of no more than a few millihartrees compared to exact results. At equilibrium geometries in the molecules we studied, this accuracy is roughly comparable to that of coupled cluster with single and double excitations and with noniterative triples [CCSD(T)]. For stretched bonds in H(2)O, our method exhibits a better overall accuracy and a more uniform behavior than CCSD(T).