Controlling the shape and topology of two-component colloidal membranes.

Changes in the geometry and topology of self-assembled membranes underlie diverse processes across cellular biology and engineering. Similar to lipid bilayers, monolayer colloidal membranes have in-plane fluid-like dynamics and out-of-plane bending elasticity. Their open edges and micrometer-length scale provide a tractable system to study the equilibrium energetics and dynamic pathways of membrane assembly and reconfiguration. Here, we find that doping colloidal membranes with short miscible rods transforms disk-shaped membranes into saddle-shaped surfaces with complex edge structures. The saddle-shaped membranes are well approximated by Enneper's minimal surfaces. Theoretical modeling demonstrates that their formation is driven by increasing the positive Gaussian modulus, which in turn, is controlled by the fraction of short rods. Further coalescence of saddle-shaped surfaces leads to diverse topologically distinct structures, including shapes similar to catenoids, trinoids, four-noids, and higher-order structures. At long timescales, we observe the formation of a system-spanning, sponge-like phase. The unique features of colloidal membranes reveal the topological transformations that accompany coalescence pathways in real time. We enhance the functionality of these membranes by making their shape responsive to external stimuli. Our results demonstrate a pathway toward control of thin elastic sheets' shape and topology-a pathway driven by the emergent elasticity induced by compositional heterogeneity.

Flow states of two dimensional active gels driven by external shear.

Using a minimal hydrodynamic model, we theoretically and computationally study the Couette flow of active gels in straight and annular two-dimensional channels subject to an externally imposed shear. The gels are isotropic in the absence of externally- or activity-driven shear, but have nematic order that increases with shear rate. Using the finite element method, we determine the possible flow states for a range of activities and shear rates. Linear stability analysis of an unconfined gel in a straight channel shows that an externally imposed shear flow can stabilize an extensile fluid that would be unstable to spontaneous flow in the absence of the shear flow, and destabilize a contractile fluid that would be stable against spontaneous flow in the absence of shear flow. These results are in rough agreement with the stability boundaries between the base shear flow state and the nonlinear flow states that we find numerically for a confined active gel. For extensile fluids, we find three kinds of nonlinear flow states in the range of parameters we study: unidirectional flows, oscillatory flows, and dancing flows. To highlight the activity-driven spontaneous component of the nonlinear flows, we characterize these states by the average volumetric flow rate and the wall stress. For contractile fluids, we only find the linear shear flow and a nonlinear unidirectional flow in the range of parameters that we studied. For large magnitudes of the activity, the unidirectional contractile flow develops a boundary layer. Our analysis of annular channels shows how curvature of the streamlines in the base flow affects the transitions among flow states.

Chiral fluid membranes with orientational order and multiple edges.

We carry out Monte Carlo simulations on fluid membranes with orientational order and multiple edges in the presence and absence of external forces. The membrane resists bending and has an edge tension, the orientational order couples with the membrane surface normal through a cost for tilting, and there is a chiral liquid crystalline interaction. In the absence of external forces, a membrane initialized as a vesicle will form a disk at low chirality, with the directors forming a smectic-A phase with alignment perpendicular to the membrane surface except near the edge. At large chirality a catenoid-like shape or a trinoid-like shape is formed, depending on the number of edges in the initial vesicle. This shape change is accompanied by cholesteric ordering of the directors and multiple π walls connecting the membrane edges and wrapping around the membrane neck. If the membrane is initialized instead in a cylindrical shape and stretched by an external force, it maintains a nearly cylindrical shape but additional liquid crystalline phases appear. For large tilt coupling and low chirality, a smectic-A phase forms where the directors are normal to the surface of the membrane. For lower values of the tilt coupling, a nematic phase appears at zero chirality with the average director oriented perpendicular to the long axis of the membrane, while for nonzero chirality a cholesteric phase appears. The π walls are tilt walls at low chirality and transition to twist walls as chirality is increased. We construct a continuum model of the director field to explain this behavior.

Deformation and orientational order of chiral membranes with free edges.

Motivated by experiments on colloidal membranes composed of chiral rod-like viruses, we use Monte Carlo methods to simulate these systems and determine the phase diagram for the liquid crystalline order of the rods and the membrane shape. We generalize the Lebwohl-Lasher model for a nematic with a chiral coupling to a curved surface with edge tension and a resistance to bending, and include an energy cost for tilting of the rods relative to the local membrane normal. The membrane is represented by a triangular mesh of hard beads joined by bonds, where each bead is decorated by a director. The beads can move, the bonds can reconnect and the directors can rotate at each Monte Carlo step. When the cost of tilt is small, the membrane tends to be flat, with the rods only twisting near the edge for low chiral coupling, and remaining parallel to the normal in the interior of the membrane. At high chiral coupling, the rods twist everywhere, forming a cholesteric state. When the cost of tilt is large, the emergence of the cholesteric state at high values of the chiral coupling is accompanied by the bending of the membrane into a saddle shape. Increasing the edge tension tends to flatten the membrane. These results illustrate the geometric frustration arising from the inability of a surface normal to have twist.

Axisymmetric membranes with edges under external force: buckling, minimal surfaces, and tethers.

We use theory and numerical computation to determine the shape of an axisymmetric fluid membrane with a resistance to bending and constant area. The membrane connects two rings in the classic geometry that produces a catenoidal shape in a soap film. In our problem, we find infinitely many branches of solutions for the shape and external force as functions of the separation of the rings, analogous to the infinite family of eigenmodes for the Euler buckling of a slender rod. Special attention is paid to the catenoid, which emerges as the shape of maximal allowable separation when the area is less than a critical area equal to the planar area enclosed by the two rings. A perturbation theory argument directly relates the tension of catenoidal membranes to the stability of catenoidal soap films in this regime. When the membrane area is larger than the critical area, we find additional cylindrical tether solutions to the shape equations at large ring separation, and that arbitrarily large ring separations are possible. These results apply for the case of vanishing Gaussian curvature modulus; when the Gaussian curvature modulus is nonzero and the area is below the critical area, the force and the membrane tension diverge as the ring separation approaches its maximum value. We also examine the stability of our shapes and analytically show that catenoidal membranes have markedly different stability properties than their soap film counterparts.

Force-Induced Formation of Twisted Chiral Ribbons.

We demonstrate that an achiral stretching force transforms disk-shaped colloidal membranes composed of chiral rods into twisted ribbons with handedness opposite the preferred twist of the rods. Using an experimental technique that enforces torque-free boundary conditions we simultaneously measure the force-extension curve and the ribbon shape. An effective theory that accounts for the membrane bending energy and uses geometric properties of the edge to model the internal liquid crystalline degrees of freedom explains both the measured force-extension curve and the force-induced twisted shape.

Achiral symmetry breaking and positive Gaussian modulus lead to scalloped colloidal membranes.

In the presence of a nonadsorbing polymer, monodisperse rod-like particles assemble into colloidal membranes, which are one-rod-length-thick liquid-like monolayers of aligned rods. Unlike 3D edgeless bilayer vesicles, colloidal monolayer membranes form open structures with an exposed edge, thus presenting an opportunity to study elasticity of fluid sheets. Membranes assembled from single-component chiral rods form flat disks with uniform edge twist. In comparison, membranes composed of a mixture of rods with opposite chiralities can have the edge twist of either handedness. In this limit, disk-shaped membranes become unstable, instead forming structures with scalloped edges, where two adjacent lobes with opposite handedness are separated by a cusp-shaped point defect. Such membranes adopt a 3D configuration, with cusp defects alternatively located above and below the membrane plane. In the achiral regime, the cusp defects have repulsive interactions, but away from this limit we measure effective long-ranged attractive binding. A phenomenological model shows that the increase in the edge energy of scalloped membranes is compensated by concomitant decrease in the deformation energy due to Gaussian curvature associated with scalloped edges, demonstrating that colloidal membranes have positive Gaussian modulus. A simple excluded volume argument predicts the sign and magnitude of the Gaussian curvature modulus that is in agreement with experimental measurements. Our results provide insight into how the interplay between membrane elasticity, geometrical frustration, and achiral symmetry breaking can be used to fold colloidal membranes into 3D shapes.

Stability of the interface of an isotropic active fluid.

