应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Gromacs/Amber/NAMD » 【求助】gromacs下对大分子体系QM/MM计算的问题
21/1返回列表
查看: 1693  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

come4c

铁虫 (初入文坛)

[交流] 【求助】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应该选用何种算法？
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-03-29 10:58:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sijifengsd

银虫 (小有名气)
★
come4c(金币+1):谢谢参与
呵呵，请问您用Gromacs设置qmmm连接原子的时候虚原子怎么写入gro文件的？这个问题我很感兴趣，但不知道怎么处理，谢谢赐教
赞一下(1人) 回复此楼
2楼2011-06-07 21:46:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 come4c 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +11 wangmingqi 2026-08-28 15/750 2026-08-29 00:17 by fangyl2005
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 17:09 by winsaint
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
关闭
请填处理意见
关闭