## Spatio-temporal patterns in Growing Bacterial Suspensions: Impact of Growth dynamics

1 ∗ 1 †

Pratikshya Jena and Shradha Mishra Department of Physics , Indian Institute of Technology (BHU) Varanasi , 221005

(Dated: August 2, 2024)

The field of active matter explores the behaviors of self-propelled agents out of equilibrium, with active suspensions, such as swimming bacteria in solutions, serving as impactful models. These systems exhibit spatio-temporal patterns akin to active turbulence, driven by internal energy injection. While bacterial turbulence in dense suspensions is well-studied, the dynamics in growing bacterial suspensions are less understood. This work presents a phenomenological coarse-grained model for growing bacterial suspensions, incorporating hydrodynamic equations for bacterial density, orientation, and fluid velocity, with birth and death terms for colony growth. Starting with low density and random orientations, the model shows the development of local ordering as bacterial density increases. As density continues to rise, the model captures four distinct phases-dilute, clustered, turbulent, and trapped-based on structural patterns and dynamics, with the turbulent phase characterized by spatio-temporal vortex structures, aligning with observations in dense bacterial suspensions.

## I. INTRODUCTION

The emerging field of active matter deals with the analysis of the individual and collective properties of self-propelled agents in out of equilibrium [1-6]. Among these, the investigation of active suspensions has been an intriguing research subject for the past few years due to their unique behaviour. And the collection of swimming bacteria in a solution is a simplest realization of active suspension and has a substantial impact on the context of active matter. The bacteria bath [7-11] are interesting model to study due to the simplicity in manipulating the experimental control parameters like concentration of bacteria, activity, and swimming behaviour.

Some studies [34, 35] have focused on the role of activity and there are other studies performed and characterised the turbulence in dense bacterial suspensions [20, 33, 36-38]. However, bacterial turbulence in the growing bacteria suspension are still poorly explored although the growth phenomena is always present in natural bacterial suspensions. In recent experimental work [39], the authors studied the dynamics of growing bacteria at the oil-water interface and discovered the existence of four dynamical phases. But a theoretical model to explore the appearance of distinct phases in growing bacterial suspension is still lacking. This motivated us to present our current work.

Several research endeavors have been conducted on this model and reported a variety of novel nonequilibrium phenomena including emergence of collective motion [12-14], pattern formation [15, 16], reduction in viscosity [17], enhancement of diffusivity [18] and energy transfer [19]. Apart from these interesting properties, the occurrence of spatio-temporal patterns resembling the features of active turbulence in such systems of micro swimmers is highly unexpected and has been observed by researchers both experimentally and theoretically [20-25]. In conventional turbulence, the vortices of different sizes are developed when inertial force dominates over the viscous force and need a supply of external energy to maintain it [26-30]. Alternatively in dense bacterial suspensions, the energy injection is internal and turbulence can occur in over-damp systems. This motivated the researchers in recent years to explore the bacterial turbulence in greater details [31-33].

∗ pratikshyajena.rs.phy20@itbhu.ac.in

† smishra.phy@iitbhu.ac.in

We develop a phenomenological coarse-grained model for a growing bacterial suspension by writing the hydrodynamic equations of motion for the bacterial local density, orientation and fluid velocity. An additional birth and death term is introduced to take care of growth of colony.

We started with low density and random orientations of bacteria and with time local ordering among bacteria develops as bacteria grows with time with increasing density. As density increases further, the locally ordered state transit to the dynamical vortex state. Further at high densities system becomes heterogeneous with slower dynamics. Thus, based on the structural pattern and dynamics of the system, four distinct phases are identified; namely dilute, clustered, turbulent and trapped phase. Notably, the turbulent phase exhibits spatio-temporal vortex structures, which are typically observed in dense bacterial suspensions [20, 34, 36].

The structure of the paper is organized as follows. In section II, we provide a detailed description of our model and the numerical scheme used for the simulation. Next, section III presents a comprehensive results of the model. Lastly, we discuss the relevance of our work on the basis of results obtained, summarizing in section IV.

## II. MODEL

We consider a suspension of growing bacterial colony in two-dimensions. The periodic boundary condition is implemented in both the directions. The system is modeled using coarse-grained equations of motion for the local density ρ ( r , t ), orientation P ( r , t ) and fluid velocity v ( r , t ) at position r and time t . The hydrodynamic equation for the density ρ ( r , t ) is:

$$\frac { \partial \rho } { \partial t } + \nabla. ( \rho \upsilon ) + v _ { 0 } \nabla. ( \rho \mathbf P ) = D _ { \rho } \nabla ^ { 2 } \rho + g ( \rho ) \quad ( 1 ) \quad \sigma _ { \alpha \beta } ^ { t o t a }$$

the corresponding equation for the local orientation P ( r , t) is:

