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%.

Circumstellar envelopes (CSEs) of a variety of evolved stars have been found to contain ammonia (NH3) in amounts that exceed predictions from conventional chemical models by many orders of magnitude. The observations reported here were performed in order to better constrain the NH3 abundance in the CSEs of four, quite diverse, oxygen-rich stars using the NH3 ortho J_K = 1_0 - 0_0 ground-state line. We used the Heterodyne Instrument for the Far Infrared aboard Herschel to observe the NH3 J_K = 1_0 - 0_0 transition near 572.5 GHz, simultaneously with the ortho-H2O J_Ka,Kc = 1_1,0 -1_0,1 transition, toward VY CMa, OH 26.5+0.6, IRC+10420, and IK Tau. We conducted non-LTE radiative transfer modeling with the goal to derive the NH3 abundance in these objects' CSEs. For the latter two stars, Very Large Array imaging of NH3 radio-wavelength inversion lines were used to provide further constraints, particularly on the spatial extent of the NH3-emitting regions. Results. We find remarkably strong NH3 emission in all of our objects with the NH3 line intensities rivaling those obtained for the ground state H2O line. The NH3 abundances relative to H2 are very high and range from 2 x 10-7 to 3 x 10-6 for the objects we have studied. Our observations confirm and even deepen the circumstellar NH3 enigma. While our radiative transfer modeling does not yield satisfactory fits to the observed line profiles, it leads to abundance estimates that confirm the very high values found in earlier studies. New ways to tackle this mystery will include further Herschel observations of more NH3 lines and imaging with the Expanded Very Large Array.

