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

Davide Marenduzzo

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

ResearcherUnclaimed source record
cond-mat.softBiological Physicscond-mat.stat-mechphysics.flu-dynSubcellular ProcessesGenomicsnlin.PSphysics.comp-phBiomoleculesCell Behaviorcond-mat.mes-hallcond-mat.mtrl-sciphysics.pop-phQuantitative Methods

Catalog footprint

What is connected

41works
14topics
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

41 published item(s)

preprint2022arXiv

Active turbulence and spontaneous phase separation in inhomogeneous extensile active gels

We report numerical results for the hydrodynamics of inhomogeneous lyotropic and extensile active nematic gels. By simulating the coupled Cahn-Hilliard, Navier-Stokes, and Beris-Edwards equation for the evolution of the composition, flow and orientational order of an active nematic, we ask whether composition variations are important to determine its emergent physics. As in active gels of uniform composition, we find that increasing either activity or nematic tendency (e.g., overall active matter concentration) triggers a transition between an isotropic passive phase and an active nematic one. We show that composition inhomogeneities are important in the latter phase, where we find three types of possible dynamical regimes. First, we observe regular patterns with defects and vortices: these exist close to the passive-active transition. Second, for larger activity, or deeper in the nematic phase, we find active turbulence, as in active gels of uniform composition, but with exceedingly large composition variation. In the third regime, which is uniquely associated with inhomogeneity and occurs for large nematic tendency and low activity, we observe spontaneous microphase separation into active and passive domains. The microphase separated regime is notable in view of the absence of an explicit demixing term in the underlying free energy which we use, and we provide a theoretical analysis based on the common tangent construction which explains its existence. We hope this regime can be probed experimentally in the future.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A renormalization group study of the dynamics of active membranes: universality classes and scaling laws

Motivated by experimental observations of patterning at the leading edge of motile eukaryotic cells, we introduce a general model for the dynamics of nearly-flat fluid membranes driven from within by an ensemble of activators. We include, in particular, a kinematic coupling between activator density and membrane slope which generically arises whenever the membrane has a non-vanishing normal speed. We unveil the phase diagram of the model by means of a perturbative field-theoretical renormalization group analysis. Due to the aforementioned kinematic coupling the natural dynamical scaling is acoustic, that is the dynamical critical exponent is 1. However, as soon as the the normal velocity of the membrane is tuned to zero, the system crosses over to diffusive dynamic scaling in mean field. Distinct critical points can be reached depending on how the limit of vanishing velocity is realised: in each of them corrections to scaling due to nonlinear coupling terms must be taken into accounts. The detailed analysis of these critical points reveals novel scaling regimes wich can be accessed with perturbative methods, together with signs of strong coupling behaviour, which establishes a promising ground for further non-perturbative calculations. Our results unify several previous studies on the dynamics of active membrane, while also identifying nontrivial scaling regimes which cannot be captured by passive theories of fluctuating interfaces and are relevant for the physics of living membranes.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Cholesteric shells: two-dimensional blue fog and finite quasicrystals

We study the phase behaviour of a quasi-two dimensional cholesteric liquid crystal shell. We characterise the topological phases arising close to the isotropic-cholesteric transition, and show that they differ in a fundamental way from those observed on a flat geometry. For spherical shells, we discover two types of quasi-two dimensional topological phases: finite quasicrystals and amorphous structures, both made up by mixtures of polygonal tessellations of half-skyrmions. These structures generically emerge instead of regular double twist lattices because of geometric frustration, which disallows a regular hexagonal tiling of curved space. For toroidal shells, the variations in the local curvature of the surface stabilises heterogeneous phases where cholesteric patterns coexist with hexagonal lattices of half-skyrmions. Quasicrystals, amorphous and heterogeneous structures could be sought experimentally by self assembling cholesteric shells on the surface of emulsion droplets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Chaotic and periodical dynamics of active chiral droplets

The interplay between the chirality of many biological molecules and the energy injected at small length-scales as the result of biological processes is at the base of the life of the cells. With the aim of unveiling the connection between these two features, here we analyze by means of lattice Boltzmann simulations the behavior of an active droplet of cholesteric liquid crystal under the effect of intense active doping, within the framework of active gel theory. We find that a droplet of chiral liquid crystal, fueled by active force dipoles, develops defect loops (closed disclination lines) that pierce the interior of the droplet, leading the droplet to develop an erratic motility mode. When the droplet is fueled by in-warding active torque dipoles, three different dynamical regimes develops at varying both the thermodynamic chirality and the strength of active energy injection: a stable rotational state at low activity, an intermittent disclination dance regime, and a turbulent state where closed disclination lines formation is inhibited and new pairs of oppositely charged surface defects leads to the development of chaotic rotational motion. Finally, we show that out-warding torque dipoles are able to sustain a periodical dynamics at higher chirality characterized by the nucleation/annihilation of pairs of disclination rings.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Competition between local erasure and long-range spreading of a single biochemical mark leads to epigenetic bistability

