BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Hiroshi Noguchi

Hiroshi Noguchi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softBiological Physicsphysics.flu-dynmath.NTphysics.comp-ph

Catalog footprint

What is connected

32works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2022arXiv

Binding of anisotropic curvature-inducing proteins onto membrane tubes

Bin/Amphiphysin/Rvs superfamily proteins and other curvature-inducing proteins have anisotropic shapes and anisotropically bend biomembrane. Here, we report how the anisotropic proteins bind the membrane tube and are orientationally ordered using mean-field theory including an orientation-dependent excluded volume. The proteins exhibit a second-order or first-order nematic transition with increasing protein density depending on the radius of the membrane tube. The tube curvatures for the maximum protein binding and orientational order are different and varied by the protein density and rigidity. As the external force along the tube axis increases, a first-order transition from a large tube radius with low protein density to a small radius with high density occurs once, and subsequently, the protein orientation tilts to the tube-axis direction. When an isotropic bending energy is used for the proteins with an elliptic shape, the force-dependence curves become symmetric and the first-order transition occurs twice. This theory quantitatively reproduces the results of meshless membrane simulation for short proteins, whereas deviations are seen for long proteins owing to the formation of protein clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Excitable reaction-diffusion waves of curvature-inducing proteins on deformable membrane tubes

Living cells employ excitable reaction-diffusion waves for internal cellular functions, in which curvature-inducing proteins are often involved. However, the role of their mechanochemical coupling is not well understood. Here, we report the membrane deformation induced by the excitable reaction-diffusion waves of curvature-inducing proteins and the alternation in the waves due to the deformation, using a coarse-grained simulation of tubular membranes with a modified FitzHugh--Nagumo model. Protein-propagating waves deform tubular membranes, and large deformations induce budding and erase waves. The wave speed and shape are determined by a combination of membrane deformation and spatial distribution of the curvature-inducing protein. Waves are also undulated in the azimuthal direction depending on the condition. Rotationally symmetric waves locally deform the tubes into a symmetric shape but maintain a straight shape on average. Our simulation method can be applied to other chemical reaction models and used to investigate various biomembrane phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Membrane shape deformation induced by curvature-inducing proteins consisting of chiral crescent binding and intrinsically disordered domains

Curvature-inducing proteins containing a Bin/Amphiphysin/Rvs domain often have intrinsically disordered domains. Recent experiments have shown that these disordered chains enhance curvature sensing and generation. Here, we report on the modification of protein-membrane interactions by disordered chains using meshless membrane simulations. The protein and bound membrane are modeled together as a chiral crescent protein rod with two excluded-volume chains. As the chain length increases, the repulsion between them reduces the cluster size of the proteins. It induces spindle-shaped vesicles and a transition between arc-shaped and circular protein assemblies in a disk-shaped vesicle. For flat membranes, an intermediate chain length induces many tubules owing to the repulsion between the protein assemblies, whereas longer chains promote perpendicular elongation of tubules. Moreover, protein rods with zero rod curvature and sufficiently long chains stabilize the spherical buds. For proteins with a negative rod curvature, an intermediate chain length induces a rugged membrane with branched protein assemblies, whereas longer chains induce the formation of tubules with periodic concave-ring structures.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Binding of thermalized and active membrane curvature-inducing proteins

Membrane phase behavior induced by the binding of curvature-inducing proteins is studied by a combination of analytical and numerical approaches. In thermal equilibrium under the detailed balance between binding and unbinding, the membrane exhibits three phases: an unbound uniform flat phase (U), a bound uniform flat phase (B), and a separated/corrugated phase (SC). In the SC phase, the bound proteins form hexagonally-ordered bowl-shaped domains. The transitions between the U and SC phases and between the B and SC phases are second order and first order, respectively. At a small spontaneous curvature of the protein or high surface tension, the transition between B and SC phases becomes continuous. Moreover, a first-order transition between the U and B phases is found at zero spontaneous curvature driven by the Casimir-like interactions between rigid proteins. Furthermore, nonequilibrium dynamics is investigated by the addition of active binding and unbinding at a constant rate. The active binding and unbinding processes alter the stability of the SC phase.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On multiplier systems and theta functions of half-integral weight for the Hilbert modular group $\mathrm{SL}_2(\mathfrak{o})$

