Shear-Induced First-Order Transition in Polar Liquid Crystals

Field Theory for Mechanical Criticality in Disordered Fiber Networks

Strain-controlled criticality governs the elasticity of jamming and fiber networks. While the upper critical dimension of jamming is believed to be d_{u}=2, non-mean-field exponents are observed in numerical studies of 2D and 3D fiber networks. The origins of this remains unclear. In this study we propose a minimal mean-field model for strain-controlled criticality of fiber networks. We then extend this to a phenomenological field theory, in which non-mean-field behavior emerges as a result of the disorder in the network structure. We predict that the upper critical dimension for such systems is d_{u}=4 using a Gaussian approximation. Moreover, we identify an order parameter for the phase transition, which has been lacking for fiber networks to date.

Mechanochemical Active Feedback Generates Convergence Extension in Epithelial Tissue

Convergence extension, the simultaneous elongation of tissue along one axis while narrowing along a perpendicular axis, occurs during embryonic development. A fundamental process that contributes to shaping the organism, it happens in many different species and tissue types. Here, we present a minimal continuum model, that can be directly linked to the controlling microscopic biochemistry, which shows spontaneous convergence extension. It is comprised of a 2D viscoelastic active material with a mechanochemical active feedback mechanism coupled to a substrate via friction. Robust convergent extension behavior emerges beyond a critical value of the activity parameter and is controlled by the boundary conditions and the coupling to the substrate. Oscillations and spatial patterns emerge in this model when internal dissipation dominates over friction, as well as in the active elastic limit.

Motor-free contractility of active biopolymer networks

Contractility in animal cells is often generated by molecular motors such as myosin, which require polar substrates for their function. Motivated by recent experimental evidence of motor-independent contractility, we propose a robust motor-free mechanism that can generate contraction in biopolymer networks without the need for substrate polarity. We show that contractility is a natural consequence of active binding-unbinding of crosslinkers that breaks the principle of detailed balance, together with the asymmetric force-extension response of semiflexible biopolymers. We have extended our earlier work to discuss the motor-free contraction of viscoelastic biopolymer networks. We calculate the resulting contractile velocity using a microscopic model and show that it can be reduced to a simple coarse-grained model under certain limits. Our model may provide an explanation of recent reports of motor-independent contractility in cells. Our results also suggest a mechanism for generating contractile forces in synthetic active materials.

Simulating structured fluids with tensorial viscoelasticity

We consider an immersed elastic body that is actively driven through a structured fluid by a motor or an external force. The behavior of such a system generally cannot be solved analytically, necessitating the use of numerical methods. However, current numerical methods omit important details of the microscopic structure and dynamics of the fluid, which can modulate the magnitudes and directions of viscoelastic restoring forces. To address this issue, we develop a simulation platform for modeling viscoelastic media with tensorial elasticity. We build on the lattice Boltzmann algorithm and incorporate viscoelastic forces, elastic immersed objects, a microscopic orientation field, and coupling between viscoelasticity and the orientation field. We demonstrate our method by characterizing how the viscoelastic restoring force on a driven immersed object depends on various key parameters as well as the tensorial character of the elastic response. We find that the restoring force depends non-monotonically on the rate of diffusion of the stress and the size of the object. We further show how the restoring force depends on the relative orientation of the microscopic structure and the pulling direction. These results imply that accounting for previously neglected physical features, such as stress diffusion and the microscopic orientation field, can improve the realism of viscoelastic simulations. We discuss possible applications and extensions to the method.

Spontaneous organization and phase separation of skyrmions in chiral active matter

Skyrmions are topologically protected vortex-like excitations that hold promise for applications such as information processing and electron manipulation. Here we combine theoretical analysis and numerical simulations to show that skyrmions can spontaneously emerge in chiral active matter without external confinements or regulation. Strikingly, these activity-driven skyrmions can either self-organize into a periodic, stable square lattice consisting of half Néel skyrmions and antiskyrmions, where the in-plane flows display an antiferromagnetic vortex array, or undergo phase separation between skyrmions with different topological numbers. We identify that the emerging skyrmion dynamics stems from the competition between the chiral and polar coherence length scales dictated by the interplay of intrinsic chirality, polarity, and elasticity in the system. Our results reveal unanticipated topological excitations, self-organization, and phase separation in non-equilibrium systems and also suggest a potential way towards engineering complicated bespoke skyrmionic structures through manipulating active matter.

