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.

We present a new census of Galactic and extragalactic symbiotic stars (SySts). This compilation contains 323 known and 87 candidate SySts. Of the confirmed SySts, 257 are Galactic and 66 extragalactic. The spectral energy distributions (SEDs) of 348 sources have been constructed using 2MASS and AllWISE data. Regarding the Galactic SySts, 74% are S-types, 13% D and 3.5% D$^{\prime}$. S-types show an SED peak between 0.8 and 1.7 $\mu$m, whereas D-type show a peak at longer wavelengths between 2 and 4 $\mu$m. D$^{\prime}$-type, on the other hand, display a nearly flat profile. Gaia distances and effective temperatures are also presented. According to their Gaia distances, S-type are found to be members of both thin and thick Galactic disk populations, while S+IR- and D-types are mainly thin disk sources. Gaia temperatures show a reasonable agreement with the temperatures derived from SEDs within their uncertainties. A new census of the OVI $\lambda$6830 Raman-scattered line in SySts is also presented. From a sample of 298 SySts with available optical spectra, 55% are found to emit the line. No significant preference is found among the different types. The report of the OVI $\lambda$6830 Raman-scattered line in non-SySts is also discussed as well as the correlation between the Raman-scattered OVI line and X-ray emission. We conclude that the presence of the OVI Raman-scattered line still provides a strong criterion for identifying a source as a SySt.

In this second paper in a series of papers based on the most-up-to-date catalogue of symbiotic stars (SySts), we present a new approach for identifying and distinguishing SySts from other Halpha emitters in photometric surveys using machine learning algorithms such as classification tree, linear discriminant analysis, and K-nearest neighbour. The motivation behind of this work is to seek for possible colour indices in the regime of near- and mid-infrared covered by the 2MASS and WISE surveys. A number of diagnostic colour-colour diagrams are generated for all the known Galactic SySts and several classes of stellar objects that mimic SySts such as planetary nebulae, post-AGB, Mira, single K and M giants, cataclysmic variables, Be, AeBe, YSO, weak and classical T Tauri stars, and Wolf-Rayet. The classification tree algorithm unveils that primarily J-H, W1-W4 and Ks-W3 and secondarily H-W2, W1-W2 and W3-W4 are ideal colour indices to identify SySts. Linear discriminant analysis method is also applied to determine the linear combination of 2MASS and AllWISE magnitudes that better distinguish SySts. The probability of a source being a SySt is determined using the K-nearest neighbour method on the LDA components. By applying our classification tree model to the list of candidate SySts (Paper I), the IPHAS list of candidate SySts, and the DR2 VPHAS+ catalogue, we find 125 (72 new candidates) sources that pass our criteria while we also recover 90 per cent of the known Galactic SySts.