Let $F$ be a totally real number field and $\mathfrak{o}$ the ring of integers of $F$. We study theta functions which are Hilbert modular forms of half-integral weight for the Hilbert modular group $\mathrm{SL}_2(\mathfrak{o})$. We obtain an equivalent condition that there exists a multiplier system of half-integral weight for $\mathrm{SL}_2(\mathfrak{o})$. We determine the condition of $F$ that there exists a theta function which is a Hilbert modular form of half-integral weight for $\mathrm{SL}_2(\mathfrak{o})$. The theta function is defined by a sum on a fractional ideal $\mathfrak{a}$ of $F$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Effects of Cavitation on Karman Vortex Behind Circular-Cylinder Arrays

The effects of cavitation on the flow around a circular-cylinder array are studied by using a molecular dynamics simulation. Cavitation significantly affects on vortex shedding characteristics. As the cavitation develops,the vibration acting on the cylinders decreases and eventually disappears. The further cavitation development generates a longer vapor region next to the cylinders, and the vortex streets are formed at further positions from the cylinders. The neighboring Karman vortexes are synchronized in the antiphase in the absence of the cavitation. This synchronization is weakened by the cavitation, and an asymmetric wake mode can be induced. These findings help mechanical designs of fluid machinery that include cylinder arrays.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Construction of nuclear envelope shape by a high-genus vesicle with pore-size constraint

Nuclear pores have an approximately uniform distribution in the nuclear envelope of most living cells. Hence, the morphology of the nuclear envelope is a spherical stomatocyte with a high genus. We have investigated the morphology of high-genus vesicles under pore-size constraint using dynamically triangulated membrane simulations. Bending-energy minimization without volume or other constraints produces a circular-cage stomatocyte, where the pores are aligned in a circular line on an oblate bud. As the pore radius is reduced, the circular-pore alignment is more stabilized than a random pore distribution on a spherical bud. However, we have clarified the conditions for the formation of a spherical stomatocyte: a small perinuclear volume, osmotic pressure within nucleoplasm, and repulsion between the pores. When area-difference elasticity is taken into account, the formation of cylindrical or budded tubules from the stomatocyte and discoidal stomatocyte is found.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Membrane tubule formation by banana-shaped proteins with or without transient network structure

In living cells, membrane morphology is regulated by various proteins. Many membrane reshaping proteins contain a Bin/Amphiphysin/Rvs (BAR) domain, which consists of a banana-shaped rod. The BAR domain bends the biomembrane along the rod axis and the features of this anisotropic bending have recently been studied. Here, we report on the role of the BAR protein rods in inducing membrane tubulation, using large-scale coarse-grained simulations. We reveal that a small spontaneous side curvature perpendicular to the rod can drastically alter the tubulation dynamics at high protein density, whereas no significant difference is obtained at low density. A percolated network is intermediately formed depending on the side curvature. This network suppresses tubule protrusion, leading to the slow formation of fewer tubules. Thus, the side curvature, which is generated by protein--protein and membrane--protein interactions, plays a significant role in tubulation dynamics. We also find that positive surface tensions and the vesicle membrane curvature can stabilize this network structure by suppressing the tubulation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Non-uniqueness of local stress of three-body potentials in molecular simulations

Microscopic stress fields are widely used in molecular simulations to understand mechanical behavior. Recently, decomposition methods of multibody forces to central force pairs between the interacting particles have been proposed. Here, we introduce a force center of a three-body potential and propose different force decompositions that also satisfy the conservation of translational and angular momentum. We compare the force decompositions by stress-distribution magnitude and discuss their difference in the stress profile of a bilayer membrane using coarse-grained and atomistic molecular dynamics simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Shape deformation of lipid membranes by banana-shaped protein rods: Comparison with isotropic inclusions and membrane rupture

