eq:D1\])\], is caused by the diverging length scale below jamming, i.e. $\xi\sim|\Delta\phi|^{-1}$ \[Eq. (\[eq:xi1\])\]. As shown in Fig. \[fig:xi\](b), the linear relation (dotted line) well explains our results if the length scale $\xi$ is smaller than $10d_0$. If the system is above jamming, the length scale increases (more than $10d_0$) with the increase of $\phi$. However, the diffusivity over the shear rate $D/\dot{\gamma}$ starts to deviate from the linear relation (dotted line) and the length scale reaches the system size $L/2\simeq 44d_0$ (shaded region). We conclude that this deviation is caused by finite-size effects and further studies of different system sizes are necessary (as in Refs. [@diff_shear_md3; @diff_shear_md4]) to figure out the relation between $D/\dot{\gamma}$ and $\xi$ in this regime, which we postpone as a future work.

Discussions {#sec:disc}
===========

In this study, we have numerically investigated rheological and transport properties of soft athermal particles under shear. Employing MD simulations of two-dimensional disks, we have clarified how the rheology, self-diffusion, and size of rigid clusters vary with the control parameters,i.e. the externally imposed shear rate $\dot{\gamma}$ and packing fraction of the disks $\phi$. Our main result is the critical scaling of the diffusivity (Sec. \[sub:diff\]) and size of rigid clusters (Sec. \[sub:rigid\]), where their dependence on both $\dot{\gamma}$ and $\phi$ is reported \[Eqs. (\[eq:D1\]), (\[eq:D2\]), (\[eq:xi1\]), and (\[eq:xi2\])\]. The diffusivity has been calculated on both sides of jamming (by a single numerical protocol) to unify the understanding of self-diffusion of soft particulate systems: We found that (i) the diffusivity below jamming exhibits a crossover from the linear scaling $D\sim\dot{\gamma}$ to the power-law $D\sim\dot{\gamma}^{3/4}$. Such a crossover can be also seen in previous simulations [@diff_she
eq: D1\])\ ], is caused by the diverging length scale below jamming,   i.e.   $ \xi\sim|\Delta\phi|^{-1}$ \[Eq.   (\[eq: xi1\])\ ]. As testify in Fig.   \[fig: xi\](b), the analogue relation (dotted course) well explain our results if the length plate $ \xi$ is smaller than $ 10d_0$. If the system is above jamming, the length plate increase (more than $ 10d_0 $) with the increase of $ \phi$. However, the diffusivity over the shear pace $ D/\dot{\gamma}$ starts to deviate from the linear relation (dotted lineage) and the length scale reaches the system size $ L/2\simeq 44d_0 $ (shaded region). We reason that this deviation is caused by finite - size effects and further study of different system size are necessary (as in Refs.   [ @diff_shear_md3; @diff_shear_md4 ]) to figure out the sexual intercourse between $ D/\dot{\gamma}$ and $ \xi$ in this regime, which we postpone as a future work. 

 Discussions { # sec: disc } 
 = = = = = = = = = = = 

 In this study, we have numerically investigate rheological and transport properties of soft athermal particles under shear. Employing MD simulations of two - dimensional disks, we have clarified how the rheology, self - diffusion, and size of rigid clusters vary with the control parameters, i.e.   the externally imposed shear pace $ \dot{\gamma}$ and tamp down fraction of the disks $ \phi$. Our independent resultant role is the critical scaling of the diffusivity (Sec.   \[sub: diff\ ]) and size of rigid clusters (Sec.   \[sub: rigid\ ]), where their dependence on both $ \dot{\gamma}$ and $ \phi$ is report \[Eqs.   (\[eq: D1\ ]), (\[eq: D2\ ]), (\[eq: xi1\ ]), and (\[eq: xi2\])\ ]. The diffusivity has been calculate on both sides of jamming (by a single numeric protocol) to unify the understanding of self - diffusion of soft particulate systems: We found that (i) the diffusivity below jamming exhibits a crossover from the linear scaling $ D\sim\dot{\gamma}$ to the power - jurisprudence $ D\sim\dot{\gamma}^{3/4}$. Such a crossover can be besides watch in previous simulations [ @diff_she
eq:D1\])\], is caused by the diverglng length scale below oamming, j.e. $\xi\sim|\Ddlta\phi|^{-1}$ \[Eq. (\[eq:xi1\])\]. As shown in Fmg. \[fit:xi\](b), uke linear relation (dutted linv) well explaiis our results ih the lekyth sdwle $\ri$ is smaller thsn $10d_0$. If tha system is abmvd lamming, the length scale increases (mjre tham $10f_0$) with the inctease jf $\pgp$. Mowever, the diffusivity over the shear gate $D/\dot{\gamma}$ stsrts to deviate from the llneag relation (dotted pine) and thg legtth scale rexches the system size $M/2\simeq 44d_0$ (shaded region). We concljde tkat this deciqtilt is caused by fpnite-size effcbts and furthet studies of dlfferxnt wystem sizes are necevsary (as in Refs. [@disf_shear_md3; @dnff_shear_md4]) to figure iur the reldtiov begwetn $V/\dof{\gamma}$ anv $\xi$ in thia regime, whuch we postpone as s slyure work.

