Source author record

V. Bujarrabal

V. Bujarrabal appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

45works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

45 published item(s)

preprint2022arXiv

Chemistry of nebulae around binary post-AGB stars: A molecular survey of mm-wave lines

Context. There is a class of binary post-asymptotic giant branch (post-AGB) stars that exhibit remarkable near-infrared (NIR) excess. Such stars are surrounded by Keplerian or quasi-Keplerian disks, as well as extended outflows composed of gas escaping from the disk. This class can be subdivided into disk- and outflow-dominated sources, depending on whether it is the disk or the outflow that represents most of the nebular mass, respectively. The chemistry of this type of source has been practically unknown thus far. Methods. We focused our observations on the 1.3, 2, 3 mm bands of the 30 m IRAM telescope and on the 7 and 13 mm bands of the 40 m Yebes telescope. Our observations add up around 600 hours of telescope time. Results. We present the first single-dish molecular survey of mm-wave lines in nebulae around binary post-AGB stars. We conclude that the molecular content is relatively low in nebulae around binary post-AGB stars, as their molecular lines and abundances are especially weaker compared with AGB stars. This fact is very significant in those sources where the Keplerian disk is the dominant component of the nebula. The study of their chemistry allows us to classify nebulae around AC Her, the Red Rectangle, AI CMi, R Sct, and IRAS 20056+1834 as O-rich, while that of 89 Her is probably C-rich. The calculated abundances of the detected species other than CO are particularly low compared with AGB stars. The initial stellar mass derived from the 17O/18O ratio for the Red Rectangle and 89 Her is compatible with the central total stellar mass derived from previous mm-wave interferometric maps. The very low 12CO/13CO ratios found in binary post-AGB stars reveal a high 13CO abundance compared to AGB and other post-AGB stars.

preprint2022arXiv

Dissecting the central regions of OH 231.8 + 4.2 with ALMA: a salty rotating disk at the base of a young bipolar outflow

We present ALMA continuum and molecular line emission maps at $\sim$1 mm of OH 231.8, a well studied bipolar nebula around an AGB star. The excellent angular resolution of our maps ($\sim$20 mas) allows us to scrutinise the central nebular regions of OH 231.8, which hold the clues to unravel how this iconic object assembled its complex nebular architecture. We report, for the first time in this object and others of its kind, the discovery of a rotating circumbinary disk selectively traced by NaCl, KCl, and H$_2$O emission lines. This represents the first detection of KCl in an oxygen-rich AGB circumstellar envelope. The rotating disk, of radius $\sim$30 au, lies at the base of a young bipolar wind traced by SiO and SiS emission, which also presents signs of rotation at its base. The NaCl equatorial structure is characterised by a mean rotation velocity of $\sim$4 km s$^{-1}$ and extremely low expansion speeds, $\sim$3 km s$^{-1}$. The outflow has a predominantly expansive kinematics characterised by a constant radial velocity gradient of $\sim$65 km s$^{-1}$ arcsec$^{-1}$ at its base. Beyond $r$$\sim$350 au, the gas in the outflow continues radially flowing at a constant terminal speed of $\sim$16 km s$^{-1}$. Our continuum maps reveal a spatially resolved dust disk-like structure perpendicular to the outflow, with the NaCl, KCl and H$_2$O emission arising from the disk's surface layers. Within the disk, we also identify an unresolved point continuum source, which likely represents the central Mira-type star QX Pup enshrouded by a $\sim$3 $R_{\star}$ component of hot ($\sim$1400 K) freshly formed dust. The point source is slightly off-centered from the disk centroid, enabling us for the first time to place constraints to the orbital separation and period of the central binary system, $a$$\sim$20 au and $P_{\rm orb}$$\sim$55 yr, respectively. (abridged).

preprint2022arXiv

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.

preprint2021arXiv

The ionised and molecular mass of post-common-envelope planetary nebulae. The missing mass problem

Most planetary nebulae (PNe) show beautiful, axisymmetric morphologies despite their progenitor stars being essentially spherical. Close binarity is widely invoked to help eject an axisymmetric nebula, after a brief phase of engulfment of the secondary within the envelope of the Asymptotic Giant Branch (AGB) star, known as the common envelope (CE). The evolution of the AGB would thus be interrupted abruptly, its envelope being rapidly ejected to form the PN, which a priori would be more massive than a PN coming from the same star, were it single. We aim at testing this hypothesis by investigating the mass of a sample of 21 post-CE PNe, ~1/5th of the known total population, and comparing them to a large sample of `regular' (i.e. not known to host close binaries) PNe. We have gathered data on the ionised and molecular content of our sample and carried out new molecular observations. We derive the ionised and molecular masses of the sample by means of a systematic approach, using tabulated, dereddened H-beta fluxes for finding the ionised mass, and CO 2-1 and 3-2 observations for the molecular mass. There is a general lack of molecular content in post-CE PNe, with few exceptions. Once we derive the ionised and molecular masses, we find that post-CE PNe arising from Single-Degenerate (SD) systems are just as massive, on average, as `regular' PNe, whereas post-CE PNe arising from Double-Degenerate (DD) systems are considerably more massive, and show larger linear momenta and kinetic energy than SD systems and `regular' PNe. Reconstruction of the CE of four objects suggests that the mass of SD nebulae actually amounts to a very small fraction of the envelope of their progenitor stars. This leads to the uncomfortable question of where the rest of the envelope is and why we cannot detect it in the stars' vicinity, thus raising serious doubts on our understanding of these intriguing objects.

preprint2020arXiv

Molecular remnant of Nova 1670 (CK Vulpeculae). II. A three-dimensional view on the gas distribution and velocity field

CK Vul is the remnant of an energetic eruption known as Nova 1670 that is thought to be caused by a stellar merger. The remnant is composed of (1) a large hourglass nebula of recombining gas (of 71\arcsec size), very similar to some classical planetary and pre-planetary nebulae (PPNe), and (2) of a much smaller and cooler inner remnant prominent in millimeter-wave emission from molecules. We investigate the three-dimensional spatio-kinematic structure of both components. The analysis of the hourglass structure yields a revised distance to the object of >2.6 kpc, at least 3.7 times greater than so far assumed. At this distance, the stellar remnant has a bolometric luminosity >12 L$_{\odot}$ and is surrounded by molecular material of total mass >0.8 M$_{\odot}$ (the latter value has a large systematic uncertainty). We also analyzed the architecture of the inner molecular nebula using ALMA observations of rotational emission lines obtained at subarcsecond resolution. We find that the distribution of neutral and ionized gas in the lobes can be reproduced by several nested and incomplete shells or jets with different velocity fields and varying orientations. The analysis indicates that the molecular remnant was created in several ejection episodes, possibly involving an interacting binary system. We calculated the linear momentum ($\approx$10$^{40}$ g cm s$^{-1}$) and kinetic energy ($\approx$10$^{47}$ erg) of the CK Vul outflows and find them within the limits typical for classical PPNe. Given the similarities of the CK Vul outflows to PPNe, we suggest there may CK Vul analogs among wrongly classified PPNe with low intrinsic luminosities, especially among PPNe with post-red-giant-branch central stars.

preprint2020arXiv

Yebes 40 m radio telescope and the broad band NANOCOSMOS receivers at 7 mm and 3 mm for line surveys