Stochastic Hydrodynamics of Complex Fluids: Discretisation and Entropy Production

Many complex fluids can be described by continuum hydrodynamic field equations, to which noise must be added in order to capture thermal fluctuations. In almost all cases, the resulting coarse-grained stochastic partial differential equations carry a short-scale cutoff, which is also reflected in numerical discretisation schemes. We draw together our recent findings concerning the construction of such schemes and the interpretation of their continuum limits, focusing, for simplicity, on models with a purely diffusive scalar field, such as ‘Model B’ which describes phase separation in binary fluid mixtures. We address the requirement that the steady-state entropy production rate (EPR) must vanish for any stochastic hydrodynamic model in a thermal equilibrium. Only if this is achieved can the given discretisation scheme be relied upon to correctly calculate the nonvanishing EPR for ‘active field theories’ in which new terms are deliberately added to the fluctuating hydrodynamic equations that break detailed balance. To compute the correct probabilities of forward and time-reversed paths (whose ratio determines the EPR), we must make a careful treatment of so-called ‘spurious drift’ and other closely related terms that depend on the discretisation scheme. We show that such subtleties can arise not only in the temporal discretisation (as is well documented for stochastic ODEs with multiplicative noise) but also from spatial discretisation, even when noise is additive, as most active field theories assume. We then review how such noise can become multiplicative via off-diagonal couplings to additional fields that thermodynamically encode the underlying chemical processes responsible for activity. In this case, the spurious drift terms need careful accounting, not just to evaluate correctly the EPR but also to numerically implement the Langevin dynamics itself.

Macroscopic current generated by local division and apoptosis in a minimal model of tissue dynamics

It has been known that the motion of self-propelled particles inside an asymmetric channel can be rectified to give rise to a macroscopic and unidirectional current. In this paper, we show that division and apoptosis process alone (which are ubiquitous in all living systems), without any self-propulsion, are sufficient to give rise to a macroscopic and unidirectional current in a similar channel. More specifically, we consider a minimal computational model of two-dimensional living tissues, with two active ingredients: local particle division and apoptosis, and we found a net steady state current along the channel. This mechanism is fundamentally different from that of self-propelled particles, since this requires the system to be dense. Finally, our results might have applications in tissue engineering such as controlling tissue growth via a geometrically non-uniform substrate.

Permeation Instabilities in Active Polar Gels

We present a theory of active, permeating, polar gels, based on a two-fluid model. An active relative force between the gel components creates a steady-state current. We analyze its stability, while considering two polar coupling terms to the relative current: a permeation-deformation term, which describes network deformation by the solvent flow, and a permeation-alignment term, which describes the alignment of the polarization field by the network deformation and flow. Novel instability mechanisms emerge at finite wave vectors, suggesting the formation of periodic domains and mesophases. Our results can be used to determine the physical conditions required for various types of multicellular migration across tissues.

Shear-induced orientational ordering in an active glass former

Dense assemblies of self-propelled particles that can form solid-like states also known as active or living glasses are abundant around us, covering a broad range of length scales and timescales: from the cytoplasm to tissues, from bacterial biofilms to vehicular traffic jams, and from Janus colloids to animal herds. Being structurally disordered as well as strongly out of equilibrium, these systems show fascinating dynamical and mechanical properties. Using extensive molecular dynamics simulation and a number of distinct dynamical and mechanical order parameters, we differentiate three dynamical steady states in a sheared model active glassy system: 1) a disordered state, 2) a propulsion-induced ordered state, and 3) a shear-induced ordered state. We supplement these observations with an analytical theory based on an effective single-particle Fokker-Planck description to rationalize the existence of the shear-induced orientational ordering behavior in an active glassy system without explicit aligning interactions of, for example, Vicsek type. This ordering phenomenon occurs in the large persistence time limit and is made possible only by the applied steady shear. Using a Fokker-Planck description with parameters that can be measured independently, we make testable predictions for the joint distribution of single-particle position and orientation. These predictions match well with the joint distribution measured from direct numerical simulation. Our results are of relevance for experiments exploring the rheological response of dense active colloids and jammed active granular matter systems.

