1
|
Kříž K, van Maaren PJ, van der Spoel D. Impact of Combination Rules, Level of Theory, and Potential Function on the Modeling of Gas- and Condensed-Phase Properties of Noble Gases. J Chem Theory Comput 2024; 20:2362-2376. [PMID: 38477573 DOI: 10.1021/acs.jctc.3c01257] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/14/2024]
Abstract
The systems of noble gases are particularly instructive for molecular modeling due to the elemental nature of their interactions. They do not normally form bonds nor possess a (permanent) dipole moment, and the only forces determining their bonding/clustering stems from van der Waals forces─dispersion and Pauli repulsion, which can be modeled by empirical potential functions. Combination rules, that is, formulas to derive parameters for pair potentials of heterodimers from parameters of corresponding homodimers, have been studied at length for the Lennard-Jones 12-6 potentials but not in great detail for other, more accurate, potentials. In this work, we examine the usefulness of nine empirical potentials in their ability to reproduce quantum mechanical (QM) benchmark dissociation curves of noble gas dimers (He, Ne, Ar, Kr, and Xe homo- and heterodimers), and we systematically study the efficacy of different permutations of combination relations for each parameter of the potentials. Our QM benchmark comprises dissociation curves computed by several different coupled cluster implementations as well as symmetry-adapted perturbation theory. The two-parameter Lennard-Jones potentials were decisively outperformed by more elaborate potentials that sport a 25-30 times lower root-mean-square error (RMSE) when fitted to QM dissociation curves. Very good fits to the QM dissociation curves can be achieved with relatively inexpensive four- or even three-parameter potentials, for instance, the damped 14-7 potential (Halgren, J. Am. Chem. Soc. 1992, 114, 7827-7843), a four-parameter Buckingham potential (Werhahn et al., Chem. Phys. Lett. 2015, 619, 133-138), or the three-parameter Morse potential (Morse, Phys. Rev. 1929, 34, 57-64). Potentials for heterodimers that are generated from combination rules have an RMSE that is up to 20 times higher than potentials that are directly fitted to the QM dissociation curves. This means that the RMSE, in particular, for light atoms, is comparable in magnitude to the well-depth of the potential. Based on a systematic permutation of combination rules, we present one or more combination rules for each potential tested that yield a relatively low RMSE. Two new combination rules are introduced that perform well, one for the van der Waals radius σij as ( 1 2 ( σ i 3 + σ j 3 ) ) 1 / 3 and one for the well-depth ϵij as ( 1 2 ( ϵ i - 2 + ϵ j - 2 ) ) - 1 / 2 . The QM data and the fitted potentials were evaluated in the gas phase against experimental second virial coefficients for homo- and heterodimers, the latter of which allowed evaluation of the combination rules. The fitted models were used to perform condensed phase molecular dynamics simulations to verify the melting points, liquid densities at the melting point, and the enthalpies of vaporization produced by the models for pure substances. Subtle differences in the benchmark potentials, in particular, the well-depth, due to the level of theory used were found here to have a profound effect on the macroscopic properties of noble gases: second virial coefficients or the bulk properties in simulations. By explicitly including three-body dispersion in molecular simulations employing the best pair potential, we were able to obtain accurate melting points as well as satisfactory densities and enthalpies of vaporization.
Collapse
Affiliation(s)
- Kristian Kříž
- Department of Cell and Molecular Biology, Uppsala University, Box 596, Uppsala SE-75124, Sweden
| | - Paul J van Maaren
- Department of Cell and Molecular Biology, Uppsala University, Box 596, Uppsala SE-75124, Sweden
| | - David van der Spoel
- Department of Cell and Molecular Biology, Uppsala University, Box 596, Uppsala SE-75124, Sweden
| |
Collapse
|
2
|
Ströker P, Meier K. Vapor-liquid equilibrium and thermodynamic properties of saturated argon and krypton from Monte Carlo simulations using ab initio potentials. J Chem Phys 2024; 160:094503. [PMID: 38426525 DOI: 10.1063/5.0196466] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2024] [Accepted: 01/29/2024] [Indexed: 03/02/2024] Open
Abstract
Vapor-liquid equilibria and thermodynamic properties of saturated argon and krypton were calculated by semi-classical Monte Carlo simulations with the NpT + test particle method using ab initio potentials for the two-body and nonadditive three-body interactions. The NpT + test particle method was extended to the calculation of second-order thermodynamic properties, such as the isochoric and isobaric heat capacities or the speed of sound, of the saturated liquid and vapor by using our recently developed approach for the systematic calculation of arbitrary thermodynamic properties in the isothermal-isobaric ensemble. Generally, the results for all simulated properties agree well with experimental data and the current reference equations of state for argon and krypton. In particular, the results for the vapor pressure and for the density and speed of sound of the saturated liquid and vapor agree with the most accurate experimental data for both noble gases almost within the uncertainty of these data, a level of agreement unprecedented for many-particle simulations. This study demonstrates that the vapor-liquid equilibrium and thermodynamic properties at saturation of a pure fluid can be predicted by Monte Carlo simulations with high accuracy when the intermolecular interactions are described by state-of-the-art ab initio pair and nonadditive three-body potentials and quantum effects are accounted for.
Collapse
Affiliation(s)
- Philipp Ströker
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| | - Karsten Meier
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| |
Collapse
|
3
|
Xie Y, Glick ZL, Sherrill CD. Assessment of three-body dispersion models against coupled-cluster benchmarks for crystalline benzene, carbon dioxide, and triazine. J Chem Phys 2023; 158:094110. [PMID: 36889937 DOI: 10.1063/5.0143712] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
To study the contribution of three-body dispersion to crystal lattice energies, we compute the three-body contributions to the lattice energies for crystalline benzene, carbon dioxide, and triazine using various computational methods. We show that these contributions converge quickly as the intermolecular distances between the monomers grow. In particular, the smallest value among the three pairwise intermonomer closest-contact distances, Rmin, shows a strong correlation with the three-body contribution to the lattice energy, and, here, the largest of the closest-contact distances, Rmax, serves as a cutoff criterion to limit the number of trimers to be considered. We considered all trimers up to Rmax=15Å. The trimers with Rmin<4Å contribute 90.4%, 90.6%, and 93.9% of the total three-body contributions for crystalline benzene, carbon dioxide, and triazine, respectively, for the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] method. For trimers with Rmin>4Å, the second-order Møller-Plesset perturbation theory (MP2) supplemented with the Axilrod-Teller-Muto (ATM) three-body dispersion correction reproduces the CCSD(T) values for the cumulative three-body contributions with errors of less than 0.1 kJ mol-1. Moreover, three-body contributions are converged within 0.15 kJ mol-1 by Rmax=10Å. From these results, it appears that in molecular crystals where dispersion dominates the three-body contribution to the lattice energy, the trimers with Rmin>4Å can be computed with the MP2+ATM method to reduce the computational cost, and those with Rmax>10Å appear to be basically negligible.
Collapse
Affiliation(s)
- Yi Xie
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
4
|
Ströker P, Hellmann R, Meier K. Thermodynamic properties of krypton from Monte Carlo simulations using ab initio potentials. J Chem Phys 2022; 157:114504. [PMID: 36137797 DOI: 10.1063/5.0107851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Ten different thermodynamic properties of the noble gas krypton were calculated by Monte Carlo simulations in the isothermal-isobaric ensemble using a highly accurate ab initio pair potential, Feynman-Hibbs corrections for quantum effects, and an extended Axilrod-Teller-Muto potential to account for nonadditive three-body interactions. Fourteen state points at a liquid and a supercritical isotherm were simulated. To obtain results representative for macroscopic systems, simulations with several particle numbers were carried out and extrapolated to the thermodynamic limit. Our results agree well with experimental data from the literature, an accurate equation of state for krypton, and a recent virial equation of state (VEOS) for krypton in the region where the VEOS has converged. These results demonstrate that very good agreement between simulation and experiment can only be achieved if nonadditive three-body interactions and quantum effects are taken into account.
Collapse
Affiliation(s)
- Philipp Ströker
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| | - Robert Hellmann
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| | - Karsten Meier
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| |
Collapse
|
5
|
Ströker P, Hellmann R, Meier K. Thermodynamic properties of argon from Monte Carlo simulations using ab initio potentials. Phys Rev E 2022; 105:064129. [PMID: 35854585 DOI: 10.1103/physreve.105.064129] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2022] [Accepted: 05/16/2022] [Indexed: 06/15/2023]
Abstract
Ten different thermodynamic properties of the noble gas argon in the liquid and supercritical regions were obtained from semiclassical Monte Carlo simulations in the isothermal-isobaric ensemble using ab initio potentials for the two-body and nonadditive three-body interactions. Our results for the density and speed of sound agree with the most accurate experimental data for argon almost within the uncertainty of these data, a level of agreement unprecedented for many-particle simulations. This demonstrates the high predictive but yet unexploited power of ab initio potentials in the field of molecular modeling and simulation for thermodynamic properties of fluids.
Collapse
Affiliation(s)
- Philipp Ströker
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| | - Robert Hellmann
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| | - Karsten Meier
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| |
Collapse
|
6
|
Graham RS, Wheatley RJ. Machine learning for non-additive intermolecular potentials: quantum chemistry to first-principles predictions. Chem Commun (Camb) 2022; 58:6898-6901. [PMID: 35642644 DOI: 10.1039/d2cc01820a] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Prediction of thermophysical properties from molecular principles requires accurate potential energy surfaces (PES). We present a widely-applicable method to produce first-principles PES from quantum chemistry calculations. Our approach accurately interpolates three-body non-additive interaction data, using the machine learning technique, Gaussian Processes (GP). The GP approach needs no bespoke modification when the number or type of molecules is changed. Our method produces highly accurate interpolation from significantly fewer training points than typical approaches, meaning ab initio calculations can be performed at higher accuracy. As an exemplar we compute the PES for all three-body cross interactions for CO2-Ar mixtures. From these we calculate the CO2-Ar virial coefficients up to 5th order. The resulting virial equation of state (EoS) is convergent for densities up to the critical density. Where convergent, the EoS makes accurate first-principles predictions for a range of thermophysical properties for CO2-Ar mixtures, including the compressibility factor, speed of sound and Joule-Thomson coefficient. Our method has great potential to make wide-ranging first-principles predictions for mixtures of comparably sized molecules. Such predictions can replace the need for expensive, laborious and repetitive experiments and inform the continuum models required for applications.
Collapse
Affiliation(s)
- Richard S Graham
- School of Mathematical Sciences, University of Nottingham, Nottingham NG7 2RD, UK.
| | | |
Collapse
|
7
|
Nitzke I, Pohl S, Thol M, Span R, Vrabec J. How well does the Tang-Toennies potential represent the thermodynamic properties of argon? Mol Phys 2022. [DOI: 10.1080/00268976.2022.2078240] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Isabel Nitzke
- Thermodynamics and Process Engineering, Technische Universtität Berlin, Berlin, Germany
| | - Sven Pohl
- Thermodynamics, Ruhr-Universität Bochum, Bochum, Germany
| | - Monika Thol
- Thermodynamics, Ruhr-Universität Bochum, Bochum, Germany
| | - Roland Span
- Thermodynamics, Ruhr-Universität Bochum, Bochum, Germany
| | - Jadran Vrabec
- Thermodynamics and Process Engineering, Technische Universtität Berlin, Berlin, Germany
| |
Collapse
|
8
|
Alboul L, Lishchuk SV. Bulk viscosity of gaseous argon from molecular dynamics simulations. Phys Rev E 2022; 105:054135. [PMID: 35706273 DOI: 10.1103/physreve.105.054135] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2021] [Accepted: 05/06/2022] [Indexed: 06/15/2023]
Abstract
The bulk viscosity of dilute argon gas is calculated using molecular dynamics simulations in the temperature range 150-500 K and is found to be proportional to density squared in the investigated range of densities 0.001-1 kg m^{-3}. A comparison of the results obtained using Lennard-Jones and Tang-Toennies models of pair interaction potential reveals that the value of the bulk viscosity coefficient is sensitive to the choice of the pair interaction model. The inclusion of the Axilrod-Teller-Muto three-body interaction in the model does not noticeably affect the values of the bulk viscosity in dilute states, contrary to the previously investigated case of dense fluids.
Collapse
Affiliation(s)
- Lyuba Alboul
- Industry & Innovation Research Institute, Sheffield Hallam University, Sheffield S1 1WB, United Kingdom
| | - Sergey V Lishchuk
- Thermodynamics and Process Engineering, Technische Universität Berlin, 10587 Berlin, Germany
| |
Collapse
|
9
|
Abstract
A broad range of approaches to many-body dispersion are discussed, including empirical approaches with multiple fitted parameters, augmented density functional-based approaches, symmetry adapted perturbation theory, and a supermolecule approach based on coupled cluster theory. Differing definitions of "body" are considered, specifically atom-based vs molecule-based approaches.
Collapse
Affiliation(s)
- Peng Xu
- Department of Chemistry, Iowa State University, Ames, Iowa 50014, United States
| | - Melisa Alkan
- Department of Chemistry, Iowa State University, Ames, Iowa 50014, United States
| | - Mark S Gordon
- Department of Chemistry, Iowa State University, Ames, Iowa 50014, United States
| |
Collapse
|
10
|
Hartl B, Sharma S, Brügner O, Mertens SFL, Walter M, Kahl G. Reliable Computational Prediction of the Supramolecular Ordering of Complex Molecules under Electrochemical Conditions. J Chem Theory Comput 2020; 16:5227-5243. [PMID: 32536160 PMCID: PMC7426907 DOI: 10.1021/acs.jctc.9b01251] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2019] [Indexed: 11/28/2022]
Abstract
We propose a computationally lean, two-stage approach that reliably predicts self-assembly behavior of complex charged molecules on metallic surfaces under electrochemical conditions. Stage one uses ab initio simulations to provide reference data for the energies (evaluated for archetypical configurations) to fit the parameters of a conceptually much simpler and computationally less expensive force field of the molecules: classical, spherical particles, representing the respective atomic entities; a flat and perfectly conducting wall represents the metallic surface. Stage two feeds the energies that emerge from this force field into highly efficient and reliable optimization techniques to identify via energy minimization the ordered ground-state configurations of the molecules. We demonstrate the power of our approach by successfully reproducing, on a semiquantitative level, the intricate supramolecular ordering observed experimentally for PQP+ and ClO4- molecules at an Au(111)-electrolyte interface, including the formation of open-porous, self-host-guest, and stratified bilayer phases as a function of the electric field at the solid-liquid interface. We also discuss the role of the perchlorate ions in the self-assembly process, whose positions could not be identified in the related experimental investigations.
Collapse
Affiliation(s)
- Benedikt Hartl
- Institute
for Theoretical Physics and Center for Computational Materials Science
(CMS), TU Wien, 1040 Wien, Austria
| | - Shubham Sharma
- FIT
Freiburg Centre for Interactive Materials and Bioinspired Technologies, Georges-Köhler-Allee 105, 79110 Freiburg, Germany
| | - Oliver Brügner
- FIT
Freiburg Centre for Interactive Materials and Bioinspired Technologies, Georges-Köhler-Allee 105, 79110 Freiburg, Germany
| | - Stijn F. L. Mertens
- Department
of Chemistry, Lancaster University, Lancaster LA1 4YB, United Kingdom
- Institute
of Applied Physics, TU Wien, 1040 Wien, Austria
| | - Michael Walter
- FIT
Freiburg Centre for Interactive Materials and Bioinspired Technologies, Georges-Köhler-Allee 105, 79110 Freiburg, Germany
- Cluster
of Excellence livMatS@FIT—Freiburg Center for Interactive Materials
and Bioinspired Technologies, University
of Freiburg, Georges-Köhler-Allee 105, D-79110 Freiburg, Germany
- MikroTribologie
Centrum μTC, Fraunhofer IWM, Wöhlerstrasse 11, 79108 Freiburg, Germany
| | - Gerhard Kahl
- Institute
for Theoretical Physics and Center for Computational Materials Science
(CMS), TU Wien, 1040 Wien, Austria
| |
Collapse
|
11
|
Pattnaik P, Raghunathan S, Kalluri T, Bhimalapuram P, Jawahar CV, Priyakumar UD. Machine Learning for Accurate Force Calculations in Molecular Dynamics Simulations. J Phys Chem A 2020; 124:6954-6967. [DOI: 10.1021/acs.jpca.0c03926] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Punyaslok Pattnaik
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| | - Shampa Raghunathan
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| | - Tarun Kalluri
- Center for Visual Information Technology, KCIS, International Institute of Information Technology, Hyderabad 500 032, India
| | - Prabhakar Bhimalapuram
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| | - C. V. Jawahar
- Center for Visual Information Technology, KCIS, International Institute of Information Technology, Hyderabad 500 032, India
| | - U. Deva Priyakumar
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India
| |
Collapse
|
12
|
Márquez-Mijares M, Roncero O, Villarreal P, González-Lezana T. Theoretical methods for the rotation–vibration spectra of triatomic molecules: distributed Gaussian functions compared with hyperspherical coordinates. INT REV PHYS CHEM 2018. [DOI: 10.1080/0144235x.2018.1514187] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Maykel Márquez-Mijares
- Instituto Superior de Tecnologías y Ciencias Aplicadas InsTec, Havana University, Quinta de los Molinos, Havana, Cuba
| | - Octavio Roncero
- Instituto de Física Fundamental, Consejo Superior de Investigaciones Científicas IFF-CSIC, Madrid, Spain
| | - Pablo Villarreal
- Instituto de Física Fundamental, Consejo Superior de Investigaciones Científicas IFF-CSIC, Madrid, Spain
| | - Tomás González-Lezana
- Instituto de Física Fundamental, Consejo Superior de Investigaciones Científicas IFF-CSIC, Madrid, Spain
| |
Collapse
|
13
|
Pérez-Jiménez AJ, Brémond E, Adamo C, Sancho-García JC. Communication: Accurate description of interaction energies and three-body effects in weakly bound molecular complexes by PBE-QIDH models. J Chem Phys 2018; 149:041101. [PMID: 30068200 DOI: 10.1063/1.5042153] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We apply a recently developed parameter-free double-hybrid density functional belonging to the quadratic-integrand double-hybrid model to calculate association energies (ΔE) and three-body effects (Δ3E) arising from intermolecular interactions in weakly bound supramolecular complexes (i.e., the dataset 3B-69). The model behaves very accurately for trimer association energies and is found to outperform widely used density functional approximations while approaching the accuracy of more costly ab initio methods for three-body effects. The results are further improved when we add some specific corrections for the remaining dispersion interactions, D3(BJ) or VV10 for two-body effects and Axilrod-Teller-Muto for three-body effects, leading to marginal deviations (less than 1 kcal/mol for ΔE and around 0.03-0.04 kcal/mol for Δ3E) with respect to benchmark results.
Collapse
Affiliation(s)
- A J Pérez-Jiménez
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| | - E Brémond
- Université Paris Diderot, Sorbonne Paris Cité, ITODYS, UMR CNRS 7086, F-75013 Paris, France
| | - C Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris, F-75005 Paris, France
| | - J C Sancho-García
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| |
Collapse
|
14
|
Baumketner A, Stelmakh A, Cai W. Cluster Crystals Stabilized by Hydrophobic and Electrostatic Interactions. J Phys Chem B 2018; 122:2669-2682. [PMID: 29432018 DOI: 10.1021/acs.jpcb.7b11662] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Cluster crystals are crystalline materials in which each site is occupied by multiple identical particles, atoms, colloids, or polymers. There are two classes of systems that make cluster crystals. One is composed of particles that interact via potentials that are bound at the origin and thus are able to penetrate each other. The other consists of non-interpenetrating particles whose interaction potential diverges at the origin. The goal of this work is to find which systems of the second class can make cluster crystals that are stable at room temperature. First, the general properties of the required potentials are established using an analytical model and Monte Carlo simulations. Next, we ask how such potentials can be constructed by combining hydrophobic attraction and electrostatic repulsion. A colloid model with a hard-sphere core and a repulsive wall is introduced to mimic the hydrophobic interaction. Charge is added to create long-range repulsion. A search in the parameter space of the colloid size, counterion type, and charge configuration uncovers several models for which effective colloid-colloid interaction, determined in explicit solvent as a potential of mean force, has the necessary shape. For the effective potential, cluster crystals are confirmed as low free-energy configurations in replica-exchange molecular dynamics simulations, which also generate the respective transition temperature. The model that exhibits a transition above room temperature is further studied in explicit solvent. Simulations on a 10 ns time scale show that crystalline conformations are stable below the target temperature but disintegrate rapidly above it, supporting the idea that hydrophobic and electrostatic interactions are sufficient to induce an assembly of cluster crystals. Finally, we discuss which physical systems are good candidates for experimental observations of cluster crystals.
Collapse
Affiliation(s)
- A Baumketner
- Institute for Condensed Matter Physics , NAS of Ukraine , 1 Svientsistsky Street , Lviv 79011 , Ukraine
| | - A Stelmakh
- Department of Chemistry , Ivan Franko Lviv National University , 6 Kyrylo and Mefodii Street , Lviv 79005 , Ukraine
| | - W Cai
- Department of Mathematics , Southern Methodist University , Dallas , Texas 75252 , United States
| |
Collapse
|
15
|
Tashakor S, Noorbala MR, Payvandy P, Mohammadi-Manesh H, Namazian M. Introducing a novel method based on the imperialistic competitive algorithm to determine fluorine intermolecular potential from ab initio calculations and calculation of some properties via MD simulations. MOLECULAR SIMULATION 2017. [DOI: 10.1080/08927022.2017.1366655] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
| | | | - Pedram Payvandy
- Department of Textile Engineering, Yazd University, Yazd, Iran
| | | | | |
Collapse
|
16
|
Vlasiuk M, Sadus RJ. Predicting vapor-liquid phase equilibria with augmented ab initio interatomic potentials. J Chem Phys 2017; 146:244504. [DOI: 10.1063/1.4986917] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Maryna Vlasiuk
- Centre for Molecular Simulation, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| | - Richard J. Sadus
- Centre for Molecular Simulation, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| |
Collapse
|
17
|
Misquitta AJ, Stone AJ. Ab Initio Atom–Atom Potentials Using CamCASP: Theory and Application to Many-Body Models for the Pyridine Dimer. J Chem Theory Comput 2016; 12:4184-208. [DOI: 10.1021/acs.jctc.5b01241] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Alston J. Misquitta
- School
of Physics and Astronomy, Queen Mary, University of London, London E1 4NS, United Kingdom
| | - Anthony J. Stone
- University Chemical Laboratory, Lensfield Road, Cambridge, CB2 1EW, United Kingdom
| |
Collapse
|
18
|
Abstract
Symmetry-adapted perturbation theory (SAPT) provides a unique set of advantages for parameterizing next-generation force fields from first principles. SAPT provides a direct, basis-set superposition error free estimate of molecular interaction energies, a physically intuitive energy decomposition, and a seamless transition to an asymptotic picture of intermolecular interactions. These properties have been exploited throughout the literature to develop next-generation force fields for a variety of applications, including classical molecular dynamics simulations, crystal structure prediction, and quantum dynamics/spectroscopy. This review provides a brief overview of the formalism and theory of SAPT, along with a practical discussion of the various methodologies utilized to parameterize force fields from SAPT calculations. It also highlights a number of applications of SAPT-based force fields for chemical systems of particular interest. Finally, the review ends with a brief outlook on the future opportunities and challenges that remain for next-generation force fields based on SAPT.
Collapse
Affiliation(s)
- Jesse G McDaniel
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706;
| | - J R Schmidt
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706;
| |
Collapse
|
19
|
Pai SJ, Bae YC. Ab initio potential energy surface for methane and carbon dioxide and application to vapor-liquid coexistence. J Chem Phys 2014; 141:064303. [DOI: 10.1063/1.4891983] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
20
|
Goujon F, Malfreyt P, Tildesley DJ. The gas-liquid surface tension of argon: A reconciliation between experiment and simulation. J Chem Phys 2014; 140:244710. [DOI: 10.1063/1.4885351] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
|
21
|
McDaniel JG, Schmidt JR. First-Principles Many-Body Force Fields from the Gas Phase to Liquid: A “Universal” Approach. J Phys Chem B 2014; 118:8042-53. [DOI: 10.1021/jp501128w] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jesse G. McDaniel
- Theoretical
Chemistry Institute
and Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - J. R. Schmidt
- Theoretical
Chemistry Institute
and Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| |
Collapse
|
22
|
del Río F, Díaz-Herrera E, Guzmán O, Moreno-Razo JA, Ramos JE. Analytical equation of state with three-body forces: Application to noble gases. J Chem Phys 2013; 139:184503. [DOI: 10.1063/1.4829055] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
|
23
|
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] [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
|
24
|
Hatz R, Korpinen M, Hänninen V, Halonen L. Characterization of the Dispersion Interactions and an ab Initio Study of van der Waals Potential Energy Parameters for Coinage Metal Clusters. J Phys Chem A 2012; 116:11685-93. [DOI: 10.1021/jp307448n] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Richard Hatz
- Laboratory of Physical
Chemistry, Department
of Chemistry, University of Helsinki, P.O. Box 55 (A.I. Virtasen aukio 1), FIN-00014 University of Helsinki,
Finland
| | - Markus Korpinen
- Laboratory of Physical
Chemistry, Department
of Chemistry, University of Helsinki, P.O. Box 55 (A.I. Virtasen aukio 1), FIN-00014 University of Helsinki,
Finland
| | - Vesa Hänninen
- Laboratory of Physical
Chemistry, Department
of Chemistry, University of Helsinki, P.O. Box 55 (A.I. Virtasen aukio 1), FIN-00014 University of Helsinki,
Finland
| | - Lauri Halonen
- Laboratory of Physical
Chemistry, Department
of Chemistry, University of Helsinki, P.O. Box 55 (A.I. Virtasen aukio 1), FIN-00014 University of Helsinki,
Finland
| |
Collapse
|
25
|
Wiebke J, Pahl E, Schwerdtfeger P. Sensitivity of the thermal and acoustic virial coefficients of argon to the argon interaction potential. J Chem Phys 2012; 137:064702. [DOI: 10.1063/1.4740465] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
Wiebke J, Schwerdtfeger P, Moyano G, Pahl E. An atomistic fourth-order virial equation of state for Argon from first principles calculations. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.08.024] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
27
|
Szalewicz K. Symmetry-adapted perturbation theory of intermolecular forces. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.86] [Citation(s) in RCA: 357] [Impact Index Per Article: 27.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
28
|
Meuwly M, Doll JD. Finite-temperature quantum simulations of mixed rare gas clusters. J Chem Phys 2010; 132:234315. [PMID: 20572713 DOI: 10.1063/1.3431080] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
Finite-temperature quantum Monte Carlo simulations are presented for mixed neon/argon rare gas clusters containing up to n=10 atoms. For the smallest clusters (n=3) comparison with rigorous bound state calculations and experiments shows that the present approach is accurate to within fractions of wavenumbers for energies and to within a few percent or better for rotational constants. For larger cluster sizes, for which no rigorous quantum calculations are available, comparison with experiment becomes even more favorable. In all simulations accurate pair potentials for the rare gas-rare gas interactions are employed and comparison with high-level electronic structure calculations suggest that many-body interactions play a minor role. For the largest clusters investigated (Ne(4)Ar(6)) gradual melting of the neon phase is observed while the argon-phase remains structurally intact.
Collapse
Affiliation(s)
- Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | | |
Collapse
|
29
|
Non-additive interactions of nucleobases in model dinucleotide steps occurring in B-DNA crystals. J Mol Model 2010; 16:1721-9. [PMID: 20449619 DOI: 10.1007/s00894-010-0722-8] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2009] [Accepted: 03/26/2010] [Indexed: 10/19/2022]
Abstract
Non-additivity of base-base interactions in all ten possible model dinucleotide steps were analyzed on MP2/aug-cc-pvDZ quantum chemistry level. Conformations of four nucleobases exactly matched to ones occurring in B-DNA crystals. In most of thw 162 analyzed tetramers both three- and four-body contributions are negligible except for d(GpG) steps. However, in these dinucleotides both contributions are always of opposite signs and in all cases the sum of all non-additive part of intermolecular interactions do not exceed 2.6 kcal mol(-1). This stands for less than 5% of the overall binding energy of dinucleotide steps. Also replacements of guanine with 8-oxoguanine in d(GpG) systems introduces non-additivity of the same magnitude as for canonical dinucleotides. It is observed linear relationships between values of total binding energy obtained in the tetramer basis set and estimated energy exclusively in dimers basis sets with assumption of pairwise additivities. For all analyzed dinucleotides steps there are also linear correlations between amount of non-additive contributions and magnitude of pairs interactions. Based on differences in electrostatic contribution to the total binding energy of four nucleobases and polarity of dinucleotide steps three distinct classes of dinucleotide steps were identified.
Collapse
|
30
|
Pitonák M, Neogrády P, Hobza P. Three- and four-body nonadditivities in nucleic acid tetramers: a CCSD(T) study. Phys Chem Chem Phys 2009; 12:1369-78. [PMID: 20119615 DOI: 10.1039/b919354e] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Three- and four-body nonadditivities in the uracil tetramer (in DNA-like geometry) and the GC step (in crystal geometry) were investigated at various levels of the wave-function theory: HF, MP2, MP3, L-CCD, CCSD and CCSD(T). All of the calculations were performed using the 6-31G**(0.25,0.15) basis set, whereas the HF, MP2 and the MP3 nonadditivities were, for the sake of comparison, also determined with the much larger aug-cc-pVDZ basis set. The HF and MP2 levels do not provide reliable values for many-body terms, making it necessary to go beyond the MP2 level. The benchmark CCSD(T) three- and four-body nonadditivities are reasonably well reproduced at the MP3 level, and almost quantitative agreement is obtained (fortuitously) either on the L-CCD level or as an average of the MP3 and the CCSD results. Reliable values of many-body terms (especially their higher-order correlation contributions) are obtained already when the rather small 6-31G**(0.25,0.15) basis set is used. The four-body term is much smaller when compared to the three-body terms, but it is definitely not negligible, e.g. in the case of the GC step it represents about 16% of all of the three- and four-body terms. While investigating the geometry dependence of many-body terms for the GG step at the MP3/6-31G**(0.25,0.15) level, we found that it is necessary to include at least three-body terms in the determination of optimal geometry parameters.
Collapse
Affiliation(s)
- M Pitonák
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, v. v. i. and Center for Biomolecules and Complex Molecular Systems, Flemingovo nám. 2, 166 10 Praha 6, Czech Republic
| | | | | |
Collapse
|
31
|
Beran GJO. Approximating quantum many-body intermolecular interactions in molecular clusters using classical polarizable force fields. J Chem Phys 2009; 130:164115. [PMID: 19405569 DOI: 10.1063/1.3121323] [Citation(s) in RCA: 86] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Many-body intermolecular interaction expansions provide a promising avenue for the efficient quantum mechanical treatment of molecular clusters and condensed-phase systems, but the computationally expensive three-body and higher terms are often nontrivial. When polar molecules are involved, these many-body terms are typically dominated by electrostatic induction effects, which can be approximated relatively easily. We demonstrate an accurate and inexpensive hybrid quantum/classical model in which one- and two-body interactions are computed quantum mechanically, while the many-body induction effects are approximated with a simple classical polarizable force field. Whereas typical hybrid quantum/classical models partition a system spatially into distinct quantum and classical regions, the model demonstrated here partitions based on the order in the many-body interaction series. This enables a spatially homogeneous treatment of the entire system, which could prove advantageous in studying a wide range of condensed-phase molecular systems.
Collapse
Affiliation(s)
- Gregory J O Beran
- Department of Chemistry, University of California, Riverside, California 92521, USA.
| |
Collapse
|
32
|
Márquez-Mijares M, Pérez de Tudela R, González-Lezana T, Roncero O, Miret-Artés S, Delgado-Barrio G, Villarreal P, Baccarelli I, Gianturco FA, Rubayo-Soneira J. A theoretical investigation on the spectrum of the Ar trimer for high rotational excitations. J Chem Phys 2009; 130:154301. [DOI: 10.1063/1.3115100] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
33
|
Maerzke KA, Murdachaew G, Mundy CJ, Schenter GK, Siepmann JI. Self-Consistent Polarization Density Functional Theory: Application to Argon. J Phys Chem A 2009; 113:2075-85. [DOI: 10.1021/jp808767y] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Katie A. Maerzke
- Department of Chemistry and Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota 55455, and Chemical & Materials Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352
| | - Garold Murdachaew
- Department of Chemistry and Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota 55455, and Chemical & Materials Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352
| | - Christopher J. Mundy
- Department of Chemistry and Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota 55455, and Chemical & Materials Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352
| | - Gregory K. Schenter
- Department of Chemistry and Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota 55455, and Chemical & Materials Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352
| | - J. Ilja Siepmann
- Department of Chemistry and Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota 55455, and Chemical & Materials Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352
| |
Collapse
|
34
|
Oakley MT, Wheatley RJ. Additive and nonadditive models of vapor-liquid equilibrium in CO2 from first principles. J Chem Phys 2009; 130:034110. [DOI: 10.1063/1.3059008] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
35
|
Cencek W, Patkowski K, Szalewicz K. Full-configuration-interaction calculation of three-body nonadditive contribution to helium interaction potential. J Chem Phys 2009; 131:064105. [DOI: 10.1063/1.3204319] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
36
|
Eskandari Nasrabad A. Theory and atomistic simulation of krypton fluid. J Chem Phys 2008; 129:244504. [DOI: 10.1063/1.3046564] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
37
|
Pahl E, Calvo F, Koči L, Schwerdtfeger P. Accurate Melting Temperatures for Neon and Argon from Ab Initio Monte Carlo Simulations. Angew Chem Int Ed Engl 2008; 47:8207-10. [DOI: 10.1002/anie.200802743] [Citation(s) in RCA: 58] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
38
|
Genaue Schmelztemperaturen für Neon und Argon aus Ab-initio-Monte-Carlo-Simulationen. Angew Chem Int Ed Engl 2008. [DOI: 10.1002/ange.200802743] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
39
|
Bukowski R, Szalewicz K, Groenenboom GC, van der Avoird A. Polarizable interaction potential for water from coupled cluster calculations. II. Applications to dimer spectra, virial coefficients, and simulations of liquid water. J Chem Phys 2008; 128:094314. [PMID: 18331100 DOI: 10.1063/1.2832858] [Citation(s) in RCA: 78] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The six-dimensional CC-pol interaction potential for the water dimer was used to predict properties of the dimer and of liquid water, in the latter case after being supplemented by a nonadditive potential. All the results were obtained purely from first principles, i.e., without any fitting to experimental data. Calculations of the vibration-rotation-tunneling levels of (H(2)O)(2) and (D(2)O)(2), a very sensitive test of the potential surface, gave results in good agreement with experimental high-resolution spectra. Also the virial coefficients and properties of liquid water agree well with measured values. The present model performs better than published force fields for water in a simultaneous reproduction of experimental data for dimer spectra, virials, and properties of the liquid.
Collapse
Affiliation(s)
- Robert Bukowski
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | | | | | | |
Collapse
|
40
|
Nasrabad AE. Monte Carlo simulations of thermodynamic and structural properties of Mie(14,7) fluids. J Chem Phys 2008; 128:154514. [DOI: 10.1063/1.2901164] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
Li J, Zhou Z, Sadus RJ. Role of nonadditive forces on the structure and properties of liquid water. J Chem Phys 2007; 127:154509. [PMID: 17949175 DOI: 10.1063/1.2786449] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The role of nonadditive interactions on the structure and dielectric properties of water is investigated at different temperatures using molecular dynamics. A new intermolecular potential is developed which contains an ab initio description of two-body additive interactions plus nonadditive contributions from both three-body interactions and polarization. Polarization is the main nonadditive influence, resulting in improved agreement with experiment for the radial distribution function, dielectric constant, and dipole moment. A comparison is also made with other widely used intermolecular potentials. The new potential provides a superior prediction of the dielectric constant and dipole moment. It also predicts the relative contribution of hydrogen bonding better than the SPC/E potential [Berendsen et al., J. Phys. Chem. 91, 6269 (1987)].
Collapse
Affiliation(s)
- Jianhui Li
- Centre for Molecular Simulation, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| | | | | |
Collapse
|
42
|
Cencek W, Jeziorska M, Akin-Ojo O, Szalewicz K. Three-Body Contribution to the Helium Interaction Potential. J Phys Chem A 2007; 111:11311-9. [PMID: 17595067 DOI: 10.1021/jp072106n] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Two nonadditive three-body analytic potentials for helium were obtained: one based on three-body symmetry-adapted perturbation theory (SAPT) and the other one on supermolecular coupled-cluster theory with single, double, and noniterative triple excitations [CCSD(T)]. Large basis sets were used, up to the quintuple-zeta doubly augmented size. The fitting functions contain an exponentially decaying component describing the short-range interactions and damped inverse powers expansions for the third- and fourth-order dispersion contributions. The SAPT and CCSD(T) potentials are very close to each other. The largest uncertainty of the potentials comes from the truncation of the level of theory and can be estimated to be about 10 mK or 10% at trimer's minimum configuration. The relative uncertainties for other configurations are also expected to be about 10% except for regions where the nonadditive contribution crosses zero. Such uncertainties are of the same order of magnitude as the current uncertainties of the two-body part of the potential.
Collapse
Affiliation(s)
- Wojciech Cencek
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | | | | | | |
Collapse
|
43
|
Eskandari Nasrabad A, Laghaei R. Computational studies on thermodynamic properties, effective diameters, and free volume of argon using anab initiopotential. J Chem Phys 2006; 125:084510. [PMID: 16965032 DOI: 10.1063/1.2338310] [Citation(s) in RCA: 24] [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 quantum mechanical derived ab initio interaction potential for the argon dimer was tested in molecular simulations to reproduce the thermophysical properties of the vapor-liquid phase equilibria using the Gibbs ensemble Monte Carlo simulations as well as the liquid and supercritical equation of state using the NVT Monte Carlo simulations. The ab initio interaction potential was taken from the literature. A recently developed theory [R. Laghaei et al., J. Chem. Phys. 124, 154502 (2006)] was used to compute the effective diameters of argon in fluid phases and the results were subsequently applied in the generic van der Waals theory to compute the free volume of argon. The calculated densities of the coexisting phases, the vapor pressure, and the equation of state show excellent agreement with experimental values. The effective diameters and free volumes of argon are given over a wide range of densities and temperatures. An empirical formula was used to fit the effective diameters as a function of density and temperature. The computed free volume will be used in future investigations to calculate the transport properties of argon.
Collapse
|
44
|
Wang L, Sadus RJ. Influence of two-body and three-body interatomic forces on gas, liquid, and solid phases. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2006; 74:021202. [PMID: 17025412 DOI: 10.1103/physreve.74.021202] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/06/2006] [Indexed: 05/12/2023]
Abstract
Accurate molecular dynamics simulations are reported which quantify the contributions of two- and three-body interactions in the gas, liquid, and solid phases of argon at both subcritical and supercritical conditions. The calculations use an accurate two-body potential in addition to contributions from three-body dispersion interactions from third-order triple-dipole interactions. The number dependence of three-body interactions is quantified, indicating that a system size of at least five hundred atoms is required for reliable calculations. The results indicate that, although the contribution of three-body interaction to the overall energy is small, three-body interactions significantly affect the pressure at which vapor-liquid and solid-liquid transitions are observed. In particular, three-body interactions substantially increase the pressure of the freezing point. Unlike two-body interactions, which vary with both density and temperature, for a given density, three-body interactions have a near-constant 'background' value irrespective of the temperature. Both two-body interactions and kinetic energy have an important role in vapor-liquid equilibria whereas solid-liquid equilibria are dominated by two-body interactions.
Collapse
Affiliation(s)
- Liping Wang
- Centre for Molecular Simulation, Swinburne University of Technology, PO Box 218 Hawthorn, Victoria 3122, Australia
| | | |
Collapse
|
45
|
Sadus RJ. Molecular simulation of the thermophysical properties of fluids: phase behaviour and transport properties. MOLECULAR SIMULATION 2006. [DOI: 10.1080/08927020600592977] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
46
|
Patkowski K, Murdachaew G, Fou CM, Szalewicz * K. Accurateab initiopotential for argon dimer including highly repulsive region. Mol Phys 2005. [DOI: 10.1080/00268970500130241] [Citation(s) in RCA: 93] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
47
|
Baccarelli I, Gianturco FA, González-Lezana T, Delgado-Barrio G, Miret-Artés S, Villarreal P. Bound-state energies in argon trimers via a variational expansion: The effects from many-body corrections. J Chem Phys 2005; 122:144319. [PMID: 15847534 DOI: 10.1063/1.1879972] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
In this paper we study the bound-state energies and geometries of Ar(3) for J=0, using the distributed Gaussian functions method that provides a configurational description of the different structures contributing to these states. Atom-atom potentials are employed and three-body long-range effects are also included in the computational treatment by adding to the sum of potentials the Axilrod-Teller triple-dipole correction for the whole rotationless energy spectrum. An estimate of the total number of bound states for the Ar trimer is given. With respect to previous calculations, limited to the lower-lying states, our results show slightly larger nonadditive effects and are further able to predict the full range of the bound spectrum. Changes on the geometries of a large part of the vibrationally excited states of Ar(3) when the Axilrod-Teller term is included in the molecular potential are found by the present study.
Collapse
Affiliation(s)
- I Baccarelli
- Department of Chemistry, Instituto Nationale di fisica Nucleare, University of Rome La Sapienza, 00185 Rome, Italy
| | | | | | | | | | | |
Collapse
|
48
|
Affiliation(s)
- Dongxu Li
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637
| | - Stuart A. Rice
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637
| |
Collapse
|
49
|
Garrison SL, Sandler SI. An Accurate Acetylene Intermolecular Potential for Phase Behavior Predictions from Quantum Chemistry. J Phys Chem B 2004. [DOI: 10.1021/jp048731z] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Stephen L. Garrison
- Department of Chemical Engineering, Center for Molecular and Engineering Thermodynamics, University of Delaware, Newark, Delaware, 19716
| | - Stanley I. Sandler
- Department of Chemical Engineering, Center for Molecular and Engineering Thermodynamics, University of Delaware, Newark, Delaware, 19716
| |
Collapse
|
50
|
Nasrabad AE, Laghaei R, Deiters UK. Prediction of the thermophysical properties of pure neon, pure argon, and the binary mixtures neon-argon and argon-krypton by Monte Carlo simulation usingab initiopotentials. J Chem Phys 2004; 121:6423-34. [PMID: 15446941 DOI: 10.1063/1.1783271] [Citation(s) in RCA: 74] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Gibbs ensemble Monte Carlo simulations were used to test the ability of intermolecular pair potentials derived ab initio from quantum mechanical principles, enhanced by Axilrod-Teller triple-dipole interactions, to predict the vapor-liquid phase equilibria of pure neon, pure argon, and the binary mixtures neon-argon and argon-krypton. The interaction potentials for Ne-Ne, Ar-Ar, Kr-Kr, and Ne-Ar were taken from literature; for Ar-Kr a different potential has been developed. In all cases the quantum mechanical calculations had been carried out with the coupled-cluster approach [CCSD(T) level of theory] and with correlation consistent basis sets; furthermore an extrapolation scheme had been applied to obtain the basis set limit of the interaction energies. The ab initio pair potentials as well as the thermodynamic data based on them are found to be in excellent agreement with experimental data; the only exception is neon. It is shown, however, that in this case the deviations can be quantitatively explained by quantum effects. The interaction potentials that have been developed permit quantitative predictions of high-pressure phase equilibria of noble-gas mixtures.
Collapse
Affiliation(s)
- A E Nasrabad
- Institute of Physical Chemistry, University of Cologne, Luxemburger Strasse 116, 50939 Koeln, Germany
| | | | | |
Collapse
|