## A ctive M atter I nvasion

P reprint , compiled A ugust 5, 2019

Felix Kempf , Romain Mueller , Erwin Frey , Julia M. Yeomans , and Amin Doostmohammadi 1 2 1 2 2 GLYPH&lt;3&gt;

1 Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München - Theresienstr. 37, D-80333 Munich, Germany. 2 The Rudolf Peierls Centre for Theoretical Physics - Clarendon Laboratory, Parks Road, Oxford, OX1 3PU, UK.

## A bstract

Biological active materials such as bacterial biofilms and eukaryotic cells thrive in confined micro-spaces. Here, we show through numerical simulations that confinement can serve as a mechanical guidance to achieve distinct modes of collective invasion when combined with growth dynamics and the intrinsic activity of biological materials. We assess the dynamics of the growing interface and classify these collective modes of invasion based on the activity of the constituent particles of the growing matter. While at small and moderate activities the active material grows as a coherent unit, we find that blobs of active material collectively detach from the cohort above a well-defined activity threshold. We further characterise the mechanical mechanisms underlying the crossovers between di GLYPH&lt;11&gt; erent modes of invasion and quantify their impact on the overall invasion speed.

## 1 I ntroduction

Understanding the mechanisms by which living (active) systems such as cells and bacteria invade and navigate through their surroundings is of pivotal importance in many physiological and pathological processes. Depending on the physical and chemical properties of their microenvironment active matter can show distinct modes of invasion, from single cell migration to groups of cells moving as a collective. The latter has been identified as the primary mode of cancer cell invasion [1, 2]. Similarly collective migration is the prime mode of growth and invasion by bacterial biofilms [3, 4, 5].

Among physical environmental factors, geometrical constraints have a strong influence on the dynamics of collective migration in biological matter. Various physiological processes such as bacterial filtering rely on active matter living within pores, cavities, and constrictions [6]. In vivo , cancerous cells are known to preferentially move along pre-existing tracks of least resistance, such as myelinated axons or blood vessels, when invading into healthy tissue [7, 8]. Similarly, during biofilm formation interaction between the bacterial colony and the surface it grows on is of major importance [9, 10]. The relevance of the interaction between biological matter and confining geometries in many physiological processes has prompted extensive experiments. In particular, in-vitro experiments on epithelial cells [11, 12, 13, 14, 15, 16, 17], bacteria [18, 19], and mixtures of biofilaments and molecular motors [20, 21] have shown that confinement to channel-like geometries generically alters the flow of biological matter in a significant way. Striking examples are the crossover from chaotic flows of bacteria [22, 18] and microtubule motor protein mixtures [23, 24] to polarised / movement along the long channel axis and the emergence of shear flow in fibroblast cells [17] confined within rectangular geometries. Moreover, theoretical and computational studies of active matter interacting with a periodic array of obstacles have shown a reduction in the e GLYPH&lt;11&gt; ective di GLYPH&lt;11&gt; usion coe GLYPH&lt;14&gt; cient of active particles with increasing density of obstacles [25] and modification of the active dispersion within periodic arrays in the presence of applied flows [26].

While recent research has been primarily focused on the crosstalk between the confinement and the activity of the particles in order to determine patterns of motion, a complete understanding of the patterns of collective invasion and spreading also requires a detailed investigation of growth dynamics. From the physical perspective, such interplay between activity, confinement, and growth can be accompanied by additional complexities arising from hydrodynamic interactions between growing active matter and the confinement, orientation dynamics of elongated active particles, and emergent collective phenomena such as active turbulence [27, 28, 29, 30, 31] and active topological defects [32, 33, 34, 35]. Due to this interconnection of complex physical processes, the mechanistic understanding of active matter invasion within geometrical constraints remains largely unexplored. Recent numerical analyses of a self-propelled particle model [36], neglecting orientational dynamics and hydrodynamic e GLYPH&lt;11&gt; ects, have shown that the physics of growing active matter can explain some of the observed experimental phenomena [11], such as caterpillar motion of the advancing front and the enhancement of collective migration speed in thin capillaries, by assuming a coupling between the curvature of the front and the motility of the particles at the leading edge.

Here we use a generic model of active matter that accounts for hydrodynamics, orientational e GLYPH&lt;11&gt; ects, and growth to investigate collective patterns of active matter invasion of capillaries. We identify three di GLYPH&lt;11&gt; erent regimes of invasion, each with distinct interface shapes, flow patterns, orientational ordering, and topological defect dynamics. In particular, we show that above a certain threshold in the strength of the active forces, highly dynamic deformations are formed at the interface between the invading active matter and the surrounding medium. At higher activities, we also find a second threshold beyond which blobs of active matter begin to detach from the growing active column to enhance the invasion of free space. We explain the first crossover in terms of intrinsic hydrodynamic instabilities of the bulk active matter accompanied by additional instabilities arising from the presence of the growing interface. The latter can be understood by the ability of active stresses to overcome the stabilising e GLYPH&lt;11&gt; ect of the surface tension and pinching o GLYPH&lt;11&gt; the growing interface.

Figure 1: Simulation setup (A) Schematic drawing of simulation setup. The active phase (yellow) invades from the broader reservoir (region below black dashed line) where it can grow into the narrower capillary (width d ) that is initially filled with isotropic liquid (green). The broad gray arrow illustrates the invasion direction. Di GLYPH&lt;11&gt; erent heights are marked for later discussion of observables: h max, the highest point of the active phase which is connected to the reservoir, h inv , the highest point that any patch of active phase reaches, and h min, the highest point where the capillary is filled with active fluid to the full width. (B) A representative simulation snapshot to illustrate the dynamics. Black bars mark the director field; the two green speckles inside the nematic phase are cores of + / 1 2 topological defects.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00768v1-figure-8ef4457b-b7d5-444d-947c-6f2b5c54b837.jpg)

Together, our study reveals several possible invasion patterns and collective dynamics when active matter invades a capillary constriction thus providing a framework for further experimental investigations.

## 2 M odel

We employ a two-phase model of active nematohydrodynamics [37] to explore the growth of an active layer into otherwise isotropic surroundings within a confined channel and use a hybrid lattice Boltzmann (LB) method to numerically solve the equations [38]. This choice was motivated by the generality of this class of continuum models which is due to the fact that only local conservation laws, a nematic interaction between the systems' constituents, and perpetual injection of energy at the smallest length-scale are assumed. Models of this kind reproduce a variety of non-equilibrium flows such as stable arrays of vortices or chaotic flows, together with diverse stationary and non-stationary patterns of nematic ordering and topological defects [39, 40, 33, 29, 28, 41, 42, 43]. On a phenomenological level, these are successful in modelling the collective dynamics in biological systems in cases such as microtubule motor-protein / mixtures [44, 33, 40], cellular monolayers [45, 46, 47, 35], or bacteria [48]. In addition to the bulk dynamics, the interaction with walls or obstacles of di GLYPH&lt;11&gt; erent shape can add further complexity. Theoretically, using the active nematohydrodynamics framework, a transition to spontaneous flows in channels has been predicted [49, 50], as well as more complex states with intricate interplay of defects and vortices, and a transition to active turbulence [51, 52, 53].

## 2.1 Governing equations

We consider a two-dimensional model of a two-phase system consisting of an isotropic fluid phase and an active nematic phase [37]. The orientational order in the nematic phase is characterised by the symmetric and traceless nematic tensor Q GLYPH&lt;11&gt;GLYPH&lt;12&gt; = S GLYPH&lt;16&gt; 2 n n GLYPH&lt;11&gt; GLYPH&lt;12&gt; GLYPH&lt;0&gt; GLYPH&lt;14&gt; GLYPH&lt;11&gt;GLYPH&lt;12&gt; GLYPH&lt;17&gt; [54] with S (magnitude of the nematic order) and n GLYPH&lt;11&gt; (director) indicating the magnitude and the direction of the nematic order, respectively. The relative density of the nematic phase is measured by a scalar phase field GLYPH&lt;30&gt; which is 0 in the purely isotropic and 1 in the purely nematic phase.

The free energy density of the system is given by (using the Einstein summation convention):