Yebes 40\,m radio telescope is the main and largest observing instrument at Yebes Observatory and it is devoted to Very Long Baseline Interferometry (VLBI) and single dish observations since 2010. It has been covering frequency bands between 2\,GHz and 90\,GHz in discontinuous and narrow windows in most of the cases, to match the current needs of the European VLBI Network (EVN) and the Global Millimeter VLBI Array (GMVA). Nanocosmos project, a European Union funded synergy grant, opened the possibility to increase the instantaneous frequency coverage to observe many molecular transitions with single tunnings in single dish mode. This reduces the observing time and maximises the output from the telescope. We present the technical specifications of the recently installed 31.5-50 GHz (Q band) and 72-90.5 GHz (W band) receivers along with the main characteristics of the telescope at these frequency ranges. We have observed IRC+10216, CRL 2688 and CRL 618, which harbour a rich molecular chemistry, to demonstrate the capabilities of the new instrumentation for spectral observations in single dish mode. The results show the high sensitivity of the telescope in the Q band. The spectrum of IRC+10126 offers a signal to noise ratio never seen before for this source in this band. On the other hand, the spectrum normalised by the continuum flux towards CRL\,618 in the W band demonstrates that the 40~m radio telescope produces comparable results to those from the IRAM 30~m radio telescope, although with a smaller sensitivity. The new receivers fulfil one of the main goals of Nanocosmos and open the possibility to study the spectrum of different astrophysical media with unprecedented sensitivity.

preprint2018arXiv

On the circumstellar envelopes of semi-regular long-period variables

The mass loss process along the AGB phase is crucial for the formation of circumstellar envelopes (CSEs), which in the post-AGB phase will evolve into planetary nebulae (PNe). There are still important issues that need to be further explored in this field; in particular, the formation of axially symmetric PNe from spherical CSEs. To address the problem, we have conducted high S/N IRAM 30m observations of 12CO J=1-0 and J=2-1, and 13CO J=1-0 in a volume-limited unbiased sample of semi-regular variables (SRs). We also conducted Yebes 40m SiO J=1-0 observations in 1/2 of the sample in order to complement our 12CO observations. We report a moderate correlation between mass loss rates and the 12CO(1-0)-to-12CO(2-1) intensity ratio, introducing a possible new method for determining mass loss rates of SRs with short analysis time. We also find that for several stars the SiO profiles are very similar to the 12CO profiles, a totally unexpected result unless these are non-standard envelopes.

preprint2016arXiv

28 SiO v=0 J=1-0 emission from evolved stars

Observations of 28SiO v=0 J=1-0 line emission (7-mm wavelength) from AGB stars show in some cases peculiar profiles, composed of a central intense component plus a wider plateau. Very similar profiles have been observed in CO lines from some AGB stars and most post-AGB nebulae and, in these cases, they are clearly associated with the presence of conspicuous axial symmetry and bipolar dynamics. We present systematic observations of 28SiO v=0 J=1-0 emission in 28 evolved stars, performed with the 40~m radio telescope of the IGN in Yebes, Spain. We find that the composite core plus plateau profiles are almost always present in O-rich Miras, OH/IR stars, and red supergiants. They are also found in one S-type Mira ($χ$ Cyg), as well as in two semiregular variables (X Her and RS Cnc) that are known to show axial symmetry. In the other objects, the profiles are simpler and similar to those of other molecular lines. The composite structure appears in the objects in which SiO emission is thought to come from the very inner circumstellar layers, prior to dust formation. The central spectral feature is found to be systematically composed of a number of narrow spikes, except for X Her and RS Cnc, in which it shows a smooth shape that is very similar to that observed in CO emission. These spikes show a significant (and mostly chaotic) time variation, while in all cases the smooth components remain constant within the uncertainties. The profile shape could come from the superposition of standard wide profiles and a group of weak maser spikes. Alternatively, we speculate that the very similar profiles detected in objects that are axisymmetric may be indicative of the systematic presence of a significant axial symmetry in the very inner circumstellar shells around AGB stars; the presence of such symmetry would be independent of the probable weak maser effects in the central spikes.

preprint2016arXiv

A lambda 3mm and 1mm line survey toward the yellow hypergiant IRC +10420

Our knowledge of the chemical properties of the circumstellar ejecta of the most massive evolved stars is particularly poor. We aim to study the chemical characteristics of the prototypical yellow hypergiant star, IRC +10420. For this purpose, we obtained full line surveys at 1 and 3mm atmospheric windows. We have identified 106 molecular emission lines from 22 molecular species. Approximately half of the molecules detected are N-bearing species, in particular HCN, HNC, CN, NO, NS, PN, and N2H+. We used rotational diagrams to derive the density and rotational temperature of the different molecular species detected. We introduced an iterative method that allows us to take moderate line opacities into account. We have found that IRC +10420 presents high abundances of the N-bearing molecules compared with O-rich evolved stars. This result supports the presence of a N-rich chemistry, expected for massive stars. Our analysis also suggests a decrease of the 12C/13C ratio from \gtrsim 7 to \sim 3.7 in the last 3800 years, which can be directly related to the nitrogen enrichment observed. In addition, we found that SiO emission presents a significant intensity decrease for high-J lines when compared with older observations. Radiative transfer modeling shows that this variation can be explained by a decrease in the infrared (IR) flux of the dust. The origin of this decrease might be an expansion of the dust shell or a lower stellar temperature due to the pulsation of the star.

preprint2016arXiv

A second post-AGB nebula that contains gas in rotation and in expansion: ALMA maps of IW Car

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.

preprint2016arXiv

ALMA high spatial resolution observations of the dense molecular region of NGC 6302

The mechanism behind the shaping of bipolar planetary nebulae is still poorly understood. Accurately tracing the molecule-rich equatorial regions of post-AGB stars can give valuable insight into the ejection mechanisms at work. We investigate the physical conditions, structure and velocity field of the dense molecular region of the planetary nebula NGC 6302 by means of ALMA band 7 interferometric maps. The high spatial resolution of the $^{12}$CO and $^{13}$CO J=3-2 ALMA data allows for an analysis of the geometry of the ejecta in unprecedented detail. We built a spatio-kinematical model of the molecular region with the software SHAPE and performed detailed non-LTE calculations of excitation and radiative transfer with the SHAPEMOL plug-in. We find that the molecular region consists of a massive ring out of which a system of fragments of lobe walls emerge and enclose the base of the lobes visible in the optical. The general properties of this region are in agreement with previous works, although the much greater spatial resolution of the data allows for a very detailed description. We confirm that the mass of the molecular region is 0.1 M$_{\odot}$. Additionally, we report a previously undetected component at the nebular equator, an inner, younger ring inclined $\sim$60$^\circ$ with respect to the main ring, showing a characteristic radius of 7.5$\times$10$^{16}$ cm, a mass of 2.7$\times$10$^{-3}$ M$_{\odot}$, and a counterpart in optical images of the nebula. This inner ring has the same kinematical age as the northwest optical lobes, implying it was ejected approximately at the same time, hundreds of years after the ejection of the bulk of the molecular ring-like region. We discuss a sequence of events leading to the formation of the molecular and optical nebulae, and briefly speculate on the origin of this intriguing inner ring.

