skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Acceleration of Semiempirical QM/MM Methods through Message Passage Interface (MPI), Hybrid MPI/Open Multiprocessing, and Self-Consistent Field Accelerator Implementations

Journal Article · · Journal of Chemical Theory and Computation

The strategy and implementation of scalable and efficient semiempirical (SE) QM/MM methods in CHARMM are described. The serial version of the code was first profiled to identify routines that required parallelization. Afterward, the code was parallelized and accelerated with three approaches. The first approach was the parallelization of the entire QM/MM routines, including the Fock matrix diagonalization routines, using the CHARMM message passage interface (MPI) machinery. In the second approach, two different self-consistent field (SCF) energy convergence accelerators were implemented using density and Fock matrices as targets for their extrapolations in the SCF procedure. In the third approach, the entire QM/MM and MM energy routines were accelerated by implementing the hybrid MPI/open multiprocessing (OpenMP) model in which both the task- and loop-level parallelization strategies were adopted to balance loads between different OpenMP threads. The present implementation was tested on two solvated enzyme systems (including <100 QM atoms) and an SN2 symmetric reaction in water. The MPI version exceeded existing SE QM methods in CHARMM, which include the SCC-DFTB and SQUANTUM methods, by at least 4-fold. The use of SCF convergence accelerators further accelerated the code by ~12–35% depending on the size of the QM region and the number of CPU cores used. Although the MPI version displayed good scalability, the performance was diminished for large numbers of MPI processes due to the overhead associated with MPI communications between nodes. This issue was partially overcome by the hybrid MPI/OpenMP approach which displayed a better scalability for a larger number of CPU cores (up to 64 CPUs in the tested systems).

Research Organization:
Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC)
Sponsoring Organization:
USDOE
OSTI ID:
1480127
Journal Information:
Journal of Chemical Theory and Computation, Vol. 13, Issue 8; ISSN 1549-9618
Publisher:
American Chemical SocietyCopyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 11 works
Citation information provided by
Web of Science

References (95)

Theoretical studies of enzymic reactions: Dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme journal May 1976
A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations journal July 1990
Perspective: Quantum mechanical methods in biochemistry and biophysics journal October 2016
New faster CHARMM molecular dynamics engine journal December 2013
The energy landscape of adenylate kinase during catalysis journal January 2015
Role of Protein Dynamics in Allosteric Control of the Catalytic Phosphoryl Transfer of Insulin Receptor Kinase journal September 2015
Die Berechnung optischer und elektrostatischer Gitterpotentiale journal January 1921
Particle mesh Ewald: An N ⋅log( N ) method for Ewald sums in large systems journal June 1993
Zonal methods for the parallel execution of range-limited N-body simulations journal January 2007
Classical Electrostatics for Biomolecular Simulations journal August 2013
Preface journal September 1994
OpenMP: an industry standard API for shared-memory programming journal January 1998
GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers journal September 2015
Fast Parallel Algorithms for Short-Range Molecular Dynamics journal March 1995
CHARMM: The biomolecular simulation program journal July 2009
Scalable molecular dynamics with NAMD journal January 2005
A smooth particle mesh Ewald method journal November 1995
Ewald mesh method for quantum mechanical calculations journal March 2012
Efficient particle-mesh Ewald based approach to fixed and induced dipolar interactions journal December 2000
An efficient algorithm for multipole energies and derivatives based on spherical harmonics and extensions to particle mesh Ewald journal May 2014
Multipolar Ewald Methods, 1: Theory, Accuracy, and Performance journal January 2015
Transition state stabilization and substrate strain in enzyme catalysis: ab initio QM/MM modelling of the chorismate mutase reaction journal January 2004
Direct dynamics calculations with NDDO (neglect of diatomic differential overlap) molecular orbital theory with specific reaction parameters journal June 1991
Exploring SCC-DFTB Paths for Mapping QM/MM Reaction Mechanisms journal July 2007
Reaction Path Force Matching: A New Strategy of Fitting Specific Reaction Parameters for Semiempirical Methods in Combined QM/MM Simulations journal June 2014
Parameterization of NDDO wavefunctions using genetic algorithms. An evolutionary approach to parameterizing potential energy surfaces and direct dynamics calculations for organic reactions journal February 1995
Specific Reaction Parametrization of the AM1/d Hamiltonian for Phosphoryl Transfer Reactions:  H, O, and P Atoms journal January 2007
Description of Phosphate Hydrolysis Reactions with the Self-Consistent-Charge Density-Functional-Tight-Binding (SCC-DFTB) Theory. 1. Parameterization journal November 2008
AM1 Specific Reaction Parameters for Reactions of Hydroxide Ion with Halomethanes in Complex Environments: Development and Testing journal September 2013
Quantum mechanical study of solvent effects in a prototype S N 2 reaction in solution: Cl attack on CH 3 Cl journal February 2014
AM1/d-CB1: A Semiempirical Model for QM/MM Simulations of Chemical Glycobiology Systems journal September 2014
Parameterization of DFTB3/3OB for Sulfur and Phosphorus for Chemical and Biological Applications journal March 2014
Specific parametrisation of a hybrid potential to simulate reactions in phosphatases journal January 2006
A PM3/d specific reaction parameterization for iron atom in the hydrogen abstraction catalyzed by soybean lipoxygenase-1 journal January 2007
Parameterization of semiempirical methods to treat nucleophilic attacks to biological phosphates: AM1/d parameters for phosphorus journal April 2003
Computational chemistry on Fujitsu vector–parallel processors: Development and performance of applications software journal July 2000
GPU Linear Algebra Libraries and GPGPU Programming for Accelerating MOPAC Semiempirical Quantum Chemistry Calculations
  • Maia, Julio Daniel Carvalho; Urquiza Carvalho, Gabriel Aires; Mangueira, Carlos Peixoto
  • Journal of Chemical Theory and Computation, Vol. 8, Issue 9 https://doi.org/10.1021/ct3004645