$$\text{lation} \quad f = \frac { 1 } { 2 } D \phi ^ { 2 } \left ( 1 - \phi \right ) ^ { 2 } + \frac { 1 } { 2 } C _ { L Q } \left ( \phi - \frac { 1 } { 2 } Q _ { \alpha \beta } Q _ { \alpha \beta } \right ) ^ { 2 } \\ \text{to the} \quad \text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$
$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$                                                                                                                                                                                                         \text{                                                                                                                                                                                                       
                                                                                                                                                                                                        }                                                                                                                                                                                                       .                                                                                                                                                                                                      .                                                                                                                                                                                                        -                                                                                                                                                                                                        +                                                                                                                                                                                                        <                                                                                                                                                                                                        </                                                                                                                                                                                                        >                                                                                                                                                                                                        :                                                                                                                                                                                                        |                                                                                                                                                                                                        ($$

where D CLQ , , K GLYPH&lt;30&gt; , and KQ are positive material constants. Since Q GLYPH&lt;11&gt;GLYPH&lt;12&gt; Q GLYPH&lt;11&gt;GLYPH&lt;12&gt; = 2 S 2 , the second term ensures a tight coupling between the magnitude of the nematic order S , and the phase field GLYPH&lt;30&gt; .. Together with the first term, which corresponds to a Cahn-Hilliard free energy [55, 56], this leads to well-defined interfaces betweens a nematic ( S = 1 ; GLYPH&lt;30&gt; = 1) and an isotropic ( S = 0 ; GLYPH&lt;30&gt; = 0) phase (Fig. 1). The third and fourth terms are elastic energies. The first, third, and fourth terms contribute to the surface energy and the fourth term also penalises bulk deformations in Q GLYPH&lt;11&gt;GLYPH&lt;12&gt; . The free energy then reads

$$\mathcal { F } = \int \mathfrak { d } ^ { 2 } \mathbf r \, f.$$

The order parameters evolve according to the following equations:

$$\partial _ { t } \phi + \partial _ { \beta } \left ( \phi u _ { \beta } \right ) = & \Gamma _ { \phi } \Delta \mu,$$

$$\L a m { - \quad } \left ( \partial _ { t } + u _ { \kappa } \partial _ { \kappa } \right ) Q _ { \alpha \beta } = - \, \xi \Sigma _ { \alpha \beta \kappa \lambda } E _ { \kappa \lambda } - T _ { \alpha \beta \kappa \lambda } \Omega _ { \kappa \lambda } + \Gamma _ { Q } H _ { \alpha \beta }, \quad ( 4 )$$

with

$$\mu = & \frac { \partial f } { \partial \phi } - \partial _ { \gamma } \frac { \partial f } { \partial \left ( \partial _ { \gamma } \phi \right ) }$$

$$\ddot { \text{lat} } \begin{array} { c } \bar { \cdot } _ { \Sigma } & \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar {  \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot  } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar{ \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot \ } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar \cdot \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar {\cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot } \bar { \cdot }
 \ddot \end{array} \,,$$

$$\text{ray} \quad \text{E} _ { \alpha \beta } = & \frac { 1 } { 2 } \left ( \partial _ { \beta } u _ { \alpha } + \partial _ { \alpha } u _ { \beta } \right ), \quad \Omega _ { \alpha \beta } = \frac { 1 } { 2 } \left ( \partial _ { \beta } u _ { \alpha } - \partial _ { \alpha } u _ { \beta } \right ), \quad ( 7 )$$

GLYPH&lt;6&gt;

GLYPH&lt;0&gt;

GLYPH&lt;16&gt;

GLYPH&lt;17&gt;

GLYPH&lt;11&gt;GLYPH&lt;12&gt;GLYPH&lt;20&gt;GLYPH&lt;21&gt;

=

Q

GLYPH&lt;11&gt;GLYPH&lt;12&gt;

Q

GLYPH&lt;20&gt;GLYPH&lt;21&gt;

GLYPH&lt;14&gt;

GLYPH&lt;11&gt;GLYPH&lt;20&gt;

Q

GLYPH&lt;21&gt;GLYPH&lt;12&gt;

+

GLYPH&lt;14&gt;

GLYPH&lt;21&gt;GLYPH&lt;12&gt;

$$- \left ( Q _ { \alpha \lambda } + \delta _ { \alpha \lambda } \right ) \delta _ { \kappa \beta } + \delta _ { \alpha \beta } \left ( Q _ { \kappa \lambda } + \delta _ { \kappa \lambda } \right ), \text{ \quad \ \ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text$$

$$\stackrel { 1, \leq } { \text{tion} } \quad T _ { \alpha \beta \kappa \lambda } = & Q _ { \alpha \kappa } \delta _ { \beta \lambda } - \delta _ { \alpha \kappa } Q _ { \beta \lambda }.$$

The l.h.s. of Eqs. (3) and (4) are convective derivatives with the underlying velocity field u GLYPH&lt;11&gt; . Advection and di GLYPH&lt;11&gt; usion drive the dynamics in GLYPH&lt;30&gt; ; GLYPH&lt;0&gt; GLYPH&lt;30&gt; is the corresponding di GLYPH&lt;11&gt; usion constant. GLYPH&lt;0&gt; Q is a rotational di GLYPH&lt;11&gt; usion constant which, together with the molecular field H GLYPH&lt;11&gt;GLYPH&lt;12&gt; , controls di GLYPH&lt;11&gt; usive relaxation in Q GLYPH&lt;11&gt;GLYPH&lt;12&gt; . In addition, the interplay of flow u GLYPH&lt;11&gt; and order Q GLYPH&lt;11&gt;GLYPH&lt;12&gt; is less trivial. The first and second terms on the r.h.s. of the Eq. (4) form the co-rotational derivative which accounts for the response of the orientation field to the extensional ( E GLYPH&lt;11&gt;GLYPH&lt;12&gt; ) and rotational ( GLYPH&lt;10&gt; GLYPH&lt;11&gt;GLYPH&lt;12&gt; ) components of the velocity gradients respectively. GLYPH&lt;24&gt; is the tumbling parameter which determines the relative influence of the rate of strain on the director orientation. It depends on the geometry of the active particles, for prolate ellipsoids GLYPH&lt;24&gt; &gt; 0, while for oblate ellipsoids GLYPH&lt;24&gt; &lt; 0, and for spherical particles GLYPH&lt;24&gt; = 0 [57].

The velocity field u GLYPH&lt;11&gt; obeys the Navier-Stokes equations:

$$\partial _ { t } \rho + \partial _ { \beta } \left ( \rho u _ { \beta } \right ) = & 0, & ( 1 0 ) \underset { \text{phase} } { \text{ is a } }$$

$$\rho \left ( \partial _ { t } + u _ { \beta } \partial _ { \beta } \right ) u _ { \alpha } = & \partial _ { \beta } \Pi _ { \alpha \beta }, \quad \ \ \ ( 1 1 ) \ \frac { \text{flow} } { \text{eau} }.$$

where GLYPH&lt;26&gt; is the density of the fluid and GLYPH&lt;5&gt; GLYPH&lt;11&gt;GLYPH&lt;12&gt; is the stress tensor comprising viscous stress, pressure contribution, elastic stresses, and the active stress:

$$\Pi _ { \alpha \beta } ^ { v i s c } = & 2 \rho \eta E _ { \alpha \beta }, & ( 1 2 ) \, \frac { \text{aviam} } { \text{sample} }$$

$$\Pi _ { \alpha \beta } ^ { p } & = - \, \frac { \rho } { 3 } \delta _ { \alpha \beta }, & ( 1 3 ) \, \ \text{conditi} \\ \dots &$$

$$\Pi _ { \alpha \beta } ^ { e l, 1 } = \left ( f - \mu \phi \right ) \delta _ { \alpha \beta }$$

$$- \, \frac { \partial f } { \partial \left ( \partial _ { \beta } \phi \right ) } \partial _ { \alpha } \phi - \frac { \partial f } { \partial \left ( \partial _ { \beta } Q _ { \kappa \lambda } \right ) } \partial _ { \alpha } Q _ { \kappa \lambda }, \quad \ ( 1 4 ) \, \overset { \ v a l a } { \ e m e r } ; \text{syste} ; \text{research}.$$

$$\Pi _ { \alpha \beta } ^ { e l, 2 } & = \left ( \xi \Sigma _ { \alpha \beta \kappa \lambda } + T _ { \alpha \beta \kappa \lambda } \right ) H _ { \kappa \lambda }, & ( 1 5 ) \underset { \text{perpen} } { \text{perpen} }$$

$$\Pi _ { \alpha \beta } ^ { a c t } & = - \zeta \phi Q _ { \alpha \beta }. & ( 1 6 ) \quad \text{In the 1}$$

Here, GLYPH&lt;17&gt; is the viscosity, and the elastic stresses GLYPH&lt;5&gt; el i ; GLYPH&lt;11&gt;GLYPH&lt;12&gt; describe feedback from variations in the order parameters on the fluid flow [37]. The definition of the active stress GLYPH&lt;5&gt; act GLYPH&lt;11&gt;GLYPH&lt;12&gt; is such that any gradient in Q GLYPH&lt;11&gt;GLYPH&lt;12&gt; generates a flow field and drives the system at small length-scales, with strength determined by the magnitude of the activity GLYPH&lt;16&gt; . A positive (negative) GLYPH&lt;16&gt; corresponds to an extensile (contractile) material. The dipole flow-fields generated with this ansatz correspond to those of microswimmers - 'pushers' generate extensile stresses, 'pullers' contractile stresses [58]. The active stress continuously drives the system out of thermodynamic equilibrium. It will be of major importance in the following discussions.

## 2.2 Simulation setup

In Fig. 1, we show the geometry for which we will solve the equations. This simulation setup corresponds to the experimental configurations for cell monolayers in references [11, 12, 13, 14, 15]. It consists of a wider reservoir that feeds a narrower capillary. At t = 0, the active phase is restricted to the large reservoir.

In order to generate new active material in the reservoir, we locally increase GLYPH&lt;30&gt; to values higher than 1 0 at random positions. : For every point in the active phase in the reservoir a growth event takes place with probability r GLYPH&lt;30&gt; (1 GLYPH&lt;0&gt; GLYPH&lt;30&gt;=GLYPH&lt;30&gt; c ). This means that for a given time GLYPH&lt;28&gt; g , a source term GLYPH&lt;11&gt;GLYPH&lt;30&gt; is added to the r.h.s. of Eq. (3)

in a circle of radius rg around this site. Di GLYPH&lt;11&gt; usion and convection then lead to a spreading of the active phase, causing it to rise into the thinner capillary. This local implementation leads to additional flow in the reservoir as growth events generate dipolelike flow fields and a growth pressure through the isotropic part of GLYPH&lt;5&gt; el ; 1 . The logistic-growth-like saturation in the probability r GLYPH&lt;30&gt; (1 GLYPH&lt;0&gt; GLYPH&lt;30&gt;=GLYPH&lt;30&gt; c ) prevents unbounded growth of GLYPH&lt;30&gt; . We find that the details of this implementation or even the geometry of a reservoir are not important for the qualitative dynamics in the capillary (see ESI section A y ).

In this study, we use the following parameters: D = 0 08, : CLQ = 0 15, : KQ = 0 02, : K GLYPH&lt;30&gt; = 0 08, : GLYPH&lt;0&gt; GLYPH&lt;30&gt; = 0 2, : GLYPH&lt;0&gt; Q = 0 4, : GLYPH&lt;24&gt; = 0 7, : GLYPH&lt;17&gt; = 1 6, = r = 0 001, : GLYPH&lt;28&gt; g = 10000, rg = 5, GLYPH&lt;11&gt; = 0 01, : GLYPH&lt;30&gt; c = 1 2, : unless otherwise noted. We vary GLYPH&lt;16&gt; from 0 0020 to 0 0100 to : : investigate how invasion depends on the strength of the active driving. The free energy coe GLYPH&lt;14&gt; cients are chosen such that there is a well-defined interface between the active and the passive phase. GLYPH&lt;24&gt; = 0 7 leads to alignment of the director to an external : flow in a passive nematic. Parameter fitting of the continuum equations to physical active systems remains a topic of research; therefore, we consider a generic parameter set that has been shown to reproduce the flow vortex-lattice generated by a dense assembly of endothelial cells [59, 60], and the flow fields of dividing Madin-Darby Canine Kidney cells[47]. Independent samples for parameter sets are obtained by varying the initial condition randomly. Boundaries are no-slip with respect to the fluid and von-Neumann for GLYPH&lt;30&gt; and Q GLYPH&lt;11&gt;GLYPH&lt;12&gt; , meaning that there is no preferred anchoring of the director in the absence of activity. Variation of the boundary conditions can potentially lead to the emergence of additional exotic phases of invasion, but their systematic study lies beyond the scope of this work. For the reservoir, we have periodic boundary conditions in the direction perpendicular to the capillary axis.

In the following, we characterise the strength of the activity by the dimensionless activity number A = d = ( p K Q =GLYPH&lt;16&gt; ) = d = GLYPH&lt;3&gt; GLYPH&lt;16&gt; . The activity-induced length-scale GLYPH&lt;3&gt; GLYPH&lt;16&gt; = p KQ =GLYPH&lt;16&gt; emerges from the competition between the activity driving the dynamics and the elastic resistance against deformations in the director field [28, 61, 62], while the capillary width d imposes an upper limit for hydrodynamic interactions across the capillary. The activity number A relates these two length-scales [51, 52].

## 3 R esults

## 3.1 Distinct invasion regimes

For a first characterisation of the invasion behaviour of active matter into a capillary, we will focus on how the phenomenological changes in the structure and the dynamics of the interface between the active phase and the isotropic fluid depend on the dimensionless activity A . Surprisingly, we observe that this can be categorised into three di GLYPH&lt;11&gt; erent regimes separated by two well-defined crossovers under variations of this single quantity alone.

Figure 2A shows characteristic snapshots for the three di GLYPH&lt;11&gt; erent regimes. For small activities, below a certain threshold ( A GLYPH&lt;24&gt; 16), the system is characterised by a flat interface between the isotropic and active phase that advances steadily from the reservoir. We denote this as regime I. Due to the extensile ac-

15

d

Nc

Figure 2: Phenomenological regimes for varying activity number. (A) Snapshots from simulations showing typical configurations of the di GLYPH&lt;11&gt; erent regimes. From left to right: Regime I with flat interface, regime II with a deformed interface, and regime III where clusters (red circle) start to detach from the main active phase. Corresponding ESI Movies are flows0020.mp4, flows0035.mp4, and flows0060.mp4, respectively, see ESI section C y . (B) Quantitatively, the crossover from regime I to II is characterised by a sudden increase in ( h max GLYPH&lt;0&gt; h min ) = d from 0 to 1 (see Fig. 1 for the definition of h max and h min). By contrast, within regimes I and II, this quantity is nearly constant. In regime III it increases with activity. (C) While there is only one single coherent patch of active material in regimes I and II, the defining property of the crossover from II to III is the appearance of additional clusters. Nc denotes the number of clusters in addition to the bulk active phase. Black arrows in B and C mark values for which snapshots in A were taken. Errorbars depict 1.96 SEM (standard error of the mean). To estimate the SEM, simulations with di GLYPH&lt;11&gt; erent initial noise were performed for identical parameter configurations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00768v1-figure-162d8f55-3356-4b04-919e-802c3b351e6f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00768v1-figure-783764af-c3f0-4053-b1dc-58a8015199e4.jpg)

tivity of the particles (i.e. GLYPH&lt;16&gt; &gt; 0), the active stresses generate a preferential orientation of the director parallel to the interface, an e GLYPH&lt;11&gt; ect called 'active anchoring' [37]. For low A , the director field remains homogeneous throughout the nematic phase and the hydrodynamic instabilities are suppressed.

regularly created at the boundary and move throughout the bulk of the active phase (see ESI Movie flows0060.mp4 y ).

When approaching the crossover from regime I to regime II by increasing the activity, bend deformations of the director field arise in the bulk without influencing the shape of the interface (see ESI Movies flows0020.mp4 and flows0030.mp4 y ). The defining property of regime II is the crossover to a state where the interface is deformed to an S-shape while keeping a 90 GLYPH&lt;14&gt; contact angles at both walls. As discussed in more detail in the next section, within this regime the deformed interface initially advances with the leading point on one side of the capillary wall (see ESI Movie flows0035.mp4 y ). With a further increase in activity, we observe a periodic switching of the S-shape of the interface from one side of the capillary to the other side as the system approaches the third invasion regime.

In both regimes, I and II, the active phase remains as one single coherent phase which is connected to the reservoir. This changes with the crossover to regime III at A GLYPH&lt;24&gt; 20, where a higher activity leads to the dispatching of small clusters of active material from the main body of the active phase, that protrude deeper into the capillary. In addition, the active interface and the director field are strongly deformed and highly dynamic. Defects are

## 3.2 Detailed characterisation of invasion regimes and crossovers

To understand the connection between the di GLYPH&lt;11&gt; erent interface dynamics and the overall invasion process, we further characterise the defining features of the three di GLYPH&lt;11&gt; erent regimes. Quantitatively, the crossovers between the regimes can be read o GLYPH&lt;11&gt; from two observables: (i) ( h max GLYPH&lt;0&gt; h min ) = d , which characterises deformations of the isotropic-active interface in the capillary (see Fig. 1A for definition of h max and h min), and (ii) Nc which is the number of clusters of the active phase that are in the system in addition to the main body of the active phase connected to the reservoir. As the activity number is increased, the crossover from regime I to regime II is reflected by a jump in the mean value of h max GLYPH&lt;0&gt; h min from 0 to the width d of the capillary, or equivalently a jump of ( h max GLYPH&lt;0&gt; h min ) = d from 0 to 1 (Fig. 2B), which mirrors the S-shape of the interface in regime II. Within the parameter range studied here, we did not observe any variation in the positioning of the crossover from regime I to regime II under variation of D . The second crossover, from regime II to regime III, is also marked by a further jump in ( h max GLYPH&lt;0&gt; h min ) = d , but is even more evident in Nc which becomes finite (Fig. 2C), reflecting the emergence of additional smaller clusters that de-

Figure 3: Comparison of flows across regimes. (A) Root mean-squared velocity in units of d over the active time scale GLYPH&lt;28&gt; GLYPH&lt;16&gt; = GLYPH&lt;17&gt;=GLYPH&lt;16&gt; [56] plotted against A . Crossover from regime I to II coincides with the appearance of finite flows. Strength of flows increases inside regimes II and III with a small plateau at the crossover from regime II to regime III. Lower case letters indicate the values at which examples in (B-D) were taken. (B-D) Characteristic velocity and orientation fields near the interface. (B) For A = 15 5, the system is in regime I, the interface is flat : with no flows near it. (C) At A = 16 7, in regime II, the interface : is bent and the characteristic flow that turns at the interface is visible. (D) In regime III ( A = 20 4), the interface can take vari-: ous forms, + 1 2-defects (green dot) appear in the bulk and the = flow becomes highly dynamical and complex. ESI Movies corresponding to panels B-D are flows0030.mp4, flows0035.mp4, and flows0052.mp4, respectively, see ESI section C y .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00768v1-figure-97a6b959-9713-4974-87fa-6c377145ba1b.jpg)

tach from the main active phase. It is noteworthy that a similar quantitative behaviour is observed when varying the width of the channel (see ESI Fig. 2 in ESI section B y ), indicating that the activity number is the relevant dimensionless parameter in this setup.

## 3.2.1 Properties of regime II.

The crossover from regime I to II is also accompanied by sharp changes in the flow and director fields in the capillary. This not only sheds light on the underlying mechanism, but also has consequences for the overall speed of invasion as we will see later (section 3.3). Figures 3A,B show that below the first crossover, there are no flows in the capillary ( v rms = 0). For higher activities, with a non-zero ( h max GLYPH&lt;0&gt; h min ) = d , finite flows are generated (Figs. 3A,C). This is reminiscent of the well-established spontaneous flow transitions in confined active matter [49, 38], which arise due to the formation of hydrodynamic instabilities in active nematics above a certain threshold of activity. While in those cases the active force alone is responsible for the transition, in the present situation the growth dynamics in the reservoir induces long-range e GLYPH&lt;11&gt; ects which influence the position of the crossover (for a system without reservoir, see ESI section A y ). The advection of active material along with these flows causes the characteristic S-shape of the interface in this regime. Interestingly, v rms is increasing with higher activity numbers within regime II while ( h max GLYPH&lt;0&gt; h min ) = d GLYPH&lt;24&gt; 1 throughout.

Moreover, because of the impact of the extensile activity that keeps the director aligned with the interface, the increasing bend deformation looks similar to a backwards pointing + 1 2-= topological defect at the less advanced side of the interface (Fig. 3C; right panel). Corresponding to the characteristic selfmotility of + 1 2-topological defects [63, 51, 35], the interface = on this side continues to move backwards relative to the mean interface-position stretching out the interface until the defect detaches from the interface and moves into the active phase. Concomitantly with this defect absorption into the nematic, the interface re-attaches to the wall at a higher position. Repeatedly, this leads to a periodic shape-change in the interface (see Fig. 4A and the ESI Movie flows0035.mp4 y ).

As the activity is increased further, the initial spontaneous transition to flows is followed by the generation of counter-rotating vortices along the capillary. As the active phase grows and advances through the capillary, a lattice of vortices is formed behind the interface, reminiscent of vortex-lattices in a nongrowing confined active matter [52, 64], which appear when the activity-induced vorticity length scale becomes comparable to the capillary width. Interestingly, however, in this growing active matter, the orientation of the vortex closest to the interface also determines the parity of the interface. As a consequence, a change in the orientation of the most forward vortex leads to a switching at the interface, with the S-shape of the progressive front flipping from one side of the channel to the other with respect to the capillary axis (see Fig. 4B and the ESI Movie flows0046.mp4 y ).

Within regime II, the rms-velocity v rms grows approximately linearly with the activity number A and the flow field remains structured until transitioning to active turbulence in regime III (Fig. 3D).

## 3.2.2 Properties of regime III.

While the crossover to regime II was solely controlled by the hydrodynamic instability to spontaneous flow formation within the bulk of the active matter in the capillary, the crossover from regime II to III, where active clusters appear, is strongly influenced by the competition between activity and surface tension at the interface. This interfacial e GLYPH&lt;11&gt; ect can be clearly seen in a close-up view of the interface at the onset of cluster detachment (Fig. 4C): At a su GLYPH&lt;14&gt; ciently high activity number relatively long and thin layers of active matter with an approximately uniform director field are formed at the sides of the capillary wall. As the simulation time goes by, a small deformation at the tip of the layer develops and grows until it breaks away from the active layer, forming a detached cluster (see ESI Movie flows0052.mp4 y ). The formation of clusters means an increase in total interface length which is increasingly unfavourable for higher surface tension and thus requires higher active stresses to be generated. Therefore, the breakup of the thin layer can be understood as the destabilising e GLYPH&lt;11&gt; ect of the activity that leads to the interface deformation dominating over

Figure 4: Interface dynamics in regimes II and III. (A) Characteristic periodic motion in regime II ( A = 17 9). The : + 1 2-defect = like bend-deformation on the right wall moves backwards relative to the interface until it is pinched o GLYPH&lt;11&gt; and dissolved in the nematic bulk. Simultaneously the interface relaxes back to a less stretched state and the cycle begins anew. (B) Switching in the S-shape for higher activities in regime II ( A = 19 2). When ap-: proaching the second crossover, vortices emerge in the capillary (panels on the right). When the most forward vortex changes its rotation direction, the director follows and bend deformations as well as interface shape adapt accordingly. (C) Time series of a cluster-formation event in regime III ( A = 20 4). The nematic : layer at the wall becomes unstable until finally a cluster detaches from the main active phase. In each of (A-C) the simulation time between all sequential snapshots is equal and is measured in units of the active time scale GLYPH&lt;28&gt; GLYPH&lt;16&gt; = GLYPH&lt;17&gt;=GLYPH&lt;16&gt; [65], A: GLYPH&lt;14&gt; t = 480 GLYPH&lt;28&gt; GLYPH&lt;16&gt; , B: GLYPH&lt;14&gt; t = 2760 GLYPH&lt;28&gt; GLYPH&lt;16&gt; , C: GLYPH&lt;14&gt; t = 1248 GLYPH&lt;28&gt; GLYPH&lt;16&gt; . ESI Movies corresponding to panels B and C are flows0046.mp4 and flows0052.mp4, respectively, see ESI section C y .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00768v1-figure-7012b5ef-eaeb-4679-abb9-c9db8acfd488.jpg)

the stabilising e GLYPH&lt;11&gt; ect of the surface tension, working to keep the interface straight. However, especially in this complex setup, it is not possible to decide to what degree this e GLYPH&lt;11&gt; ect is purely interfacial and whether dynamics in the bulk are important. The simultaneous appearance of + 1 2 topological defects in the ne-= matic phase (see ESI Fig. 3 y ) hints that changes in the bulk dynamics are also involved in the crossover as flow dynamics become more reminiscent of turbulent dynamics and higher active stresses are generated. Interestingly, directly after the crossover the number of additional clusters peaks. With a further increase in A the clusters are still present, however the number of clusters goes down. A close look at the interface reveals that at very high activities the active material forms extremely long and thin layers on the capillary walls and the detached clusters quickly reattach to the main body of the active nematic, leading to a smaller number of clusters on average.

## 3.3 Invasion capability

In the context of our original research question we now investigate the influence of the di GLYPH&lt;11&gt; erent regimes on the capability of the active system to claim new territories. To this end, we show in Fig. 5A the total amount of active material in the capillary at the end of the simulation (9 5 10 : 6 simulation steps) GLYPH&lt;8&gt; T : = R GLYPH&lt;30&gt; ( r ; T ) d 2 r = d 2 , which we denote as the invasion index. In the absence of flows (regime I, when the process is purely di GLYPH&lt;11&gt; usive), the invasion index GLYPH&lt;8&gt; T is almost independent of A and it is significantly lower than in regimes II and III, where activity-induced flows enhance the advective transport of active material from the reservoir to the capillary. After the crossover from regime I to II, GLYPH&lt;8&gt; T follows a linear dependence on A similar to the linear increase of vrms in this region. The crossover from regime II to III does not alter this quantity in a qualitative way, i.e. the additional clusters have no influence on the total amount of active material that invades the capillary. However, they can be shown to invade about 0 5 to 1 0 capillary widths deeper into : : the capillary than the tip of the coherent active phase connected to the reservoir. This is best illustrated by plotting the relative di GLYPH&lt;11&gt; erence of the maximum invasion height h inv (see Fig. 1A) and the maximum height of the interface h max: ( h inv GLYPH&lt;0&gt; h max) = d (Fig. 5B). Together, these results indicate that the flows in the capillary and the interfacial dynamics control the rate at which active material enters the capillary and thus determine the invasion speed.

## 4 C onclusions

In this work, we have presented a two-phase computational framework for studying the growth of active matter within an isotropic fluid medium. Our approach accounts for hydrodynamic e GLYPH&lt;11&gt; ects and the orientational dynamics of active particles, and is able to reproduce emergent active phenomena including spontaneous flow generation and active turbulence. Motivated by the prevalence of active matter growing within confined spaces and constrictions, we investigated the dynamics in capillaries of di GLYPH&lt;11&gt; erent sizes, focusing on the combined e GLYPH&lt;11&gt; ects of growth and the active stresses that are continuously exerted by active particles on their surroundings.

