Showing 1–50 of 147 results for author: Marchetti, M C

Asymmetric fluctuations and self-folding of active interfaces

Authors: Liang Zhao, Paarth Gulati, Fernando Caballero, Itamar Kolvin, Raymond Adkins, M. Cristina Marchetti, Zvonimir Dogic

Abstract: We study the structure and dynamics of the interface separating a passive fluid from a microtubule-based active fluid. Turbulent-like active flows power giant interfacial fluctuations, which exhibit pronounced asymmetry between regions of positive and negative curvature. Experiments, numerical simulations, and theoretical arguments reveal how the interface breaks up the spatial symmetry of the fun… ▽ More We study the structure and dynamics of the interface separating a passive fluid from a microtubule-based active fluid. Turbulent-like active flows power giant interfacial fluctuations, which exhibit pronounced asymmetry between regions of positive and negative curvature. Experiments, numerical simulations, and theoretical arguments reveal how the interface breaks up the spatial symmetry of the fundamental bend instability to generate local vortical flows that lead to asymmetric interface fluctuations. The magnitude of interface deformations increases with activity: In the high activity limit, the interface self-folds invaginating passive droplets and generating a foam-like phase, where active fluid is perforated with passive droplets. These results demonstrate how active stresses control the structure, dynamics, and break-up of soft, deformable, and reconfigurable liquid-liquid interfaces. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 26 pages, 22 figures

Nonreciprocal interactions drive capillary waves in active liquid crystals

Authors: Paarth Gulati, Fernando Caballero, Itamar Kolvin, Zhihong You, M. Cristina Marchetti

Abstract: Active liquid crystals exert nonequilibrium stresses on their surroundings through constant consumption of energy, giving rise to dynamical steady states not present in equilibrium. The paradigmatic example of an active liquid crystal is a suspension of microtubule bundles powered by kinesin motor proteins, which exhibits self-sustained spatiotemporal chaotic flows. This system has been modelled u… ▽ More Active liquid crystals exert nonequilibrium stresses on their surroundings through constant consumption of energy, giving rise to dynamical steady states not present in equilibrium. The paradigmatic example of an active liquid crystal is a suspension of microtubule bundles powered by kinesin motor proteins, which exhibits self-sustained spatiotemporal chaotic flows. This system has been modelled using continuum theories that couple the microtubule orientation to active flows. Recently the focus has shifted to the interfacial properties of mixtures of active liquid crystals and passive fluids. Active/passive interfaces have been shown to support propagating capillary waves in the absence of inertia and offer a promising route for relating experimental parameters to those of the continuum theory. In this paper we report the derivation of a minimal model that captures the linear dynamics of the interface between an active liquid crystal and a passive fluid. We show that the dynamics of the interface, although powered by active flows throughout the bulk, is qualitatively captured by equations that couple non-reciprocally interface height and nematic director at the interface. This minimal model reproduces the dynamical structure factor evaluated from numerical simulations and the qualitative form of the wave dispersion relation seen in experiments. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: 14 pages, 7 figures

Multiphase Field Model of Cells on a Substrate

Authors: Michael Chiang, Austin Hopkins, Benjamin Loewe, Davide Marenduzzo, M. Cristina Marchetti

Abstract: Multiphase field models have emerged as an important computational tool for understanding biological tissue while resolving single-cell properties. While they have successfully reproduced many experimentally observed behaviors of living tissue, the theoretical underpinnings have not been fully explored. We show that a two-dimensional version of the model, which is commonly employed to study tissue… ▽ More Multiphase field models have emerged as an important computational tool for understanding biological tissue while resolving single-cell properties. While they have successfully reproduced many experimentally observed behaviors of living tissue, the theoretical underpinnings have not been fully explored. We show that a two-dimensional version of the model, which is commonly employed to study tissue monolayers, can be derived from a three-dimensional version in the presence of a substrate. We also show how viscous forces, which arise from friction between different cells, can be included in the model. Finally, we numerically simulate a tissue monolayer, and find that intercellular friction tends to solidify the tissue. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Comments: 14 pages, 2 figures

Intercellular Friction and Motility Drive Orientational Order in Cell Monolayers

Authors: Michael Chiang, Austin Hopkins, Benjamin Loewe, M. Cristina Marchetti, Davide Marenduzzo

Abstract: Spatiotemporal patterns in multicellular systems are important to understanding tissue dynamics, for instance, during embryonic development and disease. Here, we use a multiphase field model to study numerically the behavior of a near-confluent monolayer of deformable cells with intercellular friction. Varying friction and cell motility drives a solid-liquid transition, and near the transition bou… ▽ More Spatiotemporal patterns in multicellular systems are important to understanding tissue dynamics, for instance, during embryonic development and disease. Here, we use a multiphase field model to study numerically the behavior of a near-confluent monolayer of deformable cells with intercellular friction. Varying friction and cell motility drives a solid-liquid transition, and near the transition boundary, we find the emergence of nematic order of cell deformation driven by shear-aligning cellular flows. Intercellular friction endows the monolayer with a finite viscosity, which significantly increases the spatial correlation in the flow and, concomitantly, the extent of nematic order. We also show that hexatic and nematic order are tightly coupled and propose a mechanical-geometric model for the colocalization of +1/2 nematic defects and 5-7 disclination pairs, which are the structural defects in the hexatic phase. Such topological defects coincide with regions of high cell-cell overlap, suggesting that they may mediate cellular extrusion from the monolayer, as found experimentally. Our results delineate a mechanical basis for the recent observation of nematic and hexatic order in multicellular collectives in experiments and simulations and pinpoint a generic pathway to couple topological and physical effects in these systems. △ Less

Submitted 31 October, 2023; originally announced October 2023.

Motility induced phase separation of deformable cells

Authors: Austin Hopkins, Benjamin Loewe, Michael Chiang, Davide Marenduzzo, M. Cristina Marchetti

Abstract: Using a multi-phase field model, we examine how particle deformability, which is a proxy for cell stiffness, affects motility induced phase separation (MIPS). We show that purely repulsive deformable, i.e., squishy, cells phase separate more effectively than their rigid counterparts. This can be understood as due to the fact that deformability increases the effective duration of collisions. In add… ▽ More Using a multi-phase field model, we examine how particle deformability, which is a proxy for cell stiffness, affects motility induced phase separation (MIPS). We show that purely repulsive deformable, i.e., squishy, cells phase separate more effectively than their rigid counterparts. This can be understood as due to the fact that deformability increases the effective duration of collisions. In addition, the dense regions become increasingly disordered as deformability increases. Our results contextualize the applicability of MIPS to biological systems and have implications for how cells in biological systems may self-organize △ Less

Submitted 9 August, 2023; originally announced August 2023.

Comments: 7 pages, 6 figures

Nonreciprocal pattern formation of conserved fields

Authors: Fridtjof Brauns, M. Cristina Marchetti

Abstract: In recent years, nonreciprocally coupled systems have received growing attention. Previous work has shown that the interplay of nonreciprocal coupling and Goldstone modes can drive the emergence of temporal order such as traveling waves. We show that these phenomena are generically found in a broad class of pattern-forming systems, including mass-conserving reaction--diffusion systems and viscoela… ▽ More In recent years, nonreciprocally coupled systems have received growing attention. Previous work has shown that the interplay of nonreciprocal coupling and Goldstone modes can drive the emergence of temporal order such as traveling waves. We show that these phenomena are generically found in a broad class of pattern-forming systems, including mass-conserving reaction--diffusion systems and viscoelastic active gels. All these systems share a characteristic dispersion relation that acquires a non-zero imaginary part at the edge of the band of unstable modes and exhibit a regime of propagating structures (traveling wave bands or droplets). We show that models for these systems can be mapped to a common "normal form" that can be seen as a spatially extended generalization of the FitzHugh--Nagumo model, providing a unifying dynamical-systems perspective. We show that the minimal nonreciprocal Cahn--Hilliard (NRCH) equations exhibit a surprisingly rich set of behaviors, including interrupted coarsening of traveling waves without selection of a preferred wavelength and transversal undulations of wave fronts in two dimensions. We show that the emergence of traveling waves and their speed are precisely predicted from the local dispersion relation at interfaces far away from the homogeneous steady state. Our work thus generalizes previously studied nonreciprocal phase transitions and shows that interfaces are the relevant collective excitations governing the rich dynamical patterns of conserved fields. △ Less

Submitted 21 April, 2024; v1 submitted 15 June, 2023; originally announced June 2023.

Comments: 21 pages, 12 figures (including appendices), videos: https://youtube.com/playlist?list=PLdejNN8-Lt6HyU9G9eB8Pp0-dyJvg9AJs; v2: minor changes in introduction and discussion; v3 major revision; v4 minor changes, final accepted version

Journal ref: Phys. Rev. X 14, 021014 (2024)

Vorticity phase separation and defect lattices in the isotropic phase of active liquid crystals

Authors: Fernando Caballero, Zhihong You, M. Cristina Marchetti

