1
|
Hadad RE, Roy A, Rabani E, Redmer R, Baer R. Stochastic density functional theory combined with Langevin dynamics for warm dense matter. Phys Rev E 2024; 109:065304. [PMID: 39020867 DOI: 10.1103/physreve.109.065304] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2024] [Accepted: 05/17/2024] [Indexed: 07/19/2024]
Abstract
This study overviews and extends a recently developed stochastic finite-temperature Kohn-Sham density functional theory to study warm dense matter using Langevin dynamics, specifically under periodic boundary conditions. The method's algorithmic complexity exhibits nearly linear scaling with system size and is inversely proportional to the temperature. Additionally, a linear-scaling stochastic approach is introduced to assess the Kubo-Greenwood conductivity, demonstrating exceptional stability for dc conductivity. Utilizing the developed tools, we investigate the equation of state, radial distribution, and electronic conductivity of hydrogen at a temperature of 30 000 K. As for the radial distribution functions, we reveal a transition of hydrogen from gaslike to liquidlike behavior as its density exceeds 4g/cm^{3}. As for the electronic conductivity as a function of the density, we identified a remarkable isosbestic point at frequencies around 7 eV, which may be an additional signature of a gas-liquid transition in hydrogen at 30 000 K.
Collapse
Affiliation(s)
| | | | - Eran Rabani
- Department of Chemistry, University of California, Berkeley, California 94720, USA; Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA; and The Raymond and Beverly Sackler Center of Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 69978, Israel
| | | | | |
Collapse
|
2
|
Xie H, Li ZH, Wang H, Zhang L, Wang L. Deep Variational Free Energy Approach to Dense Hydrogen. PHYSICAL REVIEW LETTERS 2023; 131:126501. [PMID: 37802941 DOI: 10.1103/physrevlett.131.126501] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/19/2022] [Revised: 08/01/2023] [Accepted: 08/14/2023] [Indexed: 10/08/2023]
Abstract
We developed a deep generative model-based variational free energy approach to the equations of state of dense hydrogen. We employ a normalizing flow network to model the proton Boltzmann distribution and a fermionic neural network to model the electron wave function at given proton positions. By jointly optimizing the two neural networks we reached a comparable variational free energy to the previous coupled electron-ion Monte Carlo calculation. The predicted equation of state of dense hydrogen under planetary conditions is denser than the findings of ab initio molecular dynamics calculation and empirical chemical model. Moreover, direct access to the entropy and free energy of dense hydrogen opens new opportunities in planetary modeling and high-pressure physics research.
Collapse
Affiliation(s)
- Hao Xie
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100190, China
| | - Zi-Hang Li
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100190, China
| | - Han Wang
- Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Fenghao East Road 2, Beijing 100094, China
| | - Linfeng Zhang
- DP Technology, Beijing 100080, China
- AI for Science Institute, Beijing 100080, China
| | - Lei Wang
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- Songshan Lake Materials Laboratory, Dongguan, Guangdong 523808, China
| |
Collapse
|
3
|
Tiihonen J, Clay RC, Krogel JT. Toward quantum Monte Carlo forces on heavier ions: Scaling properties. J Chem Phys 2021; 154:204111. [PMID: 34241166 DOI: 10.1063/5.0052266] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Quantum Monte Carlo (QMC) forces have been studied extensively in recent decades because of their importance with spectroscopic observables and geometry optimization. Here, we benchmark the accuracy and computational cost of QMC forces. The zero-variance zero-bias (ZVZB) force estimator is used in standard variational and diffusion Monte Carlo simulations with mean-field based trial wavefunctions and atomic pseudopotentials. Statistical force uncertainties are obtained with a recently developed regression technique for heavy tailed QMC data [P. Lopez Rios and G. J. Conduit, Phys. Rev. E 99, 063312 (2019)]. By considering selected atoms and dimers with elements ranging from H to Zn (1 ≤ Zeff ≤ 20), we assess the accuracy and the computational cost of ZVZB forces as the effective pseudopotential valence charge, Zeff, increases. We find that the costs of QMC energies and forces approximately follow simple power laws in Zeff. The force uncertainty grows more rapidly, leading to a best case cost scaling relationship of approximately Zeff 6.5(3) for diffusion Monte Carlo. We find that the accessible system size at fixed computational cost scales as Zeff -2, insensitive to model assumptions or the use of the "space warp" variance-reduction technique. Our results predict the practical cost of obtaining forces for a range of materials, such as transition metal oxides where QMC forces have yet to be applied, and underscore the importance of further developing force variance-reduction techniques, particularly for atoms with high Zeff.
Collapse
Affiliation(s)
- Juha Tiihonen
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Raymond C Clay
- Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Jaron T Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| |
Collapse
|
4
|
Nakano K, Attaccalite C, Barborini M, Capriotti L, Casula M, Coccia E, Dagrada M, Genovese C, Luo Y, Mazzola G, Zen A, Sorella S. TurboRVB: A many-body toolkit for ab initio electronic simulations by quantum Monte Carlo. J Chem Phys 2020; 152:204121. [DOI: 10.1063/5.0005037] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Kousuke Nakano
- International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy
- Japan Advanced Institute of Science and Technology (JAIST), Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
| | - Claudio Attaccalite
- Aix-Marseille Université, CNRS, CINaM UMR 7325, Campus de Luminy, 13288 Marseille, France
| | | | - Luca Capriotti
- New York University, Tandon School of Engineering, 6 MetroTech Center, Brooklyn, New York 11201, USA
- Department of Mathematics, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Michele Casula
- Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie (IMPMC), Sorbonne Université, CNRS UMR 7590, IRD UMR 206, MNHN, 4 Place Jussieu, 75252 Paris, France
| | - Emanuele Coccia
- Department of Chemical and Pharmaceutical Sciences, University of Trieste, Via L. Giorgieri 1, 34127 Trieste, Italy
| | - Mario Dagrada
- Forescout Technologies, John F. Kennedylaan 2, 5612AB Eindhoven, The Netherlands
| | - Claudio Genovese
- International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy
| | - Ye Luo
- Computational Science Division, Argonne National Laboratory, 9700 S. Cass Avenue, Lemont, Illinois 60439, USA
- Argonne Leadership Computing Facility, Argonne National Laboratory, 9700 S. Cass Avenue, Lemont, Illinois 60439, USA
| | | | - Andrea Zen
- Department of Earth Sciences, University College London, Gower Street, London WC1E 6BT, United Kingdom
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - Sandro Sorella
- International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy
| |
Collapse
|
5
|
Kühne TD, Iannuzzi M, Del Ben M, Rybkin VV, Seewald P, Stein F, Laino T, Khaliullin RZ, Schütt O, Schiffmann F, Golze D, Wilhelm J, Chulkov S, Bani-Hashemian MH, Weber V, Borštnik U, Taillefumier M, Jakobovits AS, Lazzaro A, Pabst H, Müller T, Schade R, Guidon M, Andermatt S, Holmberg N, Schenter GK, Hehn A, Bussy A, Belleflamme F, Tabacchi G, Glöß A, Lass M, Bethune I, Mundy CJ, Plessl C, Watkins M, VandeVondele J, Krack M, Hutter J. CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations. J Chem Phys 2020; 152:194103. [PMID: 33687235 DOI: 10.1063/5.0007045] [Citation(s) in RCA: 922] [Impact Index Per Article: 230.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.
Collapse
Affiliation(s)
- Thomas D Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Marcella Iannuzzi
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Mauro Del Ben
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Vladimir V Rybkin
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Patrick Seewald
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Frederick Stein
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Teodoro Laino
- IBM Research Europe, CH-8803 Rüschlikon, Switzerland
| | - Rustam Z Khaliullin
- Department of Chemistry, McGill University, CH-801 Sherbrooke St. West, Montreal, Quebec H3A 0B8, Canada
| | - Ole Schütt
- Department of Materials, ETH Zürich, CH-8092 Zürich, Switzerland
| | | | - Dorothea Golze
- Department of Applied Physics, Aalto University, Otakaari 1, FI-02150 Espoo, Finland
| | - Jan Wilhelm
- Institute of Theoretical Physics, University of Regensburg, Universitätsstraße 31, D-93053 Regensburg, Germany
| | - Sergey Chulkov
- School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln, United Kingdom
| | | | - Valéry Weber
- IBM Research Europe, CH-8803 Rüschlikon, Switzerland
| | | | | | | | | | - Hans Pabst
- Intel Extreme Computing, Software and Systems, Zürich, Switzerland
| | - Tiziano Müller
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Robert Schade
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Manuel Guidon
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Samuel Andermatt
- Integrated Systems Laboratory, ETH Zürich, CH-8092 Zürich, Switzerland
| | - Nico Holmberg
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, 00076 Aalto, Finland
| | - Gregory K Schenter
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Anna Hehn
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Augustin Bussy
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Fabian Belleflamme
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Gloria Tabacchi
- Department of Science and High Technology, University of Insubria and INSTM, via Valleggio 9, I-22100 Como, Italy
| | - Andreas Glöß
- BASF SE, Carl-Bosch-Straße 38, D-67056 Ludwigshafen am Rhein, Germany
| | - Michael Lass
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Iain Bethune
- Hartree Centre, Science and Technology Facilities Council, Sci-Tech Daresbury, Warrington WA4 4AD, United Kingdom
| | - Christopher J Mundy
- Physical Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, USA
| | - Christian Plessl
- Department of Computer Science and Paderborn Center for Parallel Computing, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Matt Watkins
- School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln, United Kingdom
| | - Joost VandeVondele
- Swiss National Supercomputing Centre (CSCS), ETH Zürich, Zürich, Switzerland
| | - Matthias Krack
- Laboratory for Scientific Computing and Modelling, Paul Scherrer Institute, CH-5232 Villigen PSI, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| |
Collapse
|
6
|
Arnon E, Rabani E, Neuhauser D, Baer R. Efficient Langevin dynamics for “noisy” forces. J Chem Phys 2020; 152:161103. [DOI: 10.1063/5.0004954] [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)
- Eitam Arnon
- Fritz Haber Research Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Eran Rabani
- Department of Chemistry, University of California and Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA and The Raymond and Beverly Sackler Center for Computational Molecular and Materials Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| | - Daniel Neuhauser
- Department of Chemistry and Biochemistry, University of California at Los Angeles, Los Angeles, California 90095, USA
| | - Roi Baer
- Fritz Haber Research Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| |
Collapse
|
7
|
Nakano K, Maezono R, Sorella S. All-Electron Quantum Monte Carlo with Jastrow Single Determinant Ansatz: Application to the Sodium Dimer. J Chem Theory Comput 2019; 15:4044-4055. [DOI: 10.1021/acs.jctc.9b00295] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kousuke Nakano
- International School for Advanced Studies (SISSA), Via Bonomea 265, 34136, Trieste, Italy
- School of Information Science, Japan Advanced Institute of Science and Technology (JAIST), Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
| | - Ryo Maezono
- School of Information Science, Japan Advanced Institute of Science and Technology (JAIST), Asahidai 1-1, Nomi, Ishikawa 923-1292, Japan
- Computational Engineering Applications Unit, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Sandro Sorella
- International School for Advanced Studies (SISSA), Via Bonomea 265, 34136, Trieste, Italy
| |
Collapse
|
8
|
Wang Z, Chern GW, Batista CD, Barros K. Gradient-based stochastic estimation of the density matrix. J Chem Phys 2018. [DOI: 10.1063/1.5017741] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Zhentao Wang
- Department of Physics and Astronomy, The University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Gia-Wei Chern
- Department of Physics, University of Virginia, Charlottesville, Virginia 22904, USA
| | - Cristian D. Batista
- Department of Physics and Astronomy, The University of Tennessee, Knoxville, Tennessee 37996, USA
- Quantum Condensed Matter Division and Shull-Wollan Center, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Kipton Barros
- Theoretical Division and CNLS, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
9
|
Mouhat F, Sorella S, Vuilleumier R, Saitta AM, Casula M. Fully Quantum Description of the Zundel Ion: Combining Variational Quantum Monte Carlo with Path Integral Langevin Dynamics. J Chem Theory Comput 2017; 13:2400-2417. [DOI: 10.1021/acs.jctc.7b00017] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Félix Mouhat
- IMPMC, Sorbonne Universités, Université Pierre et Marie Curie Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, 4 place Jussieu, 75252 Paris, France
| | - Sandro Sorella
- International
School for Advanced Studies (SISSA), Via Bonomea 26, 34136 Trieste, Italy
- INFM Democritos
National Simulation Center, 34151 Trieste, Italy
| | - Rodolphe Vuilleumier
- PASTEUR,
Département de chimie, École normale supérieure,
UPMC Univ Paris 06, CNRS, PSL Research University, 75005 Paris, France
- Sorbonne Universités, UPMC Univ Paris 06, École normale supérieure, CNRS, Processus d’activation
sélective par transfert d’énergie uni-électronique
ou radiatif (PASTEUR), 75005 Paris, France
| | - Antonino Marco Saitta
- IMPMC, Sorbonne Universités, Université Pierre et Marie Curie Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, 4 place Jussieu, 75252 Paris, France
| | - Michele Casula
- IMPMC, Sorbonne Universités, Université Pierre et Marie Curie Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, 4 place Jussieu, 75252 Paris, France
| |
Collapse
|
10
|
Mazzola G, Sorella S. Accelerating ab initio Molecular Dynamics and Probing the Weak Dispersive Forces in Dense Liquid Hydrogen. PHYSICAL REVIEW LETTERS 2017; 118:015703. [PMID: 28106448 DOI: 10.1103/physrevlett.118.015703] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/06/2016] [Indexed: 06/06/2023]
Abstract
We propose an ab initio molecular dynamics method, capable of dramatically reducing the autocorrelation time required for the simulation of classical and quantum particles at finite temperatures. The method is based on an efficient implementation of a first order Langevin dynamics modified by means of a suitable, position dependent acceleration matrix S. Here, we apply this technique to both Lennard-Jones models, to demonstrate the accuracy and speeding-up of the sampling, and within a quantum Monte Carlo based wave function approach, for determining the phase diagram of high-pressure hydrogen with simulations much longer than the autocorrelation time. With the proposed method, we are able to equilibrate in a few hundred steps even close to the liquid-liquid phase transition (LLT). Within our approach, we find that the LLT transition is consistent with recent density functionals predicting a much larger transition pressure when the long range dispersive forces are taken into account.
Collapse
Affiliation(s)
| | - Sandro Sorella
- International School for Advanced Studies (SISSA) Via Beirut 2, 4 34014 Trieste, Italy and INFM Democritos National Simulation Center, Via Beirut 2-4, I-34014 Trieste, Italy
| |
Collapse
|
11
|
Kawasaki A, Sugino O. Four-body correlation embedded in antisymmetrized geminal power wave function. J Chem Phys 2016; 145:244110. [DOI: 10.1063/1.4972991] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
12
|
Zen A, Coccia E, Gozem S, Olivucci M, Guidoni L. Quantum Monte Carlo Treatment of the Charge Transfer and Diradical Electronic Character in a Retinal Chromophore Minimal Model. J Chem Theory Comput 2016; 11:992-1005. [PMID: 25821414 PMCID: PMC4357234 DOI: 10.1021/ct501122z] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2014] [Indexed: 01/22/2023]
Abstract
![]()
The
penta-2,4-dieniminium cation (PSB3) displays similar ground
state and first excited state potential energy features as those of
the retinal protonated Schiff base (RPSB) chromophore in rhodopsin.
Recently, PSB3 has been used to benchmark several electronic structure
methods, including highly correlated multireference wave function
approaches, highlighting the necessity to accurately describe the
electronic correlation in order to obtain reliable properties even
along the ground state (thermal) isomerization paths. In this work,
we apply two quantum Monte Carlo approaches, the variational Monte
Carlo and the lattice regularized diffusion Monte Carlo, to study
the energetics and electronic properties of PSB3 along representative
minimum energy paths and scans related to its thermal cis–trans isomerization. Quantum Monte Carlo
is used in combination with the Jastrow antisymmetrized geminal power
ansatz, which guarantees an accurate and balanced description of the
static electronic correlation thanks to the multiconfigurational nature
of the antisymmetrized geminal power term, and of the dynamical correlation,
due to the presence of the Jastrow factor explicitly depending on
electron–electron distances. Along the two ground state isomerization
minimum energy paths of PSB3, CASSCF calculations yield wave functions
having either charge transfer or diradical character in proximity
of the two transition state configurations. Here, we observe that
at the quantum Monte Carlo level of theory, only the transition state
with charge transfer character can be located. The conical intersection,
which becomes highly sloped, is observed only if the path connecting
the two original CASSCF transition states is extended beyond the diradical
one, namely by increasing the bond-length-alternation (BLA). These
findings are in good agreement with the results obtained by MRCISD+Q
calculations, and they demonstrate the importance of having an accurate
description of the static and dynamical correlation when studying
isomerization and transition states of conjugated systems.
Collapse
|
13
|
Sorella S, Devaux N, Dagrada M, Mazzola G, Casula M. Geminal embedding scheme for optimal atomic basis set construction in correlated calculations. J Chem Phys 2015; 143:244112. [DOI: 10.1063/1.4938089] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- S. Sorella
- International School for Advanced Studies (SISSA), Via Beirut 2-4, 34014 Trieste, Italy and INFM Democritos National Simulation Center, Trieste, Italy
| | - N. Devaux
- Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, Case 115, 4 Place Jussieu, 75252 Paris Cedex 05, France
| | - M. Dagrada
- Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, Case 115, 4 Place Jussieu, 75252 Paris Cedex 05, France
| | - G. Mazzola
- Theoretische Physik, ETH Zurich, 8093 Zurich, Switzerland
| | - M. Casula
- CNRS and Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, Case 115, 4 Place Jussieu, 75252 Paris Cedex 05, France
| |
Collapse
|
14
|
Barborini M, Guidoni L. Ground State Geometries of Polyacetylene Chains from Many-Particle Quantum Mechanics. J Chem Theory Comput 2015; 11:4109-18. [PMID: 26405437 PMCID: PMC4570661 DOI: 10.1021/acs.jctc.5b00427] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2015] [Indexed: 11/28/2022]
Abstract
Due to the crucial role played by electron correlation, the accurate determination of ground state geometries of π-conjugated molecules is still a challenge for many quantum chemistry methods. Because of the high parallelism of the algorithms and their explicit treatment of electron correlation effects, Quantum Monte Carlo calculations can offer an accurate and reliable description of the electronic states and of the geometries of such systems, competing with traditional quantum chemistry approaches. Here, we report the structural properties of polyacetylene chains H-(C₂H₂)(N)-H up to N = 12 acetylene units, by means of Variational Monte Carlo (VMC) calculations based on the multi-determinant Jastrow Antisymmetrized Geminal Power (JAGP) wave function. This compact ansatz can provide for such systems an accurate description of the dynamical electronic correlation as recently detailed for the 1,3-butadiene molecule [J. Chem. Theory Comput. 2015 11 (2), 508-517]. The calculated Bond Length Alternation (BLA), namely the difference between the single and double carbon bonds, extrapolates, for N → ∞, to a value of 0.0910(7) Å, compatible with the experimental data. An accurate analysis was able to distinguish between the influence of the multi-determinantal AGP expansion and of the Jastrow factor on the geometrical properties of the fragments. Our size-extensive and self-interaction-free results provide new and accurate ab initio references for the structures of the ground state of polyenes.
Collapse
Affiliation(s)
- Matteo Barborini
- Dipartimento
di Ingegneria e Scienze dell’Informazione e Matematica, Università degli studi dell’Aquila, Via Vetoio (Coppito), 67100 L’Aquila, Italy
- Dipartimento
di Scienze Fisiche e Chimiche, Università
degli studi dell’Aquila, Via Vetoio (Coppito), 67100 L’Aquila, Italy
| | - Leonardo Guidoni
- Dipartimento
di Scienze Fisiche e Chimiche, Università
degli studi dell’Aquila, Via Vetoio (Coppito), 67100 L’Aquila, Italy
| |
Collapse
|
15
|
Dupuy N, Bouaouli S, Mauri F, Sorella S, Casula M. Vertical and adiabatic excitations in anthracene from quantum Monte Carlo: Constrained energy minimization for structural and electronic excited-state properties in the JAGP ansatz. J Chem Phys 2015; 142:214109. [DOI: 10.1063/1.4922048] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Nicolas Dupuy
- Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, case 115, 4 place Jussieu, 75252 Paris Cedex 05, France
| | - Samira Bouaouli
- Laboratoire de Chimie Théorique, Université Pierre et Marie Curie, case 115, 4 place Jussieu, 75252 Paris Cedex 05, France
| | - Francesco Mauri
- CNRS and Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, case 115, 4 place Jussieu, 75252 Paris Cedex 05, France
| | - Sandro Sorella
- International School for Advanced Studies (SISSA), Via Beirut 2-4, 34014 Trieste, Italy and INFM Democritos National Simulation Center, Trieste, Italy
| | - Michele Casula
- CNRS and Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, case 115, 4 place Jussieu, 75252 Paris Cedex 05, France
| |
Collapse
|
16
|
Zen A, Luo Y, Mazzola G, Guidoni L, Sorella S. Ab initio molecular dynamics simulation of liquid water by quantum Monte Carlo. J Chem Phys 2015; 142:144111. [PMID: 25877566 DOI: 10.1063/1.4917171] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Although liquid water is ubiquitous in chemical reactions at roots of life and climate on the earth, the prediction of its properties by high-level ab initio molecular dynamics simulations still represents a formidable task for quantum chemistry. In this article, we present a room temperature simulation of liquid water based on the potential energy surface obtained by a many-body wave function through quantum Monte Carlo (QMC) methods. The simulated properties are in good agreement with recent neutron scattering and X-ray experiments, particularly concerning the position of the oxygen-oxygen peak in the radial distribution function, at variance of previous density functional theory attempts. Given the excellent performances of QMC on large scale supercomputers, this work opens new perspectives for predictive and reliable ab initio simulations of complex chemical systems.
Collapse
Affiliation(s)
- Andrea Zen
- Dipartimento di Fisica, “La Sapienza” - Università di Roma, piazzale Aldo Moro 5, 00185 Rome, Italy
- London Centre for Nanotechnology, University College London, London WC1E 6BT, United Kingdom
| | - Ye Luo
- SISSA–International School for Advanced Studies, Via Bonomea 26, 34136 Trieste, Italy
- Democritos Simulation Center CNR–IOM Istituto Officina dei Materiali, 34151 Trieste, Italy
| | - Guglielmo Mazzola
- SISSA–International School for Advanced Studies, Via Bonomea 26, 34136 Trieste, Italy
- Democritos Simulation Center CNR–IOM Istituto Officina dei Materiali, 34151 Trieste, Italy
| | - Leonardo Guidoni
- Dipartimento di Fisica, “La Sapienza” - Università di Roma, piazzale Aldo Moro 5, 00185 Rome, Italy
- Dipartimento di Scienze Fisiche e Chimiche, Università degli Studi dell’ Aquila, via Vetoio, 67100 L’ Aquila, Italy
| | - Sandro Sorella
- SISSA–International School for Advanced Studies, Via Bonomea 26, 34136 Trieste, Italy
- Democritos Simulation Center CNR–IOM Istituto Officina dei Materiali, 34151 Trieste, Italy
| |
Collapse
|
17
|
Mazzola G, Sorella S. Distinct metallization and atomization transitions in dense liquid hydrogen. PHYSICAL REVIEW LETTERS 2015; 114:105701. [PMID: 25815949 DOI: 10.1103/physrevlett.114.105701] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/20/2014] [Indexed: 06/04/2023]
Abstract
We perform molecular dynamics simulations driven by accurate quantum Monte Carlo forces on dense liquid hydrogen. There is a recent report of a complete atomization transition between a mixed molecular-atomic liquid and a completely dissociated fluid in an almost unaccessible pressure range [Nat. Commun. 5, 3487 (2014)]. Here, instead, we identify a different transition between the fully molecular liquid and the mixed-atomic fluid at ∼400 GPa, i.e., in a much more interesting pressure range. We provide numerical evidence supporting the metallic behavior of this intermediate phase. Therefore, we predict that the metallization at finite temperature occurs in this partially dissociated molecular fluid, well before the complete atomization of the liquid. At high temperature this first-order transition becomes a crossover, in very good agreement with the experimental observation. Several systematic tests supporting the quality of our large scale calculations are also reported.
Collapse
Affiliation(s)
- Guglielmo Mazzola
- International School for Advanced Studies (SISSA), and CRS Democritos, CNR-INFM, -Via Bonomea 265, I-34136 Trieste, Italy
| | - Sandro Sorella
- International School for Advanced Studies (SISSA), and CRS Democritos, CNR-INFM, -Via Bonomea 265, I-34136 Trieste, Italy
| |
Collapse
|