$$& \frac { \partial P } { \partial t } + \lambda _ { 1 } ( P. \nabla ) P + \lambda _ { 2 } \nabla | P | ^ { 2 } + ( v. \nabla ) P = D _ { r } \left ( \frac { \rho } { \rho _ { c } } - 1 \right ) P \begin{array} { c c } \quad \ v y \ u n \\ \quad \text{with } t \\ \text{super} \\ - \beta | P | ^ { 2 } P - \omega _ { \alpha \beta } P _ { \beta } - v _ { 1 } v _ { \alpha \beta } P _ { \beta } - \sigma _ { 1 } \nabla \rho + D _ { T } \nabla ^ { 2 } P \end{array} \begin{array} { c } \quad \text{$\upsilon$ with } t \\ \quad \text{$\text{$is the} \\ \text{position} \\ \text{$position} \\ \text{$-tensil} \\ \text{$tutorial} \end{array} \end{array}$$

and the equation for the fluid velocity v ( r , t) is

$$\frac { \partial v } { \partial t } + v. \nabla v = \eta \nabla ^ { 2 } v - \nabla p + \nabla. \sigma _ { \alpha \beta } ^ { t o t a l } \quad \ ( 3 ) \quad \stackrel { \ell \times } { \text{ field} }$$

The density equation 1 is a continuity equation with additional term g ρ ( ) = ρ (1 -ρ ρ 0 ) which incorporates birth and death of the bacteria. This type of term was first used by Malthus, [40] and consequently, the flocks with birth and death are also known as Malthusian flocks and later the same term is used in other models of polar flock with birth and death [4, 6]. We fix the ρ 0 = 10, such that the birth rate is much higher than the death rate. The second term ∇ . ( ρ v ) in equation 1 refers the convection with response to the fluid flow and the third term v 0 ∇ . ( ρ P ) is the convection due to the bacteria's own orientation. The v 0 is the self-propulsion speed of bacteria. The first term on the right hand side is the diffusion term, which tries to homogenise the density of bacteria in the suspension.

The equation 2 presents the dynamics of orientation of bacteria and is quite similar to the Toner-Tu equation [41] for flocking with few additional terms due to the coupling of orientation field with background fluid velocity. The second and third terms on left hand side λ 1 ( P . ∇ ) P and λ 2 ∇| P | 2 are the convective nonlinearities present due to absence of Galilean invariance in these systems [41, 42]. The last term on the left hand side is the convective non-linearity due to fluid velocity. The first two terms and the last term D T ∇ 2 P on the right hand side are the deterministic part of Time dependent Ginzburg Landau (TDGL) model [43] for orientable rod type bacteria's. In addition two terms on the right hand side ω αβ = ( 1 2 ∂ α v β -∂ β v α ) and v αβ = 1 2 ( ∂ α v β + ∂ β v α ) represent, how the flow affects the local orientation of bacteria and are called the vorticity and strain-rate tensors respectively. The term, ∇ ρ is the inertial coupling in local orientation of bacteria in the suspension.

The equation for the fluid velocity equation 3 contains the terms present in Navier-Stokes equation with the additional force term due to the active and passive stresses in the suspension. The origin of the term ∇ .σ total αβ is due to the stresses experienced by the fluid including the internal stress due to viscosity and also the external stress by the suspended bacteria. In this case the total stress σ total αβ = σ αβ + σ a αβ , σ αβ and σ a αβ are the passive and active contributions respectively. The constitutive equation for stress tensor is formulated as:

$$( 1 ) \quad & \sigma _ { \alpha \beta } ^ { t o t a l } = 2 \eta _ { 1 } v _ { \alpha \beta } + \frac { 1 } { 2 } \nu ( P _ { \alpha } h _ { \beta } + P _ { \beta } h _ { \alpha } - \frac { d } { 2 } ( P _ { \gamma } h _ { \gamma } \delta _ { \alpha \beta } ) ) + \zeta q _ { \alpha \beta } \\ \intertext { P ( r, } \quad & \tau _ { \mu, - \ e \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots - \ e \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots +\dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots + \dots +\dots - \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \ dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e \dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ edots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \e Dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + - \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots - \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dts + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dts - \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e ds + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e ds - \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e d: + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e dots + \ e d$$

The first term corresponds to the internal stress caused by the viscosity, the second term is the passive stress with the flow coupling coefficient ν attributable to any suspended passive particles. Additionally, the term ζq αβ is the active stresses where q αβ = P P α β -1 d δ αβ . The positive and negative values of ζ correspond to the ex -tensile and contractile stresses respectively in most bacterial suspension [1]. For ζ = 0, the hydrodynamic coupling is purely passive. The equation 3 is solved using the stream function approach [44] and taking the 'curl' ' × ' on both sides to get the equation for the vorticity field ω = ∇× v ,

$$\frac { \partial \omega } { \partial t } + ( v. \nabla ) \omega = \eta \nabla ^ { 2 } \omega + \nabla \times ( \nabla. \sigma _ { \alpha \beta } ^ { t } ) \quad ( 5 )$$

Further we get the Poisson's equation for a scalar field ψ using

$$v = ( \partial _ { y } \psi - \partial _ { x } \psi )$$

$$\nabla ^ { 2 } \psi = - \omega$$

The intrinsic length and time scales in the system is chosen by persistent length l 0 = v 0 D r = 0 4 and . τ = 1 /D r = 2 respectively. Later all the lengths and times can be defined in terms of l 0 and τ , such that after writing the rescaled equations, all terms are dimensionless. To solve the system we perform the numerical integration of equations 1,2, 5, 6 and 7 using Euler's scheme with grid size ∆ x = 0 5 = 1 25 . . l 0 and ∆ t = 0 0025 . τ . The above choice of grid size lead to roughly each vertices of the square grid consist of approximately 75 -100 E. Coli bacteria (keeping typical size of bacteria roughly 10 µm ) and the reorientation rate τ is approximately after 100 steps of persistent motion (with typical speed of bacteria of the order of 5 µm/s [45]).

The system is initialized with random, isotropic and homogeneous state with initial number density ρ = ⟨ ρ ⟩ ± 0 05, . where initial mean density ⟨ ρ ⟩ = 0 1. . The initial state for orientation is random with two components of P = ( P , P x y ) and scalar field ψ randomly chosen between ± 0 1 . and ± 0 05 respectively. . The coefficients

FIG. 1: The plots (a-d) showcase the real snapshots of bacterial density at different mean densities of bacteria in the system. The color on the heatmap presents the values of local density at each lattice point. The plots (e-h) illustrate the zoomed portion (square box) of the above snapshots in sequence. The black arrows in (e-h) present the local bacterial orientation at that lattice point. The length of the arrows are proportional to the magnitude of local P . For better visualization, we have drawn the arrows by coarse-graining the lattice points in the cell of size 20 × 20.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-bac083ad-844f-4ee1-9934-466da3ebdd65.jpg)

FIG. 2: The plot (a-b) showcase the fluctuation in rms value of global order parameter P rmsf ( ⟨ ρ ⟩ ) and fluid velocity v rmsf ( ⟨ ρ ⟩ ) vs. mean density ⟨ ρ ⟩ in sequence.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-2a6b7b02-8eb7-4ccc-ad8a-b37a0cfda94a.jpg)

( D ,λ ,λ , β ρ 1 2 and D T ) of different terms are set to 1. The strength of active coupling ζ = -2 v 0 and v 0 is fixed to 0 2. . The results are also tested for the extensile case ζ &gt; 0 and found qualitatively same as the contractile case. The critical density ρ c beyond which mean field theory suggest ordering among the particles is fixed to 0 5. . The coefficient of inertial term σ 1 is set to 1 for most of the detailed results. To see the effect of inertia later it is varied from 1 to 100. Larger value of σ favours the stronger alignment of orientation P in the direction of gradient of density. The stability of the system is checked for the set of chosen parameter before start of the simulation.

The simulation is carried out for difrerent system sizes i.e, K = 256 512 600 , , for total simulation time t =

3000 = 1500 τ . All the real snapshots in the Results section are generated for K = 256. The one simulation step is counted when equations 1, 2, 5, 6 are updated for all lattice points. Statistical averaging is performed over 100 independent realisations for better statistics.

## III. RESULTS

Starting with the random homogeneous initial state and low mean density we first analysed the snapshots of the system at different times. Since, the density of bacteria is growing with time, all the measurements are done with respect to mean density. In Fig 1(a-d), we showcase the snapshots of the local

FIG. 3: The plots (a-b) depict the vorticity field Ω( r ) in dilute, clustered, turbulent and trapped phases respectively in the system for ⟨ ρ ⟩ = 0 12 2 94 0 75 10 00 in sequence. These snapshots are taken at a particular value . , . , . , . of global mean density as marked by a blue tick in the Fig. 2. The color of the heatmap shows the magnitude of vorticity in all the phases. The plots (e-h) showcase the zoomed picture of the square box in the above plots (a-b) in sequence. Details of the arrows are same as in Fig. 1.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-9f9eab00-db64-416e-859c-60a24a04c34c.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-2a061e04-05d5-45d9-8d74-da1f231e66b5.jpg)

