picasso.nonthermal
: Non-thermal pressure support#
|
Non-thermal pressure fraction computed from the Shaw+10 model |
|
Non-thermal pressure fraction computed from an adapted version of the Nelson+14 model (see Notes) |
|
Generic expression for non-thermal pressure fraction: a power law evolution with radius, plus a constant plateau (see Notes) |
- picasso.nonthermal.f_nt_shaw10(r_R500, z, alpha_0=0.18, beta=0.5, n_nt=0.8)[source]#
Non-thermal pressure fraction computed from the Shaw+10 model
- Parameters:
r_R500 (array-like) – Radii normalized to R500c
alpha_0 (float, optional) – Non-thermal pressure fraction at z=0 and r=R500c, by default 0.18
beta (float, optional) – Redshift evolution rate, by default 0.5
n_nt (float, optional) – Power-law radial dependence, by default 0.8
- Returns:
Non-thermal pressure fraction for each radius
- Return type:
array-like
- picasso.nonthermal.f_nt_nelson14(r_RDelta, A, B, C)[source]#
Non-thermal pressure fraction computed from an adapted version of the Nelson+14 model (see Notes)
- Parameters:
r_RDelta (array-like) – Radii normalized to RDelta
A (float) – \(A\) parameter values
B (float) – \(B\) parameter values
C (float) – \(C\) parameter values
- Returns:
Non-thermal pressure fraction for each radius
- Return type:
array-like
Notes
The model is computed as:
\[f_{nt} = 1 - A \times \left\{ 1 + \exp \left[ -\left( \frac{r}{B \times R_\Delta} \right)^{C} \right] \right\}. \]This is a modified version of the original model (see eq. 7 in Nelson+14), where the radius is expressed in units of RDelta instead of R200m.
- picasso.nonthermal.f_nt_generic(r_RDelta, a, b, c)[source]#
Generic expression for non-thermal pressure fraction: a power law evolution with radius, plus a constant plateau (see Notes)
- Parameters:
r_RDelta (array-like) – Radii normalized to RDelta
A (float) – Non-thermal pressure fraction in the cluster center
B (float) – Non-thermal pressure fraction at r=RDelta
C (float) – Power law radial dependence of f_nt
- Returns:
Non-thermal pressure fraction for each radius
- Return type:
array-like
Notes
The model is computed as:
\[f_{nt} = A + (B-A) \left(\frac{r}{R_\Delta}\right)^C \]