LAMMPS Performance Benchmark on VSC-1 and VSC-2

Size: px
Start display at page:

Download "LAMMPS Performance Benchmark on VSC-1 and VSC-2"

Transcription

1 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, 2012

2 Objectives LAMMPS features LAMMPS benchmark tests Parallel performance issues on VSC-1 and VSC-2 Memory performance VSC-1 VSC-2 comparison

3 LAMMPS features Large-scale Atomic/Molecular Massively Parallel Simulator classical molecular dynamics code for modeling o atomic o polymeric o biological o metallic o granular and coarse-grained o hybrid and mesoscopic systems implemented various force fields ensembles, constrains, boundary conditions, integrators multi-replica models (parallel tempering, ) pre- and post-processing specialized features (e.g. GC MC, peridynamics)

4 LAMMPS features runs on a single processor or in parallel runs efficiently in parallel using MPI technique highly portable C++ easily modified and extended runs on various platforms including GPU (CUDA & OpenCL) open-source code, distributed under GNU Public License developed at Sandia National Laboratories core group: S. Plimpton, A. Thomson & P. Crozier

5 LAMMPS in parallel spatial decomposition techniques to partition the simulation domain into small 3d sub domains, one of which is assigned to each processor for computational efficiency LAMMPS uses neighbor lists to keep track of nearby particles processors communicate and store "ghost" atom information for atoms that border their sub domain CPU δt {Comm} = communication {Bond} {Neigh} {Pair} {Other + Output} Compute bonded terms If needed update the neighbor list Compute short- and longrange interaction terms for energy/forces Loop: i = 1 to N Collect forces, time integration and update positions), adjust T/p, print/write output

6 AMMPS on VSC-1 and VSC-2 VSC-1 (Intel X5550) VSC-2 (Opteron 6132HE) LAMMPS version 20Feb oct2011 C++ compiler opt. flags icc (11.1) -O icc (12.1.2) -O MPI Openmpi Openmpi FFTW fftw fftw-3.3

7 LAMMPS benchmark tests In LAMMPS standard distribution 5 benchmark tests Short-range forces modeled with a cut off distance Chain: Polymer chain melt (coarse-grained, FENE/LJ potential) LJ: Lennard-Jones liquid (LJ pot.) EAM: EAM metallic solid (EAM pot.) Chute: granular chute flow (granular pot.) Long-range forces Rhodopsin: solvated rhodopsin protein (CHARMM) All five tests can be run as Fixed-size (default particles, 100 time steps) Scaled-size the size increases with increasing number of cores

8 LAMMPS benchmark tests our test based on our VSC project: molecular dynamics study of wetting of mineral surfaces (SiO2, clays, FeOOH) FF: CLAYFF for minerals and SPC/E for water θ E = E + E + E CLAYFF = total Coul LJ bond stretch e qiq j σ ij σ ij = + 4ε ij + k r r 4πε 0 i j rij i j r ij r ij ( ) 1 ij 0 2

9 LAMMPS benchmark tests Description Chain LJ Rhodo Our test bead-spring polymer melt of 100-mer chains Atomic fluid Rhodopsin protein in solvated lipid bilayer Kaolinite layer with water droplet of 500 water molecules # of atoms k/64k/128k 6940 FF FENE/LJ LJ CHARMM CLAYFF/SPC/E cutoff 2^(1/6) σ 2.5 σ 10 Å 12 / 40 Å Long-rage N/A N/A PPPM Ewald Ensemble NVE NVE NpT NVT T / p 1 kbt/ε 0 kbt/ε 300 K / 1 atm 300 K Time step τ τ 2 fs 1 fs Run time τ 5000 τ 200 ps 200 ps Output 100 steps at the end 50 steps 50 steps

10 AMMPS VSC-1 scalability 250 Ideal scaling Test benchmark (6940 atoms) Chain benchmark (32000 atoms) LJ benchmark (32000 atoms) 200 Rhodo benchmark (32000 atoms) Rhodo benchmark (64000 atoms) Rhodo benchmark ( atoms) Scalability # of cores

11 AMMPS VSC-2 scalability Scalability Ideal scaling Test benchmark (6940 atoms) Chain benchmark (32000 atoms) LJ benchmark (32000 atoms) Rhodo benchmark (32000 atoms) Rhodo benchmark (64000 atoms) Rhodo benchmark ( atoms) # of cores

