ALMA ’ s Acute View of pPNe : Through the Magnifying Glass . . . and What We Found There

We present recent Atacama Large Millimeter/submillimeter Array (ALMA)-based studies of circumstellar envelopes (CSEs) around Asymptotic Giant Branch (AGB) stars and pre-Planetary Nebulae (pPNe). In only a few years of operation, ALMA is revolutionising the field of AGB-to-PN research by providing unprecedentedly detailed information on the complex nebular architecture (at large but also on small scales down to a few ∼10 AU from the centre), dynamics and chemistry of the outflows/envelopes of low-to-intermediate mass stars in their late stages of the evolution. Here, we focus on continuum and molecular line mapping studies with high angular resolution and sensitivity of some objects that are key to understanding the complex PN-shaping process. In particular, we offer (i) a brief summary of ALMA observations of rotating disks in post-AGB objects and (ii) report on ALMA observations of OH 231.8+4.2 providing the most detailed and accurate description of the global nebular structure and kinematics of this iconic object to date.


Introduction
Due to its unique capabilities, the Atacama Large Millimeter/submillimeter Array (ALMA) has an immense potential to make great advances, and to answer major questions, in the field of Asymptotic Giant Branch (AGB)-to-Planetary Nebulae (PNe) evolution.The first ALMA images of an AGB star, R Sculptoris [1], showed a spiral structure inscribed in the circumstellar envelope (CSE), and, for the first time, astronomers could get full three-dimensional information about this outrageous structure, created by a hidden companion star orbiting the red giant.These ALMA maps were the prelude for future discoveries and marked the beginning of a new era in AGB/post-AGB/pre-Planetary Nebulae (pPNe)/PNe research.
At the time of giving this presentation, i.e., during the Asymmetric Planetary Nebulae (APN) VII meeting, ALMA was in cycle five, offering already excellent angular resolution over a broad range of frequencies (from ∼100 to ∼870 GHz).At 1 mm the angular resolution was below 20 milliarcsec, which represents an improvement of more than a factor of 10 with respect to other mm/sub-mm-interferometers.This resolution enables us (i) to obtain a close up view of the central regions of pPNe and to look for relevant structures that are directly related to, or take part in, the wind collimation process and (ii) to gather a very detailed characterization of large-scale structures (and to discover new ones) that can help us to build up the nebular shaping history of these objects.In this presentation, we show some examples of these two types of studies.In Section 2, a very brief summary of ALMA observations of rotating disks in post-AGB objects is offered.In Section 3, we present (a small part of) our ALMA observations of OH 231.8+4.2, a key object to understand PN-shaping.

