.. _parameters: Parameters and Data Files ========================= Parameters ---------- For all on/off integer flags, 0 is off and 1 is on. .. c:var:: int use_grackle Flag to activate the grackle machinery. Default: 0. .. c:var:: int with_radiative_cooling Flag to include radiative cooling and actually update the thermal energy during the chemistry solver. If off, the chemistry species will still be updated. The most common reason to set this to off is to iterate the chemistry network to an equilibrium state. Default: 1. .. c:var:: int primordial_chemistry Flag to control which primordial chemistry network is used. Default: 0. - 0: no chemistry network. Radiative cooling for primordial species is solved by interpolating from lookup tables calculated with Cloudy. A simplified set of functions are available (though not required) for use in this mode. For more information, see :ref:`tabulated-mode`. - 1: 6-species atomic H and He. Active species: H, H\ :sup:`+`, He, He\ :sup:`+`, \ :sup:`++`, e\ :sup:`-`. - 2: 9-species network including atomic species above and species for molecular hydrogen formation. This network includes formation from the H\ :sup:`-` and H\ :sub:`2`\ :sup:`+` channels, three-body formation (H+H+H and H+H+H\ :sub:`2`), H\ :sub:`2` rotational transitions, chemical heating, and collision-induced emission (optional). Active species: above + H\ :sup:`-`, H\ :sub:`2`, H\ :sub:`2`\ :sup:`+`. - 3: 12-species network include all above plus HD rotation cooling. Active species: above + D, D\ :sup:`+`, HD. .. note:: In order to make use of the non-equilibrium chemistry network (:c:data:`primordial_chemistry` options 1-3), you must add and advect baryon fields for each of the species used by that particular option. .. c:var:: int h2_on_dust Flag to enable H\ :sub:`2` formation on dust grains, dust cooling, and dust-gas heat transfer follow `Omukai (2000) `_. This assumes that the dust to gas ratio scales with the metallicity. Default: 0. .. c:var:: int metal_cooling Flag to enable metal cooling using the Cloudy tables. If enabled, the cooling table to be used must be specified with the :c:data:`grackle_data_file` parameter. Default: 0. .. note:: In order to use the metal cooling, you must add and advect a metal density field. .. c:var:: int cmb_temperature_floor Flag to enable an effective CMB temperature floor. This is implemented by subtracting the value of the cooling rate at T\ :sub:`CMB` from the total cooling rate. Default: 1. .. c:var:: int UVbackground Flag to enable a UV background. If enabled, the cooling table to be used must be specified with the :c:data:`grackle_data_file` parameter. Default: 0. .. c:var:: char* grackle_data_file Path to the data file containing the metal cooling and UV background tables. Default: "". .. c:var:: float Gamma The ratio of specific heats for an ideal gas. A direct calculation for the molecular component is used if :c:data:`primordial_chemistry` > 1. Default: 5/3. .. c:var:: int three_body_rate Flag to control which three-body H\ :sub:`2` formation rate is used. - 0: `Abel, Bryan & Norman (2002) `_ - 1: `Palla, Salpeter & Stahler (1983) `_ - 2: `Cohen & Westberg (1983) `_ - 3: `Flower & Harris (2007) `_ - 4: `Glover (2008) `_ - 5: `Forrey (2013) `_. The first five options are discussed in `Turk et. al. (2011) `_. Default: 0. .. c:var:: int cie_cooling Flag to enable H\ :sub:`2` collision-induced emission cooling from `Ripamonti & Abel (2004) `_. Default: 0. .. c:var:: int h2_optical_depth_approximation Flag to enable H\ :sub:`2` cooling attenuation from `Ripamonti & Abel (2004) `_. Default: 0. .. c:var:: int photoelectric_heating Flag to enable a spatially uniform heating term approximating photo-electric heating from dust from `Tasker & Bryan (2008) `_. Default: 0. .. c:var:: int photoelectric_heating_rate If :c:data:`photoelectric_heating` enabled, the heating rate in units of erg cm\ :sup:`-3` s\ :sup:`-1`. Default: 8.5e-26. .. c:var:: int Compton_xray_heating Flag to enable Compton heating from an X-ray background following `Madau & Efstathiou (1999) `_. Default: 0. .. c:var:: float LWbackground_intensity Intensity of a constant Lyman-Werner H\ :sub:`2` photo-dissociating radiation field in units of 10\ :sup:`-21` erg s\ :sup:`-1` cm\ :sup:`-2` Hz\ :sup:`-1` sr\ :sup:`-1`. Default: 0. .. c:var:: int LWbackground_sawtooth_suppression Flag to enable suppression of Lyman-Werner flux due to Lyman-series absorption (giving a sawtooth pattern), taken from `Haiman & Abel, & Rees (2000) `_. Default: 0. .. c:var:: int omp_nthreads Sets the number of OpenMP threads. If not set, this will be set to the maximum number of threads possible, as determined by the system or as configured by setting the ``OMP_NUM_THREADS`` environment variable. Note, Grackle must be compiled with OpenMP support enabled. See :ref:`openmp`. Data Files ---------- All data files are located in the **input** directory in the source. The first three files contain the heating and cooling rates for both primordial and metal species as well as the UV background photo-heating and photo-ionization rates. For all three files, the valid density and temperature range is given below. Extrapolation is performed when outside of the data range. The metal cooling rates are stored for solar metallicity and scaled linearly with the metallicity of the gas. Valid range: - number density: -10 < log\ :sub:`10` (n\ :sub:`H` / cm\ :sup:`-3`) < 4 - temperature: the temperature range is 1 < log\ :sub:`10` (T / K) < 9. Data files: - **CloudyData_noUVB.h5** - cooling rates for collisional ionization equilibrium. - **CloudyData_UVB=FG2011.h5** - heating and cooling rates and UV background rates from the work of `Faucher-Giguere et. al. (2009) `_, updated in 2011. The maxmimum redshift is 10.6. Above that, collisional ionization equilibrium is assumed. - **CloudyData_UVB=HM2012.h5** - heating and cooling rates and UV background rates from the work of `Haardt & Madau (2012) `_. The maximum redshift is 15.13. Above that, collisional ionization equilibrium is assumed. The final file includes only metal cooling rates under collisional ionization equilibrium, i.e., no incident radiation field. This table extends to higher densities and also varies in metallicity rather than scaling proportional to the solar value. This captures the thermalization of metal coolants occuring at high densities, making this table more appropriate for simulations of collapsing gas-clouds. Valid range: - number density: -6 < log\ :sub:`10` (n\ :sub:`H` / cm\ :sup:`-3`) < 12 - metallicity: -6 < log\ :sub:`10` (Z / Z\ :sub:`sun`) < 1 - temperature: the temperature range is 1 < log\ :sub:`10` (T / K) < 8. Data file: - **cloudy_metals_2008_3D.h5** - collisional ionization equilibrium, metal cooling rates only.