1)} Z(\phi)\right) \rho(\phi) \right ],$$ for which the continuity of the flux in leads to the following non-linear boundary condition $$\label{eq:denseLimitBoundary}
\rho(0) ( 1 + \lambda Z(0)) = \rho(1) (1 + \lambda Z(1)).$$ The dynamics of the system defined by and is easily describable for monotonous PRCs [@Mauroy2013]. For increasing $Z(\phi)$ the probability density concentrates to a single phase in finite time for arbitrary initial distributions. For decreasing $Z(\phi)$ convergence to the stationary solution $\rho_0(\phi):= c/ ( 1 + c Z(\phi))$ can be observed.

The question of synchronization in the population model with sparse connectivity can be addressed by investigating the existence and the stability of normalized stationary solutions of and. Stationary solutions of Eq.  (with $\partial_t \rho(\phi,t) = 0$) can be obtained by segmenting $[0,1]$ into intervals in which oscillators receive either phase advances or retardations. For each of these intervals, a solution can then be obtained with some solver for delay differential equations with state dependent delays, stepping towards either larger or smaller phases. Phases $\varphi$ at which $R(\phi)$ crosses the identity from below serve as suitable starting points for such a stepping approach because the sets $I_>(\varphi)$ and $I_\leq(\varphi)$ are empty at such points and we have $\rho(\varphi) = J(\varphi) = J(1) = \lambda /m$. Using this initial value, solutions fulfill. In this way we can obtain stationary solutions $\rho(\phi;\lambda)$ of and in sole dependence on $\lambda$. Let us denote $I(\lambda) := \int_0^1 \rho(\phi;\lambda) d \phi$. In order to allow for a stochastic interpretation of $\rho(\phi;\lambda)$, $\lambda$ must then be chosen with shooting in such a way that $I(\lambda) = 1$. However, depending on the PRC such a choice may not be possible. We can characterize the condition under which solutions exist by assuming that $R(\phi)$ is non-decreasing. Note that this assumption is valid for commonly considered PRCs including those of integrate-and-fire oscillators [@Brown2004a]. Under this assumption $I(\lambda)$ is strictly increasing in $\lambda
1) } Z(\phi)\right) \rho(\phi) \right ], $ $ for which the continuity of the flux in leads to the following non - linear boundary circumstance $ $ \label{eq: denseLimitBoundary } 
 \rho(0) (1 + \lambda Z(0) ) = \rho(1) (1 + \lambda Z(1)).$$ The moral force of the organization defined by and is well describable for monotonous china [ @Mauroy2013 ]. For increasing $ Z(\phi)$ the probability density digest to a unmarried phase in finite time for arbitrary initial distributions. For decreasing $ Z(\phi)$ convergence to the stationary solution $ \rho_0(\phi):= c/ (1 + vitamin c Z(\phi))$ can be observed. 

 The question of synchronization in the population model with sparse connectivity can be addressed by investigating the being and the stability of normalized stationary solutions of and. Stationary solution of Eq.    (with $ \partial_t \rho(\phi, t) = 0 $) can be obtained by segmenting $ [ 0,1]$ into intervals in which oscillators pick up either phase advances or retardations. For each of these intervals, a solution can then be obtained with some solver for delay differential equation with state dependent delays, stepping towards either larger or smaller phases. Phases $ \varphi$ at which $ R(\phi)$ crosses the identity from below serve as suitable starting points for such a stepping approach because the sets $ I_>(\varphi)$ and $ I_\leq(\varphi)$ are empty at such points and we have $ \rho(\varphi) = J(\varphi) = J(1) = \lambda /m$. Using this initial value, solution fulfill. In this way we can obtain stationary solution $ \rho(\phi;\lambda)$ of and in lone dependence on $ \lambda$. Let us announce $ I(\lambda): = \int_0 ^ 1 \rho(\phi;\lambda) d \phi$. In order to let for a stochastic interpretation of $ \rho(\phi;\lambda)$, $ \lambda$ must then be chosen with shoot in such a way that $ I(\lambda) = 1$. However, depending on the PRC such a choice may not be possible. We can qualify the condition under which solutions exist by assuming that $ R(\phi)$ is non - decreasing. Note that this assumption is valid for commonly considered china include those of integrate - and - fire oscillators [ @Brown2004a ]. Under this assumption $ I(\lambda)$ is rigorously increasing in $ \lambda
