Development of a Flexible‐Monomer Two‐Body Carbon Dioxide Potential and Its Application to Clusters up to (CO2)13

A flexible‐monomer two‐body potential energy function was developed that approaches the high level CCSD(T)/CBS potential energy surface (PES) of carbon dioxide (CO2) systems. This function was generated by fitting the electronic energies of unique CO2 monomers and dimers to permutationally invariant polynomials. More than 200,000 CO2 configurations were used to train the potential function. Comparisons of the PESs of six orientations of flexible CO2 dimers were evaluated to demonstrate the accuracy of the potential. Furthermore, the potential function was used to determine the minimum energy structures of CO2 clusters containing as many as 13 molecules. For isomers of (CO2)3, the potential demonstrated energetic agreement with the M06‐2X functional and structural agreement of the B2PLYP‐D functional at substantially reduced computational costs. A separate function, fit to MP2/aug‐cc‐pVDZ reference energies, was developed to directly compare the two‐body potential to the ab initio MP2 level of theory. © 2017 Wiley Periodicals, Inc.

[1]  R. Fletcher,et al.  Practical Methods of Optimization: Fletcher/Practical Methods of Optimization , 2000 .

[2]  R. Georges,et al.  Infrared spectroscopy of (CO2)N nanoparticles (30 , 2003 .

[3]  Zhen Xie,et al.  Permutationally Invariant Polynomial Basis for Molecular Energy Surface Fitting via Monomial Symmetrization. , 2010, Journal of chemical theory and computation.

[4]  V. Babin,et al.  Development of a "First Principles" Water Potential with Flexible Monomers. II: Trimer Potential Energy Surface, Third Virial Coefficient, and Small Clusters. , 2014, Journal of chemical theory and computation.

[5]  N. Moazzen-Ahmadi,et al.  Spectroscopy of dimers, trimers and larger clusters of linear molecules , 2013 .

[6]  A. Wodtke,et al.  Electron attachment time-of-flight mass spectrometry reveals fcc bulk-like packing in CO2 aggregates , 2002 .

[7]  N. Moazzen-Ahmadi,et al.  High resolution infrared spectroscopy of carbon dioxide clusters up to (CO2)13. , 2011, The Journal of chemical physics.

[8]  M. Head‐Gordon,et al.  A fifth-order perturbation comparison of electron correlation theories , 1989 .

[9]  H. Lüthi,et al.  HIGH-LEVEL AB INITIO COMPUTATIONS OF STRUCTURES AND RELATIVE ENERGIES OF TWO ISOMERS OF THE CO2 TRIMER , 1999 .

[10]  M. Klein,et al.  Ab initio molecular dynamics study of supercritical carbon dioxide including dispersion corrections. , 2009, The Journal of chemical physics.

[11]  M. Hochlaf,et al.  Explicit correlation treatment of the potential energy surface of CO2 dimer. , 2014, The Journal of chemical physics.

[12]  M J Elrod,et al.  Many-body effects in intermolecular forces. , 1994, Chemical reviews.

[13]  Kazuo Kitaura,et al.  Exploring chemistry with the fragment molecular orbital method. , 2012, Physical chemistry chemical physics : PCCP.

[14]  S. R. Gadre,et al.  Quantum chemical investigations on molecular clusters. , 2014, Chemical reviews.

[15]  S. R. Gadre,et al.  Structures, energetics and vibrational spectra of CO2 clusters through molecular tailoring and cluster building algorithm. , 2012, Physical chemistry chemical physics : PCCP.

[16]  K. Jose,et al.  An ab initio investigation on (CO2)n and CO2(Ar)m clusters: geometries and IR spectra. , 2008, The Journal of chemical physics.

[17]  Dominic A Sirianni,et al.  Comparison of Explicitly Correlated Methods for Computing High-Accuracy Benchmark Energies for Noncovalent Interactions. , 2017, Journal of chemical theory and computation.

[18]  B. Cabral,et al.  A first principles approach to the electronic properties of liquid and supercritical CO2. , 2015, The Journal of chemical physics.

[19]  Kuang Yu,et al.  Many-body effects are essential in a physically motivated CO2 force field. , 2012, The Journal of chemical physics.

[20]  D. Truhlar,et al.  The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals , 2008 .

[21]  G. Giffin,et al.  Effect of high pressure CO2 on the structure of PMMA: a FT-IR study. , 2011, The journal of physical chemistry. B.

[22]  Pavel Hobza,et al.  Describing Noncovalent Interactions beyond the Common Approximations: How Accurate Is the "Gold Standard," CCSD(T) at the Complete Basis Set Limit? , 2013, Journal of chemical theory and computation.

[23]  Andrea Lombardi,et al.  A full dimensional grid empowered simulation of the CO2 + CO2 processes , 2012, J. Comput. Chem..

[24]  Trygve Helgaker,et al.  Basis-set convergence in correlated calculations on Ne, N2, and H2O , 1998 .

[25]  H. Huber,et al.  The CO2 Trimer: Between Gas and Fluid Phase , 2001 .

[26]  Volodymyr Babin,et al.  Development of a "First Principles" Water Potential with Flexible Monomers: Dimer Potential Energy Surface, VRT Spectrum, and Second Virial Coefficient. , 2014, Journal of chemical theory and computation.

[27]  Q. Gong,et al.  Identifying isomers of carbon-dioxide clusters by laser-driven Coulomb explosion , 2014 .

[28]  D. Truhlar,et al.  Exploring the Limit of Accuracy of the Global Hybrid Meta Density Functional for Main-Group Thermochemistry, Kinetics, and Noncovalent Interactions. , 2008, Journal of chemical theory and computation.

[29]  Water pair and three-body potential of spectroscopic quality from Ab initio calculations , 2000, Physical review letters.

[30]  Volodymyr Babin,et al.  Development of a "First-Principles" Water Potential with Flexible Monomers. III. Liquid Phase Properties. , 2014, Journal of chemical theory and computation.

[31]  Andreas W Götz,et al.  On the representation of many-body interactions in water. , 2015, The Journal of chemical physics.

[32]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[33]  Teresa Head-Gordon,et al.  Convergence of the Many-Body Expansion for Energy and Forces for Classical Polarizable Models in the Condensed Phase. , 2016, Journal of chemical theory and computation.

[34]  N. Moazzen-Ahmadi,et al.  Spectroscopic identification of carbon dioxide clusters: (CO2)6 to (CO2)13. , 2011, Physical chemistry chemical physics : PCCP.

[35]  Harry Partridge,et al.  The determination of an accurate isotope dependent potential energy surface for water from extensive ab initio calculations and experimental data , 1997 .

[36]  T. Seward,et al.  Thermodynamic properties of carbon dioxide clusters by M06-2X and dispersion-corrected B2PLYP-D theory , 2013 .

[37]  H. Liu,et al.  Finite temperature properties of (CO 2 ) n clusters , 2003 .

[38]  Ryan P A Bettens,et al.  Energy-Based Molecular Fragmentation Methods. , 2015, Chemical reviews.

[39]  Y. Negishi,et al.  Structural evolution in (CO2)n clusters (n<103) as studied by mass spectrometry , 2002 .

[40]  Betsy M. Rice,et al.  Intermolecular potential of carbon dioxide dimer from symmetry-adapted perturbation theory , 1999 .

[41]  George Jackson,et al.  SAFT-γ force field for the simulation of molecular fluids. 1. A single-site coarse grained model of carbon dioxide. , 2011, The journal of physical chemistry. B.

[42]  N. Moazzen-Ahmadi,et al.  Five intermolecular vibrations of the CO2 dimer observed via infrared combination bands. , 2016, The Journal of chemical physics.

[43]  A. McKellar,et al.  High-resolution infrared spectroscopy of carbon dioxide dimers, trimers, and larger clusters , 2010 .

[44]  Soumya Ganguly Neogi,et al.  Structural and spectroscopic studies of carbon dioxide clusters: a combined genetic algorithm and DFT based study , 2014, Structural Chemistry.

[45]  Kamy Sepehrnoori,et al.  CO2 injection for enhanced oil recovery in Bakken tight oil reservoirs , 2015 .

[46]  Bin Jiang,et al.  Potential energy surfaces from high fidelity fitting of ab initio points: the permutation invariant polynomial - neural network approach , 2016 .

[47]  Hiroshi Takeuchi,et al.  Geometry optimization of carbon dioxide clusters (CO2)n for 4 < or = n < or = 40. , 2008, The journal of physical chemistry. A.

[48]  Krzysztof Szalewicz,et al.  Predictions of the Properties of Water from First Principles , 2007, Science.

[49]  Rita Prosmiti,et al.  i-TTM Model for Ab Initio-Based Ion-Water Interaction Potentials. 1. Halide-Water Potential Energy Functions. , 2016, The journal of physical chemistry. B.

[50]  S. Grimme,et al.  Double-hybrid density functionals with long-range dispersion corrections: higher accuracy and extended applicability. , 2007, Physical chemistry chemical physics : PCCP.

[51]  A. Tikhonov,et al.  Numerical Methods for the Solution of Ill-Posed Problems , 1995 .

[52]  G. Beran,et al.  Predicting finite-temperature properties of crystalline carbon dioxide from first principles with quantitative accuracy , 2015, Chemical science.

[53]  Hans Peter Lüthi,et al.  Interaction energies of van der Waals and hydrogen bonded systems calculated using density functional theory: Assessing the PW91 model , 2001 .

[54]  R. Georges,et al.  IR signature of (CO2)N clusters: size, shape and structural effects. , 2005, Physical chemistry chemical physics : PCCP.

[55]  K. J. Børve,et al.  Size of free neutral CO2 clusters from carbon 1s ionization energies. , 2011, The journal of physical chemistry. A.

[56]  Gábor Csányi,et al.  Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions , 2016, Chemical reviews.

[57]  Sotiris S. Xantheas,et al.  Development of transferable interaction models for water. IV. A flexible, all-atom polarizable potential (TTM2-F) based on geometry dependent charges derived from an ab initio monomer dipole moment surface , 2002 .

[58]  D. Nesbitt,et al.  Geometric isomerism in clusters: High resolution infrared spectroscopy of a noncyclic CO2 trimer , 1996 .

[59]  J. Light,et al.  Vibrations of the carbon dioxide dimer , 2000 .

[60]  Shridhar R Gadre,et al.  Molecular tailoring approach: a route for ab initio treatment of large clusters. , 2014, Accounts of chemical research.

[61]  Spencer R Pruitt,et al.  Fragmentation methods: a route to accurate calculations on large systems. , 2012, Chemical reviews.

[62]  Joel M. Bowman,et al.  Permutationally invariant potential energy surfaces in high dimensionality , 2009 .

[63]  Marat Valiev,et al.  Fast electron correlation methods for molecular clusters in the ground and excited states , 2005 .

[64]  Ian R. McDonald,et al.  Interaction site models for carbon dioxide , 1981 .

[65]  Jeanette M. Sperhac,et al.  High‐resolution infrared diode laser spectroscopy of (CO2)3: Vibrationally averaged structures, resonant dipole vibrational shifts, and tests of CO2–CO2 pair potentials , 1995 .

[66]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[67]  Trygve Helgaker,et al.  Basis set convergence of the interaction energy of hydrogen-bonded complexes , 1999 .

[68]  G. T. Fraser,et al.  Structure and vibrational dynamics of the CO2 dimer from the sub‐Doppler infrared spectrum of the 2.7 μm Fermi diad , 1988 .

[69]  Stefan Bachu,et al.  CO2 storage in geological media: Role, means, status and barriers to deployment , 2008 .