Disdussiogs {#sec:disc}
===========

In this study, we have numericalny jnvestigated rheologicao and transport propetties of sjft athermal particles under shear. Employing MD shmulaviuns on twu-eilensional disks, we have clarified how the rhejmoby, self-diffusion, and size og gibyd clusters vxry wicg fhe control paramehers,i.e. tre exrernally ympoxed shear rate $\dot{\gamma}$ and packing frabtiob of the disks $\phi$. Our main rzsult os thr critical scaling of tke difrusivity (Sef. \[sub:diff\]) zvd size of rigid cllstess (Sec. \[sub:rigid\]), where their dependenre on both $\dut{\gakma}$ anq $\phi$ is rfportcg \[Eqs. (\[eq:D1\]), (\[eq:D2\]), (\[eq:xi1\]), wnd (\[ee:xh2\])\]. The diffksivity has been calculated on upth sides of jdmmpng (by a finglc numerical proeocol) to unify the unberstavding of svlf-diffusmon of soft [articulate sfdtems: We fouid that (i) the difdusivith below jamming exhibits a crossovwr from the linear scxming $D\sim\dot{\gamna}$ uo the power-law $C\sio\doe{\gwmka}^{3/4}$. Vuch a crossmver cav be auso seen in orevoous simulations [@difx_she
eq:D1\])\], is caused by the diverging length jamming, $\xi\sim|\Delta\phi|^{-1}$ \[Eq. As shown in (dotted well explains our if the length $\xi$ is smaller than $10d_0$. If system is above jamming, the length scale increases (more than $10d_0$) with the of $\phi$. However, the diffusivity over the shear rate $D/\dot{\gamma}$ starts to deviate the relation line) the length scale reaches the system size $L/2\simeq 44d_0$ (shaded region). We conclude that this deviation caused by finite-size effects and further studies of system sizes are necessary in Refs. [@diff_shear_md3; @diff_shear_md4]) to out relation between and in regime, which we as a future work. Discussions {#sec:disc} =========== In this study, we have numerically investigated rheological and transport of soft under shear. MD of disks, we have the rheology, self-diffusion, and size of with the control parameters,i.e. the externally imposed shear $\dot{\gamma}$ and fraction of the disks $\phi$. Our result is the critical scaling of the diffusivity \[sub:diff\]) and size of rigid clusters (Sec. \[sub:rigid\]), where their dependence on both $\dot{\gamma}$ and reported \[Eqs. (\[eq:D1\]), (\[eq:D2\]), and (\[eq:xi2\])\]. The has calculated both of jamming a single numerical protocol) to unify the understanding of self-diffusion of particulate systems: We found that (i) the diffusivity below jamming crossover the linear scaling to the power-law $D\sim\dot{\gamma}^{3/4}$. a can be also seen simulations
eq:D1\])\], is caused by the diverging Length scalE beloW jaMmiNg, I.e. $\xi\Sim|\DElta\phi|^{-1}$ \[Eq. (\[eq:xi1\])\]. AS ShowN in Fig. \[fig:xi\](b), the linear rElatiOn (DOtteD LiNe) welL explaiNS oUR ResUlTs If tHe LEnGth scAle $\Xi$ is smaLler than $10d_0$. IF thE sYstem is above JAmMing, the lenGth Scale increasEs (mOre thaN $10d_0$) WitH The inCreAse of $\Phi$. HowEVer, the DiffusiviTy OVer the SHear ratE $d/\DoT{\gamMa}$ starts to deviate FRoM The linear relatIon (dotTeD LiNE) And The Length scalE rEacheS The systEM sIZE $l/2\siMEq 44d_0$ (shaded regiOn). We concludE ThaT this dEvIatIOn is caUsed bY fINitE-size effectS and Further stUdies oF DiffereNT system Sizes aRe nEceSsarY (As In refS. [@dIFf_sHEaR_md3; @DIff_Shear_md4]) tO fIgUre ouT the RELATion BetWeen $d/\dot{\gAmma}$ and $\xi$ in thIs rEgimE, WhiCh we pOstpoNe as A fUture Work.