A Quest for Rotating Disks in Post-AGB Objects
Rotating disks are invoked by most wind collimation theories, e.g., [2][3][4], and are suspected to be present in a particular class of post-AGB stars with near-IR emission excess (see [5] for a review).The near-IR emission excess indicates hot dust located in a stable disk-like structure close to the star (at ≈10 14 -10 15 cm).In contrast to most pPNe, these objects show narrow CO emission profiles consistent with rotating disks [6].The prototype of this class of post-AGB objects is the Red Rectangle, which was the first source in which a Keplerian disk was directly confirmed and spatially resolved with interferometric CO emission observations [7].Modelling of these data enabled a detailed characterization of the structure and kinematics of the disk, as well as an estimate of the mass of the central binary system of ∼1.5 M .
More recently, several molecular transitions in the Red Rectangle have been mapped with ALMA with better angular resolution and sensitivity [8,9].This has enabled not only an improved characterization of the disk, but has also revealed the presence of a tenuous X-shaped wind that emanates from the disk at the low velocity of ∼3-10 km s −1 .The mass-loss rate of the wind, which may be material photoevaporating from the disk, implies a lifetime of the disk of ∼5000-10,000 yr.
The quest for rotating disks continues, and more post-AGB objects with near-IR excess are being observed in CO emissions.IW Car [10] and IRAS 08544-4133 [11] are the most recent examples of post-AGB stars in which rotating disks and winds emerging from the disk have been spatially resolved with ALMA.Rotating disks around post-AGB stars had been mapped earlier, also with the IRAM interferometer, e.g., 89 Her [12] and AC Her [13].
The properties of the disks and the X-winds (when detected) are similar in all post-AGB targets (Table 1), although a very low-number statistics still prevails.We emphasize that, given the large dimensions of the disks mapped to date (with typical outer radii of the order of 1000 AU), these structures do not represent the compact accretion disks (around the compact companion of the mass-lossing star) postulated by most wind-collimation theories.The rotating disks reported as of today are relatively large structures around the central binary system, typically with a total mass of M ∼ 1-2 M .The role of these circumbinary disks in the wind collimation and/or, more generally, in the PN-shaping process, remains unclear.Recently, the first rotating disk in an AGB star, L 2 Pup, has been detected and spatially resolved with ALMA [14,15].In this case, the disk is smaller and less massive than the post-AGB stars shown before, but it is also tentatively circumbinary.
Rotating disks in post-AGB objects have typical masses of M disk ≈ 10 −3 -10 −2 M , the Red Rectangle being the most massive detected so far, and excitation temperatures of ≈100 K.For objects where a rotating disk and an X-wind component are simultaneously detected, the mass of the latter is only a small fraction (10%) of the mass of the disk.These X-winds are very different from the massive (∼0.1-1M ) and high-momentum bipolar outflows common to most pPNe [16].It is paradoxical that none of the pPNe with massive, high-momentum bipolar outflows have shown clear evidence for rotation at their cores.Therefore, ironically, rotating disk are only found in objects where a powerful jet-launching engine does not need to be invoked.This is a mystery that ALMA will hopefully solve in the coming years, as more rotating disks and pPNe are observed.Characterization with ALMA of the nebular structure and dynamics of 'standard' pNe, i.e., with dominant bipolar outflows carrying a significant amount of mass and linear momentum, is essential to understand how they were formed.We show ALMA observations of OH 231.8+4.2, an object in a brief but key evolutionary stage of the AGB-to-PN transition.
OH 231.8+4.2 is an outstanding pPN-like nebula (Figure 1) around a Mira-type pulsating star (QX Pup).This uncommon situation makes OH 231.8+4.2 of special interest, since it demonstrates that asymmetries can develop (and become very prominent) while the central star is still on the AGB.OH 231.8+4.2 has been extensively studied at many wavelengths by several authors e.g., [17][18][19][20][21][22][23][24] (and references therein).The light from the central source inside OH 231.8's core is highly obscured by dust along the line-of-sight and is only seen indirectly scattered by the dust in the lobe walls, indicating the presence of an A 0-type main-sequence companion to the mass-losing AGB star, of spectral type M8-10 III [21].It is believed that the bulk of the nebular mass, which is in the form of molecular gas (∼1 M ), has been ejected by the AGB star at a very high rate, Ṁ ≈ 10 −4 M yr −1 [17].Before ALMA came into play, two major large-scale components had been identified: (i) An equatorial waist expanding at low velocity, <30 km s −1 and (ii) a highly collimated bipolar outflow with expansion velocities that increase linearly with the distance from the centre up to ∼400 km s −1 (deprojected).The linear velocity gradient in the outflow could have resulted from the sudden interaction between collimated fast winds (CFWs) on the ambient AGB material ∼800 yr ago.
We present a summary of our recent observations of OH 231.8+4.2 with ALMA.These data have provided us with a high-definition image of the dominant nebular components described above, but also, more importantly, have enabled us (1) to identify the position of the mass-losing star relative to the nebula and (2) to discover a number of outflow components previously unknown.We also report the first detection of Na 37 Cl and CH 3 OH in OH 231.8+4.2, with CH 3 OH also being a first detection around an AGB star.A comprehensive study of this dataset is reported by [25].

Observations
OH 231.8+4.2 was observed with the ALMA 12 m array as part of project 2015.1.00256.S on July 2016.We mapped the continuum emission and a series of spectral lines in band seven (∼294-345 GHz).Our project has two major science goals: SG1 is a five-point mosaic of 12 CO/ 13 CO (3-2) and other molecules covering a ∼19 ×54 area along the bipolar outflow, and SG2 is a single-pointing observation towards the centre to map a selection of lines, e.g., SiO and CS, amongst others.The resolution of these observations is ∼0.2-0. 3. Full observational details are given in [25].

