[gmx-users] seperating bilayers (again)
Alan Dodd
anoddlad at yahoo.com
Thu Apr 13 13:15:45 CEST 2006
15, 20, even 30 nm all result in the same effect, just
at different rates. I didn't try anything larger than
that - several of my systems would be dangerously
close to being under 2x the cutoff radius across if I
went even slightly higher.
--- Eric Jakobsson <jake at ncsa.uiuc.edu> wrote:
> What is your van der Waals cut-off distance? From
> our experience
> years ago, we sometimes saw this if our van der
> Waals cut-off was too
> short; i.e., less than 15 angstroms. That seems
> long, since it falls
> off as sixth power, but because the long-range van
> der Waals is all
> attractive, there is no shielding (unlike the
> electrostatic force).
>
> At 11:03 AM 4/6/2006, you wrote:
> >Hello all,
> >I've been having some trouble recently with the
> >leaflets in my bilayer gently drifting apart over a
> >period of a few hundred picoseconds. Increasing
> the
> >vdW radius hasn't helped (up to 2nm), although it
> did
> >delay the seperation. Likewise, increasing tau_p
> from
> >5 up to 10. The really bizarre thing is that I've
> >been simulating peptides in bilayers using 3.2.1 on
> >our cluster for about 18 months now, using the
> >lipid-modified version of ffgmx and the parameters
> >listed below, and it's only since I've started
> using
> >3.3 on a new machine that I've had any trouble.
> Any
> >suggestions on what might be going wrong?
> >
> >title = Yo
> >cpp = /lib/cpp
> >constraints = all-bonds
> >integrator = md
> >dt = 0.0015 ; ps !
> >nsteps = 5000000 ; total 7500 ps.
> >nstcomm = 1
> >nstxout = 1000
> >nstvout = 1000
> >nstfout = 0
> >nstlog = 1000
> >nstenergy = 1000
> >nstlist = 10
> >ns_type = grid
> >coulombtype = PME
> >rlist = 1.5
> >rcoulomb = 1.5
> >rvdw = 1.5
> >pme_order = 8
> >fourierspacing = 0.2
> >optimize_fft = yes
> >; Berendsen temperature coupling is on in two
> groups
> >Tcoupl = berendsen
> >tc-grps = DOPC SOL Protein
> Cl
> >tau_t = 0.1 0.1 0.1
> 0.1
> >ref_t = 300 300 300
> 300
> >; Energy monitoring
> >energygrps = DOPC SOL Protein
> >; Isotropic pressure coupling is now on
> >Pcoupl = berendsen
> >Pcoupltype = surface-tension
> >tau_p = 5
> >compressibility = 4.5e-5 4.5e-5
> >ref_p = 410 1.0
> >; Generate velocites is off at 300 K.
> >gen_vel = no
> >gen_temp = 300.0
> >gen_seed = 173529
> >
> >
> >__________________________________________________
> >Do You Yahoo!?
> >Tired of spam? Yahoo! Mail has the best spam
> protection around
> >http://mail.yahoo.com
> >_______________________________________________
> >gmx-users mailing list gmx-users at gromacs.org
> >http://www.gromacs.org/mailman/listinfo/gmx-users
> >Please don't post (un)subscribe requests to the
> list. Use the
> >www interface or send it to
> gmx-users-request at gromacs.org.
> >Can't post? Read
> http://www.gromacs.org/mailing_lists/users.php
>
> ---------------------------------
> Eric Jakobsson, Ph.D.
> Professor, Department of Molecular and Integrative
> Physiology, and of
> Biochemistry, and of the Center for Biophysics and
> Computational Biology
> Senior Research Scientist, National Center for
> Supercomputing Applications
> Professor, Beckman Institute for Advanced Science
> and Technology
> 3261 Beckman Institute, mc251
> University of Illinois, Urbana, IL 61801
> ph. 217-244-2896 fax 217-244-2909
>
>
>
>
> _______________________________________________
> gmx-users mailing list gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the
> list. Use the
> www interface or send it to
> gmx-users-request at gromacs.org.
> Can't post? Read
> http://www.gromacs.org/mailing_lists/users.php
>
__________________________________________________
Do You Yahoo!?
Tired of spam? Yahoo! Mail has the best spam protection around
http://mail.yahoo.com
More information about the gromacs.org_gmx-users
mailing list