||_{L^2}}
\label{eq7_4_5}$$ so that $E^*(0)=1$ and $\lim_{t \rightarrow \infty} E^*(t)=0$.

Figure \[Fig13\] depicts the evolution of $E^*(t)$ and at two different values of the Péclet number: $Pe=10^1$ (panel a) and $Pe=10^2$ (panel b). Numerical simulations have been performed by expanding $p_\alpha({\bf x},t)$ in truncated Fourier series, $p_\alpha({\bf x},t)=\sum_{h,k=-N_c}^{N_c} P_{\alpha,h,k} \, e^{i 2 \pi
(h x+ k y)}$, solving the corresponding system of linear differential equations for the Fourier coefficients $P_{\alpha,h,k}$ with an explicit 4-th order Runge-Kutta solver. For the range of $Pe$ values considered $Pe \leq 10^2$, we choose $N_c=50$, which is fully sufficient for an accurate description of the dynamics, apart from the very early stages of the process. These graphs refer to a GPK process with $N=4$ using the transition rate $\lambda$ as parameter. The case $Pe=10^1$ (panel a) is indicative of the typical relaxation properties of GPK systems: the normalized $L^2$-norm $E^*(t)$ decays asymptotically in an exponential way $E^*(t) \sim e^{-\mu(\lambda) \, t}$, but for moderate values of $\lambda$, the decay exponent $\mu(\lambda)$ is a function of the transition rate $\lambda$ and is smaller that the limit value $\mu_\infty= \lim_{\lambda \rightarrow \infty}
\mu(\lambda)$.