Aims: By performing a detailed radiative transfer analysis, we determine fractional abundances of circumstellar H2O in the envelopes around six M-type asymptotic giant branch stars. The models are also used to predict H2O spectral line emission for the upcoming Herschel/HIFI mission. Methods: We use Infrared space observatory long wavelength spectrometer spectra to constrain the circumstellar fractional abundance distribution of ortho-H2O, using a non-local thermal equilibrium, and non-local, radiative transfer code based on the accelerated lambda iteration formalism. The mass-loss rates and kinetic temperature structures for the sample stars are determined through radiative transfer modelling of CO line emission based on the Monte-Carlo method. The density and temperature profiles of the circumstellar dust grains are determined through spectral energy distribution modelling using the publicly available code Dusty. Results: The determined ortho-H2O abundances lie between 1e-4 and 1.5e-3 relative to H2, with the exception of WX Psc, which has a much lower estimated ortho-H2O abundance of only 2e-6, possibly indicating H_2O adsorption onto dust grains or recent mass-loss-rate modulations. The estimated abundances are uncertain by, at best, a factor of a few. Conclusions: The high water abundance found for the majority of the sources suggests that either the `normal' chemical processes are very effective in producing H2O, or else non-local thermal equilibrium atmospheric chemistry, grain surface reactions, or a release of H_2O (e.g. from icy bodies like Kuiper belt objects) play a role. We provide predictions for ortho-H2O lines in the spectral window of Herschel/HIFI.

The main aim is to derive reliable mass-loss rates and circumstellar SiO abundances for a sample of 40 S-type AGB stars based on new multi-transitional CO and SiO radio line observations. In addition, the results are compared to previous results for M-type AGB stars and carbon stars to look for trends with chemical type. The circumstellar envelopes are assumed to be spherically symmetric and formed by a constant mass-loss rate. The mass-loss rates are estimated from fitting the CO observations using a non-local, non-LTE radiative transfer code. Once the physical properties of the circumstellar envelopes are determined, the same radiative transfer code is used to model the observed SiO lines in order to derive circumstellar abundances and the sizes of the SiO line-emitting regions. We have estimated mass-loss rates of 40 S-type AGB stars and find that the derived mass-loss rates have a distribution that resembles those previously derived for similar samples of M-type AGB stars and carbon stars. The estimated mass-loss rates also correlate well with the corresponding expansion velocity. In all, this indicates that the mass loss is driven by the same mechanism in all three chemical types of AGB stars. In addition, we have estimated the circumstellar fractional abundance of SiO relative to H2 in 26 of the sample S-type AGB stars. The derived SiO abundances are, on average, about an order of magnitude higher than predicted by stellar atmosphere thermal equilibrium chemistry, indicating that non-equilibrium chemical processes determines the abundance of SiO in the circumstellar envelope. Moreover, a comparison with the results for M-type AGB stars and carbon stars show that for a certain mass-loss rate, the circumstellar SiO abundance seems independent (although with a large scatter) of the C/O-ratio.

Aims: Spectrally resolved circumstellar H2O(1_10 - 1_01) lines have been obtained towards three M-type AGB stars using the Odin satellite. This provides additional strong constrains on the properties of circumstellar H2O and the circumstellar envelope. Methods: ISO and Odin satellite H2O line data are used as constraints for radiative transfer models. Special consideration is taken to the spectrally resolved Odin line profiles, and the effect of excitation to the first excited vibrational states of the stretching modes (nu1=1 and nu3=1) on the derived abundances is estimated. A non-local, radiative transfer code based on the ALI formalism is used. Results: The H2O abundance estimates are in agreement with previous estimates. The inclusion of the Odin data sets stronger constraints on the size of the H2O envelope. The H2O(1_10 - 1_01) line profiles require a significant reduction in expansion velocity compared to the terminal gas expansion velocity determined in models of CO radio line emission, indicating that the H2O emission lines probe a region where the wind is still being accelerated. Including the nu3=1 state significantly lowers the estimated abundances for the low-mass-loss-rate objects. This shows the importance of detailed modelling, in particular the details of the infrared spectrum in the range 3 to 6 micron, to estimate accurate circumstellar H2O abundances. Conclusions: Spectrally resolved circumstellar H2O emission lines are important probes of the physics and chemistry in the inner regions of circumstellar envelopes around asymptotic giant branch stars. Predictions for H2O emission lines in the spectral range of the upcoming Herschel/HIFI mission indicate that these observations will be very important in this context.

In the recent literature there has been some doubt as to the reliability of CO multi-transitional line observations as a mass-loss-rate estimator for AGB stars. Mass-loss rates for 10 intermediate- to high-mass-loss-rate AGB stars are derived using a detailed non-LTE, non-local radiative transfer code based on the Monte-Carlo method to model the CO radio line intensities. The circumstellar envelopes are assumed to be spherically symmetric and formed by constant mass-loss rates. The energy balance is solved self-consistently and the effects of dust on the radiation field and thermal balance are included. An independent estimate of the mass-loss rate is also obtained from the combination of dust radiative transfer modelling with a dynamical model of the gas and dust particles. We find that the CO radio line intensities and shapes are successfully reproduced for the majority of our objects assuming a constant mass-loss rate. Moreover, the CO line intensities are only weakly dependent on the adopted micro-turbulent velocity, in contrast to recent claims in the literature. The two methods used in the present work to derive mass-loss-rates are consistent within a factor of ~3 for intermediate- to high-mass-loss-rate objects, indicating that this is a lower limit to the uncertainty in present mass-loss-rate estimates. We find a tentative trend with chemistry. Mass-loss rates from the dust/dynamical model are systematically higher than those from the CO model for the carbon stars and vice versa for the M-type stars. This could be ascribed to a discrepancy in the adopted CO/H_2-abundance ratio, but we caution that the sample is small and systematic errors cannot be excluded.

H13CN J=8-7 sub-millimetre line emission produced in the circumstellar envelope around the extreme carbon star IRC+10216 has been imaged at sub-arcsecond angular resolution using the SMA. Supplemented by a detailed excitation analysis the average fractional abundance of H13CN in the inner wind (< 5E15 cm) is estimated to be about 4E-7, translating into a total HCN fractional abundance of 2E-5 using the isotopic ratio 12C/13C=50. Multi-transitional single-dish observations further requires the H13CN fractional abundance to remain more or less constant in the envelope out to a radius of about 4E16 cm, where the HCN molecules are effectively destroyed, most probably, by photodissociation. The large amount of HCN present in the inner wind provides effective line cooling that can dominate over that generated from CO line emission. It is also shown that great care needs to be taken in the radiative transfer modelling where non-local, and non-LTE, effects are important and where the radiation field from thermal dust grains plays a major role in exciting the HCN molecules. The amount of HCN present in the circumstellar envelope around IRC+10216 is consistent with predicted photospheric values based on equilibrium chemical models and indicates that any non-equilibrium chemistry occurring in the extended pulsating atmosphere has no drastic net effect on the fractional abundance of HCN molecules that enters the outer envelope. It further suggests that few HCN molecules are incorporated into dust grains.

New SiS multi-transition (sub-)millimetre line observations of a sample of AGB stars with varying photospheric C/O-ratios and mass-loss rates are presented. A combination of low- and high-energy lines are important in constraining the circumstellar distribution of SiS molecules. A detailed radiative transfer modelling of the observed SiS line emission is performed, including the effect of thermal dust grains in the excitation analysis. We find that the circumstellar fractional abundance of SiS in these environments has a strong dependence on the photospheric C/O-ratio as expected from chemical models. The carbon stars (C/O>1) have a mean fractional abundance of 3.1E-6, about an order of magnitude higher than found for the M-type AGB stars (C/O<1) where the mean value is 2.7E-7. These numbers are in reasonable agreement with photospheric LTE chemical models. SiS appears to behave similar to SiO in terms of photodissociation in the outer part of the circumstellar envelope. In contrast to previous results for the related molecule SiO, there is no strong correlation of the fractional abundance with density in the CSE, as would be the case if freeze-out onto dust grains were important. However, possible time-variability of the line emission in the lower J transitions and the sensitivity of the line emission to abundance gradients in the inner part of the CSE may mask a correlation with the density of the wind. There are indications that the SiS fractional abundance could be significantly higher closer to the star which, at least in the case of M-type AGB stars, would require non-equilibrium chemical processes.

The S-stars have been suggested to be a brief transitional phase as stars evolve from oxygen-rich M-type stars into carbon stars, through the dredge up of carbon from He-shell burning. As possible transition objects, S-stars might help achieve a deeper understanding of the chemical evolution as a star ascends the AGB, as well as shed more light on the mass-loss mechanism. We have initiated a large survey of 40 S-stars to observe line emission in common molecules such as CO, SiO, HCN, CS and SiS. Detailed radiative transfer modelling of multi-transition CO radio line observations towards a sample of 40 S-stars shows that the mass-loss rate distribution of S-stars is consistent with those found for M-type AGB stars and carbon stars. Initial results from modelling of the circumstellar SiO emission are also presented.

We have performed high spatial resolution observations of SiO line emission for a sample of 11 AGB stars using the ATCA, VLA and SMA interferometers. Detailed radiative transfer modelling suggests that there are steep chemical gradients of SiO in their circumstellar envelopes. The emerging picture is one where the radial SiO abundance distribution starts at an initial high abundance, in the case of M-stars consistent with LTE chemistry, that drastically decreases at a radius of ~1E15 cm. This is consistent with a scenario where SiO freezes out onto dust grains. The region of the wind with low abundance is much more extended, typically ~1E16 cm, and limited by photodissociation. The surpisingly high SiO abundances found in carbon stars requires non-equilibrium chemical processes.

In the radio regime the mass-loss rate of AGB stars is best probed using molecular (and atomic) line emission arising in the CSE formed by the stellar wind. The numerical modelling of the circumstellar emission where intricate interplays between physical and chemical processes take place, is a challenge. The derived mass-loss rates depend crucially on the assumptions in the circumstellar model, of which some can be constrained if enough observational data exist. Therefore, a reliable mass-loss-rate determination for an individual star requires, in addition to a detailed radiative transfer analysis, good observational constraints in the form of multi-line observations and radial brightness distributions. Of the methods used to estimate mass-loss rates from galactic AGB stars those based on radiative transfer modelling of CO line emission are most commonly used and possibly also the most accurate. Typically, CO multi-transitional observations can constrain the mass-loss rate to better than 50%, within the adopted circumstellar model. Comparison with complementary methods, such as estimates based on dust radiative transfer modelling coupled with a dynamical model, are consistent within a factor of three.

The large quantity and high quality of modern radio and infrared line observations require efficient modeling techniques to infer physical and chemical parameters such as temperature, density, and molecular abundances. We present a computer program to calculate the intensities of atomic and molecular lines produced in a uniform medium, based on statistical equilibrium calculations involving collisional and radiative processes and including radiation from background sources. Optical depth effects are treated with an escape probability method. The program is available on the World Wide Web at http://www.sron.rug.nl/~vdtak/radex/index.shtml . The program makes use of molecular data files maintained in the Leiden Atomic and Molecular Database (LAMDA), which will continue to be improved and expanded. The performance of the program is compared with more approximate and with more sophisticated methods. An Appendix provides diagnostic plots to estimate physical parameters from line intensity ratios of commonly observed molecules. This program should form an important tool in analyzing observations from current and future radio and infrared telescopes.

This paper presents a large spectral line and continuum survey of 8 deeply embedded, low-mass protostellar cores using the Submillimeter Array. Each source was observed in high excitation lines of some of the most common molecular species, CO, HCO+, CS, SO, H2CO, CH3OH and SiO. Line emission from 11 species originating from warm and dense gas have been imaged at high angular resolution (1-3"; typically 200-600 AU) together with continuum emission at 230 GHz (1.3 mm) and 345 GHz (0.8 mm). Compact continuum emission is observed for all sources which likely originates in marginally optically thick circumstellar disks, with typical lower limits to their masses of 0.1 M_sun (1-10% of the masses of their envelopes) and having a dust opacity law with beta approximately 1. Prominent outflows are present in CO 2-1 observations in all sources: the most diffuse outflows are found in the sources with the lowest ratios of disk-to-envelope mass, and it is suggested that these sources are in a phase where accretion of matter from the envelope has almost finished and the remainder of the envelope material is being dispersed by the outflows. Other characteristic dynamical signatures are found with inverse P Cygni profiles indicative of infalling motions seen in the 13CO 2-1 lines toward NGC1333-IRAS4A and -IRAS4B. Outflow-induced shocks are present on all scales in the protostellar environments and are most clearly traced by the emission of CH3OH in NGC1333-IRAS4A and -IRAS4B. These observations suggest that the emission of CH3OH and H2CO from these proposed "hot corinos" are related to the shocks caused by the protostellar outflows. Only one source, NGC1333-IRAS2A, has evidence for hot, compact CH3OH emission coincident with the embedded protostar.

We have used a detailed non-LTE radiative transfer code to model new APEX CO(J=3-2) data, and existing CO radio line data, on a sample of 40 AGB S-stars. The derived mass-loss-rate distribution has a median value of 2E-7 Msun/yr, and resembles values obtained for similar samples of M-stars and carbon stars. Possibly, there is a scarcity of high-mass-loss-rate (>=1E-5 Msun/yr) S-stars. The distribution of envelope gas expansion velocities is similar to that of the M-stars, the median is 7.5 km/s, while the carbon stars, in general, have higher gas expansion velocities. The mass-loss rate correlates well with the gas expansion velocity, in accordance with results for M-stars and carbon stars.

This paper presents new APEX submillimetre molecular line observations of three low-mass protostars, IRS7A, IRS7B, and IRAS32, in the R Coronae Australis molecular cloud complex. The molecular excitation analysis is performed using a statistical equilibrium radiative transfer code. The derived beam averaged fractional abundances vary by less than a factor of two among the three sources, except those of H2CO and CH3OH, which show differences of about an order of magnitude. The molecular abundances are similar to those typically found in other star-forming regions in the Galaxy, such as the $\rho$~Oph and Perseus molecular clouds. There is a marked difference in the kinetic temperatures derived for the protobinary source IRS7 from H2CO (40-60 K) and CH3OH (20 K), possibly indicating a difference in origin of the emission from these two molecules.

New interferometric observations of SiO J=5-4 circumstellar line emission around the carbon star IRC+10216, using the Submillimeter Array, are presented. Complemented by multi-transition single-dish observations, including infrared observations of ro-vibrational transitions, detailed radiative transfer modelling suggests that the fractional abundance of SiO in the inner part of the envelope, between approximately 3-8 stellar radii, is as high as 1.5E-6. This is more than an order of magnitude higher than predicted by equilibrium stellar atmosphere chemistry in a carbon-rich environment and indicative of the importance of non-LTE chemical processes. In addition to the compact component, a spatially more extended (2.4E16 cm) low-fractional-abundance (1.7E-7) region is required to fit the observations. This suggests that the majority of the SiO molecules are effectively accreted onto dust grains in the inner wind while the remaining gas-phase molecules are eventually photodissociated at larger distances. Evidence of departure from a smooth wind is found in the observed visibilities, indicative of density variations of a factor 2 to 5 on an angular scale corresponding to a time scale of about 200 years. Additionally, constraints on the velocity structure of the wind are obtained.

New SiO multi-transition millimetre line observations of a sample of carbon stars, including J=8-7 observations with the APEX telescope, are used to probe the role of non-equilibrium chemistry and the influence of grains in circumstellar envelopes of carbon stars. A detailed radiative transfer modelling, including the effect of dust emission in the excitation analysis, of the observed SiO line emission is performed. A combination of low- and high-energy lines are important in constraining the abundance distribution. It is found that the fractional abundance of SiO in these C-rich environments can be several orders of magnitude higher than predicted by equilibrium stellar atmosphere chemistry. In fact, the SiO abundance distribution of carbon stars closely mimic that of M-type (O-rich) AGB stars. A possible explanation for this behaviour is a shock-induced chemistry, but also the influence of dust grains, both as a source for depletion as well as production of SiO, needs to be further investigated. As observed for M-type AGB stars, a clear trend that the SiO fractional abundance decreases as the mass-loss rate of the star increases is found for the carbon stars. This indicates that SiO is accreted onto dust grains in the circumstellar envelopes.

Mid-infrared (23-35 micron) emission from the deeply embedded "Class 0" protostar IRAS16293-2422 is detected with the Spitzer Space Telescope infrared spectrograph. A detailed radiative transfer model reproducing the full spectral energy distribution (SED) from 23 micron to 1.3 mm requires a large inner cavity of radius 600 AU in the envelope to avoid quenching the emission from the central sources. This is consistent with a previous suggestion based on high angular resolution millimeter interferometric data. An alternative interpretation using a 2D model of the envelope with an outflow cavity can reproduce the SED but not the interferometer visibilities. The cavity size is comparable to the centrifugal radius of the envelope and therefore appears to be a natural consequence of the rotation of the protostellar core, which has also caused the fragmentation leading to the central protostellar binary. With a large cavity such as required by the data, the average temperature at a given radius does not increase above 60-80 K and although hot spots with higher temperatures may be present close to each protostar, these constitute a small fraction of the material in the inner envelope. The proposed cavity will also have consequences for the interpretation of molecular line data, especially of complex species probing high temperatures in the inner regions of the envelope.

The Submillimeter Array has opened a new window to study the innermost warm and dense regions of the envelopes and disks around deeply embedded protostars. This paper presents high-angular resolution (< 2'') submillimeter observations of the class 0 young stellar object NGC1333-IRAS2A. Dust continuum emission and lines of complex organic molecules such as CH3OCH3 and CH3OCHO, high excitation CH3OH transitions, deuterated methanol CH3OD as well as lines of CO, HCN, H13CN, SO and SO2 are detected on < 200 AU scales. The observations are interpreted using detailed radiative transfer models of the physical and chemical structure, consistent with both single-dish and interferometer data. The continuum emission is explained by an extended envelope and a compact but resolved component, presumably a circumstellar disk with a diameter of 200-300 AU and a mass of a few times 0.01-0.1 M_sun. If related to the rotation of the envelope, then the size of this disk suggests a centrifugal barrier of 200-300 AU, which implies that the temperature in the envelope does not increase above 100 K. Its large size also suggests that the build-up of disks proceeds rapidly throughout the early protostellar stages. The smaller (< 100 AU) disks found around other deeply embedded protostars may be a result of tidal truncation. The high-resolution observations of SO can be explained with a simple constant abundance, ~1E-9, constrained through single-dish observations, whereas those of H13CN and the organic species require high abundances, increased by one to two orders of magnitude, or an additional compact source of emission at small scales. The compact molecular emission could originate in a hot core region of the inner envelope, but a more likely reservoir is the circumstellar disk.

This paper presents the third in a series of single-dish studies of molecular abundances in the envelopes around a large sample of 18 low-mass pre- and protostellar objects. It focuses on typical grain mantle products and organic molecules, including H2CO, CH3OH and CH3CN. With a few exceptions, all H2CO lines can be fit by constant abundances of 7e(-11)-8e(-9) throughout the envelopes if ortho- and para lines are considered independently. The current observational dataset does not require a large H2CO abundance enhancement in the inner warm regions, but this can also not be ruled out. Through comparison of the H2CO abundances of the entire sample, the H2CO ortho-para ratio is constrained to be 1.6+/-0.3 consistent with thermalization on grains at temperatures of 10-15 K. The H2CO abundances can be related to the empirical chemical network established on the basis of our previously reported survey of other species and is found to be closely correlated with that of the nitrogen-bearing molecules. These correlations reflect the freeze-out of molecules at low temperatures and high densities. An improved fit to the data is obtained with a `drop' abundance structure in which the abundance decrease when the temperature is lower than the evaporation temperature and the density high enough so that the timescale for depletion is less than the lifetime of the core. The derived abundance structures for CH3OH are consistent with a scenario where CH3OH probes the action of compact outflows on the envelopes, which is further supported by comparison to high frequency, high excitation CS J=10-9 and HDO line profiles. The extent to which the outflow dominates the abundance enhancements compared with the passively heated inner envelope depends on the filling factors of the two components in the observing beam.

