Dynamical patterns in nematic active matter on a sphere

Individual particle persistence antagonizes global ordering in populations of nematically aligning self-propelled particles

The transition from individual to collective motion plays a significant role in many biological processes. While the implications of different types of particle-particle interactions for the emergence of particular modes of collective motion have been well studied, it is unclear how particular types of individual migration patterns influence collective motion. Here, motivated by swarming bacteria Myxococcus xanthus, we investigate the combined effects of the individual pattern of migration and particle-particle interactions on the emergence of collective migration. We analyze the effects of a feature of individual pattern migration, the persistence of motion, on the collective properties of the system that emerge from interactions among individuals, particularly when nematic velocity alignment interaction mediates collective dynamics. We find, through computer simulations and mathematical analysis, that an initially disordered migratory state can become globally ordered by increasing either the particle-particle alignment interaction strength or the persistence of individual migration. In contrast, we find that persistence prevents the emergence of global nematic order when both persistence and nematic alignment are comparatively high. We conclude that behavior at the population level not only depends on interactions between individuals but also on their own intrinsic behavior.

Active smectics on a sphere

The dynamics of active smectic liquid crystals confined on a spherical surface is explored through an active phase field crystal model. Starting from an initially randomly perturbed isotropic phase, several types of topological defects are spontaneously formed, and then annihilate during a coarsening process until a steady state is achieved. The coarsening process is highly complex involving several scaling laws of defect densities as a function of time where different dynamical exponents can be identified. In general the exponent for the final stage towards the steady state is significantly larger than that in the passive and in the planar case, i.e. the coarsening is getting accelerated both by activity and by the topological and geometrical properties of the sphere. A defect type characteristic for this active system is a rotating spiral of evolving smectic layering lines. On a sphere this defect type also determines the steady state. Our results can in principle be confirmed by dense systems of synthetic or biological active particles.

Symmetry Breaking of Self-Propelled Topological Defects in Thin-Film Active Chiral Nematics

Active nematics represent a range of dense active matter systems which can engender spontaneous flows and self-propelled topological defects. Two-dimensional (2D) active nematic theory and simulation have been successful in explaining many quasi-2D experiments in which self-propelled +1/2 defects are observed to move along their symmetry axis. However, many active liquid crystals are essentially chiral nematic, but their twist mode becomes irrelevant under the 2D assumption. Here, we use theory and simulation to examine a three-dimensional active chiral nematic confined to a thin film, thus forming a quasi-2D system. We predict that the self-propelled +1/2 disclination in a curved thin film can break its mirror symmetry by moving circularly. Our prediction is confirmed by hydrodynamic simulations of thin spherical-shell and thin cylindrical-shell systems. In the spherical-shell confinement, the four emerged +1/2 disclinations exhibit rich dynamics as a function of activity and chirality. As such, we have proposed a new symmetry-breaking scenario in which self-propelled defects in quasi-2D active nematics can acquire an active angular velocity, greatly enriching their dynamics for finer control and emerging applications.

Effect of confinement and topology: 2-TIPS vs. MIPS

2-TIPS (two temperature induced phase separation) refers to the phase separation phenomenon observed in mixtures of active and passive particles which are modelled using scalar activity. The active particles are connected to a thermostat at high temperature while the passive particles are connected to the thermostat at low temperature and the relative temperature difference between "hot" and "cold" particles is taken as the measure of the activity χ of the non-equilibrium system. The study of such binary mixtures of hot and cold particles under various kinds of confinement is an important problem in many physical and biological processes. The nature and extent of phase separation are heavily influenced by the geometry of confinement, activity, and density of the non-equilibrium binary mixture. Investigating such 3D binary mixtures confined by parallel walls, we observe that the active and passive particles phase separate, but the extent of phase separation is reduced compared to bulk phase separation at high densities and enhanced at low densities. However, when the binary mixture of active and passive particles is confined inside a spherical cavity, the phase separation is radial for small radii of the confining sphere and the extent of phase separation is higher compared to their bulk counterparts. Confinement leads to interesting properties in the passive (cold) region like enhanced layering and high compression in the direction parallel to the confining wall. In 2D, both the bulk and confined systems of the binary mixture show a significant decrement in the extent of phase separation at higher densities. This observation is attributed to the trapping of active particles inside the passive cluster, which increases with density. Thus the 2D systems show structures more akin to dense-dilute phase co-existence, which is observed in motility induced phase separation in 2D active systems. The binary mixture constrained on the spherical surface also shows similar phase co-existence. Our analyses reveal that the coexistent densities observed in 2-TIPS on the spherical surface agree with the findings of previous studies on MIPS in active systems on a sphere.

