Rotating and Expanding Gas in Binary Post-AGB Stars

There is a class of binary post-AGB stars (binary system including a post-AGB star) that are surrounded by Keplerian disks and outflows resulting from gas escaping from the disk. To date, there are seven sources that have been studied in detail through interferometric millimeter-wave maps of CO lines (ALMA/NOEMA). For the cases of the Red Rectangle, IW Carinae, IRAS 08544-4431, and AC Herculis, it is found that around greater than 85% of the total nebular mass is located in the disk with Keplerian dynamics. The remainder of the nebular mass is located in an expanding component. This outflow is probably a disk wind consisting of material escaping from the rotating disk. These sources are the disk-dominated nebulae. On the contrary, our maps and modeling of 89 Herculis, IRAS 19125+0343, and R Scuti, which allowed us to study their morphology, kinematics, and mass distribution, suggest that, in these sources, the outflow clearly is the dominant component of the nebula (around 75% of the total nebular mass), resulting in a new subclass of nebulae around binary post-AGB stars: the outflow-dominated sources.Besides CO, the chemistry of this type of source has been practically unknown thus far. We also present a very deep single-dish radio molecular survey in the 1.3, 2, 3, 7, and 13 mm bands (around 600 h of telescope time). Our results and detections allow us to classify our sources as O- or C-rich. We also conclude that the calculated abundances of the detected molecular species other than CO are particularly low, compared with AGB stars. This fact is very significant in those sources where the rotating disk is the dominant component of the nebula.


Introduction
There is a type of post-AGB star characterized by their spectral energy distributions (SEDs), which shows a near-infrared (NIR) excess indicating the presence of hot dust close to the the stellar system [3,4]. Their IR spectra reveal the presence of highly processed dust grains, so the dust might be located in stable structures [5][6][7]. All the above suggests the presence of circumbinary disks. Their disk-like shape has been confirmed by interferometric IR data see, e.g., [8,9]). Their radial velocity curves reveal that the post-AGB stars are part of a binary system (see, e.g., [4]). The systematic detection of binary systems in these objects strongly suggests that the angular momentum of the disks comes from the stellar system.
Observations of 12 CO and 13 CO in the J = 2 − 1 and J = 1 − 0 lines (230.538 and 220.398 GHz, respectively) have been well analyzed in sources with such NIR excess [10].
There are two types of CO line profiles: (a) narrow CO line profiles characteristic of rotating disks and weak wings, which implies that most of the nebular mass is contained in the disk (with Keplerian dynamics), and (b) composite CO line profiles including a narrow component, which very probably represents emission from the rotating disk, and strong wings, which represents emission from the outflow, which could dominate the nebula [10]. These types of line profiles are also found in young stars surrounded by a rotating disk made of remnants of interstellar medium (ISM) and those expected from disk-emission modeling (see, e.g., [10,11]). These results indicate that the CO emission lines of our sources come from disks with Keplerian or quasi-Keplerian rotation.
The study of the chemistry of this class of binary post-AGB stars, together with the very detailed kinematic analysis of Keplerian disks and outflows around these sources, is based on published articles (see [1,2]). This paper is organized as follows. Technical information of our observations is given in Section 2. In Section 3, we present mm-wave interferometric maps and models of the most representative cases of a disk-dominated nebula (AC Herculis), an outflow-dominated nebula (R Scuti), and an intermediate case in between the disk-and the outflow-dominated nebula (89 Herculis). We present the first molecular survey in this kind of object in Section 4, together with discussions about molecular intensities and chemistry. Finally, we summarize our conclusions in Section 5.

Observations
We show interferometric maps of our sources using the NOEMA interferometer. Observations of the 12 CO J = 2 − 1 rotational transition were carried out towards AC Herculis, R Scuti, and 89 Herculis. Observations of the 13 CO J = 2 − 1 rotational transition were also obtained for 89 Herculis.
Our single-dish observations were performed using the 30 m IRAM telescope (Granada, Spain) and the 40 m Yebes telescope (Guadalajara, Spain). We observed at the 1.3, 2, 3, 7, and 13 mm bands. Our observations required a total telescope time of ∼600 h distributed over the two telescopes and for several projects to observe the nebulae around the next binary post-AGB stars: AC Her, the Red Rectangle, HD 52961, IRAS 19157−0257, IRAS 18123+0511, IRAS 19125+0343, AI CMi, IRAS 20056+1834, and R Sct.

NOEMA Observations
In this section, we present the results directly obtained from the observations for AC Her, R Sct, and 89 Her [1].
We show our NOEMA maps per velocity channel and position-velocity (PV) diagrams along the equatorial rotating disk and along the axis of the nebula (see Figures 1, 2, 3, 4, 5, and 6).