The nature of the mechanism responsible for producing the spectacular, geometrically thin, spherical shells found around some carbon stars has been an enigma for some time. Based on extensive radiative transfer modelling of both CO line emission and dust continuum radiation for all objects with known detached molecular shells, we present compelling evidence that these shells show clear signs of interaction with a surrounding medium. The derived masses of the shells increase with radial distance from the central star while their velocities decrease. A simple model for interacting winds indicate that the mass-loss rate producing the faster moving wind has to be almost two orders of magnitudes higher (~10^-5 solar masses per year) than the slower AGB wind (a few 10^-7 solar masses per year) preceding this violent event. At the same time, the present-day mass-loss rates are very low indicating that the epoch of high mass-loss rate was relatively short, on the order of a few hundred years. This, together with the number of sources exhibiting this phenomenon, suggests a connection with He-shell flashes (thermal pulses). We report the detection of a detached molecular shell around the carbon star DR Ser, as revealed from new single-dish CO (sub-)millimetre line observations. The properties of the shell are similar to those characterising the young shell around U Cam.

Radiative transfer models of multi-transition observations are used to determine molecular abundances as functions of position in pre- and protostellar cores. The data require a "drop" abundance profile with radius, with high abundances in the outermost regions probed by low excitation 3 mm lines, and much lower abundances at intermediate zones probed by higher frequency lines. The results are illustrated by detailed analysis of CO and HCO+ lines for a subset of objects. We propose a scenario in which the molecules are frozen out in a region of the envelope where the temperature is low enough (<= 40 K) to prevent immediate desorption, but where the density is high enough (> 10^4-10^5 cm^{-3}) that the freeze-out timescales are shorter than the lifetime of the core. The size of the freeze-out zone is thereby a record of the thermal and dynamical evolution of the cores. Fits to CO data for a sample of 16 objects indicate that the size of the freeze-out zone decreases significantly between Class 0 and I objects, explaining the variations in, for example, CO abundances with envelope masses. However, the corresponding timescales are 10^{5+/-0.5} years, with no significant difference between Class 0 and I objects. These timescales suggest that the dense pre-stellar phase with heavy depletions lasts only a short time, of order 10^5 yr, in agreement with recent chemical-dynamical models.

