The H II Region/PDR Connection: Self-Consistent Calculations of Physical Conditions in Star-Forming Regions

We have performed a series of calculations designed to reproduce infrared diagnostics used to determine physical conditions in star forming regions. We self-consistently calculate the thermal and chemical structure of an H II region and photodissociation region (PDR) that are in pressure equilibrium. This differs from previous work, which used separate calculations for each gas phase. Our calculations span a wide range of stellar temperatures, gas densities, and ionization parameters. We describe improvements made to the spectral synthesis code Cloudy that made these calculations possible. These include the addition of a molecular network with ~1000 reactions involving 68 molecular species and improved treatment of the grain physics. Data from the Spitzer First Look Survey, along with other archives, are used to derive important physical characteristics of the H II region and PDR. These include stellar temperatures, electron densities, ionization parameters, UV radiation flux, and PDR density. Finally, we calculate the contribution of the H II region to PDR emission line diagnostics, which allows for a more accurate determination of physical conditions in the PDR.Comment: 60 pages, 35 figures, to be published in the Astrophysical Journal. Version with full resolution is available at http://www.pa.uky.edu/~nicholas/hii_pdr_high_res.pd

A Herschel study of Planetary Nebulae

We present Herschel PACS and SPIRE images of the dust shells around the planetary nebulae NGC 650, NGC 6853, and NGC 6720, as well as images showing the dust temperature in their shells. The latter shows a rich structure, which indicates that internal extinction in the UV is important despite the highly evolved status of the nebulae.Comment: 2 pages, 1 figure, 2012, proceedings IAU Symposium 283 Planetary Nebulae: An Eye to the Futur

Accurate determination of the free-free Gaunt factor. II - relativistic Gaunt factors

When modelling an ionised plasma, all spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum. Until now no data set exists that would meet these needs completely. We have therefore produced a table of relativistic Gaunt factors over a much wider range of parameter space than has ever been produced before. We present tables of the thermally averaged Gaunt factor covering the range log10(gamma^2) = -6 to 10 and log10(u) = -16 to 13 for all atomic numbers Z = 1 through 36. The data were calculated using the relativistic Bethe-Heitler-Elwert (BHE) approximation and were subsequently merged with accurate non-relativistic results in those parts of the parameter space where the BHE approximation is not valid. These data will be incorporated in the next major release of the spectral synthesis code Cloudy. We also produced tables of the frequency integrated Gaunt factor covering the parameter space log10(gamma^2) = -6 to 10 for all values of Z between 1 and 36. All the data presented in this paper are available online.Comment: 8 pages, 8 figures, 2 table

Radiative cooling in collisionally and photo ionized plasmas

We discuss recent improvements in the calculation of the radiative cooling in both collisionally and photo ionized plasmas. We are extending the spectral simulation code Cloudy so that as much as possible of the underlying atomic data is taken from external databases, some created by others, some developed by the Cloudy team. This paper focuses on recent changes in the treatment of many stages of ionization of iron, and discusses its extensions to other elements. The H-like and He-like ions are treated in the iso-electronic approach described previously. Fe II is a special case treated with a large model atom. Here we focus on Fe III through Fe XXIV, ions which are important contributors to the radiative cooling of hot, 1e5 to 1e7 K, plasmas and for X-ray spectroscopy. We use the Chianti atomic database to greatly expand the number of transitions in the cooling function. Chianti only includes lines that have atomic data computed by sophisticated methods. This limits the line list to lower excitation, longer wavelength, transitions. We had previously included lines from the Opacity Project database, which tends to include higher energy, shorter wavelength, transitions. These were combined with various forms of the g-bar approximation, a highly approximate method of estimating collision rates. For several iron ions the two databases are almost entirely complementary. We adopt a hybrid approach in which we use Chianti where possible, supplemented by lines from the Opacity Project for shorter wavelength transitions. The total cooling including the lightest thirty elements differs significantly from some previous calculations