aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
-rw-r--r--fg21sim/configs/config.spec10
-rw-r--r--fg21sim/extragalactic/clusters/halo.py5
2 files changed, 5 insertions, 10 deletions
diff --git a/fg21sim/configs/config.spec b/fg21sim/configs/config.spec
index 1ba8dd3..adc4134 100644
--- a/fg21sim/configs/config.spec
+++ b/fg21sim/configs/config.spec
@@ -374,19 +374,15 @@ stream = option("stderr", "stdout", "", default="stderr")
# Giant radio halos
#
[[halos]]
- # A custom parameter to tune the turbulent acceleration efficiency.
- # NOTE: The smaller this parameter, the shorter the acceleration
- # timescale, therefore more efficient acceleration.
+ # A custom factor to tune the turbulent acceleration efficiency.
+ # NOTE: This parameter incorporates the efficiency factor describing
+ # the effectiveness of the ICM plasma instabilities.
f_acc = float(default=1.0, min=0.1, max=10)
# The factor that is multiplied to the turbulence injection radius
# to derive the radio halo radius.
f_radius = float(default=1.0, min=0.5, max=10)
- # An efficiency factor describing the effectiveness of plasma
- # instabilities (e.g., due to spatial or temporal intermittency).
- zeta_ins = float(default=0.1, min=0.1, max=1.0)
-
# The fraction of merger energy transferred into the turbulence.
eta_turb = float(default=0.1, min=0.1, max=0.5)
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index 81b9c0d..4e87f3f 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -151,7 +151,6 @@ class RadioHalo1M:
self.configs = configs
self.f_acc = configs.getn(comp+"/f_acc")
self.f_radius = configs.getn(comp+"/f_radius")
- self.zeta_ins = configs.getn(comp+"/zeta_ins")
self.eta_turb = configs.getn(comp+"/eta_turb")
self.eta_e = configs.getn(comp+"/eta_e")
self.x_cr = configs.getn(comp+"/x_cr")
@@ -334,9 +333,9 @@ class RadioHalo1M:
k_L = 2 * np.pi / L_turb
cs = helper.speed_sound(self.kT(t_merger)) # [km/s]
v_t = self._velocity_turb(t_merger) # [km/s]
- tau = self.x_cr * cs**3 / (8*k_L * self.zeta_ins * v_t**4)
+ tau = self.x_cr * cs**3 / (8*k_L * v_t**4)
tau *= AUC.s2Gyr * AUC.kpc2km # [s kpc/km] -> [Gyr]
- tau *= self.f_acc # custom tune parameter
+ tau /= self.f_acc # tune factor (folded with "zeta_ins")
# Impose the maximum acceleration timescale
if tau > tau_max: