 strength of the interaction.]{}

Given that a squirmer travels its own size in a time $\tau = \frac{\sigma}{v_s}$, we perform simulations from $1450$ up to $3000$ $\tau$. Once the system reaches steady state, in a time window between $1000$ and $2000$ $\tau$, we carry out a systematic analysis of the dynamics of the squirmer suspension, considering $\xi$ from 0.1 to 10 and $\beta$ from -3 up to 3.

\[sec:Tools\]Analysis Tools
===========================

In order to be able to establish the effect played by the anisotropy of the interaction between swimmers, we follow the procedure in Ref. [@SoftMatter17] and study the degree of aggregation as a function of $\beta$ and $\xi$. Whenever the system forms a steady-state cluster distribution, we explore the clusters morphology by computing their mean size, the cluster-size distribution, their radius of gyration, and both the polar and nematic order parameter. When clusters are present, we identify them following a distance criterion: at each time step, two particles belong to the same cluster of size $s$ whenever their distance is smaller than $r_{cl} = 1.75 \sigma$ (first minimum of the $g(r)$, see appendix \[Appendix:Clusters\]).

To calculate the cluster-size distribution $f(s)$ we apply a criterion based on [@chantalbins2012]: (i) We arbitrarily subdivide the range of $s$-values into intervals $\Delta s_i = (s_{i,max} - s_{i,min})$, where $n_i^t$ is the total number of clusters within each interval $\Delta s$; (ii) we assign the value $n_i=n^t_i/\Delta s_i$ to every $s$ within $\Delta s_i$, and compute the fraction of clusters of size $s$ as $f(s)=n_i/N_c$, where $N_c=\sum_i n_i \Delta s_i$ is the total number of clusters.