We study the linear stability of an isotropic active fluid in three different geometries: a film of active fluid on a rigid substrate, a cylindrical thread of fluid, and a spherical fluid droplet. The active fluid is modeled by the hydrodynamic theory of an active nematic liquid crystal in the isotropic phase. In each geometry, we calculate the growth rate of sinusoidal modes of deformation of the interface. There are two distinct branches of growth rates; at long wavelength, one corresponds to the deformation of the interface, and one corresponds to the evolution of the liquid crystalline degrees of freedom. The passive cases of the film and the spherical droplet are always stable. For these geometries, a sufficiently large activity leads to instability. Activity also leads to propagating damped or growing modes. The passive cylindrical thread is unstable for perturbations with wavelength longer than the circumference. A sufficiently large activity can make any wavelength unstable, and again leads to propagating damped or growing modes. Our calculations are carried out for the case of zero Frank elasticity. While Frank elasticity is a stabilizing mechanism as it penalizes distortions of the order parameter tensor, we show that it has a small effect on the instabilities considered here.

Enhancement of Microorganism Swimming Speed in Active Matter.

We study a swimming undulating sheet in the isotropic phase of an active nematic liquid crystal. Activity changes the effective shear viscosity, reducing it to zero at a critical value of activity. Expanding in the sheet amplitude, we find that the correction to the swimming speed due to activity is inversely proportional to the effective shear viscosity. Our perturbative calculation becomes invalid near the critical value of activity; using numerical methods to probe this regime, we find that activity enhances the swimming speed by an order of magnitude compared to the passive case.

Helical motion of the cell body enhances Caulobacter crescentus motility.

We resolve the 3D trajectory and the orientation of individual cells for extended times, using a digital tracking technique combined with 3D reconstructions. We have used this technique to study the motility of the uniflagellated bacterium Caulobacter crescentus and have found that each cell displays two distinct modes of motility, depending on the sense of rotation of the flagellar motor. In the forward mode, when the flagellum pushes the cell, the cell body is tilted with respect to the direction of motion, and it precesses, tracing out a helical trajectory. In the reverse mode, when the flagellum pulls the cell, the precession is smaller and the cell has a lower translation distance per rotation period and thus a lower motility. Using resistive force theory, we show how the helical motion of the cell body generates thrust and can explain the direction-dependent changes in swimming motility. The source of the cell body precession is believed to be associated with the flexibility of the hook that connects the flagellum to the cell body.

Minimal model for transient swimming in a liquid crystal.

When a microorganism begins swimming from rest in a Newtonian fluid such as water, it rapidly attains its steady-state swimming speed since changes in the velocity field spread quickly when the Reynolds number is small. However, swimming microorganisms are commonly found or studied in complex fluids. Because these fluids have long relaxation times, the time to attain the steady-state swimming speed can also be long. In this article we study the swimming startup problem in the simplest liquid crystalline fluid: a two-dimensional hexatic liquid crystal film. We study the dependence of startup time on anchoring strength and Ericksen number, which is the ratio of viscous to elastic stresses. For strong anchoring, the fluid flow starts up immediately but the liquid crystal field and swimming velocity attain their sinusoidal steady-state values after a time proportional to the relaxation time of the liquid crystal. When the Ericksen number is high, the behavior is the same as in the strong-anchoring case for any anchoring strength. We also find that the startup time increases with the ratio of the rotational viscosity to the shear viscosity, and then ultimately saturates once the rotational viscosity is much greater than the shear viscosity.

Bacterial rheotaxis.

The motility of organisms is often directed in response to environmental stimuli. Rheotaxis is the directed movement resulting from fluid velocity gradients, long studied in fish, aquatic invertebrates, and spermatozoa. Using carefully controlled microfluidic flows, we show that rheotaxis also occurs in bacteria. Excellent quantitative agreement between experiments with Bacillus subtilis and a mathematical model reveals that bacterial rheotaxis is a purely physical phenomenon, in contrast to fish rheotaxis but in the same way as sperm rheotaxis. This previously unrecognized bacterial taxis results from a subtle interplay between velocity gradients and the helical shape of flagella, which together generate a torque that alters a bacterium's swimming direction. Because this torque is independent of the presence of a nearby surface, bacterial rheotaxis is not limited to the immediate neighborhood of liquid-solid interfaces, but also takes place in the bulk fluid. We predict that rheotaxis occurs in a wide range of bacterial habitats, from the natural environment to the human body, and can interfere with chemotaxis, suggesting that the fitness benefit conferred by bacterial motility may be sharply reduced in some hydrodynamic conditions.