Abstract: We use numerical simulations and linear stability analysis to study the dynamics of an active liquid crystal film on a substrate in the regime where the passive system would be isotropic. Extensile activity builds up local orientational order and destabilizes the quiescent isotropic state above a critical activity value, eventually resulting in spatiotemporal chaotic dynamics akin to the one obser… ▽ More We use numerical simulations and linear stability analysis to study the dynamics of an active liquid crystal film on a substrate in the regime where the passive system would be isotropic. Extensile activity builds up local orientational order and destabilizes the quiescent isotropic state above a critical activity value, eventually resulting in spatiotemporal chaotic dynamics akin to the one observed ubiquitously in the nematic state. Here we show that tuning substrate friction yields a variety of emergent structures at intermediate activity, including lattices of flow vortices with associated regular arrangements of topological defects and a new state where flow vortices trap pairs of $+1/2$ defect that chase each other tail. These chiral units spontaneously pick the sense of rotation and organize in a hexagonal lattice, surrounded by a diffuse flow of opposite rotation to maintain zero net vorticity. The length scale of these emergent structures is set by the screening length $l_η=\sqrt{η/Γ}$ of the flow, controlled by the shear viscosity $η$ and the substrate friction $Γ$, and can be captured by simple mode selection of the vortical flows. We demonstrate that the emergence of coherent structures can be interpreted as a phase separation of vorticity, where friction plays a role akin to that of birth/death processes in breaking conservation of the phase separating species and selecting a characteristic scale for the patterns. Our work shows that friction provides an experimentally accessible tuning parameter for designing controlled active flows. △ Less

Submitted 22 September, 2023; v1 submitted 7 June, 2023; originally announced June 2023.

Finite Elasticity of the Vertex Model and its Role in Rigidity of Curved Cellular Tissues

Authors: Arthur Hernandez, Michael F. Staddon, Michael Moshe, M. Cristina Marchetti

Abstract: Using a mean field approach and simulation, we study the non-linear mechanical response of the vertex model (VM) of biological tissue under compression and dilation. The VM is known to exhibit a transition between rigid and fluid-like, or floppy, states driven by geometric incompatibility. Target perimeter and area set a target shape which may not be geometrically achievable, thereby engendering f… ▽ More Using a mean field approach and simulation, we study the non-linear mechanical response of the vertex model (VM) of biological tissue under compression and dilation. The VM is known to exhibit a transition between rigid and fluid-like, or floppy, states driven by geometric incompatibility. Target perimeter and area set a target shape which may not be geometrically achievable, thereby engendering frustration. Previously, an asymmetry in the linear elastic response was identified at the rigidity transition between compression and dilation. Here we show and characterize how the asymmetry extends away from the transition point for finite strains. Under finite compression, an initially solid VM can totally relax perimeter tension, and thereby have reduced bulk and shear modulus. Conversely, an initially floppy VM under dilation can rigidify and have a higher bulk and shear modulus. These observations imply that re-scaling of cell area shifts the transition between rigid and floppy states. Based on this insight, we calculate the re-scaling of cell area engendered by intrinsic curvature and write a prediction for the rigidity transition in the presence of curvature. The shift of the rigidity transition in the presence of curvature for the VM provides a new metric for predicting tissue rigidity from image data for curved tissues in a manner analogous to the flat case. △ Less

Submitted 10 March, 2023; originally announced March 2023.

Comments: 10 pages, 3 figures

Patterning of morphogenetic anisotropy fields

Authors: Zihang Wang, M. Cristina Marchetti, Fridtjof Brauns

Abstract: Orientational order, encoded in anisotropic fields, plays an important role during the development of an organism. A striking example of this is the freshwater polyp Hydra, where topological defects in the muscle fiber orientation have been shown to localize to key features of the body plan. This body plan is organized by morphogen concentration gradients, raising the question how muscle fiber ori… ▽ More Orientational order, encoded in anisotropic fields, plays an important role during the development of an organism. A striking example of this is the freshwater polyp Hydra, where topological defects in the muscle fiber orientation have been shown to localize to key features of the body plan. This body plan is organized by morphogen concentration gradients, raising the question how muscle fiber orientation, morphogen gradients and body shape interact. Here, we introduce a minimal model that couples nematic orientational order to the gradient of a morphogen field. We show that on a planar surface alignment to a radial concentration gradient can induce unbinding of topological defects, as observed during budding and tentacle formation in Hydra, and stabilize aster/vortex-like defects, as observed at a Hydra's mouth. On curved surfaces mimicking the morphologies of Hydra in various stages of development -- from spheroid to adult -- our model reproduces the experimentally observed reorganization of orientational order. Our results suggest how gradient alignment and curvature effects may work together to control orientational order during development and lays the foundations for future modeling efforts that will include the tissue mechanics that drive shape deformations. △ Less

Submitted 22 March, 2023; v1 submitted 23 December, 2022; originally announced December 2022.

Comments: v2: Published version with minor revisions. v1: 12 pages, 7 figures (4 main, 3 suppl.). Movies: https://youtube.com/playlist?list=PLdejNN8-Lt6HJ52YfhyegpI0PGoRq8pSb

Journal ref: PNAS 120 (13), e2220167120 (2023)

Design rules for controlling active topological defects

Authors: Suraj Shankar, Luca V. D. Scharrer, Mark J. Bowick, M. Cristina Marchetti

Abstract: Topological defects play a central role in the physics of many materials, including magnets, superconductors and liquid crystals. In active fluids, defects become autonomous particles that spontaneously propel from internal active stresses and drive chaotic flows stirring the fluid. The intimate connection between defect textures and active flow suggests that properties of active materials can be… ▽ More Topological defects play a central role in the physics of many materials, including magnets, superconductors and liquid crystals. In active fluids, defects become autonomous particles that spontaneously propel from internal active stresses and drive chaotic flows stirring the fluid. The intimate connection between defect textures and active flow suggests that properties of active materials can be engineered by controlling defects, but design principles for their spatiotemporal control remain elusive. Here we propose a symmetry-based additive strategy for using elementary activity patterns, as active topological tweezers, to create, move and braid such defects. By combining theory and simulations, we demonstrate how, at the collective level, spatial activity gradients act like electric fields which, when strong enough, induce an inverted topological polarization of defects, akin to a negative susceptibility dielectric. We harness this feature in a dynamic setting to collectively pattern and transport interacting active defects. Our work establishes an additive framework to sculpt flows and manipulate active defects in both space and time, paving the way to design programmable active and living materials for transport, memory and logic. △ Less

Submitted 2 May, 2024; v1 submitted 1 December, 2022; originally announced December 2022.

Comments: 11 pages (including Methods), 5 figures. Changed title and format, final version

Report number: 2212.00666

Journal ref: Proc. Nat. Acad. Sci. 121 (21) e2400933121 (2024)

The role of non-affine deformations in the elastic behavior of the cellular vertex model

Authors: Michael F. Staddon, Arthur Hernandez, Mark J. Bowick, Michael Moshe, M. Cristina Marchetti

Abstract: The vertex model of epithelia describes the apical surface of a tissue as a tiling of polygonal cells, with a mechanical energy governed by deviations in cell shape from preferred, or target, area, $A_0$, and perimeter, $P_0$. The model exhibits a rigidity transition driven by geometric incompatibility as tuned by the target shape index, $p_0 = P_0 / \sqrt{A_0}$. For… ▽ More The vertex model of epithelia describes the apical surface of a tissue as a tiling of polygonal cells, with a mechanical energy governed by deviations in cell shape from preferred, or target, area, $A_0$, and perimeter, $P_0$. The model exhibits a rigidity transition driven by geometric incompatibility as tuned by the target shape index, $p_0 = P_0 / \sqrt{A_0}$. For $p_0 > p_*(6) = \sqrt{8 \sqrt{3}} \approx 3.72$, with $p_*(6)$ the perimeter of a regular hexagon of unit area, a cell can simultaneously attain both the preferred area and preferred perimeter. As a result, the tissue is in a mechanically soft compatible state, with zero shear and Young's moduli. For $p_0 < p_*(6)$, it is geometrically impossible for any cell to realize the preferred area and perimeter simultaneously, and the tissue is in an incompatible rigid solid state. Using a mean-field approach, we present a complete analytical calculation of the linear elastic moduli of an ordered vertex model. We analyze a relaxation step that includes non-affine deformations, leading to a softer response than previously reported. The origin of the vanishing shear and Young's moduli in the compatible state is the presence of zero-energy deformations of cell shape. The bulk modulus exhibits a jump discontinuity at the transition and can be lower in the rigid state than in the fluid-like state. The Poisson's ratio can become negative which lowers the bulk and Young's moduli. Our work provides a unified treatment of linear elasticity for the vertex model and demonstrates that this linear response is protocol-dependent. △ Less

Submitted 11 April, 2023; v1 submitted 30 November, 2022; originally announced November 2022.

Journal ref: Soft Matter, 2023

Molecular-scale substrate anisotropy and crowding drive long-range nematic order of cell monolayers

Authors: Yimin Luo, Mengyang Gu, Minwook Park, Xinyi Fang, Younghoon Kwon, Juan Manuel Urueña, Javier Read de Alaniz, Matthew E. Helgeson, M. Cristina Marchetti, Megan T. Valentine

Abstract: The ability of cells to reorganize in response to external stimuli is important in areas ranging from morphogenesis to tissue engineering. Elongated cells can co-align due to steric effects, forming states with local order. We show that molecular-scale substrate anisotropy can direct cell organization, resulting in the emergence of nematic order on tissue scales. To quantitatively examine the diso… ▽ More The ability of cells to reorganize in response to external stimuli is important in areas ranging from morphogenesis to tissue engineering. Elongated cells can co-align due to steric effects, forming states with local order. We show that molecular-scale substrate anisotropy can direct cell organization, resulting in the emergence of nematic order on tissue scales. To quantitatively examine the disorder-order transition, we developed a high-throughput imaging platform to analyze velocity and orientational correlations for several thousand cells over days. The establishment of global, seemingly long-ranged order is facilitated by enhanced cell division along the substrate's nematic axis, and associated extensile stresses that restructure the cells' actomyosin networks. Our work, which connects to a class of systems known as active dry nematics, provides a new understanding of the dynamics of cellular remodeling and organization in weakly interacting cell collectives. This enables data-driven discovery of cell-cell interactions and points to strategies for tissue engineering. △ Less

Submitted 24 October, 2022; originally announced October 2022.

Comments: 29 pages, 7 figures

Constitutive model for the rheology of biological tissue

Authors: Suzanne M. Fielding, James O. Cochran, Junxiang Huang, Dapeng Bi, M. Cristina Marchetti

Abstract: The rheology of biological tissue is key to processes such as embryo development, wound healing and cancer metastasis. Vertex models of confluent tissue monolayers have uncovered a spontaneous liquid-solid transition tuned by cell shape; and a shear-induced solidification transition of an initially liquid-like tissue. Alongside this jamming/unjamming behaviour, biological tissue also displays an i… ▽ More The rheology of biological tissue is key to processes such as embryo development, wound healing and cancer metastasis. Vertex models of confluent tissue monolayers have uncovered a spontaneous liquid-solid transition tuned by cell shape; and a shear-induced solidification transition of an initially liquid-like tissue. Alongside this jamming/unjamming behaviour, biological tissue also displays an inherent viscoelasticity, with a slow time and rate dependent mechanics. With this motivation, we combine simulations and continuum theory to examine the rheology of the vertex model in nonlinear shear across a full range of shear rates from quastistatic to fast, elucidating its nonlinear stress-strain curves after the inception of shear of finite rate, and its steady state flow curves of stress as a function of strain rate. We formulate a rheological constitutive model that couples cell shape to flow and captures both the tissue solid-liquid transition and its rich linear and nonlinear rheology. △ Less

Submitted 29 September, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

Comments: 5 pages, 3 figures; plus supplemental material

Defect self-propulsion in active nematic films with spatially-varying activity

Authors: Jonas Rønning, M. Cristina Marchetti, Luiza Angheluta

Abstract: We study the dynamics of topological defects in active nematic films with spatially-varying activity and consider two setups: i) a constant activity gradient, and ii) a sharp jump in activity. A constant gradient of extensile (contractile) activity endows the comet-like $+1/2$ defect with a finite vorticity that drives the defect to align its nose in the direction of decreasing (increasing) gradie… ▽ More We study the dynamics of topological defects in active nematic films with spatially-varying activity and consider two setups: i) a constant activity gradient, and ii) a sharp jump in activity. A constant gradient of extensile (contractile) activity endows the comet-like $+1/2$ defect with a finite vorticity that drives the defect to align its nose in the direction of decreasing (increasing) gradient. A constant gradient does not, however, affect the known self-propulsion of the $+1/2$ defect and has no effect on the $-1/2$ that remains a non-motile particle. A sharp jump in activity acts like a wall that traps the defects, affecting the translational and rotational motion of both charges. The $+1/2$ defect slows down as it approaches the interface and the net vorticity tends to reorient the defect polarization so that it becomes perpendicular to the interface. The $-1/2$ defect acquires a self-propulsion towards the activity interface, while the vorticity-induced active torque tends to align the defect to a preferred orientation. This effective attraction of the negative defects to the wall is consistent with the observation of an accumulation of negative topological charge at both active/passive interfaces and physical boundaries. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Comments: 14 pages, 8 figures

Controlling liquid-liquid phase behavior with an active fluid

Authors: Alexandra M. Tayar, Fernando Caballero, Trevor Anderberg, Omar A. Saleh, M. Cristina Marchetti, Zvonimir Dogic

Abstract: Demixing of binary liquids is a ubiquitous transition, which is explained using a well-established thermodynamic formalism that requires equality of intensive thermodynamics parameters across the phase boundaries. Demixing transitions also occur when binary fluid mixtures are driven away from equilibrium, for example, by external shear flow. Predicting demixing transition under non-equilibrium non… ▽ More Demixing of binary liquids is a ubiquitous transition, which is explained using a well-established thermodynamic formalism that requires equality of intensive thermodynamics parameters across the phase boundaries. Demixing transitions also occur when binary fluid mixtures are driven away from equilibrium, for example, by external shear flow. Predicting demixing transition under non-equilibrium non-potential conditions remains, however, a challenge. We drive liquid-liquid phase separation (LLPS) of attractive DNA nanostar molecules away from equilibrium using an internally driven microtubule-based active fluid. Activity lowers the critical temperature and narrows the coexistence concentrations, but only when there are mechanical bonds between the liquid droplets and the reconfiguring active fluid. Similar behaviors are observed in numerical simulations, suggesting that activity suppression of liquid-liquid phase separation is a generic feature of active LLPS. Our work describes a platform for building soft active materials with feedback control while also providing insight into cell biology, where phase separation emerged as a ubiquitous self-organizational principle. △ Less

Submitted 31 August, 2022; v1 submitted 26 August, 2022; originally announced August 2022.

Comments: 9 pages, 6 figures, supplementary information

Dynamics of active liquid interfaces

Authors: Raymond Adkins, Itamar Kolvin, Zhihong You, Sven Witthaus, M. Cristina Marchetti, Zvonimir Dogic

Abstract: Controlling interfaces of phase separating fluid mixtures is key to creating diverse functional soft materials. Traditionally, this is accomplished with surface-modifying chemical agents. Using experiment and theory, we study how mechanical activity shapes soft interfaces that separate an active and a passive fluid. Chaotic flows in the active fluid give rise to giant interfacial fluctuations and… ▽ More Controlling interfaces of phase separating fluid mixtures is key to creating diverse functional soft materials. Traditionally, this is accomplished with surface-modifying chemical agents. Using experiment and theory, we study how mechanical activity shapes soft interfaces that separate an active and a passive fluid. Chaotic flows in the active fluid give rise to giant interfacial fluctuations and non-inertial propagating active waves. At high activities, stresses disrupt interface continuity and drive droplet generation, producing an emulsion-like active state comprised of finite-sized droplets. When in contact with a solid boundary, active interfaces exhibit non-equilibrium wetting transitions, wherein the fluid climbs the wall against gravity. These results demonstrate the promise of mechanically driven interfaces for creating a new class of soft active matter. △ Less

Submitted 11 August, 2022; originally announced August 2022.

Comments: 62 pages, 5 figures, 13 supplementary figures

Journal ref: Science, 11 Aug 2022, Vol 377, Issue 6607 pp. 768-772

Active nematic defects in compressible and incompressible flows

Authors: Supavit Pokawanvit, Zhitao Chen, Zhihong You, Luiza Angheluta, M. Cristina Marchetti, Mark J. Bowick

Abstract: We study the dynamics of active nematic films on a substrate driven by active flows with or without the incompressible constraint.Through simulations and theoretical analysis, we show that arch patterns are stable in the compressible case, whilst they become unstable under the incompressibility constraint. For compressible flows at high enough activity, stable arches organize themselves into a sme… ▽ More We study the dynamics of active nematic films on a substrate driven by active flows with or without the incompressible constraint.Through simulations and theoretical analysis, we show that arch patterns are stable in the compressible case, whilst they become unstable under the incompressibility constraint. For compressible flows at high enough activity, stable arches organize themselves into a smectic-like pattern, which induce an associated global polar ordering of $+1/2$ nematic defects. By contrast, divergence-free flows give rise to a local nematic order of the $+1/2$ defects, consisting of anti-aligned pairs of neighboring defects, as established in previous studies. △ Less

Submitted 28 October, 2022; v1 submitted 27 June, 2022; originally announced June 2022.

Comments: 10 pages, 4 figures

Activity-suppressed phase separation

Authors: Fernando Caballero, M. Cristina Marchetti

Abstract: We use a continuum model to examine the effect of activity on a phase separating mixture of an extensile active nematic and a passive fluid. We highlight the distinct role of previously considered interfacial active stresses and bulk active stresses that couple to liquid crystalline degrees of freedom. Interfacial active stresses can arrest phase separation, as previously demonstrated. Bulk extens… ▽ More We use a continuum model to examine the effect of activity on a phase separating mixture of an extensile active nematic and a passive fluid. We highlight the distinct role of previously considered interfacial active stresses and bulk active stresses that couple to liquid crystalline degrees of freedom. Interfacial active stresses can arrest phase separation, as previously demonstrated. Bulk extensile active stresses can additionally strongly suppress phase separation by sustained self-stirring of the fluid, substantially reducing the size of the coexistence region in the temperature/concentration plane relative to that of the passive system. The phase separated state is a dynamical emulsion-like steady state of continuously splitting and merging droplets, as suggested by recent experiments. Using scaling analysis and simulations, we identify various regimes for the dependence of droplet size on activity. These results can provide a criterion for identifying the mechanisms responsible for arresting phase separation in experiments. △ Less

Submitted 22 December, 2022; v1 submitted 25 June, 2022; originally announced June 2022.

Journal ref: Phys. Rev. Lett. 129, 268002 (2022)

Boundaries control active channel flows

Authors: Paarth Gulati, Suraj Shankar, M. Cristina Marchetti

Abstract: Boundary conditions dictate how fluids, including liquid crystals, flow when pumped through a channel. Can boundary conditions also be used to control internally driven active fluids that generate flows spontaneously? By using numerical simulations and stability analysis we explore how surface anchoring of active agents at the boundaries and substrate drag can be used to rectify coherent flow of a… ▽ More Boundary conditions dictate how fluids, including liquid crystals, flow when pumped through a channel. Can boundary conditions also be used to control internally driven active fluids that generate flows spontaneously? By using numerical simulations and stability analysis we explore how surface anchoring of active agents at the boundaries and substrate drag can be used to rectify coherent flow of an active polar fluid in a 2D channel. Upon increasing activity, a succession of dynamical states is obtained, from laminar flow to vortex arrays to eventual turbulence, that are controlled by the interplay between the hydrodynamic screening length and the extrapolation length quantifying the anchoring strength of the orientational order parameter. We highlight the key role of symmetry in both flow and order and show that coherent laminar flow with net throughput is only possible for weak anchoring and intermediate activity. Our work demonstrates the possibility of controlling the nature and properties of active flows in a channel simply by patterning the confining boundaries. △ Less

Submitted 17 May, 2022; originally announced May 2022.

Comments: 12 pages, 9 figures

Yield Stress and Compliance in Active Cell Monolayers

Authors: Austin Hopkins, Michael Chiang, Benjamin Loewe, Davide Marenduzzo, M. Cristina Marchetti

Abstract: The rheology of biological tissue plays an important role in many processes, from organ formation to cancer invasion. Here, we use a multi-phase field model of motile cells to simulate active microrheology within a tissue monolayer. When unperturbed, the tissue exhibits a transition between a solid-like state and a fluid-like state tuned by cell motility and deformability - the ratio of the energe… ▽ More The rheology of biological tissue plays an important role in many processes, from organ formation to cancer invasion. Here, we use a multi-phase field model of motile cells to simulate active microrheology within a tissue monolayer. When unperturbed, the tissue exhibits a transition between a solid-like state and a fluid-like state tuned by cell motility and deformability - the ratio of the energetic costs of steric cell-cell repulsion and cell surface tension. When perturbed, solid tissues exhibit yield-stress behavior, with a threshold force for the onset of motion of a probe particle that vanishes upon approaching the solid-to-liquid transition. This onset of motion is qualitatively different in the low and high deformability regimes. At high deformability, the tissue is amorphous when solid, it responds compliantly to deformations, and the probe transition to motion is smooth. At low deformability, the monolayer is more ordered translationally and stiffer, and the onset of motion appears discontinuous. Our results suggest that cellular or nanoparticle transport in different types of tissues can be fundamentally different, and point to ways in which it can be controlled. △ Less

Submitted 3 March, 2022; originally announced March 2022.

Comments: Main text (6 pages, 4 figures) plus SM (8 pages, 9 figures)

Flow around topological defects in active nematic films

Authors: Jonas Rønning, M. Cristina Marchetti, Mark J. Bowick, Luiza Angheluta

Abstract: We study the active flow around isolated defects and the self-propulsion velocity of $+1/2$ defects in an active nematic film with both viscous dissipation (with viscosity $η$) and frictional damping $Γ$ with a substrate. The interplay between these two dissipation mechanisms is controlled by the hydrodynamic dissipation length $\ell_d=\sqrt{η/Γ}$ that screens the flows. For an isolated defect, in… ▽ More We study the active flow around isolated defects and the self-propulsion velocity of $+1/2$ defects in an active nematic film with both viscous dissipation (with viscosity $η$) and frictional damping $Γ$ with a substrate. The interplay between these two dissipation mechanisms is controlled by the hydrodynamic dissipation length $\ell_d=\sqrt{η/Γ}$ that screens the flows. For an isolated defect, in the absence of screening from other defects, the size of the vortical flows around the defect is controlled by the system size $R$. In the presence of friction that leads to a finite value of $\ell_d$, the vorticity field decays to zero on the lengthscales larger than $\ell_d$. We show that the self-propulsion velocity of $+1/2$ defects grows with $R$ in small systems where $R<\ell_d$, while in the infinite system limit or when $R\gg \ell_d$, it approaches a constant value determined by $\ell_d$. △ Less

Submitted 27 January, 2022; v1 submitted 16 November, 2021; originally announced November 2021.

Comments: 21 pages, 8 figures

Journal ref: Proc. R. Soc. A. 478.2257 (2022): 20210879

Anomalous elasticity of cellular tissue vertex model

Authors: Arthur Hernandez, Michael F. Staddon, Mark J. Bowick, M. Cristina Marchetti, Michael Moshe

Abstract: Vertex Models, as used to describe cellular tissue, have an energy controlled by deviations of each cell area and perimeter from target values. The constrained nonlinear relation between area and perimeter leads to new mechanical response. Here we provide a mean-field treatment of a highly simplified model: a uniform network of regular polygons with no topological rearrangements. Since all polygon… ▽ More Vertex Models, as used to describe cellular tissue, have an energy controlled by deviations of each cell area and perimeter from target values. The constrained nonlinear relation between area and perimeter leads to new mechanical response. Here we provide a mean-field treatment of a highly simplified model: a uniform network of regular polygons with no topological rearrangements. Since all polygons deform in the same way, we only need to analyze the ground states and the response to deformations of a single polygon (cell). The model exhibits the known transition between a fluid/compatible state, where the cell can accommodate both target area and perimeter, and a rigid/incompatible state. %The rigid solid-like state has a single gapped ground state. We calculate and measure the mechanical resistance to various deformation protocols and discover that at the onset of rigidity, where a single zero-energy ground-state exists, %We show that in the incompatible state, where a single frustrated ground-state exists, linear elasticity fails to describe the mechanical response to even infinitesimal deformations. In particular we identify a breakdown of reciprocity expressed via different moduli for compressive and tensile loads, implying non-analyticity of the energy functional. We give a pictorial representation in configuration space that reveals that the complex elastic response of the Vertex Model arises from the presence of two distinct sets of reference states (associated with target area and target perimeter). △ Less

Submitted 20 April, 2022; v1 submitted 21 September, 2021; originally announced September 2021.

Journal ref: Phys. Rev. E 105, 064611 (2022)

Shear-driven solidification and nonlinear elasticity in epithelial tissues

Authors: Junxiang Huang, James O. Cochran, Suzanne M. Fielding, M. Cristina Marchetti, Dapeng Bi

Abstract: Biological processes, from morphogenesis to tumor invasion, spontaneously generate shear stresses inside living tissue. The mechanisms that govern the transmission of mechanical forces in epithelia and the collective response of the tissue to bulk shear deformations remain, however, poorly understood. Using a minimal cell-based computational model, we investigate the constitutive relation of confl… ▽ More Biological processes, from morphogenesis to tumor invasion, spontaneously generate shear stresses inside living tissue. The mechanisms that govern the transmission of mechanical forces in epithelia and the collective response of the tissue to bulk shear deformations remain, however, poorly understood. Using a minimal cell-based computational model, we investigate the constitutive relation of confluent tissues under simple shear deformation. We show that an initially undeformed fluidlike tissue acquires finite rigidity above a critical applied strain. This is akin to the shear-driven rigidity observed in other soft matter systems. Interestingly, shear-driven rigidity can be understood by a critical scaling analysis in the vicinity of the second order critical point that governs the liquid-solid transition of the undeformed system. We further show that a solidlike tissue responds linearly only to small strains and but then switches to a nonlinear response at larger stains, with substantial stiffening. Finally, we propose a mean-field formulation for cells under shear that offers a simple physical explanation of shear-driven rigidity and nonlinear response in a tissue. △ Less

Submitted 27 November, 2022; v1 submitted 21 September, 2021; originally announced September 2021.

Journal ref: Phys. Rev. Lett. 128, 178001 (2022)

Symmetry, Thermodynamics and Topology in Active Matter

Authors: Mark J. Bowick, Nikta Fakhri, M. Cristina Marchetti, Sriram Ramaswamy

Abstract: This article summarizes some of the open questions in the field of active matter that have emerged during Active20, a nine-week program held at the Kavli Institute for Theoretical Physics (KITP) in Spring 2020. The article does not provide a review of the field, but rather a personal view of the authors, informed by contributions of all participants, on new directions in active matter research. Th… ▽ More This article summarizes some of the open questions in the field of active matter that have emerged during Active20, a nine-week program held at the Kavli Institute for Theoretical Physics (KITP) in Spring 2020. The article does not provide a review of the field, but rather a personal view of the authors, informed by contributions of all participants, on new directions in active matter research. The topics highlighted include: the ubiquitous occurrence of spontaneous flows and active turbulence and the theoretical and experimental challenges associated with controlling and harnessing such flows; the role of motile topological defects in ordered states of active matter and their possible biological relevance; the emergence of non-reciprocal effective interactions and the role of chirality in active systems and their intriguing connections to non-Hermitian quantum mechanics; the progress towards a formulation of the thermodynamics of active systems thanks to the feedback between theory and experiments; the impact of the active matter framework on our understanding of the emergent mechanics of biological tissue. These seemingly diverse phenomena all stem from the defining property of active matter - assemblies of self-driven entities that individually break time-reversal symmetry and collectively organize in a rich variety of nonequilibrium states. △ Less

Submitted 7 August, 2021; v1 submitted 1 July, 2021; originally announced July 2021.

Comments: 23 pages, 10 figures

Journal ref: Phys. Rev. X12, 010501 (2022)

Alignment interactions drive structural transitions in biological tissues

Authors: Matteo Paoluzzi, Luca Angelani, Giorgio Gosti, M Cristina Marchetti, Ignacio Pagonabarraga, Giancarlo Ruocco

Abstract: Experimental evidence shows that there is a feedback between cell shape and cell motion. How this feedback impacts the collective behavior of dense cell monolayers remains an open question. We investigate the effect of a feedback that tends to align the cell crawling direction with cell elongation in a biological tissue model. We find that the alignment interaction promotes nematic patterns in the… ▽ More Experimental evidence shows that there is a feedback between cell shape and cell motion. How this feedback impacts the collective behavior of dense cell monolayers remains an open question. We investigate the effect of a feedback that tends to align the cell crawling direction with cell elongation in a biological tissue model. We find that the alignment interaction promotes nematic patterns in the fluid phase that eventually undergo a non-equilibrium phase transition into a quasi-hexagonal solid. Meanwhile, highly asymmetric cells do not undergo the liquid-to-solid transition for any value of the alignment coupling. In this regime, the dynamics of cell centers and shape fluctuation show features typical of glassy systems. △ Less

