1
|
Masumian E, Boese AD. Benchmarking Swaths of Intermolecular Interaction Components with Symmetry-Adapted Perturbation Theory. J Chem Theory Comput 2024; 20:30-48. [PMID: 38117939 PMCID: PMC10782453 DOI: 10.1021/acs.jctc.3c00801] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2023] [Revised: 11/19/2023] [Accepted: 11/22/2023] [Indexed: 12/22/2023]
Abstract
A benchmark database for interaction energy components of various noncovalent interactions (NCIs) along their dissociation curve is one of the essential needs in theoretical chemistry, especially for the development of force fields and machine-learning methods. We utilize DFT-SAPT or SAPT(DFT) as one of the most accurate methods to generate an extensive stock of the energy components, including dispersion energies extrapolated to the complete basis set limit (CBS). Precise analyses of the created data, and benchmarking the total interaction energies against the best available CCSD(T)/CBS values, reveal different aspects of the methodology and the nature of NCIs. For example, error cancellation effects between the S2 approximation and nonexact xc-potentials occur, and large charge transfer energies in some systems, including heavy atoms, can explain the lower accuracy of DFT-SAPT. This method is perfect for neutral complexes containing light nonmetals, while other systems with heavier atoms should be treated carefully. In the last part, a representative data set for all NCIs is extracted from the original data.
Collapse
Affiliation(s)
- Ehsan Masumian
- Physical and Theoretical Chemistry,
Department of Chemistry, University of Graz, 8010 Graz, Austria
| | - A. Daniel Boese
- Physical and Theoretical Chemistry,
Department of Chemistry, University of Graz, 8010 Graz, Austria
| |
Collapse
|
2
|
Garcia J, Podeszwa R, Szalewicz K. SAPT codes for calculations of intermolecular interaction energies. J Chem Phys 2020; 152:184109. [PMID: 32414261 DOI: 10.1063/5.0005093] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
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
|
3
|
Patkowski K. Recent developments in symmetry‐adapted perturbation theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1452] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Konrad Patkowski
- Department of Chemistry and Biochemistry Auburn University Auburn Alabama
| |
Collapse
|
4
|
Yokogawa D. Isotropic Site-Site Dispersion Potential Determined from Localized Frequency-Dependent Density Susceptibility. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2019. [DOI: 10.1246/bcsj.20190142] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Daisuke Yokogawa
- Graduate School of Arts and Science, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| |
Collapse
|
5
|
Misquitta AJ, Stone AJ. ISA-Pol: distributed polarizabilities and dispersion models from a basis-space implementation of the iterated stockholder atoms procedure. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2371-4] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
6
|
Analysis of polarization in hydrogen bonded complexes: An asymptotic projection approach. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.02.034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
7
|
Van Vleet MJ, Misquitta AJ, Schmidt JR. New Angles on Standard Force Fields: Toward a General Approach for Treating Atomic-Level Anisotropy. J Chem Theory Comput 2018; 14:739-758. [DOI: 10.1021/acs.jctc.7b00851] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Mary J. Van Vleet
- Theoretical
Chemistry Institute and Department of Chemistry, University of Wisconsin−Madison, Madison, Wisconsin 53706, United States
| | - Alston J. Misquitta
- Department
of Physics and Astronomy, Queen Mary University of London, London E1 4NS, United Kingdom
| | - J. R. Schmidt
- Theoretical
Chemistry Institute and Department of Chemistry, University of Wisconsin−Madison, Madison, Wisconsin 53706, United States
| |
Collapse
|
8
|
Aina AA, Misquitta AJ, Price SL. From dimers to the solid-state: Distributed intermolecular force-fields for pyridine. J Chem Phys 2017; 147:161722. [DOI: 10.1063/1.4999789] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alexander A. Aina
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Alston J. Misquitta
- School of Physics and Astronomy, Queen Mary, University of London, London E1 4NS, United Kingdom
| | - Sarah L. Price
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| |
Collapse
|
9
|
LYNDEN-BELL RUTHM, STONE ANTHONYJ. A model with charges and polarizability for CS2 in an ionic liquid. J CHEM SCI 2017. [DOI: 10.1007/s12039-017-1243-0] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
10
|
Nocito D, Beran GJO. Averaged Condensed Phase Model for Simulating Molecules in Complex Environments. J Chem Theory Comput 2017; 13:1117-1129. [PMID: 28170251 DOI: 10.1021/acs.jctc.6b00890] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The need for configurational sampling dramatically increases the cost of combined quantum mechanics/molecular mechanics (QM/MM) simulations of chemical processes in solution. We developed an averaged condensed phase environment (ACPE) model that constructs an effective polarizable environment directly from explicitly sampled molecular dynamics configurations via the K-means++ algorithm and a mathematically rigorous translation of the molecular mechanics parameters. The model captures detailed heterogeneous features in the environment that may be difficult to describe using a conventional polarizable continuum model. Instead of performing repeated QM/MM calculations for each new configuration of the environment, the ACPE approach allows one to perform a single QM calculation on an averaged configuration. Here, we demonstrate the model by computing electronic excitation energies for several small molecules in solution. The ACPE model predicts the excitation energies in excellent agreement with conventional configurational averaging yet with orders of magnitude of reduction in the computational cost.
Collapse
Affiliation(s)
- Dominique Nocito
- Department of Chemistry, University of California , Riverside, California 92521, United States
| | - Gregory J O Beran
- Department of Chemistry, University of California , Riverside, California 92521, United States
| |
Collapse
|
11
|
Wang H, Yang W. Determining polarizable force fields with electrostatic potentials from quantum mechanical linear response theory. J Chem Phys 2017; 144:224107. [PMID: 27305996 DOI: 10.1063/1.4953558] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We developed a new method to calculate the atomic polarizabilities by fitting to the electrostatic potentials (ESPs) obtained from quantum mechanical (QM) calculations within the linear response theory. This parallels the conventional approach of fitting atomic charges based on electrostatic potentials from the electron density. Our ESP fitting is combined with the induced dipole model under the perturbation of uniform external electric fields of all orientations. QM calculations for the linear response to the external electric fields are used as input, fully consistent with the induced dipole model, which itself is a linear response model. The orientation of the uniform external electric fields is integrated in all directions. The integration of orientation and QM linear response calculations together makes the fitting results independent of the orientations and magnitudes of the uniform external electric fields applied. Another advantage of our method is that QM calculation is only needed once, in contrast to the conventional approach, where many QM calculations are needed for many different applied electric fields. The molecular polarizabilities obtained from our method show comparable accuracy with those from fitting directly to the experimental or theoretical molecular polarizabilities. Since ESP is directly fitted, atomic polarizabilities obtained from our method are expected to reproduce the electrostatic interactions better. Our method was used to calculate both transferable atomic polarizabilities for polarizable molecular mechanics' force fields and nontransferable molecule-specific atomic polarizabilities.
Collapse
Affiliation(s)
- Hao Wang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
12
|
Huang J, Mei Y, König G, Simmonett AC, Pickard FC, Wu Q, Wang LP, MacKerell AD, Brooks BR, Shao Y. An Estimation of Hybrid Quantum Mechanical Molecular Mechanical Polarization Energies for Small Molecules Using Polarizable Force-Field Approaches. J Chem Theory Comput 2017; 13:679-695. [PMID: 28081366 DOI: 10.1021/acs.jctc.6b01125] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In this work, we report two polarizable molecular mechanics (polMM) force field models for estimating the polarization energy in hybrid quantum mechanical molecular mechanical (QM/MM) calculations. These two models, named the potential of atomic charges (PAC) and potential of atomic dipoles (PAD), are formulated from the ab initio quantum mechanical (QM) response kernels for the prediction of the QM density response to an external molecular mechanical (MM) environment (as described by external point charges). The PAC model is similar to fluctuating charge (FQ) models because the energy depends on external electrostatic potential values at QM atomic sites; the PAD energy depends on external electrostatic field values at QM atomic sites, resembling induced dipole (ID) models. To demonstrate their uses, we apply the PAC and PAD models to 12 small molecules, which are solvated by TIP3P water. The PAC model reproduces the QM/MM polarization energy with a R2 value of 0.71 for aniline (in 10,000 TIP3P water configurations) and 0.87 or higher for other 11 solute molecules, while the PAD model has a much better performance with R2 values of 0.98 or higher. The PAC model reproduces reference QM/MM hydration free energies for 12 solute molecules with a RMSD of 0.59 kcal/mol. The PAD model is even more accurate, with a much smaller RMSD of 0.12 kcal/mol, with respect to the reference. This suggests that polarization effects, including both local charge distortion and intramolecular charge transfer, can be well captured by induced dipole type models with proper parametrization.
Collapse
Affiliation(s)
- Jing Huang
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , 20 Penn Street, Baltimore, Maryland 21201, United States.,Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute , 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University , Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Gerhard König
- Max-Planck-Institut für Kohlenforschung , 45470 Mülheim an der Ruhr, NRW Germany, EU
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute , 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Frank C Pickard
- Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute , 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Qin Wu
- Center for Functional Nanomaterials, Brookhaven National Laboratory , Upton, New York 11973, United States
| | - Lee-Ping Wang
- Department of Chemistry, University of California , 1 Shields Avenue, Davis, California 95616, United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , 20 Penn Street, Baltimore, Maryland 21201, United States
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute , 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Yihan Shao
- Q-Chem Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, United States.,Department of Chemistry and Biochemistry, University of Oklahoma , Norman, Oklahoma 73019, United States
| |
Collapse
|
13
|
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
|
14
|
Van Vleet MJ, Misquitta AJ, Stone AJ, Schmidt JR. Beyond Born-Mayer: Improved Models for Short-Range Repulsion in ab Initio Force Fields. J Chem Theory Comput 2016; 12:3851-70. [PMID: 27337546 DOI: 10.1021/acs.jctc.6b00209] [Citation(s) in RCA: 57] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Short-range repulsion within intermolecular force fields is conventionally described by either Lennard-Jones (A/r(12)) or Born-Mayer (A exp(-Br)) forms. Despite their widespread use, these simple functional forms are often unable to describe the interaction energy accurately over a broad range of intermolecular distances, thus creating challenges in the development of ab initio force fields and potentially leading to decreased accuracy and transferability. Herein, we derive a novel short-range functional form based on a simple Slater-like model of overlapping atomic densities and an iterated stockholder atom (ISA) partitioning of the molecular electron density. We demonstrate that this Slater-ISA methodology yields a more accurate, transferable, and robust description of the short-range interactions at minimal additional computational cost compared to standard Lennard-Jones or Born-Mayer approaches. Finally, we show how this methodology can be adapted to yield the standard Born-Mayer functional form while still retaining many of the advantages of the Slater-ISA approach.
Collapse
Affiliation(s)
- Mary J Van Vleet
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison , Madison, Wisconsin 53706, United States
| | - Alston J Misquitta
- Department of Physics and Astronomy, Queen Mary University of London , London E1 4NS, United Kingdom
| | - Anthony J Stone
- Department of Chemistry, University of Cambridge , Cambridge CB2 1EW, United Kingdom
| | - J R Schmidt
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison , Madison, Wisconsin 53706, United States
| |
Collapse
|
15
|
McDaniel JG, Choi E, Son CY, Schmidt JR, Yethiraj A. Ab Initio Force Fields for Imidazolium-Based Ionic Liquids. J Phys Chem B 2016; 120:7024-36. [DOI: 10.1021/acs.jpcb.6b05328] [Citation(s) in RCA: 53] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jesse G. McDaniel
- Department of Chemistry and ‡Department of Physics, University of Wisconsin—Madison, Madison, Wisconsin 53706, United States
| | - Eunsong Choi
- Department of Chemistry and ‡Department of Physics, University of Wisconsin—Madison, Madison, Wisconsin 53706, United States
| | - Chang Yun Son
- Department of Chemistry and ‡Department of Physics, University of Wisconsin—Madison, Madison, Wisconsin 53706, United States
| | - J. R. Schmidt
- Department of Chemistry and ‡Department of Physics, University of Wisconsin—Madison, Madison, Wisconsin 53706, United States
| | - Arun Yethiraj
- Department of Chemistry and ‡Department of Physics, University of Wisconsin—Madison, Madison, Wisconsin 53706, United States
| |
Collapse
|
16
|
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
|
17
|
McDaniel JG, Choi E, Son CY, Schmidt JR, Yethiraj A. Conformational and Dynamic Properties of Poly(ethylene oxide) in an Ionic Liquid: Development and Implementation of a First-Principles Force Field. J Phys Chem B 2016; 120:231-43. [DOI: 10.1021/acs.jpcb.5b10065] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jesse G. McDaniel
- Department
of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, United States
| | - Eunsong Choi
- Department
of Physics, University of Wisconsin, Madison, Wisconsin 53706, United States
| | - Chang-Yun Son
- Department
of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, United States
| | - J. R. Schmidt
- Department
of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, United States
| | - Arun Yethiraj
- Department
of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, United States
| |
Collapse
|
18
|
Misquitta AJ, Stone AJ. Accurate Induction Energies for Small Organic Molecules: 1. Theory. J Chem Theory Comput 2015; 4:7-18. [PMID: 26619975 DOI: 10.1021/ct700104t] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
The induction energy often plays a very important role in determining the structure and properties of clusters of organic molecules, but only in recent years has an effort been made to include this energy in such calculations, notably in the field of organic crystal structure prediction. In this paper and the following one in this issue we provide ab initio methods suitable for the accurate inclusion of the induction energy for molecules containing as many as 30 atoms or so. These techniques are based on Symmetry-Adapted Perturbation Theory using Density Functional Theory [SAPT(DFT)] and use distributed polarizabilities computed using the recently developed density-fitting algorithm with constrained refinement. With this approach we are able to obtain induction models of varying complexity and study the effects of overlap and related numerical issues. Basis set effects on the exact and asymptotic induction energies are investigated, and the roles of higher-order induction energies and many-body effects are explored.
Collapse
Affiliation(s)
- Alston J Misquitta
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, U.K., and University College London, 20 Gordon Street, London WC1H 0AJ, U.K
| | - Anthony J Stone
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, U.K., and University College London, 20 Gordon Street, London WC1H 0AJ, U.K
| |
Collapse
|
19
|
Millot C, Chaumont A, Engler E, Wipff G. Distributed polarizability models for imidazolium-based ionic liquids. J Phys Chem A 2014; 118:8842-51. [PMID: 25133873 DOI: 10.1021/jp505539y] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Quantum chemical calculations are used to derive distributed polarizability models sufficiently accurate and compact to be used in classical molecular dynamics simulations of imidazolium-based room temperature ionic liquids. Two distributed polarizability models are fitted to reproduce the induction energy of three imidazolium cations (1,3-dimethyl-, 1-ethyl-3-methyl-, and 1-butyl-3-methylimidazolium) and four anions (tetrafluoroborate, hexafluorophosphate, nitrate, and thiocyanate) polarized by a point charge located successively on a grid of surrounding points. The first model includes charge-flow polarizabilities between first-neighbor atoms and isotropic dipolar polarizability on all atoms (except H), while the second model includes anisotropic dipolar polarizabilities on all atoms (except H). For the imidazolium cations, particular attention is given to the transferability of the distributed polarizability sets. The molecular polarizability and its anisotropy rebuilt by the distributed models are found to be in good agreement with the exact ab initio values for the three cations and 23 additional conformers of 1-ethyl-3-methyl-, 1-butyl-3-methyl-, 1-pentyl-3-methyl-, and 1-hexyl-3-methylimidazolium cations.
Collapse
Affiliation(s)
- Claude Millot
- Université de Lorraine, CNRS, SRSMC , UMR 7565, Equipes TMS/ReSolve, Faculté des Sciences et Technologies, Boulevard des Aiguillettes, BP 70239, Vandoeuvre-lès-Nancy F-54506, France
| | | | | | | |
Collapse
|
20
|
Song HJ, Zhang YG, Li H, Zhou T, Huang FL. All-atom, non-empirical, and tailor-made force field for α-RDX from first principles. RSC Adv 2014. [DOI: 10.1039/c4ra07195f] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
|
21
|
Choi E, McDaniel JG, Schmidt JR, Yethiraj A. First-Principles, Physically Motivated Force Field for the Ionic Liquid [BMIM][BF4]. J Phys Chem Lett 2014; 5:2670-2674. [PMID: 26277961 DOI: 10.1021/jz5010945] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Molecular simulations play an important role in establishing structure-property relations in complex fluids such as room-temperature ionic liquids. Classical force fields are the starting point when large systems or long times are of interest. These force fields must be not only accurate but also transferable. In this work, we report a physically motivated force field for the ionic liquid 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM][BF4]) based on symmetry-adapted perturbation theory. The predictions (from molecular dynamics simulations) of the liquid density, enthalpy of vaporization, diffusion coefficients, viscosity, and conductivity are in excellent agreement with experiment, with no adjustable parameters. The explicit energy decomposition inherent in the force field enables a quantitative analysis of the important physical interactions in these systems. We find that polarization is crucial and there is little evidence of charge transfer. We also argue that the often used procedure of scaling down charges in molecular simulations of ionic liquids is unphysical for [BMIM][BF4]. Because all intermolecular interactions in the force field are parametrized from first-principles, we anticipate good transferability to other ionic liquid systems and physical conditions.
Collapse
Affiliation(s)
- Eunsong Choi
- †Department of Physics, University of Wisconsin, 1150 University Avenue, Madison, Wisconsin 53706, United States
| | - Jesse G McDaniel
- ‡Department of Chemistry, University of Wisconsin, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - J R Schmidt
- ‡Department of Chemistry, University of Wisconsin, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Arun Yethiraj
- ‡Department of Chemistry, University of Wisconsin, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
22
|
Krishtal A, Van Alsenoy C, Geerlings P. Evaluating interaction energies of weakly bonded systems using the Buckingham-Hirshfeld method. J Chem Phys 2014; 140:184105. [DOI: 10.1063/1.4873133] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
23
|
Xu P, Zahariev F, Gordon MS. The R–7 Dispersion Interaction in the General Effective Fragment Potential Method. J Chem Theory Comput 2014; 10:1576-87. [DOI: 10.1021/ct500017n] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Peng Xu
- Department
of Chemistry, Iowa State University, Ames, Iowa 50011, United States
| | - Federico Zahariev
- Department
of Chemistry, Iowa State University, Ames, Iowa 50011, United States
| | - Mark S. Gordon
- Department
of Chemistry, Iowa State University, Ames, Iowa 50011, United States
| |
Collapse
|
24
|
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] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
Haslmayr J, Renger T. Qualitative change of character of dispersive interaction with intermolecular distance. J Chem Phys 2013; 139:044103. [DOI: 10.1063/1.4813501] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
Rob F, Szalewicz K. Distributed molecular polarisabilities and asymptotic intermolecular interaction energies†. Mol Phys 2013. [DOI: 10.1080/00268976.2013.808770] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
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
|
27
|
|
28
|
El Kerdawy A, Murray JS, Politzer P, Bleiziffer P, Heßelmann A, Görling A, Clark T. Directional Noncovalent Interactions: Repulsion and Dispersion. J Chem Theory Comput 2013; 9:2264-75. [DOI: 10.1021/ct400185f] [Citation(s) in RCA: 61] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ahmed El Kerdawy
- Computer-Chemie-Centrum, Department
Chemie und Pharmazie, Friedrich-Alexander-Universität Erlangen-Nürnberg,
Nägelsbachstrasse 25, 91052 Erlangen, Germany
| | - Jane S. Murray
- CleveTheoComp, 1951 W. 26th
Street, Suite 409, Cleveland, Ohio 44113, United States
- Department of Chemistry, University
of New Orleans, New Orleans, Louisiana 70148, United States
| | - Peter Politzer
- CleveTheoComp, 1951 W. 26th
Street, Suite 409, Cleveland, Ohio 44113, United States
- Department of Chemistry, University
of New Orleans, New Orleans, Louisiana 70148, United States
| | - Patrick Bleiziffer
- Chair of Theoretical Chemistry,
Department Chemie und Pharmazie, Friedrich-Alexander-Universität
Erlangen-Nürnberg, Egerlandstr. 3, 91058 Erlangen, Germany
| | - Andreas Heßelmann
- Chair of Theoretical Chemistry,
Department Chemie und Pharmazie, Friedrich-Alexander-Universität
Erlangen-Nürnberg, Egerlandstr. 3, 91058 Erlangen, Germany
| | - Andreas Görling
- Chair of Theoretical Chemistry,
Department Chemie und Pharmazie, Friedrich-Alexander-Universität
Erlangen-Nürnberg, Egerlandstr. 3, 91058 Erlangen, Germany
- Interdisciplinary Center for
Molecular Materials, Department Chemie und Pharmazie, Friedrich-Alexander-Universität
Erlangen-Nürnberg, Nägelsbachstrasse 25, 91052 Erlangen,
Germany
| | - Timothy Clark
- Computer-Chemie-Centrum, Department
Chemie und Pharmazie, Friedrich-Alexander-Universität Erlangen-Nürnberg,
Nägelsbachstrasse 25, 91052 Erlangen, Germany
- Interdisciplinary Center for
Molecular Materials, Department Chemie und Pharmazie, Friedrich-Alexander-Universität
Erlangen-Nürnberg, Nägelsbachstrasse 25, 91052 Erlangen,
Germany
- Centre for Molecular Design,
University
of Portsmouth, King Henry Building, King Henry I Street, Portsmouth,
PO1 2DY, United Kingdom
| |
Collapse
|
29
|
McDaniel JG, Schmidt JR. Physically-motivated force fields from symmetry-adapted perturbation theory. J Phys Chem A 2013; 117:2053-66. [PMID: 23343200 DOI: 10.1021/jp3108182] [Citation(s) in RCA: 95] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a general methodology for generating accurate and transferable ab initio force fields, employing the framework of symmetry adapted perturbation theory (SAPT). The resulting force fields are "physically-motivated" in that they contain separate, explicit terms to account for the various fundamental intermolecular interactions, such as exchange, electrostatics, induction, and dispersion, with each term parametrized to a corresponding term in the SAPT energy decomposition. Crucially, the resulting force fields are largely compatible with existing, standard simulation packages, requiring only minimal modifications. We present several novel parametrization techniques that yield robust, physically meaningful atomic parameters that are transferable between molecular environments. We demonstrate the accuracy and generality of our method by validating against experimental second virial coefficients for a variety of small molecules. We then show that the resulting atomic parameters can be combined using physically motivated ansatzes to accurately predict arbitrary heteromolecular interaction energies, with example applications including prediction of gas adsorption in functionalized metal-organic framework materials.
Collapse
Affiliation(s)
- Jesse G McDaniel
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | | |
Collapse
|
30
|
Marenich AV, Cramer CJ, Truhlar DG. Reduced and quenched polarizabilities of interior atoms in molecules. Chem Sci 2013. [DOI: 10.1039/c3sc50242b] [Citation(s) in RCA: 51] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
|
31
|
Polarization effects in protein–ligand calculations extend farther than the actual induction energy. Theor Chem Acc 2012. [DOI: 10.1007/s00214-012-1159-1] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
32
|
Totton TS, Misquitta AJ, Kraft M. Assessing the Polycyclic Aromatic Hydrocarbon Anisotropic Potential with Application to the Exfoliation Energy of Graphite. J Phys Chem A 2011; 115:13684-93. [DOI: 10.1021/jp208088s] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Tim S. Totton
- Department of Chemical Engineering and Biotechnology, University of Cambridge, New Museums Site, Pembroke Street, Cambridge CB2 3RA, United Kingdom
| | - Alston J. Misquitta
- Department of Physics, Cavendish Laboratory, University of Cambridge, J J Thomson Avenue, Cambridge, CB3 0HE, United Kingdom
| | - Markus Kraft
- Department of Chemical Engineering and Biotechnology, University of Cambridge, New Museums Site, Pembroke Street, Cambridge CB2 3RA, United Kingdom
| |
Collapse
|
33
|
Elking DM, Perera L, Duke R, Darden T, Pedersen LG. A finite field method for calculating molecular polarizability tensors for arbitrary multipole rank. J Comput Chem 2011; 32:3283-95. [PMID: 21915883 DOI: 10.1002/jcc.21914] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2011] [Revised: 07/06/2011] [Accepted: 07/25/2011] [Indexed: 11/07/2022]
Abstract
A finite field method for calculating spherical tensor molecular polarizability tensors α(lm;l'm') = ∂Δ(lm)/∂ϕ(l'm')* by numerical derivatives of induced molecular multipole Δ(lm) with respect to gradients of electrostatic potential ϕ(l'm')* is described for arbitrary multipole ranks l and l'. Interconversion formulae for transforming multipole moments and polarizability tensors between spherical and traceless Cartesian tensor conventions are derived. As an example, molecular polarizability tensors up to the hexadecapole-hexadecapole level are calculated for water using the following ab initio methods: Hartree-Fock (HF), Becke three-parameter Lee-Yang-Parr exchange-correlation functional (B3LYP), Møller-Plesset perturbation theory up to second order (MP2), and Coupled Cluster theory with single and double excitations (CCSD). In addition, intermolecular electrostatic and polarization energies calculated by molecular multipoles and polarizability tensors are compared with ab initio reference values calculated by the Reduced Variation Space method for several randomly oriented small molecule dimers separated by a large distance. It is discussed how higher order molecular polarizability tensors can be used as a tool for testing and developing new polarization models for future force fields.
Collapse
Affiliation(s)
- Dennis M Elking
- University of North Carolina, Department of Chemistry, Chapel Hill, North Carolina 27599, USA
| | | | | | | | | |
Collapse
|
34
|
Liu RF, Ángyán JG, Dobson JF. Dispersion interaction in hydrogen-chain models. J Chem Phys 2011; 134:114106. [PMID: 21428606 DOI: 10.1063/1.3563596] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We have investigated the dispersion interaction in hydrogen chain models via density functional theory-based symmetry-adapted perturbation theory using the asymptotically corrected PBE0 energy functional. The quasimetallic and the insulating prototype systems were chosen to be hydrogen chains with equally and alternately spaced H(2) units, respectively. The dependence of the dispersion energy on the chain length for quasimetallic and insulating cases has been determined for two chains arranged either in pointing or in parallel geometries. The results are compared with those previously calculated from a continuum coupled-plasmon approach [Phys. Rev. B 77, 075436 (2008)]. The interaction energy has also been modeled by pairwise summations over short fragments of the chains, demonstrating the failure of the additivity principle for the quasimetallic case, while confirming that the additivity is a qualitatively reasonable hypothesis for the insulating case.
Collapse
Affiliation(s)
- Ru-Fen Liu
- CRM2, Institut Jean Barriol, Nancy University and CNRS, 54506 Vandoeuvre-lès-Nancy, France.
| | | | | |
Collapse
|
35
|
Söderhjelm P, Kongsted J, Ryde U. Conformational Dependence of Isotropic Polarizabilities. J Chem Theory Comput 2011; 7:1404-14. [DOI: 10.1021/ct100714e] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Pär Söderhjelm
- Department of Chemistry and Applied Biosciences, Computational Science, ETH Zürich, Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
| | - Jacob Kongsted
- Department of Physics and Chemistry, University of Southern Denmark, Campusvej 55, 5230 Odense M, Denmark
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
36
|
Luque FJ, Dehez F, Chipot C, Orozco M. Polarization effects in molecular interactions. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.32] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- F. Javier Luque
- Departament de Fisicoquímica, Facultat de Farmàcia, and Institut de Biomedicina (IBUB), Universitat de Barcelona, Barcelona, Spain
| | - François Dehez
- Equipe de dynamique des assemblages membranaires, Nancy University, Vandoeuvre‐lès‐Nancy Cedex, France
| | - Christophe Chipot
- Equipe de dynamique des assemblages membranaires, Nancy University, Vandoeuvre‐lès‐Nancy Cedex, France
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana‐Champaign Urbana, Illinois, USA
| | - Modesto Orozco
- Joint IRB‐BSC Program in Computational Biology, Institut de Recerca Biomèdica, Barcelona, Spain
- Department of Life Sciences, Barcelona Supercomputing Centre, Barcelona, Spain
- Departament de Bioquímica, Facultat de Biologia, Universitat de Barcelona, Barcelona, Spain
| |
Collapse
|
37
|
Rajchel Ł, Zuchowski PS, Szcześniak MM, Chałasiński G. Density functional theory approach to noncovalent interactions via monomer polarization and Pauli blockade. PHYSICAL REVIEW LETTERS 2010; 104:163001. [PMID: 20482044 DOI: 10.1103/physrevlett.104.163001] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2009] [Indexed: 05/29/2023]
Abstract
We propose a "DFT+dispersion" treatment which avoids double counting of dispersion terms by deriving the dispersion-free density functional theory (DFT) interaction energy and combining it with DFT-based dispersion. The formalism involves self-consistent polarization of DFT monomers restrained by the exclusion principle via the Pauli-blockade technique. Any exchange-correlation potential can be used within monomers, but only the exchange operates between them. The applications to rare-gas dimers, ion-rare-gas interactions, and hydrogen bonds demonstrate that the interaction energies agree with benchmark values.
Collapse
Affiliation(s)
- Łukasz Rajchel
- Department of Chemistry, Oakland University, Rochester, Michigan 48309-4477, USA.
| | | | | | | |
Collapse
|
38
|
Sato T, Nakai H. Density functional method including weak interactions: Dispersion coefficients based on the local response approximation. J Chem Phys 2010; 131:224104. [PMID: 20001021 DOI: 10.1063/1.3269802] [Citation(s) in RCA: 154] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
A new method to calculate the atom-atom dispersion coefficients in a molecule is proposed for the use in density functional theory with dispersion (DFT-D) correction. The method is based on the local response approximation due to Dobson and Dinte [Phys. Rev. Lett. 76, 1780 (1996)], with modified dielectric model recently proposed by Vydrov and van Voorhis [J. Chem. Phys. 130, 104105 (2009)]. The local response model is used to calculate the distributed multipole polarizabilities of atoms in a molecule, from which the dispersion coefficients are obtained by an explicit frequency integral of the Casimir-Polder type. Thus obtained atomic polarizabilities are also used in the damping function for the short-range singularity. Unlike empirical DFT-D methods, the local response dispersion (LRD) method is able to calculate the dispersion energy from the ground-state electron density only. It is applicable to any geometry, free from physical constants such as van der Waals radii or atomic polarizabilities, and computationally very efficient. The LRD method combined with the long-range corrected DFT functional (LC-BOP) is applied to calculations of S22 weakly bound complex set [Phys. Chem. Chem. Phys. 8, 1985 (2006)]. Binding energies obtained by the LC-BOP+LRD agree remarkably well with ab initio references.
Collapse
Affiliation(s)
- Takeshi Sato
- Research Institute for Science and Engineering, Waseda University, Tokyo 169-8555, Japan
| | | |
Collapse
|
39
|
Totton TS, Misquitta AJ, Kraft M. A First Principles Development of a General Anisotropic Potential for Polycyclic Aromatic Hydrocarbons. J Chem Theory Comput 2010; 6:683-95. [PMID: 26613299 DOI: 10.1021/ct9004883] [Citation(s) in RCA: 57] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Standard empirical atom-atom potentials are shown to be unable to describe the binding of polycyclic aromatic hydrocarbon (PAH) molecules in the variety of configurations seen in clusters. The main reason for this inadequacy is the lack of anisotropy in these potentials. We have constructed an anisotropic atom-atom intermolecular potential for the benzene molecule from first principles using a symmetry-adapted perturbation theory based on density functional theory (SAPT(DFT)), interaction energy calculations and the Williams-Stone-Misquitta method for obtaining molecular properties in distributed form. Using this potential as a starting point, we have constructed a transferable anisotropic potential to model intermolecular interactions between PAHs. This new potential has been shown to accurately model interaction energies for a variety of dimer configurations for four different PAH molecules, including certain configurations which are poorly modeled with current isotropic potentials. It is intended that this potential will form the basis for further work on the aggregation of PAHs.
Collapse
Affiliation(s)
- Tim S Totton
- Department of Chemical Engineering and Biotechnology, University of Cambridge, New Museums Site, Pembroke Street, Cambridge CB2 3RA, United Kingdom, and Department of Physics,Cavendish Laboratory, University of Cambridge, J J Thomson Avenue, Cambridge, CB3 0HE, United Kingdom
| | - Alston J Misquitta
- Department of Chemical Engineering and Biotechnology, University of Cambridge, New Museums Site, Pembroke Street, Cambridge CB2 3RA, United Kingdom, and Department of Physics,Cavendish Laboratory, University of Cambridge, J J Thomson Avenue, Cambridge, CB3 0HE, United Kingdom
| | - Markus Kraft
- Department of Chemical Engineering and Biotechnology, University of Cambridge, New Museums Site, Pembroke Street, Cambridge CB2 3RA, United Kingdom, and Department of Physics,Cavendish Laboratory, University of Cambridge, J J Thomson Avenue, Cambridge, CB3 0HE, United Kingdom
| |
Collapse
|
40
|
Söderhjelm P, Kongsted J, Genheden S, Ryde U. Estimates of ligand-binding affinities supported by quantum mechanical methods. Interdiscip Sci 2010; 2:21-37. [DOI: 10.1007/s12539-010-0083-0] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2009] [Revised: 10/29/2009] [Accepted: 11/11/2009] [Indexed: 12/01/2022]
|
41
|
|
42
|
Misquitta AJ, Welch GW, Stone AJ, Price SL. A first principles prediction of the crystal structure of C6Br2ClFH2. Chem Phys Lett 2008. [DOI: 10.1016/j.cplett.2008.02.113] [Citation(s) in RCA: 73] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
43
|
Welch GWA, Karamertzanis PG, Misquitta AJ, Stone AJ, Price SL. Is the Induction Energy Important for Modeling Organic Crystals? J Chem Theory Comput 2008; 4:522-32. [DOI: 10.1021/ct700270d] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Gareth W. A. Welch
- Christopher Ingold Laboratory, Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, U.K., and University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, U.K
| | - Panagiotis G. Karamertzanis
- Christopher Ingold Laboratory, Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, U.K., and University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, U.K
| | - Alston J. Misquitta
- Christopher Ingold Laboratory, Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, U.K., and University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, U.K
| | - Anthony J. Stone
- Christopher Ingold Laboratory, Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, U.K., and University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, U.K
| | - Sarah L. Price
- Christopher Ingold Laboratory, Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, U.K., and University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, U.K
| |
Collapse
|
44
|
Söderhjelm P, Öhrn A, Ryde U, Karlström G. Accuracy of typical approximations in classical models of intermolecular polarization. J Chem Phys 2008; 128:014102. [DOI: 10.1063/1.2814240] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
|
45
|
Misquitta AJ, Stone AJ, Price SL. Accurate Induction Energies for Small Organic Molecules. 2. Development and Testing of Distributed Polarizability Models against SAPT(DFT) Energies. J Chem Theory Comput 2007; 4:19-32. [DOI: 10.1021/ct700105f] [Citation(s) in RCA: 71] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Alston J. Misquitta
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, U.K., and University College London, 20 Gordon Street, London WC1H 0AJ, U.K
| | - Anthony J. Stone
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, U.K., and University College London, 20 Gordon Street, London WC1H 0AJ, U.K
| | - Sarah L. Price
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, U.K., and University College London, 20 Gordon Street, London WC1H 0AJ, U.K
| |
Collapse
|
46
|
Lillestolen TC, Wheatley RJ. First-principles calculation of local atomic polarizabilities. J Phys Chem A 2007; 111:11141-6. [PMID: 17918806 DOI: 10.1021/jp073151y] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Common methods of determining atomic polarizabilities suffer from the inclusion of nonlocal effects such as charge polarization. A new method is described for determining fully ab initio atomic polarizabilities based on calculating the response of atomic multipoles to the local electrostatic potential. The localized atomic polarizabilities are then used to calculate induction energies that are compared to ab initio induction energies to test their usefulness in practical applications. These polarizabilities are shown to be an improvement over the corresponding molecular polarizabilities, in terms of both absolute accuracy and the convergence of the multipolar induction series. The transferability of localized polarizabilities for the alkane series is also discussed.
Collapse
Affiliation(s)
- T C Lillestolen
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, U.K.
| | | |
Collapse
|
47
|
Calculating intermolecular potentials with SIMPER: the water–nitrogen and water–oxygen interactions, dispersion energy coefficients, and preliminary results for larger molecules. INT REV PHYS CHEM 2007. [DOI: 10.1080/01442350701371539] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
48
|
Wheatley RJ, Tulegenov AS, Bichoutskaia * E. Intermolecular potentials from supermolecule and monomer calculations. INT REV PHYS CHEM 2007. [DOI: 10.1080/014423504200207772] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
| | - Akyl S. Tulegenov
- a School of Chemistry , University of Nottingham , Nottingham NG7 2RD, UK
| | | |
Collapse
|
49
|
|
50
|
Söderhjelm P, Krogh JW, Karlström G, Ryde U, Lindh R. Accuracy of distributed multipoles and polarizabilities: Comparison between the LoProp and MpProp models. J Comput Chem 2007; 28:1083-90. [PMID: 17279548 DOI: 10.1002/jcc.20632] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Localized multipole moments up to the fifth moment as well as localized dipole polarizabilities are calculated with the MpProp and the newly developed LoProp methods for a total of 20 molecules, predominantly derived from amino acids. A comparison of electrostatic potentials calculated from the multipole expansion obtained by the two methods with ab initio results shows that both methods reproduce the electrostatic interaction with an elementary charge with a mean absolute error of approximately 1.5 kJ/mol at contact distance and less than 0.1 kJ/mol at distances 2 A further out when terms up to the octupole moments are included. The polarizabilities are tested with homogenous electric fields and are found to have similar accuracy. The MpProp method gives better multipole moments unless diffuse basis sets are used, whereas LoProp gives better polarizabilities.
Collapse
Affiliation(s)
- P Söderhjelm
- Department of Theoretical Chemistry, Lund University, Chemical Center, P.O.B. 124, Lund 22100, Sweden
| | | | | | | | | |
Collapse
|