The gas that is present in the interstellar medium is usually very far removed from (local) thermodynamic equilibrium, and in some cases may also not be in a steady-state equilibrium with its surroundings. The physics of this material is complex and one needs a sophisticated numerical code to study it. For this purpose the open-source photoionization code Cloudy was created. It models the physical state of the gas and predicts the spectrum that it emits. Cloudy is continually being developed to improve the treatment of the microphysical processes and the database of fundamental data that it uses. In this paper we will discuss how we are developing the code to improve our high-density predictions by implementing better collisional-radiative models for all ions. We will also briefly discuss the experimental mode in Cloudy to model gas that is not in steady-state equilibrium and present a preliminary model of recombining gas in a planetary nebula that is on the cooling track. We finish with a short discussion of how we are speeding up the code by using parallelization.

The 180-day Space Telescope and Optical Reverberation Mapping campaign on NGC 5548 discovered an anomalous period, the broad-line region (BLR) holiday, in which the emission lines decorrelated from the continuum variations. This is important since the correlation between the continuum-flux variations and the emission-line response is the basic assumption for black hole (BH) mass determinations through reverberation mapping. During the BLR holiday, the high-ionization intrinsic absorption lines also decorrelated from the continuum as a result of variable covering factor of the line of sight (LOS) obscurer. The emission lines are not confined to the LOS, so this does not explain the BLR holiday. If the LOS obscurer is a disk wind, its streamlines must extend down to the plane of the disk and the base of the wind would lie between the BH and the BLR, forming an equatorial obscurer. This obscurer can be transparent to ionizing radiation, or can be translucent, blocking only parts of the SED, depending on its density. An emission-line holiday is produced if the wind density increases only slightly above its transparent state. Both obscurers are parts of the same wind, so they can have associated behavior in a way that explains both holidays. A very dense wind would block nearly all ionizing radiation, producing a Seyfert 2 and possibly providing a contributor to the changing-look AGN phenomenon. Disk winds are very common and we propose that the equatorial obscurers are too, but mostly in a transparent state.

Fundamental atomic transition parameters, such as oscillator strengths and wavelengths, play a key role in modelling and understanding the chemical composition of stars in the universe. Despite the significant work under way to produce these parameters for many ions, uncertainties in these parameters remain large and can limit the accuracy of chemical abundance determinations. The Belgian repository of fundamental atomic data and stellar spectra (BRASS) aims to provide a large systematic and homogeneous quality assessment of the atomic data available for quantitative spectroscopy. BRASS shall compare synthetic spectra against extremely high quality observed spectra, at a resolution of ~85000 and signal-noise ratios of ~1000, for around 20 bright BAFGK spectral type stars, in order to evaluate the atomic data available for over a thousand potentially useful spectral lines. A large-scale homogeneous selection of atomic lines is performed by synthesising theoretical spectra of literature atomic lines, for FGK-type stars including the Sun, resulting in a selection of 1091 theoretically deep and unblended lines, in the wavelength range 4200-6800~\AA, which may be suitable for quality assessment. Astrophysical log(gf) values are determined for the 1091 transitions using two commonly employed methods. The agreement of these log(gf) values are used to select well-behaving lines for quality assessment. 845 atomic lines were found to be suitable for quality assessment, of which 408 were found to be robust against systematic differences between analysis methods. Around 54% of the quality-assessed lines were found to have at least one literature log(gf) value in agreement with our derived values, though the remaining values can disagree by as much as 0.5 dex. Only 38% of FeI lines were found to have sufficiently accurate log(gf) values, increasing to ~70-75% for the remaining Fe-group lines.

The flux variations in the emission lines in active galactic nuclei (AGNs) are driven by variations in the ionizing continuum flux --which are usually reflected in the observable UV-optical continuum. The "Reverberation mapping" technique measures the delay between line and continuum variations to determine the size of the line emitting region, this is the basis for measurements of the central black hole mass in AGNs. The Space Telescope and Optical Reverberation Mapping Project (AGN STORM) on NGC 5548 in 2014 is the most intensive multi-wavelength AGN monitoring campaign ever. For most of the campaign, the emission-line variations followed changes in the continuum with a time lag, as expected. However, the lines varied independently of the observed UV-optical continuum during a 60 -- 70 day "holiday." To understand this remarkable phenomenon, we study the intrinsic absorption lines present in NGC 5548. We identify a novel cycle that reproduces the absorption line variability and thus identify the physics that allows the holiday to occur. In our model, variations in this obscurer's line-of-sight covering factor modify the soft X-ray continuum. This leads to changes in the ionization of helium gas in the broad-line region. Ionizing radiation produced by recombining helium then affects the ionization of other species as observed during the AGN STORM holiday. It is likely that any other model which selectively changes the soft X-ray part of the continuum during the holiday can also explain the anomalous emission line behavior observed.

At intermediate to high densities, electron (de-)excitation collisions are the dominant process for populating or depopulating high Rydberg states. In particular, the accurate knowledge of the energy changing ($n$-changing) collisional rates is determinant for predicting the radio recombination spectra of gaseous nebula. The different datasets present in the literature come either from impact parameter calculations or semi-empirical fits and the rate coefficients agree within a factor of two. We show in this paper that these uncertainties cause errors lower than 5% in the emission of radio recombination lines (RRL) of most ionized plasmas of typical nebulae. However, in special circumstances where the transitions between Rydberg levels are amplified by maser effects, the errors can increase up to 20%. We present simulations of the optical depth and H$n\alpha$ line emission of Active Galactic Nuclei (AGN) Broad Line Regions (BLRs) and the Orion Nebula Blister to showcase our findings.

Radio continuum observations trace thermal emission of ionized plasma in planetary nebulae and bring useful information on nebular geometries. A model of homogeneous sphere or shell cannot fit the nebular spectra and brightness temperatures. Two alternative models have been proposed in the literature: the first one consists of two homogeneous components, while the other one is a model of a shell with a significant radial density gradient. On the other side, prolate ellipsoidal shell model can successfully fit the surface brightness distribution of selected objects. We verify the existing models using data collected in radio surveys covering wide range of frequencies. In about 50% cases, density gradient can be excluded, and none of the remaining objects could be confirmed. None of the observed planetary nebulae show the spectral index of 0.6 in the optically thick part of the spectrum, which is a value predicted for a shell containing strong radial density gradient. Radio spectra can be fitted with a model of prolate ellipsoidal shell, but also by a shell containing temperature variations in planetary nebulae. At least eight planetary nebulae show two component spectra, with one compact component showing much higher optical thickness than the other one. Unexpectedly, a group of planetary nebulae with lowest surface brightness show non-negligible optical thickness. Their emission comes from compact and dense structures, comprising only small part of total nebular mass.

The bipolar nebula Menzel 3 (Mz 3) was observed as part of the \textit{Herschel} Planetary Nebula Survey (\textit{HerPlaNS}), which used the PACS and SPIRE instruments aboard the \textit{Herschel Space Observatory} to study a sample of planetary nebulae (PNe). In this paper, one of the series describing \textit{HerPlaNS} results, we report the detection of H I recombination lines (HRLs) in the spectrum of Mz 3. Inspection of the spectrum reveals the presence of 12 HRLs in the 55 to 680 $\mu$m range covered by the PACS and SPIRE instruments (H11$\alpha$ to H21$\alpha$ and H14$\beta$). The presence of HRLs in this range is unusual for PNe and has not been reported in Mz 3 before. Our analysis indicates that the HRLs we observed are enhanced by laser effect occurring in the core of Mz 3. Our arguments for this are: (i) the available Mz 3 optical to submillimetre HRL $\alpha$ line intensity ratios are not well reproduced by the spontaneous emission of optically thin ionized gas, as would be typical for nebular gas in PNe; (ii) the compact core of Mz 3 is responsible for a large fraction of the Herschel HRLs emission; (iii) the line intensity ratios for Mz 3 are very similar to those in the core emission of the well known star MWC 349A, where laser effect is responsible for the enhancement of HRLs in the Herschel wavelength range; (iv) the physical characteristics relevant to cause laser effect in the core of MWC 349A are very similar to those in the core of Mz 3.

We perform a comprehensive analysis of the planetary nebula (PN) NGC6781 to investigate the physical conditions of each of its ionized, atomic, and molecular gas and dust components and the object's evolution, based on panchromatic observational data ranging from UV to radio. Empirical nebular elemental abundances, compared with theoretical predictions via nucleosynthesis models of asymptotic giant branch (AGB) stars, indicate that the progenitor is a solar-metallicity, 2.25-3.0 Msun initial-mass star. We derive the best-fit distance of 0.46 kpc by fitting the stellar luminosity (as a function of the distance and effective temperature of the central star) with the adopted post-AGB evolutionary tracks. Our excitation energy diagram analysis indicate high excitation temperatures in the photodissociation region (PDR) beyond the ionized part of the nebula, suggesting extra heating by shock interactions between the slow AGB wind and the fast PN wind. Through iterative fitting using the Cloudy code with empirically-derived constraints, we find the best-fit dusty photoionization model of the object that would inclusively reproduce all of the adopted panchromatic observational data. The estimated total gas mass (0.41 Msun) corresponds to the mass ejected during the last AGB thermal pulse event predicted for a 2.5 Msun initial-mass star. A significant fraction of the total mass (about 70 percent) is found to exist in the PDR, demonstrating the critical importance of the PDR in PNe that are generally recognized as the hallmark of ionized/H+ regions.

We describe the 2017 release of the spectral synthesis code Cloudy. A major development since the previous release has been exporting the atomic data into external data files. This greatly simplifies updates and maintenance of the data. Many large datasets have been incorporated with the result that we can now predict well over an order of magnitude more emission lines when all databases are fully used. The use of such large datasets is not realistic for most calculations due to the time and memory needs, and we describe the limited subset of data we use by default. Despite the fact that we now predict significantly more lines than the previous Cloudy release, this version is faster because of optimization of memory access patterns and other tuning. The size and use of the databases can easily be adjusted in the command-line interface. We give examples of the accuracy limits using small models, and the performance requirements of large complete models. We summarize several advances in the H- and He-like iso-electronic sequences. We use our complete collisional-radiative models of the ionization of these one and two-electron ions to establish the highest density for which the coronal or interstellar medium (ISM) approximation works, and the lowest density where Saha or local thermodynamic equilibrium can be assumed. The coronal approximation fails at surprisingly low densities for collisional ionization equilibrium but is valid to higher densities for photoionized gas clouds. Many other improvements to the physics have been made and are described. These include the treatment of isotropic continuum sources such as the cosmic microwave background (CMB) in the reported output, and the ability to follow the evolution of cooling non-equilibrium clouds.