We find three distinct invasion patterns as the activity is varied: (i) at small activities invasion is completely controlled by the growth dynamics and the interface between the active phase and the surroundings remains stable with no, or slight, deformations. (ii) Increasing activity beyond a given threshold results in spontaneous flow generation which significantly enhances the invasion within the capillary and is accompanied by deformations of the interface. (iii) At yet higher strengths of activity, a second crossover appears, where active clusters begin to detach from the main body of the invading active phase, further enhancing the invasion of the surrounding space.

Our analysis shows that the crossovers between the various invasion regimes are controlled by di GLYPH&lt;11&gt; erent mechanisms. The

Figure 5: Propensity of the system to invade the capillary. (A) Invasion index GLYPH&lt;8&gt; T : = R GLYPH&lt;30&gt; ( r ; T ) d 2 r = d 2 defined as the amount of active material inside the capillary after 9 5 10 : 6 simulation steps plotted versus activity number A . In regime I, the invasion shows little dependence on A until the appearance of flows at the crossover from regime I to II leads to a sharp increase. (B) Relative di GLYPH&lt;11&gt; erence of the maximum invasion height and the maximum height of the interface ( h inv GLYPH&lt;0&gt; h max) = d against A . In regimes I and II, this di GLYPH&lt;11&gt; erence is trivially zero. In regime III it shows an increase of up to one capillary-width meaning the clusters protrude significantly deeper into the capillary than the main body of the growing active matter.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00768v1-figure-6d50335c-e648-4e9e-8245-bbb5cb079cc6.jpg)

crossover between regimes I and II is governed by the wellestablished spontaneous flow formation in confined active matter [49, 50], which is due to a hydrodynamic instability of the active nematic. We show that the onset of flow within the bulk is accompanied by the deformation of the interface and can even result in a periodic switching of the active protrusions from one side of the capillary to the other, as flow vortices start to form behind the advancing front. When activity becomes stronger, turbulent flows and defects are observed in the bulk of the active material and, dependent on its surface tension, the system also becomes able to rip apart the active-to-passive interface marking the third invasion regime.

Our results thus highlight the significant di GLYPH&lt;11&gt; erences in invasion patterns when the activity of growing matter is accounted for. From a biological perspective, this could be linked to the invasive behaviour of biological matter such as growing colonies of cells or bacteria. It suggests that qualitative changes in the spread of these organisms into confined spaces can be caused by changes in the availability and conversion rate of chemical energy into mechanical stresses. One could conjecture that more invasive cell lines or bacterial strands could regulate their invasion of pores or cavities by tuning their strength of motion. For example, inducing the formation of additional clusters might be beneficial for cell-lines that aim at aggressive invasion of surrounding tissue while, for bacteria, it might be beneficial to remain coherent in order to profit from the benefits of cooperative behaviour.

Considering possible experimental realisations, our model predictions apply to active growing systems, where hydrodynamic interactions and nematic orientational order play an important role. Such nematic ordering has been reported in various biological systems including bacterial colonies [48, 66], cultures of amoeboid cells [67], fibroblast cells [46, 16], human bronchial cells [68], neural progenitor stem cells [69], and Madin-Darby Kanine Kidney (MDCK) epithelial monolayers [35, 70]. Investigating the degree to which experiment and theory can be matched by the adaptation of parameters and boundary conditions in the model or the reasons for deviations is a powerful way of unraveling the role of mechanics in determining the behaviour of such active biological matter. It would also be interesting to compare to models that resolve individual particles such as models of self-propelled particles [36], phase field approaches [71, 72, 73], or cellular Potts models [74]. In contrast to these, our model by construction cannot resolve phenomena or implement mechanisms that act on the single particle-level, for example dynamics of cell size or shape, but due to its generality could contribute to a deeper understanding of the generic fundamental mechanisms that underlie the invasion of microscopic biological systems.

Finally, several improvements can be envisaged for the current model. In the context of growing tissues, studies have highlighted the importance of leader cells at the progression front [75, 76]. Within particle-based models such an e GLYPH&lt;11&gt; ect is shown to be captured by introducing a curvature-dependent motility for the cells at the interface [77, 36]. Within our proposed framework, a similar e GLYPH&lt;11&gt; ect could be modelled by introducing a curvature dependence to the activity coe GLYPH&lt;14&gt; cient. In addition, the natural environment of active material is often more complicated than the isotropic fluid considered here and active invasion happens through viscoelastic media. For example, cells invade the interconnected networks of collagen matrices and bacteria secrete their own extracellular matrices. The continuum framework presented here could be extended to include viscoelasticity by introducing additional order parameters representing polymer conformation.

## C onflicts of interest

There are no conflicts to declare.

## A cknowledgements

E.F. and F.K. would like to thank Sophia Scha GLYPH&lt;11&gt; er, Matthias Zorn, and Joachim Rädler for stimulating discussions. A.D, F.K., and J.M.Y. would like to thank Kristian Thijssen for helpful discussions. This research was supported by the German Excellence Initiative via the program 'Nanosystems Initiative Munich' (NIM) and the Deutsche Forschungsgemeinschaft (DFG) via project B03 within the Collaborative Research Center (SFB 1032) 'nanoagents'. R.M. was supported by grant P2EZP2\_165261 of the Swiss National Science Foundation. A.D. was supported by the Royal Commission for Exhibition of 1851 Research Fellowship.

The preprint template has been taken from https://github.com/ brenhinkeller/preprint-template.tex/

## R eferences

- [1] A. G. Clark and D. M. Vignjevic. Modes of cancer cell invasion and the role of the microenvironment. Curr. Opin. Cell Biol. , 36:13-22, 2015.
- [2] Ralitza Staneva, Fatima El Marjou, Jorge Barbazan, Denis Krndija, Sophie Richon, Andrew Clark, and Danijela Matic Vignjevic. Cancer cells in the tumor core exhibit spatially coordinated migration patterns. J. Cell Sci. , 132(6):jcs.220277, 2019. ISSN 0021-9533. doi: 10.1242 jcs.220277. / URL http://jcs.biologists. org/content/early/2019/02/11/jcs.220277 .
- [3] Yilin Wu. Collective motion of bacteria in two dimensions. Quant. Biol. , 3:199-205, 2015.
- [4] Catherine Even, Christian Marlière, Jean-Marc Ghigo, Jean-Marc Allain, Alba Marcellan, and Eric Raspaud. Recent advances in studying single bacteria and biofilm mechanics. Adv. Colloid Interface Sci. , 247:573-588, 2017.
- [5] Raimo Hartmann, Praveen K Singh, Philip Pearce, Rachel Mok, Boya Song, Francisco Díaz-Pascual, Jörn Dunkel, and Knut Drescher. Emergence of three-dimensional order and structure in growing biofilms. Nat. Phys. , 15:251-256, 2018.
- [6] Chunping Yang, Hong Chen, Guangming Zeng, Guanlong Yu, and Shenglian Luo. Biomass accumulation and control strategies in gas biofiltration. Biotechnol. Adv. , 28(4):531 -540, 2010. ISSN 07349750. doi: https: // doi.org 10.1016 j.biotechadv.2010. / / 04.002. URL http://www.sciencedirect.com/ science/article/pii/S0734975010000418 .
- [7] Pavlo G. Gritsenko, Olga Ilina, and Peter Friedl. Interstitial guidance of cancer invasion. The Journal of Pathology , 226(2):185-199, 2011. doi: 10.1002 path. / 3031. URL https://onlinelibrary.wiley.com/ doi/abs/10.1002/path.3031 .
- [8] Bettina Weigelin, Gert-Jan Bakker, and Peter Friedl. Intravital third harmonic generation microscopy of collective melanoma cell invasion. IntraVital , 1(1):32-43, 2012. doi: 10.4161 intv.21223. URL / https://doi.org/10.4161/ intv.21223 .
- [9] Luanne Hall-Stoodley, J. William Costerton, and Paul Stoodley. Bacterial biofilms: from the natural environment to infectious diseases. Nat. Rev. Microbiol. , 2:95108, Feb 2004. URL http://dx.doi.org/10.1038/ nrmicro821 . Review Article.
- [10] Jacinta C Conrad and Ryan Poling-Skutvik. Confined flow: Consequences and implications for bacteria and biofilms. Annu. Rev. Chem. Biomol. Eng. , 9:175-200, 2018.
- [11] Sri Ram Krishna Vedula, Man Chun Leong, Tan Lei Lai, Pascal Hersen, Alexandre J. Kabla, Chwee Teck Lim, and Benoît Ladoux. Emerging modes of collective cell migration induced by geometrical constraints. Proc. Natl. Acad. Sci. U. S. A. , 109(32):12974-12979, 2012. ISSN 0027-8424. doi: 10.1073 pnas.1119313109. URL / http://www.pnas.org/content/109/32/12974 .
- [12] Anna Kristina Marel, Matthias Zorn, Christoph Klingner, Roland Wedlich-Söldner, Erwin Frey, and Joachim O.
- Rädler. Flow and di GLYPH&lt;11&gt; usion in channel-guided cell migration. Biophys. J. , 107(5):1054-1064, 2014. ISSN 15420086. doi: 10.1016 j.bpj.2014.07.017. /
- [13] Anna-Kristina Marel, Nils Podewitz, Matthias Zorn, Joachim Oskar Rädler, and Jens Elgeti. Alignment of cell division axes in directed epithelial cell migration. New J. Phys. , 16(11):115005, 2014. URL http://stacks. iop.org/1367-2630/16/i=11/a=115005 .
- [14] Yongliang Yang, Nima Jamilpour, Baoyin Yao, Zachary S. Dean, Reza Riahi, and Pak Kin Wong. Probing leader cells in endothelial collective migration by plasma lithography geometric confinement. Sci. Rep. , 6:srep22707, Mar 2016. URL http://dx.doi.org/10.1038/srep22707 . Article.
- [15] Estelle Gauquelin, Sham Tlili, Cyprien Gay, Grégoire Peyret, René-Marc Mège, Marc A. Fardin, and Benoît Ladoux. Influence of proliferation on the motions of epithelial monolayers invading adherent strips. Soft Matter , 15:2798-2810, 2019. doi: 10.1039 C9SM00105K. URL / http://dx.doi.org/10.1039/C9SM00105K .
- [16] Guillaume Duclos, Christoph Erlenkämper, Jean-François Joanny, and Pascal Silberzan. Topological defects in confined populations of spindle-shaped cells. Nat. Phys. , 13: 58-62, Sep 2017. URL http://dx.doi.org/10.1038/ nphys3876 .
- [17] G. Duclos, C. Blanch-Mercader, V. Yashunsky, G. Salbreux, J.-F. Joanny, J. Prost, and P. Silberzan. Spontaneous shear flow in confined cellular nematics. Nat. Phys. , 14(7):728-732, 2018. ISSN 1745-2481. doi: 10.1038 s41567-018-0099-7. URL / https://doi.org/ 10.1038/s41567-018-0099-7 .
- [18] H Wioland, E Lushi, and R E Goldstein. Directed collective motion of bacteria under channel confinement. New J. Phys. , 18(7):075002, 2016. URL http://stacks.iop. org/1367-2630/18/i=7/a=075002 .
- [19] Jacinta C. Conrad and Ryan Poling-Skutvik. Confined flow: Consequences and implications for bacteria and biofilms. Annu. Rev. Chem. Biomol. Eng. , 9(1):175-200, 2018. doi: 10.1146 annurev-chembioeng-060817-084006. / URL https://doi.org/10.1146/ annurev-chembioeng-060817-084006 . PMID: 29561646.
- [20] Kazuya Suzuki, Makito Miyazaki, Jun Takagi, Takeshi Itabashi, and Shin'ichi Ishiwata. Spatial confinement of active microtubule networks induces large-scale rotational cytoplasmic flow. Proc. Natl. Acad. Sci. U. S. A. , 114(11): 2922-2927, 2017. ISSN 0027-8424. doi: 10.1073 pnas. / 1616001114. URL http://www.pnas.org/content/ 114/11/2922 .
- [21] Kun-Ta Wu, Jean Bernard Hishamunda, Daniel T. N. Chen, Stephen J. DeCamp, Ya-Wen Chang, Alberto Fernández-Nieves, Seth Fraden, and Zvonimir Dogic. Transition from turbulent to coherent flows in confined three-dimensional active fluids. Science , 355(6331): eaal1979, 2017. ISSN 0036-8075. doi: 10.1126 science. / aal1979. URL http://science.sciencemag.org/ content/355/6331/eaal1979 .

