SANDIA REPORT
SAND 2006-7085
Unlimited Release
Printed November 2006
Substructured Multibody Molecular
Dynamics
Paul S. Crozier, Andrei I. Drăgănescu, Gary S. Grest, Ahmed E. Ismail,
Richard B. Lehoucq, Rudranarayan M. Mukherjee, Steven J. Plimpton,
Mark J. Stevens, and Thomas B. Woolf
Prepared by
Sandia National Laboratories
Albuquerque, New Mexico 87185 and Livermore, California 94550
Sandia is a multiprogram laboratory operated by Sandia Corporation,
a Lockheed Martin Company, for the United States Department of Energy’s
National Nuclear Security Administration under Contract DE-AC04-94AL85000.
Approved for public release; further dissemination unlimited.
2
Issued by Sandia National Laboratories, operated for the United States Department of Energy by
Sandia Corporation.
NOTICE: This report was prepared as an account of work sponsored by an agency of the United
States Government. Neither the United States Government, nor any agency thereof, nor any of
their employees, nor any of their contractors, subcontractors, or their employees, make any
warranty, express or implied, or assume any legal liability or responsibility for the accuracy,
completeness, or usefulness of any information, apparatus, product, or process disclosed, or
represent that its use would not infringe privately owned rights. Reference herein to any specific
commercial product, process, or service by trade name, trademark, manufacturer, or otherwise,
does not necessarily constitute or imply its endorsement, recommendation, or favoring by the
United States Government, any agency thereof, or any of their contractors or subcontractors. The
views and opinions expressed herein do not necessarily state or reflect those of the United States
Government, any agency thereof, or any of their contractors.
Printed in the United States of America. This report has been reproduced directly from the best
available copy.
Available to DOE and DOE contractors from
U.S. Department of Energy
Office of Scientific and Technical Information
P.O. Box 62
Oak Ridge, TN 37831
Telephone: (865) 576-8401
Facsimile: (865) 576-5728
E-Mail: reports@adonis.osti.gov
Online ordering: http://www.osti.gov/bridge
Available to the public from
U.S. Department of Commerce
National Technical Information Service
5285 Port Royal Rd.
Springfield, VA 22161
Telephone: (800) 553-6847
Facsimile: (703) 605-6900
E-Mail: orders@ntis.fedworld.gov
Online order: http://www.ntis.gov/help/ordermethods.asp?loc=7-4-0#online
3
SAND 2006-7085
Unlimited Release
Printed November 2006
Substructured Multibody Molecular Dynamics
Paul S. Crozier
Multiscale Computational Materials Methods
Andrei I. Drăgănescu and Richard B. Lehoucq
Computational Math/Algorithms
Gary S. Grest
Surface & Interface Sciences
Ahmed E. Ismail
Perf. Assmt. & Dec. Anal.
Steven J. Plimpton
Computational Biology
Mark J. Stevens
Biomolecular Analys & Imaging
Sandia National Laboratories
P.O. Box 5800
Albuquerque, NM 87185
Rudranarayan M. Mukherjee
Department of M.A.N.E.
Rensselaer Polytechnic Institute
Troy, New York 12180
Thomas B. Woolf
Johns Hopkins University
School of Medicine
Baltimore, Maryland 21205
Abstract
We have enhanced our parallel molecular dynamics (MD) simulation software LAMMPS (Large-scale
Atomic/Molecular Massively Parallel Simulator, lammps.sandia.gov
) to include many new features for
accelerated simulation including articulated rigid body dynamics via coupling to the Rensselaer
Polytechnic Institute code POEMS (Parallelizable Open-source Efficient Multibody Software). We use
new features of the LAMMPS software package to investigate rhodopsin photoisomerization, and water
model surface tension and capillary waves at the vapor-liquid interface. Finally, we motivate the recipes
of MD for practitioners and researchers in numerical analysis and computational mechanics.
4
5
Contents
Contents 5
Chapter 1: Substructured molecular dynamics simulations using multibody dynamics
algorithms through LAMMPS-POEMS coupling 7
Abstract 7
1. Introduction 7
2. Modeling Approach 8
2.1 Kinematic Model 8
2.2. Generating the rigid body properties 10
3. Algorithm Overview 10
3.1. Mathematical Preliminaries 11
4. ()On Forward Dynamics Analysis 12
4.1. Recursive Kinematic Relationships 12
4.2. Triangularization 13
4.3. Back-Substitutions 14
5. Time Integration 14
6. Applications and Results 15
6.1. Water Box 15
6.2. Alanine Dipeptide 15
6.3. DNA simulations 16
6.4. Box of alkanes 16
6.5. C-Terminal of Ribosomal 16
6.6. C-Terminal of Rubisco 16
6.7. Rhodopsin 17
7. Software Development 17
8. Conclusions 19
Acknowledgments 19
References 20
Chapter 2: How a small change in retinal leads to G-protein activation: initial events
suggested by molecular dynamics calculations 21
Abstract 21
Introduction 21
Simulation method 22
Results 23
C11-C12 dihedral transition 23
Retinal – Ala 169 gap narrows 24
Helix tilt and kink angle transitions 24
PSB counterion switch 26
C3 loop 27
Large-scale structural changes 27
Discussion 27
Conclusions 33
Table I 35
Figures 36
References 46