| 查看: 1433 | 回复: 4 | ||||
[交流]
AMBER力场 improper 参数 在NAMD中不能被读入（0 IMPROPER）
|
||||
|
我用amber力场（ff12SB, ff99SB两个都试过了）在NAMD中跑MD，在NAMD输出log文件中，显示 0 IMPROPER，如下： PARM file in AMBER 7 format Warning: Skipping ATOMIC_NUMBER in parm file while seeking MASS. Warning: Encounter 10-12 H-bond term Warning: Found 87876 H-H bonds. Info: LONG-RANGE LJ: APPLYING ANALYTICAL CORRECTIONS TO ENERGY AND PRESSURE Info: LONG-RANGE LJ: AVERAGE A AND B COEFFICIENTS -1.95755e+07 AND -18374.1 Info: SUMMARY OF PARAMETERS: Info: 104 BONDS Info: 219 ANGLES Info: 139 DIHEDRAL Info: 0 IMPROPER Info: 0 CROSSTERM Info: 0 VDW Info: 325 VDW_PAIRS Info: 0 NBTHOLE_PAIRS Info: READING EXTRA BONDS FILE restraints/step6.1_equilibration.dihe Info: READ 524 EXTRA DIHEDRALS Info: READ 524 EXTRA IMPROPERS Info: TIME FOR READING PDB FILE: 9.53674e-07 Info: Info: **************************** Info: STRUCTURE SUMMARY: Info: 381982 ATOMS Info: 381366 BONDS Info: 221289 ANGLES Info: 455628 DIHEDRALS Info: 0 IMPROPERS Info: 0 CROSSTERMS Info: 914447 EXCLUSIONS Info: 330613 RIGID BONDS Info: 815333 DEGREES OF FREEDOM Info: 139245 HYDROGEN GROUPS Info: 4 ATOMS IN LARGEST HYDROGEN GROUP Info: 139245 MIGRATION GROUPS Info: 4 ATOMS IN LARGEST MIGRATION GROUP Info: TOTAL MASS = 2.33708e+06 amu Info: TOTAL CHARGE = -0.00574125 e Info: MASS DENSITY = 0.994389 g/cm^3 Info: ATOM DENSITY = 0.0978737 atoms/A^3 Info: ***************************** 如上所示，在parameter summary 和 structure summary中，IMPROPER 一项为 0 ！！ 我试过两个力场ff99SB，ff12SB，及 AmberTool1.2、 AmberTool15两个版本，用来生成.prmtop （leap 显示存在 IMPROPER） 结果NAMD读出的IMPROPER都是0 ，虽然NAMD可以正常运行！ 经测试，发现：若用CHARMM力场在NAMD跑，确实给出了相应的IMPROPER数 请问 有没有人 注意到这个奇葩的现象呢？  1. 难道NAMD读AMBER力场时已将IMPROPER纳入到DIHEDRAL了？ 2. 还是压根就是程序本身的bug，果真如此，NAMD跑出的结果就是不对了？ 求解 求解～  |
回复此楼» 收录本帖的淘帖专辑推荐
 科研之宝，所向披靡 |