12 AMMPS VSC-1 and 2 speed-up VSC-1 Rhodo 32k VSC-1 Rhodo 64k VSC-1 Rhodo 128k time [s] VSC-2 Rhodo 32k VSC-2 Rhodo 64k VSC-2 Rhodo 128k # of cores

13 MMPS VSC-1 and VSC-2 run time [s] Test Chain LJ k 32k # cores VSC-1 VSC-2 ratè VSC-1 VSC-2 rate VSC-1 VSC-2 rate Aver.: 0.47 Aver.: 0.72 Aver.: 0.84 Rhodo 32k 64k 128k # cores VSC-1 VSC-2 rate VSC-1 VSC-2 rate VSC-1 VSC-2 rate Aver.: 0.51 Aver.: 0.51 Aver.: 0.54

14 MPS VSC-1 and VSC-2 parallel efficiency Efficiency / % efficiency >80% for optimal sub-domain size of >1000 particles VSC-1 Rhodo 32k VSC-1 Rhodo 128k VSC-2 Rhodo 32k VSC-2 Rhodo 128k ln2(n) cores

15 MPS VSC-1 and VSC-2 parallel efficiency 120 VSC-1 vs VSC Efficiency / % VSC-1 Rhodo 32k VSC-1 Rhodo 128k VSC-2 Rhodo 32k VSC-2 Rhodo 128k ln2(n) cores ln2(n) cores

16 MMPS - memory per processor [MB] Test Chair LJ Rhodopsin k 32k 32k 64k 128k # cores VSC-1 VSC-2 VSC-1 VSC-2 VSC-1 VSC-2 VSC- 1 VSC-2 VSC-1 VSC-2 VSC-1 VSC Rhodopsin 32 k 64 k 128k

17 CONCLUSIONS LAMMPS performance boost depends on model: better parallel efficiency for long-range models LAMMPS shows a good parallel performance, however o parallel efficiency of LAMMPS varies from the size of the benchmark data o performance advantage extends as cluster size increases o LAMMPS scales better to more processors for larger systems LAMMPS runs faster on VSC-1 than on VSC-2 (~0.5 for long-range models, for short-range models)

Efficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers

Efficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers Efficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers Jaewoon Jung (RIKEN, RIKEN AICS) Yuji Sugita (RIKEN, RIKEN AICS, RIKEN QBiC, RIKEN ithes) Molecular Dynamics

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 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

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

RWTH Aachen University

RWTH Aachen University IPCC @ RWTH Aachen University Optimization of multibody and long-range solvers in LAMMPS Rodrigo Canales William McDoniel Markus Höhnerbach Ahmed E. Ismail Paolo Bientinesi IPCC Showcase November 2016

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

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

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

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

Computational Chemistry - MD Simulations

Computational Chemistry - MD Simulations Computational Chemistry - MD Simulations P. Ojeda-May pedro.ojeda-may@umu.se Department of Chemistry/HPC2N, Umeå University, 901 87, Sweden. May 2, 2017 Table of contents 1 Basics on MD simulations Accelerated

More information

Coarse-Grained Models!

Coarse-Grained Models! Coarse-Grained Models! Large and complex molecules (e.g. long polymers) can not be simulated on the all-atom level! Requires coarse-graining of the model! Coarse-grained models are usually also particles

More information

The Fast Multipole Method in molecular dynamics

The Fast Multipole Method in molecular dynamics The Fast Multipole Method in molecular dynamics Berk Hess KTH Royal Institute of Technology, Stockholm, Sweden ADAC6 workshop Zurich, 20-06-2018 Slide BioExcel Slide Molecular Dynamics of biomolecules

More information

Accelerating Three-Body Potentials using GPUs NVIDIA Tesla K20X

Accelerating Three-Body Potentials using GPUs NVIDIA Tesla K20X Using a Hybrid Cray Supercomputer to Model Non-Icing Surfaces for Cold- Climate Wind Turbines Accelerating Three-Body Potentials using GPUs NVIDIA Tesla K20X GE Global Research Masako Yamada Opportunity

More information

Kobe-Brown Simulation Summer School 2015 Project: DPD Simulation of a Membrane

