| 查看: 5104 | 回复: 9 | ||||
[求助]
关于模拟离子液体力场的选择（opls-aa)
|
|
各位虫子们，大家模拟离子液体一般用什么力场呀？ 我选择用opls_aa力场，因为这个力场中有专门的离子液体的力场。我用的lammps软件。在模拟中，问题好多，跑着跑着，原子就会missing。我的初始位置，分子间隔的挺远，大约10埃，我dump force后发现，力都好大，不清楚啥个原因。 我想换个力场试试，不知哪个力场好点。 谢谢大家了。 |
回复此楼» 收录本帖的淘帖专辑推荐
 分子模拟要术 | 分子动力学-初级问题资源 |
» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
国自然面上复盘~欢迎讨论
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
基金不中，共勉
已经有5人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有3人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
系统查不到
已经有10人回复
-
怎么查啊
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教关于materials studio建模和OPLS力场
已经有7人回复
- 【求助】请问如何下载OPLS-AA力场
已经有7人回复
043114076
木虫 (正式写手)
- 应助: 49 (小学生)
- 金币: 3364.9
- 散金: 52
- 红花: 23
- 帖子: 788
- 在线: 369.5小时
- 虫号: 570592
- 注册: 2008-06-07
- 性别: GG
- 专业: 理论和计算化学
2楼2013-01-27 18:56:18
cfmzxf84
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 5003.4
- 散金: 822
- 红花: 1
- 帖子: 908
- 在线: 328.6小时
- 虫号: 814536
- 注册: 2009-07-23
- 专业: 化学环境污染与健康
3楼2013-01-27 20:44:33
|
是用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
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
6楼2013-01-28 22:32:27
Bessel
木虫 (正式写手)
- 模拟EPI: 2
- 应助: 47 (小学生)
- 金币: 1506.5
- 散金: 4755
- 红花: 8
- 帖子: 609
- 在线: 234.1小时
- 虫号: 370174
- 注册: 2007-05-13
- 专业: 凝聚态物性I:结构、力学和
7楼2013-01-29 02:58:58
|
我在nve系综下跑的minimize。我的参数是 minimize 1.0e-4 1.0e-6 100 1000 min_style hftn min_modify dmax 0.2 我的结果是：（我从log里面拷出来的） minimize 1.0e-4 1.0e-6 100 1000 WARNING: Resetting reneighboring criteria during minimization (min.cpp:167) Ewald initialization ... G vector = 0.227629 vectors: actual 1d max = 2748 10 4630 Memory usage per processor = 24.4557 Mbytes Step Temp E_pair E_mol TotEng Press 0 0 161.97025 59.620071 221.59032 -139.69835 52 0 139.99619 32.197295 172.19348 -13.155359 Loop time of 0.369925 on 48 procs for 52 steps with 72 atoms Minimization stats: Stopping criterion = energy tolerance Energy initial, next-to-last, final = 221.590316114 172.205685403 172.193484683 Force two-norm initial, final = 334.999 8.36533 Force max component initial, final = 89.5407 2.13135 Final line search alpha, max atom move = 0.000551252 0.00117491 Iterations, force evaluations = 52 255 Pair time (%) = 0.000742118 (0.200613) Bond time (%) = 0.000933076 (0.252234) Kspce time (%) = 0.182548 (49.3473) Neigh time (%) = 0 (0) Comm time (%) = 0.176163 (47.6214) Outpt time (%) = 0 (0) Other time (%) = 0.00953816 (2.57841) Nlocal: 1.5 ave 21 max 0 min Histogram: 44 0 0 0 0 0 1 1 0 2 Nghost: 38.5625 ave 72 max 0 min Histogram: 12 0 0 3 6 3 7 2 1 14 Neighs: 27.2708 ave 471 max 0 min Histogram: 44 0 0 0 2 0 0 0 1 1 Total # of neighbors = 1309 Ave neighs/atom = 18.1806 Ave special neighs/atom = 10.2222 Neighbor list builds = 0 Dangerous builds = 0 min_style hftn min_modify dmax 0.2 真是太谢谢你了。。。 今天我又这样跑了跑：就是在nvt系综下，timestep从0.1开始，每100000步增加0.1，结果到了1后竟然没有错误提示了，但是每个原子的受力力还在几十多左右。 |
8楼2013-01-29 09:46:56
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
9楼2013-02-08 13:03:14
10楼2013-02-15 10:53:51
投票:
