Phase transitions in a fluid surface model with a deficit angle term

Nambu-Goto model is investigated by using the canonical Monte Carlo simulation technique on dynamically triangulated surfaces of spherical topology. We find that the model has four distinct phases; crumpled, branched-polymer, linear, and tubular. The linear phase and the tubular phase appear to be separated by a first-order transition. It is also found that there is no long-range two-dimensional order in the model. In fact, no smooth surface can be seen in the whole region of the curvature modulus α, which is the coefficient of the deficit angle term in the Hamiltonian. The bending energy, which is not included in the Hamiltonian, remains large even at sufficiently large α in the tubular phase. On the other hand, the surface is spontaneously compactified into a one-dimensional smooth curve in the linear phase; one of the two degrees of freedom shrinks, and the other degree of freedom remains along the curve. Moreover, we find that the rotational symmetry of the model is spontaneously broken in the tubular phase just as in the same model on the fixed connectivity surfaces.     

Phase transition of compartmentalized surface models

Two types of surface models have been investigated by Monte Carlo simulations on triangulated spheres with compartmentalized domains. Both models are found to undergo a first-order collapsing transition and a first-order surface fluctuation transition. The first model is a fluid surface one. The vertices can freely diffuse only inside the compartments, and they are prohibited from the free diffusion over the surface due to the domain boundaries. The second is a skeleton model. The surface shape of the skeleton model is maintained only by the domain boundaries, which are linear chains with rigid junctions. Therefore, we can conclude that the first-order transitions occur independent of whether the shape of surface is mechanically maintained by the skeleton (=the domain boundary) or by the surface itself.     

Shape transformation transitions in a model of fixed-connectivity surfaces supported by skeletons

A compartmentalized surface model of Nambu and Goto is studied on triangulated spherical surfaces by using the canonical Monte Carlo simulation technique. One-dimensional bending energy is defined on the skeletons and at the junctions, and the mechanical strength of the surface is supplied by the one-dimensional bending energy defined on the skeletons and junctions. The compartment size is characterized by the total number L^′ of bonds between the two-neighboring junctions and is assumed to have values in the range from L^′ = 2 to L^′ = 8 in the simulations, while that of the previously reported model is characterized by L^′ = 1, where all vertices of the triangulated surface are the junctions. Therefore, the model in this paper is considered to be an extension of the previous model in the sense that the previous model is obtained from the model in this paper in the limit of L^′↦1. The model in this paper is identical to the Nambu-Goto surface model without curvature energies in the limit of L^′↦∞ and hence is expected to be ill-defined at sufficiently large L^′. One remarkable result obtained in this paper is that the model has a well-defined smooth phase even at relatively large L^′ just as the previous model of L^′↦ 1. It is also remarkable that the fluctuations of surface in the smooth phase are crucially dependent on L^′; we can see no surface fluctuation when L^′≤ 2, while relatively large fluctuations are seen when L^′≥ 3.     

Phase structure of intrinsic curvature models on dynamically triangulated disk with fixed boundary length

A first-order phase transition is found in two types of intrinsic curvature models defined on dynamically triangulated surfaces of disk topology. The intrinsic curvature energy is included in the Hamiltonian. The smooth phase is separated from a non-smooth phase by the transition. The crumpled phase, which is different from the non-smooth phase, also appears at sufficiently small curvature coefficient α. The phase structure of the model on the disk is identical to that of the spherical surface model, which was investigated by us and reported previously. Thus, we found that the phase structure of the fluid surface model with intrinsic curvature is independent of whether the surface is closed or open.     

Phase structure of a surface model with many fine holes

We study the phase structure of a surface model by using the canonical Monte Carlo simulation technique on triangulated, fixed connectivity, and spherical surfaces with many fine holes. The size of a hole is assumed to be of the order of lattice spacing (or bond length) and hence can be negligible compared to the surface size in the thermodynamic limit. We observe in the numerical data that the model undergoes a first-order collapsing transition between the smooth phase and the collapsed phase. Moreover the Hasudorff dimension H remains in the physical bound, i.e., H < 3 not only in the smooth phase but also in the collapsed phase at the transition point. The second observation is that the collapsing transition is accompanied by a continuous transition of surface fluctuations. This second result distinguishes the model in this paper and the previous one with many holes, whose size is of the order of the surface size, because the previous surface model with large-sized holes has only the collapsing transition and no transition of surface fluctuations.     

Non-trivial effect of the in-plane shear elasticity on the phase transitions of fixed-connectivity meshwork models

We numerically study the phase structure of two types of triangulated spherical surface models, which includes an in-plane shear energy in the Hamiltonian, and we found that the phase structure of the models is considerably influenced by the presence of the in-plane shear elasticity. The models undergo a first-order collapsing transition and a first-order (or second-order) transition of surface fluctuations; the latter transition was reported to be of second-order in the first model without the in-plane shear energy. This leads us to conclude that the in-plane elasticity strengthens the transition of surface fluctuations. We also found that the in-plane elasticity decreases the variety of phases in the second model without the in-plane energy. The Hamiltonian of the first model is given by a linear combination of the Gaussian bond potential, a one-dimensional bending energy, and the in-plane shear energy. The second model is obtained from the first model by replacing the Gaussian bond potential with the Nambu-Goto potential, which is defined by the summation over the area of triangles.     

