Use Au-Si MEAM Potential in LAMMPS

Use Au-Si MEAM Potential in LAMMPS

Yanming wang, Seunghwa Ryu and Wei Cai

June, 2012

This tutorial explains how to use Au-Si MEAM potential in LAMMPS on a cluster like mc2. It starts with the installation of LAMMPS on mc2, then talks about the compiling of the Au-Si MEAM potential. This tutorial also includes a test case to make sure the result from LAMMPS with this MEAM potential is consistent with that from MD++.

---

Installing LAMMPS on mc2

LAMMPS is a classical molecular dynamics code for large-scale atomic/molecular massively parallel simulation. On mc2, you can get LAMMPS from our SVN server using the following commands.

mkdir -p ~/Codes/LAMMPS.svn
cd ~/Codes/LAMMPS.svn
svn co https://micro.stanford.edu/svn/LAMMPS .

Typing the following SVN command from within the "LAMMPS.svn" directory, you can update your LAMMPS package.

svn update

Compiling MEAM potential in LAMMPS

Preparing for parallel run

First, you need to download the fftw.tar.gz and upload this archive file to your account on mc2. In your home directory, make a new folder named usr, then type the following command to unzip the fftw library in your usr folder.

cd ~usr/
tar -zxvf fftw.tar.gz

ALso add the following lines at the end of your ~/.bash_profile file to clarify the compiler.

module load intel/intel-11
module load mvapich/1.2rc1-intel-11-dell-gen2

Preparing for compiling MEAM potential

Au-Si MEAM potential is a potential developed by Dr. Seunghwa Ryu and Prof. Wei Caifor gold-silicon binary system. To compile this potential in LAMMPS, you need to open the Makefile file under the LAMMPS.svn directory, finding the following line

cd lib/meam;  make -f Makefile.$(libFC.$(SYS)) F90FLAGS="-O -D_EROSE_FUNCTION=1"

and replace it with

cd lib/meam;  make -f Makefile.$(libFC.$(SYS)) F90FLAGS="-O -D_EROSE_FUNCTION=3"

Next, in the ~/.bash_profile file, add the following line to specify the system you are using.

export TARGET=mc2

Finally, go to LAMMPS.svn directory and type

make

This will make lmp_mc2 in LAMMPS.svn/src/ directory.

Test case:Comparsion with MD++

This section illustrates how to compare the simulation result from LAMMPS using the MEAM potential with that from MD++. A successful test will show MD++ and LAMMPS report the same potential energy and forces from a given atoms' configuration.

MD++ test case

Supposing you've known how to run MD++ code from our manuals and downloaded the MD++ package from our SVN server, you can compile the MEAM potential typing the command,

make meam-lammps build=R SYS=mc2