1
|
Spencer RJ, Zhanserkeev AA, Yang EL, Steele RP. The Near-Sightedness of Many-Body Interactions in Anharmonic Vibrational Couplings. J Am Chem Soc 2024; 146:15376-15392. [PMID: 38771156 DOI: 10.1021/jacs.4c03198] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
Couplings between vibrational motions are driven by electronic interactions, and these couplings carry special significance in vibrational energy transfer, multidimensional spectroscopy experiments, and simulations of vibrational spectra. In this investigation, the many-body contributions to these couplings are analyzed computationally in the context of clathrate-like alkali metal cation hydrates, including Cs+(H2O)20, Rb+(H2O)20, and K+(H2O)20, using both analytic and quantum-chemistry potential energy surfaces. Although the harmonic spectra and one-dimensional anharmonic spectra depend strongly on these many-body interactions, the mode-pair couplings were, perhaps surprisingly, found to be dominated by one-body effects, even in cases of couplings to low-frequency modes that involved the motion of multiple water molecules. The origin of this effect was traced mainly to geometric distortion within water monomers and cancellation of many-body effects in differential couplings, and the effect was also shown to be agnostic to the identity of the ion. These outcomes provide new understanding of vibrational couplings and suggest the possibility of improved computational methods for the simulation of infrared and Raman spectra.
Collapse
Affiliation(s)
- Ryan J Spencer
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Asylbek A Zhanserkeev
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Emily L Yang
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Ryan P Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| |
Collapse
|
2
|
Sattasathuchana T, Xu P, Bertoni C, Kim YL, Leang SS, Pham BQ, Gordon MS. The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems. J Chem Theory Comput 2024; 20:2445-2461. [PMID: 38450638 DOI: 10.1021/acs.jctc.3c01309] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/08/2024]
Abstract
The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree-Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. The accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, Rcut. Rcut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small Rcut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.
Collapse
Affiliation(s)
- Tosaporn Sattasathuchana
- Department of Chemistry, Iowa State University and Ames National Laboratory, Ames, Iowa 50011, United States
| | - Peng Xu
- Department of Chemistry, Iowa State University and Ames National Laboratory, Ames, Iowa 50011, United States
| | - Colleen Bertoni
- Argonne Leadership Computing Facility, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Yu Lim Kim
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Sarom S Leang
- EP Analytics, Inc., 9909 Mira Mesa Blvd Ste. 230, San Diego, California 92131, United States
| | - Buu Q Pham
- Department of Chemistry, Iowa State University and Ames National Laboratory, Ames, Iowa 50011, United States
| | - Mark S Gordon
- Department of Chemistry, Iowa State University and Ames National Laboratory, Ames, Iowa 50011, United States
| |
Collapse
|
3
|
Zahariev F, Xu P, Westheimer BM, Webb S, Galvez Vallejo J, Tiwari A, Sundriyal V, Sosonkina M, Shen J, Schoendorff G, Schlinsog M, Sattasathuchana T, Ruedenberg K, Roskop LB, Rendell AP, Poole D, Piecuch P, Pham BQ, Mironov V, Mato J, Leonard S, Leang SS, Ivanic J, Hayes J, Harville T, Gururangan K, Guidez E, Gerasimov IS, Friedl C, Ferreras KN, Elliott G, Datta D, Cruz DDA, Carrington L, Bertoni C, Barca GMJ, Alkan M, Gordon MS. The General Atomic and Molecular Electronic Structure System (GAMESS): Novel Methods on Novel Architectures. J Chem Theory Comput 2023; 19:7031-7055. [PMID: 37793073 DOI: 10.1021/acs.jctc.3c00379] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/06/2023]
Abstract
The primary focus of GAMESS over the last 5 years has been the development of new high-performance codes that are able to take effective and efficient advantage of the most advanced computer architectures, both CPU and accelerators. These efforts include employing density fitting and fragmentation methods to reduce the high scaling of well-correlated (e.g., coupled-cluster) methods as well as developing novel codes that can take optimal advantage of graphical processing units and other modern accelerators. Because accurate wave functions can be very complex, an important new functionality in GAMESS is the quasi-atomic orbital analysis, an unbiased approach to the understanding of covalent bonds embedded in the wave function. Best practices for the maintenance and distribution of GAMESS are also discussed.
Collapse
Affiliation(s)
- Federico Zahariev
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Peng Xu
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Bryce M Westheimer
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Simon Webb
- VeraChem LLC, 12850 Middlebrook Road, Suite 205, Germantown, Maryland 20874-5244, United States
| | - Jorge Galvez Vallejo
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
- Research School of Computer Science, Australian National University, Canberra, ACT 2601, Australia
| | - Ananta Tiwari
- EP Analytics, Inc., 9909 Mira Mesa Boulevard, Suite 230, San Diego, California 92131, United States
| | - Vaibhav Sundriyal
- Department of Computational Modeling and Simulation Engineering, Old Dominion University, Norfolk, Virginia 23529, United States
| | - Masha Sosonkina
- Department of Computational Modeling and Simulation Engineering, Old Dominion University, Norfolk, Virginia 23529, United States
| | - Jun Shen
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - George Schoendorff
- Propellants Branch, Rocket Propulsion Division, Aerospace Systems Directorate, Air Force Research Laboratory, AFRL/RQRP, Edwards Air Force Base, California 93524, United States
| | - Megan Schlinsog
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Tosaporn Sattasathuchana
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Klaus Ruedenberg
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Luke B Roskop
- Hewlett-Packard Enterprise, 2131 Lindau Lane #1000, Bloomington, Minnesota 55425, United States
| | | | - David Poole
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
- School of Chemistry & Biochemistry, Georgia Institute of Technology, Athens, Georgia 30332, United States
| | - Piotr Piecuch
- Department of Chemistry and Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, United States
| | - Buu Q Pham
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Vladimir Mironov
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Joani Mato
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, United States
| | - Sam Leonard
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Sarom S Leang
- EP Analytics, Inc., 9909 Mira Mesa Boulevard, Suite 230, San Diego, California 92131, United States
| | - Joe Ivanic
- Advanced Biomedical Computational Science, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21702, United States
| | - Jackson Hayes
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Taylor Harville
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Karthik Gururangan
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Emilie Guidez
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217, United States
| | - Igor S Gerasimov
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Christian Friedl
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Katherine N Ferreras
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - George Elliott
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Dipayan Datta
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Daniel Del Angel Cruz
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Laura Carrington
- EP Analytics, Inc., 9909 Mira Mesa Boulevard, Suite 230, San Diego, California 92131, United States
| | - Colleen Bertoni
- Argonne Leadership Computing Facility, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Giuseppe M J Barca
- Research School of Computer Science, Australian National University, Canberra, ACT 2601, Australia
| | - Melisa Alkan
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| | - Mark S Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50014, United States
| |
Collapse
|
4
|
Barcza B, Szirmai Á, Tajti A, Stanton JF, Szalay PG. Benchmarking Aspects of Ab Initio Fragment Models for Accurate Excimer Potential Energy Surfaces. J Chem Theory Comput 2023; 19:3580-3600. [PMID: 37236166 PMCID: PMC10694823 DOI: 10.1021/acs.jctc.3c00104] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2023] [Indexed: 05/28/2023]
Abstract
While Coupled-Cluster methods have been proven to provide an accurate description of excited electronic states, the scaling of the computational costs with the system size limits the degree for which these methods can be applied. In this work different aspects of fragment-based approaches are studied on noncovalently bound molecular complexes with interacting chromophores of the fragments, such as π-stacked nucleobases. The interaction of the fragments is considered at two distinct steps. First, the states localized on the fragments are described in the presence of the other fragment(s); for this we test two approaches. One method is founded on QM/MM principles, only including the electrostatic interaction between the fragments in the electronic structure calculation with Pauli repulsion and dispersion effects added separately. The other model, a Projection-based Embedding (PbE) using the Huzinaga equation, includes both electrostatic and Pauli repulsion and only needs to be augmented by dispersion interactions. In both schemes the extended Effective Fragment Potential (EFP2) method of Gordon et al. was found to provide an adequate correction for the missing terms. In the second step, the interaction of the localized chromophores is modeled for a proper description of the excitonic coupling. Here the inclusion of purely electrostatic contributions appears to be sufficient: it is found that the Coulomb part of the coupling provides accurate splitting of the energies of interacting chromophores that are separated by more than 4 Å.
Collapse
Affiliation(s)
- Bónis Barcza
- Laboratory
of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1117 Budapest, Hungary
- György
Hevesy Doctoral School, Institute of Chemistry, ELTE Eötvös Loránd University, H-1117 Budapest, Hungary
| | - Ádám
B. Szirmai
- Laboratory
of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1117 Budapest, Hungary
- György
Hevesy Doctoral School, Institute of Chemistry, ELTE Eötvös Loránd University, H-1117 Budapest, Hungary
| | - Attila Tajti
- Laboratory
of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1117 Budapest, Hungary
| | - John F. Stanton
- Quantum
Theory Project, Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Péter G. Szalay
- Laboratory
of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1117 Budapest, Hungary
| |
Collapse
|
5
|
Barcza B, Szirmai ÁB, Szántó KJ, Tajti A, Szalay PG. Comparison of approximate intermolecular potentials for ab initio fragment calculations on medium sized N-heterocycles. J Comput Chem 2022; 43:1079-1093. [PMID: 35478353 PMCID: PMC9321956 DOI: 10.1002/jcc.26866] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Revised: 03/26/2022] [Accepted: 03/29/2022] [Indexed: 01/15/2023]
Abstract
The ground state intermolecular potential of bimolecular complexes of N‐heterocycles is analyzed for the impact of individual terms in the interaction energy as provided by various, conceptually different theories. Novel combinations with several formulations of the electrostatic, Pauli repulsion, and dispersion contributions are tested at both short‐ and long‐distance sides of the potential energy surface, for various alignments of the pyrrole dimer as well as the cytosine–uracil complex. The integration of a DFT/CCSD density embedding scheme, with dispersion terms from the effective fragment potential (EFP) method is found to provide good agreement with a reference CCSD(T) potential overall; simultaneously, a quantum mechanics/molecular mechanics approach using CHELPG atomic point charges for the electrostatic interaction, augmented by EFP dispersion and Pauli repulsion, comes also close to the reference result. Both schemes have the advantage of not relying on predefined force fields; rather, the interaction parameters can be determined for the system under study, thus being excellent candidates for ab initio modeling.
Collapse
Affiliation(s)
- Bónis Barcza
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| | - Ádám B Szirmai
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| | - Katalin J Szántó
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| | - Attila Tajti
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| | - Péter G Szalay
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| |
Collapse
|
6
|
Khazhiev S, Khusainov M, Khalikov R, Kataev V, Tyumkina T, Mescheryakova E, Khalilov L, Kuznetsov V. Structure and conformational analysis of 5,5-bis(bromomethyl)-2-trichloromethyl-1,3-dioxane by XRD, NMR and computer simulation. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2021.132326] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
7
|
Carter-Fenk K, Herbert JM. Appraisal of dispersion damping functions for the effective fragment potential method. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2055504] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Kevin Carter-Fenk
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH, USA
| |
Collapse
|
8
|
Loco D, Lagardère L, Adjoua O, Piquemal JP. Atomistic Polarizable Embeddings: Energy, Dynamics, Spectroscopy, and Reactivity. Acc Chem Res 2021; 54:2812-2822. [PMID: 33961401 PMCID: PMC8264944 DOI: 10.1021/acs.accounts.0c00662] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Indexed: 12/20/2022]
Abstract
The computational modeling of realistic extended systems, relevant in, e.g., Chemistry and Biophysics, is a fundamental problem of paramount importance in contemporary research. Enzymatic catalysis and photoinduced processes in pigment-protein complexes are typical problems targeted by computer-aided approaches, to complement experiments as interpretative tools at a molecular scale. The daunting complexity of this task lies in between the opposite stringent requirements of results' reliability for structural/dynamical properties and related intermolecular interactions, and a mandatory principle of realism in the modeling strategy. Therefore, in practice, a truly realistic computational model of a biologically relevant system can easily fail to meet the accuracy requirement, in order to balance the excessive computational cost necessary to reach the desired precision.To address such an "accuracy vs reality" dualistic requirement, mixed quantum mechanics/classical mechanics approaches within Atomistic (i.e., preserving the discrete particle configuration) Polarizable Embeddings (QM/APEs) methods have been proposed over the years. In this Account, we review recent developments in the design and application of general QM/APE methods, targeting situations where a local intrinsically quantum behavior is coupled to a large molecular system (i.e., an environment), often involving processes with different dynamical time scales, in order to avoid brute-force, unpractical quantum chemistry calculations on the complete system.In the first place, our interest is devoted to the available APEs models presently implemented in computational software, highlighting the quantum chemistry methods that can be used to treat the QM subsystem. We review the coupling strategy between the QM subsystem and the APE, which requires to examine the way the QM/MM mutual interactions are accounted for and how the polarization of the classical environment is considered with respect to (wrt) the quantum variables. Because of the need of reliable molecular and macromolecular structures, a pivotal aspect to address here is the handling of the system dynamics (i.e., gradients wrt nuclear positions are required), especially for large molecular assemblies composed by an overwhelming number of atoms, exploring many conformations on a complex energy landscape.Alongside, we highlight our views on the necessary steps to take toward more accurate general-purposes and transferable explicit embeddings. The main objective to achieve here is to design a more physically grounded multiscale approach. To do so, one should apply advanced new generation classical models to account for refined induction effects that are able to (i) improve the quality of QM/MM interaction energies; (ii) enhance transferability by avoiding the compulsory partial (or total) reparameterization of the classical model. Moreover, the extension of recent developments originating from the field of advanced classical molecular dynamics (MD) to the realm of QM/APE methods is a key direction to improve both speed and efficiency for the phase space exploration of systems of growing size and complexity.Lastly, we point out specific research topics where an advanced QM/APE dynamics can certainly shed some light. For example, we discuss chemical reactions in "harsh" environments and the case of spectroscopic theoretical modeling where the inclusion of refined environment effects is often mandatory.
Collapse
Affiliation(s)
- Daniele Loco
- Laboratoire
de Chimie Théorique, Sorbonne Université,
UMR 7616 CNRS, 75005 Paris, France
| | - Louis Lagardère
- Laboratoire
de Chimie Théorique, Sorbonne Université,
UMR 7616 CNRS, 75005 Paris, France
- Intitut
Parisien de Chimie Physique et Théorique, Sorbonne Université, FR 2622 CNRS, 75005 Paris, France
| | - Olivier Adjoua
- Laboratoire
de Chimie Théorique, Sorbonne Université,
UMR 7616 CNRS, 75005 Paris, France
| | - Jean-Philip Piquemal
- Laboratoire
de Chimie Théorique, Sorbonne Université,
UMR 7616 CNRS, 75005 Paris, France
- Institut
Universitaire de France, F-75005 Paris, France
- Department
of Biomedical Engineering, The University
of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
9
|
Kataev VA, Mesheryakova SA, Mesheryakova ES, Tyumkina TV, Khalilov LM, Lazarev VV, Kuznetsov VV. Synthesis, Structure, and Conformational Analysis of N-(2,4-Dichlorophenyl)-2-[6-methyl-2,4-dioxo-3-(thietan-3-yl)-1,2,3,4-tetrahydropyrimidine-1-yl]acetamide. RUSS J GEN CHEM+ 2021. [DOI: 10.1134/s1070363221050054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
10
|
Xu P, Sattasathuchana T, Guidez E, Webb SP, Montgomery K, Yasini H, Pedreira IFM, Gordon MS. Computation of host-guest binding free energies with a new quantum mechanics based mining minima algorithm. J Chem Phys 2021; 154:104122. [PMID: 33722015 DOI: 10.1063/5.0040759] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new method called QM-VM2 is presented that efficiently combines statistical mechanics with quantum mechanical (QM) energy potentials in order to calculate noncovalent binding free energies of host-guest systems. QM-VM2 efficiently couples the use of semi-empirical QM (SEQM) energies and geometry optimizations with an underlying molecular mechanics (MM) based conformational search, to find low SEQM energy minima, and allows for processing of these minima at higher levels of ab initio QM theory. A progressive geometry optimization scheme is introduced as a means to increase conformational sampling efficiency. The newly implemented QM-VM2 is used to compute the binding free energies of the host molecule cucurbit[7]uril and a set of 15 guest molecules. The results are presented along with comparisons to experimentally determined binding affinities. For the full set of 15 host-guest complexes, which have a range of formal charges from +1 to +3, SEQM-VM2 based binding free energies show poor correlation with experiment, whereas for the ten +1 complexes only, a significant correlation (R2 = 0.8) is achieved. SEQM-VM2 generation of conformers followed by single-point ab initio QM calculations at the dispersion corrected restricted Hartree-Fock-D3(BJ) and TPSS-D3(BJ) levels of theory, as post-processing corrections, yields a reasonable correlation with experiment for the full set of host-guest complexes (R2 = 0.6 and R2 = 0.7, respectively) and an excellent correlation for the +1 formal charge set (R2 = 1.0 and R2 = 0.9, respectively), as long as a sufficiently large basis set (triple-zeta quality) is employed. The importance of the inclusion of configurational entropy, even at the MM level, for the achievement of good correlation with experiment was demonstrated by comparing the calculated ΔE values with experiment and finding a considerably poorer correlation with experiment than for the calculated free energy ΔE - TΔS. For the complete set of host-guest systems with the range of formal charges, it was observed that the deviation of the predicted binding free energy from experiment correlates somewhat with the net charge of the systems. This observation leads to a simple empirical interpolation scheme to improve the linear regression of the full set.
Collapse
Affiliation(s)
- Peng Xu
- Department of Chemistry, Iowa State University, Ames, Iowa 50014, USA
| | | | - Emilie Guidez
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80204, USA
| | - Simon P Webb
- VeraChem LLC, 12850 Middlebrook Rd. Ste 205, Germantown, Maryland 20874-5244, USA
| | | | - Hussna Yasini
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80204, USA
| | - Iara F M Pedreira
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80204, USA
| | - Mark S Gordon
- Department of Chemistry, Iowa State University, Ames, Iowa 50014, USA
| |
Collapse
|
11
|
Błasiak B, Bednarska JD, Chołuj M, Góra RW, Bartkowiak W. Ab initio effective one-electron potential operators: Applications for charge-transfer energy in effective fragment potentials. J Comput Chem 2021; 42:398-411. [PMID: 33349929 DOI: 10.1002/jcc.26462] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2020] [Accepted: 11/18/2020] [Indexed: 12/19/2022]
Abstract
The concept of effective one-electron potentials (EOPs) has proven to be extremely useful in efficient description of electronic structure of chemical systems, especially extended molecular aggregates such as interacting molecules in condensed phases. Here, a general method for EOP-based elimination of electron repulsion integrals is presented, that is tuned toward the fragment-based calculation methodologies such as the second generation of the effective fragment potentials (EFP2) method. Two general types of the EOP operator matrix elements are distinguished and treated either via the distributed multipole expansion or the extended density fitting (DF) schemes developed in this work. The EOP technique is then applied to reduce the high computational costs of the effective fragment charge-transfer (CT) terms being the bottleneck of EFP2 potentials. The alternative EOP-based CT energy model is proposed, derived within the framework of intermolecular perturbation theory with Hartree-Fock noninteracting reference wavefunctions, compatible with the original EFP2 formulation. It is found that the computational cost of the EFP2 total interaction energy calculation can be reduced by up to 38 times when using the EOP-based formulation of CT energy, as compared to the original EFP2 scheme, without compromising the accuracy for a wide range of weakly interacting neutral and ionic molecular fragments. The proposed model can thus be used routinely within the EFP2 framework.
Collapse
Affiliation(s)
- Bartosz Błasiak
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland
| | - Joanna D Bednarska
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland
| | - Marta Chołuj
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland
| | - Robert W Góra
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland
| | - Wojciech Bartkowiak
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland
| |
Collapse
|
12
|
Valiakhmetova OY, Kuznetsov VV. Conformational Analysis of
2-Isopropyl-5-methoxy-5-methyl-1,3,2-dioxaborinane in Chloroform Solution: Effect of
“Magic” Solvent Molecule. RUSSIAN JOURNAL OF ORGANIC CHEMISTRY 2021. [DOI: 10.1134/s1070428021010036] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
|
13
|
Ghiasi R, Zandiyeh Z. Theoretical study of the influence of solvent polarity on the 31P and 13C NMR parameters of the Ru(PH3)4(η2-benzyne) complex. INORG CHEM COMMUN 2021. [DOI: 10.1016/j.inoche.2020.108412] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
14
|
Kuznetsov V. Stereochemistry of Simple Molecules inside Nanotubes and Fullerenes: Unusual Behavior of Usual Systems. Molecules 2020; 25:molecules25102437. [PMID: 32456128 PMCID: PMC7287839 DOI: 10.3390/molecules25102437] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2020] [Revised: 05/19/2020] [Accepted: 05/21/2020] [Indexed: 12/13/2022] Open
Abstract
Over the past three decades, carbon nanotubes and fullerenes have become remarkable objects for starting the implementation of new models and technologies in different branches of science. To a great extent, this is defined by the unique electronic and spatial properties of nanocavities due to the ramified π-electron systems. This provides an opportunity for the formation of endohedral complexes containing non-covalently bonded atoms or molecules inside fullerenes and nanotubes. The guest species are exposed to the force field of the nanocavity, which can be described as a combination of electronic and steric requirements. Its action significantly changes conformational properties of even relatively simple molecules, including ethane and its analogs, as well as compounds with C-O, C-S, B-B, B-O, B-N, N-N, Al-Al, Si-Si and Ge-Ge bonds. Besides that, the cavity of the host molecule dramatically alters the stereochemical characteristics of cyclic and heterocyclic systems, affects the energy of pyramidal nitrogen inversion in amines, changes the relative stability of cis and trans isomers and, in the case of chiral nanotubes, strongly influences the properties of R- and S- enantiomers. The present review aims at primary compilation of such unusual stereochemical effects and initial evaluation of the nature of the force field inside nanotubes and fullerenes.
Collapse
Affiliation(s)
- Valerij Kuznetsov
- Ufa State Aviation Technical University, K. Marksa, 12, Ufa 450008, Russia;
- Ufa State Petroleum Technological University, Kosmonavtov, 1, Ufa 450062, Russia
| |
Collapse
|
15
|
Barca GMJ, Bertoni C, Carrington L, Datta D, De Silva N, Deustua JE, Fedorov DG, Gour JR, Gunina AO, Guidez E, Harville T, Irle S, Ivanic J, Kowalski K, Leang SS, Li H, Li W, Lutz JJ, Magoulas I, Mato J, Mironov V, Nakata H, Pham BQ, Piecuch P, Poole D, Pruitt SR, Rendell AP, Roskop LB, Ruedenberg K, Sattasathuchana T, Schmidt MW, Shen J, Slipchenko L, Sosonkina M, Sundriyal V, Tiwari A, Galvez Vallejo JL, Westheimer B, Włoch M, Xu P, Zahariev F, Gordon MS. Recent developments in the general atomic and molecular electronic structure system. J Chem Phys 2020; 152:154102. [PMID: 32321259 DOI: 10.1063/5.0005188] [Citation(s) in RCA: 494] [Impact Index Per Article: 123.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods such as the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree-Fock, and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as it is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite such as GAMESS is maintained and developed is considered. Future developments are briefly summarized.
Collapse
Affiliation(s)
- Giuseppe M J Barca
- Research School of Computer Science, Australian National University, Canberra, ACT 2601, Australia
| | - Colleen Bertoni
- Argonne Leadership Computing Facility, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Laura Carrington
- EP Analytics, 12121 Scripps Summit Dr. Ste. 130, San Diego, California 92131, USA
| | - Dipayan Datta
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Nuwan De Silva
- Department of Physical and Biological Sciences, Western New England University, Springfield, Massachusetts 01119, USA
| | - J Emiliano Deustua
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Umezono 1-1-1, Tsukuba 305-8568, Japan
| | - Jeffrey R Gour
- Microsoft, 15590 NE 31st St., Redmond, Washington 98052, USA
| | - Anastasia O Gunina
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Emilie Guidez
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217, USA
| | - Taylor Harville
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Stephan Irle
- Computational Science and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830, USA
| | - Joe Ivanic
- Advanced Biomedical Computational Science, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21702, USA
| | - Karol Kowalski
- Physical Sciences Division, Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box 999, Richland, Washington 99352, USA
| | - Sarom S Leang
- EP Analytics, 12121 Scripps Summit Dr. Ste. 130, San Diego, California 92131, USA
| | - Hui Li
- Department of Chemistry, University of Nebraska, Lincoln, Nebraska 68588, USA
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, People's Republic of China
| | - Jesse J Lutz
- Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Ilias Magoulas
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Joani Mato
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Vladimir Mironov
- Department of Chemistry, Lomonosov Moscow State University, Leninskie Gory 1/3, Moscow 119991, Russian Federation
| | - Hiroya Nakata
- Kyocera Corporation, Research Institute for Advanced Materials and Devices, 3-5-3 Hikaridai Seika-cho, Souraku-gun, Kyoto 619-0237, Japan
| | - Buu Q Pham
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - David Poole
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Spencer R Pruitt
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Alistair P Rendell
- Research School of Computer Science, Australian National University, Canberra, ACT 2601, Australia
| | - Luke B Roskop
- Cray Inc., a Hewlett Packard Enterprise Company, 2131 Lindau Ln #1000, Bloomington, Minnesota 55425, USA
| | - Klaus Ruedenberg
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | | | - Michael W Schmidt
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Jun Shen
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Lyudmila Slipchenko
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Masha Sosonkina
- Department of Computational Modeling and Simulation Engineering, Old Dominion University, Norfolk, Virginia 23529, USA
| | - Vaibhav Sundriyal
- Department of Computational Modeling and Simulation Engineering, Old Dominion University, Norfolk, Virginia 23529, USA
| | - Ananta Tiwari
- EP Analytics, 12121 Scripps Summit Dr. Ste. 130, San Diego, California 92131, USA
| | - Jorge L Galvez Vallejo
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Bryce Westheimer
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Marta Włoch
- 530 Charlesina Dr., Rochester, Michigan 48306, USA
| | - Peng Xu
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Federico Zahariev
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Mark S Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| |
Collapse
|
16
|
Tran AL, Guidez EB. Quantum Mechanical Modeling of the Interactions between Noble Metal (Ag and Au) Nanoclusters and Water with the Effective Fragment Potential Method. ACS OMEGA 2020; 5:7446-7455. [PMID: 32280887 PMCID: PMC7144145 DOI: 10.1021/acsomega.0c00132] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/10/2020] [Accepted: 03/13/2020] [Indexed: 06/11/2023]
Abstract
Explicit solvent interactions can significantly alter the physical and chemical properties of noble metal (e.g., gold and silver) nanoclusters. In order to compute these solvent interactions at a reasonable computational cost, a quantum mechanical (QM)/molecular mechanics (MM) approach, where the metal nanocluster is treated with full QM and the water molecules are treated with a MM force field, can be used. However, classical MM force fields were typically parameterized using molecules containing main group elements as the reference. The accuracy of noble metal-solvent interactions obtained with these force fields therefore remains unpredictable. The effective fragment potential (EFP) force field, designed to model explicitly solvated systems, represents an attractive method to simulate solvated noble metal nanoclusters because it is derived from first principles and contains few or no fitted parameters, depending on implementation. At the density functional theory-optimized geometries, good correlation is obtained between the nanocluster-water interaction energies computed with EFP and those computed with the reference coupled cluster singles, doubles, and perturbative triples method. It is shown that the EFP method gives qualitatively accurate interaction energies at medium-large intermolecular distances for various molecular configurations. In order to achieve higher quantitative accuracy, the first solvation shell should be treated with full QM, if possible. EFP is therefore a promising method for the QM modeling of explicitly solvated silver and gold nanoclusters.
Collapse
|