Commit message (Collapse) | Author | Age | Files | Lines | ||
---|---|---|---|---|---|---|
... | ||||||
* | 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 | |
| | ||||||
* | clusters/halo: Add "kT_obs" and "tback_merger" properties | Aaron LI | 2017-10-19 | 2 | -2/+19 | |
| | ||||||
* | clusters/halo: Use "kT_cluster()" instead of "mass_to_kT()" | Aaron LI | 2017-10-19 | 2 | -4/+4 | |
| | ||||||
* | clusters/halo: Add kT_cluster() to calculate the cluster ICM temperature | Aaron LI | 2017-10-19 | 1 | -0/+25 | |
| | ||||||
* | clusters/helper: Add kT_virial() to calculate the virial temperature | Aaron LI | 2017-10-19 | 1 | -0/+36 | |
| | ||||||
* | clusters/halo: rewrite power/flux/Tb calc with k-correction considered | Aaron LI | 2017-10-19 | 3 | -130/+98 | |
| | ||||||
* | clusters/halo: Fix load halos data; fix save halos catalog | Aaron LI | 2017-10-19 | 1 | -27/+28 | |
| | | | | | The "template" item should be removed before convert the halos data into a Pandas DataFrame to be saved. | |||||
* | clusters/halo: Force to use cluster catalog if use_dump_halos_data=True | Aaron LI | 2017-10-19 | 1 | -0/+5 | |
| | ||||||
* | clusters/halo: Support use the previously dumped halos raw data | Aaron LI | 2017-10-19 | 1 | -1/+11 | |
| | | | | | | Add option "use_dump_halos_data" to control whether to just use the previously dumped halos data, therefore, the radio emissions at additional frequencies can be simply calculated. | |||||
* | cluster/halos: Split out "_calc_halos_emission()" method | Aaron LI | 2017-10-19 | 1 | -11/+32 | |
| | ||||||
* | clusters/halo: Add set_electron_spectrum() method; update doc a bit | Aaron LI | 2017-10-19 | 1 | -1/+26 | |
| | ||||||
* | clusters/halo: Add and use attribute "self.halos_data_dumpfile" | Aaron LI | 2017-10-19 | 1 | -2/+4 | |
| | ||||||
* | clusters/halo: Rewrite halos catalog and raw data save | Aaron LI | 2017-10-19 | 1 | -15/+39 | |
| | | | | | Replace option "halos_dumpfile" with "halos_catalog_outfile", and use option "dump_halos_data" to control whether to dump the raw data. | |||||
* | clusters: "catalog_outfile" is required | Aaron LI | 2017-10-19 | 1 | -3/+1 | |
| | ||||||
* | clusters/halo: Improve time step determination a bit | Aaron LI | 2017-10-19 | 2 | -3/+6 | |
| | ||||||
* | clusters/halo: Increase the time step, but require a minimum of 20 steps | Aaron LI | 2017-10-18 | 1 | -0/+7 | |
| | ||||||
* | clusters: Update radio halo radius estimation | Aaron LI | 2017-10-18 | 2 | -23/+25 | |
| | | | | | | | | | 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: More regular shape for radio halos | Aaron LI | 2017-10-16 | 1 | -6/+5 | |
| | ||||||
* | clusters: Update the documents/comments a b it | Aaron LI | 2017-10-14 | 2 | -11/+6 | |
| | ||||||
* | clusters/solver.py: Add FokkerPlanckTests with 3 cases | Aaron LI | 2017-10-12 | 1 | -9/+158 | |
| | | | | | | | | The FokkerPlanckSolver is validated with all the 3 test cases! References: * Park & Petrosian 1996, ApJS, 103, 255 * Donnert & Brunetti 2014, MNRAS, 443, 3564 | |||||
* | clusters/solver.py: significantly improve fix_boundary() | Aaron LI | 2017-10-12 | 1 | -10/+30 | |
| | | | | Also suggest that ``buffer_np`` be specified to 5%-10% of ``x_np``. | |||||
* | clusters/solver.py: Small tweaks | Aaron LI | 2017-10-07 | 1 | -8/+6 | |
| | ||||||
* | clusters/solver.py: Adjust parameter tc,uc orders | Aaron LI | 2017-10-06 | 1 | -10/+7 | |
| | ||||||
* | clusters/solver.py: Stop calc. at "tstop" by adding "tstep" parameter | Aaron LI | 2017-10-06 | 1 | -12/+19 | |
| | ||||||
* | clusters/solver.py: Update comments/references/descriptions etc. | Aaron LI | 2017-10-06 | 1 | -25/+31 | |
| | ||||||
* | clusters/psformalism.py: Rewrite sample_z_m() method | Aaron LI | 2017-10-04 | 2 | -38/+54 | |
| | | | | | Also update clusters/main.py to calculate the halo mass distributions before sampling the mass and redshifts for clusters. | |||||
* | clusters: Use "hmf" to calculate halo mass functions/distributions | Aaron LI | 2017-10-04 | 1 | -69/+136 | |
| | | | | | | | | | | | | * New dependency "hmf" (halo mass functions) module * Calculate halo mass distributions/functions (dndlnm) with respect to masses and redshifts, instead of use the previous data file ("ps_data") * New section "[extragalactic][psformalism]" in configurations * New functions to write and read the dndlnm data TODO: * update the method to sample (mass, redshift) for clusters from the dndlnm data | |||||
* | clusters/psformalism.py: Improve catalog comments | Aaron LI | 2017-10-02 | 1 | -4/+7 | |
| | ||||||
* | clusters/emission: Assume electron pitch angle be pi/2 | Aaron LI | 2017-09-20 | 1 | -23/+17 | |
| | | | | And some cleanup and small changes | |||||
* | clusters/emission: Rewrite synchrotron kernel function | Aaron LI | 2017-09-20 | 1 | -9/+38 | |
| | | | | | | | | | | | Add asymptotic functions to calculate the values beyond the interpolation bounds (e.g., <1e-3 and >10), otherwise, the calculated synchrotron emissivity is overestimated at the higher frequencies. By rewrite this synchrotron kernel function, the calculated results is consistent with the theoretical/analytical results, e.g., the synchrotron radiation of a population of electrons of power-law index n is also a power-law with index (n-1)/2. | |||||
* | clusters: Add option "use_output_catalog" to continue halo simulation | Aaron LI | 2017-09-03 | 1 | -4/+19 | |
| | ||||||
* | clusters: Add config option "boost" | Aaron LI | 2017-09-01 | 1 | -2/+5 | |
| | | | | | | | NOTE: The cluster number (especially the bright ones) may be rather small within a sky patch. This ``boost`` config option increase the expected cluster number by a specified times, for better testing purpose. | |||||
* | Update "Fnu_to_Tb_fast()" to "Fnu_to_Tb()" | Aaron LI | 2017-09-01 | 6 | -9/+9 | |
| | ||||||
* | clusters/psformalism: Remove config "f_darkmatter" | Aaron LI | 2017-08-29 | 1 | -3/+2 | |
| | | | | Use "1 - baryon_fraction = 1 - Ob0/Om0" to replace "f_darkmatter". | |||||
* | clusters/halo.py: Update injection rate calculation | Aaron LI | 2017-08-26 | 1 | -10/+5 | |
| | | | | | The electrons are assumed to be injected throughout the whole cluster ICM/volume. | |||||
* | Return the list of filepath to the written sky files | Aaron LI | 2017-08-26 | 1 | -0/+8 | |
| | ||||||
* | clusters/main.py: Update sky header; update comments etc. | Aaron LI | 2017-08-26 | 1 | -5/+7 | |
| | ||||||
* | configs: Remove option "save"; all enabled components will be saved | Aaron LI | 2017-08-26 | 1 | -1/+0 | |
| |