The results from a single-dish molecular line survey of a set of 18 deeply embedded young stellar objects are summarized. More than 40 lines from 16 different species were observed with the JCMT, Onsala, IRAM 30m and SEST telescopes. The multi-transition data are analyzed using a temperature and density structure derived from models of the dust continuum emission. For the outer envelope (>300 AU), the data indicate a `drop' abundance profile for many species, with normal abundances in the outer- and innermost regions and highly depleted abundances in an intermediate zone. This zone is bounded at the outer than the lifetime of the core, and at the inner edge by the evaporation temperature of the species involved. In the innermost envelope (<300 AU), all ices evaporate resulting in jumps in the abundances of complex organic molecules such as CH3OH. A key project for Herschel will be to survey gas-phase water in these objects, whose abundance shows extreme variations with temperature. ALMA wil be able to directly image the chemical variations throughout the envelope and zoom in on the inner hot core and protoplanetary disks on scales on tens of AU.

Atomic and molecular data for the transitions of a number of astrophysically interesting species are summarized, including energy levels, statistical weights, Einstein A-coefficients and collisional rate coefficients. Available collisional data from quantum chemical calculations and experiments are extrapolated to higher energies. These data, which are made publically available through the WWW at http://www.strw.leidenuniv.nl/~moldata, are essential input for non-LTE line radiative transfer programs. An online version of a computer program for performing statistical equilibrium calculations is also made available as part of the database. Comparisons of calculated emission lines using different sets of collisional rate coefficients are presented. This database should form an important tool in analyzing observations from current and future (sub)millimetre and infrared telescopes.

Two bipolar proto-planetary nebulae, IRAS16594-4656 and IRAS17150-3224, have been detected in various molecular lines, namely CO, ^{13}CO, HCN and CN, and remain undetected in several other species. CO(J=2-1) and CO(J=3-2) line profiles are compared to new spectra of similar PPN candidates, previously undetected in CO(J=2-1): CPD-53^{o}5736, IRAS17106-3046, IRAS17245-3951 and IRAS17441-2411. CO(J=2-1) maps of IRAS16594-4656 and IRAS17150-3224 show that both PPNe have little separation between blue, centre and red-shifted emission, and also that the CO(J=2-1) emission is of a similar size to the telescope beam. Fractional abundances of all detected molecules (except CO) are calculated using the results of CO line modelling and a simple photodissociation model. For those species not detected, upper limits are derived. Comparisons between these fractional abundances and those of other PPNe show that IRAS16594-4656 and IRAS17150-3224 are quite under-abundant when compared to molecule-rich sources like CRL618, CRL2688 and OH231.8+4.2. As a reason for this deficit, the difference in circumstellar envelope/torus density between the molecule-rich sources and the molecule-poor sources is proposed, and supported by a chemical model which follows the transition of a circumstellar envelope through the AGB phase and into the PPN phase of evolution. The model includes the effects of UV radiation, cosmic rays and also X-rays. Finally, the post-AGB ages of these two objects (200-400 yr) are estimated using CN/HCN and HCN/CO ratios and both ages are found to be in agreement with previous figures cited in the literature, IRAS17150-3224 being the younger of the two PPNe.

High angular resolution Australia Telescope Compact Array (ATCA) observations of SiO `thermal' millimetre line emission towards the two oxygen-rich, low mass loss rate AGB stars R Dor and L2 Pup are presented. In both cases the emission is resolved with an overall spherical symmetry. Detailed radiative transfer modelling of the SiO line emission has been performed, and the comparison between observations and models are conducted in the visibility plane, maximizing the sensitivity. The excitation analysis suggests that the abundance of SiO is as high as 4E-5 in the inner part of the wind, close to the predicted values from stellar atmosphere models. Beyond a radius of ~ 1E15 cm the SiO abundance is significantly lower, about 3E-6, until it decreases strongly at a radius of about 3E15 cm. This is consistent with a scenario where SiO first freezes out onto dust grains, and then eventually becomes photodissociated by the interstellar UV-radiation field. In these low expansion velocity sources the turbulent broadening of the lines plays an important role in the line formation. Micro-turbulent velocity widths in the range 1.1-1.5 km/s result in a very good reproduction of the observed line shapes even if the gas expansion velocity is kept constant. This, combined with the fact that the SiO and CO lines are well fitted using the same gas expansion velocity (to within 5-10%), suggest that the envelope acceleration occurs close to the stellar photosphere, within 20-30 stellar radii.

