| 查看: 1951 | 回复: 0 | ||
Bessel木虫 (正式写手)
|
[求助]
请问fix nvt 和 fix move 不能同时用吗？
|
|
初学LAMMPS，小弟有个问题一直搞不定。 我要模拟的系统的初始状态如图一所示，两组分混合物，下面绿色的原子固定不动，做为边界墙， 我用的命令是 fix 2 lo move linear 0 0 0 上层红色的原子做NVT系综， 命令是 fix 1 hi nvt temp 3 3 200 结果模拟一开始，系统就变成如图二所示，红色的球全没有了。 我查看了dump出来的文档，红球的坐标全变成了 -nan -nan -nan 但是系统没有报错，一直运行的循环结束。 还有我尝试把绿色的原子的Z方向运动不做限制（这才是我真实想做的系统），命令如下 fix 2 lo move linear 0 0 NULL 结果也是红球全不见了。 完整的in文件如下 units lj atom_style atomic lattice fcc 0.8442 region box block 0 5 0 5 0 5 create_box 2 box region hi block 0 5 0 5 2.5 5 region lo block 0 5 0 5 0 2.5 create_atoms 1 region hi create_atoms 2 region lo mass * 1.0 velocity all create 3.0 825577 dist gaussian pair_style lj/cut 2.5 pair_coeff * * 1.0 1.0 2.5 neighbor 2.0 nsq neigh_modify every 20 delay 0 check no group hi type 1 group lo type 2 fix 1 hi nvt temp 3 3 200 fix 2 lo move linear 0 0 0 #fix 2 lo setforce 0 0 0 dump 1 all atom 100 dump.friction dump_modify 1 sort id dump 2 all image 200 *.jpg type type timestep 0.005 thermo 1000 run 5000 clear 请前辈指点一下，两组分的混合物，如何fix呢。 我尝试把红球的命令变成 fix 1 hi nve, 红球没跑丢，但是这不是我想要的系综。 谢谢！  图一  图二 |
回复此楼» 猜你喜欢
面上项目没有好文章就没希望了吗？
已经有15人回复
-
困死了
已经有8人回复
-
不知道还有没有招博士的学校了
已经有5人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有8人回复
-
材料博士申请
已经有5人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有9人回复
-
还有课题组有博士名额吗
已经有6人回复
-
关于水星近日点进动成因的质疑 与实证分析
已经有10人回复
-
博士申请
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
L-丙氨酸甲酯盐酸盐是不是在空气中不能稳定存在吗？那么容易吸水吗？
已经有4人回复
- 【共享】Infix PDF Editor 简体中文版 – 像编辑普通文本文件一样编辑PDF
已经有301人回复
- XRD能用于算晶体尺寸吗？咋有同学说算的是层间距？能不能给个算尺寸的具体例子
已经有5人回复
- fix deform 指令拉伸求教
已经有10人回复
- 关于fix nvt里的Tdamp
已经有7人回复
- 中药复方水提液十分粘稠，又不能稀释，请问如何在不稀释的情况下降低黏性
已经有34人回复
- 求解读,请问今年是不是非211,非985的高校的博一学生不能申请联合培养了
已经有12人回复
- 【求助完结】请问442和419有啥区别啊？419没有博士学位能申请吗？
已经有4人回复
- 请问Endnote里面online search中ebsco不能用是什么原因？
已经有9人回复
- 【求助】请问氧原子有两对未共用电子对，可以同时作为两个氢键的受体吗
已经有7人回复
- 【求助】求助fix ave/correlate的源文件
已经有6人回复
- 【求助】请问《环境污染治理技术与设备》和《环境工程学报》是同一个期刊吗？谢谢
已经有33人回复
- 【讨论】请问同一个学校不同专业考试科目一样过线不能调剂吗？
已经有21人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人