Source author record

D. Marenduzzo

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

43works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

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

43 published item(s)

preprint2020arXiv

A nonequilibrium strategy for fast target search on the genome

Vital biological processes such as genome repair require fast and efficient binding of selected proteins to specific target sites on DNA. Here we propose an active target search mechanism based on "chromophoresis", the dynamics of DNA-binding proteins up or down gradients in the density of epigenetic marks, or colours (biochemical tags on the genome). We focus on a set of proteins that deposit marks from which they are repelled---a case which is only encountered away from thermodynamic equilibrium. For suitable ranges of kinetic parameter values, chromophoretic proteins can perform unidirectional motion and are optimally redistributed along the genome. Importantly, they can also locally unravel a region of the genome which is collapsed due to self-interactions and "dive" deep into its core, for a striking enhancement of the efficiency of target search on such an inaccessible substrate. We discuss the potential relevance of chromophoresis for the location of DNA lesions.

preprint2020arXiv

Electrostatic inactivation of RNA viruses at air-water and liquid-liquid interfaces

Understanding the interactions between viruses and surfaces or interfaces is important, as they provide the principles underpinning the cleaning and disinfection of contaminated surfaces. Yet, the physics of such interactions is currently poorly understood. For instance, there are longstanding experimental observations suggesting that the presence of air-water interfaces can generically inactivate and kill viruses, yet the mechanism underlying this phenomenon remains unknown. Here we use theory and simulations to show that electrostatics provides one such mechanism, and that this is very general. Thus, we predict that the free energy of an RNA virus should increase by several thousands of $k_BT$ as the virion breaches an air-water interface. We also show that the fate of a virus approaching a generic liquid-liquid interface depends strongly on the detailed balance between interfacial and electrostatic forces, which can be tuned, for instance, by choosing different media to contact a virus-laden respiratory droplet. We propose that these results can be used to design effective strategies for surface disinfection. Intriguingly, tunability requires electrostatic and interfacial forces to scale similarly with viral size, which naturally occurs when charges are arranged in a double-shell distribution as in RNA viruses like influenza and all coronaviruses.

preprint2016arXiv

A 1-dimensional statistical mechanics model for nucleosome positioning on genomic DNA

The first level of folding of DNA in eukaryotes is provided by the so called '10 nm chromatin fibre', where DNA wraps around histone proteins (approx. 10 nm in size) to form nucleosomes, which go on to create a zig zagging bead on a string structure. In this work we present a 1 dimensional statistical mechanics model to study nucleosome positioning within one such 10 nm fibre. We focus on the case of genomic sheep DNA, and we start from effective potentials valid at infinite dilution and determined from high resolution in vitro salt dialysis experiments. We study positioning within a polynucleosome chain, and compare the results for genomic DNA to that obtained in the simplest case of homogeneous DNA, where the problem can be mapped to a Tonks gas. First, we consider the simple, analytically solvable, case where nucleosomes are assumed to be point like. Then, we perform numerical simulations to gauge the effect of their finite size on the nucleosomal distribution probabilities. Finally we compare nucleosome distributions and simulated nuclease digestion patterns for the two cases (homogeneous and sheep DNA), thereby providing testable predictions of the effect of sequence on experimentally observable quantities in experiments on polynucleosome chromatin fibres reconstituted in vitro.

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.

preprint2016arXiv

Shear dynamics of an inverted nematic emulsion

Here we study theoretically the dynamics of a 2D and a 3D isotropic droplet in a nematic liquid crystal under a shear flow. We find a large repertoire of possible nonequilibrium steady states as a function of the shear rate and of the anchoring of the nematic director field at the droplet surface. We first discuss homeotropic anchoring. For weak anchoring, we recover the typical behaviour of a sheared isotropic droplet in a binary fluid, which rotates, stretches and can be broken by the applied flow. For intermediate anchoring, new possibilities arise due to elastic effects in the nematic fluid. We find that in this regime the 2D droplet can tilt and move in the flow, or tumble incessantly at the centre of the channel. For sufficiently strong anchoring, finally, one or both of the topological defects which form close to the surface of the isotropic droplet in equilibrium detach from it and get dragged deep into the nematic state by the flow. In 3D, instead, the Saturn ring associated with normal anchoring disclination line can be deformed and shifted downstream by the flow, but remains always localized in proximity of the droplet, at least for the parameter range we explored. Tangential anchoring in 2D leads to a different dynamic response, as the boojum defects characteristic of this situation can unbind from the droplet under a weaker shear with respect to the normal anchoring case. Our results should stimulate further experiments with inverted liquid crystal emulsions under shear, as most of the predictions can be testable in principle by monitoring the evolution of liquid crystalline orientation patterns or by tracking the position and shape of the droplet over time.

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.

preprint2015arXiv

Binding of bivalent transcription factors to active and inactive regions folds human chromosomes into loops, rosettes and domains

