In this proceeding, we present a short review of the fascinating nebulosities of symbiotic binary R Aquarii. The R Aquarii system, comprising the central binary and surrounding nebular material, has been the subject of near-continuous study since its discovery, with a few hundred papers listed in ADS. As such, it is impossible to provide here the comprehensive review that R Aquarii deserves, instead we chose to focus on the nebulosities -- covering both our own research and other relevant results from the literature.

We are studying a class of binary post-AGB stars that seem to be systematically surrounded by equatorial disks and slow outflows. Although the rotating dynamics had only been well identified in three cases, the study of such structures is thought to be fundamental to the understanding of the formation of nebulae around evolved stars. We present ALMA maps of 12CO and 13CO J=3-2 lines in one of these sources, IRAS08544-4431. We analyzed the data by means of nebula models, which account for the expectedly composite source and can reproduce the data. From our modeling, we estimated the main nebula parameters, including the structure and dynamics and the density and temperature distributions. We discuss the uncertainties of the derived values and, in particular, their dependence on the distance. Our observations reveal the presence of an equatorial disk in rotation; a low-velocity outflow is also found, probably formed of gas expelled from the disk. The main characteristics of our observations and modeling of IRAS08544-4431 are similar to those of better studied objects, confirming our interpretation. The disk rotation indicates a total central mass of about 1.8 Mo, for a distance of 1100 pc. The disk is found to be relatively extended and has a typical diameter of ~ 4 10^16 cm. The total nebular mass is ~ 2 10^-2 Mo, of which ~ 90% corresponds to the disk. Assuming that the outflow is due to mass loss from the disk, we derive a disk lifetime of ~ 10000 yr. The disk angular momentum is found to be comparable to that of the binary system at present. Assuming that the disk angular momentum was transferred from the binary system, as expected, the high values of the disk angular momentum in this and other similar disks suggest that the size of the stellar orbits has significantly decreased as a consequence of disk formation.

R Aquarii is a symbiotic binary surrounded by a large and complex nebula with a prominent curved jet. It is one of the closest known symbiotic systems, and therefore offers a unique opportunity to study the central regions of these systems and the formation and evolution of astrophysical jets. We studied the evolution of the central jet and outer nebula of R Aqr taking advantage of a long term monitoring campaign of optical imaging, as well as of high-resolution integral field spectroscopy. Narrow-band images acquired over a period of more than 21 years are compared in order to study the expansion and evolution of all components of the R Aqr nebula. The magnification method is used to derive the kinematic ages of the features that appear to expand radially. Integral field spectroscopy of the OIII 5007A emission is used to study the velocity structure of the central regions of the jet. New extended features, further out than the previously known hourglass nebula, are detected. The kinematic distance to R Aqr is calculated to be 178 pc using the expansion of the large hourglass nebula. This nebula of R Aqr is found to be roughly 650 years old, while the inner regions have ages ranging from 125 to 290 years. The outer nebula is found to be well described by a ballistic expansion, while for most components of the jet strong deviations from such behaviour are found. We find that the Northern jet is mostly red-shifted while its Southern part is blue-shifted, apparently at odds with findings from previous studies but almost certainly a consequence of the complex nature of the jet and variations in ionisation and illumination between observations.

The Galactic object MWC 137 was suggested to belong to the group of B[e] supergiants. However, with its large-scale optical bipolar ring nebula and the high velocity jet and knots, it is a rather atypical representative of this class. We performed multi-wavelength observations spreading from the optical to the radio regime. Based on optical imaging and long-slit spectroscopic data we found that the northern parts of the large-scale nebula are predominantly blue-shifted, while the southern regions appear mostly red-shifted. We developed a geometrical model consisting of two double-cones. While various observational features can be approximated with such a scenario, the observed velocity pattern is more complex. Using near-infrared integral-field unit spectroscopy we studied the hot molecular gas in the close vicinity of the star. The emission from the hot CO gas arises in a small-scale disk revolving around the star on Keplerian orbits. While the disk itself cannot be spatially resolved, its emission is reflected by dust arranged in arc-like structures and clumps surrounding MWC 137 on small scales. In the radio regime we mapped the cold molecular gas in the outskirts of the optical nebula. We found that large amounts of cool molecular gas and warm dust embrace the optical nebula in the east, south and west. No cold gas or dust were detected in the north and north-western regions. Despite the new insights on the nebula kinematics gained from our studies, the real formation scenario of the large-scale nebula remains an open issue.

