k-1} +1$ for $k>1$. Note that the zeroth moment of the force distribution is the participation ratio $\psi$ shown in Fig. \[fig:3\]b. The mass fraction of the tensional backbone that appears at the critical strain is given by the participation ratio or zeroth moment at $\gamma_c$ [@hansen_universality_1989; @bunde_fractals_1995]. In plotting the mass of the tensional structure at the critical strain versus system size $W$, we find that the fractal dimension of this backbone appears to be the same as the euclidean dimension of $2$ (see Fig. \[fig:A8\] in the Appendix).

Of particular interest are the macroscopic properties of fiber networks such as stiffness $K$ near the transition. As we approach the critical point, we find that $K$ shows a finite discontinuity $K_c$, in agreement with prior work [@vermeulen_geometry_2017; @merkel_minimal-length_2019]. Figure \[fig:2\] shows the behavior of one random realization of a diluted triangular network very close to its critical strain $\gamma - \gamma_c \simeq 10^{-4}$. In order to find the sample-specific critical point $\gamma_c(W,i)$ for a network with size $W$, we use the bisection method [@merkel_minimal-length_2019]. By performing an initial step-wise shearing simulation for every random sample, we first find a strain value $\gamma_{R,i}$ at which the network becomes rigid, i.e., the shear stress calculated from Eq. \[eq:3\] reaches a threshold value. Here we use $10^{-9}$ for the stress threshold. Our results, however, are insensitive to the choice of the threshold value as long as we use a sufficiently small value. The prior strain value to $\gamma_{R,i}$ is considered as the nearest floppy point $\gamma_{F,i}$. Modifying the bracket \[$\gamma_{F,i}$, $\gamma_{R,i}$\] in at least $20$ bisection steps, we are able to accurately identify the critical point for every random sample $i$. After identifying the critical point, the network is sheared in a step-wise manner from $\gamma_c(W,i)$. Therefore, the final ensemble averages of a specific system size are taken
k-1 } +1 $ for $ k>1$. Note that the zeroth moment of the force distribution is the participation proportion $ \psi$ usher in Fig.   \[fig:3\]b. The mass fraction of the tensional backbone that appear at the critical tune is given by the participation proportion or zeroth moment at $ \gamma_c$ [ @hansen_universality_1989; @bunde_fractals_1995 ]. In plotting the mass of the tensional social organization at the critical breed versus system size $ W$, we witness that the fractal dimension of this backbone appears to be the same as the euclidean dimension of $ 2 $ (attend Fig.   \[fig: A8\ ] in the Appendix). 

 Of particular interest are the macroscopic properties of character networks such as stiffness $ K$ near the transition. As we approach the critical point, we detect that $ K$ shows a finite discontinuity $ K_c$, in agreement with prior work [ @vermeulen_geometry_2017; @merkel_minimal - length_2019 ]. Figure   \[fig:2\ ] shows the demeanor of one random realization of a diluted triangular network very close to its critical strain $ \gamma - \gamma_c \simeq 10^{-4}$. In order to find the sample - specific critical distributor point $ \gamma_c(W, i)$ for a network with size $ W$, we use the bisection method [ @merkel_minimal - length_2019 ]. By performing an initial step - wise shearing simulation for every random sample, we first find a stress value $ \gamma_{R, i}$ at which the network becomes fixed, i.e., the shear tension account from Eq.   \[eq:3\ ] reaches a threshold value. Here we use $ 10^{-9}$ for the stress threshold. Our results, however, are insensitive to the choice of the brink value as long as we use a sufficiently small value. The prior strain value to $ \gamma_{R, i}$ is considered as the approximate floppy point $ \gamma_{F, i}$. Modifying the bracket \[$\gamma_{F, i}$, $ \gamma_{R, i}$\ ] in at least $ 20 $ bisection steps, we are able to accurately identify the critical distributor point for every random sample $ i$. After identify the critical point, the network is sheared in a step - wise manner from $ \gamma_c(W, i)$. consequently, the final ensemble averages of a specific system size are taken
