aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/extragalactic/clusters/helper.py
diff options
context:
space:
mode:
Diffstat (limited to 'fg21sim/extragalactic/clusters/helper.py')
-rw-r--r--fg21sim/extragalactic/clusters/helper.py16
1 files changed, 14 insertions, 2 deletions
diff --git a/fg21sim/extragalactic/clusters/helper.py b/fg21sim/extragalactic/clusters/helper.py
index 18e7133..c4fdecc 100644
--- a/fg21sim/extragalactic/clusters/helper.py
+++ b/fg21sim/extragalactic/clusters/helper.py
@@ -10,6 +10,10 @@ References
Arnaud, Pointecouteau & Pratt 2005, A&A, 441, 893;
http://adsabs.harvard.edu/abs/2005A%26A...441..893
+.. [beck2005]
+ Beck & Krause 2005, AN, 326, 414
+ http://adsabs.harvard.edu/abs/2005AN....326..414B
+
.. [cassano2005]
Cassano & Brunetti 2005, MNRAS, 357, 1313
http://adsabs.harvard.edu/abs/2005MNRAS.357.1313C
@@ -353,7 +357,7 @@ def density_energy_electron(n_e, gamma):
def magnetic_field(
mass,
z=0.0,
- eta_b=CONFIGS.getn("extragalactic/clusters/eta_b"),
+ eta_b=CONFIGS.getn("extragalactic/halos/x_cr"),
kT_out=CONFIGS.getn("extragalactic/clusters/kT_out"),
):
"""
@@ -366,6 +370,14 @@ def magnetic_field(
Magnetic field energy density: u_B = B^2 / (8π),
where "B" in units of [G], then "u_B" has unit of [erg/cm^3].
+ NOTE
+ ----
+ Magnetic fields and cosmic rays are strongly coupled and exchange
+ energy. Therefore equipartition between them is assumed, i.e.,
+ X_cr (= ε_cr / ε_th) = η_b (= ε_b / ε_th)
+
+ Reference: [beck2005],App.A
+
Returns
-------
B : float
@@ -380,7 +392,7 @@ def magnetic_field(
def plasma_beta(
mass,
z=0.0,
- eta_b=CONFIGS.getn("extragalactic/clusters/eta_b"),
+ eta_b=CONFIGS.getn("extragalactic/halos/x_cr"),
kT_out=CONFIGS.getn("extragalactic/clusters/kT_out"),
):
"""