Submitted 23 September, 2021; v1 submitted 1 July, 2021; originally announced July 2021.

The role of fluid flow in the dynamics of active nematic defects

Authors: Luiza Angheluta, Zhitao Chen, M. Cristina Marchetti, Mark J. Bowick

Abstract: We adapt the Halperin-Mazenko formalism to analyze two-dimensional active nematics coupled to a generic fluid flow. The governing hydrodynamic equations lead to evolution laws for nematic topological defects and their corresponding density fields. We find that $\pm 1/2$ defects are propelled by the local fluid flow and by the nematic orientation coupled with the flow shear rate. In the overdamped… ▽ More We adapt the Halperin-Mazenko formalism to analyze two-dimensional active nematics coupled to a generic fluid flow. The governing hydrodynamic equations lead to evolution laws for nematic topological defects and their corresponding density fields. We find that $\pm 1/2$ defects are propelled by the local fluid flow and by the nematic orientation coupled with the flow shear rate. In the overdamped and compressible limit, we recover the previously obtained active self-propulsion of the +1/2 defects. Non-local hydrodynamic effects are primarily significant for incompressible flows, for which it is not possible to eliminate the fluid velocity in favor of the local defect polarization alone. For the case of two defects with opposite charge, the non-local hydrodynamic interaction is mediated by non-reciprocal pressure-gradient forces. Finally, we derive continuum equations for a defect gas coupled to an arbitrary (compressible or incompressible) fluid flow. △ Less

Submitted 5 December, 2020; originally announced December 2020.

Journal ref: New Journal of Physics 23.3 (2021): 033009

Poisson-bracket formulation of the dynamics of fluids of deformable particles

Authors: Arthur Hernandez, M. Cristina Marchetti

Abstract: Using the Poisson bracket method, we derive continuum equations for a fluid of deformable particles in two dimensions. Particle shape is quantified in terms of two continuum fields: an anisotropy density field that captures the deformations of individual particles from regular shapes and a shape tensor density field that quantifies both particle elongation and nematic alignment of elongated shapes… ▽ More Using the Poisson bracket method, we derive continuum equations for a fluid of deformable particles in two dimensions. Particle shape is quantified in terms of two continuum fields: an anisotropy density field that captures the deformations of individual particles from regular shapes and a shape tensor density field that quantifies both particle elongation and nematic alignment of elongated shapes. We explicitly consider the example of a dense biological tissue as described by the Vertex model energy, where cell shape has been proposed as a structural order parameter for a liquid-solid transition. The hydrodynamic model of biological tissue proposed here captures the coupling of cell shape to flow, and provides a starting point for modeling the rheology of dense tissue. △ Less

Submitted 22 November, 2020; originally announced November 2020.

Comments: 8 pages, 1 figure

Journal ref: Phys. Rev. E 103, 032612 (2021)

Fluctuations can induce local nematic order and extensile stress in monolayers of motile cells

Authors: Farzan Vafa, Mark J. Bowick, Boris I. Shraiman, M. Cristina Marchetti

Abstract: Recent experiments in various cell types have shown that two-dimensional tissues often display local nematic order, with evidence of extensile stresses manifest in the dynamics of topological defects. Using a mesoscopic model where tissue flow is generated by fluctuating traction forces coupled to the nematic order parameter, we show that the resulting tissue dynamics can spontaneously produce loc… ▽ More Recent experiments in various cell types have shown that two-dimensional tissues often display local nematic order, with evidence of extensile stresses manifest in the dynamics of topological defects. Using a mesoscopic model where tissue flow is generated by fluctuating traction forces coupled to the nematic order parameter, we show that the resulting tissue dynamics can spontaneously produce local nematic order and an extensile internal stress. A key element of the model is the assumption that in the presence of local nematic alignment, cells preferentially crawl along the nematic axis, resulting in anisotropy of fluctuations. Our work shows that activity can drive either extensile or contractile stresses in tissue, depending on the relative strength of the contractility of the cortical cytoskeleton and tractions by cells on the extracellular matrix. △ Less

Submitted 25 February, 2021; v1 submitted 3 November, 2020; originally announced November 2020.

Comments: 6 pages, 1 figure. v2: minor edits and clarifying comments

Journal ref: Soft Matter, 2021, 17, 3068-3073

Topological active matter

Authors: Suraj Shankar, Anton Souslov, Mark J. Bowick, M. Cristina Marchetti, Vincenzo Vitelli

Abstract: Active matter encompasses different nonequilibrium systems in which individual constituents convert energy into non-conservative forces or motion at the microscale. This review provides an elementary introduction to the role of topology in active matter through experimentally relevant examples. Here, the focus lies on topological defects and topologically protected edge modes with an emphasis on t… ▽ More Active matter encompasses different nonequilibrium systems in which individual constituents convert energy into non-conservative forces or motion at the microscale. This review provides an elementary introduction to the role of topology in active matter through experimentally relevant examples. Here, the focus lies on topological defects and topologically protected edge modes with an emphasis on the distinctive properties they acquire in active media. These paradigmatic examples represent two physically distinct classes of phenomena whose robustness can be traced to a common mathematical origin: the presence of topological invariants. These invariants are typically integer numbers that cannot be changed by continuous deformations of the relevant order parameters or physical parameters of the underlying medium. We first explain the mechanisms whereby topological defects self propel and proliferate in active nematics, leading to collective states which can be manipulated by geometry and patterning. Possible implications for active microfluidics and biological tissues are presented. We then illustrate how the propagation of waves in active fluids and solids is affected by the presence of topological invariants characterizing their dispersion relations. We discuss the relevance of these ideas for the design of robotic metamaterials and the properties of active granular and colloidal systems. Open theoretical and experimental challenges are presented as future research prospects. △ Less

Submitted 19 July, 2021; v1 submitted 1 October, 2020; originally announced October 2020.

Comments: 23 pages, 3 figures, 4 boxes

Journal ref: Nat Rev Phys 4, 380-398 (2022)

Viscoelastic control of spatiotemporal order in bacterial active matter

Authors: Song Liu, Suraj Shankar, M. Cristina Marchetti, Yilin Wu

Abstract: Active matter consists of units that generate mechanical work by consuming energy. Examples include living systems, such as assemblies of bacteria and biological tissues, biopolymers driven by molecular motors, and suspensions of synthetic self-propelled particles. A central question in the field is to understand and control the self-organization of active assemblies in space and time. Most active… ▽ More Active matter consists of units that generate mechanical work by consuming energy. Examples include living systems, such as assemblies of bacteria and biological tissues, biopolymers driven by molecular motors, and suspensions of synthetic self-propelled particles. A central question in the field is to understand and control the self-organization of active assemblies in space and time. Most active systems exhibit either spatial order mediated by interactions that coordinate the spatial structure and the motion of active agents or the temporal synchronization of individual oscillatory dynamics. The simultaneous control of spatial and temporal organization is more challenging and generally requires complex interactions, such as reaction-diffusion hierarchies or genetically engineered cellular circuits. Here, we report a novel and simple means to simultaneously control the spatial and temporal self-organization of bacterial active matter. By confining an active bacterial suspension and manipulating a single macroscopic parameter, namely the viscoelasticity of the suspending fluid, we have found that the bacterial fluid first self-organizes in space into a millimeter-scale rotating vortex; then displays temporal organization as the giant vortex switches its global chirality periodically with tunable frequency, reminiscent of a torsional pendulum - a self-driven one. Combining experiments with an active matter model, we explain this striking behavior in terms of the interplay between active forcing and viscoelastic stress relaxation. Our findings advance the understanding of bacterial behavior in complex fluids, and demonstrate experimentally for the first time that rheological properties can be harnessed to control active matter flows. Coupled with actuation, our tunable self-oscillating bacterial vortex may be used as a "clock" for locomotion of soft robots and microfluidic pumping. △ Less

Submitted 13 November, 2020; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 29 pages with 4 figure including Methods and Extended Data (with 10 figures); SI available upon request (accepted version)

Journal ref: Nature 590, 80-84 (2021)

Multi-defect Dynamics in Active Nematics

Authors: Farzan Vafa, Mark J. Bowick, M. Cristina Marchetti, Boris I. Shraiman

Abstract: Recent experiments and numerical studies have drawn attention to the dynamics of active nematics. Two-dimensional active nematics flow spontaneously and exhibit spatiotemporal chaotic flows with proliferation of topological defects in the nematic texture. It has been proposed that the dynamics of active nematics can be understood in terms of the dynamics of interacting defects, propelled by active… ▽ More Recent experiments and numerical studies have drawn attention to the dynamics of active nematics. Two-dimensional active nematics flow spontaneously and exhibit spatiotemporal chaotic flows with proliferation of topological defects in the nematic texture. It has been proposed that the dynamics of active nematics can be understood in terms of the dynamics of interacting defects, propelled by active stress. Previous work has derived effective equations of motion for individual defects as quasi-particles moving in the mean field generated by other defects, but an effective theory governing multi-defect dynamics has remained out of reach. In this paper, we examine the dynamics of 2D active nematics in the limit of strong order and overdamped compressible flow. The activity-induced defect dynamics is formulated as a perturbation of the manifold of quasi-static nematic textures explicitly parameterized by defect positions. This makes it possible to derive a set of coupled ordinary differential equations governing defect (and therefore texture) dynamics. Interestingly, because of the non-orthogonality of textures associated with individual defects, their motion is coupled through a position dependent ``collective mobility" matrix. In addition to the familiar active self-propulsion of the $+1/2$ defect, we obtain new collective effects of activity that can be interpreted in terms of non-central and non-reciprocal interactions between defects. △ Less

