Mesoscopic models of lipid bilayers and bilayers with embedded ...

Mesoscopic models of lipid bilayers and bilayers with embedded ... Mesoscopic models of lipid bilayers and bilayers with embedded ...

science.uva.nl
from science.uva.nl More from this publisher
30.08.2013 Views

26 Surface tension in lipid bilayers tion 1/No from a given pair (i, j), where No is the total number of slabs that intersect the line of centers between particle i and j. Periodic boundary conditions are always taken into account in this calculation. The integral in equation 3.5 now becomes Ns Ns γ = [PN(k) − PL(k)]∆zs = γ(k) (3.12) k=1 where ∆zs = Lz/Ns is the (uniform) width of the slabs and γ(k) is the local surface tension in slab k. To illustrate the application of this method with a very simple example, we plot in figure 3.1 the distribution of the local surface tension γ(z) across an oil/water interface (the interface is perpendicular to the z-axis), calculated from a DPD simulation of oil- and water-like phase separated particles. In the same figure we also plot the oil, water and total densities profiles to better characterize the shape of the pressure profile. Note that, because of periodic boundary conditions, the interfaces in the simulation box are actually two. The maximum of the surface tension is at the 0.6 0.4 0.2 0.0 −0.2 γ(z) oil water k=1 0 2 4 6 8 10 Z Figure 3.1: Surface tension profile γ(z) (solid black), and oil (dashed black), water (dot-dashed black), and total (solid gray) density profiles across the interfaces (perpendicular to the z-axis) between phase separated oil and water. Because of periodic boundary conditions, there are two interfaces. Note that the densities have been rescaled and shifted for better comparison with the surface tension (hence the scale in the graph ordinate refers to the surface tension only). oil/water interfaces where there are large repulsive forces between oil and water, and the density is lower than in the bulk phases. Moving away from the interface (in both directions since the system is fully symmetric), as a consequence of the repulsive forces at the interface, there is a thin region of compressed fluid where the surface tension becomes (slightly) negative. The surface tension becomes again positive oil

3.3 Constant surface tension ensemble 27 where a second density minimum occurs and then goes to zero in the bulk phases. In the literature, when describing the distribution of pressures in a lipid bilayer, the quantity that is usually reported is the difference between the lateral and the normal pressure. Hence, we find it convenient to define here the pressure profile π(z) as π(z) = [PL(z) − PN(z)]. (3.13) 3.3 Constant surface tension ensemble In this section we introduce a Monte-Carlo (MC) scheme to impose a constant surface tension in the presence of a planar interface [75]. Consider a system with constant number of particles N, constant temperature T, and constant volume V, in which an interface of area A is present. The interface gives an additional term in the energy of the system, i.e. the energy associated with the creation of the interface, which is expressed by the surface tension γ times the area of the interface. The work done on the system by compressing or stretching the interface by dA, is given by [53] dW = γdA. The partition function for such a system can be written as Q = 1 Λ3N dr N! V N exp −β(U(r N ) − γA) . (3.14) where U denotes the potential energy, γ the surface tension, A the area of the interface, and β = 1/kBT. L⊥ z ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ ¡ x L|| A y L|| Figure 3.2: Schematic representation of a simulation box for a system with a flat interface parallel to the xy-plane. The area of the interface is A = L 2 and the box dimension perpendicular to the interface (z axis) is L⊥. Consider a simulation box (see figure 3.2) with edges Lx = Ly = L parallel to

26 Surface tension in <strong>lipid</strong> <strong>bilayers</strong><br />

tion 1/No from a given pair (i, j), where No is the total number <strong>of</strong> slabs that intersect<br />

the line <strong>of</strong> centers between particle i <strong>and</strong> j. Periodic boundary conditions are always<br />

taken into account in this calculation. The integral in equation 3.5 now becomes<br />

Ns<br />

Ns<br />

<br />

<br />

γ = [PN(k) − PL(k)]∆zs = γ(k) (3.12)<br />

k=1<br />

where ∆zs = Lz/Ns is the (uniform) width <strong>of</strong> the slabs <strong>and</strong> γ(k) is the local surface<br />

tension in slab k.<br />

To illustrate the application <strong>of</strong> this method <strong>with</strong> a very simple example, we plot in<br />

figure 3.1 the distribution <strong>of</strong> the local surface tension γ(z) across an oil/water interface<br />

(the interface is perpendicular to the z-axis), calculated from a DPD simulation<br />

<strong>of</strong> oil- <strong>and</strong> water-like phase separated particles. In the same figure we also plot the<br />

oil, water <strong>and</strong> total densities pr<strong>of</strong>iles to better characterize the shape <strong>of</strong> the pressure<br />

pr<strong>of</strong>ile. Note that, because <strong>of</strong> periodic boundary conditions, the interfaces in<br />

the simulation box are actually two. The maximum <strong>of</strong> the surface tension is at the<br />

0.6<br />

0.4<br />

0.2<br />

0.0<br />

−0.2<br />

γ(z)<br />

oil water<br />

k=1<br />

0 2 4 6 8 10<br />

Z<br />

Figure 3.1: Surface tension pr<strong>of</strong>ile γ(z) (solid black), <strong>and</strong> oil (dashed black), water (dot-dashed<br />

black), <strong>and</strong> total (solid gray) density pr<strong>of</strong>iles across the interfaces (perpendicular to the z-axis)<br />

between phase separated oil <strong>and</strong> water. Because <strong>of</strong> periodic boundary conditions, there are<br />

two interfaces. Note that the densities have been rescaled <strong>and</strong> shifted for better comparison<br />

<strong>with</strong> the surface tension (hence the scale in the graph ordinate refers to the surface tension<br />

only).<br />

oil/water interfaces where there are large repulsive forces between oil <strong>and</strong> water, <strong>and</strong><br />

the density is lower than in the bulk phases. Moving away from the interface (in<br />

both directions since the system is fully symmetric), as a consequence <strong>of</strong> the repulsive<br />

forces at the interface, there is a thin region <strong>of</strong> compressed fluid where the surface<br />

tension becomes (slightly) negative. The surface tension becomes again positive<br />

oil

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!