Evidence for two extremes of ciliary motor response in a single swimming microorganism.

Because arrays of motile cilia drive fluids for a range of processes, the versatile mechano-chemical mechanism coordinating them has been under scrutiny. The protist Paramecium presents opportunities to compare how groups of cilia perform two distinct functions, swimming propulsion and nutrient uptake. We present how the body cilia responsible for propulsion and the oral-groove cilia responsible for nutrient uptake respond to changes in their mechanical environment accomplished by varying the fluid viscosity over a factor of 7. Analysis with a phenomenological model of trajectories of swimmers made neutrally buoyant with magnetic forces combined with high-speed imaging of ciliary beating reveal that the body cilia exert a nearly constant propulsive force primarily by reducing their beat frequency as viscosity increases. By contrast, the oral-groove cilia beat at a nearly constant frequency. The existence of two extremes of motor response in a unicellular organism prompts unique investigations of factors controlling ciliary beating.

Mutual diffusion of inclusions in freely suspended smectic liquid crystal films.

We study experimentally and theoretically the hydrodynamic interaction of pairs of circular inclusions in two-dimensional, fluid smectic membranes suspended in air. By analyzing their Brownian motion, we find that the radial mutual mobilities of identical inclusions are independent of their size but that the angular coupling becomes strongly size dependent when their radius exceeds a characteristic hydrodynamic length. These observations are described well for arbitrary inclusion separations by a model that generalizes the Levine-MacKintosh theory of point-force response functions and uses a boundary-element approach to calculate the mobility matrix for inclusions of finite extent.

Force-free swimming of a model helical flagellum in viscoelastic fluids.

We precisely measure the force-free swimming speed of a rotating helix in viscous and viscoelastic fluids. The fluids are highly viscous to replicate the low Reynolds number environment of microorganisms. The helix, a macroscopic scale model for the bacterial flagellar filament, is rigid and rotated at a constant rate while simultaneously translated along its axis. By adjusting the translation speed to make the net hydrodynamic force vanish, we measure the force-free swimming speed as a function of helix rotation rate, helix geometry, and fluid properties. We compare our measurements of the force-free swimming speed of a helix in a high-molecular weight silicone oil with predictions for the swimming speed in a Newtonian fluid, calculated using slender-body theories and a boundary-element method. The excellent agreement between theory and experiment in the Newtonian case verifies the high accuracy of our experiments. For the viscoelastic fluid, we use a polymer solution of polyisobutylene dissolved in polybutene. This solution is a Boger fluid, a viscoselastic fluid with a shear-rate-independent viscosity. The elasticity is dominated by a single relaxation time. When the relaxation time is short compared to the rotation period, the viscoelastic swimming speed is close to the viscous swimming speed. As the relaxation time increases, the viscoelastic swimming speed increases relative to the viscous speed, reaching a peak when the relaxation time is comparable to the rotation period. As the relaxation time is further increased, the viscoelastic swimming speed decreases and eventually falls below the viscous swimming speed.

Locomotion of helical bodies in viscoelastic fluids: enhanced swimming at large helical amplitudes.

The motion of a rotating helical body in a viscoelastic fluid is considered. In the case of force-free swimming, the introduction of viscoelasticity can either enhance or retard the swimming speed and locomotive efficiency, depending on the body geometry, fluid properties, and the body rotation rate. Numerical solutions of the Oldroyd-B equations show how previous theoretical predictions break down with increasing helical radius or with decreasing filament thickness. Helices of large pitch angle show an increase in swimming speed to a local maximum at a Deborah number of order unity. The numerical results show how the small-amplitude theoretical calculations connect smoothly to the large-amplitude experimental measurements.

Effects of viscogens on RNA transcription inside reovirus particles.

