1
|
Mills JD, Ben-Nun M, Rollin K, Bromley MWJ, Li J, Hinde RJ, Winstead CL, Sheehy JA, Boatz JA, Langhoff PW. Atomic Spectral Methods for Ab Initio Molecular Electronic Energy Surfaces: Transitioning From Small-Molecule to Biomolecular-Suitable Approaches. J Phys Chem B 2016; 120:8321-37. [PMID: 27232159 DOI: 10.1021/acs.jpcb.6b02021] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Continuing attention has addressed incorportation of the electronically dynamical attributes of biomolecules in the largely static first-generation molecular-mechanical force fields commonly employed in molecular-dynamics simulations. We describe here a universal quantum-mechanical approach to calculations of the electronic energy surfaces of both small molecules and large aggregates on a common basis which can include such electronic attributes, and which also seems well-suited to adaptation in ab initio molecular-dynamics applications. In contrast to the more familiar orbital-product-based methodologies employed in traditional small-molecule computational quantum chemistry, the present approach is based on an "ex-post-facto" method in which Hamiltonian matrices are evaluated prior to wave function antisymmetrization, implemented here in the support of a Hilbert space of orthonormal products of many-electron atomic spectral eigenstates familiar from the van der Waals theory of long-range interactions. The general theory in its various forms incorporates the early semiempirical atoms- and diatomics-in-molecules approaches of Moffitt, Ellison, Tully, Kuntz, and others in a comprehensive mathematical setting, and generalizes the developments of Eisenschitz, London, Claverie, and others addressing electron permutation symmetry adaptation issues, completing these early attempts to treat van der Waals and chemical forces on a common basis. Exact expressions are obtained for molecular Hamiltonian matrices and for associated energy eigenvalues as sums of separate atomic and interaction-energy terms, similar in this respect to the forms of classical force fields. The latter representation is seen to also provide a long-missing general definition of the energies of individual atoms and of their interactions within molecules and matter free from subjective additional constraints. A computer code suite is described for calculations of the many-electron atomic eigenspectra and the pairwise-atomic Hamiltonian matrices required for practical applications. These matrices can be retained as functions of scalar atomic-pair separations and employed in assembling aggregate Hamiltonian matrices, with Wigner rotation matrices providing analytical representations of their angular degrees of freedom. In this way, ab initio potential energy surfaces are obtained in the complete absence of repeated evaluations and transformations of the one- and two-electron integrals at different molecular geometries required in most ab inito molecular calculations, with large Hamiltonian matrix assembly simplified and explicit diagonalizations avoided employing partitioning and Brillouin-Wigner or Rayleigh-Schrödinger perturbation theory. Illustrative applications of the important components of the formalism, selected aspects of the scaling of the approach, and aspects of "on-the-fly" interfaces with Monte Carlo and molecular-dynamics methods are described in anticipation of subsequent applications to biomolecules and other large aggregates.
Collapse
Affiliation(s)
- Jeffrey D Mills
- Air Force Research Laboratory , 10 East Saturn Boulevard, Edwards AFB, California 93524-7680, United States
| | - Michal Ben-Nun
- Predictive Science, Inc. , 9990 Mesa Rim Road #170, San Diego, California 92121, United States
| | - Kyle Rollin
- Northrup Grumman Corporation , 1 Rancho Carmel Drive, San Diego, California 92128, United States
| | - Michael W J Bromley
- School of Mathematics and Physics, The University of Queensland , Brisbane, Queensland 4072, Australia
| | - Jiabo Li
- Accelrys Inc. , 10188 Telesis Court #100, San Diego, California 92121-4779, United States
| | - Robert J Hinde
- Department of Chemistry, University of Tennessee , Knoxville, Tennessee 37996-1600, United States
| | - Carl L Winstead
- A.A. Noyes Laboratory of Chemical Physics, California Institute of Technology , Pasadena, California 91125, United States
| | - Jeffrey A Sheehy
- NASA Headquarters , 300 E Street SW, Suite 5R30, Washington, DC 201546, United States
| | - Jerry A Boatz
- Air Force Research Laboratory , 10 East Saturn Boulevard, Edwards AFB, California 93524-7680, United States
| | - Peter W Langhoff
- Department of Chemistry and Biochemistry, University of California, San Diego , 9500 Gilman Drive, MS 0365, La Jolla, California 92093-0365, United States
| |
Collapse
|
2
|
NAUMKIN FY, WALES DJ. Structure and properties of Nen+;clusters from a diatomics-in-molecules approach. Mol Phys 2010. [DOI: 10.1080/002689798168970] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
|
3
|
|
4
|
Albu TV, Espinosa-García J, Truhlar DG. Computational Chemistry of Polyatomic Reaction Kinetics and Dynamics: The Quest for an Accurate CH5 Potential Energy Surface. Chem Rev 2007; 107:5101-32. [DOI: 10.1021/cr078026x] [Citation(s) in RCA: 58] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Titus V. Albu
- Department of Chemistry, Box 5055, Tennessee Technological University, Cookeville, Tennessee 38505
| | | | - Donald G. Truhlar
- Department of Chemistry and Supercomputer Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
5
|
The Calculation of Potential Energy Surfaces for Excited States. ADVANCES IN CHEMICAL PHYSICS 2007. [DOI: 10.1002/9780470142745.ch8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register]
|
6
|
Cooper DL, Gerratt J, Raimondi M. Modern Valence Bond Theory. ADVANCES IN CHEMICAL PHYSICS 2007. [DOI: 10.1002/9780470142943.ch6] [Citation(s) in RCA: 89] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
|
7
|
Model Potential Energy Surfaces for Inelastic and Charge-Transfer Processes in Ion-Molecule Collisions. ADVANCES IN CHEMICAL PHYSICS 2007. [DOI: 10.1002/9780470141403.ch3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/28/2023]
|
8
|
Tishchenko O, Truhlar DG. Optimizing the Performance of the Multiconfiguration Molecular Mechanics Method. J Phys Chem A 2006; 110:13530-6. [PMID: 17165880 DOI: 10.1021/jp0640833] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
Multiconfiguration molecular mechanics (MCMM) is a general algorithm for constructing potential energy surfaces for reactive systems (Kim, Y.; Corchado, J. C.; Villà, J.; Xing, J.; Truhlar, D. G. J. Chem. Phys. 2000, 112, 2718). This paper illustrates how the performance of the MCMM method can be improved by adopting improved molecular mechanics parameters. We carry out calculations of reaction rate constants using variational transition state theory with optimized multidimensional tunneling on the MCMM PESs for three hydrogen transfer reactions, and we compare the results to direct dynamics. We find that the MCMM method with as little as one electronic structure Hessian can describe the dynamically important regions of the ground-electronic state PES, including the corner-cutting-tunneling region of the reaction swath, with practical accuracy.
Collapse
Affiliation(s)
- Oksana Tishchenko
- Chemistry Department and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | | |
Collapse
|
9
|
Chakraborty A, Zhao Y, Lin H, Truhlar DG. Combined valence bond-molecular mechanics potential-energy surface and direct dynamics study of rate constants and kinetic isotope effects for the H+C2H6 reaction. J Chem Phys 2006; 124:044315. [PMID: 16460170 DOI: 10.1063/1.2132276] [Citation(s) in RCA: 42] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This article presents a multifaceted study of the reaction H+C(2)H(6)-->H(2)+C(2)H(5) and three of its deuterium-substituted isotopologs. First we present high-level electronic structure calculations by the W1, G3SX, MCG3-MPWB, CBS-APNO, and MC-QCISD/3 methods that lead to a best estimate of the barrier height of 11.8+/-0.5 kcal/mol. Then we obtain a specific reaction parameter for the MPW density functional in order that it reproduces the best estimate of the barrier height; this yields the MPW54 functional. The MPW54 functional, as well as the MPW60 functional that was previously parametrized for the H+CH(4) reaction, is used with canonical variational theory with small-curvature tunneling to calculate the rate constants for all four ethane reactions from 200 to 2000 K. The final MPW54 calculations are based on curvilinear-coordinate generalized-normal-mode analysis along the reaction path, and they include scaled frequencies and an anharmonic C-C bond torsion. They agree with experiment within 31% for 467-826 K except for a 38% deviation at 748 K; the results for the isotopologs are predictions since these rate constants have never been measured. The kinetic isotope effects (KIEs) are analyzed to reveal the contributions from subsets of vibrational partition functions and from tunneling, which conspire to yield a nonmonotonic temperature dependence for one of the KIEs. The stationary points and reaction-path potential of the MPW54 potential-energy surface are then used to parametrize a new kind of analytical potential-energy surface that combines a semiempirical valence bond formalism for the reactive part of the molecule with a standard molecular mechanics force field for the rest; this may be considered to be either an extension of molecular mechanics to treat a reactive potential-energy surface or a new kind of combined quantum-mechanical/molecular mechanical (QM/MM) method in which the QM part is semiempirical valence bond theory; that is, the new potential-energy surface is a combined valence bond molecular mechanics (CVBMM) surface. Rate constants calculated with the CVBMM surface agree with the MPW54 rate constants within 12% for 534-2000 K and within 23% for 200-491 K. The full CVBMM potential-energy surface is now available for use in variety of dynamics calculations, and it provides a prototype for developing CVBMM potential-energy surfaces for other reactions.
Collapse
Affiliation(s)
- Arindam Chakraborty
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | | | | | | |
Collapse
|
10
|
Lin H, Pu J, Albu TV, Truhlar DG. Efficient Molecular Mechanics for Chemical Reactions: Multiconfiguration Molecular Mechanics Using Partial Electronic Structure Hessians. J Phys Chem A 2004. [DOI: 10.1021/jp049972+] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
11
|
Varandas AJC. A realistic multi-sheeted potential energy surface for NO2(2A′) from the double many-body expansion method and a novel multiple energy-switching scheme. J Chem Phys 2003. [DOI: 10.1063/1.1586911] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
12
|
Affiliation(s)
- Donald G. Truhlar
- Department of Chemistry and Supercomputer Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, Minnesota 55455−043l
| |
Collapse
|
13
|
Albu TV, Corchado JC, Truhlar DG. Molecular Mechanics for Chemical Reactions: A Standard Strategy for Using Multiconfiguration Molecular Mechanics for Variational Transition State Theory with Optimized Multidimensional Tunneling. J Phys Chem A 2001. [DOI: 10.1021/jp011951h] [Citation(s) in RCA: 83] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Titus V. Albu
- Department of Chemistry and Supercomputer Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - José C. Corchado
- Department of Chemistry and Supercomputer Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputer Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
14
|
Varandas AJC. Multivalued Potential Energy Surfaces for Dynamics Studies. LECTURE NOTES IN CHEMISTRY 2000. [DOI: 10.1007/978-3-642-57051-3_3] [Citation(s) in RCA: 21] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
15
|
Naumkin F. DIM models for RgX2− systems: suppressed influence of spin-orbit coupling and induced multipole effects for the Ar–I2− interaction. Chem Phys 1999. [DOI: 10.1016/s0301-0104(98)00356-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
16
|
Naumkin FY, McCourt FRW. Contributions of the two conformers to the microwave spectrum and scattering cross-section of the He–Cl2 van der Waals system, evaluated from an ab initio potential energy surface. J Chem Phys 1998. [DOI: 10.1063/1.476419] [Citation(s) in RCA: 23] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
|
17
|
Single- versus multi-state DIM model for RgX2 systems: On the influence of spin–orbit coupling on Ar–I2 potentials. Chem Phys 1998. [DOI: 10.1016/s0301-0104(97)00298-x] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
18
|
Vijay A, Billing GD. Semiclassical wave packet calculations on ion–molecule reactions: Studies on B+(3Pu)+H2 reaction. J Chem Phys 1997. [DOI: 10.1063/1.474655] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
19
|
|
20
|
Lepetit B, Launay J, Le Dourneuf M. Quantum study of electronically non-adiabatic collinear reactions. II. Influence of spin-orbit transitions on the F + HH reaction. Chem Phys 1986. [DOI: 10.1016/0301-0104(86)87043-4] [Citation(s) in RCA: 24] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
21
|
Perturbation non-hermitean DIM scheme for calculations of the width of autoionizing states of collision systems. Chem Phys 1986. [DOI: 10.1016/0301-0104(86)80166-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
22
|
|
23
|
Kuntz P, Polák R. A sensitivity analysis of a diatomics-in-molecules model of the 1A′ states of H2O. Chem Phys 1985. [DOI: 10.1016/0301-0104(85)80181-6] [Citation(s) in RCA: 19] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
24
|
Koshi M, Itoh H, Matsui H. Vibrational energy transfer from highly excited anharmonic oscillators: Quasiclassical Monte Carlo trajectory study of Br2–Ar and Br2–Br system. J Chem Phys 1985. [DOI: 10.1063/1.448662] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
|
26
|
Wu AJA. Generalized diatomics-in-molecules theory. Mol Phys 1984. [DOI: 10.1080/00268978400101021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
27
|
Last I, Baer M. Three‐dimensional DIM‐3C potential energy surfaces for the reactions H+XY and X+HY (X,Y=F, Cl, Br, I). J Chem Phys 1984. [DOI: 10.1063/1.447077] [Citation(s) in RCA: 35] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Ljudkovski S, Nemukhin A, Stepanov N. Diatomics-in-molecules studies of the Li + H2O interaction. ACTA ACUST UNITED AC 1983. [DOI: 10.1016/0166-1280(83)80188-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
29
|
Vojtík J, Paidarová I. Hermitean DIM model for potential energy surfaces relevant to chemiionization processes. Chem Phys Lett 1983. [DOI: 10.1016/0009-2614(83)80536-3] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
30
|
|
31
|
Ellison FO. Scaled diatomics‐in‐molecules (SDIM) theory. I. General theory and preliminary applications to H3and H4. J Chem Phys 1983. [DOI: 10.1063/1.445370] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
32
|
Schneider F, Zülicke L, Polák R, Vojtík J. Diatomics-in-molecules potential energy surfaces of (HeH2)+: An extended basis set calculation. Chem Phys 1983. [DOI: 10.1016/0301-0104(83)85040-x] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
33
|
Group-theoretical approach in the method of diatomics in molecules and the secular equation for clusters of alkali metals. THEOR EXP CHEM+ 1983. [DOI: 10.1007/bf00522414] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
34
|
Last I. Semi-empirical potential energy surfaces for the linear HXY (X, Y = F, Cl, Br, I) systems. Chem Phys 1982. [DOI: 10.1016/0301-0104(82)88146-9] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
35
|
Whitlock PA, Muckerman JT, Fisher ER. Quasiclassical trajectory investigation of the reaction O(1D)+H2. J Chem Phys 1982. [DOI: 10.1063/1.443572] [Citation(s) in RCA: 83] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
36
|
Kuntz PJ, Schreiber JL. A systematic procedure for extracting fragment matrices for the method of diatomics‐in‐molecules fromabinitiocalculations on diatomics. J Chem Phys 1982. [DOI: 10.1063/1.443488] [Citation(s) in RCA: 31] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
37
|
Garrett BC, Truhlar DG, Magnuson AW. New semiempirical method of modeling potential energy surfaces for generalized TST and application to the kinetic isotope effects in the Cl–H–H system. J Chem Phys 1982. [DOI: 10.1063/1.443304] [Citation(s) in RCA: 42] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
38
|
Last (Lyast) I, Baer M. Semiempirical three‐dimensional potential energy surfaces suitable for both reaction channels of the XH2 system (X = F, Cl). J Chem Phys 1981. [DOI: 10.1063/1.441835] [Citation(s) in RCA: 42] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
39
|
|
40
|
Isaacson AD, Muckerman JT. Valence bond diatomics‐in‐molecules (DIM) treatment of collinear interactions of Group IIaand IIbmetal atoms with hydrogen halides: Application to the CaHCl system. J Chem Phys 1980. [DOI: 10.1063/1.440309] [Citation(s) in RCA: 22] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
Roach AC, Gimzewski E. A semiempirical valence bond π‐electron theory related to the method of diatomics‐in‐molecules. J Chem Phys 1980. [DOI: 10.1063/1.440241] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|