Kobe-Brown Simulation Summer School 2015 Project: DPD Simulation of a Membrane Kobe-Brown Simulation Summer School 2015 Project: DPD Simulation of a Membrane Clark Bowman Karen Larson Yuji Funaki Ross Parker Tae Woo Kim Week 1: Introduction to Molecular Dynamics (MD) Computer simulation

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

A Nobel Prize for Molecular Dynamics and QM/MM What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential

More information

Supporting information

Supporting information Electronic Supplementary Material ESI for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2018 Supporting information Understanding three-body contributions to coarse-grained force

More information

Supplementary Information for: Controlling Cellular Uptake of Nanoparticles with ph-sensitive Polymers

Supplementary Information for: Controlling Cellular Uptake of Nanoparticles with ph-sensitive Polymers Supplementary Information for: Controlling Cellular Uptake of Nanoparticles with ph-sensitive Polymers Hong-ming Ding 1 & Yu-qiang Ma 1,2, 1 National Laboratory of Solid State Microstructures and Department

More information

Supporting Information

Supporting Information Projection of atomistic simulation data for the dynamics of entangled polymers onto the tube theory: Calculation of the segment survival probability function and comparison with modern tube models Pavlos

More information

Introduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012

Introduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012 Introduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012 K. Kremer Max Planck Institute for Polymer Research, Mainz Overview Simulations, general considerations

More information

FENZI: GPU-enabled Molecular Dynamics Simulations of Large Membrane Regions based on the CHARMM force field and PME

FENZI: GPU-enabled Molecular Dynamics Simulations of Large Membrane Regions based on the CHARMM force field and PME 211 IEEE International Parallel & Distributed Processing Symposium : GPU-enabled Molecular Dynamics Simulations of Large Membrane Regions based on the force field and PME Narayan Ganesan, Michela Taufer

More information

Structuring of hydrophobic and hydrophilic polymers at interfaces Stephen Donaldson ChE 210D Final Project Abstract

Structuring of hydrophobic and hydrophilic polymers at interfaces Stephen Donaldson ChE 210D Final Project Abstract Structuring of hydrophobic and hydrophilic polymers at interfaces Stephen Donaldson ChE 210D Final Project Abstract In this work, a simplified Lennard-Jones (LJ) sphere model is used to simulate the aggregation,

More information

Supplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization

Supplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization Supplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization Richard H. Gee * Naida Lacevic and Laurence E. Fried University of California Lawrence

More information

MD simulation of methane in nanochannels

MD simulation of methane in nanochannels MD simulation of methane in nanochannels COCIM, Arica, Chile M. Horsch, M. Heitzig, and J. Vrabec University of Stuttgart November 6, 2008 Scope and structure Molecular model for graphite and the fluid-wall

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

Multiscale simulations of complex fluid rheology

Multiscale simulations of complex fluid rheology Multiscale simulations of complex fluid rheology Michael P. Howard, Athanassios Z. Panagiotopoulos Department of Chemical and Biological Engineering, Princeton University Arash Nikoubashman Institute of

More information

LAMMPS for Molecular Dynamics Simulation

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

More information

Protein-mineral interactions: molecular dynamics simulations capture importance of variations in mineral surface composition and structure

Protein-mineral interactions: molecular dynamics simulations capture importance of variations in mineral surface composition and structure Supplemental Information for: Protein-mineral interactions: molecular dynamics simulations capture importance of variations in mineral surface composition and structure Amity Andersen, *,1 Patrick N. Reardon,

More information

Molecular Dynamics Simulation of a Nanoconfined Water Film

Molecular Dynamics Simulation of a Nanoconfined Water Film Molecular Dynamics Simulation of a Nanoconfined Water Film Kyle Lindquist, Shu-Han Chao May 7, 2013 1 Introduction The behavior of water confined in nano-scale environment is of interest in many applications.

More information

Advanced Molecular Dynamics

Advanced Molecular Dynamics Advanced Molecular Dynamics Introduction May 2, 2017 Who am I? I am an associate professor at Theoretical Physics Topics I work on: Algorithms for (parallel) molecular simulations including GPU acceleration

More information

Mechanical Properties of Tetra-Polyethylene and Tetra-Polyethylene Oxide Diamond Networks via Molecular Dynamics Simulations

Mechanical Properties of Tetra-Polyethylene and Tetra-Polyethylene Oxide Diamond Networks via Molecular Dynamics Simulations Supplemental Information Mechanical Properties of Tetra-Polyethylene and Tetra-Polyethylene Oxide Diamond Networks via Molecular Dynamics Simulations Endian Wang and Fernando A. Escobedo Table S1 Lennard-Jones

