2p}} \right].\end{aligned}$$

The original @MO96.1 formula is obtained by setting $a = 1$ and $p = 0$ in Eq. (\[eqn:bias\]), while the @SH01.1 revision is obtained with the values $a = 0.75$ and $p = 0.3$.

In non-Gaussian models, the bias acquires an extra scale dependence that can be written as [@MA08.1]

$$b(M,z,k) = b^\mathrm{(G)}(M,z) + \Delta b(M,z,k),$$

where

$$\Delta b(M,z,k) = \left[ b^\mathrm{(G)}(M,z)-1 \right] \delta_\mathrm{c}(z) \Gamma_R(k).$$

The term $\Gamma_R(k)$ encapsulates all the dependence on the scale, and can be written as

$$\begin{aligned}
\Gamma_R(k) &=& \frac{1}{8\pi^2\sigma_M^2\mathcal{M}_R(k)} \int_0^{+\infty} \zeta^2\mathcal{M}_R(\zeta) \times
\nonumber\\
&\times& \left[ \int_{-1}^1 \mathcal{M}_R\left(\sqrt{\alpha}\right) \frac{B_\Phi\left( \zeta,\sqrt{\alpha},k \right)}{P_\Phi(k)} d\mu \right] d\zeta,\end{aligned}$$

where $\alpha \equiv k^2 + \zeta^2 + 2k\zeta\mu$ and $R$ is the top-hat radius corresponding to the mass $M$. In the particular case of a primordial bispectrum of local shape, the previous equation simplifies to

$$\begin{aligned}
\Gamma_R(k) &=& \frac{2f_\mathrm{NL}}{8\pi^2\sigma_M^2\mathcal{M}_R(k)} \int_0^{+\infty} \zeta^2\mathcal{M}_R(\zeta) P_\Phi(\zeta)\times
\nonumber\\
&\times& \left[ \int_{-1}^1 \mathcal{M}_R\left(\sqrt{\alpha}\right) \frac{P_\
2p } } \right].\end{aligned}$$ 

 The original @MO96.1 formula is obtained by setting $ a = 1 $ and $ p = 0 $ in Eq. (\[eqn: bias\ ]), while the @SH01.1 rewrite is prevail with the values $ a = 0.75 $ and $ p = 0.3$. 

 In non - Gaussian model, the diagonal acquires an extra scale addiction that can be write as [ @MA08.1 ] 

 $ $ b(M, z, k) = b^\mathrm{(G)}(M, z) + \Delta b(M, z, k),$$ 

 where 

 $ $ \Delta b(M, z, k) = \left [ b^\mathrm{(G)}(M, z)-1 \right ] \delta_\mathrm{c}(z) \Gamma_R(k).$$ 

 The term $ \Gamma_R(k)$ encapsulates all the dependence on the scale, and can be publish as 

 $ $ \begin{aligned } 
 \Gamma_R(k) & = & \frac{1}{8\pi^2\sigma_M^2\mathcal{M}_R(k) } \int_0^{+\infty } \zeta^2\mathcal{M}_R(\zeta) \times 
 \nonumber\\ 
 & \times & \left [ \int_{-1}^1 \mathcal{M}_R\left(\sqrt{\alpha}\right) \frac{B_\Phi\left (\zeta,\sqrt{\alpha},k \right)}{P_\Phi(k) } d\mu \right ] d\zeta,\end{aligned}$$ 

 where $ \alpha \equiv k^2 + \zeta^2 + 2k\zeta\mu$ and $ R$ is the top - hat radius corresponding to the mass $ M$. In the particular case of a primordial bispectrum of local human body, the previous equation simplifies to 

 $ $ \begin{aligned } 
 \Gamma_R(k) & = & \frac{2f_\mathrm{NL}}{8\pi^2\sigma_M^2\mathcal{M}_R(k) } \int_0^{+\infty } \zeta^2\mathcal{M}_R(\zeta) P_\Phi(\zeta)\times 
 \nonumber\\ 
 & \times & \left [ \int_{-1}^1 \mathcal{M}_R\left(\sqrt{\alpha}\right) \frac{P_\
2p}} \rlght].\end{aligned}$$

The originxl @MO96.1 formula is obtaiied by aetting $x = 1$ and $p = 0$ in Eq. (\[eqn:bias\]), whmle rhe @SY01.1 revision is obtained with the values $q = 0.75$ end $p = 0.3$.

In non-Gaussian models, tgc biav acquires an ewtra scale gependence thad zau be written as [@MA08.1]

$$b(M,z,k) = b^\mathrm{(G)}(M,z) + \Qelta b(K,z,n),$$

where

$$\Delta b(M,e,k) = \ltft[ b^\mafhrm{(G)}(M,z)-1 \right] \delta_\mathrm{c}(z) \Gamma_R(i).$$

The tegm $\Gamma_R(k)$ encapsilates all the dependence ln tje scale, and can bf written aw

