diff options
Diffstat (limited to 'fg21sim/extragalactic/clusters/helper.py')
| -rw-r--r-- | fg21sim/extragalactic/clusters/helper.py | 16 |
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"), ): """ |