The dsRNA genome of mammalian reovirus (MRV), like the dsDNA genomes of herpesviruses and many bacteriophages, is packed inside its icosahedral capsid in liquid-crystalline form, with concentrations near or more than 400 mg/ml. Viscosity in such environments must be high, but the relevance of viscosity for the macromolecular processes occurring there remains poorly characterized. Here, we describe the use of simple viscogens, glycerol and sucrose, to examine their effects on RNA transcription inside MRV core particles. Transcription inside MRV cores was strongly inhibited by these agents and to a greater extent than either predicted by theory or exhibited by a nonencapsidated transcriptase, suggesting that RNA transcription inside MRV cores is unusually sensitive to viscogen effects. The elongation phase of transcription was found to be a primary target of this inhibition. Similar results were obtained with particles of a second dsRNA virus, rhesus rotavirus, from a divergent taxonomic subfamily. Polymeric viscogens such as polyethylene glycol also inhibited RNA transcription inside MRV cores, but in a size-limited manner, suggesting that diffusion through channels in the MRV core is required for their activity. Modeling of the data suggested that the inherent intracapsid viscosity of both reo- and rotavirus is indeed high, two to three times the viscosity of water. The capacity for quantitative comparisons of intracapsid viscosity and effects of viscogens on macromolecular processes in confined spaces should be similarly informative in other systems.

Shapes of fluid membranes with chiral edges.

We carry out Monte Carlo simulations of a colloidal fluid membrane with a free edge and composed of chiral rodlike viruses. The membrane is modeled by a triangular mesh of beads connected by bonds in which the bonds and beads are free to move at each Monte Carlo step. Since the constituent viruses are experimentally observed to twist only near the membrane edge, we use an effective energy that favors a particular sign of the geodesic torsion of the edge. The effective energy also includes the membrane bending stiffness, edge bending stiffness, and edge tension. We find three classes of membrane shapes resulting from the competition of the various terms in the free energy: branched shapes, chiral disks, and vesicles. Increasing the edge bending stiffness smooths the membrane edge, leading to correlations among the membrane normals at different points along the edge. The normalized power spectrum for edge displacements shows a peak with increasing preferred geodesic torsion. We also consider membrane shapes under an external force by fixing the distance between two ends of the membrane and finding the shape for increasing values of the distance between the two ends. As the distance increases, the membrane twists into a ribbon, with the force eventually reaching a plateau.

Topological structure and dynamics of three-dimensional active nematics.

Topological structures are effective descriptors of the nonequilibrium dynamics of diverse many-body systems. For example, motile, point-like topological defects capture the salient features of two-dimensional active liquid crystals composed of energy-consuming anisotropic units. We dispersed force-generating microtubule bundles in a passive colloidal liquid crystal to form a three-dimensional active nematic. Light-sheet microscopy revealed the temporal evolution of the millimeter-scale structure of these active nematics with single-bundle resolution. The primary topological excitations are extended, charge-neutral disclination loops that undergo complex dynamics and recombination events. Our work suggests a framework for analyzing the nonequilibrium dynamics of bulk anisotropic systems as diverse as driven complex fluids, active metamaterials, biological tissues, and collections of robots or organisms.

Role of slip between a probe particle and a gel in microrheology.

In the technique of microrheology, macroscopic rheological parameters as well as information about local structure are deduced from the behavior of microscopic probe particles under thermal or active forcing. Microrheology requires knowledge of the relation between macroscopic parameters and the force felt by a particle in response to displacements. We investigate this response function for a spherical particle using the two-fluid model, in which the gel is represented by a polymer network coupled to a surrounding solvent via a drag force. We obtain an analytic solution for the response function in the limit of small volume fraction of the polymer network, and neglecting inertial effects. We use no-slip boundary conditions for the solvent at the surface of the sphere. The boundary condition for the network at the surface of the sphere is a kinetic friction law, for which the tangential stress of the network is proportional to relative velocity of the network and the sphere. This boundary condition encompasses both no-slip and frictionless boundary conditions as limits. Far from the sphere there is no relative motion between the solvent and network due to the coupling between them. However, the different boundary conditions on the solvent and network tend to produce different far-field motions. We show that the far-field motion and the force on the sphere are controlled by the solvent boundary conditions at high frequency and by the network boundary conditions at low frequency. At low frequencies compression of the network can also affect the force on the sphere. We find the crossover frequencies at which the effects of sliding of the sphere past the polymer network and compression of the gel become important. The effects of sliding alone can lead to an underestimation of moduli by up to 33%, while the effects of compression alone can lead to an underestimation of moduli by up to 20%, and the effects of sliding and compression combined can lead to an underestimation of moduli by up to 43%.