Interface Roughening in Nonequilibrium Phase-Separated Systems

Interfaces of phase-separated systems roughen in time due to capillary waves. Because of fluxes in the bulk, their dynamics is nonlocal in real space and is not described by the Edwards-Wilkinson or Kardar-Parisi-Zhang (KPZ) equations, nor their conserved counterparts. We show that, in the absence of detailed balance, the phase-separated interface is described by a new universality class that we term |q|KPZ. We compute the associated scaling exponents via one-loop renormalization group and corroborate the results by numerical integration of the |q|KPZ equation. Deriving the effective interface dynamics from a minimal field theory of active phase separation, we finally argue that the |q|KPZ universality class generically describes liquid-vapor interfaces in two- and three-dimensional active systems.

Classical Nucleation Theory for Active Fluid Phase Separation

Classical Nucleation Theory (CNT), linking rare nucleation events to the free-energy landscape of a growing nucleus, is central to understanding phase-change kinetics in passive fluids. Nucleation in nonequilibrium systems is much harder to describe because there is no free energy, but instead a dynamics-dependent quasipotential that typically must be found numerically. Here we extend CNT to a class of active phase-separating systems governed by a minimal field-theoretic model (Active Model B+). In the small noise and supersaturation limits that CNT assumes, we compute analytically the quasipotential, and hence, nucleation barrier, for liquid-vapor phase separation. Crucial to our results, detailed balance, although broken microscopically by activity, is restored along the instanton trajectory, which in CNT involves the nuclear radius as the sole reaction coordinate.

The Royal Society RAMP modelling initiative

Normally, science proceeds following a well-established set of principles. Studies are done with an emphasis on correctness, are submitted to a journal editor who evaluates their relevance, and then undergo anonymous peer review by experts before publication in a journal and acceptance by the scientific community via the open literature. This process is slow, but its accuracy has served all fields of science well. In an emergency situation, different priorities come to the fore. Research and review need to be conducted quickly, and the target audience consists of policymakers. Scientists must jostle for the attention of non-specialists without sacrificing rigour, and must deal not only with peer assessment but also with media scrutiny by journalists who may have agendas other than ensuring scientific correctness. Here, we describe how the Royal Society coordinated efforts of diverse scientists to help model the coronavirus epidemic. This article is part of the theme issue 'Technical challenges of modelling real-life epidemics and examples of overcoming these'.

Erratum: Capillary Interfacial Tension in Active Phase Separation [Phys. Rev. Lett. 127, 068001 (2021)]

This corrects the article DOI: 10.1103/PhysRevLett.127.068001.

Capillary Interfacial Tension in Active Phase Separation

In passive fluid-fluid phase separation, a single interfacial tension sets both the capillary fluctuations of the interface and the rate of Ostwald ripening. We show that these phenomena are governed by two different tensions in active systems, and compute the capillary tension σ_{cw} which sets the relaxation rate of interfacial fluctuations in accordance with capillary wave theory. We discover that strong enough activity can cause negative σ_{cw}. In this regime, depending on the global composition, the system self-organizes, either into a microphase-separated state in which coalescence is highly inhibited, or into an "active foam" state. Our results are obtained for Active Model B+, a minimal continuum model which, although generic, admits significant analytical progress.

Constitutive Model for Time-Dependent Flows of Shear-Thickening Suspensions

We develop a tensorial constitutive model for dense, shear-thickening particle suspensions subjected to time-dependent flow. Our model combines a recently proposed evolution equation for the suspension microstructure in rate-independent materials with ideas developed previously to explain the steady flow of shear-thickening ones, whereby friction proliferates among compressive contacts at large particle stresses. We apply our model to shear reversal, and find good qualitative agreement with particle-level, discrete-element simulations whose results we also present.

Hydrodynamically Interrupted Droplet Growth in Scalar Active Matter

Suspensions of spherical active particles often show microphase separation. At a continuum level, coupling their scalar density to fluid flow, there are two distinct explanations. Each involves an effective interfacial tension: the first mechanical (causing flow) and the second diffusive (causing Ostwald ripening). Here we show how the negative mechanical tension of contractile swimmers creates, via a self-shearing instability, a steady-state life cycle of droplet growth interrupted by division whose scaling behavior we predict. When the diffusive tension is also negative, this is replaced by an arrested regime (mechanistically distinct, but with similar scaling) where division of small droplets is prevented by reverse Ostwald ripening.

