24小时热门版块排行榜

返回列表

leila3409

捐助贵宾 (正式写手)

应助: 0 (幼儿园)
金币: 616.2
帖子: 310
在线: 71.5小时
虫号: 985885

[交流] 【求助】大家帮我看看这多出来的10个原子是从哪里来的？

我用LAMMPS建了一个SIGMA5晶界，用了下面的语句
lattice bcc ${latparam}
region whole block 0.000000 9.02928935 -9.02928935  9.02928935 0.000000 28.55312  units box
create_box 2 whole
region upper block INF INF 0.000000 9.02928935 INF INF units box
lattice bcc ${latparam} orient x  3 1 0 orient y  1 -3 0 orient z  0 0 -1
create_atoms 1 region upper
region lower block INF INF -9.02928935 0.000000 INF INF units box
lattice bcc ${latparam} orient x 3 1 0 orient y  -1 3 0  orient z  0 0 1
create_atoms 1 region lower
group upper region upper
group lower region lower
delete_atoms overlap ${overlapdist} lower upper

话说我看了图，好像已经成功了。问题我看不懂计算过程中的一些结果，请高手指点。

结果中显示
create_atoms 1 region upper
Created 200 atoms
create_atoms 1 region lower
Created 200 atoms

group upper region upper
210 atoms in group upper
group lower region lower
210 atoms in group lower

delete_atoms overlap 0.275 lower upper
System init for delete_atoms ...
Deleted 20 atoms, new total = 380

我想请问，我明明用upper的region做成一个group，这多出来的10个原子是哪里来的？谢谢!

回复此楼

» 猜你喜欢

自荐读博已经有5人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
上海工程技术大学【激光智能制造】课题组招收硕士已经有5人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
需要合成515-64-0，50g，能接单的留言已经有4人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
带资进组求博导收留已经有10人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复

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

大家帮我看看哪个单位更有发展前途已经有13人回复
大家帮我看看这个东西叫什么名字已经有5人回复
大家帮我看看这两个吸附脱附曲线主要是什么类型的孔已经有10人回复
同一个平板上挑的两个菌测序，序列差很多，大家帮我看看啊，有比对图！谢谢啦已经有18人回复
大家帮我看看这个图是啥意思谢谢啊已经有4人回复
大家来帮我看看这两个INCAR 文件有什么区别哦。。为啥计算结果不一样呢已经有4人回复
马上要选导师了，大家帮我看看选哪个方向的比较好，懂物理的进。。。已经有15人回复
求求大家帮我看看这两个物质。能详细的帮我说说嘛。已经有9人回复
请大侠帮看一下74个原子ENCUT测试应该选哪个值？已经有9人回复
新手，请高手帮忙，计算小分子中某个碳原子上的电子云密度已经有15人回复
请大家帮我看看这段回复是什么意思？已经有5人回复
大家快来看看啊，Ishii-Venturello 催化剂中有四个氧原子怎么是三价的？已经有27人回复
大家帮我看看这个图是怎么画出来的已经有12人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-03-22 16:12:28

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

amynihao

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 361.1
帖子: 443
在线: 230.7小时
虫号: 941269

我看不懂，请高手来回答吧，帮顶。

赞一下

回复此楼

2楼2011-03-23 20:40:20

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

御剑江湖

荣誉版主 (著名写手)

★ ★ ★
ghcacj(金币+3): 谢谢 2011-03-25 13:27:15
leila3409(金币+1): 5 2011-03-28 07:56:00
leila3409(金币+4): 明明回复5个，可能写错地方了，补上 2011-03-28 13:08:08

楼主没有贴出整个in文件，只能猜想：
1、初始结构没有定义好没？
2. 系综选择是否正确？
3、固定边界没，周期性边界条件。
4、势函数的选择。
可以用可视化软件看轨迹，比如VMD，atomeye等，然后可以直观地知道到底发生了什么。

赞一下(1人)

回复此楼

3楼2011-03-24 15:51:24

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

leila3409

捐助贵宾 (正式写手)

应助: 0 (幼儿园)
金币: 616.2
帖子: 310
在线: 71.5小时
虫号: 985885

★ ★
ghcacj(金币+2): 谢谢分享解决问题的方法 2011-03-28 12:10:35

引用回帖:

Originally posted by 御剑江湖 at 2011-03-24 15:51:24:
楼主没有贴出整个in文件，只能猜想：
1、初始结构没有定义好没？
2. 系综选择是否正确？
3、固定边界没，周期性边界条件。
4、势函数的选择。
可以用可视化软件看轨迹，比如VMD，atomeye等，然后可以直观地知 ...

呵呵，谢谢，这个问题已经解决了。
create_atoms 在创建晶胞时会考虑周期性。而group upper region upper不考虑周期性，所以会多出十个原子来。

赞一下(1人)

回复此楼

4楼2011-03-28 07:55:50

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