| Commit message (Collapse) | Author | Age | Files | Lines |
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
| |
This method simulates the radio halo properties for each cluster with a
recent major merger.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
|
| |
* Methods "power()", "flux()" already moved to "halo.py";
* Method "brightness()" will also be moved over there;
* Clean up thus unused parameters "z" and "radius", as well as imports.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
| |
Calculate the synchrotron power (a.k.a. specific luminosity) and flux
density from the emissivity for the halo.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
| |
This method invokes ``SynchrotronEmission`` to calculate the synchrotron
emissivity at specified frequencies for the derived electron spectrum.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
| |
Update method "_calc_ratio_a()" to "_calc_factor_a()"
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
| |
This makes the simulation agree better with observations,
e.g., <Jsyn> ~ 1e-42 [erg/s/cm^3/Hz]
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
|
|
|
|
| |
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 <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
|
|
|
| |
Use 2D grid of (discrete) samples to optimize integration speed,
avoiding the more complicated integration w.r.t. functions.
WARNING:
Current calculation results seems wrong!
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
| |
Use interpolation to optimize the speed as well as to vectorize the
function to ease calling.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
|
|
|
| |
This configuration option "beta_turb" has different meaning as
"eta_turb" in [Cassano & Brunetti 2005], though they can both be an
indicator for the turbulence acceleration efficiency.
Rename this configuration option to make the difference more clear.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
|
|
|
| |
Zero or negative diffusion coefficient leads to unstable or wrong
results due to numerical algorithm/scheme adopted to solve the
Fokker-Planck equation.
Also add a NOTE to the FokkerPlanckSolver class.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
| |
Also update copyright year/email and minor cleanups
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
| |
Update to use the global CONFIGS and COSMO
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
| |
Import the global "configs" in helper.py to make parameters lean.
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
| |
Still WIP...
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
| |
References: Komatsu et al. 2011, ApJS, 192, 18; Tab.(1)
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
| |
Signed-off-by: Aaron LI <aly@aaronly.me>
|
|
|
|
|
|
|
| |
Several methods/functions have been migrated into "helper.py", while
other methods/functions are obsolete.
Signed-off-by: Aaron LI <aly@aaronly.me>
|