| 查看: 872 | 回复: 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 |
回复此楼» 猜你喜欢
面上合作单位盖章
已经有5人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有8人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
青基如何破局
已经有7人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何在lammps里加入新的势函数形式？
已经有6人回复
- 各位高手，怎么在windows下编译LAMMPS？
已经有4人回复
- lammps的时间步长问题
已经有16人回复
- 薄膜在lammps条件下的谐振
已经有5人回复
- LAMMPS Error:lost atoms
已经有7人回复
- 寻求lammps的应用 欢迎大家畅所欲言
已经有8人回复
- 【求助】LAMMPS、NAMD和CHARMM谁最适合生物大分子模拟？
已经有8人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
三无人员面上一把中散金币
+2/1898
-
发点金币玩
+5/560
-
面上得中
+1/465
-
大连理工大学王铁胜副教授招收2027年入学博士生
+1/191
-
中国基，散金币！！！！！！
+5/144
-
北京化工大学何静教授团队诚聘青年英才和博士后
+1/93
-
东华大学材料科学与工程学院丁奇老师招收2026级硕士研究生
+1/88
-
欢迎报考中科院苏州纳米所的氢能、燃料电池方向的博士生！（长期有效）
+1/55
-
欢迎保送和统考中科院苏州纳米所的氢能、燃料电池方向的研究生！（长期有效）
+1/54
-
中国科学院苏州纳米技术与纳米仿生研究所 仿生自适应智能系统实验室招收联培生
+1/39
-
现在生物医药相关有什么兼职可以做吗？
+1/35
-
国家蛋白质科学中心（北京），招收分析化学、药学和蛋白质组代谢组学博士研究生
+1/17
-
国家蛋白质科学中心（北京），招收分析化学、药学和蛋白质组代谢组学博士研究生
+1/16
-
2027材料申博自荐C9本国科大硕
+1/15
-
新审查指南堵死了"一案双申"的老路，实用新型和发明只能二选一了？
+1/10
-
浙江大学医学院黄浩杰教授团队招聘4-5名博士后
+1/6
-
项目助理招聘启事
+1/5
-
南京大学 智能驱动与感知材料实验室 诚招申请考核博士生/科研助理/博士后
+1/4
-
南方医科大学深圳医院陈烨/余涛教授课题组联合培养博士后招聘
+1/2
-
细胞治疗全流程关键检测技术与研究策略
+1/2
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
投票:
