| 查看: 1076 | 回复: 1 | |||
[求助]
NAMD，完全一样的conf文件，但结果差异显著
|
|
HELP！ 如图所示，用NAMD，在三个温度下（不同颜色），分别跑了5次，但是组内差异甚至超过了组间差异。 一模一样的模拟，跑出来的结果为什么会差异这么大呢？ 力场不合适吗？但换成了36，并且用同一个能量最小化的结构作为初始结构，也是类似的结果。 万分感谢！ ############################################################# ## JOB DESCRIPTION ## ############################################################# # NAMD Configuration for Equilibration (NPT) of # X in a Water Box # 2017-02-18 ############################################################# ## SELF-DEFINED VARIABLE ## ############################################################# set STRUCT X_wb_10 set LAST_RUN X_wb set THIS_RUN X_md1_45_1 set FTS 0 ;# firsttimestep set TEMP 318 ; # temperature set NMD 2500000 ;# 5ns ############################################################# ## ADJUSTABLE PARAMETERS ## ############################################################# structure ../common/$STRUCT.psf coordinates ../common/$STRUCT.pdb firsttimestep $FTS ############################################################# ## SIMULATION PARAMETERS ## ############################################################# ## Input paraTypeCharmm on parameters ../common/par_all27_prot_lipid.inp temperature $TEMP ## Force-Field Parameters exclude scaled1-4 1-4scaling 1.0 cutoff 12.0 switching on switchdist 10.0 ;# cutoff - 2.0 pairlistdist 14.0 ;# cutoff + 2.0 ## Integrator Parameters timestep 2.0 ;# 2fs/step rigidBonds all ;# needed for 2fs steps nonbondedFreq 1 fullElectFrequency 2 stepspercycle 10 ## Constant Temperature Control tCouple on tCoupleTemp $TEMP ## Periodic Boundary Conditions (PBC) cellBasisVector1 71.44 0.0 0.0 cellBasisVector2 0.0 62.92 0.0 cellBasisVector3 0.0 0.0 79.15 cellOrigin 24.71 15.60 6.49 wrapAll on ## PME (for full-system periodic electrostatics) PME yes PMEGridSpacing 1.0 # manual grid definition # PMEGridSizeX 65 # PMEGridSizeY 54 # PMEGridSizeZ 70 ## Constant Pressure Control (variable volume) langevinPiston on langevinPistonTarget 1.01325 ;# in bar -> 1 atm langevinPistonPeriod 500.0 langevinPistonDecay 250.0 langevinPistonTemp $TEMP useGroupPressure yes ;# needed for rigidBonds useFlexibleCell no useConstantArea no ## Output outputName $THIS_RUN restartfreq 5000 ;# 2500steps = every 5ps dcdfreq 2500 xstFreq 2500 outputEnergies 2500 outputPressure 2500 ############################################################# ## EXTRA PARAMETERS ## ############################################################# ############################################################# ## EXECUTION SCRIPT ## ############################################################# # Minimization minimize 50000 reinitvels $TEMP run $NMD  捕获.PNG |
回复此楼» 猜你喜欢
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
青基如何破局
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有7人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
hclwonder
新虫 (正式写手)
- 应助: 9 (幼儿园)
- 金币: 2564.1
- 红花: 28
- 帖子: 371
- 在线: 68.2小时
- 虫号: 2837657
- 注册: 2013-11-30
- 性别: GG
- 专业: 药物设计与药物信息
2楼2017-03-23 11:13:41
