Re: "atoms moving too fast" error in bilayer simulation

From: Peter Freddolino (petefred_at_ks.uiuc.edu)
Date: Fri Feb 02 2007 - 09:29:44 CST

Hi Himanshu,
is this using a home-built lipid bilayer? The bilayers that come with
VMD's membrane plugin should not have any bad contacts...

What is the last gradient tolerance reported by NAMD during the
minimization process? If it's much larger than, say, 1, you may just not
be minimizing the system for long enough to get rid of all the clashes,
so the system explodes as soon as you allow it to move freely. Have you
looked to see if the bad contacts appear to be eliminated at the end of
minimization?

Peter

Himanshu Khandelia wrote:
> Hi,
>
> I am using VMD's membrane plugin to construct a membrane, solvating it,
> and then trying to equilibrate it. The membrane generated by VMD does have
> a lot of bad contacts, but I thought NAMD would handle them.
>
> I keep getting the "atoms moving too fast" error as soon as I start
> heating the bilayer after minimzation. I have tried the following, but
> nothing works so far:
>
> - heating in very small increments (~ 3 degrees every 3000 steps)
> - larger cutoffs (16 A)
> - with/without constraints on head groups
> - larger box sizes (than suggested by the minmax of water)
>
> Any suggestions into what might be going wrong ? The config file is
> attached with this email,
>
> Thank you !
>
> -Himanshu
> ------------------------------------------------------------------------
>
> set MOL pope-water
>
> set temp 500
>
> set xbox 70
> set ybox 70
> set zbox 94
>
> structure ${MOL}.psf
> coordinates ${MOL}.pdb
> temperature 0
> parameters par_all27_prot_lipid.prm
> paraTypeCharmm on
>
> outputEnergies 100
> outputTiming 100
> xstFreq 100
> dcdFreq 500
> wrapAll on
>
> timestep 1
> nonBondedFreq 2
> fullElectFrequency 4
> stepsPerCycle 20
>
> exclude scaled1-4
> 1-4scaling 1.0
> switching on
> switchDist 14
> cutoff 16
> pairlistdist 16
>
> cellBasisVector1 $xbox 00.00 00.00
> cellBasisVector2 00.00 $ybox 00.00
> cellBasisVector3 00.00 00.00 $zbox
> cellOrigin -0.58 -0.78 0.056
>
> Pme on
> PmeGridsizeX 60
> PmeGridsizeY 60
> PmeGridsizeZ 96
>
> #########################################
> # CONSTRAINTS
> #########################################
> constraints on
> consRef ${MOL}-cons.pdb
> consKFile ${MOL}-cons.pdb
> consKCol B
>
> #########################################
> # PRESSURE AND TEMPERATURE CONTRON
> #########################################
>
> langevin on
> langevinDamping 5
> langevinTemp $temp
> langevinHydrogen no
>
> useflexiblecell yes
> langevinPiston on
> langevinPistonTarget 1.01325
> langevinPistonPeriod 100
> langevinPistonDecay 50
> langevinPistonTemp $temp
>
> binaryoutput off
> outputname ${MOL}-equil
>
> # GET INIT. ENERGY
> minimize 0
>
> ###########
> # MINIMIZE
> ###########
>
> minimize 10000
> output ${mol}-cons
>
> #######
> # HEAT
> #######
>
> set tem 1
> while { $tem <= 500 } {
> langevinTemp $tem
> run 3000
> set tem [expr $tem + 3]
> }
>

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:45:59 CST