The mechanism through which cells determine their fate is intimately related to the spreading of certain biochemical (so-called epigenetic) marks along their genome. The mechanisms behind mark spreading and maintenance are not yet fully understood, and current models often assume a long-range infection-like process for the dynamics of marks, due to the polymeric nature of the chromatin fibre which allows looping between distant sites. While these existing models typically consider antagonising marks, here we propose a qualitatively different scenario which analyses the spreading of a single mark. We define a 1D stochastic model in which mark spreading/infection occurs as a long-range process whereas mark erasure/recovery is a local process, with an enhanced rate at boundaries of infected domains. In the limiting case where our model exhibits absorbing states, we find a first-order-like transition separating the marked/infected phase from the unmarked/recovered phase. This suggests that our model, in this limit, belongs to the long-range compact directed percolation universality class. The abrupt nature of the transition is retained in a more biophysically realistic situation when a basal infection/recovery rate is introduced (thereby removing absorbing states). Close to the transition there is a range of bistability where both the marked/infected and unmarked/recovered states are metastable and long lived, which provides a possible avenue for controlling fate decisions in cells. Increasing the basal infection/recovery rate, we find a second transition between a coherent (marked or unmarked) phase, and a mixed, or random, one.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Interaction between Nearly Hard Colloidal Spheres at an Oil-Water Interface

We show that the interaction potential between sterically stabilized, nearly hard-sphere [poly(methylmethacrylate)-poly(lauryl methacrylate) (PMMA-PLMA)] colloids at a water-oil interface has a negligible unscreened-dipole contribution, suggesting that models previously developed for charged particles at liquid interfaces are not necessarily applicable to sterically stabilized particles. Interparticle potentials, $U(r)$, are extracted from radial distribution functions [$g(r)$, measured by fluorescence microscopy] via Ornstein-Zernike inversion and via a reverse Monte Carlo scheme. The results are then validated by particle tracking in a blinking optical trap. Using a Bayesian model comparison, we find that our PMMA-PLMA data is better described by a screened monopole only rather than a functional form having a screened monopole plus an unscreened dipole term. We postulate that the long range repulsion we observe arises mainly through interactions between neutral holes on a charged interface, i.e., the charge of the liquid interface cannot, in general, be ignored. In agreement with this interpretation, we find that the interaction can be tuned by varying salt concentration in the aqueous phase. Inspired by recent theoretical work on point charges at dielectric interfaces, which we explain is relevant here, we show that a screened $\frac{1}{r^2}$ term can also be used to fit our data. Finally, we present measurements for poly(methyl methacrylate)-poly(12-hydroxystearic acid) (PMMA-PHSA) particles at a water-oil interface. These suggest that, for PMMA-PHSA particles, there is an additional contribution to the interaction potential. This is in line with our optical-tweezer measurements for PMMA-PHSA colloids in bulk oil, which indicate that they are slightly charged.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pattern formation in active model C with anchoring: bands, aster networks, and foams

We study the dynamics of pattern formation in a minimal model for active mixtures made of microtubules and molecular motors. We monitor the evolution of the (conserved) microtubule density and of the (non-conserved) nematic order parameter, focusing on the effects of an "anchoring" term that provides a direct coupling between the preferred microtubule direction and their density gradient. The key control parameter is the ratio between activity and elasticity. When elasticity dominates, the interplay between activity and anchoring leads to formation of banded structures that can undergo additional bending or rotational instabilities. When activity dominates, the nature of anchoring instead gives rise to a range of active cellular solids, including aster-like networks, disordered foams and spindle-like patterns. We speculate that the introduced "active model C" with anchoring is a minimal model to describe pattern formation in a biomimetic analogue of the microtubule cytoskeleton.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rotation and propulsion in 3d active chiral droplets

Chirality is a recurrent theme in the study of biological systems, in which active processes are driven by the internal conversion of chemical energy into work. Bacterial flagella, acto-myosin filaments and microtubule bundles are active systems which are also intrinsically chiral. Despite some exploratory attempt to capture the relations between chirality and motility, no intrinsically chiral system has ever been analyzed so far. To address this gap in knowledge, here we study the effects of internal active forces and torques on a three-dimensional droplet of cholesteric liquid crystal (CLC) embedded in an isotropic liquid. We consider tangential anchoring of the liquid crystal director at the droplet surface. Contrary to what happens in nematics, where moderate extensile activity leads to droplet rotation, cholesteric active droplets exhibit a lot more complex and variegated behaviors. We find that extensile force dipole activity stabilises complex defect configurations whose orbiting dynamics couples to thermodynamic chirality to propel screw-like droplet motion. Instead, dipolar torque activity may either tighten or unwind the cholesteric helix, and if tuned can power rotations with an oscillatory angular velocity of zero mean.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Soft matter science and the COVID-19 pandemic