Fractional and scaled Brownian motion on the sphere: The effects of long-time correlations on navigation strategies

We analyze fractional Brownian motion and scaled Brownian motion on the two-dimensional sphere S^{2}. We find that the intrinsic long-time correlations that characterize fractional Brownian motion collude with the specific dynamics (navigation strategies) carried out on the surface giving rise to rich transport properties. We focus our study on two classes of navigation strategies: one induced by a specific set of coordinates chosen for S^{2} (we have chosen the spherical ones in the present analysis), for which we find that contrary to what occurs in the absence of such long-time correlations, nonequilibrium stationary distributions are attained. These results resemble those reported in confined flat spaces in one and two dimensions [Guggenberger et al. New J. Phys. 21, 022002 (2019)1367-263010.1088/1367-2630/ab075f; Vojta et al. Phys. Rev. E 102, 032108 (2020)2470-004510.1103/PhysRevE.102.032108]; however, in the case analyzed here, there are no boundaries that affect the motion on the sphere. In contrast, when the navigation strategy chosen corresponds to a frame of reference moving with the particle (a Frenet-Serret reference system), then the equilibrium distribution on the sphere is recovered in the long-time limit. For both navigation strategies, the relaxation times toward the stationary distribution depend on the particular value of the Hurst parameter. We also show that on S^{2}, scaled Brownian motion, distinguished by a time-dependent diffusion coefficient with a power-scaling, is independent of the navigation strategy finding a good agreement between the analytical calculations obtained from the solution of a time-dependent diffusion equation on S^{2}, and the numerical results obtained from our numerical method to generate ensemble of trajectories.

Phase behavior and dynamics in a colloid-polymer mixture under spherical confinement

From studies via molecular dynamics simulations, we report results on structure and dynamics in mixtures of active colloids and passive polymers that are confined inside a spherical container with a repulsive boundary. All interactions in the fully passive limit are chosen in such a way that in equilibrium coexistence between colloid-rich and polymer-rich phases occurs. For most part of the studies the chosen compositions give rise to Janus-like structure: nearly one side of the sphere is occupied by the colloids and the rest by the polymers. This partially wet situation mimics approximately a neutral wall in the fully passive scenario. Following the introduction of a velocity-aligning activity to the colloids, the shape of the polymer-rich domain changes to that of an ellipsoid, around the long axis of which the colloid-rich domain attains a macroscopic angular momentum. In the steady state, the orientation of this axis evolves via diffusion, magnitude of which depends upon the strength of activity, but only weakly.

Tuneable defect-curvature coupling and topological transitions in active shells

Recent experimental observations have suggested that topological defects can facilitate the creation of sharp features in developing embryos. Whereas these observations echo established knowledge about the interplay between geometry and topology in two-dimensional passive liquid crystals, the role of activity has mostly remained unexplored. In this article we focus on deformable shells consisting of either polar or nematic active liquid crystals and demonstrate that activity renders the mechanical coupling between defects and curvature much more involved and versatile than previously thought. Using a combination of linear stability analysis and three-dimensional computational fluid dynamics, we demonstrate that such a coupling can in fact be tuned, depending on the type of liquid crystal order, the specific structure of the defect (i.e. asters or vortices) and the nature of the active forces. In polar systems, this can drive a spectacular transition from spherical to toroidal topology, in the presence of large extensile activity. Our analysis strengthens the idea that defects could serve as topological morphogens and provides a number of predictions that could be tested in in vitro studies, for instance in the context of organoids.

Activity-induced polar patterns of filaments gliding on a sphere

Active matter systems feature the ability to form collective patterns as observed in a plethora of living systems, from schools of fish to swimming bacteria. While many of these systems move in a wide, three-dimensional environment, several biological systems are confined by a curved topology. The role played by a non-Euclidean geometry on the self-organization of active systems is not yet fully understood, and few experimental systems are available to study it. Here, we introduce an experimental setup in which actin filaments glide on the inner surface of a spherical lipid vesicle, thus embedding them in a curved geometry. We show that filaments self-assemble into polar, elongated structures and that, when these match the size of the spherical geometry, both confinement and topological constraints become relevant for the emergent patterns, leading to the formation of polar vortices and jammed states. These results experimentally demonstrate that activity-induced complex patterns can be shaped by spherical confinement and topology.

Vesicle shape transformations driven by confined active filaments