density with the local orientation at mean densities ⟨ ρ ⟩ = 0 12 2 94 0 75 10 00 respectively. The presence of . , . , . , . g ρ ( ) term makes the mean density of bacteria to grow with time. In Fig. 1(a), when the density of bacteria is low ⟨ ρ ⟩ = 0 1 . the bacteria are sparsely dispersed with almost homogeneous density in the system with negligible local clustering and local orientation. Then with time as density grows, the particles starts to form small clusters and small local orientation among the bacteria develops as shown in Fig. 1(b) for ⟨ ρ ⟩ = 2 94. . As the density grows further vortex type pattern starts to develop in the orientation field of particles, as shown in Fig. 1(c) at ⟨ ρ ⟩ = 7 5 and as density increases further . system becomes heterogeneous with random orientation as shown in Fig. 1(d) for ⟨ ρ ⟩ = 10 00. . The varying spatial distribution of local density at different times are depicted by the heatmaps and local orientation of the bacteria is shown by small black arrows. Here, the length of the arrows are proportional to the magnitude of local P . For clarity arrows are drawn with averaging the local P (r) in a cell size of 20 × 20. In the bottom panel of the Fig. 1(e-h) we show the zoomed plots of local density shown in Fig. 1(a-d) respectively.

varies non-monotonically with mean density ⟨ ρ ⟩ .

The above snapshots are observed at four distinct times. Next, we focus on the properties of Bacterial suspension with continuous increase of density. The Fig. 2(a-b) illustrates the root mean square fluctuations of the order parameter P rmsf ( ⟨ ρ ⟩ ) and velocity v rmsf ( ⟨ ρ ⟩ ) vs. ⟨ ρ ⟩ respectively. The root mean square fluctuations in orientation P rmsf and velocity v rmsf is obtained by calculating the fluctuations in magnitude of global P and v , where global P and v is obtained by averaging the local orientation and velocity over all the points in space. The plots suggest that both the P rmsf ( ⟨ ρ ⟩ ) and v rmsf ( ⟨ ρ ⟩ )

Observing such non-monotonic density dependence of root mean square fluctuations and spatio-temporal snapshots of local density shown in Fig. 1(a-h), we identified the four distinct phases in the system refereed as (i) dilute, (ii) clustered, (iii) turbulent and (iv) trapped in chronological order. The vertical lines in Fig. 2(a-b) are drawn to mark the distinct phases in the system. The four snapshots in Fig. 1(a-d) are at four different densities in four regions marked by the blue short vertical lines in Fig. 2(a-b). In the dilute phase, the local density of bacteria is low (Fig. 1(a,e)), and there is almost no ordering (random orientation) and zero P rmsf ( ⟨ ρ ⟩ ) and v rmsf ( ⟨ ρ ⟩ ). In contrast, in the clustered phase, low and high-density contrasts begin to develop due to the formation of clusters (Fig. 1(b-f)), the local fluctuation in P and v starts to grow. On further increment of density, in the turbulent phase both P rmsf ( ⟨ ρ ⟩ ) and v rmsf ( ⟨ ρ ⟩ ) attain a maxima and local density becomes more inhomogeneous as shown in Fig. 1(c) and (g). Later, we will also observe that in this phase swirling structures develop in the system, and finally, the P rmsf ( ⟨ ρ ⟩ ) and v rmsf ( ⟨ ρ ⟩ ) starts to decreases as system becomes heterogeneous in the trapped phase. We named this phase as trapped due to the weaker dynamics present in the system, which we will discuss later in the paper.

The four phases we mentioned above can be compared with the recent experiments on growing bacteria on the oil-water interface [39]. The plots in Fig 2(a-b) are obtained for K = 256. We also investigated the same variables for system size of K = 512 and the results are consistent. Till now, we have identified the four phases based on the static properties of the suspension using the local density, orientational ordering and velocity

