A cell multipole based domain decomposition algorithm for molecular dynamics simulation of systems of arbitrary shape

A domain decomposition algorithm for molecular dynamics simulation of atomic and molecular systems with arbitrary shape and non-periodic boundary conditions is described. The molecular dynamics program uses cell multipole method for efficient calculation of long range electrostatic interactions and a multiple time step method to facilitate bigger time steps. The system is enclosed in a cube and the cube is divided into a hierarchy of cells. The deepest level cells are assigned to processors such that each processor has contiguous cells and static load balancing is achieved by redistributing the cells so that each processor has approximately same number of atoms. The resulting domains have irregular shape and may have more than 26 neighbors. Atoms constituting bond angles and torsion angles may straddle more than two processors. An efficient strategy is devised for initial assignment and subsequent reassignment of such multiple-atom potentials to processors. At each step, computation is overlapped with communication greatly reducing the effect of communication overhead on parallel performance. The algorithm is tested on a spherical cluster of water molecules, a hexasaccharide and an enzyme both solvated by a spherical cluster of water molecules. In each case a spherical boundary containing oxygen atoms with only repulsive interactions is used to prevent evaporation of water molecules. The algorithm shows excellent parallel efficiency even for small number of cells/atoms per processor.

[1]  P. Jollès,et al.  Chitin and chitinases. , 1999 .

[2]  M. Karplus,et al.  SIMULATIONS OF MACROMOLECULES BY MULTIPLE TIME-STEP METHODS , 1995 .

[3]  A. Leach.  Molecular Modelling: Principles and Applications , 1996 .

[4]  Ming-Jing Hwang,et al.  Derivation of Class II Force Fields. 2. Derivation and Characterization of a Class II Force Field, CFF93, for the Alkyl Functional Group and Alkane Molecules , 1994 .

[5]  William Smith,et al.  Molecular dynamics on hypercube parallel computers , 1991 .

[6]  Arieh Warshel,et al.  A local reaction field method for fast evaluation of long‐range electrostatic interactions in molecular simulations , 1992 .

[7]  H. G. Petersen,et al.  An algorithm for the simulation of condensed matter which grows as the 3/2 power of the number of particles , 1988 .

[8]  Masuhiro Mikami,et al.  Computationally efficient canonical molecular dynamics simulations by using a multiple time‐step integrator algorithm combined with the particle mesh Ewald method and with the fast multipole method , 2000 .

[9]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[10]  B. Montgomery Pettitt,et al.  Structural and energetic effects of truncating long ranged interactions in ionic and polar fluids , 1985 .

[11]  B. Berne,et al.  Large scale simulation of macromolecules in solution: Combining the periodic fast multipole method with multiple time step integrators , 1997 .

[12]  K S Wilson,et al.  Crystal structure of a bacterial chitinase at 2.3 A resolution. , 1994, Structure.

[13]  Bernard Pettitt,et al.  Peptides in ionic solutions: A comparison of the Ewald and switching function techniques , 1991 .

[14]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[15]  Tjerk P. Straatsma,et al.  NWChem: Exploiting parallelism in molecular simulations , 2000 .

[16]  H. Schreiber,et al.  Molecular dynamics studies of solvated polypeptides: Why the cut-off scheme does not work , 1992 .

[17]  T. Darden,et al.  A Multipole-Based Algorithm for Efficient Calculation of Forces and Potentials in Macroscopic Period , 1996 .

[18]  B. Montgomery Pettitt,et al.  Effects of truncating long-range interactions in aqueous ionic solution simulations , 1988 .

[19]  J. Mccammon,et al.  Molecular dynamics with stochastic boundary conditions , 1982 .

[20]  Mark R. Wilson,et al.  Replicated data and domain decomposition molecular dynamics techniques for simulation of anisotropic potentials , 1997 .

[21]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[22]  J. W. Causey,et al.  Accelerated molecular dynamics simulation with the parallel fast multipole algorithm , 1992 .