preprint2016arXiv

Herschel/HIFI observations of the circumstellar ammonia lines in IRC+10216

New high-resolution far-infrared (FIR) observations of both ortho- and para-NH3 transitions toward IRC+10216 were obtained with Herschel, with the goal of determining the ammonia abundance and constraining the distribution of NH3 in the envelope of IRC+10216. We used the Heterodyne Instrument for the Far Infrared (HIFI) on board Herschel to observe all rotational transitions up to the J=3 level (three ortho- and six para-NH3 lines). We conducted non-LTE multilevel radiative transfer modelling, including the effects of near-infrared (NIR) radiative pumping through vibrational transitions. We found that NIR pumping is of key importance for understanding the excitation of rotational levels of NH3. The derived NH3 abundances relative to molecular hydrogen were (2.8+-0.5)x10^{-8} for ortho-NH3 and (3.2^{+0.7}_{-0.6})x10^{-8} for para-NH3, consistent with an ortho/para ratio of 1. These values are in a rough agreement with abundances derived from the inversion transitions, as well as with the total abundance of NH3 inferred from the MIR absorption lines. To explain the observed rotational transitions, ammonia must be formed near to the central star at a radius close to the end of the wind acceleration region, but no larger than about 20 stellar radii (1 sigma confidence level).

preprint2016arXiv

Imaging the dust sublimation front of a circumbinary disk

We present the first near-IR milli-arcsecond-scale image of a post-AGB binary that is surrounded by hot circumbinary dust. A very rich interferometric data set in six spectral channels was acquired of IRAS08544-4431 with the new RAPID camera on the PIONIER beam combiner at the Very Large Telescope Interferometer (VLTI). A broadband image in the \textit{H} band was reconstructed by combining the data of all spectral channels using the SPARCO method. We spatially separate all the building blocks of the IRAS08544-4431 system in our milliarcsecond-resolution image. Our dissection reveals a dust sublimation front that is strikingly similar to that expected in early-stage protoplanetary disks, as well as an unexpected flux signal of $\sim$4\% from the secondary star. The energy output from this companion indicates the presence of a compact circum-companion accretion disk, which is likely the origin of the fast outflow detected in H$α$. Our image provides the most detailed view into the heart of a dusty circumstellar disk to date. Our results demonstrate that binary evolution processes and circumstellar disk evolution can be studied in detail in space and over time.

preprint2016arXiv

The millimeter IRAM-30m line survey toward IK Tau

We carried out a millimeter wavelength line survey between 79 and 356 GHz with the IRAM-30m telescope to investigate the physical and chemical properties of the molecular envelope of the oxygen-rich AGB star IK Tau. We analysed the molecular lines detected using the population diagram technique to derive rotational temperatures and column densities. Additionally, we conducted a radiative transfer analysis of the SO$_2$ lines detected. For the first time in this source we detected rotational lines in the ground vibrational state of HCO$^+$, NS, NO, and H$_2$CO, as well as several isotopologues of molecules previously identified. We also detected several rotational lines in vibrationally excited states of SiS and SiO isotopologues, and rotational lines of H$_2$O ($ν_{\mathrm{2}}$=2). We have also increased the number of rotational lines detected of molecules that were previously identified, enabling a detailed study of the molecular abundances and excitation temperatures. IK Tau displays a rich chemistry for an oxygen-rich circumstellar envelope. We highlight the detection of NS and H$_2$CO with fractional abundances of f(NS)$\sim$10$^{-8}$ and f(H$_2$CO)$\sim$[10$^{-7}$--10$^{-8}$]. Most of the molecules display rotational temperatures between 15 and 40 K. NaCl and SiS isotopologues display rotational temperatures higher than the average (~65 K). In the case of SO$_2$ a warm component with $T_{\mathrm{rot}}$$\sim$290 K is also detected. This SO$_2$ warm component is probably arising from the inner regions of the envelope (at $\sim$8$R_{*}$) where SO$_2$ has a fractional abundance of f(SO$_2$)$\sim$10$^{-6}$. This result should be considered for future investigation of the main formation channels of this, and other, parent species in the inner winds of O-rich AGB stars, which at present are not well reproduced by current chemistry models.

preprint2015arXiv

Detection of Keplerian dynamics in a disk around the post-AGB star AC Her

So far, only one rotating disk has been clearly identified and studied in AGB or post-AGB objects (in the Red Rectangle), by means of observations with high spectral and spatial resolution. However, disks are thought to play a key role in the late stellar evolution and are suspected to surround many evolved stars. We aim to extend our knowledge on these structures. We present interferometric observations of CO J=2-1 emission from the nebula surrounding the post-AGB star AC Her, a source belonging to a class of objects that share properties with the Red Rectangle and show hints of Keplerian disks. We clearly detect the Keplerian dynamics of a second disk orbiting an evolved star. Its main properties (size, temperature, central mass) are derived from direct interpretation of the data and model fitting. With this we confirm that there are disks orbiting the stars of this relatively wide class of post-AGB objects

preprint2015arXiv

Molecular ions in the O-rich evolved star OH231.8+4.2: HCO$^+$,H$^{13}$CO$^+$ and first detection of SO$^+$, N$_2$H$^+$, and H$_3$O$^+$

OH 231.8+4.2, a bipolar outflow around a Mira-type variable star, displays a unique molecular richness amongst circumstellar envelopes (CSEs) around O-rich AGB and post-AGB stars. We report line observations of the HCO+ and H13CO+ molecular ions and the first detection of SO+, N2H+, and (tentatively) H3O+ in this source. SO+ and H3O+ have not been detected before in CSEs around evolved stars. These data have been obtained as part of a full mm-wave and far-IR spectral line survey carried out with the IRAM 30 m radio telescope and with Herschel/HIFI. Except for H3O+, all the molecular ions detected in this work display emission lines with broad profiles (FWHM 50-90 km/s), which indicates that these ions are abundant in the fast bipolar outflow of OH 231.8. The narrow profile (FWHM 14 km/s) and high critical densities (>1e6cm-3 ) of the H3O+ transitions observed are consistent with this ion arising from denser, inner (and presumably warmer) layers of the fossil remnant of the slow AGB CSE at the core of the nebula. From rotational diagram analysis, we deduce excitation temperatures of Tex 10-20 K for all ions except for H3O+, which is most consistent with Tex 100 K. Although uncertain, the higher excitation temperature suspected for H3O+ is similar to that recently found for H2O and a few other molecules, which selectively trace a previously unidentified, warm nebular component.The column densities of the molecular ions reported here are in the range Ntot [1-8]x1e13 cm-2, leading to beam-averaged fractional abundances relative to H2 of X(HCO+) 1e-8, X(H13CO+) 2e-9, X(SO+) 4e-9, X(N2H+) 2e-9, and X(H3O+) 7e-9 cm-2. We have performed chemical kinetics models to investigate the formation of these ions in OH 231.8 as the result of standard gas phase reactions initiated by cosmic-ray and UV-photon ionization. (abridged).

preprint2015arXiv

New N-bearing species towards OH231.8+4.2: HNCO, HNCS, HC$_3$N and NO

