New analytical progress in the theory of vesicles under linear flow

Vesicles are becoming a quite popular model for the study of red blood cells (RBCs). This is a free boundary problem which is rather difficult to handle theoretically. Quantitative computational approaches constitute also a challenge. In addition, with numerical studies, it is not easy to scan within a reasonable time the whole parameter space. Therefore, having quantitative analytical results is an essential advance that provides deeper understanding of observed features and can be used to accompany and possibly guide further numerical development. In this paper shape evolution equations for a vesicle in a shear flow are derived analytically with precision being cubic (which is quadratic in previous theories) with regard to the deformation of the vesicle relative to a spherical shape. The phase diagram distinguishing regions of parameters where different types of motion (tank-treading, tumbling and vacillating-breathing) are manifested is presented. This theory reveals unsuspected features: including higher order terms and harmonics (even if they are not directly excited by the shear flow) is necessary, whatever the shape is close to a sphere. Not only does this theory cure a quite large quantitative discrepancy between previous theories and recent experiments and numerical studies, but also it reveals a new phenomenon: the VB mode band in parameter space, which is believed to saturate after a moderate shear rate, exhibits a striking widening beyond a critical shear rate. The widening results from excitation of fourth order harmonic. The obtained phase diagram is in a remarkably good agreement with recent three dimensional numerical simulations based on the boundary integral formulation. Comparison of our results with experiments is systematically made.Comment: a tex file and 6 figure

Bulk and wetting phenomena in a colloidal mixture of hard spheres and platelets

Density functional theory is used to study binary colloidal fluids consisting of hard spheres and thin platelets in their bulk and near a planar hard wall. This system exhibits liquid-liquid coexistence of a phase that is rich in spheres (poor in platelets) and a phase that is poor in spheres (rich in platelets). For the mixture near a planar hard wall, we find that the phase rich in spheres wets the wall completely upon approaching the liquid demixing binodal from the sphere-poor phase, provided the concentration of the platelets is smaller than a threshold value which marks a first-order wetting transition at coexistence. No layering transitions are found in contrast to recent studies on binary mixtures of spheres and non-adsorbing polymers or thin hard rods.Comment: 6 pages, 4 figure

New Approach on the General Shape Equation of Axisymmetric Vesicles

The general Helfrich shape equation determined by minimizing the curvature free energy describes the equilibrium shapes of the axisymmetric lipid bilayer vesicles in different conditions. It is a non-linear differential equation with variable coefficients. In this letter, by analyzing the unique property of the solution, we change this shape equation into a system of the two differential equations. One of them is a linear differential equation. This equation system contains all of the known rigorous solutions of the general shape equation. And the more general constraint conditions are found for the solution of the general shape equation.Comment: 8 pages, LaTex, submit to Mod. Phys. Lett.

Hydrodynamic lift of vesicles under shear flow in microgravity

The dynamics of a vesicle suspension in a shear flow between parallel plates has been investigated under microgravity conditions, where vesicles are only submitted to hydrodynamic effects such as lift forces due to the presence of walls and drag forces. The temporal evolution of the spatial distribution of the vesicles has been recorded thanks to digital holographic microscopy, during parabolic flights and under normal gravity conditions. The collected data demonstrates that vesicles are pushed away from the walls with a lift velocity proportional to $\dot{\gamma} R^3/z^2$ where $\dot{\gamma}$ is the shear rate, $R$ the vesicle radius and $z$ its distance from the wall. This scaling as well as the dependence of the lift velocity upon vesicle aspect ratio are consistent with theoretical predictions by Olla [J. Phys. II France {\bf 7}, 1533--1540 (1997)].Comment: 6 pages, 8 figure

Two-dimensional Vesicle dynamics under shear flow: effect of confinement

Dynamics of a single vesicle under shear flow between two parallel plates is studied using two-dimensional lattice-Boltzmann simulations. We first present how we adapted the lattice-Boltzmann method to simulate vesicle dynamics, using an approach known from the immersed boundary method. The fluid flow is computed on an Eulerian regular fixed mesh while the location of the vesicle membrane is tracked by a Lagrangian moving mesh. As benchmarking tests, the known vesicle equilibrium shapes in a fluid at rest are found and the dynamical behavior of a vesicle under simple shear flow is being reproduced. Further, we focus on investigating the effect of the confinement on the dynamics, a question that has received little attention so far. In particular, we study how the vesicle steady inclination angle in the tank-treading regime depends on the degree of confinement. The influence of the confinement on the effective viscosity of the composite fluid is also analyzed. At a given reduced volume (the swelling degree) of a vesicle we find that both the inclination angle, and the membrane tank-treading velocity decrease with increasing confinement. At sufficiently large degree of confinement the tank-treading velocity exhibits a non-monotonous dependence on the reduced volume and the effective viscosity shows a nonlinear behavior.Comment: 12 pages, 8 figure

Gravity currents from a dam-break in a rotating channel

