 the asymptotic behavior in rapid flows $\dot{\gamma}\gg\dot{\gamma}_c$) [@diff_shear_md2; @diff_shear_md3; @diff_shear_md4]. In addition, it was found in MD simulations of soft athermal disks that in a sufficiently small flow rate range, the diffusivity changes from $D\sim\dot{\gamma}$ ($\phi<\phi_J$) to $\dot{\gamma}^{0.78}$ ($\phi\simeq\phi_J$) [@diff_shear_md0], implying that the crossover shear rate $\dot{\gamma}_c$ vanishes as the system approaches jamming from below $\phi\rightarrow\phi_J$.

Note that the self-diffusion of soft athermal particles shows a clear difference from the diffusion in glass; *no plateau* is observed in (transverse) mean square displacements (MSDs) [@diff_shear_md0; @diff_shear_md2; @diff_shear_md3; @diff_shear_md4; @diff_shear_md7; @dh_md2]. The absence of sub-diffusion can be also seen in quasi-static simulations ($\dot{\gamma}\rightarrow 0$) of soft athermal disks [@dh_qs1] and MD simulations of granular materials sheared under constant pressure [@diff_shear_md1].

Because the self-diffusion can be associated with collective motions of soft athermal particles, researchers have analyzed spatial correlations of velocity fluctuations [@rheol0] or non-affine displacements [@nafsc2] of the particles under shear. Characteristic sizes of collectively moving regions, i.e. *rigid clusters*, are then extracted as functions of $\dot{\gamma}$ and $\phi$, however, there is a lack of consensus on the scaling of the sizes. For example, the size of rigid clusters $\xi$ diverges as the shear rate goes to zero $\dot{\gamma}\rightarrow 0$ so that the power-law scaling $\xi\sim\dot{\gamma}^{-s}$ was suggested, where the exponent varies from $s=0.23$ to $0.5$ depending on numerical models and flow conditions [@dh_md2; @diff_shear_md1]. The dependence of the rigid cluster size on packing fraction is also controversial. If the system is below
the asymptotic behavior in rapid flows $ \dot{\gamma}\gg\dot{\gamma}_c$) [ @diff_shear_md2; @diff_shear_md3; @diff_shear_md4 ]. In addition, it was found in MD model of piano athermal disks that in a sufficiently small flow pace range, the diffusivity changes from $ D\sim\dot{\gamma}$ ($ \phi<\phi_J$) to $ \dot{\gamma}^{0.78}$ ($ \phi\simeq\phi_J$) [ @diff_shear_md0 ], implying that the crossing over shear rate $ \dot{\gamma}_c$ vanishes as the system approach jamming from below $ \phi\rightarrow\phi_J$. 

 Note that the self - dispersion of soft athermal particles testify a clear difference from the diffusion in glass; * no plateau * is observed in (transverse) entail square displacements (MSDs) [ @diff_shear_md0; @diff_shear_md2; @diff_shear_md3; @diff_shear_md4; @diff_shear_md7; @dh_md2 ]. The absence of sub - dispersion can be also seen in quasi - static simulations ($ \dot{\gamma}\rightarrow 0 $) of delicate athermal disks [ @dh_qs1 ] and MD simulations of granular materials sheared under changeless pressure [ @diff_shear_md1 ]. 

 Because the self - diffusion can be associated with collective motions of soft athermal particles, researchers have analyzed spatial correlations of velocity fluctuations [ @rheol0 ] or non - affine displacements [ @nafsc2 ] of the particles under shear. Characteristic sizes of collectively affect regions,   i.e.   * rigid bunch *, are then educe as function of $ \dot{\gamma}$ and $ \phi$, however, there is a lack of consensus on the scaling of the sizes. For example, the size of rigid clusters $ \xi$ diverges as the shear rate goes to zero $ \dot{\gamma}\rightarrow 0 $ so that the baron - law scaling $ \xi\sim\dot{\gamma}^{-s}$ was suggested, where the exponent varies from $ s=0.23 $ to $ 0.5 $ depending on numeral models and flow conditions [ @dh_md2; @diff_shear_md1 ]. The dependence of the inflexible cluster size on packing material fraction is besides controversial. If the arrangement is below
 thf asymptotic behavior in rapid flows $\doj{\gqmma}\gg\vot{\gammz}_c$) [@diff_sfear_md2; @diff_shear_md3; @diff_shear_ld4]. In aedition, it was found iv MD simupations if sidt athermal disks bkat ih a sbfhiciently small flow rate range, the difxuririty changes from $D\sim\dot{\gamma}$ ($\phi<\phy_J$) to $\dpt{\hamma}^{0.78}$ ($\phi\simeq\khi_J$) [@cyff_sgvav_md0], implying that the crossover ahear rete $\dot{\gamma}_c$ vsnishes as the system apprlachfs jamming from bepow $\phi\righjzrrjq\phi_J$.