Dynamic clustering and re-dispersion in concentrated colloid-active gel composites

We study the dynamics of quasi-two-dimensional concentrated suspensions of colloidal particles in active gels by computer simulations. Remarkably, we find that activity induces a dynamic clustering of colloids even in the absence of any preferential anchoring of the active nematic director at the particle surface. When such an anchoring is present, active stresses instead compete with elastic forces and re-disperse the aggregates observed in passive colloid-liquid crystal composites. Our quasi-two-dimensional "inverse" dispersions of passive particles in active fluids (as opposed to the more common "direct" suspensions of active particles in passive fluids) provide a promising route towards the self-assembly of new soft materials.

Dynamic Vorticity Banding in Discontinuously Shear Thickening Suspensions

It has recently been argued that steady-state vorticity bands cannot arise in shear thickening suspensions because the normal stress imbalance across the interface between the bands will set up particle migrations. In this Letter, we develop a simple continuum model that couples shear thickening to particle migration. We show by linear stability analysis that homogeneous flow is unstable towards vorticity banding, as expected, in the regime of negative constitutive slope. In full nonlinear computations, we show, however, that the resulting vorticity bands are unsteady, with spatiotemporal patterns governed by stress-concentration coupling. We furthermore show that these dynamical bands also arise in direct particle simulations, in good agreement with the continuum model.

Mixtures of Blue Phase Liquid Crystal with Simple Liquids: Elastic Emulsions and Cubic Fluid Cylinders

We numerically investigate the behavior of a phase-separating mixture of a blue phase I liquid crystal with an isotropic fluid. The resulting morphology is primarily controlled by an inverse capillary number, χ, setting the balance between interfacial and elastic forces. When χ and the concentration of the isotropic component are both low, the blue phase disclination lattice templates a cubic array of fluid cylinders. For larger χ, the isotropic phase arranges primarily into liquid emulsion droplets which coarsen very slowly, rewiring the blue phase disclination lines into an amorphous elastic network. Our blue phase-simple fluid composites can be externally manipulated: an electric field can trigger a morphological transition between cubic fluid cylinder phases with different topologies.

Soft interfacial materials: from fundamentals to formulation

This article is part of the themed issue ‘Soft interfacial materials: from fundamentals to formulation’.

Viscoelastic and elastomeric active matter: Linear instability and nonlinear dynamics

We consider a continuum model of active viscoelastic matter, whereby an active nematic liquid crystal is coupled to a minimal model of polymer dynamics with a viscoelastic relaxation time τ(C). To explore the resulting interplay between active and polymeric dynamics, we first generalize a linear stability analysis (from earlier studies without polymer) to derive criteria for the onset of spontaneous heterogeneous flows (strain rate) and/or deformations (strain). We find two modes of instability. The first is a viscous mode, associated with strain rate perturbations. It dominates for relatively small values of τ(C) and is a simple generalization of the instability known previously without polymer. The second is an elastomeric mode, associated with strain perturbations, which dominates at large τ(C) and persists even as τ(C)→∞. We explore the dynamical states to which these instabilities lead by means of direct numerical simulations. These reveal oscillatory shear-banded states in one dimension and activity-driven turbulence in two dimensions even in the elastomeric limit τ(C)→∞. Adding polymer can also have calming effects, increasing the net throughput of spontaneous flow along a channel in a type of drag reduction. The effect of including strong antagonistic coupling between the nematic and polymer is examined numerically, revealing a rich array of spontaneously flowing states.

Particle-size effects in the formation of bicontinuous Pickering emulsions

We demonstrate that the formation of bicontinuous emulsions stabilized by interfacial particles (bijels) is more robust when nanoparticles rather than microparticles are used. Emulsification via spinodal demixing in the presence of nearly neutrally wetting particles is induced by rapid heating. Using confocal microscopy, we show that nanospheres allow successful bijel formation at heating rates two orders of magnitude slower than is possible with microspheres. In order to explain our results, we introduce the concept of mechanical leeway, i.e., nanoparticles benefit from a smaller driving force towards disruptive curvature. Finally, we suggest that leeway mechanisms may benefit any formulation in which challenges arise due to tight restrictions on a pivotal parameter, but where the restrictions can be relaxed by rationally changing the value of a more accessible parameter.

Active viscoelastic matter: from bacterial drag reduction to turbulent solids

