P\_[nores]{}(x,t|x\_0=0,0)=\
&&(-),\
ł[eq:parabolic-resetting-Pxt-x00]{} and P\_[res]{}(t|x\_0=0) =. ł[eq:parabolic-resetting-Pt-x00]{} Equation (\[eq:parabolic-resetting-Pt-x00\]) may be put in the scaling form P\_[res]{}(t|x\_0=0)=[G]{}(t), with ${\cal G}(y) = \tanh(y)^{3/2}/\sqrt{\sinh(y)}$. Equation (\[eq:parabolic-resetting-Pt-x00\]) yields the mean first-reset time $\langle t \rangle_{\rm res}$ for $x_0=0$ to be given by t \_[res]{}=, ł[eq:tav-x00]{} where $\Gamma$ is the Gamma function. Equations (\[eq:parabolic-resetting-Pt-x00\]) and (\[eq:tav-x00\]) yield the stationary spatial distribution on using Eq. (\[eq:Pxstat-final\]): &&P\_[st]{}(x|x\_0=0)=\
&&\_0\^(-).\
&&=()\^[1/4]{}K\_[1/4]{}(), ł[eq:stationary-px-parabolic-resetting]{} where $K_n(x)$ is the $n-$th order modified Bessel function of the second kind. Equation (\[eq:stationary-px-parabolic-resetting\]) implies that the stationary distribution is symmetric around $x=0$, which is expected since the resetting rate is symmetric around $x_0=0$. The stationary distribution (\[eq:stationary-px-parabolic-resetting\]) may be put in the scaling form P\_[st]{}(x|x\_0=0)=(), where the scaling function is given by ${\cal R}(y)=(y^2/2)^{1/4}K_{1/4}(y^2/2)$.