Continuum Maps
Maps of the continuum emission at four different frequencies (294, 304, 330 and 344 GHz) have been created using line-free channels.Our ∼294 GHz-continuum map is shown in the left panel of Figure 2 as an example.The surface brightness distribution, which is very similar at all frequencies, appears as an extended, incomplete hourglass-like structure, roughly oriented along PA ∼ 21 • (i.e., the symmetry axis of the large-scale nebula).The continuum emission distribution is clumpy.The peak of the continuum maps is attained at one of these clumps, referred to as "clump S", which has coordinates R.A. = 07 h 42 m 16. s 915 and Dec. = −14 • 42 50 06 (J2000).As we will show below, clump S enshrouds the central star QX Pup.We deduce a deconvolved radius of ∼40-70 AU for clump S, which is barely resolved.Note that clump S does not lie on the equatorial plane of the hourglass, but is clearly displaced along the axis towards the south by ∼0. 6.The total continuum flux with ALMA is consistent with that measured with single-dish telescopes (see right panel of Figure 2).Except for clump S, the spectral index of the continuum shows no significant deviations from a S ν ∝ ν 3.3  power-law across the different regions.This implies a dominant contribution by optically thin ∼75 K dust with an emissivity index of α ∼1.3.At clump S, the continuum obeys a S ν ∝ ν 2.1 dependence consistent with a population of large (>100 µm-sized) grains in the vicinity of QX Pup.North is up.Two faint, bubble-like features (fish-bowls) with molecular counterparts in our ALMA data are indicated-see Section 3.3.(Right) ALMA order-zero moment maps of 12 CO (2-1) integrated over the full width of the line profile adapted from [17].The mass and expansion velocity of the low-velocity core and the fast, bipolar lobes are indicated.The distance to OH 231.8+4.2 is d∼1500 pc.The lobes are tilted with respect to the plane of the sky by i∼35 • , with the North lobe being the closest.

Molecular Line Maps: Nebular Components
We have mapped a large number of molecular transitions, covering a range of upper-level energies from E u ∼ 20 to 1800 K and including many species (CO, CS, SO, SO 2 , SiO, SiS, H 3 O + , CH 3 OH, etc.).Different lines are found to selectively (or exclusively) trace different components of the outflow, which facilitates disentangling the complex nebular architecture of OH 231.8+4.2.In the following, we describe different structures traced by ALMA, starting from the most compact ones (within ≈10 15 cm from the mass-losing star) and ending with the most extended zones (reaching out to ≈10 17 cm).
We have detected several high-excitation lines, namely, SiS (v = 1, J = 17 − 16), 30 SiO (v = 1, J = 7 − 6), and Na 37 Cl (v = 0, J = 26 − 25), with compact emission arising entirely from clump S.These lines are known to be produced in the warm inner ( < ∼ 10 15 cm) layers of the winds of evolved mass-losing stars, where the stellar wind may have not reached its terminal expansion velocity.The emission is indeed confined to a small emitting volume of radius ∼60 AU (∼24R ), comparable to that of the continuum-emitting core.The narrow profiles imply low expansion velocities of V exp ∼ 5-7 km s −1 in these layers, for which we deduce a kinematic age of about 50 yr.In addition to disclosing the locus of QX Pup inside clump S, we have discovered a compact (∼1 × 4 ) bipolar outflow that emerges from the stellar vicinity.This outflow is exclusively sampled in our data by the SiO molecule, which is a well-known shock tracer (top panels of Figure 3).Our maps show a pair of flame-shaped lobes expanding at moderate velocities and oriented similarly to the large-scale nebula.The radial (or line-of-sight) velocity increases abruptly from the centre to two bright and compact regions located at offsets δy ∼ ±0.11 about clump S, where the maximum radial velocities (∼35 km s −1 ) and full line widths (∼23 km s −1 ) are observed.These SiO-bright regions may represent bow-shocks, denoting recent ( < ∼ 50-80 yr old) bipolar mass ejections.The lower velocities at the intermediate-to-outer regions of the flame-shaped SiO-lobes imply larger kinematic ages of t kin ∼ 400-500 yr.
Adjacent to the SiO-outflow, we identify a small-scale hourglass-shaped structure (mini-hourglass) that is probably made of compressed material formed as the SiO-outflow penetrates the dense, central regions of the nebula (bottom panels of Figure 3).The lobes and the equatorial waist of the mini-hourglass both are radially expanding with a constant velocity gradient (V exp ∝ r).The dimensions and kinematics of the mini-waist are consistent with a torus-like structure with an outer radius of about 500 AU orthogonal to, and coeval with, the SiO-outflow.The expansion velocity at the inner edge of the mini-waist (∼150 AU) is extremely low, at V exp ∼ 1 km s −1 .The ∼8 × 4 -sized hourglass nebula traced by the dust thermal continuum emission maps is the major emission component of most of the species mapped by us with ALMA.In Figure 4, we offer a condensed view of CS (7−6), one of the dense gas tracers observed by us, where this component is clearly seen.The large-hourglass is characterized by a dominant expansive kinematics described by a radial velocity gradient, ∇v ∼ 6.0-6.5 km s −1 arcsec −1 , that is sustained from the base to the tips of the lobes, where the highest velocities are observed.The sharp outer boundary of the large-waist is now accurately determined (at a radius of ∼2700 AU) as well as its internal velocity field, indicating that the waist was shaped nearly simultaneously with the high-velocity (HV) lobes, about 800-900 yr ago.
We report the discovery of two large (∼8 × 6 ) faint, bubble-like structures surrounding the central parts of the nebula (fish-bowls, see Figure 5).These are the molecular counterparts of the two faint, rounded structures seen in the HST/NIR images (Figure 1).The north and the south fish-bowls have similar dimensions, but are oriented differently.The observed spatio-kinematic distribution is consistent with the fish-bowls being a pair of hollow, thin-walled ellipsoids (probably made of swept up ambient material) radially expanding and moving away from each other.They are relatively old structures although apparently slightly (∼100-200 yr) younger than the large-scale waist and the HV-lobes.The origin of these enigmatic structures is unknown, but it may be linked to a 'critical' time period or instant of the evolution/development of the HV-lobes.

