Author
George H. Booth
Other affiliations: University of Cambridge, Princeton University, University College London
Bio: George H. Booth is an academic researcher from King's College London. The author has contributed to research in topics: Quantum Monte Carlo & Full configuration interaction. The author has an hindex of 41, co-authored 93 publications receiving 5337 citations. Previous affiliations of George H. Booth include University of Cambridge & Princeton University.
Papers
More filters
••
TL;DR: The capabilities and design philosophy of the current version of the PySCF package are document, which is as efficient as the best existing C or Fortran‐based quantum chemistry programs.
Abstract: Python-based simulations of chemistry framework (PySCF) is a general-purpose electronic structure platform designed from the ground up to emphasize code simplicity, so as to facilitate new method development and enable flexible computational workflows. The package provides a wide range of tools to support simulations of finite-size systems, extended systems with periodic boundary conditions, low-dimensional periodic systems, and custom Hamiltonians, using mean-field and post-mean-field methods with standard Gaussian basis functions. To ensure ease of extensibility, PySCF uses the Python language to implement almost all of its features, while computationally critical paths are implemented with heavily optimized C routines. Using this combined Python/C implementation, the package is as efficient as the best existing C or Fortran-based quantum chemistry programs. In this paper, we document the capabilities and design philosophy of the current version of the PySCF package. WIREs Comput Mol Sci 2018, 8:e1340. doi: 10.1002/wcms.1340
This article is categorized under:
Structure and Mechanism > Computational Materials Science
Electronic Structure Theory > Ab Initio Electronic Structure Methods
Software > Quantum Chemistry
1,042 citations
••
TL;DR: A new quantum Monte Carlo method for the simulation of correlated many-electron systems in full configuration-interaction (Slater determinant) spaces is developed, designed to simulate the underlying imaginary-time Schrödinger equation of the interacting Hamiltonian.
Abstract: We have developed a new quantum Monte Carlo method for the simulation of correlated many-electron systems in full configuration-interaction (Slater determinant) spaces. The new method is a population dynamics of a set of walkers, and is designed to simulate the underlying imaginary-time Schrodinger equation of the interacting Hamiltonian. The walkers (which carry a positive or negative sign) inhabit Slater determinant space, and evolve according to a simple set of rules which include spawning, death and annihilation processes. We show that this method is capable of converging onto the full configuration-interaction (FCI) energy and wave function of the problem, without any a priori information regarding the nodal structure of the wave function being provided. Walker annihilation is shown to play a key role. The pattern of walker growth exhibits a characteristic plateau once a critical (system-dependent) number of walkers has been reached. At this point, the correlation energy can be measured using two independent methods—a projection formula and a energy shift; agreement between these provides a strong measure of confidence in the accuracy of the computed correlation energies. We have verified the method by performing calculations on systems for which FCI calculations already exist. In addition, we report on a number of new systems, including CO, O2, CH4, and NaH—with FCI spaces ranging from 109 to 1014, whose FCI energies we compute using modest computational resources.
704 citations
••
TL;DR: The application of an exact technique, full configuration interaction quantum Monte Carlo to a variety of real solids, providing reference many-electron energies that are used to rigorously benchmark the standard hierarchy of quantum-chemical techniques, up to the ‘gold standard’ coupled-cluster ansatz.
Abstract: The properties of all materials arise largely from the quantum mechanics of their constituent electrons under the influence of the electric field of the nuclei. The solution of the underlying many-electron Schrodinger equation is a ‘non-polynomial hard’ problem, owing to the complex interplay of kinetic energy, electron–electron repulsion and the Pauli exclusion principle. The dominant computational method for describing such systems has been density functional theory. Quantum-chemical methods—based on an explicit ansatz for the many-electron wavefunctions and, hence, potentially more accurate—have not been fully explored in the solid state owing to their computational complexity, which ranges from strongly exponential to high-order polynomial in system size. Here we report the application of an exact technique, full configuration interaction quantum Monte Carlo to a variety of real solids, providing reference many-electron energies that are used to rigorously benchmark the standard hierarchy of quantum-chemical techniques, up to the ‘gold standard’ coupled-cluster ansatz, including single, double and perturbative triple particle–hole excitation operators. We show the errors in cohesive energies predicted by this method to be small, indicating the potential of this computationally polynomial scaling technique to tackle current solid-state problems. Recent developments that reduce the computational cost and scaling of wavefunction-based quantum-chemical techniques open the way to the successful application of such techniques to a variety of real-world solids. Computational descriptions of solid-state materials are currently dominated by methods based on density functional theory. An attractive and potentially more accurate approach would be to adopt the wavefunction-based methods of quantum chemistry, although these have not received as much attention because of the computational complexities involved. Now George Booth and colleagues show how recent developments that serve to reduce the computational cost and scaling of such quantum-chemical techniques open the way to their successful application to a variety of real-world solids.
537 citations
••
TL;DR: It is shown that rather small initiator spaces and numbers of walkers can converge with submilliHartree accuracy to the known full configuration-interaction (FCI) energy (in the cc-pVDZ basis), in both the equilibrium geometry and the multiconfigurational stretched case.
Abstract: We provide a very simple adaptation of our recently published quantum Monte Carlo algorithm in full configuration-interaction (Slater determinant) spaces which dramatically reduces the number of walkers required to achieve convergence. A survival criterion is imposed for newly spawned walkers. We define a set of initiator determinants such that progeny of walkers spawned from such determinants onto unoccupied determinants are able to survive, while the progeny of walkers not in this set can survive only if they are spawned onto determinants which are already occupied. The set of initiators is originally defined to be all determinants constructable from a subset of orbitals, in analogy with complete-active spaces. This set is dynamically updated so that if a noninitiator determinant reaches an occupation larger than a preset limit, it becomes an initiator. The new algorithm allows sign-coherent sampling of the FCI space to be achieved with relatively few walkers. Using the N2 molecule as an illustration, w...
377 citations
••
California Institute of Technology1, Ohio State University2, Chinese Academy of Sciences3, University of Cambridge4, Max Planck Society5, King's College London6, University of Florida7, University of Bristol8, Google9, Free University of Berlin10, University of Minnesota11, University of Notre Dame12, University of Helsinki13, Beijing Normal University14, Johns Hopkins University15, Amgen16, IBM17, University of Colorado Boulder18, Delft University of Technology19, Old Dominion University20, Columbia University21, Yale University22, Interdisciplinary Center for Scientific Computing23, University of Pittsburgh24, University of Illinois at Urbana–Champaign25, Dartmouth College26, University of Hong Kong27, University of California, Irvine28
TL;DR: PySCF as mentioned in this paper is a Python-based general-purpose electronic structure platform that supports first-principles simulations of molecules and solids as well as accelerates the development of new methodology and complex computational workflows.
Abstract: PySCF is a Python-based general-purpose electronic structure platform that supports first-principles simulations of molecules and solids as well as accelerates the development of new methodology and complex computational workflows. This paper explains the design and philosophy behind PySCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PySCF as a development environment. We then summarize the capabilities of PySCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PySCF across the domains of quantum chemistry, materials science, machine learning, and quantum information science.
374 citations
Cited by
More filters
•
28,685 citations
••
TL;DR: A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided in this paper, covering approximately the last seven years, including developments in density functional theory and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces.
Abstract: A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Moller–Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr_2 dimer, exploring zeolite-catalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.
2,396 citations
01 Jan 2015
TL;DR: Detailed benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller–Plesset methods for intermolecular interactions, and tests of the accuracy of implicit solvation models are provided.
Abstract: A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller–Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr2 dimer, exploring zeolite-catalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.
1,919 citations
••
1,908 citations
01 Feb 1995
TL;DR: In this paper, the unpolarized absorption and circular dichroism spectra of the fundamental vibrational transitions of the chiral molecule, 4-methyl-2-oxetanone, are calculated ab initio using DFT, MP2, and SCF methodologies and a 5S4P2D/3S2P (TZ2P) basis set.
Abstract: : The unpolarized absorption and circular dichroism spectra of the fundamental vibrational transitions of the chiral molecule, 4-methyl-2-oxetanone, are calculated ab initio. Harmonic force fields are obtained using Density Functional Theory (DFT), MP2, and SCF methodologies and a 5S4P2D/3S2P (TZ2P) basis set. DFT calculations use the Local Spin Density Approximation (LSDA), BLYP, and Becke3LYP (B3LYP) density functionals. Mid-IR spectra predicted using LSDA, BLYP, and B3LYP force fields are of significantly different quality, the B3LYP force field yielding spectra in clearly superior, and overall excellent, agreement with experiment. The MP2 force field yields spectra in slightly worse agreement with experiment than the B3LYP force field. The SCF force field yields spectra in poor agreement with experiment.The basis set dependence of B3LYP force fields is also explored: the 6-31G* and TZ2P basis sets give very similar results while the 3-21G basis set yields spectra in substantially worse agreements with experiment. jg
1,652 citations