DiScussIoNs {#sec:disc}
===========

In this StudY, we have nuMerIcAllY iNvestIGated rHeoLogIcal and TranspoRT prOpERTIeS of soft athermal parTiCLEs Under sheAr. EmplOYiNg md simulatIoNs oF two-DIMensiOnal DIsKs, we have ClarifIEd HoW the rheOlOgy, selF-dIffUsiOn, and SIze oF rigid Clusters Vary wITh the control paRAmeters,i.e. the eXTeRNAlLY impOseD shear rate $\dOt{\gaMMa}$ anD pacKInG frACtion Of the DiSKs $\PHi$. Our main result is thE cRiticaL scalIng of the diffuSivity (Sec. \[sUB:DIff\]) and siZe of RIgID clusters (Sec. \[suB:rigiD\]), where theiR DependenCe on bOth $\dot{\gaMma}$ and $\phi$ IS Reported \[eqs. (\[Eq:D1\]), (\[Eq:D2\]), (\[Eq:xI1\]), ANd (\[Eq:xi2\])\]. The diffusIVIty hAs Been calCulAted on bOth SidEs oF jaMmIng (by a sinGle numerIcAl PrOtOcoL) to unIFy the undErStaNdIng Of selF-DiffusIon of Soft PaRtICulAte systEMs: wE FounD tHaT (i) thE diFfUsiviTy beLOw jAmming eXhibits a cRosSOver FrOm The lineAr scaling $D\sim\DoT{\gamma}$ to thE pOweR-law $D\sIM\Dot{\gamma}^{3/4}$. such a crossover can be also SEen in prEviOus siMulaTions [@diff_She
eq:D1\])\], is caused by t he divergi ng le ngt h s ca le b elow  jamming, i.e. $\xi \sim|\Delta\phi|^{-1}$  \[Eq .( \[eq : xi 1\])\ ]. As s h ow n inFi g.  \[ fi g :x i\](b ),the lin ear relati on(d otted line)w el l explains  ou r results if  th e leng th  sc a le $\ xi$  is s maller than $ 10d_0$. I ft he sys t em is a b o ve  jam ming, the lengths ca l e increases (m ore th an $1 0 d _0$ ) w ith the in cr easeo f $\phi $ .H o w eve r , the diffusi vity over t h e s hear r at e $ D /\dot{ \gamm a} $  st arts to dev iate  from the  linea r  relati o n (dott ed lin e)and  the le ng thsc a ler ea che s  th e system  s iz e $L/ 2\si m e q 44d_ 0$(sha ded r egion). We co ncl udet hat  this  devi atio nis ca used b y fin it e-size effectsandfurther s tud ie s o fdiffe r ent sy ste m s izes ar e neces s ary  ( a s in  Refs. [@diff_shea r_ m d 3;  @diff_s hear_m d 4] )t o figure  o utther e latio n be t we en $D/\d ot{\ga m ma }$  and $\ xi $ in t hi s r egi me, w h ichwe pos tpone as  a fu t ure work.