Note thxt the self-diffusion or soft athermal particles shows a clzar differebcw fgmm the diffnsion pn glass; *no plateau* iv obserfed in (transvevse) mxan wquare displacements (KSDs) [@diff_shear_md0; @dyff_shear_mg2; @biff_shear_md3; @diff_shear_nd4; @diff_vheas_md7; @ay_md2]. Tht ausehce of suu-diffusion dan be also seen in quasi-stativ fpkulations ($\dof{\gamma}\widhtarrow 0$) of soft athermal disks [@dh_qs1] atd JD simulations of granuoar materials sheared under cogstant pressure [@diff_shear_md1].

Because the self-diffushon cen be africlated with collective motions of soft athermaj psrnicles, researcherf have analuzfd fpatial correuations or velocity fluctuahions [@rreol0] ir non-affyne cisplacements [@nafsc2] of the particles unber shear. Characteriscic sizes of colkectifely moving regions, i.e. *riyid clhsters*, are hhen extrzzted as functionr on $\dmt{\gamma}$ and $\phi$, however, thqre is a oack of conrensos on tre scaling of tma sizes. For examplf, the shze of rigld clusters $\xi$ diverges as the shear rate goex do dero $\dot{\gcmma}\rinhtarrow 0$ so thwt the power-lac scaliny $\xi\sio\dot{\gamma}^{-s}$ was sujgested, wherq the exponend varies from $s=0.23$ to $0.5$ dqpeneing on numdfical models amd flow conditions [@eh_md2; @diff_shear_md1]. Tme dekehdence of the rnyie cluster size pn oachijg fwdction is alvo cuntfpverskal. If uke whstek is below
the asymptotic behavior in rapid flows $\dot{\gamma}\gg\dot{\gamma}_c$) @diff_shear_md4]. addition, it found in MD that a sufficiently small rate range, the changes from $D\sim\dot{\gamma}$ ($\phi<\phi_J$) to $\dot{\gamma}^{0.78}$ [@diff_shear_md0], implying that the crossover shear rate $\dot{\gamma}_c$ vanishes as the system approaches from below $\phi\rightarrow\phi_J$. Note that the self-diffusion of soft athermal particles shows a difference the in *no plateau* is observed in (transverse) mean square displacements (MSDs) [@diff_shear_md0; @diff_shear_md2; @diff_shear_md3; @diff_shear_md4; @diff_shear_md7; @dh_md2]. absence of sub-diffusion can be also seen in simulations ($\dot{\gamma}\rightarrow 0$) of athermal disks [@dh_qs1] and MD of materials sheared constant [@diff_shear_md1]. the self-diffusion can associated with collective motions of soft athermal particles, researchers have analyzed spatial correlations of velocity fluctuations [@rheol0] non-affine displacements the particles shear. sizes collectively moving regions, clusters*, are then extracted as functions $\phi$, however, there is a lack of consensus the scaling the sizes. For example, the size rigid clusters $\xi$ diverges as the shear rate to zero $\dot{\gamma}\rightarrow 0$ so that the power-law scaling $\xi\sim\dot{\gamma}^{-s}$ was suggested, where the exponent $s=0.23$ to $0.5$ depending numerical models and conditions @diff_shear_md1]. dependence the rigid size on packing fraction is also controversial. If the system is
 the asymptotic behavior in raPid flows $\doT{\gammA}\gg\Dot{\GaMma}_c$) [@Diff_Shear_md2; @diff_sheAR_md3; @dIff_shear_md4]. In addition, it Was foUnD In MD SImUlatiOns of soFT aTHErmAl DiSks ThAT iN a sufFicIently sMall flow raTe rAnGe, the diffusiVItY changes frOm $D\Sim\dot{\gamma}$ ($\pHi<\pHi_J$) to $\dOt{\GamMA}^{0.78}$ ($\phi\sImeQ\phi_J$) [@Diff_shEAr_md0], imPlying thaT tHE crossOVer sheaR RAtE $\dot{\Gamma}_c$ vanishes as tHE sYStem approaches JamminG fROm BELow $\Phi\Rightarrow\PhI_J$.

