BHJet
BHJet copied to clipboard
Semi analytical black hole jet model, jetset group, Amsterdam
When I the current master branch, adjust the compilation commands in `MakeBHJet`, build and then run `bhwrap.x`, I get the following output: ``` Unphysical pair content: -88.1873 pairs per proton....
Currently the radiative term assumes we are always in the Thomson regime. It may be worth using the approach here https://ui.adsabs.harvard.edu/abs/2002ApJ...568L..81D/abstract to include a KN term (depending on the photon...
https://ui.adsabs.harvard.edu/abs/2006A%26A...451..739K/abstract have a slightly better treatment to handle the transition between the two regimes, compared to what we have now
Especially for the mixed class, this would help with avoiding having a large spike when transitioning from thermal to non-thermal particles for large values of f_nth.
Have methods for both BH scale units and cgs units
Per title, we can only do low redshift blazars right now
Ideally we would like the option of treating a hollow sheath, rather than a full cylinder. That requires changing the radiative calculations. For IC it's messy because of the RT...
Include pair annihilation/creation/gamma ray absorption, Wein spectra, Bremsstrahlung
Include pair annihilation/production and gamma/gamma annihilation production
Currently Urad takes uses a dodgy pi r^2 factor to renormalize from luminosity to energy density. Instead, change to consider different geometries (and check with re-normalization to find flux at...