Shape transitions of fluid vesicles and red blood cells in capillary flows

A low-dimensional model for the red blood cell

The red blood cell (RBC) is an important determinant of the rheological properties of blood because of its predominant number density, special mechanical properties and dynamics. Here, we develop a new low-dimensional RBC model based on dissipative particle dynamics (DPD). The model is constructed as a closed-torus-like ring of 10 colloidal particles connected by wormlike chain springs combined with bending resistance. Each colloidal particle is represented by a single DPD particle with a repulsive core. The model is able to capture the essential mechanical properties of RBCs, and allows for economical exploration of the rheology of RBC suspensions. Specifically, we find that the linear and non-linear elastic deformations of healthy and malaria-infected cells match those obtained in optical tweezers experiments. Through simulations of some key features of blood flow in vessels, i.e., the cell-free layer (CFL), the Fahraeus effect and the Fahraeus-Lindqvist effect, we verify that the new model captures the essential shear flow properties of real blood, except for capillaries of sizes comparable to the cell diameter. Finally, we investigate the influence of a geometrical constriction in the flow on the enhancement of the downstream CFL. Our results are in agreement with recent experiments.

A multiscale red blood cell model with accurate mechanics, rheology, and dynamics

Red blood cells (RBCs) have highly deformable viscoelastic membranes exhibiting complex rheological response and rich hydrodynamic behavior governed by special elastic and bending properties and by the external/internal fluid and membrane viscosities. We present a multiscale RBC model that is able to predict RBC mechanics, rheology, and dynamics in agreement with experiments. Based on an analytic theory, the modeled membrane properties can be uniquely related to the experimentally established RBC macroscopic properties without any adjustment of parameters. The RBC linear and nonlinear elastic deformations match those obtained in optical-tweezers experiments. The rheological properties of the membrane are compared with those obtained in optical magnetic twisting cytometry, membrane thermal fluctuations, and creep followed by cell recovery. The dynamics of RBCs in shear and Poiseuille flows is tested against experiments and theoretical predictions, and the applicability of the latter is discussed. Our findings clearly indicate that a purely elastic model for the membrane cannot accurately represent the RBC's rheological properties and its dynamics, and therefore accurate modeling of a viscoelastic membrane is necessary.

One-particle-thick, solvent-free, coarse-grained model for biological and biomimetic fluid membranes

Biological membranes are involved in numerous intriguing biophysical and biological cellular phenomena of different length scales, ranging from nanoscale raft formation, vesiculation, to microscale shape transformations. With extended length and time scales as compared to atomistic simulations, solvent-free coarse-grained membrane models have been exploited in mesoscopic membrane simulations. In this study, we present a one-particle-thick fluid membrane model, where each particle represents a cluster of lipid molecules. The model features an anisotropic interparticle pair potential with the interaction strength weighed by the relative particle orientations. With the anisotropic pair potential, particles can robustly self-assemble into fluid membranes with experimentally relevant bending rigidity. Despite its simple mathematical form, the model is highly tunable. Three potential parameters separately and effectively control diffusivity, bending rigidity, and spontaneous curvature of the model membrane. As demonstrated by selected examples, our model can naturally simulate dynamics of phase separation in multicomponent membranes and the topological change of fluid vesicles.

Systematic coarse-graining of spectrin-level red blood cell models

We present a rigorous procedure to derive coarse-grained red blood cell (RBC) models, which yield accurate mechanical response. Based on a semi-analytic theory the linear and nonlinear elastic properties of healthy and infected RBCs in malaria can be matched with those obtained in optical tweezers stretching experiments. The present analysis predicts correctly the membrane Young's modulus in contrast to about 50% error in predictions by previous models. In addition, we develop a stress-free model which avoids a number of pitfalls of existing RBC models, such as non-smooth or poorly controlled equilibrium shape and dependence of the mechanical properties on the initial triangulation quality. Here we employ dissipative particle dynamics for the implementation but the proposed model is general and suitable for use in many existing continuum and particle-based numerical methods.

Dynamic modes of microcapsules in steady shear flow: effects of bending and shear elasticities