More information

Parallelization of Molecular Dynamics (with focus on Gromacs) SeSE 2014 p.1/29

Parallelization of Molecular Dynamics (with focus on Gromacs) SeSE 2014 p.1/29 Parallelization of Molecular Dynamics (with focus on Gromacs) SeSE 2014 p.1/29 Outline A few words on MD applications and the GROMACS package The main work in an MD simulation Parallelization Stream computing

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

Molecular dynamics simulation of Aquaporin-1. 4 nm

Molecular dynamics simulation of Aquaporin-1. 4 nm Molecular dynamics simulation of Aquaporin-1 4 nm Molecular Dynamics Simulations Schrödinger equation i~@ t (r, R) =H (r, R) Born-Oppenheimer approximation H e e(r; R) =E e (R) e(r; R) Nucleic motion described

More information

Diffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study

Diffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study Diffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study Julia Deitz, Yeneneh Yimer, and Mesfin Tsige Department of Polymer Science University of Akron

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

Interface Resistance and Thermal Transport in Nano-Confined Liquids

Interface Resistance and Thermal Transport in Nano-Confined Liquids 1 Interface Resistance and Thermal Transport in Nano-Confined Liquids Murat Barisik and Ali Beskok CONTENTS 1.1 Introduction...1 1.2 Onset of Continuum Behavior...2 1.3 Boundary Treatment Effects on Interface

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

Force Fields for Classical Molecular Dynamics simulations of Biomolecules. Emad Tajkhorshid

Force Fields for Classical Molecular Dynamics simulations of Biomolecules. Emad Tajkhorshid Force Fields for Classical Molecular Dynamics simulations of Biomolecules Emad Tajkhorshid Theoretical and Computational Biophysics Group, Beckman Institute Departments of Biochemistry and Pharmacology,

More information

MARTINI simulation details

MARTINI simulation details S1 Appendix MARTINI simulation details MARTINI simulation initialization and equilibration In this section, we describe the initialization of simulations from Main Text section Residue-based coarsegrained

More information

PDLAMMPS - made easy

PDLAMMPS - made easy PDLAMMPS - made easy R. Rahman 1, J. T. Foster 1, and S. J. Plimpton 2 1 The University of Texas at San Antonio 2 Sandia National Laboratory February 12, 2014 1 Peridynamic theory of solids The peridynamic

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

Lecture 11: Potential Energy Functions

Lecture 11: Potential Energy Functions Lecture 11: Potential Energy Functions Dr. Ronald M. Levy ronlevy@temple.edu Originally contributed by Lauren Wickstrom (2011) Microscopic/Macroscopic Connection The connection between microscopic interactions

More information

Charge equilibration

Charge equilibration Charge equilibration Taylor expansion of energy of atom A @E E A (Q) =E A0 + Q A + 1 @Q A 0 2 Q2 A @ 2 E @Q 2 A 0 +... The corresponding energy of cation/anion and neutral atom E A (+1) = E A0 + @E @Q

More information

MOLECULAR DYNAMIC SIMULATION OF WATER VAPOR INTERACTION WITH VARIOUS TYPES OF PORES USING HYBRID COMPUTING STRUCTURES

MOLECULAR DYNAMIC SIMULATION OF WATER VAPOR INTERACTION WITH VARIOUS TYPES OF PORES USING HYBRID COMPUTING STRUCTURES MOLECULAR DYNAMIC SIMULATION OF WATER VAPOR INTERACTION WITH VARIOUS TYPES OF PORES USING HYBRID COMPUTING STRUCTURES V.V. Korenkov 1,3, a, E.G. Nikonov 1, b, M. Popovičová 2, с 1 Joint Institute for Nuclear

More information

Reactive potentials and applications

Reactive potentials and applications 1.021, 3.021, 10.333, 22.00 Introduction to Modeling and Simulation Spring 2011 Part I Continuum and particle methods Reactive potentials and applications Lecture 8 Markus J. Buehler Laboratory for Atomistic

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

Extending Parallel Scalability of LAMMPS and Multiscale Reactive Molecular Simulations