FIG. 4: (Color online) The semi-log-y plot illustrates P (Ω) of vorticity field Ω in dilute, clustered, turbulent and trapped phases. The legends represent different phases with different line styles. In inset, we plot the log-log plot of the same showing the power-law decay of distribution in the turbulent phase for a range of Ω.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-1d49a79b-c810-48ce-8686-3119be5bdf57.jpg)

fluctuations. Next, we characterise these phases based on other structural and dynamical properties of the suspension.

We calculated the vorticity field ω ( r ) at the specific densities discussed before in Fig. 3(a-d). The corresponding zoom plots are shown in bottom panel Fig. 3(e-h). The colors represent the normalised vorticity Ω( r ) = ω ( r ) /ω ¯, where ¯ is the average vorticity in the system. ω The local orientation of the bacteria P ( r ) is shown by arrows. The arrows are drawn by local coarse-graining of the P ( r ) as in Fig. 1. In Fig. 3(a-b) and (d) and (e-f, h), the vorticity field's magnitude is small, and the bacterial orientations appear quite random, as shown by the black arrows. This low magnitude of vorticity indicates a low level of swirling motion. Fig 1(c, g) clearly depicts formation of vortices in the system (the magnitude of Ω( ) increased substantially). r Local inhomogeneity in Ω is developed with random orientation of bacteria and defect type structures at the core of high and low vorticity field as can be seen by dark spots in Fig. 3(g). The animation of the vorticity field of the system starting from low to high mean density ( ⟨ ρ ⟩ ) is given in SM.

The presence of large scale vortex in bacterial suspension is one of the key characteristics of bacterial turbulence. Therefore, we characterise the statistics of the four phases by measuring the probability distribution function of the magnitude of local vorticity P (Ω). In Fig. 4, we depict the P (Ω) vs. Ω for dilute, clustered, turbulent and trapped phase on semi-logy scale and the inset shows the same plot in log -log scale. The nature of distribution is very different in turbulent phase in comparison to other three phases. The P (Ω) decay as a power law with 1 / Ω for a range of Ω, with an exponential tail at larger Ω. The appearance of power law distribution suggests the finite probability of having the vortices of large magnitude in turbulent phase, but such vortices are less probable in other three phases.

FIG. 5: (Color online) The plots (a-c) showcase the spatial correlation functions for the various i.e, bacterial density C ρ ( r ), bacterial orientation C p ( r ), fluid velocity C v ( r ) respectively in four different regimes i.e, dilute, clustered, turbulent and trapped. The plot (d) illustrates the characteristic length ( L v ( ⟨ ρ ⟩ )) vs. mean density ( ⟨ ρ ⟩ ). The legends represent the system sizes taken into consideration.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-79178148-efe3-4e2f-a6ec-26c0a5d1343c.jpg)

Thus far, our observation is either on the local spatiotemporal pattern or the global order parameters of the system. We further, quantify the characteristics of the four phases by calculating the spatial and temporal correlations of bacterial density, orientation and fluid velocity. The two-point spatial correlation functions are defined by C f ( r )= ⟨ δf ( r 0 ) · δf ( r + r 0 ) ⟩ r 0 where, f = ρ , P and v . The δf is the fluctuations from their respective mean values. The angular bracket represents the spherical averaging and average over 100 independent realizations. In Fig 5(a-c), we present the spatial correlation functions of these three fields in sequence. The four plots in each window are for four different phases as marked by the different symbols. We observe that the orientation and velocity correlation increases from dilute to turbulent phase and subsequently decreases in the trapped phase as shown in Fig. 5(b-c). But density correlation is maximum in the clustered phase and lowest in trapped phase Fig. 5(a). This implies, the density patterns formed in respective fields are more correlated in the dilute to clustered and the onset of turbulent phase but, the correlation decreases in the trapped phase. Although in the turbulent phase, the flow is chaotic and has highly unpredictable motion at various scales, at the onset of turbulence, correlated structures are formed. That leads to increase in spatial correlations initially but decay at latter times due to the destruction in ordered patterns. In the trapped phase, the bacteria motion is hindered as the system has become densely packed. This suppresses the ordering

and as a result, the correlation decreases. Unlike the other three phases, in the turbulent phase the spatial correlation for density C ρ ( r ) shows a large dip in the negative side, which represents the presence of high density of bacteria at the core of the vortices and low density around it.

Further, we calculate the characteristic length for the fluid velocity L v ( ⟨ ρ ⟩ ) of the system to examine the average size of the structures or patterns generated in the system as the growth of bacteria progresses. It defined by the distance by which the the two-point correlation function C v ( r ) decay to 0.5 of its first value. Fig 5 (d) depict the correlation length L v ( ⟨ ρ ⟩ ) vs. ⟨ ρ ⟩ . From the Figure, it is evident that as soon as the system enters the turbulent region, the length L v ( ⟨ ρ ⟩ ) has an increment initially because of the formation of vortices like ordered patterns. At later times, the turbulent phase becomes more chaotic, disrupting the formation of ordered domains. This leads to decay of correlations as well as the characteristic length. This kind of sharp decay suggests random and chaotic nature of turbulent phase in comparison to other phase. As mentioned earlier, in trapped phase the local trapping of the bacteria leads to decrement of the length L v ( ⟨ ρ ⟩ ) due to very high density of particles. To investigate the effect of system sizes, we plot the correlation lengths for all the variables for system K = 64 128 256 as depicted , , in the Figure. We observed the results are consistent with different system sizes.

We further characterise the local correlations among the three fields in the system by calculating the cross-correlations. In Fig. A1(a-c) we show the cross correlation of density-orientation C ρP , density-fluid velocity C ρv and bacterial orientation and fluid velocity C Pv respectively. The details of the three correlations are given in Appendix A. The three fields are uncorrelated in dilute and clustered phase, with zero cross-correlations, whereas as soon as the system transit to the turbulent phase local correlations start to pick non-zero values. The local ordering in the turbulent phase appears at the cost of density inhomogeneity, with high ordered region having lower density and vice versa as can be seen by big negative dip in C ρP in Fig. A1(a), which is also see by bigger negative dip in the spatial correlation C ρ ( r ) in Fig. 5(a). The C ρv shows a positive peak in the turbulent phase suggesting the strong local correlation among the bacteria density and fluid velocity. Similar trend also found in the correlation of bacterial density with vorticity. Additionally the local orientation and fluid velocity also show similar trends as density and fluid velocity with smaller peak in the turbulent phase. As we move towards the trapped phase, C ρP and C ρv approaches to zero, because density starts to become homogeneous whereas the local orientation and fluid velocity shows an anti-correlation, due to the local trapping of bacteria due to high global density and fluid velocity.