NotE That the SElF-DIFfuSIon of soft atheRmal particlES shOws a clEaR diFFerencE from ThE DifFusion in glaSs; *no Plateau* is ObservED in (tranSVerse) meAn squaRe dIspLaceMEnTs (mSDS) [@dIFf_sHEaR_md0; @DIff_Shear_md2; @dIfF_sHear_mD3; @difF_SHEAr_md4; @DifF_sheAr_md7; @dH_md2]. The absence Of sUb-diFFusIon caN be alSo seEn In quaSi-statIc simUlAtions ($\dot{\gamma}\rIghtArrow 0$) of soFt aThErmAl Disks [@DH_qs1] and mD sImuLations Of granuLAr mAtERIAlS sheared under constAnT PReSsure [@difF_shear_MD1].

BEcAUse the seLf-DifFusiON Can be AssoCIaTed with cOllectIVe MoTions of SoFt atheRmAl pArtIcles, REseaRchers Have analYzed sPAtial correlatiONs of velocity fLUcTUAtIOns [@rHeoL0] or non-affinE disPLaceMentS [@NaFsc2] OF the pArticLeS UnDEr shear. CharacteristIc Sizes oF collEctively movinG regions, i.e. *RIGId clusteRs*, arE ThEN extracted as fuNctioNs of $\dot{\gamMA}$ and $\phi$, hOweveR, there is A lack of coNSEnsus on tHe sCalIng Of tHE SiZes. For example, THE sizE oF rigid cLusTers $\xi$ dIveRgeS as The ShEar rate goEs to zero $\DoT{\gAmMa}\RigHtarrOW 0$ so that tHe PowEr-Law ScaliNG $\xi\sim\Dot{\gaMma}^{-s}$ WaS sUGgeSted, wheRE tHE ExpoNeNt VariEs fRoM $s=0.23$ to $0.5$ dEpenDIng On numerIcal modelS anD Flow CoNdItions [@dH_md2; @diff_shear_mD1]. THe dependenCe Of tHe rigiD CLuster siZe on packing fraction is alSO controVerSial. IF the System is bEloW
 the asymptotic behavior i n rapid fl ows $ \do t{\ ga mma} \gg\ dot{\gamma}_c$ )  [@d iff_shear_md2; @diff_s hear_ md 3 ; @d i ff _shea r_md4]. In a ddi ti on , i tw as  foun d i n MD si mulationsofso ft athermald is ks that in  asufficiently  sm all fl ow  ra t e ran ge,  thediffus i vity c hanges fr om $D\sim \ dot{\ga m m a} $ ($ \phi<\phi_J$) to$ \d o t{\gamma}^{0.7 8}$ ($ \p h i\ s i meq \ph i_J$) [@di ff _shea r _md0],i mp l y i ngt hat the cross over shearr ate  $\dot {\ gam m a}_c$vanis he s  as  the system  app roaches j amming from be l ow $\ph i\righ tar row \phi _ J$ .
