| 查看: 743 | 回复: 3 | ||
[求助]
液体表面模拟 charmm/namd 盒子无限变长或真空层消失
|
|
建了一个乙腈溶液的模型： ————————PBC box像这样 | vaccum | | solution | | vaccum | 现在用纯乙腈测试，发现在charmm中做NPAT时溶液层比较稳定，真空层却不是持续增长就是减小到消失。NAMD也差不多 所以怀疑是压力控制的问题 用NAMD的时候总是出现 Periodic cell has become too small for original patch grid! charmm的input 就不全贴了，只写几个关键的 set tk 298. crys defi tetr 30. 30. 116.7 90. 90. 90. crys buil nope 0 imag byre sele all end scalar mass stat !calc pmass = int ( ?stot / 50.0 ) set pmass 100 calc tmass = int( ?stot / 50.0 ) * 5 dyna leap cpt restart - ! timestep 0.001 nstep 8000000 nprint 2500 iprfq 2000 - ! iunrea 51 iunwri 41 iuncrd 31 nsavcc 2500 - ! hoover reft @tk tmass @tmass tbat @tk firs @tk tstr @tk - ! pcons pmzz @pmass pmxx 0.0 pmyy 0.0 pref 1.0 - ! pint pref 1.0 pgamma 20.0 - ! ihtfrq 0 ieqfrq 0 ntrfrq 50 - ! iasors 1 iasvel 1 iscvel 0 ichecw 0 - ! ixtf 1000 echeck 100. NAMD 我测试的不多，就先不放了 希望有小伙伴能告知下有什么比较好的做这种体系的参数（1000左右或者几百个分子，体系在30*30*60 anstrong以内），或者能看出我参数里的问题 万分感谢！ |
回复此楼» 猜你喜欢
博士延得我，科研能力直往上蹿
已经有7人回复
-
退学或坚持读
已经有27人回复
-
面上基金申报没有其他的参与者成吗
已经有5人回复
-
有70后还继续奋斗在职场上的吗？
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
-
多组分精馏求助
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2018-04-21 02:48:15
3楼2018-04-22 12:13:22
4楼2018-04-23 15:50:47