Despite the many studies in the last decades, the low-ionisation structures (LISs) of planetary nebulae (PNe) still hold several mysteries. Recent imaging surveys have demonstrated that LISs are composed of molecular gas. Here we report H$_2$ emission in the LISs of NGC 7009 and NGC 6543 by means of very deep narrow-band H$_2$ images taken with NIRI@Gemini. The surface brightness of the H2 1-0 S(1) line is estimated to be (0.46-2.9)$\times$10$^{-4}$ erg s$^{-1}$ cm$^{-2}$ sr$^{-1}$ in NGC 7009 and (0.29-0.48)$\times$10$^{-4}$ erg s$^{-1}$ cm$^{-2}$ sr$^{-1}$ in NGC 6543, with signal-to-noise ratios of 10-42 and 3-4, respectively. These findings provide further confirmation of hidden H$_2$ gas in LISs. The emission is discussed in terms of the recent proposed diagnostic diagram R(H$_2$)=H$_2$ 1-0 S(1)/H$_2$ 2-1 S(1) versus R(Br$\gamma$)=H$_2$ 1-0 S(1)/Br$\gamma$, which was suggested to trace the mechanism responsible for the H$_2$ excitation. Comparing our observations to shock and ultraviolet (UV) molecular excitation models, as well as a number of observations compiled from the literature showed that we cannot conclude for either UV or shocks as the mechanism behind the molecular emission.

We present a new approach to study planetary nebulae using integral field spectroscopy. VLT@VIMOS datacube of the planetary nebula Abell 14 is analysed in three different ways by extracting: (i) the integrated spectrum, (ii) 1-dimensional simulated long slit spectra for different position angles and (iii) spaxel-by-spaxel spectra. These data are used to built emission-line diagnostic diagrams and explore the ionization structure and excitation mechanisms combining data from 1- and 3- dimensional photoionization models. The integrated and 1D simulated spectra are suitable for developing diagnostic diagrams, while the spaxel spectra can lead to misinterpretation of the observations. We find that the emission-line ratios of Abell 14 are consistent with UV photo-ionized emission, however there are some pieces of evidence of an additional thermal mechanism. The chemical abundances confirm its previous classification as a Type I planetary nebula, without spatial variation. We find, though, variation in the ionization correction factors (ICFs) as a function of the slit position angle. The star at the geometric centre of Abell 14 has an A5 spectral type with an effective temperature of Teff = 7909$\pm$135 K and surface gravity log(g) = 1.4$\pm$0.1 cm s$^{-2}$. Hence, this star cannot be responsible for the ionization state of the nebula. Gaia parallaxes of this star yield distances between 3.6 and 4.5 kpc in good agreement with the distance derived from a 3-dimensional photoionization modelling of Abell 14, indicating the presence of a binary system at the centre of the planetary nebula.

Since the start of ALMA observatory operation, new and important chemistry of infrared cold core was revealed. Molecular transitions at millimeter range are being used to identify and to characterize these sources. We have investigated the 231 GHz ALMA archive observations of the infrared dark cloud region C9, focusing on the brighter source that we called as IRDC-C9 Main. We report the existence of two sub-structures on the continuum map of this source: a compact bright spot with high chemistry diversity that we labelled as core, and a weaker and extended one, that we labelled as tail. In the core, we have identified lines of the molecules OCS(19-18), $^{13}$CS(5-4) and CH$_{3}$CH$_{2}$CN, several lines of CH$_{3}$CHO and the k-ladder emission of $^{13}$CH$_{3}$CN.We report two different temperature regions: while the rotation diagram of CH$_{3}$CHO indicates a temperature of 25 K, the rotation diagram of $^{13}$CH$_{3}$CN indicates a warmer phase at temperature of $\sim450$K. In the tail, only the OCS(19-18) and $^{13}$CS(5-4) lines were detected. We used the $Nautilus$ and the \textsc{Radex} codes to estimate the column densities and the abundances. The existence of hot gas in the core of IRDC-C9 Main suggests the presence of a protostar, which is not present in the tail.

