arXiv:cond-mat/9804010 · Physical review. B, Condensed matter

Quantum Monte Carlo calculations of the one-body density matrix and excitation energies of silicon

Paul R. C. Kent, Randolph Q. Hood, Mike D. Towler, Richard J. Needs, Gunaretnam Rajagopal

Published Jun 15, 1998• 44 citations

Abstract

Quantum Monte Carlo (QMC) techniques are used to calculate the one-body density matrix and excitation energies for the valence electrons of bulk silicon. The one-body density matrix and energies are obtained from a Slater-Jastrow wave function with a determinant of local-density approximation (LDA) orbitals. The QMC density matrix evaluated in a basis of LDA orbitals is strongly diagonally dominant. The natural orbitals obtained by diagonalizing the QMC density matrix resemble the LDA orbitals very closely. Replacing the determinant of LDA orbitals in the wave function by a determinant of natural orbitals makes no significant difference to the quality of the wave function's nodal surface, leaving the diffusion Monte Carlo energy unchanged. The extended Koopmans's theorem for correlated wave functions is used to calculate excitation energies for silicon, which are in reasonable agreement with the available experimental data. A diagonal approximation to the theorem, evaluated in the basis of LDA orbitals, works quite well for both the quasihole and quasielectron states. We have found that this approximation has an advantageous scaling with system size, allowing more efficient studies of larger systems.

DOI: 10.1103/physrevb.57.15293

Related papers

Phase-free quantum Monte Carlo method: random walks using general basis sets

Henry Krakauer

Fermion quantum Monte Carlo (QMC) methods that work in a general basis space may be more effective for some problems than the traditional diffusion QMC method. Projection of the ground state energy is achieved by random walks in the space of Slater determinants whose one‐particle orbitals are expressed in terms of the chosen basis set. A complication is that the determinants will in general acquire complex phases. This a consequence of ground state Monte Carlo projection using the Hubbard‐Stratonovich transformation of the two‐body interaction. To control the resulting “sign” decay, we describe a method we have recently introduced for the propagation of phaseless determinants. The approximation relies on importance sampling with a trial wave function. The approximation has features in common with diffusion MC fixed node and lattice‐model constrained path methods. Using a planewave basis and non‐local pseudopotentials, we apply the method to Si atom, dimer, and 2, 16, 54 atom (216 electrons) bulk supercells. Single Slater determinant wave functions from density functional theory calculations were used as |ΨT〉 with no additional optimization. The calculated binding energy of Si2 and cohesive energy of bulk Si are in excellent agreement with experiments and are comparable to the best existing theoretical results.

AIP conference proceedingsJan 1, 20030 citations

Linear-scaling quantum Monte Carlo calculations

Andrew Williamson, Randolph Q. Hood, Jeffrey C. Grossman

Abstract This study demonstrates linear scaling with system size in fixednode diffusion QMC calculations with pseudopotentials for carbon fullerenes and for hydrogenated silicon clusters with nearly 1000 valence electrons. The range covered is C20 to C18o and SiH4 to Si211H140, The calculations were carried out with standard pseudopotentials (leaving four electrons per C or Si atom) and Slater Jastrow trial functions with orbitals from LDA density functional calculations. The key elements leading to linear scaling were constructions of sparse determinants and simplification of orbital functions. The LDA orbitals were obtained with a plane-wave basis set which was transformed to yield Wannier functions with maximum localization. These functions were then fitted to cubic spline expressions with a limited number of grid points and cut-off distances matching those of the plane-wave orbitals. The computational effort for evaluating trial functions of this type was clearly seen to behave in a linear fashion. Exploration of the effects of variation of cut-off radius revealed that energies, sensitive only to node structures, increased rapidly for cut-offs less than 7 bohr but were unaffected for cutoffs greater than 7 bohr. The accuracies of the calculations were checked in other ways. These included comparison calculations for SiH3 and Si5H12 using Gaussian, original plane-wave, and Wannier function basis sets which gave close agreement for energies.

Jun 18, 20070 citations

A sparse algorithm for the evaluation of the local energy in quantum Monte Carlo

Alán Aspuru‐Guzik, Romelia Salomón–Ferrer, Brian M. Austin +1 more

A new algorithm is presented for the sparse representation and evaluation of Slater determinants in the quantum Monte Carlo (QMC) method. The approach, combined with the use of localized orbitals in a Slater-type orbital basis set, significantly extends the size molecule that can be treated with the QMC method. Application of the algorithm to systems containing up to 390 electrons confirms that the cost of evaluating the Slater determinant scales linearly with system size.

Journal of Computational ChemistryMar 10, 200527 citations

Electronic structure quantum Monte Carlo