Much of the science underpinning the global response to the COVID-19 pandemic lies in the soft matter domain. Coronaviruses are composite particles with a core of nucleic acids complexed to proteins surrounded by a protein-studded lipid bilayer shell. A dominant route for transmission is via air-borne aerosols and droplets. Viral interaction with polymeric body fluids, particularly mucus, and cell membranes control their infectivity, while their interaction with skin and artificial surfaces underpins cleaning and disinfection and the efficacy of masks and other personal protective equipment. The global response to COVID-19 has highlighted gaps in the soft matter knowledge base. We survey these gaps, especially as pertaining to the transmission of the disease, and suggest questions that can (and need to) be tackled, both in response to COVID-19 and to better prepare for future viral pandemics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Solid-Liquid Transition of Deformable and Overlapping Active Particles

Experiments and theory have shown that cell monolayers and epithelial tissues exhibit solid-liquid and glass-liquid transitions. These transitions are biologically relevant to our understanding of embryonic development, wound healing, and cancer. Current models typically consider purely two-dimensional monolayers with no overlaps between neighboring cells. In reality, overlaps are important, and they may be precursors of cell extrusion -- a key biophysical process to maintain homeostasis in epithelial tissues. Here, we use a multi-phase field model to study the solid-liquid transition in a confluent monolayer of deformable cells which can overlap. When cells overlap rather than deform, we find that the melting transition changes from continuous to discontinuous, and that there is an intermittent regime close to the transition, where solid and liquid states alternate over time. By studying the dynamics of $5$- and $7$-fold disclinations in the hexagonal lattice formed by the cell centers, we observe that these correlate with spatial fluctuations in the cellular overlap, and that cell extrusion tends to initiate near $5$-fold disclinations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Swimming suppresses correlations in dilute suspensions of pusher microorganisms

Active matter exhibits various forms of non-equilibrium states in the absence of external forcing, including macroscopic steady-state currents. Such states are often too complex to be modelled from first principles and our understanding of their physics relies heavily on minimal models. These have mostly been studied in the case of "dry" active matter, where particle dynamics are dominated by friction with their surroundings. Significantly less is known about systems with long-range hydrodynamic interactions that belong to "wet" active matter. Dilute suspensions of motile bacteria, modelled as self-propelled dipolar particles interacting solely through long-ranged hydrodynamic fields, are arguably the most studied example from this class of active systems. Their phenomenology is well-established: at sufficiently high density of bacteria, there appear large-scale vortices and jets comprising many individual organisms, forming a chaotic state commonly known as bacterial turbulence. As revealed by computer simulations, below the onset of collective motion, the suspension exhibits very strong correlations between individual microswimmers stemming from the long-ranged nature of dipolar fields. Here we demonstrate that this phenomenology is captured by the minimal model of microswimmers. We develop a kinetic theory that goes beyond the commonly used mean-field assumption, and explicitly takes into account such correlations. Notably, these can be computed exactly within our theory. We calculate the fluid velocity variance, spatial and temporal correlation functions, the fluid velocity spectrum, and the enhanced diffusivity of tracer particles. We find that correlations are suppressed by particle self-propulsion, although the mean-field behaviour is not restored even in the limit of very fast swimming.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Debye-Hückel potential at an interface between two media

Electrostatic interactions between point charges embedded into interfaces separating dielectric media are omnipresent in soft matter systems and often control their stability. Such interactions are typically complicated and do not resemble their bulk counterparts. For instance, the electrostatic potential of a point charge at an air-water interface falls off as $r^{-3}$, where $r$ is the distance from the charge, exhibiting a dipolar behaviour. This behaviour is often assumed to be generic, and is widely referred to when interpreting experimental results. Here we explicitly calculate the in-plane potential of a point charge at an interface between two electrolyte solutions with different dielectric permittivities and Debye screening lengths. We show that the asymptotic behaviour of this potential is neither a dipole, which characterises the potential at air-water interfaces, nor a screened monopole, which describes the bulk behaviour in a single electrolyte solution. By considering the same problem in arbitrary dimensions, we find that the physics behind this difference can be traced to the asymmetric propagation of the interaction in the two media. Our results are relevant, for instance, to understand the physics of charged colloidal particles trapped at oil-water interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Kinetic roughening in active interfaces

The essential features of many interfaces driven out of equilibrium are described by the same equation---the Kardar-Parisi-Zhang (KPZ) equation. How do living interfaces, such as the cell membrane, fit into this picture? In an endeavour to answer such a question, we proposed in [F. Cagnetta, M. R. Evans, D. Marenduzzo, PRL 120, 258001 (2018)] an idealised model for the membrane of a moving cell. Here we discuss how the addition of simple ingredients inspired by the dynamics of the membrane of moving cells affects common kinetic roughening theories such as the KPZ and Edwards-Wilkinson equations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Magnetic Polymer Models for Epigenetics-Driven Chromosome Folding