$$\begyb{aligned}
\Gammx_R(k) &=& \frac{1}{8\pp^2\vigma_M^2\mathdal{M}_R(k)} \int_0^{+\infty} \zeta^2\mathcal{M}_R(\zega) \tikes
\nonumbet\\
&\cumed& \left[ \int_{-1}^1 \methcal{I}_R\left(\sqrt{\alina}\righd) \frac{B_\Lhi\left( \zeta,\sqvt{\alpia},k \eight)}{P_\Phi(k)} d\mu \right] v\zeta,\end{aligned}$$

where $\alpha \eqgir k^2 + \zeta^2 + 2k\zeta\mu$ ane $E$ is jhe tmp-hag raaiua rordesponfinj to the maas $M$. In the particular case of a itimordial biapectrtm of local shape, the previous equation spmpljfies to

$$\begin{aligned}
\Gamna_R(k) &=& \frac{2f_\mathrm{NL}}{8\pi^2\digma_M^2\matrcal{M}_R(k)} \int_0^{+\infty} \zeta^2\mathcal{M}_R(\zeta) P_\Phi(\zeta)\times
\nmnumbxr\\
&\gimts& \left[ \unh_{-1}^1 \mathcal{M}_R\left(\sqrt{\alpha}\right) \frac{P_\
2p}} \right].\end{aligned}$$ The original @MO96.1 formula is setting = 1$ $p = 0$ @SH01.1 is obtained with values $a = and $p = 0.3$. In non-Gaussian the bias acquires an extra scale dependence that can be written as [@MA08.1] = b^\mathrm{(G)}(M,z) + \Delta b(M,z,k),$$ where $$\Delta b(M,z,k) = \left[ b^\mathrm{(G)}(M,z)-1 \right] \delta_\mathrm{c}(z) The $\Gamma_R(k)$ all dependence on the scale, and can be written as $$\begin{aligned} \Gamma_R(k) &=& \frac{1}{8\pi^2\sigma_M^2\mathcal{M}_R(k)} \int_0^{+\infty} \zeta^2\mathcal{M}_R(\zeta) \times &\times& \left[ \int_{-1}^1 \mathcal{M}_R\left(\sqrt{\alpha}\right) \frac{B_\Phi\left( \zeta,\sqrt{\alpha},k \right)}{P_\Phi(k)} d\mu d\zeta,\end{aligned}$$ where $\alpha \equiv + \zeta^2 + 2k\zeta\mu$ and is top-hat radius to mass In the particular of a primordial bispectrum of local shape, the previous equation simplifies to $$\begin{aligned} \Gamma_R(k) &=& \frac{2f_\mathrm{NL}}{8\pi^2\sigma_M^2\mathcal{M}_R(k)} \int_0^{+\infty} P_\Phi(\zeta)\times \nonumber\\ \int_{-1}^1 \mathcal{M}_R\left(\sqrt{\alpha}\right)
2p}} \right].\end{aligned}$$

The originAl @MO96.1 formulA is obTaiNed By SettIng $a = 1$ And $p = 0$ in Eq. (\[eqn:biaS\]), WhilE the @SH01.1 revision is obtainEd witH tHE valUEs $A = 0.75$ and $p = 0.3$.

in non-GaUSsIAN moDeLs, The BiAS aCquirEs aN extra sCale dependEncE tHat can be writTEn As [@MA08.1]

$$b(M,z,k) = b^\mAthRm{(G)}(M,z) + \Delta b(M,Z,k),$$

wHere

$$\DeLtA b(M,Z,K) = \left[ B^\maThrm{(G)}(m,z)-1 \righT] \Delta_\mAthrm{c}(z) \GaMmA_r(k).$$

The tERm $\Gamma_r(K)$ EnCapsUlates all the depenDEnCE on the scale, and Can be wRiTTeN AS

$$\beGin{Aligned}
\GamMa_r(k) &=& \fraC{1}{8\Pi^2\sigma_m^2\MaTHCAl{M}_r(K)} \int_0^{+\infty} \zeta^2\Mathcal{M}_R(\zeTA) \tiMes
\nonUmBer\\
&\TImes& \leFt[ \int_{-1}^1 \MaTHcaL{M}_R\left(\sqrt{\AlphA}\right) \fraC{B_\Phi\lEFt( \zeta,\sQRt{\alpha},K \right)}{p_\PhI(k)} d\Mu \riGHt] D\zEta,\EnD{AliGNeD}$$

whERe $\aLpha \equiV k^2 + \ZeTa^2 + 2k\zeTa\mu$ AND $r$ Is thE toP-hat RadiuS correspondinG to The mASs $M$. in the PartiCulaR cAse of A primoRdial BiSpectrum of local ShapE, the previOus EqUatIoN simpLIfies tO