Further we investigate the dynamics of different phases.

FIG. 6: (Color online) The plot demonstrates the autocorrelation functions of fluid velocity ( C v ( )) in dilute, t turbulent and trapped phase. We compare the autocorrelation of simple fluid with the bacterial suspension (magenta dashed line). The inset shows the plot in log-log scale demonstrating the 1 /t decay of autocorrelation in case of simple fluid which is very different from the bacterial suspension.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-01ac5074-a051-48b4-a2e4-1c9ead966c55.jpg)

FIG. 7: The plot showcases the speed ( U ) of the passive probe vs. ⟨ ρ ⟩ .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-b8616e0a-e99a-45f0-abd6-db692fc3e964.jpg)

Therefore, we calculated the fluid velocity autocorrelation function. Fig 6 depict the auto-correlation function for fluid velocity C v ( ) t sequentially in dilute, turbulent and trapped phases. The details of the calculation of auto-correlation functions in given in Appendix B. It can be evident from the plot that the fluid velocity is weakly correlated in dilute phase, higher correlations are developed in turbulent phase and further the correlations suppresses in the trapped phase due to less dynamics present in the system. Unlike the dilute phase, the auto-correlation shows a strong negative values in the turbulent and trapped phase suggesting the rotation is present in the system, which is highest in the turbulent phase and suppresses in the trapped phase. Alternatively, we present a comparison of our results with the simple Navier-Stokes fluid and find an algebraic decay with 1 /t for the C v ( ) as shown t by the solid line in Fig. 6(c) [46]. That represents that the presence of active Bacteria makes the fluid

FIG. 8: (Color online) The plot presents the energy spectrum E q ( ) vs. q in different phases. The q dependence of E q ( ) is shown by different slopes in different regimes suggesting multiple scaling regimes are present in the system in different phases.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-b8183632-642e-464e-9a74-a0efb2bf4eb0.jpg)

FIG. 9: (Color online) The plot depicts the density fluctuation (¯ σ/N ¯ ) vs. ¯ N in dilute, clustered and turbulent phases. The large number fluctuations in all phases for small scale are shown by the lines of slopes 0 5. . We consider system size K= 600.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-dd46fffb-579f-45a8-a3b5-cd4c879f6df9.jpg)

turbulent, otherwise the fluid itself is in the inertial regime.

We further study the dynamics of different phases placing a probe particle in the system.

## The dynamics of different phases using a probe particle:-

The dynamics of the probe particle is governed by an over-damped equation with three contributions to its velocity coming from: local density gradient of bacteria, fluid velocity and local orientation of bacteria. The dynamics of probe is affected by the bacterial suspension, but the probe itself is passive and cannot influence the flow of surrounding fluid. Hence, the probe acts like a tracer to study the dynamics of the suspension [47, 48]. We started with low density of bacteria in the suspension and placed 1000 probe particles and characterised the speed of probe particle as the density of bacteria in the suspension increases. The details of the governing equations of motion for the probe particles and calculation of speed of particles is given in Appendix C. In the following Fig. 7 we show the mean speed U of probe particles vs. ⟨ ρ ⟩ . As the system goes from a dilute to a clustered phase and then to a turbulent phase, the speed increases, peaking during the turbulent phase, before decreasing in the trapped phase. This behavior mirrors the system's dynamics and is reflected in the plot of fluid velocity root mean square fluctuations v rmsf , as illustrated in Fig. 2. Similar kind of observation is also found in a recent experimental work [49, 50], where the authors have studied the dynamics of tracer in bacterial suspensions at different densities.

To further quantify the system's energetic, we examine the energy spectrum. So far, it is evident that the turbulent phase is highly dynamic, exhibiting persistent spatio-temporal patterns. The energy spectrum of the fluid velocity is a standard approach to characterise the turbulence in normal fluid. We calculate the energy spectrum of fluid defined by E q ( ) = ∫ &lt; v ( r ) · v ( r + R ) &gt; exp( -i q · R ) d R . This is employed in order to quantify the spread of energy cascade through the system at different length scales. Fig. 8 showcases the same for the different phases in sequence as described above for σ 1 = 1. The comparison of energy spectrum for other values of σ 1 = 10 50 100 is , , shown in Appendix E. For the dilute phase, the positive slope at lower q with E q ( ) ∼ q 2 indicates increasing energy towards intermediate q and the steep decay shows the quick dissipation of energy at small scales (large q ). The form of the spectrum remains almost similar for the clustered phase as well. The E q ( ) in the turbulent phase, increases at lower q showing a positive slope of value 5 ( E q ( ) ∼ q 5 ), which is higher than that of dilute and clustered phase. That indicates the rate of increase of energy is more in turbulent phase. The peak of energy density shifts towards lower q value showing the transfer of energy towards large scale from small scales. This is an indicative of inverse energy cascade in 2 -dimensional active turbulence [51-54]. Here, the source of energy is intrinsic and comes from the active nature of the bacteria in the system. In this phase, we find E q ( ) shows an intermediate range of q values where E q ( ) ∼ q -2 . This intermediate q range shows the range of different sizes of vortices, through which energy is dissipated and further energy dissipation happens sharply as we go for larger q values. Despite the inverse cascade, the energy dissipates at small scales also due to viscous dissipation or other mechanism of dissipation showing a negative slope of E q ( ) ∼ q -7 2 . . The length is different energy scales, where the inverse cascade plays a role can be tuned by tuning the inertia in the system. The bigger range of length scale is found when inertial term is small (small σ 1 ) and it is suppressed as we increase σ 1 . For higher value of σ 1 when system is in turbulent phase as shown in the Fig. E2, the peak of the energy spectrum shifts towards even lower value of q , representing the accumulation of energy at