In active matter systems, deformable boundaries provide a mechanism to organize internal active stresses. To study a minimal model of such a system, we perform particle-based simulations of an elastic vesicle containing a collection of polar active filaments. The interplay between the active stress organization due to interparticle interactions and that due to the deformability of the confinement leads to a variety of filament spatiotemporal organizations that have not been observed in bulk systems or under rigid confinement, including highly-aligned rings and caps. In turn, these filament assemblies drive dramatic and tunable transformations of the vesicle shape and its dynamics. We present simple scaling models that reveal the mechanisms underlying these emergent behaviors and yield design principles for engineering active materials with targeted shape dynamics.

Alignment interactions drive structural transitions in biological tissues

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 nonequilibrium phase transition into a quasihexagonal 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.

Dynamics of active nematic defects on the surface of a sphere

A nematic liquid crystal confined to the surface of a sphere exhibits topological defects of total charge +2 due to the topological constraint. In equilibrium, the nematic field forms four +1/2 defects, located at the corners of a regular tetrahedron inscribed within the sphere, since this minimizes the Frank elastic energy. If additionally the individual nematogens exhibit self-driven directional motion, the resulting active system creates large-scale flow that drives it out of equilibrium. In particular, the defects now follow complex dynamic trajectories which, depending on the strength of the active forcing, can be periodic (for weak forcing) or chaotic (for strong forcing). In this paper we derive an effective particle theory for this system, in which the topological defects are the degrees of freedom, whose exact equations of motion we subsequently determine. Numerical solutions of these equations confirm previously observed characteristics of their dynamics and clarify the role played by the time dependence of their global rotation. We also show that Onsager's variational principle offers an exceptionally transparent way to derive these dynamical equations, and we explain the defect mobility at the hydrodynamics level.

Chemically Active Particles: From One to Few on the Way to Many

Chemically active particles suspended in a liquid solution can achieve self-motility by locally changing the chemical composition of the solution via catalytic reactions at their surfaces. They operate intrinsically out of equilibrium, continuously extracting free energy from the environment to power the dissipative self-motility. The effective interactions involving active particles are, in general, nonreciprocal and anisotropic, even if the particles have simple shapes (e.g., Janus spheres). Accordingly, for chemically active particles a very rich behavior of collective motion and self-assembly may be expected to emerge, including phenomena such as microphase separation in the form of kinetically stable, finite-sized aggregates. Here, I succinctly review a number of recent experimental studies that demonstrate the self-assembly of structures, involving chemically active Janus particles, which exhibit various patterns of motion. These examples illustrate concepts such as "motors made out of motors" (as suggestively named by Fischer [Fischer, P. Nat. Phys. 2018, 14, 1072]). The dynamics of assembly and structure formation observed in these systems can provide benchmark, in-depth testing of the current understanding of motion and effective interactions produced by chemical activity. Finally, one notes that these significant achievements are likely just the beginning of the field. Recently reported particles endowed with time-dependent chemical activity or switchable reaction mechanisms open the way for exciting developments, such as periodic reshaping of self-assembled structures based on man-made internal clocks.

Newton-Cartan submanifolds and fluid membranes

We develop the geometric description of submanifolds in Newton-Cartan spacetime. This provides the necessary starting point for a covariant spacetime formulation of Galilean-invariant hydrodynamics on curved surfaces. We argue that this is the natural geometrical framework to study fluid membranes in thermal equilibrium and their dynamics out of equilibrium. A simple model of fluid membranes that only depends on the surface tension is presented and, extracting the resulting stresses, we show that perturbations away from equilibrium yield the standard result for the dispersion of elastic waves. We also find a generalization of the Canham-Helfrich bending energy for lipid vesicles that takes into account the requirements of thermal equilibrium.

Binary mixtures of active and passive particles on a sphere

We study the cooperation and segregation dynamics of binary mixtures of active and passive particles on a sphere. According to the competition between rotational diffusion and polar alignment, we find three distinct phases: a mixed phase and two different demixed phases. When rotational diffusion dominates the dynamics, the demixing is due to the aggregation of passive particles, where active and passive particles respectively occupy two hemispheres. When polar alignment is dominated, the demixing is caused by the aggregation of active particles, where active particles occupy the equator of the sphere and passive particles occupy the two poles of the sphere. In this case, there exist a circulating band cluster and two cambered surface clusters, which is a purely curvature-driven effect with no equivalent in the planar model. When rotational diffusion and polar alignment are comparable, particles are completely mixed. Our findings are relevant to the experimental pursuit of segregation dynamics of binary mixtures on curved surfaces.

A theoretical phase diagram for an active nematic on a spherical surface