» 猜你喜欢
260求调剂
已经有5人回复
-
急需调剂
已经有10人回复
-
申博/考博
已经有4人回复
-
化工学硕294分，求导师收留
已经有37人回复
-
求调剂
已经有11人回复
-
一志愿华中农业071010，320求调剂
已经有19人回复
-
304求调剂
已经有7人回复
-
求博导｜生物质基多孔碳/超级电容方向，已有相关成果，寻能源材料/碳材料方向老师
已经有3人回复
-
二苯甲酮酸类衍生物
已经有6人回复
-
接受任何调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
» 抢金币啦！回帖就可以得到:
-
散金祈福
+1/484
-
厦门华侨大学就读的泰国排名第一世界排名221位朱拉隆功大学学位正式招生了
+1/87
-
长江大学化工学院化工学硕、专硕最后一批调剂（理学、工学均可）
+1/83
-
工学门类（前两位代码08）满足国家考研A区工科门类分数线要求的同学看过来
+1/51
-
中国林科院博士招生/读博
+2/46
-
森林工程课题组招调剂学术型硕士研究生，机不可失，期待您的加入。
+1/37
-
山东航空学院安全工程专业085702第三轮接收研究生调剂，欢迎咨询
+1/35
-
瑞典 Linköping University 诚聘 氮化物基无机材料方向 博士后1名
+1/27
-
佛山大学动物科技学院2026年硕士研究生招生调剂公告（第三批）
+1/20
-
中科院合肥物质研究院招聘水系锌电博士后（长期有效）诚聘海内外优秀博士后
+1/18
-
南京林业大学-国家级青年人才团队 诚聘青年教授、博士后等青年才俊
+1/8
-
大连大学—六盘水师范学院联培硕士调剂
+1/8
-
赵松方教授联合山大于伟泳杰青、延世大Jong-Hyun Ahn院士招博士后 、博士研究生
+1/7
-
招07理学调剂生，剩余2个名额
+1/3
-
鲁东大学水产学院杨建敏教授课题组，2026年水产学硕接收调剂3-4人
+1/2
-
九江学院2026年硕士研究生调剂（086002制药工程）
+1/2
-
Advanced Synthesis & Catalysis 决定是否送审一般多久？
+1/2
-
大连海洋大学联合培养硕士调剂招生
+1/2
-
安徽工程大学化学与环境工程学院2026级硕士研究生还有部分调剂名额-资源与环境专硕
+1/1
-
华南理工大学微电子学院博士招生（射频模拟混合信号芯片、天线超表面/AI方向）
+1/1
★ ★ ★
vip_zzc(金币+2): 谢谢参与
月只蓝: 金币+2, 鼓励交流！ 2014-02-18 18:27:42
vip_zzc(金币+2): 谢谢参与
月只蓝: 金币+2, 鼓励交流！ 2014-02-18 18:27:42
|
本帖内容被屏蔽 |
3楼2014-02-18 17:28:47
4楼2014-02-18 20:31:52
★
月只蓝: 金币+1, 鼓励交流！ 2014-02-19 11:01:55
月只蓝: 金币+1, 鼓励交流！ 2014-02-19 11:01:55
|
grep FLAG .prmtop 后 %FLAG TITLE %FLAG POINTERS %FLAG ATOM_NAME %FLAG CHARGE %FLAG ATOMIC_NUMBER %FLAG MASS %FLAG ATOM_TYPE_INDEX %FLAG NUMBER_EXCLUDED_ATOMS %FLAG NONBONDED_PARM_INDEX %FLAG RESIDUE_LABEL %FLAG RESIDUE_POINTER %FLAG BOND_FORCE_CONSTANT %FLAG BOND_EQUIL_VALUE %FLAG ANGLE_FORCE_CONSTANT %FLAG ANGLE_EQUIL_VALUE %FLAG DIHEDRAL_FORCE_CONSTANT %FLAG DIHEDRAL_PERIODICITY %FLAG DIHEDRAL_PHASE %FLAG SCEE_SCALE_FACTOR %FLAG SCNB_SCALE_FACTOR %FLAG SOLTY %FLAG LENNARD_JONES_ACOEF %FLAG LENNARD_JONES_BCOEF %FLAG BONDS_INC_HYDROGEN %FLAG BONDS_WITHOUT_HYDROGEN %FLAG ANGLES_INC_HYDROGEN %FLAG ANGLES_WITHOUT_HYDROGEN %FLAG DIHEDRALS_INC_HYDROGEN %FLAG DIHEDRALS_WITHOUT_HYDROGEN %FLAG EXCLUDED_ATOMS_LIST %FLAG HBOND_ACOEF %FLAG HBOND_BCOEF %FLAG HBCUT %FLAG AMBER_ATOM_TYPE %FLAG TREE_CHAIN_CLASSIFICATION %FLAG JOIN_ARRAY %FLAG IROTAT %FLAG RADIUS_SET %FLAG RADII %FLAG SCREEN %FLAG IPOL 我们会发现在amber .prmtop文件中，并没有IMPROPER这一标识，说明AMBER应该是将 IMPROPER 和 PROPER 整合在一起了，NAMD根据FLAG读取amber力场参数，故只能得到 DIHEDRAL，我认为 NAMD .log 输出中给出的 DIHEDRAL = PROPER + IMPROPER （当然，这是对amber FF作为输入而言） 转 http://simulation.haotui.com/thread-45212-1-1.html |
5楼2014-02-18 21:22:52
简单回复
2014-02-18 15:55
回复
vip_zzc(金币+2): 谢谢参与
1 [ 发自小木虫客户端 ]