OH231.8+4.2 is a well studied O-rich CSE around an intermediate-mass evolved star that displays bipolar molecular outflows accelerated up to 400 km/s. OH231 also presents an exceptional molecular richness probably due to shock-induced chemical processes. We report the first detection in this source of HNCO, HNCS, HC3N and NO, with the IRAM-30m telescope in a sensitive mm-wavelength survey towards this target. HNCO and HNCS are also first detections in CSEs. The observed line profiles show that the emission arises in the massive central component of the envelope and at the base of the fast lobes. The NO profiles are broader than those of HNCO, HNCS, and HC3N, and most importantly, broader than the line profiles of 13CO. This indicates that the NO abundance is enhanced in the fast lobes relative to the slow, central parts. From LTE and non-LTE excitation analysis, we estimate beam-average rotational temperatures of 15-30 K (and, maybe, up to 55 K for HC3N) and fractional abundances of X(HNCO)=[0.8-1]E-7, X(HNCS)=[0.9-1]E-8, X(HC3N)=[5-7]E-9 and X(NO)=[1-2]E-6. NO is, therefore, amongst the most abundant N-bearing species in OH231. We have performed thermodynamical chemical equilibrium and chemical kinetics models to investigate the formation of these N-bearing species in OH231. The model underestimates the observed abundances for HNCO, HNCS, and HC3N by several orders of magnitude, which indicates that these molecules can hardly be products of standard UV-photon and/or cosmic-ray induced chemistry in OH231, and that other processes (e.g. shocks) play a major role in their formation. For NO, the model abundance is compatible with the observed average value. The new detections presented in this work corroborate the particularly rich chemistry of OH231, which is likely profoundly influenced by shock-induced processes, as proposed in earlier works.

preprint2015arXiv

New observations and models of circumstellar CO line emission of AGB stars in the Herschel SUCCESS programme

CONTEXT: Asymptotic giant branch (AGB) stars are in one of the latest evolutionary stages of low to intermediate-mass stars. Their vigorous mass loss has a significant effect on the stellar evolution, and is a significant source of heavy elements and dust grains for the interstellar medium. The mass-loss rate can be well traced by carbon monoxide (CO) line emission. AIMS: We present new Herschel HIFI and IRAM 30m telescope CO line data for a sample of 53 galactic AGB stars. The lines cover a fairly large range of excitation energy from the $J=1\to0$ line to the $J=9\to8$ line, and even the $J=14\to13$ line in a few cases. We perform radiative transfer modelling for 38 of these sources to estimate their mass-loss rates. METHODS: We used a radiative transfer code based on the Monte Carlo method to model the CO line emission. We assume spherically symmetric circumstellar envelopes that are formed by a constant mass-loss rate through a smoothly accelerating wind. RESULTS: We find models that are consistent across a broad range of CO lines for most of the stars in our sample, i.e., a large number of the circumstellar envelopes can be described with a constant mass-loss rate. We also find that an accelerating wind is required to fit, in particular, the higher-J lines and that a velocity law will have a significant effect on the model line intensities. The results cover a wide range of mass-loss rates ($\sim 10^{-8}$ to $2\times 10^{-5}~\mathrm{M}_\odot~\mathrm{ yr}^{-1}$) and gas expansion velocities (2 to $21.5$ km s$^{-1}$), and include M-, S-, and C-type AGB stars. Our results generally agree with those of earlier studies, although we tend to find slightly lower mass-loss rates by about 40%, on average. We also present "bonus" lines detected during our CO observations.

preprint2015arXiv

The Spanish Square Kilometre Array White Book

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

preprint2015arXiv

Very Large Telescope observations of Gomez's Hamburger: Insights into a young protoplanet candidate

Planets are thought to form in the gas and dust disks around young stars. In particular, it has been proposed that giant planets can form through the gravitational instability of massive extended disks around intermediate-mass stars. However, we still lack direct observations to constrain this mechanism. We have spatially resolved the 8.6 and 11.2 $μ$m emission of a massive protoplanetary disk seen edge on around an A star, Gomez's Hamburger (GoHam), using VISIR at the Very Large Telescope. A compact region situated at a projected distance of $350\pm50$ AU south of the central star is found to have a reduced emission.This asymmetry is fully consistent with the presence of a cold density structure, or clump, identified in earlier CO observations, and we derive physical characteristics consistent with those observations: a mass of 0.8-11.4 Jupiter masses (for a dust-to-gas mass ratio of 0.01), a radius of about 10$^2$ astronomical units, and a local density of about $10^{7}$ cm$^{-3}$. Based on this evidence, we argue that this clump, which we call GoHam b, is a promising candidate for a young protoplanet formed by gravitational instability that might be representative of the precursors of massive planets observed around A stars, such as HR 8799 or Beta pictoris. More detailed studies at high angular resolution are needed to better constrain the physical properties of this object to confirm this proposal.

preprint2014arXiv

Future mmVLBI Research with ALMA: A European vision

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

preprint2014arXiv

SHAPEMOL: a 3-D code for calculating CO line emission in planetary and protoplanetary nebulae. Detailed model fitting of the complex nebula NGC 6302

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.

preprint2014arXiv

SiO masers from AGB stars in the vibrationally excited v=1,v=2, and v=3 states

The v=1 and v=2 J=1-0 (43 GHz), and v=1 J=2-1 (86 GHz) SiO masers are intense in AGB stars and have been mapped using VLBI showing ring-like distributions. Those of the v=1, v=2 J=1-0 masers are similar, but the spots are rarely coincident, while the v=1 J=2-1 maser arises from a well separated region farther out. These relative locations can be explained by models tools that include the overlap of two IR lines of SiO and H2O. The v=3 J=1-0 line is not directly affected by any line overlap and its spot structure and position, relative to the other lines, is a good test to the standard pumping models. We present single-dish and simultaneous VLBI observations of the v=1, v=2, and v=3 J=1-0 maser transitions of 28SiO in several AGB stars. The spatial distribution of the SiO maser emission in the v=3 J=1-0 transition from AGB stars is systematically composed of a series of spots that occupy a ring-like structure. The overall ring structure is extremely similar to that found in the other 43 GHz transitions and is very different from the structure of the v=1 J=2-1 maser. The positions of the individual spots of the different 43 GHz lines are, however, very rarely coincident, which in general is separated by about 0.3 AU (between 1 and 5 mas). These results are very difficult to reconcile with standard pumping models, which predict that the masers of rotational transitions within a given vibrational state require very similar excitation conditions, while the transitions of different vibrational states should appear in different positions. However, models including line overlap tend to predict v=1, v=2, v=3 J=1-0 population inversion to occur under very similar conditions, while the requirements for v=1 J=2-1 appear clearly different, and are compatible with the observational results.

preprint2014arXiv

The Chandra Planetary Nebula Survey (ChanPlaNS). II. X-ray Emission from Compact Planetary Nebulae

