A CRITERION OF COLLECTIVE BEHAVIOR OF BACTERIA

It was established in the previous works that hydrodynamic interactions between the swimmers can lead to collective motion. Its implicit evidences were confirmed by reduction in the effective viscosity. We propose a new quantitative criterion to detect such a collective behavior. Our criterion is based on a new computationally effective RVE (representative volume element) theory based on the basic statistic moments (e-sums or generalized Eisenstein-Rayleigh sums). The criterion can be applied to various two-phase dispersed media (biological systems, composites etc). The locations of bacteria are modeled by short segments having a small width randomly embedded in medium without overlapping. We compute the e-sums of the simulated disordered sets and of the observed experimental locations of Bacillus subtilis. The obtained results show a difference between these two sets that demonstrates the collective motion of bacteria.

1. Introduction.Experimental and theoretical models have been recently developed to examine fundamental aspects of collective motion exhibited by various biological systems.Following the seminal papers [18]- [17] and works cited therein we suggest that hydrodynamic interactions between the swimmers lead to collective motion when every bacterium interacts with other ones through the viscous environment.The implicit evidences of collective motion were confirmed by reduction in the effective viscosity [16].The theoretical investigations of collective motion were based on the considering the motion in the framework of mechanical dynamical systems [13]- [20].While the above presented theoretical and experimental results are related via the effective viscosity of suspension of swimmers, these relation lack simple and direct comparisons of viscosity for experimentally observed sets of bacteria and for simulated ensembles.
In the present paper, we propose a new quantitative criterion of collective behavior.The locations of bacteria are modeled by short segments having a small width randomly embedded in medium without overlapping.First, we theoretically simulate locations of particles called below by disordered sets of bacteria (DB sets for shortness) subjected to local viscous stresses and randomly reacted on hydrodynamic interactions.Second, we calculate the basic statistic moments of the constructed DB sets in terms of the generalized Eisenstein-Rayleigh sums (e-sums for shortness) introduced in [6].Further, we construct the e-sums (2) for the observed experimental locations (31 film frames) of Bacillus subtilis in a very thin liquid film [18].As the final step, we compare the e-sums of DB sets and of the experimental locations.The obtained results show a difference between these two sets of e-sums that demonstrates the collective motion of bacteria.We do not explain reasons for the collective motion and refer to [18]- [17], [13]- [20].
Our criterion of the collective behavior is based on a new RVE (representative volume element) theory proposed in [6], on the invariance of the effective transport properties on the conformal mappings [4] and on the algebraic dependence of the viscous lattice sums on the e-sums [6].The main object of this theory is a set of e-sums systematically investigated in [10], [8], [2], [3], [4], [15].According to the new RVE theory, the set of the e-sums (see (2) below and general formulas in [6]) determines macroscopically equivalent cells, i.e. cells having the same effective properties (conductivity, viscosity etc).The RVE is chosen as the minimal size cell from all the equivalent cells.The necessary justification of this theory can be found in the above cited works.It is worth noting that the usage of the esums implicitly takes into account high order correlation functions without their hard direct computations.In particular, we make simple implementation of the RVE theory [6] which demonstrates its numerical advantages in comparison with expensive computations based on the traditional statistical notion of the RVE.
It is worth noting that the criterion can be applied to various two-phase dispersed media (biological systems, composites etc).
2. Random location of segments.In the present section, we discuss theoretical simulations of elements embedded in a medium.Bacteria are modeled by short segments having non-overlapping thin δ-security coatings.Random locations of segments are generated in the following way.The plane geometry is considered as the complex plane C of the complex variable z = x + iy with standard designations accepted in complex analysis where i denotes the imaginary unit, Re and Im the real and imaginary parts, the bar stands for the complex conjugation.
The complex numbers 1 and i can be considered as the fundamental translation vectors on the complex plane C. We introduce the (0, 0)-cell as the square The square lattice Q consists of the cells where m 1 , m 2 run over integers.
The centers b k are considered as random variables distributed in such a way that the segments generate a set of uniformly distributed non-overlapping segments.Theoretically, this distribution can be introduced as the distribution of the variable b ).The separation parameter δ is taken equal to l 4 .This value refers to the minimum distance between the centers of the bacteria.The choice of the constant 4 is based on the experimentally observed widths of bacteria.It should be noted that the segments Γ k belong to the cell Q (0,0) in the double periodic topology when the opposite sides of Q (0,0) are glued by pairs.The random variable b can be statistically realised for large N by the Monte Carlo method to get numerical results.The following constructive procedure to generate random locations of segments is used.Let a random point b 1 is taken in accordance with the uniform distribution in Q (0,0) .Next, a random angle α 1 ∈ [−π, π] is chosen.Hence, the pair of points z 1 = b 1 ± l 2 e iα1 is constructed.This is equivalent to the construction of the first segment Γ 1 determined by the point b 1 and angle α 1 .At the next step, we take a random point b 2 uniformly distributed in Q\H 1 , where Further, a random angle α 2 is selected and we check whether the segments with the ends z 2 = b 2 ± l 2 e iα2 and z 1 do not intersect.If it is true, we have the second random segment determined by the point b 2 and angle α 2 .If the segments intersect, we take a random point b 2 again in Q\H 1 and randomly select a new random angle α 2 .In the same way, we take the next point 2 e iα N with the random angle α N .Introduce the density segments associated to the conformally invariant conductivity (capacity) [4] where l is the length of the segment and N the number of segments per representative cell.The algorithm described above generates a probability distribution U depending on the density.This distribution models the DB sets.In the limit case = 0, the distribution U 0 becomes the well known Poisson distribution of points in the square.
In order to study the distributions U for various densities we will use the following e-sums [9] ( Here, E p denotes the Eisenstein functions of order p (see Appendix).The values of the basic sums e 2 , e 22 , e 33 and e 44 will be estimated using the Monte Carlo method.

3.
Computation of e-sums for DB sets.The theoretical probabilistic distributions corresponding to disordered locations of bacteria are modeled in the previous section.Now, we propose an effective computational tool to properly describe U .Let M denote the number of simulated realizations of the unit cell (Monte Carlo experiments) with the random value b.The parameters M and N must be chosen sufficiently large in order to obtain the stable averaged value of the e-sums.Theoretically, it is possible to consider only one experiment (M = 1) and to take a huge number N .But computations are less expensive if M is large and N can be not so huge as in the case M = 1.
First, we determine the minimal N for which the distribution is characterized in that the angles of the segments will be uniformly distributed over the interval [−π, π].It is assumed that the segments directions are uniformly distributed if the following inequalities are fulfilled 1 N

N n=1
Re exp(iα n ) ≤ 0.15 and 1 These conditions are satisfied for N = 500 (see Fig. 2).In order to estimate M 0 = M when the computations become stable we consider the dependence of the mean e 44 on M for M ∈ [1,1500] experiments.All the sums e 2 and e pp (p = 2, 3, 4) are estimated, but e 44 is characterized by major volatility.
Here, e 44 is equal to the mean value of (e 44 ) m calculated for b generated in the m-th numerical experiment, more precisely, The results are shown in Fig. 3.One can observe that errors do not exceed 2% for M > 700.Therefore, the computations demonstrate that we have to simulate at The values e 2 , e 22 , e 33 , e 44 are computed for the distributions U for various densities .The sum e 2 must be equal to π [9] for ideal macroscopically isotropic random locations of segments.This yields the first criterion equation for the macroscopic isotropy of structures.The sums e pp (p = 2, 3, 4) describe high order basic terms of the distributions U .The next terms could describe more precisely U for higher densities [8].The average e-sums are selected in Table 1 for the densities changing from 0.05 to 0.6 with the step 0.05.The algorithm works too slowly for higher densities.The results from Table 1 can be extended by expensive computations to higher densities.
Table 1 contains the fundamental parameters of the uniform non-overlapping distribution U of segments on the plane.The simulated theoretical distribution describes DB sets when bacteria may affect each other but these interactions are local and do not yield the collective behavior.
4. Collective motion of bacteria.In the present section, we pay attention to experimental results partially presented in [18].The images of Bacillus subtilis in 31 frames are used in computations.One of the typical frame is displayed in Fig. 4. We use algorithms of image processing and analysis to determine number, centers, angles of inclinations and length of bacteria.The density of bacteria is calculated by formula (1) and oscillates around value 0.15.The results of the image processing and analysis are applied to computation of the values of the experimental e-sums.The results are selected in Table 2.
Table 2.The e-sums for 31 film frames of Bacillus subtilis.The first column contains the number of the film frame, the second column contains the number of bacteria N detected in the frame.The next columns show basic sums.

no. N
Re[e 2 ] e 22 e 33 e 44 1 In order to compare the distributions of DB sets with the distribution of bacteria we have made theoretical calculation of 31 samples.Because of the average number of bacteria in the frame is about 2050, the minimum distance between the centers of the bacteria is 1  4 of their length and density of bacteria equals about 0.15, the theoretical calculations have been carried out for the following parameters N = 2050, = 0.15 and δ = l 4 .The length l = 0.017108 is normalized to the normalized unit area of the cell.The results of the calculations are shown in Table 3.  2 and 3 for the observed and theoretical distributions of bacteria, we can see that values of the corresponding e-sums differ.Therefore, the observed locations of bacteria are not fitted with the simulated disordered locations.We can conclude that the behavior of the bacteria is not disordered, hence, we suggest that it is collective.An exception is the values of e 2 which is close to π ≈ 3.14 for the DB sets and the observed locations (see Table 4).This demonstrates averaged isotropy of the bacteria motion for the considered data.
In order to see that the theoretical and observed distributions are essentially different, we compare averaged basic sums (see Table 4) from Table 2 and Table 3. Table 4. Comparison of the averaged e-sums for the observed bacteria locations with the e-sums computed for the DB sets ( = 0.15) from Table 2 and Table 3 It is worth noting that changes of the parameters (δ, M etc) do not essentially impact onto deviations of the e-sums of DB sets.Moreover, the DB sets have significantly higher values e pp (p = 2, 3, 4) than the observed ones (c.f., the bold line in Table 1 and the averaged data from Table 4).
The above analysis of collective behavior can be considered as the first application of the RVE theory [6] which will be extended in the future.In particular, we plan to study dynamical parameters of the bacteria distributions in time using the data displayed in Fig. 5. Appendix.Following [6], [7] we present constructive formulae for the Eisenstain-Rayleigh sums S m the Eisenstain functions E m corresponding to square lattice Q (see section 2).The Eisenstein-Rayleigh lattice sums S m can be easily calculated through the rapidly convergent series (43) and recurrent formulae (44) from paper [7].

Figure 2 .
Figure 2. The real (circles) and imaginary (crosses) parts of the averaged directions for N = 500 and for the total number of distributions M = 1500 ( = 0.25).All absolute values do not exceed 0.15.

Figure 5 .
Figure 5.The values of e 44 for subsequent frames of the film.

Table 1 .
The averaged e-sums for various densities.

Table 3 .
The e-sums calculated for 31 samples of DB sets.The parameters of distribution are N = 2050, = 0.15 and δ = l 4 . .