| 查看: 1615 | 回复: 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跑出的结果就是不对了？ 求解 求解～  |
回复此楼» 收录本帖的淘帖专辑推荐
 科研之宝，所向披靡 |
» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
上海工程技术大学激光智能制造课题组招收博士研究生
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
» 抢金币啦！回帖就可以得到:
-
[2027 考博] 上海・｜激光焊接 / 增材制造 / 机器人智能装备 招博士
+1/95
-
美国肯塔基大学Dr. Sheng Tong博士后招聘（空间可控的基因组编辑与肿瘤免疫治疗）
+1/94
-
985高校杰青团队招收推免研究生
+1/94
-
【推免】】中科院国家级人才团队招收27级推免硕士生 （材料+机械）
+1/72
-
天津大学 刘浩业教授团队 招收27级博士生
+1/37
-
山东征女友，坐标济南
+1/36
-
上海大学“生物有机电子材料及器件”团队2027年推免硕士/直博研究生招生
+1/21
-
深圳大学2027级光电信息工程/物理学研究生推免
+1/14
-
同样的抗体，为什么你的mIHC背景总比文献高？从固定到成像的全套排查手册
+1/13
-
Biofunctional Materials纳米材料赋能肿瘤新型治疗专题征稿（Scopus、DOAJ收录）
+1/12
-
【生物材料与生物打印专题征稿】Scopus期刊，免APC
+1/5
-
无人机方向 博士后与青年教师招聘公告
+1/4
-
哈工大赵玉伟教授招聘优秀硕博士生及博后
+1/4
-
口腔特刊客座招募客编Biofunctional Materials牙科生物功能材料专题，Scopus期刊
+1/3
-
药物毒理肾损伤相关生物标志物蛋白因子定量检测技术解析
+1/3
-
声场数值模拟与COMSOL有限元水声信号处理：物理信息神经网络PINN与深度学习目标识别
+1/3
-
东南大学招收博士研究生
+1/2
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/2
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/1
-
水利工程博士点师资博士后招聘
+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 [ 发自小木虫客户端 ]