- [22] Hugo Wioland, Francis G. Woodhouse, Jörn Dunkel, John O. Kessler, and Raymond E. Goldstein. Confinement stabilizes a bacterial suspension into a spiral vortex. Phys. Rev. Lett. , 110:268102, Jun 2013. doi: 10.1103 PhysRevLett.110.268102. URL / https://link. aps.org/doi/10.1103/PhysRevLett.110.268102 .
- [23] Jérôme Hardoüin, Rian Hughes, Amin Doostmohammadi, Justine Laurent, Teresa Lopez-Leon, Julia M Yeomans, Jordi Ignés-Mullol, and Francesc Sagués. Reconfigurable flows and defect landscape of confined active nematics. page arXiv:1903.01787, 2019.
- [24] Achini Opathalage, Michael M. Norton, Michael P. N. Juniper, Blake Langeslay, S. Ali Aghvami, Seth Fraden, and Zvonimir Dogic. Self-organized dynamics and the transition to turbulence of confined active nematics. Proc. Natl. Acad. Sci. U. S. A. , 116(11):4788-4797, 2019. ISSN 0027-8424. doi: 10.1073 pnas.1816733116. URL / https: //www.pnas.org/content/116/11/4788 .
- [25] Marvin Brun-Cosme-Bruny, Eric Bertin, Benoît Coasne, Philippe Peyla, and Salima Rafaï. E ective di GLYPH&lt;11&gt; GLYPH&lt;11&gt; usivity of microswimmers in a crowded environment. J. Chem. Phys. , 150(10):104901, 2019. doi: 10.1063 1.5081507. / URL https://doi.org/10.1063/1.5081507 .
- [26] Roberto Alonso-Matilla, Brato Chakrabarti, and David Saintillan. Transport and dispersion of active particles in periodic porous media. Phys. Rev. Fluids , 4: 043101, Apr 2019. doi: 10.1103 PhysRevFluids.4. / 043101. URL https://link.aps.org/doi/10.1103/ PhysRevFluids.4.043101 .
- [27] Henricus H. Wensink, Jörn Dunkel, Sebastian Heidenreich, Knut Drescher, Raymond E. Goldstein, Hartmut Löwen, and Julia M. Yeomans. Meso-scale turbulence in living fluids. Proc. Natl. Acad. Sci. U. S. A. , 109(36): 14308-14313, 2012. ISSN 0027-8424. doi: 10.1073 pnas. / 1202032109. URL https://www.pnas.org/content/ 109/36/14308 .
- [28] Sumesh P. Thampi, Ramin Golestanian, and Julia M. Yeomans. Vorticity, defects and correlations in active turbulence. Philos. Trans. R. Soc., A , 372 (2029):20130366, 2014. doi: 10.1098 rsta.2013.0366. / URL https://royalsocietypublishing.org/doi/ abs/10.1098/rsta.2013.0366 .
- [29] Luca Giomi. Geometry and topology of turbulence in active nematics. Phys. Rev. X , 5:031003, Jul 2015. doi: 10.1103 PhysRevX.5.031003. URL / https://link.aps. org/doi/10.1103/PhysRevX.5.031003 .
- [30] Martin James, Wouter J. T. Bos, and Michael Wilczek. Turbulence and turbulent pattern formation in a minimal model for active fluids. Phys. Rev. Fluids , 3: 061101, Jun 2018. doi: 10.1103 PhysRevFluids.3. / 061101. URL https://link.aps.org/doi/10.1103/ PhysRevFluids.3.061101 .
- [31] Vasil Bratanov, Frank Jenko, and Erwin Frey. New class of turbulence in active fluids. Proc. Natl. Acad. Sci. U. S. A. , 112(49):15048-15053, 2015. ISSN 0027-8424. doi: 10. 1073 pnas.1509304112. URL / https://www.pnas.org/ content/112/49/15048 .
- [32] J. Elgeti, M. E. Cates, and D. Marenduzzo. Defect hydrodynamics in 2d polar active fluids. Soft Matter , 7: 3177-3185, 2011. doi: 10.1039 C0SM01097A. / URL http://dx.doi.org/10.1039/C0SM01097A .
- [33] Luca Giomi, Mark J. Bowick, Xu Ma, and M. Cristina Marchetti. Defect annihilation and proliferation in active nematics. Phys. Rev. Lett. , 110:228101, May 2013. doi: 10.1103 PhysRevLett.110.228101. URL / https://link. aps.org/doi/10.1103/PhysRevLett.110.228101 .
- [34] Luca Giomi, Mark J Bowick, Prashant Mishra, Rastko Sknepnek, and M Cristina Marchetti. Defect dynamics in active nematics. Philos. Trans. R. Soc., A , 372 (2029):20130365, 2014. doi: 10.1098 rsta.2013.0365. / URL https://royalsocietypublishing.org/doi/ abs/10.1098/rsta.2013.0365 .
- [35] Thuan Beng Saw, Amin Doostmohammadi, Vincent Nier, Leyla Kocgozlu, Sumesh Thampi, Yusuke Toyama, Philippe Marcq, Chwee Teck Lim, Julia M. Yeomans, and Benoit Ladoux. Topological defects in epithelia govern cell death and extrusion. Nature , 544:212-216, Apr 2017. URL http://dx.doi.org/10.1038/nature21718 .
- [36] Victoria Tarle, Estelle Gauquelin, S R K Vedula, Joseph D'Alessandro, C T Lim, Benoit Ladoux, and Nir S Gov. Modeling collective cell migration in geometric confinement. Phys. Biol. , 14(3):035001, 2017. URL http: //stacks.iop.org/1478-3975/14/i=3/a=035001 .
- [37] Matthew L. Blow, Sumesh P. Thampi, and Julia M. Yeomans. Biphasic, lyotropic, active nematics. Phys. Rev. Lett. , 113:248303, Dec 2014. doi: 10.1103 PhysRevLett.113. / 248303. URL http://link.aps.org/doi/10.1103/ PhysRevLett.113.248303 .
- [38] D. Marenduzzo, E. Orlandini, M. E. Cates, and J. M. Yeomans. Steady-state hydrodynamic instabilities of active liquid crystals: Hybrid lattice boltzmann simulations. Phys. Rev. E , 76:031921, Sep 2007. doi: 10.1103 PhysRevE.76. / 031921. URL https://link.aps.org/doi/10.1103/ PhysRevE.76.031921 .
- [39] R. Aditi Simha and Sriram Ramaswamy. Hydrodynamic fluctuations and instabilities in ordered suspensions of self-propelled particles. Phys. Rev. Lett. , 89:058101, Jul 2002. doi: 10.1103 PhysRevLett.89. / 058101. URL https://link.aps.org/doi/10.1103/ PhysRevLett.89.058101 .
- [40] Sumesh P. Thampi, Ramin Golestanian, and Julia M. Yeomans. Velocity correlations in an active nematic. Phys. Rev. Lett. , 111:118101, Sep 2013. doi: 10.1103 PhysRevLett. / 111.118101. URL https://link.aps.org/doi/10. 1103/PhysRevLett.111.118101 .
- [41] Stephen J. DeCamp, Gabriel S. Redner, Aparna Baskaran, Michael F. Hagan, and Zvonimir Dogic. Orientational order of motile defects in active nematics. Nat. Mater. , 14: 1110-1115, Aug 2015. URL http://dx.doi.org/10. 1038/nmat4387 .
- [42] Dario Cortese, Jens Eggers, and Tanniemola B. Liverpool. Pair creation, motion, and annihilation of topological defects in two-dimensional nematic liquid crystals. Phys. Rev. E , 97:022704, Feb 2018. doi: 10.1103 PhysRevE.97. /

- 022704. URL https://link.aps.org/doi/10.1103/ PhysRevE.97.022704 .
- [43] Amin Doostmohammadi, Jordi Ignés-Mullol, Julia M. Yeomans, and Francesc Sagués. Active nematics. Nat. Commun. , 9(1):3246, 2018. ISSN 2041-1723. doi: 10.1038 s41467-018-05666-8. URL / https://doi.org/ 10.1038/s41467-018-05666-8 .
- [44] Tim Sanchez, Daniel T. N. Chen, Stephen J. DeCamp, Michael Heymann, and Zvonimir Dogic. Spontaneous motion in hierarchically assembled active matter. Nature , 491:431-434, Nov 2012. URL http://dx.doi.org/10. 1038/nature11591 .
- [45] Thomas Bittig, Ortrud Wartlick, Anna Kicheva, Marcos González-Gaitán, and Frank Jülicher. Dynamics of anisotropic tissue growth. New J. Phys. , 10(6):063001, 2008. URL http://stacks.iop.org/1367-2630/10/ i=6/a=063001 .
- [46] G. Duclos, S. Garcia, H. G. Yevick, and P. Silberzan. Perfect nematic order in confined monolayers of spindleshaped cells. Soft Matter , 10:2346-2353, 2014. doi: 10.1039 C3SM52323C. URL / http://dx.doi.org/10. 1039/C3SM52323C .
- [47] Amin Doostmohammadi, Sumesh P. Thampi, Thuan B. Saw, Chwee T. Lim, Benoit Ladoux, and Julia M. Yeomans. Celebrating soft matter's 10th anniversary: Cell division: a source of active stress in cellular monolayers. Soft Matter , 11:7328-7336, 2015. doi: 10.1039 C5SM01382H. URL / http://dx.doi.org/10.1039/C5SM01382H .
- [48] Dmitri Volfson, Scott Cookson, Je GLYPH&lt;11&gt; Hasty, and Lev S. Tsimring. Biomechanical ordering of dense cell populations. Proc. Natl. Acad. Sci. U. S. A. , 105(40):1534615351, 2008. ISSN 0027-8424. doi: 10.1073 pnas. / 0706805105. URL http://www.pnas.org/content/ 105/40/15346 .
- [49] R. Voituriez, J. F. Joanny, and J. Prost. Spontaneous flow transition in active polar gels. EPL , 70(3):404, 2005. URL http://stacks.iop.org/0295-5075/70/ i=3/a=404 .
- [50] S. A. Edwards and J. M. Yeomans. Spontaneous flow states in active nematics: A unified picture. EPL , 85(1):18008, 2009. URL http://stacks.iop.org/0295-5075/85/ i=1/a=18008 .
- [51] Tyler N. Shendruk, Amin Doostmohammadi, Kristian Thijssen, and Julia M. Yeomans. Dancing disclinations in confined active nematics. Soft Matter , 13:3853-3862, 2017. doi: 10.1039 C6SM02310J. URL / http://dx.doi.org/ 10.1039/C6SM02310J .
- [52] Amin Doostmohammadi, Tyler N. Shendruk, Kristian Thijssen, and Julia M. Yeomans. Onset of mesoscale turbulence in active nematics. Nat. Commun. , 8:15326, 2017. URL http://dx.doi.org/10.1038/ ncomms15326 . Article.
- [53] Michael M. Norton, Arvind Baskaran, Achini Opathalage, Blake Langeslay, Seth Fraden, Aparna Baskaran, and Michael F. Hagan. Insensitivity of active nematic liquid crystal dynamics to topological constraints. Phys. Rev. E , 97:012702, Jan 2018. doi: 10.1103 PhysRevE.97. /
- 012702. URL https://link.aps.org/doi/10.1103/ PhysRevE.97.012702 .
- [54] Pierre G. de Gennes and J. Prost. The Physics of Liquid Crystals . Oxford University Press, Oxford, 1995.
- [55] Paul M. Chaikin and Taylor C. Lubensky. Principles of Condensed Matter Physics . Cambridge University Press, Cambridge, 2000.
- [56] E. Orlandini, M. R. Swift, and J. M. Yeomans. A lattice boltzmann model of binary-fluid mixtures. EPL , 32(6):463, 1995. URL http://stacks.iop.org/0295-5075/32/ i=6/a=001 .
- [57] R. G. Larson. The structure and rheology of complex fluids . Oxford University Press, New York (N.Y.), 1999.
- [58] Sriram Ramaswamy. The mechanics and statistics of active matter. Annu. Rev. Condens. Matter Phys. , 1:323-345, 2010.
- [59] Ninna S. Rossen, Jens M. Tarp, Joachim Mathiesen, Mogens H. Jensen, and Lene B. Oddershede. Long-range ordered vorticity patterns in living tissue induced by cell division. Nat. Commun. , 5:5720, Dec 2014. URL http://dx.doi.org/10.1038/ncomms6720 . Article.
- [60] Amin Doostmohammadi, Michael F Adamer, Sumesh P Thampi, and Julia M Yeomans. Stabilization of active matter by flow-vortex lattices and defect ordering. Nat. Commun. , 7:10557, 2016.
- [61] Ewan J. Hemingway, Prashant Mishra, M. Cristina Marchetti, and Suzanne M. Fielding. Correlation lengths in hydrodynamic models of active nematics. Soft Matter , 12:7943-7952, 2016. doi: 10.1039 C6SM00812G. URL / http://dx.doi.org/10.1039/C6SM00812G .
- [62] Pau Guillamat, Jordi Ignés-Mullol, and Francesc Sagués. Taming active turbulence with patterned soft interfaces. Nat. Commun. , 8:564, 2017.
- [63] Elias Putzig and Aparna Baskaran. Phase separation and emergent structures in an active nematic fluid. Phys. Rev. E , 90:042304, Oct 2014. doi: 10.1103 PhysRevE.90. / 042304. URL https://link.aps.org/doi/10.1103/ PhysRevE.90.042304 .
- [64] Maxime Theillard, Roberto Alonso-Matilla, and David Saintillan. Geometric control of active collective motion. Soft Matter , 13:363-375, 2017.
- [65] Sumesh P. Thampi, Amin Doostmohammadi, Ramin Golestanian, and Julia M. Yeomans. Intrinsic free energy in active nematics. EPL , 112(2):28004, 2015. URL http: //stacks.iop.org/0295-5075/112/i=2/a=28004 .
- [66] He Li, Xia-qing Shi, Mingji Huang, Xiao Chen, Minfeng Xiao, Chenli Liu, Hugues Chaté, and H. P. Zhang. Datadriven quantitative modeling of bacterial active nematics. Proc. Natl. Acad. Sci. U. S. A. , 116(3):777-785, 2019. ISSN 0027-8424. doi: 10.1073 pnas.1812570116. URL / https://www.pnas.org/content/116/3/777 .
- [67] H Gruler, U Dewald, and M Eberhardt. Nematic liquid crystals formed by living amoeboid cells. Euro. Phys. J. B , 11:187-192, 1999.
- [68] C Blanch-Mercader, V Yashunsky, S Garcia, G Duclos, L Giomi, and P Silberzan. Turbulent dynamics of epithelial cell cultures. Phys. Rev. Lett. , 120:208101, 2018.