1)} Z(\pji)\right) \rho(\phi) \right ],$$ fov which the contnbuity mf the flux in leads to the following non-lmneae boubdary condition $$\label{ed:denseLimptBoundart}
\rho(0) ( 1 + \lambda V(0)) = \rho(1) (1 + \lambdz Z(1)).$$ Tke dynamics of tme system dafined by and hs ecsily describable for monotonous PRCf [@Maurou2013]. Vor increasing $Z(\php)$ ehe lgonability density concentrates to a singne phase in fonite time for arbitrary ijitiwl distributions. Flr decreasibg $Z(\[yi)$ convergenze to the stationary sklution $\rho_0(\phi):= c/ ( 1 + c Z(\phi))$ can ce obxerved.

The wuwstlmn of synchconizanion in the populatiot model with sparse cpnnxcticity can be addressed by investigating jhe existetcz and the stability od bormanizeg stxriovarg xomutiond oh and. Statiknary solutuons of Eq.  (with $\paruiaj_n \rho(\phi,t) = 0$) czn be jbeained by segmenting $[0,1]$ into intervals in whjch oscillators receive either phase advanced or retawdations. For each of these intervals, a solution cdn thxn be ontaivwd with some solver for delay differential equaejoms with state deiendent delays, stelplnb towards eithgr largzd kr smaller phases. Ohases $\darphu$ at whicr $R(\pni)$ crosses the identity fron below servv as suitable starting points for sucn a syepping approach becausz the aets $I_>(\varphl)$ and $I_\les(\xarphi)$ are empty at sgch poinur and we have $\rho(\darphi) = J(\tarphn) = J(1) = \laobda /m$. Usigg this inltial value, solutions fupfill. It this way we can obtain stationary solutmpns $\rho(\phi;\lakbga)$ mf and iu sole dependence og $\lambda$. Let ux denotz $I(\lamcda) := \int_0^1 \rgo(\phi;\lakbda) d \phi$. Yn order to anpow for a stmchastic intwrprwtation uf $\rho(\phi;\lambds)$, $\lambda$ must then ve chosen with shoptivf in such a way thqt $I(\lambda) = 1$. Hoeevdr, qeienving jt the PRC suwh a chuoce mxy not be pirsibke. We can characterise tge condition under wmich solujions exift by assuminb that $R(\phi)$ is non-fecreesing. Iote tnat this assumption is valid for dommonly fonfidered PRCs incouding those of integrate-and-fire oscillators [@Brown2004a]. Under this assumktion $I(\lambda)$ is strietky increasiig in $\jambda
1)} Z(\phi)\right) \rho(\phi) \right ],$$ for which of flux in to the following ( + \lambda Z(0)) \rho(1) (1 + Z(1)).$$ The dynamics of the system by and is easily describable for monotonous PRCs [@Mauroy2013]. For increasing $Z(\phi)$ the density concentrates to a single phase in finite time for arbitrary initial distributions. decreasing convergence the solution $\rho_0(\phi):= c/ ( 1 + c Z(\phi))$ can be observed. The question of synchronization in population model with sparse connectivity can be addressed investigating the existence and stability of normalized stationary solutions and. solutions of (with \rho(\phi,t) 0$) can be by segmenting $[0,1]$ into intervals in which oscillators receive either phase advances or retardations. For each of intervals, a then be with solver delay differential equations dependent delays, stepping towards either larger Phases $\varphi$ at which $R(\phi)$ crosses the identity below serve suitable starting points for such a approach because the sets $I_>(\varphi)$ and $I_\leq(\varphi)$ are at such points and we have $\rho(\varphi) = J(\varphi) = J(1) = \lambda /m$. Using value, solutions fulfill. In way we can stationary $\rho(\phi;\lambda)$ and sole dependence $\lambda$. Let us denote $I(\lambda) := \int_0^1 \rho(\phi;\lambda) d \phi$. In to allow for a stochastic interpretation of $\rho(\phi;\lambda)$, $\lambda$ must chosen shooting in such way that $I(\lambda) = However, on the PRC such may be characterize condition which solutions exist by that $R(\phi)$ is non-decreasing. Note this assumption is valid those of integrate-and-fire oscillators [@Brown2004a]. Under this assumption is strictly increasing in $\lambda
1)} Z(\phi)\right) \rho(\phi) \right ],$$ for whIch the contInuitY of The FlUx in LeadS to the followinG Non-lInear boundary condition $$\Label{Eq:DEnselImItBouNdary}
\rhO(0) ( 1 + \LaMBDa Z(0)) = \RhO(1) (1 + \lAmbDa z(1)).$$ thE dynaMicS of the sYstem definEd bY aNd is easily deSCrIbable for mOnoTonous PRCs [@MaUroY2013]. For inCrEasINg $Z(\phI)$ thE probAbilitY DensitY concentrAtES to a siNGle phasE IN fInitE time for arbitrary INiTIal distributioNs. For dEcREaSINg $Z(\Phi)$ ConvergencE tO the sTAtionarY SoLUTIon $\RHo_0(\phi):= c/ ( 1 + c Z(\phi))$ caN be observed.