We compare the results of the semi-classical (SC) and quantum-mechanical (QM) formalisms for angular-momentum changing transitions in Rydberg atom collisions given by Vrinceanu & Flannery, J. Phys. B 34, L1 (2001), and Vrinceanu, Onofrio & Sadeghpour, ApJ 747, 56 (2012), with those of the SC formalism using a modified Monte Carlo realization. We find that this revised SC formalism agrees well with the QM results. This provides further evidence that the rates derived from the QM treatment are appropriate to be used when modelling recombination through Rydberg cascades, an important process in understanding the state of material in the early universe. The rates for $\Delta\ell=\pm1$ derived from the QM formalism diverge when integrated to sufficiently large impact parameter, $b$. Further to the empirical limits to the $b$ integration suggested by Pengelly & Seaton, MNRAS 127, 165 (1964), we suggest that the fundamental issue causing this divergence in the theory is that it does not fully cater for the finite time taken for such distant collisions to complete.

Accurate rates for energy-degenerate l-changing collisions are needed to determine cosmological abundances and recombination. There are now several competing theories for the treatment of this process, and it is not possible to test these experimentally. We show that the H I two-photon continuum produced by astrophysical nebulae is strongly affected by l-changing collisions. We perform an analysis of the different underlying atomic processes and simulate the recombination and two-photon spectrum of a nebula containing H and He. We provide an extended set of effective recombination coefficients and updated l-changing 2s-2p transition rates using several competing theories. In principle, accurate astronomical observations could determine which theory is correct.