Epigenetics is a driving force of important and ubiquitous phenomena in nature such as cell differentiation or even metamorphosis. Oppositely to its widespread role, understanding the biophysical principles that allow epigenetics to control and rewire gene regulatory networks remains an open challenge. In this work we study the effects of epigenetic modifications on the spatial folding of chromosomes -- and hence on the expression of the underlying genes -- by mapping the problem to a class of models known as magnetic polymers. In this work we show that a first-order phase transition underlies the simultaneous spreading of certain epigenetic marks and the conformational collapse of a chromosome. Further, we describe Brownian Dynamics simulations of the model in which the topology of the polymer and thermal fluctuations are fully taken into account and that confirms our mean-field predictions. Extending our models to allow for non-equilibrium terms yields new stable phases which qualitatively agrees with observations in vivo. Our results show that statistical mechanics techniques applied to models of magnetic polymers can be successfully exploited to rationalize the outcomes of experiments designed to probe the interplay between a dynamic epigenetic landscape and chromatin organization.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Nonequilibrium Theory of Epigenomic Microphase Separation in the Cell Nucleus

Understanding the spatial organisation of the genome in the cell nucleus is one of the current grand challenges in biophysics. Certain biochemical -- or epigenetic -- marks that are deposited along the genome are thought to play an important, yet poorly understood, role in determining genome organisation and cell identity. The physical principles underlying the interplay between epigenetic dynamics and genome folding remain elusive. Here we propose and study a theory that assumes a coupling between epigenetic mark and genome densities, and which can be applied at the scale of the whole nucleus. We show that equilibrium models are not compatible with experiments and a qualitative agreement is recovered by accounting for non-equilibrium processes which can stabilise microphase separated epigenomic domains. We finally discuss the potential biophysical origin of these terms.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Lamellar ordering, droplet formation and phase inversion in exotic active emulsions

We study numerically the behaviour of a mixture of a passive isotropic fluid and an active polar gel, in the presence of a surfactant favouring emulsification. Focussing on parameters for which the underlying free energy favours the lamellar phase in the passive limit, we show that the interplay between nonequilibrium and thermodynamic forces creates a range of multifarious exotic emulsions. When the active component is contractile (e.g., an actomyosin solution), moderate activity enhances the efficiency of lamellar ordering, whereas strong activity favours the creation of passive droplets within an active matrix. For extensile activity (occurring, e.g., in microtubule-motor suspensions), instead, we observe an emulsion of spontaneously rotating droplets of different size. By tuning the overall composition, we can create high internal phase emulsions, which undergo sudden phase inversion when activity is switched off. Therefore, we find that activity provides a single control parameter to design composite materials with a strikingly rich range of morphologies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Polymer Model with Epigenetic Recolouring Reveals a Pathway for the de novo Establishment and 3D Organisation of Chromatin Domains

One of the most important problems in development is how epigenetic domains can be first established, and then maintained, within cells. To address this question, we propose a framework which couples 3D chromatin folding dynamics, to a "recolouring" process modelling the writing of epigenetic marks. Because many intra-chromatin interactions are mediated by bridging proteins, we consider a "two-state" model with self-attractive interactions between two epigenetic marks which are alike (either active or inactive). This model displays a first-order-like transition between a swollen, epigenetically disordered, phase, and a compact, epigenetically coherent, chromatin globule. If the self-attraction strength exceeds a threshold, the chromatin dynamics becomes glassy, and the corresponding interaction network freezes. By modifying the epigenetic read-write process according to more biologically-inspired assumptions, our polymer model with recolouring recapitulates the ultrasensitive response of epigenetic switches to perturbations, and accounts for multi-domain conformations, strikingly similar to the topologically-associating-domains observed in eukaryotic chromosomes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Chromosome-wide simulations uncover folding pathway and 3D organization of interphase chromosomes

Three-dimensional interphase organization of metazoan genomes has been linked to cellular identity. However, the principles governing 3D interphase genome architecture and its faithful transmission through disruptive events of cell-cycle, like mitosis, are not fully understood. By using Brownian dynamics simulations of Drosophila chromosome 3R up to time-scales of minutes, we show that chromatin binding profile of Polycomb-repressive-complex-1 robustly predicts a sub-set of topologically associated domains (TADs), and inclusion of other factors recapitulates the profile of all TADs, as observed experimentally. Our simulations show that chromosome 3R attains interphase organization from mitotic state by a two-step process in which formation of local TADs is followed by long-range interactions. Our model also explains statistical features and tracks the assembly kinetics of polycomb subnuclear clusters. In conclusion, our approach can be used to predict structural and kinetic features of 3D chromosome folding and its associated proteins in biological relevant genomic and time scales.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hydrodynamic oscillations and variable swimming speed in squirmers close to repulsive walls

We present a lattice Boltzmann study of the hydrodynamics of a fully resolved squirmer, radius R, confined in a slab of fluid between two no-slip walls. We show that the coupling between hydrodynamics and short-range repulsive interactions between the swimmer and the surface can lead to hydrodynamic trapping of both pushers and pullers at the wall, and to hydrodynamic oscillations in the case of a pusher. We further show that a pusher moves significantly faster when close to a surface than in the bulk, whereas a puller undergoes a transition between fast motion and a dynamical standstill according to the range of the repulsive interaction. Our results critically require near-field hydrodynamics; they further suggest that it should be possible to control density and speed of squirmers at a surface by tuning the range of steric and electrostatic swimmer-wall interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pattern Formation in Chemically Interacting Active Rotors with Self-Propulsion

