Re: Membrane Simulation

From: Jim Phillips (jim_at_ks.uiuc.edu)
Date: Wed Oct 31 2018 - 12:15:59 CDT

Hi Alice,

To deal with the layer of vacuum in your system you probably need a new
"useConstantHeight" option that is similar to useConstantArea. The
implementation of useConstantArea is pretty trivial if you search for it
in the source code, so you could potentially add this yourself. Then you
could control surface tension without worrying about the z dimension.

I suspect that if you tried to run a water-only simulation that the
surface tension of the water would collapse the periodic cell until the
vacuum was gone, so having the water-vacuum interface in your simulation
is going to provide a false measurement of the surface tension of the
membrane. I'd suggest running two membranes with water between them to
avoid this problem.

Jim

On Tue, 30 Oct 2018, Alice Zheng wrote:

> Dear all,
>
> I am trying to simulate a membrane bilayer on top of water. In the
> experiment, there is helium gas on top of the bilayer. If I use PME with
> periodic boundary conditions, I have to leave vacuum on top of water. When
> I use NVT for a larger system, the membrane becomes corrugated due to the
> lack of a surface tension force applied to the membrane. Changing
> SurfaceTensionTarget does not prevent the membrane from being corrugated.
>
> Is there a way to apply a surface tension force while using NVT? Or is
> there a way to use NPT for the system? Are there other ensembles that are
> appropriate for this type of system with bilayer lipid membrane situated
> between a aqueous and a gas phase?
>
> I also tried using MSM and tried to eliminate the periodicity in the z
> direction (perpendicular to the membrane). However, I ended up with an
> error message saying that a atom is outside the MSM grid when I did energy
> minimization. The periodicity in x and y could also not be displayed in VMD.
>
> Is it possible to have only 2D periodicity using PME? (Whenever I commented
> out cellBasisVector3, I ended up with the error message that "PME requires
> boundary conditions.")
>
> Any help would be greatly appreciated.
>
> Thank you,
> Alice
>

This archive was generated by hypermail 2.1.6 : Mon Dec 31 2018 - 23:21:29 CST