Use MEAM in MD++: Difference between revisions

From Micro and Nano Mechanics Group
Jump to navigation Jump to search
(New page: <P ALIGN="CENTER"> <FONT SIZE="+3" color="darkred"><STRONG> Use MEAM in MD++</STRONG></font></P> <DIV> <P ALIGN="CENTER"><STRONG>Keonwook Kang, Seunghwa Ryu and Wei Cai</STRONG></P> </DIV>...)
 
No edit summary
Line 11: Line 11:


Since then MEAM potentials for many elements, as well as binary and ternary alloy systems have been developed. The MEAM potentials allows us to quickly set up a model for a wide range of materials systems. MEAM has also been implemented in [http://lammps.sandia.gov/ LAMMPS], whose subroutines have been used in the MEAM implementation in MD++.
Since then MEAM potentials for many elements, as well as binary and ternary alloy systems have been developed. The MEAM potentials allows us to quickly set up a model for a wide range of materials systems. MEAM has also been implemented in [http://lammps.sandia.gov/ LAMMPS], whose subroutines have been used in the MEAM implementation in MD++.

Through some examples, we will describe how to use MEAM in MD++ to model a binary system, as well as allowing the simulation cell to change its shape.
<HR>
<HR>


Line 16: Line 18:


First, you need to install MD++ on your computer by following the instructions on [[MD++_Manuals | MD++ Manuals]].
First, you need to install MD++ on your computer by following the instructions on [[MD++_Manuals | MD++ Manuals]].

Second, copy file [[media:meam_test.tcl.txt‎ | meam_test.tcl]] to your input file directory of MD++. You can do so by the following commands (assuming you have installed MD++ in <tt>~/Codes/MD++</tt>).

<pre>
export MDPP=~/Codes/MD++
mkdir -p ${MDPP}/scripts/work/meam
cd ${MDPP}/scripts/work/meam
wget http://micro.stanford.edu/mediawiki-1.11.0/images/Meam_test.tcl.txt -O meam_test.tcl
</pre>


==Compile executable file==

There are several MEAM implementations in MD++,
<pre>
meam-lammps - MEAM (taken from lammps)
meam-baskes - MEAM (taken from Baskes's code dynamo)
meam - MEAM (directly implemented in MD++)
</pre>

<tt>meam-lammps</tt> is the one that is most extensively tested. We will use <tt>meam-lammps</tt> exclusively in the following discussions.

You can compile <tt>meam-lammps</tt> by typing
<pre>
cd ${MDPP}
make meam-lammps build=R SYS=su-ahpcrc
</pre>

(<tt>su-ahpcrc</tt> is one of our computing clusters at Stanford.)

If compilation is successful, this will create binary file <tt>meam-lammps_su-ahpcrc</tt> file in your <tt>${MDPP}/bin</tt> directory.


==Relaxation with fixed box==



==Relaxation at constant stress==



==MD in NVT ensemble==



==MD in N<math>\sigma</math>T ensemble==

Revision as of 17:34, 22 July 2009

Use MEAM in MD++

Keonwook Kang, Seunghwa Ryu and Wei Cai

This tutorial describes how to use the modified embedded-atom method (MEAM) potential in MD++. The MEAM potential was originally developed by Dr. Mike Baskes,

Modified embedded-atom potentials for cubic materials and impurities, Physical Review B, 46, 2727 (1992).

Since then MEAM potentials for many elements, as well as binary and ternary alloy systems have been developed. The MEAM potentials allows us to quickly set up a model for a wide range of materials systems. MEAM has also been implemented in LAMMPS, whose subroutines have been used in the MEAM implementation in MD++.

Through some examples, we will describe how to use MEAM in MD++ to model a binary system, as well as allowing the simulation cell to change its shape.


Download files

First, you need to install MD++ on your computer by following the instructions on MD++ Manuals.

Second, copy file meam_test.tcl to your input file directory of MD++. You can do so by the following commands (assuming you have installed MD++ in ~/Codes/MD++).

 export MDPP=~/Codes/MD++
 mkdir -p ${MDPP}/scripts/work/meam
 cd ${MDPP}/scripts/work/meam
 wget http://micro.stanford.edu/mediawiki-1.11.0/images/Meam_test.tcl.txt -O meam_test.tcl


Compile executable file

There are several MEAM implementations in MD++,

meam-lammps         - MEAM (taken from lammps)
meam-baskes         - MEAM (taken from Baskes's code dynamo)
meam                - MEAM (directly implemented in MD++)

meam-lammps is the one that is most extensively tested. We will use meam-lammps exclusively in the following discussions.

You can compile meam-lammps by typing

 cd ${MDPP}
 make meam-lammps build=R SYS=su-ahpcrc

(su-ahpcrc is one of our computing clusters at Stanford.)

If compilation is successful, this will create binary file meam-lammps_su-ahpcrc file in your ${MDPP}/bin directory.


Relaxation with fixed box

Relaxation at constant stress

MD in NVT ensemble

MD in NT ensemble