Post-asymptotic giant branch (post-AGB) binary stars are evolved systems that host circumbinary discs formed through mass loss during late stage binary interactions. Their structural, morphological, kinematic, and chemical similarities to planet-forming discs suggest that these systems may act as sites of “second generation” planet formation. In this study, we assess whether the disc instability mechanism – a proposed pathway for rapid, giant planet formation in some protoplanetary discs - can operate in post-AGB discs; motivated by their short lifetimes (104−5years). Using the Toomre criterion under well motivated assumptions for disc structure and size, mass, and thermal properties, we assess the conditions for gravitational instability. We first benchmark our analytical framework using well studied protoplanetary disc systems (including HL Tauri, Elias 2-27, GQ Lupi) before applying the same analysis to observed post-AGB discs. We find that post-AGB discs are generally gravitationally stable at present, due primarily to their low masses. Using viscous disk theory, we find that the discs were stable against collapse even in the past, when their masses were potentially higher. In contrast, several protoplanetary discs analysed in the same way show that they likely experienced gravitationally unstable phases early on. We also find that higher viscosity parameters (α ∼ 10−2) are better aligned with expected post-AGB disc lifetimes. Finally, we revisit the planet formation scenario proposed for the post-common envelope system NN Ser, first carried out by Schleicher and Dreizler and we show that gravitational instability could be feasible under specific, high disc mass assumptions. Overall, our results provide the first systematic theoretical assessment of gravitational instability in post-AGB discs, demonstrating that this mechanism is unlikely to dominate second generation planet formation in these systems, and underscoring the need to explore alternative pathways - such as core accretion - in future studies.

Post-asymptotic giant branch (Post-AGB) and post-red giant branch (post-RGB) binaries with stable circumbinary discs provide key insights into late stellar and disc evolution, revealing how binary interactions shape disc structure and stellar surface composition. A defining trait of such systems is the observed underabundance of refractory elements in the stellar photosphere relative to volatile elements – photospheric chemical depletion – resulting from the star accreting volatile-rich circumstellar gas. In this study, we investigated the link between photospheric depletion and disc evolution by focusing on post-AGB/post-RGB binaries with low infrared excess (hereafter ‘faint disc’ targets). We analysed high-resolution optical spectra from HERMES/Mercator and UVES/VLT for 6 Galactic and 2 LMC targets. Using E-iSpec, we homogeneously derived atmospheric parameters and chemical abundances of 29 elements from carbon to europium, and included NLTE corrections for 15 elements from carbon to barium that we calculated using pySME and pre-computed grids of departure coefficients. All targets exhibit ‘saturated’ depletion patterns, which we characterised using two-piece linear fits defined by three parameters: initial metallicity ([M/H]$_0$), turn-off temperature ($T_\textrm{turn-off}$), and depletion scale ($\nabla_\textrm{ 100 K}$). Among several findings, we highlight the bimodal distribution of $T_\textrm{turn-off}$ in faint disc targets, which allows classification into two subgroups analogous to full discs with continuous, optically thick dust ($T_\textrm{turn-off} \gt 1 100$ K), and transition discs with inner clearing ($T_\textrm{turn-off} \lt 1 100$ K). Our results imply that faint disc targets likely represent the final stages of disc dissipation, highlighting the diversity of depletion profiles, the complexity of disc-binary interactions, and the need to understand the rarity and evolution of faint disc systems.

We propose that certain white dwarf (WD) planets, such as WD 1856+534 b, may form out of material from a stellar companion that tidally disrupts from common envelope evolution with the WD progenitor star. The disrupted companion shreds into an accretion disc, out of which a gas giant protoplanet forms due to gravitational instability. To explore this scenario, we make use of detailed stellar evolution models consistent with WD 1856+534. The minimum mass companion that produces a gravitationally unstable disc after tidal disruption is $\sim$$0.15\,\mathrm{M_\odot}$. In this scenario, WD 1856+534 b might have formed at or close to its present separation, in contrast to other proposed scenarios where it would have migrated in from a much larger separation. Planet formation from tidal disruption is a new channel for producing second-generation planets around WDs.

Asymptotic giant branch (AGB) stars play a significant role in our understanding of the origin of the elements. They contribute to the abundances of C, N, and approximately 50% of the abundances of the elements heavier than iron. An aspect often neglected in studies of AGB stars is the impact of a stellar companion on AGB stellar evolution and nucleosynthesis. In this study, we update the stellar abundances of AGB stars in the binary population synthesis code binary_c and calibrate our treatment of the third dredge-up using observations of Galactic carbon stars. We model stellar populations of low- to intermediate-mass stars at solar-metallicity and examine the stellar wind contributions to C, N, O, Sr, Ba, and Pb yields at binary fractions between 0 and 1. For a stellar population with a binary fraction of 0.7, we find $\sim$20–25% less C and s-process elements ejected than from a population composed of only single stars, and we find little change in the N and O yields. We also compare our models with observed abundances from Ba stars and find our models can reproduce most Ba star abundances, but our population estimates a higher frequency of Ba stars with a surface [Ce/Y] > $+0.2\,$dex. Our models also predict the rare existence of Ba stars with masses $ \gt 10\,\textrm{M}_{\odot}$.

