From fd20ab688759f3dfc3f7796ac0b1118e21ee5b22 Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Fri, 20 Oct 2017 10:22:24 +0800 Subject: clusters/halo: Determine magnetic field by fraction of its energy density * Rewrite "magnetic_field()" function to calculate the mean magnetic field strength within ICM according its energy density fraction w.r.t. the ICM thermal energy density. * Remove config options "b_mean" and "b_index", which are replaced with the option "eta_b", the assumed magnetic energy density fraction w.r.t. the ICM thermal energy density. --- fg21sim/extragalactic/clusters/halo.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'fg21sim/extragalactic/clusters/halo.py') diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index e34d91c..e6fe152 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -229,7 +229,7 @@ class RadioHalo: Unit: [uG] """ - return helper.magnetic_field(self.M_obs) + return helper.magnetic_field(mass=self.M_obs, z=self.z_obs) @property def kT_merger(self): @@ -716,6 +716,6 @@ class RadioHalo: """ z = COSMO.redshift(t) mass = self._mass(t) # [Msun] - B = helper.magnetic_field(mass) # [uG] + B = helper.magnetic_field(mass=mass, z=z) # [uG] loss = -4.32e-4 * gamma**2 * ((B/3.25)**2 + (1+z)**4) return loss -- cgit v1.2.2