A paradigm for internally driven matter is the active nematic liquid crystal, whereby the equations of a conventional nematic are supplemented by a minimal active stress that violates time-reversal symmetry. In practice, active fluids may have not only liquid-crystalline but also viscoelastic polymer degrees of freedom. Here we explore the resulting interplay by coupling an active nematic to a minimal model of polymer rheology. We find that adding a polymer can greatly increase the complexity of spontaneous flow, but can also have calming effects, thereby increasing the net throughput of spontaneous flow along a pipe (a "drag-reduction" effect). Remarkably, active turbulence can also arise after switching on activity in a sufficiently soft elastomeric solid.

A minimal physical model captures the shapes of crawling cells

Cell motility in higher organisms (eukaryotes) is crucial to biological functions ranging from wound healing to immune response, and also implicated in diseases such as cancer. For cells crawling on hard surfaces, significant insights into motility have been gained from experiments replicating such motion in vitro. Such experiments show that crawling uses a combination of actin treadmilling (polymerization), which pushes the front of a cell forward, and myosin-induced stress (contractility), which retracts the rear. Here we present a simplified physical model of a crawling cell, consisting of a droplet of active polar fluid with contractility throughout, but treadmilling connected to a thin layer near the supporting wall. The model shows a variety of shapes and/or motility regimes, some closely resembling cases seen experimentally. Our work strongly supports the view that cellular motility exploits autonomous physical mechanisms whose operation does not need continuous regulatory effort.

Filling an emulsion drop with motile bacteria

We have measured the spatial distribution of motile Escherichia coli inside spherical water droplets emulsified in oil. At low cell concentrations, the cell density peaks at the water-oil interface; at increasing concentration, the bulk of each droplet fills up uniformly while the surface peak remains. Simulations and theory show that the bulk density results from a "traffic" of cells leaving the surface layer, increasingly due to cell-cell scattering as the surface coverage rises above ∼10%. Our findings show similarities with the physics of a rarefied gas in a spherical cavity with attractive walls.

Spontaneous motility of passive emulsion droplets in polar active gels

We study by computer simulations the dynamics of a droplet of passive, isotropic fluid, embedded in a polar active gel. The latter represents a fluid of active force dipoles, which exert either contractile or extensile stresses on their surroundings, modelling for instance a suspension of cytoskeletal filaments and molecular motors. When the polarisation of the active gel is anchored normal to the droplet at its surface, the nematic elasticity of the active gel drives the formation of a hedgehog defect; this defect then drives an active flow which propels the droplet forward. In an extensile gel, motility can occur even with tangential anchoring, which is compatible with a defect-free polarisation pattern. In this case, upon increasing activity the droplet first rotates uniformly, and then undergoes a discontinuous nonequilibrium transition into a translationally motile state, powered by bending deformations in the surrounding active medium.

Discontinuous shear thickening without inertia in dense non-Brownian suspensions

A consensus is emerging that discontinuous shear thickening (DST) in dense suspensions marks a transition from a flow state where particles remain well separated by lubrication layers, to one dominated by frictional contacts. We show here that reasonable assumptions about contact proliferation predict two distinct types of DST in the absence of inertia. The first occurs at densities above the jamming point of frictional particles; here, the thickened state is completely jammed and (unless particles deform) cannot flow without inhomogeneity or fracture. The second regime shows strain-rate hysteresis and arises at somewhat lower densities, where the thickened phase flows smoothly. DST is predicted to arise when finite-range repulsions defer contact formation until a characteristic stress level is exceeded.

Chemotactic clusters in confined run-and-tumble bacteria: a numerical investigation

We present a simulation study of pattern formation in an ensemble of chemotactic run-and-tumble bacteria, focussing on the effect of spatial confinement, either within traps or inside a maze. These geometries are inspired by previous experiments probing pattern formation in chemotactic strains of E. coli under these conditions. Our main result is that a microscopic model of chemotactic run-and-tumble particles which themselves secrete a chemoattractant is able to reproduce the main experimental observations, namely the formation of bacterial aggregates within traps and in dead ends of a maze. Our simulations also demonstrate that stochasticity plays a key role and leads to a hysteretic response when the chemotactic sensitivity is varied. We compare the results of run-and-tumble particles with simulations performed with a simplified version of the model where the active particles are smooth swimmers which respond to chemotactic gradients by rotating towards the source of chemoattractant. This class of models leads again to aggregation, but with quantitative and qualitative differences in, for instance, the size and shape of clusters.