Submitted 6 July, 2020; originally announced July 2020.

Comments: 14 pages, 3 figures

How does the extracellular matrix affect the rigidity of an embedded spheroid?

Authors: Amanda Parker, M. Cristina Marchetti, M. Lisa Manning, J. M. Schwarz

Abstract: Cellularized tissue and polymer networks can both transition from floppy to rigid as a function of their control parameters, and, yet, the two systems often mechanically interact, which may affect their respective rigidities. To study this interaction, we consider a vertex model with interfacial tension (a spheroid) embedded in a spring network in two dimensions. We identify two regimes with diffe… ▽ More Cellularized tissue and polymer networks can both transition from floppy to rigid as a function of their control parameters, and, yet, the two systems often mechanically interact, which may affect their respective rigidities. To study this interaction, we consider a vertex model with interfacial tension (a spheroid) embedded in a spring network in two dimensions. We identify two regimes with different global spheroid shapes, governed by the pressure resulting from competition between interfacial tension and tension in the network. In the first regime, the tissue remains compact, while in the second, a cavitation-like instability leads to the emergence of gaps at the tissue-network interface. Intriguingly, compression of the tissue promotes fluidization, while tension promotes cellular alignment and rigidification, with the mechanisms driving rigidification differing on either side of the instability. △ Less

Submitted 29 June, 2020; originally announced June 2020.

Comments: Main text: 5 pages, 3 figures. Supplementary Material: 13 pages, 14 figures

Nonreciprocity as a generic route to traveling states

Authors: Zhihong You, Aparna Baskaran, M. Cristina Marchetti

Abstract: We examine a non-reciprocally coupled dynamical model of a mixture of two diffusing species. We demonstrate that nonreciprocity, which is encoded in the model via antagonistic cross diffusivities, provides a generic mechanism for the emergence of traveling patterns in purely diffusive systems with conservative dynamics. In the absence of non-reciprocity, the binary fluid mixture undergoes a phase… ▽ More We examine a non-reciprocally coupled dynamical model of a mixture of two diffusing species. We demonstrate that nonreciprocity, which is encoded in the model via antagonistic cross diffusivities, provides a generic mechanism for the emergence of traveling patterns in purely diffusive systems with conservative dynamics. In the absence of non-reciprocity, the binary fluid mixture undergoes a phase transition from a homogeneous mixed state to a demixed state with spatially separated regions rich in one of the two components. Above a critical value of the parameter tuning non-reciprocity, the static demixed pattern acquires a finite velocity, resulting in a state that breaks both spatial and time translational symmetry, as well as the reflection parity of the static pattern. We elucidate the generic nature of the transition to traveling patterns using a minimal model that can be studied analytically. Our work has direct relevance to nonequilibrium assembly in mixtures of chemically interacting colloids that are known to exhibit non-reciprocal effective interactions, as well as to mixtures of active and passive agents where traveling states of the type predicted here have been observed in simulations. It also provides insight on transitions to traveling and oscillatory states seen in a broad range of nonreciprocal systems with non-conservative dynamics, from reaction-diffusion and prey-predators models to multispecies mixtures of microorganisms with antagonistic interactions. △ Less

Submitted 20 August, 2020; v1 submitted 15 May, 2020; originally announced May 2020.

Comments: 8 pages, 3 figures

Journal ref: Proc. Natl. Acad. Sci. U.S.A. 117(33), 19767 (2020)

Emergence of dynamic vortex glasses in disordered polar active fluids

Authors: Amélie Chardac, Suraj Shankar, M. Cristina Marchetti, Denis Bartolo

Abstract: In equilibrium, disorder conspires with topological defects to redefine the ordered states of matter in systems as diverse as crystals, superconductors and liquid crystals. Far from equilibrium, however, the consequences of quenched disorder on active condensed matter remain virtually uncharted. Here, we reveal a state of strongly disordered active matter with no counterparts in equilibrium: a dyn… ▽ More In equilibrium, disorder conspires with topological defects to redefine the ordered states of matter in systems as diverse as crystals, superconductors and liquid crystals. Far from equilibrium, however, the consequences of quenched disorder on active condensed matter remain virtually uncharted. Here, we reveal a state of strongly disordered active matter with no counterparts in equilibrium: a dynamical vortex glass. Combining high-content microfluidic experiments and theory, we show how colloidal flocks collectively cruise through disorder without relaxing the topological singularities of their flows. The resulting state is highly dynamical but the flow patterns, shaped by a finite density of frozen vortices, are stationary and exponentially degenerated. Quenched isotropic disorder acts as a random gauge field turning active liquids into dynamical vortex glasses. We argue that this robust mechanism should shape the collective dynamics of a broad class of disordered active matter, from synthetic active nematics to collections of living cells exploring heterogeneous media. △ Less

Submitted 7 September, 2020; v1 submitted 28 February, 2020; originally announced February 2020.

Comments: 9 pages, 3 figures (corrected a Grant number)

Journal ref: PNAS March 9, 2021 118 (10) e2018218118

Shape and size changes of adherent elastic epithelia

Authors: Benjamin Loewe, Francesco Serafin, Suraj Shankar, Mark J. Bowick, M. Cristina Marchetti

Abstract: Epithelial tissues play a fundamental role in various morphogenetic events during development and early embryogenesis. Although epithelial monolayers are often modeled as two-dimensional (2D) elastic surfaces, they distinguish themselves from conventional thin elastic plates in three important ways - the presence of an apical-basal polarity, spatial variability of cellular thickness, and their non… ▽ More Epithelial tissues play a fundamental role in various morphogenetic events during development and early embryogenesis. Although epithelial monolayers are often modeled as two-dimensional (2D) elastic surfaces, they distinguish themselves from conventional thin elastic plates in three important ways - the presence of an apical-basal polarity, spatial variability of cellular thickness, and their nonequilibrium active nature. Here, we develop a minimal continuum model of a planar epithelial tissue as an active elastic material that incorporates all these features. We start from a full three-dimensional (3D)description of the tissue and derive an effective 2D model that captures, through the curvature of the apical surface, both the apical-basal asymmetry and the spatial geometry of the tissue. Crucially, variations of active stresses across the apical-basal axis lead to active torques that can drive curvature transitions. By identifying four distinct sources of activity, we find that bulk active stresses arising from actomyosin contractility and growth compete with boundary active tensions due to localized actomyosin cables and lamellipodial activity to generate the various states spanning the morphospace of a planar epithelium. Our treatment hence unifies 3D shape deformations through the coupled mechanics of apical curvature change and in-plane expansion/contraction of substrate-adhered tissues. Finally, we discuss the implications of our results for some biologically relevant processes such as tissue folding at the onset of lumen formation. △ Less

Submitted 11 June, 2020; v1 submitted 10 February, 2020; originally announced February 2020.

Comments: 12 pages, 8 figures (including Appendix), published version including more calculations without changing results

Journal ref: Soft Matter, 16, 5282-5293, 2020

Information and motility exchange in collectives of active particles

Authors: Matteo Paoluzzi, Marco Leoni, M Cristina Marchetti

Abstract: We examine the interplay of motility and information exchange in a model of run-and-tumble active particles where the particle's motility is encoded as a bit of information that can be exchanged upon contact according to the rules of AND and OR logic gates in a circuit. Motile AND particles become non-motile upon contact with a non-motile particle. Conversely, motile OR particles remain motile upo… ▽ More We examine the interplay of motility and information exchange in a model of run-and-tumble active particles where the particle's motility is encoded as a bit of information that can be exchanged upon contact according to the rules of AND and OR logic gates in a circuit. Motile AND particles become non-motile upon contact with a non-motile particle. Conversely, motile OR particles remain motile upon collision with their non-motile counterparts. AND particles that have become non-motile additionally "reawaken", i.e., recover their motility, at a fixed rate $μ$, as in the SIS (Susceptible, Infected, Susceptible) model of epidemic spreading, where an infected agent can become healthy again, but keeps no memory of the recent infection, hence it is susceptible to a renewed infection. For $μ=0$, both AND and OR particles relax irreversibly to absorbing states of all non-motile or all motile particles, respectively. The relaxation kinetics is, however, faster for OR particles that remain active throughout the process. At finite $μ$, the AND dynamics is controlled by the interplay between reawakening and collision rates. The system evolves to a state of all motile particles (an absorbing state in the language of absorbing phase transitions) for $μ>μ_c$ and to a mixed state with coexisting motile and non-motile particles (an active state in the language of absorbing phase transitions) for $μ<μ_c$. The final state exhibits a rich structure controlled by motility-induced aggregation. Our work can be relevant to biochemical signaling in motile bacteria, the spreading of epidemics and of social consensus, as well as light-controlled organization of active colloids. △ Less

Submitted 24 June, 2020; v1 submitted 4 February, 2020; originally announced February 2020.

Journal ref: Soft Matter, 2020, Advance Article

Solid-Liquid Transition of Deformable and Overlapping Active Particles

Authors: Benjamin Loewe, Michael Chiang, Davide Marenduzzo, M. Cristina Marchetti

Abstract: 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 t… ▽ More 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. △ Less

Submitted 19 May, 2020; v1 submitted 22 December, 2019; originally announced December 2019.

Comments: Figures have been updated and the text expanded. Main text (6 pages, 4 figures) plus SM (12 pages, 16 figures). Supplementary movies can be found at https://drive.google.com/open?id=1yfMf0aFqFlgjsGGjE0Zr3zo_eiIcREB9

Journal ref: Phys. Rev. Lett. 125, 038003 (2020)

Hydrodynamics of Active Defects: from order to chaos to defect ordering

Authors: Suraj Shankar, M. Cristina Marchetti