Systems combining rod-shaped objects with self-generated motion such as suspensions of microtubules or growing bacterial colonies are commonly modeled as active nematics - nematic liquid crystals with an additional active stress term. Confining a 2D active nematic to the surface of a sphere generates novel behaviour as the four +1/2 nematic defects which are produced by the spherical geometry move round each other in an intricate dance. Here, these defects are modeled as point particles experiencing elastic forces from defect position and orientation, and self-propulsion due to activity. This model exhibits four qualitatively distinct types of trajectory state: two which are consistent with previous experimental and simulated trajectories; and two others, which are apparently novel and in regions of parameter space that may not yet have been explored. This work also explains a discrepancy between some previous point-particle models and the trajectories seen in experiments and simulations: this was due to a failure to fully account for the spherical geometry in the point-particle models.

Three-Dimensional Active Defect Loops

We describe the flows and morphological dynamics of topological defect lines and loops in three-dimensional active nematics and show, using theory and numerical modeling, that they are governed by the local profile of the orientational order surrounding the defects. Analyzing a continuous span of defect loop profiles, ranging from radial and tangential twist to wedge ±1/2 profiles, we show that the distinct geometries can drive material flow perpendicular or along the local defect loop segment, whose variation around a closed loop can lead to net loop motion, elongation, or compression of shape, or buckling of the loops. We demonstrate a correlation between local curvature and the local orientational profile of the defect loop, indicating dynamic coupling between geometry and topology. To address the general formation of defect loops in three dimensions, we show their creation via bend instability from different initial elastic distortions.

Spontaneous helical flows in active nematics lying on a cylindrical surface

Within the framework of the two-dimensional Ericksen-Leslie model, we explore the effect of geometric confinement on the spontaneous flow of active nematic gels. The nematic particles are assumed to flow on a cylindrical surface, while a degenerate tangential anchoring is enforced. Using the linear approximation of the motion equations, we show that there is a close interplay among extrinsic curvature, flow, director alignment, and activity. We find that the extrinsic curvature promotes the director alignment parallel to the cylindrical axis and is responsible for raising the critical threshold with respect to the flat case. Our analysis reveals a very rich scenario where the key quantities are the activity coefficient, the tumbling parameter, and the anisotropic viscosity ratio. Thus, solutions can exhibit a double periodicity in both the azimuthal and axial variables. As a consequence, the velocity field can make a finite angle with the cylinder axis and the active flow winds on the surface with a helical pattern, while the director oscillates around the cylinder generators. Our results can be validated on thin layers of nematic gels placed between two concentric cylinders and suggest which material properties are most suited for the design of active microfluidic devices.

Multicomponent flow on curved surfaces: A vielbein lattice Boltzmann approach

We develop and implement a finite difference lattice Boltzmann scheme to study multicomponent flows on curved surfaces, coupling the continuity and Navier-Stokes equations with the Cahn-Hilliard equation to track the evolution of the binary fluid interfaces. The standard lattice Boltzmann method relies on regular Cartesian grids, which makes it generally unsuitable to study flow problems on curved surfaces. To alleviate this limitation, we use a vielbein formalism to write the Boltzmann equation on an arbitrary geometry, and solve the evolution of the fluid distribution functions using a finite difference method. Focusing on the torus geometry as an example of a curved surface, we demonstrate drift motions of fluid droplets and stripes embedded on the surface of the torus. Interestingly, they migrate in opposite directions: fluid droplets to the outer side while fluid stripes to the inner side of the torus. For the latter we demonstrate that the global minimum configuration is unique for small stripe widths, but it becomes bistable for large stripe widths. Our simulations are also in agreement with analytical predictions for the Laplace pressure of the fluid stripes, and their damped oscillatory motion as they approach equilibrium configurations, capturing the corresponding decay timescale and oscillation frequency. Finally, we simulate the coarsening dynamics of phase separating binary fluids in the hydrodynamics and diffusive regimes for tori of various shapes, and compare the results against those for a flat two-dimensional surface. Our finite difference lattice Boltzmann scheme can be extended to other surfaces and coupled to other dynamical equations, opening up a vast range of applications involving complex flows on curved geometries.

Topology and Morphology of Self-Deforming Active Shells

We present a generic framework for modeling three-dimensional deformable shells of active matter that captures the orientational dynamics of the active particles and hydrodynamic interactions on the shell and with the surrounding environment. We find that the cross talk between the self-induced flows of active particles and dynamic reshaping of the shell can result in conformations that are tunable by varying the form and magnitude of active stresses. We further demonstrate and explain how self-induced topological defects in the active layer can direct the morphodynamics of the shell. These findings are relevant to understanding morphological changes during organ development and the design of bioinspired materials that are capable of self-organization.