We use oxygen and argon abundances for planetary nebulae (PNe) with low internal extinction (progenitor ages of > 4.5 Gyr) and high extinction (progenitor ages < 2.5 Gyr), as well as those of the HII regions, to constrain the chemical enrichment and star formation efficiency in the thin and thicker discs of M31. As the argon element is produced in larger fractions by Type Ia supernovae compared to oxygen, we found that the mean log(O/Ar) values of PNe as a function of their argon abundances [12 + log(Ar/H)] trace the interstellar medium (ISM) conditions at the time of birth of the M31 disc PN progenitors. Thus, the chemical enrichment and star formation effciency information encoded in the [α/Fe] versus [Fe/H] distribution of stars is also imprinted in the oxygen-to-argon abundance ratio log(O/Ar) versus argon abundance for the nebular emissions of the different stellar evolution phases. The chemical evolution model that reproduces the mean log(O/Ar) values as a function of argon abundance for the high- and low-extinction PNe requires a second infall of metal-poorer gas during a gas-rich (wet) satellite merger for the M31 disc region within 14 kpc. A strong starburst is ongoing in the intermediate radial range (14 < RGC< 18 kpc). In the outer M31 disc (RGC > 18 kpc), the log(O/Ar) versus argon abundance distribution of the younger high-extinction PNe indicates that they too were formed in a burst, though mostly from the metal-poorer gas. In M31, the thin disc is younger and less radially extended, formed stars at a higher star formation effciency, and had a faster chemical enrichment timescale than the more extended thicker disc. Both the thin and thicker discs in M31 reach similar high argon abundances (12 + log(Ar/H))∼ 6.7. The chemical and structural properties of the thin and thicker discs in M31 are thus remarkably different from those determined for the Milky Way thin and thick discs.

Bipolar or multipolar lobes in pre-planetary nebulae (pPNe) often exhibit intertwined outer whorled patterns, resulting from stellar wind matter accumulation during the asymptotic giant branch (AGB) phase. These structures are likely triggered by stellar or substellar companions. We regard that CW Leonis currently stands at a critical transition moment, providing a vivid illustration of the progression from an AGB star in a binary system to a pPN. We have found that CW Leonis has shown significant enhancements in its optical and near-infrared light curves over the past two decades, with the recent Hubble Space Telescope image finally revealing the long-awaited central star. Utilizing an eccentric-orbit binary model, we can reproduce the position-angle dependence of the expansion velocity in the whorled pattern around CW Leonis, suggesting a nearly face-on orbital inclination. Its contradiction to the features in the innermost circumstellar envelope, corresponding to a nearly edge-on inclination, may imply the presence of an additional companion. Our updated theoretical framework explores the complexity of the whorled pattern. Further identifying and monitoring phase-transition candidates at the tip of the AGB will provide valuable insights into the AGB-pPN transition and the role of companions in shaping the morphological evolution of these stellar objects.

As AGB stars evolve to planetary nebulae (PNe), the geometry of the ejected mass transforms from nearly spherical to extremely aspherical. The mechanisms governing this transformation are plausibly linked to binarity and the associated production of disks and jets during the transitional (post-AGB) evolutionary stage. We are carrying out an ALMA survey of a representative sample of bipolar/multipolar post-AGB objects to obtain high angular-resolution (0″.11) observations of the CO(3–2) and 6–5 emission and study the collimated outflows and central disks. We present highlights from our ongoing survey – e.g., the presence of bipolar or multipolar high-velocity outflows, dense toroidal waists, and in one case, a circular ring around the central bipolar nebula. We will use radiative transfer modeling to derive accurate outflow momenta, masses, and mass-loss rates for our sample, thereby constraining different classes of binary PN-shaping models.

Dynamics in the atmospheres of central stars of cool protoplanetary nebulae can contribute to shaping of the subsequent planetary nebulae by giving rise to stellar wind. We have investigated variation in high-resolution spectra of three relatively cool post-asymptotic giant branch stars that are surrounded by non-spherical nebulae. Spectra of HD 161796 show variability in blue wings of weak and medium strength absorption lines while red wings remain unchanged. This could be caused by warm variable outflow. An episode of infall of matter is revealed in spectra of IRAS 22272+5435. Also, episodes of molecular lines in emission are detected for this star. The emissions appear to be related to the star’s brightness changes. Qualitatively the same molecular line variability is observed in the spectrum of IRAS Z02229+6208. Possibly, the observed spectral variations are a consequence of similar dynamics as in the atmospheres of asymptotic giant branch stars.

