Dynamics of a multicomponent vesicle in shear flow

Hydrodynamics of a multicomponent vesicle under strong confinement

We numerically investigate the hydrodynamics and membrane dynamics of a multicomponent vesicle in two strongly confined geometries. This serves as a simplified model for red blood cells undergoing large deformations while traversing narrow constrictions. We propose a new parameterization for the bending modulus that remains positive for all lipid phase parameter values. For a multicomponent vesicle passing through a stenosis, we establish connections between various properties: lipid phase coarsening, size and flow profile of the lubrication layers, excess pressure, and the tank-treading velocity of the membrane. For a multicomponent vesicle passing through a contracting channel, we find that the lipid always phase separates so that the vesicle is stiffer in the front as it passes through the constriction. For both cases of confinement we find that lipid coarsening is arrested under strong confinement, and resumes at a high rate upon relief from extreme confinement. The results may be useful for efficient sorting lipid domains using microfluidic flows by controlled release of vesicles passing through strong confinement.

Lift at low Reynolds number

Lift forces are widespread in hydrodynamics. These are typically observed for big and fast objects and are often associated with a combination of fluid inertia (i.e. large Reynolds numbers) and specific symmetry-breaking mechanisms. In contrast, the properties of viscosity-dominated (i.e. low Reynolds numbers) flows make it more difficult for such lift forces to emerge. However, the inclusion of boundary effects qualitatively changes this picture. Indeed, in the context of soft and biological matter, recent studies have revealed the emergence of novel lift forces generated by boundary softness, flow gradients and/or surface charges. The aim of the present review is to gather and analyse this corpus of literature, in order to identify and unify the questioning within the associated communities, and pave the way towards future research.

Phase field modeling and computation of vesicle growth or shrinkage

We present a phase field model for vesicle growth or shrinkage induced by an osmotic pressure due to a chemical potential gradient. The model consists of an Allen-Cahn equation describing the evolution of the phase field parameter that describes the shape of the vesicle and a Cahn-Hilliard-type equation describing the evolution of the ionic fluid. We establish conditions for vesicle growth or shrinkage via a common tangent construction using free energy curves. During the membrane deformation, the model ensures total mass conservation of the ionic fluid, and we weakly enforce a surface area constraint of the vesicle. We develop a stable numerical scheme and an efficient nonlinear multigrid solver to evolve the phase and concentration fields, and we use this to evolve the fields to near equilibrium for 2D vesicles. Convergence tests confirm an [Formula: see text] accuracy for our scheme and near-optimal convergence for our multigrid solver. Numerical results reveal that the diffuse interface model captures the main features of cell shape dynamics: for a growing vesicle, there exist circle-like equilibrium shapes if the concentration difference across the membrane and the initial osmotic pressure are large enough; while for a shrinking vesicle, there exists a rich collection of finger-like equilibrium morphologies.

Three-dimensional numerical study on wrinkling of vesicles in elongation flow based on the immersed boundary method

We study the wrinkling dynamics of three-dimensional vesicles in a time-dependent elongation flow by utilizing an immersed boundary method. For a quasispherical vesicle, our numerical results well match the predictions of perturbation analysis, where similar exponential relationships between wrinkles' characteristic wavelength and the flow strength are observed. Using the same parameters as in the experiments by Kantsler et al. [V. Kantsler et al., Phys. Rev. Lett. 99, 178102 (2007)0031-900710.1103/PhysRevLett.99.178102], our simulations of an elongated vesicle are in good agreement with their results. In addition, we get rich three-dimensional morphological details, which are favorable to comprehend the two-dimensional snapshots. This morphological information helps identify wrinkle patterns. We analyze the morphological evolution of wrinkles using spherical harmonics. We find discrepancies in elongated vesicle dynamics between simulations and perturbation analysis, highlighting the importance of the nonlinear effects. Finally, we investigate the unevenly distributed local surface tension, which largely determines the position of wrinkles excited on the vesicle membrane.

Protein drift-diffusion dynamics and phase separation in curved cell membranes and dendritic spines: Hybrid discrete-continuum methods

