jianmu
/
lammps
forked from lijiext/lammps
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Remove bomh bond-order-scaling coefficient 

git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@2638 f3b2605a-c512-4ea7-a41b-209d697bcdaa
This commit is contained in:
athomps 2009-03-12 22:57:44 +00:00
parent b960f4a534
commit 5086d7f399
1 changed files with 6 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
lib/reax/reax_connect.F
View File

@ -311,21 +311,20 @@ c$$$ call mdsav(1,qfile(nprob))
ehulp=zero
ehulpp=zero
ehulppp=zero
bomh=bom(itype)
if (rapt(ih1).gt.zero.and.rapt(ih2).gt.zero) then
rhulp2=dis/rob2(ih1,ih2)
rh2p=rhulp2**ptp(itype)
ehulpp=bomh*exp(pdp(itype)*rh2p)
ehulpp=exp(pdp(itype)*rh2p)
end if
if (vnq(ih1).gt.zero.and.vnq(ih2).gt.zero) then
rhulp3=dis/rob3(ih1,ih2)
rh2pp=rhulp3**popi(itype)
ehulppp=bomh*exp(pdo(itype)*rh2pp)
ehulppp=exp(pdo(itype)*rh2pp)
end if
if (rat(ih1).gt.zero.and.rat(ih2).gt.zero) then
rh2=rhulp**bop2(itype)
ehulp=bomh*(1.0+cutoff)*exp(bop1(itype)*rh2)
ehulp=(1.0+cutoff)*exp(bop1(itype)*rh2)
end if
bor=ehulp+ehulpp+ehulppp
@ -724,21 +723,20 @@ c$$$ call mdsav(1,qfile(nprob))
ehulp=zero
ehulpp=zero
ehulppp=zero
bomh=bom(itype)
if (rapt(ih1).gt.zero.and.rapt(ih2).gt.zero) then
rhulp2=dis/rob2(ih1,ih2)
rh2p=rhulp2**ptp(itype)
ehulpp=bomh*exp(pdp(itype)*rh2p)
ehulpp=exp(pdp(itype)*rh2p)
end if
if (vnq(ih1).gt.zero.and.vnq(ih2).gt.zero) then
rhulp3=dis/rob3(ih1,ih2)
rh2pp=rhulp3**popi(itype)
ehulppp=bomh*exp(pdo(itype)*rh2pp)
ehulppp=exp(pdo(itype)*rh2pp)
end if
if (rat(ih1).gt.zero.and.rat(ih2).gt.zero) then
rh2=rhulp**bop2(itype)
ehulp=bomh*(1.0+cutoff)*exp(bop1(itype)*rh2)
ehulp=(1.0+cutoff)*exp(bop1(itype)*rh2)
end if
bor=ehulp+ehulpp+ehulppp