Biophysicists are modeling conformations of interphase chromosomes, often basing the strengths of interactions between segments distant on the genetic map on contact frequencies determined experimentally. Here, instead, we develop a fitting-free, minimal model: bivalent red and green "transcription factors" bind to cognate sites in runs of beads ("chromatin") to form molecular bridges stabilizing loops. In the absence of additional explicit forces, molecular dynamic simulations reveal that bound "factors' spontaneously cluster -- red with red, green with green, but rarely red with green -- to give structures reminiscent of transcription factories. Binding of just two transcription factors (or proteins) to active and inactive regions of human chromosomes yields rosettes, topological domains, and contact maps much like those seen experimentally. This emergent "bridging-induced attraction" proves to be a robust, simple, and generic force able to organize interphase chromosomes at all scales.

preprint2015arXiv

Field-controlled columnar and planar patterning of cholesteric colloids

We study how dispersions of colloidal particles in a cholesteric liquid crystal behave under a time-dependent electric field. By controlling the amplitude and shape of the applied field wave, we show that the system can be reproducibly driven out of equilibrium through different kinetic pathways and navigated through a glassy-like free energy landscape encompassing many competing metastable equilibria. Such states range from simple Saturn rings to complex structures featuring amorphous defect networks, or stacks of disclination loops. A non-equilibrium electric field can also trigger the alignment of particles into columnar arrays, through defect-mediated force impulses, or their repositioning within a plane. Our results are promising in terms of providing new avenues towards controlled patterning and self-assembly of soft colloid-liquid crystal composite materials.

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.

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.

preprint2014arXiv

Models for twistable elastic polymers in Brownian dynamics, and their implementation for LAMMPS

An elastic rod model for semi-flexible polymers is presented. Theory for a continuum rod is reviewed, and it is shown that a popular discretised model used in numerical simulations gives the correct continuum limit. Correlation functions relating to both bending and twisting of the rod are derived for both continuous and discrete cases, and results are compared with numerical simulations. Finally, two possible implementations of the discretised model in the multi-purpose molecular dynamics software package LAMMPS are described.

preprint2014arXiv

Threading Dynamics of Ring Polymers in a Gel

We perform large scale three-dimensional molecular dynamics simulations of unlinked and unknotted ring polymers diffusing through a background gel, here a three-dimensional cubic lattice. Taking advantage of this architecture, we propose a new method to unambiguously identify and quantify inter-ring threadings (penetrations) and to relate these to the dynamics of the ring polymers. We find that both the number and the persistence time of the threadings increase with the length of the chains, ultimately leading to a percolating network of inter-ring penetrations. We discuss the implications of these findings for the possible emergence of a topological jammed state of very long rings.

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.

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.

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.

preprint2013arXiv

Mechanically driven growth of quasi-two dimensional microbial colonies

We study colonies of non-motile, rod-shaped bacteria growing on solid substrates. In our model, bacteria interact purely mechanically, by pushing each other away as they grow, and consume a diffusing nutrient. We show that mechanical interactions control the velocity and shape of the advancing front, which leads to features that cannot be captured by established Fisher-Kolmogorov models. In particular, we find that the velocity depends on the elastic modulus of bacteria or their stickiness to the surface. Interestingly, we predict that the radius of an incompressible, strictly two-dimensional colony cannot grow linearly in time. Importantly, mechanical interactions can also account for the nonequilibrium transition between circular and branching colonies, often observed in the lab.

preprint2013arXiv

Spiral and never-settling patterns in active suspensions

We present a combined numerical and analytical study of pattern formation in an active system where particles align, possess a density-dependent motility, and are subject to a logistic reaction. This is a model for suspensions of reproducing bacteria, but it can also represent, in the ordered phase, actomyosin gels in vitro or in vivo. In the disordered phase, we find that motility suppression and growth compete to yield stable or blinking patterns, which, when dense enough, acquire internal orientational ordering, to yield asters or spirals. In the ordered phase, the reaction term leads to previously unobserved never-settling patterns which can provide a simple framework to understand the formation of motile and spiral patterns in actin.

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.

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.

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.

preprint2012arXiv

Facilitated diffusion on confined DNA

In living cells, proteins combine 3D bulk diffusion and 1D sliding along the DNA to reach a target faster. This process is known as facilitated diffusion, and we investigate its dynamics in the physiologically relevant case of confined DNA. The confining geometry and DNA elasticity are key parameters: we find that facilitated diffusion is most efficient inside an isotropic volume, and on a flexible polymer. By considering the typical copy numbers of proteins in vivo, we show that the speedup due to sliding becomes insensitive to fine tuning of parameters, rendering facilitated diffusion a robust mechanism to speed up intracellular diffusion-limited reactions. The parameter range we focus on is relevant for in vitro systems and for facilitated diffusion on yeast chromatin.

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.

preprint2012arXiv

Pattern formation in self-propelled particles with density-dependent motility

We study the behaviour of interacting self-propelled particles, whose self-propulsion speed decreases with their local density. By combining direct simulations of the microscopic model with an analysis of the hydrodynamic equations obtained by explicitly coarse graining the model, we show that interactions lead generically to the formation of a host of patterns, including moving clumps, active lanes and asters. This general mechanism could explain many of the patterns seen in recent experiments and simulations.

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.

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.

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.

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.