Intracellular facilitated diffusion: searchers, crowders, and blockers

In bacteria, regulatory proteins search for a specific DNA-binding target via "facilitated diffusion": a series of rounds of three-dimensional diffusion in the cytoplasm, and one-dimensional (1D) linear diffusion along the DNA contour. Using large scale Brownian dynamics simulations we find that each of these steps is affected differently by crowding proteins, which can either be bound to the DNA acting as a road block to the 1D diffusion, or freely diffusing in the cytoplasm. Macromolecular crowding can strongly affect mechanistic features such as the balance between three-dimensional and 1D diffusion, but leads to surprising robustness of the total search time.

Colloidal templating at a cholesteric-oil interface: assembly guided by an array of disclination lines

We simulate colloids (radius R ~ 1 μm) trapped at the interface between a cholesteric liquid crystal and an immiscible oil at which the helical order (pitch p) in the bulk conflicts with the orientation induced at the interface, stabilizing an ordered array of disclinations. For a weak anchoring strength W of the director field at the colloidal surface, this creates a template, favoring particle positions either on top of or midway between defect lines, depending on α=R/p. For small α, optical microscopy experiments confirm this picture, but for larger α no templating is seen. This may stem from the emergence at moderate W of a rugged energy landscape associated with defect reconnections.

Facilitated diffusion on mobile DNA: configurational traps and sequence heterogeneity

We present Brownian dynamics simulations of the facilitated diffusion of a protein, modeled as a sphere with a binding site on its surface, along DNA, modeled as a semiflexible polymer. We consider both the effect of DNA organization in three dimensions and of sequence heterogeneity. We find that in a network of DNA loops, which are thought to be present in bacterial DNA, the search process is very sensitive to the spatial location of the target within such loops. Therefore, specific genes might be repressed or promoted by changing the local topology of the genome. On the other hand, sequence heterogeneity creates traps which normally slow down facilitated diffusion. When suitably positioned, though, these traps can, surprisingly, render the search process much more efficient.

Bulk rheology and microrheology of active fluids

We simulate macroscopic shear experiments in active nematics and compare them with microrheology simulations where a spherical probe particle is dragged through an active fluid. In both cases we define an effective viscosity: in the case of bulk shear simulations this is the ratio between shear stress and shear rate, whereas in the microrheology case it involves the ratio between the friction coefficient and the particle size. We show that this effective viscosity, rather than being solely a property of the active fluid, is affected by the way chosen to measure it, and strongly depends on details such as the anchoring conditions at the probe surface and on both the system size and the size of the probe particle.

First-principles constitutive equation for suspension rheology

We provide a detailed derivation of a recently developed first-principles approach to calculating averages in systems of interacting, spherical Brownian particles under time-dependent flow. Although we restrict ourselves to flows which are both homogeneous and incompressible, the time dependence and geometry (e.g., shear and extension) are arbitrary. The approximations formulated within mode-coupling theory are particularly suited to dense colloidal suspensions and capture the slow relaxation arising from particle interactions and the resulting glass transition to an amorphous solid. The delicate interplay between slow structural relaxation and time-dependent external flow in colloidal suspensions thus may be studied within a fully tensorial theory.

Confined cubic blue phases under shear

We study the behaviour of confined cubic blue phases under shear flow via lattice Boltzmann simulations. We focus on the two experimentally observed phases, blue phase I and blue phase II. The disclination network of blue phase II continuously breaks and reforms under shear, leading to an oscillatory stress response in time. The oscillations are only regular for very thin samples. For thicker samples, the shear leads to a 'stick-slip' motion of part of the network along the vorticity direction. Blue phase I responds very differently: its defect network undergoes seemingly chaotic rearrangements under shear, irrespective of system size.

Colloids in active fluids: anomalous microrheology and negative drag

We simulate an experiment in which a colloidal probe is pulled through an active nematic fluid. We find that the drag on the particle is non-stokesian (not proportional to its radius). Strikingly, a large enough particle in contractile fluid (such as an actomyosin gel) can show negative viscous drag in steady state: the particle moves in the opposite direction to the externally applied force. We explain this, and the qualitative trends seen in our simulations, in terms of the disruption of orientational order around the probe particle and the resulting modifications to the active stress.

Size limit for particle-stabilized emulsion droplets under gravity

