The study of planetary nebulae started more than a century ago. Since then the understanding of these exciting objects has advanced extraordinarily. I present a personal selection of topics to reflect on some developments of PNe research.

Determining the demographics of the Galactic planetary nebula (PN) population is an important goal to further our understanding of this intriguing phase of stellar evolution. The Galactic population has more than doubled in number over the last 15 years, particularly from narrowband Hα surveys along the plane. In this review I will summarise these results, with emphasis on the time interval since the last IAU Symposium. These primarily optical surveys are not without their limitations and new surveys for PNe in the infrared similarly face a number of challenges. I will discuss the need for multi-wavelength approaches to discovery and analysis. The desire to have accurate volume-limited samples of Galactic PNe at our disposal is emphasised, which will be impacted with new data from the Gaia satellite mission. We need robust surveys of PNe and their central stars, especially volume-limited surveys, in order to clarify and quantify their evolutionary pathways.

This review presents the latest advances in the nebular studies of post-AGB objects. Post-AGB stars are great tools to test nucleosynthesis and evolution models for stars of low and intermediate masses, and the evolution of dust in harsh environment. I will present the newly discovered class of post-RGB stars, formed via binary interaction on the RGB. Binary systems can also lead to the formation of two class of aspherical post-AGB, the Proto-Planetary Nebulae and the naked post-AGBs (dusty RV Taus, a.k.a. Van Winckel’s stars).

We present infrared two color diagrams (2CDs) for a large sample of AGB stars, post-AGB stars, and planetary nebulae (PNe) to study their spectral evolution in infrared bands. We discuss possible evolutionary tracks from AGB stars to PNe on the 2CDs. We use catalogs from the available literature for the sample of AGB stars, post-AGB stars, and PNe in our Galaxy.

UWISH2 is the first unbiassed imaging survey of v =1 − 0 S(1) molecular hydrogen emission (λ = 2.122~μm) in the northern Galactic Plane. Here we discuss 284 extended emission line objects which we consider to be candidate planetary or pre-planetary nebulae. Some are clearly associated with known PN, but the majority (60%) have no previous detection. We have classified the objects according to morphology and find 53% are bipolar with half of these being new detections. The remaining objects are mostly elliptical/round (35%), fainter than the median sample flux (4.4 × 10−17 W m−2) and previously undetected (74%).

The HASH (Hong Kong/ AAO/ Strasbourg/ Hα) planetary nebula research platform is a unique data repository with a graphical interface and SQL capability that offers the community powerful, new ways to undertake Galactic PN studies. HASH currently contains multi-wavelength images, spectra, positions, sizes, morphologies and other data whenever available for 2401 true, 447 likely, and 692 possible Galactic PNe, for a total of 3540 objects. An additional 620 Galactic post-AGB stars, pre-PNe, and PPN candidates are included. All objects were classified and evaluated following the precepts and procedures established and developed by our group over the last 15 years. The complete database contains over 6,700 Galactic objects including the many mimics and related phenomena previously mistaken or confused with PNe. Curation and updating currently occurs on a weekly basis to keep the repository as up to date as possible until the official release of HASH v1 planned in the near future.

The understanding of astronomical nebulae is based on observational data (images, spectra, 3D data-cubes) and theoretical models. In this review, I present my very biased view on photoionization modeling of planetary nebulae, focusing on 1D multi-component models, on 3D models and on big database of models.

Nearly 50 years ago, in the proceedings of the first IAU symposium on planetary nebulae, Lawrence H. Aller and Stanley J. Czyzak said that “the problem of determination of the chemical compositions of planetary and other gaseous nebulae constitutes one of the most exasperating problems in astrophysics”. Although the situation has greatly improved over the years, many important problems are still open and new questions have arrived to the field, which still is an active field of study. Here I will review some of the main aspects related to the determination of gaseous abundances in PNe and some relevant results derived in the last five years, since the last IAU symposium on PNe.

In Planetary Nebulae (PNe) and HII regions ionic abundances can be derived by using collisionally excited lines (CELs) or recombination lines (ORLs). Such abundances do not coincide for the same ion and usually abundances from ORLs are larger than those from CELs by factors of 2 or larger. The origin of the discrepancy, known as the Abundance Discrepancy Factor is an open problem in astrophysics of gaseous nebulae. It has been attributed to temperature fluctuations in the plasma, tiny metal-rich inclusions embedded in the H-rich plasma, gas inhomogeneities or other processes. In this work we analyze the kinematical behavior of CELs and ORLs in two PNe ionized by [WC] stars, finding that kinematics of ORLS is incompatible with the kinematics of CELs. In particular the expansion velocities from CELs and ORLs for the same ion are different, indicating that ORLs seem to be produced in zones nearer the central star than CELs. This is in agreement with results found by other authors for individual PNe.

