Molecular dynamics simulations of materials
 Version 12
 by Alejandro Strachan
 Version 13
 by Michael Anderson
Deletions or items before changed
Additions or items after changed
1  == What is MD? ==  

2  
3  Molecular dynamics involves the numerical solution of the classical equations of motion (F=ma) for every single atom in material. The result is a very detailed description of the temporal evolution of the material: we obtain the position, velocity and force of every atom as a function of time. The total force acting on an atom comes from its interaction with other atoms and external fields and an accurate model to describe atomic interactions is critical in any MD simulation.  
4  If you are interested in MD read on; you will find more information, additional resources, simulation codes, and online tool that enables running live MD simulations using simply your webbrowser.  
5  
6  '''Advantages'''  
7  
8    As compared with other materials simulation techniques MD exhibits several advantages: The only inputs to the simulation are the interatomic potential used to compute forces on atoms and the initial configuration; all other properties and processes are derived. For example, you cannot dialin the melting temperature, or the character of the slip planes that will be activated during deformation, how a dislocation might look, nucleate or move, or how phonons interact with each other or with other defects. Consequently, phenomena like size effects in the thermal or mechanical response of materials are naturally captured by MD making a power technique in nanoscience and nanotechnology.

+  As compared with other materials simulation techniques, MD exhibits several advantages: The only inputs to the simulation are the interatomic potential used to compute forces on atoms and the initial configuration; all other properties and processes are derived. For example, you cannot dialin the melting temperature, or the character of the slip planes that will be activated during deformation, how a dislocation might look, nucleate or move, or how phonons interact with each other or with other defects. Consequently, phenomena like size effects in the thermal or mechanical response of materials are naturally captured by MD making a power technique in nanoscience and nanotechnology.

9  
10  '''Approximations'''  
11  
12  MD is not without approximations and it is important to understand them in order to design meaningful simulations and correctly interpret their results. There are two key approximations:  
13  
14    # Atomic interactions are now known exactly and an accurate model to describe them is critical in MD simulations of realistic materials. Atomic forces can be obtained from ab initio electronic structure calculations (such as density functional theory); these methods provide accurate results (not exact) but are computationally intensive. Most MD simulation use of interatomic potentials (or force fields) that describe atomic interactions in terms of functions parameterized to describe specific materials using experimental or ab initio data. These potentials are less accurate and less transferable than ab initio calculations but are computationally less intensive enabling largescale simulations (up to billions of atoms in current petaFLOP supercomputers).

+  # Atomic interactions are now known exactly, and an accurate model to describe them is critical in MD simulations of realistic materials. Atomic forces can be obtained from ''ab initio'' electronic structure calculations (such as density functional theory); these methods provide accurate results (not exact) but are computationally intensive. Most MD simulation use of interatomic potentials (or force fields) that describe atomic interactions in terms of functions parameterized to describe specific materials using experimental or ''ab initio'' data. These potentials are less accurate and less transferable than ''ab initio'' calculations but are computationally less intensive enabling largescale simulations (up to billions of atoms in current petaFLOP supercomputers).

15  # Being based on classical mechanics, MD does not capture quantum effects on ionic dynamics. These effects originate when atomic vibrations around a minimum lead confinement of position and velocity that violates Heisenberg’s uncertainty principle are important when simulations are performed at temperatures below or near the Debye temperature of the material.  
16  
17  == Learn more about MD at nanoHUB.org ==  
18  
19  * ''An Introduction to Molecular Dynamics.'' A lecture series by Alejandro Strachan. Video, voicedover presentations and notes: http://nanohub.org/resources/5838  
20  * ''Overview of Computational Nanoscience: a UC Berkeley Course.'' By Jeffrey C Grossman, Elif Ertekin. Class notes in pdf format: https://nanohub.org/resources/3944  
21  * ''Short Course on Molecular Dynamics Simulation.'' By Ashlie Martini. Class notes in pdf format: https://nanohub.org/resources/7570  
22  
23  == Online MD simulations at nanoHUB.org ==  
24  
25  The '''nanoMATERIALS simulation tool''' https://nanohub.org/resources/matsimtk enables running live MD simulations in nanoHUB.org. A simpletouse graphical interface allows users to set up the initial structure, simulation conditions (e.g. thermodynamic ensemble, temperature, strain rate), and output options; advanced options enable users to compute thermal conductivity using nonequilibrium simulations.  
26  
27  Learn more about the '''nanoMATERIALS simulation tool''' and run simulations via the following tutorials:  
28  
29  * ''Running MD on the nanoHUB: The nanoMATERIALS Simulation Toolkit.'' http://nanohub.org/resources/5843  
30  * ''Materials strength: does size matter? nanoMATERIALS simulation toolkit tutorial.'' https://nanohub.org/resources/2322  
31  
32  == MD in research and education ==  
33  
34  * ''First Principlesbased Atomistic and Mesoscale Modeling of Materials.'' https://nanohub.org/resources/434  
35  * ''Learning Module: Atomic picture of plastic deformation in metals.'' Teach and learn about the atomic level mechanisms that govern plastic deformation is metals: https://nanohub.org/topics/LearningModulePlasticityMD  
36  
37  == Other resources ==  
38  
39  '''Some of my favorite books and publications on MD'''  
40    * 
+  * Allen, M. P. and D. J. Tildesley. 1989. ''Computer Simulation of Liquids''. Oxford (UK): Oxford University Press.

41    * ''Understanding Molecular Simulation: From Algorithms to Applications'' 
+  * Frenkel, Daan and Berend Smit. 2001. ''Understanding Molecular Simulation: From Algorithms to Applications''. 2nd ed. San Diego: Academic Press..

42  
43  '''Downloadable, free of charge, MD codes'''  
44  * LAMMPS, the Largescale Atomic/Molecular Massively Parallel Simulator from Sandia National Labs, is available at: http://lammps.sandia.gov/ .  
45  * NAMD, from the University of Illinois at Urbana Champain, available at: http://www.ks.uiuc.edu/Research/namd/  
46  * GROMACS, GROningen MAchine for Chemical Simulations, is available at: http://www.gromacs.org/  
47  '''Related codes and resources'''  
48  * VMD, molecular visualization, http://www.ks.uiuc.edu/Research/vmd/ 