| 查看: 1829 | 回复: 1 | ||
| 【悬赏金币】回答本帖问题，作者来围观啊将赠送您 30 个金币 | ||
来围观啊新虫 (小有名气)
|
[求助]
GROMCS出现错误
|
|
|
求助各位大佬，在用gromacs的时候出现了以下的错误该如何解决呢？体系已经做了能量优化，这个错误是在做nvt预平衡的时候出现的。谢谢各位大佬！！ Step 1325, time 1.325 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 0.035172, max 0.281169 (between atoms 142 and 145) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 143 144 31.2 0.1005 0.0989 0.1000 Step 1331, time 1.331 (ps) LINCS WARNING relative constraint deviation after LINCS: rms 816.149048, max 9026.031250 (between atoms 143 and 144) bonds that rotated more than 30 degrees: atom 1 atom 2 angle previous, current, constraint length 103 127 48.7 0.1453 0.1471 0.1435 114 115 49.4 0.1566 0.0854 0.1435 114 124 98.6 0.1649 0.1101 0.1520 117 118 38.8 0.1835 0.3072 0.1520 117 141 106.8 0.1585 0.1745 0.1435 118 119 45.0 0.1661 0.0786 0.1435 118 126 67.7 0.1792 0.0449 0.1520 119 120 37.4 0.1027 0.1192 0.1000 121 124 35.8 0.1524 0.2155 0.1520 124 125 38.1 0.1444 0.1945 0.1435 125 126 30.7 0.1466 0.1670 0.1435 128 129 176.2 0.2776 5.3298 0.1435 128 131 69.7 0.5636 4.4779 0.1520 128 138 79.8 0.3747 3.8723 0.1520 129 130 93.1 0.0182 2.2967 0.1000 131 132 56.9 0.3678 1.8625 0.1520 131 155 118.5 1.1285 4.6649 0.1435 132 133 105.0 0.1806 15.4852 0.1435 132 140 48.7 0.0869 0.9637 0.1520 133 134 98.4 0.1397 15.0138 0.1000 135 136 87.5 0.1530 1.4493 0.1435 135 138 82.6 0.1740 2.9280 0.1520 136 137 92.0 0.1226 0.7336 0.1000 138 139 177.3 0.1450 16.4410 0.1435 139 140 59.1 0.1621 16.2713 0.1435 140 141 95.8 0.1300 0.8640 0.1435 142 145 99.8 4.3105 38.1381 0.1520 142 152 149.1 4.6295 6.3105 0.1520 145 146 31.5 2.1491 40.9124 0.1520 145 171 49.5 2.5563 41.1024 0.1435 146 147 93.5 3.6932 3.2399 0.1435 146 154 33.9 2.4452 7.8218 0.1520 147 148 89.6 1.6534 1.7318 0.1000 149 150 92.5 0.8919 0.8930 0.1435 149 152 86.2 4.6661 10.1095 0.1520 150 151 92.4 0.4329 1.0297 0.1000 152 153 64.5 5.4649 13.5725 0.1435 153 154 156.7 4.2225 6.0705 0.1435 154 155 31.8 3.5694 6.1695 0.1435 157 158 127.9 0.1700 0.2188 0.1000 159 162 34.4 0.7675 1.7566 0.1520 160 161 58.5 0.3207 0.2597 0.1000 162 163 62.7 0.5360 0.5178 0.1435 162 170 84.9 2.2265 4.8269 0.1520 163 164 74.8 0.6339 0.2147 0.1000 166 167 55.6 0.1750 0.0975 0.1000 168 169 69.1 0.4028 1.2414 0.1435 169 170 85.7 2.0814 3.6421 0.1435 170 171 124.3 1.6610 4.1818 0.1435 170 171 124.3 1.6610 4.1818 0.1435 170 171 124.3 1.6610 4.1818 0.1435 step 1331: Water molecule starting at atom 4837 can not be settled. Check for bad contacts and/or reduce the timestep if appropriate. Back Off! I just backed up step1331b.pdb to ./#step1331b.pdb.1# Back Off! I just backed up step1331c.pdb to ./#step1331c.pdb.1# Wrote pdb files with previous and current coordinates 段错误 ******************************************************************************************************************************************** 以下为minim.mdp的参数 ; minim.mdp - used as input into grompp to generate em.tpr ; Parameters describing what to do, when to stop and what to save integrator = steep ; Algorithm (steep = steepest descent minimization) emtol = 1000.0 ; Stop minimization when the maximum force < 1000.0 kJ/mol/nm emstep = 0.01 ; Energy step size nsteps = 50000 ; Maximum number of (minimization) steps to perform ; Parameters describing how to find the neighbors of each atom and how to calculate the interactions nstlist = 1 ; Frequency to update the neighbor list and long range forces ns_type = grid ; Method to determine neighbor list (simple, grid) rlist = 1.4 ; Cut-off for making neighbor list (short range forces) coulombtype = PME ; Treatment of long range electrostatic interactions rcoulomb = 1.4 ; Short-range electrostatic cut-off rvdw = 1.4 ; Short-range Van der Waals cut-off pbc = xyz ; Periodic Boundary Conditions (yes/no) ******************************************************************************************************************************************* 以下为nvt restraint的预平衡参数nvt.mdp： define = -DPOSRES ; position restrain the protein integrator = sd dt = 0.001 ; ps nsteps =100000 ; = tot nstcomm = 10000 nstxout = 10000 nstvout = 10000 nstfout = 10000 nstlog = 10000 nstenergy = 1000 nstlist = 5 pbc = xyz ns_type = grid coulombtype = PME rcoulomb = 1.2 pme_order = 5 ; cubic interpolation fourierspacing = 0.175 ; grid spacing for FFT vdwtype = cut-off rvdw = 1.2 rlist = 1.2 DispCorr = EnerPres ;optimize_fft = yes ;epsilon_r = 80 ;Couplig Temp tcoupl = v-rescale tau_t = 0.2 tc-grps = System ref_t = 298 ;Costrain Bond constraints = all-bonds constraint_algorithm = LINCS ;Coupling P ;Pcoupl = berendsen ;Pcoupltype = isotropic ;tau_p = 0.5 ;compressibility = 4.5e-5 ;ref_p=1 ;Generate velocity gen_vel = yes gen_temp = 298.0 gen_seed = 333529@smutao@smutao@oxox6085@oxox6085 |
回复此楼» 猜你喜欢
召集青年编委
已经有7人回复
-
Chemical Engineering Journal投稿3周了，一直显示With editor状态。这是送审了吗？
已经有10人回复
-
0854－0855调剂
已经有10人回复
-
南京大学-广州大学联合招聘博士后 欢迎广大优秀人才！！！
已经有11人回复
-
申博求助
已经有8人回复
-
Nature一直在编辑手里，考虑好几天了，是悬了吗
已经有7人回复
-
博士生招生！！
已经有5人回复
-
国社科项目，你们学校都限额申报吗？
已经有8人回复
-
企业博后是否能申请CSC博士后项目？谢谢
已经有6人回复
-
博士白读了
已经有40人回复
hainu_carl
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4930.8
- 散金: 154
- 红花: 1
- 帖子: 710
- 在线: 36.7小时
- 虫号: 5502821
- 注册: 2017-01-20
- 专业: 食品科学
2楼2020-09-17 07:07:57
10