We develop methods for investigating protein drift-diffusion dynamics in heterogeneous cell membranes and the roles played by geometry, diffusion, chemical kinetics, and phase separation. Our hybrid stochastic numerical methods combine discrete particle descriptions with continuum-level models for tracking the individual protein drift-diffusion dynamics when coupled to continuum fields. We show how our approaches can be used to investigate phenomena motivated by protein kinetics within dendritic spines. The spine geometry is hypothesized to play an important biological role regulating synaptic strength, protein kinetics, and self-assembly of clusters. We perform simulation studies for model spine geometries varying the neck size to investigate how phase-separation and protein organization is influenced by different shapes. We also show how our methods can be used to study the roles of geometry in reaction-diffusion systems including Turing instabilities. Our methods provide general approaches for investigating protein kinetics and drift-diffusion dynamics within curved membrane structures.

Mode-coupling approach to near-cuspidal patterns in planar fluid flows

We investigate the evolution of the interface separating two Newtonian fluids of different viscosities in two-dimensional Stokes flow driven by suction or injection. A second-order, mode-coupling theory is used to explore key morphological aspects of the emerging interfacial patterns in the stage of the flow that bridges the purely linear and fully nonlinear regimes. In the linear regime, our analysis reveals that an injection-driven expanding interface is stable, while a contracting motion driven by suction is unstable. Moreover, we find that the linear growth rate associated with this suction-driven instability is independent of the viscosity contrast between the fluids. However, second-order results tell a different story, and show that the viscosity contrast is crucial in determining the morphology of the interface. Our theoretical description is applicable to the entire range of viscosity contrasts, and provides insights on the formation of near-cusp pattern-forming structures. Reproduction of fully nonlinear, n-fold symmetric near-cuspidal shapes previously obtained through conformal mapping techniques substantiates the validity of our mode-coupling approach.

Electrohydrodynamics of Vesicles and Capsules

Giant unilamellar vesicles (GUVs) made up of phospholipid bilayer membranes (liposomes) and elastic capsules with a cross-linked, polymerized membrane, have emerged as biomimetic alternatives to investigating biological cells such as leukocytes and erythrocytes. This feature article looks at the similarities and differences in the electrohydrodynamics (EHD) of vesicles and capsules under electric fields that determines their electromechanical response. The physics of EHD is illustrated through several examples such as the electrodeformation of single and compound, spherical and cylindrical, and charged and uncharged vesicles in uniform and nonuniform electric fields, and the relevance and challenges are discussed. Both small and large deformation results are discussed. The use of EHD in understanding complex interfacial kinetics in capsules and the synthesis of nonspherical capsules using electric fields are also presented. Finally, the review looks at the large electrodeformation of water-in-water capsules and the relevance of constitutive laws in their response.

Equilibrium structure and deformation response of 2D kinetoplast sheets

The considerable interest in two-dimensional (2D) materials and complex molecular topologies calls for a robust experimental system for single-molecule studies. In this work, we study the equilibrium properties and deformation response of a complex DNA structure called a kinetoplast, a 2D network of thousands of linked rings akin to molecular chainmail. Examined in good solvent conditions, kinetoplasts appear as a wrinkled hemispherical sheet. The conformation of each kinetoplast is dictated by its network topology, giving it a unique shape, which undergoes small-amplitude thermal fluctuations at subsecond timescales, with a wide separation between fluctuation and diffusion timescales. They deform elastically when weakly confined and swell to their equilibrium dimensions when the confinement is released. We hope that, in the same way that linear DNA became a canonical model system on the first investigations of its polymer-like behavior, kinetoplasts can serve that role for 2D and catenated polymer systems.

Three-dimensional multicomponent vesicles: dynamics and influence of material properties

In this work, the nonlinear hydrodynamics of a three-dimensional multicomponent vesicle in shear flow are explored. Using a volume- and area-conserving projection method coupled to a gradient-augmented level set and surface phase field method, the dynamics are systematically studied as a function of the membrane bending rigidity difference between the components, the speed of diffusion compared to the underlying shear flow, and the strength of the phase domain energy compared to the bending energy. Using a pre-segregated vesicle, three dynamics are observed: stationary phase, phase-treading, and a new dynamic called vertical banding. These regimes are very sensitive to the strength of the domain line energy, as the vertical banding regime is not observed when the line energy is larger than the bending energy. The findings demonstrate that a complete understanding of multicomponent vesicle dynamics requires that the full three-dimensional system be modeled, and show the complexity obtained when considering heterogeneous material properties.