We demonstrate that emulsion droplets stabilized by interfacial particles become unstable beyond a size threshold set by gravity. This holds not only for colloids but also for supracolloidal glass beads, using which we directly observe the ejection of particles near the droplet base. The number of particles acting together in these ejection events decreases with time until a stable acornlike configuration is reached. Stability occurs when the weight of all remaining particles is less than the interfacial binding force of one particle. We also show the importance of the curvature of the droplet surface in promoting particle ejection.

Diffusive transport without detailed balance in motile bacteria: does microbiology need statistical physics?

Microbiology is the science of microbes, particularly bacteria. Many bacteria are motile: they are capable of self-propulsion. Among these, a significant class execute so-called run-and-tumble motion: they follow a fairly straight path for a certain distance, then abruptly change direction before repeating the process. This dynamics has something in common with Brownian motion (it is diffusive at large scales), and also something in contrast. Specifically, motility parameters such as the run speed and tumble rate depend on the local environment and hence can vary in space. When they do so, even if a steady state is reached, this is not generally invariant under time reversal: the principle of detailed balance, which restores the microscopic time-reversal symmetry of systems in thermal equilibrium, is mesoscopically absent in motile bacteria. This lack of detailed balance (allowed by the flux of chemical energy that drives motility) creates pitfalls for the unwary modeller. Here I review some statistical-mechanical models for bacterial motility, presenting them as a paradigm for exploring diffusion without detailed balance. I also discuss the extent to which statistical physics is useful in understanding real or potential microbiological experiments.

Cytoplasmic streaming in plant cells: the role of wall slip

We present a computer simulation study, via lattice Boltzmann simulations, of a microscopic model for cytoplasmic streaming in algal cells such as those of Chara corallina. We modelled myosin motors tracking along actin lanes as spheres undergoing directed motion along fixed lines. The sphere dimension takes into account the fact that motors drag vesicles or other organelles, and, unlike previous work, we model the boundary close to which the motors move as walls with a finite slip layer. By using realistic parameter values for actin lane and myosin density, as well as for endoplasmic and vacuole viscosity and the slip layer close to the wall, we find that this simplified view, which does not rely on any coupling between motors, cytoplasm and vacuole other than that provided by viscous Stokes flow, is enough to account for the observed magnitude of streaming velocities in intracellular fluid in living plant cells.

Simple model for the deformation-induced relaxation of glassy polymers

Glassy polymers show "strain hardening": at constant extensional load, their flow first accelerates, then arrests. Recent experiments have found this to be accompanied by a striking and unexplained dip in the segmental relaxation time. Here we explain such behavior by combining a minimal model of flow-induced liquefaction of a glass with a description of the stress carried by strained polymers, creating a nonfactorable interplay between aging and strain-induced rejuvenation. Under constant load, liquefaction of segmental motion permits strong flow that creates polymer-borne stress. This slows the deformation enough for the segmental modes to revitrify, causing strain hardening.

Self-assembly and nonlinear dynamics of dimeric colloidal rotors in cholesterics

We study by simulation the physics of two colloidal particles in a cholesteric liquid crystal with tangential order parameter alignment at the particle surface. The effective force between the pair is attractive at short range and favors assembly of colloid dimers at specific orientations relative to the local director field. When pulled through the fluid by a constant force along the helical axis, we find that such a dimer rotates, either continuously or stepwise with phase-slip events. These cases are separated by a sharp dynamical transition and lead, respectively, to a constant or an ever-increasing phase lag between the dimer orientation and the local nematic director.

Modelling thermal fluctuations in non-ideal fluids with the lattice Boltzmann method

Recently, we proposed a theoretical framework to include thermal fluctuations into the Lattice Boltzmann (LB) method for non-ideal fluids. Here, we apply a variant thereof to a certain class of force-based non-ideal fluid LB models. We find that ideal-gas-like noise is an exact result of the fluctuation-dissipation theorem in the hydrodynamic regime. It is shown that satisfactory equilibration of the density and fluid momentum can be obtained in a simulation over a wide range of length scales.

Crystallization mechanism of hard sphere glasses

