Added placeholder for neural network SNAP potential

# DATE: 2017-02-20 CONTRIBUTOR: Mitchell Wood mitwood@sandia.gov CITATION: Wood, M. A. and Thompson, A. P. "Quantum-Accurate Molecular Dynamics Potential for Tungsten" arXiv:1702.07042 [physics.comp-ph]

#

# Definition of SNAP+ZBL potential.

variable zblcutinner equal 4

variable zblcutouter equal 4.8

variable zblz equal 74

# Specify hybrid with SNAP and ZBL

pair_style hybrid/overlay &

zbl ${zblcutinner} ${zblcutouter} &

snap

pair_coeff 1 1 zbl ${zblz} ${zblz}

pair_coeff * * snap W_2940_2017_2.snapcoeff W_2940_2017_2.snapparam W

#Nomenclature on the snap files are Element_DakotaID_Year_Month

# Demonstrate SNAP Ta potential

# Initialize simulation

variable nsteps index 100

variable nrep equal 4

variable a equal 3.1803

units		metal

# generate the box and atom positions using a BCC lattice

variable nx equal ${nrep}

variable ny equal ${nrep}

variable nz equal ${nrep}

boundary	p p p

lattice         bcc $a

region		box block 0 ${nx} 0 ${ny} 0 ${nz}

create_box	1 box

create_atoms	1 box

mass 1 183.84

# choose potential

include W.nnsnap

# Setup output

thermo		10

thermo_modify norm yes

# Set up NVE run

timestep 0.5e-3

neighbor 1.0 bin

neigh_modify once no every 1 delay 0 check yes

# Run MD

velocity all create 300.0 4928459

fix 1 all nve

run             ${nsteps}