| 1 | 1/1 | 返回列表 |
| 查看: 692 | 回復: 0 | ||
水兒CoCo木蟲 (小有名氣)
|
[求助]
OPLSAA力場MD模擬時mdp文件
|
|
我采用以上力場進行模擬,但是運行到MD最后一步時終止了,而且還出現了Fatal error: Too many LINCS warnings (1000) 我看網上說是md.mdp參數設置不對,以下是我的文件,希望大神幫我看看哪里要改的,應該是該cut-off參數,但是我不知道這個對不對,謝謝大家啦 ![]() ![]() ![]() title = OPLS MD ; Run parameters integrator = md ; leap-frog integrator nsteps = 5000000 ; 2 * 5000000 = 10000 ps, 10 ns dt = 0.002 ; 2 fs ; Output control nstxout = 1000 ; save coordinates every 2 ps nstvout = 1000 ; save velocities every 2 ps nstxtcout = 1000 ; xtc compressed trajectory output every 2 ps nstenergy = 1000 ; save energies every 2 ps nstlog = 1000 ; update log file every 2 ps energygrps = Protein LC3 ; Bond parameters continuation = yes ; Restarting after NPT 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 cells nstlist = 5 ; 10 fs rlist = 1.0 ; short-range neighborlist cutoff (in nm) rcoulomb = 1.0 ; short-range electrostatic cutoff (in nm) rvdw = 1.0 ; short-range van der Waals cutoff (in nm) ; Electrostatics 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 ; modified Berendsen thermostat tc-grps = Protein_LC3 Water_and_ions ; two coupling groups - more accurate tau_t = 0.1 0.1 ; time constant, in ps ref_t = 310 310 ; reference temperature, one for each group, in K ; Pressure coupling is on pcoupl = Parrinello-Rahman ; Pressure coupling on in NPT pcoupltype = isotropic ; uniform scaling of box vectors tau_p = 2.0 ; time constant, in ps ref_p = 1.0 ; reference pressure, in bar compressibility = 4.5e-5 ; isothermal compressibility of water, 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 |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
孅華 2026-03-22 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 0856材料專碩353求調劑 +5 | NIFFFfff 2026-03-20 | 5/250 |
|
|
[考研] 一志愿南航材料專317分求調劑 +4 | 炸呀炸呀炸薯條 2026-03-23 | 4/200 |
|
|
[考研] 材料/農業(yè)專業(yè),07/08開頭均可,過線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|
|
[考研] 328求調劑 +4 | LHHL66 2026-03-23 | 4/200 |
|
|
[考研] 北科281學碩材料求調劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
[考研] 317求調劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
[考研] 275求調劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
[考研] 考研調劑 +3 | 呼呼?~+123456 2026-03-21 | 3/150 |
|
|
[考研] 0703化學調劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 一志愿深大,0703化學,總分302,求調劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 313求調劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 336求調劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[考研] 308求調劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研] 中南大學化學學碩337求調劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 290求調劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 353求調劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
[考研] 求調劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 生物學調劑招人。! +3 | 山海天嵐 2026-03-17 | 4/200 |
|