| 查看: 769 | 回复: 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以内），或者能看出我参数里的问题 万分感谢！ |
回复此楼» 猜你喜欢
0703调剂，一志愿天津大学319分
已经有14人回复
-
085600材料与化工301分求调剂院校
已经有16人回复
-
一志愿211，0703化学305分求调剂
已经有4人回复
-
机械专硕274求调剂，不挑专业学校
已经有8人回复
-
312求调剂
已经有13人回复
-
生物与医药求调剂
已经有7人回复
-
284求调剂
已经有5人回复
-
一志愿北交大材料工程总分358求调剂
已经有10人回复
-
一志愿南科大生物学297分，求调剂推荐
已经有4人回复
-
材料调剂
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2018-04-21 02:48:15
3楼2018-04-22 12:13:22
4楼2018-04-23 15:50:47