The assembly of curved protein rods on fluid membranes is studied using implicit-solvent meshless membrane simulations. As the rod curvature increases, the rods on a membrane tube assemble along the azimuthal direction first and subsequently along the longitudinal direction. Here, we show that both transition curvatures decrease with increasing rod stiffness. For comparison, curvature-inducing isotropic inclusions are also simulated. When the isotropic inclusions have the same bending rigidity as the other membrane regions, the inclusions are uniformly distributed on the membrane tubes and vesicles even for large spontaneous curvature of the inclusions. However, the isotropic inclusions with much larger bending rigidity induce shape deformation and are concentrated on the region of a preferred curvature. For high rod density, high rod stiffness, and/or low line tension of the membrane edge, the rod assembly induces vesicle rupture, resulting in the formation of a high-genus vesicle. A gradual change in the curvature suppresses this rupture. Hence, large stress, compared to the edge tension, induced by the rod assembly is the key factor determining rupture. For rod curvature with the opposite sign to the vesicle curvature, membrane rupture induces inversion of the membrane, leading to division into multiple vesicles as well as formation of a high-genus vesicle.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Formation of polyhedral vesicles and polygonal membrane tubes induced by banana-shaped proteins

The shape transformations of fluid membranes induced by curved protein rods are studied using meshless membrane simulations. The rod assembly at low rod density induces a flat membrane tube and oblate vesicle. It is found that the polyhedral shapes are stabilized at high rod densities. The discrete shape transition between triangular and buckled discoidal tubes is obtained and their curvature energies are analyzed by a simple geometric model. For vesicles, triangular hosohedron and elliptic-disk shapes are formed in equilibrium, whereas tetrahedral and triangular prism shapes are obtained as metastable states.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Monte Carlo study of the frame, fluctuation and internal tensions of fluctuating membranes with fixed area

Three types of surface tensions can be defined for lipid membranes: the internal tension, $σ$, conjugated to the real membrane area in the Hamiltonian, the mechanical frame tension, $τ$, conjugated to the projected area, and the "fluctuation tension", $r$, obtained from the fluctuation spectrum of the membrane height. We investigate these surface tensions by means of a Monge gauge lattice Monte Carlo simulation involving the exact, nonlinear, Helfrich Hamiltonian and a measure correction for excess entropy of the Monge gauge. Our results for the relation between $σ$ and $τ$ agrees well with the theoretical prediction of [J.-B. Fournier and C. Barbetta, Phys. Rev. Lett., 2008, 100, 078103] based on a Gaussian approximation. This provides a valuable knowledge of~$τ$ in the standard Gaussian models where the tension is controlled by $σ$. However, contrary to the conjecture in the above paper, we find that $r$ exhibits no significant difference from $τ$ over more than five decades of tension. Our results appear to be valid in the thermodynamic limit and are robust to changing the ensemble in which the membrane area is controlled.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Shape transitions of high-genus fluid vesicles

The morphologies of genus-2 to -8 fluid vesicles are studied by using dynamically triangulated membrane simulations with area-difference elasticity. It is revealed that the alignments of the membrane pores alter the vesicle shapes and the types of shape transitions for the genus $g \ge 3$. At a high reduced volume, a stomatocyte with a circular alignment of $g+1$ pores continuously transforms into a discocyte with a line of $g$ pores with increasing intrinsic area difference. In contrast, at a low volume, a stomatocyte transforms into a ($g+1$)-hedral shape and subsequently exhibits a discrete phase transition to a discocyte.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Morphological changes of amphiphilic molecular assemblies induced by chemical reaction