Extending Parallel Scalability of LAMMPS and Multiscale Reactive Molecular Simulations September 18, 2012 Extending Parallel Scalability of LAMMPS and Multiscale Reactive Molecular Simulations Yuxing Peng, Chris Knight, Philip Blood, Lonnie Crosby, and Gregory A. Voth Outline Proton Solvation

More information

Bioengineering 215. An Introduction to Molecular Dynamics for Biomolecules

Bioengineering 215. An Introduction to Molecular Dynamics for Biomolecules Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD

More information

Crystal Polymorphism in Hydrophobically Nanoconfined Water. O. Vilanova, G. Franzese Universitat de Barcelona

Crystal Polymorphism in Hydrophobically Nanoconfined Water. O. Vilanova, G. Franzese Universitat de Barcelona Crystal Polymorphism in Hydrophobically Nanoconfined Water O. Vilanova, G. Franzese Universitat de Barcelona Water s phase diagram Phase diagram of bulk water presents up to 15 different phases Typical

More information

FORSCHUNGSZENTRUM JÜLICH GmbH Zentralinstitut für Angewandte Mathematik D Jülich, Tel. (02461)

FORSCHUNGSZENTRUM JÜLICH GmbH Zentralinstitut für Angewandte Mathematik D Jülich, Tel. (02461) FORSCHUNGSZENTRUM JÜLICH GmbH Zentralinstitut für Angewandte Mathematik D-52425 Jülich, Tel. (02461) 61-6402 Interner Bericht Thermalization of Large Fluid Systems consisting of Lennard-Jones Mitures Matthias

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

Principles and Applications of Molecular Dynamics Simulations with NAMD

Principles and Applications of Molecular Dynamics Simulations with NAMD Principles and Applications of Molecular Dynamics Simulations with NAMD Nov. 14, 2016 Computational Microscope NCSA supercomputer JC Gumbart Assistant Professor of Physics Georgia Institute of Technology

More information

Computer simulation methods (2) Dr. Vania Calandrini

Computer simulation methods (2) Dr. Vania Calandrini Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing

More information

Flora Tsourtou & Vlasis G. Mavrantzas. HPC for Innovation: when Science meets Industry, May, Barcelona, Spain

Flora Tsourtou & Vlasis G. Mavrantzas. HPC for Innovation: when Science meets Industry, May, Barcelona, Spain Optimized, Large-scale Monte Carlo and Molecular Dynamics Algorithms for modelling the self-organization of two classes of materials: semifluorinated alkanes & semiconducting polymers based on thiophenes

More information

The Computational Microscope

The Computational Microscope The Computational Microscope Computational microscope views at atomic resolution... Rs SER RER C E M RER N GA L... how living cells maintain health and battle disease John Stone Our Microscope is Made

