| 查看: 5194 | 回复: 9 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
关于模拟离子液体力场的选择（opls-aa)
|
||||
|
各位虫子们，大家模拟离子液体一般用什么力场呀？ 我选择用opls_aa力场，因为这个力场中有专门的离子液体的力场。我用的lammps软件。在模拟中，问题好多，跑着跑着，原子就会missing。我的初始位置，分子间隔的挺远，大约10埃，我dump force后发现，力都好大，不清楚啥个原因。 我想换个力场试试，不知哪个力场好点。 谢谢大家了。 |
回复此楼» 收录本帖的淘帖专辑推荐
 分子模拟要术 | 分子动力学-初级问题资源 |
» 猜你喜欢
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有8人回复
-
极限拉扯，如何抉择
已经有7人回复
-
课题组会帮你通过考核么？
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教关于materials studio建模和OPLS力场
已经有7人回复
- 【求助】请问如何下载OPLS-AA力场
已经有7人回复
Bessel
木虫 (正式写手)
- 模拟EPI: 2
- 应助: 47 (小学生)
- 金币: 1506.5
- 散金: 4755
- 红花: 8
- 帖子: 609
- 在线: 234.1小时
- 虫号: 370174
- 注册: 2007-05-13
- 专业: 凝聚态物性I:结构、力学和
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
wup123: 金币+5, ★★★很有帮助 2013-01-28 11:02:36
感谢参与，应助指数 +1
wup123: 金币+5, ★★★很有帮助 2013-01-28 11:02:36
|
有原子missing的问题，可能是你的初始构型做的不好，原子相互作用力过大。这和具体是什么力场无关，最好在run命令之前加上minimize命令，可能解决你的问题。 具体的看这里：http://lammps.sandia.gov/doc/minimize.html 还有你的in文件里： fix 2 all langevin 300.0 300.0 0.2 48279 fix 3 all nve （文献上说可以先布朗运动下，所以我也设置了了一下） 更正一下，这个不是真的让体系作布朗运动，而是利用郎之万方程来控温的。但是这个不是问题的关键。 |
5楼2013-01-28 04:00:23
043114076
木虫 (正式写手)
- 应助: 49 (小学生)
- 金币: 3364.9
- 散金: 52
- 红花: 23
- 帖子: 788
- 在线: 369.5小时
- 虫号: 570592
- 注册: 2008-06-07
- 性别: GG
- 专业: 理论和计算化学
2楼2013-01-27 18:56:18
|
是用Pauda他们组的力场吗？二面角只有三个参数v1,v2,v3。我现在是刚开始研究这个，这个问题已经困了我好几天了，一直没找出啥问题来。希望你花一点点时间看看我的输入文件吧。 LAMMPS input file units real boundary p p p atom_style full #uncharged molecules, with bond, angle, dihedral, improper dihedral bond_style harmonic angle_style harmonic dihedral_style opls pair_style lj/cut/coul/long 12 pair_modify mix geometric tail yes kspace_style ewald 1.0e-4 #read_restart restart.mix.* read_data data.il #read_restart restart.mix.* #fix 1 all shake 0.0001 20 10 b 3 neighbor 2.0 bin neigh_modify delay 3 every 1 check yes #delay 3 fix 2 all langevin 300.0 300.0 0.2 48279 fix 3 all nve （文献上说可以先布朗运动下，所以我也设置了了一下） timestep 0.1 dump 1 all custom 20000 dump.atom.*.lammpstrj xu yu zu dump 2 all custom 20000 dump.force.* id mol type element fx fy fz dump_modify 1 flush yes dump_modify 2 flush yes thermo_style custom step temp press etotal vol lx ly lz thermo 1000 thermo_modify flush yes run 1000000 我是做的（mmpy）（tsfi)这个分子，就两对分子。我真的不知道哪里有问题。。。 |
4楼2013-01-28 03:07:34
6楼2013-01-28 22:32:27