The discrepancy between abundances computed using optical recombination lines (ORLs) and collisionally excited lines (CELs) is a major, unresolved problem with significant implications for the determination of chemical abundances throughout the Universe. In planetary nebulae (PNe), the most common explanation for the discrepancy is that two different gas phases coexist: a hot component with standard metallicity, and a much colder plasma enhanced in heavy elements. This dual nature is not predicted by mass loss theories, and direct observational support for it is still weak. In this work, we present our recent findings that demonstrate that the largest abundance discrepancies are associated with close binary central stars. OSIRIS-GTC tunable filter imaging of the faint O ii ORLs and MUSE-VLT deep 2D spectrophotometry confirm that O ii ORL emission is more centrally concentrated than that of [Oiii] CELs and, therefore, that the abundance discrepancy may be closely linked to binary evolution.

Recent work (Corradi et al. 2015; Jones et al. 2016) has shown that the phenomenon of extreme abundance discrepancies, where recombination line abundances exceed collisionally excited line abundances by factors of 10 or more, seem to be strongly associated with planetary nebulae with close binary central stars. To further investigate, we have obtained spectra of a sample of nebulae with known close binary central stars, using FORS2 on the VLT, and we have discovered several new extreme abundance discrepancy objects. We did not find any non-extreme discrepancies, suggesting that a very high fraction of nebulae with close binary central stars also have an extreme abundance discrepancy.

Neutron(n)-capture elements are produced by s-process nucleosynthesis in low- and intermediate-mass AGB stars, and therefore can be enriched in planetary nebulae (PNe). In the last ten years, n-capture elements have been detected in more than 100 PNe in the Milky Way and nearby galaxies. In some objects, several different n-capture elements have been detected, providing valuable constraints to models of AGB nucleosynthesis and evolution. These detections have motivated theoretical and experimental investigations of the atomic data needed to derive accurate n-capture element abundances. In this review, I discuss the methods and results of these atomic data studies, and their application to abundance determinations in PNe.

We have developed a method for determining elemental Fe-group abundances in planetary nebulae using an infrared emission line of Zn, the least refractory Fe-group species. Many planetary nebulae, particularly those of the Milky Way’s thick disk and bulge, display subsolar [Fe/H] (as inferred from Zn) although their abundances of α elements such as O, S, and Ar are nearly solar. We discuss the implications for determining enhancements of species synthesized by the progenitor star during the AGB (e.g., s-process products), and for galactic chemical evolution in view of the metallicity dependence of AGB nucleosynthetic yields.

Planetary nebulae retain the signature of the nucleosynthesis and mixing events that occurred during the previous AGB phase. Observational signatures complement observations of AGB and post-AGB stars and their binary companions. The abundances of the elements heavier than iron such as Kr and Xe in planetary nebulae can be used to complement abundances of Sr/Y/Zr and Ba/La/Ce in AGB stars, respectively, to determine the operation of the slow neutron-capture process (the s process) in AGB stars. Additionally, observations of the Rb abundance in Type I planetary nebulae may allow us to infer the initial mass of the central star. Several noble gas components present in meteoritic stardust silicon carbide (SiC) grains are associated with implantation into the dust grains in the high-energy environment connected to the fast winds from the central stars during the planetary nebulae phase.

We compare recent precise/reliable nebular abundances - as derived from high-quality optical spectra and the most recent ICFs - in a sample of Galactic planetary nebulae (PNe) with nucleosynthesis predictions (HeCNOCl) from asymptotic giant branch (AGB) ATON models in the metallicity range Z⊙/4 < Z < 2Z⊙. According to the infrared dust features, the sample is divided among carbon-, oxygen-, and double-dust chemistry (CC, OC, and DC, respectively), providing an independent proxy for the nature of the PNe progenitors. Our AGB models, with diffusive overshooting from all the convective borders, nicely reproduce the O overabundances observed in CC PNe, indicating that they evolve from low-Z low-mass (∼1 −3 M⊙) AGB stars. This indicates that O is not always a good indicator of the original ISM metallicity and that the O production by low-mass stars should be considered in galactic-evolution models. The lowest metallicity OC PNe evolve from low-mass (∼1 M⊙) O-rich AGBs, while the higher metallicity ones (all with uncertain dust classifications) display a chemical pattern similar to the DC PNe. In agreement with the recent literature, the DC PNe mostly descend from high-mass (M > 3.5 M⊙) solar/supersolar metallicity AGBs that experience hot bottom burning (HBB), but other formation channels in low-mass AGBs like extra mixing, stellar rotation, binary interaction, or He pre-enrichment cannot be disregarded until more accurate C/O ratios can be obtained. Two DC PNe show the imprint of advanced CNO processing and deep second dredge-up, suggesting progenitors masses close to the limit to evolve as core collapse supernovae (above 6 M⊙). Their actual C/O ratios, if confirmed, indicate contamination from the third dredge-up, rejecting the hypothesis that the chemical composition of such high-metallicity massive AGBs is modified exclusively by HBB.

