Search

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 non-equilibrium systems is much harder to describe because there is no free energy, but instead a dynamics-dependent quasi-potential 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 quasi-potential, and hence nucleation barrier, for liquid-vapor phase separation. Crucially 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.

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., Comment: 17 pages, 5 figures, Suppl. Movie 1 can be found at https://www2.ph.ed.ac.uk/~dmarendu/ManyColloidsMovie.html

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., Comment: includes supplemental text and link to movies; added note on higher order terms in active stress; to appear in Phys. Rev. Lett

At the surfaces of autophoretic colloids, slip velocities arise from local chemical gradients that are many-body functions of particle configuration and activity. For rapid chemical diffusion, coupled with slip-induced hydrodynamic interactions, we deduce the chemohydrodynamic forces and torques between colloids. Near a no-slip wall, the forces can be expressed as gradients of a non-equilibrium potential which, by tuning the type of activity, can be varied from repulsive to attractive. When this potential has a barrier, we find arrested phase separation with a mean cluster size set by competing chemical and hydrodynamic forces. These are controlled in turn by the monopolar and dipolar contributions to the active chemical surface fluxes., Comment: to appear in the J. Chem. Phys

We investigate numerically the behaviour 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, $\chi$, setting the balance between interfacial and elastic forces. When $\chi$ and the concentration of the isotropic component are both low, the blue phase disclination lattice templates a cubic array of fluid cylinders. For larger $\chi$, 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., Comment: 7 pages, 4 figures; contains Supplementary Material as Appendix. Supplementary Movies Movies can be found at https://www2.ph.ed.ac.uk/~dmarendu/BPMixtureSM.html

Colloidal heat engines extract power out of a fluctuating bath by manipulating a confined tracer. Considering a self-propelled tracer surrounded by a bath of passive colloids, we optimize the engine performances based on the maximum available power. Our approach relies on an adiabatic mean-field treatment of the bath particles which reduces the many-body description into an effective tracer dynamics. It leads us to reveal that, when operated at constant activity, an engine can only produce less maximum power than its passive counterpart. In contrast, the output power of an isothermal engine, operating with cyclic variations of the self-propulsion without any passive equivalent, exhibits an optimum in terms of confinement and activity. Direct numerical simulations of the microscopic dynamics support the validity of these results even beyond the mean-field regime, with potential relevance to the design of experimental engines., Comment: 7 pages, 2 figures

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 $\tau_C$. To explore the resulting interplay between active and polymeric dynamics, we first generalise 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 $\tau_C$ and is a simple generalisation of the instability known previously without polymer. The second is an elastomeric mode, associated with strain perturbations, which dominates at large $\tau_C$ and persists even as $\tau_C\to\infty$. We explore the novel dynamical states to which these instabilities lead by means of direct numerical simulations. These reveal oscillatory shear-banded states in 1D, and activity-driven turbulence in 2D even in the elastomeric limit $\tau_C\to\infty$. Adding polymer can also have calming effects, increasing the net throughput of spontaneous flow along a channel in a new type of "drag-reduction". Finally the effect of including strong, antagonistic coupling between nematic and polymer is examined numerically, revealing a rich array of spontaneously flowing states., Comment: 25 pages, 21 figures; v2 fixed rendering issue with fig 4, updated to published version

Self-propelled colloidal objects, such as motile bacteria or synthetic microswimmers, have microscopically irreversible individual dynamics - a feature they share with all living systems. The incoherent behaviour of individual swimmers can then be harnessed (or "rectified") by microfluidic devices that create systematic motions impossible in equilibrium. Examples include flow of rotor particles round a circuit, steady rotation of a gear wheel in a bacterial bath, and pumping of bacteria between chambers by "funnel gates". Here we present a computational proof-of-concept study, showing that such active rectification devices might be created directly from an unstructured "primordial soup" of motile particles, solely by using spatially modulated illumination to control their local propulsion speed. Alongside both microscopic irreversibility and speed modulation, our mechanism requires spatial symmetry breaking, such as a chevron light pattern, and strong interactions between particles, such as volume exclusion causing a collisional slow-down at high density. These four factors create a many-body rectification mechanism that generically differs from one-body microfluidic antecedents. Our work suggests that standard spatial-light-modulator technology might allow the programmable, light-induced self-assembly of active rectification devices from an unstructured particle bath.