We present results from the most recent set of observations obtained as part of the Chandra X-ray observatory Planetary Nebula Survey (ChanPlaNS), the first comprehensive X-ray survey of planetary nebulae (PNe) in the solar neighborhood (i.e., within ~1.5 kpc of the Sun). The survey is designed to place constraints on the frequency of appearance and range of X-ray spectral characteristics of X-ray-emitting PN central stars and the evolutionary timescales of wind-shock-heated bubbles within PNe. ChanPlaNS began with a combined Cycle 12 and archive Chandra survey of 35 PNe. ChanPlaNS continued via a Chandra Cycle 14 Large Program which targeted all (24) remaining known compact (R_neb <~ 0.4 pc), young PNe that lie within ~1.5 kpc. Results from these Cycle 14 observations include first-time X-ray detections of hot bubbles within NGC 1501, 3918, 6153, and 6369, and point sources in HbDs 1, NGC 6337, and Sp 1. The addition of the Cycle 14 results brings the overall ChanPlaNS diffuse X-ray detection rate to ~27% and the point source detection rate to ~36%. It has become clearer that diffuse X-ray emission is associated with young (<~5x10^3 yr), and likewise compact (R_neb<~0.15 pc), PNe with closed structures and high central electron densities (n_e>~1000 cm^-3), and rarely associated with PNe that show H_2 emission and/or pronounced butterfly structures. Hb 5 is one such exception of a PN with a butterfly structure that hosts diffuse X-ray emission. Additionally, of the five new diffuse X-ray detections, two host [WR]-type CSPNe, NGC 1501 and NGC 6369, supporting the hypothesis that PNe with central stars of [WR]-type are likely to display diffuse X-ray emission.

preprint2014arXiv

The Chandra Planetary Nebulae Survey (ChanPlaNS): III. X-ray Emission from the Central Stars of Planetary Nebulae

We present X-ray spectral analysis of 20 point-like X-ray sources detected in Chandra Planetary Nebula Survey (ChanPlaNS) observations of 59 planetary nebulae (PNe) in the solar neighborhood. Most of these 20 detections are associated with luminous central stars within relatively young, compact nebulae. The vast majority of these point-like X-ray-emitting sources at PN cores display relatively "hard" ($\geq0.5$~keV) X-ray emission components that are unlikely to be due to photospheric emission from the hot central stars (CSPN). Instead, we demonstrate that these sources are well modeled by optically-thin thermal plasmas. From the plasma properties, we identify two classes of CSPN X-ray emission: (1) high-temperature plasmas with X-ray luminosities, $L_{\rm X}$, that appear uncorrelated with the CSPN bolometric luminosity, $L_{\rm bol}$; and (2) lower-temperature plasmas with $L_{\rm X}/L_{\rm bol}\sim10^{-7}$. We suggest these two classes correspond to the physical processes of magnetically active binary companions and self-shocking stellar winds, respectively. In many cases this conclusion is supported by corroborative multiwavelength evidence for the wind and binary properties of the PN central stars. By thus honing in on the origins of X-ray emission from PN central stars, we enhance the ability of CSPN X-ray sources to constrain models of PN shaping that invoke wind interactions and binarity.

preprint2014arXiv

The wind of W Hya as seen by Herschel - I. The CO envelope

Asymptotic giant branch (AGB) stars lose their envelopes by means of a stellar wind whose driving mechanism is not understood well. Characterizing the composition and thermal and dynamical structure of the outflow provides constraints that are essential for understanding AGB evolution, including the rate of mass loss and isotopic ratios. We modeled the CO emission from the wind of the low mass-loss rate oxygen-rich AGB star W Hya using data obtained by the HIFI, PACS, and SPIRE instruments onboard the Herschel Space Observatory and ground-based telescopes. $^{12}$CO and $^{13}$CO lines are used to constrain the intrinsic $^{12}$C/$^{13}$C ratio from resolved HIFI lines. The acceleration of the outflow up to about 5.5 km/s is quite slow and can be represented by a beta-type velocity law with index 5. Beyond this point, acceleration up the terminal velocity of 7 km/s is faster. Using the J=10-9, 9-8, and 6-5 transitions, we find an intrinsic $^{12}$C/$^{13}$C ratio of $18\pm10$ for W Hya, where the error bar is mostly due to uncertainties in the $^{12}$CO abundance and the stellar flux around 4.6 $μ$m. To match the low-excitation CO lines, these molecules need to be photo-dissociated at about 500 stellar radii. The radial dust emission intensity profile measured by PACS images at 70 $μ$m shows substantially stronger emission than our model predicts beyond 20 arcsec. The initial slow acceleration of the wind implies inefficient wind driving in the lower part of the envelope. The final injection of momentum in the wind might be the result of an increase in the opacity thanks to the late condensation of dust species. The derived intrinsic isotopologue ratio for W Hya is consistent with values set by the first dredge-up and suggestive of an initial mass of 2 M$_\odot$ or more. However, the uncertainty in the main-sequence mass derived based on this isotopologic ratio is large.

preprint2014arXiv

The wind of W Hya as seen by Herschel. II. The molecular envelope of W Hya

The evolution of low- and intermediate-mass stars on the asymptotic giant branch (AGB) is mainly controlled by the rate at which these stars lose mass in a stellar wind. Understanding the driving mechanism and strength of the stellar winds of AGB stars and the processes enriching their surfaces with products of nucleosynthesis are paramount to constraining AGB evolution and predicting the chemical evolution of galaxies. In a previous paper we have constrained the structure of the outflowing envelope of W Hya using spectral lines of the $^{12}$CO molecule. Here we broaden this study by modelling an extensive set of H$_{2}$O and $^{28}$SiO lines observed by the three instruments on board Herschel using a state-of-the-art molecular excitation and radiative transfer code. The oxygen isotopic ratios and the $^{28}$SiO abundance profile can be connected to the initial stellar mass and to crucial aspects of dust formation at the base of the stellar wind, respectively. The modelling of H$_{2}$O and $^{28}$SiO confirms the properties of the envelope model of W Hya derived from $^{12}$CO lines. We find an H$_2$O ortho-to-para ratio of 2.5\,$^{+2.5}_{-1.0}$, consistent with what is expected for an AGB wind. The O$^{16}$/O$^{17}$ ratio indicates that W Hya has an initial mass of about 1.5 M$_\odot$. Although the ortho- and para-H$_{2}$O lines observed by HIFI appear to trace gas of slightly different physical properties, a turbulence velocity of $0.7\pm0.1$ km s$^{-1}$ fits the HIFI lines of both spin isomers and those of $^{28}$SiO well. The ortho- and para-H$_2^{16}$O and $^{28}$SiO abundances relative to H$_{2}$ are $(6^{+3}_{-2}) \times 10^{-4}$, $(3^{+2}_{-1}) \times 10^{-4}$, and $(3.3\pm 0.8)\times 10^{-5}$, respectively. Assuming a solar silicon-to-carbon ratio, the $^{28}$SiO line emission model is consistent with about one-third of the silicon atoms being locked up in dust particles.

preprint2013arXiv

ALMA observations of the Red Rectangle, a preliminary analysis

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.

preprint2013arXiv

Extended rotating disks around post-AGB stars

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.

preprint2013arXiv

HIFISTARS Herschel/HIFI observations of VY Canis Majoris. Molecular-line inventory of the envelope around the largest known star