We present detailed gas-phase chemical models for the envelope of the low-mass star-forming region IRAS 16293-2422. By considering both time- and space-dependent chemistry, these models are used to study both the physical structure proposed by Schoier et al. (2002), as well as the chemical evolution of this region. A new feature of our study is the use of a detailed, self-consistent radiative transfer model to translate the model abundances into line strengths and compare them directly with observations of a total of 76 transitions for 18 chemical species, and their isotopes. The model can reproduce many of the line strengths observed within 50%. The best fit is for times in the range of 3e3 - 3e4 yrs, and requires only minor modifications to our model for the high-mass star-forming region AFGL 2591. The ionization rate for the source may be higher than previously expected -- either due to an enhanced cosmic-ray ionization rate, or, more probably, to the presence of X-ray induced ionization from the center. A significant fraction of the CO is found to desorb in the temperature range of 15-40 K; below this temperature, \~90% or more of the CO is frozen out. The inability of the model to explain the HCS+, C2H, and OCS abundances suggests the importance of further laboratory studies of basic reaction rates. Finally, predictions of the abundances and spatial distributions of other species which could be observed by future facilities (e.g., Herschel-HIFI, SOFIA, millimeter arrays) are provided.

High angular resolution H2CO 218 GHz line observations have been carried out toward the low-mass protostars IRAS 16293-2422 and L1448-C using the Owens Valley Millimeter Array at ~2" resolution. Simultaneous 1.37 mm continuum data reveal extended emission which is compared with that predicted by model envelopes constrained from single-dish data. For L1448-C the model density structure works well down to the 400 AU scale to which the interferometer is sensitive. For IRAS 16293-2422, a known proto-binary object, the interferometer observations indicate that the binary has cleared much of the material in the inner part of the envelope, out to the binary separation of ~800 AU. The morphology and velocity structure of the H2CO array data have been used to investigate whether the abundance enhancements inferred from single-dish modelling are due to thermal evaporation of ices or due to liberation of the ice mantles by shocks in the inner envelope. A scenario in which the H2CO abundance drops in the cold dense part of the envelope where CO is frozen out but is undepleted in the outermost region provides good fits to the single-dish and interferometer data on short baselines for both sources. Emission on the longer baselines is best reproduced if the H2CO abundance is increased by about an order of magnitude from ~1E-10 to ~1E-9 in the inner parts of the envelope due to thermal evaporation when the temperature exceeds ~50K. Other scenarios, including weak outflow-envelope interactions and photon heating of the envelope, are discussed and predictions for future generation interferometers are presented, illustrating their potential in distinguishing these competing scenarios.

This paper presents the first substantial study of the chemistry of the envelopes around a sample of 18 low-mass pre- and protostellar objects for which physical properties have previously been derived from radiative transfer modeling of their dust continuum emission. Single-dish line observations of 24 transitions of 9 molecular species (not counting isotopes) including HCO+, N2H+, CS, SO, SO2, HCN, HNC, HC3N and CN are reported. The line intensities are used to constrain the molecular abundances by comparison to Monte Carlo radiative transfer modeling of the line strengths. An empirical chemical network is constructed on the basis of correlations between the abundances of various species. For example, it is seen that the HCO+ and CO abundances are linearly correlated, both increasing with decreasing envelope mass. Species such as CS, SO and HCN show no trend with envelope mass. In particular no trend is seen between ``evolutionary stage'' of the objects and the abundances of the main sulfur- or nitrogen-containing species. Among the nitrogen-bearing species abundances of CN, HNC and HC3N are found to be closely correlated, which can be understood from considerations of the chemical network. The CS/SO abundance ratio is found to correlate with the abundances of CN and HC3N, which may reflect a dependence on the atomic carbon abundance. An anti-correlation is found between the deuteration of HCO+ and HCN, reflecting different temperature dependences for gas-phase deuteration mechanisms. The abundances are compared to other protostellar environments. In particular it is found that the abundances in the cold outer envelope of the previously studied class 0 protostar IRAS16293-2422 are in good agreement with the average abundances for the presented sample of class 0 objects.

We present a survey of the formaldehyde emission in a sample of eight Class 0 protostars obtained with the IRAM and JCMT telescopes. The data have been analyzed with three different methods with increasing level of sophistication. We first analyze the observed emission in the LTE approximation, and derive rotational temperatures between 11 and 40 K, and column densities between 1 and 20 x 10^13 cm^-2. Second, we use a LVG code and derive larger kinetic temperatures, between 30 and 90 K, consistent with subthermally populated levels and densities from 1 to 6 x 10^5 cm^-3. The column densities from the LVG modeling are within a factor of 10 with respect to those derived in the LTE approximation. Finally, we analyze the observations based upon detailed models for the envelopes surrounding the protostars, using temperature and density profiles previously derived from continuum observations. We approximate the formaldehyde abundance across the envelope with a jump function, the jump occurring when the dust temperature reaches 100 K, the evaporation temperature of the grain mantles. The observed formaldehyde emission is well reproduced only if there is a jump, more than two orders of magnitude, in four sources. In the remaining four sources the data are consistent with a formaldehyde abundance jump, but the evidence is more marginal (~2 sigma). The inferred inner H2CO abundance varies between 1 x 10^-8 and 6 x 10^-6. We discuss the implications of these jumps for our understanding of the origin and evolution of ices in low mass star forming regions. Finally, we give predictions for the submillimeter H2CO lines, which are particularly sensitive to the abundance jumps.

