HII Regions Library
Introduction¶
Star formation in galaxies is traced by the presence of emission lines, produced by the ionized gas surrounding young star clusters. In particular, hydrogen (H) recombination lines intensity is a direct measure of the number of massive stars still alive. Thus, emission lines are important tracers of the current star formation rate (SFR). These regions are often embedded into dusty environments, and the effects of dust may be significant. Therefore, it is fundamental to model the nebular emission within a context that takes into account star-light reprocessing by dust. For this reason, during my PhD, I implemented the modelling of nebular emission by HII regions into GRASIL that was originally developed to simulate the spectral energy distribution (SED) of dusty galaxies. The implementation is based on a pre-computed library of HII region emission models.
You can find more details on this work in Panuzzo et al. 2003 and in my PhD Thesis.
Organization of the library¶
The library of HII region emission models provides the emission line luminosities for a given ionizing stellar population. It can be used in other spectrophotometric codes. You can freely download the library and use it in your code. Please acknowledge if you use this library for a publication, citing Panuzzo et al. 2003.
The library provides the emission line luminosities as a function of:
Hydrogen density nH
Filling factor e
Gas metallicity Zgas
Number of ionizing photons, per second, for hydrogen QH
The ratio QHe/QH
The ratio QO/QHe
It is composed by plain ASCII files, one for each value of hydrogen density (in number of ions per cm3) and filling factor. In the following table you can find the library’s files.
| density | filling factor | name | density | filling factor | name |
|---|---|---|---|---|---|
| 10 | 1.0 | hii_d10f10.dat.gz | 30 | 1.0 | hii_d30f10.dat.gz |
| 10 | 0.1 | hii_d10f01.dat.gz | 30 | 0.1 | hii_d30f01.dat.gz |
| 10 | 0.01 | hii_d10f001.dat.gz | 30 | 0.01 | hii_d30f001.dat.gz |
| 100 | 1.0 | hii_d100f10.dat.gz | 300 | 1.0 | hii_d300f10.dat.gz |
| 100 | 0.1 | hii_d100f01.dat.gz | 300 | 0.1 | hii_d300f01.dat.gz |
| 100 | 0.01 | hii_d100f001.dat.gz | 300 | 0.01 | hii_d300f001.dat.gz |
| 1000 | 1.0 | hii_d1000f10.dat.gz | 3000 | 0.1 | hii_d3000f01.dat.gz |
| 1000 | 0.1 | hii_d1000f01.dat.gz | 3000 | 0.01 | hii |
| 1000 | 0.01 | hii | 3000 | 0.001 | hii |
| 10000 | 0.1 | hii | |||
| 10000 | 0.01 | hii | |||
| 10000 | 0.001 | hii |
Parameter space¶
The HII region models are computed for different values of metallicity, namely: Zgas=0.0008, 0.004, 0.008, 0.015, 0.02, 0.03, 0.04, 0.05, where 0.02 means solar metallicity.
The space of Q parameters is divided into two grids, identified by a grid index. In the first grid, with index number equal 1, QHe/QH is > 0.05; in the second (grid index = 2) QHe/QH is < 0.05.
In grid number 1, QHe/QH=0.05+0.05625*i with i=0,...,8; QO/QHe=0.13+0.1*j with j=0,...,5.
In grid 2, QHe/QH=0.05, 0.025, 0.01, 0.005, 0.0025, 0.0001, 0.; while QO/QHe=0.+0.1*j with j=0,...,4.
In both the grids, QH=1046+0.2*k with k=0,...,30.
Format¶
The first line is an header; each following line correspond to an HII region model. This is an example of the first two lines of a library file (namely hii_d10f01.dat):
#density= 10.00 filling= 0.100 1)name 2)numgrid 3)zeta 4)lum 5)shi 6)sx 7)shei 8)so 9-135)intfluxline
bZ01H01E01O01 1 .0008 38.622 -2.0000 -2.4590 -3.3010 -4.1871 664505961859289591245900558898587045865964200631....Description byte by byte¶
| Characters | Format | Description |
|---|---|---|
| 1-13 | A13 | Name of the model: bytes 1-13 are a string that identify the name of the input file for CLOUDY. |
| 15 | I1 | Grid index |
| 17-21 | F5.4 | Metallicity |
| 23-28 | F6.3 | Log of the luminosity of analytical spectrum in erg/s |
| 30-37 | F8.4 | Log(QH/1050 photons/second) |
| 39-46 | F8.4 | Log(QX/1050): QX is the number of photons per second with a wavelength shorter than 750 Å |
| 48-55 | F8.4 | Log(QHe/1050) |
| 57-64 | F8.4 | Log(QO/1050) |
| 66-636 | 114I5 | Line luminosity |
Bytes from 66 to 636 correspond to the flux of 114 emission lines, an integer of 5 digits for each emission line. The line luminosity can be obtained as Luminosity=exp(intLine/1000.-50.) where intLine is the integer read from the library. I wrote a FORTRAN90 module that provide a subroutine to read these libraries.
Computed lines¶
In the following table there is the list of the computed lines. The order is the same that you can find in the library files.
Ly α 1216Å
Ly β 1025Å
Ly γ 972Å
Ly δ 949Å
Ly ε 937Å
Ly ζ 930Å
Ly η 926Å
Ly θ 922Å
H α 6563Å
H β 4861Å
H γ 4340Å
H δ 4102Å
H ε 3970Å
H ζ 3889Å
H η 3835Å
H θ 3798Å
Pa α 18752Å
Pa β 12819Å
Pa γ 10939Å
Pa δ 10050Å
Pa ε 9546Å
Pa ζ 9229Å
Pa η 9015Å
Pa θ 8863Å
Br α 40515Å
Br β 26254Å
Br γ 21657Å
Br δ 19447Å
Br ε 18175Å
Br ζ 17363Å
Br η 16808Å
Br θ 16408Å
Pf α 74585Å
Pf β 46529Å
Pf γ 37398Å
Pf δ 32964Å
Pf ε 30386Å
Pf ζ 28724Å
Pf η 27577Å
Pf θ 26746Å
Hu α 12.4µm
Hu β 75011Å
Hu γ 59071Å
Hu δ 51277Å
Hu ε 46716Å
Hu ζ 43756Å
Hu η 41700Å
Hu θ 40201Å
HeI 4472Å
HeI 5877Å
HeI 6680Å
HeI 10833Å
HeI 3889Å
HeI 7065Å
[CI]9850Å
[CI]8727Å
[CI]4621Å
[CI]610µm
[CI]369µm
[CII]157.7µm
CII]2326Å
[NI]5200Å
[NI]3467Å
[NI]10405Å
[NII]6585Å
[NII]6549Å
[NII]5756Å
[NII]122µm
[NII]205µm
NII]2141Å
[NIII]57µm
[OI]6302Å
[OI]6365Å
[OI]5578Å
[OI]63µm
[OI]145µm
[OII]3727Å
[OII]7327Å
[OII]2471Å
OIII]1663Å
[OIII]5007Å
[OIII]4960Å
[OIII]4364Å
[OIII]2321Å
[OIII]88µm
[OIII]52µm
[NeII]12.8µm
[NeIII]15.5µm
[NeIII]36µm
[NeIII]3870Å
[NeIII]3968Å
[NeIII]3343Å
[NeIII]1815Å
MgII 2800Å
[SiII]35µm
[SII]10331Å
[SII]6732Å
[SII]6717Å
[SII]4070Å
[SII]4078Å
[SIII]18.7µm
[SIII]33.5µm
[SIII]9533Å
[SIII]9071Å
[SIII]6314Å
[SIII]3723Å
[SIV]10.5µm
[ArII]7µm
[ArIII]7138Å
[ArIII]7753Å
[ArIII]5193Å
[ArIII]3110Å
[ArIII]22µm
[ArIII]9µm
Problems?!?¶
Send me an email at pasquale
- Panuzzo, P., Bressan, A., Granato, G. L., Silva, L., & Danese, L. (2003). Dust and nebular emission: I. Models for normal galaxies. Astronomy & Astrophysics, 409(1), 99–114. 10.1051/0004-6361:20031094