In supercooled liquids, vitrification generally suppresses crystallization. Yet some glasses can still crystallize despite the arrest of diffusive motion. This ill-understood process may limit the stability of glasses, but its microscopic mechanism is not yet known. Here we present extensive computer simulations addressing the crystallization of monodisperse hard-sphere glasses at constant volume (as in a colloid experiment). Multiple crystalline patches appear without particles having to diffuse more than one diameter. As these patches grow, the mobility in neighboring areas is enhanced, creating dynamic heterogeneity with positive feedback. The future crystallization pattern cannot be predicted from the coordinates alone: Crystallization proceeds by a sequence of stochastic micronucleation events, correlated in space by emergent dynamic heterogeneity.

Crystallization and aging in hard-sphere glasses

We report new results from our programme of molecular dynamics simulation of hard-sphere systems, focusing on crystallization and glass formation at high concentrations. First we consider a much larger system than hitherto, N = 86 400 equal-sized particles. The results are similar to those obtained with a smaller system, studied previously, showing conventional nucleation and growth of crystals at concentrations near melting and crossing over to a spinodal-like regime at higher concentrations where the free energy barrier to nucleation appears to be negligible. Second, we investigate the dependence on the initial state of the system. We have devised a Monte Carlo 'constrained aging' method to move the particles in such a way that crystallization is discouraged. After a period of such aging, the standard molecular dynamics programme is run. For a system of N = 3200, we find that constrained aging encourages caging of the particles and slows crystallization somewhat. Nevertheless, both aged and unaged systems crystallize at volume fraction φ = 0.61 whereas neither system shows full crystallization in the duration of the simulation at φ = 0.62, a concentration still significantly below that of random close packing.

Nonlinear dynamics and rheology of active fluids: simulations in two dimensions

We report simulations of a continuum model for (apolar, flow aligning) active fluids in two dimensions. Both free and anchored boundary conditions are considered, at parallel confining walls that are either static or moving at fixed relative velocity. We focus on extensile materials and find that steady shear bands, previously shown to arise ubiquitously in one dimension for the active nematic phase at small (or indeed zero) shear rate, are generally replaced in two dimensions by more complex flow patterns that can be stationary, oscillatory, or apparently chaotic. The consequences of these flow patterns for time-averaged steady-state rheology are examined.

Structure of blue phase III of cholesteric liquid crystals

We report large scale simulations of the blue phases of cholesteric liquid crystals. Our results suggest a structure for blue phase III, the blue fog, which has been the subject of a long debate in liquid crystal physics. We propose that blue phase III is an amorphous network of disclination lines, which is thermodynamically and kinetically stabilized over crystalline blue phases at intermediate chiralities. This amorphous network becomes ordered under an applied electric field, as seen in experiments.

Shear banding and flow-concentration coupling in colloidal glasses

We report experiments on hard-sphere colloidal glasses that show a type of shear banding hitherto unobserved in soft glasses. We present a scenario that relates this to an instability due to shear-concentration coupling, a mechanism previously thought unimportant in these materials. Below a characteristic shear rate γ(c) we observe increasingly nonlinear and localized velocity profiles. We attribute this to very slight concentration gradients in the unstable flow regime. A simple model accounts for both the observed increase of γ(c) with concentration, and the fluctuations in the flow.

Thermal fluctuations in the lattice Boltzmann method for nonideal fluids

We introduce thermal fluctuations in the lattice Boltzmann method for nonideal fluids. A fluctuation-dissipation theorem is derived within the Langevin framework and applied to a specific lattice Boltzmann model that approximates the linearized fluctuating Navier-Stokes equations for fluids based on square-gradient free-energy functionals. The obtained thermal noise is shown to ensure equilibration of all degrees of freedom in a simulation to high accuracy. Furthermore, we demonstrate that satisfactory results for most practical applications of fluctuating hydrodynamics can already be achieved using thermal noise derived in the long-wavelength limit.

Colloids in cholesterics: size-dependent defects and non-stokesian microrheology

We simulate a colloidal particle (radius R) in a cholesteric liquid crystal (pitch p) with tangential order parameter alignment at the particle surface. The local defect structure evolves from a dipolar pair of surface defects (boojums) at small R/p to a pair of twisted disclination lines wrapping around the particle at larger values. On dragging the colloid with small velocity v through the medium along the cholesteric helix axis (an active microrheology measurement), we find a hydrodynamic drag force that scales linearly with v but superlinearly with R-in striking violation of Stokes' law, as generally used to interpret such measurements.

Run-and-tumble particles with hydrodynamics: sedimentation, trapping, and upstream swimming

