| 查看: 2785 | 回复: 1 | |||
[求助]
NAMD中NVT和NPT系综问题
|
|
系统做完优化之后，开始跑分子动力学时整个系统（水箱+研究体系）突然缩小。关于类似的问题，有人说可以先跑NVT，然后再跑NPT。 不知道是否有人了解这类问题，跪求指点。 我是在NAMD中利用CHARMM立场做的计算。 输入的conf文件如下： ############################################################# ## JOB DESCRIPTION ## ############################################################# # Minimization and Equilibration of # the whole complex in a Water Box ############################################################# ## ADJUSTABLE PARAMETERS ## ############################################################# structure ionized.psf coordinates ionized.pdb set temperature 300 set outputname whole_out firsttimestep 0 ############################################################# ## SIMULATION PARAMETERS ## ############################################################# # Input paraTypeCharmm on parameters par_all27_prot_lipid_na.inp temperature $temperature # Force-Field Parameters exclude scaled1-4 1-4scaling 1.0 cutoff 14.0 switching on switchdist 10.0 pairlistdist 16.0 # Integrator Parameters timestep 2 ;# 2fs/step rigidBonds all ;# needed for 2fs steps nonbondedFreq 1 fullElectFrequency 2 stepspercycle 10 # Constant Temperature Control langevin on ;# do langevin dynamics langevinDamping 1 ;# damping coefficient (gamma) of 1/ps langevinTemp $temperature langevinHydrogen off ;# don't couple langevin bath to hydrogens # Periodic Boundary Conditions cellBasisVector1 152.1 0.0 0.0 cellBasisVector2 0.0 158.9 0.0 cellBasisVector3 0.0 0.0 152.6 cellOrigin -52.6 -16.9 -29.8 wrapAll on # PME (for full-system periodic electrostatics) PME yes PMEGridSpacing 1.0 #manual grid definition #PMEGridSizeX 45 #PMEGridSizeY 45 #PMEGridSizeZ 48 # Constant Pressure Control (variable volume) useGroupPressure yes ;# needed for rigidBonds useFlexibleCell yes useConstantArea no langevinPiston on langevinPistonTarget 1.01325 ;# in bar -> 1 atm langevinPistonPeriod 100.0 langevinPistonDecay 50.0 langevinPistonTemp $temperature # Output outputName $outputname restartfreq 500 ;# 500steps = every 1ps dcdfreq 250 xstFreq 250 outputEnergies 100 outputPressure 100 ############################################################# ## EXTRA PARAMETERS ## ############################################################# ############################################################# ## EXECUTION SCRIPT ## ############################################################# # Minimization minimize 20000 for { set TEMP 0 } { $TEMP < 300 } { incr TEMP 30 } { reinitvels $TEMP langevinTemp $TEMP run 5000 } reinitvels $temperature langevinTemp 300 langevinPistonTemp 300 run 1000000 ;# 2ns 对于NVT和NPT的conf文件的写法，也需要牛人指教，不胜感激！ |
回复此楼» 本帖已获得的红花（最新10朵）
» 猜你喜欢
307求调剂
已经有3人回复
-
085600材料与化工专硕329 求调剂
已经有14人回复
-
材料考研求调剂总分280
已经有17人回复
-
085100建筑学 寻求跨专业调剂 一志愿南大294分 校级省级国家级奖项若干 踏实肯干
已经有4人回复
-
265求调剂
已经有15人回复
-
化工求调剂！
已经有20人回复
-
266求调剂
已经有22人回复
-
化工调剂303分，过四级
已经有29人回复
-
生物学363调剂求助
已经有4人回复
-
22408 一志愿双一流人工智能300分 四六级，数据分析国奖
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
送红花一朵 |
本帖内容被屏蔽 |
2楼2020-12-14 23:19:46