Abstract: Topological defects play a prominent role in the physics of two-dimensional materials. When driven out of equilibrium in active nematics, disclinations can acquire spontaneous self-propulsion and drive self-sustained flows upon proliferation. Here we construct a general hydrodynamic theory for a two-dimensional active nematic interrupted by a large number of such defects. Our equations describe th… ▽ More Topological defects play a prominent role in the physics of two-dimensional materials. When driven out of equilibrium in active nematics, disclinations can acquire spontaneous self-propulsion and drive self-sustained flows upon proliferation. Here we construct a general hydrodynamic theory for a two-dimensional active nematic interrupted by a large number of such defects. Our equations describe the flows and spatio-temporal defect chaos characterizing active turbulence, even close to the defect unbinding transition. At high activity, nonequilibrium torques combined with many-body screening cause the active disclinations to spontaneously break rotational symmetry forming a collectively moving defect ordered polar liquid. By recognizing defects as the relevant quasiparticle excitations, we construct a comprehensive phase diagram for two-dimensional active nematics. Using our hydrodynamic approach, we additionally show that activity gradients can act like "electric fields", driving the sorting of topological charge. This demonstrates the versatility of our continuum model and its relevance for quantifying the use of spatially inhomogeneous activity for controlling active flows and for the fabrication of active devices with targeted transport capabilities. △ Less

Submitted 23 October, 2019; v1 submitted 4 July, 2019; originally announced July 2019.

Comments: 18 pages, 7 figures, additional explanation provided with results unchanged

Journal ref: Phys. Rev. X 9, 041047 (2019)

Condensate formation and multiscale dynamics in two-dimensional active suspensions

Authors: Moritz Linkmann, M. Cristina Marchetti, Guido Boffetta, Bruno Eckhardt

Abstract: The collective effects of microswimmers in active suspensions result in active turbulence, a spatiotemporally chaotic dynamics at mesoscale, which is characterized by the presence of vortices and jets at scales much larger than the characteristic size of the individual active constituents. To describe this dynamics, Navier-Stokes-based one-fluid models driven by small-scale forces have been propos… ▽ More The collective effects of microswimmers in active suspensions result in active turbulence, a spatiotemporally chaotic dynamics at mesoscale, which is characterized by the presence of vortices and jets at scales much larger than the characteristic size of the individual active constituents. To describe this dynamics, Navier-Stokes-based one-fluid models driven by small-scale forces have been proposed. Here, we provide a justification of such models for the case of dense suspensions in two dimensions (2d). We subsequently carry out an in-depth numerical study of the properties of one-fluid models as a function of the active driving in view of possible transition scenarios from active turbulence to large-scale pattern, referred to as condensate, formation induced by the classical inverse energy cascade in Newtonian 2d turbulence. Using a one-fluid model it was recently shown (Linkmann et al., Phys. Rev. Lett. (in press)) that two-dimensional active suspensions support two non-equilibrium steady states, one with a condensate and one without, which are separated by a subcritical transition. Here, we report further details on this transition such as hysteresis and discuss a low-dimensional model that describes the main features of the transition through nonlocal-in-scale coupling between the small-scale driving and the condensate. △ Less

Submitted 15 May, 2019; originally announced May 2019.

Journal ref: Phys. Rev. E 101, 022609 (2020)

Glassy Dynamics in Models of Confluent Tissue with Mitosis and Apoptosis

Authors: Michael Czajkowski, Daniel M. Sussman, M. Cristina Marchetti, M. Lisa Manning

Abstract: Recent work on particle-based models of tissues has suggested that any finite rate of cell division and cell death is sufficient to fluidize an epithelial tissue. At the same time, experimental evidence has indicated the existence of glassy dynamics in some epithelial layers despite continued cell cycling. To address this discrepancy, we quantify the role of cell birth and death on glassy states i… ▽ More Recent work on particle-based models of tissues has suggested that any finite rate of cell division and cell death is sufficient to fluidize an epithelial tissue. At the same time, experimental evidence has indicated the existence of glassy dynamics in some epithelial layers despite continued cell cycling. To address this discrepancy, we quantify the role of cell birth and death on glassy states in confluent tissues using simulations of an active vertex model that includes cell motility, cell division, and cell death. Our simulation data is consistent with a simple ansatz in which the rate of cell-life cycling and the rate of relaxation of the tissue in the absence of cell cycling contribute independently and additively to the overall rate of cell motion. Specifically, we find that a glass-like regime with caging behavior indicated by subdiffusive cell displacements can be achieved in systems with sufficiently low rates of cell cycling. △ Less

Submitted 5 May, 2019; originally announced May 2019.

Comments: Main text (10 pages, 5 figures) and Appendices (7 pages, 7 figures)

Small-scale demixing in confluent biological tissues

Authors: Preeti Sahu, Daniel M. Sussman, Matthias Rubsam, Aaron F. Mertz, Valerie Horsley, Eric R. Dufresne, Carien M. Niessen, M. Cristina Marchetti, M. Lisa Manning, J. M. Schwarz

Abstract: Surface tension governed by differential adhesion can drive fluid particle mixtures to sort into separate regions, i.e., demix. Does the same phenomenon occur in confluent biological tissues? We begin to answer this question for epithelial monolayers with a combination of theory via a vertex model and experiments on keratinocyte monolayers. Vertex models are distinct from particle models in that t… ▽ More Surface tension governed by differential adhesion can drive fluid particle mixtures to sort into separate regions, i.e., demix. Does the same phenomenon occur in confluent biological tissues? We begin to answer this question for epithelial monolayers with a combination of theory via a vertex model and experiments on keratinocyte monolayers. Vertex models are distinct from particle models in that the interactions between the cells are shape-based, as opposed to distance-dependent. We investigate whether a disparity in cell shape or size alone is sufficient to drive demixing in bidisperse vertex model fluid mixtures. Surprisingly, we observe that both types of bidisperse systems robustly mix on large lengthscales. On the other hand, shape disparity generates slight demixing over a few cell diameters, a phenomenon we term micro-demixing. This result can be understood by examining the differential energy barriers for neighbor exchanges (T1 transitions). Experiments with mixtures of wild-type and E-cadherin-deficient keratinocytes on a substrate are consistent with the predicted phenomenon of micro-demixing, which biology may exploit to create subtle patterning. The robustness of mixing at large scales, however, suggests that despite some differences in cell shape and size, progenitor cells can readily mix throughout a developing tissue until acquiring means of recognizing cells of different types. △ Less

Submitted 20 February, 2020; v1 submitted 2 May, 2019; originally announced May 2019.

Comments: 32 pages, 17 figures

Swimmer suspensions on substrates: anomalous stability and long-range order

Authors: Ananyo Maitra, Pragya Srivastava, M. Cristina Marchetti, Sriram Ramaswamy, Martin Lenz

Abstract: We present a comprehensive theory of the dynamics and fluctuations of a two-dimensional suspension of polar active particles in an incompressible fluid confined to a substrate. We show that, depending on the sign of a single parameter, a state with polar orientational order is anomalously stable (or anomalously unstable), with a nonzero relaxation (or growth) rate for angular fluctuations at zero… ▽ More We present a comprehensive theory of the dynamics and fluctuations of a two-dimensional suspension of polar active particles in an incompressible fluid confined to a substrate. We show that, depending on the sign of a single parameter, a state with polar orientational order is anomalously stable (or anomalously unstable), with a nonzero relaxation (or growth) rate for angular fluctuations at zero wavenumber. This screening of the broken-symmetry mode in the stable state does lead to conventional rather than giant number fluctuations as argued by Bricard et al., Nature ${\bf 503}$, 95 (2013), but their bend instability in a splay-stable flock does not exist and the polar phase has long-range order in two dimensions. Our theory also describes confined three-dimensional thin-film suspensions of active polar particles as well as dense compressible active polar rods, and predicts a flocking transition without a banding instability △ Less

Submitted 4 January, 2019; originally announced January 2019.

Journal ref: Phys. Rev. Lett. 124, 028002 (2020)

Fractal aggregation of active particles

Authors: Matteo Paoluzzi, Marco Leoni, M Cristina Marchetti

Abstract: We study active run-and-tumble particles with an additional two-state internal variable characterizing their motile or non-motile state. Motile particles change irreversibly into non-motile ones upon collision with a non-motile particle. The system evolves towards an absorbing state where all particles are non-motile. We initialize the system with one non-motile particles in a bath of motile ones… ▽ More We study active run-and-tumble particles with an additional two-state internal variable characterizing their motile or non-motile state. Motile particles change irreversibly into non-motile ones upon collision with a non-motile particle. The system evolves towards an absorbing state where all particles are non-motile. We initialize the system with one non-motile particles in a bath of motile ones and study numerically the kinetics of relaxation to absorbing state and its structure as function of the density of the initial bath of motile particles and of their tumbling rate. We find a crossover from fractal aggregates at low density to homogeneous ones at high density. The persistence of single-particle dynamics as quantified by the tumbling rate pushes this crossover to higher density and can be used to tune the porosity of the aggregate. At the lowest density the fractal dimension of the aggregate approaches that obtained in single-particle diffusion limited aggregation. Our results could be exploited for the design of structures of desired porosity. The model is a first step towards the study of the collective dynamics of active particles that can exchange biological information. △ Less

Submitted 24 November, 2018; v1 submitted 8 July, 2018; originally announced July 2018.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. E 98, 052603 (2018)

Phase transition to large scale coherent structures in 2d active matter turbulence

Authors: Moritz Linkmann, Guido Boffetta, M. Cristina Marchetti, Bruno Eckhardt

