aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/extragalactic
diff options
context:
space:
mode:
authorAaron LI <aly@aaronly.me>2017-10-19 21:47:19 +0800
committerAaron LI <aly@aaronly.me>2017-10-19 21:47:19 +0800
commit6938f435b70bb380a438e400f4f1cbdcd92e2faa (patch)
tree725fc413d4ba6f1631939dc6d85ae98261df5638 /fg21sim/extragalactic
parentd680f1714673e93f330b2c6f0d2be5fb87d340ff (diff)
downloadfg21sim-6938f435b70bb380a438e400f4f1cbdcd92e2faa.tar.bz2
clusters/halo: Add kT_cluster() to calculate the cluster ICM temperature
Diffstat (limited to 'fg21sim/extragalactic')
-rw-r--r--fg21sim/extragalactic/clusters/helper.py25
1 files changed, 25 insertions, 0 deletions
diff --git a/fg21sim/extragalactic/clusters/helper.py b/fg21sim/extragalactic/clusters/helper.py
index 20ad9d9..daf806b 100644
--- a/fg21sim/extragalactic/clusters/helper.py
+++ b/fg21sim/extragalactic/clusters/helper.py
@@ -139,6 +139,31 @@ def kT_virial(mass, z=0.0, radius=None):
return kT
+def kT_cluster(mass, z=0.0, radius=None):
+ """
+ Calculate the temperature of a cluster ICM.
+
+ NOTE
+ ----
+ When a cluster forms, there are accretion shocks forms around
+ the cluster (near the virial radius) which can heat the gas,
+ therefore the ICM has a higher temperature than the virial
+ temperature, which can be estimated as:
+ T_icm = T_vir + 1.5*T_1 ~ T_vir + 0.8 [keV]
+
+ Reference: Ref.[fujita2003],Eq.(49)
+
+ Returns
+ -------
+ kT_icm : float
+ The temperature of the cluster ICM.
+ Unit: [keV]
+ """
+ kT_vir = kT_virial(mass=mass, z=z, radius=radius)
+ kT_icm = kT_vir + 0.8
+ return kT_icm
+
+
def mass_to_kT(mass, z=0.0):
"""
Calculate the cluster ICM temperature from its mass using the