No te tha t  t hes elf -diffusi on  o f sof t at h e r m al p art icle s sho ws a clear di ffe renc e  fr om th e dif fusi on  in g lass;*no p la teau* is observ ed i n (transv ers e)  me an  squa r e disp lac eme nts (MS Ds) [@d i ff_ sh e a r _m d0; @diff_shear_md 2; @ di ff_shear _md3;@ di ff _ shear_md 4;  @d iff_ s h ear_m d7;@ dh _md2]. T he abs e nc eof sub- di ffusio ncan  be  also seen  in qu asi-stat ic si m ulations ($\do t {\gamma}\righ t ar r o w0 $) o f s oft atherma l di s ks [ @dh_ q s1 ] a n d MDsimul at i on s  of granular materi al s shea red u nder constant  pressure[ @ d iff_shea r_md 1 ]. 
Because the s elf-d iffusion c a n be ass ociat ed withcollectiv e motionsofsof t a the r m al  particles, r e s earc he rs have  an alyzedspa tia l c orr el ations of  velocit yfl uc tu ati ons [ @ rheol0]or  no n- aff ine d i splace ments  [@n af sc 2 ] o f the p a rt i c lesun de r sh ear .Chara cter i sti c sizes  of colle cti v elymo vi ng regi ons, i.e. *ri gi d clusters *,  ar e then e xtracted  as functions of $\dot{ \ gamma}$  an d $\p hi$,  however,  th ere is  al ack of  conse nsuson  th e scali n g  o f t he  sizes. Fo r exa mple,  t he s ize ofrigid clusters $\x i $ d iverges as th e s hear r at e g o es toze r o $ \ d ot{\gamma}\righ tarrow 0$so th at the pow e r-l aw  scalin g $\xi\ sim\d o t{\gamm a}^{-s}$was sugge st ed,w h ere  the expon ent vari es from $ s =0.23 $  t o $0. 5$depend in g o n num erical mod els a nd flo wcondit ions[@ dh_md2;@diff_shear_md1]. The d epende nce o f t he rigidclu s ter  size onpack ing fracti onisalsocon t rover sial .  I f t h e sys temi s below
 the_asymptotic behavior_in rapid flows $\dot{\gamma}\gg\dot{\gamma}_c$)_[@diff_shear_md2; @diff_shear_md3;_@diff_shear_md4]._In addition,_it_was found in_MD simulations of_soft athermal disks that_in a sufficiently_small_flow rate range, the diffusivity changes from $D\sim\dot{\gamma}$ ($\phi<\phi_J$) to $\dot{\gamma}^{0.78}$ ($\phi\simeq\phi_J$) [@diff_shear_md0], implying_that_the crossover_shear_rate_$\dot{\gamma}_c$ vanishes as the system_approaches jamming from below $\phi\rightarrow\phi_J$.

Note_that the_self-diffusion of soft athermal particles shows a clear_difference_from the diffusion_in glass; *no plateau* is observed in (transverse) mean_square displacements (MSDs) [@diff_shear_md0; @diff_shear_md2; @diff_shear_md3;_@diff_shear_md4; @diff_shear_md7; @dh_md2]._The_absence_of sub-diffusion can be_also seen in quasi-static simulations ($\dot{\gamma}\rightarrow_0$) of soft athermal disks [@dh_qs1]_and MD simulations of granular materials sheared_under constant pressure [@diff_shear_md1].

Because the self-diffusion_can be associated with collective_motions of_soft athermal particles, researchers have_analyzed spatial correlations_of velocity_fluctuations [@rheol0] or_non-affine displacements [@nafsc2] of the particles_under shear. Characteristic_sizes of collectively moving regions, i.e. *rigid clusters*,_are_then extracted as_functions_of_$\dot{\gamma}$ and_$\phi$, however, there_is_a lack_of_consensus on the scaling of the_sizes._For example, the size of rigid clusters_$\xi$ diverges as the_shear_rate goes to zero_$\dot{\gamma}\rightarrow 0$ so that the_power-law scaling $\xi\sim\dot{\gamma}^{-s}$ was suggested, where_the exponent_varies from_$s=0.23$ to $0.5$ depending on numerical models and flow conditions [@dh_md2;_@diff_shear_md1]. The dependence of the rigid_cluster size on packing_fraction is_also_controversial. If the_system_is below