![**Theory versus simulation results for the stationary spatial distribution of a free Brownian particle undergoing parabolic resetting.** The points denote simulation results, while the line stands for the exact analytical expression given in Eq. (\[eq:stationary-px-parabolic-resetting\]). The numerical results were obtained from $10^4$ independent simulations of
P\_[nores]{}(x, t|x\_0=0,0)=\ 
 & & (-),\ 
 ł[eq: parabolic - resetting - Pxt - x00 ] { } and P\_[res]{}(t|x\_0=0) =. ł[eq: parabolic - resetting - Pt - x00 ] { } Equation (\[eq: parabolic - resetting - platinum - x00\ ]) may be place in the scaling form P\_[res]{}(t|x\_0=0)=[G]{}(t), with $ { \cal G}(y) = \tanh(y)^{3/2}/\sqrt{\sinh(y)}$. Equation   (\[eq: parabolic - resetting - platinum - x00\ ]) yields the beggarly first - reset prison term $ \langle triiodothyronine \rangle_{\rm res}$ for $ x_0=0 $ to be given by t \_[res]{}=, ł[eq: tav - x00 ] { } where $ \Gamma$ is the Gamma affair. Equations (\[eq: parabolic - resetting - platinum - x00\ ]) and (\[eq: tav - x00\ ]) yield the stationary spatial distribution on using Eq.   (\[eq: Pxstat - final\ ] ): & & P\_[st]{}(x|x\_0=0)=\ 
 & & \_0\^(-).\ 
 & & =() \^[1/4]{}K\_[1/4 ] { } (), ł[eq: stationary - px - parabolic - resetting ] { } where $ K_n(x)$ is the $ n-$th order modified Bessel affair of the second kind. Equation   (\[eq: stationary - px - parabolic - resetting\ ]) implies that the stationary distribution is symmetric around $ x=0 $, which is expected since the resetting pace is symmetric around $ x_0=0$. The stationary distribution (\[eq: stationary - px - parabolic - resetting\ ]) may be put in the scaling shape P\_[st]{}(x|x\_0=0)= (), where the scaling function is given by $ { \cal R}(y)=(y^2/2)^{1/4}K_{1/4}(y^2/2)$. 

! [ * * hypothesis versus simulation solution for the stationary spatial distribution of a free Brownian atom undergoing parabolic resetting. * * The points denote simulation results, while the cable stands for the exact analytical expression given in Eq. (\[eq: stationary - px - parabolic - resetting\ ]). The numeric results were obtained from $ 10 ^ 4 $ independent simulations of
P\_[noges]{}(x,t|x\_0=0,0)=\
&&(-),\
ł[eq:parabolic-resettikg-Pxt-x00]{} and P\_[res]{}(t|r\_0=0) =. ł[eq:pacabolic-desettine-Pt-x00]{} Equation (\[eq:parabolic-resxttibg-Pt-x00\]) may be put in the scauing form P\_[res]{}(t|x\_0=0)=[G]{}(r), wiuh ${\cal G}(y) = \tanh(y)^{3/2}/\sqrt{\sinh(y)}$. Equatjln (\[eq:'acabolic-resettinn-Pt-x00\]) yields the mean firsd-rdszt time $\langle t \rangle_{\rm res}$ for $x_0=0$ eo be govfn by t \_[res]{}=, ł[eq:jav-x00]{} erere $\Gamma$ is the Gamma function. Equatjons (\[eq:karabolic-resetting-Lt-x00\]) and (\[eq:tav-x00\]) yield the shatilnary spatial distgibution on usigt Eq. (\[eq:Pxstat-winal\]): &&P\_[st]{}(x|q\_0=0)=\
&&\_0\^(-).\
&&=()\^[1/4]{}N\_[1/4]{}(), ł[eq:statiohary-px-parabolic-resetting]{} where $Y_n(x)$ ix the $n-$th ireer kodified Bewsel sunction of bne secmnd kinc. Equation (\[eq:ststimnaey-px-parabolic-resettinj\]) implies that the sjationary giatribution is symnerric droutd $x=0$, qhizh js espectef smnce the reaetting ratw is symmetric aroumd $q_0=0$. The stationzry diftwibution (\[eq:stationary-px-parabolic-resettitg\]) jay be put in the scalibg form P\_[st]{}(x|x\_0=0)=(), where tje scalind function is given by ${\cal R}(y)=(y^2/2)^{1/4}K_{1/4}(y^2/2)$.

![**Theory versus simglatikv rtsmjgw vor the stationary spatial distribution of a sdet Bgownian particle mndergoing paraboloc rrfetting.** The puints bsnkte simulation resklts, whyle tye line suands for the exact analytical ezpression ginen un Eq. (\[eq:stationary-'x-parabolic-rzsettimg\]). Thr numerical results werz obtajned from $10^4$ lndependehg simulations of
P\_[nores]{}(x,t|x\_0=0,0)=\ &&(-),\ ł[eq:parabolic-resetting-Pxt-x00]{} and P\_[res]{}(t|x\_0=0) =. ł[eq:parabolic-resetting-Pt-x00]{} may put in scaling form P\_[res]{}(t|x\_0=0)=[G]{}(t), Equation yields the mean time $\langle t res}$ for $x_0=0$ to be given t \_[res]{}=, ł[eq:tav-x00]{} where $\Gamma$ is the Gamma function. Equations (\[eq:parabolic-resetting-Pt-x00\]) and (\[eq:tav-x00\]) the stationary spatial distribution on using Eq. (\[eq:Pxstat-final\]): &&P\_[st]{}(x|x\_0=0)=\ &&\_0\^(-).\ &&=()\^[1/4]{}K\_[1/4]{}(), ł[eq:stationary-px-parabolic-resetting]{} where is $n-$th modified function of the second kind. Equation (\[eq:stationary-px-parabolic-resetting\]) implies that the stationary distribution is symmetric around $x=0$, is expected since the resetting rate is symmetric $x_0=0$. The stationary distribution may be put in the form where the function given ${\cal R}(y)=(y^2/2)^{1/4}K_{1/4}(y^2/2)$. ![**Theory simulation results for the stationary spatial distribution of a free Brownian particle undergoing parabolic resetting.** The points simulation results, line stands the analytical given in Eq. numerical results were obtained from $10^4$
P\_[nores]{}(x,t|x\_0=0,0)=\
&&(-),\
ł[eq:parabolic-reseTting-Pxt-x00]{} aNd P\_[reS]{}(t|x\_0=0) =. Ł[eq:PaRaboLic-rEsetting-Pt-x00]{} EquATion (\[Eq:parabolic-resetting-Pt-X00\]) may bE pUT in tHE sCalinG form P\_[rES]{}(t|X\_0=0)=[g]{}(T), wiTh ${\CaL G}(y) = \TaNH(y)^{3/2}/\Sqrt{\sInh(Y)}$. EquatiOn (\[eq:paraboLic-ReSetting-Pt-x00\]) yiELdS the mean fiRst-Reset time $\lanGle T \ranglE_{\rM reS}$ For $x_0=0$ tO be Given By t \_[res]{}=, Ł[Eq:tav-x00]{} Where $\GammA$ iS The GamMA functiON. eqUatiOns (\[eq:parabolic-resETtINg-Pt-x00\]) and (\[eq:tav-x00\]) Yield tHe STaTIOnaRy sPatial distRiButioN On using eQ. (\[eQ:pXStaT-Final\]): &&P\_[st]{}(x|x\_0=0)=\
&&\_0\^(-).\
&&=()\^[1/4]{}K\_[1/4]{}(), ł[eQ:stationary-PX-paRaboliC-rEseTTing]{} whEre $K_n(X)$ iS The $N-$th order modIfieD Bessel fuNction OF the secONd kind. EQuatioN (\[eq:StaTionARy-Px-ParAbOLic-REsEttINg\]) iMplies thAt ThE statIonaRY DIStriButIon iS symmEtric around $x=0$, wHicH is eXPecTed siNce thE resEtTing rAte is sYmmetRiC around $x_0=0$. The statIonaRy distribUtiOn (\[Eq:sTaTionaRY-px-parAboLic-ResettiNg\]) may be PUt iN tHE SCaLing form P\_[st]{}(x|x\_0=0)=(), where ThE SCaLing funcTion is GIvEn BY ${\cal R}(y)=(y^2/2)^{1/4}K_{1/4}(Y^2/2)$.

![**THeoRy veRSUs simUlatIOn Results fOr the sTAtIoNary spaTiAl distRiButIon Of a frEE BroWnian pArticle uNdergOIng parabolic reSEtting.** The poinTS dENOtE SimuLatIon results, wHile THe liNe stANdS foR The exAct anAlYTiCAl expression given in eq. (\[Eq:statIonarY-px-parabolic-rEsetting\]). ThE NUMerical rEsulTS wERe obtained from $10^4$ IndepEndent simuLAtions of
P\_[nores]{}(x,t|x\_0=0,0) =\
&&(-),\ 
ł[eq :pa rab ol ic-r eset ting-Pxt-x00]{ }  and  P\_[res]{}(t|x\_0=0)=. ł[ eq : para b ol ic-re setting - Pt - x 00] {}  E qua ti o n(\[eq :pa rabolic -resetting -Pt -x 00\]) may be pu t in the s cal ing form P\_ [re s]{}(t |x \_0 = 0)=[G ]{} (t),with $ { \cal G }(y) = \t an h (y)^{3 / 2}/\sqr t { \s inh( y)}$. Equation (\ [ eq : parabolic-rese tting- Pt - x0 0 \ ])yie lds the me an  firs t -resett im e $ \la n gle t \rangle _{\rm res}$ for  $x_0= 0$  to be giv en by  t \_[ res]{}=, ł[ eq:t av-x00]{}  where $\Gamma $  is the  Gamma  fu nct ion. Eq ua tio ns (\[ e q: par a bol ic-reset ti ng -Pt-x 00\] ) a n d (\ [eq :tav -x00\ ]) yield thesta tion a ryspati al di stri bu tionon usi ng Eq .(\[eq:Pxstat-fi nal\ ]): &&P\_ [st ]{ }(x |x \_0=0 ) =\
&&\ _0\ ^(- ).\
&&= ()\^[1/ 4 ]{} K\ _ [ 1 /4 ]{}(), ł[eq:statio na r y -p x-parabo lic-re s et ti n g]{} whe re  $K _n(x ) $  is t he $ n -$ th order  modif i ed  B essel f un ctionof  th e s econd kind . Equa tion (\[ eq:st a tionary-px-par a bolic-resetti n g\ ] )  i m plie s t hat the sta tion a ry d istr i bu tio n  is s ymmet ri c  a r ound $x=0$, which i sexpect ed si nce the reset ting ratei s symmetri c ar o un d  $x_0=0$. Thestati onary dist r ibution(\[eq :station ary-px-pa r a bolic-re set tin g\] ) m a y  b e put in thes c alin gform P\ _[s t]{}(x| x\_ 0=0 )=( ),wh ere the s caling f un ct io nisgiven by ${\ca lR}( y) =(y ^2/2) ^ {1/4}K _{1/4 }(y^ 2/ 2) $ .
![**The o ry v ersu ssi mula tio nresul ts f o r t he stat ionary sp ati a l di st ri butionof a free Bro wn ian partic le  un dergoi n g  parabol ic resetting.** The poi n ts deno tesimul atio n results , w hile t hel ine st ands f or th eexa c t  anal y t ic alex pression g i v enin Eq .(\[e q:stati onary-px-parabolic - res etting\]). Th e n umer i c al  re s ul t s w er e  ob t a ined from $10^4 $ independ en t  s imulations of
P\_[nores]{}(x,t|x\_0=0,0)=\
&&(-),\
ł[eq:parabolic-resetting-Pxt-x00]{} and_P\_[res]{}(t|x\_0=0) =._ł[eq:parabolic-resetting-Pt-x00]{} Equation (\[eq:parabolic-resetting-Pt-x00\]) may_be put_in_the scaling_form_P\_[res]{}(t|x\_0=0)=[G]{}(t), with ${\cal_G}(y) = \tanh(y)^{3/2}/\sqrt{\sinh(y)}$._Equation (\[eq:parabolic-resetting-Pt-x00\]) yields the mean_first-reset time $\langle_t_\rangle_{\rm res}$ for $x_0=0$ to be given by t \_[res]{}=, ł[eq:tav-x00]{} where $\Gamma$ is_the_Gamma function._Equations_(\[eq:parabolic-resetting-Pt-x00\])_and (\[eq:tav-x00\]) yield the stationary_spatial distribution on using Eq. (\[eq:Pxstat-final\]):_&&P\_[st]{}(x|x\_0=0)=\
&&\_0\^(-).\
&&=()\^[1/4]{}K\_[1/4]{}(), ł[eq:stationary-px-parabolic-resetting]{}_where $K_n(x)$ is the $n-$th order modified Bessel_function_of the second_kind. Equation (\[eq:stationary-px-parabolic-resetting\]) implies that the stationary distribution is symmetric_around $x=0$, which is expected since_the resetting rate_is_symmetric_around $x_0=0$. The stationary_distribution (\[eq:stationary-px-parabolic-resetting\]) may be put in_the scaling form P\_[st]{}(x|x\_0=0)=(), where the_scaling function is given by ${\cal R}(y)=(y^2/2)^{1/4}K_{1/4}(y^2/2)$.

![**Theory_versus simulation results for the stationary_spatial distribution of a free_Brownian particle_undergoing parabolic resetting.** The points_denote simulation results,_while the_line stands for_the exact analytical expression given in_Eq. (\[eq:stationary-px-parabolic-resetting\]). The_numerical results were obtained from $10^4$_independent_simulations of