- Open Access
Towards the Solution of the Many-Electron Problem in Real Materials: Equation of State of the Hydrogen Chain with State-of-the-Art Many-Body Methods
Phys. Rev. X 7, 031059 – Published 28 September, 2017Erratum Phys. Rev. X 11, 029901 (2021)
DOI: https://doi.org/10.1103/PhysRevX.7.031059
Abstract
We present numerical results for the equation of state of an infinite chain of hydrogen atoms. A variety of modern many-body methods are employed, with exhaustive cross-checks and validation. Approaches for reaching the continuous space limit and the thermodynamic limit are investigated, proposed, and tested. The detailed comparisons provide a benchmark for assessing the current state of the art in many-body computation, and for the development of new methods. The ground-state energy per atom in the linear chain is accurately determined versus bond length, with a confidence bound given on all uncertainties.
Physics Subject Headings (PhySH)
Corrections
24 March, 2021
Erratum
Popular Summary
One of the grand challenges of modern science is to understand and predict the quantum-mechanical behavior of a large ensemble of interacting electrons. Physical and chemical properties of materials and molecules are often the result of a delicate balance between competing behaviors in a many-electron system. Accurate computations are essential for predicting outcomes, but performing such calculations in a straightforward manner requires unattainable computational expense. Researchers have therefore come up with a variety of theoretical and numerical techniques to work around this hurdle. We present a comprehensive benchmark study of quantum many-body computational methods for addressing this challenge.
Our work focuses on two key aspects in treating real materials: the presence of long-ranged Coulomb interactions (the electrostatic force felt among electrons) and the need to study the continuum and thermodynamic limits. We characterize the relative accuracy and capabilities of 16 methods for performing many-electron calculations, which provides a survey of the state-of-the-art to guide applications. A large amount of data is produced that will be useful in benchmarking other existing and future electronic structure methods. Combining the strengths of complementary methods, we determine the equation of state (energy per atom versus interatomic spacing) of an infinite linear chain of hydrogen atoms in the continuum limit to high accuracy.
Future work can extend these benchmarks to more complex materials. Any progress in addressing the computational challenge of many-electron systems will be fundamental to “materials genome” initiatives, which attempt to discover and design new materials for a range of applications.
Article Text
References (108)
- P. A. M. Dirac, Quantum Mechanics of Many-Electron Systems, Proc. R. Soc. A 123, 714 (1929).
- J. P. F. LeBlanc et al., Solutions of the Two-Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms, Phys. Rev. X 5, 041041 (2015).
- B.-X. Zheng, C.-M. Chung, P. Corboz, G. Ehlere, M.-P. Qin, R. M. Noack, H. Shi, S. R. White, S. Zhang, and G. K.-L. Chan, Stripe Order in the Underdoped Region of the Two-Dimensional Hubbard Model, arXiv:1701.00054.
- J. Hachmann, W. Cardoen, and G. K.-L. Chan, Multireference Correlation in Long Molecules with the Quadratic Scaling Density Matrix Renormalization Group, J. Chem. Phys. 125, 144101 (2006).
- W. A. Al-Saidi, S. Zhang, and H. Krakauer, Bond Breaking with Auxiliary-Field Quantum Monte Carlo, J. Chem. Phys. 127, 144101 (2007).
- T. Tsuchimochi and G. E. Scuseria, Strong Correlations via Constrained-Pairing Mean-Field Theory, J. Chem. Phys. 131, 121102 (2009).
- A. V. Sinitskiy, L. Greenman, and D. A. Mazziotti, Strong Correlation in Hydrogen Chains and Lattices Using the Variational Two-Electron Reduced Density Matrix Method, J. Chem. Phys. 133, 014104 (2010).
- D. A. Mazziotti, Large-Scale Semidefinite Programming for Many-Electron Quantum Mechanics, Phys. Rev. Lett. 106, 083001 (2011).
- N. Lin, C. A. Marianetti, A. J. Millis, and D. R. Reichman, Dynamical Mean-Field Theory for Quantum Chemistry, Phys. Rev. Lett. 106, 096402 (2011).
- L. Stella, C. Attaccalite, S. Sorella, and A. Rubio, Strong Electronic Correlation in the Hydrogen Chain: A Variational Monte Carlo Study, Phys. Rev. B 84, 245117 (2011).
- M. Motta et al., Github repository, https://github.com/simonsfoundation/hydrogen-benchmark-PRX/.
- W. Kutzelnigg, R12-Dependent Terms in the Wave Function as Closed Sums of Partial Wave Amplitudes for Large l, Theor. Chim. Acta 68, 445 (1985).
- L. Kong and E. F. Valeev, Perturbative Correction for the Basis Set Incompleteness Error of Complete-Active-Space Self-Consistent Field, J. Chem. Phys. 133, 174126 (2010).
- R. Blankenbecler, D. J. Scalapino, and R. L. Sugar, Monte Carlo Calculations of Coupled Boson-Fermion Systems. I, Phys. Rev. D 24, 2278 (1981).
- S. Zhang, J. Carlson, and J. E. Gubernatis, Constrained Path Monte Carlo Method for Fermion Ground States, Phys. Rev. B 55, 7464 (1997).
- S. Zhang and H. Krakauer, Quantum Monte Carlo Method Using Phase-Free Random Walks with Slater Determinants, Phys. Rev. Lett. 90, 136401 (2003).
- K. Van Houcke, E. Kozik, N. Prokof’ev, and B. Svistunov, Diagrammatic Monte Carlo, Phys. Procedia 6, 95 (2010).
- S. A. Kulagin, N. Prokof’ev, O. A. Starykh, B. Svistunov, and C. N. Varney, Bold Diagrammatic Monte Carlo Method Applied to Fermionized Frustrated Spins, Phys. Rev. Lett. 110, 070601 (2013).
- I. S. Tupitsyn, A. S. Mishchenko, N. Nagaosa, and N. Prokof’ev, Coulomb and Electron-Phonon Interactions in Metals, Phys. Rev. B 94, 155145 (2016).
- G. Knizia and Garnet Kin-Lic Chan, Density Matrix Embedding: A Simple Alternative to Dynamical Mean-Field Theory, Phys. Rev. Lett. 109, 186404 (2012).
- G. Knizia and G. K. -L. Chan, Density Matrix Embedding: A Strong-Coupling Quantum Embedding Theory, J. Chem. Theory Comput. 9, 1428 (2013).
- S. Wouters, C. A. Jimenez-Hoyos, Q. Sun, and G. K. -L. Chan, A Practical Guide to Density Matrix Embedding Theory in Quantum Chemistry, J. Chem. Theory Comput. 12, 2706 (2016).
- S. R. White, Density Matrix Formulation for Quantum Renormalization Groups, Phys. Rev. Lett. 69, 2863 (1992).
- S. R. White and R. L. Martin, Ab Initio Quantum Chemistry Using the Density Matrix Renormalization Group, J. Chem. Phys. 110, 4127 (1999).
- G. K. -L. Chan and M. Head-Gordon, Highly Correlated Calculations with a Polynomial Cost Algorithm: A Study of the Density Matrix Renormalization Group, J. Chem. Phys. 116, 4462 (2002).
- G. K. -L. Chan, A. Keselman, N. Nakatani, Z. Li, and S. R. White, Matrix Product Operators, Matrix Product States, and Ab Initio Density Matrix Renormalization Group Algorithms, J. Chem. Phys. 145, 014102 (2016).
- J. J. Phillips and D. Zgid, Communication: The Description of Strong Correlation within Self-Consistent Green’s Function Second-Order Perturbation Theory, J. Chem. Phys. 140, 241101 (2014).
- A. A. Rusakov and D. Zgid, Self-Consistent Second-Order Green’s Function Perturbation Theory for Periodic Systems, J. Chem. Phys. 144, 054106 (2016).
- J. J. Phillips, A. A. Kananenka, and D. Zgid, Fractional Charge and Spin Errors in Self-Consistent Green’s Function Theory, J. Chem. Phys. 142, 194108 (2015).
- A. R. Welden, A. A. Rusakov, and D. Zgid, Exploring Connections between Statistical Mechanics and Green’s Functions for Realistic Systems. Temperature Dependent Electronic Entropy and Internal Energy from a Self-Consistent Second-Order Green’s Function, J. Chem. Phys. 145, 204106 (2016).
- N. E. Dahlen and U. von Barth, Variational Energy Functionals Tested on Atoms, Phys. Rev. B 69, 195102 (2004).
- N. E. Dahlen and U. von Barth, Variational Second-Order Møller-Plesset Theory Based on the Luttinger-Ward Functional, J. Chem. Phys. 120, 6826 (2004).
- M. Casula, C. Filippi, and S. Sorella, Diffusion Monte Carlo Method with Lattice Regularization, Phys. Rev. Lett. 95, 100201 (2005).
- P. Reynolds, D. Ceperley, B. Alder, and J. Lester, Fixed-Node Quantum Monte Carlo for Molecules, J. Chem. Phys. 77, 5593 (1982).
- H.-J. Werner and P. J. Knowles, An Efficient Internally Contracted Multiconfiguration-Reference Configuration Interaction Method, J. Chem. Phys. 89, 5803 (1988).
- P. J. Knowles and H.-J. Werner, An Efficient Method for the Evaluation of Coupling Coefficients in Configuration Interaction Calculations, Chem. Phys. Lett. 145, 514 (1988).
- C. Angeli, M. Pastore, and R. Cimiraglia, New Perspectives in Multireference Perturbation Theory: The N-Electron Valence State Approach, Theor. Chem. Acc. 117, 743 (2007).
- J. Paldus and X. Li, A Critical Assessment of Coupled Cluster Method in Quantum Chemistry, Adv. Chem. Phys. 110, 1 (1999).
- R. J. Bartlett and M. Musial, Coupled-Cluster Theory in Quantum Chemistry, Rev. Mod. Phys. 79, 291 (2007).
- I. Shavitt and R. J. Bartlett, Many-Body Methods in Chemistry and Physics (Cambridge University Press, New York, 2009).
- C. C. J. Roothaan, Self-Consistent Field Theory for Open Shells of Electronic Systems, Rev. Mod. Phys. 32, 179 (1960).
- E. Stoudenmire and S. R. White, Sliced Basis Density Matrix Renormalization Group for Electronic Structure, Phys. Rev. Lett. 119, 046401 (2017).
- C.-O. Almbladh, U. V. Barth, and R. V. Leeuwen, Variational Total Energies from - and -Derivable Theories, Int. J. Mod. Phys. B 13, 535 (1999).
- A. Stan, N. E. Dahlen, and R. van Leeuwen, Fully Self-Consistent GW Calculations for Atoms and Molecules, Europhys. Lett. 76, 298 (2006).
- A. Stan, N. E. Dahlen, and R. van Leeuwen, Levels of Self-Consistency in the GW Approximation, J. Chem. Phys. 130, 114105 (2009).
- A. A. Kananenka, E. Gull, and D. Zgid, Systematically Improvable Multiscale Solver for Correlated Electron Systems, Phys. Rev. B 91, 121111 (2015).
- T. N. Lan, A. A. Kananenka, and D. Zgid, Communication: Towards Ab Initio Self-Energy Embedding Theory in Quantum Chemistry, J. Chem. Phys. 143, 241102 (2015).
- T. Nguyen Lan, A. A. Kananenka, and D. Zgid, Rigorous Ab Initio Quantum Embedding for Quantum Chemistry Using Green’s Function Theory: Screened Interaction, Nonlocal Self-Energy Relaxation, Orbital Basis, and Chemical Accuracy, J. Chem. Theory Comput. 12, 4856 (2016).
- D. Zgid and E. Gull, Finite Temperature Quantum Embedding Theories for Correlated Systems, New J. Phys. 19, 023047 (2017).
- T. N. Lan and D. Zgid, Generalized Self-Energy Embedding Theory, arXiv:1703.06981.
- J. A. Pople and R. K. Nesbet, Self-Consistent Orbitals for Radicals, J. Chem. Phys. 22, 571 (1954).
- W. L. McMillan, Ground State of Liquid , Phys. Rev. 138, A442 (1965).
- D. Ceperley, G. V. Chester, and M. H. Kalos, Monte Carlo Simulation of a Many-Fermion Study, Phys. Rev. B 16, 3081 (1977).
- B.-X. Zheng, J. S. Kretchmer, H. Shi, S. Zhang, and Garnet Kin-Lic Chan, Cluster Size Convergence of the Density Matrix Embedding Theory and Its Dynamical Cluster Formulation: A Study with an Auxiliary-Field Quantum Monte Carlo Solver, Phys. Rev. B 95, 045103 (2017).
- T. H. Dunning, Jr., Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. The Atoms Boron through Neon and Hydrogen, J. Chem. Phys. 90, 1007 (1989).
- D. Feller, Application of Systematic Sequences of Wave Functions to the Water Dimer, J. Chem. Phys. 96, 6104 (1992).
- T. Helgaker, W. Klopper, H. Koch, and J. Noga, Basis-Set Convergence of Correlated Calculations on Water, J. Chem. Phys. 106, 9639 (1997).
- E. Stoudenmire and S. R. White, Studying Two-Dimensional Systems with the Density Matrix Renormalization Group, Annu. Rev. Condens. Matter Phys. 3, 111 (2012).
- S. Chiesa, D. M. Ceperley, R. M. Martin, and M. Holzmann, Finite-Size Error in Many-Body Simulations with Long-Range Interactions, Phys. Rev. Lett. 97, 076404 (2006).
- G. D. Purvis III and R. J. Bartlett, A Full Coupled-Cluster Singles and Doubles Model: The Inclusion of Disconnected Triples, J. Chem. Phys. 76, 1910 (1982).
- G. E. Scuseria, A. C. Scheiner, T. J. Lee, J. E. Rice, and H. F. Schaefer III, The Closed-Shell Coupled Cluster Single and Double Excitation (CCSD) Model for the Description of Electron Correlation. A Comparison with Configuration Interaction (CISD) Results, J. Chem. Phys. 86, 2881 (1987).
- J. A. Pople, M. Head-Gordon, and K. Rghavachari, Quadratic Configuration Interaction. A General Technique for Determining Electron Correlation Energies, J. Chem. Phys. 87, 5968 (1987).
- M. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria et al., Gaussian 09, Revision D.01, Gaussian, Inc., 2013.
- H. Koch, O. Christiansen, R. Kobayashi, P. Jorgensen, and T. Helgaker, A Direct Atomic Orbital Driven Implementation of the Coupled Cluster Singles and Doubles (CCSD) Model, Chem. Phys. Lett. 228, 233 (1994).
- S. Wouters and D. Van Neck, The Density Matrix Renormalization Group for Ab Initio Quantum Chemistry, Eur. Phys. J. D 68, 272 (2014).
- R. Olivares-Amaya, W. Hu, N. Nakatani, S. Sharma, J. Yang, and G. K.-L. Chan, The Ab-Initio Density Matrix Renormalization Group in Practice, J. Chem. Phys. 142, 034102 (2015).
- G. Barcza, Ö. Legeza, K. H. Marti, and M. Reiher, Quantum-Information Analysis of Electronic States of Different Molecular Structures, Phys. Rev. A 83, 012508 (2011).
- H.-J. Werner, P. J. Knowles, G. Knizia, F. R. Manby, and M. Schütz, MOLPRO: A General-Purpose Quantum Chemistry Program Package, Comput. Mol. Sci. 2, 242 (2012).
- T. Shiozaki, G. Knizia, and H.-J. Werner, Explicitly Correlated Multireference Configuration Interaction: MRCI-F12, J. Chem. Phys. 134, 034113 (2011).
- T. Shiozaki and H.-J. Werner, Multireference Explicitly Correlated F12 Theories, Mol. Phys. 111, 607 (2013).
- T. Shiozaki and H.-J. Werner, Explicitly Correlated Multireference Configuration Interaction with Multiple Reference Functions: Avoided Crossings and Conical Intersections, J. Chem. Phys. 134, 184104 (2011).
- K. G. Dyall, The Choice of a Zeroth-Order Hamiltonian for Second-Order Perturbation Theory with a Complete Active Space Self-Consistent-Field Reference Function, J. Chem. Phys. 102, 4909 (1995).
- S. Guo, M. A. Watson, W. Hu, Q. Sun, and G. K. -L. Chan, N-Electron Valence State Perturbation Theory Based on a Density Matrix Renormalization Group Reference Function, with Applications to the Chromium Dimer and a Trimer Model of Poly (p-Phenylenevinylene), J. Chem. Theory Comput. 12, no. 1583 (2016).
- S. Sharma and G. Chan, Communication: A Flexible Multi-Reference Perturbation Theory by Minimizing the Hylleraas Functional with Matrix Product States, J. Chem. Phys. 141, 111101 (2014).
- S. Sharma, G. Knizia, S. Guo, and A. Alavi, Combining Internally Contracted States and Matrix Product States to Perform Multireference Perturbation Theory, J. Chem. Theory Comput. 13, 488 (2017).
- R. L. Stratonovich, On a Method of Calculating Quantum Distribution Functions, Sov. Phys. Dokl. 2, 416 (1957).
- J. Hubbard, Calculation of Partition Functions, Phys. Rev. Lett. 3, 77 (1959).
- S. Zhang, Auxiliary-Field Quantum Monte Carlo for Correlated Electron Systems, in Emergent Phenomena in Correlated Matter, edited by E. Pavarini,, E. Koch, and U. Schollwock (Forschungszentrum Jülich, Julich, Germany, 2013), Vol. 3.
- M. Qin, H. Shi, and S. Zhang, Coupling Quantum Monte Carlo and Independent-Particle Calculations: Self-Consistent Constraint for the Sign Problem Based on the Density or the Density Matrix, Phys. Rev. B 94, 235119 (2016).
- W. Purwanto, H. Krakauer, and S. Zhang, Pressure-Induced Diamond to -Tin Transition in Bulk Silicon: A Quantum Monte Carlo Study, Phys. Rev. B 80, 214116 (2009).
- Y. Virgus, W. Purwanto, H. Krakauer, and S. Zhang, Stability, Energetics, and Magnetic States of Cobalt Adatoms on Graphene, Phys. Rev. Lett. 113, 175502 (2014).
- W. Purwanto, S. Zhang, and H. Krakauer, Excited State Calculations Using Phaseless Auxiliary-Field Quantum Monte Carlo: Potential Energy Curves of Low-Lying Singlet States, J. Chem. Phys. 130, 094107 (2009).
- F. Ma, S. Zhang, and H. Krakauer, Auxiliary-Field Quantum Monte Carlo Calculations with Multiple-Projector Pseudopotentials, Phys. Rev. B 95, 165103 (2017).
- W. A. Al-Saidi, S. Zhang, and H. Krakauer, Auxiliary-Field Quantum Monte Carlo Calculations of Molecular Systems with a Gaussian Basis, J. Chem. Phys. 124, 224101 (2006).
- F. Ma, W. Purwanto, S. Zhang, and H. Krakauer, Quantum Monte Carlo Calculations in Solids with Downfolded Hamiltonians, Phys. Rev. Lett. 114, 226401 (2015).
- W. Purwanto, H. Krakauer, Y. Virgus, and S. Zhang, Assessing Weak Hydrogen Binding on Centers: An Accurate Many-Body Study with Large Basis Sets, J. Chem. Phys. 135, 164105 (2011).
- M. Casula, S. Moroni, S. Sorella, and C. Filippi, Size-Consistent Variational Approaches to Nonlocal Pseudopotentials: Standard and Lattice Regularized Diffusion Monte Carlo Methods Revisited, J. Chem. Phys. 132, 154113 (2010).
- S. Sorella, TurboRVB Quantum Monte Carlo Package, http://people.sissa.it/~sorella/web/.
- R. Dovesi, R. Orlando, A. Erba, C. M. Zicovich-Wilson, B. Civalleri, S. Casassa, L. Maschio, M. Ferrabone, M. D. L. Pierre, P. D’Arco, Y. Noel, M. Causa, M. Rerat, and B. Kirtman, CRYSTAL14: A Program for the Ab Initio Investigation of Crystalline Solids, Int. J. Quantum Chem. 114, 1287 (2014).
- A. Zen, Y. Luo, S. Sorella, and L. Guidoni, Molecular Properties by Quantum Monte Carlo: An Investigation on the Role of the Wave Function Ansatz and the Basis Set in the Water Molecule, J. Chem. Theory Comput. 9, 4332 (2013).
- M. Marchi, S. Azadi, M. Casula, and S. Sorella, Resonating Valence Bond Wave Function with Molecular Orbitals: Application to First-Row Molecules, J. Chem. Phys. 131, 154116 (2009).
- G. Knizia, Intrinsic Atomic Orbitals: An Unbiased Bridge between Quantum Theory and Chemical Concepts, J. Chem. Theory Comput. 9, 4834 (2013).
- S. Wouters, C. A. Jiménez-Hoyos, Q. Sun, and G. K.-L. Chan, A Practical Guide to Density Matrix Embedding Theory in Quantum Chemistry, J. Chem. Theory Comput. 12, 2706 (2016).
- L. Hedin, New Method for Calculating the One-Particle Green’s Function with Application to the Electron-Gas Problem, Phys. Rev. 139, A796 (1965).
- A. A. Rusakov, J. J. Phillips, and D. Zgid, Local Hamiltonians for Quantitative Green’s Function Embedding Methods, J. Chem. Phys. 141, 194105 (2014).
- A. A. Kananenka, A. R. Welden, T. N. Lan, E. Gull, and D. Zgid, Efficient Temperature-Dependent Green’s Function Methods for Realistic Systems: Using Cubic Spline Interpolation to Approximate Matsubara Green’s Functions, J. Chem. Theory Comput. 12, 2250 (2016).
- L. Hedin, New Method for Calculating the One-Particle Green’s Function with Application to the Electron-Gas Problem, Phys. Rev. 139, A796 (1965).
- J. M. Luttinger and J. C. Ward, Ground-State Energy of a Many-Fermion System. II, Phys. Rev. 118, 1417 (1960).
- G. Baym, Self-Consistent Approximations in Many-Body Systems, Phys. Rev. 127, 1391 (1962).
- A. A. Kananenka, J. J. Phillips, and D. Zgid, Efficient Temperature-Dependent Green’s Functions Methods for Realistic Systems: Compact Grids for Orthogonal Polynomial Transforms, J. Chem. Theory Comput. 12, 564 (2016).
- F. Caruso, P. Rinke, X. Ren, A. Rubio, and M. Scheffler, Self-Consistent GW: All-Electron Implementation with Localized Basis Functions, Phys. Rev. B 88, 075105 (2013).
- A. Gaenko, A. E. Antipov, G. Carcassi, T. Chen, X. Chen, Q. Dong, L. Gamper, J. Gukelberger, R. Igarashi, S. Iskakov, M. Könz, J. P. F. LeBlanc, R. Levy, P. N. Ma, J. E. Paki, H. Shinaoka, S. Todo, M. Troyer, and E. Gull, Updated Core Libraries of the ALPS Project, arXiv:1609.03930.
- E. F. Valeev, LIBINT: A Library for the Evaluation of Molecular Integrals of Many-Body Operators over Gaussian Functions (version 2.1.0), http://libint.valeyev.net/.
- N. Prokof’ev and B. Svistunov, Bold Diagrammatic Monte Carlo Technique: When the Sign Problem is Welcome, Phys. Rev. Lett. 99, 250201 (2007).
- R. Rossi, Determinant Diagrammatic Monte Carlo in the Thermodynamic Limit, Phys. Rev. Lett. 119, 045701 (2017).
- R. Rossi, N. V. Prokof’ev, B. Svistunov, K. Van Houcke, and F. Werner, Polynomial Complexity Despite the Fermionic Sign, Europhys. Lett. 118, 10004 (2017).
- I. S. Tupitsyn and N. V. Prokof’ev, Stability of Dirac Liquids with Strong Coulomb Interaction, Phys. Rev. Lett. 118, 026403 (2017).
For and we perform fragment-size extrapolations on , which is not converged with respect to the full system size. Our extrapolated result includes a full-system correction determined from calculations on larger systems, with a fragment of size [22].