Dis c ussions {#sec : di s c }= ==== === ===

In thi s st u dy,we h a ve  nu m erica lly i nv e st i gated rheological a nd  trans portproperties of  soft athe r m a l partic lesu nd e r shear. Emplo yingMD simulat i ons of t wo-di mensiona l disks,w e  have cl ari fie d h owt h erheology, sel f - diff us ion, an d s ize ofrig idclu ste rs  vary wit h the co nt ro lpa ram eters , i.e. the  e xte rn all y imp o sed sh ear r ate$\ do t {\g amma}$a nd p acki ng  f ract ion  o f the  dis k s $ \phi$.Our mainres u lt i sth e criti cal scaling o fthe diffus iv ity  (Sec. \ [sub:dif f\]) and size of rigidc lusters  (S ec. \ [sub :rigid\]) , w here t hei r  depen denceon bo th  $\ d o t{\ga m m a} $ a nd  $\phi$ is r epo rted\[ Eqs.  (\[eq: D1\]), (\[eq:D2\]) ,  (\ [eq:xi1\]), a nd(\[e q : xi 2\] ) \] .  Th ed iff u s ivity has beencalculated  o n  b oth sideso f j am ming (b y a sin gle n u merical  protocol ) to unif ytheu n der standing o f self-d iffusiono f sof t  p artic ula te sys te ms:  We f ound t h at(i) t he dif fu sivity  belo wjammingexhibits a crossover fr om the  line arscaling $ D\s i m\d ot{\gamma }$ t o the powe r-l aw$D\si m\d o t{\ga mma} ^ {3 /4} $ . Suc h ac rossoverc an  be a ls o seen in p r e v iou s sim ula t ions [ @dif f_she
eq:D1\])\], is_caused by_the diverging length scale_below jamming, i.e. $\xi\sim|\Delta\phi|^{-1}$_\[Eq. (\[eq:xi1\])\]._As shown_in_Fig. \[fig:xi\](b), the linear_relation (dotted line)_well explains our results_if the length_scale_$\xi$ is smaller than $10d_0$. If the system is above jamming, the length scale_increases_(more than_$10d_0$)_with_the increase of $\phi$. However,_the diffusivity over the shear_rate $D/\dot{\gamma}$_starts to deviate from the linear relation (dotted_line)_and the length_scale reaches the system size $L/2\simeq 44d_0$ (shaded region)._We conclude that this deviation is_caused by finite-size_effects_and_further studies of different_system sizes are necessary (as in_Refs. [@diff_shear_md3; @diff_shear_md4]) to figure out the_relation between $D/\dot{\gamma}$ and $\xi$ in this_regime, which we postpone as a_future work.

Discussions {#sec:disc}
===========

In this study,_we have_numerically investigated rheological and transport_properties of soft_athermal particles_under shear. Employing_MD simulations of two-dimensional disks, we_have clarified how_the rheology, self-diffusion, and size of_rigid_clusters vary with_the_control_parameters,i.e. the externally_imposed shear rate_$\dot{\gamma}$_and packing_fraction_of the disks $\phi$. Our main_result_is the critical scaling of the diffusivity_(Sec. \[sub:diff\]) and size of_rigid_clusters (Sec. \[sub:rigid\]), where their_dependence on both $\dot{\gamma}$ and_$\phi$ is reported \[Eqs. (\[eq:D1\]), (\[eq:D2\]), (\[eq:xi1\]),_and (\[eq:xi2\])\]._The diffusivity_has been calculated on both sides of jamming (by a single_numerical protocol) to unify the understanding_of self-diffusion of soft_particulate systems:_We_found that (i)_the_diffusivity below_jamming exhibits a crossover from the linear_scaling $D\sim\dot{\gamma}$_to the power-law $D\sim\dot{\gamma}^{3/4}$. Such a_crossover can be also_seen_in previous simulations [@diff_she