large scales from small scales is favoured by the inertia. Similarly from the plot of the energy spectrum in the trapped phase, we get two different scaling regimes having different slopes. The peaks are at higher q values (smaller scales) in comparison to turbulent phase. That indicates the trapping of bacteria in this phase, the probability of energy transfer from small scales to large scale is less. Different from the turbulent phase, here for a range of intermediate q values energy spectrum E q ( ) increases with q as q 1 5 . , representing the energy input on relatively larger q values due to high density of bacteria. Further, the peaks of energy spectrum E q ( ) shifts towards lower q values with increasing value of σ 1 favouring the accumulation of energy at the larger scales.

Density fluctuation is one of the interesting quantity in active system. To observe the effect in growing bacteria suspension, we measure the same in the system. We calculate the density fluctuations of bacteria in cells of different sizes. The details of the density fluctuation calculation in given in Appendix D. In Fig 9 we show the plot of density fluctuation ¯ σ √ ¯ N vs. ¯ , N where ¯ N is the mean number of particles in different cells of the same size. For dilute and clustered suspension the density fluctuation is larger for small ¯ , N ¯ σ ∝ ¯ N and attains the equilibrium limit for larger ¯ , N ¯ σ ∝ √ ¯ . N Hence, at larger distances the density is homogeneous in dilute and clustered phase. As system approaches to the turbulent phase, the density fluctuations are much larger at small scales suggesting large density inhomogeneity in the system and further the density fluctuation suppresses at larger ¯ N (¯ σ/ √ ¯ N decays with ¯ ) showing the N local trapping of bacteria at the core of the vortices can be seen in the snapshot shown in Fig. 1(c) and (g). In Fig. E2 (b), we show ¯ σ/ √ ¯ N vs. ¯ for large inertia terms N ( σ 1 = 10 50 100) when system is in turbulent phase. As , , σ 1 increases, the dip in ¯ σ/ √ ¯ N shifts towards larger ¯ , N showing the large scale structures as can be observed in the snapshots in Fig. E1 (b-d). Additionally, similar information can be analysed from the energy spectrum E q ( ) in Fig. E2 (a).

## IV. DISCUSSION

We developed a coarse-grained model to study the statistical and dynamical properties of the growing bacterial suspension by allowing the growth of bacteria in the suspension with time introduced by a birth and death term in the hydrodynamic equations of motion. This model provides a realistic description of growing bacterial dynamics in suspension. In this work, we present a comprehensive results of the growing bacterial suspensions.

The properties of the suspension is characterised by local density of bacteria, local orientation and fluid velocity. Staring from low density of bacteria, with time the density of bacteria grows and system transits from dilute to clustered and further to turbulent phase, and at very high density we obtain trapped phase. It is observed that local orientation and fluid velocity develop correlation in the turbulent phase with vortex type patterns. The turbulent phase has specific spatio-temporal patterns resembling the bacterial turbulence in the previous studies [20, 34, 36]. Additionally, the turbulent phase lead to the enhanced dynamics with large speed of probe particle. The energy spectrum shows the presence of dissipation at different length scales due to presence of vortices in the turbulent phase and absent in other phases. Further, we found the turbulent nature of the suspension is suppressed, if we increase the effect of inertia σ 1 and instead wave like patterns are developed [55, 56].

Our study identifies four phases-dilute, clustered, turbulent, and trapped-in the growing bacterial suspension and these phases overlaps with the phase found in recent study of bacterial suspension on oil-water interface [39]. Unlike the previous study (theoretical model on dense bacterial suspension), the density of bacteria in our study continuously grow in time, which makes it more closer to many natural bacterial suspension.

The present model and our results can be useful to understand the dynamical and statistical properties of growing bacterial suspension in natural systems. In the current study, the suspension evolves in the clean environment, it is interesting to explore the dynamics of the suspension in the presence of random impurities.

## ACKNOWLEDGMENTS

P.J. and S. M. thank Dr. Sivasurender Chandran for useful discussions. P.J. gratefully acknowledge the DST INSPIRE fellowship for funding this project. The support and the resources provided by PARAM Shivay Facility under the National Supercomputing Mission, Government of India at the Indian Institute of Technology, Varanasi are gratefully acknowledged by all authors. S.M. thanks DST-SERB India, ECR/2017/000659, CRG/2021/006945 and MTR/2021/000438 for financial support. P.J. and S.M. also thank the Centre for Computing and Information Services at IIT (BHU), Varanasi.

## Appendix A: Cross-correlation

The cross correlation function is defined as C AB ( ⟨ ρ ⟩ ) = ⟨ δA ( r ) δB ( r ) ⟩ r . To obtain the crosscorrelation function, we calculate the product δA ( r ) δB ( r ) on each lattice point, where δA and δB are the fluctuations in the desired quantities. Furthermore, we take the average of this quantity over all the lattice points. Additionally, for better statistical average we take averaged data of 50 ensembles and plot C AB ( ⟨ ρ ⟩ ) against the global mean density ( ⟨ ρ ⟩ ).

FIG. A1: The plots (a-c) represent the cross-correlations C ρp ( ⟨ ρ ⟩ ), C ρv ( ⟨ ρ ⟩ ) and C pv ( ⟨ ρ ⟩ ) vs. ⟨ ρ ⟩ respectively in the system. This plot shows the correlation among different fields in different phases.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-22d02594-1355-42ca-90a6-2ea4d456a727.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-5c2d7f4f-746a-440a-ae62-d2f3564aa20d.jpg)