Michal Bajdich, Luboš Mitáš

Electronic structure quantum Monte Carlo Quantum Monte Carlo (QMC) is an advanced simulation methodology for studies of manybody quantum systems. The QMC approaches combine analytical insights with stochastic computational techniques for efficient solution of several classes of important many-body problems such as the stationary Schrödinger equation. QMC methods of various flavors have been applied to a great variety of systems spanning continuous and lattice quantum models, molecular and condensed systems, BEC-BCS ultracold condensates, nuclei, etc. In this review, we focus on the electronic structure QMC, i.e., methods relevant for systems described by the electron-ion Hamiltonians. Some of the key QMC achievements include direct treatment of electron correlation, accuracy in predicting energy differences and favorable scaling in the system size. Calculations of atoms, molecules, clusters and solids have demonstrated QMC applicability to real systems with hundreds of electrons while providing 90-95% of the correlation energy and energy differences typically within a few percent of experiments. Advances in accuracy beyond these limits are hampered by the so-called fixed-node approximation which is used to circumvent the notorious fermion sign problem. Many-body nodes of fermion states and their properties have therefore become one of the important topics for further progress in predictive power and efficiency of QMC calculations. Some of our recent results on the wave function nodes and related nodal domain topologies will be briefly reviewed. This includes analysis of few-electron systems and descriptions of exact and approximate nodes using transformations and projections of the highly-dimensional nodal hypersurfaces into the 3D space. Studies of fermion nodes offer new insights into topological properties of eigenstates such as explicit demonstrations that generic fermionic ground states exhibit the minimal number of two nodal domains. Recently proposed trial wave functions based on Pfaffians with pairing orbitals are presented and their nodal properties are tested in calculations of first row atoms and molecules. Finally, backflow "dressed" coordinates are introduced as another possibility for capturing correlation effects and for decreasing the fixed-node bias.

arXiv:1008.2369Acta Physica Slovaca Reviews and TutorialsApr 1, 200933 citations

Quantum Monte Carlo calculations of energy gaps from first principles

R. J. Hunt, Marcin Szyniszewski, Genki I. Prayogo +2 more

We review the use of continuum quantum Monte Carlo (QMC) methods for the calculation of energy gaps from first principles, and present a broad set of excited-state calculations carried out with the variational and fixed-node diffusion QMC methods on atoms, molecules, and solids. We propose a finite-size-error correction scheme for bulk energy gaps calculated in finite cells subject to periodic boundary conditions. We show that finite-size effects are qualitatively different in two-dimensional materials, demonstrating the effect in a QMC calculation of the band gap and exciton binding energy of monolayer phosphorene. We investigate the fixed-node errors in diffusion Monte Carlo gaps evaluated with Slater-Jastrow trial wave functions by examining the effects of backflow transformations, and also by considering the formation of restricted multideterminant expansions for excited-state wave functions. For several molecules, we examine the importance of structural relaxation in the excited state in determining excited-state energies. We study the feasibility of using variational Monte Carlo with backflow correlations to obtain accurate excited-state energies at reduced computational cost, finding that this approach can be valid. We find that diffusion Monte Carlo gap calculations can be performed with much larger time steps than are typically required to converge the total energy, at significantly diminished computational expense, but that in order to alleviate fixed-node errors in calculations on solids the inclusion of backflow correlations is sometimes necessary.

arXiv:1806.04750Physical review. B./Physical review. BAug 13, 201850 citations

Quantum Monte Carlo study of correlated electronic systems

Li Chen

Describing correlated electron systems has been a major challenge in computational condensed-matter physics. Quantum Monte Carlo, a powerful computational tool for the study of correlated systems, solves electron correlation problems explicitly. It has been taken as a benchmark method for understanding the correlated systems. Instead of making approximations to Hamiltonian, QMC methods work with the wave functions, and the computational cost scales well with the system size. With the development of parallel computing, QMC calculations on large systems are becoming more and more feasible. We have investigated two correlated systems with highly accurate fixed node QMC techniques. The first system is a correlated hydrogen model system near the metal to insulator transition. We have successfully identified the transition point by calculating spin and charge properties and analyzing the low energy Hilbert space. The second one is a strongly correlated Fe/O system. Calculations on the Fe atoms, O atoms, and FeO molecules are conducted with multiple highly accurate many-body techniques. The source of errors has been disentangled by comparing the results of the many body techniques with the experimental results. For the Fe and O atoms, the calculated properties coincide well with previous experimental results. For the basis-based techniques, the performance is mainly limited by the basis set. The calculated equilibrium bond length, excitation energy and vibrational frequency of the FeO molecules are also in close agreement with the known values from previous experiments.

Nov 30, 20170 citations