We demonstrate that active rotations in chemically signalling particles, such as autochemotactic {\it E. coli} close to walls, create a route for pattern formation based on a nonlinear yet deterministic instability mechanism. For slow rotations, we find a transient persistence of the uniform state, followed by a sudden formation of clusters contingent on locking of the average propulsion direction by chemotaxis. These clusters coarsen, which results in phase separation into a dense and a dilute region. Faster rotations arrest phase separation leading to a global travelling wave of rotors with synchronized roto-translational motion. Our results elucidate the physics resulting from the competition of two generic paradigms in active matter, chemotaxis and active rotations, and show that the latter provides a tool to design programmable self-assembly of active matter, for example to control coarsening.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pattern formation in polymerising actin flocks: spirals, spots and waves without nonlinear chemistry

We propose a model solely based on actin treadmilling and polymerisation which describes many characteristic states of actin wave formation: spots, spirals and travelling waves. In our model, as in experiments on cell recovering motility following actin depolymerisation, we choose an isotropic low density initial condition; polymerisation of actin filaments then raises the density towards the Onsager threshold where they align. We show that this alignment, in turn, destabilizes the isotropic phase and generically induces transient actin spots or spirals as part of the dynamical pathway towards a polarized phase which can either be uniform or consist of a series of actin-wave trains (flocks). Our results uncover a universal route to actin wave formation in the absence of any system specific nonlinear biochemistry, and it may help understand the mechanism underlying the observation of actin spots and waves in vivo. They also suggest a minimal setup to design similar patterns in vitro.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Predicting the three-dimensional folding of cis-regulatory regions in mammalian genomes using bioinformatic data and polymer models

The three-dimensional organisation of chromosomes can be probed using methods such as Capture-C. However it is unclear how such population level data relates to the organisation within a single cell, and the mechanisms leading to the observed interactions are still largely obscure. We present a polymer modelling scheme based on the assumption that chromosome architecture is maintained by protein bridges which form chromatin loops. To test the model we perform FISH experiments and also compare with Capture-C data. Starting merely from the locations of protein binding sites, our model accurately predicts the experimentally observed chromatin interactions, revealing a population of 3D conformations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Active Model H: Scalar Active Matter in a Momentum-Conserving Fluid

We present a continuum theory of self-propelled particles, without alignment interactions, in a momentum-conserving solvent. To address phase separation we introduce a scalar concentration field $ϕ$ with advective-diffusive dynamics. Activity creates a contribution $Σ_{ij}=-ζ((\partial_iϕ)(\partial_jϕ)-(\nablaϕ)^{2}δ_{ij}/d)$ to the deviatoric stress, where $ζ$ is odd under time reversal and $d$ is the number of spatial dimensions; this causes an effective interfacial tension contribution that is negative for contractile swimmers. We predict that domain growth then ceases at a length scale where diffusive coarsening is balanced by active stretching of interfaces, and confirm this numerically. Thus the interplay of activity and hydrodynamics is highly nontrivial, even without alignment interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Activity-induced phase separation and self-assembly in mixtures of active and passive particles

We investigate the phase behavior and kinetics of a monodisperse mixture of active (\textit{i.e.}, self-propelled) and passive isometric Brownian particles through Brownian dynamics simulations and theory. As in a purely active system, motility of the active component triggers phase separation into a dense and a dilute phase; in the dense phase we further find active-passive segregation, with "rafts" of passive particles in a "sea" of active particles. We find that phase separation from an initially disordered mixture can occur with as little as 15 percent of the particles being active. Finally, we show that a system prepared in a suitable fully segregated initial state reproducibly self-assembles an active "corona" which triggers crystallization of the passive core by initiating a compression wave. Our findings are relevant to the experimental pursuit of directed self-assembly using active particles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Clustering and pattern formation in chemorepulsive active colloids

We demonstrate that migration away from self-produced chemicals (chemorepulsion) generates a generic route to clustering and pattern formation among self-propelled colloids. The clustering instability can be caused either by anisotropic chemical production, or by a delayed orientational response to changes of the chemical environment. In each case, chemorepulsion creates clusters of a self-limiting area which grows linearly with self-propulsion speed. This agrees with recent observations of dynamic clusters in Janus colloids (albeit not yet known to be chemorepulsive). More generally, our results could inform design principles for the self-assembly of chemorepulsive synthetic swimmers and/or bacteria into nonequilibrium patterns.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Is the Kinetoplast DNA a Percolating Network of Linked Rings at its Critical Point?