Emergent topological phenomena in active polymeric fluids

Polymeric fluids show a wealth of topological phenomena, from entanglement and reptation at microscales to orientational ordering and defect production at macroscales, which can be explained by statistical-mechanical theories. In the presence of activity, the latter must be augmented by forces that cause spontaneous chain motion and fluid flow. Here, using such augmented Langevin equations, we study active polymeric solutions and melts composed of chains of hydrodynamically interacting stresslets. In a spherical volume, contractile chains are unstable and self-knot into entangled melts at both low and high densities. Extensile chains in the same geometry form an unentangled reptating state at low densities and an entangled, coherently moving, non-reptating state at high densities. On a spherical surface, contractile chains show transitions, with increasing areal density, between isotropic, orientationally ordered and micro-phase separated states. Extensile chains in the same geometry show a transition between isotropic and nematic states. In both cases, defects in orientationally ordered states are produced athermally and without conserving topological charge. Our work reproduces the phenomenology of several recent experiments, highlights the importance of hydrodynamic interactions in active polymer fluids, and suggests non-equilibrium kinetic routes to topological structures that are otherwise difficult to obtain in equilibrium.

The interplay between activity and filament flexibility determines the emergent properties of active nematics

Active nematics are microscopically driven liquid crystals that exhibit dynamical steady states characterized by the creation and annihilation of topological defects. Motivated by differences between previous simulations of active nematics based on rigid rods and experimental realizations based on semiflexible biopolymer filaments, we describe a large-scale simulation study of a particle-based computational model that explicitly incorporates filament semiflexibility. We find that energy injected into the system at the particle scale preferentially excites bend deformations, reducing the apparent filament bend modulus. The emergent characteristics of the active nematic depend on activity and flexibility only through this activity-renormalized bend 'modulus', demonstrating that apparent values of material parameters, such as the Frank 'constants', depend on activity. Thus, phenomenological parameters within continuum hydrodynamic descriptions of active nematics must account for this dependence. Further, we present a systematic way to estimate these parameters from observations of deformation fields and defect shapes in experimental or simulation data.

Extensile motor activity drives coherent motions in a model of interphase chromatin

The 3D spatiotemporal organization of the human genome inside the cell nucleus remains a major open question in cellular biology. In the time between two cell divisions, chromatin-the functional form of DNA in cells-fills the nucleus in its uncondensed polymeric form. Recent in vivo imaging experiments reveal that the chromatin moves coherently, having displacements with long-ranged correlations on the scale of micrometers and lasting for seconds. To elucidate the mechanism(s) behind these motions, we develop a coarse-grained active polymer model where chromatin is represented as a confined flexible chain acted upon by molecular motors that drive fluid flows by exerting dipolar forces on the system. Numerical simulations of this model account for steric and hydrodynamic interactions as well as internal chain mechanics. These demonstrate that coherent motions emerge in systems involving extensile dipoles and are accompanied by large-scale chain reconfigurations and nematic ordering. Comparisons with experiments show good qualitative agreement and support the hypothesis that self-organizing long-ranged hydrodynamic couplings between chromatin-associated active motor proteins are responsible for the observed coherent dynamics.

Active motion on curved surfaces

A theoretical analysis of active motion on curved surfaces is presented in terms of a generalization of the telegrapher equation. Such a generalized equation is explicitly derived as the polar approximation of the hierarchy of equations obtained from the corresponding Fokker-Planck equation of active particles diffusing on curved surfaces. The general solution to the generalized telegrapher equation is given for a pulse with vanishing current as initial data. Expressions for the probability density and the mean squared geodesic displacement are given in the limit of weak curvature. As an explicit example of the formulated theory, the case of active motion on the sphere is presented, where oscillations observed in the mean squared geodesic displacement are explained.

Active crystals on a sphere

Two-dimensional crystals on curved manifolds exhibit nontrivial defect structures. Here we consider "active crystals" on a sphere, which are composed of self-propelled colloidal particles. Our work is based on a phase-field-crystal-type model that involves a density and a polarization field on the sphere. Depending on the strength of the self-propulsion, three different types of crystals are found: a static crystal, a self-spinning "vortex-vortex" crystal containing two vortical poles of the local velocity, and a self-translating "source-sink" crystal with a source pole where crystallization occurs and a sink pole where the active crystal melts. These different crystalline states as well as their defects are studied theoretically here and can in principle be confirmed in experiments.