From 9631363cf5d3c9d6c37979e0d7a79f1f4c7d849f Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Wed, 26 Jul 2017 23:18:25 +0800 Subject: clusters/emission.py: Fix units/dimensions problem MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit NOTE: Since we just substitute ``p`` (electron momentum) for ``γ`` (electron Lorentz factor), i.e., N = int_p n(p) dp = int_γ n(γ) dγ , therefore the coefficient ``sqrt(3) e^3 B / (m c^2)`` does NOT need modification (e.g., multiply an extra ``m c``)! Signed-off-by: Aaron LI --- fg21sim/extragalactic/clusters/emission.py | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/fg21sim/extragalactic/clusters/emission.py b/fg21sim/extragalactic/clusters/emission.py index 9e097d8..2454644 100644 --- a/fg21sim/extragalactic/clusters/emission.py +++ b/fg21sim/extragalactic/clusters/emission.py @@ -185,7 +185,8 @@ class SynchrotronEmission: # Integrate over energy ``gamma`` in logarithmic grid j_nu = integrate.simps(s1d*self.gamma, np.log(self.gamma)) - coef = np.sqrt(3) * AC.e**3 * self.B / AC.c + B_gauss = self.B * 1e-6 # [uG] -> [G] + coef = np.sqrt(3) * AC.e**3 * B_gauss / AU.mec2 j_nu *= coef return j_nu -- cgit v1.2.2