FIG. E1: The plots (a-d) showcase the real snapshots of density field in the turbulent phase for σ 1 = 1 10 50 100. , , , The colors on the heatmap present the magnitude of local density ( ρ ( r )).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-6fb8f497-d2b2-4b9c-88ba-bead341573ec.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-01d84241-ef54-458d-94fb-52846c135ca7.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-3d85bce3-5da1-4cbe-9498-41b63659f29b.jpg)

## Appendix B: Auto-correlation

## Appendix D: Density fluctuation

We define the auto-correlation C v ( ) = t &lt; δ v ( t 0 ) δ v ( t + t 0 ) &gt; . And δv = ( ) v t -¯ , v where ¯ is the mean value v of v ( ). t Here, t 0 is taken as different reference points within the time span of each phase. By observing the snapshots for vorticity field, we decided different time span to obtain the auto-correlation in different phases. Next, we calculate the auto-correlation function by averaging the quantity over all the lattice points and time in a particular phase.

## Appendix C: Dynamics of probe particle

To examine the dynamical properties of different phases, we use probe particles. We update the position r = ( x,y ) by using the over-damped Langevin equation. The position update of the tracer is

$$\frac { d { \mathbf r } } { d t } = - \nabla \rho + \rho P + v \quad \quad \ \ ( C 1 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \quad \frac { ( C 1 ) } { d t } \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Here, the forces on the right hand side are local density gradient of bacteria, fluid velocity and local orientation of bacteria. Further, to calculate the speed of the tag particle, we observed the trajectory of the particle by tracking it's co-ordinates. Using the co-ordinates in every 10 unit time differences, we calculate the mean speed of the particle in each phase during it's motion. Further, the speed is averaged over the time span of each phase.

We calculate the normalized density fluctuation ¯, σ by defining ¯ = σ &lt;ρ 2 n -ρ &gt; 2 ρ nn . To calculate the quantity numerically, we divide the system in to n multiple small cells and assume these cells as n ensembles. Here, ρ n is the mean density in each ensembles in the corresponding phase and ρ nn is the average of the mean density in the corresponding phase. ¯ N is obtained by ρ n ρ nn . Finally, we plot the density fluctuation ¯ σ √ ¯ N vs. ¯ N in Fig 9.

## Appendix E: Effect of Inertia

Until now, the results presented here are for fixed inertial coupling. Now, by increasing the strength of this term we increase the response of P with respect to δρ . To elucidate the significance of inertia in the turbulent phase, we examine the growth of bacteria by varying the coefficient σ 1 from 1 to 100. In Fig E1, we showcase the real snapshots of density field in the turbulent phase by fixing σ 1 = 1, 10, 50, 100 respectively. For all cases, we make sure that the system is in turbulent phase by characterizing P rmsf and V rmsf . For lower values of σ 1 = 1, 10, vortices are generated in the system distinctively, but with increasing value, the vortices are not sharp and wave like structures are developed instead and move outward from their central location. This type of concentrated wave formation with the higher value of σ 1 also seen in the work of [55, 56].

The Fig. E2 (a) represent the effect of σ 1 on the en-

FIG. E2: (Color online) The plots (a-b) represent the energy spectrum E q ( ) vs. q and density fluctuation ¯ σ/ √ ¯ N in turbulent phase respectively. The legends show different values of σ 1 = 1 50 100. The arrows represent the , , shifting of maximum of E q ( ) with increasing σ 1 .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00403v1-figure-aa11ac13-aa2b-4197-bf0f-6fa388baaa94.jpg)

ergy spectrum in the Fourier space in turbulent phase. It shows the peaks in the energy spectrum shift towards lower q values in case of higher values of σ 1 in turbulent phase. Similarly, Fig. E2 (b) ¯ σ / 1 √ ¯ N vs. ¯ N plot in turbulent phase shows the shifting of the bending towards large N .

- [1] M. C. Marchetti, J.-F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. A. Simha, Reviews of modern physics 85 , 1143 (2013).
- [2] S. Ramaswamy, Annu. Rev. Condens. Matter Phys. 1 , 323 (2010).
- [3] S. Ramaswamy, Journal of Statistical Mechanics: Theory and Experiment 2017 , 054002 (2017).
- [4] P. Jena and S. Mishra, Physics of Fluids 35 (2023).
- [5] V. Semwal, J. Joshi, S. Dikshit, and S. Mishra, Physica A: Statistical Mechanics and its Applications 634 , 129435 (2024).
- [6] P. K. Mishra and S. Mishra, Physics of Fluids 34 (2022).
- [7] H. C. Berg, E. coli in Motion (Springer, 2004).
- [8] G. Gr´ egoire and H. Chat´ e, Physical review letters 92 , 025702 (2004).
- [9] G. Soni, B. J. Ali, Y. Hatwalne, and G. Shivashankar, Biophysical journal 84 , 2634 (2003).
- [10] X.-L. Wu and A. Libchaber, Physical Review Letters 86 , 557 (2001).
- [11] D. Saintillan and M. J. Shelley, Physical Review Letters 100 , 178103 (2008).
- [12] A. Sokolov and I. S. Aranson, Physical review letters 109 , 248109 (2012).
- [13] S. D. Ryan, A. Sokolov, L. Berlyand, and I. S. Aranson, New journal of physics 15 , 105021 (2013).
- [14] J. Gachelin, A. Rousselet, A. Lindner, and E. Clement, New Journal of Physics 16 , 025003 (2014).
- [15] A. Baskaran and M. C. Marchetti, Proceedings of the National Academy of Sciences 106 , 15567 (2009).
- [16] M. E. Cates, D. Marenduzzo, I. Pagonabarraga, and J. Tailleur, Proceedings of the National Academy of Sciences 107 , 11715 (2010).
- [17] A. Sokolov and I. S. Aranson, Physical review letters 103 , 148101 (2009).
- [18] X.-L. Wu and A. Libchaber, Physical review letters 84 , 3017 (2000).
- [19] A. Sokolov, M. M. Apodaca, B. A. Grzybowski, and I. S. Aranson, Proceedings of the National Academy of Sciences 107 , 969 (2010).
- [20] H. H. Wensink, J. Dunkel, S. Heidenreich, K. Drescher,
- R. E. Goldstein, H. L¨wen, o and J. M. Yeomans, Proceedings of the national academy of sciences 109 , 14308 (2012).
- [21] S. Mukherjee, R. K. Singh, M. James, and S. S. Ray, Nature Physics 19 , 891 (2023).
- [22] B. Zhang, P. Leishangthem, Y. Ding, and X. Xu, Proceedings of the National Academy of Sciences 118 , e2100145118 (2021).
- [23] I. S. Aranson, Reports on Progress in Physics 85 , 076601 (2022).
- [24] R. Alert, J. Casademunt, and J.-F. Joanny, Annual Review of Condensed Matter Physics 13 , 143 (2022).
- [25] I. Maryshev, A. B. Goryachev, D. Marenduzzo, and A. Morozov, Soft Matter 15 , 6038 (2019).
- [26] K. R. Sreenivasan, Rev. Mod. Phys. 71 , S383 (1999).
- [27] J. Lumley and P. Blossey, Annual Review of Fluid Mechanics 30 , 311 (1998).
- [28] R. H. Kraichnan, The Physics of Fluids 9 , 1728 (1966).

