应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps 离子液体，离子中相邻隔的挺远的原子成键了，不知什么原因
51/1返回列表
查看: 1763  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wup123

铁虫 (小有名气)

[求助] lammps 离子液体，离子中相邻隔的挺远的原子成键了，不知什么原因

我在进行模拟的时候，首先用两个分子（四个离子）测试了，结构非常好，没啥异常。但是增多原子后，发现其中一个阴离子（tfsi-）中F原子和相隔非常远的O原子（相隔五个原子）成键了，从而使整个离子变的乱起八糟。
   我检查了力场参数，采用的是opls—aa势，和文献中一致，电荷也守恒。
  想问问各位虫子，这个是什么原因呢？该怎么解决呢？
  也不知怎么搞的，我遇到好多诡异的问题。。。
回复此楼

» 本帖已获得的红花（最新10朵）

凌枫霜叶

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-02-01 11:32:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wup123

铁虫 (小有名气)
离子结构图和结果图如下， fsi中的F原子和emim中的Hoverlap了，导致能量变为负数了。
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 结构1.jpg
    • 2013-02-04 11:30:25, 11.79 K
  • 附件 2 : result.jpg
    • 2013-02-04 11:30:57, 21.26 K
5楼2013-02-04 11:32:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
把你的结构图，in文件， log文件，坐标文件都贴出来。 否则怎么看呀。
赞一下 回复此楼
2楼2013-02-04 10:27:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wup123

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by ddll_slsl at 2013-02-04 10:27:01
把你的结构图，in文件， log文件，坐标文件都贴出来。 否则怎么看呀。

好的，谢谢，马上贴上。
赞一下 回复此楼
3楼2013-02-04 11:15:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wup123

铁虫 (小有名气)
我模拟的体系是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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 数学博导 +5 学术霸王 2024-09-25 8/400 2024-09-29 07:23 by 半生梦君
[教师之家] 东南大学电气学院电力电子方向招收硕士博士 +5 蜡笔小鑫1989 2024-09-26 8/400 2024-09-29 06:55 by willbuilder
[硕博家园] 希望有师兄师姐帮忙引荐申博 +6 zhanyaqian 2024-09-25 11/550 2024-09-28 22:24 by 自强2019
[教师之家] 开心！国庆前发个小财~ +15 zjjxzh 2024-09-26 15/750 2024-09-28 19:54 by bear2007
[论文投稿] 求推荐医学类杂志 35+4 wshxtim1 2024-09-25 4/200 2024-09-28 13:56 by Wormaciae
[基金申请] 山东省基金公示结束 +6 ishow19 2024-09-24 14/700 2024-09-28 12:20 by ishow19
[教师之家] 西湖大学教授：我可以自由地选择讲课的方式, 讲授的角度 +4 zju2000 2024-09-22 4/200 2024-09-28 09:40 by bio-polymer
[基金申请] 广东省基金 +5 yingyueabc 2024-09-22 5/250 2024-09-27 22:57 by FengjiangLiu
[论文投稿] Journal of colloid and interface science期刊 10+4 1821588220 2024-09-24 10/500 2024-09-27 20:37 by 1821588220
[基金申请] 请问大家的计划书填写列表中状态更新了吗？ +7 Laker610 2024-09-25 9/450 2024-09-27 16:36 by 田田hj
[教师之家] 压力山大 +17 暮色恋伊人 2024-09-24 25/1250 2024-09-27 16:10 by 领先于123
[硕博家园] 当前读博士还有性价比吗？ +8 苏东坡二世 2024-09-23 8/400 2024-09-27 11:39 by 半生梦君
[考博] 电催化，电池方向学生自荐 30+3 Freelancedd 2024-09-24 10/500 2024-09-27 10:27 by 莱茵润色
[论文投稿] 期刊论文发表了但查询到sci未收录 50+5 yibuxiao 2024-09-25 11/550 2024-09-27 10:24 by 莱茵润色
[考博] 有机转材料申博难度 30+3 昕散 2024-09-22 10/500 2024-09-26 22:42 by 824282658
[有机交流] 请问胺的盐酸盐中氯化氢的氢会在核磁氢谱中出峰吗？ +3 rommel1975 2024-09-25 4/200 2024-09-26 13:07 by 091602
[基金申请] 请教2024后期资助到哪个阶段了 +3 拾光者5566 2024-09-23 8/400 2024-09-26 09:27 by 拾光者5566
[论文投稿] laser physics期刊投稿 5+3 mengxiangcz 2024-09-23 6/300 2024-09-25 15:08 by mengxiangcz
[论文投稿] 爱斯维尔旗下Postharvest Biology and Technology，采用什么系统查查 20+4 sdsdsxh 2024-09-22 9/450 2024-09-25 07:33 by sdsdsxh
[有机交流] 做聚酰亚胺一步法合成，机械搅拌怎么密封不漏气 +3 哼哈EN 2024-09-24 5/250 2024-09-24 14:16 by mrzhl1986
信息提示
关闭
请填处理意见
关闭