aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/extragalactic/clusters/halo.py
Commit message (Collapse)AuthorAgeFilesLines
* clusters/halo: Rework the turbulence acceleration timescale calculationAaron LI2017-12-311-47/+39
| | | | | | | | The new method adopted to determine the turbulence acceleration timescale follows [miniati2015] and [pinzke2017], in which the merger-induced turbulence cascades into small scales (and thus accelerate electrons) through the TTD (transit-time damping) mechanism by scattering mainly with the relativistic particles (cosmic-ray protons and electrons).
* clusters/halo: calc turbulence velocity dispersion based on mergerAaron LI2017-12-311-0/+38
| | | | | Also update the description of option "extragalactic/halos/eta_turb" and change its default value to 0.1
* Minor cleanupsAaron LI2017-12-311-3/+7
|
* clusters/halo: Estimate turbulence acceleration w.r.t. injection scaleAaron LI2017-12-301-15/+15
| | | | The turbulence injection scale ~ halo size
* cluster/halo: Estimate halo radius ~ 0.33*R_vir (main cluster)Aaron LI2017-12-301-1/+1
|
* clusters/halo: Tune the maximum acceleration timescale to 10 GyrAaron LI2017-11-161-6/+17
|
* clusters/halo: improve advection calculationAaron LI2017-11-161-10/+10
|
* clusters/configs: Update f_acc default to 1.0; also update descriptionsAaron LI2017-10-261-1/+1
|
* clusters/halo: Calc turbulence injection scale w.r.t. sub-cluster radiusAaron LI2017-10-261-6/+8
|
* clusters/solver: Update buffer_np description and handlingAaron LI2017-10-251-0/+2
|
* clusters/halo: Fix fp_advection() w.r.t. initial electron spectrum calc.Aaron LI2017-10-251-2/+4
|
* clusters/halo: Adjust initial electron spectrum calculationAaron LI2017-10-251-4/+6
|
* clusters/halo: Add method to derive the initial electron spectrumAaron LI2017-10-251-11/+51
| | | | | Also add the option "time_init" to control how long a period is used to derive the initial electron spectrum.
* clusters/halo: Fix passing the "configs" objectAaron LI2017-10-251-3/+5
|
* clusters/halo: Always pass the using configsAaron LI2017-10-241-3/+6
| | | | | The user-using `configs` object may be different to the global `CONFIGS`, which may be forgotten to load the up-to-date user configurations.
* clusters/halo: Change zbegin/zend to tstart/tsop (use cosmic time)Aaron LI2017-10-241-19/+20
|
* clusters/halo: Add method "_magnetic_field()"Aaron LI2017-10-241-2/+23
|
* clusters/halo: Rename property "magnetic_field" to "B_obs"Aaron LI2017-10-241-3/+3
|
* clusters/halo: Allow specify magnetic field in calc_emissivity()Aaron LI2017-10-241-9/+14
|
* clusters/halo: Use "lru_cache()", also improves some performanceAaron LI2017-10-231-12/+22
|
* clusters/halo: Add option "kT_out" to tune cluster temperatureAaron LI2017-10-221-3/+1
|
* 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>