1
|
Stachowiak M, Grabowska E, Wang XG, Carrington T, Szalewicz K, Jankowski P. Theory cracks old data: Rovibrational energy levels of orthoH 2-CO derived from experiment. Sci Adv 2024; 10:eadj8632. [PMID: 38394212 PMCID: PMC10889352 DOI: 10.1126/sciadv.adj8632] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/20/2023] [Accepted: 01/22/2024] [Indexed: 02/25/2024]
Abstract
Measurements of rovibrational spectra of clusters provide physical insight only if spectral lines can be assigned to pairs of quantum states, and further insight is obtained if one can deduce the quantitative energy-level pattern. Both steps can be so difficult that some measured spectra remain unassigned, one example is orthoH2-CO. To extend the scope of spectroscopic insights, we propose to use theoretical information in interpretation of spectra. We first performed high accuracy, full-dimensional calculations of the orthoH2-CO spectrum, at the highest practically achievable levels of electronic structure theory and quantum nuclear dynamics. Then, an iterative, theory-guided method developed here allowed us to fully interpret the spectrum of orthoH2-CO, extending the range of van der Waals clusters for which spectroscopy can provide physical insights.
Collapse
Affiliation(s)
- Marcin Stachowiak
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Gagarina 7, 87-100 Toruń, Poland
| | - Ewelina Grabowska
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Gagarina 7, 87-100 Toruń, Poland
| | - Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, USA
| | - Piotr Jankowski
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Gagarina 7, 87-100 Toruń, Poland
| |
Collapse
|
2
|
Szalewicz K, Jeziorski B. Physical mechanisms of intermolecular interactions from symmetry-adapted perturbation theory. J Mol Model 2022; 28:273. [PMID: 36006512 DOI: 10.1007/s00894-022-05190-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Accepted: 05/12/2022] [Indexed: 10/15/2022]
Abstract
Symmetry-adapted perturbation theory (SAPT) is a method for computational studies of noncovalent interactions between molecules. This method will be discussed here from the perspective of establishing the paradigm for understanding mechanisms of intermolecular interactions. SAPT interaction energies are obtained as sums of several contributions. Each contribution possesses a clear physical interpretation as it results from some specific physical process. It also exhibits a specific dependence on the intermolecular separation R. The four major contributions are the electrostatic, induction, dispersion, and exchange energies, each due to a different mechanism, valid at any R. In addition, at large R, SAPT interaction energies are seamlessly connected with the corresponding terms in the asymptotic multipole expansion of interaction energy in inverse powers of R. Since such expansion explicitly depends on monomers' multipole moments and polarizabilities, this connection provides additional insights by rigorously relating interaction energies to monomers' properties.
Collapse
Affiliation(s)
- Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, DE, 19716, USA.
| | - Bogumił Jeziorski
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02093, Warsaw, Poland
| |
Collapse
|
3
|
Wiscons RA, Nikhar R, Szalewicz K, Matzger AJ. Factors influencing hydrogen peroxide versus water inclusion in molecular crystals. Phys Chem Chem Phys 2022; 24:11206-11212. [PMID: 35481469 DOI: 10.1039/d1cp05765k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Hydrate formation is often unavoidable during crystallization, leading to performance degradation of pharmaceuticals and energetics. In some cases, water molecules trapped within crystal lattices can be substituted for hydrogen peroxide, improving the solubility of drugs and detonation performance of explosives. The present work compares hydrates and hydrogen peroxide solvates in two ways: (1) analyzing structural motifs present in crystal structures accessed from the Cambridge Structural Database and (2) developing potential energy surfaces for water and hydrogen peroxide interacting with functional groups of interest at geometries relevant to the solid state. By elucidating fundamental differences in local interactions that can be formed with molecules of hydrogen peroxide and/or water, the analyses presented here provide a foundation for the design and selection of candidate molecules for the formation of hydrogen peroxide solvates.
Collapse
Affiliation(s)
- Ren A Wiscons
- Department of Chemistry and the Macromolecular Science and Engineering Program, University of Michigan, 930 North University of Ave, Ann Arbor, Michigan 48109-1055, USA.
| | - Rahul Nikhar
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware, 19716, USA. szalewic.@udel.edu
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware, 19716, USA. szalewic.@udel.edu
| | - Adam J Matzger
- Department of Chemistry and the Macromolecular Science and Engineering Program, University of Michigan, 930 North University of Ave, Ann Arbor, Michigan 48109-1055, USA.
| |
Collapse
|
4
|
Jankowski P, Grabowska E, Szalewicz K. On the role of coupled-clusters' full triple and perturbative quadruple excitations on rovibrational spectra of van der Waals complexes. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1955989] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Piotr Jankowski
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Toruń, Poland
| | - Ewelina Grabowska
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Toruń, Poland
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, DE, USA
| |
Collapse
|
5
|
Jedwabny W, Dyguda-Kazimierowicz E, Pernal K, Szalewicz K, Patkowski K. Extension of an Atom-Atom Dispersion Function to Halogen Bonds and Its Use for Rational Design of Drugs and Biocatalysts. J Phys Chem A 2021; 125:1787-1799. [PMID: 33620223 PMCID: PMC8028329 DOI: 10.1021/acs.jpca.0c11347] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
![]()
A dispersion function Das in the form
of a damped atom–atom asymptotic expansion fitted to ab initio dispersion energies from symmetry-adapted perturbation
theory was improved and extended to systems containing heavier halogen
atoms. To illustrate its performance, the revised Das function was implemented in the multipole first-order
electrostatic and second-order dispersion (MED) scoring model. The
extension has allowed applications to a much larger set of biocomplexes
than it was possible with the original Das. A reasonable correlation between MED and experimentally determined
inhibitory activities was achieved in a number of test cases, including
structures featuring nonphysically shortened intermonomer distances,
which constitute a particular challenge for binding strength predictions.
Since the MED model is also computationally efficient, it can be used
for reliable and rapid assessment of the ligand affinity or multidimensional
scanning of amino acid side-chain conformations in the process of
rational design of novel drugs or biocatalysts.
Collapse
Affiliation(s)
- Wiktoria Jedwabny
- Department of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| | - Edyta Dyguda-Kazimierowicz
- Department of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| | - Katarzyna Pernal
- Institute of Physics, Łódź University of Technology, Wólczańska 219, 90-924 Łódź, Poland
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| |
Collapse
|
6
|
Abstract
Symmetry-adapted perturbation theory (SAPT) is a method for calculations of intermolecular (noncovalent) interaction energies. The set of SAPT codes that is described here, the current version named SAPT2020, includes virtually all variants of SAPT developed so far, among them two-body SAPT based on perturbative, coupled cluster, and density functional theory descriptions of monomers, three-body SAPT, and two-body SAPT for some classes of open-shell monomers. The properties of systems governed by noncovalent interactions can be predicted only if potential energy surfaces (force fields) are available. SAPT is the preferred approach for generating such surfaces since it is seamlessly connected to the asymptotic expansion of interaction energy. SAPT2020 includes codes for automatic development of such surfaces, enabling generation of complete dimer surfaces with a rigid monomer approximation for dimers containing about one hundred atoms. These codes can also be used to obtain surfaces including internal degrees of freedom of monomers.
Collapse
Affiliation(s)
- Javier Garcia
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Rafał Podeszwa
- Institute of Chemistry, University of Silesia at Katowice, Szkolna 9, Katowice, Poland
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
7
|
Abstract
A method is developed for automatic generation of nonreactive intermolecular two-body potential energy surfaces (PESs) including intramonomer degrees of freedom. This method, called flex-autoPES, is an extension of the autoPES method developed earlier, which assumes rigid monomers. In both cases, the whole PES development proceeds without any human intervention. The functional form used is a sum of products of site-site functions (both atomic and off-atomic sites can be used). The leading terms with sites involving different monomers are of physically motivated form. The long-range part of a PES is computed from monomer properties without using any dimer information. The close-range part is fitted to dimer interaction energies computed using electronic structure methods. Virtually any method can be used in such calculations, but the use of symmetry-adapted perturbation theory provides a seamless connection to the long-range part of the PES. The performance of the flex-autoPES code was tested by developing a full-dimensional PES for the water dimer and PESs including only some soft intramonomer degrees of freedom for the ethylene glycol dimer and for the ethylene glycol-water dimer. In the case of the water dimer, the root-mean-square error (RMSE) of the PES from the data points with negative total energies is 0.03 kcal/mol, and we expect this PES to be more accurate than any previously published PES of this type. For the ethylene glycol dimer and the ethylene glycol-water dimers, the analogous RMSEs are 0.25 and 0.1 kcal/mol, respectively.
Collapse
Affiliation(s)
- Michael P Metz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
8
|
Metz MP, Szalewicz K. A statistically guided grid generation method and its application to intermolecular potential energy surfaces. J Chem Phys 2020; 152:134111. [PMID: 32268757 DOI: 10.1063/1.5141777] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a method for the generation of points in space needed to create training data for fitting of nonlinear parametric models. This method uses statistical information extracted from an initial fit on a sparse grid to select optimal grid points in an iterative manner and is, therefore, called the iterative variance minimizing grid approach. We demonstrate the method in the case of six-dimensional intermolecular potential energy surfaces (PESs) fitted to ab initio computed interaction energies. The number of required grid points is reduced by roughly a factor of two in comparison to alternative systematic sampling methods. The method is not limited to fitting PESs and can be applied to any cases of fitting parametric models where data points may be chosen freely but are expensive to obtain.
Collapse
Affiliation(s)
- Michael P Metz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
9
|
Abstract
The Hartree-Fock plus dispersion plus first-order correlation (HFDc(1)) method consists in augmenting the HF interaction energy by the correlation part of the first-order interaction energy and the second-order dispersion and exchange-dispersion energies. All of the augmentation terms are computed using the symmetry-adapted perturbation theory based on density functional theory description of monomers [SAPT(DFT)]; thus, HFDc(1) is a fully ab initio method. A partly empirical version of this method, HFDasc(1), uses a damped asymptotic expansion for the dispersion plus exchange-dispersion term fitted to SAPT(DFT) ab initio values. The HFDc(1) interaction energies for dimers in the S22, S66, S66x8, NCCE31, IonHB, and UD-ARL benchmark data sets are more accurate than those given by most ab initio methods with comparable costs. HFDc(1) can be used routinely for dimers with nearly 200 atoms, such as included in the S12L benchmark set, giving results comparable to those obtained by the most expensive methods applicable.
Collapse
Affiliation(s)
- Javier Garcia
- Department of Physics and Astronomy , University of Delaware , Newark , Delaware 19716 , United States
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy , University of Delaware , Newark , Delaware 19716 , United States
| |
Collapse
|
10
|
Andrés J, Ayers PW, Boto RA, Carbó-Dorca R, Chermette H, Cioslowski J, Contreras-García J, Cooper DL, Frenking G, Gatti C, Heidar-Zadeh F, Joubert L, Martín Pendás Á, Matito E, Mayer I, Misquitta AJ, Mo Y, Pilmé J, Popelier PLA, Rahm M, Ramos-Cordoba E, Salvador P, Schwarz WHE, Shahbazian S, Silvi B, Solà M, Szalewicz K, Tognetti V, Weinhold F, Zins ÉL. Nine questions on energy decomposition analysis. J Comput Chem 2019; 40:2248-2283. [PMID: 31251411 DOI: 10.1002/jcc.26003] [Citation(s) in RCA: 79] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2019] [Accepted: 05/16/2019] [Indexed: 01/05/2023]
Abstract
The paper collects the answers of the authors to the following questions: Is the lack of precision in the definition of many chemical concepts one of the reasons for the coexistence of many partition schemes? Does the adoption of a given partition scheme imply a set of more precise definitions of the underlying chemical concepts? How can one use the results of a partition scheme to improve the clarity of definitions of concepts? Are partition schemes subject to scientific Darwinism? If so, what is the influence of a community's sociological pressure in the "natural selection" process? To what extent does/can/should investigated systems influence the choice of a particular partition scheme? Do we need more focused chemical validation of Energy Decomposition Analysis (EDA) methodology and descriptors/terms in general? Is there any interest in developing common benchmarks and test sets for cross-validation of methods? Is it possible to contemplate a unified partition scheme (let us call it the "standard model" of partitioning), that is proper for all applications in chemistry, in the foreseeable future or even in principle? In the end, science is about experiments and the real world. Can one, therefore, use any experiment or experimental data be used to favor one partition scheme over another? © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Juan Andrés
- Departament de Ciències Experimentals Universitat Jaume I, 12080, Castelló, Spain
| | - Paul W Ayers
- Department of Chemistry and Chemical Biology, McMaster University, 1280 Main Street West, L8S 4M1, Hamilton, Ontario, Canada
| | | | - Ramon Carbó-Dorca
- Institut de Química Computational i Catàlisi, Universitat de Girona, C/M Aurelia Capmany 69, 17003, Girona, Spain
| | - Henry Chermette
- Université Lyon 1 et UMR CNRS 5280 Institut Sciences Analytiques, Université de Lyon, 69622, Paris, France
| | - Jerzy Cioslowski
- Institute of Physics, University of Szczecin, Wielkopolska, 15, 70-451, Szczecin, Poland
| | | | - David L Cooper
- Department of Chemistry, University of Liverpool, Liverpool, L69 7ZD, United Kingdom
| | - Gernot Frenking
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerweinstr. 4, 35032, Marburg, Germany
| | - Carlo Gatti
- CNR-ISTM Istituto di Scienze e Tecnologie Molecolari, via Golgi 19, 20133, Milan, Italy and Istituto Lombardo Accademia di Scienze e Lettere, via Brera 28, 20121, Milan, Italy
| | - Farnaz Heidar-Zadeh
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg and Department of Chemistry, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Laurent Joubert
- COBRA UMR 6014 & FR 3038, INSA Rouen, CNRS, Université de Rouen Normandie, Mont-St-Aignan, France
| | - Ángel Martín Pendás
- Departamento de Química Física y Analítica, Universidad de Oviedo, 33006, Oviedo, Spain
| | - Eduard Matito
- Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU), and Donostia International Physics Center (DIPC), P.K. 1072, 20080, Donostia, Euskadi, Spain.,IKERBASQUE, Basque Foundation for Science, 48011, Bilbao, Euskadi, Spain
| | - István Mayer
- Institute of Organic Chemistry, Research Centre for Natural Sciences, Hungarian Academy of Sciences, Budapest, 1117, Hungary
| | - Alston J Misquitta
- School of Physics and Astronomy, Queen Mary University of London, Mile End Road, London, E1 4NS, United Kingdom
| | - Yirong Mo
- Chemistry Department, Western Michigan University, Kalamazoo, Michigan, 49008
| | - Julien Pilmé
- Sorbonne Université, CNRS, LCT, UMR 7616, 4 place Jussieu, 75005, Paris, France
| | - Paul L A Popelier
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester, M1 7DN, United Kingdom.,School of Chemistry, University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom
| | - Martin Rahm
- Department of Chemistry and Chemical Engineering, Chalmers University of Technology, 412 96, Gothenburg, Sweden
| | - Eloy Ramos-Cordoba
- Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU), and Donostia International Physics Center (DIPC), P.K. 1072, 20080, Donostia, Euskadi, Spain
| | - Pedro Salvador
- Institut de Química Computacional i Catàlisi, Universitat de Girona, C/M Aurelia Capmany 69, 17003, Girona, Spain
| | - W H Eugen Schwarz
- Theoretical Chemistry Center at Tsinghua University, Beijing, 100084, China.,Physical and Theoretical Chemistry Laboratory, Faculty of Science and Engineering, University of Siegen, Siegen, 57068, Germany
| | - Shant Shahbazian
- Department of Physics, Shahid Beheshti University, P.O. Box 19395-4716, G. C., Evin, 19839, Tehran, Iran
| | - Bernard Silvi
- Sorbonne Université, CNRS, LCT, UMR 7616, 4 place Jussieu, 75005, Paris, France
| | - Miquel Solà
- Institut de Química Computacional i Catàlisi, Universitat de Girona, C/M Aurelia Capmany 69, 17003, Girona, Spain
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware
| | - Vincent Tognetti
- COBRA UMR 6014 & FR 3038, INSA Rouen, CNRS, Université de Rouen Normandie, Mont-St-Aignan, France
| | - Frank Weinhold
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin, 53706
| | - Émilie-Laure Zins
- Sorbonne Université, UPMC Univ. Paris 06, MONARIS, UMR 8233, Université Pierre et Marie Curie, 4 Place Jussieu, Case Courrier 49, 75252, Paris, France
| |
Collapse
|
11
|
Abstract
A simple nonlocal functional for calculation of dispersion energies is proposed. Compared to a similar formula used earlier, we introduced a regularization to remove its singularities and used a dynamic polarizability density similar to those in the so-called van der Waals density functionals. The performance of the new functional is tested on dispersion energies for a set of representative dimers, and it is found that it is significantly more accurate than published nonlocal functionals.
Collapse
Affiliation(s)
- Muhammad Shahbaz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
12
|
Abstract
Molecular dynamics (MD) simulations of methane-water mixtures were performed using ab initio force fields for the CH4-H2O, H2O-H2O, and CH4-CH4 interactions. Both methane and water molecules were polarizable. From these calculations, the potential of mean force (PMF) between two methane molecules was extracted. Our results are compared with PMFs from a density-functional-theory (DFT) based Born-Oppenheimer type MD (BOMD) simulation, from a Monte Carlo (MC) simulation with ab initio-based force fields, and from MD simulations with empirical force fields. Our PMF is qualitatively similar to that obtained from the simulations with empirical force fields but differs significantly from those resulting from the DFT-BOMD and MC simulations. The depth of the PMF global minimum obtained in the present work is in a much better agreement with the experimental estimate than the result of the DFT-BOMD simulation, possibly due to the inability of DFT to describe the dispersion interactions and the lack of extensive sampling in the BOMD simulations. Our work indicates that, for a pair of methane molecules, there are configurations where the solvent increases the attraction between the solutes, but there are also conformations in which the solvent causes a weak net repulsion. On average, the methane molecules are more likely to be in the configuration where they are separated by a water molecule than in the one in which they are in contact even though the minimum of the PMF at the latter configuration is deeper than that at the former. Finally, we found that the water structure around methane solutes does not show a greater tetrahedral ordering than in neat bulk water.
Collapse
Affiliation(s)
- Omololu Akin-Ojo
- ICTP East Africa Institute for Fundamental Research, University of Rwanda, Rwanda
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
13
|
Shahbaz M, Szalewicz K. Evaluation of methods for obtaining dispersion energies used in density functional calculations of intermolecular interactions. Theor Chem Acc 2019. [DOI: 10.1007/s00214-019-2414-5] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
14
|
Metz MP, Szalewicz K, Sarka J, Tóbiás R, Császár AG, Mátyus E. Molecular dimers of methane clathrates: ab initio potential energy surfaces and variational vibrational states. Phys Chem Chem Phys 2019; 21:13504-13525. [DOI: 10.1039/c9cp00993k] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Motivated by the energetic and environmental relevance of methane clathrates, highly accurate ab initio potential energy surfaces (PESs) have been developed for the three possible dimers of the methane and water molecules: (H2O)2, CH4·H2O, and (CH4)2.
Collapse
Affiliation(s)
- Michael P. Metz
- Department of Physics and Astronomy
- University of Delaware
- Newark
- USA
| | | | - János Sarka
- Department of Chemistry and Biochemistry
- Texas Tech University
- Lubbock
- USA
| | - Roland Tóbiás
- Institute of Chemistry
- ELTE Eötvös Loránd University
- Budapest
- Hungary
- MTA-ELTE Complex Chemical Systems Research Group
| | - Attila G. Császár
- Institute of Chemistry
- ELTE Eötvös Loránd University
- Budapest
- Hungary
- MTA-ELTE Complex Chemical Systems Research Group
| | - Edit Mátyus
- Institute of Chemistry
- ELTE Eötvös Loránd University
- Budapest
- Hungary
| |
Collapse
|
15
|
Babikov D, Benoit D, Bowman J, Burd T, Clary D, Donovan R, Fischer I, Gianturco F, Hochlaf M, Kar S, Kirrander A, Leone S, Malcomson T, Manthe U, McCoy AB, Petersen J, Richardson J, Slavíček P, Stoecklin T, Szalewicz K, van der Avoird A, Wester R, Worth G, Zehnacker-Rentien A. Quantum dynamics of isolated molecules: general discussion. Faraday Discuss 2018; 212:281-306. [PMID: 30511069 DOI: 10.1039/c8fd90052c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
16
|
Bacic Z, Benoit D, Biczysko M, Bowman J, Bradforth S, Burd T, Chambaud G, Clary D, Crépin C, Dracinsky M, Felker P, Fischer I, Gianturco F, Hochlaf M, Kouril K, Kratochvilova I, Liu C, McCoy A, Miyazaki J, Mouhib H, Richardson J, Slaviček P, Stoecklin T, Szalewicz K, van der Avoird A, Zehnacker-Rentien A. Molecules in confinement in clusters, quantum solvents and matrices: general discussion. Faraday Discuss 2018; 212:569-601. [PMID: 30520925 DOI: 10.1039/c8fd90053a] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
17
|
Garberoglio G, Jankowski P, Szalewicz K, Harvey AH. Fully quantum calculation of the second and third virial coefficients of water and its isotopologues from ab initio potentials. Faraday Discuss 2018; 212:467-497. [PMID: 30302450 PMCID: PMC6561489 DOI: 10.1039/c8fd00092a] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Path-Integral Monte Carlo methods were applied to calculate the second, B(T), and the third, C(T), virial coefficients for water. A fully quantum approach and state-of-the-art flexible-monomer pair and three-body potentials were used. Flexible-monomer potentials allow calculations for any isotopologue; we performed calculations for both H2O and D2O. For B(T) of H2O, the quantum effect contributes 25% of the value at 300 K and is not entirely negligible even at 1000 K, in accordance with recent literature findings. The effect of monomer flexibility, while not as large as some claims in the literature, is significant compared to the experimental uncertainty. It is of opposite sign to the quantum effect, smaller in magnitude than the latter below 500 K, and varies from 2% at 300 K to 10% at 700 K. When monomer flexibility is accounted for, results from the CCpol-8sf pair potential are in excellent agreement with the available experimental data and provide reliable B(T) values at temperatures outside the range of experimental data. The flexible-monomer MB-pol pair potential yields B(T) values that are slightly too high compared to experiment. For C(T), our calculations confirm earlier findings that the use of three-body potential is necessary for meaningful predictions. However, due to various uncertainties of the potentials used, especially the three-body ones, we were not able to establish benchmark values of C(T), although our results are in qualitative agreement with available experimental data. The quantum effect, never before included for water, reduces the magnitude of the classical value for H2O by a factor of 2.5 at 300 K and is not entirely negligible even at 1000 K.
Collapse
Affiliation(s)
- Giovanni Garberoglio
- European Centre for Theoretical Studies in Nuclear Physics and Related Areas (FBK-ECT*), strada delle Tabarelle 286, I-38123 Trento, Italy. and Trento Institute for Fundamental Physics and Applications (TIFPA-INFN), via Sommarive 18, I-38213 Trento, Italy
| | - Piotr Jankowski
- Faculty of Chemistry, Nicolaus Copernicus University in Torun, Gagarina 7, PL-87-100 Torun, Poland.
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA.
| | - Allan H Harvey
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305-3337, USA.
| |
Collapse
|
18
|
Bacic Z, Benoit D, Besemer M, Bowman J, Bradforth S, Clary D, Donovan R, Fischer I, Gianturco F, Hochlaf M, Houston P, Knowles P, Leone S, Linguerri R, Manthe U, McCoy AB, Petersen J, Richardson J, Shan X, Slavíček P, Stoecklin T, Szalewicz K, van der Avoird A, Wester R, Worth G, Zehnacker-Rentien A. Precise characterisation of isolated molecules: general discussion. Faraday Discuss 2018; 212:137-155. [PMID: 30507986 DOI: 10.1039/c8fd90050g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
19
|
Shahbaz M, Szalewicz K. Do Semilocal Density-Functional Approximations Recover Dispersion Energies at Small Intermonomer Separations? Phys Rev Lett 2018; 121:113402. [PMID: 30265106 DOI: 10.1103/physrevlett.121.113402] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2018] [Revised: 06/12/2018] [Indexed: 06/08/2023]
Abstract
The methods that add dispersion energies to interaction energies computed using density-functional theory (DFT), known as DFT+D methods, taper off the dispersion energies at distances near van der Waals minima and smaller based on an assumption that DFT starts to reproduce the dispersion energies there. We show that this assumption is not correct as the alleged contribution behaves unphysically and originates to a large extent from nonexchange-correlation terms. Thus, dispersion functions correct DFT in this region for deficiencies unrelated to dispersion interactions.
Collapse
Affiliation(s)
- Muhammad Shahbaz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
20
|
Affiliation(s)
- Anthony J. Stone
- University
Chemical Laboratory, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| | - Krzysztof Szalewicz
- Department
of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
21
|
Ban L, Bowman J, Bradforth S, Chambaud G, Dracinsky M, Fischer I, Góra R, Hochlaf M, Janicki M, Kirrander A, McCoy AB, Petersen J, Richardson J, Slavicek P, Szalewicz K, Zehnacker-Rentien A. Molecules in confinement in liquid solvents: general discussion. Faraday Discuss 2018; 212:383-397. [DOI: 10.1039/c8fd90051e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
22
|
Przybytek M, Cencek W, Jeziorski B, Szalewicz K. Pair Potential with Submillikelvin Uncertainties and Nonadiabatic Treatment of the Halo State of the Helium Dimer. Phys Rev Lett 2017; 119:123401. [PMID: 29341636 DOI: 10.1103/physrevlett.119.123401] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/18/2017] [Indexed: 06/07/2023]
Abstract
The pair potential for helium is computed with accuracy improved by an order of magnitude relative to the best previous determination. For the well region, its uncertainties are now below 1 millikelvin. The main improvement is due to the use of explicitly correlated wave functions at the nonrelativistic Born-Oppenheimer (BO) level of theory. The diagonal BO and the relativistic corrections are obtained from large full configuration interaction calculations. Nonadiabatic perturbation theory is used to predict the properties of the halo state of the helium dimer. Its binding energy and the average value of the interatomic distance are found to be 138.9(5) neV and 47.13(8) Å. The binding energy agrees with its first experimental determination of 151.9(13.3) neV [Zeller et al., Proc. Natl. Acad. Sci. U.S.A. 113, 14651 (2016)PNASA60027-842410.1073/pnas.1610688113].
Collapse
Affiliation(s)
- Michał Przybytek
- Department of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Wojciech Cencek
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Bogumił Jeziorski
- Department of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
23
|
Garberoglio G, Jankowski P, Szalewicz K, Harvey AH. All-dimensional H 2-CO potential: Validation with fully quantum second virial coefficients. J Chem Phys 2017; 146:054304. [PMID: 28178790 DOI: 10.1063/1.4974993] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We use a new high-accuracy all-dimensional potential to compute the cross second virial coefficient B12(T) between molecular hydrogen and carbon monoxide. The path-integral method is used to fully account for quantum effects. Values are calculated from 10 K to 2000 K and the uncertainty of the potential is propagated into uncertainties of B12. Our calculated B12(T) are in excellent agreement with most of the limited experimental data available, but cover a much wider range of temperatures and have lower uncertainties. Similar to recently reported findings from scattering calculations, we find that the reduced-dimensionality potential obtained by averaging over the rovibrational motion of the monomers gives results that are a good approximation to those obtained when flexibility is fully taken into account. Also, the four-dimensional approximation with monomers taken at their vibrationally averaged bond lengths works well. This finding is important, since full-dimensional potentials are difficult to develop even for triatomic monomers and are not currently possible to obtain for larger molecules. Likewise, most types of accurate quantum mechanical calculations, e.g., spectral or scattering, are severely limited in the number of dimensions that can be handled.
Collapse
Affiliation(s)
- Giovanni Garberoglio
- European Centre for Theoretical Studies in Nuclear Physics and Related Areas (ECT*-FBK), via Sommarive 18, I-38123 Trento, Italy
| | - Piotr Jankowski
- Faculty of Chemistry, Nicolaus Copernicus University in Torun, Gagarina 7, PL-87-100 Torun, Poland
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Allan H Harvey
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, 325 Broadway, Boulder, Colorado 80305-3337, USA
| |
Collapse
|
24
|
Affiliation(s)
- Michael P. Metz
- Department of Physics and
Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Konrad Piszczatowski
- Department of Physics and
Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Krzysztof Szalewicz
- Department of Physics and
Astronomy, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
25
|
Taylor DE, Ángyán JG, Galli G, Zhang C, Gygi F, Hirao K, Song JW, Rahul K, Anatole von Lilienfeld O, Podeszwa R, Bulik IW, Henderson TM, Scuseria GE, Toulouse J, Peverati R, Truhlar DG, Szalewicz K. Blind test of density-functional-based methods on intermolecular interaction energies. J Chem Phys 2016; 145:124105. [DOI: 10.1063/1.4961095] [Citation(s) in RCA: 81] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- DeCarlos E. Taylor
- U.S. Army Research Laboratory,
Aberdeen Proving Ground, Aberdeen, Maryland 21005-5069, USA
| | - János G. Ángyán
- CNRS, CRM2, UMR 7036,
Vandœuvre-lès-Nancy F-54506, France and Université de Lorraine, CRM2, UMR 7036, Vandœuvre-lès-Nancy F-54506,
France
| | - Giulia Galli
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637,
USA
| | - Cui Zhang
- Department of Chemistry, University of California Davis, Davis, California 95616,
USA
| | - Francois Gygi
- Department of Computer Science,
University of California, Davis, California 95616, USA
| | - Kimihiko Hirao
- Computational Chemistry Unit, RIKEN Advanced Institute for Computational Science, Kobe, Hyogo 6500047,
Japan
| | - Jong Won Song
- Computational Chemistry Unit, RIKEN Advanced Institute for Computational Science, Kobe, Hyogo 6500047,
Japan
| | - Kar Rahul
- Computational Chemistry Unit, RIKEN Advanced Institute for Computational Science, Kobe, Hyogo 6500047,
Japan
| | - O. Anatole von Lilienfeld
- General Chemistry (ALGC), Free University Brussels (VUB), Pleinlaan 2, 1050 Brussel,
Belgium and Institute of Physical Chemistry and National Center
for Computational Design and Discovery of Novel Materials (MARVEL), Department of
Chemistry, University of Basel, 4056 Basel,
Switzerland
| | - Rafał Podeszwa
- Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice, Poland
| | | | | | | | - Julien Toulouse
- Laboratoire de Chimie Théorique, Sorbonne Universités,
Université Pierre et Marie Curie, CNRS, F-75005 Paris,
France
| | - Roberto Peverati
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455,
USA
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901,
USA
| | - Donald G. Truhlar
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455,
USA
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy,
University of Delaware, Newark, Delaware 19716,
USA
| |
Collapse
|
26
|
Reilly AM, Cooper RI, Adjiman CS, Bhattacharya S, Boese AD, Brandenburg JG, Bygrave PJ, Bylsma R, Campbell JE, Car R, Case DH, Chadha R, Cole JC, Cosburn K, Cuppen HM, Curtis F, Day GM, DiStasio Jr RA, Dzyabchenko A, van Eijck BP, Elking DM, van den Ende JA, Facelli JC, Ferraro MB, Fusti-Molnar L, Gatsiou CA, Gee TS, de Gelder R, Ghiringhelli LM, Goto H, Grimme S, Guo R, Hofmann DWM, Hoja J, Hylton RK, Iuzzolino L, Jankiewicz W, de Jong DT, Kendrick J, de Klerk NJJ, Ko HY, Kuleshova LN, Li X, Lohani S, Leusen FJJ, Lund AM, Lv J, Ma Y, Marom N, Masunov AE, McCabe P, McMahon DP, Meekes H, Metz MP, Misquitta AJ, Mohamed S, Monserrat B, Needs RJ, Neumann MA, Nyman J, Obata S, Oberhofer H, Oganov AR, Orendt AM, Pagola GI, Pantelides CC, Pickard CJ, Podeszwa R, Price LS, Price SL, Pulido A, Read MG, Reuter K, Schneider E, Schober C, Shields GP, Singh P, Sugden IJ, Szalewicz K, Taylor CR, Tkatchenko A, Tuckerman ME, Vacarro F, Vasileiadis M, Vazquez-Mayagoitia A, Vogt L, Wang Y, Watson RE, de Wijs GA, Yang J, Zhu Q, Groom CR. Report on the sixth blind test of organic crystal structure prediction methods. Acta Crystallogr B Struct Sci Cryst Eng Mater 2016; 72:439-59. [PMID: 27484368 PMCID: PMC4971545 DOI: 10.1107/s2052520616007447] [Citation(s) in RCA: 294] [Impact Index Per Article: 36.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2016] [Accepted: 05/04/2016] [Indexed: 05/05/2023]
Abstract
The sixth blind test of organic crystal structure prediction (CSP) methods has been held, with five target systems: a small nearly rigid molecule, a polymorphic former drug candidate, a chloride salt hydrate, a co-crystal and a bulky flexible molecule. This blind test has seen substantial growth in the number of participants, with the broad range of prediction methods giving a unique insight into the state of the art in the field. Significant progress has been seen in treating flexible molecules, usage of hierarchical approaches to ranking structures, the application of density-functional approximations, and the establishment of new workflows and `best practices' for performing CSP calculations. All of the targets, apart from a single potentially disordered Z' = 2 polymorph of the drug candidate, were predicted by at least one submission. Despite many remaining challenges, it is clear that CSP methods are becoming more applicable to a wider range of real systems, including salts, hydrates and larger flexible molecules. The results also highlight the potential for CSP calculations to complement and augment experimental studies of organic solid forms.
Collapse
Affiliation(s)
- Anthony M. Reilly
- The Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, England
| | - Richard I. Cooper
- Chemical Crystallography, Chemistry Research Laboratory, Mansfield Road, Oxford OX1 3TA, England
| | - Claire S. Adjiman
- Department of Chemical Engineering, Centre for Process Systems Engineering, Imperial College London, London SW7 2AZ, England
| | - Saswata Bhattacharya
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195, Berlin, Germany
| | - A. Daniel Boese
- Department of Chemistry, Institute of Physical and Theoretical Chemistry, University of Graz, Heinrichstraße 28/IV, 8010 Graz, Austria
| | - Jan Gerit Brandenburg
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Peter J. Bygrave
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Rita Bylsma
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Josh E. Campbell
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, NJ 08544, USA
| | - David H. Case
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Renu Chadha
- University Institute of Pharmaceutical Sciences, Panjab University, Chandigarh, India
| | - Jason C. Cole
- The Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, England
| | - Katherine Cosburn
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA
- Department of Physics, University of Toronto, Toronto, Canada M5S 1A7
| | - Herma M. Cuppen
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Farren Curtis
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA
- Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Graeme M. Day
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Robert A. DiStasio Jr
- Department of Chemistry, Princeton University, Princeton, NJ 08544, USA
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, NY 14853, USA
| | | | | | - Dennis M. Elking
- OpenEye Scientific Software, 9 Bisbee Court, Suite D, Santa Fe, NM 87508, USA
| | - Joost A. van den Ende
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Julio C. Facelli
- Center for High Performance Computing, University of Utah, 155 South 1452 East Room 405, Salt Lake City, UT 84112-0190, USA
- Department of Biomedical Informatics, University of Utah, 155 South 1452 East Room 405, Salt Lake City, UT 84112-0190, USA
| | - Marta B. Ferraro
- Departamento de Física and Ifiba (CONICET) Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. I (1428), Buenos Aires, Argentina
| | - Laszlo Fusti-Molnar
- OpenEye Scientific Software, 9 Bisbee Court, Suite D, Santa Fe, NM 87508, USA
| | - Christina-Anna Gatsiou
- Department of Chemical Engineering, Centre for Process Systems Engineering, Imperial College London, London SW7 2AZ, England
| | - Thomas S. Gee
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - René de Gelder
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Luca M. Ghiringhelli
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195, Berlin, Germany
| | - Hitoshi Goto
- Educational Programs on Advanced Simulation Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi 441-8580, Japan
- Department of Computer Science and Engineering, Graduate School of Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi 441-8580, Japan
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Rui Guo
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, England
| | - Detlef W. M. Hofmann
- CRS4, Parco Scientifico e Tecnologico, POLARIS, Edificio 1, 09010 PULA, Italy
- FlexCryst, Schleifweg 23, 91080 Uttenreuth, Germany
| | - Johannes Hoja
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195, Berlin, Germany
| | - Rebecca K. Hylton
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, England
| | - Luca Iuzzolino
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, England
| | - Wojciech Jankiewicz
- Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice, Poland
| | - Daniël T. de Jong
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - John Kendrick
- Faculty of Life Sciences, University of Bradford, Richmond Road, Bradford BD7 1DP, England
| | - Niek J. J. de Klerk
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Hsin-Yu Ko
- Department of Chemistry, Princeton University, Princeton, NJ 08544, USA
| | | | - Xiayue Li
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA
- Argonne Leadership Computing Facility, Argonne National Laboratory, Lemont, IL 60439, USA
| | - Sanjaya Lohani
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA
| | - Frank J. J. Leusen
- Faculty of Life Sciences, University of Bradford, Richmond Road, Bradford BD7 1DP, England
| | - Albert M. Lund
- OpenEye Scientific Software, 9 Bisbee Court, Suite D, Santa Fe, NM 87508, USA
- Department of Chemistry, University of Utah, 155 South 1452 East Room 405, Salt Lake City, UT 84112-0190, USA
| | - Jian Lv
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, People’s Republic of China
| | - Yanming Ma
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, People’s Republic of China
| | - Noa Marom
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA
- Department of Materials Science and Engineering and Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Artëm E. Masunov
- NanoScience Technology Center, University of Central Florida, 12424 Research Parkway, PAV400, Orlando, FL 32826, USA
- Department of Chemistry, University of Central Florida, 4111 Libra Drive PSB225, Orlando, FL 32816, USA
- Department of Physics, University of Central Florida, 4111 Libra Drive PSB430, Orlando, FL 32816, USA
- Department of Condensed Matter Physics, National Research Nuclear University MEPhI, Kashirskoye shosse 31, Moscow 115409, Russia
| | - Patrick McCabe
- The Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, England
| | - David P. McMahon
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Hugo Meekes
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Michael P. Metz
- Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, USA
| | - Alston J. Misquitta
- School of Physics and Astronomy, Queen Mary University of London, London E1 4NS, England
| | | | - Bartomeu Monserrat
- Cavendish Laboratory, 19, J. J. Thomson Avenue, Cambridge CB3 0HE, England
- Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08854-8019, USA
| | - Richard J. Needs
- Cavendish Laboratory, 19, J. J. Thomson Avenue, Cambridge CB3 0HE, England
| | | | - Jonas Nyman
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Shigeaki Obata
- Educational Programs on Advanced Simulation Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi 441-8580, Japan
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Artem R. Oganov
- Department of Geosciences, Center for Materials by Design, and Institute for Advanced Computational Science, SUNY Stony Brook, NY 11794-2100, USA
- Skolkovo Institute of Science and Technology, Skolkovo Innovation Centers, Bldg. 3, Moscow Region, 143026, Russia
- Moscow Institute of Physics and Technology, 9 Institutskiy Lane, Dolgoprudny City, Moscow Region 141700, Russia
- International Center for Materials Discovery, School of Materials Science and Engineering, Northwestern Polytechnical University, Xi’an 710072, China
| | - Anita M. Orendt
- Center for High Performance Computing, University of Utah, 155 South 1452 East Room 405, Salt Lake City, UT 84112-0190, USA
| | - Gabriel I. Pagola
- Departamento de Física and Ifiba (CONICET) Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. I (1428), Buenos Aires, Argentina
| | - Constantinos C. Pantelides
- Department of Chemical Engineering, Centre for Process Systems Engineering, Imperial College London, London SW7 2AZ, England
| | - Chris J. Pickard
- Department of Materials Science and Metallurgy, University of Cambridge, 27 Charles Babbage Road, Cambridge CB3 0FS, England
- Department of Physics and Astronomy, University College London, Gower St., London WC1E 6BT, England
| | - Rafal Podeszwa
- Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice, Poland
| | - Louise S. Price
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, England
| | - Sarah L. Price
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, England
| | - Angeles Pulido
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Murray G. Read
- The Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, England
| | - Karsten Reuter
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Elia Schneider
- Department of Chemistry, New York University, New York, NY 10003, USA
| | - Christoph Schober
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Gregory P. Shields
- The Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, England
| | - Pawanpreet Singh
- University Institute of Pharmaceutical Sciences, Panjab University, Chandigarh, India
| | - Isaac J. Sugden
- Department of Chemical Engineering, Centre for Process Systems Engineering, Imperial College London, London SW7 2AZ, England
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, USA
| | | | - Alexandre Tkatchenko
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195, Berlin, Germany
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg
| | - Mark E. Tuckerman
- Department of Chemistry, New York University, New York, NY 10003, USA
- Courant Institute of Mathematical Sciences, New York University, New York, NY 10012, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
| | - Francesca Vacarro
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA
- Department of Chemistry, Loyola University, New Orleans, LA 70118, USA
| | - Manolis Vasileiadis
- Department of Chemical Engineering, Centre for Process Systems Engineering, Imperial College London, London SW7 2AZ, England
| | | | - Leslie Vogt
- Department of Chemistry, New York University, New York, NY 10003, USA
| | - Yanchao Wang
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, People’s Republic of China
| | - Rona E. Watson
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, England
| | - Gilles A. de Wijs
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Jack Yang
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Qiang Zhu
- Department of Geosciences, Center for Materials by Design, and Institute for Advanced Computational Science, SUNY Stony Brook, NY 11794-2100, USA
| | - Colin R. Groom
- The Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, England
| |
Collapse
|
27
|
Loboda O, Ingrosso F, Ruiz-López MF, Szalewicz K, Millot C. Geometry-dependent distributed polarizability models for the water molecule. J Chem Phys 2016; 144:034304. [PMID: 26801031 DOI: 10.1063/1.4939519] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Geometry-dependent distributed polarizability models have been constructed by fits to ab initio calculations at the coupled cluster level of theory with up to noniterative triple excitations in an augmented triple-zeta quality basis set for the water molecule in the field of a point charge. The investigated models include (i) charge-flow polarizabilities between chemically bonded atoms, (ii) isotropic or anisotropic dipolar polarizabilities on oxygen atom or on all atoms, and (iii) combinations of models (i) and (ii). For each model, the polarizability parameters have been optimized to reproduce the induction energy of a water molecule polarized by a point charge successively occupying a grid of points surrounding the molecule. The quality of the models is ascertained by examining their ability to reproduce these induction energies as well as the molecular dipolar and quadrupolar polarizabilities. The geometry dependence of the distributed polarizability models has been explored by changing bond lengths and HOH angle to generate 125 molecular structures (reduced to 75 symmetry-unique ones). For each considered model, the distributed polarizability components have been fitted as a function of the geometry by a Taylor expansion in monomer coordinate displacements up to the sum of powers equal to 4.
Collapse
Affiliation(s)
- Oleksandr Loboda
- Université de Lorraine, SRSMC UMR 7565, Vandoeuvre-les-Nancy F-54506, France
| | - Francesca Ingrosso
- Université de Lorraine, SRSMC UMR 7565, Vandoeuvre-les-Nancy F-54506, France
| | - Manuel F Ruiz-López
- Université de Lorraine, SRSMC UMR 7565, Vandoeuvre-les-Nancy F-54506, France
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Claude Millot
- Université de Lorraine, SRSMC UMR 7565, Vandoeuvre-les-Nancy F-54506, France
| |
Collapse
|
28
|
Podeszwa R, Bukowski R, Szalewicz K. Density-Fitting Method in Symmetry-Adapted Perturbation Theory Based on Kohn-Sham Description of Monomers. J Chem Theory Comput 2015; 2:400-12. [PMID: 26626527 DOI: 10.1021/ct050304h] [Citation(s) in RCA: 144] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We present a new implementation of symmetry-adapted perturbation theory of intermolecular interactions based on Kohn-Sham description of monomers. With density-fitting of molecular integrals, the scaling of the computational cost of the method is reduced from the sixth to the fifth power of the system size. Computational requirements of some operations scaling as the fifth power have also been significantly reduced. The new method allows an accurate treatment of molecules consisting of as many as a few dozen of atoms, using both nonhybrid and hybrid density functionals.
Collapse
Affiliation(s)
- Rafał Podeszwa
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716
| | - Robert Bukowski
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716
| |
Collapse
|
29
|
Piszczatowski K, Puchalski M, Komasa J, Jeziorski B, Szalewicz K. Frequency-dependent polarizability of helium including relativistic effects with nuclear recoil terms. Phys Rev Lett 2015; 114:173004. [PMID: 25978230 DOI: 10.1103/physrevlett.114.173004] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/04/2015] [Indexed: 06/04/2023]
Abstract
Future metrology standards will be partly based on physical quantities computed from first principles rather than measured. In particular, a new pressure standard can be established if the dynamic polarizability of helium can be determined from theory with an uncertainty smaller than 0.2 ppm. We present calculations of the frequency-dependent part of this quantity including relativistic effects with full account of leading nuclear recoil terms and using highly optimized explicitly correlated basis sets. A particular emphasis is put on uncertainty estimates. At the He-Ne laser wavelength of 632.9908 nm, the computed polarizability value of 1.39181141 a.u. has uncertainty of 0.1 ppm that is 2 orders of magnitude smaller than those of the most accurate polarizability measurements. We also obtained an accurate expansion of the helium refractive index in powers of density.
Collapse
Affiliation(s)
- Konrad Piszczatowski
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Mariusz Puchalski
- Faculty of Chemistry, Adam Mickiewicz University, Umultowska 89B, 61-614 Poznań, Poland
| | - Jacek Komasa
- Faculty of Chemistry, Adam Mickiewicz University, Umultowska 89B, 61-614 Poznań, Poland
| | - Bogumił Jeziorski
- Department of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
30
|
Jankowski P, Murdachaew G, Bukowski R, Akin-Ojo O, Leforestier C, Szalewicz K. Ab initio water pair potential with flexible monomers. J Phys Chem A 2015; 119:2940-64. [PMID: 25687650 DOI: 10.1021/jp512847z] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A potential energy surface for the water dimer with explicit dependence on monomer coordinates is presented. The surface was fitted to a set of previously published interaction energies computed on a grid of over a quarter million points in the 12-dimensional configurational space using symmetry-adapted perturbation theory and coupled-cluster methods. The present fit removes small errors in published fits, and its accuracy is critically evaluated. The minimum and saddle-point structures of the potential surface were found to be very close to predictions from direct ab initio optimizations. The computed second virial coefficients agreed well with experimental values. At low temperatures, the effects of monomer flexibility in the virial coefficients were found to be much smaller than the quantum effects.
Collapse
Affiliation(s)
- Piotr Jankowski
- †Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States.,§Faculty of Chemistry, Nicolaus Copernicus University in Torun, Gagarina 7, 87-100 Torun, Poland
| | - Garold Murdachaew
- †Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Robert Bukowski
- †Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Omololu Akin-Ojo
- †Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Claude Leforestier
- ‡Institut Charles Gerhardt (CTMM)-UMR 5253, CC 15.01, Université Montpellier II-CNRS, Montpellier, Cedex 05, France 34095
| | - Krzysztof Szalewicz
- †Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
31
|
Abstract
CONSPECTUS: Until recently, it had been impossible to predict structures of molecular crystals just from the knowledge of the chemical formula for the constituent molecule(s). A solution of this problem has been achieved using intermolecular force fields computed from first principles. These fields were developed by calculating interaction energies of molecular dimers and trimers using an ab initio method called symmetry-adapted perturbation theory (SAPT) based on density-functional theory (DFT) description of monomers [SAPT(DFT)]. For clusters containing up to a dozen or so atoms, interaction energies computed using SAPT(DFT) are comparable in accuracy to the results of the best wave function-based methods, whereas the former approach can be applied to systems an order of magnitude larger than the latter. In fact, for monomers with a couple dozen atoms, SAPT(DFT) is about equally time-consuming as the supermolecular DFT approach. To develop a force field, SAPT(DFT) calculations are performed for a large number of dimer and possibly also trimer configurations (grid points in intermolecular coordinates), and the interaction energies are then fitted by analytic functions. The resulting force fields can be used to determine crystal structures and properties by applying them in molecular packing, lattice energy minimization, and molecular dynamics calculations. In this way, some of the first successful determinations of crystal structures were achieved from first principles, with crystal densities and lattice parameters agreeing with experimental values to within about 1%. Crystal properties obtained using similar procedures but empirical force fields fitted to crystal data have typical errors of several percent due to low sensitivity of empirical fits to interactions beyond those of the nearest neighbors. The first-principles approach has additional advantages over the empirical approach for notional crystals and cocrystals since empirical force fields can only be extrapolated to such cases. As an alternative to applying SAPT(DFT) in crystal structure calculations, one can use supermolecular DFT interaction energies combined with scaled dispersion energies computed from simple atom-atom functions, that is, use the so-called DFT+D approach. Whereas the standard DFT methods fail for intermolecular interactions, DFT+D performs reasonably well since the dispersion correction is used not only to provide the missing dispersion contribution but also to fix other deficiencies of DFT. The latter cancellation of errors is unphysical and can be avoided by applying the so-called dispersionless density functional, dlDF. In this case, the dispersion energies are added without any scaling. The dlDF+D method is also one of the best performing DFT+D methods. The SAPT(DFT)-based approach has been applied so far only to crystals with rigid monomers. It can be extended to partly flexible monomers, that is, to monomers with only a few internal coordinates allowed to vary. However, the costs will increase relative to rigid monomer cases since the number of grid points increases exponentially with the number of dimensions. One way around this problem is to construct force fields with approximate couplings between inter- and intramonomer degrees of freedom. Another way is to calculate interaction energies (and possibly forces) "on the fly", i.e., in each step of lattice energy minimization procedure. Such an approach would be prohibitively expensive if it replaced analytic force fields at all stages of the crystal predictions procedure, but it can be used to optimize a few dozen candidate structures determined by other methods.
Collapse
Affiliation(s)
- Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
32
|
Góra U, Cencek W, Podeszwa R, van der Avoird A, Szalewicz K. Predictions for water clusters from a first-principles two- and three-body force field. J Chem Phys 2014; 140:194101. [PMID: 24852524 DOI: 10.1063/1.4875097] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new rigid-monomer three-body potential has been developed for water by fitting it to more than 70 thousand trimer interaction energies computed ab initio using coupled-cluster methods and augmented triple-zeta-quality basis sets. This potential was used together with a modified form of a previously developed two-body potential and with a polarization model of four- and higher-body interactions to predict the energetics of the water trimer, hexamer, and 24-mer. Despite using the rigid-monomer approximation, these predictions agree better with flexible-monomer benchmarks than published results obtained with flexible-monomer force fields. An unexpected finding of our work is that simple polarization models predict four-body interactions to within a few percent, whereas for three-body interactions these models are known to have errors on the order of 50%.
Collapse
Affiliation(s)
- Urszula Góra
- Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice, Poland
| | - Wojciech Cencek
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Rafał Podeszwa
- Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice, Poland
| | - Ad van der Avoird
- Theoretical Chemistry, Institute for Molecules and Materials, Radboud University Nijmegen, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
33
|
Garberoglio G, Jankowski P, Szalewicz K, Harvey AH. Path-integral calculation of the second virial coefficient including intramolecular flexibility effects. J Chem Phys 2014; 141:044119. [DOI: 10.1063/1.4890440] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Giovanni Garberoglio
- Interdisciplinary Laboratory for Computational Science (LISC), FBK-CMM and University of Trento, via Sommarive 18, I-38123 Povo (TN), Italy
| | - Piotr Jankowski
- Department of Quantum Chemistry, Faculty of Chemistry, Nicolaus Copernicus University, Gagarina 7, PL-87-100 Toruń, Poland
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Allan H. Harvey
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, 325 Broadway, Boulder, Colorado 80305-3337, USA
| |
Collapse
|
34
|
van der Avoird A, Podeszwa R, Ensing B, Szalewicz K. Comment on “Communication: Benzene dimer—The free energy landscape” [J. Chem. Phys. 139, 201102 (2013)]. J Chem Phys 2014; 140:227101. [DOI: 10.1063/1.4882015] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
35
|
Cencek W, Szalewicz K. Erratum: “On asymptotic behavior of density functional theory” [J. Chem. Phys. 139, 024104 (2013)]. J Chem Phys 2014. [DOI: 10.1063/1.4870418] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
36
|
Rob F, Misquitta AJ, Podeszwa R, Szalewicz K. Localized overlap algorithm for unexpanded dispersion energies. J Chem Phys 2014; 140:114304. [DOI: 10.1063/1.4867969] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
37
|
Hoja J, Sax AF, Szalewicz K. Is Electrostatics Sufficient to Describe Hydrogen-Bonding Interactions? Chemistry 2014; 20:2292-300. [DOI: 10.1002/chem.201303528] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2013] [Revised: 11/25/2013] [Indexed: 11/10/2022]
|
38
|
Affiliation(s)
- Fazle Rob
- a Department of Physics and Astronomy , University of Delaware , Newark , DE , 19716 , USA
| | - Krzysztof Szalewicz
- a Department of Physics and Astronomy , University of Delaware , Newark , DE , 19716 , USA
| |
Collapse
|
39
|
|
40
|
Cencek W, Garberoglio G, Harvey AH, McLinden MO, Szalewicz K. Three-Body Nonadditive Potential for Argon with Estimated Uncertainties and Third Virial Coefficient. J Phys Chem A 2013; 117:7542-52. [DOI: 10.1021/jp4018579] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Wojciech Cencek
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Giovanni Garberoglio
- Interdisciplinary Laboratory for Computational Science
(LISC), FBK-CMM and University of Trento, via Sommarive 18, I-38123 Povo (TN), Italy
| | - Allan H. Harvey
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, 325
Broadway, Boulder, Colorado 80305, United States
| | - Mark O. McLinden
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, 325
Broadway, Boulder, Colorado 80305, United States
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
41
|
|
42
|
Jankowski P, Surin LA, Potapov A, Schlemmer S, McKellar ARW, Szalewicz K. A comprehensive experimental and theoretical study of H2−CO spectra. J Chem Phys 2013; 138:084307. [PMID: 23464151 DOI: 10.1063/1.4791712] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Piotr Jankowski
- Department of Quantum Chemistry, Faculty of Chemistry, Nicolaus Copernicus University, Gagarina 7, 87-100 Toruń, Poland
| | | | | | | | | | | |
Collapse
|
43
|
Affiliation(s)
- Omololu Akin-Ojo
- Theoretical Physics Department, African University of Science and Technology, Km 10 Airport Road, Galadimawa, Abuja, Nigeria
| | | |
Collapse
|
44
|
Abstract
We report the definition and testing of a new ab initio 12-dimensional potential for the water dimer with flexible monomers. Using our recent accurate CCpol-8s rigid water pair potential [W. Cencek, K. Szalewicz, C. Leforestier, R. van Harrevelt, and A. van der Avoird, Phys. Chem. Chem. Phys. 10, 4716 (2008)] as a reference for the undistorted monomers' geometries, a distortion correction has been added, which was taken from a former flexible-monomer ab initio potential. This correction allows us to retrieve the correct binding energy D(e)=21.0 kJ mol(-1), and leads to an equilibrium geometry in close agreement with the one obtained from benchmark calculations. The kinetic energy operator describing the flexible-monomer water dimer has been expressed in terms of Radau coordinates for each monomer and a recent general cluster polyspherical formulation describing their relative motions. Within this formulation, an adiabatic scheme has been invoked in order to decouple fast (intramolecular) modes and slow (intermolecular) ones. Different levels of approximation were tested, which differ in the way in which the residual potential coupling between the intramolecular modes located on different monomers and the dependence of the monomer rotational constants on the dimer geometry are handled. Accurate calculations of the vibration-rotation-tunneling levels of (H(2)O)(2) and (D(2)O)(2) were performed, which show the best agreement with experiments achieved so far for any water potential. Intramolecular excitations of the two monomers were calculated within two limiting cases, to account for the lack of non-adiabatic coupling between intramolecular modes due to the intermolecular motion. In the first model, the excitation was assumed to stay either on the donor or the acceptor molecule, and to hop between the two moieties upon donor-acceptor interchange. In the second model, the excitation remains on the same molecule whatever is the dimer geometry. Marginal frequency differences, less than 2 cm(-1), were obtained for all modes, and the resulting infrared shifts are in good agreement with experiments.
Collapse
Affiliation(s)
- Claude Leforestier
- Institut Charles Gerhardt (CTMM)-UMR 5253, CC 1501, Université Montpellier II-CNRS, 34095 Montpellier Cedex 05, France.
| | | | | |
Collapse
|
45
|
Garberoglio G, Jankowski P, Szalewicz K, Harvey AH. Second virial coefficients of H2 and its isotopologues from a six-dimensional potential. J Chem Phys 2012; 137:154308. [DOI: 10.1063/1.4757565] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
46
|
Cencek W, Przybytek M, Komasa J, Mehl JB, Jeziorski B, Szalewicz K. Effects of adiabatic, relativistic, and quantum electrodynamics interactions on the pair potential and thermophysical properties of helium. J Chem Phys 2012; 136:224303. [DOI: 10.1063/1.4712218] [Citation(s) in RCA: 209] [Impact Index Per Article: 17.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
|
47
|
Jankowski P, McKellar ARW, Szalewicz K. Theory Untangles the High-Resolution Infrared Spectrum of the
ortho
-H
2
-CO van der Waals Complex. Science 2012; 336:1147-50. [DOI: 10.1126/science.1221000] [Citation(s) in RCA: 68] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Affiliation(s)
- Piotr Jankowski
- Faculty of Chemistry, Nicolaus Copernicus University, Gagarina 7, 87-100 Toruń, Poland
| | - A. R. W. McKellar
- Steacie Institute for Molecular Sciences, National Research Council, Ottawa, Ontario K1A 0R6, Canada
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, USA
| |
Collapse
|
48
|
Przybytek M, Jeziorski B, Cencek W, Komasa J, Mehl JB, Szalewicz K. Onset of Casimir-Polder retardation in a long-range molecular quantum state. Phys Rev Lett 2012; 108:183201. [PMID: 22681072 DOI: 10.1103/physrevlett.108.183201] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/18/2012] [Indexed: 06/01/2023]
Abstract
Two 4He atoms form a diatomic molecule with a significant vibrational wave function amplitude at interatomic separations R>100 Å, where the retardation switches the London R(-6) decay of the potential to the Casimir-Polder R(-7) form. It has been assumed that this effect of retardation on the long-range part of the potential is responsible for the 2 Å (4%) increase of the bond length <R> of 4He2. We show that <R> is, unexpectedly, insensitive to the potential at R>20 Å and its increase is due to quantum electrodynamics effects computed by us from expressions valid at short R--beyond the validity range of Casimir-Polder theory--that seamlessly extend this theory to distances relevant for properties of long molecules.
Collapse
Affiliation(s)
- Michał Przybytek
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | | | | | | | | | | |
Collapse
|
49
|
Podeszwa R, Cencek W, Szalewicz K. Efficient Calculations of Dispersion Energies for Nanoscale Systems from Coupled Density Response Functions. J Chem Theory Comput 2012; 8:1963-9. [DOI: 10.1021/ct300200m] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rafał Podeszwa
- Institute of Chemistry, University
of Silesia, Szkolna 9, 40-006 Katowice, Poland
- Department of Physics and Astronomy,
University of Delaware, Newark, Delaware 19716, United States
| | - Wojciech Cencek
- Department of Physics and Astronomy,
University of Delaware, Newark, Delaware 19716, United States
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy,
University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
50
|
Podeszwa R, Szalewicz K. Communication: Density functional theory overcomes the failure of predicting intermolecular interaction energies. J Chem Phys 2012; 136:161102. [DOI: 10.1063/1.4707166] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Rafał Podeszwa
- Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice, Poland
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Krzysztof Szalewicz
- Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice, Poland
| |
Collapse
|