This paper presents a detailed study of the chemistry in the outflow associated with the low-mass protostar NGC1333-IRAS2A down to 3" (650 AU) scales. Millimeter-wavelength aperture-synthesis observations from the OVRO and BIMA interferometers and (sub)millimeter single-dish observations from the Onsala 20m telescope and CSO are presented. The interaction of the highly collimated protostellar outflow with a molecular condensation ~15000 AU from the central protostar is clearly traced by molecular species such as HCN, SiO, SO, CS, and CH3OH. Especially SiO traces a narrow high velocity component at the interface between the outflow and the molecular condensation. Multi-transition single-dish observations are used to distinguish the chemistry of the shock from that of the molecular condensation and to address the physical conditions therein. Statistical equilibrium calculations reveal temperatures of 20 and 70 K for the quiescent and shocked components, respectively, and densities near 10^6 cm^{-3}. Significant abundance enhancements of two to four orders of magnitude are found in the shocked region for molecules such as CH3OH, SiO and the sulfur-bearing molecules. HCO+ is seen only in the aftermath of the shock consistent with models where it is destroyed through release of H2O from grain mantles in the shock. N2H+ shows narrow lines, not affected by the outflow but rather probing the ambient cloud. Differences in abundances of HCN, H2CO and CS are seen between different outflow regions and are suggested to be related to differences in the atomic carbon abundance. Compared to the warm inner parts of protostellar envelopes, higher abundances of in particular CH3OH and SiO are found in the outflows, which may be related to density differences between the regions.

This paper investigates small-scale structures of dense gas and dust around the low-mass protostellar binary NGC1333-IRAS2 using millimeter-wavelength aperture-synthesis observations from the OVRO and BIMA interferometers. The detected 3 mm continuum emission is consistent with models of the envelope around IRAS2A, based on previously reported submillimeter continuum images, down to the 3", or 500 AU, resolution of the interferometer data. Our data constrain the contribution of an unresolved point source to 22 mJy. Within the accuracy of the parameters describing the envelope model, the point source flux has an uncertainty by up to 25%. We interpret this point source as a cold disk of mass \gtrsim 0.3 M_\odot. The same envelope model also reproduces aperture-synthesis line observations of the optically thin isotopic species C34S and H13CO+. The more optically thick main isotope lines show a variety of components in the protostellar environment: N2H+ is closely correlated with dust concentrations as seen at submillimeter wavelengths and is particularly strong toward the starless core IRAS2C. We hypothesize that N2H+ is destroyed through reactions with CO that is released from icy grains near the protostellar sources IRAS2A and B. CS, HCO+, and HCN have complex line shapes apparently affected by both outflow and infall. In addition to the east-west jet from IRAS2A, a north-south velocity gradient near this source indicates a second, perpendicular outflow. This suggests the presence of a binary companion within 0.3" (65 AU) from IRAS2A as driving source of this outflow. Alternative explanations of the velocity gradient such as rotation in a circumstellar envelope or a single, wide-angle outflow are less likely.

We present radio-interferometric observations of HCN J=1-0 line emission from the carbon star R Scl, obtained with the interim 3-mm receivers of the Australia Telescope Compact Array. The emission is resolved into a central source with a Gaussian FWHM of ~1 arcsec, which we identify as the present mass loss envelope. Using a simple photodissociation model and constraints from single-dish HCN spectra, we argue that the present mass-loss rate is low, ~2 x 10^{-7} Msol/yr, supporting the idea that R Scl had to experience a brief episode of intense mass loss in order to produce the detached CO shell at ~10 arcsec radius inferred from single-dish observations. Detailed radiative transfer modelling yields an abundance of HCN relative to H_2, of ~10^{-5} in the present-day wind. There appears to be a discrepancy between model results obtained with higher transition single-dish data included and those from the J=1-0 interferometer data alone, in that the interferometer data suggest a smaller envelope size and larger HCN abundance than the single-dish data. The lack of HCN in the detached shell is consistent with the rapid photodissociation of HCN into CN as it expands away from the star.

