Data from Dannen et al. (2019)
Photoionization calculations of the radiation force due to spectral lines in AGNs
(Accepted to ApJ)
arXiv posting
XSTAR CALCULATIONS
Here we make publicly available the results of our photoionization calculations for force multipliers
and associated heating/cooling rates for the Type 1 and 2 AGN SEDs from
Mehdipour et al. (2015).
A tarball containing all of the data is available for download
here.
These tables provide the essential data needed to perform many types of state-of-the-art radiation-hydrodynamical calcuations.
For example, these tables now make it possible to model line-driven winds with self-consistent heating and
cooling rates to also account for thermal driving - a task that has yet to be performed!
Alternatively, these tables can be used together to self-consistently model
cloud formation and acceleration in AGN using the methods developed by
Proga & Waters (2015).
If used in combination with MHD simulations, these force multipliers and heating/cooling rates can be used
to compute the first ever self-consistent models of magnetothermal line-driven winds or
cloud formation and acceleration in the presence of magnetic fields.
As a final example, they can be used to compute models of clumpy accretion flows following the methods of
Moscibrodzka & Proga (2013),
as both of these AGN equilibrium curves (see below) feature large, thermally unstable regions.
FORCE MULTIPLIER TABLES
AGN1_Fmult.dat
AGN2_Fmult.dat
These tables are formatted to be easy to parse (see the python example below), and a C++ interface to simulation codes is provided.
The first entry in each table is N_xi, the number of photoionization parameter values. The remainder of the first row contains the N_xi values of log10(xi). The remainder of the first column is all the log10 values of the optical depth parameter, t.
The entries corresponding to a given (t,xi) pair are the values of log10(M), where M is the force multiplier.
HEATING/COOLING TABLES
AGN1_heating.dat
AGN1_cooling.dat
AGN1_net.dat
AGN2_heating.dat
AGN2_cooling.dat
AGN2_net.dat
For each SED (AGN1 or AGN2), there are three tables with labels 'heating', 'cooling', and 'net' (heating - cooling).
The first entry in each table is N_xi, the number of photoionization parameter values. The remainder of the first row contains the N_xi values of xi (in cgs units this time). The remainder of the first column is all of the temperature values (in Kelvin).
The entries corresponding to a given (T,xi) pair are the rates in units of erg cm^3 s^-1.
Radiative equilibrium curves
AGN1_eq_curve.dat
AGN2_eq_curve.dat
These are the equlibrium curves (where heating balances cooling) corresponding to the Type 1/2 AGN SEDs.
These curves can be extracted from the tables above, so we provide them separately as a means to verify the table parser.
SUPPORTIVE CODE
Simple python reader
Here
is a sample script that parses the tables and generates a contour plot of the equlibirum curves and cooling rates.
Hydro interface/interpolation routine
VegasTables is a sophisticated routine written in C++ that is intended for use with hydrodynamical codes.
This code first parses a table and then performs either bilinear or bicubic interpolation to evaluate the rates
at any given (xi,T) or the force multipliers at any given (xi,t).
We employed an earlier version of this interface (which was applied to
Athena++)
in the hydrodynamical calculations published by
Dyda et al. (2017),
where the bilinear interpolation method is documented (see the Appendix).
The routines in vegas_tables.h now perform bicubic interpolation using the GSL library, as bicubic
interpolation is found to be about 5x more accurate than bilinear interpolation.
See the documentation within vegas_tables.h for example usage.
This is a general routine that works for any ordered tabular data. The data being 'ordered' is important,
as it permits immediate indexing of the entry locations, obviating the need for a table search.
Citing these tables
Please cite Dannen et al. (2019) if you publish results that utilized any of these tables.
Acknowledgements
These calculations were performed on the Cherry Creek cluster at UNLV.
Support for Program number HST-AR-14579.001-A was provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555.
Work at LANL was done under the auspices of the National
Nuclear Security Administration of the US Department
of Energy at Los Alamos National Laboratory under
Contract No. DE-AC52-06NA25396. TW is partially supported by
LANL LDRD Exploratory Research Grant 20170317ER.
(Accepted to ApJ)
arXiv posting