In this work we present a computational study of the Kinetoplast genome, modelled as a large number of semiflexible unknotted loops, which are allowed to link with each other. As the DNA density increases, the systems shows a percolation transition between a gas of unlinked rings and a network of linked loops which spans the whole system. Close to the percolation transition, we find that the mean valency of the network, i.e. the average number of loops which are linked to any one loop, is around 3 as found experimentally for the Kinetoplast DNA. Even more importantly, by simulating the digestion of the network by a restriction enzyme, we show that the distribution of oligomers, i.e. structures formed by a few loops which remain linked after digestion, quantitatively matches experimental data obtained from gel electrophoresis, provided that the density is, once again, close to the percolation transition. With respect to previous work, our analysis builds on a reduced number of assumptions, yet can still fully explain the experimental data. Our findings suggest that the Kinetoplast DNA can be viewed as a network of linked loops positioned very close to the percolation transition, and we discuss the possible biological implications of this remarkable fact.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Motility-induced phase separation and coarsening in active matter

Active systems, or active matter, are self-driven systems which live, or function, far from equilibrium - a paradigmatic example which we focus on here is provided by a suspension of self-motile particles. Active systems are far from equilibrium because their microscopic constituents constantly consume energy from the environment in order to do work, for instance to propel themselves. The nonequilibrium nature of active matter leads to a variety of non-trivial intriguing phenomena. An important one which has recently been the subject of intense interest among biological and soft matter physicists is that of the so-called "motility-induced phase separation", whereby self-propelled particles accumulate into clusters in the absence of any explicit attractive interactions between them. Here we review the physics of motility-induced phase separation, and discuss this phenomenon within the framework of the classic physics of phase separation and coarsening. We also discuss theories for bacterial colonies where coarsening may be arrested. Most of this work will focus on the case of run-and-tumble and active Brownian particles in the absence of solvent-mediated hydrodynamic interactions - we will briefly discuss at the end their role, which is not currently fully understood in this context.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Topological Patterns in Two-dimensional Gel Electrophoresis of DNA Knots

Gel electrophoresis is a powerful experimental method to probe the topology of DNA and other biopolymers. While there is a large body of experimental work which allows us to accurately separate different topoisomers of a molecule, a full theoretical understanding of these experiments has not yet been achieved. Here we show that the mobility of DNA knots depends crucially and subtly on the physical properties of the gel, and in particular on the presence of dangling ends. The topological interactions between these and DNA molecules can be described in terms of an "entanglement number", and yield a non-monotonic mobility at moderate fields. Consequently, in two-dimensional electrophoresis, gel bands display a characteristic arc pattern; this turns into a straight line when the density of dangling ends vanishes. We also provide a novel framework to accurately predict the shape of such arcs as a function of molecule length and topological complexity, which may be used to inform future experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Active polar fluid flow in finite droplets

We present a continuum level analytical model of a droplet of active contractile fluid consisting of filaments and motors. We calculate the steady state flows that result from a splayed polarisation of the filaments. We account for the interaction with an arbitrary external medium by imposing a viscous friction at the fixed droplet boundary. We then show that the droplet has non-zero force dipole and quadrupole moments, the latter of which is essential for self-propelled motion of the droplet at low Reynolds' number. Therefore, this calculation describes a simple mechanism for the motility of a droplet of active contractile fluid embedded in a 3D environment, which is relevant to cell migration in confinement (for example, embedded within a gel or tissue). Our analytical results predict how the system depends on various parameters such as the effective friction coefficient, the phenomenological activity parameter and the splay of the imposed polarisation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Activity-induced clustering in model dumbbell swimmers: The role of hydrodynamic interactions

Using a fluid-particle dynamics approach, we numerically study the effects of hydrodynamic interactions on the collective dynamics of active suspensions within a simple model for bacterial motility: each microorganism is modeled as a stroke-averaged dumbbell swimmer with prescribed dipolar force pairs. Using both simulations and qualitative arguments, we show that, when the separation between swimmers is comparable to their size, the swimmers' motions are strongly affected by activity-induced hydrodynamic forces. To further understand these effects, we investigate semidilute suspensions of swimmers in the presence of thermal fluctuations. A direct comparison between simulations with and without hydrodynamic interactions shows these to enhance the dynamic clustering at a relatively small volume fraction; with our chosen model the key ingredient for this clustering behavior is hydrodynamic trapping of one swimmer by another, induced by the active forces. Furthermore, the density dependence of the motility (of both the translational and rotational motions) exhibits distinctly different behaviors with and without hydrodynamic interactions; we argue that this is linked to the clustering tendency. Our study illustrates the fact that hydrodynamic interactions not only affect kinetic pathways in active suspensions, but also cause major changes in their steady state properties.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamics of Self-Threading Ring Polymers in a Gel

