| 查看: 1505 | 回复: 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呢？ |
回复此楼» 猜你喜欢
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用卡尔费休库仑法测定易结晶的化合物
已经有11人回复
- 分子动力学模拟长程库仑势的处理问题
已经有8人回复
- 求助关于控制电位电解库仑法的参数设置
已经有5人回复
cmkk
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1616.8
- 散金: 50
- 红花: 1
- 帖子: 25
- 在线: 21.4小时
- 虫号: 721639
- 注册: 2009-03-13
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
2楼2012-05-31 17:56:37
100