Hi all,
I would like to calculate the energy terms for the airebo-m potential for a subset of atoms. This is what I try in my script:
group mid2800 id 201:3000
compute en mid2800 pair airebo/morse
variable REBO equal c_en[1]
variable LJ equal c_en[2]
variable TORSION equal c_en[3]
compute energ all pair airebo/morse
variable REBO2 equal c_energ[1]
variable LJ2 equal c_energ[2]
variable TORSION2 equal c_energ[3]
There is a total of 3200 atoms in a straight carbon nanotube, and I made this group with the middle 2800. Then when I calculate the energy terms as described in the documentation for airebo, for the two different groups, the energies end up the same:
Step Time Temp c_px PotEng Fmax v_REBO v_LJ v_TORSION v_REBO2 v_LJ2 v_TORSION2 v_REBO2
3 0 10 7.0574824e+08 -23384.286 1.1477917 -24439.676 -36.67048 1092.0596 -24439.676 -36.67048 1092.0596 -24439.676
100 0.0485 10.566012 -6.697875e+08 -23385.109 1.2894551 -24439.851 -38.058452 1092.8004 -24439.851 -38.058452 1092.8004 -24439.851
I attached the files I used to test this. I'm using LAMMPS 16 March 2018. If anyone knots how to calculate the airebo energy terms for a subset of atoms, or what I am missing here, I'd really appreciate the help.
Thanks,
Kyle Rego
3200nanotube.data (201 KB)
knotscript.txt (2.01 KB)