aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/extragalactic/clusters/halo.py
Commit message (Collapse)AuthorAgeFilesLines
* clusters/halo: Rewrite turbulence acceleration & diffusion coefficient!Aaron LI2017-10-211-57/+90
| | | | | | | | | | | | The previous calculation of the turbulence acceleration efficiency has significant problems, which is too low, i.e., the acceleration timescale is too large (~10 Gyr). However, it is reported that turbulent acceleration has a timescale ~100 Myr (0.1 Gyr). I believe this problem is due to the over-simplification to the formula proposed by [cassano2005]. Rewrite the turbulence acceleration and diffusion coefficient calculations, adopting the method proposed by [brunetti2016]. To this end, two new options "f_lturb" and "f_acc" are introduced to tune the results.
* clusters/halos: Add more cluster and halo propertiesAaron LI2017-10-211-7/+30
| | | | | | * Rvir0, Rvir_main, Rvir_sub * kT_main, kT_sub (replace kT_merger) * radius -> Rhalo, angular_radius -> Rhalo_angular, B -> B0
* clusters/halo: Determine magnetic field by fraction of its energy densityAaron LI2017-10-201-2/+2
| | | | | | | | | | * Rewrite "magnetic_field()" function to calculate the mean magnetic field strength within ICM according its energy density fraction w.r.t. the ICM thermal energy density. * Remove config options "b_mean" and "b_index", which are replaced with the option "eta_b", the assumed magnetic energy density fraction w.r.t. the ICM thermal energy density.
* clusters/halo: Add "kT_obs" and "tback_merger" propertiesAaron LI2017-10-191-1/+16
|
* clusters/halo: Use "kT_cluster()" instead of "mass_to_kT()"Aaron LI2017-10-191-1/+1
|
* clusters/halo: rewrite power/flux/Tb calc with k-correction consideredAaron LI2017-10-191-25/+94
|
* clusters/halo: Add set_electron_spectrum() method; update doc a bitAaron LI2017-10-191-1/+26
|
* clusters/halo: Improve time step determination a bitAaron LI2017-10-191-2/+5
|
* clusters/halo: Increase the time step, but require a minimum of 20 stepsAaron LI2017-10-181-0/+7
|
* clusters: Update radio halo radius estimationAaron LI2017-10-181-4/+8
| | | | | | | | | The simulated radio halo is assumed to have a radius of the falling sub-cluster; while previously it is assumed to be 1/4 of the virial radius of the merged cluster. The new estimation can agree better with the currently observed radio halos, which generally have a angular diameter size of ~2-7 [arcmin].
* clusters: Update the documents/comments a b itAaron LI2017-10-141-2/+2
|
* clusters/halo.py: Update injection rate calculationAaron LI2017-08-261-10/+5
| | | | | The electrons are assumed to be injected throughout the whole cluster ICM/volume.
* clusters/halo.py: Rewrite electron acceleration coefficient calcAaron LI2017-08-141-26/+51
| | | | | | | | | | * Adopt the electron acceleration coefficient formula from [cassano2005] * Rename method "_tau_acceleration()" to "_chi_acceleration()", and rewrite * Add property "kT_merger" * Also save "kT_merger" and "chi" into halos data Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo.py: Simplify two coefficients usageAaron LI2017-08-141-7/+5
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters: Rename property "injection_rate" and save into halos dataAaron LI2017-08-141-37/+37
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters: Implement calc_{power,flux,brightness_mean} in helperAaron LI2017-08-111-40/+10
| | | | | | Help development & debug ... Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo: "volume" use unit "kpc^3" instead of "cm^3"Aaron LI2017-08-111-3/+6
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* Fix a bug in clusters/halo.py; Add a comment to utils/io/_create_dir()Aaron LI2017-08-061-1/+1
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo.py: Finished brief simulation descriptionAaron LI2017-08-011-3/+13
|
* clusters: Fix several small issues on wrong variable namesAaron LI2017-07-311-2/+2
| | | | | | Also update logging message for RadioHalo a bit. Signed-off-by: Aaron LI <aly@aaronly.me>
* Add `fg21sim.share` module to hold globally shared "CONFIGS" & "COSMO"Aaron LI2017-07-311-2/+1
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo.py: Add description for "time_crossing" propertyAaron LI2017-07-311-0/+6
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* cluster/halo.py: Add method "calc_brightness_mean()"Aaron LI2017-07-301-0/+48
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo.py: Add property "angular_radius" (unit: [arcsec])Aaron LI2017-07-301-0/+10
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo.py: Add methods "calc_power()" and "calc_flux()"Aaron LI2017-07-291-0/+50
| | | | | | | 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>
* clusters/halo.py: Add method "calc_emissivity()"Aaron LI2017-07-291-6/+37
| | | | | | | 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>
* clusters/halo.py: Add property "volume" (unit: [cm^3])Aaron LI2017-07-291-0/+9
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo.py: Add property "magnetic_field"Aaron LI2017-07-261-0/+11
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo: Change config "eta_turb" to "beta_turb"Aaron LI2017-07-251-3/+3
| | | | | | | | | | 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>
* clusters: Accept 1D numpy array and calculate values for all gamma'sAaron LI2017-07-231-17/+16
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo.py: Constrain tau_acc to avoid zero diffusionAaron LI2017-07-231-2/+15
| | | | | | | | | | 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>
* Use all uppercase to identify the global CONFIGS and COSMOAaron LI2017-07-221-15/+11
| | | | | | Update to use the global CONFIGS and COSMO Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters: Import global "configs" to simplify parametersAaron LI2017-07-221-52/+54
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo.py: Implement "_injection_rate" propertyAaron LI2017-07-221-6/+56
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo.py: update parameters and methodsAaron LI2017-07-221-46/+119
| | | | | | Still WIP... Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo.py: More cleanups with minor updatesAaron LI2017-07-221-70/+97
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo.py: Significant cleanupsAaron LI2017-07-221-551/+2
| | | | | | | Several methods/functions have been migrated into "helper.py", while other methods/functions are obsolete. Signed-off-by: Aaron LI <aly@aaronly.me>
* Use the global "cosmo" instance for simplificationAaron LI2017-07-211-26/+20
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* halo.py: Fix wrong cosmo.age() usageAaron LI2017-06-251-1/+1
|
* halo.py: Add argument "n0_e" to method "calc_electron_spectrum()"Aaron LI2017-06-241-4/+10
|
* halo.py: Split out property "_chi_data"Aaron LI2017-06-241-31/+31
|
* halo.py: Add property "pgrid" as same as it is in slover.pyAaron LI2017-06-241-0/+10
|
* halo.py: use age at cluster redshift z0 instead of present dayAaron LI2017-06-241-6/+12
|
* halo.py: Add property "zgrid" as well as config option "zbinsize"Aaron LI2017-06-241-0/+8
|
* clusters/halo.py: calculate electron density default end at z0Aaron LI2017-06-231-10/+12
|
* halo.py: calculate magnetic field using a scaling relationAaron LI2017-06-231-8/+22
|
* halo.py + formation.py: Add argument z0Aaron LI2017-06-231-3/+7
|
* configs: extragalactic/halos: change "radius_halo" to "radius"Aaron LI2017-06-011-5/+5
|
* clusters: Move units.py and cosmology.py to utilsAaron LI2017-06-011-2/+4
|
* clusters/halo.py: Add zbegin & zend parametersAaron LI2017-06-011-3/+37
| | | | | * Add "zbegin" and "zend" parameters to "calculate_electron_spectrum()" * Add more detail docstrings