Abstract: The collective motion of microswimmers in suspensions induce patterns of vortices on scales that are much larger than the characteristic size of a microswimmer, attaining a state called bacterial turbulence. Hydrodynamic turbulence acts on even larger scales and is dominated by inertial transport of energy. Using an established modification of the Navier-Stokes equation that accounts for the small… ▽ More The collective motion of microswimmers in suspensions induce patterns of vortices on scales that are much larger than the characteristic size of a microswimmer, attaining a state called bacterial turbulence. Hydrodynamic turbulence acts on even larger scales and is dominated by inertial transport of energy. Using an established modification of the Navier-Stokes equation that accounts for the small scale forcing of hydrodynamic flow by microswimmers, we study the properties of a dense supensions of microswimmers in two dimensions, where the conservation of enstrophy can drive an inverse cascade through which energy is accumulated on the largest scales. We find that the dynamical and statistical properties of the flow show a sharp transition to the formation of vortices at the largest length scale. The results show that 2d bacterial and hydrodynamic turbulence are separated by a subcritical phase transition. △ Less

Submitted 1 August, 2019; v1 submitted 23 June, 2018; originally announced June 2018.

Comments: postprint version

Journal ref: Phys. Rev. Lett. 122, 214503 (2019)

Continuum models of collective cell migration

Authors: Shiladitya Banerjee, M. Cristina Marchetti

Abstract: Collective cell migration plays a central role in tissue development, morphogenesis, wound repair and cancer progression. With the growing realization that physical forces mediate cell motility in development and physiology, a key biological question is how cells integrate molecular activities for force generation on multicellular scales. In this review we discuss recent advances in modeling colle… ▽ More Collective cell migration plays a central role in tissue development, morphogenesis, wound repair and cancer progression. With the growing realization that physical forces mediate cell motility in development and physiology, a key biological question is how cells integrate molecular activities for force generation on multicellular scales. In this review we discuss recent advances in modeling collective cell migration using quantitative tools and approaches rooted in soft matter physics. We focus on theoretical models of cell aggregates as continuous active media, where the feedback between mechanical forces and regulatory biochemistry gives rise to rich collective dynamical behavior. This class of models provides a powerful predictive framework for the physiological dynamics that underlies many developmental processes, where cells need to collectively migrate like a viscous fluid to reach a target region, and then stiffen to support mechanical stresses and maintain tissue cohesion. △ Less

Submitted 5 December, 2018; v1 submitted 16 May, 2018; originally announced May 2018.

Comments: Review article; Figures reproduced with permission from relevant author/publisher

Curvature-dependent tension and tangential flows at the interface of motility-induced phases

Authors: Adam Patch, Daniel Sussman, David Yllanes, M. Cristina Marchetti

Abstract: Purely repulsive active particles spontaneously undergo motility-induced phase separation (MIPS) into condensed and dilute phases. Remarkably, the mechanical tension measured along the interface between these phases is negative. In equilibrium this would imply an unstable interface that wants to expand, but these out-of-equilibrium systems display long-time stability and have intrinsically stiff b… ▽ More Purely repulsive active particles spontaneously undergo motility-induced phase separation (MIPS) into condensed and dilute phases. Remarkably, the mechanical tension measured along the interface between these phases is negative. In equilibrium this would imply an unstable interface that wants to expand, but these out-of-equilibrium systems display long-time stability and have intrinsically stiff boundaries. Here, we study this phenomenon in detail using active Brownian particle simulations and a novel frame of reference. By shifting from the global (or laboratory) frame to a local frame that follows the dynamics of the phase boundary, we observe correlations between the local curvature of the interface and the measured value of the tension. Importantly, our analysis reveals the presence of sustained local tangential motion of particles within a surface layer in both the gas and the dense regions. The combined tangential current in the gas and self-shearing of the surface of the dense phase suggest a stiffening interface that redirects particles along itself to heal local fluctuations. These currents restore the otherwise wildly fluctuating interface through an out-of-equilibrium Marangoni effect. We discuss the implications of our observations on phenomenological models of interfacial dynamics. △ Less

Submitted 18 April, 2018; originally announced April 2018.

Comments: 10 pages, 9 figures

Journal ref: Soft Matter 14 (2018) 7435-7445

Defect unbinding in active nematics

Authors: Suraj Shankar, Sriram Ramaswamy, M. Cristina Marchetti, Mark J. Bowick

Abstract: We formulate the statistical dynamics of topological defects in the active nematic phase, formed in two dimensions by a collection of self-driven particles on a substrate. An important consequence of the non-equilibrium drive is the spontaneous motility of strength +1/2 disclinations. Starting from the hydrodynamic equations of active nematics, we derive an interacting particle description of defe… ▽ More We formulate the statistical dynamics of topological defects in the active nematic phase, formed in two dimensions by a collection of self-driven particles on a substrate. An important consequence of the non-equilibrium drive is the spontaneous motility of strength +1/2 disclinations. Starting from the hydrodynamic equations of active nematics, we derive an interacting particle description of defects that includes active torques. We show that activity, within perturbation theory, lowers the defect-unbinding transition temperature, determining a critical line in the temperature-activity plane that separates the quasi-long-range ordered (nematic) and disordered (isotropic) phases. Below a critical activity, defects remain bound as rotational noise decorrelates the directed dynamics of +1/2 defects, stabilizing the quasi-long-range ordered nematic state. This activity threshold vanishes at low temperature, leading to a re-entrant transition. At large enough activity, active forces always exceed thermal ones and the perturbative result fails, suggesting that in this regime activity will always disorder the system. Crucially, rotational diffusion being a two-dimensional phenomenon, defect unbinding cannot be described by a simplified one-dimensional model. △ Less

Submitted 7 September, 2018; v1 submitted 17 April, 2018; originally announced April 2018.

Comments: 15 pages (including SI), 4 figures. Significant technical improvements without changing the results

Journal ref: Phys. Rev. Lett. 121, 108002 (2018)

Hidden entropy production and work fluctuations in an ideal active gas

Authors: Suraj Shankar, M. Cristina Marchetti

Abstract: Collections of self-propelled particles that move persistently by continuously consuming free energy are a paradigmatic example of active matter. In these systems, unlike Brownian "hot colloids", the breakdown of detailed balance yields a continuous production of entropy at steady state, even for an ideal active gas. We quantify the irreversibility for a non-interacting active particle in two dime… ▽ More Collections of self-propelled particles that move persistently by continuously consuming free energy are a paradigmatic example of active matter. In these systems, unlike Brownian "hot colloids", the breakdown of detailed balance yields a continuous production of entropy at steady state, even for an ideal active gas. We quantify the irreversibility for a non-interacting active particle in two dimensions by treating both conjugated and time-reversed dynamics. By starting with underdamped dynamics, we identify a hidden rate of entropy production required to maintain persistence and prevent the rapidly relaxing momenta from thermalizing, even in the limit of very large friction. Additionally, comparing two popular models of self-propulsion with identical dissipation on average, we find that the fluctuations and large deviations in work done are markedly different, providing thermodynamic insight into the varying extents to which macroscopically similar active matter systems may depart from equilibrium. △ Less

Submitted 1 September, 2018; v1 submitted 9 April, 2018; originally announced April 2018.

Comments: 9 pages including SI, 3 figures (published version)

Journal ref: Phys. Rev. E 98, 020604 (2018)

Anomalous glassy dynamics in simple models of dense biological tissue

Authors: Daniel M. Sussman, M. Paoluzzi, M. Cristina Marchetti, M. Lisa Manning

Abstract: In order to understand the mechanisms for glassy dynamics in biological tissues and shed light on those in non-biological materials, we study the low-temperature disordered phase of 2D vertex-like models. Recently it has been noted that vertex models have quite unusual behavior in the zero-temperature limit, with rigidity transitions that are controlled by residual stresses and therefore exhibit v… ▽ More In order to understand the mechanisms for glassy dynamics in biological tissues and shed light on those in non-biological materials, we study the low-temperature disordered phase of 2D vertex-like models. Recently it has been noted that vertex models have quite unusual behavior in the zero-temperature limit, with rigidity transitions that are controlled by residual stresses and therefore exhibit very different scaling and phenomenology compared to particulate systems. Here we investigate the finite-temperature phase of two-dimensional Voronoi and Vertex models, and show that they have highly unusual, sub-Arrhenius scaling of dynamics with temperature. We connect the anomalous glassy dynamics to features of the potential energy landscape associated with zero-temperature inherent states. △ Less

Submitted 15 December, 2017; originally announced December 2017.

Comments: 8 pages, 8 figures

Stability from activity

Authors: Ananyo Maitra, Pragya Srivastava, M. Cristina Marchetti, Juho Lintuvuori, Sriram Ramaswamy, Martin Lenz

Abstract: Suspensions of actively driven anisotropic objects exhibit distinctively nonequilibrium behaviors, and current theories predict that they are incapable of sustaining orientational order at high activity. By contrast, here we show that nematic suspensions on a substrate can display order at arbitrarily high activity due to a previously unreported, potentially stabilizing active force. The resulting… ▽ More Suspensions of actively driven anisotropic objects exhibit distinctively nonequilibrium behaviors, and current theories predict that they are incapable of sustaining orientational order at high activity. By contrast, here we show that nematic suspensions on a substrate can display order at arbitrarily high activity due to a previously unreported, potentially stabilizing active force. The resulting nonequilibrium ordered phase displays robust giant number fluctuations that cannot be suppressed even by an incompressible solvent. Our results apply to virtually all experimental assays used to investigate the active nematic ordering of self-propelled colloids, bacterial suspensions and the cytoskeleton, and have testable implications in interpreting their nonequilibrium behaviors △ Less

Submitted 7 November, 2017; originally announced November 2017.