preprint2011arXiv

Bistable defect structures in blue phase devices

Blue phases (BPs) are liquid crystals made up by networks of defects, or disclination lines. While existing phase diagrams show a striking variety of competing metastable topologies for these networks, very little is known as to how to kinetically reach a target structure, or how to switch from one to the other, which is of paramount importance for devices. We theoretically identify two confined blue phase I systems in which by applying an appropriate series of electric field it is possible to select one of two bistable defect patterns. Our results may be used to realise new generation and fast switching energy-saving bistable devices in ultrathin surface treated BPI wafers.

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.

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.

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.

preprint2011arXiv

Switching and defect dynamics in multistable liquid crystal devices

We investigate the switching dynamics of multistable nematic liquid crystal devices. In particular we identify a remarkably simple 2-dimensional (2D) device which exploits hybrid alignment at the surfaces to yield a bistable response. We also consider a 3-dimensional (3D) tristable nematic device with patterned anchoring, recently implemented in practice, and discuss how the director and disclination patterns change during switching.

preprint2011arXiv

Switching dynamics in cholesteric blue phases

Blue phases are networks of disclination lines, which occur in cholesteric liquid crystals near the transition to the isotropic phase. They have recently been used for the new generation of fast switching liquid crystal displays. Here we study numerically the steady states and switching hydrodynamics of blue phase I (BPI) and blue phase II (BPII) cells subjected to an electric field. When the field is on, there are three regimes: for very weak fields (and strong anchoring at the boundaries) the blue phases are almost unaffected, for intermediate fields the disclinations twist (for BPI) and unzip (for BPII), whereas for very large voltages the network dissolves in the bulk of the cell. Interestingly, we find that a BPII cell can recover its original structure when the field is switched off, whereas a BPI cell is found to be trapped more easily into metastable configurations. The kinetic pathways followed during switching on and off entails dramatic reorganisation of the disclination networks. We also discuss the effect of changing the director field anchoring at the boundary planes and of varying the direction of the applied field.

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.

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.

preprint2010arXiv

Non-equilibrium dynamics of an active colloidal "chucker"

We report Monte Carlo simulations of the dynamics of a "chucker": a colloidal particle which emits smaller solute particles from its surface, isotropically and at a constant rate k_c. We find that the diffusion constant of the chucker increases for small k_c, as recently predicted theoretically. At large k_c the chucker diffuses more slowly due to crowding effects. We compare our simulation results to those of a "point particle" Langevin dynamics scheme in which the solute concentration field is calculated analytically, and in which hydrodynamic effects can be included albeit in an approximate way. By simulating the dragging of a chucker, we obtain an estimate of its apparent mobility coefficient which violates the fluctuation-dissipation theorem. We also characterise the probability density profile for a chucker which sediments onto a surface which either repels or absorbs the solute particles, and find that the steady state distributions are very different in the two cases. Our simulations are inspired by the biological example of exopolysaccharide-producing bacteria, as well as by recent experimental, simulation and theoretical work on phoretic colloidal "swimmers".

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.

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.

preprint2010arXiv

The role of noise and advection in absorbing state phase transitions

We study the effect of advection and noise on the field theory for directed percolation (DP). We show that even a very small advective velocity is enough to change the universality class of the dynamic phase transition. When the noise is taken to be proportional to the square root of the population density, we find an additional nonequilibrium 'spinodal' line separating a region where an exponentially decreasing density is metastable, from another one in which it is unstable. If the noise is instead linear in the density, the phase diagram changes dramatically both quantitatively and qualitatively, and the spinodal line becomes a true phase boundary. We briefly discuss possible applications of our results to microbial sedimentation and population dynamics in rivers.

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.

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.

preprint2005arXiv

Continuum model for polymers with finite thickness

We consider the continuum limit of a recently-introduced model for discretized thick polymers, or tubes. We address both analytically and numerically how the polymer thickness influences the decay of tangent-tangent correlations and find how the persistence length scales with the thickness and the torsional rigidity of the tube centerline. At variance with the worm-like chain model, the phase diagram that we obtain for a continuous tube is richer; in particular, for a given polymer thickness there exists a threshold value for the centerline torsional rigidity separating a simple exponential decay of the tangent-tangent correlation from an oscillatory one.

preprint2005arXiv

Physics of thick polymers

We present the results of analytic calculations and numerical simulations of the behaviour of a new class of chain molecules which we call thick polymers. The concept of the thickness of such a polymer, viewed as a tube, is encapsulated by a special three body interaction and impacts on the behaviour both locally and non-locally. When thick polymers undergo compaction due to an attractive self-interaction, we find a new type of phase transition between a compact phase and a swollen phase at zero temperature on increasing the thickness. In the vicinity of this transition, short tubes form space filling helices and sheets as observed in protein native state structures. Upon increasing the chain length, or the number of chains, we numerically find a crossover from secondary structure motifs to a quite distinct class of structures akin to the semi-crystalline phase of polymers or amyloid fibers in polypeptides.