aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorAaron LI <aly@aaronly.me>2017-07-23 14:01:07 +0800
committerAaron LI <aly@aaronly.me>2017-07-23 14:01:07 +0800
commite727de48313acacb55649893851674bb38fc3c72 (patch)
treec8f587c95e08a19dd0272290ec1390f47d785240
parent4443b6e6d39f586d881b354eead261e7545757b3 (diff)
downloadfg21sim-e727de48313acacb55649893851674bb38fc3c72.tar.bz2
clusters/helper.py: Fix and replace R_halo vs. R_vir relation
Signed-off-by: Aaron LI <aly@aaronly.me>
-rw-r--r--fg21sim/extragalactic/clusters/helper.py23
1 files changed, 12 insertions, 11 deletions
diff --git a/fg21sim/extragalactic/clusters/helper.py b/fg21sim/extragalactic/clusters/helper.py
index 618c237..1367374 100644
--- a/fg21sim/extragalactic/clusters/helper.py
+++ b/fg21sim/extragalactic/clusters/helper.py
@@ -21,6 +21,10 @@ References
.. [cassano2012]
Cassano et al. 2012, A&A, 548, A100
http://adsabs.harvard.edu/abs/2012A%26A...548A.100C
+
+.. [zandanel2014]
+ Zandanel, Pfrommer & Prada 2014, MNRAS, 438, 124
+ http://adsabs.harvard.edu/abs/2014MNRAS.438..124Z
"""
@@ -64,8 +68,13 @@ def radius_halo(mass, z=0.0):
"""
Calculate the radius of (giant) radio halo for a cluster.
- The halo radius is derived from the virial radius using the scaling
- relation in [cassano2007]_.
+ The halo radius is assumed to linearly scale with the virial radius,
+ and is estimated by:
+ R_halo = R_vir / 4
+ * halo radius is ~3-6 times smaller than the virial radius;
+ Ref.[cassano2007],Sec.(1)
+ * halo half radius is ~R500/4, therefore, R_halo ~ R_vir/4;
+ Ref.[zandanel2014],Sec.(6.2)
Parameters
----------
@@ -81,17 +90,9 @@ def radius_halo(mass, z=0.0):
R_halo : float
Radius of the (expected) giant radio halo
Unit: [kpc]
-
- References
- ----------
- Ref.[cassano2007],Fig.(11)
"""
- # slope = 2.63 + np.random.normal(scale=0.5)
- slope = 2.63
- # intercept = 2.3 + np.random.normal(scale=0.05)
- intercept = 2.3
R_vir = radius_virial(mass=mass, z=z) # [kpc]
- R_halo = 10 ** (slope * np.log10(R_vir) + intercept)
+ R_halo = R_vir / 4.0 # [kpc]
return R_halo