To gain insight into the physical conditions and kinematics of the warm (100-1000 K) gas around the red hyper-giant VY CMa, we performed sensitive high spectral resolution observations of molecular lines in the sub-mm/FIR using the HIFI instrument of the Herschel Space Observatory. We observed CO, H2O, and other molecular species, sampling excitation energies from a few tens to a few thousand K. These observations are part of the Herschel Guaranteed Time Key Program HIFISTARS. We detected the J=6-5, J=10-9, and J=16-15 lines of 12CO and 13CO at about 100, 300, and 750K above the ground state (and the 13CO J=9-8 line). These lines are crucial for improving the modelling of the internal layers of the envelope around VY CMa. We also detected 27 lines of H2O and its isotopomers, and 96 lines of species such as NH3, SiO, SO, SO2 HCN, OH and others, some of them originating from vibrationally excited levels. Three lines were not unambiguously assigned. Our observations confirm that VY CMa's envelope must consist of two or more detached components. The molecular excitation in the outer layers is significantly lower than in the inner ones, resulting in strong self-absorbed profiles in molecular lines that are optically thick in this outer envelope, for instance, low-lying lines of H2O. Except for the most abundant species, CO and H2O, most of the molecular emission detected at these sub-mm/FIR wavelengths arise from the central parts of the envelope. The spectrum of VY CMa is very prominent in vibrationally excited lines, which are caused by the strong IR pumping present in the central regions. Compared with envelopes of other massive evolved stars, VY CMa's emission is particularly strong in these vibrationally excited lines, as well as in the emission from less abundant species such as H13CN, SO, and NH3.

preprint2013arXiv

The molecular envelope of CRL 618: A new model based on Herschel/HIFI observations

We study the physical properties and molecular excitation of the different warm gas components found in the protoplanetary nebula CRL 618. We revise our previous Herschel/HIFI observations, which consist of several 12CO and 13CO lines in the far-infrared/sub-mm band. These data have been re-analyzed in detail by improving calibration, the signal-to-noise-ratio, and baseline substraction. We identify the contributions of the different nebular components to the line profiles. We have used a spatio-kinematical model to better constrain the temperature, density, and kinematics of the molecular components probed by the improved CO observations. The 12CO and 13CO J=16-15, J=10-9, and J=6-5 transitions are detected in this source. The line profiles present a composite structure showing spectacular wings in some cases, which become dominant as the energy level increases. Our analysis of the high-energy CO emission with the already known low-energy J=2-1 and J=1-0 lines confirms that the high-velocity component, or fast bipolar outflow, is hotter than previously estimated with a typical temperature of ~300 K. This component may then be an example of a very recent acceleration of the gas by shocks that has not yet cooled down. We also find that the dense central core is characterized by a very low expansion velocity, ~5 km/s, and a strong velocity gradient. We conclude that this component is very likely to be the unaltered circumstellar layers that are lost in the last AGB phase, where the ejection velocity is particularly low. The physical properties of the other two nebular components, the diffuse halo and the double empty shell, more or less agrees with the estimations derived in previous models.

preprint2013arXiv

Warm gas in the rotating disk of the Red Rectangle: accurate models of molecular line emission

We aim to study the excitation conditions of the molecular gas in the rotating disk of the Red Rectangle, the only post-Asymptotic-Giant-Branch object in which the existence of an equatorial rotating disk has been demonstrated. For this purpose, we developed a complex numerical code that accurately treats radiative transfer in 2-D, adapted to the study of molecular lines from rotating disks. We present far-infrared Herschel/HIFI observations of the 12CO and 13CO J=6-5, J=10-9, and J=16-15 transitions in the Red Rectangle. We also present our code in detail and discuss the accuracy of its predictions, from comparison with well-tested codes. Theoretical line profiles are compared with the empirical data to deduce the physical conditions in the disk by means of model fitting. We conclude that our code is very efficient and produces reliable results. The comparison of the theoretical predictions with our observations reveals that the temperature of the Red Rectangle disk is typically ~ 100-150 K, about twice as high as previously deduced from mm-wave observations of lower-J lines. We discuss the relevance of these new temperature estimates for understanding the thermodynamics and dynamics of this prototype object, as well as for interpreting observations of other rarely studied post-AGB disks. Despite our sophisticated treatment of the line formation, our model cannot explain the relatively strong line-wing emission for intermediate-J transitions. We argue that a model including a rotating disk only cannot reproduce these data and suggest that there is an additional extended (probably bipolar) structure expanding at about 7--15 km/s.

preprint2012arXiv

Herschel/HIFI observations of red supergiants and yellow hypergiants: I. Molecular inventory

Red supergiant stars (RSGs) and yellow hypergiant stars (YHGs) are believed to be the high-mass counterparts of stars in the AGB and early post-AGB phases. We study the mass-loss in the post main-sequence evolution of massive stars, through the properties of their envelopes in the intermediate and warm gas layers. These are the regions where the acceleration of the gas takes place and the most recent mass-loss episodes can be seen. We used the HIFI instrument on-board the Herschel Space Observatory to observe sub-mm and FIR transitions of CO, water, and their isotopologues in a sample of two RSGs (NML Cyg and Betelgeuse) and two YHGs (IRC+10420 and AFGL 2343) stars. We present an inventory of the detected lines and analyse the information revealed by their spectral profiles. On the basis of the results presented in an earlier study, we model the CO and 13CO emission in IRC+10420 and compare it to a set of lines ranging from the mm, to the FIR. Red supergiants have stronger high-excitation lines than the YHGs, indicating that they harbour dense and hot inner shells contributing to these transitions. Consequently, these high-J lines in RSGs originate from acceleration layers that have not yet reached the circumstellar terminal velocity and have narrower profiles than their flat-topped lower-J counterparts. The YHGs tend to lack this inner component, in line with the picture of detached, hollow envelopes derived from studies at longer wavelengths. NH3 is only detected in two sources (NML Cyg, IRC+10420), which are also observed to be the strongest water-line emitters of the studied sample. In contrast, OH is detected in all sources and does not seem to correlate with the water line intensities. We show that the IRC+10420 model derived solely from mm low-J CO transitions is capable of reproducing the high-J transitions when the temperature in the inner shell is simply lowered by about 30%.

preprint2012arXiv

Modeling the physical and excitation conditions of the molecular envelope of NGC 7027

The link between the shaping of bipolar planetary nebulae and the mass ejection activity of their central stars is still poorly understood. Appropriately characterizing the evolution of the shells ejected during the late stages of evolution is vital to gain insight into the mechanism of nebular shaping. Herschel/HIFI provides an invaluable tool by opening a new window from which to probe warm molecular gas (~50-1000 K). We present a radiative-transfer, spatio-kinematic modeling of the molecular envelope of the young planetary nebula NGC 7027 in high- and low-J 12CO and 13CO transitions observed by Herschel/HIFI and IRAM 30-m, and discuss the structure and dynamics of the molecular envelope. We have developed a code which, used along with the existing SHAPE software, implements spatio-kinematic modeling with accurate non-LTE calculations of line excitation and radiative transfer in molecular species. We have used this code to build a relatively simple "russian doll" model of the molecular envelope of NGC 7027. The model nebula consists of four nested, mildly bipolar shells plus a pair of high-velocity blobs. The innermost shell is the thinnest and shows a significant jump in physical conditions (temperature, density, abundance and velocity) with respect to the adjacent shell. This is a clear indication of a shock front in the system. Each of the high-velocity blobs is divided into two sections with different physical conditions. The presence of H2O in NGC 7027, a C-rich nebula, is likely due to photo-induced chemistry from the hot central star. The computed molecular mass of the nebula is 1.3 Msun, compatible with estimates from previous works.