[29]

J. Singh, M. Rudman, and H. Blackburn, Physical Re-

view Fluids

3

,

094607 (2018).

- [30] D. Buaria and K. R. Sreenivasan, Physical Review Fluids 5 , 092601 (2020).
- [31] K. V. Kiran, A. Gupta, A. K. Verma, and R. Pandit, Physical Review Fluids 8 , 023102 (2023).
- [32] E. Secchi, R. Rusconi, S. Buzzaccaro, M. M. Salek, S. Smriga, R. Piazza, and R. Stocker, Journal of the Royal Society Interface 13 , 20160175 (2016).
- [33] C. W. Wolgemuth, Biophysical journal 95 , 1564 (2008).
- [34] J. Dunkel, S. Heidenreich, K. Drescher, H. H. Wensink, M. B¨ ar, and R. E. Goldstein, Physical review letters 110 , 228102 (2013).
- [35] I. S. Aranson, A. Sokolov, J. O. Kessler, and R. E. Goldstein, Physical Review E 75 , 040901 (2007).
- [36] H. Wioland, F. G. Woodhouse, J. Dunkel, J. O. Kessler, and R. E. Goldstein, Physical review letters 110 , 268102 (2013).
- [37] E. Lushi and C. S. Peskin, Computers &amp; Structures 122 , 239 (2013).
- [38] M. James, W. J. T. Bos, and M. Wilczek, Phys. Rev. Fluids 3 , 061101 (2018).

- [39] B. Langeslay and G. Juarez, arXiv preprint arXiv:2306.15102 (2023).
- [40] T. Malthus, in British Politics and the Environment in the Long Nineteenth Century (Routledge, 2023) pp. 7784.
- [41] J. Toner and Y. Tu, Physical review E 58 , 4828 (1998).
- [42] S. Mishra, A. Baskaran, and M. C. Marchetti, Physical Review E 81 , 061916 (2010).
- [43] Q. Tang and S. Wang, Physica D: Nonlinear Phenomena 88 , 139 (1995).
- [44] M. H. Kobayashi and J. M. C. Pereira, International journal for numerical methods in engineering 62 , 1950 (2005).
- [45] C. Bechinger, R. Di Leonardo, H. L¨wen, C. Reichhardt, o G. Volpe, and G. Volpe, Reviews of modern physics 88 , 045006 (2016).
- [46] M. Ferrario, A. Fionino, and G. Ciccotti, Physica A: Statistical Mechanics and its Applications 240 , 268 (1997).
- [47] A. P. Petroff, X.-L. Wu, and A. Libchaber, Physical
- review letters 114 , 158102 (2015).
- [48] R. Chelakkot, M. F. Hagan, and A. Gopinath, Soft matter 17 , 1091 (2021).
- [49] D. Gautam, H. Meena, S. Matheshwaran, and S. Chandran, Phys. Rev. E 110 , L012601 (2024).
- [50] C. Xie, Y. Liu, H. Luo, and G. Jing, Micromachines 13 , 746 (2022).
- [51] A. Alexakis and L. Biferale, Physics Reports 767 , 1 (2018).
- [52] G. Boffetta and S. Musacchio, Physics of Fluids 13 , 1060 (2001).
- [53] Z. Xiao, M. Wan, S. Chen, and G. L. Eyink, Journal of Fluid Mechanics 619 , 1 (2009).
- [54] G. Boffetta, A. Celani, and M. Vergassola, Phys. Rev. E 61 , R29 (2000).
- [55] P. Jain, N. Rana, S. Ramaswamy, and P. Perlekar, arXiv preprint arXiv:2401.11927 (2024).
- [56] R. Chatterjee, N. Rana, R. A. Simha, P. Perlekar, and S. Ramaswamy, Physical Review X 11 , 031063 (2021).
Pratikshya Jena Shradha Mishra
Spatio-temporal patterns in Growing Bacterial Suspensions: Impact of Growth dynamics
The field of active matter explores the behaviors of self propelled agents
out of equilibrium, with active suspensions, such as swimming bacteria in
solutions, serving as impactful models. These systems exhibit spatio-temporal
patterns akin to active turbulence, driven by internal energy injection. While
bacterial turbulence in dense suspensions is well studied, the dynamics in
growing bacterial suspensions are less understood. This work presents a
phenomenological coarse-grained model for growing bacterial suspensions,
incorporating hydrodynamic equations for bacterial density, orientation, and
fluid velocity, with birth and death terms for colony growth. Starting with low
density and random orientations, the model shows the development of local
ordering as bacterial density increases. As density continues to rise, the
model captures four distinct phases; dilute, clustered, turbulent, and trapped
based on structural patterns and dynamics, with the turbulent phase
characterized by spatio-temporal vortex structures, aligning with observations
in dense bacterial suspensions.