k-1} +1$ vor $k>1$. Note that the zerouh moment of the yirce dmstribufion is ghe participation ratio $\psi$ dhiwn ib Fig. \[fig:3\]b. The mass fraztion of nhe tensiinal vackbone tizt appears at bhe cxivical strain is given by dhe participathov xatio or zeroth moment at $\gamma_c$ [@hanfen_univrrdality_1989; @bunde_frwctakf_1995]. In ilitting the mass of the tensiohal strlcture at the criyical strain versus system sizf $W$, we find that tje fractal eimegwion of this backbone appears to bg the same as the euclidean dimevsion of $2$ (see Fug. \[dig:W8\] in the Appxndix).

Os particular interesd are tne macroscopic prmpeeties of fiber networns such as stiffnefs $K$ near tke transition. As we appeoach the crigucau pkiit, se finf tiat $K$ shows a finite duscontinuity $K_c$, in sgwvrment with pdior wjrh [@vermeulen_geometry_2017; @merkel_minimal-length_2019]. Fifure \[fig:2\] shows the behavuor of one random reapization jf a diluted triangular network very close to its critmcxl wtvain $\talma - \gamma_c \simeq 10^{-4}$. In order to find the samplq-aptcinic critical poikt $\gamma_c(W,i)$ for a mehwptk with size $W$, we usz tge bisection methof [@merkej_mininal-length_2019]. By lerforming an initial step-wuse shearing wimulation for evexy random saople, we forst find a strain valuz $\gammz_{R,i}$ at whicj the netsurk becomes rigia, i.v., tha shear stress calculated srom Eq. \[eq:3\] reaehes a tfresnold vwlue. Here ae usc $10^{-9}$ for the stress tjreshllg. Our resupts, however, are insensitive to vie choice of jhe thgeshold vclue ax long as we tse a sufficieutly smajl vauue. The prpor straii value to $\gwmma_{R,i}$ is convldered as thx nearest floppy point $\gxoma_{F,i}$. Modifyinb the brabktt \[$\gamma_{F,i}$, $\gamma_{R,i}$\] in at leaxt $20$ gisection steps, we are able to acvurxtejy iventisf the criticdl puing for dvery randon samlle $i$. After identifyhng fhe critical point, tme networj is shewred in a stel-wise manner from $\hamma_r(W,i)$. Thxreforr, tre final ensemble averages of z specifif snstem size arq tajen
k-1} +1$ for $k>1$. Note that the of force distribution the participation ratio The fraction of the backbone that appears the critical strain is given by participation ratio or zeroth moment at $\gamma_c$ [@hansen_universality_1989; @bunde_fractals_1995]. In plotting the mass the tensional structure at the critical strain versus system size $W$, we find the dimension this appears to be the same as the euclidean dimension of $2$ (see Fig. \[fig:A8\] in the Of particular interest are the macroscopic properties of networks such as stiffness near the transition. As we the point, we that shows finite discontinuity $K_c$, agreement with prior work [@vermeulen_geometry_2017; @merkel_minimal-length_2019]. Figure \[fig:2\] shows the behavior of one random realization of a triangular network to its strain - \simeq 10^{-4}$. In find the sample-specific critical point $\gamma_c(W,i)$ with size $W$, we use the bisection method By performing initial step-wise shearing simulation for every sample, we first find a strain value $\gamma_{R,i}$ which the network becomes rigid, i.e., the shear stress calculated from Eq. \[eq:3\] reaches a Here we use $10^{-9}$ the stress threshold. results, are to choice of threshold value as long as we use a sufficiently small value. prior strain value to $\gamma_{R,i}$ is considered as the nearest $\gamma_{F,i}$. the bracket \[$\gamma_{F,i}$, in at least $20$ steps, are able to accurately critical for $i$. identifying critical point, the network sheared in a step-wise manner $\gamma_c(W,i)$. Therefore, the final system size are taken
k-1} +1$ for $k>1$. Note that the zeroth momEnt of the foRce diStrIbuTiOn is The pArticipation raTIo $\psI$ shown in Fig. \[fig:3\]b. The mass FractIoN Of thE TeNsionAl backbONe THAt aPpEaRs aT tHE cRiticAl sTrain is Given by the ParTiCipation ratiO Or Zeroth momeNt aT $\gamma_c$ [@hanseN_unIversaLiTy_1989; @bUNde_frActAls_1995]. In PlottiNG the maSs of the teNsIOnal stRUcture aT THe CritIcal strain versus sYStEM size $W$, we find thAt the fRaCTaL DImeNsiOn of this baCkBone aPPears to BE tHE SAme AS the euclidean Dimension of $2$ (SEe FIg. \[fig:A8\] In The aPpendiX).

Of paRtICulAr interest aRe thE macroscoPic proPErties oF Fiber neTworks SucH as StifFNeSs $k$ neAr THe tRAnSitIOn. AS we approAcH tHe criTicaL POINt, we FinD thaT $K$ shoWs a finite discOntInuiTY $K_c$, In agrEemenT witH pRior wOrk [@verMeuleN_gEometry_2017; @merkel_miNimaL-length_2019]. FiGurE \[fIg:2\] sHoWs the BEhavioR of One Random rEalizatIOn oF a DILUtEd triangular networK vERY cLose to itS critiCAl StRAin $\gamma - \GaMma_C \simEQ 10^{-4}$. in ordEr to FInD the sampLe-specIFiC cRitical PoInt $\gamMa_C(W,i)$ For A netwORk wiTh size $w$, we use thE biseCTion method [@merkEL_minimal-lengtH_2019]. by PERfORminG an Initial step-Wise SHearIng sIMuLatIOn for Every RaNDoM Sample, we first find a sTrAin valUe $\gamMa_{R,i}$ at which thE network beCOMEs rigid, i.E., the SHeAR stress calculaTed frOm Eq. \[eq:3\] reacHEs a thresHold vAlue. Here We use $10^{-9}$ for tHE Stress thResHolD. OuR reSULtS, however, are inSENsitIvE to the cHoiCe of the ThrEshOld ValUe As long as wE use a sufFiCiEnTlY smAll vaLUe. The priOr StrAiN vaLue to $\GAmma_{R,i}$ Is conSideReD aS The Nearest FLoPPY poiNt $\GaMma_{F,I}$. MoDiFying The bRAckEt \[$\gamma_{f,i}$, $\gamma_{R,i}$\] In aT LeasT $20$ bIsEction sTeps, we are able To Accurately IdEntIfy the CRItical poInt for every random sample $I$. after idEntIfyinG the Critical pOinT, the neTwoRK is sheAred in A step-WiSe mANNer frOM $\GaMma_C(W,I)$. Therefore, THE fiNal enSeMble AverageS of a specific system SIze Are taken
k-1} +1$ for $k>1$. Note t hat the ze rothmom ent  o f th e fo rce distributi o n is  the participation rat io $\ ps i $ sh o wn  in F ig. \[f i g: 3 \ ]b.  T he  ma ss fr actio n o f the t ensional b ack bo ne that appe a rs  at the cr iti cal strain i s g iven b ythe parti cip ation  ratio or zer oth momen ta t $\ga m ma_c$ [ @ h an sen_ universality_1989 ;  @ b unde_fractals_ 1995].  I n  p l o tti ngthe mass o fthe t e nsional st r u c tur e  at the criti cal strainv ers us sys te m s i ze $W$ , wefi n d t hat the fra ctal  dimensio n of t h is back b one app ears t o b e t he s a me  a s t he euc l id ean dim ension o f$2 $ (se e Fi g . \ [fig :A8 \] i n the  Appendix).
Ofpart i cul ar in teres t ar ethe m acrosc opicpr operties of fib er n etworks s uch  a s s ti ffnes s  $K$ n ear  th e trans ition.A s w ea p p ro ach the critical p oi n t ,we findthat $ K $sh o ws a fin it e d isco n t inuit y $K _ c$ , in agr eement wi th  priorwo rk [@v er meu len _geom e try_ 2017;@merkel_ minim a l-length_2019] .  Figure \[fig : 2\ ] sh o ws t hebehavior of  one rand om r e al iza t ion o f a d il u te d  triangular network  v ery cl ose t o its critica l strain $ \ g a mma - \g amma _ c\ simeq 10^{-4}$ . Inorder to f i nd the s ample -specifi c critica l point $\ gam ma_ c(W ,i) $ fo r a network w i t h si ze  $W$, w e u se thebis ect ion  me th od [@merk el_minim al -l en gt h_2 019]. By perfo rm ing  a n i nitia l  step- wiseshea ri ng sim ulation fo r ever yra ndom  sa mp le, w e fi r stfind astrain va lue $\ga mm a_ {R,i}$at which thene twork beco me s r igid,i . e., theshear stress calculated from Eq . \ [eq:3 \] r eaches athr eshold  va l ue. He re weuse $ 10 ^{- 9 } $ for t he  st re ss thresho l d . O ur re su lts,  howeve r, are insensitive tothe choice of  th e th r e sh old va l ueas lon g as we use a suf ficientlysm a ll  value. Th e  pr io r strai n value  to $ \ gamma_{ R,i}$ isconsidere das t h e  ne arest flop py point  $\gamma_ { F,i}$ .  M odify ing  the b ra cke t \[$ \gamma _ {F, i}$,$\gamm a_ {R,i}$ \] in  a t least$20$ bisection steps, w e areabletoaccuratel y i d ent ify the c riti cal pointfor  ev ery r and o m sam ple$ i$ . A f ter i dent i fying the cr iti c a lpoint, then e t wor k isshe a red in  a s tep-wise manner f r om $\gamma_c(W ,i)$ . The ref o re,th e final ensemb leav e r ages ofaspecific sy stem siz ea re ta ken
k-1} +1$_for $k>1$._Note that the zeroth_moment of_the_force distribution_is_the participation ratio_$\psi$ shown in_Fig. \[fig:3\]b. The mass fraction_of the tensional_backbone_that appears at the critical strain is given by the participation ratio or zeroth_moment_at $\gamma_c$_[@hansen_universality_1989;_@bunde_fractals_1995]._In plotting the mass of_the tensional structure at the_critical strain_versus system size $W$, we find that the_fractal_dimension of this_backbone appears to be the same as the euclidean_dimension of $2$ (see Fig. \[fig:A8\] in_the Appendix).