Phase structure of a spherical surface model on fixed connectivity meshes

An elastic surface model is investigated by using the canonical Monte Carlo simulation technique on triangulated spherical meshes. The model undergoes a first-order collapsing transition and a continuous surface fluctuation transition. The shape of surfaces is maintained by a one-dimensional bending energy, which is defined on the mesh, and no two-dimensional bending energy is included in the Hamiltonian.     

Nucleation-and-growth problem in model lipid membranes undergoing subgel phase transitions is a problem of time scale

In this Rapid Note, we show that the problem of growth of molecular superlattice in a fully hydrated dipalmitoylphosphatidylcholine (DPPC) membrane during the gel-to-subgel phase transformation process is a problem of time scale. There are, in fact, two time scales. The first is an “integrated” or, in some sense, stagnant time scale, that reflects the well-known isotropic growth effect in the d-dimensional space, but assigns the problem to be still in a category of Debye relaxation kinetics. The fraction of old (parent) phase does not suit the Paley-Wiener criterion for relaxation functions, and the time behavior is exclusively due to the geometrical characteristics of the kinetic process. The second (multi-instantaneous) time scale, in turn, is recognised to be a “broken” (fractional time derivative) or memory-feeling (dynamic) scale, which carries some very essential physics of the phenomenon under study, and classifies the problem to be of non-Debye (viz., stretched exponential) nature. It may, in principle, contain all the important effects, like small scale coexistence, presence of collisions between domains, with possible annihilation and creation of domain boundaries, and/or a headgroup packing, hydration against lipid mobility behavior, and finally, a multitude of quasi-crystalline states. It turns out, that within the range of validity of the dynamic scale approximation proposed, the criterion for relaxation functions is very well fulfilled. Received 30 November 1998     

Elasticity and shape equation of a liquid membrane

The density of the elastic energy of a deformed membrane in a liquid state is calculated. The thermodynamic equilibrium of its different parts is taken into account. The shape equation of a closed membrane is deduced. The quantity which keeps its value, when the variations of the energy of the system are calculated, is not the area of the deformed membrane, but its area in the flat tension free state. Because of this, additional terms appear in the second variation around the stable state. The case of a lipid bilayer and its fluctuations is examined for both free and blocked exchange of molecules between the monolayers, comprising the bilayer. Received 4 February 2002 / Received in final form 15 April 2002 Published online 2 October 2002 RID="a" ID="a"e-mail: bivas@issp.bas.bg     

Pre-roughening dynamics on Si(100) stepped surfaces: a study by molecular dynamics

The thermal evolution of steps on Si(100) is well studied and experiment indicates that at temperatures below the roughening transition (i.e. T? 1000 K) the displacements of atoms at the step-edge are the basic factor of this evolution. However the evaluation of the nature and participants of these displacements is beyond experimental observations and a theoretical approach is therefore needed. The problem addressed by this study is the identification of the properties of atomic motions of step-edge atoms and this investigation is performed applying an isothermal Molecular Dynamics simulation method to simple stepped configurations on Si(100). The calculations describe the functional dependence of the motions of step-edge atoms on the step type, size and temperature and on the nature of the interatomic forces. Possible mechanisms of kink formations are suggested. Received 15 February 2002 Published online 13 August 2002     

Adhesion of fluid vesicles at chemically structured substrates

The adhesion of fluid vesicles at chemically structured substrates is studied theoretically via Monte Carlo simulations. The substrate surface is planar and repels the vesicle membrane apart from a single surface domain γ , which strongly attracts this membrane. If the vesicle is larger than the attractive γ domain, the spreading of the vesicle onto the substrate is restricted by the size of this surface domain. Once the contact line of the adhering vesicle has reached the boundaries of the γ domain, further deflation of the vesicle leads to a regime of low membrane tension with pronounced shape fluctuations, which are now governed by the bending rigidity. For a circular γ domain and a small bending rigidity, the membrane oscillates strongly around an average spherical cap shape. If such a vesicle is deflated, the contact area increases or decreases with increasing osmotic pressure, depending on the relative size of the vesicle and the circular γ domain. The lateral localization of the vesicle's center of mass by such a domain is optimal for a certain domain radius, which is found to be rather independent of adhesion strength and bending rigidity. For vesicles adhering to stripe-shaped surface domains, the width of the contact area perpendicular to the stripe varies nonmonotonically with the adhesion strength.     

Vesicle self-reproduction: The involvement of membrane hydraulic and solute permeabilities

