BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

M. E. Cates

M. E. Cates appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechBiological Physicscond-mat.mtrl-sciphysics.flu-dyncond-mat.mes-hallphysics.comp-phCell BehaviorSubcellular Processescond-matcond-mat.dis-nncond-mat.othermath-phmath.MPnlin.CGOther Quantitative BiologyPopulations and EvolutionQuantitative Methods

Catalog footprint

What is connected

56works
18topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

56 published item(s)

preprint2016arXiv

Light-induced self-assembly of active rectification devices

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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 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 $τ_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 $τ_C$ and persists even as $τ_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 $τ_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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Active Brownian Particles and Run-and-Tumble Particles: a Comparative Study

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(ρ) < 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 $ρ$. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pressure and Phase Equilibria in Interacting Active Brownian Spheres

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 $ρ$. Our exact result relates $P$, in closed form, to bulk correlators and shows that (i) $P(ρ)$ 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(ρ)$ 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$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pressure is not a state function for generic active fluids

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Self-assembly of colloid-cholesteric composites provides a possible route to switchable optical materials

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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 $\sim 10\%$. Our findings show similarities with the physics of a rarefied gas in a spherical cavity with attractive walls.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

We simulate colloids (radius $R \sim 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 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 $α= 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Flexoelectric switching in cholesteric blue phases

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Colloids in active fluids: Anomalous micro-rheology 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Complex Fluids: The Physics of Emulsions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Facilitated diffusion on mobile DNA: configurational traps and sequence heterogeneity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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, 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modeling the Relaxation of Polymer Glasses under Shear and Elongational Loads

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phase separation and rotor self-assembly in active particle suspensions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Rheology of Lamellar Liquid Crystals in Two and Three Dimensions: A Simulation Study

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Schematic Mode Coupling Theory of Glass Rheology: Single and Double Step Strains

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spontaneous symmetry breaking in active droplets provides a generic route to motility

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

When are active Brownian particles and run-and-tumble particles equivalent? Consequences for motility-induced phase separation

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 $ρ$ 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(ρ)$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A 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 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bijels Containing Magnetic Particles: A Simulation Study

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Colloids in liquid crystals: a lattice Boltzmann study

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dilatancy in the flow and fracture of stretched colloidal suspensions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Langevin theory of fluctuations in the discrete Boltzmann equation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lattice Models of Nonequilibrium Bacterial Dynamics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Migration of chemotactic bacteria in soft agar: role of gel concentration

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the structure of blue phase III

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Arrested phase separation in reproducing bacteria: a generic route to pattern formation?

We present a generic mechanism by which reproducing microorganisms, with a diffusivity that depends on the local population density, can form stable patterns. It is known that a decrease of swimming speed with density can promote separation into bulk phases of two coexisting densities; this is opposed by the logistic law for birth and death which allows only a single uniform density to be stable. The result of this contest is an arrested nonequilibrium phase separation in which dense droplets or rings become separated by less dense regions, with a characteristic steady-state length scale. Cell division mainly occurs in the dilute regions and cell death in the dense ones, with a continuous flux between these sustained by the diffusivity gradient. We formulate a mathematical model of this in a case involving run-and-tumble bacteria, and make connections with a wider class of mechanisms for density-dependent motility. No chemotaxis is assumed in the model, yet it predicts the formation of patterns strikingly similar to those believed to result from chemotactic behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hard discs under steady shear: comparison of Brownian dynamics simulations and mode coupling theory

Brownian dynamics simulations of bidisperse hard discs moving in two dimensions in a given steady and homogeneous shear flow are presented close to and above the glasstransition density. The stationary structure functions and stresses of shear-melted glass are compared quantitatively to parameter-free numerical calculations of monodisperse hard discs using mode coupling theory within the integration through transients framework. Theory qualitatively explains the properties of the yielding glass but quantitatively overestimatesthe shear-driven stresses and structural anisotropies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Lattice Boltzmann simulations of liquid crystalline fluids: active gels and blue phases

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ordering dynamics of blue phases entails kinetic stabilization of amorphous networks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Shear banding and flow-concentration coupling in colloidal glasses

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γ_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γ_c$ with concentration, and the fluctuations observed in the flow.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Thermal fluctuations in the lattice Boltzmann method for non-ideal fluids

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Zero-range processes with saturated condensation: the steady state and dynamics

We study a class of zero-range processes in which the real-space condensation phenomenon does not occur and is replaced by a saturated condensation: that is, an extensive number of finite-size "condensates" in the steady state. We determine the conditions under which this occurs, and investigate the dynamics of relaxation to the steady state. We identify two stages: a rapid initial growth of condensates followed by a slow process of activated evaporation and condensation. We analyze these nonequilibrium dynamics with a combination of meanfield approximations, first-passage time calculations and a fluctuation-dissipation type approach.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Domain growth in cholesteric blue phases: hybrid lattice Boltzmann simulations

Here we review a hybrid lattice Boltzmann algorithm to solve the equations of motion of cholesteric liquid crystals. The method consists in coupling a lattice Boltzmann solver for the Navier-Stokes equation to a finite difference method to solve the dynamical equations governing the evolution of the liquid crystalline order parameter. We apply this method to study the growth of cholesteric blue phase domains, within a cholesteric phase. We focus on the growth of blue phase II and on a thin slab geometry in which the domain wall is flat. Our results show that, depending on the chirality, the growing blue phase is either BPII with no or few defects, or another structure with hexagonal ordering. We hope that our simulations will spur further experimental investigations on quenches in micron-size blue phase samples. The computational size that our hybrid lattice Boltzmann scheme can handle suggest that large scale simulations of new generation of blue phase liquid crystal device are within reach.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hard Spheres: Crystallization and Glass Formation

Motivated by old experiments on colloidal suspensions, we report molecular dynamics simulations of assemblies of hard spheres, addressing crystallization and glass formation. The simulations cover wide ranges of polydispersity s (standard deviation of the particle size distribution divided by its mean) and particle concentration. No crystallization is observed for s > 0.07. For 0.02 < s < 0.07, we find that increasing the polydispersity at a given concentration slows down crystal nucleation. The main effect here is that polydispersity reduces the supersaturation since it tends to stabilise the fluid but to destabilise the crystal. At a given polydispersity (< 0.07) we find three regimes of nucleation: standard nucleation and growth at concentrations in and slightly above the coexistence region; "spinodal nucleation", where the free energy barrier to nucleation appears to be negligible, at intermediate concentrations; and, at the highest concentrations, a new mechanism, still to be fully understood, which only requires small re-arrangement of the particle positions. The cross-over between the second and third regimes occurs at a concentration, around 58% by volume, where the colloid experiments show a marked change in the nature of the crystals formed and the particle dynamics indicate an "ideal" glass transition.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Competition between glass transition and liquid-gas separation in attracting colloids

We present simulation results addressing the phenomena of colloidal gelation induced by attractive interactions. The liquid-gas transition is prevented by the glass arrest at high enough attraction strength, resulting in a colloidal gel. The dynamics of the system is controlled by the glass, with little effect of the liquid-gas transition. When the system separates in a liquid and vapor phases, even if the denser phase enters the non-ergodic region, the vapor phase enables the structural relaxation of the system as a whole.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Instability and spatiotemporal rheochaos in a shear-thickening fluid model

We model a shear-thickening fluid that combines a tendency to form inhomogeneous, shear-banded flows with a slow relaxational dynamics for fluid microstructure. The interplay between these factors gives rich dynamics, with periodic regimes (oscillating bands, travelling bands, and more complex oscillations) and spatiotemporal rheochaos. These phenomena, arising from constitutive nonlinearity not inertia, can occur even when the steady-state flow curve is monotonic. Our model also shows rheochaos in a low-dimensional truncation where sharply defined shear bands cannot form.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Jamming, hysteresis and oscillation in scalar models for shear thickening

We investigate shear thickening and jamming within the framework of a family of spatially homogeneous, scalar rheological models. These are based on the `soft glassy rheology' model of Sollich et al. [Phys. Rev. Lett. 78, 2020 (1997)], but with an effective temperature x that is a decreasing function of either the global stress σor the local strain l. For appropiate x=x(σ), it is shown that the flow curves include a region of negative slope, around which the stress exhibits hysteresis under a cyclically varying imposed strain rate \gd. A subclass of these x(σ) have flow curves that touch the \gd=0 axis for a finite range of stresses; imposing a stress from this range {\em jams} the system, in the sense that the strain γcreeps only logarithmically with time t, γ(t)\sim\ln t. These same systems may produce a finite asymptotic yield stress under an imposed strain, in a manner that depends on the entire stress history of the sample, a phenomenon we refer to as history--dependent jamming. In contrast, when x=x(l) the flow curves are always monotonic, but we show that some x(l) generate an oscillatory strain response for a range of steady imposed stresses. Similar spontaneous oscillations are observed in a simplified model with fewer degrees of freedom. We discuss this result in relation to the temporal instabilities observed in rheological experiments and stick--slip behaviour found in other contexts, and comment on the possible relationship with `delay differential equations' that are known to produce oscillations and chaos.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Rheological instability in a simple shear thickening model

We study the strain response to steady imposed stress in a spatially homogeneous, scalar model for shear thickening, in which the local rate of yielding Γ(l) of mesoscopic `elastic elements' is not monotonic in the local strain l. Despite this, the macroscopic, steady-state flow curve (stress vs. strain rate) is monotonic. However, for a broad class of Γ(l), the response to steady stress is not in fact steady flow, but spontaneous oscillation. We discuss this finding in relation to other theoretical and experimental flow instabilities. Within the parameter ranges we studied, the model does not exhibit rheo-chaos.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Jamming and Stress Propagation in Particulate Matter

We present simple models of particulate materials whose mechanical integrity arises from a jamming process. We argue that such media are generically "fragile", that is, they are unable to support certain types of incremental loading without plastic rearrangement. In such models, fragility is naturally linked to the marginal stability of force chain networks (granular skeletons) within the material. Fragile matter exhibits novel mechanical responses that may be relevant to both jammed colloids and cohesionless assemblies of poured, rigid grains.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Stress Propagation in Sand

We describe a new continuum approach to the modelling of stress propagation in static granular media, focussing on the conical sandpile created from a point source. We argue that the stress continuity equations should be closed by means of scale-free, local constitutive relations between different components of the stress tensor, encoding the construction history of the pile: this history determines the organization of the grains, and thereby the local relationship between stresses. Our preferred model FPA (Fixed Principle Axes) assumes that the eigendirections (but not the eigenvalues) of the stress tensor are determined forever when a material element is first buried. Stresses propagate along a nested set of archlike structures within the medium; the results are in good quantitative agreement with published experimental data. The FPA model is one of a larger class, called OSL (Oriented Stress Linearity) models, in which the direction of the characteristics for stress propagation are fixed at burial. We speculate on the connection between these characteristics and the stress paths observed microscopically.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Role of Metastable States in Phase Ordering Dynamics

We show that the rate of separation of two phases of different densities (e.g. gas and solid) can be radically altered by the presence of a metastable intermediate phase (e.g. liquid). Within a Cahn-Hilliard theory we study the growth in one dimension of a solid droplet from a supersaturated gas. A moving interface between solid and gas phases (say) can, for sufficient (transient) supersaturation, unbind into two interfaces separated by a slab of metastable liquid phase. We investigate the criteria for unbinding, and show that it may strongly impede the growth of the solid phase.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Stress Propagation and Arching in Static Sandpiles

We present a new approach to the modelling of stress propagation in static granular media, focussing on the conical sandpile constructed from a point source. We view the medium as consisting of cohesionless hard particles held up by static frictional forces; these are subject to microscopic indeterminacy which corresponds macroscopically to the fact that the equations of stress continuity are incomplete -- no strain variable can be defined. We propose that in general the continuity equations should be closed by means of a constitutive relation (or relations) between different components of the (mesoscopically averaged) stress tensor. The primary constitutive relation relates radial and vertical shear and normal stresses (in two dimensions, this is all one needs). We argue that the constitutive relation(s) should be local, and should encode the construction history of the pile: this history determines the organization of the grains at a mesoscopic scale, and thereby the local relationship between stresses. To the accuracy of published experiments, the pattern of stresses beneath a pile shows a scaling between piles of different heights (RSF scaling) which severely limits the form the constitutive relation can take ...

0 citations0 reviewsNo code linkedOpen access