1
|
Kessler J, Calcavecchia F, Kühne TD. Artificial Neural Networks as Trial Wave Functions for Quantum Monte Carlo. ADVANCED THEORY AND SIMULATIONS 2021. [DOI: 10.1002/adts.202000269] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Jan Kessler
- Dynamics of Condensed Matter and Center for Sustainable Systems Design Chair of Theoretical Chemistry University of Paderborn Warburger Str. 100 D‐33098 Paderborn Germany
| | - Francesco Calcavecchia
- Dynamics of Condensed Matter and Center for Sustainable Systems Design Chair of Theoretical Chemistry University of Paderborn Warburger Str. 100 D‐33098 Paderborn Germany
| | - Thomas D. Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design Chair of Theoretical Chemistry University of Paderborn Warburger Str. 100 D‐33098 Paderborn Germany
- Paderborn Center for Parallel Computing and Institute for Lightweight Design University of Paderborn Warburger Str. 100 D‐33098 Paderborn Germany
| |
Collapse
|
2
|
Azadi S, Booth GH, Kühne TD. Equation of state of atomic solid hydrogen by stochastic many-body wave function methods. J Chem Phys 2020; 153:204107. [DOI: 10.1063/5.0026499] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sam Azadi
- Department of Physics, King’s College London, Strand, WC2R 2LS London, United Kingdom
| | - George H. Booth
- Department of Physics, King’s College London, Strand, WC2R 2LS London, United Kingdom
| | - Thomas D. Kühne
- Department of Chemistry, Paderborn Center for Parallel Computing, Paderborn University, 33098 Paderborn, Germany
| |
Collapse
|
3
|
Sahoo SK, Heske J, Azadi S, Zhang Z, Tarakina NV, Oschatz M, Khaliullin RZ, Antonietti M, Kühne TD. On the Possibility of Helium Adsorption in Nitrogen Doped Graphitic Materials. Sci Rep 2020; 10:5832. [PMID: 32242048 PMCID: PMC7118168 DOI: 10.1038/s41598-020-62638-z] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2020] [Accepted: 03/17/2020] [Indexed: 11/29/2022] Open
Abstract
The potassium salt of polyheptazine imide (K-PHI) is a promising photocatalyst for various chemical reactions. From powder X-ray diffraction data an idealized structural model of K-PHI has been derived. Using atomic coordinates of this model we defined an energetically optimized K-PHI structure, in which the K ions are present in the pore and between the PHI-planes. The distance between the anion framework and K+ resembles a frustrated Lewis pair-like structure, which we denote as frustrated Coulomb pair that results in an interesting adsorption environment for otherwise non-adsorbing, non-polar gas molecules. We demonstrate that even helium (He) gas molecules, which are known to have the lowest boiling point and the lowest intermolecular interactions, can be adsorbed in this polarized environment with an adsorption energy of - 4.6 kJ mol-1 per He atom. The interaction between He atoms and K-PHI is partially originating from charge transfer, as disclosed by our energy decomposition analysis based on absolutely localized molecular orbitals. Due to very small charge transfer interactions, He gas adsorption saturates at 8 at%, which however can be subject to further improvement by cation variation.
Collapse
Affiliation(s)
- Sudhir K Sahoo
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Warburger Str. 100, D-33098, Paderborn, Germany
| | - Julian Heske
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Warburger Str. 100, D-33098, Paderborn, Germany
- Department of Colloid Chemistry, Max Planck Institute of Colloids and Interfaces, Am Mühlenberg 1, D-14476, Potsdam, Germany
| | - Sam Azadi
- Department of Physics, King's College London, Strand, London, WC2R 2L, United Kingdom
- Department of Physics, Imperial College London, Exhibition Road, London, SW7 2AZ, United Kingdom
| | - Zhenzhe Zhang
- Department of Chemistry, McGill University, 801 Sherbrooke Str. West, Montreal, Quebec, H3A 0B8, Canada
| | - Nadezda V Tarakina
- Department of Colloid Chemistry, Max Planck Institute of Colloids and Interfaces, Am Mühlenberg 1, D-14476, Potsdam, Germany
| | - Martin Oschatz
- Department of Colloid Chemistry, Max Planck Institute of Colloids and Interfaces, Am Mühlenberg 1, D-14476, Potsdam, Germany
- University of Potsdam, Institute of Chemistry, Karl-Liebknecht-Str. 24-25, D-14476, Potsdam, Germany
| | - Rustam Z Khaliullin
- Department of Chemistry, McGill University, 801 Sherbrooke Str. West, Montreal, Quebec, H3A 0B8, Canada
| | - Markus Antonietti
- Department of Colloid Chemistry, Max Planck Institute of Colloids and Interfaces, Am Mühlenberg 1, D-14476, Potsdam, Germany
| | - Thomas D Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Warburger Str. 100, D-33098, Paderborn, Germany.
- Paderborn Center for Parallel Computing and Institute for Lightweight Design, University of Paderborn, Warburger Str. 100, D-33098, Paderborn, Germany.
| |
Collapse
|
4
|
Chin SA. Solving fermion problems without solving the sign problem: Symmetry-breaking wave functions from similarity-transformed propagators for solving two-dimensional quantum dots. Phys Rev E 2020; 101:043304. [PMID: 32422780 DOI: 10.1103/physreve.101.043304] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2019] [Accepted: 03/17/2020] [Indexed: 06/11/2023]
Abstract
It is well known that the use of the primitive second-order propagator in path-integral Monte Carlo calculations of many-fermion systems leads to the sign problem. This work will show that by using the similarity-transformed Fokker-Planck propagator, it is possible to solve for the ground state of a large quantum dot, with up to 100 polarized electrons, without solving the sign problem. These similarity-transformed propagators naturally produce rotational symmetry-breaking ground-state wave functions previously used in the study of quantum dots and quantum Hall effects. However, instead of localizing the electrons at positions that minimize the potential energy, this derivation shows that they should be located at positions that maximize the bosonic ground-state wave function. Further improvements in the energy can be obtained by using these as initial wave functions in a ground-state path-integral Monte Carlo calculation with second- and fourth-order propagators.
Collapse
Affiliation(s)
- Siu A Chin
- Department of Physics and Astronomy, Texas A&M University, College Station, Texas 77843, USA
| |
Collapse
|
5
|
Calcavecchia F, Kühne TD. Metal-Insulator Transition of Solid Hydrogen by the Antisymmetric Shadow Wave Function. ACTA ACUST UNITED AC 2018. [DOI: 10.1515/zna-2018-0180] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Abstract
We revisit the pressure-induced molecular-atomic metal-insulator transition of solid hydrogen by means of variational quantum Monte Carlo simulations based on the antisymmetric shadow wave function. For the purpose of facilitating the study of the electronic structure of large-scale fermionic systems, the shadow wave function formalism is extended by a series of technical advancements as implemented in our HswfQMC code. Among others, these improvements include a revised optimization method for the employed shadow wave function and an enhanced treatment of periodic systems with long-range interactions. It is found that the superior accuracy of the antisymmetric shadow wave function results in a significantly increased transition pressure with respect to previous theoretical estimates.
Collapse
Affiliation(s)
- Francesco Calcavecchia
- LPMMC, UMR 5493 of CNRS , Université Grenoble Alpes , 38042 Grenoble , France
- Institute of Physics , Johannes Gutenberg University , Staudingerweg 7 , D-55128 Mainz , Germany
- Graduate School of Excellence Materials Science in Mainz , Staudingerweg 9 , D-55128 Mainz , Germany
| | - Thomas D. Kühne
- Dynamics of Condensed Matter, Department of Chemistry , University of Paderborn , Warburger Str. 100 , D-33098 Paderborn , Germany
- Paderborn Center for Parallel Computing , University of Paderborn , Warburger Str. 100 , D-33098 Paderborn , Germany
- Center for Sustainable Systems Design , University of Paderborn , Warburger Str. 100 , D-33098 Paderborn , Germany
- Institute for Lightweight Design , University of Paderborn , Warburger Str. 100 , D-33098 Paderborn , Germany
| |
Collapse
|
6
|
Ruggeri M, Moroni S, Holzmann M. Nonlinear Network Description for Many-Body Quantum Systems in Continuous Space. PHYSICAL REVIEW LETTERS 2018; 120:205302. [PMID: 29864292 DOI: 10.1103/physrevlett.120.205302] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/08/2017] [Indexed: 06/08/2023]
Abstract
We show that the recently introduced iterative backflow wave function can be interpreted as a general neural network in continuum space with nonlinear functions in the hidden units. Using this wave function in variational Monte Carlo simulations of liquid ^{4}He in two and three dimensions, we typically find a tenfold increase in accuracy over currently used wave functions. Furthermore, subsequent stages of the iteration procedure define a set of increasingly good wave functions, each with its own variational energy and variance of the local energy: extrapolation to zero variance gives energies in close agreement with the exact values. For two dimensional ^{4}He, we also show that the iterative backflow wave function can describe both the liquid and the solid phase with the same functional form-a feature shared with the shadow wave function, but now joined by much higher accuracy. We also achieve significant progress for liquid ^{3}He in three dimensions, improving previous variational and fixed-node energies.
Collapse
Affiliation(s)
- Michele Ruggeri
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Saverio Moroni
- DEMOCRITOS National Simulation Center, Istituto Officina dei Materiali del CNR and SISSA, Via Bonomea 265, I-34136 Trieste, Italy
| | - Markus Holzmann
- Univ. Grenoble Alpes, CNRS, LPMMC, 3800 Grenoble, France
- Institut Laue Langevin, BP 156, F-38042 Grenoble Cedex 9, France
| |
Collapse
|
7
|
Azadi S, Kühne TD. High-pressure hydrogen sulfide by diffusion quantum Monte Carlo. J Chem Phys 2017; 146:084503. [DOI: 10.1063/1.4976836] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sam Azadi
- Department of Materials Science, Royal School of Mines, Thomas Young Center, London Centre for Nanotechnology, Imperial College London, London SW7 2AZ, United Kingdom
| | - Thomas D. Kühne
- Dynamics of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany and Paderborn Center for Parallel Computing and Institute for Lightweight Design with Hybrid Systems, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| |
Collapse
|
8
|
Calcavecchia F, Holzmann M. Fermion sign problem in imaginary-time projection continuum quantum Monte Carlo with local interaction. Phys Rev E 2016; 93:043321. [PMID: 27176442 DOI: 10.1103/physreve.93.043321] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2016] [Indexed: 06/05/2023]
Abstract
We use the shadow wave function formalism as a convenient model to study the fermion sign problem affecting all projector quantum Monte Carlo methods in continuum space. We demonstrate that the efficiency of imaginary-time projection algorithms decays exponentially with increasing number of particles and/or imaginary-time propagation. Moreover, we derive an analytical expression that connects the localization of the system with the magnitude of the sign problem, illustrating this behavior through numerical results. Finally, we discuss the computational complexity of the fermion sign problem and methods for alleviating its severity.
Collapse
Affiliation(s)
- Francesco Calcavecchia
- LPMMC, UMR 5493 of CNRS, Université Grenoble Alpes, 38042 Grenoble, France; Institute of Physics, Johannes Gutenberg University, Staudingerweg 7, D-55128 Mainz, Germany; and Graduate School of Excellence Materials Science in Mainz, Staudingerweg 9, D-55128 Mainz, Germany
| | - Markus Holzmann
- LPMMC, UMR 5493 of CNRS, Université Grenoble Alpes, 38042 Grenoble, France and Institut Laue Langevin, BP 156, F-38042 Grenoble Cedex 9, France
| |
Collapse
|