Wrinkling of vesicles during transient dynamics in elongational flow

Determinative scrolling and folding of membranes through shrinking channels

Flat membranes ubiquitously transform into mysterious complex shapes in nature and artificial worlds. Behind the complexity, clear determinative deformation modes have been continuously found to serve as basic application rules but remain unfulfilled. Here, we decipher two elemental deformation modes of thin membranes, spontaneous scrolling and folding as passing through shrinking channels. We validate that these two modes rule the deformation of membranes of a wide thickness range from micrometer to atomic scale. Their occurrence and the determinative fold number quantitatively correlate with the Föppl-von Kármán number and shrinkage ratio. The unveiled determinative deformation modes can guide fabricating foldable designer microrobots and delicate structures of two-dimensional sheets and provide another mechanical principle beyond genetic determinism in biological morphogens.

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.

Nonlinear multimode buckling dynamics examined with semiflexible paramagnetic filaments

We present the contractile buckling dynamics of superparamagnetic filaments using experimental, theoretical, and simulation approaches. Under the influence of an orthogonal magnetic field, flexible magnetic filaments exhibit higher-order buckling dynamics that can be identified as occurring in three stages: initiation, development, and decay. Unlike initiation and decay stages where the balance between magnetic interactions and elastic forces is dominant, in the development stage, the influence of hydrodynamic drag results in transient buckling dynamics that is nonlinear along the filament contour. The inhomogeneous temporal evolution of the buckling wavelength is analyzed and the contractions under various conditions are compared.

Wrinkling dynamics of fluctuating vesicles in time-dependent viscous flow

We study the fully nonlinear, nonlocal dynamics of two-dimensional vesicles in a time-dependent, incompressible viscous flow at finite temperature. We focus on a transient instability that can be observed when the direction of applied flow is suddenly reversed, which induces compressive forces on the vesicle interface, and small-scale interface perturbations known as wrinkles develop. These wrinkles are driven by regions of negative elastic tension on the membrane. Using a stochastic immersed boundary method with a biophysically motivated choice of thermal fluctuations, we investigate the wrinkling dynamics numerically. Different from deterministic wrinkling dynamics, thermal fluctuations lead to symmetry-breaking wrinkling patterns by exciting higher order modes. This leads to more rapid and more realistic wrinkling dynamics. Our results are in excellent agreement with the experimental data by Kantsler et al. [Kantsler et al., Phys. Rev. Lett., 2007, 99, 17802]. We compare the nonlinear simulation results with perturbation theory, modified to account for thermal fluctuations. The strength of the applied flow strongly influences the most unstable wavelength characterizing the wrinkles, and there are significant differences between the results from perturbation theory and the fully nonlinear simulations, which suggests that the perturbation theory misses important nonlinear interactions. Strikingly, we find that thermal fluctuations actually have the ability to attenuate variability of the characteristic wavelength of wrinkling by exciting a wider range of modes than the deterministic case, which makes the evolution less constrained and enables the most unstable wavelength to emerge more readily. We further find that thermal noise helps prevent the vesicle from rotating if it is misaligned with the direction of the applied extensional flow.

Effects of thermal noise on the transitional dynamics of an inextensible elastic filament in stagnation flow

We investigate the dynamics of a single inextensible elastic filament subject to anisotropic friction in a viscous stagnation-point flow, by employing both a continuum model represented by Langevin type stochastic partial differential equations (SPDEs) and a dissipative particle dynamics (DPD) method. Unlike previous works, the filament is free to rotate and the tension along the filament is determined by the local inextensible constraint. The kinematics of the filament is recorded and studied with normal modes analysis. The results show that the filament displays an instability induced by negative tension, which is analogous to Euler buckling of a beam. Symmetry breaking of normal modes dynamics and stretch-coil transitions are observed above the threshold of the buckling instability point. Furthermore, both temporal and spatial noise are amplified resulting from the interaction of thermal fluctuations and nonlinear filament dynamics. Specifically, the spatial noise is amplified with even normal modes being excited due to symmetry breaking, while the temporal noise is amplified with increasing time correlation length and variance.

