24小时热门版块排行榜    

查看: 1422  |  回复: 1

cmkk

木虫 (初入文坛)

[求助] lammps用coul/cut和coul/long 算库仑势能不一致的问题

units                  metal #eV A ps
atom_style             charge
boundary               p p p
newton         on

lattice sc             10
region         mybox   block 0 10 0 10 0 10 units lattice
create_box             1 mybox
create_atoms 1         region mybox
mass 1 12
group 1 type 1
set group 1 charge  0.5

用lammps产生的1000个原子的simple cubic网格,所有原子都设置的0.5e的电量,可是当我用pair_style coul/long 来算时,得到的PotEng 是 -510.71783eV, 而用coul/cut得到的 1079.9734eV.
这俩方法得到的势能确实会有不同,可是为什么一个正一个负呢??
全部原子电量都相同的话,为什么会出现负值的PotEng呢?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cmkk

木虫 (初入文坛)

没人搭理我吗?
2楼2012-05-31 17:56:37
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cmkk 的主题更新
信息提示
请填处理意见