Active Brownian particles (ABPs) and Run-and-Tumble particles (RTPs) both self-propel at fixed speed $v$ along a body-axis ${\bf u}$ that reorients either through slow angular diffusion (ABPs) or sudden complete randomisation (RTPs). We compare the physics of these two model systems both at microscopic and macroscopic scales. Using exact results for their steady-state distribution in the presence of external potentials, we show that they both admit the same effective equilibrium regime perturbatively that breaks down for stronger external potentials, in a model-dependent way. In the presence of collisional repulsions such particles slow down at high density: their propulsive effort is unchanged, but their average speed along ${\bf u}$ becomes $v(\rho) < v$. A fruitful avenue is then to construct a mean-field description in which particles are ghost-like and have no collisions, but swim at a variable speed $v$ that is an explicit function or functional of the density $\rho$. We give numerical evidence that the recently shown equivalence of the fluctuating hydrodynamics of ABPs and RTPs in this case, which we detail here, extends to microscopic models of ABPs and RTPs interacting with repulsive forces., Comment: 32 pages, 6 figures

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., Comment: 11 pages, 8 figures

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., Comment: 37 pages, 11 figures, 1 table

We derive from first principles the mechanical pressure $P$, defined as the force per unit area on a bounding wall, in a system of spherical, overdamped, active Brownian particles at density $\rho$. Our exact result relates $P$, in closed form, to bulk correlators and shows that (i) $P(\rho)$ is a state function, independent of the particle-wall interaction; (ii) interactions contribute two terms to $P$, one encoding the slow-down that drives motility-induced phase separation, and the other a direct contribution well known for passive systems; (iii) $P(\rho)$ is equal in coexisting phases. We discuss the consequences of these results for the motility-induced phase separation of active Brownian particles, and show that the densities at coexistence do not satisfy a Maxwell construction on $P$., Comment: 6p (main text) + 4p (SI)

Pressure is the mechanical force per unit area that a confined system exerts on its container. In thermal equilibrium, it depends only on bulk properties (density, temperature, etc.) through an equation of state. Here we show that in a wide class of active systems the pressure depends on the precise interactions between the active particles and the confining walls. In general, therefore, active fluids have no equation of state, their mechanical pressures exhibit anomalous properties that defy the familiar thermodynamic reasoning that holds in equilibrium. The pressure remains a function of state, however, in some specific and well-studied active models that tacitly restrict the character of the particle-wall and/or particle-particle interactions., Comment: 8 pages + 9 SI pages, Nature Physics (2015)

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 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., Comment: 11 pages, 4 figures (supplementary information included); v2 has updated PACS numbers, and a note added with one additional citation; v3 fixes a typo in Eq. 6 and a minor correction to top phase diagram in Fig. 1

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 $\sim 10\%$. Our findings show similarities with the physics of a rarefied gas in a spherical cavity with attractive walls., Comment: 5 pages, 4 figures, Supporting Information (5 pages, 5 figures)

