24小时热门版块排行榜    

CyRhmU.jpeg
查看: 4481  |  回复: 22
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ysshj

金虫 (小有名气)

[交流] 求助:lammps非周期边界已有6人参与

向大家请教一下,lammps里的非周期边界f,s,m是什么意思?物理意义是什么?手册上只是形象的描述,不算解释吧。
lammps里其他的命令是怎么利用非周期边界条件的?只是用于原子丢失的判定吗?
假如我设置一个很大的盒子,原子只占里面的很小部分,肯定不会丢失,这时f,s,m下是一样的结果吗?
回复此楼
一个人要像一支队伍。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ysshj

金虫 (小有名气)

★ ★
月只蓝: 金币+2, 鼓励有价值的讨论 2014-02-15 16:52:26
置顶。。。
我来说明一下。
1、f,s边界下盒子实际尺寸(例如我的例子中的zlo,zhi)是不同的。
2、原子编号由data文件确定。这点请大家指正。
3、不同盒子尺寸下,原子信息存储顺序不一样。这点请大家指正。
可能的证据:
手册:For example, if your atoms extend from 0 to 50,you should not specify the box bounds as -10000 and 10000. This is because LAMMPS uses the specified box size to layout the 3d grid of processors.
试算:dump 123 all custom 1 vx.txt id type vx,看下顺序。
4、猜想:langevin温控中产生的随机力序列是按存储的原子顺序,而不是原子编号施加。这点请大家指正。cpp文件我没太看明白。
于是得到进一步猜想:
不同边界条件或盒子尺寸下,langevin温控的随机力施加情况是不一样的,因此计算结果不一样,这体现了随机性。
我的尝试:
我把fix_langevin.cpp产生随机数的部分注释掉后编译计算程序,不同边界条件或盒子尺寸下计算结果是一样的。NVT系综(采用默认Nose-Hoover温控)下计算结果也是一样的。
请大家指点。
一个人要像一支队伍。
19楼2014-02-15 16:48:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答

ysshj

金虫 (小有名气)

同行们发表下看法吧,否则我这问题显得特没趣,甚至lammps的非周期边界特没趣。
一个人要像一支队伍。
2楼2013-12-26 10:04:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dyc_2008

捐助贵宾 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
2楼: Originally posted by ysshj at 2013-12-26 10:04:39
同行们发表下看法吧,否则我这问题显得特没趣,甚至lammps的非周期边界特没趣。

真的很没趣  测测就知道了
一布福田,二修菩提
3楼2013-12-27 12:48:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ysshj

金虫 (小有名气)


月只蓝: 金币+1, 鼓励交流! 2013-12-31 16:21:15
引用回帖:
3楼: Originally posted by dyc_2008 at 2013-12-27 12:48:56
真的很没趣  测测就知道了...

终于有人回复了,感激涕零。
假如我设置一个很大的盒子,原子只占里面的很小部分,肯定不会丢失,这时f,s,m下是不一样的结果。我给Lammps/Mailing Lists发过信息,得到的回复意思应该是结果会一样
大家如果用到了非周期边界,应该都考虑过这三种的差别吧。我不想给大家什么暗示,只是想知道下大家的结果如何,最后怎么定的边界。
一个人要像一支队伍。
4楼2013-12-27 15:46:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见