From: Axel Kohlmeyer (akohlmey_at_gmail.com)
Date: Wed Oct 28 2009 - 10:55:16 CDT
On Wed, 2009-10-28 at 08:57 +0800, pxq 庞雪芹 wrote:
> namd-lHello everyone,
>
> I am trying to do a MD simulation of POPC lipid membrane according to the reference. However, I am getting a strange results here.
according to _what_ reference?
> When doing the npt, the water moleculars shrink allow the Z axes while the surface of POPC lipid change a little.
this doesn't sound strange at all. this is equilibration.
> When doing a nvt,the water moleculars get an abnormal surface ( it seams the eight angles were cutoff ) still the surface of popc lipid again change a little.
same thing here. your system seem to have a "safety vacuum" between the
lipids and the water and the bubble just moves into the corner(s). look
at the system with periodic images and you should see.
> Besides, I made this water box wtih VMD and simulate it with NAMD, FF Amber 99SB and GAFF.
>
> Following are the scripts of the simulation, attantched are the pictures of both npt and nvt output coordinate.
there are no attachments.
> Waiting for your replay and suggestions is greatly welcomed.
please note replay is for sports, reply for email.
cheers,
axel.
>
> ###########################################################################
> ## JOB DESCRIPTION NPT SIMULATION ##
> ###########################################################################
>
> #output filed
>
> # preparation for MD simulation
> set M 0
> set N 1-p
> set temperature 300
> set begin 0
>
> #############################################################
> ## Required parameters for each simulation ##
> #############################################################
>
> #amber force field parameters
> amber yes
> parmfile 3emlpopccww-cio2.prmtop
> ambercoor 3emlpopccww-cio2.inpcrd
> outputname gpcr$N
> readexclusions no
>
> # frequencies
> outputEnergies 1000
> restartfreq 1000
>
>
> # Force-Field Parameters
> exclude scaled1-4
> 1-4scaling 0.833333
> #according to amber
> cutoff 12.
> switching on
> switchdist 10.
> pairlistdist 13.5
>
> # Periodic Boundary conditions
> if {1} {
> #change according to volume change with lipid membrane
> cellBasisVector1 108. 0. 0.
> cellBasisVector2 0. 122. 0.
> cellBasisVector3 0. 0 113.
> cellOrigin 0.1 -0.3 4.3
> }
> wrapWater on
> wrapAll on
>
>
> #PME (for full-system periodic electrostatics)
> if {1} {
> PME yes
> PMEGridSizeX 120
> PMEGridSizeY 128
> PMEGridSizeZ 120
> #PMEGridSpacing 1
> }
>
> ###################################################
> #################for MD simulations################
> ###################################################
>
>
> if {1} {
>
> #
> timestep 1
> rigidBonds all ;# needed for 2fs steps
> nonbondedFreq 1
> fullElectFrequency 2
> stepspercycle 10
> rigidDieOnError on
>
> #
> firsttimestep $begin
>
> #
> temperature $temperature
>
> #
> dcdfile gpcr$N.dcd
> dcdfreq 1000
> }
>
>
> # Continuing a job from the restart files
> if {0} {
> set inputname gpcr$M
> binCoordinates $inputname.restart.coor
> binVelocities $inputname.restart.vel ;
> extendedSystem $inputname.restart.xsc
> }
>
>
> # Constant Temperature Control
> if {1} {
> langevin on ;
> langevinDamping 5 ;
> langevinTemp $temperature
> langevinHydrogen no ;
> }
>
>
> # Constant Pressure Control (variable volume)
> if {1} {
> useGroupPressure yes ;
> useFlexibleCell yes ;
> useConstantRatio yes
> useConstantArea no ;
>
> langevinPiston on
> langevinPistonTarget 1.01325 ;
> langevinPistonPeriod 200.
> langevinPistonDecay 100.
> langevinPistonTemp $temperature
> SurfaceTensionTarget 60.0 ;
> }
>
> #############################################################
> ## EXECUTION SCRIPT ##
> #############################################################
>
> # Minimization
> if {1} {
> minimize 10000
> reinitvels $temperature
> }
>
> run 2000000 ;# 50ps
>
> >
> >
> >
> >
> >
>
> ###########################################################################
> ## JOB DESCRIPTION NPT SIMULATION ##
> ###########################################################################
>
> #output filed
>
> # preparation for MD simulation
> set M 0
> set N 1-p
> set temperature 300
> set begin 0
>
> #############################################################
> ## Required parameters for each simulation ##
> #############################################################
>
> #amber force field parameters
> amber yes
> parmfile 3emlpopccww-cio2.prmtop
> ambercoor 3emlpopccww-cio2.inpcrd
> outputname gpcr$N
> readexclusions no
>
> # frequencies
> outputEnergies 1000
> restartfreq 1000
>
>
> # Force-Field Parameters
> exclude scaled1-4
> 1-4scaling 0.833333
> #according to amber
> cutoff 12.
> switching on
> switchdist 10.
> pairlistdist 13.5
>
> # Periodic Boundary conditions
> if {1} {
> #change according to volume change with lipid membrane
> cellBasisVector1 108. 0. 0.
> cellBasisVector2 0. 122. 0.
> cellBasisVector3 0. 0 113.
> cellOrigin 0.1 -0.3 4.3
> }
> wrapWater on
> wrapAll on
>
>
> #PME (for full-system periodic electrostatics)
> if {1} {
> PME yes
> PMEGridSizeX 120
> PMEGridSizeY 128
> PMEGridSizeZ 120
> #PMEGridSpacing 1
> }
>
> ###################################################
> #################for MD simulations################
> ###################################################
>
>
> if {1} {
>
> #
> timestep 1
> rigidBonds all ;# needed for 2fs steps
> nonbondedFreq 1
> fullElectFrequency 2
> stepspercycle 10
> rigidDieOnError on
>
> #
> firsttimestep $begin
>
> #
> temperature $temperature
>
> #
> dcdfile gpcr$N.dcd
> dcdfreq 1000
> }
>
>
> # Continuing a job from the restart files
> if {0} {
> set inputname gpcr$M
> binCoordinates $inputname.restart.coor
> binVelocities $inputname.restart.vel ;
> extendedSystem $inputname.restart.xsc
> }
>
>
> # Constant Temperature Control
> if {1} {
> langevin on ;
> langevinDamping 5 ;
> langevinTemp $temperature
> langevinHydrogen no ;
> }
>
>
> # Constant Pressure Control (variable volume)
> if {0} {
> useGroupPressure yes ;
> useFlexibleCell yes ;
> useConstantRatio yes
> useConstantArea no ;
>
> langevinPiston on
> langevinPistonTarget 1.01325 ;
> langevinPistonPeriod 200.
> langevinPistonDecay 100.
> langevinPistonTemp $temperature
> SurfaceTensionTarget 60.0 ;
> }
>
> #############################################################
> ## EXECUTION SCRIPT ##
> #############################################################
>
> # Minimization
> if {1} {
> minimize 10000
> reinitvels $temperature
> }
>
> run 2000000 ;# 50ps
>
>
>
> ===================================================================================================
> Pang Xueqin
> State Key Laboratory of Molecular Reaction Dynamics
> Dalian Institute of Chemical Physics
> Chinese Academy of Sciences.
> Tel: +86-411-84379352 Fax: +86-411-84675584
> ===================================================================================================
>
>
>
>
>
-- Dr. Axel Kohlmeyer akohlmey_at_gmail.com Institute for Computational Molecular Science College of Science and Technology Temple University, Philadelphia PA, USA.
This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:53:25 CST