| 查看: 1418 | 回复: 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跑出的结果就是不对了？ 求解 求解～  |
回复此楼» 收录本帖的淘帖专辑推荐
 科研之宝，所向披靡 |
» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
» 抢金币啦！回帖就可以得到:
青岛大学/康复大学导师课题组2026年博士研究生
+3/368
-
长江大学化工学院化工学硕、专硕最后一批调剂（理学、工学均可）
+1/94
-
招统计机器学习理论硕士及博士生｜南京大学人工智能学院
+1/84
-
华南理工大学微电子学院招博士研究生
+1/84
-
上海应用技术大学刘敏课题组招生---欢迎材料，物理和化学相关背景
+1/41
-
fluent 化学反应部分预混、三维全尺寸计算慢与压力小问题请教
+1/40
-
广州医科大学生物技术与工程（代码086001）专硕调剂
+5/35
-
温州大学邵银林课题组收调剂专硕研究生
+1/30
-
北京理工大学（珠海）招收人工智能方向26级博士
+1/20
-
河南科技学院——复试调剂—线上复试
+1/19
-
沈阳工业大学石油化工学院传质强化研究团队招收硕士研究生调剂
+1/16
-
粤港澳大湾区量子中心与南方科技大学联合招收2026级博士生（材料、凝聚态物理方向）
+1/16
-
澳门科技大学-系统超分子化学与功能材料课题组 博士生招生简章
+1/10
-
二次调剂名额13名，华北理工大学矿物加工工程方向接收调剂研究生
+1/9
-
湖南理工大学调剂系统已开。 化学、化工、材料化工 调剂指标充足，录取率高，欢迎调剂
+1/7
-
江西水利电力大学土木水利第二轮调剂系统已开
+1/4
-
07专业方向还有调剂名额
+1/2
-
鲁东大学水产学院杨建敏教授课题组，2026年水产学硕接收调剂3-4人
+1/2
-
武汉工程大学化学与环境工程学院招收调剂学生
+1/1
-
海南大学袁禹亮老师招收第二批博士，名额1个，2026年9月份入学（电解（海）水制氢）
+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 [ 发自小木虫客户端 ]