Fluid vesicles in flow

We review the dynamical behavior of giant fluid vesicles in various types of external hydrodynamic flow. The interplay between stresses arising from membrane elasticity, hydrodynamic flows, and the ever present thermal fluctuations leads to a rich phenomenology. In linear flows with both rotational and elongational components, the properties of the tank-treading and tumbling motions are now well described by theoretical and numerical models. At the transition between these two regimes, strong shape deformations and amplification of thermal fluctuations generate a new regime called trembling. In this regime, the vesicle orientation oscillates quasi-periodically around the flow direction while asymmetric deformations occur. For strong enough flows, small-wavelength deformations like wrinkles are observed, similar to what happens in a suddenly reversed elongational flow. In steady elongational flow, vesicles with large excess areas deform into dumbbells at large flow rates and pearling occurs for even stronger flows. In capillary flows with parabolic flow profile, single vesicles migrate towards the center of the channel, where they adopt symmetric shapes, for two reasons. First, walls exert a hydrodynamic lift force which pushes them away. Second, shear stresses are minimal at the tip of the flow. However, symmetry is broken for vesicles with large excess areas, which flow off-center and deform asymmetrically. In suspensions, hydrodynamic interactions between vesicles add up to these two effects, making it challenging to deduce rheological properties from the dynamics of individual vesicles. Further investigations of vesicles and similar objects and their suspensions in steady or time-dependent flow will shed light on phenomena such as blood flow.

Membrane compression in tumbling and vacillating-breathing regimes for quasispherical vesicles

We derive some analytical results of a well-known model for quasispherical vesicles in a linear shear flow at low deformability. Attention is focussed on the oscillatory regimes: the tumbling (TB) mode, vacillating-breathing (VB) mode, and the transition from vacillating-breathing to tumbling, depending on a control parameter Γ. It is shown that, during the VB-to-TB transition (Γ=1), the vesicle momentarily attains its maximal extension in the vorticity direction and transits through a circular profile in the shear plane for which the radius is exactly determined. In addition, we provide an explicit analytical expression for the effective membrane tension for different types of motions. We find a critical bending number below which the membrane undergoes compression at each instant and show that, during the VB-to-TB transition, a fourth-order membrane deformation is possible.

Noisy nonlinear dynamics of vesicles in flow

We present a model for the dynamics of fluid vesicles in linear flow which consistently includes thermal fluctuations and nonlinear coupling between different modes. At the transition between tank treading and tumbling, we predict a trembling motion which is at odds with the known deterministic motions and for which thermal noise is strongly amplified. In particular, highly asymmetric shapes are observed even though the deterministic flow only allows for axisymmetric ones. Our results explain quantitatively recent experimental observations [Levant and Steinberg, Phys. Rev. Lett. 109, 268103 (2012)].

Amplification of thermal noise by vesicle dynamics

A novel noise amplification mechanism resulting from the interaction of thermal fluctuations and nonlinear vesicle dynamics is reported. It is observed in a time-dependent vesicle state called trembling (TR). High spatial resolution and very long time series of TR compared to the vesicle period allow us to quantitatively analyze the generation and amplification of spatial and temporal modes of the vesicle shape perturbations. During a compression part of each TR cycle, a vesicle finds itself on the edge of the wrinkling instability, where thermally excited spatial modes are amplified.

Fluctuations, dynamics, and the stretch-coil transition of single actin filaments in extensional flows

Semiflexible polymers subject to hydrodynamic forcing play an important role in cytoskeletal motions in the cell, particularly when filaments guide molecular motors whose motions create flows. Near hyperbolic stagnation points, filaments experience a competition between bending elasticity and tension and are predicted to display suppressed thermal fluctuations in the extensional regime and a buckling instability under compression. Using a microfluidic cross-flow geometry, we verify these predictions in detail, including a fluctuation-rounded stretch-coil transition of actin filaments.