We aim to study the presence of both rotation and expansion in post-AGB nebulae, in particular around IW Car, a binary post-AGB star that was suspected to be surrounded by a Keplerian disk. We obtained high-quality ALMA observations of 12CO and 13CO J=3-2 lines in IW Car. The maps were analyzed by means of a simplified model of CO emission, based on those used for similar objects. Our observations clearly show the presence of gas components in rotation, in an equatorial disk, and expansion, which shows an hourglass-like structure with a symmetry axis perpendicular to the rotation plane and is probably formed of material extracted from the disk. Our modeling can reproduce the observations and shows moderate uncertainties. The rotation velocity corresponds to a central stellar mass of approximately 1 Mo. We also derive the total mass of the molecule-rich nebula, found to be of ~ 4 10e-3 Mo; the outflow is approximately eight times less massive than the disk. From the kinematical age of the outflow and the mass values derived for both components, we infer a (future) lifetime of the disk of approximately 5000--10000 yr.

The peculiar emission-line star MWC 137 with its extended optical nebula was recently classified as B[e] supergiant. To study the spatial distribution of its circumstellar molecular gas on small and large scales, we obtained near-infrared and radio observations using SINFONI and APEX, respectively. We find that the hot CO gas is arranged in moving clumpy ring and shell structures close to the star, while a cold CO envelope is encircling the borders of the optical nebula from the south to the west.

We present new high-quality ALMA observations of the Red Rectangle (a well known post-AGB object) in C17O J=6-5 and H13CN J=4-3 line emission and results from a new reduction of already published 13CO J=3-2 data. A detailed model fitting of all the molecular line data, including previous maps and single-dish spectra, was performed using a sophisticated code. These observations and the corresponding modeling allowed us to deepen the analysis of the nebular properties. We also stress the uncertainties in the model fitting. We confirm the presence of a rotating equatorial disk and an outflow, which is mainly formed of gas leaving the disk. The mass of the disk is ~ 0.01 Mo, and that of the CO-rich outflow is ~ 10 times smaller. High temperatures of ~ 100 K are derived for most components. From comparison of the mass values, we roughly estimate the lifetime of the rotating disk, which is found to be of about 10000 yr. Taking data of a few other post-AGB composite nebulae into account, we find that the lifetimes of disks around post-AGB stars typically range between 5000 and more than 20000 yr. The angular momentum of the disk is found to be high, ~ 9 Mo AU km/s, which is comparable to that of the stellar system at present. Our observations of H13CN show a particularly wide velocity dispersion and indicate that this molecule is only abundant in the inner Keplerian disk, at ~ 60 AU from the stellar system. We suggest that HCN is formed in a dense photodissociation region (PDR) due to the UV excess known to be produced by the stellar system, following chemical mechanisms that are well established for interstellar medium PDRs and disks orbiting young stars. We further suggest that this UV excess could lead to the efficient formation and excitation of PAHs and other C-bearing macromolecules, whose emission is very intense in the optical counterpart.