Author Posting. © Cambridge University Press, 2005. This article is posted here by permission of Cambridge University Press for personal use, not for redistribution. The definitive version was published in Journal of Fluid Mechanics 536 (2005): 253-283, doi:10.1017/S0022112005004544.The generation of a gravity current by the release of a semi-infinite region of buoyant fluid of depth $H$ overlying a deeper, denser and quiescent lower layer in a rotating channel of width $w$ is considered. Previous studies have focused on the characteristics of the gravity current head region and produced relations for the gravity current speed $c_{b}$ and width $w_b$ as a functions of the local current depth along the wall $h_b$, reduced gravity $g^\prime$, and Coriolis frequency $f$. Here, the dam-break problem is solved analytically by the method of characteristics assuming reduced-gravity flow, uniform potential vorticity and a semigeostrophic balance. The solution makes use of a local gravity current speed relation $c_{b} \,{=}\, c_b(h_b,\ldots)$ and a continuity constraint at the head to close the problem. The initial value solution links the local gravity current properties to the initiating dam-break conditions. The flow downstream of the dam consists of a rarefaction joined to a uniform gravity current with width $w_b$ (${\le}\, w$) and depth on the right-hand wall of $h_b$, terminated at the head moving at speed $c_b$. The solution gives $h_b$, $c_b$, $w_b$ and the transport of the boundary current as functions of $w/L_R$, where $L_R \,{=}\, \sqrt{g^\prime H}/f$ is the deformation radius. The semigeostrophic solution compares favourably with numerical solutions of a single-layer shallow-water model that internally develops a leading bore. Existing laboratory experiments are re-analysed and some new experiments are undertaken. Comparisons are also made with a three-dimensional shallow-water model. These show that lateral boundary friction is the primary reason for differences between the experiments and the semigeostrophic theory. The wall no-slip condition is identified as the primary cause of the experimentally observed decrease in gravity current speed with time. A model for the viscous decay is developed and shown to agree with both experimental and numerical model data.This work was supported by NSF Grants OCE-0095059 and OCE-0132903

Concise theory of chiral lipid membranes

A theory of chiral lipid membranes is proposed on the basis of a concise free energy density which includes the contributions of the bending and the surface tension of membranes, as well as the chirality and orientational variation of tilting molecules. This theory is consistent with the previous experiments [J.M. Schnur \textit{et al.}, Science \textbf{264}, 945 (1994); M.S. Spector \textit{et al.}, Langmuir \textbf{14}, 3493 (1998); Y. Zhao, \textit{et al.}, Proc. Natl. Acad. Sci. USA \textbf{102}, 7438 (2005)] on self-assembled chiral lipid membranes of DC$_{8,9}$PC. A torus with the ratio between its two generated radii larger than $\sqrt{2}$ is predicted from the Euler-Lagrange equations. It is found that tubules with helically modulated tilting state are not admitted by the Euler-Lagrange equations, and that they are less energetically favorable than helical ripples in tubules. The pitch angles of helical ripples are theoretically estimated to be about 0$^\circ$ and 35$^\circ$, which are close to the most frequent values 5$^\circ$ and 28$^\circ$ observed in the experiment [N. Mahajan \textit{et al.}, Langmuir \textbf{22}, 1973 (2006)]. Additionally, the present theory can explain twisted ribbons of achiral cationic amphiphiles interacting with chiral tartrate counterions. The ratio between the width and pitch of twisted ribbons is predicted to be proportional to the relative concentration difference of left- and right-handed enantiomers in the low relative concentration difference region, which is in good agreement with the experiment [R. Oda \textit{et al.}, Nature (London) \textbf{399}, 566 (1999)].Comment: 14 pages, 7 figure

Effective free energy for pinned membranes

We consider membranes adhered through specific receptor-ligand bonds. Thermal undulations of the membrane induce effective interactions between adhesion sites. We derive an upper bound to the free energy that is independent of interaction details. To lowest order in a systematic expansion we obtain two-body interactions which allow to map the free energy onto a lattice gas with constant density. The induced interactions alone are not strong enough to lead to a condensation of individual adhesion sites. A measure of the thermal roughness is shown to depend on the inverse square root of the density of adhesion sites, which is in good agreement with previous computer simulations.Comment: to appear as a Rapid Communication in Phys. Rev.

n-atic Order and Continuous Shape Changes of Deformable Surfaces of Genus Zero

We consider in mean-field theory the continuous development below a second-order phase transition of $n$-atic tangent plane order on a deformable surface of genus zero with order parameter $\psi = \langle e^{i n \theta} \rangle$. Tangent plane order expels Gaussian curvature. In addition, the total vorticity of orientational order on a surface of genus zero is two. Thus, the ordered phase of an $n$-atic on such a surface will have $2n$ vortices of strength $1/n$, $2n$ zeros in its order parameter, and a nonspherical equilibrium shape. Our calculations are based on a phenomenological model with a gauge-like coupling between $\psi$ and curvature, and our analysis follows closely the Abrikosov treatment of a type II superconductor just below $H_{c2}$.Comment: REVTEX, 12 page

Observation of an Efimov resonance in an ultracold mixture of atoms and weakly bound dimers

We discuss our recent observation of an atom-dimer Efimov resonance in an ultracold mixture of Cs atoms and Cs_2 Feshbach molecules [Nature Phys. 5, 227 (2009)]. We review our experimental procedure and present additional data involving a non-universal g-wave dimer state, to contrast our previous results on the universal s-wave dimer. We resolve a seeming discrepancy when quantitatively comparing our experimental findings with theoretical results from effective field theory.Comment: Conference Proceeding ICPEAC 2009 Kalamazoo, to appear in Journal of Physics: Conference Serie