We present a detailed analysis of deep VLT/X-Shooter observations of the planetary nebula Tc 1. We calculate gas temperature, density, extinction, and abundances for several species from the empirical analysis of the total line fluxes. In addition, a spatially resolved analysis of the most intense lines provides the distribution of such quantities across the nebula. The new data reveal that several lines exhibit a double peak spectral profile consistent with the blue- and red-shifted components of an expanding spherical shell. The study of such components allowed us to construct for the first time a three-dimensional morphological model, which reveals that Tc 1 is a slightly elongated spheroid with an equatorial density enhancement seen almost pole on. A few bright lines present extended wings (with velocities up to a few hundred km/s), but the mechanism producing them is not clear. We constructed photoionization models for the main shell of Tc 1. The models predict the central star temperature and luminosity, as well as the nebular density and abundances similar to previous studies. Our models indicate that Tc 1 is located at a distance of approximately 2 kpc. We report the first detection of the [Kr III] 6825 A emission line, from which we determine the Krypton abundance. Our model indicates that the main shell of Tc 1 is matter bounded; leaking H ionizing photons may explain the ionization of its faint AGB-remnant halo.

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 performed a detailed spectroscopic analysis of the fullerene C60-containing planetary nebula (PN) Lin49 in the Small Magellanic Cloud using XSHOOTER at the ESO VLT and the Spitzer/IRS instruments. We derived nebular abundances for nine elements. We used TLUSTY to derive photospheric parameters for the central star. Lin49 is C-rich and metal-deficient PN (Z~0.0006). The nebular abundances are in good agreement with Asymptotic Giant Branch nucleosynthesis models for stars with initial mass 1.25 Msun and metallicity Z = 0.001. Using the TLUSTY synthetic spectrum of the central star to define the heating and ionising source, we constructed the photoionisation model with CLOUDY that matches the observed spectral energy distribution (SED) and the line fluxes in the UV to far-IR wavelength ranges simultaneously. We could not fit the ~1-5 um SED using a model with 0.005-0.1 um-sized graphite grains and a constant hydrogen density shell owing to the prominent near-IR excess, while at other wavelengths the model fits the observed values reasonably well. We argue that the near-IR excess might indicate either (1) the presence of very small particles in the form of small carbon clusters, small graphite sheets, or fullerene precursors, or (2) the presence of a high-density structure surrounding the central star. We found that SMC C60 PNe show a near-IR excess component to lesser or greater degree. This suggests that these C60 PNe might maintain a structure nearby their central star.

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 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 study the kinematics and excitation mechanisms of H2 and [Fe II] lines in a sample of 67 emission-line galaxies with Infrared Telescope Facility SpeX near-infrared (NIR, 0.8-2.4 micrometers) spectroscopy together with new photoionisation models. H2 emission lines are systematically narrower than narrow-line region (NLR) lines, suggesting that the two are, very likely, kinematically disconnected. The new models and emission-line ratios show that the thermal excitation plays an important role not only in active galactic nuclei but also in star forming galaxies. The importance of the thermal excitation in star forming galaxies may be associated with the presence of supernova remnants close to the region emitting H2 lines. This hypothesis is further supported by the similarity between the vibrational and rotational temperatures of H2. We confirm that the diagram involving the line ratios H2 2.121/Br_gamma and [Fe II] 1.257/Pa_beta is an efficient tool for separating emission-line objects according to their dominant type of activity. We suggest new limits to the line ratios in order to discriminate between the different types of nuclear activity.

Detections of molecular lines, mainly from H2$ and CO, reveal molecular material in planetary nebulae. Observations of a variety of molecules suggest that the molecular composition in these objects differs from that found in interstellar clouds or in circumstellar envelopes. The success of the models, which are mostly devoted to explain molecular densities in specific planetary nebulae, is still partial, however. The present study aims at identifying the influence of stellar and nebular properties on the molecular composition of planetary nebulae by means of chemical models. A comparison of theoretical results with those derived from the observations may provide clues to the conditions that favor the presence of a particular molecule. A self-consistent photoionization numerical code was adapted to simulate cold molecular regions beyond the ionized zone. The code was used to obtain a grid of models and the resulting column densities are compared with those inferred from observations. Our models show that the inclusion of an incident flux of X-rays is required to explain the molecular composition derived for planetary nebulae. We also obtain a more accurate relation for the N(CO)/N(H2) ratio in these objects. Molecular masses obtained by previous works in the literature were then recalculated, showing that these masses can be underestimated by up to three orders of magnitude. We conclude that the problem of the missing mass in planetary nebulae can be solved by a more accurate calculation of the molecular mass.

