search

cosimoNigro / agnpy

Modelling jetted Active Galactic Nuclei radiative processes with python
https://agnpy.readthedocs.io/en/latest/
BSD 3-Clause "New" or "Revised" License
45 stars 32 forks source link

Mismatch of agnpy's absorption with Finke (2016) reference: DT #66

Closed cosimoNigro closed 3 years ago

cosimoNigro commented 3 years ago

As suggested by @jsitarek in issue #50, I am opening a different issue for each of the absorption crosschecks. This one regards the absorption on the photon field of the DT. I obtained the opacity vs energy computed at several distances by Finke - before I was using values I had fetched with webplotdigitizer from the paper's figures.

Here again the same problem of the BLR in issue #65 repeats: for small distances, within the disk, there is agreement, at larger distances r=10^2 R(Ly alpha) there is mismatch. On the other hand, at very large distances we have the consistency check from the approximation with the point-source behind the jet.

So I am a bit confused.

cosimoNigro commented 3 years ago

After removing the dependency on the Blob from the absorption (in PR #76), the case of the dust torus at all distances has improved. In the agnpy comparison now mu_s=1, before mu_s was taken from the blob. Here some examples, the final check compares against a point source approximating the disk in case of a very large from the BH. tau_dt_comprison_r_1e-1_R_Ly_alpha_figure_14_finke_2016 tau_dt_comprison_r_1e0_R_Ly_alpha_figure_14_finke_2016 tau_dt_comprison_r_1e1_R_Ly_alpha_figure_14_finke_2016 tau_dt_comprison_r_1e2_R_Ly_alpha_figure_14_finke_2016 tau_dt_point_source_comparison

I would close this, if it's ok with you @jsitarek. I think we have a decent agreement for this case now.

jsitarek commented 3 years ago

the agreement looks very nice, closing

jsitarek commented 3 years ago

sorry for reopening it, but I found something strange in DT absorption while investigating the BLR one. finke_dt_1 03 finke_dt_10 finke_dt_100

if r << R_dt the reference line and the agnpy line are basicaly shifted by a factor very close to 1+z. This is something that we were discussing before, but from the paper it seems that the curves are given in energy at the galaxy, and actually for BLR absorption from the same plot there is a very nice agreeement.

if the emission region is far there starts to be a disagreement between agnpy and Finke's paper. However I'm not sure of the actual parameters, table3 gives parameters for extended dust torus, while the plot is for point like. I used T=1000K (which is probably fine), and R_dt = 3.5e19 cm (which is the outer DT radius in the table for the extended DT). @cosimoNigro which parameters have you used in your comparison plots?

jsitarek commented 3 years ago

Following the discussion that we had today with @cosimoNigro I updated the plots with the proper R_dt= 1.6e19 cm (using formula 96 from the Finke'16 paper). If the Finke's points are scaled up by a factor of 2 (an issue in that paper that we discussed earlier), and are corrected by redshift of then I get a perfect agreement like Cosimo in his plots. We can close again the issue. However I think it is worth to write to Justin Finke that the BLR and DT absorption points in his Fig 14 seem to be computed w.r.t. different energies (BLR in the source frame, DT in the observer frame).

cosimoNigro commented 3 years ago

@jsitarek can we close after #89?

jsitarek commented 3 years ago

yes, closing...