preprint2012arXiv

The Chandra X-ray Survey of Planetary Nebulae (ChanPlaNS): Probing Binarity, Magnetic Fields, and Wind Collisions

We present an overview of the initial results from the Chandra Planetary Nebula Survey (ChanPlaNS), the first systematic (volume-limited) Chandra X-ray Observatory survey of planetary nebulae (PNe) in the solar neighborhood. The first phase of ChanPlaNS targeted 21 mostly high-excitation PNe within ~1.5 kpc of Earth, yielding 4 detections of diffuse X-ray emission and 9 detections of X-ray-luminous point sources at the central stars (CSPNe) of these objects. Combining these results with those obtained from Chandra archival data for all (14) other PNe within ~1.5 kpc that have been observed to date, we find an overall X-ray detection rate of ~70%. Roughly 50% of the PNe observed by Chandra harbor X-ray-luminous CSPNe, while soft, diffuse X-ray emission tracing shocks formed by energetic wind collisions is detected in ~30%; five objects display both diffuse and point-like emission components. The presence of X-ray sources appears correlated with PN density structure, in that molecule-poor, elliptical nebulae are more likely to display X-ray emission (either point-like or diffuse) than molecule-rich, bipolar or Ring-like nebulae. All but one of the X-ray point sources detected at CSPNe display X-ray spectra that are harder than expected from hot (~100 kK) central star photospheres, possibly indicating a high frequency of binary companions to CSPNe. Other potential explanations include self-shocking winds or PN mass fallback. Most PNe detected as diffuse X-ray sources are elliptical nebulae that display a nested shell/halo structure and bright ansae; the diffuse X-ray emission regions are confined within inner, sharp-rimmed shells. All sample PNe that display diffuse X-ray emission have inner shell dynamical ages <~5x10^3 yr, placing firm constraints on the timescale for strong shocks due to wind interactions in PNe.

preprint2012arXiv

Two short mass-loss events that unveil the binary heart of Minkowski's Butterfly Nebula

Studying the appearance and properties of bipolar winds is critical to understand the stellar evolution from the AGB to the planetary nebula (PN) phase. Many uncertainties exist regarding the presence and role of binary stellar systems, mainly due to the deficit of conclusive observational evidences. We investigate the extended equatorial distribution around the early bipolar planetary nebula M 2-9 ("Minkowski's Butterfly Nebula") to gather new information on the mechanism of the axial ejections. Interferometric millimeter observations of molecular emission provide the most comprehensive view of the equatorial mass distribution and kinematics in early PNe. Here we present subarcsecond angular-resolution observations of the 12CO J=2-1 line and continuum emission with the Plateau de Bure interferometer. The data reveal two ring-shaped and eccentric structures at the equatorial basis of the two coaxial optical lobes. The two rings were formed during short mass-loss episodes (~ 40 yr), separated by ~ 500 yr. Their positional and dynamical imprints provide evidence of the presence of a binary stellar system at the center, which yields critical information on its orbital characteristics, including a mass estimate for the secondary of ~< 0.2 \ms. The presence of a stellar system with a modest-mass companion at the center of such an elongated bipolar PN strongly supports the binary-based models, because these are more easily able to explain the frequent axisymmetric ejections in PNe.

preprint2011arXiv

Herschel/HIFI observations of molecular emission in protoplanetary nebulae and young planetary nebulae

We performed Herschel/HIFI observations of intermediate-excitation molecular lines in the far-infrared/submillimeter range in a sample of ten protoplanetary nebulae and young planetary nebulae. The high spectral resolution provided by HIFI yields accurate measurements of the line profiles. The observation of these high-energy transitions allows an accurate study of the excitation conditions, particularly in the warm gas, which cannot be properly studied from the low-energy lines. We have detected FIR/sub-mm lines of several molecules, in particular of 12CO, 13CO, and H2O. Emission from other species, like NH3, OH, H2^{18}O, HCN, SiO, etc, has been also detected. Wide profiles showing sometimes spectacular line wings have been found. We have mainly studied the excitation properties of the high-velocity emission, which is known to come from fast bipolar outflows. From comparison with general theoretical predictions, we find that CRL 618 shows a particularly warm fast wind, with characteristic kinetic temperature Tk >~ 200 K. In contrast, the fast winds in OH 231.8+4.2 and NGC 6302 are cold, Tk ~ 30 K. Other nebulae, like CRL 2688, show intermediate temperatures, with characteristic values around 100 K. We also discuss how the complex structure of the nebulae can affect our estimates, considering two-component models. We argue that the differences in temperature in the different nebulae can be due to cooling after the gas acceleration (that is probably due to shocks); for instance, CRL 618 is a case of very recent acceleration, less than ~ 100 yr ago, while the fast gas in OH 231.8+4.2 was accelerated ~ 1000 yr ago. We also find indications that the densest gas tends to be cooler, which may be explained by the expected increase of the radiative cooling efficiency with the density.

preprint2010arXiv

CO observations of symbiotic stellar systems

We have performed mm-wave observations with the IRAM 30m telescope of the 12CO J=2-1 and J=1-0, 13CO J=2-1 and J=1-0, and SiO J=5-4 transitions in the symbiotic stars R Aqr, CH Cyg, and HM Sge. The data were analyzed by means of a simple analytical description of the general properties of molecular emission from the inner shells around the cool star. Numerical calculations of the expected line profiles, taking into account the level population and radiative transfer under such conditions, were also performed. Weak emission of 12CO J=2-1 and J=1-0 was detected in R Aqr and CH Cyg; a good line profile of 12CO J=2-1 in R Aqr was obtained. The intensities and profile shapes of the detected lines are compatible with emission coming from a very small shell around the Mira-type star, with a radius comparable to or slightly smaller than the distance to the hot dwarf companion, 10$^{14}$ - 2 10$^{14}$ cm. We argue that other possible explanations are improbable. This region probably shows properties similar to those characteristic of the inner shells around standard AGB stars: outwards expansion at about 5 - 25 km/s, with a significant acceleration of the gas, temperatures decreasing with radius between about 1000 and 500 K, and densities ~ 10$^9$ - 3 10$^8$ cm$^{-3}$. Our model calculations are able to explain the asymmetric line shape observed in 12CO J=2-1 from R Aqr, in which the relatively weaker red part of the profile would result from selfabsorption by the outer layers (in the presence of a velocity increase and a temperature decrease with radius). The mass-loss rates are somewhat larger than in standard AGB stars, as often happens for symbiotic systems. In R Aqr, we find that the total mass of the CO emitting region is ~ 2 - 3 10$^{-5}$ Mo, corresponding to M' ~ 5 10$^{-6}$ - 10$^{-5}$ Mo/yr, and compatible with results obtained from dust emission. Taking into account other existing data on molecular emission, we suggest that the small extent of the molecule-rich gas in symbiotic systems is mainly due to molecule photodissociation by the radiation of the hot dwarf star.

preprint2010arXiv

Fast, gusty winds blowing from the core of the pre-planetary nebula M2-56