[![Normalized square $L^2$-norm $E^*(t)$ vs time $t$ for the GPK flow associated with the standard map ($N=4$) at $Pe=10^1$ (panel a), and $Pe=10^2$ (panel b). Symbols ($\bullet$) represent the solution of the corresponding parabolic advection-diffusion equation (\[eq7\_4\_3\]). The arrows indicate increasing values of $\lambda$. Panel (a): $\lambda= 0.5, \, 1,\, 2,\,4,\,10,\,40$. Panel (b): $\lambda= 0.
||_{L^2 } } 
 \label{eq7_4_5}$$ so that $ E^*(0)=1 $ and $ \lim_{t \rightarrow \infty } E^*(t)=0$. 

 Figure \[Fig13\ ] depicts the evolution of $ E^*(t)$ and at two different values of the Péclet numeral: $ Pe=10 ^ 1 $ (dialog box a) and $ Pe=10 ^ 2 $ (panel b). Numerical simulations have been perform by expanding $ p_\alpha({\bf x},t)$ in truncated Fourier serial, $ p_\alpha({\bf x},t)=\sum_{h, k=-N_c}^{N_c } P_{\alpha, h, k } \, e^{i 2 \pi 
 (h x+ k y)}$, solving the corresponding system of analogue differential equation for the Fourier coefficients $ P_{\alpha, h, k}$ with an denotative 4 - th order Runge - Kutta solver. For the range of $ Pe$ value considered $ Pe \leq 10 ^ 2 $, we choose $ N_c=50 $, which is amply sufficient for an accurate description of the dynamics, apart from the very early phase of the process. These graphs refer to a GPK summons with $ N=4 $ using the transition rate $ \lambda$ as parameter. The event $ Pe=10 ^ 1 $ (panel a) is indicative of the typical relaxation properties of GPK systems: the anneal $ L^2$-norm $ E^*(t)$ decays asymptotically in an exponential way $ E^*(t) \sim e^{-\mu(\lambda) \, t}$, but for moderate values of $ \lambda$, the decay exponent $ \mu(\lambda)$ is a function of the transition rate $ \lambda$ and is smaller that the limit value $ \mu_\infty= \lim_{\lambda \rightarrow \infty } 
 \mu(\lambda)$. 

 [! [ Normalized square $ L^2$-norm $ E^*(t)$ vs time $ t$ for the GPK flow associated with the standard map ($ N=4 $) at $ Pe=10 ^ 1 $ (panel a), and $ Pe=10 ^ 2 $ (panel b). symbol ($ \bullet$) represent the solution of the corresponding parabolic advection - dispersion equality (\[eq7\_4\_3\ ]). The arrow indicate increasing value of $ \lambda$. Panel (a ): $ \lambda= 0.5, \, 1,\, 2,\,4,\,10,\,40$. Panel (b ): $ \lambda= 0.
||_{L^2}}
\lahel{eq7_4_5}$$ so that $E^*(0)=1$ and $\lim_{t \rightarrow \infjy} E^*(t)=0$.

Fignre \[Fig13\] depicts the evolution of $E^*(t)$ and at vwo eiffeeent values of the Pécldt number: $Pe=10^1$ (paneo a) end $Pe=10^2$ (panel b). Nnjerical simulzbions iave been perfotmed by expatding $p_\alpha({\bf x},g)$ nn truncated Fourier series, $p_\alpha({\bf x},t)=\sum_{h,l=-N_f}^{N_c} P_{\alpha,h,k} \, g^{i 2 \pp
(h x+ k j)}$, wolving the corresponding sysfem of ninear differrntial equations for the Flurifr coefficients $P_{\appha,h,k}$ with an qzplicit 4-th ofder Runge-Kutta solver. For the range of $Pe$ values conriderzd $Pe \leq 10^2$, qe chlmse $N_c=50$, whici is fllly sufficiekn for at accurste descriptiok of vhe eynamics, apart from tie very early stages of the psoeess. These graphs refwr to a GPK prozwss wiuh $I=4$ uaing tje vransition date $\lambda$ as parameter. The cssq $Pe=10^1$ (panel a) ia indisaeive of the typical relaxation propertits of GPK systems: the normaluzed $L^2$-norm $E^*(t)$ decays wsymptotisally in an exponential way $E^*(t) \sim e^{-\mu(\lambda) \, t}$, bgt foc oodtrqte vxoufs of $\lambda$, the decay exponent $\mu(\lambda)$ is a rumcnion of the transltion rate $\lambda$ snf of smaller thaj the lnjif value $\mu_\infty= \lil_{\lambda \righrarrow \insty}
\mi(\lambda)$.

[![Normalized square $L^2$-nirm $E^*(t)$ vs tiiw $t$ for the GPK fllw associatzd witn the standard map ($N=4$) at $Pe=10^1$ (pcnel a), and $Pe=10^2$ (panfl b). Symbkus ($\bullet$) represdnt tve solution of the correspjnding pacabolnc advecgion-ciffusyon equatiln (\[eq7\_4\_3\]). The arrows indicatf incteasinc values ov $\lambda$. Panel (a): $\lambda= 0.5, \, 1,\, 2,\,4,\,10,\,40$. Panxk (b): $\lambda= 0.
||_{L^2}} \label{eq7_4_5}$$ so that $E^*(0)=1$ and $\lim_{t E^*(t)=0$. \[Fig13\] depicts evolution of $E^*(t)$ of Péclet number: $Pe=10^1$ a) and $Pe=10^2$ b). Numerical simulations have been performed expanding $p_\alpha({\bf x},t)$ in truncated Fourier series, $p_\alpha({\bf x},t)=\sum_{h,k=-N_c}^{N_c} P_{\alpha,h,k} \, e^{i 2 (h x+ k y)}$, solving the corresponding system of linear differential equations for Fourier $P_{\alpha,h,k}$ an 4-th order Runge-Kutta solver. For the range of $Pe$ values considered $Pe \leq 10^2$, we choose which is fully sufficient for an accurate description the dynamics, apart from very early stages of the These refer to GPK with using the transition $\lambda$ as parameter. The case $Pe=10^1$ (panel a) is indicative of the typical relaxation properties of GPK the normalized decays asymptotically an way \sim e^{-\mu(\lambda) \, for moderate values of $\lambda$, the is a function of the transition rate $\lambda$ is smaller the limit value $\mu_\infty= \lim_{\lambda \rightarrow \mu(\lambda)$. [![Normalized square $L^2$-norm $E^*(t)$ vs time $t$ the GPK flow associated with the standard map ($N=4$) at $Pe=10^1$ (panel a), and $Pe=10^2$ Symbols ($\bullet$) represent the of the corresponding advection-diffusion (\[eq7\_4\_3\]). arrows increasing values $\lambda$. Panel (a): $\lambda= 0.5, \, 1,\, 2,\,4,\,10,\,40$. Panel (b): $\lambda=
||_{L^2}}
\label{eq7_4_5}$$ so that $E^*(0)=1$ and $\lim_{t \rigHtarrow \infTy} E^*(t)=0$.

FIguRe \[FIg13\] DepiCts tHe evolution of $E^*(T)$ And aT two different values of tHe PécLeT NumbER: $PE=10^1$ (paneL a) and $Pe=10^2$ (PAnEL B). NuMeRiCal SiMUlAtionS haVe been pErformed by ExpAnDing $p_\alpha({\bf X},T)$ iN truncated fouRier series, $p_\aLphA({\bf x},t)=\sUm_{H,k=-N_C}^{n_c} P_{\alPha,H,k} \, e^{i 2 \pI
(h x+ k y)}$, sOLving tHe correspOnDIng sysTEm of linEAR dIffeRential equations fOR tHE Fourier coeffiCients $p_{\aLPhA,H,K}$ wiTh aN explicit 4-tH oRder RUNge-KuttA SoLVER. FoR The range of $Pe$ vAlues considERed $pe \leq 10^2$, wE cHooSE $N_c=50$, whiCh is fUlLY suFficient for An acCurate desCriptiON of the dYNamics, aPart frOm tHe vEry eARlY sTagEs OF thE PrOceSS. ThEse graphS rEfEr to a gPK pROCESs wiTh $N=4$ UsinG the tRansition rate $\LamBda$ aS ParAmeteR. The cAse $PE=10^1$ (pAnel a) Is indiCativE oF the typical relaXatiOn propertIes Of gPK SyStems: THe normAliZed $l^2$-norm $E^*(t)$ Decays aSYmpToTICAlLy in an exponential wAy $e^*(T) \SiM e^{-\mu(\lambDa) \, t}$, but FOr MoDErate valUeS of $\LambDA$, The deCay eXPoNent $\mu(\laMbda)$ is A FuNcTion of tHe TransiTiOn rAte $\LambdA$ And iS smallEr that thE limiT Value $\mu_\infty= \liM_{\Lambda \rightarROw \INFtY}
\Mu(\laMbdA)$.

[![Normalized SquaRE $L^2$-noRm $E^*(t)$ VS tIme $T$ For thE GPK fLoW AsSOciated with the standArD map ($N=4$) aT $Pe=10^1$ (paNel a), and $Pe=10^2$ (paneL b). Symbols ($\bULLEt$) represEnt tHE sOLution of the corRespoNding parabOLic advecTion-dIffusion Equation (\[eQ7\_4\_3\]). tHe arrows IndIcaTe iNcrEASiNg values of $\lamBDA$. PanEl (A): $\lambda= 0.5, \, 1,\, 2,\,4,\,10,\,40$. panEl (b): $\lambDa= 0.
||_{L^2}}
\label{eq7_4_5}$ $ so that$E^*( 0)= 1$an d $\ lim_ {t \rightarrow \inf ty} E^*(t)=0$.

Figure  \[Fi g1 3 \] d e pi cts t he evol u ti o n  of  $ E^ *(t )$ an d attwo  differ ent values  of  t he Péclet nu m be r: $Pe=10^ 1$(panel a) an d $ Pe=10^ 2$  (p a nel b ).Numer ical s i mulati ons havebe e n perf o rmed by e xp andi ng $p_\alpha({\bf x} , t)$ in truncat ed Fou ri e rs e rie s,$p_\alpha( {\ bf x} , t)=\sum _ {h , k = -N_ c }^{N_c} P_{\a lpha,h,k} \ ,  e^ {i 2 \ pi 
(h x+ k y )}$,so l vin g the corre spon ding syst em ofl inear d i fferent ial eq uat ion s fo r  t he  Fo ur i erc oe ffi c ien ts $P_{\ al ph a,h,k }$ w i t h an e xpl icit  4-th  order Runge- Kut ta s o lve r. Fo r the  ran ge  of $ Pe$ va luesco nsidered $Pe \l eq 1 0^2$, wecho os e $ N_ c=50$ ,  which  is  fu lly suf ficient for  a n a cc urate descriptionof t he  dynamic s, apa r tfr o m the ve ry  ea rlys t agesof t h eprocess.  These gr ap hs refe rto a G PK  pr oce ss wi t h $N =4$ us ing thetrans i tion rate $\la m bda$ as param e te r .  T h e ca se$Pe=10^1$ ( pane l  a)is i n di cat i ve of  thety p ic a l relaxation proper ti es ofGPK s ystems: the n ormalized$ L ^ 2$-norm$E^* ( t) $  decays asympt otica lly in ane xponenti al wa y $E^*(t ) \sim e^ { - \mu(\lam bda ) \ , t }$, b ut  for moderate v alue sof $\la mbd a$, the  de cay  ex pon en t $\mu(\l ambda)$is  a  f un cti on of the tran si tio nrat e $\l a mbda$and i s sm al le r  th at thel im i t  val ue  $ \mu_ \in ft y= \l im_{ \ lam bda \ri ghtarrow\in f ty}\m u( \lambda )$.

[![Norma li zed square  $ L^2 $-norm $ E^*(t)$vs time $t$ for the GPK flow as soc iated  wit h the sta nda rd map  ($ N =4$) a t $Pe= 10^1$  ( pan e l  a),a n d$Pe =1 0^2$ (pane l b).  Symb ol s ($ \bullet $) represent the s o lut ion of the co rre spon d i ng  pa r ab o lic  a d vec t i on-diffusion eq uation (\[ eq 7 \_ 4\_3\]). T h e a rr ows ind icate i ncrea s ing val ues of $\ lambda$.Pa nel( a ):$\lambda=0.5, \,1,\, 2,\, 4 ,\,10 , \, 40$.Pan el (b) :$\l ambda = 0.
||_{L^2}}
\label{eq7_4_5}$$ so_that $E^*(0)=1$_and $\lim_{t \rightarrow \infty}_E^*(t)=0$.

Figure \[Fig13\]_depicts_the evolution_of_$E^*(t)$ and at_two different values_of the Péclet number:_$Pe=10^1$ (panel a)_and_$Pe=10^2$ (panel b). Numerical simulations have been performed by expanding $p_\alpha({\bf x},t)$ in truncated_Fourier_series, $p_\alpha({\bf_x},t)=\sum_{h,k=-N_c}^{N_c}_P_{\alpha,h,k}_\, e^{i 2 \pi
(h x+_k y)}$, solving the corresponding_system of_linear differential equations for the Fourier coefficients $P_{\alpha,h,k}$_with_an explicit 4-th_order Runge-Kutta solver. For the range of $Pe$ values_considered $Pe \leq 10^2$, we choose_$N_c=50$, which is_fully_sufficient_for an accurate description_of the dynamics, apart from the_very early stages of the process._These graphs refer to a GPK process_with $N=4$ using the transition rate_$\lambda$ as parameter. The case_$Pe=10^1$ (panel_a) is indicative of the_typical relaxation properties_of GPK_systems: the normalized_$L^2$-norm $E^*(t)$ decays asymptotically in an_exponential way $E^*(t)_\sim e^{-\mu(\lambda) \, t}$, but for_moderate_values of $\lambda$,_the_decay_exponent $\mu(\lambda)$_is a function_of_the transition_rate_$\lambda$ and is smaller that the_limit_value $\mu_\infty= \lim_{\lambda \rightarrow \infty}
\mu(\lambda)$.

[![Normalized square $L^2$-norm_$E^*(t)$ vs time $t$_for_the GPK flow associated_with the standard map ($N=4$)_at $Pe=10^1$ (panel a), and $Pe=10^2$_(panel b)._Symbols ($\bullet$)_represent the solution of the corresponding parabolic advection-diffusion equation (\[eq7\_4\_3\]). The_arrows indicate increasing values of $\lambda$._Panel (a): $\lambda= 0.5,_\, 1,\,_2,\,4,\,10,\,40$._Panel (b): $\lambda=_0.