A family of unidentified infrared emission (UIE) bands is widely observed in planetary nebulae. We suggest that the carriers of the UIE bands are mixed aromatic/aliphatic organic nanoparticles (MAONs) synthesized over thousand-year time scales in the nebulae. The possible chemical pathways of synthesis is discussed. These organics are ejected into the interstellar medium and could have enriched the primordial Solar System, leading to the reservoir of complex organics in comets, asteroids, planetary satellites, and interplanetary dust particles.

We present a radial velocity and model atmosphere analysis of both components of the spectroscopic binary central star of NGC 1514, based on high-resolution, high-signal-to-noise-ratio spectrograms taken with the CFH and Subaru telescopes at Maunakea, Hawaii. Together with the Gaia parallax and other data from the literature, all this information permits to determine the basic stellar parameters (Teff, L, log g, masses) of both binary components. This allows us to empirically test the theoretical post-AGB mass-luminosity relation.

The Red Rectangle is a nebula surrounding the post-AGB star HD 44179. It is the prototype of a particular class of nebulae associated with post-AGB binaries characterised by the presence of stable circumbinary disks in (quasi-)Keplerian rotation. Here we present the results of new high-resolution (0.″ 02″05) ALMA observations of continuum and line emissions at 0.9 mm. The continuum maps are analysed through a simple model of dust emission, which can reproduce the observational data. We find that most dust emission in the Red Rectangle is concentrated in the central regions of the rotating disk and that the settlement of dust grains onto the equatorial plane is very significant, particularly in comparison with the much larger scale height displayed by the gas distribution. The diameter of the dust-emitting region is about 250 au, with a total width of about 50 au. This region coincides with the warm PDR where certain molecules (like HCN), CI, and CII are presumably formed, as well as probably PAHs. From the spectral index, we confirm the presence in the disk of large grains, with a typical radius of about 150 μm, which supports the long-lived hypothesis for this structure. We also confirm the existence of a compact ionised wind at the centre of the nebula, probably emerging from the accretion disk around the companion, for which we derive an extent of about 10 au and a total flux of 8 mJy. We also briefly present the results on molecular lines of 12CO, 13CO, and other less abundant species.

Even after decades of usage as an extragalactic standard candle, the universal bright end of the planetary nebula luminosity function (PNLF) still lacks a solid theoretical explanation. Until now, models have modeled planetary nebulae (PNe) from artificial stellar populations, without an underlying cosmological star formation history. We present PICS (PNe In Cosmological Simulations), a novel method of modeling PNe in cosmological simulations, through which PN populations for the first time naturally occur within galaxies of diverse evolutionary pathways. We find that only by using realistic stellar populations and their metallicities is it possible to reproduce the bright end of the PNLF for all galaxy types. In particular, the dependence of stellar lifetimes on metallicity has to be accounted for to produce bright PNe in metal-rich populations. Finally, PICS reproduces the statistically complete part of the PNLF observed around the Sun, down to six orders of magnitude below the bright end.

The modelling of the evolutionary phases beyond the asymptotic giant branch attracts the interest of the astrophysical community because it allows the determination of the properties of progenitor stars and to deduce the efficiency of the mechanisms able to alter the surface chemistry of the stars evolving through the asymptotic giant branch. This has been possible since improvements in the modelling of these phases, which allow a reliable determination of the luminosity with which stars evolve after the termination the asymptotic giant branch evolution.

The surface chemistry of post-asymptotic giant branch stars and planetary nebulae is shown to be tightly correlated to the various processes taking place during the asymptotic giant branch evolution. The possibility of using the observed infrared excess of these evolved stars to derive information on the dust formation process during the previous evolutionary phases is also discussed.

Ongoing improvements of sub-mm- and mm-range interferometers and single-dish radiotelescopes are progressively allowing the detailed study of planetary nebulae (PNe) in molecular species other than 12CO and 13CO. We are implementing a new set of tables for extending the capabilities of the morpho-kinematical modelling tool SHAPE+shapemol, so radiative transfer in molecular species beyond 12CO and 13CO, namely C17O, C18O, HCN, HNC, CS, SiO, HCO+, and N2H+, are enabled under the Large Velocity Gradient approximation with the ease of use of SHAPE. We present preliminary results on the simultaneous analysis of a plethora of IRAM-30m and HERSCHEL/HIFI spectra, and NOEMA maps of different species in the pre-PN nebula M 1-92, which show interesting features such as a previously undetected pair of polar, turbulent, high-temperature blobs, or a 17O/18O isotopic ratio of 1.7, which indicates the AGB should have turned C-rich, as opposed to the apparent nature of its O-rich nebula.