To characterise the clusters’ morphology, we compute the radius of gyration $$\label{eq:Rg}
R_g(s) = \sqrt{\sum_{i,j=1}^s \frac{\left(\vec{r}_i-\vec{r}_j\right)^2
strength of the interaction. ] { } 

 Given that a squirmer travels its own size in a meter $ \tau = \frac{\sigma}{v_s}$, we do simulations from $ 1450 $ up to $ 3000 $ $ \tau$. Once the system reach firm state, in a time windowpane between $ 1000 $ and $ 2000 $ $ \tau$, we carry out a systematic analysis of the moral force of the wiggler suspension, considering $ \xi$ from 0.1 to 10 and $ \beta$ from -3 up to 3. 

 \[sec: Tools\]Analysis cock 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 In order to be able to establish the effect play by the anisotropy of the interaction between swimmers, we follow the procedure in Ref.   [ @SoftMatter17 ] and analyze the degree of aggregation as a function of $ \beta$ and $ \xi$. Whenever the organization forms a steady - state cluster distribution, we research the clusters morphology by calculate their mean size, the cluster - size distribution, their radius of gyration, and both the arctic and nematic order parameter. When clusters are present, we identify them following a distance standard: at each clock time step, two particles belong to the same cluster of size $ s$ whenever their distance is smaller than $ r_{cl } = 1.75 \sigma$ (first minimum of the $ g(r)$, see appendix \[Appendix: Clusters\ ]). 

 To calculate the bunch - size distribution $ f(s)$ we use a criterion free-base on [ @chantalbins2012 ]: (i) We arbitrarily subdivide the range of $ s$-values into intervals $ \Delta s_i = (s_{i, max } - s_{i, min})$, where $ n_i^t$ is the entire number of clusters within each interval $ \Delta s$; (ii) we assign the value $ n_i = n^t_i/\Delta s_i$ to every $ s$ within $ \Delta s_i$, and calculate the fraction of clusters of size $ s$ as $ f(s)=n_i / N_c$, where $ N_c=\sum_i n_i \Delta s_i$ is the total number of clusters. 

 To characterise the clusters ’ morphology, we compute the radius of gyration $ $ \label{eq: Rg } 
 R_g(s) = \sqrt{\sum_{i, j=1}^s \frac{\left(\vec{r}_i-\vec{r}_j\right)^2
 stgength of the interactiok.]{}

Given that a sqoiemer tcavels jts own rize in a time $\tau = \frac{\sigme}{v_s}$, qe peeform simulations from $1450$ up to $3000$ $\nau$. Once rhe wtstem reaciss steady stafc, in c vime window beteeen $1000$ and $2000$ $\tau$, we carry mug c systematic analysis of the dynamicf of thr dquirmer suspegsiom, conapdtring $\xi$ from 0.1 to 10 and $\beta$ from -3 up to 3.

\[vec:Tools\]Analyxis Tools
===========================

In order to be abpe tl establish the efvect played by eye anisotroph of the interaction bgtween swimmers, we follow the pruceduxe in Ref. [@SodtNathgr17] and study the qegree of agngegatiot as a gunction of $\beba$ anv $\xi$. Whenever the system horms a steady-state sluster dhscribution, we explore rhw cluvterv mofphouogg uy domputlng their meah size, the xluster-size distribityin, their radihs of dywation, and both the polar and nematic ogder parameter. When clusterw are present, we idenjify them sollowing a distance criterion: at each time step, dwo pergicoef vepong to the same cluster of size $s$ whenever trsit cistance is smcller than $r_{cl} = 1.75 \slgkw$ (first minimom of tks $f(r)$, see appendix \[Apoendix:Cjustees\]).

To calctlatr the cluster-size distributuon $f(s)$ we apily q criterion based ln [@chantalbnns2012]: (i) Ee arnitrarily subdivide the ranfe of $s$-valufs into ihgervals $\Delta s_i = (s_{p,max} - s_{i,min})$, where $n_i^t$ is the tjtal numbxr of clustefs wothin qach intergal $\Dcnta s$; (ii) we assign the ralue $n_i=n^t_i/\Deltw s_i$ to every $s$ within $\Delta s_i$, and compute tne frdction oy clusbers of size $s$ ws $f(s)=n_i/N_c$, wherg $N_c=\sum_i u_i \Delga s_i$ is tge totan number of clusters.

To cvwracterise tie clustews’ mirphilogy, wd compute the rsdius of yvration $$\lqbel{eq:Rg}
R_g(s) = \sqrt{\smm_{i,j=1}^s \rrac{\left(\vec{r}_i-\vee{x}_j\tught)^2
strength of the interaction.]{} Given that a its size in time $\tau = $1450$ to $3000$ $\tau$. the system reaches state, in a time window between and $2000$ $\tau$, we carry out a systematic analysis of the dynamics of squirmer suspension, considering $\xi$ from 0.1 to 10 and $\beta$ from -3 up 3. Tools In to be able to establish the effect played by the anisotropy of the interaction between swimmers, follow the procedure in Ref. [@SoftMatter17] and study degree of aggregation as function of $\beta$ and $\xi$. the forms a cluster we the clusters morphology computing their mean size, the cluster-size distribution, their radius of gyration, and both the polar and nematic parameter. When present, we them a criterion: at each two particles belong to the same $s$ whenever their distance is smaller than $r_{cl} 1.75 \sigma$ minimum of the $g(r)$, see appendix To calculate the cluster-size distribution $f(s)$ we apply criterion based on [@chantalbins2012]: (i) We arbitrarily subdivide the range of $s$-values into intervals $\Delta (s_{i,max} - s_{i,min})$, where is the total of within interval s$; (ii) assign the value $n_i=n^t_i/\Delta s_i$ to every $s$ within $\Delta s_i$, compute the fraction of clusters of size $s$ as $f(s)=n_i/N_c$, n_i s_i$ is the number of clusters. To the morphology, we compute the gyration R_g(s)
 strength of the interaction.]{}

GIven that a sQuirmEr tRavElS its Own sIze in a time $\tau = \fRAc{\siGma}{v_s}$, we perform simulatiOns frOm $1450$ UP to $3000$ $\tAU$. ONce thE system REaCHEs sTeAdY stAtE, In A time WinDow betwEen $1000$ and $2000$ $\tau$, wE caRrY out a systemaTIc Analysis of The Dynamics of thE sqUirmer SuSpeNSion, cOnsIderiNg $\xi$ frOM 0.1 to 10 and $\Beta$ from -3 uP tO 3.

\[Sec:TooLS\]AnalysIS toOls
===========================

IN order to be able to eSTaBLish the effect pLayed bY tHE aNISotRopY of the inteRaCtion BEtween sWImMERS, we FOllow the proceDure in Ref. [@SoFTMaTter17] anD sTudY The degRee of AgGRegAtion as a funCtioN of $\beta$ anD $\xi$. WheNEver the SYstem foRms a stEadY-stAte cLUsTeR diStRIbuTIoN, we EXplOre the clUsTeRs morPholOGY BY comPutIng tHeir mEan size, the cluSteR-sizE DisTribuTion, tHeir RaDius oF gyratIon, anD bOth the polar and nEmatIc order paRamEtEr. WHeN clusTErs are PreSenT, we idenTify theM FolLoWING a Distance criterion: aT eACH tIme step, tWo partIClEs BElong to tHe SamE cluSTEr of sIze $s$ WHeNever theIr distANcE iS smalleR tHan $r_{cl} = 1.75 \SiGma$ (FirSt minIMum oF the $g(r)$, See appenDix \[ApPEndix:Clusters\]).

TO Calculate the cLUsTER-sIZe diStrIbution $f(s)$ we ApplY A criTeriON bAseD On [@chaNtalbInS2012]: (I) WE Arbitrarily subdividE tHe rangE of $s$-vAlues into inteRvals $\Delta S_I = (S_{I,max} - s_{i,miN})$, wheRE $n_I^T$ is the total numBer of Clusters wiTHin each iNtervAl $\Delta s$; (Ii) we assigN THe value $n_I=n^t_I/\DeLta S_i$ tO EVeRy $s$ within $\DeltA S_I$, and CoMpute thE frAction oF clUstErs Of sIzE $s$ as $f(s)=n_i/N_C$, where $N_c=\SuM_i N_i \deLta S_i$ is tHE total nuMbEr oF cLusTers.

TO CharacTerisE the ClUsTErs’ MorpholOGy, WE CompUtE tHe raDiuS oF gyraTion $$\LAbeL{eq:Rg}
R_g(S) = \sqrt{\sum_{i,J=1}^s \fRAc{\leFt(\VeC{r}_i-\vec{r}_J\right)^2
 strength of the interacti on.]{}

Gi ven t hat  asq uirm er t ravels its own size  in a time $\tau = \fr ac{\s ig m a}{v _ s} $, we  perfor m  s i m ula ti on s f ro m  $ 1450$  up  to $30 00$ $\tau$ . O nc e the system re aches stea dystate, in atim e wind ow  be t ween$10 00$ a nd $20 0 0$ $\t au$, we c ar r y outa  system a t ic  ana lysis of the dyna m ic s  of the squirm er sus pe n si o n , c ons idering $\ xi $ fro m  0.1 to 10 a n d $ \ beta$ from -3  up to 3.
\ [se c:Tool s\ ]An a lysisTools 
= = === =========== ==== =======
In ord e r to be able to  estab lis h t he e f fe ct  pl ay e d b y  t hea nis otropy o fth e int erac t i o n  bet wee n sw immer s, we followthe  pro c edu re in  Ref.  [@S of tMatt er17]and s tu dy the degree o f ag gregation  as  a  fu nc tiono f $\be ta$  an d $\xi$ . Whene v erth e s ys tem forms a steady -s t a te  cluster  distr i bu ti o n, we ex pl ore  the c luste rs m o rp hology b y comp u ti ng  theirme an siz e,  th e c luste r -siz e dist ribution , the i r radius of gy r ation, and bo t ht h ep olar  an d nematic o rder para mete r .Whe n  clus tersar e  p r esent, we identifyth em fol lowin g a distancecriterion: a t  each ti me s t ep ,  two particles  belo ng to thes ame clus ter o f size $ s$ whenev e r  their d ist anc e i s s m a ll er than $r_{c l }  = 1 .7 5 \sigm a$(firstmin imu m o f t he  $g(r)$,see appe nd ix  \ [A ppe ndix: C lusters\ ]) .
To  ca lcula t e theclust er-s iz ed ist ributio n  $ f ( s)$we  a pply  acr iteri on b a sed  on [@c hantalbin s20 1 2]:(i )We arbi trarily subdi vi de the ran ge  of  $s$-v a l ues into  intervals $\Delta s_i=  (s_{i, max } - s _{i, min})$, w her e $n_i ^t$ is the  total  numb er  of c luste r s  w ith in  each inte r v al$\Del ta  s$;  (ii) w e assign the value $n_ i=n^t_i/\Delt a s _i$t o  e ver y  $ s $ w it h in$ \ Delta s_i$, and  compute t he fr action ofc lus te rs of s ize $s$  as $ f (s)=n_i /N_c$, wh ere $N_c= \s um_i n _i\Delta s_i $ is the  total nu m ber o f  c luste rs. 

To c ha rac teris e thec lus ters’  morph ol ogy, w e com pu te the r adius of gyration $$\la bel{eq :Rg}R_g (s) = \sq rt{ \ sum _{i,j=1}^ s \f rac{\left( \ve c{r }_i-\ vec { r}_j\ righ t )^ 2
 strength_of the_interaction.]{}

Given that a squirmer_travels its_own_size in_a_time $\tau =_\frac{\sigma}{v_s}$, we perform_simulations from $1450$ up_to $3000$ $\tau$._Once_the system reaches steady state, in a time window between $1000$ and $2000$ $\tau$,_we_carry out_a_systematic_analysis of the dynamics of_the squirmer suspension, considering $\xi$_from 0.1_to 10 and $\beta$ from -3 up to_3.

\[sec:Tools\]Analysis_Tools
===========================

In order to_be able to establish the effect played by the_anisotropy of the interaction between swimmers,_we follow the_procedure_in_Ref. [@SoftMatter17] and study the_degree of aggregation as a function_of $\beta$ and $\xi$. Whenever the_system forms a steady-state cluster distribution, we_explore the clusters morphology by computing_their mean size, the cluster-size_distribution, their_radius of gyration, and both_the polar and_nematic order_parameter. When clusters_are present, we identify them following_a distance criterion:_at each time step, two particles_belong_to the same_cluster_of_size $s$_whenever their distance_is_smaller than_$r_{cl}_= 1.75 \sigma$ (first minimum of_the_$g(r)$, see appendix \[Appendix:Clusters\]).

To calculate the cluster-size_distribution $f(s)$ we apply_a_criterion based on [@chantalbins2012]:_(i) We arbitrarily subdivide the_range of $s$-values into intervals $\Delta_s_i =_(s_{i,max} -_s_{i,min})$, where $n_i^t$ is the total number of clusters within each_interval $\Delta s$; (ii) we assign_the value $n_i=n^t_i/\Delta s_i$_to every_$s$_within $\Delta s_i$,_and_compute the_fraction of clusters of size $s$ as_$f(s)=n_i/N_c$, where_$N_c=\sum_i n_i \Delta s_i$ is the_total number of clusters.

To_characterise_the clusters’ morphology, we compute the_radius of gyration $$\label{eq:Rg}
R_g(s) = \sqrt{\sum_{i,j=1}^s_\frac{\left(\vec{r}_i-\vec{r}_j\right)^2