The spatial structure of the emission lines and continuum over the 50 arcsecond extent of the nearby, O-rich, planetary nebula NGC 7009 (Saturn Nebula) have been observed with the MUSE integral field spectrograph on the ESO Very Large Telescope. Science Verification data, in <0.6 arcsecond seeing, have been reduced and analysed as images over the wavelength range 4750-9350A. Emission line maps over the bright shells are presented, from neutral to the highest ionization available (He II and [Mn V]). For collisionally excited lines (CELs), maps of electron temperature (T_e from [N II] and [S III]) and electron density (N_e from [S II] and [Cl III]) are available and for optical recombination lines (ORLs) temperature (from the Paschen jump and ratio of He I lines) and density (from high Paschen lines). These estimates are compared: for the first time, maps of the differences in CEL and ORL T_e's have been derived, and correspondingly a map of t^2 between a CEL and ORL temperature, showing considerable detail. Total abundances of He and O were formed, the latter using three ionization correction factors. However the map of He/H is not flat, departing by ~2% from a constant value, with remnants corresponding to ionization structures. Ionization correction factor methods are compared for O abundance, but none delivers a flat map. An integrated spectrum over an area of 2340 square arcseconds was also formed and compared to 1D photoionization models. The spatial variation of a range of nebular parameters illustrates the complexity of the ionized media in NGC 7009. These MUSE data are very rich with detections of many lines over areas of hundreds of square arcseconds and follow-on studies are indicated. (Abridged)

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 present a new updated catalogue of Galactic and extragalactic symbiotic stars (SySts). Since the last catalogue of SySts (Belczynski et al. 2000), the number of known SySts has significantly increased. Our new catalogue contains 316 known and 82 candidates SySts. Of the confirmed Systs 252 are located in our Galaxy and 64 in nearby galaxies. This reflects an increase of $\sim$50% in the population of Galactic SySts and $\sim$400% in the population of extragalactic SySts. The spectral energy distribution (SED) of 334 (known and candidates) SySts have been constructed using the 2MASS and WISE data. These SEDs are used to provide a robust reclassification in scheme of S- (74%), D- (15%) and D'-types (2.5%). The SEDs of S- and D-type peak between 0.8 and 1.6$\mu$m and between 1.6 and 4$\mu$m, respectively, whereas those of D'-type exhibit a plateau profile. Moreover, we provide the first compilation of SySts that exhibit the OVI Raman-scattered line at 6830\AA. Our analysis shows that 55% of the Galactic SySts exhibit that line in their spectrum, whereas this percentage is different from galaxy to galaxy.

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 have mapped the CO(3-2) line emission around the Mira AB system at 0.5 resolution using the Atacama Large Millimeter/submillimeter Array (ALMA). The CO map shows amazing complexity. The circumstellar gas has been shaped by different dynamical actors during the evolution of the system and several morphological components can be identified. The companion is marginally resolved in continuum emission and is currently at 0.487$\pm$0.006 separation. In the main line component, centered on the stellar velocity, spiral arcs around Mira A are found. The spiral appears to be relatively flat and oriented in the orbital plane. An accretion wake behind the companion is clearly visible and the projected arc separation is of order 5''. In the blue wing of the line emission, offset from the main line, several large ($\sim$5-10''), opposing arcs are found. We tentatively suggest that this structure is created by the wind of Mira B blowing a bubble in the expanding envelope of Mira A.

Detached shells are believed to be created during a thermal pulse, and constrain the time scales and physical properties of one of the main drivers of late stellar evolution. We aim at determining the morphology of the detached dust shells around the carbon AGB stars R Scl and V644 Sco, and compare this to observations of the detached gas shells. We observe the polarised, dust-scattered stellar light around these stars using the PolCor instrument mounted on the ESO 3.6m telescope. Observations were done with a coronographic mask to block out the direct stellar light. The polarised images clearly show the detached shells. Using a dust radiative transfer code to model the dust-scattered polarised light, we constrain the radii and widths of the shells to 19.5 arcsec and 9.4 arcsec for the detached dust shells around R Scl and V644 Sco, respectively. Both shells have an overall spherical symmetry and widths of approx. 2 arcsec. For R Scl we can compare the observed dust emission directly with high spatial-resolution maps of CO(3-2) emission from the shell observed with ALMA. We find that the dust and gas coincide almost exactly, indicating a common evolution. The data presented here for R Scl are the most detailed observations of the entire dusty detached shell to date. For V644 Sco these are the first direct measurements of the detached shell. Also here we find that the dust most likely coincides with the gas shell. The observations are consistent with a scenario where the detached shells are created during a thermal pulse. The determined radii and widths will constrain hydrodynamical models describing the pre-pulse mass loss, the thermal pulse, and post-pulse evolution of the star.

Optical integral-field spectroscopy was used to investigate the planetary nebula NGC 3242. We analysed the main morphological components of this source, including its knots, but not the halo. In addition to revealing the properties ofthe physical and chemical nature of this nebula, we also provided reliable spatially resolved constraints that can be used for future photoionisation modelling of the nebula. The latter is ultimately necessary to obtain a fully self-consistent 3D picture of the physical and chemical properties of the object. The observations were obtained with the VIMOS instrument attached to VLT-UT3. Maps and values for specific morphological zones for the detected emission-lines were obtained and analysed with routines developed by the authors to derive physical and chemical conditions of the ionised gas in a 2D fashion. We obtained spatially resolved maps and mean values of the electron densities, temperatures, and chemical abundances, for specific morphological structures in NGC 3242. These results show the pixel-to-pixel variations of the the small- and large-scale structures of the source. These diagnostic maps provide information free from the biases introduced by traditional single long-slit observations. In general, our results are consistent with a uniform abundance distribution for the object, whether we look at abundance maps or integrated fluxes from specified morphological structures. The results indicate that special care should be taken with the calibration of the data and that only data with extremely good signal-to-noise ratio and spectral coverage should be used to ensure the detection of possible spatial variations.

