| 查看: 845 | 回复: 2 | |||
[交流]
【求助】Gaussian core model 的LAMMPS模拟
|
|
大家好 我用LAMMPS模拟Gaussian core model ，发现模拟后粒子都集中在一个非常小的区域，都重叠在一起了。不知道是为什么。请各位指点一下 下面是我的LAMMPS模拟的in 文件 #GCM AT 0.04 units lj boundary p p p atom_style atomic lattice fcc 3.419952 region box block 0 8.0 0 8.0 0 8.0 create_box 1 box create_atoms 1 box pair_style gauss 3.2 pair_coeff * * 1.0 0.9 mass 1 1.0 timestep 0.05 thermo 100 neighbor 0.5 bin neigh_modify every 4 delay 0 check yes velocity all create 0.04 4928459 rot yes dist gaussian fix 1 all nvt temp 0.04 0.04 1.5 drag 0.3 thermo_style custom step temp press vol etotal log log_1.gcm run 500000 dump 1 all custom 100 dump_1.gcm id x y z #fix 3 all msd 100 msd_1.gcm run 10000 unfix 1 clear next x jump in.eam |
回复此楼» 猜你喜欢
filecode
已经有8人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有4人回复
-
filecode
已经有5人回复
-
求各位大神看下
已经有31人回复
-
咨询面上基金
已经有7人回复
-
影响面上的因素
已经有11人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有11人回复
-
一个有机合成实验室都需要哪些设备？
已经有11人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有30人回复
-
好消息？这个有何含义？？？
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何在lammps里加入新的势函数形式？
已经有6人回复
- 各位高手，怎么在windows下编译LAMMPS？
已经有4人回复
- lammps的时间步长问题
已经有16人回复
- 薄膜在lammps条件下的谐振
已经有5人回复
- LAMMPS Error:lost atoms
已经有7人回复
- 寻求lammps的应用 欢迎大家畅所欲言
已经有8人回复
- 【求助】LAMMPS、NAMD和CHARMM谁最适合生物大分子模拟？
已经有8人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
祈福
+5/660
-
一边期待，一边忐忑：发榜前的科研人正在想什么？
+5/230
-
1,3,2-Dioxaphospholane, 2-(2,2,2-trifluoroethoxy)-
+1/91
-
三(三甲基硅基)硼酸酯（CAS：4325-85-3）
+1/90
-
1,3,2-Dioxaphospholane, 2-(2,2,3,3,3-pentafluoropropoxy)-
+1/90
-
乙二醇二甲基环硅醚（CAS：57915-66-9）
+1/90
-
三(三甲基硅基)甘油醚（CAS：6787-10-6 ）
+1/89
-
香港城市大学郑星课题组诚聘博士后 (生态学、环境科学、植物学、大气科学背景)
+1/85
-
上海交通大学诚聘有机化学方向博士生/科研助理
+1/33
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
+2/14
-
希望大家都高中！！
+1/10
-
中山大学可穿戴电路系统方向招2027届考核制博士生
+1/9
-
有合作举办scopus会议的吗？
+1/8
-
瑞典林雪平大学 高峰课题组 诚聘 氮化物基无机材料方向 博士后1名
+1/7
-
高弯曲强度玻璃制品的特殊加工工艺
+1/4
-
北京师范大学（珠海校区）赵亚然副教授招收2027年电催化硕士/博士研究生
+1/3
-
2026年University of South Carolina继续招收博士生
+1/2
-
【香港城市大学】范战西教授课题组现招收博士生（全奖, 2027年秋季入学）
+1/2
-
致每一位等待国自然结果的科研人：我们想听听你的心声
+1/1
-
国自然发榜在即：此刻的你，是什么心情？
+1/1
2楼2010-11-26 12:54:25
★ ★ ★ ★
cxz314(金币+5):非常感谢！我在核实一下 2010-11-27 19:57:37
ghcacj(金币+4):谢谢 2010-11-29 18:23:05
cxz314(金币+5):非常感谢！我在核实一下 2010-11-27 19:57:37
ghcacj(金币+4):谢谢 2010-11-29 18:23:05
|
cxz314(金币+1):谢谢！是这样的Gaussian core model 的约化单位和LJ势的约化单位一样的 这个。。。我有点不明白，pair_style gauss的势函数形式是 E=Aexp(-Br^2) ，这样的势，要用约化单位么？那你的sigma，epsilon之类是如何设置的？都是1么？那你的 lattice fcc 3.419952，还有region, 都是约化单位么？建议你再核实一下。不过我不是很了解这个势，也可能没有问题。 |
3楼2010-11-27 17:39:21