24小时热门版块排行榜

返回列表

bigsmallg

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 111
红花: 1
帖子: 36
在线: 21.6小时
虫号: 4198276
注册: 2015-11-05
专业: 力学

[求助] 关于lammps中老外的compte命令

关于lammps的一个问题，求助大神们。
最近在做关于晶界能的计算，并查阅了一下老外的做法，其中有关于compute这个命令的写法。
老外的：
compute peratom all pe/atom
compute eatoms all reduce sum c_peratom
我的：
compute 1 all pe

我就想问，他那个eatoms是为嘛要这样搞，直接像我这样不行吗？并不是一个老外这样做，好几个人的命令都是这样子。而且，他最后thermo了pe，pe和eatoms算的是一样的。
结果如下：
Step PotEng eatoms
   0 -14908.862 -14908.862
   1 -14909.227 -14909.227
   2 -14909.345 -14909.345
   3 -14909.396 -14909.396
   4 -14909.461 -14909.461
   5 -14909.475 -14909.475
   6 -14909.482 -14909.482
   7 -14909.488 -14909.488
   8    -14909.5    -14909.5
   9 -14909.507 -14909.507
   10 -14909.517 -14909.517
   11 -14909.525 -14909.525
   12 -14909.53 -14909.53
   13 -14909.536 -14909.536
   14 -14909.539 -14909.539
   15 -14909.54 -14909.54
   16 -14909.54 -14909.54
这只是复制了一点log文件的内容，之后也是同样相等的结果。

而且就算根据定义来说，总能减去单个原子势能，是个什么意思？