Modern instrumentation in radioastronomy constitutes a valuable tool for studying the Universe: ALMA has reached unprecedented sensitivities and spatial resolution, while Herschel/HIFI has opened a new window for probing molecular warm gas (~50-1000 K). On the other hand, the software SHAPE has emerged in the past few years as a standard tool for determining the morphology and velocity field of different kinds of gaseous emission nebulae via spatio-kinematical modelling. SHAPE implements radiative transfer solving, but it is only available for atomic species and not for molecules. Being aware of the growing importance of the development of tools for simplifying the analyses of molecular data, we introduce shapemol, a complement to SHAPE, with which we intend to fill the so-far under-developed molecular niche. shapemol enables user-friendly, spatio-kinematic modelling with accurate non-LTE calculations of excitation and radiative transfer in CO lines. It allows radiative transfer solving in the 12CO and 13CO J=1-0 to J=17-16 lines, but its implementation permits easily extending the code to different molecular species. shapemol allows easily generating synthetic maps and line profiles to match against interferometric or single-dish observations. We fully describe shapemol and discuss its limitations and the sources of uncertainty to be expected in the final synthetic profiles or maps. As an example of the power and versatility of shapemol, we build a model of the molecular envelope of the planetary nebula NGC 6302 and compare it with 12CO and 13CO J=2-1 interferometric maps from SMA and high-J transitions from Herschel/HIFI. We find the molecular envelope to have a complex, broken ring-like structure with an inner, hotter region and several 'fingers' and high-velocity blobs, emerging outwards from the plane of the ring. We derive a mass of 0.11 Msun for the molecular envelope.

It is now clear that a binary formation pathway is responsible for a significant fraction of planetary nebulae, and this increased sample of known binaries means that we are now in a position to begin to constrain their influence on the formation and evolution of their host nebulae. Here, we will review current understanding of how binarity influences the resulting nebulae, based on observations and modelling of both the central binary systems and the planetary nebulae themselves. We will also highlight the most important test-cases which have proved the most interesting in studying the evolution of binaries into and through the planetary nebula phase.

We aim to study equatorial disks in rotation and axial outflows in post-AGB objects, as to disclose the formation and shaping mechanisms in planetary nebulae. So far, both disks and outflows had not been observed simultaneously. We have obtained high-quality ALMA observations of 12CO and 13CO J=3-2 and 12CO J=6-5 line emission in the Red Rectangle, the only post-AGB/protoplanetary object in which a disk in rotation has been mapped up to date. These observations provide an unprecedented description of the complex structure of this source. Together with an equatorial disk in rotation, we find a low-velocity outflow that occupies more or less the region placed between the disk and the optical X-shaped nebula. From our observations and preliminary modeling of the data, we confirm the previously known properties of the disk and obtain a first description of the structure, dynamics, and physical conditions of the outflow.

There is a group of binary post-AGB stars that show a conspicuous NIR excess, usually assumed to arise from hot dust in very compact possibly rotating disks. These stars are surrounded by significantly fainter nebulae than the "standard", well studied protoplanetary and planetary nebulae (PPNe, PNe). We present high-sensitivity mm-wave observations of CO lines in 24 objects of this type. CO emission is detected in most observed sources and the line profiles show that the emissions very probably come from disks in rotation. We derive typical values of the disk mass between 1e-3 and 1e-2 Mo, about two orders of magnitude smaller than the (total) masses of "standard" PPNe. The high-detection rate (upper limits being in fact not very significant) clearly confirm that the NIR excess of these stars arises from compact disks in rotation, very probably the inner parts of those found here. Low-velocity outflows are also found in about eight objects, with moderate expansion velocities of ~ 10 km/s, to be compared with the velocities of about 100 km/s often found in "standard" PPNe. Except for two sources with complex profiles, the outflowing gas in our objects represents a minor nebular component. Our simple estimates of the disk typical sizes yields values ~ 0.5 - 1 arcsec, i.e. between 5e15 and 3e16 cm. Estimates of the linear momenta carried by the outflows, which can only be performed in a few well studied objects, also yield moderate values, compared with the linear momenta that can be released by the stellar radiation pressure (contrary, again, to the case of the very massive and fast bipolar outflows in "standard" PPNe, that are strongly overluminous). The mass and dynamics of nebulae around various classes of post-AGB stars differ very significantly, and we can expect the formation of PNe with very different properties.