Shape transformations of amphiphilic molecular assemblies induced by chemical reaction are studied using coarse-grained molecular simulations. A binding reaction between hydrophilic and hydrophobic molecules is considered. It is found that the reaction induces transformation of an oil droplet to a tubular vesicle via bicelles and vesicles with discoidal arms. The discoidal arms close into vesicles, which are subsequently fused into the tubular vesicle. Under the chemical reaction, the bicelle-to-vesicle transition occurs at smaller sizes than in the absence of the hydrophobic molecules. It is revealed that the enhancement of this transition is due to embedded hydrophobic particles that reduce the membrane bending rigidity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Shape transformations of toroidal vesicles

Morphologies of genus-1 and 2 toroidal vesicles are studied numerically by dynamically triangulated membrane models and experimentally by confocal laser microscopy. Our simulation results reproduce shape transformations observed in our experiments well. At large reduced volumes of the genus-1 vesicles, obtained vesicle shapes agree with the previous theoretical prediction, in which axisymmetric shapes are assumed: double-necked stomatocyte, discoidal toroid, and circular toroid. However, for small reduced volumes, it is revealed that a non-axisymmetric discoidal toroid and handled discocyte exist in thermal equilibrium in the parameter range, in which the previous theory predicts axisymmetric discoidal shapes. Polygonal toroidal vesicles and subsequent budding transitions are also found. The entropy caused by shape fluctuations slightly modifies the stability of the vesicle shapes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Two- or three-step assembly of banana-shaped proteins coupled with shape transformation of lipid membranes

BAR superfamily proteins have a banana-shaped domain that causes the local bending of lipid membranes. We study as to how such a local anisotropic curvature induces effective interaction between proteins and changes the global shape of vesicles and membrane tubes using meshless membrane simulations. The proteins are modeled as banana-shaped rods strongly adhered to the membrane. Our study reveals that the rods assemble via two continuous directional phase separations unlike a conventional two-dimensional phase separation. As the rod curvature increases, in the membrane tube the rods assemble along the azimuthal direction and subsequently along the longitudinal direction accompanied by shape transformation of the tube. In the vesicle, in the addition to these two assembly processes, further increase in the rod curvature induces tubular scaffold formation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Confinement-induced shape transitions in multilamellar vesicles

Morphologies of a vesicle confined in a spherical vesicle were explored experimentally by fast confocal laser microscopy and numerically by a dynamically-triangulated membrane model with area-difference elasticity. The confinement was found to induce several novel shapes of the inner vesicles, that had been never observed in unilamellar vesicles: double and quadruple stomatocytes, slit vesicle, and vesicles of two or three compartments with various shapes. The simulations reproduced the experimental results very well and some of the shape transitions can be understood by a simple theoretical model for axisymmetric shapes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Entropy-driven aggregation in multilamellar membranes

Membrane-fluctuation-induced attraction between ligand--receptor sites binding neighboring membranes is studied using meshless membrane simulations and the Weil--Farago two-dimensional lattice model. For the adhesion sites binding two membranes, this entropic interaction is too weak by itself for the adhesion sites to form a large stable domain. However, it is found that this attraction is enhanced sufficiently to induce large domains either when the sites bind three or more neighboring membranes together or have anchors that harden surrounding membranes. The latter effect is understood by the Asakura--Oosawa type of effective potential in the depletion theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Structure formation of surfactant membranes under shear flow

Shear-flow-induced structure formation in surfactant-water mixtures is investigated numerically using a meshless-membrane model in combination with a particle-based hydrodynamics simulation approach for the solvent. At low shear rates, uni-lamellar vesicles and planar lamellae structures are formed at small and large membrane volume fractions, respectively. At high shear rates, lamellar states exhibit an undulation instability, leading to rolled or cylindrical membrane shapes oriented in the flow direction. The spatial symmetry and structure factor of this rolled state agree with those of intermediate states during lamellar-to-onion transition measured by time-resolved scatting experiments. Structural evolution in time exhibits a moderate dependence on the initial condition.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Structure formation in binary mixtures of lipids and detergents: Self-assembly and vesicle division