AC Herculis
The 12 CO J = 2 − 1 mm-wave interferometric maps are presented in Figure 1. We can see in the left panel of Figure 2 the PV diagram along the equatorial direction, and it very clearly shows the characteristic signature of rotation with Keplerian dynamics. On the contrary, the analysis of the PV diagram along the nebula axis direction helps us to detect the presence of an axially outflowing component (see Figure 2 right). A theoretical PV diagram along the axis direction in the presence of a disk with Keplerian dynamics might show emission with a form close to a rhombus, with equal or very similar emission in all four quadrants of the PV diagram. Nevertheless, we do not see equal emission in the four quadrants: we see slightly inclined emission at central velocities at around ±1 . This fact can be explained by the existence of an expanding component that surrounds the rotating disk.

R Scuti
We present combined NOEMA maps and 30 m maps of R Sct in 12 CO J = 2 − 1 emission in Figure 3 and PV diagrams in Figure 4. In both figures, we clearly see two components: an intense inner region and an extended component of around ∼40 surrounding the inner region. This extended and expanding component contains most of the total nebular mass (see Section 3.4.2). The PV diagram along the equatorial direction shows an intense central clump in the innermost region of the nebula that may represent the unresolved rotating disk (see Figure 4 left). The velocity dispersion from the inner (and unresolved) central condensation is similar to other post-AGB nebulae with disks, including a significant lack of blueshifted emission (see, e.g., [12]). The PV diagram along the nebula axis reveals the structure of the nebula (see Figure 4 right): we clearly see two large cavities at approximately ±10 . We see this type of structure in other pPNe, such as M 2−56 [13] or M 1−92 [14].

89 Her
We present NOEMA maps and PV diagrams of 89 Her of 13 CO J = 2 − 1 emission (and 12 CO J = 2 − 1; see [1] for further details) in Figures 5 and 6. We see an intense central clump and an extended hourglass-shaped structure surrounding this central clump. For a distance of 1 kpc, the size of the hourglass-like structure is, at least, 10,000 AU.

Models
Our models consist of a disk with present Keplerian dynamics and an extended and expanding component escaping from the rotating disk and surrounding it. The outflowing component can present different shapes, such as an hourglass, an ellipsoid, etc. We assume LTE populations, which is a reasonable assumption for low-J rotational levels of CO transitions. We consider potential laws for the density (n) and rotational temperature (T). Additionally, we also consider Keplerian dynamics in the rotating disk (V rot K ) and radial expansion in the extended component (V exp ). We must highlight that our code produces results that can be quantitatively compared to observations.

AC Her
Our proposed model for the structure of AC Her (Figure 7; see also Figures 1 and 2) is very similar to the one found for the Red Rectangle, IRAS 08544−4431, and IW Car (see [12,15,16]). The total mass of the nebula is 8.3 × 10 −4 M . Our model predicts that the mass of the outflow must be ≤5% of the total mass. Thus, AC Her is clearly a binary post-AGB star surrounded by a disk-dominated nebula, due to the fact that the mass of the Keplerian disk is, at least, 19 times larger than that of the outflow. The Keplerian rotation velocity field of the disk is compatible with a central total stellar mass of ∼1 M .

R Sct
Our proposed model for the structure of R Sct is presented (Figure 7; see also Figures 3 and 4). There are slight differences between the PV diagrams, but all of them are accounted for in our uncertainties. We note that our best-fit model cannot be very different from other models. The nature of R Sct is not yet clear, but our interferometric maps firmly suggest that this source is also a binary post-AGB star surrounded by a disk with Keplerian dynamics and by a high-mass extended and expanding component. The mass of the nebula is found to be ∼ 3.2 × 10 −2 M and approximately ∼25% of the nebular material would be placed in the rotating disk. This fact, together with the large size of the outflow, allows us to classify this source as an outflowdominated post-AGB nebula. The disk with Keplerian dynamics is compatible with a central stellar mass of 1.7 M . Figure 7. Density distribution and structure of our best-fit model for the disk and outflow of AC Her, 89 her, and R Sct as seen from the observer. All models are shown with the same scale so that we can see their relative sizes.

89 Her
The total mass of the nebula around 89 Her is 1.4 × 10 −2 M and our proposed model predicts that the mass of the hourglass must be ∼50% of the total mass (Figure 7; see also Figures 5 and 6). Thus, this source is in between the disk-and outflow-dominated sources. We find that the disk with Keplerian dynamics is compatible with a central stellar mass of 1.7 M .