We study of the dynamics of ring polymers confined to diffuse in a background gel at low concentrations. We do this in order to probe the inter-play between topology and dynamics in ring polymers. We develop an algorithm that takes into account the possibility that the rings hinder their own motion by passing through themselves, i.e. "self-threading". Our results suggest that the number of self-threadings scales extensively with the length of the rings and that this is substantially independent of the details of the model. The slowing down of the rings' dynamics is found to be related to the fraction of segments that can contribute to the motion. Our results give a novel perspective on the motion of ring polymers in gel, for which a complete theory is still lacking, and may help us to understand the irreversible trapping of ring polymers in gel electrophoresis experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase behaviour of active Brownian particles: The role of dimensionality

Recently, there has been much interest in activity-induced phase separations in concentrated suspensions of "active Brownian particles" (ABPs), self-propelled spherical particles whose direction of motion relaxes through thermal rotational diffusion. To date, almost all these studies have been restricted to 2 dimensions. In this work we study activity-induced phase separation in 3D and compare the results with previous and new 2D simulations. To this end, we performed state-of-the-art Brownian dynamics simulations of up to 40 million ABPs -- such very large system sizes are unavoidable to evade finite size effects in 3D. Our results confirm the picture established for 2D systems in which an activity-induced phase separation occurs, with strong analogies to equilibrium gas-liquid spinodal decomposition, in spite of the purely non-equilibrium nature of the driving force behind the phase separation. However, we also find important differences between the 2D and 3D cases. Firstly, the shape and position of the phase boundaries is markedly different for the two cases. Secondly, for the 3D coarsening kinetics we find that the domain size grows in time according to the classical diffusive $t^{1/3}$ law, in contrast to the nonstandard subdiffusive exponent observed in 2D.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scalar $ϕ^4$ field theory for active-particle phase separation