Colloidal particles dispersed in liquid crystals can form new materials with tunable elastic and electro-optic properties. In a periodic `blue phase' host, particles should template into colloidal crystals with potential uses in photonics, metamaterials, and transformational optics. Here we show by computer simulation that colloid/cholesteric mixtures can give rise to regular crystals, glasses, percolating gels, isolated clusters, twisted rings and undulating colloidal ropes. This structure can be tuned via particle concentration, and by varying the surface interactions of the cholesteric host with both the particles and confining walls. Many of these new materials are metastable: two or more structures can arise under identical thermodynamic conditions. The observed structure depends not only on the formulation protocol, but also on the history of an applied electric field. This new class of soft materials should thus be relevant to design of switchable, multistable devices for optical technologies such as smart glass and e-paper., Comment: Manuscript with 3 figures plus supporting text and figures

We simulate colloids (radius $R \sim 1\mu$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 weak anchoring strength W of the director field at the colloidal surface, this creates a template, favoring particle positions eitheron top of or midway between defect lines, depending on $\alpha = R/p$. For small $\alpha$, optical microscopy experiments confirm this picture, but for larger $\alpha$ no templating is seen. This may stem from the emergence at moderate W of a rugged energy landscape associated with defect reconnections., Comment: 5 pages, 4 figures

In bacteria, regulatory proteins search for a specific DNA binding target via "facilitated diffusion": a series of rounds of 3D diffusion in the cytoplasm, and 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 3D and 1D diffusion, but leads to surprising robustness of the total search time.

We present computer simulations of the response of a flexoelectric blue phase network, either in bulk or under confinement, to an applied field. We find a transition in the bulk between the blue phase I disclination network and a parallel array of disclinations along the direction of the applied field. Upon switching off the field, the system is unable to reconstruct the original blue phase but gets stuck in a metastable phase. Blue phase II is comparatively much less affected by the field. In confined samples, the anchoring at the walls and the geometry of the device lead to the stabilisation of further structures, including field-aligned disclination loops, splayed nematic patterns, and yet more metastable states. Our results are relevant to the understanding of the switching dynamics for a class of new, "superstable", blue phases which are composed of bimesogenic liquid crystals, as these materials combine anomalously large flexoelectric coefficients, and low or near-zero dielectric anisotropy., Comment: 15 pages, 9 figures

We present Brownian dynamics simulations of the facilitated diffusion of a protein, modelled as a sphere with a binding site on its surface, along DNA, modelled as a semi-flexible polymer. We consider both the effect of DNA organisation in 3D, and of sequence heterogeneity. We find that in a network of DNA loops, as 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.

We explore a generic mechanism whereby a droplet of active matter acquires motility by the spontaneous breakdown of a discrete symmetry. The model we study offers a simple representation of a "cell extract" comprising, e.g., a droplet of actomyosin solution. (Such extracts are used experimentally to model the cytoskeleton.) Actomyosin is an active gel whose polarity describes the mean sense of alignment of actin fibres. In the absence of polymerization and depolymerization processes ('treadmilling'), the gel's dynamics arises solely from the contractile motion of myosin motors; this should be unchanged when polarity is inverted. Our results suggest that motility can arise in the absence of treadmilling, by spontaneous symmetry breaking (SSB) of polarity inversion symmetry. Adapting our model to wall-bound cells in two dimensions, we find that as wall friction is reduced, treadmilling-induced motility falls but SSB-mediated motility rises. The latter might therefore be crucial in three dimensions where frictional forces are likely to be modest. At a supra-cellular level, the same generic mechanism can impart motility to aggregates of non-motile but active bacteria; we show that SSB in this (extensile) case leads generically to rotational as well as translational motion., Comment: 13 pages, 8 figures, 1 table

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., Comment: 12 pages, 10 figures

These lectures start with the mean field theory for a symmetric binary fluid mixture, addressing interfacial tension, the stress tensor, and the equations of motion (Model H). We then consider the phase separation kinetics of such a mixture: coalescence, Ostwald ripening, its prevention by trapped species, coarsening of bicontinuous states, and the role of shear flow. The third topic addressed is the stabilization of emulsions by using surfactants to reduce or even eliminate the interfacial tension between phases; the physics of bending energy, which becomes relevant in the latter case, is then presented briefly. The final topic is the creation of long-lived metastable emulsions by adsorption of colloidal particles or nano-particles at the fluid-fluid interface; alongside spherical droplets, these methods can be used to create a range of unconventional structures with potentially interesting properties that are only now being explored., Comment: Lecture notes for les Houches 2012 Summer School on Soft Interfaces

Glassy polymers show strain hardening: at constant extensional load, their flow first accelerates, then arrests. Recent experiments under such loading have found this to be accompanied by a striking dip in the segmental relaxation time. This can be explained by a minimal nonfactorable model combining flow-induced melting of a glass with the buildup of stress carried by strained polymers. Within this model, liquefaction of segmental motion permits strong flow that creates polymer-borne stress, slowing the deformation enough for the segmental (or solvent) modes to then re-vitrify. Here we present new results for the corresponding behavior under step-stress shear loading, to which very similar physics applies. To explain the unloading behavior in the extensional case requires introduction of a crinkle factor describing a rapid loss of segmental ordering. We discuss in more detail here the physics of this, which we argue involves non-entropic contributions to the polymer stress, and which might lead to some important differences between shear and elongation. We also discuss some fundamental and possibly testable issues concerning the physical meaning of entropic elasticity in vitrified polymers. Finally we present new results for the startup of steady shear flow, addressing the possible role of transient shear banding.

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.

Active Brownian particles (ABPs, such as self-phoretic colloids) swim at fixed speed $v$ along a body-axis ${\bf u}$ that rotates by slow angular diffusion. Run-and-tumble particles (RTPs, such as motile bacteria) swim with constant $\u$ until a random tumble event suddenly decorrelates the orientation. We show that when the motility parameters depend on density $\rho$ but not on ${\bf u}$, the coarse-grained fluctuating hydrodynamics of interacting ABPs and RTPs can be mapped onto each other and are thus strictly equivalent. In both cases, a steeply enough decreasing $v(\rho)$ causes phase separation in dimensions $d=2,3$, even when no attractive forces act between the particles. This points to a generic role for motility-induced phase separation in active matter. However, we show that the ABP/RTP equivalence does not automatically extend to the more general case of $\u$-dependent motilities.

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, 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 may thus be studied within a fully tensorial theory.

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., Comment: 5 pages, 3 figures

Adding a non-adsorbing polymer to passive colloids induces an attraction between the particles via the `depletion' mechanism. High enough polymer concentrations lead to phase separation. We combine experiments, theory and simulations to demonstrate that using active colloids (such as motile bacteria) dramatically changes the physics of such mixtures. First, significantly stronger inter-particle attraction is needed to cause phase separation. Secondly, the finite size aggregates formed at lower inter-particle attraction show unidirectional rotation. These micro-rotors demonstrate the self assembly of functional structures using active particles. The angular speed of the rotating clusters scales approximately as the inverse of their size, which may be understood theoretically by assuming that the torques exerted by the outermost bacteria in a cluster add up randomly. Our simulations suggest that both the suppression of phase separation and the self assembly of rotors are generic features of aggregating swimmers, and should therefore occur in a variety of biological and synthetic active particle systems., Comment: Main text: 6 pages, 5 figures. Supplementary information: 5 pages, 4 figures. Supplementary movies available from httP://www.pnas.org/lookup/suppl/doi:10.1073/pnas.1116334109/-/DCSupplemental