An extensive radiative transfer analysis of circumstellar SiO `thermal' radio line emission from a large sample of M-type AGB stars has been performed. The sample contains 18 irregulars of type Lb (IRV), 7 and 34 semiregulars of type SRa and SRb (SRV), respectively, and 12 Miras. New observational material, which contains data in several SiO rotational lines in the ground vibrational state, is presented for the IRV/SRVs. The detection rate was about 60%. SiO fractional abundances have been determined through radiative transfer modelling. The abundance distribution of the IRV/SRV sample has a median value of 8E-6. The high mass loss rate Miras have a bimodal abundance distribution, a low abundance group (on average 4E-7) and a high abundance group (on average 7E-6). The derived SiO abundances are in all cases well below the abundance expected from stellar atmosphere chemistry, on average by a factor of ten. In addition, there is a trend of decreasing SiO abundance with increasing mass loss rate. This is interpreted in terms of depletion of SiO molecules by the formation of silicate grains in the circumstellar envelopes, with an efficiency which increases with the mass loss rate. The SiO and CO radio line profiles differ in shape. In general, the dominating part of the SiO profiles are narrower than the CO profiles, but they have low-intensity wings which cover the full velocity range of the CO profile. This is interpreted (as has been done also by others) as due to the influence of gas acceleration in the region which produces a significant fraction of the SiO line emission. Finally, a number of sources which have peculiar CO line profiles are discussed from the point of view of their SiO line properties.

We present a millimetre molecular line survey of seven high mass loss rate carbon stars in both the northern and southern skies. A total of 196 emission lines (47 transitions) from 24 molecular species were detected. The observed CO emission is used to determine mass-loss rates and the physical structure of the circumstellar envelope, such as the density and temperature structure, using a detailed radiative transfer analysis. This enables abundances for the remaining molecular species to be determined. The derived abundances generally vary between the sources by no more than a factor of five indicating that circumstellar envelopes around carbon stars with high mass-loss rates have similar chemical compositions. However, there are some notable exceptions. The most striking difference between the abundances are reflecting the spread in the 12C/13C-ratio of about an order of magnitude between the sample stars, mainly reflecting the results of nucleosynthesis. The abundance of SiO also shows a variation of more than an order of magnitude between the sources and is on the average more than an order of magnitude more abundant than predicted from photospheric chemistry in thermal equilibrium. The over abundance of SiO is consistent with dynamical modelling of the stellar atmosphere and the inner parts of the wind where a pulsation-driven shock has passed. This scenario is possibly further substantiated by the relatively low amount of CS present in the envelopes. The chemistry occurring in the outer envelope is consistent with current photochemical models.

General results from a 3-5 micron spectroscopic survey of nearby low-mass young stellar objects are presented. L and M-band spectra have been obtained of \~50 low mass embedded young stars using the ISAAC spectrometer mounted on UT1-Antu at Paranal Observatory. For the first time, a consistent census of the CO, H2O ices and the minor ice species CH3OH and OCN- and warm CO gas present around young stars is obtained, using large number statistics and resolving powers of up to R=10000. The molecular structure of circumstellar CO ices, the depletion of gaseous CO onto grains in protoplanetary disks, the presence of hot gas in the inner parts of circumstellar disks and in outflows and infalls are studied. Furthermore, the importance of scattering effects for the interpretation of the spectra have been addressed.

High angular resolution H2CO line observations have been carried out for the low mass proto-binary star IRAS 16293-2422 using the OVRO millimeter array. Simultaneous continuum observations reveal that the source most probably contains two accretion disks, less than 250 AU in size, through which matter is fed onto the stars. The observations suggest that the binary has cleared most of the material in the inner part of the envelope, out to the binary separation. The H2CO interferometer observations indicate the presence of large scale rotation roughly perpendicular to the large scale CO outflow associated with this source. The H2CO emission is dominated by the compact dense and hot (T>100 K) gas close to the positions of the continuum peaks.

A study of the chemical structure of the envelopes around a sample protostars is introduced. Physical models for the envelopes derived using 1D radiative transfer modeling of their dust continuum emission are used as input for Monte Carlo modeling of single-dish line observations in order to establish the chemical inventories for the different sources. CO and HCO+ abundances are found to be correlated with envelope mass: the sources with the most massive envelopes show the lowest abundances, supporting the idea that the depletion of these molecules is most efficient in the colder, denser environments. Other molecules like CS and HCN do not show a similar trend. Deuterated species like DCO+ and DCN also show signs of correlations with envelope mass and temperature: the deuterium fractionation of HCO+ is most prominent for the coldest and most massive envelopes in the sample and at the same time anti-correlated with the fractionation of HCN. This puts constraints on the low temperature gas-phase deuterium chemistry.

Molecular line data for the rotational transitions of a number of astrophysically interesting species have been collected and are made publically available through the www. These data, including energy levels, statistical weights, Einstein A-coefficients and collisional rate coefficients, are essential input for non-LTE molecular radiative transfer programs. A computer program for performing statistical equilibrium calculations is made publically available for use as well. This database should form an important tool in analyzing observations from current and future infrared and (sub)millimetre telescopes.

Comparison is made between a number of independent computer programs for radiative transfer in molecular rotational lines. The test models are spherically symmetric circumstellar envelopes with a given density and temperature profile. The first two test models have a simple power law density distribution, constant temperature and a fictive 2-level molecule, while the other two test models consist of an inside-out collapsing envelope observed in rotational transitions of HCO+. For the 2-level molecule test problems all codes agree well to within 0.2%, comparable to the accuracy of the individual codes, for low optical depth and up to 2% for high optical depths (tau=4800). The problem of the collapsing cloud in HCO+ has a larger spread in results, ranging up to 12% for the J=4 population. The spread is largest at the radius where the transition from collisional to radiative excitation occurs. The resulting line profiles for the HCO+ J=4-3 transition agree to within 10%, i.e., within the calibration accuracy of most current telescopes. The comparison project and the results described in this paper provide a benchmark for future code development, and give an indication of the typical accuracy of present day calculations of molecular line transfer.

We present a 4.5-4.85 micron R=5000 spectrum of the low mass class I young stellar object GSS 30 IRS1 (L=25L_sun) in the rho Ophiuchus core, observed with the infrared spectrometer (ISAAC) on the Very Large Telescope (VLT-UT1). Strong line emission from the ro-vibrational transitions of 12CO and 13CO is detected. In total more than 40 distinct lines are seen in the covered region. The line emission is spatially extended and detected up to 2" = 320 AU from the central source but is spectrally unresolved (Delta v < 30 km/s). This is the first time strong emission in the fundamental ro-vibrational band from CO has been observed from an embedded young stellar object. The line fluxes were modeled using a 1-dimensional full radiative transfer code, which shows that the emission is fully consistent with a gas in LTE at a single well constrained temperature (T=515+/-5 K). Furthermore, the ratios between lines from the two detected isotopic species of CO show that the 12CO lines must be optically thick. However, this is inconsistent with the observed spatial extent of the emission, since this implies such low CO column densities that the lines are optically thin. A likely solution to the discrepancy is that the lines are emitted by a smaller more dense region and then scattered in the bipolar cavity present around the central star. This gives a rough estimate of the total molecular gas mass of 1-100 M_earth and a physical extent of ~20-100 AU. We propose that the most likely origin of the line emission is post-shocked gas in a dense dissociative accretion shock from the inner 10-50 AU of a circumstellar disk. The presence of a shock capable of dissociating molecules in the disk will have implications for the chemical evolution of disks around young low mass stars.

We have determined mass loss rates and gas expansion velocities for a sample of 69 M-type irregular (IRV; 22 objects) and semiregular (SRV; 47 objects) AGB-variables using a radiative transfer code to model their circumstellar CO radio line emission. We believe that this sample is representative for the mass losing stars of this type. The (molecular hydrogen) mass loss rate distribution has a median value of 2.0E-7 solar masses per year. M-type IRVs and SRVs with a mass loss rate in excess of 5E-7 solar masses per year must be very rare, and among these mass losing stars the number of sources with mass loss rates below a few 10E-8 solar masses per year must be small. We find no significant difference between the IRVs and the SRVs in terms of their mass loss characteristics. Among the SRVs the mass loss rate shows no dependence on the period. Likewise the mass loss rate shows no correlation with the stellar temperature. The gas expansion velocity distribution has a median of 7.0 km/s. The mass loss rate and the gas expansion velocity correlate well, a result in line with theoretical predictions for an optically thin, dust-driven wind. In general, the model produces line profiles which acceptably fit the observed ones. We have compared the results of this M-star sample with a similar C-star sample analysed in the same way. The mass loss rate characteristics are very similar for the two samples. On the contrary, the gas expansion velocity distributions are clearly different. In particular, the number of low-velocity sources is much higher in the M-star sample. We found no example of the sharply double-peaked CO line profile, which is evidence of a large, detached CO-shell, among the M-stars. About 10% of the C-stars show this phenomenon.

We have modelled the circumstellar envelope of {\it o} Ceti (Mira) using new observational constraints. These are obtained from photospheric light scattered in near-IR vibrational-rotational lines of circumstellar CO molecules at 4.6 micron: absolute fluxes, the radial dependence of the scattered intensity, and two line ratios. Further observational constraints are provided by ISO observations of far-IR emission lines from highly excited rotational states of the ground vibrational state of CO, and radio observations of lines from rotational levels of low excitation of CO. A code based on the Monte-Carlo technique is used to model the circumstellar line emission. We find that it is possible to model the radio and ISO fluxes, as well as the highly asymmetric radio-line profiles, reasonably well with a spherically symmetric and smooth stellar wind model. However, it is not possible to reproduce the observed NIR line fluxes consistently with a `standard model' of the stellar wind. This is probably due to incorrectly specified conditions of the inner regions of the wind model, since the stellar flux needs to be larger than what is obtained from the standard model at the point of scattering, i.e., the intermediate regions at approximately 100-400 stellar radii (2"-7") away from the star. Thus, the optical depth in the vibrational-rotational lines from the star to the point of scattering has to be decreased. This can be accomplished in several ways. For instance, the gas close to the star (within approximately 2") could be in such a form that light is able to pass through, either due to the medium being clumpy or by the matter being in radial structures (which, further out, developes into more smooth or shell-like structures).