Self-assembly dynamics in binary surfactant mixtures and structure changes of lipid vesicles induced by detergent solution are studied using coarse-grained molecular simulations. Disk-shaped micelles, the bicelles, are stabilized by detergents surrounding the rim of a bilayer disk of lipids. The self-assembled bicelles are considerably smaller than bicelles formed from vesicle rupture, and their size is determined by the concentrations of lipids and detergents and the interactions between the two species. The detergent-adsorption induces spontaneous curvature of the vesicle bilayer and results in vesicle division into two vesicles or vesicle rupture into worm-like micelles. The division occurs mainly via the inverse pathway of the modified stalk model. For large spontaneous curvature of the monolayers of the detergents, a pore is often opened, thereby leading to vesicle division or worm-like micelle formation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Structure formation in binary mixtures of surfactants: vesicle opening-up to bicelles and octopus-like micelles

Micelle formation in binary mixtures of surfactants is studied using a coarse-grained molecular simulation. When a vesicle composed of lipid and detergent types of molecules is ruptured, a disk-shaped micelle, the bicelle, is typically formed. It is found that cup-shaped vesicles and bicelles connected with worm-like micelles are also formed depending on the surfactant ratio and critical micelle concentration. The obtained octopus shape of micelles agree with those observed in the cryo-TEM images reported in [S. Jain and F. S. Bates, Macromol. 37, 1511 (2004).]. Two types of connection structures between the worm-like micelles and the bicelles are revealed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Anisotropic surface tension of buckled fluid membrane

Solid sheets and fluid membranes exhibit buckling under lateral compression. Here, it is revealed that fluid membranes have anisotropic buckling surface tension contrary to solid sheets. Surprisingly, the surface tension perpendicular to the buckling direction shows stronger dependence than that parallel to it. Our theoretical predictions are supported by numerical simulations of a meshless membrane model. This anisotropic tension can be used to measure the membrane bending rigidity. It is also found phase synchronization occurs between multilayered buckled membranes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Line tension of branching junctions of bilayer membranes

Branching of bilayer membranes appear in the inverted hexagonal phase as well as in metastable states of the lamellar phase such as membrane fusion intermediates. A method for estimating the line tension of the branching junction is proposed for molecular simulations. The line tension is calculated from the pressure tensor of equiangularly branched membranes. The simulation results agree very well with the theoretical prediction of Hamm and Kozlov's tilt model. The transition between the lamellar and inverted hexagonal phases is also investigated using the tilt model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Measurement of the bending rigidity and spontaneous curvature of fluid membranes in simulations

Several numerical methods for measuring the bending rigidity and the spontaneous curvature of fluid membranes are studied using two types of meshless membrane models. The bending rigidity is estimated from the thermal undulations of planar and tubular membranes and the axial force of tubular membranes. We found a large dependence of its estimate value from the thermal undulation analysis on the upper-cutoff frequency q_{cut} of the least squares fit. The inverse power-spectrum fit with an extrapolation to q_{cut} -> 0 yields the smallest estimation error among the investigated methods. The spontaneous curvature is estimated from the axial force of tubular membranes and the average curvature of bent membrane stripes. The results of these methods show good agreement with each other.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamic Modes of Microcapsules in Steady Shear Flow: Effects of Bending and Shear Elasticities

