| 查看: 802 | 回复: 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 |
回复此楼» 猜你喜欢
反应求助
已经有6人回复
-
基于自然哲学类比的风化壳型稀土矿
已经有16人回复
-
看《给阿ma的情书》有感
已经有10人回复
-
2026年WR青拔进展
已经有7人回复
-
职称论文投稿
已经有11人回复
-
中！中！中！
已经有10人回复
-
会评什么时候开始？
已经有4人回复
-
国自然申请五篇代表作大比拼，感觉这个是最重要的
已经有12人回复
-
无聊看看时间戳打发时间
已经有4人回复
-
咨询
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何在lammps里加入新的势函数形式？
已经有6人回复
- 各位高手，怎么在windows下编译LAMMPS？
已经有4人回复
- lammps的时间步长问题
已经有16人回复
- 薄膜在lammps条件下的谐振
已经有5人回复
- LAMMPS Error:lost atoms
已经有7人回复
- 寻求lammps的应用 欢迎大家畅所欲言
已经有8人回复
- 【求助】LAMMPS、NAMD和CHARMM谁最适合生物大分子模拟？
已经有8人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
代友征（上海）-江浙沪独女生，等一个年下的他
+3/342
-
征一人，此生平安欢喜度日常
+1/69
-
坐标北京，诚征男友，找到会删帖。试试运气！说不定就有了！
+1/43
-
温州医科大学康复医学院招聘博后！！
+1/32
-
天津大学 合成生物与生物制造学院 招27级推免硕士研究生及博士研究生
+1/26
-
沉默
+2/20
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/10
-
北理工集成电路杰青团队 | 诚招科助理
+1/10
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/10
-
送什么好呢
+1/9
-
东北师范大学荒漠与草地生态学方向诚聘博士后青年才俊
+1/8
-
国际普刊/正规ISSN/DOI/一周见刊
+1/7
-
西湖大学工学院薄膜电子实验室科研助理及博后招聘
+1/6
-
墨尔本大学(26年QS19)招全奖博士/CSC博士(补齐全奖)等-材料/生物电子/器官芯片等
+1/6
-
中国科学院生态环境研究中心环境工程与公卫方向招（联合培养）博士和硕士研究生
+1/5
-
密苏里大学生物材料合成生物学博士后招聘
+1/4
-
广东医科大学点击化学药物团队招有机、药化、药学等方向博士后（40w 或 55w年薪）
+1/2
-
厦门大学航空航天学院智能制造课题组招2026年硕士、博士生若干名
+1/2
-
北京理工大学团簇研究团队招聘博士后 （化学+材料）
+1/1
-
【急招】2026秋季入学全额奖学金博士招生，北师香港浸会大学，大湾区大学合作培养
+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