aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/extragalactic/clusters/halo.py
Commit message (Collapse)AuthorAgeFilesLines
* clusters/halo: improve turbulence velocity dispersion calculationAaron LI2018-01-011-9/+9
| | | | | | | Use the merged total mass instead of the main cluster mass only as the denominator to determine the turbulence velocity dispersion, which reduces the scatter and thus the acceleration results will be less sensitive to the (mostly unknown) merging turbulence properties.
* clusters/halo: Add back "mach_turbulence" attributeAaron LI2017-12-311-3/+14
|
* 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>