Mode coupling theory (MCT) has had notable successes in addressing the rheology of hard-sphere colloidal glasses, and also soft colloidal glasses such as star-polymers. Here, we explore the properties of a recently developed MCT-based schematic constitutive equation under idealized experimental protocols involving single and then double step strains. We find strong deviations from expectations based on factorable, BKZ-type constitutive models. Specifically, a nonvanishing stress remains long after the application of two equal and opposite step strains; this residual stress is a signature of plastic deformation. We also discuss the distinction between hypothetically instantaneous step strains and fast ramps. These are not generally equivalent in our MCT approach, with the latter more likely to capture the physics of experimental `step' strains. The distinction points to the different role played by reversible anelastic, and irreversible plastic rearrangements.

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., Comment: 18 pages (incl. appendix), 10 figures, accepted in J R Soc Interface

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., Comment: 4 + 5 pages

We present large scale computer simulations of the nonlinear bulk rheology of lamellar phases (smectic liquid crystals) at moderate to large values of the shear rate (Peclet numbers 10-100), in both two and three dimensions. In two dimensions we find that modest shear rates align the system and stabilise an almost regular lamellar phase, but high shear rates induce the nucleation and proliferation of defects, which in steady state is balanced by the annihilation of defects of opposite sign. The critical shear rate at onset of this second regime is controlled by thermodynamic and kinetic parameters; we offer a scaling analysis that relates the critical shear rate to a critical "capillary number" involving those variables. Within the defect proliferation regime, the defects may be partially annealed by slowly decreasing the applied shear rate; this causes marked memory effects, and history-dependent rheology. Simulations in three dimensions show instead shear-induced ordering even at the highest shear rates studied here. This suggests that the critical shear rate shifts markedly upward on increasing dimensionality. This may in part reflect the reduced constraints on defect motion, allowing them to find and annihilate each other more easily. Residual edge defects in the 3D aligned state mostly point along the flow velocity, an orientation impossible in two dimensions., Comment: 18 pages, 12 figures

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 disinclination 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., Comment: 10 pages, 5 figures, to be published in J. Phys. Condens. Matter (2012)

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 stabilised over crystalline blue phases at intermediate chiralities}. This amorphous network becomes ordered under an applied electric field, as seen in experiments., Comment: 9 pages, 6 figures

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 non-factorable 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 re-vitrify, causing strain hardening., Comment: 9 pages, 5 figures

We propose a hybrid lattice Boltzmann algorithm to simulate the hydrodynamics of colloidal particles inside a liquid crystalline host. To validate our algorithm, we study the static and the microrheology of a colloid in a nematic, with tangential anchoring of the director field at the particle surface, and we confirm theories and experiments showing that the drag force in a nematic is markedly anisotropic. We then apply our method to consider the case of a colloid inside a cholesteric, and with normal anchoring at the surface. We show that by tuning the ratio between particle size and cholesteric pitch it is possible to control the defect configuration around the particle, and to stabilise novel figure-of-eight or highly twisted loops close to the colloid surface., Comment: 7 pages, 5 figures

We study the migration of chemotactic wild-type Escherichia coli populations in semisolid (soft) agar in the concentration range C = 0.15-0.5% (w/v). For C < 0.35%, expanding bacterial colonies display characteristic chemotactic rings. At C = 0.35%, however, bacteria migrate as broad circular bands rather than sharp rings. These are growth/diffusion waves arising because of suppression of chemotaxis by the agar and have not been previously reported experimentally to our knowledge. For C = 0.4-0.5%, expanding colonies do not span the depth of the agar and develop pronounced front instabilities. The migration front speed is weakly dependent on agar concentration at C < 0.25%, but decreases sharply above this value. We discuss these observations in terms of an extended Keller-Segel model for which we derived novel transport parameter expressions accounting for perturbations of the chemotactic response by collisions with the agar. The model makes it possible to fit the observed front speed decay in the range C = 0.15-0.35%, and its solutions qualitatively reproduce the observed transition from chemotactic to growth/diffusion bands. We discuss the implications of our results for the study of bacteria in porous media and for the design of improved bacteriological chemotaxis assays., Comment: 28 pages, 5 figures. Published online at http://www.sciencedirect.com/science/article/pii/S0006349511007211

Bicontinuous, interfacially jammed emulsion gels (bijels) represent a class of soft solid materials in which interpenetrating domains of two immiscible fluids are stabilized by an interfacial colloidal monolayer. Such structures can be formed by arrested spinodal decomposition from an initially single-phase colloidal suspension. Here we explore by lattice Boltzmann simulation the possible effects of using magnetic colloids in bijels. This may allow additional control over the structure, during or after formation, by application of a magnetic field or field gradient. These effects are modest for typical parameters based on the magnetic nanoparticles used in conventional ferrofluids, although significantly larger particles might be appropriate here. Field gradient effects, which are cumulative across a sample, could then allow a route for controlled breakdown of bijels as they do for particle-stabilized droplet emulsions.

Concentrated particulate suspensions, commonplace in the pharmaceutical, cosmetic and food industries, display intriguing rheology. In particular, the dramatic increase in viscosity with strain rate (shear thickening and jamming) which is often observed at high volume fractions, is of strong practical and fundamental importance. Yet manufacture of these products and their subsequent dispensing often involves flow geometries substantially different from that of simple shear flow experiments. Here we show that the elongation and breakage of a filament of a colloidal fluid under tensile loading is closely related to the jamming transition seen in its shear rheology. However, the modified flow geometry reveals important additional effects. Using a model system with nearly hard-core interactions, we provide evidence of surprisingly strong viscoelasticity in such a colloidal fluid under tension. With high speed photography we also directly observe dilatancy and granulation effects, which lead to fracture above a critical elongation rate., Comment: 22pages, 4 figures

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., Comment: 4 pages, 4 figures

The discrete Boltzmann equation for both the ideal and a non-ideal fluid is extended by adding Langevin noise terms in order to incorporate the effects of thermal fluctuations. After casting the fluctuating discrete Boltzmann equation in a form appropriate to the Onsager-Machlup theory of linear fluctuations, the statistical properties of the noise are determined by invoking a fluctuation-dissipation theorem at the kinetic level. By integrating the fluctuating discrete Boltzmann equation, the fluctuating lattice Boltzmann equation is obtained, which provides an efficient way to solve the equations of fluctuating hydrodynamics for ideal and non-ideal fluids. Application of the framework to a generic force-based non-ideal fluid model leads to ideal gas-type thermal noise. Simulation results indicate proper thermalization of all degrees of freedom., Comment: 24 pages, 2 figures, v2: added corrections and clarifications; as published in J. Stat. Mech

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 1D for the active nematic phase at small (or indeed zero) shear rate, are generally replaced in 2D 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., Comment: 16 pages, 15 eps figures, submitted to Phys. Rev. E

We study a model of self propelled particles exhibiting run and tumble dynamics on lattice. This non-Brownian diffusion is characterised by a random walk with a finite persistence length between changes of direction, and is inspired by the motion of bacteria such as E. coli. By defining a class of models with multiple species of particle and transmutation between species we can recreate such dynamics. These models admit exact analytical results whilst also forming a counterpart to previous continuum models of run and tumble dynamics. We solve the externally driven non-interacting and zero-range versions of the model exactly and utilise a field theoretic approach to derive the continuum fluctuating hydrodynamics for more general interactions. We make contact with prior approaches to run and tumble dynamics off lattice and determine the steady state and linear stability for a class of crowding interactions, where the jump rate decreases as density increases. In addition to its interest from the perspective of nonequilibrium statistical mechanics, this lattice model constitutes and efficient tool to simulate a class of interacting run and tumble models relevant to bacterial motion, so long as certain conditions (that we derive) are met., Comment: 33 pages, 12 figures

We report experiments on hard sphere colloidal glasses that reveal a type of shear banding hitherto unobserved in soft glasses. We present a scenario that relates this to an instability arising from shear-concentration coupling, a mechanism previously thought unimportant in this class of materials. Below a characteristic shear rate $\dot\gamma_c$ we observe increasingly non-linear velocity profiles and strongly localized flows. We attribute this trend to very slight concentration gradients (likely to evade direct detection) arising in the unstable flow regime. A simple model accounts for both the observed increase of $\dot\gamma_c$ with concentration, and the fluctuations observed in the flow., Comment: 4 pages, 4 figures, accepted for publication in Phys. Rev. Lett

Lattice Boltzmann simulations have become a method of choice to solve the hydrodynamic equations of motion of a number of complex fluids. Here we review some recent applications of lattice Boltzmann to study the hydrodynamics of liquid crystalline materials. In particular, we focus on the study of (a) the exotic blue phases of cholesteric liquid crystals, and (b) active gels - a model system for actin plus myosin solutions or bacterial suspensions. In both cases lattice Boltzmann studies have proved useful to provide new insights into these complex materials., Comment: 16 pages, 3 figures

We introduce thermal fluctuations in the lattice Boltzmann method for non-ideal 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., Comment: 15 pages, 5 figures

The cubic blue phases of liquid crystals are fascinating and technologically promising examples of hierarchically structured soft materials, comprising ordered networks of defect lines (disclinations) within a liquid crystalline matrix. We present the first large-scale simulations of their domain growth, starting from a blue phase nucleus within a supercooled isotropic or cholesteric background. The nucleated phase is thermodynamically stable; one expects its slow orderly growth, creating a bulk cubic. Instead, we find that the strong propensity to form disclinations drives the rapid disorderly growth of a metastable amorphous defect network. During this process the original nucleus is destroyed; re-emergence of the stable phase may therefore require a second nucleation step. Our findings suggest that blue phases exhibit hierarchical behavior in their ordering dynamics, to match that in their structure., Comment: 11 pages, 5 figures, 2 supplementary figures, 2 supplementary tables, accepted by PNAS

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.

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 2 and 3 dimensional, are stabilized by a field, but we significantly refine the phase boundaries, which were previously estimated by means of a semi-analytical 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., Comment: submitted to Physical Review E, 7 pages (excluding figures), 12 figures