Conclusions
Our ALMA observations of OH 231.8+4.2 unveil a series of new substructures that point to a nebular formation/shaping history significantly more complex than previously thought, and in particular, indicative of multiple non-spherical mass ejections.The origin of the bipolar ejections that led to the formation of the outstanding nebular architecture of OH 231.8+4.2 remains unknown, but the presence of present-day bipolar ejections indicate that the CFW engine is still active at its core.
There are notable differences in the velocities and overall symmetry between the large-scale ∼800 yr old CO-outflow and the most current ( < ∼ 80 yr old) bipolar ejections traced by SiO.These differences could indicate that the jet-launching mechanism itself has changed, or that the initial conditions under which the jets formed (presumably, involving disk-mediated mass transfer from the AGB to the companion star that blows the jets) have been modified in such a way that they now produce slower and more symmetric bipolar ejections than in the past.These changes may be intimately linked to major adjustments in the binary system configuration, which may have occurred after the powerful ejections that resulted in the formation of the massive and fast CO-outflow.
The position of QX Pup off-centre from the waist of the large-hourglass is one of the most puzzling discoveries from our ALMA data.Perhaps the combination of orbital motion and recoil of the binary system after strong asymmetrical mass ejections could explain, at least partially, its mysterious location.The orbital parameters of the central binary system are very poorly known.A loose upper limit to the orbital separation of a < 150 AU is deduced from the relative separation between clump S and the centroid of the SiO-outflow (presumably collimated by the companion).
New ∼0.05-resolution observations with ALMA by our team (2017.1.00706.S) will enable us to determine the precise location of the SiO outflow engine relative to QX Pup and the scale at which the wind collimation begins, bringing us to the closest we have ever been to watching live the launch of a bipolar outflow from a mass-losing star and to understanding how OH 231.8+4.2 assembled its complex nebular architecture.

3 .
Nebular Architecture and Dynamics of pPNe: The ALMA View of OH 231.8+4.2

Figure 1 .
Figure 1.(Left) False color HST image of the reflection nebulosity (red-to-yellow) and Hα-emitting lobes (blue) of OH 231.8+4.2 (Credit: ESA/Hubble & NASA.Acknowledgement: Judy Schmidt).North is up.Two faint, bubble-like features (fish-bowls) with molecular counterparts in our ALMA data are indicated-see Section 3.3.(Right) ALMA order-zero moment maps of 12 CO (2-1) integrated over the full width of the line profile adapted from [17].The mass and expansion velocity of the low-velocity core and the fast, bipolar lobes are indicated.The distance to OH 231.8+4.2 is d∼1500 pc.The lobes are tilted with respect to the plane of the sky by i∼35 • , with the North lobe being the closest.

Figure 2 .
Figure 2. (Left) ALMA 294 GHz-continuum emission map of OH 231.8+4.2.Contour level spacing is 0.8 mJy/beam.The clean beam (HPBW = 0. 31 × 0. 25, PA = −84.5 • ) is plotted at the bottom-left corner.The compact region at the centre where the continuum emission peaks is referred to as clump S and it is a dust (and gas) region enshrouding the AGB star QX Pup.(Right) Spectral energy distribution.Pink symbols are our ALMA continuum measurements (squares = integrated flux; stars = peak-flux surface density at clump S) and the dotted and dashed lines are S ν ∝ ν 3.3 and S ν ∝ ν 2.1 fitted to those.

Figure 4 .
Figure 4. Overview of the CS (7−6) emission maps of OH 231.8+4.2 showing the integrated intensity map, atop the HST/F110W image (left), the one-dimensional spectral profile (middle), and the axial PV diagram (right).Different nebular components are indicated (namely, the large-hourglass, the HV-lobes, and the mini-hourglass).See text for details.

Figure 5 .
Figure 5. Velocity-channel maps of 12 CO (3−2) in the velocity range where two faint, elongated bubble-like features (dubbed the fish-bowls) are discovered in the central regions of OH 231.8.