| 查看: 1703 | 回复: 1 | |||
[交流]
【求助】gromacs下对大分子体系QM/MM计算的问题
|
|
我用的是gromacs4.0.3+MOPAC7.0 模拟对象是MHC/VIRAL PEPTIDE/T-CELL RECEPTOR三者复合体(1AO7)，五条链，一共11000个原子，加上水和离子后270000左右的原子。 主要是计算MHC和VIRAL PEPTIDE的结合能，在这两者上分别取若干氨基酸用QM计算，其余用MM计算。整个体系能量最小化到1000一下，开始模拟是有如下错误， RHF CALCULATION, NO. OF DOUBLY OCCUPIED LEVELS = 35 keywords are: PRECISE GEO-OK CHARGE=0 GRAD MMOK ANALYT AM1 Back Off! I just backed up qmmm.trr to ./#qmmm.trr.1# Back Off! I just backed up qmmm.xtc to ./#qmmm.xtc.1# Back Off! I just backed up qmmm.edr to ./#qmmm.edr.1# WARNING: This run will generate roughly 3863 Mb of data starting mdrun 'HLA-A 0201' 5000 steps, 5.0 ps. counting the number of NH-C=O groups step 0 Step 44, time 0.044 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 0.000010, max 0.000487 (between atoms 2254 and 2257) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 3941 3942 34.1 0.0945 0.0945 0.0945 Step 45, time 0.045 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 0.000013, max 0.000668 (between atoms 2254 and 2257) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 3941 3942 41.1 0.0945 0.0945 0.0945 Step 46, time 0.046 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 0.000017, max 0.000853 (between atoms 2254 and 2256) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 3941 3942 49.3 0.0945 0.0945 0.0945 Step 47, time 0.047 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 0.000023, max 0.001332 (between atoms 2254 and 2256) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 3941 3942 58.9 0.0945 0.0945 0.0945 ALL CONVERGERS ARE NOW FORCED ON SHIFT=10, PULAY ON, CAMP-KING ON AND ITERATION COUNTER RESET MAXIMUM SIZE: 9 EMAT 0.100000E+01-0.100000E+01 -0.100000E+01 0.000000E+00 EMAT Infinity 0.000000E+00 0.000000E+00 0.000000E+00 Step 48, time 0.048 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 0.000034, max 0.001928 (between atoms 2254 and 2256) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 3941 3942 70.1 0.0945 0.0945 0.0945 Step 49, time 0.049 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 0.000059, max 0.002527 (between atoms 3941 and 3942) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 3941 3942 90.0 0.0945 0.0947 0.0945 2254 2256 30.9 0.1012 0.1012 0.1010 t = 0.049 ps: Water molecule starting at atom 272720 can not be settled. Check for bad contacts and/or reduce the timestep. Back Off! I just backed up step49b.pdb to ./#step49b.pdb.1# Back Off! I just backed up step49c.pdb to ./#step49c.pdb.1# Wrote pdb files with previous and current coordinates Step 50, time 0.05 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 0.000399, max 0.028944 (between atoms 3941 and 3942) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 3941 3942 90.0 0.0947 0.0972 0.0945 2254 2256 34.1 0.1012 0.1012 0.1010 t = 0.050 ps: Water molecule starting at atom 272852 can not be settled. Check for bad contacts and/or reduce the timestep. Back Off! I just backed up step50b.pdb to ./#step50b.pdb.1# Back Off! I just backed up step50c.pdb to ./#step50c.pdb.1# Wrote pdb files with previous and current coordinates Step 51, time 0.051 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 0.000548, max 0.039885 (between atoms 3941 and 3942) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 3941 3942 90.0 0.0972 0.0983 0.0945 2254 2256 35.2 0.1012 0.1012 0.1010 t = 0.051 ps: Water molecule starting at atom 272720 can not be settled. Check for bad contacts and/or reduce the timestep. Back Off! I just backed up step51b.pdb to ./#step51b.pdb.1# Back Off! I just backed up step51c.pdb to ./#step51c.pdb.1# Wrote pdb files with previous and current coordinates Step 52, time 0.052 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 0.000232, max 0.016477 (between atoms 3941 and 3942) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 3941 3942 90.0 0.0983 0.0961 0.0945 2254 2256 33.8 0.1011 0.1010 0.1010 ALL CONVERGERS ARE NOW FORCED ON SHIFT=10, PULAY ON, CAMP-KING ON AND ITERATION COUNTER RESET MAXIMUM SIZE: 2 EMAT 0.100000E+01-0.100000E+01 -0.100000E+01 0.000000E+00 EMAT Infinity 0.000000E+00 0.000000E+00 0.000000E+00 Step 53, time 0.053 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 0.000045, max 0.002539 (between atoms 2254 and 2257) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 3941 3942 84.1 0.0960 0.0945 0.0945 2254 2256 30.3 0.1010 0.1009 0.1010 t = 0.053 ps: Water molecule starting at atom 272720 can not be settled. Check for bad contacts and/or reduce the timestep. Back Off! I just backed up step53b.pdb to ./#step53b.pdb.1# Back Off! I just backed up step53c.pdb to ./#step53c.pdb.1# Wrote pdb files with previous and current coordinates ALL CONVERGERS ARE NOW FORCED ON SHIFT=10, PULAY ON, CAMP-KING ON AND ITERATION COUNTER RESET """""""""""""UNABLE TO ACHIEVE SELF-CONSISTENCE DELTAE= -0.1994E-03 DELTAP= NaN ------------------------------------------------------------------------------- PRECISE GEO-OK CHARGE=0 GRAD MMOK ANALYT AM1 A FAILURE HAS OCCURRED, TREAT RESULTS WITH CAUTION!! ++++----**** FAILED TO ACHIEVE SCF. ****----++++ AM1 CALCULATION VERSION 7.00 Sun Mar 27 22:07:56 2011 FOR SOME REASON THE SCF CALCULATION FAILED. THE RESULTS WOULD BE MEANINGLESS, SO WILL NOT BE PRINTED. TRY TO FIND THE REASON FOR THE FAILURE BY USING "PL". CHECK YOUR GEOMETRY AND ALSO TRY USING SHIFT OR PULAY. ATOM CHEMICAL BOND LENGTH BOND ANGLE TWIST ANGLE NUMBER SYMBOL (ANGSTROMS) (DEGREES) (DEGREES) (I) NA:I NB:NA:I NC:NB:NA:I NA NB NC 1 N 2 H 0.00000 0 3 C 0.00000 0.00000 0 0 4 H 0.00000 0.00000 0.00000 0 0 0 5 C 0.00000 0.00000 0.00000 0 0 0 6 H 0.00000 0.00000 0.00000 0 0 0 7 H 0.00000 0.00000 0.00000 0 0 0 8 H 0.00000 0.00000 0.00000 0 0 0 9 C 0.00000 0.00000 0.00000 0 0 0 10 O 0.00000 0.00000 0.00000 0 0 0 11 H 0.00000 0.00000 0.00000 0 0 0 12 H 0.00000 0.00000 0.00000 0 0 0 13 N 0.00000 0.00000 0.00000 0 0 0 14 C 0.00000 0.00000 0.00000 0 0 0 15 H 0.00000 0.00000 0.00000 0 0 0 16 C 0.00000 0.00000 0.00000 0 0 0 17 H 0.00000 0.00000 0.00000 0 0 0 18 H 0.00000 0.00000 0.00000 0 0 0 19 C 0.00000 0.00000 0.00000 0 0 0 20 H 0.00000 0.00000 0.00000 0 0 0 21 H 0.00000 0.00000 0.00000 0 0 0 22 C 0.00000 0.00000 0.00000 0 0 0 23 H 0.00000 0.00000 0.00000 0 0 0 24 H 0.00000 0.00000 0.00000 0 0 0 25 C 0.00000 0.00000 0.00000 0 0 0 26 O 0.00000 0.00000 0.00000 0 0 0 27 H 0.00000 0.00000 0.00000 0 0 0 28 H 0.00000 0.00000 0.00000 0 0 0 另外求教对大体系能量最小化应该怎么设置比较合理？emtol设到多少比较合适？constraints应该选用何种算法？ |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
求教各位大神：2026教育部人文社科青年基金项目何时公示呀？
已经有7人回复
-
浙江师范大学是怎么坑我的
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 氨基酸在线衍生测定过程中遇到的一些问题
已经有6人回复
- 本人查阅文献，遇到一个问题，请教哪位给我解答哈，在下将不胜感激！
已经有18人回复
- 陕西省榆林市大能源体系建设2008-2018年重大项目规划表
已经有34人回复
- 面试时应对跳槽原因问题的诀窍
已经有7人回复
- 关于熵变的计算，有人知道怎么计算大的体系的熵变吗？比如将近一千个残基的体系？
已经有11人回复
- 分子模拟时间越长越好吗？
已经有44人回复
- AMBER中qm/mm的问题
已经有5人回复
- 求助各位，关于初学分子模拟
已经有9人回复
- 【求助】gromacs的又一个问题
已经有22人回复
- 【求助】psf文件的生成
已经有19人回复
» 抢金币啦！回帖就可以得到:
-
华南理工大学 光电材料方向 王竣可教授课题组招聘博士后、博士生及科研助理
+5/435
-
坐标北京，诚征女友，奔结婚，非诚勿扰！
+5/250
-
材料理论模拟或理论计算的作用有多大
+1/197
-
中国科学院上海有机化学研究所左智伟课题组招收2027年入学博士研究生
+1/89
-
山东征女友，坐标济南
+1/38
-
石河子大学211双一流学科分子筛等多孔材料与生物质资源化利用方向招2027级推免生
+1/32
-
211诚招材化相关专业硕士、博士、博后（长期有效）推免-保研-调剂-申博
+1/31
-
浙江师范大学引进人才被解聘，拒绝支付50万绩效并被索赔211万
+1/29
-
27博士申请自荐
+1/21
-
深圳大学应用技术学院招聘凝聚态物理博士后
+1/16
-
中科院西北研究院招收遥感、地信相关专业联培研究生及科研助理！
+1/11
-
欢迎推免生/准推免生报名（约80名每年）--长期有效--物理/电子等所有理工专业
+1/9
-
无人机方向 博士后与青年教师招聘公告
+1/8
-
【特刊征稿】BM临床诊断与治疗生物功能材料专题（Scopus、DOAJ 收录）
+1/5
-
【牙科生物材料特刊征稿】对牙科方向友好（Scopus、DOAJ收录的正规国际期刊）
+1/5
-
北理工集成电路杰青团队 | 诚招科助理
+1/4
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/4
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/4
-
口腔特刊客座招募客编Biofunctional Materials牙科生物功能材料专题，Scopus期刊
+1/4
-
哈工大（深圳）招收工程热物理专业2027入学博士生
+1/1
2楼2011-06-07 21:46:58