Determination of the 3He isotope is important to many fields of astrophysics, including stellar evolution, chemical evolution, and cosmology. The isotope is produced in stars which evolve through the planetary nebula phase. Planetary nebulae are the final evolutionary phase of low- and intermediate-mass stars, where the extensive mass lost by the star on the asymptotic giant branch is ionised by the emerging white dwarf. This ejecta quickly disperses and merges with the surrounding ISM. 3He abundances in planetary nebulae have been derived from the hyperfine transition of the ionised 3He, 3He+, at the radio rest frequency 8.665 GHz. 3He abundances in PNe can help test models of the chemical evolution of the Galaxy. Many hours have been put into trying to detect this line, using telescopes like the Effelsberg 100m dish of the Max Planck Institute for Radio Astronomy, the National Radio Astronomy Observatory (NRAO) 140-foot telescope, the NRAO Very Large Array, the Arecibo antenna, the Green Bank Telescope, and only just recently, the Deep Space Station 63 antenna from the Madrid Deep Space Communications Complex.

X-ray emission from planetary nebulae (PNe) provides unique insight on the formation and evolution of PNe. Past observations and the ongoing Chandra Planetary Nebulae Survey (ChanPlaNS) provide a consensus on the two types of X-ray emission detected from PNe: extended and compact point-like sources. Extended X-ray emission arises from a shocked “hot bubble” plasma that resides within the nebular shell. Cooler than expected hot bubble plasma temperatures spurred a number of potential solutions with one emerging as the likely dominate process. The origin of X-ray emission from compact sources at the location of the central star is less clear. These sources might arise from one or combinations of the following processes: self-shocking stellar winds, spun-up binary companions, and/or accretion, perhaps from mass transfer, PN fallback, or debris disks. In the discovery phase, X-ray studies of PNe have mainly focused on the origin of the various emission processes. New directions incorporate multi-wavelength observations to study the influence of X-ray emission on the rest of the electromagnetic spectrum.

We developed a model for wind-blown bubbles with temperature and density profiles based on self-similar solutions including thermal conduction. We constructed also heat-conduction bubbles with chemical discontinuities. The X-ray emission is computed using the well-documented CHIANTI code (v6.0.1). These bubble models are used to (re)analyse the high-resolution X-ray spectrum of the hot bubble of BD+30°3639, and they appeared to be much superior to constant temperature approaches.

We found for the X-ray emission of BD+30°3639 that temperature-sensitive and abundance-sensitive line ratios computed on the basis of heat-conducting wind-blown bubbles and with abundances as found in the stellar photosphere/wind can only be reconciled with the observations if the hot bubble of BD+30°3639 is chemically stratified, i.e. if it contains also a small mass fraction (≃ 3 %) of hydrogen-rich matter immediately behind the conduction front. Neon appears to be strongly enriched, with a mass fraction of at least about 0.06.

The fast stellar winds can blow bubbles in the circumstellar material ejected from previous phases of stellar evolution. These are found at different scales, from planetary nebulae (PNe) around stars evolving to the white dwarf stage, to Wolf-Rayet (WR) bubbles and up to large-scale bubbles around massive star clusters. In all cases, the fast stellar wind is shock-heated and a hot bubble is produced. Processes of mass evaporation and mixing of nebular material and heat conduction occurring at the mixing layer between the hot bubble and the optical nebula are key to determine the thermal structure of these bubbles and their evolution. In this contribution we review our current understanding of the X-ray observations of hot bubbles in PNe and present the first spatially-resolved study of a mixing layer in a PN.

Infrared spectra from the Spitzer Space Telescope trace the evolution of carbon-rich dust from the asymptotic giant branch (AGB) to young planetary nebulae (PNe). On the AGB, amorphous carbon dominates the dust, but SiC and MgS also appear. In more evolved systems with warmer central stars, the spectra reveal the unidentified 21 μm feature, features from aliphatic hydrocarbons, and spectra from polycyclic aromatic hydrocarbons (PAHs), often with shifted feature positions indicative of the presence of aliphatics. More evolved systems with hot central stars show more typical PAH spectra, along with fullerenes and/or an emission feature known as the big-11 feature at ~11 μm. This features arises from a combination of SiC and PAHs, and it is usually accompanied by a shoulder at 18 μm, which while unidentified might be from cool silicate grains. The strong emission from MgS and SiC in young PNe probably arises from coatings on carbonaceous grains.