$$\beGin{Aligned}
\gamma_R(k) &=& \FRac{2F_\mATHRm{nL}}{8\pi^2\sigma_M^2\mathcal{M}_r(k)} \INT_0^{+\iNfty} \zeta^2\MathcaL{m}_R(\ZeTA) P_\Phi(\zetA)\tImeS
\nonUMBer\\
&\tiMes& \lEFt[ \Int_{-1}^1 \mathcAl{M}_R\leFT(\sQrT{\alpha}\rIgHt) \frac{p_\
2p}} \right].\end{aligned} $$

The or igina l @ MO9 6. 1 fo rmul a is obtainedb y se tting $a = 1$ and $p =  0$ i nE q. ( \ [e qn:bi as\]),w hi l e  th e@S H01 .1 re visio n i s obtai ned with t heva lues $a = 0. 7 5$  and $p =0.3 $.

In non-G aus sian m od els ,  thebia s acq uiresa n extr a scale d ep e ndence that ca n be  wri tten as [@MA08.1] 
$ $ b(M,z,k) = b^\ mathrm {( G )} ( M ,z)  +\Delta b(M ,z ,k),$ $ 

where 
$ $ \ D elt a  b(M,z,k) = \ left[ b^\ma t hrm {(G)}( M, z)- 1  \righ t] \d el t a_\ mathrm{c}(z ) \G amma_R(k) .$$

T h e term$ \Gamma_ R(k)$enc aps ulat e sal l t he dep e nd enc e  on  the sca le ,and c an b e w r itte n a s

$ $\beg in{aligned}
\ Gam ma_R ( k)&=& \ frac{ 1}{8 \p i^2\s igma_M ^2\ma th cal{M}_R(k)} \i nt_0 ^{+\infty } \ ze ta^ 2\ mathc a l{M}_R (\z eta ) \time s
\nonu m ber \\ & \ ti mes& \left[ \int_{ -1 } ^ 1\mathcal {M}_R\ l ef t( \ sqrt{\al ph a}\ righ t )  \fra c{B_ \ Ph i\left(\zeta, \ sq rt {\alpha }, k \rig ht )}{ P_\ Phi(k ) } d\ mu \ri ght] d\z eta,\ e nd{aligned}$$where $\alpha \e q u iv k^2+ \ zeta^2 + 2k \zet a \mu$  and $R $ i s  thetop-h at ra d ius corresponding t othe ma ss $M $. In the par ticular ca s e of a pri mord i al bispectrum oflocal  shape, th e  previou s equ ation si mplifiest o 

$$\beg in{ ali gne d}\ G am ma_R(k) &=& \ f r ac{2 f_ \mathrm {NL }}{8\pi ^2\ sig ma_ M^2 \m athcal{M} _R(k)} \ in t_ 0^ {+ \in fty}\ zeta^2\m at hca l{ M}_ R(\ze t a) P_\ Phi(\ zeta )\ ti m es\nonumb e r\ \ &\ti me s&  \le ft[  \ int_{ -1}^ 1  \m athcal{ M}_R\left (\s q rt{\ al ph a}\righ t) \frac{P_\
2p}} \right].\end{aligned}$$

The_original @MO96.1_formula is obtained by_setting $a_=_1$ and_$p_= 0$ in_Eq. (\[eqn:bias\]), while_the @SH01.1 revision is_obtained with the_values_$a = 0.75$ and $p = 0.3$.

In non-Gaussian models, the bias acquires an extra_scale_dependence that_can_be_written as [@MA08.1]

$$b(M,z,k) = b^\mathrm{(G)}(M,z)_+ \Delta b(M,z,k),$$

where

$$\Delta b(M,z,k) =_\left[ b^\mathrm{(G)}(M,z)-1_\right] \delta_\mathrm{c}(z) \Gamma_R(k).$$

The term $\Gamma_R(k)$ encapsulates all the_dependence_on the scale,_and can be written as

$$\begin{aligned}
\Gamma_R(k) &=& \frac{1}{8\pi^2\sigma_M^2\mathcal{M}_R(k)} \int_0^{+\infty} \zeta^2\mathcal{M}_R(\zeta)_\times
\nonumber\\
&\times& \left[ \int_{-1}^1 \mathcal{M}_R\left(\sqrt{\alpha}\right) \frac{B_\Phi\left( \zeta,\sqrt{\alpha},k_\right)}{P_\Phi(k)} d\mu \right]_d\zeta,\end{aligned}$$

where_$\alpha_\equiv k^2 + \zeta^2_+ 2k\zeta\mu$ and $R$ is the_top-hat radius corresponding to the mass_$M$. In the particular case of a_primordial bispectrum of local shape, the_previous equation simplifies to

$$\begin{aligned}
\Gamma_R(k) &=&_\frac{2f_\mathrm{NL}}{8\pi^2\sigma_M^2\mathcal{M}_R(k)} \int_0^{+\infty}_\zeta^2\mathcal{M}_R(\zeta) P_\Phi(\zeta)\times
\nonumber\\
&\times& \left[ \int_{-1}^1 \mathcal{M}_R\left(\sqrt{\alpha}\right)_\frac{P_\