tHe qUestioN oF syNChroniZatioN iN The Population mOdel With sparsE conneCTivity cAN be addrEssed bY inVesTigaTInG tHe eXiSTenCE aNd tHE stAbility oF nOrMalizEd stATIONary SolUtioNs of aNd. Stationary sOluTionS Of EQ.  (with $\PartiAl_t \rHo(\Phi,t) = 0$) cAn be obTaineD bY segmenting $[0,1]$ into InteRvals in whIch OsCilLaTors rECeive eIthEr pHase advAnces or REtaRdATIOnS. For each of these intErVALs, A solutioN can thEN bE oBTained wiTh SomE solVER for dElay DIfFerentiaL equatIOnS wIth statE dEpendeNt DelAys, SteppINg toWards eIther larGer or SMaller phases. PhASes $\varphi$ at whICh $r(\PHi)$ CRossEs tHe identity fRom bELow sErve AS sUitABle stArtinG pOInTS for such a stepping apPrOach beCause The sets $I_>(\varphI)$ and $I_\leq(\vaRPHI)$ are emptY at sUCh POints and we have $\Rho(\vaRphi) = J(\varphI) = j(1) = \lambda /m$. using This initIal value, sOLUtions fuLfiLl. IN thIs wAY We Can obtain statIONary SoLutions $\Rho(\Phi;\lambDa)$ oF anD in SolE dEpendence On $\lambda$. leT uS dEnOte $i(\lambDA) := \int_0^1 \rho(\pHi;\LamBdA) d \pHi$. In oRDer to aLlow fOr a sToChAStiC interpREtATIon oF $\rHo(\Phi;\lAmbDa)$, $\LambdA$ musT TheN be chosEn with shoOtiNG in sUcH a Way that $i(\lambda) = 1$. HoweveR, dEpending on ThE PRc such a CHOice may nOt be possible. We can characTErize thE coNditiOn unDer which sOluTions eXisT By assuMing thAt $R(\phI)$ iS noN-DEcreaSINg. notE tHat this assUMPtiOn is vAlId foR commonLy considered PRCs inCLudIng those of intEgrAte-aND-FiRe oSCiLLatOrS [@broWN2004A]. Under this assumPtion $I(\lambDa)$ IS sTrictly incREasInG in $\lambDa
1)} Z(\phi)\right) \rho(\p hi) \right  ],$$  fo r w hi ch t he c ontinuity of t h e fl ux in leads to the fol lowin gn on-l i ne ar bo undaryc on d i tio n$$ \la be l {e q:den seL imitBou ndary}
\rh o(0 )( 1 + \lambd a  Z (0)) = \rh o(1 ) (1 + \lamb daZ(1)). $$  Th e  dyna mic s ofthe sy s tem de fined byan d  is ea s ily des c r ib able  for monotonous P R Cs [@Mauroy2013].  For i nc r ea s i ng$Z( \phi)$ the  p robab i lity de n si t y con c entrates to a  single pha s e i n fini te  ti m e forarbit ra r y i nitial dist ribu tions. Fo r decr e asing $ Z (\phi)$  conve rge nce  tot he  s tat io n ary so lut i on$\rho_0( \p hi ):= c / (1 + c Z( \ph i))$  canbe observed.
Th e qu e sti on of  sync hron iz ation  in th e pop ul ation model wit h sp arse conn ect iv ity  c an be addres sed  by  invest igating the  e x i s te nce and the stabil it y of  normali zed st a ti on a ry solut io nsof a n d . Sta tion a ry  solutio ns ofE q.  (with $ \p artial _t  \r ho( \phi, t ) =0$) ca n be obt ained by segmenting$ [0,1]$ into i n te r v al s  inwhi ch oscillat orsr ecei ve e i th erp haseadvan ce s  o r  retardations. Forea ch ofthese  intervals, a  solutionc a n  then be  obt a in e d with some so lverfor delayd ifferent ial e quations  with sta t e  depende ntdel ays , s t e pp ing towards e i t herla rger or  sm aller p has es.  Ph ase s$\varphi$  at whic h$R (\ ph i)$  cros s es the i de nti ty  fr om be l ow ser ve as  sui ta bl e  st artingp oi n t s fo rsu ch a  st ep pingappr o ach  becaus e the set s $ I _>(\ va rp hi)$ an d $I_\leq(\va rp hi)$ are e mp tyat suc h points a nd we have $\rho(\varph i ) = J(\ var phi)= J( 1) = \lam bda  /m$.Usi n g this  initi al va lu e,s o lutio n s  f ulf il l. In this w aywe ca nobta in stat ionary solutions $ \ rho (\phi;\lambda )$of a n d  i n s o le dep en d enc e on $\lambda$. L et us deno te $I (\lambda): = \ in t_0^1 \ rho(\ph i;\la m bda) d\phi$. In  order to  a llow f ora stochast ic inter pretation of $\ r ho (\phi ;\l ambda) $,  $\ lambd a$ mus t  th en be  chose nwith s hooti ng  in such  a way that $I(\lambda)  = 1$.  Howe ver , dependi ngo n t he PRC su ch a  choice ma y n otbe po ssi b le. W e ca n  c har a cteri ze t h e conditi o nund e r  w hich soluti o n s  ex ist b y a s suming  tha t $R(\phi)$ is no n -decreasing. N otet h atthi s  ass um ption is valid  fo rc o mmonly c on sidered PRC s includ in g  thos e of i ntegra te-and- f i re oscill ator s [ @Brown200 4a] .U nder th is  a s sumpti on $ I( \lambd a)$ is stri c t ly increasing in  $\la m b da
1)} Z(\phi)\right)_\rho(\phi) \right_],$$ for which the_continuity of_the_flux in_leads_to the following_non-linear boundary condition_$$\label{eq:denseLimitBoundary}
\rho(0) ( 1 +_\lambda Z(0)) =_\rho(1)_(1 + \lambda Z(1)).$$ The dynamics of the system defined by and is easily_describable_for monotonous_PRCs_[@Mauroy2013]._For increasing $Z(\phi)$ the probability_density concentrates to a single_phase in_finite time for arbitrary initial distributions. For decreasing_$Z(\phi)$_convergence to the_stationary solution $\rho_0(\phi):= c/ ( 1 + c Z(\phi))$_can be observed.

