Arxiv:1611.08305V1 [Astro-Ph.GA] 24 Nov 2016
Total Page:16
File Type:pdf, Size:1020Kb
Draft version November 28, 2016 Typeset using LATEX twocolumn style in AASTeX61 NEBULAR CONTINUUM AND LINE EMISSION IN STELLAR POPULATION SYNTHESIS MODELS Nell Byler,1 Julianne J. Dalcanton,1 Charlie Conroy,2 and Benjamin D. Johnson2 1Department of Astronomy, University of Washington, Box 351580, Seattle, WA 98195, USA 2Department of Astronomy, Harvard University, Cambridge, MA, USA (Received November 16, 2016) Submitted to ApJ ABSTRACT Accounting for nebular emission when modeling galaxy spectral energy distributions (SEDs) is important, as both line and continuum emission can contribute significantly to the total observed flux. In this work, we present a new nebular emission model integrated within the Flexible Stellar Population Synthesis code that computes the total line and continuum emission for complex stellar populations using the photoionization code Cloudy. The self-consistent coupling of the nebular emission to the matched ionizing spectrum produces emission line intensities that correctly scale with the stellar population as a function of age and metallicity. This more complete model of galaxy SEDs will improve estimates of global gas properties derived with diagnostic diagrams, star formation rates based on Hα, and stellar masses derived from NIR broadband photometry. Our models agree well with results from other photoionization models and are able to reproduce observed emission from H ii regions and star-forming galaxies. Our models show improved agreement with the observed H ii regions in the Ne iii/O ii plane and show satisfactory agreement with He ii emission from z = 2 galaxies when including rotating stellar models. Models including post-asymptotic giant branch stars are able to reproduce line ratios consistent with low-ionization emission regions (LIERs). Keywords: Galaxies | galaxies: abundances | galaxies: ISM | galaxies: star formation arXiv:1611.08305v1 [astro-ph.GA] 24 Nov 2016 Corresponding author: Nell Byler [email protected] 2 Byler et al. 1. INTRODUCTION For Starburst-99, line emission is computed from a The light emerging from galaxies is a complex com- normalized library of stellar UV spectra, with absolute bination of emission from stars and gas, processed by fluxes derived from the stellar SED. While these analytic any intervening dust. To fully model this emission, it prescriptions are computationally efficient, they cannot is necessary to include the effects of each of these com- account for temporal or chemical evolution. ponents. For star-forming galaxies, the UV and optical For a more detailed analysis of the stellar and nebular flux is dominated by the light produced by young, lu- energy distributions, population synthesis models can minous stars and the surrounding networks of ionized be coupled with photoionization models. Photoioniza- gas. The latter produces nebular emission, which can tion models have proven to be essential in interpreting contribute as much as 20-60% of broadband fluxes and the emission-line properties of H ii regions in terms of which is responsible for nearly all the optical emission the properties of the stars and gas (e.g., Dopita et al. lines present in the spectra of star-forming galaxies (An- 2000) and the nebular emission from galaxies in terms ders & Fritze-v. Alvensleben 2003; Reines et al. 2010). of macroscopic star formation parameters (e.g., Brinch- Emission lines are routinely used at low and high red- mann et al. 2004). Popular photoionization codes in- shift to measure key physical properties of entire galax- clude Cloudy (Ferland et al. 2013) and Mappings- ies, such as the star formation rate (SFR) and metallicity III (Groves et al. 2004), both of which compute the (e.g., Tremonti et al. 2004; Kewley & Ellison 2008). full radiative transfer through a gas cloud and predict Nebular emission is comprised of two components: (1) the resultant emission spectrum. Most implementations the nebular continuum, which is a continuous emission model the total nebular emission as the sum of emission spectrum that consists of free-free (bremsstrahlung), from multiple H ii regions with different ages and phys- free-bound (recombination continuum), and two-photon ical properties, though at the cost of being much more emission; and (2) nebular line emission, which is primar- computationally expensive(e.g., Kewley et al. 2001; Moy ily produced by radiative recombination processes and et al. 2001; Charlot & Longhetti 2001; Dopita et al. emission from forbidden and fine structure line transi- 2006). We present a nebular emission model within the Flex- tions. The strength of emission from these two compo- 1 nents depends on both the ionizing radiation field and ible Stellar Population Synthesis code (FSPS , Conroy the metallicity of the gas. The amount of nebular emis- et al. 2009) that adopts the best features of the realis- sion thus varies from galaxy to galaxy, and can evolve tic photoionization models and implements them within with cosmic time. a more flexible stellar population framework. We cou- Stellar Population Synthesis (SPS) models used to in- ple the ionizing spectrum and stellar metallicity with terpret galaxy observations account for the light emit- the gas phase metallicity to self-consistently compute ted by stars and the reprocessing of that light by dust, the total line and continuum emission. Following the but only a handful of current SPS codes include nebu- process of Charlot & Longhetti(2001), we use simple lar emission, despite its important effects on the output stellar populations (SSPs) as the ionizing source for the spectrum (see reviews in Walcher et al. 2011; Conroy gas clouds using the photoionization code Cloudy. The 2013). The contribution from nebular emission can be resultant nebular line and continuum emission are em- calculated with varying levels of sophistication. The bedded within FSPS as pre-computed tables. We can simplest approach computes the line and continuum then combine the results from the SSPs to produce self- emission analytically, using the number of photons in consistent spectra for arbitrary star formation histories the Lyman continuum to calculate the strength of emis- (SFHs). This strategy maintains the flexibility of FSPS sion as a function of wavelength. The most complex without the need to rerun the computationally expensive and accurate approach uses photoionization models to photoionization models for each output stellar popula- compute the transfer of ionizing photons exactly. tion. The SPS codes Pegase´ (Fioc & Rocca-Volmerange The self-consistent implementation of the nebular 1999), PopStar (Moll´aet al. 2009), and Starburst- model is two-fold: First, the nebular line and contin- 99 (Leitherer et al. 1999) all include analytic prescrip- uum emission predicted by Cloudy is added to the tions for nebular emission, where it is assumed that all same spectrum that was used to produce the emission stellar photons with energies greater than 13.6 eV are lines, and is thus directly tied to the ionizing contin- converted into nebular emission. The nebular contin- uum of each SSP as a function of age, metallicity, and uum contribution is calculated based on emission coef- ionization parameter. Second, by linking the stellar ficients for free-free, bound-free, and two-photon transi- metallicity in FSPS with the gas-phase abundances in tions. Hydrogenic line intensities are similarly computed Cloudy we couple the metallicity dependent changes by translating the number of photons in the Lyman con- in the ionizing EUV spectral shape to the changes in tinuum directly into a line intensity based on Case-B ap- the gas coolants, which is reflected in both the temper- proximations. For Pegase´ , the line intensities for other elements are included based on empirical line strengths. 1 available on GitHub https://github.com/cconroy20/fsps 3 ature and ionization structure of the nebula. This is a constant throughout the gas) to those run with an iso- particularly important feature, as it allows us to more baric density law (which dictates that the gas pressure accurately model the simultaneous signature of stars, follows Pgas = ntot · k · Te) and found that our results gas, and dust in the integrated spectra of galaxies, as did not change by more than a few percent. Differences discussed below. between Mappings-III and Cloudy models are there- In x 2 we introduce the nebular emission model and fore unlikely to be due solely to the different application our means of coupling FSPS and Cloudy. We first dis- of density and pressure laws. cuss broad trends in the ionizing spectra in x 3, before moving on to discuss the results of the nebular model 2.1.2. Gas Chemical Content in x 4. We begin x 4 by discussing the ionization struc- We adopt the gas phase abundances specified by Do- ture of individual H ii regions as a function of age and pita et al.(2000), which are based on the solar abun- metallicity (x 4.1). We then use these results to discuss dances from Anders & Grevesse(1989). We assume that the origin of common optical and NIR lines and their the gas phase metallicity scales with the metallicity of variation in line strength (x 4.2). In x 5 we generate the stellar population, given that the metallicity of the model grids of line ratios and showcase their ability to newly-formed stars should be identical to the metallic- reproduce lines observed in H ii regions and star-forming ity of the gas cloud from which the stars formed. Metal galaxies. We consider several complicating features in abundances are solar-scaled, with the exception of nitro- x 6, including evolutionary tracks and dust, followed by gen, which has a known secondary nucleosynthetic con- our conclusions in x 7. tribution. For the scaling of nitrogen with metallicity, we follow the piecewise relationship between nitrogen 2. METHODS AND IMPLEMENTATION and oxygen specified by Dopita et al.(2000). In the following sections we discuss the parametriza- Elements like carbon and nitrogen are observed to be tion of the nebular emission model and our means of heavily depleted onto dust grains in H ii regions.