- [69] Kyogo Kawaguchi, Ryoichiro Kageyama, and Masaki Sano. Topological defects control collective dynamics in neural progenitor cell cultures. Nature , 545:327331, Apr 2017. URL http://dx.doi.org/10.1038/ nature22321 .
- [70] Thuan Beng Saw, Wang Xi, Benoit Ladoux, and Chwee Teck Lim. Biological tissues as active nematic liquid crystals. Adv. Mater. , 30(47):1802579, 2018. doi: 10. 1002 adma.201802579. URL / https://onlinelibrary. wiley.com/doi/abs/10.1002/adma.201802579 .
- [71] Danying Shao, Wouter-Jan Rappel, and Herbert Levine. Computational model for cell morphodynamics. Phys. Rev. Lett. , 105:108104, Sep 2010. doi: 10.1103 PhysRevLett. / 105.108104. URL https://link.aps.org/doi/10. 1103/PhysRevLett.105.108104 .
- [72] Igor S. Aranson. Physical Models of Cell Motility . Springer International Publishing, Switzerland, 2016.
- [73] Romain Mueller, Julia M. Yeomans, and Amin Doostmohammadi. Emergence of active nematic behavior in monolayers of isotropic cells. Phys. Rev. Lett. , 122: 048004, Feb 2019. doi: 10.1103 PhysRevLett.122. / 048004. URL https://link.aps.org/doi/10.1103/ PhysRevLett.122.048004 .
- [74] Florian Thuero GLYPH&lt;11&gt; , Andriy Goychuk, Matthias Reiter, and Erwin Frey. Bridging the gap between single cell migration and collective dynamics. page bioRxiv 548677, 2019. / doi: 10.1101 548677. / URL https://www.biorxiv. org/content/early/2019/02/13/548677 .
- [75] Mathieu Poujade, Erwan Grasland-Mongrain, A Hertzog, J Jouanneau, Philippe Chavrier, Benoît Ladoux, Axel Buguin, and Pascal Silberzan. Collective migration of an epithelial monolayer in response to a model wound. Proc. Natl. Acad. Sci. U. S. A. , 104:15988-15993, 2007.
- [76] Xavier Trepat, Michael R Wasserman, Thomas E Angelini, Emil Millet, David A Weitz, James P Butler, and Je GLYPH&lt;11&gt; rey J Fredberg. Physical forces during collective cell migration. Nat. Phys. , 5(6):426, 2009.
- [77] Shirley Mark, Roie Shlomovitz, Nir S Gov, Mathieu Poujade, Erwan Grasland-Mongrain, and Pascal Silberzan. Physical model of the dynamic instability in an expanding cell culture. Biophys. J. , 98:361-370, 2010.

## S upplementary M aterials

## A R ole of reservoir and growth

In this section we examine the importance of the cell division and the presence of a reservoir on the dynamics reported in the main text. ESI Figure 1 shows simulation results for two systems that di GLYPH&lt;11&gt; er from the one considered in the main text, by the absence of a wider reservoir, and by the absence of cell division, respectively. We see that growth is an essential factor to create the phenomena reported in the main text as in the absence of growth the behaviour is qualitatively di GLYPH&lt;11&gt; erent. If the reservoir is absent, meaning that there is only a capillary of uniform width where growth takes place in the lower region, we in contrast observe the first crossover at the same values for the activity A (ESI Fig. 1A), together with the same change in the invasion speed (ESI Fig. 1B). The second crossover is however significantly shifted to higher values of A (see ESI Fig. 1C,D), indicating that the reservoir has long-range e GLYPH&lt;11&gt; ects on the dynamics in the capillary. These, however, do not change the qualitative picture.

## B U niversality for different channel widths

As described in the main text, the activity number A characterises the ratio of the channel width d to the intrinsic length scale of the active phase. In ESI Figure 2 the observables ( h max GLYPH&lt;0&gt; h min ) = d and Nc are plotted against A for varying channel width d . Both crossovers happen at the same values for A , independent of d . In addition, both observables are of the same order of magnitude and within regime II ( h max GLYPH&lt;0&gt; h min ) = d GLYPH&lt;24&gt; 1. These results indicate that the activity number is the relevant dimensionless parameter in this setup and that ( h max GLYPH&lt;0&gt; h min ) = d is a meaningful observable.

## C D escription of M ovies

All movies show the time-evolution of the velocity field on the left hand side and the orientation field on the right hand side, in the same way as panels B-D in Figure 3. We choose di GLYPH&lt;11&gt; erent values of the activity GLYPH&lt;16&gt; to display the qualitatively di GLYPH&lt;11&gt; erent phenomena. There are the following files:

Higher quality movies can be found on https: // doi.org 10.1039 c9sm01210a. / /

flows0020.mp4 GLYPH&lt;16&gt; = 0 0020. System is in regime I, director ho-: mogeneous and parallel to the interface. Corresponds to Figure 2A, leftmost panel.

flows0030.mp4 GLYPH&lt;16&gt; = 0 0030. : System is in regime I, but the director starts to show deviations from the perfectly homogeneous configuration, especially close to the reservoir. Corresponds to Figure 3B.

flows0035.mp4 GLYPH&lt;16&gt; = 0 0035. : System is in regime 2, showing characteristic flow- and director-fields. At later times the periodic dynamics described in section 3.2.1 can be observed. Corresponds to Figures 2A (middle panel) and 3C.

flows0046.mp4 GLYPH&lt;16&gt; = 0 0046. : System is the higher activity region of regime II. The video shows two examples of the switching behaviour described in section 3.2.2 and Figure 4B.

flows0052.mp4 GLYPH&lt;16&gt; = 0 0052. : System is in regime III. This example shows the typical birth of small active clusters, corresponding to the situations shown in Figures 3D and 4C.

flows0060.mp4 GLYPH&lt;16&gt; = 0 0060. : System is deep in regime III. The dynamics is more turbulent than in the video for GLYPH&lt;16&gt; = 0 0052; there are more defects and clusters. Corre-: sponds to the example shown in Figure 2A on the right panel.

## D S upplementary figures

Supplementary Figure 1: Comparison of di GLYPH&lt;11&gt; erent modifications of the setup. (A) ( h max GLYPH&lt;0&gt; h min ) = d as a function of A . The crossover from regime I to regime II takes place at the same values as in the unmodified setup for a system without reservoir but with growth (blue circles), but is absent for a system lacking growth (red crosses). Vertical lines mark the crossover points in the unmodified setup. (B) Total amount of active material in the system at the end of the simulation relative to the starting value. In the non-growing system, the amount of active material is trivially constant over time. The system without reservoir shows the same qualitative dynamics as the unmodified setup. Data in this figure has been taken from one sample per value for A only. (C) The number of clusters Nc as a function of A . The crossover from regime II to regime III takes place at higher values as in the unmodified setup for a system with reservoir but without growth, which also explains the longer linear range in subfigure B. Clusters are also present for a system lacking growth, but it takes longer for them to appear. (D) The number of + 1 2 defects = N + as a function of A . The appearance of topological defects, for both setups, coincides with the appearance of additional clusters (compare subfigure C).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00768v1-figure-ce3b5684-7cab-4646-b6cd-97d08ac5cceb.jpg)

Supplementary Figure 2: Universality for di GLYPH&lt;11&gt; erent channel widths. ( h max GLYPH&lt;0&gt; h min ) = d (panel A) and Nc (panel B) against A . Both crossovers happen at the same values for both channel widths.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00768v1-figure-f9061d57-c9af-4721-99ca-3f04361c301c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00768v1-figure-8e2c2c53-590e-4017-9901-05c4a4698230.jpg)

Supplementary Figure 3: Transition between regimes II and III. Comparison of number of clusters Nc (A) and number of + 1 2 = topological defects N + in the system (B). Appearance of clusters and defects coincide around the second crossover.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00768v1-figure-ff9c32ac-2dfa-4785-88e4-6d7b770fee04.jpg)

Supplementary Figure 4: Graphical abstract. Biological materials such as bacterial biofilms and eukaryotic cells combine their intrinsic activity with growth dynamics to create distinct patterns of motion for invading confined spaces.
Felix Kempf Romain Mueller Erwin Frey Julia M. Yeomans Amin Doostmohammadi
Active Matter Invasion
Biological active materials such as bacterial biofilms and eukaryotic cells
thrive in confined micro-spaces. Here, we show through numerical simulations
that confinement can serve as a mechanical guidance to achieve distinct modes
of collective invasion when combined with growth dynamics and the intrinsic
activity of biological materials. We assess the dynamics of the growing
interface and classify these collective modes of invasion based on the activity
of the constituent particles of the growing matter. While at small and moderate
activities the active material grows as a coherent unit, we find that blobs of
active material collectively detach from the cohort above a well-defined
activity threshold. We further characterise the mechanical mechanisms
underlying the crossovers between different modes of invasion and quantify
their impact on the overall invasion speed.