The dynamics of microcapsules in steady shear flow were studied using a theoretical approach based on three variables: the Taylor deformation parameter α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 at low viscosity ηin of internal fluid (αD and θ relaxes to constant values), (ii) tumbling (TB) at high ηin (θ rotates), and (iii) swinging (SW) at middle ηin and high shear rates γ (θ 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.

Polymer-functionalized nanoparticles: from stealth viruses to biocompatible quantum dots

In this article, we focus on nanoparticles that have been functionalized by polymers. We draw our examples from nanoparticle systems that have found biomedical and therapeutic applications. Our aim is to highlight the physical principles that might explain why these systems have been found to be successful in biomedical applications and to highlight other physical properties that might lead to new applications. We consider viruses, gold nanoparticles, magnetic nanoparticles and quantum dots, focussing attention on the ways in which functionalization by polymers has been used to alter the physical characteristics of the particular nanoparticle to improve its function as a possible therapy. In the case of viral vectors, polymer functionalization tunes the biocompatibility, suppressing the binding of antibodies and conferring the nanoparticle with stealth properties. By contrast, the inorganic nanoparticles comprise materials in a form that is not normally encountered in the human body, and polymer functionalization is necessary to ensure biocompatibility.

Why do red blood cells have asymmetric shapes even in a symmetric flow?

Understanding why red blood cells (RBCs) move with an asymmetric shape (slipperlike shape) in small blood vessels is a long-standing puzzle in blood circulatory research. By considering a vesicle (a model system for RBCs), we discovered that the slipper shape results from a loss in stability of the symmetric shape. It is shown that the adoption of a slipper shape causes a significant decrease in the velocity difference between the cell and the imposed flow, thus providing higher flow efficiency for RBCs. Higher membrane rigidity leads to a dramatic change in the slipper morphology, thus offering a potential diagnostic tool for cell pathologies.

A particle dynamic model of red blood cell aggregation kinetics

To elucidate the relationship between microscopic red blood cell (RBC) interactions and macroscopic rheological behavior, we propose a two-dimensional particle model capable of mimicking the main characteristics of RBC aggregation kinetics. The mechanical model of RBCs sheared in Couette flow is based on Newton law. We assumed a hydrodynamic force to move particles, a force to describe aggregation and an elasticity force. The role of molecular mass and concentration of neutral polymers on aggregation [Neu, B., and H. J. Meiselman. Biophys. J. 83:2482-2490, 2002] could be mimicked. Specifically, it was shown that for any shear rate (SR), the mean aggregate size (MAS) grew with time until it reached a constant value, which is consistent with in vitro experiments. It was also demonstrated that we could mimic the modal relationship between MAS and SR and the occurrence of maximum aggregation at about 0.1 s(-1). As anticipated, simulations indicated that an increase in aggregation force augmented MAS. Further, augmentation of the depletion layer thickness influenced MAS only for SR close to zero, which is a new finding. To conclude, our contribution reveals that the aggregation force intensity and SR influence the steady state MAS, and that the depletion and layer thickness affect the aggregation speed.

Swinging and synchronized rotations of red blood cells in simple shear flow

The dynamics of red blood cells (RBCs) in simple shear flow was studied using a theoretical approach based on three variables: a shape parameter, the inclination angle theta, and phase angle phi of the membrane rotation. At high shear rate and low viscosity contrast of internal fluid, RBCs exhibit tank-treading motion, where phi rotates with swinging oscillation of shape and theta . At low shear rate, tumbling motion occurs and theta rotates. In the middle region between these two phases, it is found that synchronized rotation of phi and theta with integer ratios of the frequencies occurs in addition to intermittent rotation. These dynamics are robust to the modification of the potential of the RBC shape and membrane rotation. Our results agree well with recent experiments.

Elastic capsules in shear flow: analytical solutions for constant and time-dependent shear rates

We investigate the dynamics of microcapsules in linear shear flow within a reduced model with two degrees of freedom. In previous work for steady shear flow, the dynamic phases of this model, i.e. swinging, tumbling and intermittent behaviour, have been identified using numerical methods. In this paper, we integrate the equations of motion in the quasi-spherical limit analytically for time-constant and time-dependent shear flow using matched asymptotic expansions. Using this method, we find analytical expressions for the mean tumbling rate in general time-dependent shear flow. The capsule dynamics is studied in more detail when the inverse shear rate is harmonically modulated around a constant mean value for which a dynamic phase diagram is constructed. By a judicious choice of both modulation frequency and phase, tumbling motion can be induced even if the mean shear rate corresponds to the swinging regime. We derive expressions for the amplitude and width of the resonance peaks as a function of the modulation frequency.

Spatially resolved simulations of membrane reactions and dynamics: multipolar reaction DPD

Biophysical chemistry of mesoscale systems and quantitative modeling in systems biology now require a simulation methodology unifying chemical reaction kinetics with essential collective physics. This will enable the study of the collective dynamics of complex chemical and structural systems in a spatially resolved manner with a combinatorially complex variety of different system constituents. In order to allow a direct link-up with experimental data (e.g. high-throughput fluorescence images) the simulations must be constructed locally, i.e. mesoscale phenomena have to emerge from local composition and interactions that can be extracted from experimental data. Under suitable conditions, the simulation of such local interactions must lead to processes such as vesicle budding, transport of membrane-bounded compartments and protein sorting, all of which result from a sophisticated interplay between chemical and mechanical processes and require the link-up of different length scales. In this work, we show that introducing multipolar interactions between particles in dissipative particle dynamics (DPD) leads to extended membrane structures emerging in a self-organized manner and exhibiting the necessary mechanical stability for transport, correct scaling behavior, and membrane fluidity so as to provide a two-dimensional self-organizing dynamic reaction environment for kinetic studies in the context of cell biology.

The Deformation Behavior of Multiple Red Blood Cells in a Capillary Vessel

The deformation of multiple red blood cells in a capillary flow was studied numerically. The immersed boundary method was used for the fluid red blood cells interaction. The membrane of the red blood cell was modeled as a hyperelastic thin shell. The numerical results show that the apparent viscosity in the capillary flow is more sensitive to the change of shear coefficient of the membrane than the bending coefficient and surface dilation coefficient, and the increase in the shear coefficient results in an increase in the pressure drop in the blood flow in capillary vessels in order to sustain the same flux rate of red blood cells.

Modeling morphological instabilities in lipid membranes with anchored amphiphilic polymers

Anchoring molecules, like amphiphilic polymers, are able to dynamically regulate membrane morphology. Such molecules insert their hydrophobic groups into the bilayer, generating a local membrane curvature. In order to minimize the elastic energy penalty, a dynamic shape instability may occur, as in the case of the curvature-driven pearling instability or the polymer-induced tubulation of lipid vesicles. We review recent works on modeling of such instabilities by means of a mesoscopic dynamic model of the phase-field kind, which take into account the bending energy of lipid bilayers.

Vesicles in Poiseuille flow

Blood microcirculation critically depends on the migration of red cells towards the flow centerline. We identify theoretically the ratio of the inner over the outer fluid viscosities lambda as a key parameter. At low lambda, the vesicle deforms into a tank-treading ellipsoid shape far away from the flow centerline. The migration is always towards the flow centerline, unlike drops. Above a critical lambda, the vesicle tumbles or breaths and migration is suppressed. A surprising coexistence of two types of shapes at the centerline, a bulletlike and a parachutelike shape, is predicted.

Flow-induced clustering and alignment of vesicles and red blood cells in microcapillaries

The recent development of microfluidic devices allows the investigation and manipulation of individual liquid microdroplets, capsules, and cells. The collective behavior of several red blood cells (RBCs) or microcapsules in narrow capillaries determines their flow-induced morphology, arrangement, and effective viscosity. Of fundamental interest here is the relation between the flow behavior and the elasticity and deformability of these objects, their long-range hydrodynamic interactions in microchannels, and thermal membrane undulations. We study these mechanisms in an in silico model, which combines a particle-based mesoscale simulation technique for the fluid hydrodynamics with a triangulated-membrane model. The 2 essential control parameters are the volume fraction of RBCs (the tube hematocrit, H(T)), and the flow velocity. Our simulations show that already at very low H(T), the deformability of RBCs implies a flow-induced cluster formation above a threshold flow velocity. At higher H(T) values, we predict 3 distinct phases: one consisting of disordered biconcave-disk-shaped RBCs, another with parachute-shaped RBCs aligned in a single file, and a third with slipper-shaped RBCs arranged as 2 parallel interdigitated rows. The deformation-mediated clustering and the arrangements of RBCs and microcapsules are relevant for many potential applications in physics, biology, and medicine, such as blood diagnosis and cell sorting in microfluidic devices.

Diffusion in a fluid membrane with a flexible cortical cytoskeleton

We calculate the influence of a flexible network of long-chain proteins, which is anchored to a fluid membrane, on protein diffusion in this membrane. This is a model for the cortical cytoskeleton and the lipid bilayer of the red blood cell, which we apply to predict the influence of the cytoskeleton on the diffusion coefficient of a mobile band 3 protein. Using the pressure field that the cytoskeleton exerts on the membrane, from the steric repulsion between the diffusing protein and the cytoskeletal filaments, we define a potential landscape for the diffusion within the bilayer. We study the changes to the diffusion coefficient on removal of one type of anchor proteins, e.g., in several hemolytic anemias, as well as for isotropic and anisotropic stretching of the cytoskeleton. We predict an overall increase of the diffusion for a smaller number of anchor proteins and increased diffusion for anisotropic stretching in the direction of the stretch, because of the decrease in the spatial frequency as well as in the height of the potential barriers.

Relaxation dynamics of fluid membranes

We study the effect of membrane viscosity in the dynamics of liquid membranes-possibly with free or internal boundaries-driven by conservative forces (curvature elasticity and line tension) and dragged by the bulk dissipation of the ambient fluid and the friction occurring when the amphiphilic molecules move relative to each other. To this end, we formulate a continuum model which includes a form of the governing equations for a two-dimensional viscous fluid moving on a curved, time-evolving surface. The effect of membrane viscosity has received very limited attention in previous continuum studies of the dynamics of fluid membranes, although recent coarse-grained discrete simulations suggest its importance. By applying our model to the study of vesiculation and membrane fusion in a simplified geometry, we conclude that membrane viscosity plays a dominant role in the relaxation dynamics of fluid membranes of sizes comparable to those found in eukaryotic cells, and is not negligible in many large synthetic systems of current interest.

Coarse-grained red blood cell model with accurate mechanical properties, rheology and dynamics

We present a coarse-grained red blood cell (RBC) model with accurate and realistic mechanical properties, rheology and dynamics. The modeled membrane is represented by a triangular mesh which incorporates shear inplane energy, bending energy, and area and volume conservation constraints. The macroscopic membrane elastic properties are imposed through semi-analytic theory, and are matched with those obtained in optical tweezers stretching experiments. Rheological measurements characterized by time-dependent complex modulus are extracted from the membrane thermal fluctuations, and compared with those obtained from the optical magnetic twisting cytometry results. The results allow us to define a meaningful characteristic time of the membrane. The dynamics of RBCs observed in shear flow suggests that a purely elastic model for the RBC membrane is not appropriate, and therefore a viscoelastic model is required. The set of proposed analyses and numerical tests can be used as a complete model testbed in order to calibrate the modeled viscoelastic membranes to accurately represent RBCs in health and disease.

Attractive colloidal rods in shear flow

The effect of shear flow on the isotropic-nematic phase transition of attractive colloidal rods is investigated by a combination of simulations and experiments. The isotropic phase aligns with the flow, while the nematic phase undergoes a collective rotational motion which frustrates the merging of the coexisting regions. The location of binodals, spinodals, and the tumbling-to-aligning transition line in the shear-rate versus concentration plane are investigated. The phase diagrams in the shear-concentration plane for the various strengths of attractions can be mapped onto a master curve by appropriate scaling.

Fast measurements of concentration profiles inside deformable objects in microflows with reduced spatial coherence digital holography

We investigate the use of a digital holographic microscope working with partially coherent spatial illumination to study concentration profiles inside confined deformable bodies flowing in microchannels. The studied phenomenon is rapidly changing in time and requires the recording of the complete holographic information for every frame. For this purpose, we implemented one of the classical methods of off-axis digital holography: the Fourier method. Digital holography allows one to numerically investigate a volume by refocusing the different planes of depth, allowing one to locate the objects under investigation in three dimensions. Furthermore, the phase is directly related to the refractive index, thus to the concentration inside the body. Based on simple symmetry assumptions, we present an original method for determining the concentration profiles inside deformable objects in microconfined flows. Details of the optical and numerical implementation, as well as exemplative experimental results are presented.

Accurate coarse-grained modeling of red blood cells

We develop a systematic coarse-graining procedure for modeling red blood cells (RBCs) using arguments based on mean-field theory. The three-dimensional RBC membrane model takes into account the bending energy, in-plane shear energy, and constraints of fixed surface area and fixed enclosed volume. The coarse-graining procedure is general, it can be used for arbitrary level of coarse-graining and does not employ any fitting parameters. The sensitivity of the coarse-grained model is investigated and its behavior is validated against available experimental data and in dissipative particle dynamics (DPD) simulations of RBCs in capillary and shear flows.

Cellular-scale hydrodynamics

Microfluidic tools are providing many new insights into the chemical, physical and physicochemical responses of cells. Both suspension-level and single-cell measurements have been studied. We review our studies of these kinds of problems for red blood cells with particular focus on the shapes of individual cells in confined geometries, the development and use of a 'differential manometer' for evaluating the mechanical response of individual cells or other objects flowing in confined geometries, and the cross-streamline drift of cells that pass through a constriction. In particular, we show how fluid mechanical effects on suspended cells can be studied systematically in small devices, and how these features can be exploited to develop methods for characterizing physicochemical responses and possibly for the diagnosis of cellular-scale changes to environmental factors.

Toward a thermodynamically consistent picture of the phase-field model of vesicles: curvature energy

We extend our recent work on phase-field model for vesicles [D. Jamet and C. Misbah, Phys. Rev. E 76, 051907 (2007)]-where only the membrane local incompressibility was treated-to the situation where the bending forces and spontaneous curvature are included. We show how the general phase-field equations can be derived within a thermodynamic consistent picture. We analyze a general form of the bending energy, where the Helfrich bending force is treated as a special case. The dynamical evolution equation derived here for the velocity field allows one to write down a constitutive law of the composite fluid: The ambient fluid plus the membrane. This constitutive law has a viscoelastic form, the viscous part arises from the fluid, while the elastic one represents the action of the membrane. It is shown that the elastic stress tensor is not symmetric, owing to bending torque, inherent to a diffuse membrane model.

Transport coefficients of off-lattice mesoscale-hydrodynamics simulation techniques

The viscosity and self-diffusion constant of particle-based mesoscale hydrodynamic methods, multiparticle collision dynamics (MPC), and dissipative particle dynamics, are investigated, both with and without angular-momentum conservation. Analytical results are derived for fluids with an ideal-gas equation of state and a finite-time-step dynamics, and compared with simulation data. In particular, the viscosity is derived in a general form for all variants of the MPC method. In general, very good agreement between theory and simulations is obtained.

Intravital two-photon microscopy: focus on speed and time resolved imaging modalities

Initially used mainly in the neurosciences, two-photon microscopy has become a powerful tool for the analysis of immunological processes. Here, we describe currently available two-photon microscopy techniques with a focus on novel approaches that allow very high image acquisition rates compared with state-of-the-art systems. This improvement is achieved through a parallelization of the excitation process: multiple beams scan the sample simultaneously, and the fluorescence is collected with sensitive charge-coupled device (CCD)-based line or field detectors. The new technique's performance is compared with conventional single beam laser-scanning systems that detect signals by means of photomultipliers. We also discuss the use of time- and polarization-resolved fluorescence detection, especially fluorescence lifetime imaging (FLIM), which goes beyond simple detection of cells and tissue structures and allows insight into cellular physiology. We focus on the analysis of endogenous fluorophores such as NAD(P)H as a way to analyze the redox status in cells with subcellular resolution. Here, high-speed imaging setups in combination with novel ways of data analysis allow the generation of FLIM data sets almost in real time. The implications of this technology for the analysis of immune reactions and other cellular processes are discussed.

Dynamic states of cells adhering in shear flow: from slipping to rolling

Motivated by rolling adhesion of white blood cells in the vasculature, we study how cells move in linear shear flow above a wall to which they can adhere via specific receptor-ligand bonds. Our computer simulations are based on a Langevin equation accounting for hydrodynamic interactions, thermal fluctuations, and adhesive interactions. In contrast to earlier approaches, our model not only includes stochastic rules for the formation and rupture of bonds, but also fully resolves both receptor and ligand positions. We identify five different dynamic states of motion in regard to the translational and angular velocities of the cell. The transitions between the different states are mapped out in a dynamic state diagram as a function of the rates for bond formation and rupture. For example, as the cell starts to adhere under the action of bonds, its translational and angular velocities become synchronized and the dynamic state changes from slipping to rolling. We also investigate the effect of nonmolecular parameters. In particular, we find that an increase in viscosity of the medium leads to a characteristic expansion of the region of stable rolling to the expense of the region of firm adhesion, but not to the expense of the regions of free or transient motion. Our results can be used in an inverse approach to determine single bond parameters from flow chamber data on rolling adhesion.

Vesicles and red blood cells in shear flow

We describe the similarities and the specificities of the behaviour of individual soft particles, namely, drops, lipid vesicles and red blood cells subjected to a shear flow. We highlight that their motion depends in a non-trivial way on the particle mechanical properties. We detail the effect of the presence of a wall with or without wall-particle attractive interaction from a biological perspective.

Two-dimensional fluctuating vesicles in linear shear flow

The stochastic motion of a two-dimensional vesicle in linear shear flow is studied at finite temperature. In the limit of small deformations from a circle, Langevin-type equations of motion are derived, which are highly nonlinear due to the constraint of constant perimeter length. These equations are solved in the low-temperature limit and using a mean-field approach, in which the length constraint is satisfied only on average. The constraint imposes non-trivial correlations between the lowest deformation modes at low temperature. We also simulate a vesicle in a hydrodynamic solvent by using the multi-particle collision dynamics technique, both in the quasi-circular regime and for larger deformations, and compare the stationary deformation correlation functions and the time autocorrelation functions with theoretical predictions. Good agreement between theory and simulations is obtained.

Dynamics of nearly spherical vesicles in an external flow

Tank-treading, tumbling, and trembling are different types of the vesicle behavior in an external flow. We derive a dynamical equation enabling us to establish a state of nearly spherical vesicles. For a 2D external flow, the character of the vesicle dynamics is determined by two dimensionless parameters, depending on the vesicle excess area, fluid viscosities, membrane viscosity and bending modulus, strength of the flow, and ratio of the elongational and rotational components of the flow. The tank-treading to tumbling transition occurs via a saddle-node bifurcation, whereas the tank-treading to trembling transition occurs via a Hopf bifurcation. A slowdown of vesicle dynamics should be observed in a vicinity of a point separating the transition lines. We show that the slowdown can be described by a power law with two different critical exponents 1/4 and 1/2 corresponding to the slowdown of tumbling and trembling cycles.

Continuous blood cell separation by hydrophoretic filtration

We propose a new hydrophoretic method for continuous blood cell separation using a microfluidic device composed of slanted obstacles and filtration obstacles. The slanted obstacles have a larger height and gap than the particles in order to focus them to a sidewall by hydrophoresis. In the successive structure, the height and gap of the filtration obstacles with a filtration pore are set between the diameters of small and large particles, which defines the critical separation diameter. Accordingly, the particles smaller than the criterion freely pass through the gap and keep their focused position. In contrast, the particles larger than the criterion collide against the filtration obstacle and move into the filtration pore. The microfluidic device was characterized with polystyrene beads with a minimum diameter difference of 7.3%. We completely separated polystyrene microbeads of 9 and 12 microm diameter with a separation resolution of approximately 6.2. This resolution is increased by 6.4-fold compared with our previous separation method based on hydrophoresis (S. Choi and J.-K. Park, Lab Chip, 2007, 7, 890, ref. 1). In the isolation of white blood cells (WBCs) from red blood cells (RBCs), the microfluidic device isolated WBCs with 210-fold enrichment within a short filtration time of approximately 0.3 s. These results show that the device can be useful for the binary separation of a wide range of biological particles by size. The hydrophoretic filtration as a sample preparation unit offers potential for a power-free cell sorter to be integrated into disposable lab-on-a-chip devices.

Towards a thermodynamically consistent picture of the phase-field model of vesicles: local membrane incompressibility

A phase-field model for vesicles including hydrodynamics was presented in two and three dimensions [T. Biben and C. Misbah, Phys. Rev. E 67, 031908 (2003); T. Biben, K. Kassner, and C. Misbah, Phys. Rev. E 72, 041921 (2005)]. A particularly important feature for vesicles is that their membrane is locally incompressible. In these works a tension field defined everywhere in the bulk was introduced in order to fulfill local membrane inextensibility. Here we reconsider the original model by treating the phase field as a thermodynamic variable and develop a picture which is consistent with the second law of thermodynamics. This enables us to write the phase-field evolution equations in terms of a thermodynamical potential. This potential acquires, at global equilibrium, a Lyapunov functional character. The goal of this paper is twofold: (i) The first and primary goal is purely conceptual, in that we can write down a first and second principle for membranes, from which the evolution equations follow, thanks to the evaluation of the entropy production and the use of concepts of irreversible thermodynamics. (ii) Due to the monotonous character of the evolution of the functional (at global equilibrium), we expect this formulation to be more appropriate for numerical studies. The formalism developed to account for the local incompressibility of the membrane is believed to offer a systematic framework in order to include naturally other physical ingredients, as briefly discussed here and demonstrated in future works.

Fluctuations of coupled fluid and solid membranes with application to red blood cells

The fluctuation spectra and the intermembrane interaction of two membranes at a fixed average distance are investigated. Each membrane can either be a fluid or a solid membrane, and in isolation, its fluctuations are described by a bare or a wave-vector-dependent bending modulus, respectively. The membranes interact via their excluded-volume interaction; the average distance is maintained by an external, homogeneous pressure. For strong coupling, the fluctuations can be described by a single, effective membrane that combines the elastic properties. For weak coupling, the fluctuations of the individual, noninteracting membranes are recovered. The case of a composite membrane consisting of one fluid and one solid membrane can serve as a microscopic model for the plasma membrane and cytoskeleton of the red blood cell. We find that, despite the complex microstructure of bilayers and cytoskeletons in a real cell, the fluctuations with wavelengths lambda greater, similar 400 nm are well described by the fluctuations of a single, polymerized membrane (provided that there are no inhomogeneities of the microstructure). The model is applied to the fluctuation data of discocytes ("normal" red blood cells), a stomatocyte, and an echinocyte. The elastic parameters of the membrane and an effective temperature that quantifies active, metabolically driven fluctuations are extracted from the experiments.

Relevance of angular momentum conservation in mesoscale hydrodynamics simulations

The angular momentum is conserved in fluids with a few exceptions such as ferrofluids. However, it can be violated locally in fluid simulations to reduce computational costs. The effects of this violation are investigated using a particle-based simulation method, multiparticle collision dynamics, which can switch on or off angular-momentum conservation. To this end, we study circular Couette flows between concentric and eccentric cylinders, where nonphysical torques due to the lack of the angular-momentum conservation are found whereas the velocity field is not affected. In addition, in simulations of fluids with different viscosities in contact and star polymers in solvent, incorrect angular velocities occur. These results quantitatively agree with the theoretical predictions based on the macroscopic stress tensor.

Biological Soft Materials

With one or two exceptions, biological materials are "soft", meaning that they combine viscous and elastic elements. This mechanical behavior results from self-assembled supramolecular structures that are stabilized by noncovalent interactions. It is an ongoing and profound challenge to understand the self-organization of biological materials. In many cases, concepts can be imported from soft-matter physics and chemistry, which have traditionally focused on materials such as colloids, polymers, surfactants, and liquid crystals. Using these ideas, it is possible to gain a new perspective on phenomena as diverse as DNA condensation, protein and peptide fibrillization, lipid partitioning in rafts, vesicle fusion and budding, and others, as discussed in this selective review of recent highlights from the literature.

Tank-tread frequency of the red cell membrane: dependence on the viscosity of the suspending medium

Single human red cells were suspended in media with viscosities ranging from 12.9 to 109 mPa s and subjected to shear flow ranging from 1/s to 290/s in a rheoscope. This is a transparent cone-plate chamber adapted to a microscope. The motion of the membrane around red cells oriented in a steady-state fashion in the shear field (tank-tread motion) was videotaped. The projected length and width of the cells as well as the frequency of tank-tread motion were measured. One-thousand eight-hundred seventy-three cells of three blood donors were evaluated. The frequency increased with the mean shear rate in an almost linear fashion. The slope of this dependence increased weakly with the viscosity of the suspending medium. No correlation was found between the frequency and four morphological red cell parameters: the projected length and width of the cells as well as the ratio and the square root of the product of these quantities. The energy dissipation within the red cell membrane was estimated based on the measured parameters and compared to the energy dissipation in the undisturbed shear flow. At constant mean shear rate the rise of the energy dissipation with viscosity is slower whereas at constant viscosity the rise with the shear rate is steeper than in the undisturbed shear flow. A fit of the data collected in this work to a theoretical red cell model might allow one to determine intrinsic mechanical constants in the low deformation regime.

Shape transitions in lipid membranes and protein mediated vesicle fusion and fission

In the cell, the plasma membrane is often densely decorated by transmembrane proteins. The morphology and dynamics of the membrane are strongly influenced by the presence of proteins. In this paper, we use a coarse-grained model to explore the composite membrane-protein system and develop a simulation methodology based on thermodynamic integration to examine free energy changes during membrane shape transitions. The authors show that a critical concentration of conical membrane proteins or proteins with nonzero spontaneous curvature can drive the formation of small vesicles. The driving force of vesicle budding stems from the preference of proteins to gather in regions of high curvature. A sufficiently high concentration of proteins therefore can influence the topology of the membrane. The biological significance of our results is discussed.

Shear-driven release of a bud from a multicomponent vesicle

The authors study the response of a multicomponent budded vesicle to an imposed shear flow using dissipative particle dynamics. Under certain circumstances, phase separation in the vesicle membrane leads to the formation of a minority domain which deforms into a nearly spherical bud in order to reduce its interfacial energy. The authors show that an imposed shear force has a varying effect on the vesicle, tending either to separate the bud from the vesicle or to stretch the bud open, depending on the vesicle orientation. The authors examine the interplay of membrane bending rigidity, line tension, and shear in determining the behavior of the vesicle. With the appropriate design, vesicles can be made to release buds in a controlled manner in response to shear. The authors outline a regime in which bud release is favorable.

Coarse-grained simulation of amphiphilic self-assembly

The authors present a computer simulation study of amphiphilic self-assembly performed using a computationally efficient single-site model based on Gay-Berne [J. Chem. Phys. 74, 3316 (1981)] and Lennard-Jones particles. Molecular dynamics simulations of these systems show that free self-assembly of micellar, bilayer, and inverse micelle arrangements can be readily achieved for a single model parametrization. This self-assembly is predominantly driven by the anisotropy of the amphiphile-solvent interaction, amphiphile-amphiphile dispersive interactions being found to be of secondary importance. While amphiphile concentration is the main determinant of phase stability, molecular parameters such as head group size and interaction strength also have measurable affects on system properties.

Dynamics of vesicle self-assembly and dissolution

The dynamics of membranes is studied on the basis of a particle-based meshless surface model, which was introduced earlier [Phys. Rev. E 73, 021903 (2006)]. The model describes fluid membranes with bending energy and-in the case of membranes with boundaries-line tension. The effects of hydrodynamic interactions are investigated by comparing Brownian dynamics with a particle-based mesoscale solvent simulation (multiparticle collision dynamics). Particles self-assemble into vesicles via disk-shaped membrane patches. The time evolution of assembly is found to consist of three steps: particle assembly into discoidal clusters, aggregation of clusters into larger membrane patches, and finally vesicle formation. The time dependence of the cluster distribution and the mean cluster size is evaluated and compared with the predictions of Smoluchowski rate equations. On the other hand, when the line tension is suddenly decreased (or the temperature is increased), vesicles dissolve via pore formation in the membrane. Hydrodynamic interactions are found to speed up the dynamics in both cases. Furthermore, hydrodynamics makes vesicle more spherical in the membrane-closure process.

Particles held by springs in a linear shear flow exhibit oscillatory motion

The dynamics of small spheres that are held by linear springs in a low Reynolds number shear flow at neighboring locations is investigated. The flow elongates the beads, and the interplay of the shear gradient with the nonlinear behavior of the hydrodynamic interaction among the spheres causes in a large range of parameters a bifurcation to a surprising oscillatory bead motion. The parameter ranges wherein this bifurcation is either super- or subcritical are determined.

Meshless membrane model based on the moving least-squares method

A meshless particle-based membrane model is proposed. The particles possess no internal degree of freedom and interact via a potential, which has three different contributions: a short-range repulsive pair potential, an attractive multibody potential, and a curvature potential based on the moving least-squares method. Brownian dynamics simulations are employed to demonstrate that the particles self-assemble into a membrane and to study equilibrium properties, such as bending rigidity, surface tension, line tension, and diffusion constant. The bending rigidity and line tension are shown to depend on different potential parameters and can therefore be varied independently. The finite-size effects of nearly planar membranes are investigated. This model is well suited to study the membrane dynamics with topological changes.