Recent theories predict phase separation among orientationally disordered active particles whose propulsion speed decreases rapidly enough with density. Coarse-grained models of this process show time-reversal symmetry (detailed balance) to be restored for uniform states, but broken by gradient terms; hence detailed-balance violation is strongly coupled to interfacial phenomena. To explore the subtle generic physics resulting from such coupling we here introduce `Active Model B'. This is a scalar $ϕ^4$ field theory (or phase-field model) that minimally violates detailed balance via a leading-order square-gradient term. We find that this additional term has modest effects on coarsening dynamics, but alters the static phase diagram by creating a jump in (thermodynamic) pressure across flat interfaces. Both results are surprising, since interfacial phenomena are always strongly implicated in coarsening dynamics but are, in detailed-balance systems, irrelevant for phase equilibria.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Topology Regulation during Replication of the Kinetoplast DNA

We study theoretically the replication of Kinetoplast DNA consisting of several thousands separate mini-circles found in organisms of the class Kinetoplastida. When the cell is not actively dividing these are topologically connected in a marginally linked network of rings with only one connected component. During cell division each mini-circle is removed from the network, duplicated and then re-attached, along with its progeny. We study this process under the hypothesis that there is a coupling between the topological state of the mini-circles and the expression of genetic information encoded on them, leading to the production of Topoisomerase. This model describes a self-regulating system capable of full replication that reproduces several previous experimental findings. We find that the fixed point of the system depends on a primary free parameter of the model: the ratio between the rate of removal of mini-circles from the network (R) and their (re)attachment rate (A). The final topological state is found to be that of a marginally linked network structure in which the fraction of mini-circles linked to the largest connected component approaches unity as R/A decreases. Finally we discuss how this may suggest an evolutionary trade-off between the speed of replication and the accuracy with which a fully topologically linked state is produced.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A continuum theory of phase separation kinetics for active Brownian particles

Active Brownian particles (ABPs), when subject to purely repulsive interactions, are known to undergo activity-induced phase separation broadly resembling an equilibrium (attraction-induced) gas-liquid coexistence. Here we present an accurate continuum theory for the dynamics of phase-separating ABPs, derived by direct coarse-graining, capturing leading-order density gradient terms alongside an effective bulk free energy. Such gradient terms do not obey detailed balance; yet we find coarsening dynamics closely resembling that of equilibrium phase separation. Our continuum theory is numerically compared to large-scale direct simulations of ABPs and accurately accounts for domain growth kinetics, domain topologies and coexistence densities.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Enhanced diffusion of tracer particles in dilute bacterial suspensions

Swimming bacteria create long-range velocity fields that stir a large volume of fluid and move around passive particles dispersed in the fluid. Recent experiments and simulations have shown that long-time mean-squared displacement of passive particles in a bath of swimming bacteria exhibits diffusive behaviour with the effective diffusion coefficient significantly larger than its thermal counterpart. Comprehensive theoretical prediction of this effective diffusion coefficient and understanding of the enhancement mechanism remain a challenge. Here, we adapt the kinetic theory by Lin et al., J. Fluid Mech. 669, 167 (2011) developed for 'squirmers' to the bacterial case to quantitatively predict enhanced diffusivity of tracer particles in dilute two- and three-dimensional suspensions of swimming bacteria. We demonstrate that the effective diffusion coefficient is a product of the bacterial number density, their swimming speed, a geometric factor characterising the velocity field created by a single bacterium, and a numerical factor. We show that the numerical factor is, in fact, a rather strong function of the system parameters, most notably the run length of the bacteria, and that these dependencies have to be taken into account to quantitatively predict the enhanced diffusivity. We perform molecular-dynamics-type simulations to confirm the conclusions of the kinetic theory. Our results are in a good agreement with the values of enhanced diffusivity measured in recent two- and three-dimensional experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modelling the effect of myosin X motors on filopodia growth

We present a numerical simulation study of the dynamics of filopodial growth in the presence of active transport by myosin X motors. We employ both a microscopic agent-based model, which captures the stochasticity of the growth process, and a continuum mean-field theory which neglects fluctuations. We show that in the absence of motors, filopodia growth is overestimated by the continuum mean-field theory. Thus fluctuations slow down the growth, especially when the protrusions are driven by a small number (10 or less) of F-actin fibres, and when the force opposing growth (coming from membrane elasticity) is large enough. We also show that, with typical parameter values for eukaryotic cells, motors are unlikely to provide an actin transport mechanism which enhances filopodial size significantly, unless the G-actin concentration within the filopodium greatly exceeds that of the cytosol bulk. We explain these observations in terms of order-of-magnitude estimates of diffusion-induced and advection-induced growth of a bundle of Brownian ratchets.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rheology of Cubic Blue Phases

We study the behaviour of cubic blue phases under shear flow via lattice Boltzmann simulations. We focus on the two experimentally observed phases, Blue Phase I (BPI) and Blue Phase II (BPII). The disclination network of Blue Phase II continuously breaks and reforms under steady shear, leading to an oscillatory stress response in time. For larger shear rates, the structure breaks up into a Grandjean texture with a cholesteric helix lying along the flow gradient direction. Blue Phase I leads to a very different response. Here, oscillations are only possible for intermediate shear rates -- very slow flow causes a transition of the initially ordered structure into an amorphous network with an apparent yield stress. Larger shear rates lead to another amorphous state with different structure of the defect network. For even larger flow rates the same break-up into a Grandjean texture as for Blue Phase II is observed. At the highest imposed flow rates both cubic blue phases adopt a flow-aligned nematic state. Our results provide the first theoretical investigation of sheared blue phases in large systems, and are relevant to understanding the bulk rheology of these materials.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Sedimentation of Knotted Polymers

We investigate the sedimentation of knotted polymers by means of stochastic rotation dynamics, a molecular dynamics algorithm that takes hydrodynamics fully into account. We show that the sedimentation coefficient s, related to the terminal velocity of the knotted polymers, increases linearly with the average crossing number n_c of the corresponding ideal knot. To the best of our knowledge, this provides the first direct computational confirmation of this relation, postulated on the basis of experiments in "The effect of ionic conditions on the conformations of supercoiled DNA. I. sedimentation analysis" by Rybenkov et al., for the case of sedimentation. Such a relation was previously shown to hold with simulations for knot electrophoresis. We also show that there is an accurate linear dependence of s on the inverse of the radius of gyration R_g^-1, more specifically with the inverse of the R_g component that is perpendicular to the direction along which the polymer sediments. When the polymer sediments in a slab, the walls affect the results appreciably. However, R_g^-1 remains to a good precision linearly dependent on n_c. Therefore, R_g^-1 is a good measure of a knot's complexity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Colloids in a bacterial bath: simulations and experiments

We present a joint experimental and computational study of the effect of bacterial motion on micron-scale colloids contained in a two-dimensional suspension of Bacillus subtilis. With respect to previous work using E. coli, here we introduce a novel experimental set-up that allows us to realise a two-dimensional bacterial suspension insensitive to either evaporation or fluid flow. By analysing the mean square displacements of both bacteria and colloids, we confirm the existence of a crossover from super-diffusive behaviour at short time scales to normal diffusion at longer times. We also study the same two-dimensional system by means of numerical simulations, using a suspension of self-propelled dumbbells or the Vicsek model, which has been previously used to study the dynamics of active particles. Our numerical results obtained with both models are in broad agreement with the experimental trends, but only the dumbbell simulations can match the experimental data quantitatively. The level of agreement we find suggest that steric interactions due to collisions are important players in determining collective motion of the bacterial bath, and should complement hydrodynamic interactions in experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Polymers with spatial or topological constraints: theoretical and computational results

In this review we provide an organized summary of the theoretical and computational results which are available for polymers subject to spatial or topological constraints. Because of the interdisciplinary character of the topic, we provide an accessible, non-specialist introduction to the main topological concepts, polymer models, and theoretical/computational methods used to investigate dense and entangled polymer systems. The main body of our review deals with: (i) the effect that spatial confinement has on the equilibrium topological entanglement of one or more polymer chains and (ii) the metric and entropic properties of polymer chains with fixed topological state. These problems have important technological applications and implications for the life-sciences. Both aspects, especially the latter, are amply covered. A number of selected open problems are finally highlighted.

0 citations0 reviewsNo code linkedOpen access