We simulate by lattice Boltzmann the nonequilibrium steady states of run-and-tumble particles (inspired by a minimal model of bacteria), interacting by far-field hydrodynamics, subject to confinement. Under gravity, hydrodynamic interactions barely perturb the steady state found without them, but for particles in a harmonic trap such a state is quite changed if the run length is larger than the confinement length: a self-assembled pump is formed. Particles likewise confined in a narrow channel show a generic upstream flux in Poiseuille flow: chiral swimming is not required.

Thermodynamics of blue phases in electric fields

We present extensive numerical studies to determine the phase diagrams of cubic and hexagonal blue phases in an electric field. We confirm the earlier prediction that hexagonal phases, both two and three dimensional, are stabilized by a field, but we significantly refine the phase boundaries, which were previously estimated by means of a semianalytical approximation. In particular, our simulations show that the blue phase I-blue phase II transition at fixed chirality is largely unaffected by electric field, as observed experimentally.

Crystallization of hard-sphere glasses

We study by molecular dynamics the interplay between arrest and crystallization in hard spheres. For state points in the plane of volume fraction (0.54 <or= varphi <or= 0.63) and polydispersity (0 <or= s <or= 0.085), we delineate states that spontaneously crystallize from those that do not. For noncrystallizing (or precrystallization) samples we find isodiffusivity lines consistent with an ideal glass transition at varphi_{g} approximately 0.585, independent of s. Despite this, for s < 0.05, crystallization occurs at varphi > varphi_{g}. This happens on time scales for which the system is aging, and a diffusive regime in the mean square displacement is not reached; by those criteria, the system is a glass. Hence, contrary to a widespread assumption in the colloid literature, the occurrence of spontaneous crystallization within a bulk amorphous state does not prove that this state was an ergodic fluid rather than a glass.

First-principles constitutive equation for suspension rheology

Using mode-coupling theory, we derive a constitutive equation for the nonlinear rheology of dense colloidal suspensions under arbitrary time-dependent homogeneous flow. Generalizing previous results for simple shear, this allows the full tensorial structure of the theory to be identified. Macroscopic deformation measures, such as the Cauchy-Green tensors, thereby emerge. So does a direct relation between the stress and the distorted microstructure, illuminating the interplay of slow structural relaxation and arbitrary imposed flow. We present flow curves for steady planar and uniaxial elongation and compare these to simple shear. The resulting nonlinear Trouton ratios point to a tensorially nontrivial dynamic yield condition for colloidal glasses.

Shearing active gels close to the isotropic-nematic transition

We study numerically the rheological properties of a slab of active gel close to the isotropic-nematic transition. The flow behavior shows a strong dependence on the sample size, boundary conditions, and on the bulk constitutive curve, which, on entering the nematic phase, acquires an activity-induced discontinuity at the origin. The precursor of this within the metastable isotropic phase for contractile systems (e.g., actomyosin gels) gives a viscosity divergence; its counterpart for extensile suspensions admits instead a shear-banded flow with zero apparent viscosity.

Arrest of fluid demixing by nanoparticles: a computer simulation study

We use lattice Boltzmann simulations to investigate the formation of arrested structures upon demixing of a binary solvent containing neutrally wetting colloidal particles. Previous simulations for symmetric fluid quenches pointed to the formation of "bijels": bicontinuous interfacially jammed emulsion gels. These should be created when a glassy monolayer of particles forms at the fluid-fluid interface, arresting further demixing and rigidifying the structure. Experimental work has broadly confirmed this scenario, but it shows that bijels can also be formed in volumetrically asymmetric quenches. Here, we present new simulation results for such quenches, compare these to the symmetric case, and find a crossover to an arrested droplet phase at strong asymmetry. We then make extensive new analyses of the postarrest dynamics in our simulated bijel and droplet structures, on time scales comparable to the Brownian time for colloid motion. Our results suggest that, on these intermediate time scales, the effective activation barrier to ejection of particles from the fluid-fluid interface is smaller by at least 2 orders of magnitude than the corresponding barrier for an isolated particle on a flat interface.

Statistical mechanics of interacting run-and-tumble bacteria

We consider self-propelled particles undergoing run-and-tumble dynamics (as exhibited by E. coli) in one dimension. Building on previous analyses at drift-diffusion level for the one-particle density, we add both interactions and noise, enabling discussion of domain formation by "self-trapping," and other collective phenomena. Mapping onto detailed-balance systems is possible in certain cases.