Issue |
A&A
Volume 518, July-August 2010
Herschel: the first science highlights
|
|
---|---|---|
Article Number | L145 | |
Number of page(s) | 5 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/201014641 | |
Published online | 16 July 2010 |
Herschel: the first science highlights
LETTER TO THE EDITOR
PACS and SPIRE spectroscopy of the red supergiant VY CMa![[*]](/icons/foot_motif.png)
P. Royer1 - L. Decin1,12 - R. Wesson2 - M. J. Barlow2 - E. T. Polehampton7,9 - M. Matsuura2,15 - M. Agúndez14,16 - J. A. D. L. Blommaert1 - J. Cernicharo14 - M. Cohen4 - F. Daniel14 - P. Degroote1 - W. De Meester1 - K. Exter1 - H. Feuchtgruber3 - W. K. Gear5 - H. L. Gomez5 - M. A. T. Groenewegen11 - P. C. Hargrave5 - R. Huygen1 - P. Imhof10 - R. J. Ivison6 - C. Jean1 - F. Kerschbaum13 - S. J. Leeks7 - T. Lim7 - R. Lombaert1 - G. Olofsson8 - T. Posch13 - S. Regibo1 - G. Savini2 - B. Sibthorpe6 - B. M. Swinyard7 - B. Vandenbussche1 - C. Waelkens1 - D. K. Witherick2 - J. A. Yates2
1 - Instituut voor Sterrenkunde, Katholieke Universiteit Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
2 - Dept of Physics & Astronomy, University College London, Gower St, London WC1E 6BT, UK
3 - Max-Planck-Institut für extraterrestrische Physik,
Giessenbachstrasse, 85748, Germany
4 - Radio Astronomy Laboratory, University of California at Berkeley, CA 94720, USA
5 - School of Physics and Astronomy, Cardiff University, 5 The Parade, Cardiff, Wales CF24 3YB, UK
6 - UK Astronomy Technology Centre, Royal Observatory Edinburgh, Blackford Hill, Edinburgh EH9 3HJ, UK
7 - Space Science and Technology Department, Rutherford Appleton Laboratory, Oxfordshire, OX11 0QX, UK
8 - Dept of Astronomy, Stockholm University, AlbaNova University Center, Roslagstullsbacken 21, 10691 Stockholm, Sweden
9 - Department of Physics, University of Lethbridge, Lethbridge, Alberta, T1J 1B1, Canada
10 - Blue Sky Spectroscopy, 9/740 4 Ave S, Lethbridge, Alberta T1J 0N9, Canada
11 - Royal Observatory of Belgium, Ringlaan 3, 1180 Brussels, Belgium
12 - Sterrenkundig Instituut Anton Pannekoek, University of Amsterdam, Science Park 904, 1098 Amsterdam, The Netherlands
13 - University of Vienna, Department of Astronomy,
Türkenschanzstraße 17, 1180 Vienna, Austria
14 - Astrophysics Dept, CAB (INTA-CSIC), Crta Ajalvir km4, 28805 Torrejon de Ardoz, Madrid, Spain
15 - Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking, Surrey RH5 6NT, UK
16 - LUTH, Observatoire de Paris-Meudon, LERMA UMR CNRS 8112, 5 place Jules Janssen, 92195 Meudon Cedex, France
Received 31 March 2010 / Accepted 20 April 2010
Abstract
With a luminosity >105
and a mass-loss rate of
yr-1,
the red supergiant VY CMa truly is a spectacular object. Because
of its extreme evolutionary state, it could explode as supernova any
time. Studying its circumstellar material, into which the supernova
blast will run, provides interesting constraints on supernova
explosions and on the rich chemistry taking place in such complex
circumstellar envelopes. We have obtained spectroscopy of VY CMa
over the full wavelength range offered by the PACS and SPIRE
instruments of Herschel, i.e. 55-672 micron. The
observations show the spectral fingerprints of more than
900 spectral lines, of which more than half belong to water. In
total, we have identified 13 different molecules and some of their
isotopologues. A first analysis shows that water is abundantly present,
with an ortho-to-para ratio as low as
1.3:1, and that chemical non-equilibrium processes determine the abundance fractions in the inner envelope.
Key words: line: identification - stars: chemically peculiar - circumstellar matter - stars: mass-loss - supergiants
1 Introduction
The oxygen-rich red supergiant VY CMa is one of the most extreme objects of its kind. Its large mass-loss rate results in a very dense circumstellar envelope and an enormous infrared excess (Harwit et al. 2001). Infrared observations are invaluable for in-depth studies of the complex chemical and dynamical processes at work in its circumstellar envelope. By this means, VY CMa has been shown to harbour a rich, water-dominated chemistry (e.g. Neufeld et al. 1999; Polehampton et al. 2010). The unprecedented wavelength coverage and spectral resolution of the Herschel instruments will undoubtedly push our understanding of the mass-loss related phenomena a step further.
Table 1: Journal of observations.
In this context, VY CMa has been observed as part of the Herschel guaranteed time key programme ``Mass-loss of Evolved StarS'', MESS (Groenewegen et al., in prep.). The spectra discussed hereafter have been obtained with the PACS and SPIRE instruments of Herschel (Pilbratt et al. 2010). The instruments and their in-orbit performance are described in Poglitsch et al. (2010) and Griffin et al. (2010). The SPIRE calibration methods and accuracy are outlined by Swinyard et al. (2010).
![]() |
Figure 1: The complete PACS + SPIRE spectrum of VY CMa. The various spectral bands are presented in order of increasing wavelength. For the sake of clarity, PACS band B2A is not included, as it matches the wavelength range of PACS-B3A. The ISO-LWS spectrum is overplotted for comparison (Polehampton et al. 2010). For PACS, the spectrum of the central spaxel is displayed, after application of the point source correction (version 2.0). The strongest spectral features of the main molecules are indicated. |
Open with DEXTER |
2 Observations and data reduction
The PACS observations of VY CMa consist of 5 spectral energy distribution (SED) scans, acquired on 3 and 13 November 2009 (see Table 1 for details). All observations were performed with a non-standard version of the PACS-SED observing mode. The difference with the standard mode resides in the wavelength coverage, as our SEDs cover the complete PACS wavelength range in the four bands, in two observations. All observations were ``chop-nodded'', and contained one single nodding cycle, and one single up-down scan in wavelength.
The data were reduced with the nominal pipeline, except that we combined the spectra obtained at the two nod positions after rebinning only. The flux calibration is based on a constant detector response and the wavelength ranges affected by light leaks have been excluded from the present analysis.
The original design of the PACS SED observing mode included
Nyquist spectral sampling. However, in the early version that was used
for these observations, a few wavelength ranges were still slightly
undersampled. Consequently, rebinning the spectra at the nominal
resolution of the spectrometer returns a number of empty bins. For the
modeling, we worked at the nominal resolution, and used linear
interpolation at the wavelengths of the empty bins. For the line
identifications, we have rebinned the spectrum to a slightly degraded
resolution (half Nyquist).
The SPIRE spectrum of VY CMa was obtained on operational day (OD) 123 and consisted of 30 repetitions. Each repetition consisted of 1 forward and 1 reverse scan of the FTS. The total on-source integration time was 3996 s. The unapodized spectral resolution was 1.4 GHz (0.048 cm-1), after apodization (using extended Norton-Beer function 1.5; Naylor & Tahic 2007) this is 2.1 GHz (0.07 cm-1).
![]() |
Figure 2:
Line peak intensity, expressed in Jy, for the o-H2O (71,6-62,5) line at 66.09 |
Open with DEXTER |
3 Results
3.1 Observational results
Figure 1 shows the full PACS and SPIRE spectrum of VY CMa. More than 930 molecular emission lines have been detected, 95% of which we could already identify. More than 400 of these lines originate from water vapour, but many different molecules contribute to the spectra: the identified species and their isotopologues are 12CO, 13CO, C17O, C18O, ortho-H216O (o-H216O), para-H216O (p-H216O), o-H217O, o-H218O, p-H218O, CN, CS, SO, NH3, OH, 28SiO, 29SiO, 30SiO, SiS, HCN, SO2, H2S, [C II]. H3O+ is possibly present as well, as reported by Phillips et al. (1992) and Polehampton et al. (2010), but due to the high density of lines, detailed modeling is needed to give a definitive answer (Matsuura et al. 2010, in prep.). Both low-excitation (e.g., CO J=4-3, at 55 K) and high-excitation (e.g., o-H2O (168,9-159,6) at 5500 K) are found, allowing us to study the molecular envelope from the outermost layers down to the dust formation region close to the star. In this way, the thermodynamical structure of the complex envelope and the chemical processes in the envelope can be determined.
The flux calibrations of the two instruments match within 15%,
the different spectral ranges of a given instrument agree within 5%.
The fluxes measured by PACS and IRAS at 60 and 100 m are similar within 5 and 15%. The agreement between SPIRE and ISO-LWS is within a few percent at 195
m, and the PACS and ISO-LWS spectra coincide to
% over all wavelengths. Note that the ISO-LWS fluxes were found to be high in comparison with IRAS. Gry et al. (2003)
reported that the LWS fluxes were 20% above those of IRAS at the flux
level of VY CMa. On the other hand, the PACS flux in Fig. 1 might be slightly underestimated, due to the extension of the source.
Indeed, although it appears unresolved at first sight, in lines as well as in continuum (Fig. 2), VY CMa can not be considered as a point source for PACS. Its flux, integrated over the field, shows an excess of % in the blue and
%
in the red compared to the flux measured in the central spatial pixel
(spaxel) after application of the point source correction, for the
observations of OD 173 (shown in Fig. 1). The same comparison, carried out on the Neptune SED that was obtained on the same OD shows Neptune (
2
diameter during the observation) to be compatible with a point source,
hence confirming the extension of the emission around VY CMa. The
ratio between central spaxel and integrated flux was computed for all
raster positions. All of these ratios show comparable or larger
excesses, hence ruling out the possibility that, on OD 173, the
satellite was pointing significantly out of the peak of emission.
Additionally, these central to integrated flux ratios, as well
as dedicated measurements in the spectra, reveal that the line to
continuum ratios are constant to within 5% over the whole field of
view, at all PACS wavelengths, the only exception being the [C II] line at 157.74 m,
which has a stronger integrated flux due to its interstellar origin.
Finally, we created maps of the line ratios between high and
low-excitation transitions of water, e.g. o-H2O (71,6-62,5) 66.09
m / o-H2O (33,0-22,1) 66.43
m in the blue and o-H2O (73,4-72,5) 166.82
m / o-H2O (21,2-10,1) 179.53
m
in the red part of the PACS SED.
These maps show no structure either, supporting the interpretation that
all the emission is localized in the central part of our field. Hence,
although detailed comparison with a high quality PSF is still to be
performed, all PACS observations are consistent with an extension
of 5 to 10
,
which is compatible with the central part of the nebula as imaged in the visible, near and mid-infrared by Smith et al. (2001).
3.2 Modeling the envelope: CO, H
O, HCN and SiO
![]() |
Figure 3:
The continuum-subtracted PACS spectrum (black) between 173 and 188.5 |
Open with DEXTER |
Due to the rich water spectrum, most of the emission lines are blended. In addition, many line transitions have high optical depth. Consequently, the fractional molecular abundance can only be retrieved from a proper solution of the radiative transfer equation coupled to the statistical equilibrium equations. Although we realize that the envelope around VY CMa has a very complex structure (Smith et al. 2009), we have assumed a spherically symmetric envelope. Using this assumption, the average density over all directions of 12CO, 13CO, HCN, o-H216O, p-H216O, and 28SiO can be determined. A more detailed study will be presented in Matsuura et al. (in prep.).
We have used the non-local thermodynamic equilibrium (non-LTE) radiative transfer code GASTRoNOoM (Decin et al. 2010,2006), with line lists and collisional rates as specified by Decin et al. (2010).
The kinetic temperature is calculated from the balance of cooling and
heating, with the main contributions coming from adiabatic expansion,
grain-gas collisions, and CO and H2O ro-vibrational
transitions. The conservation of momentum delivers the velocity
structure. The model parameters are summarized in Table 2, a comparison between observed data and theoretical predictions is shown in Fig. 3. Specifically, we derive a mean mass-loss rate of
yr-1, a fractional abundance of SiO of
,
of HCN we find
and of o-H2O being
.
We derive a 12CO/13CO
isotopic ratio of 60, and an ortho-to-para water ratio (OPR) of 1.27:1.
Due to the high water optical depth, the latter remains rather
uncertain.
Table 2: Model parameters for VY CMa.
4 Discussion
HCN: Ziurys et al. (2009) derived a HCN fractional abundance [HCN/H2] of
in the spherical outflow and of
in the so-called red and blue flow for a mass-loss rate of
yr-1
and a distance of 1500 pc, being similar to our results. This
abundance is orders of magnitude higher than the thermodynamic
equilibrium value (TE) [HCN/H
]
of
(Duari et al. 1999). A similar high abundance of HCN has also been derived for oxygen-rich (O-rich) asymptotic giant branch (AGB) stars (Ziurys et al. 2009; Decin et al. 2010). Photochemical predictions for the outer envelope of O-rich AGB stars (Willacy & Millar 1997)
are clearly too low by at least an order of magnitude. Pulsation driven
shock models can account for this high concentration of HCN, yielding a
HCN fractional abundance in the inner wind (between 2.2 and 5
)
ranging from
to
(Duari et al. 1999; Cherchneff 2006).
SiO: the derived SiO abundance is consistent with the results of Ziurys et al. (2007). Both inner wind shock-induced non-equilibrium chemistry (Duari et al. 1999; Cherchneff 2006) and outer wind photochemical processes (Willacy & Millar 1997) predict almost the same abundance value.
H2O: Zubko et al. (2004) derived a water abundance [H2O/H2] of
from fitting ISO SWS and LWS observations, with an OPR of 1:1,
similar to our results. Low OPR ratios are also derived for the
hypergiant NML Cyg by Zubko et al. (2004) and for the O-rich AGB star W Hya (1:1.3 by Zubko & Elitzur (2000) and 1:1 by Barlow et al. (1996)). An OPR close to 1:1 reflects a cold spin temperature, around 15 K (Mumma et al. 1987).
It is thought that the OPR value should correspond to the physical
conditions where water was formed, i.e. the warm and dense regions of
the envelope where the chemistry is under thermodynamical equilibrium
and controlled by three body reactions. Hence, we could expect an OPR
value of 3. The large opacity of the water lines very sensibly reduces
our sensitivity to the OPR value, as also noted by Barlow et al. (1996), so that, even using H218O, we can at present not rule out the aforementioned expected value.
12CO/13CO: from low-excitation CO lines, Muller et al. (2007) derived also a 12CO/13CO ratio of 60, while Milam et al. (2009) obtained a 12C/13C value of 46 for the spherical wind. This indicates that the isotopologue ratio is (roughly) constant in the spherical wind.
5 Conclusion
Through the detection of more than 900 emission lines of various chemical species in VY CMa, PACS and SPIRE have shown their diagnostic strength to unravel the chemically complex mass-loss processes at work in evolved stars. Specifically, we have derived the HCN, SiO and water fractional abundances of the circumstellar material around VY CMa, which also shows a suprisingly low ortho-to-para water ratio, close to 1. The origin of the latter is not yet understood, but the Herschel observations of many evolved stars in the MESS programme might shed light on this topic.
AcknowledgementsPACS has been developed by a consortium of institutes led by MPE (Germany) and including UVIE (Austria); KUL, CSL, IMEC (Belgium); CEA, OAMP (France); MPIA (Germany); IFSI, OAP/AOT, OAA/CAISMI, LENS, SISSA (Italy); IAC (Spain). This development has been supported by the funding agencies BMVIT (Austria), ESA-PRODEX (Belgium), CEA/CNES (France), DLR (Germany), ASI (Italy), and CICT/MCT (Spain). SPIRE has been developed by a consortium of institutes led by Cardiff Univ. (UK) and including Univ. Lethbridge (Canada); NAOC (China); CEA, LAM (France); IFSI, Univ. Padua (Italy); IAC (Spain); Stockholm Observatory (Sweden); Imperial College London, RAL, UCL-MSSL, UKATC, Univ. Sussex (UK); Caltech, JPL, NHSC, Univ. Colorado (USA). This development has been supported by national funding agencies: CSA (Canada); NAOC (China); CEA, CNES, CNRS (France); ASI (Italy); MCINN (Spain); SNSB (Sweden); STFC (UK); and NASA (USA). MG, JB, WDM, KE, RH, CJ, SR, PR & BV acknowledge support from the Belgian Federal Science Policy Office via the PRODEX Programme of ESA. L.D. acknowledges support from the Fund for Scientific Research - Flanders (FWO). F.K. acknowledges funding by the Austrian Science Fund FWF under project numbers P18939-N16 & I163-N16. The authors thank D. Ladjal, & B. Acke for their support.
References
- Barlow, M. J., Nguyen-Q-Rieu, Truong-Bach, et al. 1996, A&A, 315, L241 Cherchneff, I. 2006, A&A, 456, 1001 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Choi, Y. K., Hirota, T., Honma, M., et al. 2008, PASJ, 60, 1007 [NASA ADS] [Google Scholar]
- Decin, L., Hony, S., de Koter, A., et al. 2006, A&A, 456, 549 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Decin, L., De Beck, E., Brünken, S., et al. 2010, A&A, 516, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Duari, D., Cherchneff, I., & Willacy, K. 1999, A&A, 341, L47 [NASA ADS] [Google Scholar]
- Griffin, M. J., et al. 2010, A&A, 518, L3 [Google Scholar]
- Gry, C., Swinyard, B., Harwood, A., et al. 2003, The ISO Handbook, Vol. III - LWS - The Long Wavelength Spectrometer, ed. Müller, T., et al. [Google Scholar]
- Harwit, M., Malfait, K., Decin, L., et al. 2001, ApJ, 557, 844 [NASA ADS] [CrossRef] [Google Scholar]
- Milam, S. N., Woolf, N. J., & Ziurys, L. M. 2009, ApJ, 690, 837 [NASA ADS] [CrossRef] [Google Scholar]
- Muller, S., Dinh-V-Trung, Lim, J., et al. 2007, ApJ, 656, 1109 [NASA ADS] [CrossRef] [Google Scholar]
- Mumma, M. J., Weaver, H. A., & Larson, H. P. 1987, A&A, 187, 419 [NASA ADS] [Google Scholar]
- Naylor, D. A., & Tahic, M. K. 2007, J. Opt. Soc. Amer. A, 24, 3644 [Google Scholar]
- Neufeld, D. A., Feuchtgruber, H., Harwit, M., et al. 1999, ApJ, 517, L147 [Google Scholar]
- Phillips, T. G., van Dishoeck, E. F., & Keene, J. 1992, ApJ, 399, 533 [NASA ADS] [CrossRef] [Google Scholar]
- Pilbratt, G. L., et al. 2010, A&A, 518, L1 [CrossRef] [EDP Sciences] [Google Scholar]
- Poglitsch, A., et al. 2010, A&A, 518, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Polehampton, E. T., Menten, K. M., van der Tak, F. F. S., & White, G. J. 2010, A&A, 510, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Smith, N., Hinkle, K. H., & Ryde, N. 2009, AJ, 137, 3558 [Google Scholar]
- Smith, N., Humphreys, R. M., Davidson, K., et al. 2001, AJ, 121, 1111 [NASA ADS] [CrossRef] [Google Scholar]
- Swinyard, B. M., et al. 2010, A&A, 518, L4 [Google Scholar]
- Willacy, K., & Millar, T. J. 1997, A&A, 324, 237 [NASA ADS] [Google Scholar]
- Ziurys, L. M., Milam, S. N., Apponi, A. J., & Woolf, N. J. 2007, Nature, 447, 1094 [Google Scholar]
- Ziurys, L. M., Tenenbaum, E. D., Pulliam, R. L., Woolf, N. J., & Milam, S. N. 2009, ApJ, 695, 1604 [NASA ADS] [CrossRef] [Google Scholar]
- Zubko, V., & Elitzur, M. 2000, ApJ, 544, L137 [NASA ADS] [CrossRef] [Google Scholar]
- Zubko, V., Li, D., Lim, T., Feuchtgruber, H., & Harwit, M. 2004, ApJ, 610, 427 [NASA ADS] [CrossRef] [Google Scholar]
Footnotes
- ... CMa
- Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.
- ... Nyquist)
- The standard PACS-SED was modified on OD 305, and the full Nyquist sampling has been implemented since then.
All Tables
Table 1: Journal of observations.
Table 2: Model parameters for VY CMa.
All Figures
![]() |
Figure 1: The complete PACS + SPIRE spectrum of VY CMa. The various spectral bands are presented in order of increasing wavelength. For the sake of clarity, PACS band B2A is not included, as it matches the wavelength range of PACS-B3A. The ISO-LWS spectrum is overplotted for comparison (Polehampton et al. 2010). For PACS, the spectrum of the central spaxel is displayed, after application of the point source correction (version 2.0). The strongest spectral features of the main molecules are indicated. |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Line peak intensity, expressed in Jy, for the o-H2O (71,6-62,5) line at 66.09 |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
The continuum-subtracted PACS spectrum (black) between 173 and 188.5 |
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.