Dynamics of red blood cells and vesicles in microchannels of oscillating width

We have studied the dynamics of red blood cells and fluid lipid vesicles in hydrodynamic flow fields created by microchannels with periodically varying channel width. For red blood cells we find a transition from a regime with oscillating tilt angle and fixed shape to a regime with oscillating shape with increasing flow velocity. We have determined the crossover to occur at a critical ratio L(y)/v(m) ≈ 2.2 × 10⁻³ s with channel width L(y) and red blood cell velocity v(m). These oscillations are superposed by shape transitions from a discocyte to a slipper shape at low velocities and a slipper to parachute transition at high flow velocities.

Stretching and relaxation of vesicles

We study the shape relaxation of spherical giant unilamellar vesicles which have been deformed far from equilibrium into ellipsoids using optical tweezers. The relaxation back to a sphere is determined by elastic constants of the vesicles, and their excess area, parameters that are obtained for each stretched vesicle from shape fluctuations in thermal equilibrium, as well as low Reynolds number fluid flow. The relaxation time could be compared favorably to a simple formula which encompasses the joint effect of membrane rigidity and fluid flow. The time constant of the stretched vesicle is slower than that of its thermal fluctuations, which agrees with a recent theory; however, it is one order of magnitude faster than predicted.

Forming transmembrane channels using end-functionalized nanotubes

Using dissipative particle dynamics (DPD) simulations, we examine the interaction between amphiphilic nanotubes and lipid bilayer membranes. The nanotubes are represented by a hydrophobic shaft that is end-functionalized with hydrophilic groups. Nanotubes that are capped by a monolayer of hydrophilic beads or also encompass hydrophilic "hairs" on just one end of the shaft are found to spontaneously penetrate and assume a transmembrane position; the process, however, depends critically on the membrane tension. On the other hand, nanotubes that include hydrophilic hairs at both ends of the hydrophobic shaft are not observed to spontaneously self-organize into the bilayer. When the membrane is stretched to form a pore, the nanotubes with two hairy ends adsorb on the edge of the pore and become localized in the membrane, thus forming a transmembrane channel. The findings from these studies provide guidelines for creating biomimetic nanotube channels that are capable of selectively transporting molecules through the membrane in response to changes in the local environment.

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. The thermal fluctuations can induce transitions between two orbits at very low shear amplitudes. 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, synthetic capsules, and lipid vesicles.

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.

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.

Dynamical regimes and hydrodynamic lift of viscous vesicles under shear

The dynamics of two-dimensional viscous vesicles in shear flow, with different fluid viscosities etain and etaout 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 etain/etaout. 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 ycm of the vesicle center of mass from the wall is well described by an effective power law ycm(-2) for intermediate distances 0.8Rp< approximately ycm< approximately 3Rp with vesicle radius Rp. The boundary-integral calculation indicates that the lift force decays asymptotically as 1/[ycm ln(ycm)] far from the wall.

Dynamics of a vesicle in general flow

An approach to quantitatively study vesicle dynamics as well as biologically-related micro-objects in a fluid flow, which is based on the combination of a dynamical trap and a control parameter, the ratio of the vorticity to the strain rate, is suggested. The flow is continuously varied between rotational, shearing, and elongational in a microfluidic 4-roll mill device, the dynamical trap, that allows scanning of the entire phase diagram of motions, i.e., tank-treading (TT), tumbling (TU), and trembling (TR), using a single vesicle even at lambda = eta(in)/eta(out) = 1, where eta(in) and eta(out) are the viscosities of the inner and outer fluids. This cannot be achieved in pure shear flow, where the transition between TT and either TU or TR is attained only at lambda>1. As a result, it is found that the vesicle dynamical states in a general are presented by the phase diagram in a space of only 2 dimensionless control parameters. The findings are in semiquantitative accord with the recent theory made for a quasi-spherical vesicle, although vesicles with large deviations from spherical shape were studied experimentally. The physics of TR is also uncovered.