| 查看: 770 | 回复: 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 |
回复此楼» 猜你喜欢
反应很差，大量原料没有反应
已经有8人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有13人回复
-
评审有感
已经有23人回复
-
提交了我也来说说感想
已经有6人回复
-
Sci. Bull. 悲剧经验
已经有7人回复
-
今年审到国自然15份，谈谈感受
已经有20人回复
-
国自然评分
已经有4人回复
-
如果工作了想读博，可以边工作边读全日制嘛？
已经有5人回复
-
如何实现卤原子转化
已经有8人回复
-
河北省自然科学基金
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何在lammps里加入新的势函数形式？
已经有6人回复
- 各位高手，怎么在windows下编译LAMMPS？
已经有4人回复
- lammps的时间步长问题
已经有16人回复
- 薄膜在lammps条件下的谐振
已经有5人回复
- LAMMPS Error:lost atoms
已经有7人回复
- 寻求lammps的应用 欢迎大家畅所欲言
已经有8人回复
- 【求助】LAMMPS、NAMD和CHARMM谁最适合生物大分子模拟？
已经有8人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
祈福，面上申请连续三年了，希望今年面上能中！！！
+1/952
-
求助：透明涂层样品-测涂层厚度。厚度范围小于100nm
+1/100
-
钯炭、铂炭、钌炭等加氢催化剂
+1/80
-
空军军医大学第一附属医院药学部科研助理招聘
+1/78
-
替投稿20天的论文祈福
+1/75
-
注册小木虫周年大礼包领取不了
+5/50
-
10 年TOP猎头｜免费岗位推荐 + 简历优化，直击大厂 offer
+1/49
-
河南省医学科学院王宁利院士科研团队2026年博士、博后、硕士招聘
+1/40
-
化学化工学院招聘博士2人-黔南民族师范学院
+2/34
-
华南师范大学化学学院能源转换与存储材料实验室电合成/电催化课题组招收2026级博士生
+1/30
-
上海交大招收 2027级“明志荣誉博士生” （化学、生物方向-津贴约12万/年）
+1/27
-
青岛农业大学国家杰青团队蔬菜病害绿色防控方向招聘博士后
+1/27
-
海南大学生物医学工程学院 （光学诊疗团队）诚邀优秀人才
+1/14
-
山西大学招收2026级申请考核制化学材料类博士研究生一名
+1/7
-
博士后招聘（高薪40万+）
+1/6
-
邀请函 | 华南理工大学集成电路技术前沿论坛（含电子信息）诚邀您参加！
+1/5
-
西门子医疗（深圳）招聘
+1/4
-
中科院力学所CAE科研助理招聘
+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