journal August 2012
Semiempirical Quantum Chemical Calculations Accelerated on a Hybrid Multicore CPU–GPU Computing Platform journal June 2012
Unified Approach for Molecular Dynamics and Density-Functional Theory journal November 1985
Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals journal June 2001
Curvy-steps approach to constraint-free extended-Lagrangian ab initio molecular dynamics, using atom-centered basis functions: Convergence toward Born–Oppenheimer trajectories journal December 2004
Acceleration of Semiempirical Quantum Mechanical Calculations by Extended Lagrangian Molecular Dynamics Approach journal July 2013
An Efficient Linear-Scaling Ewald Method for Long-Range Electrostatic Interactions in Combined QM/MM Calculations journal October 2004
Extension of MNDO to d Orbitals:  Parameters and Results for the Second-Row Elements and for the Zinc Group journal January 1996
Acceleration of Ab Initio QM/MM Calculations under Periodic Boundary Conditions by Multiscale and Multiple Time Step Approaches journal September 2014
Convergence acceleration of iterative sequences. the case of scf iteration journal July 1980
ImprovedSCF convergence acceleration journal January 1982
Fast semiempirical calculations: Fast Semiempirical Calculations journal June 1982
Extended Lagrangian Born–Oppenheimer molecular dynamics with dissipation journal June 2009
Lagrangian formulation with dissipation of Born-Oppenheimer molecular dynamics using the density-functional tight-binding method journal July 2011
Fock matrix dynamics journal March 2004
Accelerated, energy-conserving Born–Oppenheimer molecular dynamics via Fock matrix extrapolation journal January 2005
New Developments in Molecular Orbital Theory journal April 1951
Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Theory, Implementation, and Parameters journal January 2016
Electronic Population Analysis on LCAO–MO Molecular Wave Functions. I journal October 1955
CHARMM: A program for macromolecular energy, minimization, and dynamics calculations journal July 1983
The implementation of a fast and accurate QM/MM potential method in Amber journal January 2008
COSMO: a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient journal January 1993
Smooth Solvation Method for d-Orbital Semiempirical Calculations of Biological Reactions. 1. Implementation journal May 2005
Beyond the MNDO model: Methodical considerations and numerical results journal July 1993
Orthogonalization corrections for semiempirical methods journal April 2000
Ground states of molecules. 38. The MNDO method. Approximations and parameters journal June 1977
Ground states of molecules. 39. MNDO results for molecules containing hydrogen, carbon, nitrogen, and oxygen journal June 1977
Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model journal June 1985
Optimization of parameters for semiempirical methods I. Method journal March 1989
Optimization of parameters for semiempirical methods II. Applications journal March 1989
Optimization of parameters for semiempirical methods. III Extension of PM3 to Be, Mg, Zn, Ga, Ge, As, Se, Cd, In, Sn, Sb, Te, Hg, Tl, Pb, and Bi journal April 1991
A combinedab initio quantum mechanical and molecular mechanical method for carrying out simulations on complex molecular systems: Applications to the CH3Cl + Cl? exchange reaction and gas phase protonation of polyethers journal December 1986
A Generalized Hybrid Orbital (GHO) Method for the Treatment of Boundary Atoms in Combined QM/MM Calculations journal June 1998
Comparison of simple potential functions for simulating liquid water journal July 1983
All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins journal April 1998
All-atom empirical force field for nucleic acids: I. Parameter optimization based on small molecule and condensed phase macromolecular target data journal January 2000
Development and current status of the CHARMM force field for nucleic acids journal January 2000
Extending the treatment of backbone energetics in protein force fields: Limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations journal August 2004
Crystal structure of the activated insulin receptor tyrosine kinase in complex with peptide substrate and ATP analog journal September 1997
Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes journal March 1977
Parametrization and Benchmark of DFTB3 for Organic Molecules journal November 2012
Parametrization of DFTB3/3OB for Magnesium and Zinc for Chemical and Biological Applications journal August 2014
Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling journal February 1977
THE weighted histogram analysis method for free-energy calculations on biomolecules. I. The method journal October 1992
Implementation of an adaptive umbrella sampling method for the calculation of multidimensional potential of mean force of chemical reactions in solution journal September 2003
The Scalasca performance toolset architecture journal January 2010
Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties journal September 1998
A QM/MM Implementation of the Self-Consistent Charge Density Functional Tight Binding (SCC-DFTB) Method journal January 2001
A Perspective on Enzyme Catalysis journal August 2003
Connecting Protein Conformational Dynamics with Catalytic Function As Illustrated in Dihydrofolate Reductase journal January 2013
Enzymatic transition states and dynamic motion in barrier crossing journal July 2009
At the dawn of the 21st century: Is dynamics the missing link for understanding enzyme catalysis? journal January 2010
Perspectives on Electrostatics and Conformational Motions in Enzyme Catalysis journal January 2015
Role of Dynamics in Enzyme Catalysis: Substantial versus Semantic Controversies journal December 2014
Are there dynamical effects in enzyme catalysis? Some thoughts concerning the enzymatic chemical step journal September 2015
Revisiting the finite temperature string method for the calculation of reaction tubes and free energies journal May 2009
Free energy of conformational transition paths in biomolecules: The string method and its application to myosin VI journal February 2011
A Simple and Accurate Method To Calculate Free Energy Profiles and Reaction Rates from Restrained Molecular Simulations of Diffusive Processes journal May 2016
Escaping free-energy minima journal September 2002

Cited By (1)

Force matching as a stepping stone to QM/MM CB[8] host/guest binding free energies: a SAMPL6 cautionary tale journal October 2018