Molecular hydrogen emission is commonly observed in planetary nebulae. Images taken in infrared H2 emission lines show that at least part of the molecular emission is produced inside the ionised region. In the best-studied case, the Helix nebula, the H2 emission is produced inside cometary knots (CKs), comet-shaped structures believed to be clumps of dense neutral gas embedded within the ionised gas. Most of the H2 emission of the CKs seems to be produced in a thin layer between the ionised diffuse gas and the neutral material of the knot, in a mini photodissociation region (PDR). However, PDR models published so far cannot fully explain all the characteristics of the H2 emission of the CKs. In this work, we use the photoionisation code \textsc{Aangaba} to study the H2 emission of the CKs, particularly that produced in the interface H^+/H^0 of the knot, where a significant fraction of the H2 1-0S(1) emission seems to be produced. Our results show that the production of molecular hydrogen in such a region may explain several characteristics of the observed emission, particularly the high excitation temperature of the H2 infrared lines. We find that the temperature derived from H2 observations even of a single knot, will depend very strongly on the observed transitions, with much higher temperatures derived from excited levels. We also proposed that the separation between the H_alpha and NII peak emission observed in the images of CKs may be an effect of the distance of the knot from the star, since for knots farther from the central star the NII line is produced closer to the border of the CK than H_alpha.

The analysis and interpretation of the H2 line emission from planetary nebulae have been done in the literature assuming that the molecule survives only in regions where the hydrogen is neutral, as in photodissociation, neutral clumps or shocked regions. However, there is strong observational and theoretical evidence that at least part of the H2 emission is produced inside the ionized region of such objects. The aim of the present work is to calculate and analyze the infrared line emission of H2 produced inside the ionized region of planetary nebulae using a one-dimensional photoionization code. The photoionization code Aangaba was improved in order to calculate the statistical population of the H2 energy levels and the intensity of the H2 infrared emission lines in physical conditions typical of planetary nebulae. A grid of models was obtained and the results are analyzed and compared with the observational data. We show that the contribution of the ionized region to the H2 line emission can be important, particularly in the case of nebulae with high temperature central stars. This result explains why H2 emission is more frequently observed in bipolar planetary nebulae (Gatley's rule), since this kind of object typically has hotter stars. Collisional excitation plays an important role on the population of the rovibrational levels of the electronic ground state of H2. Radiative mechanisms are also important, particularly for the upper vibrational levels. Formation pumping can have minor effects on the line intensities produced by de-excitation from very high rotational levels, especially in dense and dusty environments. We included the effect of the H2 on the thermal equilibrium of the gas, concluding that H2 only contributes to the thermal equilibrium in the case of a very high temperature of the central star or a high dust-to-gas ratio, mainly through collisional de-excitation.

In the present work, we use a photoionisation code to study the H2 emission of the Helix nebula (NGC 7293) cometary knots, particularly that produced in the interface H+/H0 of the knot, where a significant fraction of the H2 1-0 S(1) emission seems to be produced. Our results show that the production of molecular hydrogen in such region may explain several characteristics of the observed emission, particularly the high excitation temperature of the H2 infrared lines.

This paper presents an analysis of the concentration of the hydrogen molecule inside the ionized region of planetary nebulae. The equations corresponding to the ionization and chemical equilibria of H, H+, H-, H2, H2+, and H3+ are coupled with the equations of ionization and thermal balance for a photoionized atomic gas. Forty different reactions related to the formation or the destruction of these species are included. The presence of dust is taken into account, since grains act as catalysts for the production of H2, as well as shield the molecules against the stellar ionizing radiation. We analyze the effect of the stellar ionizing continuum, as well as of the gas and grain properties on the calculated H2 mass. It is shown that a significant concentration of H2 can survive inside the ionized region of planetary nebulae, mostly in the inner region of the recombination zone. The total H2 to total hydrogen mass ratio inside the ionized region increases with the central star temperature, and, depending on the PN physical conditions, it can be of the order of 10^-6 or even higher. The increase of the recombination zone with the stellar temperature can account for such correlation. This can explain why the H2 emission is more frequently observed in bipolar planetary nebulae (Gatley's rule), since this kind of object has typically hotter stars. Applying our results for the planetary nebula NGC 6720, we obtain an H2 to hydrogen mass ratio similar to the value obtained from the observed H2 line emission.