[23]  A. H. Juffer,et al.  DYNAMIC SURFACE BOUNDARY-CONDITIONS - A SIMPLE BOUNDARY MODEL FOR MOLECULAR-DYNAMICS SIMULATIONS , 1993 .

[24]  Wataru Shinoda,et al.  Parallel molecular dynamics simulation: Implementation of PVM for a lipid membrane , 1999 .

[25]  J. Mccammon,et al.  Dynamics of Proteins and Nucleic Acids , 2018 .

[26]  Bruce Hendrickson,et al.  A new parallel method for molecular dynamics simulation of macromolecular systems , 1996 .

[27]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[28]  B. Berne,et al.  A Multiple-Time-Step Molecular Dynamics Algorithm for Macromolecules , 1994 .

[29]  K. Esselink,et al.  Efficient Parallel Implementation of Molecular Dynamics on a Toroidal Network. Part II. Multi-particle Potentials , 1993 .

[30]  Alan C. Belch,et al.  Molecular dynamics simulations of tips2 water restricted by a spherical hydrophobic boundary , 1985 .

[31]  Berend Smit,et al.  Efficient Parallel Implementation of Molecular Dynamics on a Toroidal Network. Part I. Parallelizing Strategy , 1993 .

[32]  Sharon Brunett,et al.  Molecular dynamics for very large systems on massively parallel computers: The MPSim program , 1997 .

[33]  Grant S. Heffelfinger,et al.  Parallel atomistic simulations , 2000 .

[34]  Benoît Roux,et al.  Dominant solvation effects from the primary shell of hydration: Approximation for molecular dynamics simulations , 1995 .

[35]  David Brown,et al.  A domain decomposition parallel processing algorithm for molecular dynamics simulations of systems of arbitrary connectivity. , 1997 .

[36]  Berend Smit,et al.  Understanding molecular simulation: from algorithms to applications , 1996 .

[37]  B. Berne,et al.  A new molecular dynamics method combining the reference system propagator algorithm with a fast multipole method for simulating proteins and other complex systems , 1995 .

[38]  O. Steinhauser,et al.  Taming cut-off induced artifacts in molecular dynamics studies of solvated polypeptides. The reaction field method. , 1992, Journal of molecular biology.

[39]  Martin Billeter,et al.  Point-centered domain decomposition for parallel molecular dynamics simulation , 2000 .

[40]  K. Schulten,et al.  Difficulties with multiple time stepping and fast multipole algorithm in molecular dynamics , 1997 .

[41]  W. Goddard,et al.  Atomic level simulations on a million particles: The cell multipole method for Coulomb and London nonbond interactions , 1992 .

[42]  H. C. Andersen.  Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

[43]  Jaroslav M. Ilnytskyi,et al.  A domain decomposition molecular dynamics program for the simulation of flexible molecules with an arbitrary topology of Lennard–Jones and/or Gay–Berne sites , 2001 .

[44]  Minoru Saito,et al.  Molecular dynamics simulations of proteins in water without the truncation of long-range Coulomb interactions , 1992 .

[45]  John Zahorjan,et al.  Parallel short-range molecular dynamics using the Ādhāra runtime system , 1997 .

[46]  Minoru Saito,et al.  Molecular dynamics simulations of proteins in solution: Artifacts caused by the cutoff approximation , 1994 .

[47]  Charles L. Brooks,et al.  The influence of long-range force truncation on the thermodynamics of aqueous ionic solutions , 1987 .

[48]  Ming-Jing Hwang,et al.  Derivation of class II force fields. I. Methodology and quantum force field for the alkyl functional group and alkane molecules , 1994, J. Comput. Chem..

[49]  Laxmikant V. Kale,et al.  NAMD2: Greater Scalability for Parallel Molecular Dynamics , 1999 .

[50]  John Zahorjan,et al.  Dynamic-domain-decomposition parallel molecular dynamics , 1997 .