| 查看: 1007 | 回复: 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 |
回复此楼» 猜你喜欢
299求调剂
已经有8人回复
-
一志愿北京理工大学本科211材料工程294求调剂
已经有6人回复
-
300求调剂，材料科学英一数二
已经有8人回复
-
招收生物学/细胞生物学调剂
已经有5人回复
-
070305高分子化学与物理 304分求调剂
已经有7人回复
-
289求调剂
已经有13人回复
-
一志愿哈尔滨工业大学材料与化工方向336分
已经有9人回复
-
081200-11408-276学硕求调剂
已经有6人回复
-
调剂求院校招收
已经有5人回复
-
调剂310
已经有8人回复
hclwonder
新虫 (正式写手)
- 应助: 9 (幼儿园)
- 金币: 2564.1
- 红花: 28
- 帖子: 371
- 在线: 68.2小时
- 虫号: 2837657
- 注册: 2013-11-30
- 性别: GG
- 专业: 药物设计与药物信息
2楼2017-03-23 11:13:41
