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 ...
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
- Page 1: Mesoscopic models of lipid bilayers
- Page 4 and 5: Promotiecommissie: Promotor: • pr
- Page 6 and 7: ii CONTENTS 5.3 Double-tail lipid b
- Page 8 and 9: 2 Introduction 1.1 The cell membran
- Page 10 and 11: 4 Introduction between the beads, a
- Page 12 and 13: 6 Introduction whether the preferre
- Page 14 and 15: 8 Simulation method for coarse-grai
- Page 16 and 17: 10 Simulation method for coarse-gra
- Page 18 and 19: 12 Simulation method for coarse-gra
- Page 20 and 21: 14 Simulation method for coarse-gra
- Page 22 and 23: 16 Simulation method for coarse-gra
- Page 24 and 25: 18 Simulation method for coarse-gra
- Page 27 and 28: III Surface tension in lipid bilaye
- Page 29 and 30: 3.2 Method of calculation of surfac
- Page 31: 3.2 Method of calculation of surfac
- Page 35 and 36: 3.3 Constant surface tension ensemb
- Page 37 and 38: 3.4 Surface tension in lipid bilaye
- Page 39 and 40: 3.4 Surface tension in lipid bilaye
- Page 41 and 42: 3.4 Surface tension in lipid bilaye
- Page 43 and 44: IV Structural characterization of l
- Page 45 and 46: 4.2 Structural quantities 39 been r
- Page 47 and 48: 4.3 Computational details 41 lipid
- Page 49 and 50: 4.4 Results and discussion 43 ρ(z)
- Page 51 and 52: 4.4 Results and discussion 45 one l
- Page 53 and 54: 4.4 Results and discussion 47 WH HT
- Page 55 and 56: 4.4 Results and discussion 49 Shill
- Page 57 and 58: 4.4 Results and discussion 51 chain
- Page 59 and 60: 4.4 Results and discussion 53 4.4.4
- Page 61: 4.4 Results and discussion 55 headg
- Page 64 and 65: 58 Phase behavior of coarse-grained
- Page 66 and 67: 60 Phase behavior of coarse-grained
- Page 68 and 69: 62 Phase behavior of coarse-grained
- Page 70 and 71: 64 Phase behavior of coarse-grained
- Page 72 and 73: 66 Phase behavior of coarse-grained
- Page 74 and 75: 68 Phase behavior of coarse-grained
- Page 76 and 77: 70 Phase behavior of coarse-grained
- Page 78 and 79: 72 Phase behavior of coarse-grained
- Page 80 and 81: 74 Phase behavior of coarse-grained
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