1
|
Schiltz C, Rappoport D, Mandelshtam VA. Implementation of the self-consistent phonons method with ab initio potentials (AI-SCP). J Chem Phys 2023; 158:2890485. [PMID: 37184023 DOI: 10.1063/5.0146682] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 04/24/2023] [Indexed: 05/16/2023] Open
Abstract
The self-consistent phonon (SCP) method allows one to include anharmonic effects when treating a many-body quantum system at thermal equilibrium. The system is then described by an effective temperature-dependent harmonic Hamiltonian, which can be used to estimate its various dynamic and static properties. In this paper, we combine SCP with ab initio (AI) potential energy evaluation in which case the numerical bottleneck of AI-SCP is the evaluation of Gaussian averages of the AI potential energy and its derivatives. These averages are computed efficiently by the quasi-Monte Carlo method utilizing low-discrepancy sequences leading to a fast convergence with respect to the number, S, of the AI energy evaluations. Moreover, a further substantial (an-order-of-magnitude) improvement in efficiency is achieved once a numerically cheap approximation of the AI potential is available. This is based on using a perturbation theory-like (the two-grid) approach in which it is the average of the difference between the AI and the approximate potential that is computed. The corresponding codes and scripts are provided.
Collapse
Affiliation(s)
- Colin Schiltz
- Department of Chemistry, University of California, Irvine, California 92697, USA
| | - Dmitrij Rappoport
- Department of Chemistry, University of California, Irvine, California 92697, USA
| | | |
Collapse
|
2
|
Parneix P, Maupin R, Attal L, Calvo F, Falvo C. Extracting vibrational anharmonicities from short driven molecular dynamics trajectories. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02740-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
3
|
Brown SE, Mandelshtam VA. Self-consistent phonons: An accurate and practical method to account for anharmonic effects in equilibrium properties of general classical or quantum many-body systems. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.06.008] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
4
|
Hermes MR, Hirata S. Stochastic many-body perturbation theory for anharmonic molecular vibrations. J Chem Phys 2015; 141:084105. [PMID: 25173003 DOI: 10.1063/1.4892614] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new quantum Monte Carlo (QMC) method for anharmonic vibrational zero-point energies and transition frequencies is developed, which combines the diagrammatic vibrational many-body perturbation theory based on the Dyson equation with Monte Carlo integration. The infinite sums of the diagrammatic and thus size-consistent first- and second-order anharmonic corrections to the energy and self-energy are expressed as sums of a few m- or 2m-dimensional integrals of wave functions and a potential energy surface (PES) (m is the vibrational degrees of freedom). Each of these integrals is computed as the integrand (including the value of the PES) divided by the value of a judiciously chosen weight function evaluated on demand at geometries distributed randomly but according to the weight function via the Metropolis algorithm. In this way, the method completely avoids cumbersome evaluation and storage of high-order force constants necessary in the original formulation of the vibrational perturbation theory; it furthermore allows even higher-order force constants essentially up to an infinite order to be taken into account in a scalable, memory-efficient algorithm. The diagrammatic contributions to the frequency-dependent self-energies that are stochastically evaluated at discrete frequencies can be reliably interpolated, allowing the self-consistent solutions to the Dyson equation to be obtained. This method, therefore, can compute directly and stochastically the transition frequencies of fundamentals and overtones as well as their relative intensities as pole strengths, without fixed-node errors that plague some QMC. It is shown that, for an identical PES, the new method reproduces the correct deterministic values of the energies and frequencies within a few cm(-1) and pole strengths within a few thousandths. With the values of a PES evaluated on the fly at random geometries, the new method captures a noticeably greater proportion of anharmonic effects.
Collapse
Affiliation(s)
- Matthew R Hermes
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| |
Collapse
|
5
|
Pereverzev A, Sewell TD. Obtaining the Hessian from the force covariance matrix: Application to crystalline explosives PETN and RDX. J Chem Phys 2015; 142:134110. [DOI: 10.1063/1.4916614] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andrey Pereverzev
- Department of Chemistry, University of Missouri-Columbia, Columbia, Missouri 65211-7600, USA
| | - Thomas D. Sewell
- Department of Chemistry, University of Missouri-Columbia, Columbia, Missouri 65211-7600, USA
| |
Collapse
|
6
|
|
7
|
Thaunay F, Dognon JP, Ohanessian G, Clavaguéra C. Vibrational mode assignment of finite temperature infrared spectra using the AMOEBA polarizable force field. Phys Chem Chem Phys 2015. [DOI: 10.1039/c5cp02270c] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The Driven Molecular Dynamics approach has been adapted and associated with the AMOEBA polarizable force field to assign and visualize vibrational modes in infrared spectra obtained by molecular dynamics simulations.
Collapse
Affiliation(s)
- Florian Thaunay
- Laboratoire de Chimie Moléculaire
- Ecole Polytechnique
- CNRS
- 91128 Palaiseau Cedex
- France
| | - Jean-Pierre Dognon
- CEA/Saclay
- DSM/IRAMIS/NIMBE
- CNRS
- UMR 3685
- Laboratoire de Chimie Moléculaire et de Catalyse pour l'Energie
| | - Gilles Ohanessian
- Laboratoire de Chimie Moléculaire
- Ecole Polytechnique
- CNRS
- 91128 Palaiseau Cedex
- France
| | - Carine Clavaguéra
- Laboratoire de Chimie Moléculaire
- Ecole Polytechnique
- CNRS
- 91128 Palaiseau Cedex
- France
| |
Collapse
|
8
|
Hermes MR, Hirata S. Stochastic algorithm for size-extensive vibrational self-consistent field methods on fully anharmonic potential energy surfaces. J Chem Phys 2014; 141:244111. [DOI: 10.1063/1.4904220] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Matthew R. Hermes
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
- CREST, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| |
Collapse
|
9
|
Calvo F, Falvo C, Parneix P. Atomistic Modeling of Vibrational Action Spectra in Polyatomic Molecules: Nuclear Quantum Effects. J Phys Chem A 2014; 118:5427-36. [DOI: 10.1021/jp5040147] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- F. Calvo
- ILM, Université Lyon I and CNRS UMR 5306, 43 Bd du 11 Novembre 1918, F69622 Villeurbanne Cedex, France
| | - C. Falvo
- Institut des Sciences
Moléculaires d’Orsay, UMR CNRS 8214, Université Paris Sud 11, Bât.
210, F91405 Orsay Cedex, France
| | - P. Parneix
- Institut des Sciences
Moléculaires d’Orsay, UMR CNRS 8214, Université Paris Sud 11, Bât.
210, F91405 Orsay Cedex, France
| |
Collapse
|
10
|
Georgescu I, Jitomirskaya S, Mandelshtam VA. On the ground state calculation of a many-body system using a self-consistent basis and quasi-Monte Carlo: An application to water hexamer. J Chem Phys 2013; 139:204104. [DOI: 10.1063/1.4829836] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
|
11
|
|
12
|
Simon A, Spiegelman F. Water clusters adsorbed on polycyclic aromatic hydrocarbons: Energetics and conformational dynamics. J Chem Phys 2013; 138:194309. [DOI: 10.1063/1.4805015] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Affiliation(s)
- Matthew R. Hermes
- Department of Chemistry, University of Illinois at Urbana−Champaign,
600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana−Champaign,
600 South Mathews Avenue, Urbana, Illinois 61801, United States
- CREST, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama
332-0012, Japan
| |
Collapse
|
14
|
Brown SE, Georgescu I, Mandelshtam VA. Self-consistent phonons revisited. II. A general and efficient method for computing free energies and vibrational spectra of molecules and clusters. J Chem Phys 2013; 138:044317. [DOI: 10.1063/1.4788977] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
15
|
Calvo F, Falvo C, Parneix P. A simple but accurate potential for the naphthalene-argon complex: Applications to collisional energy transfer and matrix isolated IR spectroscopy. J Chem Phys 2013; 138:034305. [DOI: 10.1063/1.4773469] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
16
|
Georgescu I, Mandelshtam VA. Self-consistent phonons revisited. I. The role of thermal versus quantum fluctuations on structural transitions in large Lennard-Jones clusters. J Chem Phys 2012; 137:144106. [DOI: 10.1063/1.4754819] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ionuţ Georgescu
- Chemistry Department, University of California, Irvine, California 92697, USA.
| | | |
Collapse
|
17
|
Simon A, Rapacioli M, Mascetti J, Spiegelman F. Vibrational spectroscopy and molecular dynamics of water monomers and dimers adsorbed on polycyclic aromatic hydrocarbons. Phys Chem Chem Phys 2012; 14:6771-86. [PMID: 22495405 DOI: 10.1039/c2cp40321h] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
This paper reports structures, energetics, dynamics and spectroscopy of H2O and (H2O)2 systems adsorbed on coronene (C24H12), a compact polycyclic aromatic hydrocarbon (PAH). On-the-fly Born-Oppenheimer molecular dynamics simulations are performed for temperatures T varying from 10 to 300 K, on a potential energy surface obtained within the self-consistent-charge density-functional based tight-binding (SCC-DFTB) approach. Anharmonic infrared (IR) spectra are extracted from these simulations. We first benchmark the SCC-DFTB semi-empirical hamiltonian vs. DFT (Density Functional Theory) calculations that include dispersion, on (C6H6)(H2O)1,2 small complexes. We find that charge corrections and inclusion of dispersion contributions in DFTB are necessary to obtain consistent structures, energetics and IR spectra. Using this Hamiltonian, the structures, energetics and IR features of the low-energy isomers of (C24H12)(H2O)1,2 are found to be similar to the DFT ones, with evidence for a stabilizing edge-coordination. The temperature dependence of the motions of H2O and (H2O)2 on the surface of C24H12 is analysed, revealing ultra-fast periodic motion. The water dimer starts diffusing at a higher temperature than the water monomer (150 K vs. 10 K respectively), which appears to be consistent with the binding energies. Qualitative and quantitative analyses of the effects of T on the IR spectra are performed. Anharmonic factors in particular are derived and it is shown that they can be used as signatures for the presence of PAH-water complexes. Finally, this paper lays the foundations for the studies of larger (PAH)m(H2O)n clusters, that can be treated with the efficient computational approach benchmarked in this paper.
Collapse
Affiliation(s)
- Aude Simon
- Laboratoire de Chimie et Physique Quantiques (LCPQ), Université de Toulouse [UPS] and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France.
| | | | | | | |
Collapse
|
18
|
Van-Oanh NT, Falvo C, Calvo F, Lauvergnat D, Basire M, Gaigeot MP, Parneix P. Improving anharmonic infrared spectra using semiclassically prepared molecular dynamics simulations. Phys Chem Chem Phys 2012; 14:2381-90. [DOI: 10.1039/c2cp23101h] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
19
|
Mathias G, Baer MD. Generalized Normal Coordinates for the Vibrational Analysis of Molecular Dynamics Simulations. J Chem Theory Comput 2011; 7:2028-39. [DOI: 10.1021/ct2001304] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Gerald Mathias
- Lehrstuhl für BioMolekulare Optik, Ludwig−Maximilians Universität München, Oettingenstrasse 67, 80538 München, Germany
| | - Marcel D. Baer
- Chemical and Materials Science Division, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| |
Collapse
|