| 查看: 730 | 回复: 5 | |||
[交流]
【求助】Melting of Lipid Tails后脂双层有的地方分离已有2人参与
|
Melting of Lipid Tails后脂双层有的地方分离 这是我的conf文件： ## JOB DESCRIPTION ## ############################################################# # Min. and Eq. of Aqp1 # embedded in POPE membrane, water. # Melting lipid tails. PME, Constant Volume. ############################################################# ## ADJUSTABLE PARAMETERS ## ############################################################# structure ../02-MEMBRANE/aqp1_popew.psf coordinates ../02-MEMBRANE/aqp1_popew.pdb outputName aqp1_popewmineq-01 set temperature 300 # Continuing a job from the restart files if {0} { set inputname aqp1-popew binCoordinates $inputname.restart.coor binVelocities $inputname.restart.vel ;# remove the "temperature" entry if you use this! extendedSystem $inputname.restart.xsc } firsttimestep 0 ############################################################# ## SIMULATION PARAMETERS ## ############################################################# # Input paraTypeCharmm on parameters ../par_all27_prot_lipid.prm # NOTE: Do not set the initial velocity temperature if you # have also specified a .vel restart file! temperature $temperature # Periodic Boundary Conditions # NOTE: Do not set the periodic cell basis if you have also # specified an .xsc restart file! if {1} { cellBasisVector1 105. 0. 0. cellBasisVector2 0. 105. 0. cellBasisVector3 0. 0. 73. cellOrigin -0.077429056167 0.074017807841 -0.44251784668 } wrapWater on wrapAll on # Force-Field Parameters exclude scaled1-4 1-4scaling 1.0 cutoff 12. switching on switchdist 10. pairlistdist 13.5 # Integrator Parameters timestep 2.0 ;# 2fs/step rigidBonds all ;# needed for 2fs steps nonbondedFreq 1 fullElectFrequency 2 stepspercycle 20 #PME (for full-system periodic electrostatics) if {1} { PME yes PMEGridSizeX 106 PMEGridSizeY 106 PMEGridSizeZ 74 } # Constant Temperature Control langevin on ;# do langevin dynamics langevinDamping 1 ;# damping coefficient (gamma) of 5/ps langevinTemp $temperature # Constant Pressure Control (variable volume) if {0} { useGroupPressure yes ;# needed for 2fs steps useFlexibleCell yes ;# no for water box, yes for membrane useConstantArea no ;# no for water box, yes for membrane langevinPiston on langevinPistonTarget 1.01325 ;# in bar -> 1 atm langevinPistonPeriod 200. langevinPistonDecay 50. langevinPistonTemp $temperature } restartfreq 1000 ;# 1000steps = every 2ps dcdfreq 1000 xstFreq 1000 outputEnergies 50 outputPressure 50 # Fixed Atoms Constraint (set PDB beta-column to 1) if {1} { fixedAtoms on fixedAtomsFile ../02-MEMBRANE/myfixatomfile.pdb fixedAtomsCol B fixedAtomsForces on } ############################################################# ## EXTRA PARAMETERS ## ############################################################# # Put here any custom parameters that are specific to # this job (e.g., SMD, TclForces, etc...) ############################################################# ## EXECUTION SCRIPT ## ############################################################# # Minimization if {1} { minimize 10000 reinitvels $temperature } run 500000 ;# 0.5 ns |
回复此楼» 猜你喜欢
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
2楼2011-03-18 19:50:38
3楼2011-03-18 22:57:40
wally8962
木虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 5366.7
- 红花: 2
- 帖子: 1062
- 在线: 471.1小时
- 虫号: 600024
- 注册: 2008-09-12
- 性别: GG
- 专业: 生物力学
4楼2011-03-19 10:06:27
5楼2011-03-19 11:02:10
★ ★ ★
zh1987hs(金币+3): 鼓励思考 2011-03-20 22:50:09
zh1987hs(金币+3): 鼓励思考 2011-03-20 22:50:09
|
这个问题算是暂时解决了，运行后看看吧： 1. 打开压力 2. 在运行Melting of Lipid Tails前，应执行： set sel [atomselect top "protein and backbone"] set all [atomselect top all] $all set beta 0 $sel set beta 1 $all writepdb myfixatomfile.pdb 产生的myfixatomfile.pdb放在conf文件fixedAtomsFile 后面 [ Last edited by nufang19a on 2011-3-19 at 22:56 ] |
6楼2011-03-19 21:30:33