1
|
Liu Y, Guo H. A Gaussian Process Based Δ-Machine Learning Approach to Reactive Potential Energy Surfaces. J Phys Chem A 2023; 127:8765-8772. [PMID: 37815868 DOI: 10.1021/acs.jpca.3c05318] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/12/2023]
Abstract
The Gaussian process (GP) is an efficient nonparametric machine learning (ML) method. A distinct advantage of the GP is its ability to provide an estimate of statistical uncertainties. This is particularly useful in constructing high-dimensional potential energy surfaces (PESs) from ab initio data as it offers an optimal way to add new geometries to reduce the overall error. In this work, GP is employed in the context of Δ-machine learning (Δ-ML), in which a correction PES to an inaccurate low-level PES is constructed using a small number of high-level ab initio calculations. This new method is tested in three prototypical reactive systems, namely, the H + H2 → H + H2, OH + H2 → H2O + H, and H + CH4 → H2 + CH3 reactions. The results show that the GP-based Δ-ML approach is more efficient than its direct application in constructing high-level PESs. We also compare the new method to a previously proposed neural-network-based Δ-ML approach [Liu and Li J. Phys. Chem. Lett. 2022, 13, 4729-4738]. The results indicate that the two Δ-ML methods have comparable efficiencies in constructing accurate PESs.
Collapse
Affiliation(s)
- Yang Liu
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
2
|
Clark T. How deeply should we analyze non-covalent interactions? J Mol Model 2023; 29:66. [PMID: 36757533 PMCID: PMC9911493 DOI: 10.1007/s00894-023-05460-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Accepted: 01/18/2023] [Indexed: 02/10/2023]
Abstract
CONTEXT Just how much effort and detail should we invest in analyzing interactions of the order of 5 kcal mol-1? This comment attempts to provide a conciliatory overview of what is often a contentious field and to pose some questions that I hope will eventually lead at least to some consensus. METHODS This is an opinion article without calculations or data.
Collapse
Affiliation(s)
- Timothy Clark
- Computer-Chemistry-Center, Department of Chemistry and Pharmacy, Friedrich-Alexander-University Erlangen-Nuernberg, Naegelsbachstrasse 25, 91052, Erlangen, Germany.
| |
Collapse
|
3
|
Non-covalent interactions from a Quantum Chemical Topology perspective. J Mol Model 2022; 28:276. [PMID: 36006513 PMCID: PMC9411098 DOI: 10.1007/s00894-022-05188-7] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Accepted: 02/07/2022] [Indexed: 11/12/2022]
Abstract
About half a century after its little-known beginnings, the quantum topological approach called QTAIM has grown into a widespread, but still not mainstream, methodology of interpretational quantum chemistry. Although often confused in textbooks with yet another population analysis, be it perhaps an elegant but somewhat esoteric one, QTAIM has been enriched with about a dozen other research areas sharing its main mathematical language, such as Interacting Quantum Atoms (IQA) or Electron Localisation Function (ELF), to form an overarching approach called Quantum Chemical Topology (QCT). Instead of reviewing the latter’s role in understanding non-covalent interactions, we propose a number of ideas emerging from the full consequences of the space-filling nature of topological atoms, and discuss how they (will) impact on interatomic interactions, including non-covalent ones. The architecture of a force field called FFLUX, which is based on these ideas, is outlined. A new method called Relative Energy Gradient (REG) is put forward, which is able, by computation, to detect which fragments of a given molecular assembly govern the energetic behaviour of this whole assembly. This method can offer insight into the typical balance of competing atomic energies both in covalent and non-covalent case studies. A brief discussion on so-called bond critical points is given, highlighting concerns about their meaning, mainly in the arena of non-covalent interactions.
Collapse
|
4
|
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
|
5
|
Sauza-de la Vega A, Rocha-Rinza T, Guevara-Vela JM. Cooperativity and Anticooperativity in Ion-Water Interactions: Implications for the Aqueous Solvation of Ions. Chemphyschem 2021; 22:1269-1285. [PMID: 33635563 DOI: 10.1002/cphc.202000981] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2020] [Revised: 02/24/2021] [Indexed: 01/03/2023]
Abstract
Non-additive effects in hydrogen bonds (HB) take place as a consequence of electronic charge transfers. Therefore, it is natural to expect cooperativity and anticooperativity in ion-water interactions. Nevertheless, investigations on this matter are scarce. This paper addresses the interactions of (i) the cations Li+ , Na+ , K+ , Be2+ , Mg2+ , and Ca2+ together with (ii) the anions F- , Cl- , Br- , NO3 - and SO4 2- with water clusters (H2 O)n , n=1-8, and the effects of these ions on the HBs within the complete molecular adducts. We used quantum chemical topology tools, specifically the quantum theory of atoms in molecules and the interacting quantum atoms energy partition to investigate non-additive effects among the interactions studied herein. Our results show a decrease on the interaction energy between ions and the first neighbouring water molecules with an increment of the coordination number. We also found strong cooperative effects in the interplay between HBs and ion-dipole interactions within the studied systems. Such cooperativity affects considerably the interactions among ions with their first and second solvation shells in aqueous environments. Overall, we believe this article provides valuable information about how ion-dipole contacts interact with each other and how they relate to other interactions, such as HBs, in the framework of non-additive effects in aqueous media.
Collapse
Affiliation(s)
- Arturo Sauza-de la Vega
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, Delegación Coyoacán C.P., 04510, CDMX, México
| | - Tomás Rocha-Rinza
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, Delegación Coyoacán C.P., 04510, CDMX, México
| | - José Manuel Guevara-Vela
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, Delegación Coyoacán C.P., 04510, CDMX, México
| |
Collapse
|
6
|
Di Pasquale N, Davie SJ, Popelier PLA. The accuracy of ab initio calculations without ab initio calculations for charged systems: Kriging predictions of atomistic properties for ions in aqueous solutions. J Chem Phys 2018; 148:241724. [DOI: 10.1063/1.5022174] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Affiliation(s)
- Nicodemo Di Pasquale
- Manchester Institute of Biotechnology (MIB), 131 Princess
Street, Manchester M1 7DN, United Kingdom and School of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL,
United Kingdom
| | - Stuart J. Davie
- Manchester Institute of Biotechnology (MIB), 131 Princess
Street, Manchester M1 7DN, United Kingdom and School of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL,
United Kingdom
| | - Paul L. A. Popelier
- Manchester Institute of Biotechnology (MIB), 131 Princess
Street, Manchester M1 7DN, United Kingdom and School of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL,
United Kingdom
| |
Collapse
|
7
|
Uteva E, Graham RS, Wilkinson RD, Wheatley RJ. Interpolation of intermolecular potentials using Gaussian processes. J Chem Phys 2017; 147:161706. [DOI: 10.1063/1.4986489] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
|
8
|
Geometry Optimization with Machine Trained Topological Atoms. Sci Rep 2017; 7:12817. [PMID: 28993674 PMCID: PMC5634454 DOI: 10.1038/s41598-017-12600-3] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2017] [Accepted: 09/06/2017] [Indexed: 11/19/2022] Open
Abstract
The geometry optimization of a water molecule with a novel type of energy function called FFLUX is presented, which bypasses the traditional bonded potentials. Instead, topologically-partitioned atomic energies are trained by the machine learning method kriging to predict their IQA atomic energies for a previously unseen molecular geometry. Proof-of-concept that FFLUX’s architecture is suitable for geometry optimization is rigorously demonstrated. It is found that accurate kriging models can optimize 2000 distorted geometries to within 0.28 kJ mol−1 of the corresponding ab initio energy, and 50% of those to within 0.05 kJ mol−1. Kriging models are robust enough to optimize the molecular geometry to sub-noise accuracy, when two thirds of the geometric inputs are outside the training range of that model. Finally, the individual components of the potential energy are analyzed, and chemical intuition is reflected in the independent behavior of the three energy terms \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${E}_{{\rm{intra}}}^{{\rm{A}}}$$\end{document}EintraA(intra-atomic), \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${V}_{{\rm{cl}}}^{\text{AA}\text{'}}$$\end{document}VclAA' (electrostatic) and \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${V}_{{\rm{x}}}^{\text{AA}\text{'}}$$\end{document}VxAA' (exchange), in contrast to standard force fields.
Collapse
|
9
|
Kolb B, Marshall P, Zhao B, Jiang B, Guo H. Representing Global Reactive Potential Energy Surfaces Using Gaussian Processes. J Phys Chem A 2017; 121:2552-2557. [DOI: 10.1021/acs.jpca.7b01182] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Brian Kolb
- Department
of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
- Department
of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Paul Marshall
- Department
of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Bin Zhao
- Department
of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Bin Jiang
- Department
of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Hua Guo
- Department
of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
10
|
Davie SJ, Di Pasquale N, Popelier PLA. Incorporation of local structure into kriging models for the prediction of atomistic properties in the water decamer. J Comput Chem 2016; 37:2409-22. [PMID: 27535711 PMCID: PMC5031213 DOI: 10.1002/jcc.24465] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2016] [Revised: 06/06/2016] [Accepted: 07/07/2016] [Indexed: 01/17/2023]
Abstract
Machine learning algorithms have been demonstrated to predict atomistic properties approaching the accuracy of quantum chemical calculations at significantly less computational cost. Difficulties arise, however, when attempting to apply these techniques to large systems, or systems possessing excessive conformational freedom. In this article, the machine learning method kriging is applied to predict both the intra-atomic and interatomic energies, as well as the electrostatic multipole moments, of the atoms of a water molecule at the center of a 10 water molecule (decamer) cluster. Unlike previous work, where the properties of small water clusters were predicted using a molecular local frame, and where training set inputs (features) were based on atomic index, a variety of feature definitions and coordinate frames are considered here to increase prediction accuracy. It is shown that, for a water molecule at the center of a decamer, no single method of defining features or coordinate schemes is optimal for every property. However, explicitly accounting for the structure of the first solvation shell in the definition of the features of the kriging training set, and centring the coordinate frame on the atom-of-interest will, in general, return better predictions than models that apply the standard methods of feature definition, or a molecular coordinate frame. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Stuart J Davie
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain and School of Chemistry, University of Manchester, Oxford Road, Manchester, M13 9PL, Great Britain
| | - Nicodemo Di Pasquale
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain and School of Chemistry, University of Manchester, Oxford Road, Manchester, M13 9PL, Great Britain
| | - Paul L A Popelier
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain and School of Chemistry, University of Manchester, Oxford Road, Manchester, M13 9PL, Great Britain.
| |
Collapse
|
11
|
McDonagh JL, Vincent MA, Popelier PL. Partitioning dynamic electron correlation energy: Viewing Møller-Plesset correlation energies through Interacting Quantum Atom (IQA) energy partitioning. Chem Phys Lett 2016. [DOI: 10.1016/j.cplett.2016.09.019] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
12
|
Fletcher TL, Popelier PLA. Multipolar Electrostatic Energy Prediction for all 20 Natural Amino Acids Using Kriging Machine Learning. J Chem Theory Comput 2016; 12:2742-51. [DOI: 10.1021/acs.jctc.6b00457] [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)
- Timothy L. Fletcher
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain
- School
of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| | - Paul L. A. Popelier
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain
- School
of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| |
Collapse
|
13
|
Di Pasquale N, Davie SJ, Popelier PLA. Optimization Algorithms in Optimal Predictions of Atomistic Properties by Kriging. J Chem Theory Comput 2016; 12:1499-513. [DOI: 10.1021/acs.jctc.5b00936] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Nicodemo Di Pasquale
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain
- School of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| | - Stuart J. Davie
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain
- School of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| | - Paul L. A. Popelier
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain
- School of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| |
Collapse
|
14
|
On Quantum Chemical Topology. CHALLENGES AND ADVANCES IN COMPUTATIONAL CHEMISTRY AND PHYSICS 2016. [DOI: 10.1007/978-3-319-29022-5_2] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
|
15
|
Fletcher TL, Davie SJ, Popelier PLA. Prediction of Intramolecular Polarization of Aromatic Amino Acids Using Kriging Machine Learning. J Chem Theory Comput 2015; 10:3708-19. [PMID: 26588516 DOI: 10.1021/ct500416k] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Present computing power enables novel ways of modeling polarization. Here we show that the machine learning method kriging accurately captures the way the electron density of a topological atom responds to a change in the positions of the surrounding atoms. The success of this method is demonstrated on the four aromatic amino acids histidine, phenylalanine, tryptophan, and tyrosine. A new technique of varying training set sizes to vastly reduce training times while maintaining accuracy is described and applied to each amino acid. Each amino acid has its geometry distorted via normal modes of vibration over all local energy minima in the Ramachandran map. These geometries are then used to train the kriging models. Total electrostatic energies predicted by the kriging models for previously unseen geometries are compared to the true energies, yielding mean absolute errors of 2.9, 5.1, 4.2, and 2.8 kJ mol(-1) for histidine, phenylalanine, tryptophan, and tyrosine, respectively.
Collapse
Affiliation(s)
- Timothy L Fletcher
- Manchester Institute of Biotechnology (MIB) , 131 Princess Street, Manchester, M1 7DN, Great Britain.,School of Chemistry, University of Manchester , Oxford Road, Manchester, M13 9PL, Great Britain
| | - Stuart J Davie
- Manchester Institute of Biotechnology (MIB) , 131 Princess Street, Manchester, M1 7DN, Great Britain.,School of Chemistry, University of Manchester , Oxford Road, Manchester, M13 9PL, Great Britain
| | - Paul L A Popelier
- Manchester Institute of Biotechnology (MIB) , 131 Princess Street, Manchester, M1 7DN, Great Britain.,School of Chemistry, University of Manchester , Oxford Road, Manchester, M13 9PL, Great Britain
| |
Collapse
|
16
|
Transferable kriging machine learning models for the multipolar electrostatics of helical deca-alanine. Theor Chem Acc 2015. [DOI: 10.1007/s00214-015-1739-y] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
17
|
Hughes TJ, Cardamone S, Popelier PLA. Realistic sampling of amino acid geometries for a multipolar polarizable force field. J Comput Chem 2015; 36:1844-57. [PMID: 26235784 PMCID: PMC4973712 DOI: 10.1002/jcc.24006] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2015] [Revised: 06/19/2015] [Accepted: 06/20/2015] [Indexed: 12/19/2022]
Abstract
The Quantum Chemical Topological Force Field (QCTFF) uses the machine learning method kriging to map atomic multipole moments to the coordinates of all atoms in the molecular system. It is important that kriging operates on relevant and realistic training sets of molecular geometries. Therefore, we sampled single amino acid geometries directly from protein crystal structures stored in the Protein Databank (PDB). This sampling enhances the conformational realism (in terms of dihedral angles) of the training geometries. However, these geometries can be fraught with inaccurate bond lengths and valence angles due to artefacts of the refinement process of the X-ray diffraction patterns, combined with experimentally invisible hydrogen atoms. This is why we developed a hybrid PDB/nonstationary normal modes (NM) sampling approach called PDB/NM. This method is superior over standard NM sampling, which captures only geometries optimized from the stationary points of single amino acids in the gas phase. Indeed, PDB/NM combines the sampling of relevant dihedral angles with chemically correct local geometries. Geometries sampled using PDB/NM were used to build kriging models for alanine and lysine, and their prediction accuracy was compared to models built from geometries sampled from three other sampling approaches. Bond length variation, as opposed to variation in dihedral angles, puts pressure on prediction accuracy, potentially lowering it. Hence, the larger coverage of dihedral angles of the PDB/NM method does not deteriorate the predictive accuracy of kriging models, compared to the NM sampling around local energetic minima used so far in the development of QCTFF.
Collapse
Affiliation(s)
- Timothy J Hughes
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester, M1 7DN, Great Britain
- School of Chemistry, University of Manchester, Oxford Road, Manchester, M13 9PL, Great Britain
| | - Salvatore Cardamone
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester, M1 7DN, Great Britain
- School of Chemistry, University of Manchester, Oxford Road, Manchester, M13 9PL, Great Britain
| | - Paul L A Popelier
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester, M1 7DN, Great Britain
- School of Chemistry, University of Manchester, Oxford Road, Manchester, M13 9PL, Great Britain
| |
Collapse
|
18
|
Mills MJL, Popelier PLA. Electrostatic Forces: Formulas for the First Derivatives of a Polarizable, Anisotropic Electrostatic Potential Energy Function Based on Machine Learning. J Chem Theory Comput 2014; 10:3840-56. [DOI: 10.1021/ct500565g] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Matthew J. L. Mills
- Manchester
Institute of Biotechnology (MIB), University of Manchester, 131 Princess
Street, Manchester M1 7DN, Great Britain
| | - Paul L. A. Popelier
- School
of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| |
Collapse
|
19
|
Spherical tensor multipolar electrostatics and smooth particle mesh Ewald summation: a theoretical study. J Mol Model 2014; 20:2256. [DOI: 10.1007/s00894-014-2256-y] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2013] [Accepted: 04/22/2014] [Indexed: 10/25/2022]
|
20
|
Cisneros GA, Karttunen M, Ren P, Sagui C. Correction to Classical Electrostatics for Biomolecular Simulations. Chem Rev 2014. [DOI: 10.1021/cr500124k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
21
|
Multipolar electrostatics based on the Kriging machine learning method: an application to serine. J Mol Model 2014; 20:2172. [DOI: 10.1007/s00894-014-2172-1] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2013] [Accepted: 02/07/2014] [Indexed: 10/25/2022]
|
22
|
Yuan Y, Mills MJL, Popelier PLA. Multipolar electrostatics for proteins: Atom-atom electrostatic energies in crambin. J Comput Chem 2013; 35:343-59. [DOI: 10.1002/jcc.23469] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2013] [Revised: 09/04/2013] [Accepted: 09/29/2013] [Indexed: 12/31/2022]
Affiliation(s)
- Yongna Yuan
- Manchester Institute of Biotechnology (MIB); 131 Princess Street Manchester M1 7DN United Kingdom
- School of Chemistry, University of Manchester; Oxford Road Manchester M13 9PL United Kingdom
| | - Matthew J. L. Mills
- Manchester Institute of Biotechnology (MIB); 131 Princess Street Manchester M1 7DN United Kingdom
- School of Chemistry, University of Manchester; Oxford Road Manchester M13 9PL United Kingdom
| | - Paul L. A. Popelier
- Manchester Institute of Biotechnology (MIB); 131 Princess Street Manchester M1 7DN United Kingdom
- School of Chemistry, University of Manchester; Oxford Road Manchester M13 9PL United Kingdom
| |
Collapse
|