| 查看: 720 | 回复: 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 |
回复此楼» 猜你喜欢
所感
已经有3人回复
-
要不要辞职读博？
已经有7人回复
-
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何在lammps里加入新的势函数形式？
已经有6人回复
- 各位高手，怎么在windows下编译LAMMPS？
已经有4人回复
- lammps的时间步长问题
已经有16人回复
- 薄膜在lammps条件下的谐振
已经有5人回复
- LAMMPS Error:lost atoms
已经有7人回复
- 寻求lammps的应用 欢迎大家畅所欲言
已经有8人回复
- 【求助】LAMMPS、NAMD和CHARMM谁最适合生物大分子模拟？
已经有8人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
求一个访问学者邀请函，非常非常感谢
+1/677
-
北京理工大学郑长松教授课题组诚招2026年秋季博士/硕士研究生
+3/353
-
上海科技大学物质科学与技术学院｜王平鸾课题组长期招聘（博后/博硕/科研助理）
+1/188
-
双一流南京医科大学招计算机、AI、统计、生物信息等方向26年9月入学博士
+1/177
-
原子层沉积（ALD）磁控溅射PECVD等微纳代工服务：18817872921
+1/87
-
内蒙古大学能源材料化学研究院招聘2026年博士生
+1/73
-
考核制博士自荐
+1/37
-
中科院理化技术研究所张飞龙研究员/王树涛研究员团队招生（博士/硕士）
+1/31
-
中国地质大学（北京）王琳课题组招收2026年硕转博/申请-考核博士研究生-1月8日截止
+1/28
-
SCI，计算机相关可以写
+1/23
-
征女友 @长安
+1/21
-
SCI，计算机相关可以写
+1/21
-
SCI，计算机相关可以写
+1/21
-
SCI计算机相关论文
+1/11
-
长江学者团队招聘药学/生物信息学等方向高校教师7名(地点杭州、有事业编)+博后5名
+1/8
-
重庆医科大学-药学院-新靶标教育部医药基础研究创新中心-药物化学2026年博士招生
+1/8
-
上海理工顾敏院士/李蔚团队招收2026级博士研究生 （集成光学、量子信息方向）
+1/3
-
海南大学化学院—功能分子器件团队博士后招聘
+1/1
-
香港中文大学（深圳）下一代半导体激光器概念验证中心（筹）招聘工程师
+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