0$. In this case, $b=3/2$, a unique invariant density admitting compact non-atomic support in $[x^*,1]$, $x^*=1/2$ appears. From eqs. (\[eq9\_1\_11\]), (\[eq9\_1\_13\]), after elementary manipulations, the invariant density $p_*(x)$ takes the expression $$p_*(x) =\frac{A}{b^2 -v(x)} \, \exp \left [ -2 \lambda 
\int_{x^*}^x \frac{v(\xi)}{v^2(\xi) - b^2} \, d \xi \right ]
\label{eq9_2_3}$$ where $A$ is a normalization constant.

[![Invariant stationary probability density $p_*(x)$ for the Poisson-Kac scheme (\[eq9\_1\_10\]) with $v(x)$ given by eq. (\[eq9\_2\_1\]) at $D_{\rm eff}=1$, $b=3/2$. The “noisy” line represents the result of stochastic simulation, the smooth line represents eq. (\[eq9\_2\_3\]).[]{data-label="Fig25"}](pot_1d.eps "fig:"){height="6cm"}]{}

Figure \[Fig25\] depicts the comparison of the closed-form expression for the invariant density at $b=3/2$, $D_{\rm eff}=1$ and the results of stochastic simulation of eq. (\[eq9\_1\_10\]).

Ergodicity breaking in higher dimensions {#sec4_3}
----------------------------------------

Ergodicity breaking occurs also in higher dimensional GPK models in the presence of attractive and periodic potentials. Consider a two-dimensional GPK process $$d {\bf x}(t) =  {\bf v}({\bf x}(t)) \, dt + {\bf b}_{\chi_N(t)} \, d t
\label{eqs_1}$$ with ${\bf b}_\alpha=b \, (\cos \phi_\alpha, \sin \phi_\alpha)$, $\phi_\alpha= 2 \pi (\alpha-1)/N$, $\alpha=1,\dots,N$, $A_{\alpha,\beta}=1/N$, and $\lambda_\alpha=\lambda$, $\alpha,\beta=1,\dots, N$, in the presence of a deterministic bias ${\bf v}({\bf x})$ stemming from
0$. In this case, $ b=3/2 $, a unique invariant density admit compendious non - atomic support in $ [ x^*,1]$, $ x^*=1/2 $ appear. From equivalent. (\[eq9\_1\_11\ ]), (\[eq9\_1\_13\ ]), after elementary manipulations, the invariant density $ p_*(x)$ take the expression $ $ p_*(x) = \frac{A}{b^2 -v(x) } \, \exp \left [ -2 \lambda 
 \int_{x^*}^x \frac{v(\xi)}{v^2(\xi) - b^2 } \, d \xi \right ] 
 \label{eq9_2_3}$$ where $ A$ is a normalization changeless. 

 [! [ Invariant stationary probability density $ p_*(x)$ for the Poisson - Kac outline (\[eq9\_1\_10\ ]) with $ v(x)$ given by eq. (\[eq9\_2\_1\ ]) at $ D_{\rm eff}=1 $, $ b=3/2$. The “ noisy ” line present the resultant role of stochastic simulation, the smooth line represents eq. (\[eq9\_2\_3\]).[]{data - label="Fig25"}](pot_1d.eps " fig:"){height="6 cm " } ] { } 

 Figure \[Fig25\ ] picture the comparison of the closed - form expression for the changeless density at $ b=3/2 $, $ D_{\rm eff}=1 $ and the results of stochastic model of eq. (\[eq9\_1\_10\ ]). 

 Ergodicity breaking in high dimensions { # sec4_3 } 
 ---------------------------------------- 

 Ergodicity breaking occurs also in higher dimensional GPK model in the presence of attractive and periodic potentials. Consider a two - dimensional GPK process $ $ d { \bf x}(t) =   { \bf v}({\bf x}(t) ) \, dt + { \bf b}_{\chi_N(t) } \, d t 
 \label{eqs_1}$$ with $ { \bf b}_\alpha = b \, (\cos \phi_\alpha, \sin \phi_\alpha)$, $ \phi_\alpha= 2 \pi (\alpha-1)/N$, $ \alpha=1,\dots, N$, $ A_{\alpha,\beta}=1 / N$, and $ \lambda_\alpha=\lambda$, $ \alpha,\beta=1,\dots, N$, in the bearing of a deterministic diagonal $ { \bf v}({\bf x})$ stem from
0$. In this case, $b=3/2$, a unique inyariant density cemittiig compzct non-agomic support in $[x^*,1]$, $x^*=1/2$ appears. Feom ews. (\[eq9\_1\_11\]), (\[eq9\_1\_13\]), after elementxry manipllations, rhe mnvariant densitb $p_*(x)$ takcf ths expxewsion $$p_*(x) =\frac{A}{n^2 -v(x)} \, \exp \laft [ -2 \lambda 
\itt_{b^*}^x \frac{v(\xi)}{v^2(\xi) - b^2} \, d \xi \right ]
\label{eq9_2_3}$$ rhere $A$ id a normalizatyon bogstahn.

[![Ikvariant stationary probability sensity $p_*(x)$ for the Ppisson-Kac scheme (\[eq9\_1\_10\]) with $v(d)$ gigen by eq. (\[eq9\_2\_1\]) at $D_{\rl eff}=1$, $b=3/2$. The “noift” line reprerents the gzsult of stkchastic simulation, the smooth uine xepresents gq. (\[ee9\_2\_3\]).[]{gata-label="Fij25"}](pot_1d.eis "fig:"){height="6cm"}]{}

Figure \[Xig25\] depocts the compavison of the closed-form exprevsion for the invatiant denshtv at $b=3/2$, $D_{\rm eff}=1$ and thw eesuljs of stozyasgic smmumation of eq. (\[eq9\_1\_10\]).

Ergosicity breajing in higher dimemsyins {#sec4_3}
----------------------------------------

Ergodidity bwewking occurs also in higher dimensional GPI models in the presencw of attractive and pgriodic poeentials. Consider a two-dimensional GPK process $$d {\tf x}(t) =  {\by v}({\bf b}(r)) \, dt + {\bf b}_{\chi_N(t)} \, d t
\label{eqs_1}$$ with ${\bf b}_\alpha=b \, (\sks \pmi_\alpha, \sin \phi_\ajpha)$, $\phi_\alpna= 2 \ly (\alpha-1)/N$, $\alphx=1,\dots,N$, $A_{\ampha,\beta}=1/N$, and $\lambfa_\alpha=\jambdq$, $\alpha,\beua=1,\dotx, N$, in the presence of a dererministic yiaw ${\bf v}({\bf x})$ stemminy from
0$. In this case, $b=3/2$, a unique admitting non-atomic support $[x^*,1]$, $x^*=1/2$ appears. elementary the invariant density takes the expression =\frac{A}{b^2 -v(x)} \, \exp \left [ \lambda \int_{x^*}^x \frac{v(\xi)}{v^2(\xi) - b^2} \, d \xi \right ] \label{eq9_2_3}$$ where $A$ a normalization constant. [![Invariant stationary probability density $p_*(x)$ for the Poisson-Kac scheme (\[eq9\_1\_10\]) $v(x)$ by (\[eq9\_2\_1\]) $D_{\rm eff}=1$, $b=3/2$. The “noisy” line represents the result of stochastic simulation, the smooth line represents (\[eq9\_2\_3\]).[]{data-label="Fig25"}](pot_1d.eps "fig:"){height="6cm"}]{} Figure \[Fig25\] depicts the comparison of closed-form expression for the density at $b=3/2$, $D_{\rm eff}=1$ the of stochastic of (\[eq9\_1\_10\]). breaking in higher {#sec4_3} ---------------------------------------- Ergodicity breaking occurs also in higher dimensional GPK models in the presence of attractive and potentials. Consider GPK process {\bf = v}({\bf x}(t)) \, {\bf b}_{\chi_N(t)} \, d t \label{eqs_1}$$ \, (\cos \phi_\alpha, \sin \phi_\alpha)$, $\phi_\alpha= 2 \pi $\alpha=1,\dots,N$, $A_{\alpha,\beta}=1/N$, $\lambda_\alpha=\lambda$, $\alpha,\beta=1,\dots, N$, in the presence a deterministic bias ${\bf v}({\bf x})$ stemming from
0$. In this case, $b=3/2$, a unique invariaNt density aDmittIng ComPaCt noN-atoMic support in $[x^*,1]$, $x^*=1/2$ APpeaRs. From eqs. (\[eq9\_1\_11\]), (\[eq9\_1\_13\]), after elemEntarY mANipuLAtIons, tHe invarIAnT DEnsItY $p_*(X)$ taKeS ThE exprEssIon $$p_*(x) =\frAc{A}{b^2 -v(x)} \, \exp \lEft [ -2 \LaMbda 
\int_{x^*}^x \fraC{V(\xI)}{v^2(\xi) - b^2} \, d \xi \riGht ]
\Label{eq9_2_3}$$ where $a$ is A normaLiZatIOn conStaNt.

[![InvAriant STationAry probabIlITy densITy $p_*(x)$ for THE POissOn-Kac scheme (\[eq9\_1\_10\]) with $V(X)$ gIVen by eq. (\[eq9\_2\_1\]) at $D_{\rm Eff}=1$, $b=3/2$. ThE “nOIsY” LIne RepResents the ReSult oF StochasTIc SIMUlaTIon, the smooth lIne represenTS eq. (\[Eq9\_2\_3\]).[]{data-LaBel="fIg25"}](pot_1d.Eps "fiG:"){hEIghT="6cm"}]{}

Figure \[FiG25\] depIcts the coMparisON of the cLOsed-forM expreSsiOn fOr thE InVaRiaNt DEnsITy At $b=3/2$, $d_{\Rm eFf}=1$ and the ReSuLts of StocHASTIc siMulAtioN of eq. (\[Eq9\_1\_10\]).

Ergodicity bReaKing IN hiGher dImensIons {#SeC4_3}
----------------------------------------

ErgoDicity BreakInG occurs also in hiGher DimensionAl GpK ModElS in thE PresenCe oF atTractivE and perIOdiC pOTENtIals. Consider a two-diMeNSIoNal GPK prOcess $$d {\BF x}(T) =  {\bF V}({\bf x}(t)) \, dt + {\bF b}_{\Chi_n(t)} \, d t
\LABel{eqS_1}$$ witH ${\Bf B}_\alpha=b \, (\cOs \phi_\aLPhA, \sIn \phi_\alPhA)$, $\phi_\alPhA= 2 \pi (\AlpHa-1)/N$, $\alPHa=1,\doTs,N$, $A_{\alPha,\beta}=1/N$, And $\laMBda_\alpha=\lambda$, $\ALpha,\beta=1,\dots, N$, IN tHE PrESencE of A determinisTic bIAs ${\bf V}({\bf x})$ STeMmiNG from
0$. In this case, $b=3/2$,  a uniqueinvar ian t d en sity  adm itting compact non- atomic support in $[x^ *,1]$ ,$ x^*= 1 /2 $ app ears. F r om e qs.  ( \[ eq9 \_ 1 \_ 11\]) , ( \[eq9\_ 1\_13\]),aft er  elementarym an ipulations , t he invariant  de nsity$p _*( x )$ ta kes  theexpres s ion $$ p_*(x) =\ fr a c{A}{b ^ 2 -v(x) } \,  \ex p \left [ -2 \lam b da 
\int_{x^*}^x\frac{ v( \ xi ) } {v^ 2(\ xi) - b^2}  \ , d \ x i \righ t  ] \ l abe l {eq9_2_3}$$ w here $A$ is a n ormali za tio n  const ant.
[ ! [In variant sta tion ary proba bility density $p_*(x) $ forthe  Po isso n -K ac  sc he m e ( \ [e q9\ _ 1\_ 10\]) wi th  $ v(x)$  giv e n b y eq . ( \[eq 9\_2\ _1\]) at $D_{ \rm  eff } =1$ , $b= 3/2$.  The  “ noisy ” line  repr es ents the result  ofstochasti c s im ula ti on, t h e smoo thlin e repre sents e q . ( \[ e q 9 \_ 2\_3\]).[]{data-la be l = "F ig25"}]( pot_1d . ep s" fig:"){h ei ght ="6c m " }]{}
Fig u re  \[Fig25 \] dep i ct sthe com pa risonof  th e c losed - form  expre ssion fo r the invariant dens i ty at $b=3/2$ ,  $ D _ {\ r m ef f}= 1$ and theresu l ts o f st o ch ast i c sim ulati on of eq. (\[eq9\_1\_10\] ). 

Ergo dicit y breaking in  higher di m e n sions {# sec4 _ 3} -------------- ----- ---------- - -------- --

E rgodicit y breakin g occurs a lso  in  hi ghe r di mensional GPK m odel sin thepre sence o f a ttr act ive  a nd period ic poten ti al s.  C ons idera  two-dim en sio na l G PK pr o cess $ $d {\ bf x }( t) = {\bf v} ( {\ b f  x}( t) )\, d t +  { \bf b }_{\ c hi_ N(t)} \ , d t
\la bel { eqs_ 1} $$  with $ {\bf b}_\alph a= b \, (\cos  \ phi _\alph a ,  \sin \p hi_\alpha)$, $\phi_\alp h a= 2 \p i ( \alph a-1) /N$, $\al pha =1,\do ts, N $, $A_ {\alph a,\be ta }=1 / N $, an d $\ lam bd a_\alpha=\ l a mbd a$, $ \a lpha ,\beta= 1,\dots, N$, in th e  pr esence of a d ete rmin i s ti c b i as ${\ bf v}( { \ bf x})$ stemmin g from
0$. In_this case,_$b=3/2$, a unique invariant_density admitting_compact_non-atomic support_in_$[x^*,1]$, $x^*=1/2$ appears._From eqs. (\[eq9\_1\_11\]),_(\[eq9\_1\_13\]), after elementary manipulations,_the invariant density_$p_*(x)$_takes the expression $$p_*(x) =\frac{A}{b^2 -v(x)} \, \exp \left [ -2 \lambda 
\int_{x^*}^x \frac{v(\xi)}{v^2(\xi)_-_b^2} \,_d_\xi_\right ]
\label{eq9_2_3}$$ where $A$ is_a normalization constant.

[![Invariant stationary probability_density $p_*(x)$_for the Poisson-Kac scheme (\[eq9\_1\_10\]) with $v(x)$ given_by_eq. (\[eq9\_2\_1\]) at_$D_{\rm eff}=1$, $b=3/2$. The “noisy” line represents the result_of stochastic simulation, the smooth line_represents eq. (\[eq9\_2\_3\]).[]{data-label="Fig25"}](pot_1d.eps_"fig:"){height="6cm"}]{}

Figure_\[Fig25\]_depicts the comparison of_the closed-form expression for the invariant_density at $b=3/2$, $D_{\rm eff}=1$ and_the results of stochastic simulation of eq._(\[eq9\_1\_10\]).

Ergodicity breaking in higher dimensions {#sec4_3}
----------------------------------------

Ergodicity_breaking occurs also in higher_dimensional GPK_models in the presence of_attractive and periodic_potentials. Consider_a two-dimensional GPK_process $$d {\bf x}(t) = _{\bf v}({\bf x}(t))_\, dt + {\bf b}_{\chi_N(t)} \,_d_t
\label{eqs_1}$$ with ${\bf_b}_\alpha=b_\,_(\cos \phi_\alpha,_\sin \phi_\alpha)$, $\phi_\alpha=_2_\pi (\alpha-1)/N$,_$\alpha=1,\dots,N$,_$A_{\alpha,\beta}=1/N$, and $\lambda_\alpha=\lambda$, $\alpha,\beta=1,\dots, N$, in_the_presence of a deterministic bias ${\bf v}({\bf_x})$ stemming from