We report optical long-slit spectra and direct imaging (ground-based and with HST) of the pre-planetary nebula (pPN) M2-56 obtained at different epochs. The optical nebula is composed by shock-excited material distributed in two pairs of nested lobes with different sizes and surface brightness. The compact, bright inner lobes (ILs) have an angular size of ~1.5"x1" each and display closed, bow-shaped ends. The extended, faint outer lobes (OLs), which enclose the inner ones, have an agular size of ~13"x10". Within the ILs and the OLs the velocity increases with the distance to the center, however, the ILs show expansion velocities larger than the OLs. Consistent with the large speeds reached by the ILs (of up to ~350 km/s at the tips), we have measured the expansive proper motions of the knots (~0.03 arcsec/yr) by comparing two-epoch HST images. Moreover, we have discovered remarkable changes with time in the continuum and line emission spectrum of M2-56. In 1998, we detected a burst of Halpha emission from the nebula nucleus that is interpreted as an indication of a dense, fast (~350-500 km/s) bipolar wind from the nebula's core (referred to as "F1-wind"). Such a wind has been recently ejected (after 1989) probably as a short-duration mass-loss event. Our data also reveal an optically thick, compact structure (cocoon?) and a HII region around the central star that result from further post-AGB mass-loss after the F1-wind. Recent brightening of the scattered stellar continuum as well as an increase of scattered Halpha emission along the lobes is reported, both results pointing to a decrease of the optical depth of the circumstellar material enshrouding the star. The data presented here unveil the complex post-AGB mass-loss history of this object, whose rapid evolution is driven by multiple episodes of mass outflow, not regularly spaced in time... (abridged).

preprint2010arXiv

The widespread occurence of water vapor in the circumstellar envelopes of carbon-rich AGB stars: first results from a survey with Herschel/HIFI

We report the preliminary results of a survey for water vapor in a sample of eight C stars with large mid-IR continuum fluxes: V384 Per, CIT 6, V Hya, Y CVn, IRAS 15194-5115, V Cyg, S Cep, and IRC+40540. This survey, performed using the HIFI instrument on board the Herschel Space Observatory, entailed observations of the lowest transitions of both ortho- and para-water: the 556.936 GHz 1(10)-1(01) and 1113.343 GHz 1(11)-0(00) transitions, respectively. Water vapor was unequivocally detected in all eight of the target stars. Prior to this survey, IRC+10216 was the only carbon-rich AGB star from which thermal water emissions had been discovered, in that case with the use of the Submillimeter Wave Astronomy Satellite (SWAS). Our results indicate that IRC+10216 is not unusual, except insofar as its proximity to Earth leads to a large line flux that was detectable with SWAS. The water spectral line widths are typically similar to those of CO rotational lines, arguing against the vaporization of a Kuiper belt analog (Ford & Neufeld 2001) being the general explanation for water vapor in carbon-rich AGB stars. There is no apparent correlation between the ratio of the integrated water line fluxes to the 6.3 micron continuum flux - a ratio which measures the water outflow rate - and the total mass-loss rate for the stars in our sample.

preprint2009arXiv

The physical conditions in Gomez's Hamburger (IRAS 18059-3211), a pre-MS rotating disk

We aim to study the structure, dynamics and physical conditions of Gomez's Hamburger (IRAS 18059-3211; GoHam). We confirm that GoHam essentially consists of a flaring disk in keplerian rotation around a young, probably pre-MS star. We present high resolution SMA maps of 12CO J=2-1, 13CO J=2-1, 12CO J=3-2, and C17O J=3-2, as well as data on 12CO J=6-5 and the continuum flux at these wavelengths. Spatial resolutions up to 1" are obtained. Except for the C17O data, the dynamical ranges are larger than 10. The maps are compared to a numerical model, which simulates the emission of a rotating disk with the expected general properties of such objects; a very satisfactory fitting of our maps is obtained. The meaning and reliability of our results are thoroughly discussed. Our observations allow measurement of the main properties of GoHam at scales between ~ 1" (~ 5 10^15 cm, for the assumed distance, 300 pc) and the total extent of the nebula, 14". We are able to measure the global structure of the gas-rich disk, which is found to be flaring, and its dynamics, which is clearly dominated by keplerian rotation, with a very small degree of turbulence. The combination of different lines, particularly showing different opacities, allows us to reasonably estimate the distributions of the gas temperature and density. We clearly find a significant and sharp increase in temperature at large distances from the equator, accompanied by a decrease in density of the same order. Finally, we identify and study a condensation in the southern part of the disk that has no counterparts in the rest of the nebula. This condensation is quite extended (about 5 10^15 cm), contains a significant amount of mass (roughly, ~ 6 10^-3 Mo), and seems to be associated with a detectable distortion of the global rotation kinematics.

preprint2008arXiv

Relative astrometry of the J=1-0, v=1 and v=2 SiO masers towards R Leonis Minoris using VERA

Oxygen-rich Asymptotic Giant Branch (AGB) stars are intense emitters of SiO and H$_2$O maser lines at 43 (J=1-0, v=1 and 2) and 22 GHz, respectively. VLBI observations of the maser emission provides a unique tool to sample the innermost layers of the circumstellar envelopes in AGB stars. Nevertheless, the difficulties in achieving astrometrically aligned v=1 and v=2 SiO maser maps have traditionally prevented a unique interpretation of the observations in terms of physical underlying conditions, which depend on the nature of the SiO pumping mechanism. We have carried out observations of the SiO and H$_2$O maser emission towards RLMi, using the astrometric capabilities of VERA. Due to the too-weak emission of the reference calibrator we had to develop a special method to accurately relate the coordinates for both transitions. We present relative astrometrically aligned v=1 and v=2 J=1-0 SiO maser maps, at multiple epochs, and discuss the astrophysical results. The incorporation of astrometric information into the maps of SiO masers challenges the weak points in the current theoretical models, which will need further refinements to address the observations results.

preprint2007arXiv

Arcsecond-resolution 12CO mapping of the yellow hypergiants IRC +10420 and AFGL 2343

IRC +10420 and AFGL 2343 are the unique, known yellow hypergiants (YHGs) presenting a heavy circumstellar envelope (CSE). We aim to study the morphology, exceptional kinematics, and excitation conditions of their CSEs, and the implications for mass-loss processes. We have mapped the 12CO J=2-1 and 1-0 emission in these YHGs with the IRAM Plateau de Bure interferometer and the 30m telescope. We developed LVG models in order to analyze their circumstellar characteristics. The maps show that the overall shape of both CSEs is approximately spherical, although they also reveal several aspherical features. The CSE around IRC +10420 shows a rounded extended halo surrounding a bright inner region, with both components presenting aspherical characteristics. It presents a brightness minimum at the center. The envelope around AFGL 2343 is a detached shell, showing spherical symmetry and clumpiness at a level of about 15% of the maximum brightness. The envelopes expand isotropically at about 35 km/s, about two or three times faster than typical CSEs around AGB stars. High temperatures (~ 200 K) are derived for the innermost regions in IRC +10420, while denser and cooler (~ 30 K) gas is found in AFGL 2343. The mass-loss processes in these YHGs have been found to be similar. The deduced mass-loss rates (~ 10E-4 - 10E-3 Msun/yr) are much higher than those obtained in AGB stars, and they present significant variations on time scales of ~ 1000 yr.