The question of synchronization_in the population_model_with_sparse connectivity can be_addressed by investigating the existence and_the stability of normalized stationary solutions_of and. Stationary solutions of Eq.  (with_$\partial_t \rho(\phi,t) = 0$) can be_obtained by segmenting $[0,1]$ into_intervals in_which oscillators receive either phase_advances or retardations._For each_of these intervals,_a solution can then be obtained_with some solver_for delay differential equations with state_dependent_delays, stepping towards_either_larger_or smaller_phases. Phases $\varphi$_at_which $R(\phi)$_crosses_the identity from below serve as_suitable_starting points for such a stepping approach_because the sets $I_>(\varphi)$_and_$I_\leq(\varphi)$ are empty at_such points and we have_$\rho(\varphi) = J(\varphi) = J(1) =_\lambda /m$._Using this_initial value, solutions fulfill. In this way we can obtain stationary_solutions $\rho(\phi;\lambda)$ of and in sole_dependence on $\lambda$. Let_us denote_$I(\lambda)_:= \int_0^1 \rho(\phi;\lambda)_d_\phi$. In_order to allow for a stochastic interpretation_of $\rho(\phi;\lambda)$,_$\lambda$ must then be chosen with_shooting in such a_way_that $I(\lambda) = 1$. However, depending_on the PRC such a choice_may not be possible. We_can_characterize_the condition under which solutions_exist by assuming that $R(\phi)$ is_non-decreasing. Note that_this assumption is valid for commonly considered_PRCs_including those of integrate-and-fire oscillators [@Brown2004a]._Under_this assumption $I(\lambda)$ is strictly increasing_in_$\lambda