LAMMPS for Molecular Dynamics Simulation

Size: px
Start display at page:

Download "LAMMPS for Molecular Dynamics Simulation"

Transcription

1 Gwangju Institute of Science and Technology (GIST) Super Computing & Collaboration Environment Technology Center (SCENT) LAMMPS for Molecular Dynamics Simulation Taekhee Ryu Molecular Modeling Laboratory (MML) School of Materials Science and Engineering (MSE) Gwangju Institute of Science and Technology (GIST)

2 Overview: What is LAMMPS Lesson 1: Preparing input file Lesson 2: Running LAMMPS Lesson 3: Viewing output file

3 What is LAMMPS Large-scale Atomic/Molecular Massively Parallel Simulation - By Dr. S. Plimpton ( Free open-source code (C++) Model System with only a few particles up to millions or billions (depends on computational power) Designed for Parallel applications Capable for Atomic, polyatomic, biological, metallic, granular system 3

4 LAMMPS doesn t Build Molecular systems Use a builder code or write out your own config file Can build lattices, grain boundaries, etc Assign force-field coefficients auto-magically Only hard for molecules, use another MD code Compute lots of diagnostics on-the-fly Often better left to post-processing Visualize your output Many packages do this VMD, Rasmol, AtomEye, Ensight,.. 4

5 Our research using LAMMPS Multiscale molecular modeling of Nano-device containing fullerene derivatives l Unimolecular rectifier Ø Single organic molecule between metals can function as rectif iers, one-way conductors of electrical current l Absorption model on Au(111) l Self-assembled monolayer 5

6 Our research using LAMMPS l P3HT/PCBM BHJ l Amorphous phase of PVDF Molecular dynamic simulation 6

7 Overview of LAMMPS simulation INPUT (**.cmd, in.**, data.**) RUNNING OUTPUT (log.***, dump.***) INPUT - Initial positions and velocities - Boundary conditions - Potential energy function - Etc. RUNNING - Molecular Dynamics Simulation OUTPUT - Trajectories - Forces, Energy - Temperature - Pressure - Etc. 7

8 LAMMPS Input Reads an input script (ASCⅡtext) One command per line Command name + arguments (in.***) atom_style molecular read_data water.data fix 1 all nve run Data files(typical components) (data.***) Summary of molecular connectivity Cell dimensions Force field parameters(pair, bond, angle, dihedral, improper) Atom info (molecule #, atom type #, charge, xyz coordinates) Bond info (bonds, angles, etc) Examples sub-dir has many input scripts 8

9 Example Problems in LAMMPS Crack: crack growth in a LJ crystal (2d) Flow: couette/poisseuille flow between walls (2d) Friction: rubbing of 2 irregular surfaces (2d) Indent: crystal response to spherical indenter (2d) Melt: rapid melt of 3d LJ system Micelle: self-assembly of tiny lipid molecules (2d) Min: energy minimization of LJ melt (2d) Obstacle: flow around obstacles (2d) Peptide: granular particle pour and flow (2d/3d) Shear: shear of a metal slab with void (quasi-3d) Peptide: dynamics of a small solvated peptide chain (5-mer) 9

10 Input script Full-atom model of 5-mer peptide CHARMM force-field SHAKE bond/angle constraints 2 fs timestep 10

11 Input script 1. initialization 2. Atom definition 3. Settings 4. Run a simulation 11

12 Input script - initialization 1. initialization Set parameters that need to be defined before atoms are created or read-in from a file Unit, dimension, newton, processors, boundary, atom style, atom modify 12

13 Input script - initialization Units 1. initialization sets the style of units used for a simulation Beginning of an input script lj, real, metal, si, cgs, electron 13

14 Input script - initialization atom_style 1. initialization Define what style of atoms to use in a simulation Beginning of an input script 14

15 Input script - initialization 1. initialization 2. Atom definition Boundary - set the style of boundaries for the global simulatio n box in each dimension 3. Settings Processors - Specify how processors are mapped as a 3d logic al grid to the global simulation box Dimension - Set the dimensionality of the simulation. By defaul t LAMMPS runs 3d simulations. 15

16 Input script - initialization 1. initialization 2. Atom definition pair_style set the formula uses to compute pairwise interactions 3. Settings 16

17 Input script - initialization r 1. initialization 2. Atom definition bond_style uses to compute bond interactions between pairs of atoms 3. Settings 17

18 Input script - initialization 1. initialization r r 2. Atom definition angle_style uses to compute angle interactions between triplets of atoms 3. Settings 18

19 Input script - initialization 1. initialization dihedral_style 2. Atom definition uses to compute dihedral interactions between quadruplets of atoms 3. Settings 19

20 Input script - initialization 1. initialization improper_style 2. Atom definition uses to compute dihedral interactions between quadruplets of atoms 20

21 Input script - initialization 1. initialization 2. Atom definition kspace_style Define a K-space solver for LAMMPS to use each timestep to compute long-range Coulombic interactions or long-range 1/r^ N interactions Style none or ewald or pppm or pppm/cg or The pppm style invokes a particle-particle particle-mesh solver (Hockney) which maps atom charge to a 3d mesh, uses 3d F FTs to solve Poisson's equation on the mesh, then interpolate s electric fields on the mesh points back to the atoms 21

22 Input script - atom definition 1. initialization 2. Atom definition Read them in from a data or restart file via the read data or read restart command Contain molecular topology information 3. Settings Create atoms on a lattice with no molecular topology Lattice, region, create box, create atoms 22

23 Input script Read_data - Read in a data file containing information LA MMPS needs to run a simulation (data.***) 1. initialization # of components in system Simulation box boundaries Mass 23

24 Input script Read_data - Read in a data file containing information LA MMPS needs to run a simulation (data.***) 1. initialization 24

25 Input script Read_data - Read in a data file containing information LA MMPS needs to run a simulation (data.***) 1. initialization 25

26 Input script Read_data - Read in a data file containing information LA MMPS needs to run a simulation (data.***) 1. initialization 26

27 Input script Read_data - Read in a data file containing information LA MMPS needs to run a simulation (data.***) 1. initialization 27

28 Input script Read_data - Read in a data file containing information LA MMPS needs to run a simulation (data.***) 1. initialization 28

29 Input script - atom definition read_data - Read in a data file containing information LAMMPS needs to run a simulation (data.***) 1. initialization 2. Atom definition read_restart - Read in a previously saved simulation from a restart file 29

30 Input script - atom definition Create atoms on a lattice with no molecular topology Lattice, region, create box, create atoms Lattice simply a set of points in space determined by a unit cell with basis atoms used to define a wide variety of crystallographic latt ices syntax example 30

31 Input script - atom definition Create atoms on a lattice with no molecular topology Lattice, region, create box, create atoms region This command defines a geometric region of sp ace create box This command creates a simulation box based on t he specified region create atoms This command creates atoms on a lattice, or a si ngle atom, or a random collection of atoms 31

32 Input script - settings 1. initialization 2. Atom definition 3. Settings 3. Run a simulation 32

33 Input script - settings neighbor This command sets parameters that affect the building of pairwise neighbor lists The style value selects what algorithm is used to build the list The bin style creates the list by binning which is an operation that scales linearly with N/P, the number of atoms per processor where N = total number of atoms and P = number of processors 33

34 Input script - settings Neigh_modify This command sets parameters that affect the building and use of pairwise neighbor lists The delay setting means never build a new list until at least N steps after the previous build 34

35 Input script - settings timestep Set the timestep size for subsequent molecular dynamics simulations 35

36 Input script - settings thermo_style Set the style and content for printing thermodynamic data to the screen and log file. Style multi prints a multiple-line listing of thermodyna mic info thermo_style custom etotal ke temp pe ebond eangle edihed eimp evdwl ecoul elong press 36

37 Input script - settings thermo Compute and print thermodynamic info on timesteps that are a multiple of N and at the beginning and end of a simulation 37

38 Input script - settings fix Set a fix that will be applied to a group of atoms 38

39 Input script - settings fix Set a fix that will be applied to a group of atoms 39

40 Input script - settings group Identify a collection of atoms as belonging to a group The example above adds all atoms with IDs from 1 t o 12 to the group named peptide type 40

41 Input script - settings dump snapshot of atom quantities to one or more files every N timesteps in one of several styles 41

42 Input script 1. initialization 2. Atom definition 3. Settings run Run or continue dynamics for a specified number of timesteps 3. Run a simulation 42

43 INPUT (**.cmd, in.**, data.**) Tutorial RUNNING OUTPUT (log.***, dump.***) VMD Data.***. 43

44 VMD Data.***. Tutorial 44

45 VMD Data.***. Tutorial 45

46 VMD Data.***. Tutorial 46

47 VMD Data.***. Tutorial 47

48 VMD Data.***. Tutorial 48

49 Hand-on exercise Download the files Xming ( VMD ( Putty (SSH Access Client Program) Input file (in.peptide) Topology file (data.peptide) Executable file (lammps_kigi.cmd) Analysis file (thermo_extract.c) Install the programs Xming, Putty 49

50 Hand-on exercise PuTTY SSH Remote Access Client Program putty 50

51 Hand-on exercise Connect to the GIST supercomputer (kigi) localhost:0.0 51

52 Hand-on exercise Connect to the GIST supercomputer (kigi) Login01.plsi.or.kr, Port-22 52

53 Hand-on exercise Connect to the GIST supercomputer (kigi) $ ssh kigi-lg01 à connect to the kigi $ cp lammps.tar /pwork01/edun/pedunxx $ cd /pwokr01/edun/pedunxx $ tar xvf lammps.tar $ cd lammps 53

54 Hand-on exercise Lammps_kigi.cmd -> executable file Job name Job-time Input file 54

55 Example-1 Hand-on exercise Full-atom model of 5-mer peptide CHARMM force-field SHAKE bond/angle constraints 2 fs timestep 55

56 Hand-on exercise Running LAMMPS on PLSI supercomputer $ llsubmit lammps_kigi.cmd # submit the job $ llcancel job_id # cancel the job $ llq # job status $ llstatus # system status 56

57 Output There are two basic types of LAMMPS output. 1. Thermodynamic output (log.***) 1. A list of quantities printed every few time steps to the screen and log file 2. Dump files (dump.***) 1. Snapshots of atoms and various per-atom values and are w ritten at a specified frequency. 2. Default format is simple : id, type, x, y, z 3. VMD, Rasmol, AtomEye, XMOVIE 57

58 Output Example-1 (log.lammps) 58

59 Output Analysis of thermodynamic output (log.lammps) thermo_extract.c To run this program please compile it as -> gcc -O3 thermo_extract.c -lm -o thermo_extract -p for thermo parameter to extract -s if the output format is single line, or -m is the output format is multi line -> thermo_extract p PotEng -m log.lammps >&PotEng.tot 59

60 Output Analysis of thermodynamic output (log.***) -> gnuplot -> plot PotEng u 1:2 w l 60

61 Visualizer AtomEye - Rasmol - XMOVIE Auxiliary tool distributed with LAMMPS Very fast, simple viz 2d projection of 3d systems /home01/applic/lammps/lammps-18apr10/intel/ssh/tools/x movie/xmovie scale dump.*** VMD - Windows, Linux, MacOS X 3d hi-quality viz 61

62 Additional tools Amber2lmp AMBER <-> LAMMPS Binary2txt binary dump file -> ASCⅡ test files Ch2lmp CHARMM <-> LAMMPS Chain create the chains and solvent Lmp2arc LAMMPS -> Accelrys s insight MD Lamp2cfg LAMMPS -> AtomEYE(*.cfg) Lmp2vmd LAMMPS -> VMD Msi2lmp LAMMPS -> Accelrys s insight MD Python, micelle2d, matlab, ipp, emacs, eff, eam database, data 2xmovie, creatatoms 62

63 Thank you 63

Molecular Dynamics Simulations

Molecular Dynamics Simulations Molecular Dynamics Simulations Dr. Kasra Momeni www.knanosys.com Outline LAMMPS AdHiMad Lab 2 What is LAMMPS? LMMPS = Large-scale Atomic/Molecular Massively Parallel Simulator Open-Source (http://lammps.sandia.gov)

More information

MAGICS Workshop Thermal Conductivity Hands-On Session

MAGICS Workshop Thermal Conductivity Hands-On Session MAGICS Workshop Thermal Conductivity Hands-On Session This work was supported as part of the Computational Materials Sciences Program funded by the U.S. Department of Energy, Office of Science, Basic Energy

More information

Table of Contents. Table of Contents LAMMPS trajectories into QuantumATK

Table of Contents. Table of Contents LAMMPS trajectories into QuantumATK Table of Contents Table of Contents LAMMPS trajectories into QuantumATK Preparing a LAMMPS script Running LAMMPS Importing to QuantumATK Postprocessing LAMMPS trajectories References 1 2 2 3 3 5 8 QuantumATK

More information

Molecular Dynamics Simulation of Nanometric Machining Under Realistic Cutting Conditions Using LAMMPS

Molecular Dynamics Simulation of Nanometric Machining Under Realistic Cutting Conditions Using LAMMPS Molecular Dynamics Simulation of Nanometric Machining Under Realistic Cutting Conditions Using LAMMPS Rapeepan Promyoo Thesis Presentation Advisor: Dr. Hazim El-Mounayri Department of Mechanical Engineering

More information

LAMMPS Performance Benchmark on VSC-1 and VSC-2

LAMMPS Performance Benchmark on VSC-1 and VSC-2 LAMMPS Performance Benchmark on VSC-1 and VSC-2 Daniel Tunega and Roland Šolc Institute of Soil Research, University of Natural Resources and Life Sciences VSC meeting, Neusiedl am See, February 27-28,

More information

JASS Modeling and visualization of molecular dynamic processes

JASS Modeling and visualization of molecular dynamic processes JASS 2009 Konstantin Shefov Modeling and visualization of molecular dynamic processes St Petersburg State University, Physics faculty, Department of Computational Physics Supervisor PhD Stepanova Margarita

More information

Mechanical Characterization via Full Atomistic Simulation: Applications to Nanocrystallized Ice

Mechanical Characterization via Full Atomistic Simulation: Applications to Nanocrystallized Ice Mechanical Characterization via Full Atomistic Simulation: Applications to Nanocrystallized Ice A thesis presented By Arvand M.H. Navabi to The Department of Civil and Environmental Engineering in partial

More information

Modeling and visualization of molecular dynamic processes

Modeling and visualization of molecular dynamic processes JASS 2009 Konstantin Shefov Modeling and visualization of molecular dynamic processes St-Petersburg State University Physics faculty Department of Computational Physics Supervisor PhD Stepanova Margarita

More information

Introduction to Benchmark Test for Multi-scale Computational Materials Software

Introduction to Benchmark Test for Multi-scale Computational Materials Software Introduction to Benchmark Test for Multi-scale Computational Materials Software Shun Xu*, Jian Zhang, Zhong Jin xushun@sccas.cn Computer Network Information Center Chinese Academy of Sciences (IPCC member)

More information

Winmostar tutorial LAMMPS Melting point V X-Ability Co,. Ltd. 2017/8/17

Winmostar tutorial LAMMPS Melting point V X-Ability Co,. Ltd. 2017/8/17 Winmostar tutorial LAMMPS Melting point V7.025 X-Ability Co,. Ltd. question@winmostar.com Contents Configure I. Build solid phase II. Equilibration of solid phase III. Equilibration of liquid phase IV.

More information

2008 Biowerkzeug Ltd.

2008 Biowerkzeug Ltd. 2008 Biowerkzeug Ltd. 1 Contents Summary...3 1 Simulation...4 1.1 Setup...4 1.2 Output...4 2 Settings...5 3 Analysis...9 3.1 Setup...9 3.2 Input options...9 3.3 Descriptions...10 Please note that we cannot

More information

DISCRETE ELEMENT METHOD Basic Theory and Algorithms

DISCRETE ELEMENT METHOD Basic Theory and Algorithms Introduction to Dynamical Particle Simulations DISCRETE ELEMENT METHOD Basic Theory and Algorithms Daniel Queteschiner* Christoph Kloss** *daniel.queteschiner@gmx.at **christoph.kloss@cfdem.com Discrete

More information

Hands-on : Model Potential Molecular Dynamics

Hands-on : Model Potential Molecular Dynamics Hands-on : Model Potential Molecular Dynamics OUTLINE 0. DL_POLY code introduction 0.a Input files 1. THF solvent molecule 1.a Geometry optimization 1.b NVE/NVT dynamics 2. Liquid THF 2.a Equilibration

More information

Homework 2: Structure and Defects in Crystalline and Amorphous Materials

Homework 2: Structure and Defects in Crystalline and Amorphous Materials Homework 2: Structure and Defects in Crystalline and Amorphous Materials Muhammad Ashraful Alam Network of Computational Nanotechnology Discovery Park, Purdue University. In Lectures 4-6, we have discussed

More information

Chapter 3 Getting started with molecular dynamics modeling

Chapter 3 Getting started with molecular dynamics modeling Chapter 3 Getting started with molecular dynamics modeling Abstract In this chapter we provide a quick introduction to molecular dynamics modeling. In molecular dynamics the motion of a set of atoms is

More information

Applications of Molecular Dynamics

Applications of Molecular Dynamics June 4, 0 Molecular Modeling and Simulation Applications of Molecular Dynamics Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo Tohru

More information

EAM. ReaxFF. PROBLEM B Fracture of a single crystal of silicon

EAM. ReaxFF. PROBLEM B Fracture of a single crystal of silicon PROBLEM B Fracture of a single crystal of silicon This problem set utilizes a new simulation method based on Computational Materials Design Facility (CMDF) to model fracture in a brittle material, silicon.

More information

IAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.

IAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1. IAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.978 PDF) http://web.mit.edu/mbuehler/www/teaching/iap2006/intro.htm

More information

Lab 1: Empirical Energy Methods Due: 2/14/18

Lab 1: Empirical Energy Methods Due: 2/14/18 Lab 1: Empirical Energy Methods Due: 2/14/18 General remarks on scientific scripting Scientific scripting for managing the input and output data is an important component of modern materials computations,

More information

Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.

Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar. Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar. 25, 2002 Molecular Dynamics: Introduction At physiological conditions, the

More information

Computational Modeling of Li Diffusion Using Molecular Dynamics

Computational Modeling of Li Diffusion Using Molecular Dynamics Computational Modeling of Li Diffusion Using Molecular Dynamics A paper presented to the faculty of the Department of Physics of Wake Forest University in partial fulfillment of the requirements for graduation

More information

Preparation of Entangled Polymer Melts of Various Architecture for Coarse-grained Models

Preparation of Entangled Polymer Melts of Various Architecture for Coarse-grained Models Preparation of Entangled Polymer Melts of Various Architecture for Coarse-grained Models by Yelena R. Sliozberg and Jan W. Andzelm ARL-TR-5679 September 2011 Approved for public release; distribution unlimited.

More information

1 Introduction. command intended for command prompt

1 Introduction. command intended for command prompt Guest Lecture, Smith College, CS 334, BioInformatics 21 October 2008 GROMACS, Position Restrained MD, Protein Catalytic Activity Filip Jagodzinski 1 Introduction GROMACS (GROningen MAchine for Chemistry

More information

Tutorial I: IQ MOL and Basic DFT and MP2 Calculations 1 / 30

Tutorial I: IQ MOL and Basic DFT and MP2 Calculations 1 / 30 Tutorial I: IQ MOL and Basic DFT and MP2 Calculations Q-Chem User Workshop, Denver March 21, 2015 1 / 30 2 / 30 Introduction to IQMOL DFT and MP2 Calculations 3 / 30 IQMOL and Q-CHEM IQMOL is an open-source

More information

Universal Repulsive Contribution to the. Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information

Universal Repulsive Contribution to the. Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information Universal Repulsive Contribution to the Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information B. Shadrack Jabes, Dusan Bratko, and Alenka Luzar Department of Chemistry,

More information

The Molecular Dynamics Method

The Molecular Dynamics Method The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation

More information

NEW OPEN SOURCE SOFTWARE FOR BUILDING MOLECULAR DYNAMICS SYSTEMS. Welcome to MD Studio

NEW OPEN SOURCE SOFTWARE FOR BUILDING MOLECULAR DYNAMICS SYSTEMS. Welcome to MD Studio NEW OPEN SOURCE SOFTWARE FOR BUILDING MOLECULAR DYNAMICS SYSTEMS Welcome to MD Studio Background PROBLEM SPACE Molecular Dynamics (MD) Simulations are hampered by long build times for an initial MD cell.

More information

HIP Scientific Paper Example

HIP Scientific Paper Example HIP Scientific Paper Example Molecular Dynamics of Peptides on Graphitic Surfaces A. N. Camden a, D. E. Cossey b, Y.H. Ngo b, R. J. Berry b a Wright State University, 3460 Colonel Glenn Hwy, Dayton OH

More information

Protein Structures. 11/19/2002 Lecture 24 1

Protein Structures. 11/19/2002 Lecture 24 1 Protein Structures 11/19/2002 Lecture 24 1 All 3 figures are cartoons of an amino acid residue. 11/19/2002 Lecture 24 2 Peptide bonds in chains of residues 11/19/2002 Lecture 24 3 Angles φ and ψ in the

More information

Hands-on Course in Computational Structural Biology and Molecular Simulation BIOP590C/MCB590C. Course Details

Hands-on Course in Computational Structural Biology and Molecular Simulation BIOP590C/MCB590C. Course Details Hands-on Course in Computational Structural Biology and Molecular Simulation BIOP590C/MCB590C Emad Tajkhorshid Center for Computational Biology and Biophysics Email: emad@life.uiuc.edu or tajkhors@uiuc.edu

More information

Using the MMTSB Tool Set

Using the MMTSB Tool Set Using the MMTSB Tool Set MMTSB NIH Research Resource PSC Workshop 2003 Architecture of MMTSB Tool Set MMTSB Tool Set Package o Available from: (free for academic use) http://mmtsb.scripps.edu/software/mmtsbtoolset.html

More information

LAMMPS Simulation of a Microgravity Shear Cell 299r Progress Report 2 Taiyo Wilson. Introduction:

LAMMPS Simulation of a Microgravity Shear Cell 299r Progress Report 2 Taiyo Wilson. Introduction: Introduction: After our initial simulations outlined in the first progress report, we determined that a shear cell with a bidisperse composition of particles would be a promising system to explore. However,

More information

The Molecular Dynamics Method

The Molecular Dynamics Method H-bond energy (kcal/mol) - 4.0 The Molecular Dynamics Method Fibronectin III_1, a mechanical protein that glues cells together in wound healing and in preventing tumor metastasis 0 ATPase, a molecular

More information

EBIS-PIC 2D. 2D EBIS Simulation Code. Version 0.1. Copyright ( ) January FAR-TECH, Inc Science Center Dr., Ste 150. San Diego CA 92121

EBIS-PIC 2D. 2D EBIS Simulation Code. Version 0.1. Copyright ( ) January FAR-TECH, Inc Science Center Dr., Ste 150. San Diego CA 92121 EBIS-PIC 2D 2D EBIS Simulation Code Version 0.1 Copyright ( ) January 2017 by 10350 Science Center Dr., Ste 150 San Diego CA 92121 Phone 858-455-6655 Email support@far-tech.com URL http://www.far-tech.com

More information

MOLECULAR DYNAMICS SIMULATIONS OF FRICTION FORCE VERSUS LOAD

MOLECULAR DYNAMICS SIMULATIONS OF FRICTION FORCE VERSUS LOAD MOLECULAR DYNAMICS SIMULATIONS OF FRICTION FORCE VERSUS LOAD Ana-Camelia Pirghie Stefan cel Mare University of Suceava, Department of Mechanics and Technology, camelia.pirghie@fim.usv.ro Abstract: The

More information

CPMD Tutorial Atosim/RFCT 2009/10

CPMD Tutorial Atosim/RFCT 2009/10 These exercices were inspired by the CPMD Tutorial of Axel Kohlmeyer http://www.theochem.ruhruni-bochum.de/ axel.kohlmeyer/cpmd-tutor/index.html and by other tutorials. Here is a summary of what we will

More information

Tutorial. Getting started. Sample to Insight. March 31, 2016

Tutorial. Getting started. Sample to Insight. March 31, 2016 Getting started March 31, 2016 Sample to Insight CLC bio, a QIAGEN Company Silkeborgvej 2 Prismet 8000 Aarhus C Denmark Telephone: +45 70 22 32 44 www.clcbio.com support-clcbio@qiagen.com Getting started

More information

Winmostar tutorial LAMMPS Polymer modeling V X-Ability Co,. Ltd. 2017/7/6

Winmostar tutorial LAMMPS Polymer modeling V X-Ability Co,. Ltd. 2017/7/6 Winmostar tutorial LAMMPS Polymer modeling V7.021 X-Ability Co,. Ltd. question@winmostar.com 2017/7/6 Contents Configure I. Register a monomer II. Define a polymer III. Build a simulation cell IV. Execute

More information

Introduction Molecular Structure Script Console External resources Advanced topics. JMol tutorial. Giovanni Morelli.

Introduction Molecular Structure Script Console External resources Advanced topics. JMol tutorial. Giovanni Morelli. Gen 19th, 2017 1 2 Create and edit Display and view Mesurament and labelling Surface and Orbitals 3 4 from Database Protein Enzyme Crystal Structure and Unit Cell 5 Symmetry Animation General information

More information

Table of Contents. Table of Contents Opening a band gap in silicene and bilayer graphene with an electric field

Table of Contents. Table of Contents Opening a band gap in silicene and bilayer graphene with an electric field Table of Contents Table of Contents Opening a band gap in silicene and bilayer graphene with an electric field Bilayer graphene Building a bilayer graphene structure Calculation and analysis Silicene Optimizing

More information

Large-Scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Simulations of the Molecular Crystal αrdx

Large-Scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Simulations of the Molecular Crystal αrdx Large-Scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Simulations of the Molecular Crystal αrdx by Lynn B. Munday and Peter W. Chung ARL-TR-6579 August 2013 Approved for public release; distribution

More information

Table of Contents. Table of Contents Initialize from a converged state. Introduction Examples

Table of Contents. Table of Contents Initialize from a converged state. Introduction Examples Table of Contents Table of Contents Initialize from a converged state Introduction Examples Water molecule example Looping over parameters Stepping up in bias Special example: anti-parallel spin in MTJs

More information

LysinebasedTrypsinActSite. A computer application for modeling Chymotrypsin

LysinebasedTrypsinActSite. A computer application for modeling Chymotrypsin LysinebasedTrypsinActSite A computer application for modeling Chymotrypsin Version.2 May 2006 LysTAS A computer application for modeling chymotrypsin Version.2 May 2006 Table of Contents Page. Introduction

More information

VASP Tutorial. Nodo Lee. Ph. D. Candidate

VASP Tutorial. Nodo Lee. Ph. D. Candidate VASP Tutorial July 17, 2014 @ SCENT HPC Summer School @ GIST Nodo Lee Ph. D. Candidate Molecular Modeling Laboratory (MML) School of Materials Science and Engineering (MSE) Gwangju Institute of Science

More information

APBS electrostatics in VMD - Software. APBS! >!Examples! >!Visualization! >! Contents

APBS electrostatics in VMD - Software. APBS! >!Examples! >!Visualization! >! Contents Software Search this site Home Announcements An update on mailing lists APBS 1.2.0 released APBS 1.2.1 released APBS 1.3 released New APBS 1.3 Windows Installer PDB2PQR 1.7.1 released PDB2PQR 1.8 released

More information

Exercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions

Exercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions Exercise 2: Solvating the Structure HyperChem lets you place a molecular system in a periodic box of water molecules to simulate behavior in aqueous solution, as in a biological system. In this exercise,

More information

MENA 9520 FME Modelling Tutorial 5 ( )

MENA 9520 FME Modelling Tutorial 5 ( ) MENA 9520 FME Modelling Tutorial 5 (25.02.2011) Task 5: Understanding type of bonding using charge-density plots Exercise 5.1: Visualizing charge density in Si: 1. mkdir charge 2. Copy a converged CTRL

More information

CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields

CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields Department of Chemistry and Biochemistry, Concordia University! page 1 of 6 CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields INTRODUCTION The goal of this tutorial

More information

Supporting Information Soft Nanoparticles: Nano Ionic Networks of Associated Ionic Polymers

Supporting Information Soft Nanoparticles: Nano Ionic Networks of Associated Ionic Polymers Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2016 Supporting Information Soft Nanoparticles: Nano Ionic Networks of Associated Ionic Polymers Dipak

More information

Using GPUs for faster LAMMPS particle simulations

Using GPUs for faster LAMMPS particle simulations Using GPUs for faster LAMMPS particle simulations Paul S. Crozier Exploiting New Computer Architectures in Molecular Dynamics Simulations March 23, 2011 Sandia National Laboratories is a multi-program

More information

EDACA software demonstration for Molecular Dynamics simulations of EXAFS spectra

EDACA software demonstration for Molecular Dynamics simulations of EXAFS spectra EDACA software demonstration for Molecular Dynamics simulations of EXAFS spectra Alexei Kuzmin Dr.phys., Head of the Lab Institute of Solid State Physics, University of Latvia Kengaraga street 8, LV-1063

More information

QUANTUM CHEMISTRY WITH GAUSSIAN : A VERY BRIEF INTRODUCTION (PART 2)

QUANTUM CHEMISTRY WITH GAUSSIAN : A VERY BRIEF INTRODUCTION (PART 2) QUANTUM CHEMISTRY WITH GAUSSIAN : A VERY BRIEF INTRODUCTION (PART 2) TARAS V. POGORELOV AND MIKE HALLOCK SCHOOL OF CHEMICAL SCIENCES, UIUC This tutorial continues introduction to Gaussian [2]. Here we

More information

Dock Ligands from a 2D Molecule Sketch

Dock Ligands from a 2D Molecule Sketch Dock Ligands from a 2D Molecule Sketch March 31, 2016 Sample to Insight CLC bio, a QIAGEN Company Silkeborgvej 2 Prismet 8000 Aarhus C Denmark Telephone: +45 70 22 32 44 www.clcbio.com support-clcbio@qiagen.com

More information

Introduction to molecular dynamics

Introduction to molecular dynamics 1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.

More information

Molecular dynamics simulations of Si binding and diffusion on the native and thermal silicon oxide surfaces

Molecular dynamics simulations of Si binding and diffusion on the native and thermal silicon oxide surfaces The University of Toledo The University of Toledo Digital Repository Theses and Dissertations 2012 Molecular dynamics simulations of Si binding and diffusion on the native and thermal silicon oxide surfaces

More information

Lysozyme pka example - Software. APBS! >!Examples! >!pka calculations! >! Lysozyme pka example. Background

Lysozyme pka example - Software. APBS! >!Examples! >!pka calculations! >! Lysozyme pka example. Background Software Search this site Home Announcements An update on mailing lists APBS 1.2.0 released APBS 1.2.1 released APBS 1.3 released New APBS 1.3 Windows Installer PDB2PQR 1.7.1 released PDB2PQR 1.8 released

More information

Molecular modeling with InsightII

Molecular modeling with InsightII Molecular modeling with InsightII Yuk Sham Computational Biology/Biochemistry Consultant Phone: (612) 624 7427 (Walter Library) Phone: (612) 624 0783 (VWL) Email: shamy@msi.umn.edu How to run InsightII

More information

QuantumWise. QuantumWise is now part of Synopsys

QuantumWise. QuantumWise is now part of Synopsys Table of Contents Table of Contents NiSi2 Si interface Create the NiSi2/Si device The screening region Increase the length of the central region Set-up the calculation for the undoped device Dope the device

More information

Boris Mantisi. Contents. Context 2. How to install APoGe 3. How APoGe works 4. Input File 5. Crystal File 8. How to use APoGe 10

Boris Mantisi. Contents. Context 2. How to install APoGe 3. How APoGe works 4. Input File 5. Crystal File 8. How to use APoGe 10 Boris Mantisi Contents Context 2 How to install APoGe 3 How APoGe works 4 Input File 5 Crystal File 8 How to use APoGe 10 Context Polycrystalline structure plays an important role in the macroscopic properties

More information

Potential Energy (hyper)surface

Potential Energy (hyper)surface The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = " d dx U(x) Conformation

More information

Simulation of Molecules in LAMMPS

Simulation of Molecules in LAMMPS Simulation of Molecules in LAMMPS David Keffer Department of Materials Science & Engineering University of Tennessee, Knoxville date begun: February, 016 date last updated: February 3, 016 Table of Contents

More information

Winmostar tutorial LAMMPS Polymer Annealing V X-Ability Co., Ltd. 2018/01/15

Winmostar tutorial LAMMPS Polymer Annealing V X-Ability Co., Ltd. 2018/01/15 Winmostar tutorial LAMMPS Polymer Annealing V8.007 X-Ability Co., Ltd. question@winmostar.com 2018/01/15 Summary In this tutorial we will calculate glass transition temperature from the cooling process

More information

Supporting information of Self-Assembly, Structure and Pi- Conjugation at the Interface of Poly-3-Hexylthiophene and Carbon Nanotubes

Supporting information of Self-Assembly, Structure and Pi- Conjugation at the Interface of Poly-3-Hexylthiophene and Carbon Nanotubes Supporting information of Self-Assembly, Structure and Pi- Conjugation at the Interface of Poly-3-Hexylthiophene and Carbon Nanotubes Details of the molecular dynamics simulations and computational procedures

More information

Introduction to Hartree-Fock calculations in Spartan

Introduction to Hartree-Fock calculations in Spartan EE5 in 2008 Hannes Jónsson Introduction to Hartree-Fock calculations in Spartan In this exercise, you will get to use state of the art software for carrying out calculations of wavefunctions for molecues,

More information

Molecular Dynamics Simulation of Force- Controlled Nanoindentation

Molecular Dynamics Simulation of Force- Controlled Nanoindentation University of Arkansas, Fayetteville ScholarWorks@UARK Mechanical Engineering Undergraduate Honors Theses Mechanical Engineering 12-2015 Molecular Dynamics Simulation of Force- Controlled Nanoindentation

More information

The Molecular Dynamics Simulation Process

The Molecular Dynamics Simulation Process The Molecular Dynamics Simulation Process For textbooks see: M.P. Allen and D.J. Tildesley. Computer Simulation of Liquids.Oxford University Press, New York, 1987. D. Frenkel and B. Smit. Understanding

More information

Parallel Utility for Modeling of Molecular Aggregation

Parallel Utility for Modeling of Molecular Aggregation Parallel Utility for Modeling of Molecular Aggregation Peter Spijker, Bart Markvoort and Peter Hilbers BioModeling and bioinformatics Eindhoven University of Technology Eindhoven,The Netherlands p.spijker@pumma.nl

More information

How is molecular dynamics being used in life sciences? Davide Branduardi

How is molecular dynamics being used in life sciences? Davide Branduardi How is molecular dynamics being used in life sciences? Davide Branduardi davide.branduardi@schrodinger.com Exploring molecular processes with MD Drug discovery and design Protein-protein interactions Protein-DNA

More information

Journal Name ARTICLE. ESI for: CHARMM force field parameterization protocol for selfassembling peptide amphiphiles: The Fmoc moiety

Journal Name ARTICLE. ESI for: CHARMM force field parameterization protocol for selfassembling peptide amphiphiles: The Fmoc moiety Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2016Please do not adjust margins ESI for: Received 00th January 20xx, Accepted 00th

More information

What is Classical Molecular Dynamics?

What is Classical Molecular Dynamics? What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated

More information

FlexMD Manual. ADF Program System Release 2013

FlexMD Manual. ADF Program System Release 2013 FlexMD Manual ADF Program System Release 2013 Scientific Computing & Modelling NV Vrije Universiteit, Theoretical Chemistry De Boelelaan 1083; 1081 HV Amsterdam; The Netherlands WWW: www.scm.com E-mail:

More information

Structure of the HIV-1 Capsid Assembly by a hybrid approach

Structure of the HIV-1 Capsid Assembly by a hybrid approach Structure of the HIV-1 Capsid Assembly by a hybrid approach Juan R. Perilla juan@ks.uiuc.edu Theoretical and Computational Biophysics Group University of Illinois at Urbana-Champaign Virion HIV infective

More information

Chapter 32. Computations of radial distributions functions, PMFs and diffusion constants

Chapter 32. Computations of radial distributions functions, PMFs and diffusion constants Chapter 32 Computations of radial distributions functions, PMFs and diffusion constants After studying this chapter, you will be able to perform the 1) Calculations of radial distribution functions between

More information

Molecular Dynamics Simulations

Molecular Dynamics Simulations Molecular Dynamics Simulations Dr. Kasra Momeni www.knanosys.com Outline Long-range Interactions Ewald Sum Fast Multipole Method Spherically Truncated Coulombic Potential Speeding up Calculations SPaSM

More information

All-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010

All-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010 All-atom Molecular Mechanics Trent E. Balius AMS 535 / CHE 535 09/27/2010 Outline Molecular models Molecular mechanics Force Fields Potential energy function functional form parameters and parameterization

More information

Supervisors: Committee: Investigation of H-center diffusion in sodium chloride by molecular dynamics simulation

Supervisors: Committee: Investigation of H-center diffusion in sodium chloride by molecular dynamics simulation Investigation of H-center diffusion in sodium chloride by molecular dynamics simulation by Karun Datadien 18-7-2013 Supervisors: M.Sc. V. Veligura Dipl.-Ing. Dr. G. Hlawacek Committee: M.Sc. V. Veligura

More information

ArcGIS Enterprise: What s New. Philip Heede Shannon Kalisky Melanie Summers Sam Williamson

ArcGIS Enterprise: What s New. Philip Heede Shannon Kalisky Melanie Summers Sam Williamson ArcGIS Enterprise: What s New Philip Heede Shannon Kalisky Melanie Summers Sam Williamson ArcGIS Enterprise is the new name for ArcGIS for Server What is ArcGIS Enterprise ArcGIS Enterprise is powerful

More information

Fundamentals of Computational Science

Fundamentals of Computational Science Fundamentals of Computational Science Dr. Hyrum D. Carroll August 23, 2016 Introductions Each student: Name Undergraduate school & major Masters & major Previous research (if any) Why Computational Science

More information

Molecular dynamics simulation of graphene formation on 6H SiC substrate via simulated annealing

Molecular dynamics simulation of graphene formation on 6H SiC substrate via simulated annealing Molecular dynamics simulation of graphene formation on 6H SiC substrate via simulated annealing Yoon Tiem Leong @ Min Tjun Kit School of Physics Universiti Sains Malaysia 1 Aug 2012 Single layer graphene

More information

ChemNetworks. A Complex Network Analysis Tool For Chemical Systems. User s Manual. Version 1.0. November 2013

ChemNetworks. A Complex Network Analysis Tool For Chemical Systems. User s Manual. Version 1.0. November 2013 ChemNetworks A Complex Network Analysis Tool For Chemical Systems User s Manual Version 1.0 November 2013 Dr. Abdullah Ozkanlar and Prof. Aurora E. Clark Department of Chemistry, Washington State University

More information

Basics of Electronic-Structure Theory

Basics of Electronic-Structure Theory A guideline through the tutorial 23.06.2009 For each exercise in this tutorial, a separate directory should be created in your home directory in which the input files are generated and the calculation

More information

Speed-up of ATK compared to

Speed-up of ATK compared to What s new @ Speed-up of ATK 2008.10 compared to 2008.02 System Speed-up Memory reduction Azafulleroid (molecule, 97 atoms) 1.1 15% 6x6x6 MgO (bulk, 432 atoms, Gamma point) 3.5 38% 6x6x6 MgO (k-point sampling

More information

NAMD TUTORIAL. Unix/MacOSX Version. NAMD Developers: James Phillips, David Hardy

NAMD TUTORIAL. Unix/MacOSX Version. NAMD Developers: James Phillips, David Hardy University of Illinois at Urbana Champaign NIH Center for Macromolecular Modeling and Bioinformatics Beckman Institute Computational Biophysics Workshop NAMD TUTORIAL Unix/MacOSX Version NAMD Developers:

More information

Development of Molecular Dynamics Simulation System for Large-Scale Supra-Biomolecules, PABIOS (PArallel BIOmolecular Simulator)

Development of Molecular Dynamics Simulation System for Large-Scale Supra-Biomolecules, PABIOS (PArallel BIOmolecular Simulator) Development of Molecular Dynamics Simulation System for Large-Scale Supra-Biomolecules, PABIOS (PArallel BIOmolecular Simulator) Group Representative Hisashi Ishida Japan Atomic Energy Research Institute

More information

Vibrational Spectroscopy Practical 2:

Vibrational Spectroscopy Practical 2: Vibrational Spectroscopy Practical 2: The aims of today's session are to introduce you to running larger-scale CASTEP jobs on supercomputer clusters, creation of input files for and setup of CASTEP jobs,

More information

Department of Chemical Engineering University of California, Santa Barbara Spring Exercise 3. Due: Thursday, 5/3/12

Department of Chemical Engineering University of California, Santa Barbara Spring Exercise 3. Due: Thursday, 5/3/12 Department of Chemical Engineering ChE 210D University of California, Santa Barbara Spring 2012 Exercise 3 Due: Thursday, 5/3/12 Objective: To learn how to write & compile Fortran libraries for Python,

More information

Department of Chemical Engineering University of California, Santa Barbara Spring Exercise 2. Due: Thursday, 4/19/09

Department of Chemical Engineering University of California, Santa Barbara Spring Exercise 2. Due: Thursday, 4/19/09 Department of Chemical Engineering ChE 210D University of California, Santa Barbara Spring 2012 Exercise 2 Due: Thursday, 4/19/09 Objective: To learn how to compile Fortran libraries for Python, and to

More information

Supplementary Materials

Supplementary Materials Supplementary Materials Atomistic Origin of Brittle Failure of Boron Carbide from Large Scale Reactive Dynamics Simulations; Suggestions toward Improved Ductility Qi An and William A. Goddard III * Materials

More information

Contact Angle Of A Nano-Drop On A Heterogeneous Surface

Contact Angle Of A Nano-Drop On A Heterogeneous Surface Virginia Commonwealth University VCU Scholars Compass Theses and Dissertations Graduate School 2010 Contact Angle Of A Nano-Drop On A Heterogeneous Surface John Ritchie Virginia Commonwealth University

More information

Computational Structural Biology and Molecular Simulation. Introduction to VMD Molecular Visualization and Analysis

Computational Structural Biology and Molecular Simulation. Introduction to VMD Molecular Visualization and Analysis Computational Structural Biology and Molecular Simulation Introduction to VMD Molecular Visualization and Analysis Emad Tajkhorshid Department of Biochemistry, Beckman Institute, Center for Computational

More information

Project 2. Chemistry of Transient Species in Planetary Atmospheres: Exploring the Potential Energy Surfaces of CH 2 S

Project 2. Chemistry of Transient Species in Planetary Atmospheres: Exploring the Potential Energy Surfaces of CH 2 S Chemistry 362 Spring 2018 Dr. Jean M. Standard March 21, 2018 Project 2. Chemistry of Transient Species in Planetary Atmospheres: Exploring the Potential Energy Surfaces of CH 2 S In this project, you

More information

Materials that you may find helpful when working through this exercise

Materials that you may find helpful when working through this exercise Detailed steps illustrating how to use VASP on the Suns in Fitz 177 For use in lab: 11/10/2009 (Original file by Dr. Rachel Getman, 11/18/2007. Editted for use by Dorrell McCalman 11/09/2009.) Note on

More information

High-Performance Scientific Computing

High-Performance Scientific Computing High-Performance Scientific Computing Instructor: Randy LeVeque TA: Grady Lemoine Applied Mathematics 483/583, Spring 2011 http://www.amath.washington.edu/~rjl/am583 World s fastest computers http://top500.org

More information

Accelerating Three-Body Molecular Dynamics Potentials Using NVIDIA Tesla K20X GPUs. GE Global Research Masako Yamada

Accelerating Three-Body Molecular Dynamics Potentials Using NVIDIA Tesla K20X GPUs. GE Global Research Masako Yamada Accelerating Three-Body Molecular Dynamics Potentials Using NVIDIA Tesla K20X GPUs GE Global Research Masako Yamada Overview of MD Simulations Non-Icing Surfaces for Wind Turbines Large simulations ~ 1

More information

LAMMPS for Dummies. Once you have a working machine running, download the most recent stable version of LAMMP from

LAMMPS for Dummies. Once you have a working machine running, download the most recent stable version of LAMMP from LAMMPS for Dummies G.A. Frank and I. Sticco Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Pabellón I, Ciudad Universitaria, 1428 Buenos Aires, Argentina.

More information

Topics in Soft Matter Physics Theory and Modeling Introduction to Computer Modeling & Lab Set-up

Topics in Soft Matter Physics Theory and Modeling Introduction to Computer Modeling & Lab Set-up Topics in Soft Matter Physics Theory and Modeling Introduction to Computer Modeling & Lab Set-up Topics Today VMD Tutorial Molecular Simulation Monte Carlo (MC) Method Molecular Dynamics (MD) Method Euler

More information

CGWAS 2013: Simulating Stellar Collapse

CGWAS 2013: Simulating Stellar Collapse CGWAS 2013: Simulating Stellar Collapse C. D. Ott July 22, 2013 1 Introduction In this exercise, we will be using the code GR1D by O Connor & Ott (2010) [1] to simulate stellar collapse to a neutron star

More information

2.3 Modeling Interatomic Interactions Pairwise Potentials Many-Body Potentials Studying Biomolecules: The Force

2.3 Modeling Interatomic Interactions Pairwise Potentials Many-Body Potentials Studying Biomolecules: The Force Contents 1 Introduction to Computational Meso-Bio-Nano (MBN) Science and MBN EXPLORER.... 1 1.1 Meso-Bio-Nano Science: A Novel Field of Interdisciplinary Research.... 1 1.1.1 Structure and Dynamics of

More information

Chem 253. Tutorial for Materials Studio

Chem 253. Tutorial for Materials Studio Chem 253 Tutorial for Materials Studio This tutorial is designed to introduce Materials Studio 7.0, which is a program used for modeling and simulating materials for predicting and rationalizing structure

More information

"C:\Program Files\ANSYS Inc\v190\CFX\bin\perllib\cfx5solve.pl" -batch -ccl runinput.ccl -fullname "Fluid Flow CFX_002"

C:\Program Files\ANSYS Inc\v190\CFX\bin\perllib\cfx5solve.pl -batch -ccl runinput.ccl -fullname Fluid Flow CFX_002 This run of the CFX Release 19.0 Solver started at 19:06:17 on 05 Jun 2018 by user ltval on DESKTOP-JCG0747 (intel_xeon64.sse2_winnt) using the command: "C:\Program Files\ANSYS Inc\v190\CFX\bin\perllib\cfx5solve.pl"

More information