Commit message (Collapse) | Author | Age | Files | Lines | |
---|---|---|---|---|---|
* | clusters/halo: impose the maximum acceleration timescale | Aaron LI | 2018-01-01 | 1 | -3/+9 |
| | |||||
* | clusters/mergertree: implement itermain(), lmain, imain(idx) | Aaron LI | 2018-01-01 | 1 | -0/+45 |
| | |||||
* | clusters/mergertree: Use utils.io for pickle dump and load | Aaron LI | 2018-01-01 | 1 | -10/+4 |
| | |||||
* | clusters/mergertree: import matplotlib only when used | Aaron LI | 2018-01-01 | 1 | -4/+6 |
| | | | | In this way to avoid dependency on matplotlib. Also add print messages. | ||||
* | clusters/mergertree: improve print formats for "show_mtree()" | Aaron LI | 2018-01-01 | 1 | -4/+7 |
| | |||||
* | Minor changes | Aaron LI | 2018-01-01 | 1 | -22/+22 |
| | |||||
* | clusters/formation: Move parameter "ratio_major" to recent_major_merger() | Aaron LI | 2018-01-01 | 2 | -32/+19 |
| | | | | | | * Simplify "max_merger" property to not check whether it is a major merger. * Rename "max_merger" to "maximum_merger" * Rename "simulate_mergertree()" to "simulate_mtree()" | ||||
* | clusters/halo: Change option "tau_merger" to "time_traceback" | Aaron LI | 2018-01-01 | 1 | -4/+4 |
| | |||||
* | clusters/halo: adjust radio halo size estimation | Aaron LI | 2018-01-01 | 1 | -3/+3 |
| | |||||
* | clusters/halo: improve turbulence velocity dispersion calculation | Aaron LI | 2018-01-01 | 1 | -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" attribute | Aaron LI | 2017-12-31 | 2 | -4/+15 |
| | |||||
* | clusters/halo: Rework the turbulence acceleration timescale calculation | Aaron LI | 2017-12-31 | 1 | -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 merger | Aaron LI | 2017-12-31 | 1 | -0/+38 |
| | | | | | Also update the description of option "extragalactic/halos/eta_turb" and change its default value to 0.1 | ||||
* | clusters/halo/helper: Add fmass_nfw() to calc NFW mass profile | Aaron LI | 2017-12-31 | 1 | -0/+30 |
| | |||||
* | Minor cleanups | Aaron LI | 2017-12-31 | 1 | -3/+7 |
| | |||||
* | Minor cleanups | Aaron LI | 2017-12-30 | 1 | -4/+4 |
| | |||||
* | clusters/halo: support to drop out the most powerful halos | Aaron LI | 2017-12-30 | 1 | -0/+26 |
| | | | | | The new option "extragalactic/clusters/halo_dropout" is added to specify how many halos to be dropped out. | ||||
* | clusters/halo: Add virial_speed() helper function | Aaron LI | 2017-12-30 | 1 | -2/+13 |
| | |||||
* | clusters/halo: Estimate turbulence acceleration w.r.t. injection scale | Aaron LI | 2017-12-30 | 3 | -16/+39 |
| | | | | The turbulence injection scale ~ halo size | ||||
* | cluster/halo: Estimate halo radius ~ 0.33*R_vir (main cluster) | Aaron LI | 2017-12-30 | 2 | -8/+17 |
| | |||||
* | clusters/halo: Tune the maximum acceleration timescale to 10 Gyr | Aaron LI | 2017-11-16 | 1 | -6/+17 |
| | |||||
* | clusters/halo: improve advection calculation | Aaron LI | 2017-11-16 | 1 | -10/+10 |
| | |||||
* | Update logging messages | Aaron LI | 2017-10-29 | 1 | -0/+2 |
| | |||||
* | clusters/configs: Update f_acc default to 1.0; also update descriptions | Aaron LI | 2017-10-26 | 1 | -1/+1 |
| | |||||
* | clusters/halo: Calc turbulence injection scale w.r.t. sub-cluster radius | Aaron LI | 2017-10-26 | 1 | -6/+8 |
| | |||||
* | clusters/solver: Update buffer_np description and handling | Aaron LI | 2017-10-25 | 2 | -4/+11 |
| | |||||
* | clusters/solver: Remove an empty method | Aaron LI | 2017-10-25 | 1 | -8/+0 |
| | |||||
* | clusters/halo: Fix fp_advection() w.r.t. initial electron spectrum calc. | Aaron LI | 2017-10-25 | 1 | -2/+4 |
| | |||||
* | clusters/halo: Adjust initial electron spectrum calculation | Aaron LI | 2017-10-25 | 1 | -4/+6 |
| | |||||
* | clusters/halo: Add method to derive the initial electron spectrum | Aaron LI | 2017-10-25 | 1 | -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" object | Aaron LI | 2017-10-25 | 2 | -4/+6 |
| | |||||
* | clusters/halo: Always pass the using configs | Aaron LI | 2017-10-24 | 2 | -9/+12 |
| | | | | | 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 LI | 2017-10-24 | 1 | -19/+20 |
| | |||||
* | clusters/halo: Add method "_magnetic_field()" | Aaron LI | 2017-10-24 | 1 | -2/+23 |
| | |||||
* | clusters/halo: Rename property "magnetic_field" to "B_obs" | Aaron LI | 2017-10-24 | 2 | -4/+4 |
| | |||||
* | clusters/halo: Allow specify magnetic field in calc_emissivity() | Aaron LI | 2017-10-24 | 1 | -9/+14 |
| | |||||
* | clusters/halos: Fix halo radius calculation | Aaron LI | 2017-10-23 | 1 | -3/+5 |
| | | | | Also change the parameter "nr" default to 2.0 for draw_halo() function. | ||||
* | clusters/psformalism: Use functools.lru_cache() | Aaron LI | 2017-10-23 | 1 | -0/+3 |
| | |||||
* | clusters/halo: merge halo_rprofile() & draw_halo(); Fix halo radius error | Aaron LI | 2017-10-23 | 2 | -40/+28 |
| | | | | | | | * Fix the error that misses to divide the halo radius by 3 for the e-folding radius, which causes the generated halo image has a 3x size! * Merge function "halo_rprofile()" into "draw_halo()" * Skip transform circular halo into elliptical one if "felong" is None | ||||
* | clusters/halo: Use "lru_cache()", also improves some performance | Aaron LI | 2017-10-23 | 3 | -13/+26 |
| | |||||
* | clusters/formation: Fix problem that no mergers occurred | Aaron LI | 2017-10-22 | 2 | -2/+9 |
| | |||||
* | clusters/halo: Add "use_max_merger" option to allow use max merger events | Aaron LI | 2017-10-22 | 1 | -20/+30 |
| | |||||
* | clusters/formation: Add method to found the maximum merger event | Aaron LI | 2017-10-22 | 1 | -1/+51 |
| | |||||
* | clusters/halo: Add option "kT_out" to tune cluster temperature | Aaron LI | 2017-10-22 | 2 | -5/+7 |
| | |||||
* | clusters/halo: Rewrite turbulence acceleration & diffusion coefficient! | Aaron LI | 2017-10-21 | 2 | -58/+92 |
| | | | | | | | | | | | | 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 properties | Aaron LI | 2017-10-21 | 2 | -12/+39 |
| | | | | | | * Rvir0, Rvir_main, Rvir_sub * kT_main, kT_sub (replace kT_merger) * radius -> Rhalo, angular_radius -> Rhalo_angular, B -> B0 | ||||
* | clusters: Fix a comment about "felong" | Aaron LI | 2017-10-21 | 1 | -2/+1 |
| | |||||
* | clusters/halo: Add function "speed_sound()" | Aaron LI | 2017-10-21 | 1 | -0/+28 |
| | | | | Also add "gamma" the adiabatic index of ideal gas to utils/units.py | ||||
* | clusters/halo: Determine magnetic field by fraction of its energy density | Aaron LI | 2017-10-20 | 2 | -32/+26 |
| | | | | | | | | | | * 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: Remove "mass_to_kT()", replaced by "kT_cluster()" | Aaron LI | 2017-10-20 | 1 | -39/+0 |
| |