Strategic Plan for Astronomy in the Netherlands 2011 - 2020, written by the Netherlands Committee for Astronomy (NCA), on behalf of the excellence research school in astronomy NOVA, (combining the university astronomy institutes of the universities of Amsterdam, Groningen, Leiden and Nijmegen), the NWO division of Physical Sciences, the Netherlands Institute for Radio Astronomy ASTRON and the Netherlands Institute for Space Research SRON. The Strategic plan outlines the scientific priorities for Dutch astronomy in the next decade; the instrumentation effort required to address these priorities, and the connection between astronomical instrumentation and technology development and fundamental technological R&D; the financial contours needed to realise the priorities; and the role of Dutch astronomy in education and outreach. The Strategic Plan also includes a retrospective on the achievements since the last Strategic Plan (2000) and a forward look beyond 2020.

We present deep images of dust continuum emission at 450, 800, and 850 micron of the dark cloud LDN 1689N which harbors the low-mass young stellar objects (YSOs) IRAS 16293-2422A and B (I16293A and I16293B) and the cold prestellar object I16293E. Toward the positions of I16293A and E we also obtained spectra of CO-isotopomers and deep submillimeter observations of chemically related molecules with high critical densities. To I16293A we report the detection of the HDO 1_01 - 0_00 and H2O 1_10 - 1_01 ground-state transitions as broad self-reversed emission profiles with narrow absorption, and a tentative detection of H2D+ 1_10 - 1_11. To I16293E we detect weak emission of subthermally excited HDO 1_01 - 0_00. Based on this set of submillimeter continuum and line data we model the envelopes around I16293A and E. The density and velocity structure of I16293A is fit by an inside-out collapse model, yielding a sound speed of a=0.7 km/s, an age of t=(0.6--2.5)e4 yr, and a mass of 6.1 Msun. The density in the envelope of I16293E is fit by a radial power law with index -1.0+/-0.2, a mass of 4.4 Msun, and a constant temperature of 16K. These respective models are used to study the chemistry of the envelopes of these pre- and protostellar objects. The [HDO]/[H2O] abundance ratio in the warm inner envelope of I16293A of a few times 1e-4 is comparable to that measured in comets. This supports the idea that the [HDO]/[H2O] ratio is determined in the cold prestellar core phase and conserved throughout the formation process of low-mass stars and planets.

We present submillimeter observations of rotational transitions of carbon monoxide from J = 2 -> 1 up to 7 -> 6 for a sample of Asymptotic Giant Branch stars and red supergiants. It is the first time that the high transitions J = 6 -> 5 and 7 -> 6 are included in such a study. With line radiative transfer calculations, we aim to determine the mass-loss history of these stars by fitting the CO line intensities. We find that the observed line intensities of the high transitions, including the J = 4 -> 3 transition, are significantly lower than the predicted values. We conclude that the physical structure of the outflow of Asymptotic Giant Branch stars is more complex than previously thought. In order to understand the observed line intensities and profiles, a physical structure with a variable mass-loss rate and/or a gradient in stochastic gas velocity is required. A case study of the AGB star WX Psc is performed. We find that the CO line strengths may be explained by variations in mass-loss on time scales similar to those observed in the separated arc-like structures observed around post-AGB stars. In addition, a gradient in the stochastic velocity may play a role. Until this has been sorted out fully, any mass loss determinations based upon single CO lines will remain suspect.

The nature of the tenuous disks around A stars has raised considerable controversy in the literature during the past decade. The debate whether or not the disk around beta Pictoris contains gaseous molecular hydrogen is only the most recent example. Since CO is in general a poor tracer for the gas content of these low mass disks, we discuss here detailed emission line calculations for alternative tracers like C and C+, based on recent optically thin disk models by Kamp & van Zadelhoff (2001). The [CII] 157.7 mum line was searched toward Vega and beta Pictoris -- the most prominent A stars with disks -- using ISO LWS data, and a tentative detection is reported toward the latter object. From a comparison with emission line observations as well as absorption line studies of both stars, the gas-to-dust ratio is constrained to lie between 0.5 and 9 for beta Pictoris. For Vega the [CII] observations indicate an upper limit of 0.2 M_Earth for the disk gas mass. Predicted line intensities of C+ and C are presented for a range of models and appear promising species to trace the gas content in the disks around A stars with future instrumental capabilities (SOFIA, Herschel, APEX and ALMA). Searches for CO emission should focus on the J=3-2 line.

We present observations of the v2=0 and vibrationally excited v2=1 J=9-8 rotational lines of HCN at 797 GHz toward the deeply embedded massive young stellar object GL 2591, which provide the missing link between the extended envelope traced by lower-J line emission and the small region of hot (T_ex >= 300 K), abundant HCN seen in 14 micron absorption with the Infrared Space Observatory (ISO). The line ratio yields T_ex=720^+135_-100 K and the line profiles reveal that the hot gas seen with ISO is at the velocity of the protostar, arguing against a location in the outflow or in shocks. Radiative transfer calculations using a depth-dependent density and temperature structure show that the data rule out a constant abundance throughout the envelope, but that a model with a jump of the abundance in the inner part by two orders of magnitude matches the observations. Such a jump is consistent with the sharp increase in HCN abundance at temperatures >~230 K predicted by recent chemical models in which atomic oxygen is driven into water at these temperatures. Together with the evidence for ice evaporation in this source, this result suggests that we may be witnessing the birth of a hot core. Thus, GL 2591 may represent a rare class of objects at an evolutionary stage just preceding the `hot core' stage of massive star formation.

We report the detection of the 1_{10}-1_{11} ground state transition of ortho-H_2D^+ at 372.421 GHz in emission from the young stellar object NGC 1333 IRAS 4A. Detailed excitation models with a power-law temperature and density structure yield a beam-averaged H_2D^+ abundance of 3 x 10^{-12} with an uncertainty of a factor of two. The line was not detected toward W 33A, GL 2591, and NGC 2264 IRS, in the latter source at a level which is 3-8 times lower than previous observations. The H_2D^+ data provide direct evidence in support of low-temperature chemical models in which H_2D^+ is enhanced by the reaction of H_3^+ and HD. The H_2D^+ enhancement toward NGC 1333 IRAS 4A is also reflected in the high DCO^+/HCO^+ abundance ratio. Simultaneous observations of the N_2H^+ 4-3 line show that its abundance is about 50-100 times lower in NGC 1333 IRAS 4A than in the other sources, suggesting significant depletion of N_2. The N_2H^+ data provide independent lower limits on the H_3^+ abundance which are consistent with the abundances derived from H_2D^+. The corresponding limits on the H_3^+$ column density agree with recent near-infrared absorption measurements of H_3^+ toward W 33A and GL 2591.