Modelling pattern formation through differential repulsion

Scheme of the predictions of the linear stability analysis by acting on the interspecies repulsion force

Form of the potential $ \tilde{\Phi}^{ST}(x) $ for $ R = 2 $ and $ \frac{\nu^{ST}_c}{\nu^{ST}_d} \frac{\kappa^{ST}}{2} = 1 $. The potential is repulsive on its support

Values of $ \lambda_1 $ (blue curve), $ \lambda_2 $ (orange curve) and their mean (yellow dotted line) near $ z = 0 $ (z = 0.1), plotted as functions of the scaling parameter $ s $ for $ R = 1 $, $ D^A = D^B = 1, c'^{AA} = c'^{AB} = c'^{BA} = 1, c'^{BB} = 10 $

(I): Values of $ \lambda_1(z) $ as functions of $ z $ for $ R = 1 $, $ D^A = D^B = 1, \, c^{AA} = c^{AB} = c^{BA} = 1, \, c^{BB} = 10 $ and for different values of $ s $ in the instability regime: $ s = 30 $ (blue curve), $ s = 50 $ (orange curve), $ s = 70 $ (yellow curve), and $ s = 90 $ (red curve). (II): same plots for $ \lambda_2(z) $. (III) Plot of $ z^* $ defined in (36) as a function of parameter $ s $

Microscopic simulations for Cases 1-4 for parameters described in Table 1. $ A $-cells are represented as red disks, $ B $-cells as green disks. For each subsection, the left figure is obtained starting from a homogeneous distribution of particles, the right one from a segregated initial distribution ($ B $-cells on the half-left of the domain, $ A $-cells on the right)

Macroscopic simulations for Cases 1-4 for parameters described in table 1 for the final time of simulations equal to $ T = 8000 $

(I) Simulations of the microscopic and macroscopic models for $ \kappa^{AA} = 4, \kappa^{BB} = \tilde{\kappa}^{AB} = \tilde{\kappa}^{BA} = 1 $ for which $ s^* \approx 2.1 $. Simulations of the microscopic model are performed with $ N_A = N_B = 500 $ (IA) and $ N_A = N_B = 2000 $ (IB) particles. Simulations of the macroscopic model correspond to (IC). We consider 7 values of the interspecies repulsion intensity, from left to right: for $ 2.05 = s<s^* $, for $ s^*<s = \{2.15, 2.2, 2.5, 4, 6, 10 \} $. Type B cells are represented in green, type A cells in red. (II) In (IIA-C), we show the values of the quantifiers at time equilibrium as functions of $ s $. Fig. (IA) shows the mean elongation of the green clusters, (IIB) shows the number of green clusters and (IIC) shows the overlapping amount $ Q $ described by (73). Black curves are obtained with the microscopic model for $ N_A = N_B = 500 $ (corresponding to Figures (IA)), yellow curves are for $ N_A = N_B = 2000 $ and correspond to Figures (IB) and red curves are obtained with the macroscopic model (Figures (IC)). The two bottom figures correspond to a zoom of the corresponding curves close to the transition region for $ s $

Quantifiers computed on the simulation images as functions of the logarithm of the simulation time for $ \kappa^{AA} = 4, \kappa^{BB} = \tilde{\kappa}^{AB} = \tilde{\kappa}^{BA} = 1 $ for the microscopic model with $ N_A = N_B = 500 $ (green curves), $ N_A = N_B = 2000 $ (blue curves), $ N_A = N_B = 4000 $ (yellow curves) and for the macro model (red curve). (I) Green cluster elongation, (II) Number of cell clusters and (III) Overlapping amount $ Q $. On Figure (I), we superimpose linear fits (dotted lines) for short times and large times, showing the two timescales (two slopes) of the micro model compared to the unique timescale of the macro dynamics (single slope)