First Molecular Survey in Binary Post-AGB Stars
The chemistry of this kind of binary post-AGB source with rotating disks is practically unknown. We present a very deep and wide survey of radio lines in ten of our sources: AC Her, the Red Rectangle, 89 Her, HD 52961, IRAS 19157−0257 IRAS 18123+0511, IRAS 19125+0343, AI CMi, IRAS 20056+1834, and R Sct. All of them have been observed at 7 and 13 mm, and most of them have been also observed in the 1.3, 2, and 3 mm bands; see Table 1 (see [2]).

Molecular Richness
We show in Figure ?? integrated intensity ratios between the main rare molecules (SO, SiO, SiS, CS, and HCN) and CO ( 13 CO J = 2 − 1 and 12 CO J = 1 − 0). Additionally, we also compare these molecular integrated intensities with the 12, 25, and 60 µm IR emission. The averaged of our results are represented with black horizontal lines. We compare our results with the molecular emission of AGB stars (blue and red horizontal lines represent averaged values of the molecular emission for O-and C-rich AGB stars, respectively). We note that the average of our results always presents low molecular emission in molecules other than CO. Note the large range (logarithmic scale) of intensity ratios. These low intensities are more remarkable in the disk-dominated sources, such as AC Her and the Red Rectangle (see [1,2]).  [17,18]). Our sample of nebulae around binary post-AGB stars clearly presents low molecular emission in molecules other than CO.

The Discrimination between O-and C-Rich Envelopes
Evolved stars present a O/C > 1 and O/C < 1 chemistry, so they can be classified as O-and C-rich environments, respectively. The O/C abundance ratio has important effects on the molecular abundances. The lines of O-bearing molecules are much more intense in O-rich environments than in the C-rich ones (such as SiO and SO). On the contrary, the lines of C-bearing molecules are much more intense in C-rich environments (see, e.g., [17,18]). Additionally, SiO and H 2 O maser emission is exclusive to O-rich environments (see, e.g., [19]).
We analyze the integrated intensities of pairs of molecular transitions, and this analysis is crucial to distinguish between O-and C-rich environments (see Figure 9). When an O-bearing molecule is compared with a C-bearing one, we find that the integrated intensities are larger in Orich than in C-rich sources. Our results are compared with CSEs around AGB stars, because they are prototypical environments rich in molecules. Based on the maser detection of O-bearing molecules (SiO maser emission in R Sct, AI CMi, and IRAS 20056+1834; H 2 O maser emission in R Sct and AI CMi), we classify some of our sources as O-rich. On the contrary, and based on the integrated intensity ratios, we classify the nebula around 89 Her as C-rich (see Figure 9). Therefore, the nebula around AC Her, the Red Rectangle, AI CMi, IRAS 20056+1834, and R Sct presents a O/C > 1 chemistry, while 89 Her presents a O/C < 1 environment [2].

Conclusions
There is a class of post-AGB star that is part of a binary system with a significant NIR excess that is surrounded by a disk with Keplerian dynamics and an extended and expanding component composed of gas escaping from the disk and surrounding it.
Based on our observational data and model results, we find disk-dominated sources that present ≥85% of the total nebular mass located in the Keplerian disk. This is the case of AC Herculis. We also find a subclass of these binary post-AGB stars, in which the disk contains ∼25% of the total mass of the nebula, such as R Scuti. The extended components of these outflow-dominated sources are mainly composed of cold gas. Moreover, our NOEMA maps and modeling suggest that the nebula around 89 Her is in an intermediate case between both the disk-and the outflow-dominated sources, since around 50% of the nebular mass is located in the rotating disk. See Section 3 for further details. HD 52961 and IRAS 1957−0247 would also belong to this intermediate case. However, the existence of this intermediate type is not clear, because these objects were classified as intermediate sources under high uncertainties and they could belong to either subclass: the disk-or the outflow-dominated sources. In the case of 89 Her, our new 30 m IRAM on-the-fly observations recover all the filtered flux. These maps show a larger hourglass-like structure compared to that in the NOEMA maps. According to these new maps and preliminary results, the hourglass-like structure around 89 Her could contain most of the material (Gallardo Cava et al., in prep).
We present the first survey in the search for molecules other than CO in binary post-AGB stars surrounded by Keplerian disks (see Section 4). The emission of molecules other than CO in our sources is low and this fact is especially remarkable in the disk-dominated nebulae. Additionally, and according to our analysis, we catalog the chemistry of 89 Her as C-rich. On the contrary, we find O-rich environments in AC Her, the Red Rectangle, AI CMi, IRAS 20056+1834, and R Sct.