More information

ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below

ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below Introduction In statistical physics Monte Carlo methods are considered to have started in the Manhattan project (1940

More information

Glass-Transition and Side-Chain Dynamics in Thin Films: Explaining. Dissimilar Free Surface Effects for Polystyrene and Poly(methyl methacrylate)

Glass-Transition and Side-Chain Dynamics in Thin Films: Explaining. Dissimilar Free Surface Effects for Polystyrene and Poly(methyl methacrylate) Supporting Information for Glass-Transition and Side-Chain Dynamics in Thin Films: Explaining Dissimilar Free Surface Effects for Polystyrene and Poly(methyl methacrylate) David D. Hsu, Wenjie Xia, Jake

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

A Brief Guide to All Atom/Coarse Grain Simulations 1.1

A Brief Guide to All Atom/Coarse Grain Simulations 1.1 A Brief Guide to All Atom/Coarse Grain Simulations 1.1 Contents 1. Introduction 2. AACG simulations 2.1. Capabilities and limitations 2.2. AACG commands 2.3. Visualization and analysis 2.4. AACG simulation

More information

Computational Modeling of Protein Dynamics. Yinghao Wu Department of Systems and Computational Biology Albert Einstein College of Medicine Fall 2014

Computational Modeling of Protein Dynamics. Yinghao Wu Department of Systems and Computational Biology Albert Einstein College of Medicine Fall 2014 Computational Modeling of Protein Dynamics Yinghao Wu Department of Systems and Computational Biology Albert Einstein College of Medicine Fall 2014 Outline Background of protein dynamics Basic computational

More information

Reactive Molecular Dynamics on Massively Parallel Heterogeneous Architectures

Reactive Molecular Dynamics on Massively Parallel Heterogeneous Architectures Reactive Molecular Dynamics on Massively Parallel Heterogeneous Architectures Sudhir B Kylasa, Hasan Metin Aktulga, Ananth Y Grama Abstract We present a parallel implementation of the ReaxFF force field

More information

Molecular Mechanics / ReaxFF

Molecular Mechanics / ReaxFF Molecular Dynamics simulations Lecture 09: Molecular Mechanics / ReaxFF Dr. Olli Pakarinen University of Helsinki Fall 2012 Lecture notes based on notes by Dr. Jani Kotakoski, 2010 CONTENTS Molecular mechanics

More information

Developing Monovalent Ion Parameters for the Optimal Point Charge (OPC) Water Model. John Dood Hope College

Developing Monovalent Ion Parameters for the Optimal Point Charge (OPC) Water Model. John Dood Hope College Developing Monovalent Ion Parameters for the Optimal Point Charge (OPC) Water Model John Dood Hope College What are MD simulations? Model and predict the structure and dynamics of large macromolecules.

More information

ABSTRACT. Xiaofei Ma, M.S., Low molecular weight dicarboxylic acids constitute a significant fraction of

ABSTRACT. Xiaofei Ma, M.S., Low molecular weight dicarboxylic acids constitute a significant fraction of ABSTRACT Title of Document: MOLECULAR DYNAMICS SIMULATION OF DICARBOXYLIC ACID COATED AQUEOUS AEROSOL: STRUCTURE AND PROCESSING OF WATER VAPOR Xiaofei Ma, M.S., 2010 Directed By: Professor Michael R. Zachariah,

More information

Distribution of chains in polymer brushes produced by a grafting from mechanism

Distribution of chains in polymer brushes produced by a grafting from mechanism SUPPLEMENTARY INFORMATION Distribution of chains in polymer brushes produced by a grafting from mechanism Andre Martinez, Jan-Michael Y. Carrillo, Andrey V. Dobrynin,, * and Douglas H. Adamson, * Department

More information

Molecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia

Molecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia Molecular Dynamics Simulations Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia 1 An Introduction to Molecular Dynamics Simulations Macroscopic properties

More information

WATER PERMEATION THROUGH GRAPHENE NANOSLIT BY MOLECULAR DYNAMICS SIMULATION

WATER PERMEATION THROUGH GRAPHENE NANOSLIT BY MOLECULAR DYNAMICS SIMULATION WATER PERMEATION THROUGH GRAPHENE NANOSLIT BY MOLECULAR DYNAMICS SIMULATION Taro Yamada 1 and Ryosuke Matsuzaki 2 1 Department of Mechanical Engineering, Tokyo University of Science, 2641 Yamazaki, Noda,

More information

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

Advantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations

Advantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations The Hilltop Review Volume 7 Issue 1 Winter 2014 Article 10 December 2014 Advantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations Tai-Hsien Wu Western Michigan University

More information

Autocorrelation Study for a Coarse-Grained Polymer Model

Autocorrelation Study for a Coarse-Grained Polymer Model Autocorrelation Study for a Coarse-Grained Polymer Model KaiK.Qi and andmichael M. Bachmann Bachmann he Center for CSP,UGA Simulational Physics he University of Georgia alk given at the 14th Leipzig Workshop

More information

Computation of non-bonded interactions: Part 1

Computation of non-bonded interactions: Part 1 Computation of non-bonded interactions: Part 1 Computation of non-local (non-bonded) interactions in the potential function E p non-local scales with the number of atoms as N. For the large molecular systems

More information

Molecular Dynamics. A very brief introduction

Molecular Dynamics. A very brief introduction Molecular Dynamics A very brief introduction Sander Pronk Dept. of Theoretical Physics KTH Royal Institute of Technology & Science For Life Laboratory Stockholm, Sweden Why computer simulations? Two primary

More information

Supporting Information for Solid-liquid Thermal Transport and its Relationship with Wettability and the Interfacial Liquid Structure

Supporting Information for Solid-liquid Thermal Transport and its Relationship with Wettability and the Interfacial Liquid Structure Supporting Information for Solid-liquid Thermal Transport and its Relationship with Wettability and the Interfacial Liquid Structure Bladimir Ramos-Alvarado, Satish Kumar, and G. P. Peterson The George

More information

Water models in classical simulations

Water models in classical simulations Water models in classical simulations Maria Fyta Institut für Computerphysik, Universität Stuttgart Stuttgart, Germany Water transparent, odorless, tasteless and ubiquitous really simple: two H atoms attached

More information

This semester. Books

This semester. Books Models mostly proteins from detailed to more abstract models Some simulation methods This semester Books None necessary for my group and Prof Rarey Molecular Modelling: Principles and Applications Leach,

More information

Efficient Molecular Dynamics on Heterogeneous Architectures in GROMACS

Efficient Molecular Dynamics on Heterogeneous Architectures in GROMACS Efficient Molecular Dynamics on Heterogeneous Architectures in GROMACS Berk Hess, Szilárd Páll KTH Royal Institute of Technology GTC 2012 GROMACS: fast, scalable, free Classical molecular dynamics package

More information

Using Molecular Dynamics to Compute Properties CHEM 430

Using Molecular Dynamics to Compute Properties CHEM 430 Using Molecular Dynamics to Compute Properties CHEM 43 Heat Capacity and Energy Fluctuations Running an MD Simulation Equilibration Phase Before data-collection and results can be analyzed the system

More information

Huge-Scale Molecular Dynamics Simulation of Multi-bubble Nuclei

Huge-Scale Molecular Dynamics Simulation of Multi-bubble Nuclei 1/20 Huge-Scale Molecular Dynamics Simulation of Multi-bubble Nuclei H. Watanabe ISSP, The M. Suzuki H. Inaoka N. Ito Kyushu University RIKEN AICS The, RIKEN AICS Outline 1. Introduction 2. Benchmark results

More information

Non-bonded interactions

Non-bonded interactions speeding up the number-crunching continued Marcus Elstner and Tomáš Kubař December 3, 2013 why care? number of individual pair-wise interactions bonded interactions proportional to N: O(N) non-bonded interactions

More information

Dynamic Scheduling for Work Agglomeration on Heterogeneous Clusters

Dynamic Scheduling for Work Agglomeration on Heterogeneous Clusters Dynamic Scheduling for Work Agglomeration on Heterogeneous Clusters Jonathan Lifflander, G. Carl Evans, Anshu Arya, Laxmikant Kale University of Illinois Urbana-Champaign May 25, 2012 Work is overdecomposed

More information

Molecular Mechanics Force Fields

Molecular Mechanics Force Fields Molecular Mechanics Force Fields Basic Premise If we want to study a protein, piece of DNA, biological membranes, polysaccharide, crystal la;ce, nanomaterials, diffusion in liquids, the number of electrons

More information

Ion-Gated Gas Separation through Porous Graphene

Ion-Gated Gas Separation through Porous Graphene Online Supporting Information for: Ion-Gated Gas Separation through Porous Graphene Ziqi Tian, Shannon M. Mahurin, Sheng Dai,*,, and De-en Jiang *, Department of Chemistry, University of California, Riverside,

More information

Virtual material Design

Virtual material Design Fraunhofer Institute for Algorithms and Scientific Computing SCAI Virtual material Design Foto: viastore systems 1 virtual Material Design Many of the challenges of the 21st century such as the development

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

Investigation of an Unusual Phase Transition Freezing on heating of liquid solution

Investigation of an Unusual Phase Transition Freezing on heating of liquid solution Investigation of an Unusual Phase Transition Freezing on heating of liquid solution Calin Gabriel Floare National Institute for R&D of Isotopic and Molecular Technologies, Cluj-Napoca, Romania Max von

More information

Enhancing Amber for Use on the Blue Waters High- Performance Computing Resource

Enhancing Amber for Use on the Blue Waters High- Performance Computing Resource Enhancing Amber for Use on the Blue Waters High- Performance Computing Resource Daniel R. Roe a, Jason M. Swails b, Romelia Salomon-Ferrer c, Adrian E. Roitberg b, and Thomas E. Cheatham III a * a. Department

More information

Folding WT villin in silico Three folding simulations reach native state within 5-8 µs

Folding WT villin in silico Three folding simulations reach native state within 5-8 µs Modeling of Cryo-EM Maps Workshop Baylor College of Medicine Klaus Schulten, U. Illinois at Urbana-Champaign Molecular Modeling Flexible Fitting 1: Introduction to Molecular Dynamics MD force field Equation

More information

ABSRACT. Department of Mechanical Engineering. Professor Peter W. Chung Department of Mechanical Engineering

ABSRACT. Department of Mechanical Engineering. Professor Peter W. Chung Department of Mechanical Engineering ABSRACT Title of thesis: WETTING OF GRAPHENE Joseph E. Andrews, Master of Science, 2016 Thesis directed by: Professor Siddhartha Das Department of Mechanical Engineering Professor Peter W. Chung Department

More information

Intermolecular Forces and Monte-Carlo Integration 열역학특수연구

Intermolecular Forces and Monte-Carlo Integration 열역학특수연구 Intermolecular Forces and Monte-Carlo Integration 열역학특수연구 2003.3.28 Source of the lecture note. J.M.Prausnitz and others, Molecular Thermodynamics of Fluid Phase Equiliria Atkins, Physical Chemistry Lecture

More information

Microscopic Picture of Aging in SiO 2 : A Computer Simulation

Microscopic Picture of Aging in SiO 2 : A Computer Simulation Microscopic Picture of Aging in SiO 2 : A Computer Simulation Katharina Vollmayr-Lee, Robin Bjorkquist, Landon M. Chambers Bucknell University & Göttingen 7 td tb r n (t) 6 R 5 4 3 2 ti t [ns] waiting

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

File ISM02. Dynamics of Soft Matter

File ISM02. Dynamics of Soft Matter File ISM02 Dynamics of Soft Matter 1 Modes of dynamics Quantum Dynamics t: fs-ps, x: 0.1 nm (seldom of importance for soft matter) Molecular Dynamics t: ps µs, x: 1 10 nm Brownian Dynamics t: ns>ps, x:

More information

Molecular Dynamics. What to choose in an integrator The Verlet algorithm Boundary Conditions in Space and time Reading Assignment: F&S Chapter 4

Molecular Dynamics. What to choose in an integrator The Verlet algorithm Boundary Conditions in Space and time Reading Assignment: F&S Chapter 4 Molecular Dynamics What to choose in an integrator The Verlet algorithm Boundary Conditions in Space and time Reading Assignment: F&S Chapter 4 MSE485/PHY466/CSE485 1 The Molecular Dynamics (MD) method

More information

Molecular Dynamics 9/6/16

Molecular Dynamics 9/6/16 Molecular Dynamics What to choose in an integrator The Verlet algorithm Boundary Conditions in Space and time Reading Assignment: Lesar Chpt 6, F&S Chpt 4 1 The Molecular Dynamics (MD) method for classical

More information

MOLECULAR DYNAMICS SIMULATIONS OF CH4 CLATHRATE HYDRATE DISSOCIATION ADJACENT TO HYDRATED SILICA SURFACES

MOLECULAR DYNAMICS SIMULATIONS OF CH4 CLATHRATE HYDRATE DISSOCIATION ADJACENT TO HYDRATED SILICA SURFACES Proceedings of the 7th International Conference on Gas Hydrates (ICGH 2011), Edinburgh, Scotland, United Kingdom, July 17-21, 2011. MOLECULAR DYNAMICS SIMULATIONS OF CH4 CLATHRATE HYDRATE DISSOCIATION

More information

A Framework for Hybrid Parallel Flow Simulations with a Trillion Cells in Complex Geometries

A Framework for Hybrid Parallel Flow Simulations with a Trillion Cells in Complex Geometries A Framework for Hybrid Parallel Flow Simulations with a Trillion Cells in Complex Geometries SC13, November 21 st 2013 Christian Godenschwager, Florian Schornbaum, Martin Bauer, Harald Köstler, Ulrich

More information

Monte Carlo strategies for first-principles simulations of elemental systems

Monte Carlo strategies for first-principles simulations of elemental systems Monte Carlo strategies for first-principles simulations of elemental systems Lev Gelb Department of Materials Science and Engineering, University of Texas at Dallas XSEDE12 Lev Gelb (UT Dallas) Monte Carlo

More information

Supporting Information for. Hydrogen Bonding Structure at Zwitterionic. Lipid/Water Interface

Supporting Information for. Hydrogen Bonding Structure at Zwitterionic. Lipid/Water Interface Supporting Information for Hydrogen Bonding Structure at Zwitterionic Lipid/Water Interface Tatsuya Ishiyama,, Daichi Terada, and Akihiro Morita,, Department of Applied Chemistry, Graduate School of Science

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