Of particular_interest_are_the macroscopic properties of_fiber networks such as stiffness $K$_near the transition. As we approach_the critical point, we find that $K$_shows a finite discontinuity $K_c$, in_agreement with prior work [@vermeulen_geometry_2017;_@merkel_minimal-length_2019]. Figure \[fig:2\]_shows the behavior of one_random realization of_a diluted_triangular network very_close to its critical strain $\gamma_- \gamma_c \simeq_10^{-4}$. In order to find the_sample-specific_critical point $\gamma_c(W,i)$_for_a_network with_size $W$, we_use_the bisection_method_[@merkel_minimal-length_2019]. By performing an initial step-wise_shearing_simulation for every random sample, we first_find a strain value_$\gamma_{R,i}$_at which the network_becomes rigid, i.e., the shear_stress calculated from Eq. \[eq:3\] reaches a_threshold value._Here we_use $10^{-9}$ for the stress threshold. Our results, however, are insensitive_to the choice of the threshold_value as long as_we use_a_sufficiently small value._The_prior strain_value to $\gamma_{R,i}$ is considered as the_nearest floppy_point $\gamma_{F,i}$. Modifying the bracket \[$\gamma_{F,i}$,_$\gamma_{R,i}$\] in at least_$20$_bisection steps, we are able to_accurately identify the critical point for_every random sample $i$. After_identifying_the_critical point, the network is_sheared in a step-wise manner from_$\gamma_c(W,i)$. Therefore, the_final ensemble averages of a specific system_size_are taken