Shear-Driven Mechanism of Temperature Gradient Formation in Microfluidic Nematic Devices: Theory and Numerical Studies

The purpose of this paper is to show some routes in describing the mechanism responsible for the formation of the temperature difference ΔT at the boundaries of the microfluidic hybrid aligned nematic (HAN) channel, initially equal to zero, if one sets up the stationary hydrodynamic flow vst or under the effect of an externally applied shear stress (SS) to the bounding surfaces. Calculations based on the nonlinear extension of the classical Ericksen–Leslie theory, supplemented by thermomechanical correction of the SS σzx and Rayleigh dissipation function while accounting for the entropy balance equation, show that the ΔT is proportional to the heat flux q across the HAN channel and grows by up to several degrees under the influence of the externally applied SS. The role of vst=ust(z)i^ with a sharp triangular profile ust(z) across the HAN in the production of the highest ΔT is also investigated.

Odd Viscosity in Active Matter: Microscopic Origin and 3D Effects

In common fluids, viscosity is associated with dissipation. However, when time-reversal symmetry is broken a new type of nondissipative "viscosity" emerges. Recent theories and experiments on classical 2D systems with active spinning particles have heightened interest in "odd viscosity," but a microscopic theory for it in active materials is still absent. Here, we present such first-principles microscopic Hamiltonian theory, valid for both 2D and 3D, showing that odd viscosity is present in any system, even at zero temperature, with globally or locally aligned spinning components. Our work substantially extends the applicability of odd viscosity into 3D fluids, and specifically to internally driven active materials, such as living matter (e.g., actomyosin gels). We find intriguing 3D effects of odd viscosity such as propagation of anisotropic bulk shear waves and breakdown of Bernoulli's principle.

Motor-Free Contractility in Active Gels

Animal cells form contractile structures to promote various functions, from cell motility to cell division. Force generation in these structures is often due to molecular motors such as myosin that require polar substrates for their function. Here, we propose a motor-free mechanism that can generate contraction in biopolymer networks without the need for polarity. This mechanism is based on active binding and unbinding of cross-linkers that breaks the principle of detailed balance, together with the asymmetric force-extension response of semiflexible biopolymers. We find that these two ingredients can generate steady state contraction via a nonthermal, ratchetlike process. We calculate the resulting force-velocity relation using both coarse-grained and microscopic models.

Soft channel formation and symmetry breaking in exotic active emulsions

We use computer simulations to study the morphology and rheological properties of a bidimensional emulsion resulting from a mixture of a passive isotropic fluid and an active contractile polar gel, in the presence of a surfactant that favours the emulsification of the two phases. By varying the intensity of the contractile activity and of an externally imposed shear flow, we find three possible morphologies. For low shear rates, a simple lamellar state is obtained. For intermediate activity and shear rate, an asymmetric state emerges, which is characterized by shear and concentration banding at the polar/isotropic interface. A further increment in the active forcing leads to the self-assembly of a soft channel where an isotropic fluid flows between two layers of active material. We characterize the stability of this state by performing a dynamical test varying the intensity of the active forcing and shear rate. Finally, we address the rheological properties of the system by measuring the effective shear viscosity, finding that this increases as active forcing is increased-so that the fluid thickens with activity.

Microfluidic control over topological states in channel-confined nematic flows

Compared to isotropic liquids, orientational order of nematic liquid crystals makes their rheological properties more involved, and thus requires fine control of the flow parameters to govern the orientational patterns. In microfluidic channels with perpendicular surface alignment, nematics discontinuously transition from perpendicular structure at low flow rates to flow-aligned structure at high flow rates. Here we show how precise tuning of the driving pressure can be used to stabilize and manipulate a previously unresearched topologically protected chiral intermediate state which arises before the homeotropic to flow-aligned transition. We characterize the mechanisms underlying the transition and construct a phenomenological model to describe the critical behaviour and the phase diagram of the observed chiral flow state, and evaluate the effect of a forced symmetry breaking by introduction of a chiral dopant. Finally, we induce transitions on demand through channel geometry, application of laser tweezers, and careful control of the flow rate.

It is interesting phenomenon that chiral order can emerge in intrinsically achiral liquid crystals. Here Čopar et al. demonstrate achiral-to-chiral transition of the nematic liquid crystals flow in microfluidic channels and their behaviour, stability, and dependence on geometric and material parameters.