The role of central star binarity in the shaping of planetary nebulae (PNe) has been the subject of much debate, with single stars believed to be incapable of producing the most highly collimated morphologies. However, observational support for binary-induced shaping has been sadly lacking. Here, we highlight the results of a continuing programme to spatio-kinematically model the morphologies of all PNe known to contain a close binary central star. Spatio-kinematical modelling is imperative for these objects, as it circumvents the degeneracy between morphology and orientation which can adversely affect determinations of morphology based on imaging alone. Furthermore, spatio-kinematical modelling accurately determines the orientation of the nebular shell, allowing the theoretically predicted perpendicular alignment, between nebular symmetry axis and binary orbital plane, to be tested. To date, every PN subjected to this investigation has displayed the predicted alignment, indicating that binarity has played an important role in the formation and evolution of these nebulae. The further results from this programme will be key, not only in determining whether binary interaction is responsible for shaping the studied PNe, but also in assessing the importance of binarity in the formation and evolution of all PNe in general.

Close binarity can play a significant role in the shaping of planetary nebulae (PNe) as the system evolves through the common-envelope phase. We present the detection of two of the shortest orbital periods among PN binary central stars. These are Hen 2-428, a bipolar PN, and V458 Vul, a recent nova surrounded by a mildly bipolar planetary nebula. The properties of the central stars of these systems, of their nebulae and their possible fate are discussed.

The kinematics, shaping, density distribution, expansion distance, and ionized mass of the nebula Hen 2-104, and the nature of its symbiotic Mira are investigated. A combination of multi-epoch HST images and VLT integral field high-resolution spectroscopy is used to study the nebular dynamics both along the line of sight and in the plane of the sky. These observations allow us to construct a 3-D spatio-kinematical model of the nebula, which together with the measurement of its apparent expansion in the plane of the sky over a period of 4 years, provides the expansion parallax for the nebula. The integral field data featuring the [S{\sc ii}] $\lambda\lambda$671.7,673.1 emission line doublet provide us with a density map of the inner lobes of the nebula, which together with the distance estimation allow us to estimate its ionized mass. We find densities ranging from n$_e$=500 to 1000 cm$^{-3}$ in the inner lobes and from 300 to 500cm$^{-3}$ in the outer lobes. We determine an expansion parallax distance of 3.3$\pm$0.9 kpc to Hen 2-104, implying an unexpectedly large ionized mass for the nebula of the order of one tenth of a solar mass.

We investigate the dynamics of the nebula around the symbiotic star Hen 2-147, determine its expansion parallax, and compare it with the distance obtained via the Period-Luminosity relation for its Mira variable. The geometry of the nebula is found to be that of a knotty annulus of ionized gas inclined to the plane of sky and expanding with a velocity of ~90 km/s. A straightforward application of the expansion parallax method provides a distance of 1.5+-0.4 kpc, which is a factor of two lower than the distance of 3.0+-0.4 kpc obtained from the Period-Luminosity relationship for the Mira (which has a pulsation period of 373 days). The discrepancy is removed if, instead of expanding matter, we are observing the expansion of a shock front in the plane of the sky. This shock interpretation is further supported by the broadening of the nebular emission lines.

The structure and kinematics of the bipolar nebula Mz 3 have been investigated by means of HST, CTIO and ESO images and spectra. At least four distinct outflows have been identified which, from the inside to the outside, are the following: a pair of bright bipolar lobes, two opposite highly collimated column-shaped outflows, a conical system of radial structure, and a very dim, previously unnoticed, low-latitude and flattened (ring-like) radial outflow. A simple Hubble-law describes the velocity field of the ballisticaly expanding lobes, columns and rays, suggesting that their shaping has being done at very early stages of evolution, in a sort of eruptive events with increasing degree of collimation and expansion ages ranging from ~600 for the inner structures to ~1600 years (per kpc to the nebula) for the largest ones.