The dynamics of microcapsules in steady shear flow was studied using a theoretical approach based on three variables: The Taylor deformation parameter $α_{\rm D}$, the inclination angle $θ$, and the phase angle $ϕ$ of the membrane rotation. It is found that the dynamic phase diagram shows a remarkable change with an increase in the ratio of the membrane shear and bending elasticities. A fluid vesicle (no shear elasticity) exhibits three dynamic modes: (i) Tank-treading (TT) at low viscosity $η_{\rm {in}}$ of internal fluid ($α_{\rm D}$ and $θ$ relaxes to constant values), (ii) Tumbling (TB) at high $η_{\rm {in}}$ ($θ$ rotates), and (iii) Swinging (SW) at middle $η_{\rm {in}}$ and high shear rate $\dotγ$ ($θ$ oscillates). All of three modes are accompanied by a membrane ($ϕ$) rotation. For microcapsules with low shear elasticity, the TB phase with no $ϕ$ rotation and the coexistence phase of SW and TB motions are induced by the energy barrier of $ϕ$ rotation. Synchronization of $ϕ$ rotation with TB rotation or SW oscillation occurs with integer ratios of rotational frequencies. At high shear elasticity, where a saddle point in the energy potential disappears, intermediate phases vanish, and either $ϕ$ or $θ$ rotation occurs. This phase behavior agrees with recent simulation results of microcapsules with low bending elasticity.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamic Modes of Red Blood Cells in Oscillatory Shear Flow

The dynamics of red blood cells (RBCs) in oscillatory shear flow was studied using differential equations of three variables: a shape parameter, the inclination angle $θ$, and phase angle $ϕ$ of the membrane rotation. In steady shear flow, three types of dynamics occur depending on the shear rate and viscosity ratio. i) tank-treading (TT): $ϕ$ rotates while the shape and $θ$ oscillate. ii) tumbling (TB): $θ$ rotates while the shape and $ϕ$ oscillate. iii) intermediate motion: both $ϕ$ and $θ$ rotate synchronously or intermittently. In oscillatory shear flow, RBCs show various dynamics based on these three motions. For a low shear frequency with zero mean shear rate, a limit-cycle oscillation occurs, based on the TT or TB rotation at a high or low shear amplitude, respectively. This TT-based oscillation well explains recent experiments. In the middle shear amplitude, RBCs show an intermittent or synchronized oscillation. As shear frequency increases, the vesicle oscillation becomes delayed with respect to the shear oscillation. At a high frequency, multiple limit-cycle oscillations coexist. For a high mean shear rate with small shear oscillation, the shape and $θ$ oscillate in the TT motion but only one attractor exists even at high shear frequencies. The measurement of these oscillatory modes is a promising tool for quantifying the viscoelasticity of RBCs and synthetic capsules.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamics of Fluid Vesicles in Flow through Structured Microchannels

The dynamics of fluid vesicles is studied under flow in microchannels, in which the width varies periodically along the channel. Three types of flow instabilities of prolate vesicles are found. For small quasi-spherical vesicles -- compared to the average channel width -- perturbation theory predicts a transition from a state with orientational oscillations of a fixed prolate shape to a state with shape oscillations of symmetrical ellipsoidal or bullet-like shapes with increasing flow velocity. Experimentally, such orientational oscillations are observed during the slow migration of a vesicle towards the centerline of the channel. For larger vesicles, mesoscale hydrodynamics simulations and experiments show similar symmetric shape oscillation at reduced volumes $V^* \gtrsim 0.9$. However, for non-spherical vesicles with $V^* \lesssim 0.9$, shapes are found with two symmetric or a single asymmetric tail.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Solvent-free coarse-grained lipid model for large-scale simulations

A coarse-grained molecular model, which consists of a spherical particle and an orientation vector, is proposed to simulate lipid membrane on a large length scale. The solvent is implicitly represented by an effective attractive interaction between particles. A bilayer structure is formed by orientation-dependent (tilt and bending) potentials. In this model, the membrane properties (bending rigidity, line tension of membrane edge, area compression modulus, lateral diffusion coefficient, and flip-flop rate) can be varied over broad ranges. The stability of the bilayer membrane is investigated via droplet-vesicle transition. The rupture of the bilayer and worm-like micelle formation can be induced by an increase in the spontaneous curvature of the monolayer membrane.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dynamical regimes and hydrodynamic lift of viscous vesicles under shear