V4334 Sgr (a.k.a. Sakurai's object) is the central star of an old planetary nebula that underwent a very late thermal pulse a few years before its discovery in 1996. We have been monitoring the evolution of the optical emission line spectrum since 2001. The goal is to improve the evolutionary models by constraining them with the temporal evolution of the central star temperature. In addition the high resolution spectral observations obtained by X-shooter and ALMA show the temporal evolution of the different morphological components.

Cosmological models can be constrained by determining primordial abundances. Accurate predictions of the He I spectrum are needed to determine the primordial helium abundance to a precision of $< 1$% in order to constrain Big Bang Nucleosynthesis models. Theoretical line emissivities at least this accurate are needed if this precision is to be achieved. In the first paper of this series, which focused on H I, we showed that differences in $l$-changing collisional rate coefficients predicted by three different theories can translate into 10% changes in predictions for H I spectra. Here we consider the more complicated case of He atoms, where low-$l$ subshells are not energy degenerate. A criterion for deciding when the energy separation between $l$ subshells is small enough to apply energy-degenerate collisional theories is given. Moreover, for certain conditions, the Bethe approximation originally proposed by Pengelly & Seaton (1964) is not sufficiently accurate. We introduce a simple modification of this theory which leads to rate coefficients which agree well with those obtained from pure quantal calculations using the approach of Vrinceanu et al. (2012). We show that the $l$-changing rate coefficients from the different theoretical approaches lead to differences of $\sim 10$% in He I emissivities in simulations of H II regions using spectral code Cloudy.

Hydrogen and helium emission lines in nebulae form by radiative recombination. This is a simple process which, in principle, can be described to very high precision. Ratios of He I and H I emission lines can be used to measure the He$^+$/H$^+$ abundance ratio to the same precision as the recombination rate coefficients. This paper investigates the controversy over the correct theory to describe dipole $l$-changing collisions ($nl\rightarrow nl'=l\pm 1$) between energy-degenerate states within an $n$-shell. The work of Pengelly & Seaton (1964) has, for half-a-century, been considered the definitive study which "solved" the problem. Recent work by Vrinceanu et al.(2012) recommended the use of rate coefficients from a semi-classical approximation which are nearly an order of magnitude smaller than those of Pengelly & Seaton (1964), with the result that significantly higher densities are needed for the $nl$ populations to come into local thermodynamic equilibrium. Here, we compare predicted H~I emissivities from the two works and find widespread differences, of up to $\approx 10$%. This far exceeds the 1% precision required to obtain the primordial He/H abundance ratio from observations so as to constrain Big Bang cosmologies. We recommend using the rate coefficients of Pengelly & Seaton (1964) for $l$-changing collisions, to describe the H recombination spectrum, based-on their quantum mechanical representation of the long-range dipole interaction.

Recent atomic physics calculations for Si II are employed within the Cloudy modelling code to analyse Hubble Space Telescope (HST) STIS ultraviolet spectra of three cool stars, Beta-Geminorum, Alpha-Centauri A and B, as well as previously published HST/GHRS observations of Alpha-Tau, plus solar quiet Sun data from the High Resolution Telescope and Spectrograph. Discrepancies found previously between theory and observation for line intensity ratios involving the 3s$^{2}$3p $^{2}$P$_{J}$--3s3p$^{2}$ $^{4}$P$_{J^{\prime}}$ intercombination multiplet of Si II at 2335 Angs are significantly reduced, as are those for ratios containing the 3s$^{2}$3p $^{2}$P$_{J}$--3s3p$^{2}$ $^{2}$D$_{J^{\prime}}$ transitions at 1816 Angs. This is primarily due to the effect of the new Si II transition probabilities. However, these atomic data are not only very different from previous calculations, but also show large disagreements with measurements, specifically those of Calamai et. al. (1993) for the intercombination lines. New measurements of transition probabilities for Si II are hence urgently required to confirm (or otherwise) the accuracy of the recently calculated values. If the new calculations are confirmed, then a long-standing discrepancy between theory and observation will have finally been resolved. However, if the older measurements are found to be correct, then the agreement between theory and observation is simply a coincidence and the existing discrepancies remain.

We describe a new atomic and molecular database we developed for use in the spectral synthesis code Cloudy. The design of Stout is driven by the data needs of Cloudy, which simulates molecular, atomic, and ionized gas with kinetic temperatures 2.8 K < T < 1e10 K and densities spanning the low to high-density limits. The radiation field between photon energies $10^{-8}$ Ry and 100 MeV is considered, along with all atoms and ions of the lightest 30 elements, and ~100 molecules. For ease of maintenance, the data are stored in a format as close as possible to the original data sources. Few data sources include the full range of data we need. We describe how we fill in the gaps in the data or extrapolate rates beyond their tabulated range. We tabulate data sources both for the atomic spectroscopic parameters and for collision data for the next release of Cloudy. This is not intended as a review of the current status of atomic data, but rather a description of the features of the database which we will build upon.

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.

Many massive galaxies at the centres of relaxed galaxy clusters and groups have vast reservoirs of cool (~10,000 K) and cold (~100 K) gas. In many low redshift brightest group and cluster galaxies this gas is lifted into the hot ISM in filamentary structures, which are long lived and are typically not forming stars. Two important questions are how far do these reservoirs cool and if cold gas is abundant what is the cause of the low star formation efficiency? Heating and excitation of the filaments from collisions and mixing of hot particles in the surrounding X-ray gas describes well the optical and near infra-red line ratios observed in the filaments. In this paper we examine the theoretical properties of dense, cold clouds emitting in the far infra-red and submillimeter through the bright lines of [C II]157 \mu m , [O I]63 \mu m and CO, exposed to these energetic ionising particles. While some emission lines may be optically thick we find this is not sufficient to model the emission line ratios. Models where the filaments are supported by thermal pressure support alone also cannot account for the cold gas line ratios but a very modest additional pressure support, either from turbulence or magnetic fields can fit the observed [O I]/[C II] line ratios by decreasing the density of the gas. This may also help stabilise the filaments against collapse leading to the low rates of star formation. Finally we make predictions for the line ratios expected from cold gas under these conditions and present diagnostic diagrams for comparison with further observations. We provide our code as an Appendix.

In our series of papers presenting the Herschel imaging of evolved planetary nebulae, we present images of the dust distribution in the Helix nebula (NGC 7293). Images at 70, 160, 250, 350, and 500 micron were obtained with the PACS and SPIRE instruments on board the Herschel satellite. The broadband maps show the dust distribution over the main Helix nebula to be clumpy and predominantly present in the barrel wall. We determined the spectral energy distribution of the main nebula in a consistent way using Herschel, IRAS, and Planck flux values. The emissivity index of 0.99 +/- 0.09, in combination with the carbon rich molecular chemistry of the nebula, indicates that the dust consists mainly of amorphous carbon. The dust excess emission from the central star disk is detected at 70 micron and the flux measurement agree with previous measurement. We present the temperature and dust column density maps. The total dust mass across the Helix nebula (without its halo) is determined to be 0.0035 solar mass at a distance of 216 pc. The temperature map shows dust temperatures between 22 and 42 K, which is similar to the kinetic temperature of the molecular gas, strengthening the fact that the dust and gas co-exist in high density clumps. Archived images are used to compare the location of the dust emission in the far infrared (Herschel) with the ionized (GALEX, Hbeta) and molecular hydrogen component. The different emission components are consistent with the Helix consisting of a thick walled barrel-like structure inclined to the line of sight. The radiation field decreases rapidly through the barrel wall.

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 micron data and improved imaging quality at 100 and 160 micron compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 micron [O I] line flux, eliminating the possibility that line contaminations distort the previously estimated dust mass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 micron flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5+-0.1 Msun of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3 Msun of amorphous carbon and 0.5 Msun of silicates, totalling 0.8 Msun of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies.

The evolution of central stars of planetary nebulae was so far documented in just a few cases. However, spectra collected a few decades ago may provide a good reference for studying the evolution of central stars using the emission line fluxes of their nebulae. We investigated evolutionary changes of the [OIII] 5007 A line flux in the spectra of planetary nebulae. We compared nebular fluxes collected during a decade or longer. We used literature data and newly obtained spectra. A grid of Cloudy models was computed using existing evolutionary models, and the models were compared with the observations. An increase of the [OIII] 5007 A line flux is frequently observed in young planetary nebulae hosting H-rich central stars. The increasing nebular excitation is the response to the increasing temperature and hardening radiation of the central stars. We did not observe any changes in the nebular fluxes in the planetary nebulae hosting late-type Wolf-Rayet (WR) central stars. This may indicate a slower temperature evolution (which may stem from a different evolutionary status) of late-[WR] stars. In young planetary nebulae with H-rich central stars, the evolution can be followed using optical spectra collected during a decade or longer. The observed evolution of H-rich central stars is consistent with the predictions of the evolutionary models provided in the literature. Late-[WR] stars possibly follow a different evolutionary path.

We announce a new facility in the spectral code CLOUDY that enables tracking the evolution of a cooling parcel of gas with time. For gas cooling from temperatures relevant to galaxy clusters, earlier calculations estimated the [Fe XIV] {\lambda}5303 / [Fe X] {\lambda}6375 luminosity ratio, a critical diagnostic of a cooling plasma, to slightly less than unity. By contrast, our calculations predict a ratio ~3. We revisit recent optical coronal line observations along the X-ray cool arc around NGC 4696 by Canning et al. (2011), which detected [Fe X] {\lambda}6375, but not [Fe XIV] {\lambda}5303. We show that these observations are not consistent with predictions of cooling flow models. Differential extinction could in principle account for the observations, but it requires extinction levels (A_V > 3.625) incompatible with previous observations. The non-detection of [Fe XIV] implies a temperature ceiling of 2.1 million K. Assuming cylindrical geometry and transonic turbulent pressure support, we estimate the gas mass at ~1 million solar masses. The coronal gas is cooling isochorically. We propose that the coronal gas has not condensed out of the intracluster medium, but instead is the conductive or mixing interface between the X-ray plume and the optical filaments. We present a number of emission lines that may be pursued to test this hypothesis and constrain the amount of intermediate temperature gas in the system.

V4334 Sgr (Sakurai's object) is an enigmatic evolved star that underwent a very late thermal pulse a few years before its discovery in 1996. It ejected a new, hydrogen-deficient nebula in the process. Emission lines from the newly ejected gas were first discovered in 1998 (He I 1083 nm) and 2001 (optical). We have monitored the optical emission spectrum since. From 2001 through 2007 the optical spectrum showed an exponential decline in flux, consistent with a shock that occurred around 1998 and started cooling soon after that. In this paper we show that since 2008 the line fluxes have been continuously rising again. Our preliminary interpretation is that this emission comes from a region close to the central star, and is excited by a second shock. This shock may have been induced by an increase in the stellar mass loss and wind velocity associated with a rise in the stellar temperature.

We present a framework for the interpretation of the far-infrared spectra of galaxies in which we have expanded the model parameters compared with previous work by varying the ionization parameter $U$, column density $N$(H), and gas density at the cloud face $n$(H$^{+}$) for a central starburst or AGN. We compare these models carried out with the $Cloudy$ spectral synthesis code to trends in line-to-total far-infrared luminosity ratios, far-infrared fine-structure line ratios, IRAS colors, and OH and H$_{2}$O column densities with trends found in the well-studied sample of ten nearby galaxies from the IRAS Bright Galaxy Sample with infrared luminosities greater than 10$^{10}$ L$_{\odot}$ and IRAS 60 micron fluxes equal to or greater than that of the nearby ULIRG Arp 220. We find that the spectral sequence extending from normal starburst-type emission line spectra to ULIRG-type absorption-dominated spectra with significant absorption from excited levels, can be best explained by simultaneously increasing the hydrogen column density, from as low as 10$^{21}$ cm$^{-2}$ to as high as 10$^{24.8}$ cm$^{-2}$ or greater, and the ionization parameter, from as low as 10$^{-4}$ to as high as 1. The starburst models best reproduce most of the sequence, while AGN models are somewhat better able to produce the high OH and H$_{2}$O column densities in Arp 220. Our results suggest that the molecular ISM in ULIRG-like, absorption-dominated systems is located close to and at least partially obscures the source of power throughout much of the far-infrared, which must be taken into account in order to properly interpret diagnostics of both their sources of power and of feedback.

Modern 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 emitted by an ionized plasma. Until now no set of data exists that would meet this need in a fully satisfactory way. We have therefore undertaken to produce a table of very accurate non-relativistic Gaunt factors over a much wider range of parameters than has ever been produced before. We first produced a table of non-averaged Gaunt factors, covering the parameter space log10(epsilon_i) = -20 to +10 and log10(w) = -30 to +25. We then continued to produce a table of thermally averaged Gaunt factors covering the parameter space log10(gamma^2) = -6 to +10 and log10(u) = -16 to +13. Finally we produced a table of the frequency integrated Gaunt factor covering the parameter space log10(gamma^2) = -6 to +10. All the data presented in this paper are available online.

We report the first detections of OH$^+$ emission in planetary nebulae (PNe). As part of an imaging and spectroscopy survey of 11 PNe in the far-IR using the PACS and SPIRE instruments aboard the Herschel Space Observatory, we performed a line survey in these PNe over the entire spectral range between 51 and 672$\mu$m to look for new detections. OH$^+$ rotational emission lines at 152.99, 290.20, 308.48, and 329.77$\mu$m were detected in the spectra of three planetary nebulae: NGC 6445, NGC 6720, and NGC 6781. Excitation temperatures and column densities derived from these lines are in the range of 27 to 47 K and 2$\times$10$^{10}$ to 4 $\times$10$^{11}$ cm$^{-2}$, respectively. In PNe, the OH+ rotational line emission appears to be produced in the photodissociation region (PDR) in these objects. The emission of OH+ is observed only in PNe with hot central stars (T$_{eff}$ > 100000 K), suggesting that high-energy photons may play a role in the OH+ formation and its line excitation in these objects, as it seems to be the case for ultraluminous galaxies.

This work follows Lykins et al. discussion of classic plasma cooling function at low density and solar metallicity. Here we focus on how the cooling function changes over a wide range of density (n_H<10^12 cm^(-3)) and metallicity (Z<30Z _sun ). We find that high densities enhance the ionization of elements such as hydrogen and helium until they reach local thermodynamic equilibrium. By charge transfer, the metallicity changes the ionization of hydrogen when it is partially ionized. We describe the total cooling function as a sum of four parts: those due to H&He, the heavy elements, electron-electron bremsstrahlung and grains. For the first 3 parts, we provide a low-density limit cooling function, a density dependence function, and a metallicity dependence function. These functions are given with numerical tables and analytical fit functions. For grain cooling, we only discuss in ISM case. We then obtain a total cooling function that depends on density, metallicity and temperature. As expected, collisional de-excitation suppresses the heavy elements cooling. Finally, we provide a function giving the electron fraction, which can be used to convert the cooling function into a cooling rate.

We analysed the planetary nebula Hen 2-260 using optical spectroscopy and photometry. We compared our observations with the data from literature to search for evolutionary changes. The nebular line fluxes were modelled with the Cloudy photoionization code to derive the stellar and nebular parameters. The planetary nebula shows a complex structure and possibly a bipolar outflow. The nebula is relatively dense and young. The central star is just starting $\rm O^+$ ionization ($\rm T_{eff} \approx 30,000 \, K$). Comparison of our observations with literature data indicates a 50% increase of the [OIII] 5007 \AA\ line flux between 2001 and 2012. We interpret it as the result of the progression of the ionization of $\rm O^{+}$. The central star evolves to higher temperatures at a rate of $\rm 45 \pm 7\,K\, yr^{-1}$. The heating rate is consistent with a final mass of $\rm 0.626 ^{+0.003}_{-0.005} \, M_{\odot}$ or $\rm 0.645 ^{+0.008}_{-0.008} \, M_{\odot}$ for two different sets of post-AGB evolutionary tracks from literature. The photometric monitoring of Hen 2-260 revealed variations on a timescale of hours or days. The variability may be caused by pulsations of the star. The temperature evolution of the central star can be traced using spectroscopic observations of the surrounding planetary nebula spanning a timescale of roughly a decade. This allows us to precisely determine the stellar mass, since the pace of the temperature evolution depends critically on the core mass. The kinematical age of the nebula is consistent with the age obtained from the evolutionary track. The final mass of the central star is close to the mass distribution peak for central stars of planetary nebulae found in other studies. The object belongs to a group of young central stars of planetary nebulae showing photometric variability.

The Helix Nebula (NGC 7293) is the closest planetary nebulae. Therefore, it is an ideal template for photochemical studies at small spatial scales in planetary nebulae. We aim to study the spatial distribution of the atomic and the molecular gas, and the structure of the photodissociation region along the western rims of the Helix Nebula as seen in the submillimeter range with Herschel. We use 5 SPIRE FTS pointing observations to make atomic and molecular spectral maps. We analyze the molecular gas by modeling the CO rotational lines using a non-local thermodynamic equilibrium (non-LTE) radiative transfer model. For the first time, we have detected extended OH+ emission in a planetary nebula. The spectra towards the Helix Nebula also show CO emission lines (from J= 4 to 8), [NII] at 1461 GHz from ionized gas, and [CI] (2-1), which together with the OH+ lines, trace extended CO photodissociation regions along the rims. The estimated OH+ column density is (1-10)x1e12 cm-2. The CH+ (1-0) line was not detected at the sensitivity of our observations. Non-LTE models of the CO excitation were used to constrain the average gas density (n(H2)=(1-5)x1e5 cm-3) and the gas temperature (Tk= 20-40 K). The SPIRE spectral-maps suggest that CO arises from dense and shielded clumps in the western rims of the Helix Nebula whereas OH+ and [CI] lines trace the diffuse gas and the UV and X-ray illuminated clumps surface where molecules reform after CO photodissociation. [NII] traces a more diffuse ionized gas component in the interclump medium.

We present 48 Herschel/PACS spectra of evolved stars in the wavelength range of 67-72 $\mu$m. This wavelength range covers the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{2-2x}\text{Fe}_{(2x)}\text{SiO}_{4}$). The width and wavelength position of this band are sensitive to the temperature and composition of the crystalline olivine. Our sample covers a wide range of objects: from high mass-loss rate AGB stars (OH/IR stars, $\dot M \ge 10^{-5}$ M$_\odot$/yr), through post-AGB stars with and without circumbinary disks, to planetary nebulae and even a few massive evolved stars. The goal of this study is to exploit the spectral properties of the 69 $\mu$m band to determine the composition and temperature of the crystalline olivine. Since the objects cover a range of evolutionary phases, we study the physical and chemical properties in this range of physical environments. We fit the 69 $\mu$m band and use its width and position to probe the composition and temperature of the crystalline olivine. For 27 sources in the sample, we detected the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{(2-2x)}\text{Fe}_{(2x)}\text{SiO}_{4}$). The 69 $\mu$m band shows that all the sources produce pure forsterite grains containing no iron in their lattice structure. The temperature of the crystalline olivine as indicated by the 69 $\mu$m band, shows that on average the temperature of the crystalline olivine is highest in the group of OH/IR stars and the post-AGB stars with confirmed Keplerian disks. The temperature is lower for the other post-AGB stars and lowest for the planetary nebulae. A couple of the detected 69 $\mu$m bands are broader than those of pure magnesium-rich crystalline olivine, which we show can be due to a temperature gradient in the circumstellar environment of these stars. continued...

This is the first of a series of investigations into far-IR characteristics of 11 planetary nebulae (PNs) under the Herschel Space Observatory Open Time 1 program, Herschel Planetary Nebula Survey (HerPlaNS). Using the HerPlaNS data set, we look into the PN energetics and variations of the physical conditions within the target nebulae. In the present work, we provide an overview of the survey, data acquisition and processing, and resulting data products. We perform (1) PACS/SPIRE broadband imaging to determine the spatial distribution of the cold dust component in the target PNs and (2) PACS/SPIRE spectral-energy-distribution (SED) and line spectroscopy to determine the spatial distribution of the gas component in the target PNs. For the case of NGC 6781, the broadband maps confirm the nearly pole-on barrel structure of the amorphous carbon-richdust shell and the surrounding halo having temperatures of 26-40 K. The PACS/SPIRE multi-position spectra show spatial variations of far-IR lines that reflect the physical stratification of the nebula. We demonstrate that spatially-resolved far-IR line diagnostics yield the (T_e, n_e) profiles, from which distributions of ionized, atomic, and molecular gases can be determined. Direct comparison of the dust and gas column mass maps constrained by the HerPlaNS data allows to construct an empirical gas-to-dust mass ratio map, which shows a range of ratios with the median of 195+-110. The present analysis yields estimates of the total mass of the shell to be 0.86 M_sun, consisting of 0.54 M_sun of ionized gas, 0.12 M_sun of atomic gas, 0.2 M_sun of molecular gas, and 4 x 10^-3 M_sun of dust grains. These estimates also suggest that the central star of about 1.5 M_sun initial mass is terminating its PN evolution onto the white dwarf cooling track.

We analyse the remnants of CK Vul (Nova Vul 1670) using optical imaging and spectroscopy. The imaging, obtained between 1991 and 2010, spans 5.6% of the life-time of the nebula. The flux of the nebula decreased during the last 2 decades. The central source still maintains the ionization of the innermost part of the nebula, but recombination proceeds in more distant parts of the nebula. Surprisingly, we discovered two stars located within 10 arcsec of the expansion centre of the radio emission that are characterized by pronounced long term variations and one star with high proper motion. The high proper motion star is a foreground object, and the two variable stars are background objects. The photometric variations of two variables are induced by a dusty cloud ejected by CK Vul and passing through the line of sight to those stars. The cloud leaves strong lithium absorption in the spectra of the stars. We discuss the nature of the object in terms of recent observations.

A variety of astronomical environments contain clouds irradiated by a combination of isotropic and beamed radiation fields. For example, molecular clouds may be irradiated by the isotropic cosmic microwave background (CMB), as well as by a nearby active galactic nucleus (AGN). These radiation fields excite atoms and molecules and produce emission in different ways. We revisit the escape probability theorem and derive a novel expression that accounts for the presence of external radiation fields. We show that when the field is isotropic the escape probability is reduced relative to that in the absence of external radiation. This is in agreement with previous results obtained under ad hoc assumptions or with the two-level system, but can be applied to complex many-level models of atoms or molecules. This treatment is in the development version of the spectral synthesis code Cloudy. We examine the spectrum of a Spitzer cloud embedded in the local interstellar radiation field, and show that about 60 percent of its emission lines are sensitive to background subtraction. We argue that this geometric approach could provide an additional tool toward understanding the complex radiation fields of starburst galaxies.

We present an analysis of the far-infrared and submillimetre molecular emission line spectrum of the luminous M-supergiant VY CMa, observed with the SPIRE and PACS spectrometers aboard the Herschel Space Observatory. Over 260 emission lines were detected in the 190-650-micron SPIRE FTS spectra, with one-third of the observed lines being attributable to H2O. Other detected species include CO, 13CO, H2^18O, SiO, HCN, SO, SO2, CS, H2S, and NH3. Our model fits to the observed 12CO and 13CO line intensities yield a 12C/13C ratio of 5.6+-1.8, consistent with measurements of this ratio for other M supergiants, but significantly lower than previously estimated for VY CMa from observations of lower-J lines. The spectral line energy distribution for twenty SiO rotational lines shows two temperature components: a hot component at 1000 K, which we attribute to the stellar atmosphere and inner wind, plus a cooler ~200 K component, which we attribute to an origin in the outer circumstellar envelope. We fit the line fluxes of 12CO, 13CO, H2O and SiO, using the SMMOL non-LTE line transfer code, with a mass-loss rate of 1.85x10^-4 Msun yr^-1 between 9 R* and 350 R*. To fit the observed line fluxes of 12CO, 13CO, H2O and SiO with SMMOL non-LTE line radiative transfer code, along with a mass-loss rate of 1.85x10^-4 Msun yr^-1. To fit the high rotational lines of CO and H2O, the model required a rather flat temperature distribution inside the dust condensation radius, attributed to the high H2O opacity. Beyond the dust condensation radius the gas temperature is fitted best by an r^-0.5 radial dependence, consistent with the coolant lines becoming optically thin. Our H2O emission line fits are consistent with an ortho:para ratio of 3 in the outflow.

As part of the Herschel Guaranteed Time Key Project MESS (Mass loss of Evolved StarS) we have imaged a sample of planetary nebulae. In this paper we present the PACS and SPIRE images of the classical bipolar planetary nebula NGC 650. We used these images to derive a temperature map of the dust. We also constructed a photoionization and dust radiative transfer model using the spectral synthesis code Cloudy. To constrain this model, we used the PACS and SPIRE fluxes and combined these with hitherto unpublished IUE and Spitzer IRS spectra as well as various other data from the literature. The temperature map combined with the photoionization model were used to study various aspects of the central star, the nebula, and in particular the dust grains in the nebula. The central star parameters are determined to be T_eff = 208 kK and L = 261 L_sol assuming a distance of 1200 pc. The stellar temperature is much higher than previously published values. We confirm that the nebula is carbon-rich with a C/O ratio of 2.1. The nebular abundances are typical for a type IIa planetary nebula. With the photoionization model we determined that the grains in the ionized nebula are large (assuming single-sized grains, they would have a radius of 0.15 micron. Most likely these large grains were inherited from the asymptotic giant branch phase. The PACS 70/160 micron temperature map shows evidence for two radiation components heating the grains. The first component is direct emission from the central star, while the second component is diffuse emission from the ionized gas (mainly Ly alpha). We show that previous suggestions that there is a photo-dissociation region surrounding the ionized region are incorrect. The neutral material resides in dense clumps inside the ionized region. These may also harbor stochastically heated very small grains in addition to the large grains.

We report spectroscopic and imaging observations of rotational transitions of cold CO and SiO in the ejecta of SN1987A, the first such emission detected in a supernova remnant. In addition to line luminosities for the CO J=1-0, 2-1, 6-5, and 7-6 transitions, we present upper limits for all other transitions up to J=13-12, collectively measured from the Atacama Large Millimeter Array (ALMA), the Atacama Pathfinder EXperiment (APEX), and the Herschel Spectral and Photometric Imaging REceiver (SPIRE). Simple models show the lines are emitted from at least 0.01 solar masses of CO at a temperature > 14 K, confined within at most 35% of a spherical volume expanding at ~ 2000 km/s. Moreover, we locate the emission within 1'' of the central debris. These observations, along with a partial observation of SiO, confirm the presence of cold molecular gas within supernova remnants and provide insight into the physical conditions and chemical processes in the ejecta. Furthermore, we demonstrate the powerful new window into supernova ejecta offered by submillimeter observations.

Planetary nebulae (PNe) are powerful tracers of evolved stellar populations. Among the 3000 known PNe in the Galaxy, about 600 are located within the 520 square-degree area covered by the VVV survey. The VVV photometric catalogue provides an important new dataset for the study of PNe, with high-resolution imaging in five near-infrared bands. Aperture photometry of known PNe in the VVV area was retrieved from source catalogues. Care was taken to minimise any confusion with field stars. The colours of the PNe we are determined for H-Ks, J-H, Z-Y, and Y-J, and compared to stars and to other types of emission line objects. Cloudy photo-ionisation models were used to predict colours for typical PNe. We present near-infrared photometry for 353 known PNe. The best separation from other objects is obtained in the H-Ks vs. J-H diagram. We calculated the emission-line contribution to the in-band flux based on a model for NGC 6720: we find that this is highest in the Z and Y bands at over 50%, lower in the J band at 40%, and lowest in the H and Ks bands at 20%. A new view of PNe in the wavelength domain of the Z and Y bands is shown. Photo-ionisation models are used to explore the observed colours in these bands. The Y band is shown to be dominated by HeI 1.083 mu and HeII 1.012 mu, and colours involving this band are very sensitive to the temperature of the ionizing star. The VVV survey represents a unique dataset for studing crowded and obscured regions in the Galactic plane. The diagnostic diagrams presented here allow one to study the properties of known PNe and to uncover objects not previously classified.

The Fe unresolved transition array (UTAs) produce prominent features in the 15-17?A wavelength range in the spectra of Active Galactic Nuclei (AGN). Here we present new calculations of the energies and oscillator strengths of inner- shell lines from Fe XIV, Fe XV, and Fe XVI. These are crucial ions since they are dominant at inflection points in the gas thermal stability curve, and UTA excitation followed by autoionization is an important ionization mechanism for these species. We incorporate these, and data reported in previous papers, into the plasma simulation code Cloudy. This updated physics is subsequently employed to reconsider the thermally stable phases in absorbing media in Active Galactic Nuclei. We show how the absorption profile of the Fe XIV UTA depends on density, due to the changing populations of levels within the ground configuration.

This is a summary of the 2013 release of the plasma simulation code Cloudy. Cloudy models the ionization, chemical, and thermal state of material that may be exposed to an external radiation field or other source of heating, and predicts observables such as emission and absorption spectra. It works in terms of elementary processes, so is not limited to any particular temperature or density regime. This paper summarizes advances made since the last major review in 1998. Much of the recent development has emphasized dusty molecular environments, improvements to the ionization / chemistry solvers, and how atomic and molecular data are used. We present two types of simulations to demonstrate the capability of the code. We consider a molecular cloud irradiated by an X-ray source such as an Active Nucleus and show how treating EUV recombination lines and the full SED affects the observed spectrum. A second example illustrates the very wide range of particle and radiation density that can be considered.

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.

The optical [N I] doublet near 5200 {\AA} is anomalously strong in a variety of emission-line objects. We compute a detailed photoionization model and use it to show that pumping by far-ultraviolet (FUV) stellar radiation previously posited as a general explanation applies to the Orion Nebula (M42) and its companion M43; but, it is unlikely to explain planetary nebulae and supernova remnants. Our models establish that the observed nearly constant equivalent width of [N I] with respect to the dust-scattered stellar continuum depends primarily on three factors: the FUV to visual-band flux ratio of the stellar population; the optical properties of the dust; and the line broadening where the pumping occurs. In contrast, the intensity ratio [N I]/H{\beta} depends primarily on the FUV to extreme-ultraviolet ratio, which varies strongly with the spectral type of the exciting star. This is consistent with the observed difference of a factor of five between M42 and M43, which are excited by an O7 and B0.5 star respectively. We derive a non-thermal broadening of order 5 km/s for the [N I] pumping zone and show that the broadening mechanism must be different from the large-scale turbulent motions that have been suggested to explain the line-widths in this H II region. A mechanism is required that operates at scales of a few astronomical units, which may be driven by thermal instabilities of neutral gas in the range 1000 to 3000 K. In an appendix, we describe how collisional and radiative processes are treated in the detailed model N I atom now included in the Cloudy plasma code.

Far ultraviolet emission has been detected from a knot of Halpha emission in the Horseshoe filament, far out in the NGC 1275 nebula. The flux detected relative to the brightness of the Halpha line in the same spatial region is very close to that expected from Hydrogen two-photon continuum emission in the particle heating model of Ferland et al. (2009) if reddening internal to the filaments is taken into account. We find no need to invoke other sources of far ultraviolet emission such as hot stars or emission lines from CIV in intermediate temperature gas to explain these data.

In this paper we give a progress report on the Herschel imaging and spectroscopic observations of planetary nebulae that are carried out as part of the MESS guaranteed time key program. We present and discuss imaging and temperature maps of NGC 6720, NGC 650, and NGC 6853, as well as PACS and SPIRE spectroscopy of NGC 7027.

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.

We have supplemented existing spectra of Barnard's Loop with high accuracy spectrophotometry of one new position. Cloudy photoionization models were calculated for a variety of ionization parameters and stellar temperatures and compared with the observations. After testing the procedure with recent observations of M43, we establish that Barnard's Loop is photoionized by four candidate ionizing stars, but agreement between the models and observations is only possible if Barnard's Loop is enhanced in heavy elements by about a factor of 1.4. Barnard's Loop is very similar in properties to the brightest components of the Orion-Eridanus Bubble and the Warm Ionized Medium (WIM). We are able to establish models that bound the range populated in low-ionization color-color diagrams (I([SII])/I(H{\alpha}) versus I([NII])/I(H{\alpha})) using only a limited range of ionization parameters and stellar temperatures. Previously established variations in the relative abundance of heavy elements render uncertain the most common method of determining electron temperatures for components of the Orion-Eridanus Bubble and the WIM based on only the I([NII])/I(H{\alpha}) ratio, although we confirm that the lowest surface brightness components of the WIM are on average of higher electron temperature. The electron temperatures for a few high surface brightness WIM components determined by direct methods are comparable to those of classical bright H II regions. In contrast, the low surface brightness HII regions studied by the Wisconsin H{\alpha} Mapper are of lower temperatures than the classical bright HII regions.

MESS (Mass-loss of Evolved StarS) is a Guaranteed Time Key Program that uses the PACS and SPIRE instruments on board the Herschel Space Observatory to observe a representative sample of evolved stars, that include asymptotic giant branch (AGB) and post-AGB stars, planetary nebulae and red supergiants, as well as luminous blue variables, Wolf-Rayet stars and supernova remnants. In total, of order 150 objects are observed in imaging and about 50 objects in spectroscopy. This paper describes the target selection and target list, and the observing strategy. Key science projects are described, and illustrated using results obtained during Herschel's science demonstration phase. Aperture photometry is given for the 70 AGB and post-AGB stars observed up to October 17, 2010, which constitutes the largest single uniform database of far-IR and sub-mm fluxes for late-type stars.

In this paper we will discuss the images of Planetary Nebulae that have recently been obtained with PACS and SPIRE on board the Herschel satellite. This comprises results for NGC 650 (the little Dumbbell nebula), NGC 6853 (the Dumbbell nebula), and NGC 7293 (the Helix nebula).

Herschel PACS and SPIRE images have been obtained of NGC 6720 (the Ring Nebula). This is an evolved planetary nebula with a central star that is currently on the cooling track, due to which the outer parts of the nebula are recombining. From the PACS and SPIRE images we conclude that there is a striking resemblance between the dust distribution and the H2 emission, which appears to be observational evidence that H2 forms on grain surfaces. We have developed a photoionization model of the nebula with the Cloudy code which we used to determine the physical conditions of the dust and investigate possible formation scenarios for the H2. We conclude that the most plausible scenario is that the H2 resides in high density knots which were formed after the recombination of the gas started when the central star entered the cooling track. Hydrodynamical instabilities due to the unusually low temperature of the recombining gas are proposed as a mechanism for forming the knots. H2 formation in the knots is expected to be substantial after the central star underwent a strong drop in luminosity about one to two thousand years ago, and may still be ongoing at this moment, depending on the density of the knots and the properties of the grains in the knots.

We investigate consequences of the discovery that Fe II emission in quasars, one of the spectroscopic signatures of "Eigenvector 1", may originate in infalling clouds. Eigenvector 1 correlates with the Eddington ratio L/L_Edd so that Fe II/Hbeta increases as L/L_Edd increases. We show that the "force multiplier", the ratio of gas opacity to electron scattering opacity, is ~ 10^3 - 10^4 in Fe II-emitting gas. Such gas would be accelerated away from the central object if the radiation force is able to act on the entire cloud. As had previously been deduced, infall requires that the clouds have large column densities so that a substantial amount of shielded gas is present. The critical column density required for infall to occur depends on L/L_Edd, establishing a link between Eigenvector 1 and the Fe II/Hbeta ratio. We see predominantly the shielded face of the infalling clouds rather than the symmetric distribution of emitters that has been assumed. The Fe II spectrum emitted by the shielded face is in good agreement with observations thus solving several long-standing mysteries in quasar emission lines.

The Orion Bar is one of the nearest and best-studied photodissociation or photon-dominated regions (PDRs). Observations reveal the presence of H2 lines from vibrationally or rotationally excited upper levels that suggest warm gas temperatures (400 to 700 K). However, standard models of PDRs are unable to reproduce such warm rotational temperatures. In this paper we attempt to explain these observations with new comprehensive models which extend from the H+ region through the Bar and include the magnetic field in the equation of state. We adopt the model parameters from our previous paper which successfully reproduced a wide variety of spectral observations across the Bar. In this model the local cosmic-ray density is enhanced above the galactic background, as is the magnetic field, and which increases the cosmic-ray heating elevating the temperature in the molecular region. The pressure is further enhanced above the gas pressure in the H+ region by the momentum transferred from the absorbed starlight. Here we investigate whether the observed H2 lines can be reproduced with standard assumptions concerning the grain photoelectric emission. We also explore the effects due to the inclusion of recently computed H2 + H2, H2 + H and H2 + He collisional rate coefficients.

The observed faintness of infrared fine-structure line emission along with the warm far-infrared (FIR) colors of ultraluminous infrared galaxies (ULIRGs) is a long-standing problem. In this work, we calculate the line and continuum properties of a cloud exposed to an Active Galactic Nucleus (AGN) and starburst spectral energy distribution (SED). We use an integrated modeling approach, predicting the spectrum of ionized, atomic, and molecular environments in pressure equilibrium. We find that the effects of high ratios of impinging ionizing radiation density to particle density (i.e. high ionization parameters, or U) can reproduce many ULIRG observational characteristics. Physically, as U increases, the fraction of UV photons absorbed by dust increases, corresponding to fewer photons available to photoionize and heat the gas, producing what is known as a "dust-bounded" nebula. We show that high U effects can explain the "[C II] deficit", the ~1 dex drop in the [C II] 158 micron /FIR ratio seen in ULIRGs when compared to starburst or normal galaxies. Additionally, by increasing U through increasing the ionizing photon flux, warmer dust and thus higher IRAS F(60)/F(100) ratios result. High U effects also predict an increase in [O I]63 micron /[C II] 158 micron and a gradual decline in [O III] 88 micron /FIR, similar to the magnitude of the trends observed, and yield a reasonable fit to [Ne V]14 micron /FIR ratio AGN observations.

[Abridged] Photoionization, whether by starlight or other sources, has difficulty in accounting for the observed spectra of the optical filaments that often surround central galaxies in large clusters. Our first paper examined whether heating by energetic particles or dissipative MHD wave can account for the observations. Here we include atomic and low-ionization regions. The model of the hydrogen atom, along with all elements of the H-like iso-electronic sequence, is now fully nl-resolved. We show how the predicted HI spectrum differs from the pure recombination case. The second update is to the rates for H^0 - H2 inelastic collisions. We now use the values computed by Wrathmall et al. The rates are often much larger. We calculate the chemistry, ionization, temperature, gas pressure, and emission-line spectrum for a wide range of gas densities and collisional heating rates. We assume that the filaments are magnetically confined and free to move along field lines so that the gas pressure is equal to that of the surrounding hot gas. A mix of clouds, some being dense and cold and others hot and tenuous, can exist. The observed spectrum will be the integrated emission from clouds with different densities and temperatures but the same pressure. We assume that the gas filling factor is given by a power law in density. The power-law index is set by matching the observed intensities of IR H2 lines relative to optical HI lines. We conclude that the filaments are heated by ionizing particles, either conducted in from surrounding regions or produced in situ by processes related to MHD waves.

Our understanding of physical processes in Photodissociation regions or Photon Dominated Regions (PDRs) largely depends on the ability of spectral synthesis codes to reproduce the observed infrared emission-line spectrum. In this paper, we explore the sensitivity of a single PDR model to microphysical details. Our calculations use the Cloudy spectral synthesis code, recently modified to include a wealth of PDR physical processes. We show how the chemical/thermal structure of a PDR, along with the calculated spectrum, changes when the treatment of physical processes such as grain physics and atomic/molecular rates are varied. We find a significant variation in the intensities of PDR emission lines, depending on different treatments of the grain physics. We also show how different combinations of the cosmic-ray ionization rate, inclusion of grain-atom/ion charge transfer, and the grain size distribution can lead to very similar results for the chemical structure. Additionally, our results show the utility of Cloudy for the spectral modeling of molecular environments.

While in the past spherodicity was assumed, and still is used in modeling of most nebulae, we know now that only a small number of planetary nebulae (PNe) are really spherical or at least nearly round. Round planetary nebulae are the minority of objects. In the case of those objects that underwent a very late helium flash (called VLTP objects or ``born-again'' PNe) it seems to be different. The first, hydrogen-rich PN, is more or less round. The ejecta from the VLTP event, in contrast, are extremely asymmetrical.

While in the past spheroidicity was assumed, and still is used in modeling of most nebulae, we know now that only a small number of planetary nebulae (PNe) are really spherical or at least nearly round. Round planetary nebulae are the minority of objects. In case of those objects that underwent a very late helium flash (called VLTP or "born-again" PNe) it seems to be different. The first, hydrogen rich PN, is more or less round. The ejecta from the VLTP event is extremely asymmetrically. Angular momentum is mostly assumed to be the main reason for the asymmetry in PNe. Thus we have to find processes either changing their behavior within a few hundred to a few thousands of years or change their properties dramatically due to the variation of the abundance. They most likely have a strong link or dependency with the abundance of the ejecta.

The optical filaments found in many cooling flows in galaxy clusters consist of low density ($\sim 10^3 \pcc$) cool ($\sim 10^3$ K) gas surrounded by significant amounts of cosmic-ray and magnetic-field energy. Their spectra show anomalously strong low-ionization and molecular emission lines when compared with galactic molecular clouds exposed to ionizing radiation such as the Orion complex. Previous studies have shown that the spectra cannot be produced by O-star photoionization. Here we calculate the physical conditions in dusty gas that is well shielded from external sources of ionizing photons and is energized either by cosmic rays or dissipative MHD waves. Strong molecular hydrogen lines, with relative intensities similar to those observed, are produced. Selection effects introduced by the microphysics produce a correlation between the \htwo line upper level energy and the population temperature. These selection effects allow a purely collisional gas to produce \htwo emission that masquerades as starlight-pumped \htwo but with intensities that are far stronger. This physics may find application to any environment where a broad range of gas densities or heating rates occur.

We present the results of a 25-year program to monitor the radio flux evolution of the planetary nebula NGC7027. We find significant evolution of the spectral flux densities. The flux density at 1465 MHz, where the nebula is optically thick, is increasing at a rate of 0.251+-0.015 % per year, caused by the expansion of the ionized nebula. At frequencies where the emission is optically thin, the spectral flux density is changing at a rate of -0.145+-0.005 % per year, caused by a decrease in the number of ionizing photons coming from the central star. A distance of 980+-100 pc is derived. By fitting interpolated models of post-AGB evolution to the observed changes, we find that over the 25-yr monitoring period, the stellar temperature has increased by 3900+-900 K and the stellar bolometric luminosity has decreased by 1.75+-0.38 %. We derive a distance-independent stellar mass of 0.655+-0.01 solar masses adopting the Bloecker stellar evolution models, or about 0.04 solar masses higher when using models of Vassiliadis & Wood which may provide a better fit. A Cloudy photoionization model is used to fit all epochs at all frequencies simultaneously. The differences between the radio flux density predictions and the observed values show some time-independent residuals of typically 1 %. A possible explanation is inaccuracies in the radio flux scale of Baars et al. We propose an adjustment to the flux density scale of the primary radio flux calibrator 3C286, based on the Cloudy model of NGC7027. We also calculate precise flux densities for NGC7027 for all standard continuum bands used at the VLA, as well as for some new 30GHz experiments.

context: The spectrum of IRAS 16594-4656 shows shock excited H_2 emission and collisionally excited emission lines such as[O I],[C I],and [Fe II]. aim: The goal is to determine the location of the H_2 and [Fe II] shock emission, to determine the shock velocities,and constrain the physical properties in the shock. methods: High resolution spectra of the H_2 1-0 S(1),H_2 2-1 S(1), [Fe II], and Pa$\beta$ emission lines were obtained with the near infrared spectrograph Phoenix on Gemini South. results: The position-velocity diagrams of H_2 1-0 S(1), H_2 2-1 S(1), and [Fe II] are presented. The H_2 and [Fe II] emission is spatially extended. The collisionally excited [O I] and [C I] optical emission lines have a similar double peaked profile compared to the extracted H_2 profile and appear to be produced in the same shock. They all indicate an expansion velocity of ~8 km/s and the presence of a neutral, very high density region with $n_{\rm e}$ about 3 x 10^6 to 5 x10^7 cm$^{-3}$. The [Fe II] emission however is single peaked. It has a gaussian FWHM of 30 km/s and a total width of 62 km/s at 1% of the peak. The Pa$\beta$ profile is even wider with a gaussian FWHM of 48 km/s and a total width of 75 km/s at 1% of the peak. conclusions: The H$_2$ emission is excited in a slow 5 to 20 km/s shock into dense material at the edge of the lobes, caused by the interaction of the AGB ejecta and the post-AGB wind. The 3D representation of the H_2 data shows a hollow structure with less H_2 emission in the equatorial region. The [Fe II] emission is not present in the lobes, but originates close to the central star in fast shocks in the post-AGB wind or in a disk. The Pa$\beta$ emission also appears to originate close to the star.

The spatial distribution of the cosmic-ray flux is important in understanding the Interstellar Medium (ISM) of the Galaxy. This distribution can be analyzed by studying different molecular species along different sight lines whose abundances are sensitive to the cosmic-ray ionization rate. Recently several groups have reported an enhanced cosmic-ray ionization rate in diffuse clouds compared to the standard value, zeta(average)=2.5e-17 s^-1, measured toward dense molecular clouds. In an earlier work we reported an enhancement in the cosmic ray rate of 20 towards HD185418. McCall et al. have reported enhancements of 48 towards zeta Persei based on the observed abundance of H$_{3}^+ while Le Petit et al. found a cosmic ray enhancement of ~10 to be consistent with their models for this same sight line. Here we revisit zeta Persei and perform a detailed calculation using a self-consistent treatment of the hydrogen chemistry, grain physics, energy and ionization balance, and excitation physics. We show that the value of zeta deduced from the H$_{3}^+$ column density in the diffuse region of the sightline depends strongly on the properties of the grains because they remove free electrons and change the hydrogen chemistry. The observations are largely consistent with a cosmic ray enhancement of 40, with several diagnostics indicating higher values. This underscores the importance of a full treatment of grain physics in studies of interstellar chemistry.

CK Vul is classified as, amongst others, the slowest known nova, a hibernating nova, or a very late thermal pulse object. Following its eruption in AD 1670, the star remained visible for 2 years. A 15-arcsec nebula was discovered in the 1980's, but the star itself has not been detected since the eruption. We here present radio images which reveal an 0.1-arcsec radio source with a flux of 1.5 mJy at 5 GHz. Deep Halpha images show a bipolar nebula with a longest extension of 70 arcsec, with the previously known compact nebula at its waist. The emission-line ratios show that the gas is shock-ionized, at velocities >100 km/s. Dust emission yields an envelope mass of ~0.05 Msun. Echelle spectra indicate outflow velocities up to 360 km/s. From a comparison of images obtained in 1991 and 2004 we find evidence for expansion of the nebula, consistent with an origin in the 1670 explosion; the measured expansion is centred on the radio source. No optical or infrared counterpart is found at the position of the radio source. The radio emission is interpreted as thermal free-free emission from gas with Tel ~ 10000 K. The radio source may be due to a remnant circumbinary disk, similar to those seen in some binary post-AGB stars. We discuss possible classifications of this unique outburst, including that of a sub-Chandrasekhar mass supernova, a nova eruption on a cool, low-mass white dwarf, or a thermal pulse induced by accretion from a circumbinary disk.

We investigate the reheating of the very late thermal pulse (VLTP) object V4334 Sgr (Sakurai's Object) using radio observations from the Very Large Array, and optical spectra obtained with the Very Large Telescope. We find a sudden rise of the radio flux at 5 and 8 GHz - from <= 90 micro-Jy and 80 +/- 30 micro-Jy in February 2005 to 320 micro-Jy and 280 micro-Jy in June 2006. Optical line emission is also evolving, but the emission lines are fading. The optical line emission and early radio flux are attributed to a fast shock (and not photoionization as was reported earlier) which occurred around 1998. The fading is due to post-shock cooling and recombination. The recent rapid increase in radio flux is evidence for the onset of photoionization of carbon starting around 2005. The current results indicate an increase in the stellar temperature to 12 kK in 2006. The mass ejected in the VLTP eruption is M_ej >= 1e-4 Msol, but could be as high as 1e-2 Msol, depending mainly on the distance and the clumping factor of the outflow. We derive a distance between 1.8 and 5 kpc. A high mass loss could expose the helium layer and yield abundances compatible with those of [WC] and PG1159 stars.

One of the puzzling findings in the study of the chemical evolution of (post-)AGB stars is why very similar stars (in terms of metallicity, spectral type, infrared properties, etc...) show a very different photospheric composition. We aim at extending the still limited sample of s-process enriched post-AGB stars, in order to obtain a statistically large enough sample that allows us to formulate conclusions concerning the 3rd dredge-up occurrence. We selected two post-AGB stars on the basis of IR colours indicative of a past history of heavy mass loss: IRAS08281-4850 and IRAS14325-6428. They are cool sources in the locus of the Planetary Nebulae (PNe) in the IRAS colour-colour diagram. Abundances of both objects were derived for the first time on the basis of high-quality UVES and EMMI spectra, using a critically compiled line list with accurate log(gf) values, together with the latest Kurucz model atmospheres. Both objects have very similar spectroscopically defined effective temperatures of 7750-8000K. They are strongly carbon and s-process enriched, with a C/O ratio of 1.9 and 1.6, and an [ls/Fe] of +1.7 and +1.2, for IRAS08281 and IRAS14325 resp. Moreover, the spectral energy distributions (SEDs) point to heavy mass-loss during the preceding AGB phase. IRAS08281 and IRAS14325 are prototypical post-AGB objects in the sense that they show strong post 3rd dredge-up chemical enrichments. The neutron irradiation has been extremely efficient, despite the only mild sub-solar metallicity. This is not conform with the recent chemical models. The existence of very similar post-AGB stars without any enrichment emphasizes our poor knowledge of the details of the AGB nucleosynthesis and dredge-up phenomena. We call for a very systematic chemical study of all cool sources in the PN region of the IRAS colour-colour diagram.

We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly.

We report on the progress of our on-going campaign to monitor the evolution of the VLTP objects V4334 Sgr and V605 Aql, as well as the suspected (V)LTP object CK Vul. V4334 Sgr does not show signs of increased ionization compared to our previous observations in 2004. We obtained the first radio detection of V605 Aql, indicating a strong increase in radio flux since 1987. We also present the first radio detection of CK Vul and discuss the expansion of the material ejected during the 1670 event.

New Planetary Nebulae (PNe) were discovered through an [O III] 5007 A emission line survey in the Galactic bulge region with l>0 deg. We detected 240 objects, including 44 new PNe. Deep Halpha+[N II] CCD images as well as low resolution spectra were obtained for the new PNe in order to study them in detail. Preliminary photo-ionization models of the new PNe with Cloudy resulted in first estimates of the physical parameters and abundances. They are compared to the abundances of Galactic PNe.

The identification of two new Planetary Nebulae in the Sagittarius Dwarf Spheroidal Galaxy (Sgr) is presented. This brings the total number to four. The first, StWr 2-21, belongs to the main body of Sgr. The second, the halo PN BoBn 1, has a location, distance and velocity in agreement with the leading tidal tail of Sgr. We estimate that 10 per cent of the Galactic halo consists of Sgr debris. The specific frequency of PNe indicates a total luminosity of Sgr, including its tidal tails, of M_V=-14.1. StWr 2-21 shows a high abundance of [O/H]=-0.23, which confirms the high-metallicity population in Sgr uncovered by Bonaficio et al. (2004). The steep metallicity--age gradient in Sgr is due to ISM removal during the Galactic plane passages, ISM reformation due to stellar mass loss, and possibly accretion of metal-enriched gas from our Galaxy. The ISM re-formation rate of Sgr, from stellar mass loss, is 5 X 10^-4 M_sun yr^-1, amounting to ~10^6 M_sun per orbital period. HST images reveal well-developed bipolar morphologies, and provide clear detections of the central stars. All three stars with deep spectra show WR-lines, suggesting that the progenitor mass and metallicity determines whether a PN central star develops a WR spectrum. One Sgr PN belongs to the class of IR-[WC] stars. Expansion velocities are determined for three nebulae. Comparison with hydrodynamical models indicates an initial density profile of rho ~ r^-3. This is evidence for increasing mass-loss rates on the AGB. Peak mass-loss rates are indicated of ~ 10^-4 M_sun yr^-1. The IR-[WC] PN, He 2-436, provides the sole direct detection of dust in a dwarf spheroidal galaxy, to date.

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.

Much of the baryonic matter in the Universe is in the form of H2 which includes most of the gas in galactic and extragalactic interstellar clouds. Molecular hydrogen plays a significant role in establishing the thermal balance in many astrophysical environments and can be important as a spectral diagnostic of the gas. Modeling and interpretation of observations of such environments requires a quantitatively complete and accurate treatment of H2. Using this micro-physical model of H2, illustrative calculations of prototypical astrophysical environments are presented. This work forms the foundation for future investigations of these and other environments where H2 is an important constituent.

Ever since the pioneering study of Spitzer, it has been widely recognized that grains play an important role in the heating and cooling of photo-ionized environments. A detailed code is necessary to model grains in a photo-ionized medium. In this paper we will use the spectral synthesis code Cloudy for this purpose. The grain model of this code has recently undergone a comprehensive upgrade. Part of this is the newly developed hybrid grain charge model, which will be described in detail. This model allows discrete charge states of very small grains to be modelled accurately while simultaneously avoiding the overhead of fully resolving the charge distribution of large grains, thus making the model both accurate and computationally efficient. A comprehensive comparison with the fully resolved charge state models of Weingartner & Draine (2001) shows that the agreement is very satisfactory for realistic size distributions. The effect of the grain size distribution on the line emission from photo-ionized regions is studied by taking standard models for an H II region and a planetary nebula and adding a dust component to the models with varying grain size distributions. A comparison of the models shows that varying the size distribution has a dramatic effect on the emitted spectrum, and affects the ionization balance as well. These results clearly demonstrate that the grain size distribution is an important parameter in photo-ionization models.

Observations of infrared emission from AGN show that grains exist over a broad range of distances from the central object, extending to the point where they are destroyed by sublimation. These ~1000 K grains produce much of the 1 micron continuum. In this region closest to the central object there must be a gaseous component associated with the hot grains. This paper employs a state of the art grain model and shows that the gas must be very hot, with temperatures in the neighborhood of 10^6 K. The dusty component has a covering factor of roughly 50% and so this region also reprocesses much of the total x-ray emission. Our explicit models of the IR through x-ray spectral energy distributions allow the x-ray component to be predicted from IR observations. We are creating a grid of such predictions and will make them available as an XSPEC add-in, allowing this spectral component to be included in quantitative modeling of AGN spectra.

In this paper we present results obtained with the new grain code in Cloudy which underline the strong effect of photo-electric heating by grains in photo-ionized regions. We study the effect that the distribution of grain sizes has on the magnitude of the effect, and show that this effect is nothing short of dramatic. This makes the grain size distribution an important parameter in modeling of photo-ionized regions such as H II regions and planetary nebulae.

We present the basic properties of the multipolar post-AGB star IRAS 16594-4656, and discuss in particular its near infrared spectrum which shows shock excited H2 and [Fe II] emission lines.

In this paper we study the near-infrared emission spectrum of IRAS 16594-4656, a bipolar post-AGB star with spectral type B7 and no observed ionization. Using optical and near-infrared photometry we determined the total extinction towards this object to be A_V=7.5 +/- 0.4 mag and derived a distance of 2.2 +/- 0.4 kpc, assuming a luminosity of 10^4 L_sol. The near-infrared spectrum shows strong H_2 emission lines and some typical metastable shock excited lines such as [Fe II] 1.257 & 1.644 micron. We determined the rotational and vibrational excitation temperatures, as well as the ortho-to-para ratio of the molecular hydrogen. Based on these we argue that the H_2 emission is mainly collisionally excited. Line ratios indicate that the H_2 emission originates in a +/- 25 km/s C-type shock. On the other hand, the metastable lines, and especially the [Fe II] emission lines, indicate the presence of a +/- 75 km/s J-type shock. Hence we postulate that the H_2 emission originates where the stellar wind (with an observed terminal velocity of +/- 126 km/s) is funneled through an equatorial density enhancement, impinging almost tangentially upon the circumstellar material. The [Fe II] emission either occurs along the walls of the bipolar lobes where the transverse shock velocity would be higher, or could originate much closer to the central star in shocks in the post-AGB wind itself, or possibly even an accretion disk. Further high resolution near-infrared spectra are currently being obtained to confirm the proposed geometry and kinematics.

The identification of spectral lines can be a tedious process requiring the interrogation of large spectroscopic databases, but it does lend itself to software algorithms that can determine the characteristics of candidate line identifications. We present here criteria used for the identification of lines and a logic developed for a line identification software package called EMILI, which uses the v2.04 Atomic Line List as the basic line database. EMILI considers all possible database transitions within the wavelength uncertainties for observed optical emission lines and computes an approximate intensity for each candidate line. It searches for other multiplet members that are expected to be seen with each candidate line, and rank orders all of the tentative line identifications for each observed line based on a set of criteria. When applied to the spectra of the Orion Nebula and the planetary nebula IC 418, EMILI's recommended line ID's agree well with those of previous traditional manual line assignments. The existence of a semi-automated procedure should give impetus to the study of very high signal-to-noise spectra, enabling the identification of previously unidentified spectral lines to be handled with ease and consistency.

Ever since the pioneering study of Spitzer, it has been widely recognized that grains play an important role in the heating and cooling of photo-ionized environments. This includes the diffuse ISM, as well as H II regions, planetary nebulae, and photo-dissociation regions. A detailed code is necessary to model grains in a photo-ionized medium since the interactions of grains with their environment include a host of microphysical processes, and their importance can only be judged by performing a complete simulation. In this paper we will use the spectral synthesis code Cloudy for this purpose. A comprehensive upgrade of the grain model has been recently incorporated in Cloudy, and certain aspects of this upgrade will be discussed. Special emphasis will be on the new grain charge model. We will consider in detail the physics of grains in both ionized and neutral environments, and will present a calculation of photo-electric heating and collisional cooling rates for a range of physical conditions and grain materials and for a range of grain sizes (including a realistic size distribution). We conclude with a brief discussion of the problems currently hampering progress in this field. The new grain model will be used to model the silicate emission in the Ney-Allen nebula, and will help us better understand the nature of the grains in that part of the Orion complex.

We have investigated the infrared properties of a sample of objects with IRAS colors typical of planetary nebulae. The selected objects have not yet evolved to the planetary nebula stage since they have no detectable radio continuum emission and they are therefore likely post-AGB stars.

We investigate the infrared properties of the IRAS sources not confirmed as PNe in Van de Steene & Pottasch (1993). We observed 28 objects in the N-band of which 20 were detected and 5 were resolved. We obtained medium resolution Br_gamma spectra and we took high resolution JHKL images of these 20 objects. We critically assessed the identification of the IRAS counterpart in the images. Of sixteen positively identified objects, seven show Br_gamma in absorption, another six objects show Br_gamma in emission on top of a continuum. In another three objects, no clear Br_gamma absorption or emission was visible. Our detection rate of emission line objects is higher compared to previous studies. None of the objects with Br_gamma in emission were detected above a detection limit of 0.6 mJy/beam at 6 cm, while they should have been easily detected if the radio flux was optically thin and Case B recombination was applicable. Their SED's are typical for post-AGB stars. We also constructed various color-color diagrams using the near-infrared and IRAS magnitudes. No distinction can be made between the objects showing Br_gamma in emission, absorption, or a flat spectrum in the near and far-infrared color-color diagrams. The near-infrared color-color diagrams show evidence for a very large range of extinction, which in part is of circumstellar origin. Near-infrared versus far-infrared color-color diagrams show trends that are consistent with the expected evolution of the circumstellar shell. This sample of post-AGB stars show a larger range in color and are generally closer to the galactic plane than the ones known so far.

This paper describes the most important differences between Cloudy 90 and Cloudy 94, the most recent upgrade of the code which has now been posted on the WWW. It also describes the upgrade of the grain model in Cloudy which is currently underway.

When measuring diameters of partially resolved sources like planetary nebulae, H II regions or galaxies, often a technique called gaussian deconvolution is used. This technique yields a gaussian diameter which subsequently has to be multiplied with a conversion factor to obtain the true angular diameter of the source. This conversion factor is a function of the Full Width at Half Maximum (FWHM) of the beam. It also depends on the intrinsic surface brightness distribution of the source. In this paper the theory behind this technique will be studied. This study will be restricted to circularly symmetric geometries and beams. First an implicit equation will be derived, from which the conversion factor for a given surface brightness distribution and beam size can be solved. Explicit expressions for the conversion factor will be derived from this equation which are valid in cases where the beam size is larger than the intrinsic size of the source. A more detailed discussion will be given for two simple geometries: a circular constant surface brightness disk and a spherical constant emissivity shell with arbitrary inner radius. The theory is subsequently used to construct a new technique for determining the FWHM of an arbitrary observed surface brightness distribution. Usually the FWHM of the source and beam are measured using gaussian fits, but second moments can also be used. The alternative use of second moments in this context is studied here for the first time and it is found that in this case the conversion factor has a different value which is independent of the beam size. In the limit for infinitely large beam sizes, the values of the conversion factors for both techniques are equal.

We present H alpha emission line spectra of 7 post-AGB stars. Typically, they have a P-Cygni type profile with a very strong emission component and a relatively narrow absorption component. These lines are formed in the fast post-AGB wind which in most cases has a velocity of around 100 km/s. Hence the emission originates close to the central star and the absorption occurs in the fast wind itself or in the region where the fast wind sweeps up the slow wind. The emission profile has very broad wings. It is still unclear what causes these wings.

The calculation of accurate collision strengths for atomic transitions has been a long standing problem in quantitative spectroscopy. Most modern calculations are based on the R-matrix method and problems pertaining to the use of this method have led to a discussion of the accuracy of these results. More in particular, based on an analysis of the spectra of NGC 3918 and NGC 6302, Clegg et al. (1987) and Oliva et al. (1996) have questioned R-matrix calculations for the infrared [Ne V] fine-structure transitions. Using improved flux measurements for the [Ne V] lines, we show that the conclusion that these collision strengths would be too high, is not correct. The discrepancies found by Clegg et al. (1987) can be explained by the inaccuracy of the [Ne V] 342.6 nm flux they adopted. The discrepancies found by Oliva et al. (1996) can be explained by the inaccuracy of the LRS flux for the [Ne V] 14.32 micron line. Based on the data presented in this paper there is no reason to assume that there are any problems with the R-matrix calculations for Ne^4+ of Lennon & Burke (1994). We show that the data are accurate at the 30 % level or better. This confirms the validity of the close coupling method.

When measuring diameters of partially resolved sources often a technique called gaussian deconvolution is used. This technique yields a gaussian diameter which subsequently has to be multiplied with a conversion factor to obtain the true angular diameter of the source. This conversion factor is a function of the FWHM of the beam or point spread function and also depends on the intrinsic surface brightness distribution of the source. In this paper conversion factors are presented for a number of simple geometries: a circular constant surface brightness disk and a spherical constant emissivity shell, using a range of values for the inner radius. Also more realistic geometries are studied, based on a spherically symmetric photo-ionization model of a planetary nebula. This enables a study of optical depth effects, a comparison between images in various emission lines and the use of power law density distributions. It is found that the conversion factor depends quite critically on the intrinsic surface brightness distribution, which is usually unknown. The uncertainty is particularly large if extended regions of low surface brightness are present in the nebula. In such cases the use of gaussian or second moment deconvolution is not recommended. As an alternative, a new algorithm is presented which allows the determination of the intrinsic FWHM of the source using only the observed surface brightness distribution and the FWHM of the beam. Tests show that this implicit deconvolution method works well in realistic conditions, even when the signal-to-noise is low, provided that the beam size is less than roughly 2/3 of the observed FWHM and the beam profile can be approximated by a gaussian.

In this paper we study new infrared spectra of the evolved planetary nebula NGC 6445 obtained with ISO. These data show that the thermal emission from the grains is very cool and has a low flux compared to H beta. A model of the ionized region is constructed, using the photo-ionization code CLOUDY 90.05. Based on this model, we show from depletions in the gas phase elements that little grain destruction can have occurred in the ionized region of NGC 6445. We also argue that dust-gas separation in the nebula is not plausible. The most likely conclusion is that grains are residing inside the ionized region of NGC 6445 and that the low temperature and flux of the grain emission are caused by the low luminosity of the central star and the low optical depth of the grains. This implies that the bulk of the silicon-bearing grains in this nebula were able to survive exposure to hard UV photons for at least several thousands of years, contradicting previously published results. A comparison between optical and infrared diagnostic line ratios gives a marginal indication for the presence of a t^2-effect in the nebula. However, the evidence is not convincing and other explanations for the discrepancy are also plausible. The off-source spectrum taken with ISO-LWS clearly shows the presence of a warm cirrus component with a temperature of 24 K as well as a very cold component with a temperature of 7 K. Since our observation encompasses only a small region of the sky, it is not clear how extended the 7 K component is and whether it contributed significantly to the FIRAS spectrum taken by COBE. Because our line of sight is in the galactic plane, the very cold component could be a starless core.

Infrared spectra of NGC 7027 and NGC 6543 ranging from 2.4 to 45 micron were obtained with the Short Wavelength Spectrometer on board the Infrared Space Observatory. A first analysis of these spectra, with the aid of photo-ionization models, is presented. We report the first detection of the Ar VI 4.53 micron and Ne VI 7.65 micron lines in the spectrum of NGC 7027. When compared with older observations it is clear that the Ar VI line and possibly also other lines have increased in strength since 1981. We argue that a likely explanation for this variability is a change in the spectral energy distribution of the central star, possibly an increase in effective temperature. However, this result needs to be confirmed by further observations. We also report a non-detection of the O IV 25.9 micron line and the first detection of the Na III 7.32 micron line in the spectrum of NGC 6543. The non-detection is not expected based on a blackbody approximation for the spectrum of the central star. The ionization threshold for O(3+) is just beyond the He II limit, and the absence of this line shows that the stellar flux drops at least by a factor 350 at the He II limit. Modeling the O IV line may prove to be a valuable test for atmosphere models.

A parameter study of the spectral evolution of a typical post-AGB star, with particular emphasis on the evolution of the infrared colours, is presented. The models are based on the latest evolutionary tracks for hydrogen burning post-AGB stars. For such tracks the evolutionary rate is very dependent on the assumed mass loss rate as a function of time. We investigate this effect by modifying the mass loss prescription. The newly calculated evolutionary rates and density distributions are used to model the spectral evolution of a post-AGB star with the photo-ionization code Cloudy, including dust in the radiative transfer. Different assumptions for the dust properties and dust formation are considered. It is shown that by varying these parameters in a reasonable way, entirely different paths are followed in the IRAS colour-colour diagram. First of all, the effects of the evolution of the central star on the expanding dust shell can not be neglected. Also the dust properties and the definition of the end of the AGB phase have an important effect. The model tracks show that objects occupying the same location in the IRAS colour-colour diagram can have a different evolutionary past, and therefore the position in the IRAS colour-colour diagram alone can not a priori give a unique determination of the evolutionary status of an object. An alternative colour-colour diagram, the K-[12] vs. [12]-[25] diagram, is presented. The tracks in this diagram seem less affected by particulars of the grain emission. This diagram may be a valuable additional tool for studying post-AGB evolution.

We have constructed photo-ionization models of five galactic bulge planetary nebulae using our automatic method which enables a fully self-consistent determination of the physical parameters of a planetary nebula. The models are constrained using the spectrum, the IRAS and radio fluxes and the angular diameter of the nebula. We also conducted a literature search for physical parameters determined with classical methods for these nebulae. Comparison of the distance independent physical parameters with published data shows that the stellar temperatures generally are in good agreement and can be considered reliable. The literature data for the electron temperature, electron density and also for the abundances show a large spread, indicating that the use of line diagnostics is not reliable and that the accuracy of these methods needs to be improved. Comparison of the various abundance determinations indicates that the uncertainty in the electron temperature is the main source of uncertainty in the abundance determination. The stellar magnitudes predicted by the photo-ionization models are in good agreement with observed values.

Radio continuum emission has been detected from a planetary nebula in the Large Magellanic Cloud: this is the first radio continuum detection for any planetary nebula outside our galaxy. The radio flux density is a factor of two lower than predicted from the \hbeta\ flux. This could be due either to a two-component morphology or to a stellar contribution to the \hbeta\ emission. We have modeled the optical and infrared spectrum using the photo-ionisation code Cloudy: a very good fit is obtained if we assume the nebula is carbon rich. The derived diameter implies an evolutionary age of the nebula of $<1000\yr$, similar to what is derived from the IRAS colours alone. The central star, which is of [WC] type, has a much higher stellar temperature than Galactic [WC] stars showing similar circumstellar IRAS colours. An explanation could be that the expansion velocity of the nebula is lower than those of its galactic counterparts. This radio detection indicates that accurate nebular luminosities could be determined from their radio emission for many LMC planetary nebulae.