(Difference between revisions)
Jump to: navigation, search
(LAMMPS: Added video link)
(LAMMPS: Added video link)
Line 34: Line 34:
** [[LAMMPS_Interstitial_Formation_Energy | How to calculate interstitial formation energy for copper]]
** [[LAMMPS_Interstitial_Formation_Energy | How to calculate interstitial formation energy for copper]]
** [[LAMMPS reactive deformation of a single polyethylene chain | Reactive Deformation of a Polymer Chain ]]
** [[LAMMPS reactive deformation of a single polyethylene chain | Reactive Deformation of a Polymer Chain ]]
** Video tutorials for LAMMPS - [ Intro]; [ Input file Part 1]
** Video tutorials for LAMMPS - [ Intro] [ Input file Part 1] [ Input file Part 2]
== Visualization ==
== Visualization ==

Revision as of 17:10, 23 March 2017




Tensile Loading of an Aluminum Single Crystal. Movie showing deformation of single crystal aluminum loaded in the <100> direction at a strain rate of 1010 s-1 and a temperature of 300 K.

The nanoscale material models are molecular dynamics codes and tools used to ascertain properties at the atomistic scale. These simulations generally use interatomic potentials, or force fields, developed using properties obtained from both electronic scale) calculations and experiments, and feed these results into higher scale models, such as dislocation dynamics at the microscale, or continuum models at the macroscale. To date, much of the research at the atomistic scale has focused on informing continuum models for multiscale modeling of metal and polymer material systems. This particular site contains production and research codes that have been developed both at CAVS and outside for performing and analyzing atomistic simulation results. The production codes have user's manuals and a theoretical manual and have been used in practice to solve complex atomistic problems at the nanoscale. The codes that are research codes have not enjoyed the wealth of application and might not have a user's manual or a theoretical manual. We caution the user that there is some risk in using the research version of the codes. Another resource for computational chemistry can be found at computational chemistry.

Finally, to garner more information about the information bridges between length scales go to the MSU Education page.


If you are just beginning with atomistic codes, we recommend that you familiarize yourself with LAMMPS, MATLAB (pre- and post-processing), and some of the visualization codes.


Here are some brief tutorials for learning to use LAMMPS.


This section shows links to visualization packages used at the atomistic scale. Of these, AtomEye, Ensight, OVITO, and VMD are most frequently used at CAVS. AtomEye, OVITO, and VMD are open source codes.

MEAM Parameter Calibration

MEAM Parameter Calibration (MPC) is a graphical MATLAB application for:

  1. interactive editing of MEAM library and parameter files,
  2. running LAMMPS with an input file containing the commands 'pair_style meam' and 'pair_coeff * * LIBRARY_FILE ELEMENTS PARAMETER_FILE ATOM_TYPES', and
  3. automatic calibration of user-specified MEAM parameters.

See the MEAM Parameter Calibration - MPC Wiki page.

Tutorial videos for using the MPC tool can be found here, here, and here.

Preprocessing & Postprocessing Codes

This section includes codes used for preprocessing and postprocessing atomistic results. This section can also include scripts used to generate initial structures for inclusion in molecular dynamics simulations. Additionally, this subsection will include examples of xyz coordinate files that can be used in conjunction with the LAMMPS read_data command to upload.

K-12 Projects

This project(s) is designed to help introduce high school students to STEM-related 'relevant' research in physics and materials science and engineering.

Material Models

Molecular Dynamics Codes

This section includes links to molecular dynamics codes. LAMMPS[1] (Large-scale Atomic/Molecular Massively Parallel Simulator) is commonly used for many molecular dynamics simulations related to metal and polymer systems at CAVS. LAMMPS' Fortran predecessor WARP can also be used for parallel molecular dynamics simulations. Last, DYNAMO is commonly used for MEAM (modified embedded atom method)[2] interatomic potential generation.

Interatomic Potentials available online

For more information on interatomic potential generation using electronic structure information, use the following links.

Atomistic Research

This section includes interatomic potential information for atomistic simulations. Embedded atom method[3] potentials can be found at the NIST Interatomic Potential website. A number of modified embedded atom method[2] potentials have been developed here at CAVS for lightweight metals and steel research. Some published and ongoing interatomic potential work at CAVS includes


ICME Overview of Mechanical Properties of the Lipid Bilayer during Traumatic Brain Injury