We present an analysis of high resolution spectra in the optical and near infrared wavelength region for cool protoplanetary nebulae with the goal to identify lines of carbon bearing molecules and some less studied neutron capture elements. The site of formation of CN and C2 lines in the spectra of IRAS 22272+5435 and IRAS Z02229+6208 appears to be an extended stellar atmosphere and inner/outer CSE. The abundance of thulium was estimated for the first time in the photosphere of IRAS22272+5435, log ε(Tm) ≃ 1.5.

A high sensitivity and high-angular resolutions infrared space telescope, the James Webb Space Telescope (JWST), allowed us to study dust and molecules in unprecedented details. This contribution highlights the first year of JWST’s scientific operation, and reports prospects of dust and molecular studies in the coming future.

We report on recent results from our successful and pioneering observational program with ALMA to study emerging ultracompact HII regions of pre-planetary nebulae (pPNe) using mm-wavelength recombination lines (mRRLs) as new optimal tracers. We focus on our study of two poster-child pPNe, namely, M 2-9 and CRL 618. We reveal the structure and kinematics of the enigmatic inner nebular regions of these objects with an unprecedented angular resolution down to 20-30 mas (∼15-30 AU linear scales). For both targets, the ionized central regions are elongated along the main symmetry axis of the large-scale nebulae, consistent with bipolar winds, and show notable axial velocity gradients with expansion velocities of up to ∼100 kms. The intensity and width of the H30α profiles are found to be time variable, denoting changes on scales of a few years of the physical properties and kinematics of the present-day post-AGB ejections. Our ongoing analysis involves 3D, non-LTE radiative transfer modeling of the mRRLs and free-free continuum emission. This approach allows us to provide an exceptionally detailed description of the physical conditions in the innermost layers of these well known pPNe.

The carbon-rich star V Hydrae is known to be in its transition from the asymptotic giant branch phase to a bipolar planetary nebula. However, the origin of its bipolar shape, as well as its long-term periodic obscuration, are not well understood. Using high-resolution spectra from the HERMES/Mercator spectrograph spanning over 10-yr monitoring, we disentangled the orbital signal from the intrinsic Mira-like variations of V Hya. The orbital solution obtained is compatible with the light-curve modulation, which we interpret as an obscuration caused by an extended dusty environment surrounding the unseen companion. Additionally, we report high-velocity phase-dependent absorption features in the sodium doublet and explain them by the presence of a conical jet ejected from the (accretion disc of the) companion. We confront our binary scenario with previous radio observations and discuss how this jet could shape the large-scale bipolar outflow of the system.

We performed numerical simulations of the common envelope (CE) interaction between two intermediate-mass asymptotic giant branch (AGB) stars and their low-mass companions. For the first time, formation and growth of dust in the envelope is calculated explicitly. We find that the first dust grains appear as early as ∼1–3 yrs after the onset of the CE, and are smaller than grains formed later. As the simulations progress, a high-opacity dusty shell forms, resulting in the CE photosphere being up to an order of magnitude larger than it would be without the inclusion of dust. At the end of the simulations, the total dust yield is ∼8.2×10−3 M⊙ (∼2.2×10−2 M⊙) for a CE with a 1.7 M⊙ (3.7 M⊙) AGB star. Dust formation does not substantially lead to more mass unbinding or substantially alter the orbital evolution.

We proposed a machine learning approach to identify and distinguish dusty stellar sources employing supervised methods and categorizing point sources, mainly evolved stars, using photometric and spectroscopic data collected over the IR sky. Spectroscopic data is typically used to identify specific infrared sources. However, our goal is to determine how well these sources can be identified using multiwavelength data. Consequently, we developed a robust training set of spectra of confirmed sources from the Large and Small Magellanic Clouds derived from SAGE-Spec Spitzer Legacy and SMC-Spec Spitzer Infrared Spectrograph (IRS) spectral catalogs. Subsequently, we applied various learning classifiers to distinguish stellar subcategories comprising young stellar objects (YSOs), C-rich asymptotic giant branch (CAGB), O-rich AGB stars (OAGB), Red supergiant (RSG), and post-AGB stars. We have classified around 700 counts of these sources. It should be highlighted that despite utilizing the limited spectroscopic data we trained, the accuracy and models’ learning curve provided outstanding results for some of the models. Therefore, the Support Vector Classifier (SVC) is the most accurate classifier for this limited dataset.