The dynamics of two-dimensional viscous vesicles in shear flow, with different fluid viscosities $η_{\rm in}$ and $η_{\rm out}$ inside and outside, respectively, is studied using mesoscale simulation techniques. Besides the well-known tank-treading and tumbling motions, an oscillatory swinging motion is observed in the simulations for large shear rate. The existence of this swinging motion requires the excitation of higher-order undulation modes (beyond elliptical deformations) in two dimensions. Keller-Skalak theory is extended to deformable two-dimensional vesicles, such that a dynamical phase diagram can be predicted for the reduced shear rate and the viscosity contrast $η_{\rm in}/η_{\rm out}$. The simulation results are found to be in good agreement with the theoretical predictions, when thermal fluctuations are incorporated in the theory. Moreover, the hydrodynamic lift force, acting on vesicles under shear close to a wall, is determined from simulations for various viscosity contrasts. For comparison, the lift force is calculated numerically in the absence of thermal fluctuations using the boundary-integral method for equal inside and outside viscosities. Both methods show that the dependence of the lift force on the distance $y_{\rm {cm}}$ of the vesicle center of mass from the wall is well described by an effective power law $y_{\rm {cm}}^{-2}$ for intermediate distances $0.8 R_{\rm p} \lesssim y_{\rm {cm}} \lesssim 3 R_{\rm p}$ with vesicle radius $R_{\rm p}$. The boundary-integral calculation indicates that the lift force decays asymptotically as $1/[y_{\rm {cm}}\ln(y_{\rm {cm}})]$ far from the wall.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dynamics of Fluid Vesicles in Oscillatory Shear Flow

The dynamics of fluid vesicles in oscillatory shear flow was studied using differential equations of two variables: the Taylor deformation parameter and inclination angle $θ$. In a steady shear flow with a low viscosity $η_{\rm {in}}$ of internal fluid, the vesicles exhibit steady tank-treading motion with a constant inclination angle $θ_0$. In the oscillatory flow with a low shear frequency, $θ$ oscillates between $\pm θ_0$ or around $θ_0$ for zero or finite mean shear rate $\dotγ_{\rm m}$, respectively. As shear frequency $f_γ$ increases, the vesicle oscillation becomes delayed with respect to the shear oscillation, and the oscillation amplitude decreases. At high $f_γ$ with $\dotγ_{\rm m}=0$, another limit-cycle oscillation between $θ_0-π$ and $-θ_0$ is found to appear. In the steady flow, $θ$ periodically rotates (tumbling) at high $η_{\rm {in}}$, and $θ$ and the vesicle shape oscillate (swinging) at middle $η_{\rm {in}}$ and high shear rate. In the oscillatory flow, the coexistence of two or more limit-cycle oscillations can occur for low $f_γ$ in these phases. For the vesicle with a fixed shape, the angle $θ$ rotates back to the original position after an oscillation period. However, it is found that a preferred angle can be induced by small thermal fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Membrane simulation models from nm to $μ$m scale

Recent developments in lipid membrane models for simulations are reviewed. To reduce computational costs, various coarse-grained molecular models have been proposed. Among them, implicit solvent (solvent-free) molecular models are relatively more coarse-grained and efficient for simulating large bilayer membranes. On a $μ$m scale, the molecular details are typically negligible and the membrane can be described as a continuous curved surface. The theoretical models for fluid and elastic membranes with mesh or meshless discretizations are presented. As examples of applications, the dynamics of vesicles in flows, vesicle formation, and membrane fusion are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Swinging and Tumbling of Fluid Vesicles in Shear Flow

The dynamics of fluid vesicles in simple shear flow is studied using mesoscale simulations of dynamically-triangulated surfaces, as well as a theoretical approach based on two variables, a shape parameter and the inclination angle, which has no adjustable parameters. We show that between the well-known tank-treading and tumbling states, a new ``swinging'' state can appear. We predict the dynamic phase diagram as a function of the shear rate, the viscosities of the membrane and the internal fluid, and the reduced vesicle volume. Our results agree well with recent experiments.

0 citations0 reviewsNo code linkedOpen access