1
|
Salazar-Lozas H, Guevara-Vela JM, Pendás ÁM, Francisco E, Rocha-Rinza T. Partition of the electronic energy of the PM7 method via the interacting quantum atoms approach. Phys Chem Chem Phys 2022; 24:19521-19530. [PMID: 35938407 DOI: 10.1039/d2cp02013k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Partitions of the electronic energy such as that provided by the Interacting Quantum Atoms (IQA) approach have given valuable insights for numerous chemical systems and processes. Unfortunately, this kind of analysis may involve the integration of scalar fields over very irregular volumes, a condition which leads to a large and often prohibitive computational effort. These circumstances have limited the use of these energy partitions to systems comprising a few tens of atoms at most. On the other hand, semiempirical methods have proved useful in the study of systems of several thousands of atoms. Therefore, the goal of this work is to carry out partitions of the semiempirical method PM7 in compliance with the IQA approach. For this purpose, we computed one- and two-atomic energetic contributions whose sum equals the PM7 electronic energy. We illustrate how one might exploit the partition of electronic energies computed via the PM7 method by considering small organic and inorganic molecules and the energetics of individual hydrogen bond interactions within several water clusters which include (H2O)30, (H2O)50 and (H2O)100. We also considered the solvation of the amphiphilic caprylate anion to exemplify how to exploit the energy partition proposed in this paper. Overall, this investigation shows how the approach put forward herein might give further insights of the interactions occurring within complex systems in physical and biological chemistry.
Collapse
Affiliation(s)
- Hugo Salazar-Lozas
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, Alcaldía Coyoacán C.P. 04510, Ciudad de México, Mexico.
| | | | - Ángel Martín Pendás
- Departamento de Química Física y Analítica, Universidad de Oviedo, Julián Claveria 8, 33006, Oviedo, Spain
| | - Evelio Francisco
- Departamento de Química Física y Analítica, Universidad de Oviedo, Julián Claveria 8, 33006, Oviedo, Spain
| | - Tomás Rocha-Rinza
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, Alcaldía Coyoacán C.P. 04510, Ciudad de México, Mexico.
| |
Collapse
|
2
|
Ballesteros F, Lao KU. Accelerating the Convergence of Self-Consistent Field Calculations Using the Many-Body Expansion. J Chem Theory Comput 2021; 18:179-191. [PMID: 34881906 DOI: 10.1021/acs.jctc.1c00765] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The balance between cost-effective and sufficiently accurate methods represents the proverbial "promised land" for quantum chemistry calculations. The burden thus falls upon theoretical and computational chemists to provide such alternatives to mitigate the issues that arise from the employ of finite computing resources. In this paper, we attempt to demonstrate the importance of the quality of the initial guess for the self-consistent field (SCF) calculation when considering cost reduction techniques. We broach this challenge by using the many body expansion (MBE) to yield high quality density matrices (DMs) which, in turn, are applied as an SCF initial guess. The MBE-DM approaches combined with purification schemes and distance-based cutoff schemes can serve as initial guesses to reduce the SCF cycles necessary for convergence or derive energy directly through one Fock build. To this end, four unique types of clusters including water clusters, fluoride anion water clusters, sodium cation water clusters, and ammonium-bisulfate salt clusters have been used to test the performance of MBE-DM where its truncation at three-body expansion, MBE(3)-DM, shows vast improvement for those four clusters with reductions in the number of SCF cycles up to 40% as compared with the traditional superposition of atomic densities (SAD) guess. Other types of typical initial guesses, superposition of atomic potentials (SAP) and basis set projection (BSP), perform much worse than MBE-DM and SAD. In addition, the MBE-DM shows consistency across an array of fragment types irrespective of charges, size, level of theory, and basis set selection. Through MBE(3)-DM with the distance cutoff and the average purification scheme, the energy can be obtained directly with a mere 3.2 mH of the mean absolute deviation (MAD) for (H2O)N=6-55 which is at least 73 times better than the energy prediction using the typical initial guesses (SAD, SAP, and BSP). The corresponding MAD per monomer is only 0.14 mH which reaches the threshold of the "dynamical accuracy". The promising results of the methods outlined in this paper not only indicate two direct routes for computational cost reduction but also lay the possible foundation for composite techniques (i.e., ab initio sampling) that make best use of near converged values as their starting point.
Collapse
Affiliation(s)
- Francisco Ballesteros
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, United States
| | - Ka Un Lao
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, United States
| |
Collapse
|
3
|
Hughes ZE, Ren E, Thacker JCR, Symons BCB, Silva AF, Popelier PLA. A FFLUX Water Model: Flexible, Polarizable and with a Multipolar Description of Electrostatics. J Comput Chem 2020; 41:619-628. [PMID: 31747059 PMCID: PMC7004022 DOI: 10.1002/jcc.26111] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2019] [Revised: 10/21/2019] [Accepted: 10/31/2019] [Indexed: 12/15/2022]
Abstract
Key to progress in molecular simulation is the development of advanced models that go beyond the limitations of traditional force fields that employ a fixed, point charge-based description of electrostatics. Taking water as an example system, the FFLUX framework is shown capable of producing models that are flexible, polarizable and have a multipolar description of the electrostatics. The kriging machine-learning methods used in FFLUX are able to reproduce the intramolecular potential energy surface and multipole moments of a single water molecule with chemical accuracy using as few as 50 training configurations. Molecular dynamics simulations of water clusters (25-216 molecules) using the new FFLUX model reveal that incorporating charge-quadrupole, dipole-dipole, and quadrupole-charge interactions into the description of the electrostatics results in significant changes to the intermolecular structuring of the water molecules. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Zak E. Hughes
- Manchester Institute of Biotechnology, The University of ManchesterManchesterM1 7DNUnited Kingdom
- Department of ChemistryThe University of ManchesterManchesterM13 9PLUnited Kingdom
- School of Chemistry and Biosciences, University of BradfordBradfordBD7 1DPUnited Kingdom
| | - Emmanuel Ren
- Manchester Institute of Biotechnology, The University of ManchesterManchesterM1 7DNUnited Kingdom
- Department of ChemistryThe University of ManchesterManchesterM13 9PLUnited Kingdom
| | - Joseph C. R. Thacker
- Manchester Institute of Biotechnology, The University of ManchesterManchesterM1 7DNUnited Kingdom
- Department of ChemistryThe University of ManchesterManchesterM13 9PLUnited Kingdom
| | - Benjamin C. B. Symons
- Manchester Institute of Biotechnology, The University of ManchesterManchesterM1 7DNUnited Kingdom
- Department of ChemistryThe University of ManchesterManchesterM13 9PLUnited Kingdom
| | - Arnaldo F. Silva
- Manchester Institute of Biotechnology, The University of ManchesterManchesterM1 7DNUnited Kingdom
- Department of ChemistryThe University of ManchesterManchesterM13 9PLUnited Kingdom
| | - Paul L. A. Popelier
- Manchester Institute of Biotechnology, The University of ManchesterManchesterM1 7DNUnited Kingdom
- Department of ChemistryThe University of ManchesterManchesterM13 9PLUnited Kingdom
| |
Collapse
|
4
|
Petković M, Nakarada Đ, Etinski M. When hydroquinone meets methoxy radical: Hydrogen abstraction reaction from the viewpoint of interacting quantum atoms. J Comput Chem 2018; 39:1868-1877. [PMID: 29799128 DOI: 10.1002/jcc.25359] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2018] [Revised: 04/16/2018] [Accepted: 04/27/2018] [Indexed: 01/01/2023]
Abstract
Interacting Quantum Atoms methodology is used for a detailed analysis of hydrogen abstraction reaction from hydroquinone by methoxy radical. Two pathways are analyzed, which differ in the orientation of the reactants at the corresponding transition states. Although the discrepancy between the two barriers amounts to only 2 kJ/mol, which implies that the two pathways are of comparable probability, the extent of intra-atomic and inter-atomic energy changes differs considerably. We thus demonstrated that Interacting Quantum Atoms procedure can be applied to unravel distinct energy transfer routes in seemingly similar mechanisms. Identification of energy components with the greatest contribution to the variation of the overall energy (intra-atomic and inter-atomic terms that involve hydroquinone's oxygen and the carbon atom covalently bound to it, the transferring hydrogen and methoxy radical's oxygen), is performed using the Relative energy gradient method. Additionally, the Interacting Quantum Fragments approach shed light on the nature of dominant interactions among selected fragments: both Coulomb and exchange-correlation contributions are of comparable importance when considering interactions of the transferring hydrogen atom with all other atoms, whereas the exchange-correlation term dominates interaction between methoxy radical's methyl group and hydroquinone's aromatic ring. This study represents one of the first applications of Interacting Quantum Fragments approach on first order saddle points. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Milena Petković
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, Belgrade, 11 158, Serbia
| | - Đura Nakarada
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, Belgrade, 11 158, Serbia
| | - Mihajlo Etinski
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, Belgrade, 11 158, Serbia
| |
Collapse
|
5
|
McDonagh JL, Silva AF, Vincent MA, Popelier PLA. Machine Learning of Dynamic Electron Correlation Energies from Topological Atoms. J Chem Theory Comput 2017; 14:216-224. [PMID: 29211469 DOI: 10.1021/acs.jctc.7b01157] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
We present an innovative method for predicting the dynamic electron correlation energy of an atom or a bond in a molecule utilizing topological atoms. Our approach uses the machine learning method Kriging (Gaussian Process Regression with a non-zero mean function) to predict these dynamic electron correlation energy contributions. The true energy values are calculated by partitioning the MP2 two-particle density-matrix via the Interacting Quantum Atoms (IQA) procedure. To our knowledge, this is the first time such energies have been predicted by a machine learning technique. We present here three important proof-of-concept cases: the water monomer, the water dimer, and the van der Waals complex H2···He. These cases represent the final step toward the design of a full IQA potential for molecular simulation. This final piece will enable us to consider situations in which dispersion is the dominant intermolecular interaction. The results from these examples suggest a new method by which dispersion potentials for molecular simulation can be generated.
Collapse
Affiliation(s)
- James L McDonagh
- Manchester Institute of Biotechnology, The University of Manchester , 131 Princess Street, Manchester M1 7DN, Great Britain
| | - Arnaldo F Silva
- Manchester Institute of Biotechnology, The University of Manchester , 131 Princess Street, Manchester M1 7DN, Great Britain
| | - Mark A Vincent
- School of Chemistry, The University of Manchester , Oxford Road, Manchester M13 9PL, Great Britain
| | - Paul L A Popelier
- Manchester Institute of Biotechnology, The University of Manchester , 131 Princess Street, Manchester M1 7DN, Great Britain.,School of Chemistry, The University of Manchester , Oxford Road, Manchester M13 9PL, Great Britain
| |
Collapse
|
6
|
Maxwell PI, Popelier PLA. Unfavorable regions in the ramachandran plot: Is it really steric hindrance? The interacting quantum atoms perspective. J Comput Chem 2017; 38:2459-2474. [PMID: 28841241 PMCID: PMC5659141 DOI: 10.1002/jcc.24904] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2017] [Revised: 07/10/2017] [Accepted: 07/18/2017] [Indexed: 01/06/2023]
Abstract
Accurate description of the intrinsic preferences of amino acids is important to consider when developing a biomolecular force field. In this study, we use a modern energy partitioning approach called Interacting Quantum Atoms to inspect the cause of the φ and ψ torsional preferences of three dipeptides (Gly, Val, and Ile). Repeating energy trends at each of the molecular, functional group, and atomic levels are observed across both (1) the three amino acids and (2) the φ/ψ scans in Ramachandran plots. At the molecular level, it is surprisingly electrostatic destabilization that causes the high-energy regions in the Ramachandran plot, not molecular steric hindrance (related to the intra-atomic energy). At the functional group and atomic levels, the importance of key peptide atoms (Oi-1 , Ci , Ni , Ni+1 ) and some sidechain hydrogen atoms (Hγ ) are identified as responsible for the destabilization seen in the energetically disfavored Ramachandran regions. Consistently, the Oi-1 atoms are particularly important for the explanation of dipeptide intrinsic behavior, where electrostatic and steric destabilization unusually complement one another. The findings suggest that, at least for these dipeptides, it is the peptide group atoms that dominate the intrinsic behavior, more so than the sidechain atoms. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Peter I. Maxwell
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain and School of Chemistry, University of Manchester, Oxford RoadManchesterGreat BritainM13 9PL
| | - 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 RoadManchesterGreat BritainM13 9PL
| |
Collapse
|
7
|
Silva AF, Vincent MA, McDonagh JL, Popelier PLA. The Transferability of Topologically Partitioned Electron Correlation Energies in Water Clusters. Chemphyschem 2017; 18:3360-3368. [PMID: 29094804 DOI: 10.1002/cphc.201700890] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2017] [Revised: 09/21/2017] [Indexed: 11/06/2022]
Abstract
The electronic effects that govern the cohesion of water clusters are complex, demanding the inclusion of N-body, Coulomb, exchange and correlation effects. Here we present a much needed quantitative study of the effect of correlation (and hence dispersion) energy on the stabilization of water clusters. For this purpose we used a topological energy partitioning method called Interacting Quantum Atoms (IQA) to partition water clusters into topological atoms, based on a MP2/6-31G(d,p) wave function, and modified versions of GAUSSIAN09 and the Quantum Chemical Topology (QCT) program MORFI. Most of the cohesion in the water clusters provided by electron correlation comes from intramolecular energy stabilization. Hydrogen bond-related interactions tend to largely cancel each other. Electron correlation energies are transferable in almost all instances within 1 kcal mol-1 . This observed transferability is very important to the further development of the QCT force field FFLUX, especially to the future modelling of liquid water.
Collapse
Affiliation(s)
- Arnaldo F Silva
- Manchester Institute of Biotechnology (MIB), the University of Manchester, 131 Princess Street, Manchester, M1 7DN, Great Britain), Fax: (+44) 161 3064511.,School of Chemistry, the University of Manchester, Oxford Road, Manchester, M13 9PL, Great Britain
| | - Mark A Vincent
- Manchester Institute of Biotechnology (MIB), the University of Manchester, 131 Princess Street, Manchester, M1 7DN, Great Britain), Fax: (+44) 161 3064511.,School of Chemistry, the University of Manchester, Oxford Road, Manchester, M13 9PL, Great Britain
| | - James L McDonagh
- Manchester Institute of Biotechnology (MIB), the University of Manchester, 131 Princess Street, Manchester, M1 7DN, Great Britain), Fax: (+44) 161 3064511.,School of Chemistry, the University of Manchester, Oxford Road, Manchester, M13 9PL, Great Britain
| | - Paul L A Popelier
- Manchester Institute of Biotechnology (MIB), the University of Manchester, 131 Princess Street, Manchester, M1 7DN, Great Britain), Fax: (+44) 161 3064511.,School of Chemistry, the University of Manchester, Oxford Road, Manchester, M13 9PL, Great Britain
| |
Collapse
|
8
|
Maxwell PI, Popelier PLA. Accurate prediction of the energetics of weakly bound complexes using the machine learning method kriging. Struct Chem 2017. [DOI: 10.1007/s11224-017-0928-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|