Conditions for self-reproduction are sought for a growing vesicle with its growth defined by an exponential increase of vesicle membrane area and by adequate flow of the solution across the membrane. In the first step of the presumed vesicle self-reproduction process, the initially spherical vesicle must double its volume in the doubling time of the membrane area and, through the appropriate shape transformations, attain the shape of two equal spheres connected by an infinitesimally thin neck. The second step involves separation of the two spheres and relies on conditions that cause the neck to be broken. In this paper we consider the first step of this self-reproduction process for a vesicle suspended in a solution whose solute can permeate the vesicle membrane. It is shown that vesicle self-reproduction occurs only for certain combinations of the values of membrane hydraulic and solute permeabilities and the external solute concentration, these quantities being related to the mechanical properties of the membrane and the membrane area doubling time. The analysis includes also the relaxation of a perturbed system towards stationary self-reproduction behavior and the case where the final shape consists of two connected spheres of different radii.     

On coupling between the orientation and the shape of a vesicle under a shear flow

A simple 2D model of deformable vesicles tumbling in a shear under flow is introduced in order to account for the main qualitative features observed experimentally as shear rates are increased. The simplicity of the model allows for a full analytical tractability while retaining the essential physical ingredients. The model reveals that the main axes of the vesicle undergo oscillations which are coupled to the vesicle orientation in the flow. The model reproduces and sheds light on the main novel features reported in recent experiments [M. Mader et al., Eur. Phys. J. E. 19, 389 (2006)], namely that both coefficients A and B that enter the Keller-Skalak equation, dψ/dt = A+Bcos(2 ψ) (ψ is the vesicle orientation angle in the shear flow), undergo a collapse upon increasing shear rate.     

Domain-induced budding in buckling membranes

We present a phase field model on buckling membranes to analyze phase separation and budding on soft membranes. By numerically integrating dynamic equations, it turns out that the formation of caps is greatly influenced by the presence of a little excess area due to the surface area constraint. When cap-shaped domains are created, domain coalescence is mainly observed not between domains with same budding directions, but between domains with opposite budding directions, because the bending energy between two domains is larger in the former case. Although we do not introduce spontaneous curvature like Helfrich model, we obtain some suggestions related to the slow dynamics of the phase separation on vesicles.     

Dynamic model and stationary shapes of fluid vesicles

A phase-field model that takes into account the bending energy of fluid vesicles is presented. The Canham-Helfrich model is derived in the sharp-interface limit. A dynamic equation for the phase-field has been solved numerically to find stationary shapes of vesicles with different topologies and the dynamic evolution towards them. The results are in agreement with those found by minimization of the Canham-Helfrich free energy. This fact shows that our phase-field model could be applied to more complex problems of instabilities.     

Effect of an electric field on a floating lipid bilayer: A neutron reflectivity study

We present here a neutron reflectivity study of the influence of an alternative electric field on a supported phospholipid double bilayer. We report for the first time a reproducible increase of the fluctuation amplitude leading to the complete unbinding of the floating bilayer. Results are in good agreement with a semi-quantitative interpretation in terms of negative electrostatic surface tension.     

Dynamics of viscous vesicles in shear flow

The dynamics of giant lipid vesicles under shear flow is experimentally investigated. Consistent with previous theoretical and numerical studies, two flow regimes are identified depending on the viscosity ratio between the interior and the exterior of the vesicle, and its reduced volume or excess surface. At low viscosity ratios, a tank-treading motion of the membrane takes place, the vesicle assuming a constant orientation with respect to the flow direction. At higher viscosity ratios, a tumbling motion is observed in which the whole vesicle rotates with a periodically modulated velocity. When the shear rate increases, this tumbling motion becomes increasingly sensitive to vesicle deformation due to the elongational component of the flow and significant deviations from simpler models are observed. A good characterization of these various flow regimes is essential for the validation of analytical and numerical models, and to relate microscopic dynamics to macroscopic rheology of suspensions of deformable particles, such as blood.     

Fluctuation spectra of free and supported membrane pairs

Fluctuation spectra of fluid compound membrane systems are calculated. The systems addressed contain two (or more) almost parallel membranes that are connected by harmonic tethers or by a continuous, harmonic confining potential. Additionally, such a compound system can be attached to a supporting substrate. We compare quasi-analytical results for tethers with analytical results for corresponding continuous models and investigate under what circumstances the discrete nature of the tethers actually influences the fluctuations. A tethered, supported membrane pair with similar bending rigidities and stiff tethers can possess a nonmonotonic fluctuation spectrum with a maximum. A nonmonotonic spectrum with a maximum and a minimum can occur for an either free or supported membrane pair of rather different bending rigidities and for stiff tethers. Typical membrane displacements are calculated for supported membrane pairs with discrete or continuous interacting potentials. Thereby an estimate of how close the constituent two membranes and the substrate typically approach each other is given. For a supported membrane pair with discrete or continuous interactions, the typical displacements of each membrane are altered with respect to a single supported membrane, where those of the membrane near the substrate are diminished and those of the membrane further away are enhanced.