| 查看: 1724 | 回复: 8 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
lammps 离子液体，离子中相邻隔的挺远的原子成键了，不知什么原因
|
||
|
我在进行模拟的时候，首先用两个分子（四个离子）测试了，结构非常好，没啥异常。但是增多原子后，发现其中一个阴离子（tfsi-）中F原子和相隔非常远的O原子（相隔五个原子）成键了，从而使整个离子变的乱起八糟。 我检查了力场参数，采用的是opls—aa势，和文献中一致，电荷也守恒。 想问问各位虫子，这个是什么原因呢？该怎么解决呢？ 也不知怎么搞的，我遇到好多诡异的问题。。。 |
回复此楼» 本帖已获得的红花（最新10朵）
凌枫霜叶» 猜你喜欢
唉
已经有8人回复
-
对象没有，还非常想读博，难以抉择
已经有66人回复
-
博后面上如何提高创新分？
已经有8人回复
-
佛山大学和深职大怎么选？
已经有13人回复
-
青拔答辩通知出来了
已经有5人回复
-
2024基金必中
已经有8人回复
-
博士该不该读？
已经有16人回复
-
求助光催化方向二区末，三区期刊
已经有8人回复
-
今天公派访问学者项目会出结果吗？
已经有6人回复
-
我投稿了Elsevier，后来申请撤稿，你不允许我撤稿，强制发表然后问我要钱
已经有12人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- lammps进行分子动力学模拟之前是不是一定要进行结构优化？
已经有4人回复
- 诚聘锂电研发人才！
已经有7人回复
- 招聘锂电研发测试人才
已经有7人回复
- 招聘锂电方面的人才
已经有5人回复
- 国内一流锂电企业招聘研发人才！
已经有18人回复
- ]ATL诚聘各类电池英才
已经有16人回复
- 求助？怎么从剑桥晶体库得到物质的结构数据
已经有5人回复
- 招聘研发人才！
已经有8人回复
- 诚聘锂电研发人才！
已经有7人回复
- 诚聘研发人才！
已经有6人回复
- 聚贤猎头
已经有5人回复
- 电池领域工作机会多个
已经有24人回复
- 锂电池相关招聘启示
已经有8人回复
- 关于模拟离子液体力场的选择（opls-aa)
已经有9人回复
- 博士工作
已经有8人回复
- 2012年3月份招聘黄河科技学院纳米功能材料研究所
已经有27人回复
- 2012黄河科技学院纳米功能材料研究所招聘材料类化学类及材料计算类博士
已经有27人回复
- 2011年9月前黄河科技学院纳米功能材料研究所招聘材料类化学类及材料计算类博士
已经有15人回复
- 2011年7月黄河科技学院纳米功能材料研究所材料类化学类及分子模拟计算类博士招聘
已经有15人回复
- 黄河科技学院纳米材料研究所今年最后一批招聘
已经有8人回复
- 急聘——黄河科技学院纳米功能材料研究所
已经有11人回复
- 2011黄河科技学院纳米功能材料研究所招聘材料类化学类及分子模拟方面博士
已经有51人回复
- 【讨论】lammps data 文件
已经有3人回复
7楼2013-02-08 09:06:00
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
2楼2013-02-04 10:27:01
3楼2013-02-04 11:15:01
|
我模拟的体系是emim阳离子，和fsi阴离子。 我的in文件如下： 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 improper_style cvff pair_style lj/cut/coul/long 12.0 pair_modify mix arithmetic tail yes kspace_style pppm 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 velocity all create 300.0 4567712 rot yes dist gaussian compute 1 all temp #fix 2 all npt temp 300.0 300.0 100.0 iso 1.0 1.0 10.0 #fix 2 all langevin 300.0 300.0 200 48279 #fix 3 all nve fix 2 all nvt temp 300.0 300.0 50.0 minimize 1.0e-4 1.0e-6 100 1000 min_style cg min_modify dmax 0.2 timestep 0.05 dump 1 all custom 20 dump.atom.*.lammpstrj id x y z dump 2 all custom 200000 dump.force.* id mol type element fx fy fz #dump 3 all dcd 2000 dump.trj.lammpstrj dump_modify 1 flush yes dump_modify 2 flush yes #dump_modify 3 flush yes thermo_style custom step c_1 pe etotal vol lx ly lz thermo 1000 thermo_modify flush yes restart 100000 restart.mix.* run 100000 timestep 0.1 #thermo 1000 run 100000 我的log文件 24000 273.6706 272.92199 317.7888 125000 50 50 50 25000 325.54712 263.24033 316.61201 125000 50 50 50 26000 313.35585 255.15618 306.52916 125000 50 50 50 27000 3516.0726 1330.2863 1906.7272 125000 50 50 50 28000 1783.7174 315.87659 608.30734 125000 50 50 50 29000 1473.6364 568.98349 810.57815 125000 50 50 50 30000 1474.526 -32.645347 209.09516 125000 50 50 50 31000 671.23014 38.678556 148.72308 125000 50 50 50 32000 464.68752 -90.201486 -14.018499 125000 50 50 50 33000 11929.712 727.72594 2683.5372 125000 50 50 50 34000 1249.3459 -505.03535 -300.2119 125000 50 50 50 35000 6864.2858 180.49005 1305.8523 125000 50 50 50 36000 418.5098 -612.91378 -544.30138 125000 50 50 50 37000 174.50306 -677.42917 -648.82034 125000 50 50 50 38000 79.074924 -680.99101 -668.0271 125000 50 50 50 39000 149.36694 -684.89448 -660.40658 125000 50 50 50 40000 34276.857 762.54228 6382.0463 125000 50 50 50 41000 405.71061 -658.47541 -591.96137 125000 50 50 50 42000 52.579346 -689.37579 -680.7557 125000 50 50 50 43000 37.702359 -697.65724 -691.47615 125000 50 50 50 44000 77.636161 -698.36347 -685.63545 125000 50 50 50 45000 1068.4743 -541.76354 -366.593 125000 50 50 50 46000 769.90494 -593.81149 -467.58977 125000 50 50 50 47000 909.93159 -769.17108 -619.99276 125000 50 50 50 Rank 19 [Mon Feb 4 02:49:12 2013] [c13-3c2s3n1] application called MPI_Abort(MPI_COMM_WORLD, 1) - process 19 ERROR on proc 19: Bond atoms 26 27 missing on proc 19 at step 47662 (neigh_bond.cpp:49) 从第3000步开始，能量突然骤减，我调出结构图来一看，是结构出现了问题 |
4楼2013-02-04 11:22:57
5