A detailed radiative transfer analysis of the observed continuum and molecular line emission toward the deeply embedded young stellar object IRAS 16293-2422 is performed. The continuum modelling is used to constrain the temperature and density distributions in the envelope, enabling quantitative estimates of various molecular abundances. The molecular excitation analysis reveals that the emission from some molecular species is well reproduced assuming a constant fractional abundance throughout the envelope. The abundances and isotope ratios are generally close to typical values found in cold molecular clouds in these cases, and there is a high degree of deuterium fractionation. There are, however, a number of notable exceptions. Lines covering a wide range of excitation conditions indicate for some molecules, e.g., H2CO, CH3OH, SO, SO2 and OCS, a drastic increase in their abundances in the warm and dense inner region of the circumstellar envelope. The location at which this increase occurs is consistent with the radius at which ices are expected to thermally evaporate off the grains. In all, there is strong evidence for the presence of a `hot core' close to the protostar, whose physical properties are similar to those detected towards most high mass protostars except for a scaling factor. However, the small scale of the hot gas and the infalling nature of the envelope lead to very different chemical time scales between low mass and high mass hot cores, such that only very rapidly produced second-generation complex molecules can be formed in IRAS 16293-2422. Alternatively, the ices may be liberated due to grain-grain collisions in turbulent shear zones where the outflow interacts with the envelope.

An extensive modelling of CO line emission from the circumstellar envelopes around a number of carbon stars is performed. By combining radio observations and infrared observations obtained by ISO the circumstellar envelope characteristics are probed over a large radial range. In the radiative transfer analysis the observational data are consistently reproduced assuming a spherically symmetric and smooth wind expanding at a constant velocity. The combined data set gives better determined envelope parameters, and puts constraints on the mass loss history of these carbon stars. The importance of dust in the excitation of CO is addressed using a radiative transfer analysis of the observed continuum emission, and it is found to have only minor effects on the derived line intensities. The analysis of the dust emission also puts further constraints on the mass loss rate history. The stars presented here are not likely to have experienced any drastic long-term mass loss rate modulations, at least less than a factor of about 5, over the past thousands of years. Only three, out of nine, carbon stars were observed long enough by ISO to allow a detection of CO far-infrared rotational lines.

We present 1D radiative transfer modelling of the envelopes of a sample of 18 low-mass protostars and pre-stellar cores with the aim of setting up realistic physical models, for use in a chemical description of the sources. The density and temperature profiles of the envelopes are constrained from their radial profiles obtained from SCUBA maps at 450 and 850 micron and from measurements of the source fluxes ranging from 60 micron to 1.3 mm. The densities of the envelopes within ~10000 AU can be described by single power-laws r^{-p} for the class 0 and I sources with p ranging from 1.3 to 1.9, with typical uncertainties of +/- 0.2. Four sources have flatter profiles, either due to asymmetries or to the presence of an outer constant density region. No significant difference is found between class 0 and I sources. The power-law fits fail for the pre-stellar cores, supporting recent results that such cores do not have a central source of heating. The derived physical models are used as input for Monte Carlo modelling of submillimeter C18O and C17O emission. It is found that class I objects typically show CO abundances close to those found in local molecular clouds, but that class 0 sources and pre-stellar cores show lower abundances by almost an order of magnitude implying that significant depletion occurs for the early phases of star formation. While the 2-1 and 3-2 isotopic lines can be fitted using a constant fractional CO abundance throughout the envelope, the 1-0 lines are significantly underestimated, possibly due to contribution of ambient molecular cloud material to the observed emission. The difference between the class 0 and I objects may be related to the properties of the CO ices.

Using a detailed radiative transfer analysis, combined with an energy balance equation for the gas, we have performed extensive modelling of circumstellar CO radio line emission from a large sample of optically bright carbon stars. We determine some of the basic parameters that characterize circumstellar envelopes (CSEs), e.g., the stellar mass loss rate, the gas expansion velocity, and the kinetic temperature structure of the gas. The derived mass loss rates span almost four orders of magnitude, from 5E-9 up to 2E-5 solar masses per year, with the median mass loss rate being 3E-7 solar masses per year. We estimate that the estimated mass loss rates are typically accurate to 50% within the adopted circumstellar model. The physical conditions prevailing in the CSEs vary considerably over such a large range of mass loss rates. Among other things, it appears that the dust-to-gas mass ratio and/or the dust properties change with the mass loss rate. We find that the mass loss rate and the gas expansion velocity are well correlated, and that both of them clearly depend on the pulsational period and (with larger scatter) the stellar luminosity. Moreover, the mass loss rate correlates weakly with the stellar effective temperature, in the sense that the cooler stars tend to have higher mass loss rates, but there seems to be no correlation with the stellar C/O-ratio. We conclude that the mass loss rate increases with increased regular pulsation and/or luminosity, and that the expansion velocity increases as an effect of increasing mass loss rate (for low mass loss rates) and luminosity.

We present observations of circumstellar millimetre-wave 13CO line emission towards a sample of 20 cool carbon stars. Using a detailed radiative transfer model we estimate the circumstellar 12CO/13CO-ratios, which we believe accurately measure the important stellar 12C/13C-ratios. For those optically bright carbon stars where it is possible, our derived 12C/13C-ratios are compared with the photospheric results, obtained with different methods. Our estimates agree well with those of Lambert et al. (1986). It is shown that a straightforward determination of the 12CO/13CO-ratio from observed line intensity ratios is often hampered by optical depth effects, and that a detailed radiative transfer analysis is needed in order to determine reliable isotope ratios.

We have modelled radio and far-infrared rotational lines of 12^CO and 13^CO obtained from the circumstellar envelope of the infrared-bright carbon star IRAS 15194-5115. Eleven rotational lines between J=1-0 and 21-20 and nine rotational lines between J=1-0 and 22-21 in the ground vibrational states of 12^CO and 13^CO, respectively, provide the observational constraints. A model of the circumstellar envelope with a constant mass-loss rate (1 10^-5Msun/yr) and 12^CO/13^CO-ratio (5.5) is consistent with our observed FIR and radio data, suggesting that the wind characteristics have not changed significantly over the past few thousand years. Thus, IRAS 15194-5115 appears to be a highly evolved AGB-star, but the carbon star properties combined with the inferred low 12^C/13^C-ratio make the evolutionary status of this star uncertain. It may have been a J-star for which the 12^C}/13^C-ratio has remained low, or it may be a star of 5 to 8 solar masses, which has recently become a carbon star due to quenching of hot bottom burning. The dust properties or the relative amount of dust in the envelope appear to be different from those in the envelope of the well known carbon star IRC+10 216.