We aim to detect the magnetic field and infer its properties around four AGB stars using H$_2$O maser observations. The sample we observed consists of the following sources: the semi-regular variable RT Vir and the Mira variables AP Lyn, IK Tau, and IRC+60370. We observed the 6$_{1,6}-5_{2,3}$ H$_2$O maser rotational transition, in full-polarization mode, to determine its linear and circular polarization. Based on the Zeeman effect, one can infer the properties of the magnetic field from the maser polarization analysis. We detected a total of 238 maser features, in three of the four observed sources. No masers were found toward AP Lyn. The observed masers are all located between 2.4 and 53.0 AU from the stars. Linear and circular polarization was found in 18 and 11 maser features, respectively. We more than doubled the number of AGB stars in which magnetic field has been detected from H$_2$O maser polarization, as our results confirm the presence of fields around IK Tau, RT Vir and IRC+60370. The strength of the field along the line of sight is found to be between 47 and 331 mG in the H$_2$O maser region. Extrapolating this result to the surface of the stars, assuming a toroidal field ($\propto$ r$^{-1}$), we find magnetic fields of 0.3-6.9 G on the stellar surfaces. If, instead of a toroidal field, we assume a poloidal field ($\propto$ r$^{-2}$), then the extrapolated magnetic field strength on the stellar surfaces are in the range between 2.2 and $\sim$115 G. Finally, if a dipole field ($\propto$ r$^{-3}$) is assumed, the field strength on the surface of the star is found to be between 15.8 and $\sim$1945 G. The magnetic energy of our sources is higher than the thermal and kinetic energy in the H$_2$O maser region of this class of objects. This leads us to conclude that, indeed, magnetic fields probably play an important role in shaping the outflows of evolved stars. (abridged)

Most of PNe are not spherical. The loss of spherical symmetry occurs somewhere between the AGB and PN phase. The cause of this change of morphology is not yet well known, but magnetic fields are one of the possible agents. Its origin remains to be determined, and potentially requires the presence of a massive companion to the AGB star. Therefore, further detections of the magnetic field around evolved stars (in particular those thought to be part of a binary system) are crucial to improve our understanding of the origin and role of magnetism on evolved stars. One such binaries is the pre-PN OH231.8, around which a magnetic field was detected in the OH maser region of the outer circumstellar envelope. We aim to detect and infer the properties of the magnetic field of this source in the water maser region. We observed the 6_{1,6}-5_{2,3} water maser rotational transition to determine its linear and circular polarization. These emissions are located within the inner regions of OH231.8 (at few tens of AU). We detected 30 water maser features, which occur in two distinct regions that are moving apart with a velocity on the sky of 2.3 mas/year. Taking into account the inclination angle of the source, this corresponds to an average separation velocity of 21 km/s. Based on the velocity gradient of the maser emission, the masers appear to be dragged along the direction of the nebula jet. Linear polarization is present in 3 of the features, and circular polarization was detected in the 2 brightest ones. We found that the strength of the magnetic field is |B_{||}|~45 mG which, when assuming a toroidal magnetic field, implies B~2.5 G on the stellar surface. The morphology of the field is not yet determined, but the high scatter found on the directions of the linear polarization vectors could indicate that the masers occur near the tangent points of a toroidal field.

In this paper we present an analysis of the physical and chemical conditions of the planetary nebula NGC 40 through spatially-resolved spectroscopic maps. We also introduce a new algorithm --2D_NEB-- based on the well-established IRAF nebular package, which was developed to enable the use of the spectroscopic maps to easily estimate the astrophysical quantities of ionised nebulae. The 2D_NEB was benchmarked, and we clearly show that it works properly, since it compares nicely with the IRAF nebular software. Using this software, we derive the maps of several physical parameters of NGC 40. From these maps, we conclude that Te[NII] shows only a slight temperature variation from region to region, with its values constrained between ~8,000 K and ~9,500 K. Electron densities, on the other hand, have a much more prominent spatial variation, as Ne[SII] values vary from ~1,000 cm^(-3) to ~3,000 cm^(-3). Maps of the chemical abundances also show significant variations. From the big picture of our work, we strongly suggest that analysis with spatial resolution be mandatory for more complete study of the physical and chemical properties of planetary nebulae.