| 2 | 1/1 | 返回列表 |
| 查看: 2826 | 回復: 1 | ||
469187892銀蟲 (小有名氣)
|
[求助]
gromacs 運行npt時出錯
|
|
Fatal error: Too many LINCS warnings (2290) If you know what you are doing you can adjust the lincs warning threshold in your mdp file or set the environment variable GMX_MAXCONSTRWARN to -1 but normally it is better to fix the problem 輸入文件是 title = NPT Equilibration for KALP15-DPPC define = -DPOSRES ; position restrain the protein ; Run parameters integrator = md ; leap-frog integrator nsteps = 500000 ; 2 * 500000 = 1000 ps (1 ns) dt = 0.002 ; 2 fs ; Output control nstxout = 1 ; save coordinates every 0.2 ps nstvout = 1 ; save velocities every 0.2 ps nstenergy = 100 ; save energies every 0.2 ps nstlog = 100 ; update log file every 0.2 ps ; Bond parameters freezegrps = HEM_2FE_4FE_3FE_r_57_r_62_r_77_r_65_r_150_r_153_r_217_r_156_r_213_r_160_r_207_r_335_r_442 freezedim = Y Y Y continuation = yes ; Restarting after NVT constraint_algorithm = lincs ; holonomic constraints constraints = all-bonds ; all bonds (even heavy atom-H bonds) constrained lincs_iter = 1 ; accuracy of LINCS lincs_order = 4 ; also related to accuracy ; Neighborsearching ns_type = grid ; search neighboring grid cels nstlist = 5 ; 10 fs rlist = 1.2 ; short-range neighborlist cutoff (in nm) rcoulomb = 1.2 ; short-range electrostatic cutoff (in nm) rvdw = 1.2 ; short-range van der Waals cutoff (in nm) ; Electrostatics nhchainlength =1 coulombtype = PME ; Particle Mesh Ewald for long-range electrostatics pme_order = 4 ; cubic interpolation fourierspacing = 0.16 ; grid spacing for FFT ; Temperature coupling is on tcoupl = v-rescale ; More accurate thermostat tc-grps = Protein_HEM_MOL_2FE_4FE_3FE DPPC SOL_CL ; three coupling groups - more accurate tau_t = 0.5 0.5 0.5 ; time constant, in ps ref_t = 50 50 50 ; reference temperature, one for each group, in K ; Pressure coupling is on pcoupl = Parrinello-Rahman ; Pressure coupling on in NPT pcoupltype = semiisotropic ; uniform scaling of x-y box vectors, independent z tau_p = 10.0 ; time constant, in ps ref_p = 1.0 1.0 ; reference pressure, x-y, z (in bar) compressibility = 4.5e-5 4.5e-5 ; isothermal compressibility, bar^-1 ; Periodic boundary conditions pbc = xyz ; 3-D PBC ; Dispersion correction DispCorr = EnerPres ; account for cut-off vdW scheme ; Velocity generation gen_vel = no ; Velocity generation is off ; COM motion removal ; These options remove motion of the protein/bilayer relative to the solvent/ions nstcalcenergy = 1 comm-mode = Linear comm-grps = Protein_HEM_MOL_2FE_4FE_3FE_DPPC SOL_CL ; Scale COM of reference coordinates refcoord_scaling = com 能幫忙找找錯誤嗎?是參數設置的不合理,還是結構沒有優(yōu)化好? 我從網上下載的DPPC膜,http://moose.bio.ucalgary.ca/ind ... ures_and_Topologies, 還需要做npt, 和md 優(yōu)化嗎? 大家做膜蛋白是先做nvt,還是npt? 問題有點多 |

| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南昌大學324求調劑 +6 | hanamiko 2026-03-27 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 315分求調劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
[考研] 0856材料化工調劑 總分330 +10 | zhubinhao 2026-03-27 | 10/500 |
|
|
[考研] 086000調劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
[考研] 272求調劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 334求調劑 +3 | 雨清天晴 2026-03-21 | 3/150 |
|
|
[考研] 316求調劑 +5 | 江辭666 2026-03-26 | 5/250 |
|
|
[考研] 081200-11408-276學碩求調劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 284求調劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
[考研] 342求調劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
[考研] 一志愿 南京郵電大學 288分 材料考研 求調劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 調劑310 +3 | 溫柔的晚安 2026-03-25 | 4/200 |
|
|
[考研] 334分 一志愿武理-080500 材料求調劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] 食品專碩 一志愿雙一流 328 +3 | xiaom99 2026-03-21 | 4/200 |
|
|
[考研] 一志愿吉大化學322求調劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
[考研] 335求調劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
[考研] 求調劑 +7 | 十三加油 2026-03-21 | 7/350 |
|
|
[考研] 284求調劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
[考研] 306求調劑 +5 | 來好運來來來 2026-03-22 | 5/250 |
|