Movie showing dislocation nucleation from a Sigma 3 asymmetric tilt grain boundary.








Coarse Grain Simulations
Polymer Atomistic Research. Movie showing deformation of an amorphous polyethylene structure with 20 chains of 1000 monomers length. The strain rate is 1010 s-1 and the temperature is 100 K[12][13].

An example of tensile deformation in amorphous polyethylene using a united atom method potential.

All Atom Simulations

Reactive molecular dynamics simulation of hydrocarbon-based polymers, such as polyethylene and polypropylene, is now possible using the recently parameterized modified embedded-atom method (MEAM) potential for hydrocarbons (C/H system).[14] For more information about the potential, please visit the following link:

Python Based Testing of Atomistic Potentials


  1. S. Plimpton, "Fast Parallel Algorithms for Short-Range Molecular Dynamics," J. Comp. Phys., 117, 1-19 (1995).
  2. 2.0 2.1 Baskes, M.I. (1992). Modified embedded-atom potentials for cubic materials and impurities. Phys. Rev. B, 46, 2727 (
  3. Murray S. Daw, Stephen M. Foiles, Michael I. Baskes,(1993) The embedded-atom method: a review of theory and applications, Materials Science Reports, Volume 9, Issues 7-8, Pages 251-310. (
  4. 4.0 4.1 4.2 Tschopp, M. A., & McDowell, D.L. (2007). Structures and energies of Sigma3 asymmetric tilt grain boundaries in Cu and Al. Philosophical Magazine, 87, 3147-3173 (
  5. 5.0 5.1 5.2 Tschopp, M. A., & McDowell, D.L. (2007). Asymmetric tilt grain boundary structure and energy in copper and aluminum. Philosophical Magazine, 87, 3871-3892 (
  6. Spearot, D.E., Tschopp, M.A., Jacob, K.I., McDowell, D.L., "Tensile strength of <100> and <110> tilt bicrystal copper interfaces," Acta Materialia 55 (2007) p. 705-714 (
  7. Tschopp, M.A., Spearot, D.E., McDowell, D.L., "Atomistic simulations of homogeneous dislocation nucleation in single crystal copper," Modelling and Simulation in Materials Science and Engineering 15 (2007) 693-709 (
  8. 8.0 8.1 8.2 Tschopp, M.A., McDowell, D.L., "Influence of single crystal orientation on homogeneous dislocation nucleation under uniaxial loading," Journal of Mechanics and Physics of Solids 56 (2008) 1806-1830. (
  9. K. Solanki, M.F. Horstemeyer, M. I. Baskes, and H. Feng, Multiscale study of dynamic void collapse in single crystals, Mechanics of Materials Volume 37, Issues 2-3, February-March 2005, Pages 317-330
  10. Tang, T., Kim, S., & Horstemeyer, M. (2010). Fatigue Crack Growth in Magnesium Single Crystals under Cyclic Loading: Molecular Dynamics Simulation. Computational Materials Science, 48, 426., 48, 426-439 (
  11. Barrett, C.D., El Kadiri, H., Tschopp, M.A. (2011). Breakdown of the Schmid Law in Homogenous and Heterogenous Nucleation Events of Slip and Twinning in Magnesium. Journal of Mechanics and Physics of Solids, in review.
  12. 12.0 12.1 Hossain, D., Tschopp, M.A., Ward, D.K., Bouvard, J.L., Wang, P., Horstemeyer, M.F., "Molecular dynamics simulations of deformation mechanisms of amorphous polyethylene," Polymer, 51 (2010) 6071-6083.
  13. 13.0 13.1 Tschopp, M.A., Ward, D.K., Bouvard, J.L., Horstemeyer, M.F., "Atomic Scale Deformation Mechanisms of Amorphous Polyethylene under Tensile Loading," TMS 2011 Conference Proceedings, accepted.
  14. S Nouranian, MA Tschopp, SR Gwaltney, MI Baskes, and MF Horstemeyer, An Interatomic Potential for Saturated Hydrocarbons Based on the Modified Embedded-Atom Method, Physical Chemistry Chemical Physics 16 (13) (2014):6233-6249.


This category has the following 3 subcategories, out of 3 total.

Pages in category "Nanoscale"

The following 75 pages are in this category, out of 75 total.



C cont.








M cont.










Personal tools

Material Models