If you add position restraints to your DPPC molecules, then you are
changing your effective order parameter. The PMF for the solute along
the normal to a restrained bilayer will have a different shape than a
PMF for the solute along the normal to an urestrained bilayer. Whether
or not this is a problem will depend on the question that you are trying
to answer.
If the bilayer is falling apart, then you certainly need to do
something. But if it is just locally ordering/disordering, then I don't
see the big problem, besides the effect that this has on convergence
(DOI: 10.1021/ct200316w).
Chris.
-- original message --
Dear Gromacs Users,
I am trying to extract the potential of mean force of a small
molecule in a DPPC bilayer. To this end, I applied the methodology
described in an online manual written by Justin Lemkul. My problem is
when I run biasing simulations of the molecule near the interface
(DPPC/water), some lipid molecules move to the water phase. This has as
a consequence a local disorder of the bilayer. Below is the parameters I
employ for the pull code:
pull = umbrella
pull_geometry = distance
pull_dim = N N Y
pull_start = yes
pull_ngroups = 1
pull_group0 = DPPC
pull_group1 = DTC
pull_rate1 = 0.0
pull_k1 = 1000
Is the position restraints on DPPC molecules a solution to my problem?
Thanks in advance
Best regards
Giovani
--
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
