HII-CHI-Mistry is a collection of python subroutines aimed at the calculation of chemical abundances and physical properties using emission line fluxes from ionised gaseous nebulae. Main website of the code is available here.
HII-CHI-Mistry was written for Python v.2.7, but since its latest versions it is also compatible with Python v.3. It requires the Python library NumPy.
HII-CHI-Mistry uses grid of models (libraries) to perform the estimations. Each version of code has its own libraries, created from photoionization models computed with Cloudy v.17.
Most recent versions of the codes allow the user to introduce a different file. However, this file introduced must be located under the corresponding folder and with the proper format (see documentation for each code).
In any case, defaults libraries must not be removed.
HII-CHI-Mistry presents different variations depending on the spectral range analysed and the analysis to be performed with the input data. These are the available options:
- HII-CHI-Mistry for optical emission lines (HCm).
- HII-CHI-Mistry for ultraviolet emission lines (HCm-uv).
- HII-CHI-Mistry for infrared emission lines (HCm-ir).
- HII-CHI-Mistry for equivalent effective temperature (HCm-teff).
The latest version available is HII-CHI-Mistry v.5.22. Details on its usage and list of changes can be found here. The code is described in the following papers:
- Version for Star-Forming Galaxies (Pérez-Montero 2014).
- Version for Seyferts 2 (Pérez-Montero et al. 2019).
- Version for Low-Ionization Active Galactic Nuclei (Pérez-Díaz et al. 2021).
The latest version available is HII-CHI-Mistry-UV v.5.0. Details on its usage and list of changes can be found here. The code is described in the following papers:
- Version for Star-Forming Galaxies (Pérez-Montero & Amorín 2017).
The latest version available is HII-CHI-Mistry-IR v.3.01. Details on its usage and list of changes can be found here. The code is described in the following papers:
- Version for Star-Forming Galaxies (Fernández-Ontiveros et al. 2021).
- Version for Active Galactic Nuclei (Pérez-Díaz et al. 2022).
The latest version available is HII-CHI-Mistry-Tefff v.5.4. Details on its usage and list of changes can be found here. The code is described in the following papers:
- Version using plane-parallel or spherical geometry (Pérez-Montero et al. 2019).
- Version using density-bounded models to estimate absorbed photons (Pérez-Montero et al. 2020).
- Latest version accepting as input [NII] emission lines, as explained in Pérez-Montero et al. (2023).
Additionally, a new release of an infrared version of this code has been published, namely HII-CHI-Mistry-Teff-IR v.2.2. Details on its usage and list of changes can be found here. The code is described in the following paper:
- Introduction to the code based on infrared emission lines (Pérez-Montero et al. (2024)).
The input file must be written in text format with a first row of labels indicating some or all of the following columns:
- 'ID': identification for each row.
- 'OII_3727' and 'eOII_3727': emission line ratio [OII] 3727/Hbeta and its error.
- 'NeIII_3868' and 'eNeIII_3868': emission line ratio [NeIII] 3868/Hbeta and its error.
- 'OIII_4363' and 'eOIII_4363': emission line ratio [OIII] 4363/Hbeta and its error.
- 'OIII_4959' and 'eOIII_4959': emission line ratio [OIII] 4959/Hbeta and its error.1
- 'OIII_5007' and 'eOIII_5007': emission line ratio [OIII] 5007/Hbeta and its error.1
- 'NII_6584' and 'eNII_6584': emission line ratio [NII] 6584/Hbeta and its error.
- 'SII_6725' and 'eSII_6725': emission line ratio [SII] 6717+6731/Hbeta and its error. It is possible to give both emission lines separated as 'SII_6717' and 'SII_6731'.
1 It is possible to use only one of the two strong nebular [OIII] emission lines.
The above emission lines must be reddening corrected. In case no errors are provided, it is advisable not using MonteCarlo iterations at all. All comments must be placed at the beginning of the document preceded by the symbol "#".
The output file is also a text format, named with the original name of the input file and followed by the extension "_hcm-output.dat". The first columns will show the information provided in the input file. In case ID column is missing, the code will automatically assign a cardinal to each row. The last seven columns show:
- Grid 'i': Index of the grid employe (further details in the .readme file of the code).
- 'O/H' and 'eO/H': estimation of the oxygen abundance 12+log(O/H) and its uncertainty.
- 'N/O' and 'eN/O': estimation of the nitrogen-to-oxygen abundance log(N/O) and its uncertainty.
- 'U' and 'eU': estimation of the ionisation parameter log(U) and its uncertainty.
The input file must be written in text format with a first row of labels indicating some or all of the following columns:
- 'ID': identification for each row.
- 'Lya_1216' and 'eLya_1216': emission line Lya HI 1216 and its error.
- 'NV_1239' and 'eNV_1239': emission line NV] 1239 and its error.
- 'CIV_1549' and 'eCIV_1549': emission line CIV 1549 and its error.
- 'HeII_1640' and 'eHeII_1640': emission line HeII 1640 and its error.
- 'OIII_1665' and 'eOIII_1665': emission line OIII] 1665 and its error.
- 'CIII_1909' and 'eCIII_1909': emission line CIII 1909 and its error.
- 'Hb_4861' and 'eHb_4861': emission line Hb HI 4861 and its error.
- 'OIII_5007' and 'eOIII_5007': emission line [OIII] 5007 and its error.
The above emission lines must be reddening corrected and in the same units. In case no errors are provided, it is advisable not using MonteCarlo iterations at all. All comments must be placed at the beginning of the document preceded by the symbol "#".
The output file is also a text format, named with the original name of the input file and followed by the extension "_hcm-uv-output.dat". The first columns will show the information provided in the input file. In case ID column is missing, the code will automatically assign a cardinal to each row. The last seven columns show:
- Grid 'i': Index of the grid employe (further details in the .readme file of the code).
- 'O/H and 'eO/H': estimation of the oxygen abundance 12+log(O/H) and its uncertainty.
- 'C/O and 'eC/O': estimation of the carbon-to-oxygen abundance log(C/O) and its uncertainty.
- 'U' and 'eU': estimation of the ionisation parameter log(U) and its uncertainty.
The input file must be written in text format with a first row of labels indicating some or all of the following columns:
- 'ID': identification for each row.
- 'HI_4m' and 'eHI_4m': emission line HI 4.05 mic and its error.
- 'HI_7m' and 'eHI_7m': emission line HI 7.46 mic and its error.
- 'SIV_10m' and 'eSIV_10m': emission line [SIV] 10.5 mic and its error.
- 'HI_12m' and 'eHI_12m': emission line HI 12.4 mic and its error.
- 'NeII_12m' and 'eNII_12m': emission line [NeII] 12.8 mic and its error.
- 'NeV_14m' and 'eNeV_14m': emission line [NeV] 14.3 mic and its error.
- 'NeIII_15m' and 'eNeIII_15m': emission line [NeIII] 15.5 mic and its error.
- 'SIII_18m' and 'eSIII_18m': emission line [SIII] 18.7 mic and its error.
- 'NeV_24m' and 'eNeV_24m': emission line [NeV] 24.2 mic and its error.
- 'OIV_26m' and 'eOIV_26m': emission line [OIV] 25.9 mic and its error.
- 'SIII_33m' and 'eSIII_33m': emission line [SIII] 33.7 mic and its error.
- 'OIII_52m' and 'eOIII_52m': emission line [OIII] 52 mic and its error.
- 'NIII_57m' and 'eNIII_57m': emission line [NIII] 57 mic and its error.
- 'OIII_88m' and 'eOIII_88m': emission line [OIII] 88 mic and its error.
- 'NII_122m' and 'eNII_122m': emission line [NII] 122 mic and its error.
- 'NII_205m' and 'eNII_205m': emission line [NII] 205 mic and its error.
The above emission lines must be given in the same units. In case no errors are provided, it is advisable not using MonteCarlo iterations at all. All comments must be placed at the beginning of the document preceded by the symbol "#".
The output file is also a text format, named with the original name of the input file and followed by the extension "_hcm-ir-output.dat". The first columns will show the information provided in the input file. In case ID column is missing, the code will automatically assign a cardinal to each row. The last seven columns show:
- Grid 'i': Index of the grid employe (further details in the .readme file of the code).
- 'O/H' and 'eO/H': estimation of the oxygen abundance 12+log(O/H) and its uncertainty.
- 'N/O' and 'eN/O': estimation of the nitrogen-to-oxygen abundance log(N/O) and its uncertainty.
- 'U' and 'eU': estimation of the ionisation parameter log(U) and its uncertainty.
The input file must be written in text format with a first row of labels indicating some or all of the following columns:
- 'ID': identification for each row.
- '12logOH' and 'e12logOH': if known, oxygen abundance 12+log(O/H) and its error.1
- 'OII_3727' and 'eOII_3727': emission line [OII] 3727 and its error.
- 'OIII_4959' and 'eOIII_3727': emission line [OIII] 4959 and its error.2
- 'OIII_5007' and 'eOIII_5007': emission line [OIII] 5007 and its error.2
- 'SII_6725' and 'eSII_6725': emission line [SII] 6717+6731 and its error. It is possible to give both emission lines separated as 'SII_6717' and 'SII_6731'.
- 'SIII_9069' and 'eSIII_9069': emission line [SIII] 9069 and its error.3
- 'SIII_9532' and 'eSIII_9532': emission line [SIII] 9532 and its error.3
- 'HeI_4471' and 'eHeI_4471': emission line HeI 4471 and its error.
- 'HeI_5876' and 'eHeI_5876': emission line HeI 5876 and its error.
- 'HeII_4686' and 'eHeII_4686': emission line HeII 4686 and its error.
- 'ArIV_4730' and 'eAIV_4740': emission line [ArIV] 4740 and its error.
- 'ArIII_7135' and 'eAIII_7135': emission line [ArIII] 7135 and its error.
- 'NII_6584' and 'eNII_6584': emission line [NII] 6584 and its error.
1 If oxygen abundances are unknown, HCm will estimate them. 2 It is possible to use only one of the two strong nebular [OIII] emission lines. 3 It is possible to use only one of the two strong nebular [SIII] emission lines.
The above emission lines must be reddening corrected. In case no errors are provided, it is advisable not using MonteCarlo iterations at all. All comments must be placed at the beginning of the document preceded by the symbol "#".
The output file is also a text format, named with the original name of the input file and followed by the extension "_hcm-teff-output.dat". The first columns will show the information provided in the input file. In case ID column is missing, the code will automatically assign a cardinal to each row. The last seven columns show:
- 'O/H' and 'eO/H': estimation or input value of the oxygen abundance 12+log(O/H) and its uncertainty.
- 'Teff' and 'eTeff': estimation of the effective temperature in K and its uncertainty.4
- 'U' and 'eU': estimation of the ionisation parameter log(U) and its uncertainty.
4 If required, columns 'Teff' and 'eTeff' are replaced by 'f_abs' and 'ef_abs' (fraction of absorbed photons).
The input file must be written in text format with a first row of labels indicating some or all of the following columns:
- 'ID': identification for each row.
- '12logOH' and 'e12logOH': if known, oxygen abundance 12+log(O/H) and its error.1
- 'ArII_7m' and 'eArII_7m': [ArII] 6.98 mic and its error
- 'ArV_7m' and 'eArV_7m': [ArV] 7.90 mic and its error
- 'ArIII_9m' and 'eArIII_9m': [ArIII] 8.99 mic and its error
- 'SIV_10m' and 'eSIV_10m': [SIV] 10.5 mic and its error
- 'NeII_12m' and 'eNeII_12m': [NeII] 12.8 mic and its error
- 'ArV_13m' and 'eArV_13m': [ArV] 13.1 mic and its error
- 'NeV_14m' and 'eNeV_14m': [NeV] 14.3 mic and its error
- 'NeIII_15m' and 'eNeIII_15m': [NeIII] 15.5 mic and its error
- 'SIII_18m' and 'eSIII_18m': [SIII] 18.7 mic and its error
- 'NeV_24m' and 'eNeV_24m': [NeV] 24.2 mic and its error
- 'OIV_25m' and 'eOIV_25m': [OIV] 25.9 mic and its error
- 'SIII_33m' and 'eSIII_33m': [SIII] 33.7 mic and ist error
- 'OIII_52m' and 'eOIII_52m': [OIII] 52 mic and its error
- 'NIII_57m' and 'eNIII_57m': [NII] 57 mic and its error
- 'OIII_88m' and 'eOIII_88m': [OIII] 88 mic and its error
- 'NII_122m' and 'eNII_122m': [NII] 122 mic and its error
- 'NII_205m' and 'eNII_205m': [NII] 205 mic and its error
The above emission lines must be reddening corrected. If no information exists about a certain column it must typed as zero or not introduced in the input file. If the error is not known or if it is not going to be taken into account in the calculations, it is advisable not using MonteCarlo iterations at all. Regarding lines the routine will only provide a solution if at least one low-to-high excitation is given (e.g. [NeII] and [NeIII] and/or [SIII] and [SIV]). If only two low-excitation or high-excitation lines are given the program will provide 0 values in the results.
The output file is also a text format, named with the original name of the input file and followed by the extension "_HCm-Teff-IR-output.dat". The first columns will show the information provided in the input file. In case ID column is missing, the code will automatically assign a cardinal to each row. The last seven columns show:
- 'O/H' and 'eO/H': estimation or input value of the oxygen abundance 12+log(O/H) and its uncertainty.
- 'Teff' and 'eTeff': estimation of the effective temperature in K and its uncertainty.
- 'U' and 'eU': estimation of the ionisation parameter log(U) and its uncertainty.
This program has been made thanks to the financial support from the Spanish AYA project Estallidos.
Further questions, comments and suggestions are welcome to:
- Enrique Pérez-Montero | Pronouns: he/his | Affiliation: IAA-CSIC | Publications: List | Mail: epm[at]iaa.es
- Borja Pérez-Díaz | Pronouns: he/his | Affiliation: IAA-CSIC | Publications: List | Mail: bperez[at]iaa.es
- Juan Antonio Fernández-Ontiveros | Pronouns: he/his | Affiliation: CEFCA | Publications: List | Mail: j.a.fernandez.ontiveros[at]gmail.com