There is currently a severe discrepancy between theoretical models of dust formation in core-collapse supernovae (CCSNe), which predict $\gtrsim 0.01$ M$_\odot$ of ejecta dust forming within $\sim 1000$ days, and observations at these epochs, which infer much lower masses. We demonstrate that, in the optically thin case, these low dust masses are robust despite significant observational and model uncertainties. For a sample of 11 well-observed CCSNe, no plausible model reaches carbon dust masses above $10^{-4}$ M$_\odot$, or silicate masses above $\sim 10^{-3}$ M$_\odot$. Optically thick models can accommodate larger dust masses, but the dust must be clumped and have a low ($<0.1$) covering fraction to avoid conflict with data at optical wavelengths. These values are insufficient to reproduce the observed infrared fluxes, and the required covering fraction varies not only between SNe but between epochs for the same object. The difficulty in reconciling large dust masses with early-time observations of CCSNe, combined with well-established detections of comparably large dust masses in supernova remnants, suggests that a mechanism for late-time dust formation is necessary.

The luminous Type IIn SN 2010jl shows strong signs of interaction between the SN ejecta and dense circumstellar material. Dust may be present in the unshocked ejecta, the cool, dense shell between the shocks in the interaction region, or in the circumstellar medium. We present and model new optical and infrared photometry and spectroscopy of SN 2010jl from 82 to 1367 days since explosion. We evaluate the photometric and spectroscopic evolution using the radiative transfer codes MOCASSIN and DAMOCLES, respectively. We propose an interaction scenario and investigate the resulting dust formation scenarios and dust masses. We find that SN 2010jl has been continuously forming dust based on the evolution of its infrared emission and optical spectra. There is evidence for pre-existing dust in the circumstellar medium as well as new dust formation in the cool, dense shell and/or ejecta. We estimate that 0.005--0.01M$_{\odot}$ of predominantly carbon dust grains has formed in SN 2010jl by ~1400d post-outburst.

Core-collapse supernovae can condense large masses of dust post-explosion. However, sputtering and grain-grain collisions during the subsequent passage of the dust through the reverse shock can potentially destroy a significant fraction of the newly formed dust before it can reach the interstellar medium. Here we show that in oxygen-rich supernova remnants like Cassiopeia A the penetration and trapping within silicate grains of the same impinging ions of oxygen, silicon and magnesium that are responsible for grain surface sputtering can significantly reduce the net loss of grain material. We model conditions representative of dusty clumps (density contrast $\chi=100$) passing through the reverse shock in the oxygen-rich Cassiopeia A remnant and find that, compared to cases where the effect is neglected, as well as facilitating the formation of grains larger than those that had originally condensed, ion trapping increases the surviving masses of silicate dust by factors of up to two to four, depending on initial grain radii. For higher density contrasts ($\chi\gtrsim180$), we find that the effect of gas accretion on the surface of dust grains surpasses ion trapping, and the survival rate increases to ${\sim}55 \%$ of the initial dust mass for $\chi=256$.

We search for far-infrared (FIR) counterparts of known supernova remnants (SNRs) in the Galactic plane (360 degrees in longitude and b = +/- 1 deg ) at 70 - 500 micron with Herschel. We detect dust signatures in 39 SNRs out of 190, made up of 13 core-collapse supernovae (CCSNe), including 4 Pulsar Wind Nebulae (PWNe), and 2 Type Ia SNe. A further 24 FIR detected SNRs have unknown types. We confirm the FIR detection of ejecta dust within G350.1-0.3, adding to the known sample of ~10 SNRs containing ejecta dust. We discover dust features at the location of a radio core at the centre of G351.2+0.1, indicating FIR emission coincident with a possible Crab-like compact object, with dust temperature and mass of Td = 45.8 K and Md = 0.18 solar mass, similar to the PWN G54.1+0.3. We show that the detection rate is higher among young SNRs. We produce dust temperature maps of 11 SNRs and mass maps of those with distance estimates, finding dust at temperatures 15 < Td < 40 K. If the dust is heated by shock interactions the shocked gas must be relatively cool and/or have a low density to explain the observed low grain temperatures.

We calculate dust spectral energy distributions (SEDs) for a range of grain sizes and compositions, using physical properties appropriate for five pulsar wind nebulae (PWNe) from which dust emission associated with the ejecta has been detected. By fitting the observed dust SED with our models, with the number of grains of different sizes as the free parameters, we are able to determine the grain size distribution and total dust mass in each PWN. We find that all five PWNe require large ($\ge 0.1 \, {\rm \mu m}$) grains to make up the majority of the dust mass, with strong evidence for the presence of micron-sized or larger grains. Only two PWNe contain non-negligible quantities of small ($<0.01 \, {\rm \mu m}$) grains. The size distributions are generally well-represented by broken power laws, although our uncertainties are too large to rule out alternative shapes. We find a total dust mass of $0.02-0.28 \, {\rm M}_\odot$ for the Crab Nebula, depending on the composition and distance from the synchrotron source, in agreement with recent estimates. For three objects in our sample, the PWN synchrotron luminosity is insufficient to power the observed dust emission, and additional collisional heating is required, either from warm, dense gas as found in the Crab Nebula, or higher temperature shocked material. For G$54.1$+$0.3$, the dust is heated by nearby OB stars rather than the PWN. Inferred dust masses vary significantly depending on the details of the assumed heating mechanism, but in all cases large mass fractions of micron-sized grains are required.

We present high angular resolution (~80 mas) ALMA continuum images of the SN 1987A system, together with CO $J$=2 $\!\rightarrow\!$ 1, $J$=6 $\!\rightarrow\!$ 5, and SiO $J$=5 $\!\rightarrow\!$ 4 to $J$=7 $\!\rightarrow\!$ 6 images, which clearly resolve the ejecta (dust continuum and molecules) and ring (synchrotron continuum) components. Dust in the ejecta is asymmetric and clumpy, and overall the dust fills the spatial void seen in H$\alpha$ images, filling that region with material from heavier elements. The dust clumps generally fill the space where CO $J$=6 $\!\rightarrow\!$ 5 is fainter, tentatively indicating that these dust clumps and CO are locationally and chemically linked. In these regions, carbonaceous dust grains might have formed after dissociation of CO. The dust grains would have cooled by radiation, and subsequent collisions of grains with gas would also cool the gas, suppressing the CO $J$=6 $\!\rightarrow\!$ 5 intensity. The data show a dust peak spatially coincident with the molecular hole seen in previous ALMA CO $J$=2 $\!\rightarrow\!$ 1 and SiO $J$=5 $\!\rightarrow\!$ 4 images. That dust peak, combined with CO and SiO line spectra, suggests that the dust and gas could be at higher temperatures than the surrounding material, though higher density cannot be totally excluded. One of the possibilities is that a compact source provides additional heat at that location. Fits to the far-infrared--millimeter spectral energy distribution give ejecta dust temperatures of 18--23K. We revise the ejecta dust mass to $\mathrm{M_{dust}} = 0.2-0.4$M$_\odot$ for carbon or silicate grains, or a maximum of $<0.7$M$_\odot$ for a mixture of grain species, using the predicted nucleosynthesis yields as an upper limit.

Observations of the ejecta of core-collapse supernovae have shown that dust grains form in over-dense gas clumps in the expanding ejecta. The clumps are later subject to the passage of the reverse shock and a significant amount of the newly formed dust material can be destroyed due to the high temperatures and high velocities in the post-shock gas. To determine dust survival rates, we have performed a set of hydrodynamic simulations using the grid-based code AstroBEAR in order to model a shock wave interacting with a clump of gas and dust. Afterwards, dust motions and dust destruction rates are computed using our newly developed external, post-processing code Paperboats, which includes gas and plasma drag, grain charging, kinematic and thermal sputtering as well as grain-grain collisions. We have determined dust survival rates for the oxygen-rich supernova remnant Cassiopeia A as a function of initial grain sizes, dust materials and clump gas densities.

The reverse shock in the ejecta of core-collapse supernovae is potentially able to destroy newly formed dust material. In order to determine dust survival rates, we have performed a set of hydrodynamic simulations using the grid-based code AstroBEAR in order to model a shock wave interacting with clumpy supernova ejecta. Dust motions and destruction rates were computed using our newly developed external, post-processing code Paperboats, which includes gas drag, grain charging, sputtering and grain-grain collisions. We have determined dust destruction rates for the oxygen-rich supernova remnant Cassiopeia A as a function of initial grain sizes and clump gas density. We found that up to 30 % of the carbon dust mass is able to survive the passage of the reverse shock if the initial grain size distribution is narrow with radii around ~10 - 50 nm for high gas densities, or with radii around ~0.5 - 1.5 ${\mu}$m for low and medium gas densities. Silicate grains with initial radii around 10 - 30 nm show survival rates of up to 40 % for medium and high density contrasts, while silicate material with micron sized distributions is mostly destroyed. For both materials, the surviving dust mass is rearranged into a new size distribution that can be approximated by two components: a power-law distribution of small grains and a log-normal distribution of grains having the same size range as the initial distribution. Our results show that grain-grain collisions and sputtering are synergistic and that grain-grain collisions can play a crucial role in determining the surviving dust budget in supernova remnants.

Eta Car is one of the most luminous and massive stars in our Galaxy and is the brightest mid-infrared (mid-IR) source in the sky, outside our solar system. Since the late 1990s the central source has dramatically brightened at ultraviolet and optical wavelengths. This might be explained by a decrease in circumstellar dust extinction. We aim to establish the mid-IR flux evolution and further our understanding of the star's ultraviolet and optical brightening. Mid-IR images from $8-20~\mu$m were obtained in 2018 with VISIR at the Very Large Telescope. Archival data from 2003 and 2005 are retrieved from the ESO Science Archive Facility and historical records are collected from publications. We present the highest angular resolution mid-IR images of $\eta$ Car to date at the corresponding wavelengths ($\geq 0.22''$). We reconstruct the mid-IR evolution of the spectral energy distribution of the spatially integrated Homunculus nebula from 1968 to 2018 and find no long-term changes. Eta Car's bolometric luminosity has been stable over the past five decades. We do not observe a long-term decrease in the mid-IR flux densities that could be associated with the brightening at ultraviolet and optical wavelengths, but circumstellar dust must be declining in our line-of-sight only. Short-term flux variations within about 25% of the mean levels could be present.

We have modelled the near-infrared to radio images of the Crab Nebula with a Bayesian SED model to simultaneously fit its synchrotron, interstellar and supernova dust emission. We infer an interstellar dust extinction map with an average $A_{\text{V}}$=1.08$\pm$0.38 mag, consistent with a small contribution (<22%) to the Crab's overall infrared emission. The Crab's supernova dust mass is estimated to be between 0.032 and 0.049 M$_{\odot}$ (for amorphous carbon grains) with an average dust temperature $T_{\text{dust}}$=41$\pm$3K, corresponding to a dust condensation efficiency of 8-12%. This revised dust mass is up to an order of magnitude lower than some previous estimates, which can be attributed to our different interstellar dust corrections, lower SPIRE flux densities, and higher dust temperature than were used in previous studies. The dust within the Crab is predominantly found in dense filaments south of the pulsar, with an average V-band dust extinction of $A_{\text{V}}$=0.20-0.39 mag, consistent with recent optical dust extinction studies. The modelled synchrotron power-law spectrum is consistent with a radio spectral index $\alpha_{\text{radio}}$=0.297$\pm$0.009 and an infrared spectral index $\alpha_{\text{IR}}$=0.429$\pm$0.021. We have identified a millimetre excess emission in the Crab's central regions, and argue that it most likely results from two distinct populations of synchrotron emitting particles. We conclude that the Crab's efficient dust condensation (8-12%) provides further evidence for a scenario where supernovae can provide substantial contributions to the interstellar dust budgets in galaxies.

In order to understand the contribution of core-collapse supernovae to the dust budget of the early universe, it is important to understand not only the mass of dust that can form in core-collapse supernovae but also the location and rate of dust formation. SN 2005ip is of particular interest since dust has been inferred to have formed in both the ejecta and the post-shock region behind the radiative reverse shock. We have collated eight optical archival spectra that span the lifetime of SN 2005ip and we additionally present a new X-shooter optical-near-IR spectrum of SN 2005ip at 4075d post-discovery. Using the Monte Carlo line transfer code DAMOCLES, we have modelled the blueshifted broad and intermediate width H$\alpha$, H$\beta$ and He I lines from 48d to 4075d post-discovery using an ejecta dust model. We find that dust in the ejecta can account for the asymmetries observed in the broad and intermediate width H$\alpha$, H$\beta$ and He I line profiles at all epochs and that it is not necessary to invoke post-shock dust formation to explain the blueshifting observed in the intermediate width post-shock lines. Using a Bayesian approach, we have determined the evolution of the ejecta dust mass in SN 2005ip over 10 years presuming an ejecta dust model, with an increasing dust mass from ~10$^{-8}$ M$_{\odot}$ at 48d to a current dust mass of $\sim$0.1 M$_{\odot}$.

We model the thermal dust emission from dust grains heated by synchrotron radiation and by particle collisions, under conditions appropriate for four different shocked and unshocked gas components of the Cassiopeia A (Cas A) supernova remnant (SNR). By fitting the resulting spectral energy distributions (SEDs) to the observed SNR dust fluxes, we determine the required mass of dust in each component. We find the observed SED can be reproduced by $\sim 0.6 \, {\rm M_\odot}$ of silicate grains, the majority of which is in the unshocked ejecta and heated by the synchrotron radiation field. Warmer dust, located in the X-ray emitting reverse shock and blastwave regions, contribute to the shorter wavelength infrared emission but make only a small fraction of the total dust mass. Carbon grains can at most make up $\sim 25 \%$ of the total dust mass. Combined with estimates for the gas masses, we obtain dust-to-gas mass ratios for each component, which suggest that the condensation efficiency in the ejecta is high, and that dust in the shocked ejecta clumps is well protected from destruction by sputtering in the reverse shock.

We search for far-infrared (FIR) counterparts of known supernova remnants (SNRs) in the Galactic plane (10 deg <| l |< 60 deg) at 70-500 micron using the Herschel Infrared Galactic Plane Survey (Hi-GAL). Of 71 sources studied, we find that 29 (41 per cent) SNRs have a clear FIR detection of dust emission associated with the SNR. Dust from 8 of these is in the central region, and 4 indicate pulsar wind nebulae (PWNe) heated ejecta dust. A further 23 have dust emission in the outer shell structures which is potentially related to swept up material. Many Galactic SNe have dust signatures but we are biased towards detecting ejecta dust in young remnants and those with a heating source (shock or PWN). We estimate the dust temperature and mass contained within three PWNe, G11.2-0.3, G21.5-0.9, and G29.7-0.3 using modified blackbody fits. To more rigorously analyse the dust properties at various temperatures and dust emissivity index beta, we use point process mapping (PPMAP). We find significant quantities of cool dust (at 20-40 K) with dust masses of Md = 0.34 +/- 0.14 solar mass, Md = 0.29 +/- 0.08 solar mass, and Md = 0.51 +/- 0.13 solar mass for G11.2-0.3, G21.5-0.9, and G29.7-0.3 respectively. We derive the dust emissivity index for the PWN ejecta dust in G21.5-0.3 to be beta = 1.4 +/- 0.5 compared to dust in the surrounding medium where beta = 1.8 +/- 0.1.

The spatial structure of the emission lines and continuum over the 50 arcsecond extent of the nearby, O-rich, planetary nebula NGC 7009 (Saturn Nebula) have been observed with the MUSE integral field spectrograph on the ESO Very Large Telescope. Science Verification data, in <0.6 arcsecond seeing, have been reduced and analysed as images over the wavelength range 4750-9350A. Emission line maps over the bright shells are presented, from neutral to the highest ionization available (He II and [Mn V]). For collisionally excited lines (CELs), maps of electron temperature (T_e from [N II] and [S III]) and electron density (N_e from [S II] and [Cl III]) are available and for optical recombination lines (ORLs) temperature (from the Paschen jump and ratio of He I lines) and density (from high Paschen lines). These estimates are compared: for the first time, maps of the differences in CEL and ORL T_e's have been derived, and correspondingly a map of t^2 between a CEL and ORL temperature, showing considerable detail. Total abundances of He and O were formed, the latter using three ionization correction factors. However the map of He/H is not flat, departing by ~2% from a constant value, with remnants corresponding to ionization structures. Ionization correction factor methods are compared for O abundance, but none delivers a flat map. An integrated spectrum over an area of 2340 square arcseconds was also formed and compared to 1D photoionization models. The spatial variation of a range of nebular parameters illustrates the complexity of the ionized media in NGC 7009. These MUSE data are very rich with detections of many lines over areas of hundreds of square arcseconds and follow-on studies are indicated. (Abridged)

The equatorial ring of Supernova (SN) 1987A has been exposed to forward shocks from the SN blast wave, and it has been suggested that these forward shocks have been causing on-going destruction of dust in the ring. We obtained SOFIA FORCAST 11.1, 19.7 and 31.5 micron photometry of SN\,1987A in 2016. Compared with Spitzer measurements 10 years earlier, the 31.5 micron flux has significantly increased. The excess at 31.5 micron appears to be related to the Herschel 70 micron excess, which was detected 5 years earlier. The dust mass needed to account for the the 31.5--70 micron excess is 3--7x10^-4 Msun, more than ten times larger than the ring dust mass (1x10^-5 Msun) estimate from the data 10-years earlier. We argue that dust grains are re-formed or grown in the post-shock regions in the ring after forward shocks have destroyed pre-existing dust grains in the ring and released refractory elements into gas. In the post-shock region, atoms can stick to surviving dust grains, and the dust mass may have increased (grain growth), or dust grains might have condensed directly from the gas. An alternative possibility is that the outer part of the expanding ejecta dust might have been heated by X-ray emission from the circumstellar ring. The future development of this excess could reveal whether grains are reformed in the post-shocked region of the ring or eject dust is heated by X-ray.

We model the molecular emission from cometary knots in planetary nebulae (PNe) using a combination of photoionization and photodissociation region (PDR) codes, for a range of central star properties and gas densities. Without the inclusion of ionizing extreme ultraviolet (EUV) radiation, our models require central star temperatures $T_*$ to be near the upper limit of the range investigated in order to match observed H$_2$ and OH$^+$ surface brightnesses consistent with observations - with the addition of EUV flux, our models reproduce observed OH$^+$ surface brightnesses for $T_* \ge 100 \, {\rm kK}$. For $T_* < 80 \, {\rm kK}$, the predicted OH$^+$ surface brightness is much lower, consistent with the non-detection of this molecule in PNe with such central star temperatures. Our predicted level of H$_2$ emission is somewhat weaker than commonly observed in PNe, which may be resolved by the inclusion of shock heating or fluorescence due to UV photons. Some of our models also predict ArH$^+$ and HeH$^+$ rotational line emission above detection thresholds, despite neither molecule having been detected in PNe, although the inclusion of photodissociation by EUV photons, which is neglected by our models, would be expected to reduce their detectability.

We have performed combined photoionization and photodissociation region (PDR) modelling of a Crab Nebula filament subjected to the synchrotron radiation from the central pulsar wind nebula, and to a high flux of charged particles; a greatly enhanced cosmic ray ionization rate over the standard interstellar value, $\zeta_0$, is required to account for the lack of detected [C I] emission in published Herschel SPIRE FTS observations of the Crab Nebula. The observed line surface brightness ratios of the OH$^+$ and ArH$^+$ transitions seen in the SPIRE FTS frequency range can only be explained with both a high cosmic ray ionization rate and a reduced ArH$^+$ dissociative recombination rate compared to that used by previous authors, although consistent with experimental upper limits. We find that the ArH$^+$/OH$^+$ line strengths and the observed H$_2$ vibration-rotation emission can be reproduced by model filaments with $n_{\rm{H}} = 2 \times 10^4$ cm$^{-3}$, $\zeta = 10^7 \zeta_0$ and visual extinctions within the range found for dusty globules in the Crab Nebula, although far-infrared emission from [O I] and [C II] is higher than the observational constraints. Models with $n_{\rm{H}} = 1900$ cm$^{-3}$ underpredict the H$_2$ surface brightness, but agree with the ArH$^+$ and OH$^+$ surface brightnesses and predict [O I] and [C II] line ratios consistent with observations. These models predict HeH$^+$ rotational emission above detection thresholds, but consideration of the formation timescale suggests that the abundance of this molecule in the Crab Nebula should be lower than the equilibrium values obtained in our analysis.

Infrared observations of the dusty, massive Homunculus Nebula around the luminous blue variable $\eta$ Carinae are crucial to characterize the mass-loss history and help constrain the mechanisms leading to the Great Eruption. We present the 2.4 - 670 $\mu$m spectral energy distribution, constructed from legacy ISO observations and new spectroscopy obtained with the {\em{Herschel Space Observatory}}. Using radiative transfer modeling, we find that the two best-fit dust models yield compositions which are consistent with CNO-processed material, with iron, pyroxene and other metal-rich silicates, corundum, and magnesium-iron sulfide in common. Spherical corundum grains are supported by the good match to a narrow 20.2 $\mu$m feature. Our preferred model contains nitrides AlN and Si$_3$N$_4$ in low abundances. Dust masses range from 0.25 to 0.44 $M_\odot$ but $M_{\rm{tot}} \ge$ 45 $M_\odot$ in both cases due to an expected high Fe gas-to-dust ratio. The bulk of dust is within a 5$"$ $\times$ 7$"$ central region. An additional compact feature is detected at 390 $\mu$m. We obtain $L_{\rm{IR}}$ = 2.96 $\times$ 10$^6$ $L_\odot$, a 25\% decline from an average of mid-IR photometric levels observed in 1971-1977. This indicates a reduction in circumstellar extinction in conjunction with an increase in visual brightness, allowing 25-40\% of optical and UV radiation to escape from the central source. We also present an analysis of $^{12}$CO and $^{13}$CO $J = 5-4$ through $9-8$ lines, showing that the abundances are consistent with expectations for CNO-processed material. The [$^{12}$C~{\sc{ii}}] line is detected in absorption, which we suspect originates in foreground material at very low excitation temperatures.

We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta morphology is likely elongated. The difference of the CO and SiO line profiles is consistent with hydrodynamic simulations, which show that Rayleigh-Taylor instabilities cause mixing of gas, with heavier elements much more disturbed, making more elongated structure. We obtained isotopologue ratios of 28SiO/29SiO>13, 28SiO/30SiO>14, and 12CO/13CO>21, with the most likely limits of 28SiO/29SiO>128, 28SiO/30SiO>189. Low 29Si and 30Si abundances in SN 1987A are consistent with nucleosynthesis models that show inefficient formation of neutron-rich isotopes in a low metallicity environment, such as the Large Magellanic Cloud. The deduced large mass of HCO+ (~5x10^-6 Msun) and small SiS mass (<6x10^-5 Msun) might be explained by some mixing of elements immediately after the explosion. The mixing might have caused some hydrogen from the envelope to sink into carbon and oxygen-rich zones after the explosion, enabling the formation of a substantial mass of HCO+. Oxygen atoms may have penetrated into silicon and sulphur zones, suppressing formation of SiS. Our ALMA observations open up a new window to investigate chemistry, dynamics and explosive-nucleosynthesis in supernovae.

Theoretical models predict that core-collapse supernovae (CCSNe) can be efficient dust producers (0.1-1.0 Msun), potentially accounting for most of the dust production in the early Universe. Observational evidence for this dust production efficiency is however currently limited to only a few CCSN remnants (e.g., SN1987A, Crab Nebula). In this paper, we revisit the dust mass produced in Cassiopeia A (Cas A), a ~330-year old O-rich Galactic supernova remnant (SNR) embedded in a dense interstellar foreground and background. We present the first spatially resolved analysis of Cas A based on Spitzer and Herschel infrared and submillimetre data at a common resolution of ~0.6 arcmin for this 5 arcmin diameter remnant following a careful removal of contaminating line emission and synchrotron radiation. We fit the dust continuum from 17 to 500 micron with a four-component interstellar medium (ISM) and supernova (SN) dust model. We find a concentration of cold dust in the unshocked ejecta of Cas A and derive a mass of 0.3-0.5 Msun of silicate grains freshly produced in the SNR, with a lower limit of >=0.1-0.2 Msun. For a mixture of 50% of silicate-type grains and 50% of carbonaceous grains, we derive a total SN dust mass between 0.4 Msun and 0.6 Msun. These dust mass estimates are higher than from most previous studies of Cas A and support the scenario of supernova dominated dust production at high redshifts. We furthermore derive an interstellar extinction map for the field around Cas A which towards Cas A gives average values of A_V=6-8 mag, up to a maximum of A_V=15 mag.

We have conducted a search for optical circumstellar absorption lines in the spectra of 16 debris disc host stars. None of the stars in our sample showed signs of emission line activity in either H$_{\alpha}$, Ca II or Na I, confirming their more evolved nature. Four stars were found to exhibit narrow absorption features near the cores of the photospheric Ca II and Na I D lines (when Na I D data were available). We analyse the characteristics of these spectral features to determine whether they are of circumstellar or interstellar origins. The strongest evidence for circumstellar gas is seen in the spectrum of HD110058, which is known to host a debris disc observed close to edge-on. This is consistent with a recent ALMA detection of molecular gas in this debris disc, which shows many similarities to the $\beta$ Pictoris system.

We present Monte Carlo line transfer models that investigate the effects of dust on the very late time emission line spectra of the core collapse supernovae SN 1980K and SN 1993J and the young core collapse supernova remnant Cassiopeia A. Their blue-shifted emission peaks, resulting from the removal by dust of redshifted photons emitted from the far sides of the remnants, and the presence of extended red emission wings are used to constrain dust compositions and radii and to determine the masses of dust in the remnants. We estimate dust masses of between 0.08 - 0.15 M$_\odot$ for SN 1993J at year 16, 0.12 - 0.30 M$_\odot$ for SN 1980K at year 30 and $\sim$1.1 M$_\odot$ for Cas A at year $\sim$330. Our models for the strong oxygen forbidden lines of Cas A require the overall modelled profiles to be shifted to the red by between 700 - 1000 km s$^{-1}$, consistent with previous estimates for the shift of the dynamical centroid of this remnant.

Using the PACS and SPIRE spectrometers on-board the Herschel Space Observatory, we obtained spectra of two red supergiants (RSGs) in the Large Magellanic Cloud (LMC). Multiple rotational CO emission lines (J=6-5 to 15-14) and 15 H2O lines were detected from IRAS 05280-6910, and one CO line was detected from WOH G64. This is the first time CO rotational lines have been detected from evolved stars in the LMC. Their CO line intensities are as strong as those of the Galactic RSG, VY CMa. Modelling the CO lines and the spectral energy distribution results in an estimated mass-loss rate for IRAS 05280-6910 of 3x10^-4 Msun per yr. The model assumes a gas-to-dust ratio and a CO-to-H2 abundance ratio is estimated from the Galactic values scaled by the LMC metallicity ([Fe/H]~-0.3), i.e., that the CO-to-dust ratio is constant for Galactic and LMC metallicities within the uncertainties of the model. The key factor determining the CO line intensities and the mass-loss rate found to be the stellar luminosity.

The young star beta Pictoris is well known for its dusty debris disk, produced through the grinding down by collisions of planetesimals, kilometre-sized bodies in orbit around the star. In addition to dust, small amounts of gas are also known to orbit the star, likely the result from vaporisation of violently colliding dust grains. The disk is seen edge on and from previous absorption spectroscopy we know that the gas is very rich in carbon relative to other elements. The oxygen content has been more difficult to assess, however, with early estimates finding very little oxygen in the gas at a C/O ratio 20x higher than the cosmic value. A C/O ratio that high is difficult to explain and would have far-reaching consequences for planet formation. Here we report on observations by the far-infrared space telescope Herschel, using PACS, of emission lines from ionised carbon and neutral oxygen. The detected emission from C+ is consistent with that previously reported being observed by the HIFI instrument on Herschel, while the emission from O is hard to explain without assuming a higher-density region in the disk, perhaps in the shape of a clump or a dense torus, required to sufficiently excite the O atoms. A possible scenario is that the C/O gas is produced by the same process responsible for the CO clump recently observed by ALMA in the disk, and that the re-distribution of the gas takes longer than previously assumed. A more detailed estimate of the C/O ratio and the mass of O will have to await better constraints on the C/O gas spatial distribution.

(Abridged) We present the first public release of high-quality data products (DR1) from Hi-GAL, the {\em Herschel} infrared Galactic Plane Survey. Hi-GAL is the keystone of a suite of continuum Galactic Plane surveys from the near-IR to the radio, and covers five wavebands at 70, 160, 250, 350 and 500 micron, encompassing the peak of the spectral energy distribution of cold dust for 8 < T < 50K. This first Hi-GAL data release covers the inner Milky Way in the longitude range 68{\deg} > l > -70{\deg} in a |b|<1{\deg} latitude strip. Photometric maps have been produced with the ROMAGAL pipeline, that optimally capitalizes on the excellent sensitivity and stability of the bolometer arrays of the {\em Herschel} PACS and SPIRE photometric cameras, to deliver images of exquisite quality and dynamical range, absolutely calibrated with {\em Planck} and {\em IRAS}, and recovering extended emission at all wavelengths and all spatial scales. The compact source catalogues have been generated with the CuTEx algorithm, specifically developed to optimize source detection and extraction in the extreme conditions of intense and spatially varying background that are found in the Galactic Plane in the thermal infrared. Hi-GAL DR1 images will be accessible via a dedicated web-based image cutout service. The DR1 Compact Source Catalogues are delivered as single-band photometric lists containing, in addition to source position, peak and integrated flux and source sizes, a variety of parameters useful to assess the quality and reliability of the extracted sources, caveats and hints to help this assessment are provided. Flux completeness limits in all bands are determined from extensive synthetic source experiments and depend on the specific line of sight along the Galactic Plane. Hi-GAL DR1 catalogues contain 123210, 308509, 280685, 160972 and 85460 compact sources in the five bands, respectively.

The large field and wavelength range of MUSE is well suited to mapping Galactic planetary nebulae (PN). The bright PN NGC 7009 was observed with MUSE on the VLT during the Science Verification of the instrument in seeing of 0.6". Emission line maps in hydrogen Balmer and Paschen lines were formed from analysis of the MUSE cubes. The measured electron temperature and density from the MUSE cube were employed to predict the theoretical hydrogen line ratios and map the extinction distribution across the nebula. After correction for the interstellar extinction to NGC 7009, the internal dust-to-gas ratio (A_V/N_H) has been mapped for the first time in a PN. The extinction map of NGC 7009 has considerable structure, broadly corresponding to the morphological features of the nebula. A large-scale feature in the extinction map, consisting of a crest and trough, occurs at the rim of the inner shell. The nature of this feature was investigated and instrumental and physical causes considered; no convincing mechanisms were identified to produce this feature, other than mass loss variations in the earlier asymptotic giant branch phase. The dust-to-gas ratio A_V/N_H increases from 0.7 times the interstellar value to >5 times from the centre towards the periphery of the ionized nebula. The integrated A_V/N_H is about 2 times the mean ISM value. It is demonstrated that extinction mapping with MUSE provides a powerful tool for studying the distribution of PN internal dust and the dust-to-gas ratio. (Abridged.)

Infrared Herschel imaging and spectroscopic observations of the nebula M1-67 around the Wolf-Rayet star WR 124 have been obtained along with optical imaging observations. The infrared images reveal a clumpy dusty nebula that extends up to 1 pc. The comparison with the optical images shows that the ionized gas nebula coincides with the dust nebula, the dust and the gas being mixed together. A photodissociation region is revealed from the infrared spectroscopic analysis. The analysis of the infrared spectrum of the nebula, where forbidden emission lines of ionized elements were detected, showed that the nebula consists of mildly processed material with the calculated abundance number ratios being N/O = 1.0 +/- 0.5 and C/O = 0.46 +/- 0.27. Based on a radiative transfer model, the dust mass of the nebula was estimated to be 0.22 Msun with a population of large grains being necessary to reproduce the observations. The comparison of the mass-loss rate and the abundance ratios to theoretical models of stellar evolution led to the conclusion that the nebular ejection took place during a RSG/YSG evolutionary phase of a central star with an initial mass of 32 Msun.

The late time optical and near-IR line profiles of many core-collapse supernovae exhibit a red-blue asymmetry as a result of greater extinction by internal dust of radiation emitted from the receding parts of the supernova ejecta. We present here a new code, DAMOCLES, that models the effects of dust on the line profiles of core-collapse supernovae in order to determine the masses of newly formed dust. As noted by Lucy et al. (1989), the presence of an extended red scattering wing in late-time line profiles can also indicate dust formation. We find that dust-affected line profiles need not necessarily be flux-biased towards to the blue, although the profile peak will always be blue-shifted. We have collated optical spectra of SN 1987A from a variety of archival sources and have modelled the evolution of the H$\alpha$ line from days 714 to 3604, as well as that of the [OI] 6300,6363A doublet between days 714 and 1478. A variety of evidence points to the presence of clumping and we find that our clumped dust models require significantly higher dust masses than smoothly distributed dust models. Our line profile fits imply day 714 dust masses of <3 $\times$ 10$^{-3}$ M$_{\odot}$ for all grain types apart from very high albedo pure magnesium silicates, for which up to 0.07M$_{\odot}$ can be accommodated. Large grain radii (>0.6$\mu$m) are generally required to fit the line profiles even at the earlier epochs. We find that a large dust mass (>0.1M$_{\odot}$) had formed by day 3604 and infer that the majority of the present dust mass must have formed after this epoch. Our findings agree with recent estimates from SED fits for the dust mass evolution of SN 1987A and support the inference that the majority of SN 1987A's dust formed many years after the initial explosion.

SN 2011ja was a bright (I = -18.3) Type II supernova occurring in the nearby edge on spiral galaxy NGC 4945. Flat-topped and multi-peaked H-alpha and H-beta spectral emission lines appear between 64 - 84 days post-explosion, indicating interaction with a disc-like circumstellar medium inclined 30-45 degrees from edge-on. After day 84 an increase in the H- and K-band flux along with heavy attenuation of the red wing of the emission lines are strong indications of early dust formation, likely located in the cool dense shell created between the forward shock of the SN ejecta and the reverse shock created as the ejecta plows into the existing CSM. Radiative transfer modeling reveals both ~1.5 x 10^-4 Msun of pre-existing dust located ~ 10^16.7 cm away and ~ 5 x 10^-5 Msun of newly formed dust. Spectral observations after 1.5 years reveal the possibility that the fading SN is located within a young (3-6 Myr) massive stellar cluster, which when combined with tentative 56Ni mass estimates of 0.2 Msun may indicate a massive (> 25 Msun) progenitor for SN 2011ja.

The interaction of ionizing and far-ultraviolet radiation with the interstellar medium is of great importance. It results in the formation of regions in which the gas is ionized, beyond which are photodissociation regions (PDRs) in which the gas transitions to its atomic and molecular form. Several numerical codes have been implemented to study these two main phases of the interstellar medium either dynamically or chemically. In this paper we present TORUS-3DPDR, a new self-consistent code for treating the chemistry of three-dimensional photoionization and photodissociation regions. It is an integrated code coupling the two codes TORUS, a hydrodynamics and Monte Carlo radiation transport code, and 3D-PDR, a photodissociation regions code. The new code uses a Monte Carlo radiative transfer scheme to account for the propagation of the ionizing radiation including the diffusive component as well as a ray-tracing scheme based on the HEALPix package in order to account for the escape probability and column density calculations. Here, we present the numerical techniques we followed and we show the capabilities of the new code in modelling three-dimensional objects including single or multiple sources. We discuss the effects introduced by the diffusive component of the UV field in determining the thermal balance of PDRs as well as the effects introduced by a multiple sources treatment of the radiation field. We find that diffuse radiation can positively contribute to the formation of CO. With this new code, three-dimensional synthetic observations for the major cooling lines are possible, for making feasible a detailed comparison between hydrodynamical simulations and observations.

Aim: The late stages of stellar evolution are mainly governed by the mass of the stars. Low- and intermediate-mass stars lose copious amounts of mass during the asymptotic giant branch (AGB) which obscure the central star making it difficult to study the stellar spectra and determine the stellar mass. In this study, we present observational data that can be used to determine lower limits to the stellar mass. Method: Spectra of nine heavily reddened AGB stars taken by the Herschel Space Observatory display numerous molecular emission lines. The strongest emission lines are due to H2O. We search for the presence of isotopologues of H2O in these objects. Result: We detected the 16O and 17O isotopologues of water in these stars, but lines due to H2^{18}O are absent. The lack of 18O is predicted by a scenario where the star has undergone hot-bottom burning which preferentially destroys 18O relative to 16O and 17O. From stellar evolution calculations, this process is thought to occur when the stellar mass is above 5 Msun for solar metallicity. Hence, observations of different isotopologues of H2O can be used to help determine the lower limit to the initial stellar mass. Conclusion: From our observations, we deduce that these extreme OH/IR stars are intermediate-mass stars with masses of >= 5 Msun. Their high mass-loss rates of ~ 1.0e-4 Msun/yr may affect the enrichment of the interstellar medium and the overall chemical evolution of our Galaxy.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

We have constructed MOCASSIN photoionization plus dust radiative transfer models for the Crab Nebula core-collapse supernova (CCSN) remnant, using either smooth or clumped mass distributions, in order to determine the chemical composition and masses of the nebular gas and dust. We computed models for several different geometries suggested for the nebular matter distribution but found that the observed gas and dust spectra are relatively insensitive to these geometries, being determined mainly by the spectrum of the pulsar wind nebula which ionizes and heats the nebula. Smooth distribution models are ruled out since they require 16-49 Msun of gas to fit the integrated optical nebular line fluxes, whereas our clumped models re quire 7.0 Msun of gas. A global gas-phase C/O ratio of 1.65 by number is derived, along with a He/H number ratio of 1.85, neither of which can be matched by current CCSN yield predictions. A carbonaceous dust composition is favoured by the observed gas-phase C/O ratio: amorphous carbon clumped model fits to the Crab's Herschel and Spitzer infrared spectral energy distribution imply the presence of 0.18-0.27 Msun of dust, corresponding to a gas to dust mass ratio of 26-39. Mixed dust chemistry models can also be accommodated, comprising 0.11-0.13 Msun of amorphous carbon and 0.39-0.47 Msun of silicates. Power-law grain size distributions with mass distributions that are weighted towards the largest grain radii are derived, favouring their longer-term survival when they eventually interact with the interstellar medium. The total mass of gas plus dust in the Crab Nebula is 7.2 +/- 0.5 Msun, consistent with a progenitor star mass of 9 Msun.

In our series of papers presenting the Herschel imaging of evolved planetary nebulae, we present images of the dust distribution in the Helix nebula (NGC 7293). Images at 70, 160, 250, 350, and 500 micron were obtained with the PACS and SPIRE instruments on board the Herschel satellite. The broadband maps show the dust distribution over the main Helix nebula to be clumpy and predominantly present in the barrel wall. We determined the spectral energy distribution of the main nebula in a consistent way using Herschel, IRAS, and Planck flux values. The emissivity index of 0.99 +/- 0.09, in combination with the carbon rich molecular chemistry of the nebula, indicates that the dust consists mainly of amorphous carbon. The dust excess emission from the central star disk is detected at 70 micron and the flux measurement agree with previous measurement. We present the temperature and dust column density maps. The total dust mass across the Helix nebula (without its halo) is determined to be 0.0035 solar mass at a distance of 216 pc. The temperature map shows dust temperatures between 22 and 42 K, which is similar to the kinetic temperature of the molecular gas, strengthening the fact that the dust and gas co-exist in high density clumps. Archived images are used to compare the location of the dust emission in the far infrared (Herschel) with the ionized (GALEX, Hbeta) and molecular hydrogen component. The different emission components are consistent with the Helix consisting of a thick walled barrel-like structure inclined to the line of sight. The radiation field decreases rapidly through the barrel wall.

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 micron data and improved imaging quality at 100 and 160 micron compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 micron [O I] line flux, eliminating the possibility that line contaminations distort the previously estimated dust mass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 micron flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5+-0.1 Msun of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3 Msun of amorphous carbon and 0.5 Msun of silicates, totalling 0.8 Msun of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies.

We present a comprehensive spectral and morphological analysis of the remnant of Supernova (SN) 1987A with the Australia Telescope Compact Array (ATCA) and the Atacama Large Millimeter/submillimeter Array (ALMA). The non-thermal and thermal components of the radio emission are investigated in images from 94 to 672 GHz ($\lambda$ 3.2 mm to 450 $\mu$m), with the assistance of a high-resolution 44 GHz synchrotron template from the ATCA, and a dust template from ALMA observations at 672 GHz. An analysis of the emission distribution over the equatorial ring in images from 44 to 345 GHz highlights a gradual decrease of the east-to-west asymmetry ratio with frequency. We attribute this to the shorter synchrotron lifetime at high frequencies. Across the transition from radio to far infrared, both the synchrotron/dust-subtracted images and the spectral energy distribution (SED) suggest additional emission beside the main synchrotron component ($S_{\nu}\propto\nu^{-0.73}$) and the thermal component originating from dust grains at $T\sim22$ K. This excess could be due to free-free flux or emission from grains of colder dust. However, a second flat-spectrum synchrotron component appears to better fit the SED, implying that the emission could be attributed to a pulsar wind nebula (PWN). The residual emission is mainly localised west of the SN site, as the spectral analysis yields $-0.4\lesssim\alpha\lesssim-0.1$ across the western regions, with $\alpha\sim0$ around the central region. If there is a PWN in the remnant interior, these data suggest that the pulsar may be offset westward from the SN position.

The discovery with the {\it Herschel Space Observatory} of bright far infrared and submm emission from the ejecta of the core collapse supernova SN\,1987A has been interpreted as indicating the presence of some 0.4--0.7\,M$_\odot$ of dust. We have constructed radiative transfer models of the ejecta to fit optical to far-infrared observations from the literature at epochs between 615 days and 24 years after the explosion, to determine when and where this unexpectedly large amount of dust formed. We find that the observations by day 1153 are consistent with the presence of 3$\times$10$^{-3}$M$_\odot$ of dust. Although this is a larger amount than has previously been considered possible at this epoch, it is still very small compared to the amount present in the remnant after 24 years, and significantly higher dust masses at the earlier epochs are firmly ruled out by the observations, indicating that the majority of the dust must have formed at very late times. By 8515-9200 days after the explosion, 0.6--0.8\,M$_\odot$ of dust is present, and dust grains with radii greater than 2\,$\mu$m are required to obtain a fit to the observed SED. This suggests that the dust mass increase at late times was caused by accretion onto and coagulation of the dust grains formed at earlier epochs. These findings provide further confirmation that core collapse supernovae can create large quantities of dust, and indicate that the reason for small dust masses being estimated in many cases is that the vast majority of the dust forms long after most supernovae have been detectable at mid-infrared wavelengths.

We report on the discovery of strong intensity variations in the high rotational lines of abundant molecular species towards the archetypical circumstellar envelope of IRC+10216. The observations have been carried out with the HIFI instrument on board \textit{Herschel}\thanks{\textit{Herschel} is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA} and with the IRAM\thanks{This work was based on observations carried out with the IRAM 30-meter telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain)} 30-m telescope. They cover several observing periods spreading over 3 years. The line intensity variations for molecules produced in the external layers of the envelope most probably result from time variations in the infrared pumping rates. We analyze the main implications this discovery has on the interpretation of molecular line emission in the envelopes of Mira-type stars. Radiative transfer calculations have to take into account both the time variability of infrared pumping and the possible variation of the dust and gas temperatures with stellar phase in order to reproduce the observation of molecular lines at different epochs. The effect of gas temperature variations with stellar phase could be particularly important for lines produced in the innermost regions of the envelope. Each layer of the circumstellar envelope sees the stellar light radiation with a different lag time (phase). Our results show that this effect must be included in the models. The sub-mm and FIR lines of AGB stars cannot anymore be considered as safe intensity calibrators.

The INT/WFC Photometric H-Alpha Survey of the Northern Galactic Plane (IPHAS) is a 1800 square degrees imaging survey covering Galactic latitudes |b| < 5 deg and longitudes l = 30 to 215 deg in the r, i and H-alpha filters using the Wide Field Camera (WFC) on the 2.5-metre Isaac Newton Telescope (INT) in La Palma. We present the first quality-controlled and globally-calibrated source catalogue derived from the survey, providing single-epoch photometry for 219 million unique sources across 92% of the footprint. The observations were carried out between 2003 and 2012 at a median seeing of 1.1 arcsec (sampled at 0.33 arcsec/pixel) and to a mean 5\sigma-depth of 21.2 (r), 20.0 (i) and 20.3 (H-alpha) in the Vega magnitude system. We explain the data reduction and quality control procedures, describe and test the global re-calibration, and detail the construction of the new catalogue. We show that the new calibration is accurate to 0.03 mag (rms) and recommend a series of quality criteria to select the most reliable data from the catalogue. Finally, we demonstrate the ability of the catalogue's unique (r-Halpha, r-i) diagram to (1) characterise stellar populations and extinction regimes towards different Galactic sightlines and (2) select H-alpha emission-line objects. IPHAS is the first survey to offer comprehensive CCD photometry of point sources across the Galactic Plane at visible wavelengths, providing the much-needed counterpart to recent infrared surveys.

We present the first results of our search for new, extended Planetary Nebulae (PNe) based on careful, systematic, visual scrutiny of the imaging data from the INT Photometric H-alpha Survey of the Northern Galactic Plane (IPHAS). The newly uncovered PNe will help to improve the census of this important population of Galactic objects that serve as key windows into the late stage evolution of low to intermediate mass stars. They will also facilitate study of the faint end of the ensemble Galactic PN luminosity function. The sensitivity and coverage of IPHAS allows PNe to be found in regions of greater extinction in the Galactic Plane and/or those PNe in a more advanced evolutionary state and at larger distances compared to the general Galactic PN population. Using a set of newly revised optical diagnostic diagrams in combination with access to a powerful, new, multi-wavelength imaging database, we have identified 159 true, likely and possible PNe for this first catalogue release. The ability of IPHAS to unveil PNe at low Galactic latitudes and towards the Galactic Anticenter, compared to previous surveys, makes this survey an ideal tool to contribute to the improvement of our knowledge of the whole Galactic PN population

We present a three dimensional map of extinction in the Northern Galactic Plane derived using photometry from the IPHAS survey. The map has fine angular ($\sim 10$ arcmin) and distance (100 pc) sampling allied to a significant depth ($\gtrsim 5$ kpc). We construct the map using a method based on a hierarchical Bayesian model as previously described by Sale (2012). In addition to mean extinction, we also measure differential extinction, which arises from the fractal nature of the ISM, and show that it will be the dominant source of uncertainty in estimates of extinction to some arbitrary position. The method applied also furnishes us with photometric estimates of the distance, extinction, effective temperature, surface gravity, and mass for $\sim 38$ million stars. Both the extinction map and the catalogue of stellar parameters are made publicly available via http://www.iphas.org/extinction .

We present a model of the photodissociation regions of NGC 4038, which is part of the Antennae galaxies. We have considered one-dimensional slabs of uniform density all having a maximum $A_V=10\,{\rm mag}$, interacting with plane-parallel radiation. The density range in our simulations spans four orders of magnitude ($100\le n\le 10^6\,{\rm cm}^{-3}$) and the UV field strength spans more than three orders of magnitude ($10\le\chi\le10^{4.5}$ multiples of the ${\rm Draine}$ field), from which we generated a grid of about 1400 simulations. We compare our results with Herschel SPIRE-FTS, CSO and ISO-LWS observations of eight CO transition lines ($J=1-0$ to $8-7$) and the [CI] 609\mu m and [OI] 146\mu m fine structure lines. We find that the molecular and atomic emission lines trace different gas components of NGC 4038, thus single emission models are insufficient to reproduce the observed values. In general, low-$J$ CO transition lines correspond to either low density regions interacting with low UV field strengths, or high density regions interacting with high UV field strengths. Higher $J$ CO transition lines are less dependent on the UV field strength and are fitted by gas with density $n\sim10^{4.5}-10^{5.2}\,{\rm cm}^{-3}$. We find that the observed fine structure line ratio of [CI] 609\mu m/[OI] 146\mu m is reproduced by clouds subject to weaker UV fields compared to the CO lines. We make estimates of the $X_{\rm CO}$ factor which relates the CO emission with the column density of molecular hydrogen, and find that it is less than the canonical Milky Way value.

For three nebulae that have early-WN Wolf-Rayet exciting stars, NGC 6888, WR 8 and Abell 48, we have obtained Herschel-PACS line scans of the [N III] 57 um and [O III] 88 micron lines, along with the 122 and 205 micron lines of [N II]. From the former two lines we have derived N$^{2+}$/O$^{2+}$ abundance ratios, equal to the overall N/O ratio under a wide range of nebular conditions. We find that all of the nebulae observed possess significant nitrogen enrichment, with derived N/O ratios greater than solar. The two nebulae with massive Wolf-Rayet exciting stars, NGC 6888 and WR8 are found to have N/O ratios that are enhanced by factors of 7 - 10 relative to the solar N/O ratio, consistent with an origin as material ejected just before the onset of the Wolf-Rayet phase. The other nebula, Abell 48, has recently been reclassified as a member of the rare class of three planetary nebulae that have early-WN central stars and are not of Peimbert Type I. We derive a nebular N/O ratio for it that is a factor of 4 enhanced relative to solar and slightly above the range of N/O values that have been measured for the other three members of its [WN] planetary nebula class.

The Helix Nebula (NGC 7293) is the closest planetary nebulae. Therefore, it is an ideal template for photochemical studies at small spatial scales in planetary nebulae. We aim to study the spatial distribution of the atomic and the molecular gas, and the structure of the photodissociation region along the western rims of the Helix Nebula as seen in the submillimeter range with Herschel. We use 5 SPIRE FTS pointing observations to make atomic and molecular spectral maps. We analyze the molecular gas by modeling the CO rotational lines using a non-local thermodynamic equilibrium (non-LTE) radiative transfer model. For the first time, we have detected extended OH+ emission in a planetary nebula. The spectra towards the Helix Nebula also show CO emission lines (from J= 4 to 8), [NII] at 1461 GHz from ionized gas, and [CI] (2-1), which together with the OH+ lines, trace extended CO photodissociation regions along the rims. The estimated OH+ column density is (1-10)x1e12 cm-2. The CH+ (1-0) line was not detected at the sensitivity of our observations. Non-LTE models of the CO excitation were used to constrain the average gas density (n(H2)=(1-5)x1e5 cm-3) and the gas temperature (Tk= 20-40 K). The SPIRE spectral-maps suggest that CO arises from dense and shielded clumps in the western rims of the Helix Nebula whereas OH+ and [CI] lines trace the diffuse gas and the UV and X-ray illuminated clumps surface where molecules reform after CO photodissociation. [NII] traces a more diffuse ionized gas component in the interclump medium.

We present 48 Herschel/PACS spectra of evolved stars in the wavelength range of 67-72 $\mu$m. This wavelength range covers the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{2-2x}\text{Fe}_{(2x)}\text{SiO}_{4}$). The width and wavelength position of this band are sensitive to the temperature and composition of the crystalline olivine. Our sample covers a wide range of objects: from high mass-loss rate AGB stars (OH/IR stars, $\dot M \ge 10^{-5}$ M$_\odot$/yr), through post-AGB stars with and without circumbinary disks, to planetary nebulae and even a few massive evolved stars. The goal of this study is to exploit the spectral properties of the 69 $\mu$m band to determine the composition and temperature of the crystalline olivine. Since the objects cover a range of evolutionary phases, we study the physical and chemical properties in this range of physical environments. We fit the 69 $\mu$m band and use its width and position to probe the composition and temperature of the crystalline olivine. For 27 sources in the sample, we detected the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{(2-2x)}\text{Fe}_{(2x)}\text{SiO}_{4}$). The 69 $\mu$m band shows that all the sources produce pure forsterite grains containing no iron in their lattice structure. The temperature of the crystalline olivine as indicated by the 69 $\mu$m band, shows that on average the temperature of the crystalline olivine is highest in the group of OH/IR stars and the post-AGB stars with confirmed Keplerian disks. The temperature is lower for the other post-AGB stars and lowest for the planetary nebulae. A couple of the detected 69 $\mu$m bands are broader than those of pure magnesium-rich crystalline olivine, which we show can be due to a temperature gradient in the circumstellar environment of these stars. continued...

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

The VST Photometric Halpha Survey of the Southern Galactic Plane and Bulge (VPHAS+) is surveying the southern Milky Way in u, g, r, i and Halpha at 1 arcsec angular resolution. Its footprint spans the Galactic latitude range -5 < b < +5 at all longitudes south of the celestial equator. Extensions around the Galactic Centre to Galactic latitudes +/-10 bring in much of the Galactic Bulge. This ESO public survey, begun on 28th December 2011, reaches down to 20th magnitude (10-sigma) and will provide single-epoch digital optical photometry for around 300 million stars. The observing strategy and data pipelining is described, and an appraisal of the segmented narrowband Halpha filter in use is presented. Using model atmospheres and library spectra, we compute main-sequence (u - g), (g - r), (r - i) and (r - Halpha) stellar colours in the Vega system. We report on a preliminary validation of the photometry using test data obtained from two pointings overlapping the Sloan Digital Sky Survey. An example of the (u - g, g - r) and (r - Halpha, r - i) diagrams for a full VPHAS+ survey field is given. Attention is drawn to the opportunities for studies of compact nebulae and nebular morphologies that arise from the image quality being achieved. The value of the u band as the means to identify planetary-nebula central stars is demonstrated by the discovery of the central star of NGC 2899 in survey data. Thanks to its excellent imaging performance, the VST/OmegaCam combination used by this survey is a perfect vehicle for automated searches for reddened early-type stars, and will allow the discovery and analysis of compact binaries, white dwarfs and transient sources.

The protonated Argon ion, $^{36}$ArH$^{+}$, has been identified recently in the Crab Nebula (Barlow et al. 2013) from Herschel spectra. Given the atmospheric opacity at the frequency of its $J$=1-0 and $J$=2-1 rotational transitions (617.5 and 1234.6 GHz, respectively), and the current lack of appropriate space observatories after the recent end of the Herschel mission, future studies on this molecule will rely on mid-infrared observations. We report on accurate wavenumber measurements of $^{36}$ArH$^{+}$ and $^{38}$ArH$^{+}$ rotation-vibration transitions in the $v$=1-0 band in the range 4.1-3.7 $\mu$m (2450-2715 cm$^{-1}$). The wavenumbers of the $R$(0) transitions of the $v$=1-0 band are 2612.50135$\pm$0.00033 and 2610.70177$\pm$0.00042 cm$^{-1}$ ($\pm3\sigma$) for $^{36}$ArH$^{+}$ and $^{38}$ArH$^{+}$, respectively. The calculated opacity for a gas thermalized at a temperature of 100 K and a linewidth of 1 km.s$^{-1}$ of the $R$(0) line is $1.6\times10^{-15}\times N$($^{36}$ArH$^+$). For column densities of $^{36}$ArH$^+$ larger than $1\times 10^{13}$ cm$^{-2}$, significant absorption by the $R$(0) line can be expected against bright mid-IR sources.

Noble gas molecules have not hitherto been detected in space. From spectra obtained with the Herschel Space Observatory, we report the detection of emission in the 617.5 GHz and 1234.6 GHz J = 1-0 and 2-1 rotational lines of {36}ArH^+ at several positions in the Crab Nebula, a supernova remnant known to contain both H2 molecules and regions of enhanced ionized argon emission. {36}Ar is believed to have originated from explosive nucleosynthesis in massive stars during core-collapse supernova events. Its detection in the Crab Nebula, the product of such a supernova event, confirms this expectation. The likely excitation mechanism for the observed {36}ArH^+ emission lines is electron collisions in partially ionized regions with electron densities of a few hundred per centimeter cubed.

Supernova (SN) explosions are crucial engines driving the evolution of galaxies by shock heating gas, increasing the metallicity, creating dust, and accelerating energetic particles. In 2012 we used the Atacama Large Millimeter/Submillimeter Array to observe SN 1987A, one of the best-observed supernovae since the invention of the telescope. We present spatially resolved images at 450um, 870um, 1.4mm, and 2.8mm, an important transition wavelength range. Longer wavelength emission is dominated by synchrotron radiation from shock-accelerated particles, shorter wavelengths by emission from the largest mass of dust measured in a supernova remnant (>0.2Msun). For the first time we show unambiguously that this dust has formed in the inner ejecta (the cold remnants of the exploded star's core). The dust emission is concentrated to the center of the remnant, so the dust has not yet been affected by the shocks. If a significant fraction survives, and if SN 1987A is typical, supernovae are important cosmological dust producers.

We present Herschel SPIRE-FTS observations of the Antennae (NGC 4038/39), a well studied, nearby ($22$ Mpc) ongoing merger between two gas rich spiral galaxies. We detect 5 CO transitions ($J=4-3$ to $J=8-7$), both [CI] transitions and the [NII]$205\mu m$ transition across the entire system, which we supplement with ground based observations of the CO $J=1-0$, $J=2-1$ and $J=3-2$ transitions, and Herschel PACS observations of [CII] and [OI]$63\mu m$. Using the CO and [CI] transitions, we perform both a LTE analysis of [CI], and a non-LTE radiative transfer analysis of CO and [CI] using the radiative transfer code RADEX along with a Bayesian likelihood analysis. We find that there are two components to the molecular gas: a cold ($T_{kin}\sim 10-30$ K) and a warm ($T_{kin} \gtrsim 100$ K) component. By comparing the warm gas mass to previously observed values, we determine a CO abundance in the warm gas of $x_{CO} \sim 5\times 10^{-5}$. If the CO abundance is the same in the warm and cold gas phases, this abundance corresponds to a CO $J=1-0$ luminosity-to-mass conversion factor of $\alpha_{CO} \sim 7 \ M_{\odot}{pc^{-2} \ (K \ km \ s^{-1})^{-1}}$ in the cold component, similar to the value for normal spiral galaxies. We estimate the cooling from H$_2$, [CII], CO and [OI]$63\mu m$ to be $\sim 0.01 L_{\odot}/M_{\odot}$. We compare PDR models to the ratio of the flux of various CO transitions, along with the ratio of the CO flux to the far-infrared flux in NGC 4038, NGC 4039 and the overlap region. We find that the densities recovered from our non-LTE analysis are consistent with a background far-ultraviolet field of strength $G_0\sim 1000$. Finally, we find that a combination of turbulent heating, due to the ongoing merger, and supernova and stellar winds are sufficient to heat the molecular gas.

Many nearby main-sequence stars have been searched for debris using the far-infrared Herschel satellite, within the DEBRIS, DUNES and Guaranteed-Time Key Projects. We discuss here 11 stars of spectral types A to M where the stellar inclination is known and can be compared to that of the spatially-resolved dust belts. The discs are found to be well aligned with the stellar equators, as in the case of the Sun's Kuiper belt, and unlike many close-in planets seen in transit surveys. The ensemble of stars here can be fitted with a star-disc tilt of ~<10 degrees. These results suggest that proposed mechanisms for tilting the star or disc in fact operate rarely. A few systems also host imaged planets, whose orbits at tens of AU are aligned with the debris discs, contrary to what might be expected in models where external perturbers induce tilts.

We present an analysis of the far-infrared and submillimetre molecular emission line spectrum of the luminous M-supergiant VY CMa, observed with the SPIRE and PACS spectrometers aboard the Herschel Space Observatory. Over 260 emission lines were detected in the 190-650-micron SPIRE FTS spectra, with one-third of the observed lines being attributable to H2O. Other detected species include CO, 13CO, H2^18O, SiO, HCN, SO, SO2, CS, H2S, and NH3. Our model fits to the observed 12CO and 13CO line intensities yield a 12C/13C ratio of 5.6+-1.8, consistent with measurements of this ratio for other M supergiants, but significantly lower than previously estimated for VY CMa from observations of lower-J lines. The spectral line energy distribution for twenty SiO rotational lines shows two temperature components: a hot component at 1000 K, which we attribute to the stellar atmosphere and inner wind, plus a cooler ~200 K component, which we attribute to an origin in the outer circumstellar envelope. We fit the line fluxes of 12CO, 13CO, H2O and SiO, using the SMMOL non-LTE line transfer code, with a mass-loss rate of 1.85x10^-4 Msun yr^-1 between 9 R* and 350 R*. To fit the observed line fluxes of 12CO, 13CO, H2O and SiO with SMMOL non-LTE line radiative transfer code, along with a mass-loss rate of 1.85x10^-4 Msun yr^-1. To fit the high rotational lines of CO and H2O, the model required a rather flat temperature distribution inside the dust condensation radius, attributed to the high H2O opacity. Beyond the dust condensation radius the gas temperature is fitted best by an r^-0.5 radial dependence, consistent with the coolant lines becoming optically thin. Our H2O emission line fits are consistent with an ortho:para ratio of 3 in the outflow.

We present H\alpha images of an ionized nebula surrounding the M2-5Ia red supergiant (RSG) W26 in the massive star cluster Westerlund 1. The nebula consists of a circumstellar shell or ring ~0.1pc in diameter and a triangular nebula ~0.2pc from the star that in high-resolution Hubble Space Telescope images shows a complex filamentary structure. The excitation mechanism of both regions is unclear since RSGs are too cool to produce ionizing photons and we consider various possibilities. The presence of the nebula, high stellar luminosity and spectral variability suggest that W26 is a highly evolved RSG experiencing extreme levels of mass-loss. As the only known example of an ionized nebula surrounding a RSG W26 deserves further attention to improve our understanding of the final evolutionary stages of massive stars.

We present new photometric data from our Herschel Key Programme, the Dwarf Galaxy Survey (DGS), dedicated to the observation of the gas and dust in 48 low-metallicity environments. They were observed with PACS and SPIRE onboard Herschel at 70,100,160,250,350, and 500 microns. We focus on a systematic comparison of the derived FIR properties (FIR luminosity, dust mass, dust temperature and emissivity index) with more metal-rich galaxies and investigate the detection of a potential submm excess. The data reduction method is adapted for each galaxy to derive the most reliable photometry from the final maps. PACS flux densities are compared with the MIPS 70 and 160 microns bands. We use colour-colour diagrams and modified blackbody fitting procedures to determine the dust properties of the DGS galaxies. We also include galaxies from the Herschel KINGFISH sample, containing more metal-rich environments, totalling 109 galaxies. The location of the DGS galaxies on Herschel colour-colour diagrams highlights the differences in global environments of low-metallicity galaxies. The dust in DGS galaxies is generally warmer than in KINGFISH galaxies (T_DGS~32 K, T_KINGFISH~23 K). The emissivity index, beta, is ~1.7 in the DGS, but metallicity does not make a strong effect on beta. The dust-to-stellar mass ratio is lower in low-metallicity galaxies: M_dust/M_star~0.02% for the DGS vs 0.1% for KINGFISH. Per unit dust mass, dwarf galaxies emit ~6 times more in the FIR than higher metallicity galaxies. Out of the 22 DGS galaxies detected at 500 micron, 41% present an excess in the submm not explained by our dust SED model. The excess mainly appears in lower metallicity galaxies (12+log(O/H) < 8.3), and the strongest excesses are detected in the most metal-poor galaxies. We stress the need for observations longwards of the Herschel wavelengths to detect any submm excess appearing beyond 500 micron.

As part of the Herschel Guaranteed Time Key Project MESS (Mass loss of Evolved StarS) we have imaged a sample of planetary nebulae. In this paper we present the PACS and SPIRE images of the classical bipolar planetary nebula NGC 650. We used these images to derive a temperature map of the dust. We also constructed a photoionization and dust radiative transfer model using the spectral synthesis code Cloudy. To constrain this model, we used the PACS and SPIRE fluxes and combined these with hitherto unpublished IUE and Spitzer IRS spectra as well as various other data from the literature. The temperature map combined with the photoionization model were used to study various aspects of the central star, the nebula, and in particular the dust grains in the nebula. The central star parameters are determined to be T_eff = 208 kK and L = 261 L_sol assuming a distance of 1200 pc. The stellar temperature is much higher than previously published values. We confirm that the nebula is carbon-rich with a C/O ratio of 2.1. The nebular abundances are typical for a type IIa planetary nebula. With the photoionization model we determined that the grains in the ionized nebula are large (assuming single-sized grains, they would have a radius of 0.15 micron. Most likely these large grains were inherited from the asymptotic giant branch phase. The PACS 70/160 micron temperature map shows evidence for two radiation components heating the grains. The first component is direct emission from the central star, while the second component is diffuse emission from the ionized gas (mainly Ly alpha). We show that previous suggestions that there is a photo-dissociation region surrounding the ionized region are incorrect. The neutral material resides in dense clumps inside the ionized region. These may also harbor stochastically heated very small grains in addition to the large grains.

We report spectroscopic and imaging observations of rotational transitions of cold CO and SiO in the ejecta of SN1987A, the first such emission detected in a supernova remnant. In addition to line luminosities for the CO J=1-0, 2-1, 6-5, and 7-6 transitions, we present upper limits for all other transitions up to J=13-12, collectively measured from the Atacama Large Millimeter Array (ALMA), the Atacama Pathfinder EXperiment (APEX), and the Herschel Spectral and Photometric Imaging REceiver (SPIRE). Simple models show the lines are emitted from at least 0.01 solar masses of CO at a temperature > 14 K, confined within at most 35% of a spherical volume expanding at ~ 2000 km/s. Moreover, we locate the emission within 1'' of the central debris. These observations, along with a partial observation of SiO, confirm the presence of cold molecular gas within supernova remnants and provide insight into the physical conditions and chemical processes in the ejecta. Furthermore, we demonstrate the powerful new window into supernova ejecta offered by submillimeter observations.

The Dwarf Galaxy Survey (DGS) program is studying low-metallicity galaxies using 230h of far-infrared (FIR) and submillimetre (submm) photometric and spectroscopic observations of the Herschel Space Observatory and draws to this a rich database of a wide range of wavelengths tracing the dust, gas and stars. This sample of 50 galaxies includes the largest metallicity range achievable in the local Universe including the lowest metallicity (Z) galaxies, 1/50 Zsun, and spans 4 orders of magnitude in star formation rates. The survey is designed to get a handle on the physics of the interstellar medium (ISM) of low metallicity dwarf galaxies, especially on their dust and gas properties and the ISM heating and cooling processes. The DGS produces PACS and SPIRE maps of low-metallicity galaxies observed at 70, 100, 160, 250, 350, and 500 mic with the highest sensitivity achievable to date in the FIR and submm. The FIR fine-structure lines, [CII] 158 mic, [OI] 63 mic, [OI] 145 mic, [OIII] 88 mic, [NIII] 57 mic and [NII] 122 and 205 mic have also been observed with the aim of studying the gas cooling in the neutral and ionized phases. The SPIRE FTS observations include many CO lines (J=4-3 to J=13-12), [NII] 205 mic and [CI] lines at 370 and 609 mic. This paper describes the sample selection and global properties of the galaxies, the observing strategy as well as the vast ancillary database available to complement the Herschel observations. The scientific potential of the full DGS survey is described with some example results included.

Aim : In order to study the history of mass loss in extreme OH/IR stars, we observed a number of these objects using CO as a tracer of the density and temperature structure of their circumstellar envelopes. Method : Combining CO observations from the Herschel Space Observatory with those from the ground, we trace mass loss rates as a function of radius in five extreme OH/IR stars. Using radiative transfer modelling, we modelled the dusty envelope as well as the CO emission. The high-rotational transitions of CO indicate that they originate in a dense superwind region close to the star while the lower transitions tend to come from a more tenuous outer wind which is a result of the mass loss since the early AGB phase. Result : The models of the circumstellar envelopes around these stars suggest that they have entered a superwind phase in the past 200 - 500 years. The low 18O/17O (~ 0.1 compared to the solar abundance ratio of ~ 5) and 12C/13C (3-30 cf. the solar value of 89) ratios derived from our study support the idea that these objects have undergone hot-bottom burning and hence that they are massive M >= 5 solar-mass AGB stars.

[Abridged] We present an integrated photometric spectral energy distribution (SED) of the Magellanic-type galaxy NGC 4449 from the far-ultraviolet (UV) to the submillimetre, including new observations acquired by the Herschel Space Observatory. We include integrated UV photometry from the Swift Ultraviolet and Optical Telescope using a measurement technique which is appropriate for extended sources with coincidence loss. In this paper, we examine the available multiwavelength data to infer a range of ages, metallicities and star formation rates for the underlying stellar populations, as well as the composition and the total mass of dust in NGC 4449. We present an iterative scheme, which allows us to build an in-depth and multicomponent representation of NGC 4449 `bottom-up', taking advantage of the broad capabilities of the photoionization and radiative transfer code MOCASSIN (MOnte CArlo SimulationS of Ionized Nebulae). We fit the observed SED, the global ionization structure and the emission line intensities, and infer a recent SFR of 0.4 Msolar/yr and a total stellar mass of approximately 1e9 Msolar emitting with a bolometric luminosity of 5.7e9 Lsolar. Our fits yield a total dust mass of 2.9e6 Msolar including 2 per cent attributed to polycyclic aromatic hydrocarbons. We deduce a dust to gas mass ratio of 1/190 within the modelled region. While we do not consider possible additional contributions from even colder dust, we note that including the extended HI envelope and the molecular gas is likely to bring the ratio down to as low as ~ 1/800.

We present the sub-millimeter spectra from 450 GHz to 1550 GHz of eleven nearby active galaxies observed with the SPIRE Fourier Transform Spectrometer (SPIRE/FTS) onboard Herschel. We detect CO transitions from J_up = 4 to 12, as well as the two [CI] fine structure lines at 492 and 809 GHz and the [NII] 461 GHz line. We used radiative transfer models to analyze the observed CO spectral line energy distributions (SLEDs). The FTS CO data were complemented with ground-based observations of the low-J CO lines. We found that the warm molecular gas traced by the mid-J CO transitions has similar physical conditions (n_H2 ~ 10^3.2 - 10^3.9 cm^-3 and T_kin ~ 300 - 800 K) in most of our galaxies. Furthermore, we found that this warm gas is likely producing the mid-IR rotational H2 emission. We could not determine the specific heating mechanism of the warm gas, however it is possibly related to the star-formation activity in these galaxies. Our modeling of the [CI] emission suggests that it is produced in cold (T_kin < 30 K) and dense (n_H2 > 10^3 cm^-3) molecular gas. Transitions of other molecules are often detected in our SPIRE/FTS spectra. The HF J=1-0 transition at 1232 GHz is detected in absorption in UGC05101 and in emission in NGC7130. In the latter, near-infrared pumping, chemical pumping, or collisional excitation with electrons are plausible excitation mechanisms likely related to the AGN of this galaxy. In some galaxies few H2O emission lines are present. Additionally, three OH+ lines at 909, 971, and 1033 GHz are identified in NGC7130.

We study for the first time the complete sub-millimeter spectra (450 GHz to 1550 GHz) of a sample of nearby active galaxies observed with the SPIRE Fourier Transform Spectrometer (SPIRE/FTS) onboard Herschel. The CO ladder (from Jup = 4 to 12) is the most prominent spectral feature in this range. These CO lines probe warm molecular gas that can be heated by ultraviolet photons, shocks, or X-rays originated in the active galactic nucleus or in young star-forming regions. In these proceedings we investigate the physical origin of the CO emission using the averaged CO spectral line energy distribution (SLED) of six Seyfert galaxies. We use a radiative transfer model assuming an isothermal homogeneous medium to estimate the molecular gas conditions. We also compare this CO SLED with the predictions of photon and X-ray dominated region (PDR and XDR) models.

We investigate a region of the Galactic plane, between 120 <= l <= 140 and -1 <= b <= +4, and uncover a population of moderately reddened (E(B-V) \sim 1) classical Be stars within and beyond the Perseus and Outer Arms. 370 candidate emission line stars (13 <= r <= 16) selected from the INT Photometric H-alpha Survey of the Northern Galactic plane (IPHAS) have been followed up spectroscopically. A subset of these, 67 stars with properties consistent with those of classical Be stars, have been observed at sufficient spectral resolution (Delta_lambda \sim 2 - 4 Angstrom) at blue wavelengths to narrow down their spectral types. We determine these to a precision estimated to be +/- 1 sub-type and then we measure reddenings via SED fitting with reference to appropriate model atmospheres. Corrections for contribution to colour excess from circumstellar discs are made using an established scaling to H-alpha emission equivalent width. Spectroscopic parallaxes are obtained after luminosity class has been constrained via estimates of distances to neighbouring A/F stars with similar reddenings. Overwhelmingly, the stars in the sample are confirmed as luminous classical Be stars at heliocentric distances ranging from 2 kpc up to \sim 12 kpc. However, the errors are presently too large to enable the cumulative distribution function with respect to distance to distinguish between models placing the stars exclusively in spiral arms, or in a smooth exponentially-declining distribution.

Using VLT/FLAMES optical integral field unit observations, we present the first spatially resolved spectroscopic study of the well known blue compact galaxy Haro11, thought to be a local analogue to high redshift Lyman Break Galaxies. Haro11 displays complex emission line profiles, consisting of narrow (FWHM<200 km/s) and broad (FWHM 200-300 km/s) components. We identify three distinct emission knots kinematically connected to one another. A chemodynamical analysis is presented, revealing that spatially resolved ionic and elemental abundances do not agree with those derived from integrated spectra across the galaxy. We conclude that this is almost certainly due to the surface-brightness-weighting of electron temperature in integrated spectra, leading to higher derived abundances. We find that the eastern knot has a low gas density, but a higher temperature (by ~4,000 K) and consequently an oxygen abundance ~0.4 dex lower than the neighbouring regions. A region of enhanced N/O is found specifically in Knot C, confirming previous studies that found anomalously high N/O ratios in this system. Maps of the Wolf-Rayet feature at 4686A reveal large WR populations (~900-1500 stars) in Knots A and B. The lack of WR stars in Knot C, combined with an age of ~7.4 Myr suggest that a recently completed WR phase may be responsible for the observed N/O excess. Conversely, the absence of N-enriched gas and strong WR emission in Knots A and B suggests we are observing these regions at an epoch where stellar ejecta has yet to cool and mix with the ISM.

Context. The interaction between stellar winds and the interstellar medium (ISM) can create complex bow shocks. The photometers on board the Herschel Space Observatory are ideally suited to studying the morphologies of these bow shocks. Aims. We aim to study the circumstellar environment and wind-ISM interaction of the nearest red supergiant, Betelgeuse. Methods. Herschel PACS images at 70, 100, and 160 micron and SPIRE images at 250, 350, and 500 micron were obtained by scanning the region around Betelgeuse. These data were complemented with ultraviolet GALEX data, near-infrared WISE data, and radio 21 cm GALFA-HI data. The observational properties of the bow shock structure were deduced from the data and compared with hydrodynamical simulations. Results. The infrared Herschel images of the environment around Betelgeuse are spectacular, showing the occurrence of multiple arcs at 6-7 arcmin from the central target and the presence of a linear bar at 9 arcmin. Remarkably, no large-scale instabilities are seen in the outer arcs and linear bar. The dust temperature in the outer arcs varies between 40 and 140 K, with the linear bar having the same colour temperature as the arcs. The inner envelope shows clear evidence of a non-homogeneous clumpy structure (beyond 15 arcsec), probably related to the giant convection cells of the outer atmosphere. The non-homogeneous distribution of the material even persists until the collision with the ISM. A strong variation in brightness of the inner clumps at a radius of 2 arcmin suggests a drastic change in mean gas and dust density some 32 000 yr ago. Using hydrodynamical simulations, we try to explain the observed morphology of the bow shock around Betelgeuse. Conclusions: [abbreviated]

Photodissociation regions (PDRs) define the transition zone between an ionized and a dark molecular region. They consist of neutral gas which interacts with far-ultraviolet radiation and are characterized by strong infrared line emission. Various numerical codes treating one-dimensional PDRs have been developed in the past, simulating the complexity of chemical reactions occurring and providing a better understanding of the structure of a PDR. In this paper we present the three-dimensional code, 3D-PDR, which can treat PDRs of arbitrary density distribution. The code solves the chemistry and the thermal balance self-consistently within a given three-dimensional cloud. It calculates the total heating and cooling functions at any point in a given PDR by adopting an escape probability method. It uses a HEALPix-based ray-tracing scheme to evaluate the attenuation of the far-ultraviolet radiation in the PDR and the propagation of the far-infrared/submm line emission out of the PDR. We present benchmarking results and apply 3D-PDR to i) a uniform-density spherical cloud interacting with a plane-parallel external radiation field, ii) a uniform-density spherical cloud interacting with a two-component external radiation field, and iii) a cometary globule interacting with a plane-parallel external radiation field. We find that the code is able to reproduce the benchmarking results of various other one-dimensional numerical codes treating PDRs. We also find that the accurate treatment of the radiation field in the fully three-dimensional treatment of PDRs can in some cases leads to different results when compared to a standard one-dimensional treatment.

Some planetary systems harbour debris disks containing planetesimals such as asteroids and comets. Collisions between such bodies produce small dust particles, the spectral features of which reveal their composition and, hence, that of their parent bodies. A measurement of the composition of olivine crystals has been done for the protoplanetary disk HD100546 and for olivine crystals in the warm inner parts of planetary systems. The latter compares well with the iron-rich olivine in asteroids (x<0.29). In the cold outskirts of the Beta Pictoris system, an analogue to the young Solar System, olivine crystals were detected but their composition remained undetermined, leaving unknown how the composition of the bulk of Solar System cometary olivine grains compares with that of extrasolar comets. Here we report the detection of the 69-micrometre-wavelength band of olivine crystals in the spectrum of Beta Pictoris. Because the disk is optically thin, we can associate the crystals with an extrasolar proto-Kuiper belt a distance of 15-45 astronomical units from the star (one astronomical unit is the Sun-Earth distance), determine their magnesium-rich composition (x=0.01+/-0.001) and show that they make up 3.6+/-1.0 per cent of the total dust mass. These values are strikingly similar to those for the dust emitted by the most primitive comets in the Solar System, even though Beta Pictoris is more massive and more luminous and has a different planetary system architecture.

Whether supernovae are major sources of dust in galaxies is a long-standing debate. We present infrared and submillimeter photometry and spectroscopy from the Herschel Space Observatory of the Crab Nebula between 51 and 670 micron as part of the Mass Loss from Evolved StarS program. We compare the emission detected with Herschel with multiwavelength data including millimeter, radio, mid-infrared and archive optical images. We carefully remove the synchrotron component using the Herschel and Planck fluxes measured in the same epoch. The contribution from line emission is removed using Herschel spectroscopy combined with Infrared Space Observatory archive data. Several forbidden lines of carbon, oxygen and nitrogen are detected where multiple velocity components are resolved, deduced to be from the nitrogen-depleted, carbon-rich ejecta. No spectral lines are detected in the SPIRE wavebands; in the PACS bands, the line contribution is 5% and 10% at 70 and 100 micron and negligible at 160 micron. After subtracting the synchrotron and line emission, the remaining far-infrared continuum can be fit with two dust components. Assuming standard interstellar silicates, the mass of the cooler component is 0.24(+0.32)(-0.08) Msolar for T = 28.1(+5.5)(-3.2)K. Amorphous carbon grains require 0.11 +/- 0.01 Msolar of dust with T = 33.8(+2.3)(-1.8) K. A single temperature modified-blackbody with 0.14Msolar and 0.08Msolar for silicate and carbon dust respectively, provides an adequate fit to the far- infrared region of the spectral energy distribution but is a poor fit at 24-500 micron. The Crab Nebula has condensed most of the relevant refractory elements into dust, suggesting the formation of dust in core-collapse supernova ejecta is efficient.

We present new Herschel/PACS images at 70, 100, and 160 micron of the well-known, nearby, carbon-rich asymptotic giant branch star IRC+10216 revealing multiple dust shells in its circumstellar envelope. For the first time, dust shells (or arcs) are detected until 320 arcsec. The almost spherical shells are non-concentric and have an angular extent between 40 deg and 200 deg. The shells have a typical width of 5 arcsec - 8 arcsec, and the shell separation varies in the range of 10 arcsec - 35 arcsec, corresponding to 500-1700 yr. Local density variations within one arc are visible. The shell/intershell density contrast is typically 4, and the arcs contain some 50% more dust mass than the smooth envelope. The observed (nested) arcs record the mass-loss history over the past 16 000 yr, but Rayleigh-Taylor and Kelvin-Helmholtz instabilities in the turbulent astropause and astrosheath will erase any signature of the mass-loss history for at least the first 200 000 yr of mass loss. Accounting for the bowshock structure, the envelope mass around IRC+10216 contains >2Msun of gas and dust mass. It is argued that the origin of the shells is related to non-isotropic mass-loss events and clumpy dust formation.

Using VLT/FLAMES optical IFU observations, we present a detailed study of UM 448, a nearby Blue Compact Galaxy (BCG) previously reported to have an anomalously high N/O abundance ratio. NTT/SuSI2 images reveal a morphology suggestive of a merger of two systems of contrasting colour, whilst our H-alpha emission maps resolve UM 448 into three separate regions that do not coincide with the stellar continuum peaks. UM 448 exhibits complex emission line profiles, with lines consisting of a narrow, central component, an underlying broad component and a third, narrow blue-shifted component. Radial velocity maps show signs of solid body rotation across UM 448, with a projected rotation axis that correlates with the continuum morphology of the galaxy. A spatially-resolved, chemodynamical analysis is presented. Whilst the eastern tail of UM 448 has electron temperatures (Te) that are typical of BCGs, we find a region within the main body of the galaxy where the narrow and broad [O III] 4363 line components trace temperatures differing by 5000K and oxygen abundances differing by 0.4 dex. We measure spatially resolved and integrated ionic and elemental abundances for O, N, S and Ne throughout UM 448, and find they do not agree, possibly due the flux-weighting of Te from the integrated spectrum. This has significant implications for abundances derived from long-slit and integrated spectra of star-forming galaxies in the nearby and distant universe. A region of enhanced N/O ratio is indeed found, extended over a ~0.6 kpc^2 region within the main body of the galaxy. Contrary to previous studies, however, we do not find evidence for a large Wolf-Rayet population, and conclude that WR stars alone cannot be responsible for producing the observed N/O excess. Instead, the location and disturbed morphology of the N-enriched region suggests that interaction-induced inflow of metal-poor gas may be responsible.

(abridged) Our goal is to describe the multi-phase ISM of the IR bright low-metallicity galaxy Haro 11, dissecting the photoionised and photodissociated gas components. We present observations of the mid- and far-IR fine-structure cooling lines obtained with the Spitzer/IRS and Herschel/PACS spectrometers. We use the spectral synthesis code Cloudy to methodically model the ionised and neutral gas from which these lines originate. We find that the mid- and far-IR lines account for ~1% of the total IR luminosity L_TIR. Haro 11 is undergoing a phase of intense star formation, as traced by the brightest line [OIII] 88um, with L_[OIII]/L_TIR ~0.3%, and high ratios of [NeIII]/[NeII] and [SIV]/[SIII]. Due to their different origins, the observed lines require a multi-phase modeling comprising: a compact HII region, dense fragmented photodissociation regions (PDRs), a diffuse extended low-ionisation/neutral gas which has a volume filling factor of at least 90%, and porous warm dust in proximity to the stellar source. For a more realistic picture of the ISM of Haro 11 we would need to model the clumpy source and gas structures. We combine these 4 model components to explain the emission of 17 spectral lines, investigate the global energy balance of the galaxy through its spectral energy distribution, and establish a phase mass inventory. While the ionic emission lines of Haro 11 essentially originate from the dense HII region component, a diffuse low-ionisation gas is needed to explain the [NeII], [NII], and [CII] line intensities. The [OIII] 88um line intensity is not fully reproduced by our model, hinting towards the possible presence of yet another low-density high-ionisation medium. The [OI] emission is consistent with a dense PDR of low covering factor, and we find no evidence for an X-ray dominated component. The PDR component accounts for only 10% of the [CII] emission.

Fine-structure mid-infrared emission lines with critical densities in the regime 10^4 to 10^6 cm^-3 can be employed to probe the outflow from Wolf-Rayet stars at radii of \sim 10^15 cm. Narrow-band mid-infrared imaging and spectroscopy of the nearest WR star to the sun, gamma Velorum is analysed for spatially resolved forbidden line emission in the Wolf-Rayet outer wind. The [S IV] 10.52 and [Ne II] 12.81 micron emission regions are found to be spatially extended, compared to unresolved continuum and He and C recombination line emission. The [S IV] and [{Ne II] emission line distributions have a high degree of azimuthal symmetry, indicating a spherically symmetric outflow. A model wind with a modest degree of clumping (clumping factor f \sim 10) provides a better match to the observations than an unclumped model. The overall line intensity distributions are consistent with a freely expanding, spherically symmetric 1/$r^2$ outflow with constant ionization fraction and modestly clumped density structure.

V1280 Sco is one of the slowest dust-forming nova ever historically observed. We performed multi-epoch high-spatial resolution observations of the circumstellar dusty environment of V1280 Sco to investigate the level of asymmetry of the ejecta We observed V1280 Sco in 2009, 2010 and 2011 using unprecedented high angular resolution techniques. We used the NACO/VLT adaptive optics system in the J, H and K bands, together with contemporaneous VISIR/VLT mid-IR imaging that resolved the dust envelope of V1280 Sco, and SINFONI/VLT observations secured in 2011. We report the discovery of a dusty hourglass-shaped bipolar nebula. The apparent size of the nebula increased from 0.30" x 0.17" in July 2009 to 0.64" x 0.42" in July 2011. The aspect ratio suggests that the source is seen at high inclination. The central source shines efficiently in the K band and represents more than 56+/-5% of the total flux in 2009, and 87+/-6% in 2011. A mean expansion rate of 0.39+/-0.03 mas per day is inferred from the VISIR observations in the direction of the major axis, which represents a projected upper limit. Assuming that the dust shell expands in that direction as fast as the low-excitation slow ejecta detected in spectroscopy, this yields a lower limit distance to V1280 Sco of 1kpc; however, the systematic errors remain large due to the complex shape and velocity field of the dusty ejecta. The dust seems to reside essentially in the polar caps and no infrared flux is detected in the equatorial regions in the latest dataset. This may imply that the mass-loss was dominantly polar.

CW Leo has been observed six times between October 2009 and June 2012 with the SPIRE instrument on board the Herschel satellite. Variability has been detected in the flux emitted by the central star with a period of 639 \pm 4 days, in good agreement with determinations in the literature. Variability is also detected in the bow shock around CW Leo that had previously been detected in the ultraviolet and Herschel PACS/SPIRE data. Although difficult to prove directly, our working hypothesis is that this variability is directly related to that of the central star. In this case, fitting a sine curve with the period fixed to 639 days results in a time-lag in the variability between bow shock and the central star of 402 \pm 37 days. The orientation of the bow shock relative to the plane of the sky is unknown (but see below). For an inclination angle of zero degrees, the observed time-lag translates into a distance to CW Leo of 130 \pm 13 pc, and for non-zero inclination angles the distance is smaller. Fitting the shape of the bow shock with an analytical model (Wilkin 1996), the effect of the inclination angle on the distance may be estimated. Making the additional assumption that the relative peculiar velocity between the interstellar medium (ISM) and CW Leo is determined entirely by the star space velocity with respect to the local standard of rest (i.e. a stationary ISM), the inclination angle is found to be (-33.3 \pm 0.8) degrees based on the observed proper motion and radial velocity. Using the Wilkin model, our current best estimate of the distance to CW Leo is 123 \pm 14 pc. For a distance of 123 pc, we derive a mean luminosity of 7790 \pm 150 Lsol (internal error).

We present nebular phase optical imaging and spectroscopy and near/mid-IR imaging of the Type II SN 2006bc. Observations reveal the central wavelength of the symmetric H$\alpha$ line profile to be red-shifted with respect to the host galaxy H$\alpha$ emission by day 325. Such an phenomenon has been argued to result from an asymmetric explosion in the iron-peak elements resulting in a larger mass of $^{56}$Ni and higher excitation of hydrogen on the far side of the SN explosion. We also observe a gradual blue-shifting of this H$\alpha$ peak which is indicative of dust formation in the ejecta. Although showing a normal peak brightness, V $\sim$ -17.2, for a core-collapse SN, 2006bc fades by $\sim$6 mag during the first 400 days suggesting either a relatively low $^{56}$Ni yield, an increase in extinction due to new dust, or both. A short duration flattening of the light curve is observed from day 416 to day 541 suggesting an optical light echo. Based on the narrow time window of this echo, we discuss implications on the location and geometry of the reflecting ISM. With our radiative transfer models, we find an upper limit of 2 x 10$^{-3}$ M$_{\odot}$ of dust around SN 2006bc. In the event that all of this dust were formed during the SN explosion, this quantity of dust is still several orders of magnitude lower than that needed to explain the large quantities of dust observed in the early universe.

Fomalhaut is a young, nearby star that is suspected to harbor an infant planetary system, interspersed with one or more belts of dusty debris. We present far-infrared images obtained with the Herschel Space Observatory with an angular resolution between 5.7 and 36.7 arcsec at wavelengths between 70 and 500 micrometer. The images show the main debris belt in great detail. Even at high spatial resolution, the belt appears smooth. The region in between the belt and the central star is not devoid of material; thermal emission is observed here as well. Also at the location of the star, excess emission is detected. We use a dynamical model together with radiative-transfer tools to derive the parameters of the debris disk. We include detailed models of the interaction of the dust grains with radiation, for both the radiation pressure and the temperature determination. Comparing these models to the spatially resolved temperature information contained in the images allows us to place strong constraints on the presence of grains that will be blown out of the system by radiation pressure. We use this to derive the dynamical parameters of the system. The appearance of the belt points towards a remarkably active system in which dust grains are produced at a very high rate by a collisional cascade in a narrow region filled with dynamically excited planetesimals. Dust particles with sizes below the blow-out size are abundantly present. The equivalent of 2000 one-km-sized comets are destroyed every day, out of a cometary reservoir amounting to 110 Earth masses. From comparison of their scattering and thermal properties, we find evidence that the dust grains are fluffy aggregates, which indicates a cometary origin. The excess emission at the location of the star may be produced by hot dust with a range of temperatures, but may also be due to gaseous free-free emission from a stellar wind.

We report optical and mid-infrared photometry of SN 1980K between 2004 and 2010, which show slow monotonic fading consistent with previous spectroscopic and photometric observations made 8 to 17 years after outburst. The slow rate-of-change over two decades suggests that this evolution may result from scattered and thermal light echoes off of extended circumstellar material. We present a semi- analytic dust radiative-transfer model that uses an empirically corrected effective optical depth to provide a fast and robust alternative to full Monte-Carlo radiative transfer modeling for homogenous dust at low to intermediate optical depths. We find that unresolved echoes from a thin circumstellar shell 14-15 lt-yr from the progenitor, and containing about 0.02 Msun of carbon-rich dust, can explain the broadband spectral and temporal evolution. The size, mass and dust composition are in good agreement with the contact discontinuity observed in scattered echoes around SN 1987A. The origin of slowly-changing high-velocity [O I] and Halpha lines is also considered. We propose an origin in shocked high-velocity metal-rich clumps of ejecta, rather than arising in the impact of ejecta on slowly-moving circumstellar material, as is the case with hot spots in SN 1987A.

The C-rich AGB star IRC+10216 undergoes strong mass loss, and quasi-periodic density enhancements in the circumstellar matter have been reported. CO is ubiquitous in the CSE, while CCH emission comes from a spatially confined shell. With the IRAM 30m telescope and Herschel/HIFI, we recently detected unexpectedly strong emission from the CCH N=4-3, 6-5, 7-6, 8-7, and 9-8 transitions, challenging the available chemical and physical models. We aim to constrain the physical properties of IRC+10216's CSE, including the effect of episodic mass loss on the observed emission. In particular, we aim to determine the excitation region and conditions of CCH and to reconcile these with interferometric maps of the N=1-0 transition. Via radiative-transfer modelling, we provide a physical description of the CSE, constrained by the SED and a sample of 20 high-resolution and 29 low-resolution CO lines. We further present detailed radiative-transfer analysis of CCH. Assuming a distance of 150pc, the SED is modelled with a stellar luminosity of 11300Lsun and a dust-mass-loss rate of 4.0\times10^{-8}Msun/yr. Based on the analysis of 20 high resolution CO observations, an average gas-mass-loss rate for the last 1000yrs of 1.5\times10^{-5}Msun/yr is derived. This gives a gas-to-dust-mass ratio of 375, typical for an AGB star. The gas kinetic temperature throughout the CSE is described by 3 powerlaws: it goes as r^{-0.58} for r<9R*, as r^{-0.40} for 9<=r<=65R*, and as r^{-1.20} for r>65R*. This model successfully describes all 49 CO lines. We show the effect of wind-density enhancements on the CCH-abundance profile, and the good agreement of the model with the CCH N=1-0 transition and with the lines observed with the 30m telescope and HIFI. We report on the importance of radiative pumping to the vibrationally excited levels of CCH and the significant effect this has on the excitation of all levels of the CCH-molecule.

The origin of interstellar dust in galaxies is poorly understood, particularly the relative contributions from supernovae and the cool stellar winds of low-intermediate mass stars. Here, we present Herschel PACS and SPIRE photometry at 70-500um of the historical young supernova remnants: Kepler and Tycho; both thought to be the remnants of Type Ia explosion events. We detect a warm dust component in Kepler's remnant with T = 82K and mass 0.0031Msun; this is spatially coincident with thermal X-ray emission optical knots and filaments, consistent with the warm dust originating in the circumstellar material swept up by the primary blast wave of the remnant. Similarly for Tycho's remnant, we detect warm dust at 90K with mass 0.0086Msun. Comparing the spatial distribution of the warm dust with X-rays from the ejecta and swept-up medium, and Ha emission arising from the post-shock edge, we show that the warm dust is swept up interstellar material. We find no evidence of a cool (25-50 K) component of dust with mass >0.07Msun as observed in core-collapse remnants of massive stars. Neither the warm or cold dust components detected here are spatially coincident with supernova ejecta material. We compare the lack of observed supernova dust with a theoretical model of dust formation in Type Ia remnants which predicts dust masses of 0.088(0.017)Msun for ejecta expanding into surrounding densities of 1(5)cm-3. The model predicts that silicon- and carbon-rich dust grains will encounter the interior edge of the observed dust emission at 400 years confirming that the majority of the warm dust originates from swept up circumstellar or interstellar grains (for Kepler and Tycho respectively). The lack of cold dust grains in the ejecta suggests that Type Ia remnants do not produce substantial quantities of iron-rich dust grains and has important consequences for the 'missing' iron mass observed in ejecta.

We use the photoionisation and dust radiative transfer code MOCASSIN to create a model of the dwarf irregular galaxy NGC 4449. The best-matching model reproduces the global optical emission line fluxes and the observed spectral energy distribution (SED) spanning wavelengths from the UV to sub-mm, and requires the bolometric luminosity of 6.25e9 Lsolar for the underlying stellar component, M_d/M_g of 1/680 and M_d of 2.2e6 Msolar.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

We present BVRIJHK band photometry of 6 core-collapse supernovae, SNe 1999bw, 2002hh, 2003gd, 2004et, 2005cs, and 2006bc measured at late epochs (>2 yrs) based on Hubble Space Telescope (HST), Gemini north, and WIYN telescopes. We also show the JHK lightcurves of a supernova impostor SN 2008S up to day 575. Of our 43 HST observations in total, 36 observations are successful in detecting the light from the SNe alone and measuring magnitudes of all the targets. HST observations show a resolved scattered light echo around SN 2003gd at day 1520 and around SN 2002hh at day 1717. Our Gemini and WIYN observations detected SNe 2002hh and 2004et, as well. Combining our data with previously published data, we show VRIJHK-band lightcurves and estimate decline magnitude rates at each band in 4 different phases. Our prior work on these lightcurves and other data indicate that dust is forming in our targets from day ~300-400, supporting SN dust formation theory. In this paper we focus on other physical properties derived from the late time light curves. We estimate 56Ni masses for our targets (0.5-14 x 10^{-2} Msun) from the bolometric lightcurve of each for days ~150-300 using SN 1987A as a standard (7.5 x 10^{-2} Msun). The flattening or sometimes increasing fluxes in the late time light curves of SNe 2002hh, 2003gd, 2004et and 2006bc indicate the presence of light echos. We estimate the circumstellar hydrogen density of the material causing the light echo and find that SN 2002hh is surrounded by relatively dense materials (n(H) >400 cm^{-3}) and SNe 2003gd and 2004et have densities more typical of the interstellar medium (~1 cm^{-3}). The 56Ni mass appears well correlated with progenitor mass with a slope of 0.31 x 10^{-2}, supporting the previous work by Maeda et al. (2010), who focus on more massive Type II SNe. The dust mass does not appear to be correlated with progenitor mass.

HII regions and planetary nebulae (PNe) both emit at radio and infrared (IR) wavelengths, and angularly small HII regions can be mistaken for PNe. This problem of classification is most severe for HII regions in an early evolutionary stage, those that are extremely distant, or those that are both young and distant. Using data from the Herschel Hi-GAL survey, as well as WISE and the Spitzer MIPSGAL and GLIMPSE surveys, we wish to establish characteristic IR colors that can be used to distinguish between HII regions and PNe. We perform aperture photometry measurements for a sample of 126 HII regions and 43 PNe at wavelengths from 8.0 micron to 500 micron. We find that HII regions and PNe have distinct IR colors. The most robust discriminating color criteria are [F12/F8] < 0.3, [F160/F12] > 1.3, and [F160/F24] > 0.8 (or alternately [F160/F22] > 0.8), where the brackets indicate the log of the flux ratio. All three of these criteria are individually satisfied by over 98% of our sample of HII regions and by about 10% of our sample of PNe. Combinations of these colors are more robust in separating the two populations; for example all HII regions and no PNe satisfy [F12/F8] < 0.4 and [F160/F22] > 0.8. When applied to objects of unknown classification, these criteria prove useful in separating the two populations. The dispersion in color is relatively small for HII regions; this suggests that any evolution in these colors with time for HII regions must be relatively modest. The spectral energy distributions (SEDs) of HII regions can be separated into "warm" and "cold" components. The "cold" component is well-fit by a grey-body of temperature 25K. The SEDs of nearly two-thirds of our sample of HII regions peak at 160 micron and one third peak at 70 micron. For PNe, 67% of the SEDs peak at 70 micron, 23% peak at either 22 micron or 24 micron, and 9% (two sources) peak at 160 micron.

In this paper we give a progress report on the Herschel imaging and spectroscopic observations of planetary nebulae that are carried out as part of the MESS guaranteed time key program. We present and discuss imaging and temperature maps of NGC 6720, NGC 650, and NGC 6853, as well as PACS and SPIRE spectroscopy of NGC 7027.

In 2007, R Coronae Borealis (R CrB) went into an historically deep and long decline. In this state, the dust acts like a natural coronagraph at visible wavelengths, allowing faint nebulosity around the star to be seen. Imaging has been obtained from 0.5 to 500 micron with Gemini/GMOS, HST/WFPC2, Spitzer/MIPS, and Herschel/SPIRE. Several of the structures around R CrB are cometary globules caused by wind from the star streaming past dense blobs. The estimated dust mass of the knots is consistent with their being responsible for the R CrB declines if they form along the line of sight to the star. In addition, there is a large diffuse shell extending up to 4 pc away from the star containing cool 25 K dust that is detected all the way out to 500 micron. The SED of R CrB can be well fit by a 150 AU disk surrounded by a very large diffuse envelope which corresponds to the size of the observed nebulosity. The total masses of the disk and envelope are 10^-4 and 2 M(Sun), respectively, assuming a gas-to-dust ratio of 100. The evidence pointing toward a white-dwarf merger or a final-helium-shell flash origin for R CrB is contradictory. The shell and the cometary knots are consistent with a fossil planetary nebula. Along with the fact that R CrB shows significant Lithium in its atmosphere, this supports the final-helium-shell flash. However, the relatively high inferred mass of R CrB and its high fluorine abundance support a white-dwarf merger.

We present an analysis of multi-epoch observations of the Type II-P supernova SN 2004et. New and archival optical spectra of SN 2004et are used to study the evolution of the Halpha and [O I] 6300A line profiles between days 259 and 646. Mid-infrared imaging was carried out between 2004 to 2010. We include Spitzer `warm' mission photometry at 3.6 and 4.5um obtained on days 1779, 1931 and 2151, along with ground-based and HST optical and near-infrared observations obtained between days 79 and 1803. Multi-wavelength light curves are presented, as well as optical-infrared spectral energy distributions (SEDs) for multiple epochs. Starting from about day 300, the optical light curves provide evidence for an increasing amount of circumstellar extinction attributable to newly formed dust, with the additional extinction reaching 0.8-1.5 magnitudes in the V-band by day 690. The overall SEDs were fitted with multiple blackbody components, in order to investigate the luminosity evolution of the supernova, and then with Monte Carlo radiative transfer models using smooth or clumpy dust distributions, in order to estimate how much new dust condensed in the ejecta. The luminosity evolution was consistent with the decay of 56Co in the ejecta up until about day 690, after which an additional emission source is required, in agreement with the findings of Kotak et al. (2009). Clumped dust density distributions consisting of 20% amorphous carbons and 80% silicates by mass were able to match the observed optical and infrared SEDs, with dust masses that increased from 8x10^{-5} Msun on day 300 to 1.5x10^{-3} Msun on day 690, still significantly lower than the values needed for core collapse supernovae to make a significant contribution to the dust enrichment of galaxies.

The extraction of chemical abundances of ionised nebulae from a limited spectral range is usually hampered by the lack of emission lines corresponding to certain ionic stages. So far, the missing emission lines have been accounted for by the ionisation correction factors (ICFs), constructed under simplistic assumptions like spherical geometry by using 1-D photoionisation modelling. In this contribution we discuss the results (Goncalves et al. 2011, in prep.) of our ongoing project to find a new set of ICFs to determine total abundances of N, O, Ne, Ar, and S, with optical spectra, in the case of non-spherical PNe. These results are based on a grid of 3-D photoionisation modelling of round, elliptical and bipolar shaped PNe, spanning the typical PN luminosities, effective temperatures and densities. We show that the additional corrections --to the widely used Kingsburgh and Barlow (1994) ICFs-- are always higher for bipolars than for ellipticals. Moreover, these additional corrections are, for bipolars, up to: 17% for oxygen, 33% for nitrogen, 40% for neon, 28% for argon and 50% for sulphur. Finally, on top of the fact that corrections change greatly with shape, they vary also greatly with the central star temperature, while the luminosity is a less important parameter.

We present Herschel PACS and SPIRE images of the dust shells around the planetary nebulae NGC 650, NGC 6853, and NGC 6720, as well as images showing the dust temperature in their shells. The latter shows a rich structure, which indicates that internal extinction in the UV is important despite the highly evolved status of the nebulae.

We have obtained spectroscopic observations of the nebulae around seven Wolf-Rayet (WR) stars, four in the Milky Way (WR 8, 16, 18 and 40) and three in the LMC (BAT99-2, -11 and -38). They were observed using the ESO NTT with the EFOSC-2 instrument, except for one nebula, NGC 3199, which was observed using the UVES echelle spectrometer on the VLT. The aims of these observations were to (a) quantify the degree of chemical enrichment in WR nebulae which had previously been suggested to have abundances reflecting nucleosynthetic processing; and (b) to attempt to detect the far-red lines of neutral carbon (e.g. [C I] 9850A) in the nebulae around WC stars. Nebular densities, temperatures and elemental abundances were derived using standard emission line diagnostics. Our high spectral resolution UVES dataset for NGC 3199 allowed the determination of line broadening temperatures using lines from several different species. These showed consistent patterns and fair agreement with the nebular temperatures derived from diagnostic line ratios. Amongst the Galactic WR nebulae, the nebulae around WR 8, 16 and 40 were found to be strongly nitrogen-enriched while NGC 3199 displayed abundances which were similar to those of Galactic H II regions. Amongst the LMC WR nebulae, none displayed an N/O ratio that significantly exceeded LMC H II region values. The far-red [C I] lines were detected from NGC 3199, with the strength of the 8727A line indicating that at least part of it was produced by C+ recombination rather than collisional excitation. The line widths of the [C I] 9824 and 9850A lines were however the same as those of collisionally excited nebular lines, indicating that those lines originated largely from ionized regions within the nebula.

We present a 3D photoionization model of the PN NGC 6302, one of the most complex objects of its kind. Our Mocassin model is composed of an extremely dense circumstellar disk and a large pair of diffuse bipolar lobes, a combination necessary to reproduce the observed spectrum. The masses of these components gives a total nebular mass of 4.7Mo. Discrepancies between our model fit and the observations are attributed to complex density inhomogeneities in the nebula. The potential to resolve such discrepancies with more complex models is confirmed by a range of models introducing small-scale structures. Compared to solar abundances He is enhanced by 50%, C is slightly subsolar, O is solar, and N is enhanced by a factor of 6. These imply a significant 3rd dredge-up coupled with hot-bottom burning CN-cycle conversion of dredged-up C to N. The central star is partly obscured by the edge-on circumstellar disk and its properties are not well constrained. Emission from a number of high-ionization `coronal' lines provides constraints on the form of the high-energy ionizing flux. Using a solar abundance stellar atmosphere we are unable to fit all of the observed line fluxes, but a substantially better fit was obtained using a 220,000K H-deficient stellar atmosphere with L*=14,300 Lo. The H-deficient nature of the central star suggests it has undergone a late thermal pulse, and fits to evolutionary tracks imply a central star mass of 0.73-0.82Mo. Timescales for these tracks suggest the object left the top of the AGB ~2100 years ago, in agreement with studies of the recent mass-loss event that formed the bipolar lobes. Based on the modelled nebular and central star masses we estimate the initial mass of the central star to be 5.5Mo, in agreement with that derived from evolutionary tracks. (Abstract truncated)

We report far-infrared and submillimeter observations of Supernova 1987A, the star that exploded on February 23, 1987 in the Large Magellanic Cloud, a galaxy located 160,000 lightyears away. The observations reveal the presence of a population of cold dust grains radiating with a temperature of ~17-23 K at a rate of about 220 solar luminosity. The intensity and spectral energy distribution of the emission suggests a dust mass of ~0.4-0.7 solar mass. The radiation must originate from the SN ejecta and requires the efficient precipitation of all refractory material into dust. Our observations imply that supernovae can produce the large dust masses detected in young galaxies at very high redshifts.

SN 2010jl was an extremely bright, Type IIn SNe which showed a significant IR excess no later than 90 days after explosion. We have obtained Spitzer 3.6 and 4.5 \mum and JHK observations of SN 2010jl \sim90 days post explosion. Little to no reddening in the host galaxy indicated that the circumstellar material lost from the progenitor must lie in a torus inclined out of the plane of the sky. The likely cause of the high mid-IR flux is the reprocessing of the initial flash of the SN by pre-existing circumstellar dust. Using a 3D Monte Carlo Radiative Transfer code, we have estimated that between 0.03-0.35 Msun of dust exists in a circumstellar torus around the SN located 6 \times 10 ^17 cm away from the SN and inclined between 60-80\cdot to the plane of the sky. On day 90, we are only seeing the illumination of approximately 5% of this torus, and expect to see an elevated IR flux from this material up until day \sim 450. It is likely this dust was created in an LBV-like mass loss event of more than 3 Msun, which is large but consistent with other LBV progenitors such as {\eta} Carinae.

Thermal images of cold dust in the Central Molecular Zone of the Milky Way, obtained with the far-infrared cameras on-board the Herschel satellite, reveal a 3x10^7 solar masses ring of dense and cold clouds orbiting the Galactic Center. Using a simple toy-model, an elliptical shape having semi-major axes of 100 and 60 parsecs is deduced. The major axis of this 100-pc ring is inclined by about 40 degrees with respect to the plane-of-the-sky and is oriented perpendicular to the major axes of the Galactic Bar. The 100-pc ring appears to trace the system of stable x_2 orbits predicted for the barred Galactic potential. Sgr A* is displaced with respect to the geometrical center of symmetry of the ring. The ring is twisted and its morphology suggests a flattening-ratio of 2 for the Galactic potential, which is in good agreement with the bulge flattening ratio derived from the 2MASS data.

In 2001, the discovery of circumstellar water vapour around the ageing carbon star IRC+10216 was announced. This detection challenged the current understanding of chemistry in old stars, since water vapour was predicted to be absent in carbon-rich stars. Several explanations for the occurrence of water vapour were postulated, including the vaporization of icy bodies (comets or dwarf planets) in orbit around the star, grain surface reactions, and photochemistry in the outer circumstellar envelope. However, the only water line detected so far from one carbon-rich evolved star can not discriminate, by itself, between the different mechanisms proposed. Here we report on the detection by the Herschel satellite of dozens of water vapour lines in the far-infrared and sub-millimetre spectrum of IRC+10216, including some high-excitation lines with energies corresponding to ~1000 K. The emission of these high-excitation water lines can only be explained if water vapour is present in the warm inner region of the envelope. A plausible explanation for the formation of warm water vapour appears to be the penetration of ultraviolet (UV) photons deep into a clumpy circumstellar envelope. This mechanism triggers also the formation of other molecules such as ammonia, whose observed abundances are much higher than hitherto predicted.

Post-AGB stars are key objects for the study of the dramatic morphological changes of low- to intermediate-mass stars on their evolution from the Asymptotic Giant Branch (AGB) towards the Planetary Nebula stage. There is growing evidences that binary interaction processes may very well have a determining role in the shaping process of many objects, but so far direct evidence is still weak. We aim at a systematic study of the dust distribution around a large sample of Post-AGB stars as a probe of the symmetry breaking in the nebulae around these systems. We used imaging in the mid-infrared to study the inner part of these evolved stars to probe direct emission from dusty structures in the core of Post-AGB stars in order to better understand their shaping mechanisms. We imaged a sample of 93 evolved stars and nebulae in the mid-infrared using VISIR/VLT, T-Recs/Gemini South and Michelle/Gemini North. We found that all the the Proto-Planetary Nebulae we resolved show a clear departure from spherical symmetry. 59 out of the 93 observed targets appear to be non resolved. The resolved targets can be divided in two categories. The nebulae with a dense central core, that are either bipolar and multipolar. The nebulae with no central core have an elliptical morphology.The dense central torus observed likely host binary systems which triggered fast outflows that shaped the nebulae.

SN 2007it is a bright, Type IIP supernova which shows indications of both pre-existing and newly formed dust. The visible photometry shows a bright late-time luminosity, powered by the 0.09 M$_{\sun}$ of $^{56}$Ni present in the ejecta. There is also a sudden drop in optical brightness after day 339, and a corresponding brightening in the IR due to new dust forming in the ejecta. CO and SiO emission, generally thought to be precursors to dust formation, may have been detected in the mid-IR photometry of SN 2007it. The optical spectra show stronger than average [O I] emission lines and weaker than average [Ca II] lines, which may indicate a 16 - 27 M$_{\sun}$ progenitor, on the higher end of expected Type IIP masses. Multi-component [O I] lines are also seen in the optical spectra, most likely caused by an asymmetric blob or a torus of oxygen core material being ejected during the SN explosion. Interaction with circumstellar material prior to day 540 may have created a cool dense shell between the forward and reverse shocks where new dust is condensing. At late times there is also a flattening of the visible lightcurve as the ejecta luminosity fades and a surrounding light echo becomes visible. Radiative transfer models of SN 2007it SEDs indicate that up to 10$^{-4}$ M$_{\sun}$ of new dust has formed in the ejecta, which is consistent with the amount of dust formed in other core collapse supernovae.

MESS (Mass-loss of Evolved StarS) is a Guaranteed Time Key Program that uses the PACS and SPIRE instruments on board the Herschel Space Observatory to observe a representative sample of evolved stars, that include asymptotic giant branch (AGB) and post-AGB stars, planetary nebulae and red supergiants, as well as luminous blue variables, Wolf-Rayet stars and supernova remnants. In total, of order 150 objects are observed in imaging and about 50 objects in spectroscopy. This paper describes the target selection and target list, and the observing strategy. Key science projects are described, and illustrated using results obtained during Herschel's science demonstration phase. Aperture photometry is given for the 70 AGB and post-AGB stars observed up to October 17, 2010, which constitutes the largest single uniform database of far-IR and sub-mm fluxes for late-type stars.

Galactic cirrus emission at far-infrared wavelengths affects many extragalactic observations. Separating this emission from that associated with extragalactic objects is both important and difficult. In this paper we discuss a particular case, the M81 group, and the identification of diffuse structures prominent in the infrared, but also detected at optical wavelengths. The origin of these structures has previously been controversial, ranging from them being the result of a past interaction between M81 and M82 or due to more local Galactic emission. We show that over of order a few arcminute scales the far-infrared (Herschel 250 &\mu&m) emission correlates spatially very well with a particular narrow velocity (2-3 km/s) component of the Galactic HI. We find no evidence that any of the far-infrared emission associated with these features actually originates in the M81 group. Thus we infer that the associated diffuse optical emission must be due to galactic light back scattered off dust in our galaxy. Ultra-violet observations pick out young stellar associations around M81, but no detectable far-infrared emission. We consider in detail one of the Galactic cirrus features, finding that the far-infrared HI relation breaks down below arc minute scales and that at smaller scales there can be quite large dust temperature variations.

In this paper we will discuss the images of Planetary Nebulae that have recently been obtained with PACS and SPIRE on board the Herschel satellite. This comprises results for NGC 650 (the little Dumbbell nebula), NGC 6853 (the Dumbbell nebula), and NGC 7293 (the Helix nebula).

IPHASXJ194359.5+170901 is a new high-excitation planetary nebula with remarkable characteristics. It consists of a knotty ring expanding at a speed of 28 km/s, and a fast collimated outflow in the form of faint lobes and caps along the direction perpendicular to the ring. The expansion speed of the polar caps is 100 km/s, and their kinematical age is twice as large as the age of the ring. Time-resolved photometry of the central star of IPHASXJ194359.5+170901 reveals a sinusoidal modulation with a period of 1.16 days. This is interpreted as evidence for binarity of the central star, the brightness variations being related to the orbital motion of an irradiated companion. This is supported by the spectrum of the central star in the visible range, which appears to be dominated by emission from the irradiated zone, consisting of a warm (6000-7000 K) continuum, narrow C III, C IV, and N III emission lines, and broader lines from a flat H I Balmer sequence in emission. IPHASXJ194359.5+170901 helps to clarify the role of (close) binaries in the formation and shaping of planetary nebulae. The output of the common-envelope evolution of the system is a strongly flattened circumstellar mass deposition, a feature that seems to be distinctive of this kind of binary system. Also, IPHASXJ194359.5+170901 is among the first post-CE PNe for which the existence of a high-velocity polar outflow has been demonstrated. Its kinematical age might indicate that the polar outflow is formed before the common-envelope phase. This points to mass transfer onto the secondary as the origin, but alternative explanations are also considered.

We present the first results of a high-spectral-resolution survey of the carbon-rich evolved star IRC+10216 that was carried out with the HIFI spectrometer onboard Herschel. This survey covers all HIFI bands, with a spectral range from 488 to 1901GHz. In this letter we focus on the band-1b spectrum, in a spectral range 554.5-636.5GHz, where we identified 130 spectral features with intensities above 0.03 K and a signal-to-noise ratio >5. Detected lines arise from HCN, SiO, SiS, CS, CO, metal-bearing species and, surprisingly, silicon dicarbide (SiC2). We identified 55 SiC2 transitions involving energy levels between 300 and 900 K. By analysing these rotational lines, we conclude that SiC2 is produced in the inner dust formation zone, with an abundance of ~2x10^-7 relative to molecular hydrogen. These SiC2 lines have been observed for the first time in space and have been used to derive an SiC2 rotational temperature of ~204 K and a source-averaged column density of ~6.4x10^15 cm^-2. Furthermore, the high quality of the HIFI data set was used to improve the spectroscopic rotational constants of SiC2.

The determination of reliable distances to Planetary Nebulae (PNe) is one of the major limitations in the study of this class of objects in the Galaxy. The availability of new photometric surveys such as IPHAS covering large portions of the sky gives us the opportunity to apply the "extinction method" to determine distances of a large number of objects. The technique is applied to a sample of 137 PNe located between -5 and 5 degrees in Galactic latitude, and between 29.52 and 215.49 degrees in longitude. The characteristics of the distance-extinction method and the main sources of errors are carefully discussed. The data on the extinction of the PNe available in the literature, complemented by new observations, allow us to determine extinction distances for 70 PNe. A comparison with statistical distance scales from different authors is presented.

We present time-resolved optical spectroscopy of V458 Vulpeculae (Nova Vul 2007 No. 1) spread over a period of 15 months starting 301 days after its discovery. Our data reveal radial velocity variations in the HeII {\lambda}5412 and HeII {\lambda}4686 emission lines. A period analysis of the radial velocity curves resulted in a period of 98.09647 \pm 0.00025 min (0.06812255 \pm 0.00000017 d) which we identify with the orbital period of the binary system. V458 Vul is therefore the planetary nebula central binary star with the shortest period known. We explore the possibility of the system being composed of a relatively massive white dwarf (M1 \gsim 1.0 Msun) accreting matter from a post-asymptotic giant branch star which produced the planetary nebula observed. In this scenario, the central binary system therefore underwent two common-envelope episodes. A combination of previous photoionisation modelling of the nebular spectra, post-asymptotic giant branch evolutionary tracks and the orbital period favour a mass of M2 \sim 0.6 Msun for the donor star. Therefore, the total mass of the system may exceed the Chandrasekhar mass, which makes V458 Vul a Type Ia supernova progenitor candidate.

Recent large scale Galactic Plane H$\alpha$ surveys allow a re-examination of the environs of Wolf-Rayet (WR) stars for the presence of a circumstellar nebula. Using the morphologies of WR nebulae known to be composed of stellar ejecta as a guide, we constructed ejecta nebula criteria similar to those of Chu (1991) and searched for likely WR ejecta nebula in the SHS H$\alpha$ survey. A new Wolf-Rayet ejecta nebula around WR 8 is found and its morphology discussed. The fraction of WR stars with ejecta type nebulae is roughly consistent between the MilkyWay (MW) and LMC at around 5-6%, with the MW sample dominated by nitrogen rich WR central stars (WN type) and the LMC stars having a higher proportion of carbon rich WR central stars (WC type). We compare our results with those of previous surveys, including those of Marston (1997) and Miller & Chu (1993), and find broad consistency. We investigate several trends in the sample: most of the clear examples of ejecta nebulae have WNh central stars; and very few ejecta nebulae have binary central stars. Finally, the possibly unique evolutionary status of the nebula around the binary star WR 71 is explored.

SPIRE, the Spectral and Photometric Imaging Receiver, is the Herschel Space Observatory's submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 {\mu}m, and an imaging Fourier transform spectrometer (FTS) covering 194-671 {\mu}m (447-1550 GHz). In this paper we describe the initial approach taken to the absolute calibration of the SPIRE instrument using a combination of the emission from the Herschel telescope itself and the modelled continuum emission from solar system objects and other astronomical targets. We present the photometric, spectroscopic and spatial accuracy that is obtainable in data processed through the "standard" pipelines. The overall photometric accuracy at this stage of the mission is estimated as 15% for the photometer and between 15 and 50% for the spectrometer. However, there remain issues with the photometric accuracy of the spectra of low flux sources in the longest wavelength part of the SPIRE spectrometer band. The spectrometer wavelength accuracy is determined to be better than 1/10th of the line FWHM. The astrometric accuracy in SPIRE maps is found to be 2 arcsec when the latest calibration data are used. The photometric calibration of the SPIRE instrument is currently determined by a combination of uncertainties in the model spectra of the astronomical standards and the data processing methods employed for map and spectrum calibration. Improvements in processing techniques and a better understanding of the instrument performance will lead to the final calibration accuracy of SPIRE being determined only by uncertainties in the models of astronomical standards.

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

The interstellar medium is enriched primarily by matter ejected from evolved low and intermediate mass stars. The outflows from these stars create a circumstellar envelope in which a rich gas-phase and dust-nucleation chemistry takes place. We observed the nearest carbon-rich evolved star, IRC+10216, using the PACS (55-210 {\mu}m) and SPIRE (194-672 {\mu}m) spectrometers on board Herschel. We find several tens of lines from SiS and SiO, including lines from the v=1 vibrational level. For SiS these transitions range up to J=124-123, corresponding to energies around 6700K, while the highest detectable transition is J=90-89 for SiO, which corresponds to an energy around 8400K. Both species trace the dust formation zone of IRC+10216, and the broad energy ranges involved in their detected transitions permit us to derive the physical properties of the gas and the particular zone in which each species has been formed. This allows us to check the accuracy of chemical thermodynamical equilibrium models and the suggested depletion of SiS and SiO due to accretion onto dust grains.

We report on the detection of anhydrous hydrochloric acid (hydrogen chlorine, HCl) in the carbon-rich star IRC+10216 using the spectroscopic facilities onboard the Herschel satellite. Lines from J=1-0 up to J=7-6 have been detected. From the observed intensities, we conclude that HCl is produced in the innermost layers of the circumstellar envelope with an abundance relative to H2 of 5x10^-8 and extends until the molecules reach its photodissociation zone. Upper limits to the column densities of AlH, MgH, CaH, CuH, KH, NaH, FeH, and other diatomic hydrides have also been obtained.

We obtained Herschel PACS and SPIRE images of the thermal emission of the debris disk around the A5V star {\beta} Pic. The disk is well resolved in the PACS filters at 70, 100, and 160 {\mu}m. The surface brightness profiles between 70 and 160 {\mu}m show no significant asymmetries along the disk, and are compatible with 90% of the emission between 70 and 160 {\mu}m originating in a region closer than 200 AU to the star. Although only marginally resolving the debris disk, the maps obtained in the SPIRE 250 - 500 {\mu}m filters provide full-disk photometry, completing the SED over a few octaves in wavelength that had been previously inaccessible. The small far-infrared spectral index ({\beta} = 0.34) indicates that the grain size distribution in the inner disk (<200AU) is inconsistent with a local collisional equilibrium. The size distribution is either modified by non-equilibrium effects, or exhibits a wavy pattern, caused by an under-abundance of impactors which have been removed by radiation pressure.

We present five band imaging of the Vega debris disc obtained using the Herschel Space Observatory. These data span a wavelength range of 70-500 um with full-width half-maximum angular resolutions of 5.6-36.9". The disc is well resolved in all bands, with the ring structure visible at 70 and 160 um. Radial profiles of the disc surface brightness are produced, and a disc radius of 11" (~ 85 AU) is determined. The disc is seen to have a smooth structure thoughout the entire wavelength range, suggesting that the disc is in a steady state, rather than being an ephemeral structure caused by the recent collision of two large planetesimals.

We present far-infrared and submillimetre spectra of three carbon-rich evolved objects, AFGL 2688, AFGL 618 and NGC 7027. The spectra were obtained with the SPIRE Fourier transform spectrometer on board the Herschel Space Observatory, and cover wavelengths from 195-670 um, a region of the electromagnetic spectrum hitherto difficult to study in detail. The far infrared spectra of these objects are rich and complex, and we measure over 150 lines in each object. Lines due to 18 different species are detected. We determine physical conditions from observations of the rotational lines of several molecules, and present initial large velocity gradient models for AFGL 618. We detect water in AFGL 2688 for the first time, and confirm its presence in AFGL 618 in both ortho and para forms. In addition, we report the detection of the J=1-0 line of CH+ in NGC 7027.

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

With a luminosity > 10^5 Lsun and a mass-loss rate of about 2.10-4 Msun/yr, 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 VYCMa over the full wavelength range offered by the PACS and SPIRE instruments of Herschel, i.e. 55 to 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.

Using the 3.5-m Herschel Space Observatory, imaging photometry of Cas A has been obtained in six bands between 70 um and 500 um with the PACS and SPIRE instruments, with angular resolutions ranging from 6 to 37". In the outer regions of the remnant the 70-um PACS image resembles the 24-um image Spitzer image, with the emission attributed to the same warm dust component, located in the reverse shock region. At longer wavelengths, the three SPIRE bands are increasingly dominated by emission from cold interstellar dust knots and filaments, particularly across the central, western and southern parts of the remnant. Nonthermal emission from the northern part of the remnant becomes prominent at 500 um. We have estimated and subtracted the contributions from the nonthermal, warm dust and cold interstellar dust components. We confirm and resolve for the first time a cool (~35 K) dust component, emitting at 70-160 um, that is located interior to the reverse shock region, with an estimated mass of 0.075 Msun.

The standard method of mapping the interstellar medium in a galaxy, by observing the molecular gas in the CO 1-0 line and the atomic gas in the 21-cm line, is largely limited with current telescopes to galaxies in the nearby universe. In this letter, we use SPIRE observations of the galaxies M99 and M100 to explore the alternative approach of mapping the interstellar medium using the continuum emission from the dust. We have compared the methods by measuring the relationship between the star-formation rate and the surface density of gas in the galaxies. We find the two methods give relationships with a similar dispersion, confirming that observing the continuum emission from the dust is a promising method of mapping the interstellar medium in galaxies.

We present Herschel SPIRE and PACS photometeric observations of the low metallicity (Z ~ 0.35 solar) nearby dwarf galaxy, NGC 1705, in six wavelength bands as part of the Dwarf Galaxy Survey guaranteed time Herschel Key Program. We confirm the presence of two dominant circumnuclear IR-bright regions surrounding the central super star cluster that had been previously noted at mid-IR wavelengths and in the sub-mm by LABOCA. On constructing a global spectral energy distribution using the SPIRE and PACS photometry, in conjunction with archival IR measurements, we note the presence of an excess at sub-mm wavelengths. This excess suggests the presence of a significant cold dust component within NGC 1705 and was modeled as an additional cold component in the SED. Although alternative explanations for the sub-mm excess beyond 350 microns, such as changes to the dust emissivity cannot be ruled out, the most likely explanation for the observed submillimetre excess is that of an additional cold dust component.

We present the observations of the starburst galaxy M82 taken with the Herschel SPIRE Fourier Transform Spectrometer. The spectrum (194-671 {\mu}m) shows a prominent CO rotational ladder from J = 4-3 to 13-12 emitted by the central region of M82. The fundamental properties of the gas are well constrained by the high J lines observed for the first time. Radiative transfer modeling of these high-S/N 12CO and 13CO lines strongly indicates a very warm molecular gas component at ~500 K and pressure of ~3x10^6 K cm^-3, in good agreement with the H_2 rotational lines measurements from Spitzer and ISO. We suggest that this warm gas is heated by dissipation of turbulence in the interstellar medium (ISM) rather than X-rays or UV flux from the straburst. This paper illustrates the promise of the SPIRE FTS for the study of the ISM of nearby galaxies.

Herschel PACS and SPIRE images have been obtained over a 30'x30' area around the well-known carbon star CW Leo (IRC +10 216). An extended structure is found in an incomplete arc of ~22' diameter, which is cospatial with the termination shock due to interaction with the interstellar medium (ISM) as defined by Sahai & Chronopoulos from ultraviolet GALEX images. Fluxes are derived in the 70, 160, 250, 350, and 550 um bands in the region where the interaction with the ISM takes place, and this can be fitted with a modified black body with a temperature of 25+-3 K. Using the published proper motion and radial velocity for the star, we derive a heliocentric space motion of 25.1 km/s. Using the PACS and SPIRE data and the analytical formula of the bow shock structure, we infer a de-projected standoff distance of the bow shock of R0 = (8.0+-0.3)x10^17 cm. We also derive a relative velocity of the star with respect to the ISM of (106.6+-8.7)/sqrt(n_ISM) km/s, where n_ISM is the number density of the local ISM.

We present the first Herschel PACS and SPIRE images of the low-metallicity galaxy NGC6822 observed from 70 to 500 mu and clearly resolve the HII regions with PACS and SPIRE. We find that the ratio 250/500 is dependent on the 24 mu surface brightness in NGC6822, which would locally link the heating processes of the coldest phases of dust in the ISM to the star formation activity. We model the SEDs of some regions HII regions and less active regions across the galaxy and find that the SEDs of HII regions show warmer ranges of dust temperatures. We derive very high dust masses when graphite is used in our model to describe carbon dust. Using amorphous carbon, instead, requires less dust mass to account for submm emission due to its lower emissivity properties. This indicates that SED models including Herschel constraints may require different dust properties than commonly used.

With appropriate spatial resolution, images of spiral galaxies in thermal infrared (~10 micron and beyond) often reveal a bright central component, distinct from the stellar bulge, superimposed on a disk with prominent spiral arms. ISO and Spitzer studies have shown that much of the scatter in the mid-infrared colors of spiral galaxies is related to changes in the relative importance of these two components, rather than to other modifications, such as the morphological type or star formation rate, that affect the properties of the galaxy as a whole. With the Herschel imaging capability from 70 to 500 micron, we revisit this two-component approach at longer wavelengths, to see if it still provides a working description of the brightness distribution of galaxies, and to determine its implications on the interpretation of global far-infrared properties of galaxies.

We use Herschel Space Observatory data to place observational constraints on the peak and Rayleigh-Jeans slope of dust emission observed at 70-500 microns in the nearby spiral galaxy M81. We find that the ratios of wave bands between 160 and 500 microns are primarily dependent on radius but that the ratio of 70 to 160 micron emission shows no clear dependence on surface brightness or radius. These results along with analyses of the spectral energy distributions imply that the 160-500 micron emission traces 15-30 K dust heated by evolved stars in the bulge and disc whereas the 70 micron emission includes dust heated by the active galactic nucleus and young stars in star forming regions.

We present Herschel-SPIRE observations at 250-500um of the giant elliptical galaxy M86 and examine the distribution of the resolved cold dust emission and its relation with other galactic tracers. The SPIRE images reveal three dust components: emission from the central region; a dust lane extending north-south; and a bright emission feature 10kpc to the south-east. We estimate that approximately 10^6 solar masses of dust is spatially coincident with atomic and ionized hydrogen, originating from stripped material from the nearby spiral NGC4438 due to recent tidal interactions with M86. The gas-to-dust ratio of the cold gas component ranges from ~20-80. We discuss the different heating mechanisms for the dust features.

Herschel PACS and SPIRE images have been obtained of NGC 6720 (the Ring Nebula). This is an evolved planetary nebula with a central star that is currently on the cooling track, due to which the outer parts of the nebula are recombining. From the PACS and SPIRE images we conclude that there is a striking resemblance between the dust distribution and the H2 emission, which appears to be observational evidence that H2 forms on grain surfaces. We have developed a photoionization model of the nebula with the Cloudy code which we used to determine the physical conditions of the dust and investigate possible formation scenarios for the H2. We conclude that the most plausible scenario is that the H2 resides in high density knots which were formed after the recombination of the gas started when the central star entered the cooling track. Hydrodynamical instabilities due to the unusually low temperature of the recombining gas are proposed as a mechanism for forming the knots. H2 formation in the knots is expected to be substantial after the central star underwent a strong drop in luminosity about one to two thousand years ago, and may still be ongoing at this moment, depending on the density of the knots and the properties of the grains in the knots.

M82 is a unique representative of a whole class of galaxies, starbursts with superwinds, in the Very Nearby Galaxy Survey with Herschel. In addition, its interaction with the M81 group has stripped a significant portion of its interstellar medium from its disk. SPIRE maps now afford better characterization of the far-infrared emission from cool dust outside the disk, and sketch a far more complete picture of its mass distribution and energetics than previously possible. They show emission coincident in projection with the starburst wind and in a large halo, much more extended than the PAH band emission seen with Spitzer. Some complex substructures coincide with the brightest PAH filaments, and others with tidal streams seen in atomic hydrogen. We subtract the far-infrared emission of the starburst and underlying disk from the maps, and derive spatially-resolved far-infrared colors for the wind and halo. We interpret the results in terms of dust mass, dust temperature, and global physical conditions. In particular, we examine variations in the dust physical properties as a function of distance from the center and the wind polar axis, and conclude that more than two thirds of the extraplanar dust has been removed by tidal interaction, and not entrained by the starburst wind.

We present a detailed analysis of the far-IR properties of the bright, lensed, z = 2.3, SMG, SMM J2135-0102, using new observations with Herschel, SCUBA-2 and the VLA. These data allow us to constrain the galaxy's SED and show that it has an intrinsic rest-frame 8-1000um luminosity, L(bol), of (2.3 +/- 0.2) x 10^12 L(sun) and a likely SFR of ~400 M(sun)/yr. The galaxy sits on the far-IR/radio correlation for far-IR-selected galaxies. At ~>70um, the SED can be described adequately by dust components with T(d) ~ 30 and 60K. Using SPIRE's Fourier Transform Spectrometer we report a detection of the [CII] 158um cooling line. If the [CII], CO and far-IR continuum arise in photo-dissociation regions, we derive a characteristic gas density, n ~ 10^3 cm^-3, and a far-UV radiation field, G_0, 10^3x stronger than the Milky Way. L([CII])/L(bol) is significantly higher than in local ULIRGs but similar to the values found in local star-forming galaxies and starburst nuclei. This is consistent with SMM J2135-0102 being powered by starburst clumps distributed across ~2 kpc, evidence that SMGs are not simply scaled-up ULIRGs. Our results show that SPIRE's FTS has the ability to measure the redshifts of distant, obscured galaxies via the blind detection of atomic cooling lines, but it will not be competitive with ground-based CO-line searches. It will, however, allow detailed study of the integrated properties of high-redshift galaxies, as well as the chemistry of their ISM, once more suitably bright candidates have been found.

SN 2007od exhibits characteristics that have rarely been seen in a Type IIP supernova (SN). Optical V band photometry reveals a very steep brightness decline between the plateau and nebular phases of ~4.5 mag, likely due to SN 2007od containing a low mass of 56Ni. The optical spectra show an evolution from normal Type IIP with broad Halpha emission, to a complex, four component Halpha emission profile exhibiting asymmetries caused by dust extinction after day 232. This is similar to the spectral evolution of the Type IIn SN 1998S, although no early-time narrow (~200 km s-1) Halpha component was present in SN 2007od. In both SNe, the intermediate-width Halpha emission components are thought to arise in the interaction between the ejecta and its circumstellar medium (CSM). SN 2007od also shows a mid-IR excess due to new dust. The evolution of the Halpha profile and the presence of the mid-IR excess provide strong evidence that SN 2007od formed new dust before day 232. Late-time observations reveal a flattening of the visible lightcurve. This flattening is a strong indication of the presence of a light echo, which likely accounts for much of the broad, underlying Halpha component seen at late-times. We believe the multi-peaked Halpha emission is consistent with the interaction of the ejecta with a circumstellar ring or torus (for the inner components at \pm1500 km s-1), and a single blob or cloud of circumstellar material out of the plane of the CSM ring (for the outer component at -5000 km s-1). The most probable location for the formation of new dust is in the cool dense shell created by the interaction between the expanding ejecta and its CSM. Monte Carlo radiative transfer modeling of the dust emission from SN 2007od implies that up to 4x 10-4Msun of new dust has formed. This is similar to the amounts of dust formed in other CCSNe such as SNe 1999em, 2004et, and 2006jc.

We present chemical and radiative transfer models for the many far-IR ortho- and para-H2O lines that were observed from the Orion-KL region in high resolution Fabry-Perot (FP) mode by the Long Wavelength Spectrometer (LWS) on board the Infrared Space Observatory (ISO). The chemistry of the region was first studied by simulating the conditions in the different known components of Orion-KL: chemical models for a hot core, a plateau and a ridge were coupled with an accelerated Lambda$-iteration (ALI) radiative transfer model to predict H2O line fluxes and profiles. Our models include the first 45 energy levels of ortho- and para-H2O. We find that lines arising from energy levels below 560 K were best reproduced by a gas of density 3x10^5 cm^-3 at a temperature of 70-90 K, expanding at a velocity of 30 km s^-1 and with a H2O/H2 abundance ratio of the order of 2 - 3 x 10^-5, similar to the abundance derived by Cernicharo et al. (2006). However, the model that best reproduces the fluxes and profiles of H2O lines arising from energy levels above 560 K has a significantly higher H2O/H2 abundance, 1 - 5 x 10^-4, originating from gas of similar density, in the Plateau region, that has been heated to 300 K, relaxing to 90-100 K. We conclude that the observed water lines do not originate from high temperature shocks.

The Herschel Reference Survey is a guaranteed time Herschel key project and will be a benchmark study of dust in the nearby universe. The survey will complement a number of other Herschel key projects including large cosmological surveys that trace dust in the distant universe. We will use Herschel to produce images of a statistically-complete sample of 323 galaxies at 250, 350 and 500 micron. The sample is volume-limited, containing sources with distances between 15 and 25 Mpc and flux limits in the K-band to minimize the selection effects associated with dust and with young high-mass stars and to introduce a selection in stellar mass. The sample spans the whole range of morphological types (ellipticals to late-type spirals) and environments (from the field to the centre of the Virgo Cluster) and as such will be useful for other purposes than our own. We plan to use the survey to investigate (i) the dust content of galaxies as a function of Hubble type, stellar mass and environment, (ii) the connection between the dust content and composition and the other phases of the interstellar medium and (iii) the origin and evolution of dust in galaxies. In this paper, we describe the goals of the survey, the details of the sample and some of the auxiliary observing programs that we have started to collect complementary data. We also use the available multi-frequency data to carry out an analysis of the statistical properties of the sample.

Far-infrared [C ii] 158 micron, [O i] 145 micron and [O i] 63 micron fine structure emission line fluxes were measured from archival Infrared Space Observatory Long Wavelength Spectrometer spectra of 46 extragalactic sources, with 28 sources providing detections in all three lines. For 12 of the sources, the contribution to the [C ii] 158 micron line flux from H ii regions could be estimated from their detected [N ii] 122 micron line fluxes. The measured [C ii]/[O i] and [O i] 63/145 micron line flux ratios were compared with those from a grid of PDR models previously computed using the UCL PDR code. Persistent offsets between the observed and modelled line ratios could be partly attributed to the effects of [O i] 63 micron self-absorption. Using the SMMOL code, we calculated model [O i] line profiles and found that the strength of the [O i] 63 micron line was reduced by 20-80%, depending on the PDR parameters. We conclude that high PDR densities and radiation field strengths, coupled with the effects of [O i] 63 micron self-absorption, are likely to provide the best match to the observed line flux ratios.

During their evolution, Asymptotic Giant Branch (AGB) stars experience a high mass-loss which leads to the formation of a Circumstellar Envelope (CSE) of dust and gas. The mass-loss process is the most important phenomenon during this evolutionary stage. In order to understand it, it is important to study the physical parameters of the CSE. The emission of the CSE in the (sub)millimetre range is dominated by the dust continuum. This means that (sub)millimetre observations are a key tool in tracing the dust and improving our knowledge of the mass-loss process. In this paper we analyse new sub-millimetre observations of 9 evolved stars in order to constrain the CSE physical parameters. The data were taken by the new APEX bolometer LABoCa. The fluxes at 870 um are derived and the extended emission is investigated. We compute the spectral energy distribution (SEDs) using a 1D radiative transfer code, DUSTY which we compared to literature data. Grain properties are calculated using both spherical grains distribution and a Continuous Distribution of Ellipsoids (CDE) and a comparison between the two is drawn. Synthetic surface brightness maps have been derived from the modelling and were compared to the LABoCa brightness maps. We detected the presence of extended emission around four stars. Physical parameters of the circumstellar envelope are derived from SED modelling, such as the dust chemical composition, the dust condensation temperature and the total mass of the envelope. It proves difficult however to fit the SED and the intensity profile simultaneously.

The INT Photometric H-Alpha Survey (IPHAS) is particular effective in tracing asymptotic giant branch (AGB) stars across the Galactic disk due to its use of broad-band r' and i' filters, and the clear separation of the dwarf and giant branches in the colour-colour diagram. Here we discuss some results from our studies of the Galactic AGB population that has been revealed by IPHAS. This includes a photometric study, our spectroscopic follow-up campaign, and the photometric technique that allows us to separate O-rich, S-type and carbon stars across the entire Galactic disk based on their IPHAS photometry alone. We also discuss the potential impact of this technique on our understanding of the dredge-up process in AGB stars, Galactic chemical evolution, and the structure and properties of the Galactic disk.

SN 2008S erupted in early 2008 in the grand design spiral galaxy NGC 6946. The progenitor was detected by Prieto et al. in Spitzer Space Telescope images taken over the four years prior to the explosion, but was not detected in deep optical images, from which they inferred a self-obscured object with a mass of about 10 Msun. We obtained Spitzer observations of SN 2008S five days after its discovery, as well as coordinated Gemini and Spitzer optical and infrared observations six months after its outburst. We have constructed radiative transfer dust models for the object before and after the outburst, using the same r^-2 density distribution of pre-existing amorphous carbon grains for all epochs and taking light-travel time effects into account for the early post-outburst epoch. We rule out silicate grains as a significant component of the dust around SN 2008S. The inner radius of the dust shell moved outwards from its pre-outburst value of 85 AU to a post-outburst value of 1250 AU, attributable to grain vaporisation by the light flash from SN 2008S. Although this caused the circumstellar extinction to decrease from Av = 15 before the outburst to 0.8 after the outburst, we estimate that less than 2% of the overall circumstellar dust mass was destroyed. The total mass-loss rate from the progenitor star is estimated to have been (0.5-1.0)x10^-4 Msun yr^-1. The derived dust mass-loss rate of 5x10^-7 Msun yr^-1 implies a total dust injection into the ISM of up to 0.01 Msun over the suggested duration of the self-obscured phase. We consider the potential contribution of objects like SN 2008S to the dust enrichment of galaxies.

In a previous paper [arXiv:0712.2391], we presented the selection criteria needed to search for symbiotic stars in IPHAS, the INT Halpha survey of the Northern Galactic plane. IPHAS gives us the opportunity to make a systematic, complete search for symbiotic stars in a magnitude-limited volume. Follow-up spectroscopy at different telescopes worldwide of a sample of sixty two symbiotic star candidates is presented. Seven out of nineteen S-type candidates observed spectroscopically are confirmed to be genuine symbiotic stars. The spectral type of their red giant components, as well as reddening and distance, were computed by modelling the spectra. Only one new D-type symbiotic system, out of forty-three candidates observed, was found. This was as expected (see discussion in our paper on the selection criteria). The object shows evidence for a high density outflow expanding at a speed larger than 65 km/s. Most of the other candidates are lightly reddened classical T Tauri stars and more highly reddened young stellar objects that may be either more massive young stars of HAeBe type or classical Be stars. In addition, a few notable objects have been found, such as three new Wolf-Rayet stars and two relatively high-luminosity evolved massive stars. We also found a helium-rich source, possibly a dense ejecta hiding a WR star, which is surrounded by a large ionized nebula.

We report on optical integral field spectroscopy of two unrelated blue compact galaxies mapped with the 13 x 13 arcsec^2 VIMOS integral field unit at a resolution of 0.33 x 0.33 arcsec^2. Continuum and background subtracted emission line maps in the light of [O III] 5007, H-alpha, and [N II] 6584 are presented. Both galaxies display signs of ongoing perturbation and/or interaction. UM 420 is resolved for the first time to be a merging system composed of two starbursting components with an 'arm-like' structure associated with the largest component. UM 462 which is a disrupted system of irregular morphology is resolved into at least four starbursting regions. Maps of the H-alpha radial velocity and FWHM are discussed. No underlying broad line region was detected from either galaxy as the emission lines are well-fitted with single Gaussian profiles only. Electron temperatures and densities as well as the abundances of helium, oxygen, nitrogen, and sulphur were computed from spectra integrated over the whole galaxies and for each area of recent star formation. Maps of the O/H ratio are presented: these galaxies show oxygen abundances that are ~20 per cent solar. No evidence of substantial abundance variations across the galaxies that would point to significant nitrogen or oxygen self-enrichment is found (<0.2 dex limit). Contrary to previous observations, this analysis does not support the classification of these BCGs as Wolf-Rayet galaxies as the characteristic broad emission line features have not been detected in our spectra. Baldwin-Phillips-Terlevich emission line ratio diagrams which were constructed on a pixel by pixel basis indicate that the optical spectra of these systems are predominantly excited by stellar photoionization.

We present a library of 139 near-IR spectra of cool asymptotic giant branch stars that will be useful for comparison with theoretical model atmosphere calculations and for modeling the integrated emission from intermediate-age stellar populations. The source list was selected from the `extremely red' region of the INT Photometric H Alpha Survey (IPHAS) colour-colour plane that is overwhelmingly dominated by very late-type stars. The spectral library also includes a large fraction of S-type and carbon stars. We present a number of spectral classification sequences highlighting the various molecular features identified and discuss a number of rare features with uncertain identifications in the literature. With its focus on particularly cool photospheres this catalogue serves as a companion to recent spectroscopic atlases of MK standards in the near-IR. Finally the relationship between IPHAS (r'-i') and (r'-H Alpha) colours and spectroscopically determined properties is discussed and a strong correlation between (r'-H Alpha) colour and the C/O abundance index for S-type and carbon stars is noted. This relation has the potential to separate O-rich, S-type and carbon stars in the Galaxy based on their photometry alone.

The UV-Excess Survey of the Northern Galactic Plane images a 10x185 degree wide band, centered on the Galactic Equator using the 2.5m Isaac Newton Telescope in four bands (U,g,r,HeI5875) down to ~21st-22nd magnitude (~20th in HeI5875). The setup and data reduction procedures are described. Simulations of the colours of main-sequence stars, giant, supergiants, DA and DB white dwarfs and AM CVn stars are made, including the effects of reddening. A first look at the data of the survey (currently 30% complete) is given.

A study of the blue compact dwarf (BCD) galaxy Mrk 996 based on high resolution optical VLT VIMOS integral field unit spectroscopy is presented. Mrk 996 displays multi-component line emission, with most line profiles consisting of a narrow, central Gaussian with an underlying broad component. The broad HI Balmer component splits into two separate broad components inside a 1".5 radius from the nucleus; these are attributed to a two-armed mini-spiral. The rotation curve of Mrk 996 derived from the H\alpha narrow component yields a total mass of 5x10^8 Msol within a radius of 3 kpc. The high excitation energy, high critical density [O III] 4363 and [N II] 5755 lines are only detected from the inner region and exist purely in broad component form, implying unusual excitation conditions. Surface brightness, radial velocity, and FWHM maps for several emission components are presented. A separate physical analysis of the broad and narrow emission line regions is undertaken. We derive an upper limit of 10,000 K for the electron temperature of the narrow line gas, together with an electron density of 170 cm^-3, typical of normal H II regions. For the broad line component, we estimate a temperature of 11,000 K and an electron density of 10^7 cm^-3. The broad line emission regions show a N/H enrichment factor of ~20 relative to the narrow line regions, but no He/H, O/H, S/H, or Ar/H enrichment is inferred. The dominant line excitation mechanism is photoionisation by the ~3000 WR stars and ~150,000 O-type stars estimated to be present in the core. This is indeed a peculiar BCD, with extremely dense zones of gas in the core, through which stellar outflows and possible shock fronts permeate contributing to the excitation of the broad line emission. [Abridged]

We report on an analysis of the gas and dust budget in the the interstellar medium (ISM) of the Large Magellanic Cloud (LMC). Recent observations from the Spitzer Space Telescope enable us to study the mid-infrared dust excess of asymptotic giant branch (AGB) stars in the LMC. This is the first time we can quantitatively assess the gas and dust input from AGB stars over a complete galaxy, fully based on observations. The integrated mass-loss rate over all intermediate and high mass-loss rate carbon-rich AGB candidates in the LMC is 8.5x10^-3 solar mass per year, up to 2.1x10^-2 solar mass per year. This number could be increased up to 2.7x10^-2 solar mass per year, if oxygen-rich stars are included. This is overall consistent with theoretical expectations, considering the star formation rate when these low- and intermediate-mass stars where formed, and the initial mass functions. AGB stars are one of the most important gas sources in the LMC, with supernovae (SNe), which produces about 2-4x10^-2 solar mass per year. At the moment, the star formation rate exceeds the gas feedback from AGB stars and SNe in the LMC, and the current star formation depends on gas already present in the ISM. This suggests that as the gas in the ISM is exhausted, the star formation rate will eventually decline in the LMC, unless gas is supplied externally. Our estimates suggest `a missing dust-mass problem' in the LMC, which is similarly found in high-z galaxies: the accumulated dust mass from AGB stars and possibly SNe over the dust life time (400--800 Myrs) is significant less than the dust mass in the ISM. Another dust source is required, possibly related to star-forming regions.

We have identified 17 A-type stars in the Galactic Plane that have mid-IR excesses at 8 micron. From the observed colors in the (r'-H_alpha)-(r'-i') plane, we first identified 23050 early A-type main sequence (MS) star candidates in the Isaac Newton Photometric H-Alpha Survey (IPHAS) point source database that are located in Spitzer GLIMPSE Galactic Plane fields. Imposing the requirement that they be detected in all seven 2MASS and IRAC bands led to a sample of 2692 candidate A-type stars with fully sampled 0.6 to 8 micron SEDs. Optical classification spectra of 18 of the IPHAS candidate A-type MS stars showed that all but one could be well fitted using main sequence A-type templates, with the other being an A-type supergiant. Out of the 2692 A-type candidates 17 (0.6%) were found to have 8-micron excesses above the expected photospheric values. Taking into account non-A-Type contamination estimates, the 8-micron excess fraction is adjusted to ~0.7%. The distances to these sources range from 0.7-2.5 kpc. Only 10 out of the 17 excess stars had been covered by Spitzer MIPSGAL survey fields, of which 5 had detectable excesses at 24 micron. For sources with excesses detected in at least two mid-IR wavelength bands, blackbody fits to the excess SEDs yielded temperatures ranging from 270 to 650 K, and bolometric luminosity ratios L_IR/L* from 2.2x10^{-3}-1.9x10^{-2}, with a mean value of 7.9x10^{-3} (these bolometric luminosities are lower limits as cold dust is not detectable by this survey). Both the presence of mid-IR excesses and the derived bolometric luminosity ratios are consistent with many of these systems being in the planet-building transition phase between the early protoplanetary disk phase and the later debris disk phase.

We present photometric analysis and follow-up spectroscopy for a population of extremely red stellar objects extracted from the point-source catalogue of the INT Photometric H-Alpha Survey (IPHAS) of the northern galactic plane. The vast majority of these objects have no previous identification. Analysis of optical, near- and mid-infrared photometry reveals that they are mostly highly-reddened asymptotic giant branch stars, with significant levels of circumstellar material. We show that the distribution of these objects traces galactic extinction, their highly reddened colours being a product of both interstellar and circumstellar reddening. This is the first time that such a large sample of evolved low-mass stars has been detected in the visual and allows optical counterparts to be associated with sources from recent infrared surveys. Follow-up spectroscopy on some of the most interesting objects in the sample has found significant numbers of S-type stars which can be clearly separated from oxygen-rich objects in the IPHAS colour-colour diagram. We show that this is due to the positions of different molecular bands relative to the narrow-band H-alpha filter used for IPHAS observations. The IPHAS (r' - H-alpha) colour offers a valuable diagnostic for identifying S-type stars. A selection method for identifying S-type stars in the galactic plane is briefly discussed and we estimate that over a thousand new objects of this type may be discovered, potentially doubling the number of known objects in this short but important evolutionary phase.

I review the prospects for studies of the advanced evolutionary stages of low-, intermediate- and high-mass stars by the JWST and concurrent facilities, with particular emphasis on how they may help elucidate the dominant contributors to the interstellar dust component of galaxies. Observations extending from the mid-infrared to the submillimeter can help quantify the heavy element and dust species inputs to galaxies from AGB stars. JWST's MIRI mid-infrared instrument will be so sensitive that observations of the dust emission from individual intergalactic AGB stars and planetary nebulae in the Virgo Cluster will be feasible. The Herschel Space Observatory will enable the last largely unexplored spectral region, the far-IR to the submillimeter, to be surveyed for new lines and dust features, while SOFIA will cover the wavelength gap between JWST and Herschel, a spectral region containing important fine structure lines, together with key water-ice and crystalline silicate bands. Spitzer has significantly increased the number of Type II supernovae that have been surveyed for early-epoch dust formation but reliable quantification of the dust contributions from massive star supernovae of Type II, Type Ib and Type Ic to low- and high-redshift galaxies should come from JWST MIRI observations, which will be able to probe a volume over 1000 times larger than Spitzer.

Nova V458 Vul erupted on 2007 August 8th and reached a visual magnitude of 8.1 a few days later. H$\alpha$ images obtained six weeks before the outburst as part of the IPHAS galactic plane survey reveal an 18th magnitude progenitor surrounded by an extended nebula. Subsequent images and spectroscopy of the nebula reveal an inner nebular knot increasing rapidly in brightness due to flash ionization by the nova event. We derive a distance of 13 kpc based on light travel time considerations, which is supported by two other distance estimation methods. The nebula has an ionized mass of 0.2 M$_{\odot}$ and a low expansion velocity: this rules it out as ejecta from a previous nova eruption, and is consistent with it being a ~14,000 year old planetary nebula, probably the product of a prior common envelope (CE) phase of evolution of the binary system. The large derived distance means that the mass of the erupting WD component of the binary is high. We identify two possible evolutionary scenarios, in at least one of which the system is massive enough to produce a Type Ia supernova on merging.

We present the 3.5m SPICA space telescope, a proposed Japanese-led JAXA-ESA mission scheduled for launch around 2017. The actively cooled (<5 K), single aperture telescope and monolithic mirror will operate from ~3.5 to ~210 um and will provide superb sensitivity in the mid- and far-IR spectral domain (better than JWST at lambda > 18 um). SPICA is one of the few space missions selected to go to the next stage of ESA's Cosmic Vision 2015-2025 selection process. In this White Paper we present the main specifications of the three instruments currently baselined for SPICA: a mid-infrared (MIR) coronagraph (~3.5 to ~27 um) with photometric and spectral capabilities (R~200), a MIR wide-field camera and high resolution spectrometer (R~30,000), and a far-infrared (FIR ~30 to ~210 um) imaging spectrometer - SAFARI - led by a European consortium. We discuss their capabilities in the context of MIR direct observations of exo-planets (EPs) and multiband photometry/high resolution spectroscopy observations of transiting exo-planets. We conclude that SPICA will be able to characterize the atmospheres of transiting exo-planets down to the super-Earth size previously detected by ground- or space-based observatories. It will also directly detect and characterize Jupiter/Neptune-size planets orbiting at larger separation from their parent star (>5-10 AU), by performing quantitative atmospheric spectroscopy and studying proto-planetary and debris disks. In addition, SPICA will be a scientific and technological precursor for future, more ambitious, IR space missions for exo-planet direct detection as it will, for example, quantify the prevalence exo-zodiacal clouds in planetary systems and test coronographic techniques, cryogenic systems and lightweight, high quality telescopes. (abridged)

The INT/WFC Photometric H-alpha Survey of the Northern Galactic Plane (IPHAS) is an imaging survey being carried out in H-alpha, r' and i' filters, with the Wide Field Camera (WFC) on the 2.5-metre Isaac Newton Telescope (INT) to a depth of r'=20 (10 sigma). The survey is aimed at revealing large scale structure in our local galaxy, and also the properties of key early and late populations making up the Milky Way. Mapping emission line objects enables a particular focus on objects in the young and old stages of stellar evolution ranging from early T-Tauri stars to late planetary nebulae. In this paper we present the IPHAS Initial Data Release, primarily a photometric catalogue of about 200 million unique objects, coupled with associated image data covering about 1600 square degrees in three passbands. We note how access to the primary data products has been implemented through use of standard virtual observatory publishing interfaces. Simple traditional web access is provided to the main IPHAS photometric catalogue, in addition to a number of common catalogues (such as 2MASS) which are of immediate relevance. Access through the AstroGrid VO Desktop opens up the full range of analysis options, and allows full integration with the wider range of data and services available through the Virtual Observatory. The IDR represents the largest dataset published primarily through VO interfaces to date, and so stands as an examplar of the future of survey data mining. Examples of data access are given, including a cross-matching of IPHAS photometry with sources in the UKIDSS Galactic Plane Survey that validates the existing calibration of the best data.

As part of the first high resolution far-IR spectral survey of the Orion KL region (Lerate et al. 2006), we observed 20 CO emission lines with Jup=16 to Jup=39 (upper levels from approx 752 K to 4294 K above the ground state). Observations were taken using the Long Wavelength Spectrometer (LWS) on board the Infrared Space Observatory (ISO), in its high resolution Fabry-Perot (FP) mode (approx 33 km s$^{-1}$). We present here an analysis of the final calibrated CO data, performed with a more sophisticated modelling technique than hitherto, including a detailed analysis of the chemistry, and discuss similarities and differences with previous results. The inclusion of chemical modelling implies that atomic and molecular abundances are time-predicted by the chemistry. This provides one of the main differences with previous studies in which chemical abundances needed to be assumed as initial condition. The chemistry of the region is studied by simulating the conditions of the different known components of the KL region: chemical models for a hot core, a plateau and a ridge are coupled with an accelerated Lambda-iteration (ALI)radiative transfer model to predict line fluxes and profiles. We conclude that the CO transitions with 18<Jup<25 mainly arise from a hot core of diameter 0.02 pc and a density of 10$^{7}$ cm$^{-3}$ rather from the plateau as previous studies had indicated.

(Abridged) Results from the first dedicated study of Galactic PNe by means of optical integral field spectroscopy with the VLT FLAMES Argus IFU are presented. Three typical Galactic-disk PNe have been mapped with the 11.5''x7.2'' Argus array: two dimensional spectral maps of NGC 5882, 6153 and 7009 with 297 spatial pixels per target were obtained at sub-arcsec resolutions and 297 spectra per target were obtained in the 396.4-507.8 nm range. Spatially resolved maps of emission lines and of nebular physical properties were produced. The abundances of helium and of doubly ionized carbon and oxygen were derived from optical recombination lines (ORLs), while those of O^2+ were also derived from the collisionally excited lines (CELs). The abundance discrepancy problem was investigated by mapping the ratio of ORL/CEL abundances for O^2+ (the abundance discrepancy factor; ADF) across the face of the PNe. The ADF varies between targets and also with position within the targets attaining values of ~40 in the case of NGC 6153. Correlations of the ADF with geometric distance from the nucleus, as well as with [O III] electron temperature, plasma ionization state and other physical properties are established. Very small values of the temperature fluctuation parameter in the plane of the sky are found in all cases. It is argued that these results provide further evidence for the existence in typical PNe of a distinct nebular component consisting of hydrogen-deficient plasma. The zones containing this posited component appear as undulations in the C II and C II ORL abundance diagnostics of about 2 spatial pixels across; any associated structures should have physical sizes of less than ~1000 astronomical units. We propose that circumstellar disks, Abell 30-type knots, or Helix-type cometary globules may be involved.

We have analysed deep optical spectra of the `born-again' planetary nebula Abell 58 and its hydrogen-deficient knot, surrounding V605 Aql, which underwent a nova-like eruption in 1919. The electron temperature we derive for the central knot varies widely depending on the diagnostic used. The [O III] nebular-to-auroral transition ratio gives a temperature of 20800 K, while the ratio of the [N II] nebular and auroral lines give Te=15200 K. The helium line ratios 5876/4471 and 6678/4471 imply temperatures of 350 K and 550 K respectively. Weakly temperature-sensitive O II recombination line ratios imply similarly low electron temperatures. Abundances derived from recombination lines are vastly higher than those found from collisionally excited lines, with the abundance discrepancy factor (adf) for O2+ reaching 89 -- the second highest known value after that found for the hydrogen deficient knots in Abell 30. The observed temperature diagnostics and abundances support the idea that, like Abell 30, the knot of Abell 58 contains some very cold ionised material. Although the central star is carbon-rich (C/O>1), the knot is found to be oxygen-rich, a situation not predicted by the single-star `born again' theory of its formation. We argue that the abundances in the ejecta observed in A 30 and A 58 have more in common with neon novae than with Sakurai's Object, which is believed to have undergone a final helium flash. In particular, the C/O ratio of less than unity and presence of substantial quantities of neon in the ejecta of both Abell 30 and Abell 58 are not predicted by very late thermal pulse models.

We present a three-dimensional photoionisation and dust radiative transfer model of NGC 6302, an extreme, high-excitation planetary nebula. We use the 3D photoionisation code Mocassin} to model the emission from the gas and dust. We have produced a good fit to the optical emission-line spectrum, from which we derived a density distribution for the nebula. A fit to the infrared coronal lines places strong constraints on the properties of the unseen ionising source. We find the best fit comes from using a 220,000 K hydrogen-deficient central star model atmosphere, indicating that the central star of this PN may have undergone a late thermal pulse. We have also fitted the overall shape of the ISO spectrum of NGC 6302 using a dust model with a shallow power-law size distribution and grains up to 1.0 micron in size. To obtain a good fit to the infrared SED the dust must be sufficiently recessed within the circumstellar disk to prevent large amounts of hot dust at short wavelengths, a region where the ISO spectrum is particularly lacking. These and other discoveries are helping to unveil many properties of this extreme object and trace it's evolutionary history.

We present new, very late-time optical photometry and spectroscopy of the interesting Type II-P supernova, SN 2002hh, in NGC 6946. Gemini/GMOS-N has been used to acquire visible spectra at six epochs between 2004 August and 2006 July, following the evolution of the SN from age 661 to 1358 days. Few optical spectra of Type II supernovae with ages greater than one year exist. In addition, g'r'i' images were acquired at all six epochs. The spectral and photometric evolution of SN 2002hh has been very unusual. Measures of the brightness of this SN, both in the R and I bands as well as in the H-alpha emission flux, show no significant fading over an interval of nearly two years. The most straightforward explanation for this behavior is that the light being measured comes not only from the SN itself but also from an echo off of nearby dust. Echoes have been detected previously around several SNe but these echoes, at their brightest, were ~8 mag below the maximum brightness of the SN. At V~21 mag, the putative echo dominates the light of SN 2002hh and is only ~4 mag below the outburst's peak brightness. There is an estimated 6 magnitudes of total extinction in V towards SN 2002hh. The proposed explanation of a differential echo/SN absorption is inconsistent with the observed BVRI colors.

We present the scientific motivation and observing plan for an upcoming detection survey for debris disks using the James Clerk Maxwell Telescope. The SCUBA-2 Unbiased Nearby Stars (SUNS) Survey will observe 500 nearby main sequence and sub-giant stars (100 of each of the A, F, G, K and M spectral classes) to the 850 micron extragalactic confusion limit to search for evidence of submillimeter excess, an indication of circumstellar material. The survey distance boundaries are 8.6, 16.5, 22, 25 and 45 pc for M, K, G, F and A stars, respectively, and all targets lie between the declinations of -40 deg to 80 deg. In this survey, no star will be rejected based on its inherent properties: binarity, presence of planetary companions, spectral type or age. This will be the first unbiased survey for debris disks since IRAS. We expect to detect ~125 debris disks, including ~50 cold disks not detectable in current shorter wavelength surveys. A substantial amount of complementary data will be required to constrain the temperatures and masses of discovered disks. High resolution studies will likely be required to resolve many of the disks. Therefore, these systems will be the focus of future observational studies using a variety of observatories to characterize their physical properties. For non-detected systems, this survey will set constraints (upper limits) on the amount of circumstellar dust, of typically 200 times the Kuiper Belt mass, but as low as 10 times the Kuiper Belt mass for the nearest stars in the sample (approximately 2 pc).

We present a study of the effects of clumping on the emergent spectral energy distribution (SED) from dusty supernova (SN) shells illuminated by a diffuse radiation source distributed throughout the medium. (...) The fully 3D radiation transport problem is solved using a Monte Carlo code, MOCASSIN, and we present a set of models aimed at investigating the sensitivity of the SEDs to various clumping parameters. We find that, contrary to the predictions of analytical prescriptions, the combination of an optical and IR observational data set is sufficient to constrain dust masses even in the case where optically thick clumps are present. Using both smoothly varying and clumped grain density distributions, we obtain new estimates for the mass of dust condensed by the Type II SN 1987A by fitting the optical and infrared spectrophotometric data of Wooden et al. (1993) at two epochs (day 615 and day 775). (...) From our numerical models we derive dust masses for SN 1987A that are comparable to previous analytic clumped graphite grain mass estimates, and at least two orders of magnitude below the 0.1-0.3 Msol that have been predicted to condense as dust grains in primordial core collapse supernova ejecta. This low condensation efficiency for SN 1987A is in contrast to the case of SN 2003gd, for which a dust condensation efficiency as large as 0.12 has recently been estimated. (Abridged)

We present the detection and modeling of more than 70 far-IR pure rotational lines of water vapor, including the 18O and 17O isotopologues, towards Orion KL. Observations were performed with the Long Wavelength Spectrometer Fabry-Perot (LWS/FP; R~6800-9700) on board the Infrared Space Observatory (ISO) between ~43 and ~197 um. The water line profiles evolve from P-Cygni type profiles (even for the H2O18 lines) to pure emission at wavelengths above ~100 um. We find that most of the water emission/absorption arises from an extended flow of gas expanding at 25+-5 kms^-1. Non-local radiative transfer models show that much of the water excitation and line profile formation is driven by the dust continuum emission. The derived beam averaged water abundance is 2-3x10^-5. The inferred gas temperature Tk=80-100 K suggests that: (i) water could have been formed in the "plateau" by gas phase neutral-neutral reactions with activation barriers if the gas was previously heated (e.g. by shocks) to >500 K and/or (ii) H2O formation in the outflow is dominated by in-situ evaporation of grain water-ice mantles and/or (iii) H2O was formed in the innermost and warmer regions (e.g. the hot core) and was swept up in ~1000 yr, the dynamical timescale of the outflow.

The Isaac Newton Telescope Photometric H-alpha Survey (IPHAS) is currently mapping the Northern Galactic plane reaching to r'=20 mag with typically 1" resolution. Hundreds of Planetary Nebulae (PNe), both point-like and resolved, are expected to be discovered. We report on the discovery of the first new PN from this survey: it is an unusual object located at a large galactocentric distance and has a very low oxygen abundance. The nebula shows an intricate morphology: there is an inner ring surrounding the central star, bright inner lobes with an enhanced waist, and very faint lobular extensions reaching up to more than 100". We classify it as a quadrupolar PN, a rather unusual class of planetary showing two pairs of misaligned lobes. From long-slit spectroscopy we derive Te[NII] =12800+-1000K, Ne = 390+-40 cm-3, and chemical abundances typical of Peimbert's Type I nebulae (He/H =0.13, N/O =1.8) with an oxygen abundance of 12 + log(O/H)=8.17+-0.15. A kinematic distance of 7.0 kpc is derived, implying an unusually large size of >4 pc for the nebula. The photometry of the central star indicates the presence of a relatively cool companion. This, and the evidence for a dense circumstellar disk and quadrupolar morphology, all of which are rare among PNe, support the hypothesis that this morphology is related to binary interaction.

We present late-time optical and mid-infrared observations of the Type-II supernova 2003gd in NGC 628. Mid-infrared excesses consistent with cooling dust in the ejecta are observed 499-678 days after outburst, and are accompanied by increasing optical extinction and growing asymmetries in the emission-line profiles. Radiative-transfer models show that up to 0.02 solar masses of dust has formed within the ejecta, beginning as early as 250 days after outburst. These observations show that dust formation in supernova ejecta can be efficient and that massive-star supernovae can be major dust producers throughout the history of the Universe.

We have carried out a high spectral resolution line survey towards the Orion Kleinmann-Low (KL) cluster from 44-188 um. The observations were taken with the Long Wavelength Spectrometer (LWS) in Fabry-Perot mode, on board the Infrared Space Observatory (ISO). A total of 152 lines are clearly detected and a further 34 features are present as possible detections. The spectrum is dominated by the molecular species H2O, OH and CO, along with [OI] and [CII] lines from PDR or shocked gas and [OIII], [NIII] lines from the foreground M42 HII region. Several isotopic species, as well as NH3, are also detected. HDO and H3O+ are tentatively detected for the first time in the far-infrared range towards Orion-KL. A basic analysis of the line observations is carried out, by comparing with previous measurements and published models and deriving rotational temperatures and column densities in the case of the molecular species. The complexity of the region requires more sophisticated models for the interpretation of all the line observations.

The assumption of spherical symmetry is not justified for the vast majority of PNe. The interpretation of spatially-resolved observations cannot rely solely on the application of 1D codes, which may yield incorrect abundances determinations resulting in misleading conclusions. The 3D photoionisation code MOCASSIN (Monte CAarlo SimulationS of ionised Nebulae) is designed to remedy these shortcomings. The 3D transfer of both primary and secondary radiation is treated self-consistently without the need of approximations. The code was benchmarked and has been applied to the study of several PNe. The current version includes a fully self-consistent radiative transfer treatment for dust grains mixed within the gas, taking into account the microphysics of dust-gas interactions within the geometry-independent Monte Carlo transfer. The new code provides an excellent tool for the self-consistent analysis of dusty ionised regions showing asymmetries and/or density and chemical inhomogeneities. Work is currently in progress to incorporate the processes that dominate the thermal balance of photo-dissociation regions (PDRs), as well as the formation and destruction processes for all the main molecular species.

We have acquired high spectral resolution observations (R=150,000) of the planetary nebulae NGC 7009 and NGC 6153, using bHROS on Gemini South. Observations of this type may provide a key to understanding why optical recombination lines (ORLs) yield systematically higher heavy element abundances for photoionized nebulae than do the classical forbidden collisionally excited lines (CELs) emitted by the same ions; NGC 7009 and NGC 6153 have notably high ORL/CEL abundance discrepancy factors (ADFs) of 5 and 10, respectively. Due to the opposite temperature dependences of ORLs and CELs, ORLs should be preferentially emitted by colder plasma. Our bHROS observations of NGC 7009 reveal that the [O III] 4363A CEL has a FWHM linewidth that is 1.5 times larger than that shown by O II ORLs in the same spectrum, despite the fact that all of these lines are emitted by the O2+ ion. The bHROS spectra of NGC 6153 also show that its O II ORLs have significantly narrower linewidths than do the [O III] 4363A and 5007A lines but, in addition, the [O III] 4363A and 5007A lines show very different velocity profiles, implying the presence of large temperature variations in the nebula.

Recent weak emission-line long-slit surveys and modelling studies of PNe have convincingly argued in favour of the existence of an unknown component in the planetary nebula plasma consisting of cold, hydrogen-deficient gas, as an explanation for the long-standing recombination-line versus forbidden-line temperature and abundance discrepancy problems. Here we describe the rationale and initial results from a detailed spectroscopic study of three Galactic PNe undertaken with the VLT FLAMES integral-field unit spectrograph, which advances our knowledge about the small-scale physical properties, chemical abundances and velocity structure of these objects across a two-dimensional field of view, and opens up for exploration an uncharted territory in the study and modelling of PNe and photoionized nebulae in general.

We present high quality optical spectroscopic observations of the planetary nebula (PN) Hf 2-2. The spectrum exhibits many prominent optical recombination lines (ORLs) from heavy element ions. Analysis of the H {\sc i} and He {\sc i} recombination spectrum yields an electron temperature of $\sim 900$ K, a factor of ten lower than given by the collisionally excited [O {\sc iii}] forbidden lines. The ionic abundances of heavy elements relative to hydrogen derived from ORLs are about a factor of 70 higher than those deduced from collisionally excited lines (CELs) from the same ions, the largest abundance discrepancy factor (adf) ever measured for a PN. By comparing the observed O {\sc ii} $\lambda$4089/$\lambda$4649 ORL ratio to theoretical value as a function of electron temperature, we show that the O {\sc ii} ORLs arise from ionized regions with an electron temperature of only $\sim 630$ K. The current observations thus provide the strongest evidence that the nebula contains another previously unknown component of cold, high metallicity gas, which is too cool to excite any significant optical or UV CELs and is thus invisible via such lines. The existence of such a plasma component in PNe provides a natural solution to the long-standing dichotomy between nebular plasma diagnostics and abundance determinations using CELs on the one hand and ORLs on the other.

As part of the first far-IR line survey towards Orion KL, we present the detection of seven new rotationally excited OH Lambda-doublets (at 48, 65, 71, 79, 98 and 115 um). Observations were performed with the Long Wavelength Spectrometer (LWS) Fabry-Perots on board the Infrared Space Observatory (ISO). In total, more than 20 resolved OH rotational lines, with upper energy levels up to 620 K, have been detected at an angular and velocity resolutions of 80$'' and 33 km s^-1 respectively. OH line profiles show a complex behavior evolving from pure absorption, P-Cygni type to pure emission. We also present a large scale 6' declination raster in the OH ^2\Pi_3/2 J=5/2^+-3/2^- and ^2\Pi_3/2 J=7/2^-5/2^+ lines (at 119.441 and 84.597 um) revealing the decrease of excitation outside the core of the cloud. From the observed profiles, mean intrinsic line widths and velocity offsets between emission and absorption line peaks we conclude that most of the excited OH arises from Orion outflow(s), i.e. the ``plateau'' component. We determine an averaged OH abundance relative to H_2 of X(OH)=(0.5-1.0)x10^-6, a kinetic temperature of 100 K and a density of n(H_2)=5x10^5 cm^-3. Even with these conditions, the OH excitation is heavily coupled with the strong dust continuum emission from the inner hot core regions and from the expanding flow itself.

We have carried out JHK polarimetric observations of eleven dusty young stars, by using the polarimeter module IRPOL2 with the near-IR camera UIST on the 3.8-m United Kingdom Infrared Telescope (UKIRT). Our sample targeted systems for which UKIRT-resolvable discs had been predicted by model fits to their spectral energy distributions. Our observations have confirmed the presence of extended polarized emission around TW Hya and around HD 169142. HD 150193 and HD 142666 show the largest polarization values among our sample, but no extended structure was resolved. By combining our observations with HST coronographic data from the literature, we derive the J- and H-band intrinsic polarization radial dependences of TW Hya's disc. We find the disc's polarizing efficiency is higher at H than at J, and we confirm that the J- and H-band percentage polarizations are reasonably constant with radius in the region between 0.9 and 1.3 arcseconds from the star. We find that the objects for which we have detected extended polarizations are those for which previous modelling has suggested the presence of flared discs, which are predicted to be brighter than flat discs and thus would be easier to detect polarimetrically.

We report the discovery of 5 massive Wolf-Rayet (WR) stars resulting from a programme of follow-up spectroscopy of candidate emission line stars in the AAO/UKST Southern Galactic Plane H-alpha survey. The 6195-6775 angstrom spectra of the stars are presented and discussed. A WC9 class is assigned to all 5 stars through comparison of their spectra with those of known late-type WC stars, bringing the known total number of Galactic WC9 stars to 44. Whilst three of the five WC9 stars exhibit near infrared (NIR) excesses characteristic of hot dust emission -- as seen in the great majority of known WC9 stars -- we find that two of the stars show no discernible evidence of such excesses. This increases the number of known WC9 stars without NIR excesses to 7. Reddenings and distances for all 5 stars are estimated.

We present deep optical spectra of the archetypal young planetary nebula (PN) NGC 7027, covering a wavelength range from 3310 to 9160 A. The observations were carried out by uniformly scanning a long slit across the entire nebular surface, thus yielding average optical spectra for the whole nebula. A total of 937 emission features are detected. The extensive line list presented here should prove valuable for future spectroscopic analyses of emission line nebulae. The optical data, together with the archival IUE and ISO spectra, are used to probe the temperature and density structures and to determine the elemental abundances from lines produced by different excitation mechanisms. The C++/H+, N++/H+, O++/H+ and Ne++/H+ ionic abundance ratios derived from optical recombination lines (ORLs) are found to be only slightly higher than those derived from collisionally excited lines (CELs). We conclude that whatever mechanism is causing the BJ/CEL temperature discrepanies and the ORL/CEL abundance discrepancies that have been observed in many PNe, it has an insignificant effect on this bright young compact PN. The properties of the central star are also discussed. Based on the integrated spectrum and using the energy-balance method, we have derived an effective temperature of 219 000 K for the ionizing star. Finally, we report the first detection in the spectrum of this bright young PN of Raman-scattered O VI features at 6830 and 7088 A, pointing to the existence of abundant neutral hydrogen around the ionized regions. (abridged)

Airborne and space-based low-resolution spectroscopy in the 1980s discovered tantalizing quantitative relationships between the gas phase C/O abundance ratio in planetary nebulae (PNe) and the fractions of total far-infrared luminosity radiated by the 7.7 and 11.3 micron bands (the C=C stretch and C-H bend, respectively), of polycyclic aromatic hydrocarbons. Only a very small sample of nebulae was studied in this context, limited by airborne observations of the 7.7 micron band, or the existence of adequate IRAS Low Resolution Spectrometer data for the 11.3 micron band. To investigate these trends further, we have expanded the sample of planetaries available for this study using Infrared Space Observatory (ISO) low-resolution spectra secured with the Short Wavelength Spectrometer (SWS) and the Long Wavelength Spectrometer (LWS). The new sample of 43 PNe, of which 17 are detected in PAH emission, addresses the range from C/O=0.2-13 with the objective of trying to delineate the pathways by which carbon dust grains might have formed in planetaries. For the 7.7-micron and 11.3- micron bands, we confirm that the ratio of band strength to total infrared luminosity is correlated with the nebular C/O ratio. Expressed in equivalent width terms, the cut-on C/O ratio for the 7.7-micron band is found to be 0.6(+0.2,-0.4), in good accord with that found from sensitive ground-based measurements of the 3.3-micron band too.

We present the first 3D Monte Carlo (MC) photoionisation code to include a fully self-consistent treatment of dust radiative transfer (RT) within a photoionised region. This is the latest development (Version 2.0) of the gas-only photoionisation code MOCASSIN (Ercolano et al., 2003a), and employs a stochastic approach to the transport of radiation, allowing both the primary and secondary components of the radiation field to be treated self-consistently, whilst accounting for the scattering of radiation by dust grains mixed with the gas, as well as the absorption and emission of radiation by both the gas and the dust components. A set of rigorous benchmark tests have been carried out for dust-only spherically symmetric geometries and 2D disk configurations. MOCASSIN's results are found to be in agreement with those obtained by well established dust-only RT codes that employ various approaches to the solution of the RT problem. A model of the dust and of the photoionised gas components of the planetary nebula (PN) NGC 3918 is also presented as a means of testing the correct functioning of the RT procedures in a case where both gas and dust opacities are present. The two components are coupled via the heating of dust grains by the absorption of both UV continuum photons and resonance line photons emitted by the gas. The MOCASSIN results show agreement with those of a 1D dust and gas model of this nebula published previously, showing the reliability of the new code, which can be applied to a variety of astrophysical environments.

Dust emission from the Type II supernova SN 2002hh in NGC 6946 has been detected at mid-infrared wavelengths by the Spitzer Space Telescope from 590 to 758 days after outburst and confirmed by higher angular resolution Gemini-N mid-IR observations. The day-600 5.8-24 micron emission can be fit by a 290-K blackbody having a luminosity of 1.6E+7 L_Sol. The minimum emitting radius of 1.1E+17 cm is too large for the emitting dust to have been formed in the supernova ejecta. Using radiative transfer models and realistic dust grain parameters, fits to the observed flux distribution could be obtained with an optically thick dust shell having a mass of 0.10-0.15 M_Sol, corresponding to a total dust+gas mass in excess of 10 M_Sol, suggesting a massive M supergiant or luminous blue variable precursor to this self-obscured object.

The core of planetary nebula NGC 6302 is filled with high-excitation photoionised gas at low expansion velocities. It represents a unique astrophysical situation in which to search for hyperfine structure (HFS) in coronal emission lines from highly ionised species. HFS is otherwise blended by thermal or velocity broadening. Spectra containing [AlVI]~3.66um 3P2-3P1, obtained with Phoenix on Gemini-South at resolving powers of up to 75000, resolve the line into five hyperfine components separated by 20 to 60km/s due to the coupling of the I=5/2 nuclear spin of 27Al with the total electronic angular momentum J. 26Al has a different nuclear spin of I=5, and a different HFS, which allows us to place a 3sigma upper limit on the 26Al/27Al abundance ratio of 1/33. We measure the HFS magnetic-dipole coupling constants for [AlVI], and provide the first estimates of the electric-quadrupole HFS coupling constants obtained through astronomical observations of an atomic transition.

We present a comparison of electron densities derived from optical forbidden line diagnostic ratios for a sample of over a hundred nebulae. We consider four density indicators, the [O II] $\lambda3729/\lambda3726$, [S II] $\lambda6716/\lambda6731$, [Cl III] $\lambda5517/\lambda5537$ and [Ar IV] $\lambda4711/\lambda4740$ doublet ratios. Except for a few H II regions for which data from the literature were used, diagnostic line ratios were derived from our own high quality spectra. For the [O II] doublet ratio, we find that our default atomic data set, consisting of transition probabilities (Aij) from Zeippen (1982} and collision strengths from Pradhan (1976), fit the observations well, although at high electron densities, the [O II]doublet ratio yields densities systematically lower than those given by the [S II] doublet ratio, suggesting that the ratio of Aij of the [O II] doublet,$A(\lambda3729)/A(\lambda3726)$, given by Zeippen (1982) may need to be revised upwards by ~6%. Our analysis also shows that the more recent calculations of [O II] A value by Zeippen (1987a) and collision strengths by McLaughlin & Bell (1998) are inconsistent with the observations at the high and low density limits, respectively, and can therefore be ruled out. We confirm the earlier result of Copetti & Writzl (2002) that the [O II] A values calculated by Wiese et al. (1996) yield electron densities systematically lower than those deduced from the [S II] doublet ratio and that the discrepancy is most likely caused by errors in the A values calculated by Wiese et al. Using our default atomic data set for [\ion{O}{ii}], we find that $N_{\rm e}([\ion{O}{ii}]) \la N_{\rm e}([\ion{S}{ii}]) \approx N_{\rm e}([\ion{Cl}{iii}])< N_{\rm e}([\ion{Ar} {iv}])$.

Deep optical spectra of the high excitation planetary nebula NGC 1501 and its W04 central star are presented. A recombination line abundance analysis of the central star's emission-line spectrum yields He:C:O mass fractions of 0.36:0.48:0.16, similar to those of PG1159 stars. A detailed empirical analysis of the nebular collisionally excited line (CEL) and optical recombination line (ORL) spectrum is presented, together with fully three-dimensional photoionisation modelling of the nebula. We found very large ORL-CEL abundance discrepancy factors (ADFs) for O^2+ (32) and Ne^2+ (33). The mean value of 5100K for the T_e derived from HeI recombination lines ratios is 6000K lower than the value of 11100K implied by the [OIII] line ratio. A three-dimensional photoionisation model of NGC 1501 was constructed using the photoionisation code MOCASSIN, based on our new spectroscopic data and using the three-dimensional electron density distribution determined from long-slit echellograms of the nebula by Ragazzoni et al. (2001). The initial models showed higher degrees of ionisation of heavy elements than indicated by observations. We investigated the importance of the missing low-temperature dielectronic recombination rates for third-row elements and have estimated upper limits to their rate coefficients. Our model heavily under-predicts the optical recombination line emission. We conclude that the presence of a hydrogen-deficient, metal-rich component is necessary to explain the observed ORL spectrum of this object. The existence of such knots could also provide a softening of the radiation field, via the removal of ionising photons by absorption in the knots, thereby helping to alleviate the over-ionisation of the heavy elements in our models.

We present the first 3D Monte Carlo photoionisation code to include a fully self-consistent treatment of dust radiative transfer (RT) within the photoionised region. This is the latest development of the recently published pure photoionisation code MOCASSIN (Ercolano et al., 2003a) and it is currently undergoing several benchmarking tests. The preliminary results of these tests are presented in these conference proceedings. The new code provides the ideal tool for the analysis of dusty ionised regions showing asymmetries and/or density and chemical inhomogeneities

[Abridged] Deep optical observations of the spectra of 12 Galactic planetary nebulae (PNe) and 3 Magellanic Cloud PNe were presented in Paper I by Tsamis et al. (2003b), who carried out an abundance analysis using the collisionally excited forbidden lines. Here, the relative intensities of faint optical recombination lines (ORLs) from ions of carbon, nitrogen and oxygen are analysed in order to derive the abundances of these ions relative to hydrogen. We define an abundance discrepancy factor (ADF) as the ratio of the abundance derived for a heavy element ion from its recombination lines to that derived for the same ion from its ultraviolet, optical or infrared collisionally excited lines (CELs). All of the PNe in our sample are found to have ADF's that exceed unity. There is no dependence of the magnitude of the ADF upon the excitation energy of the UV, optical or IR CEL transition used, indicating that classical nebular temperature fluctuations--i.e. in a chemically homogeneous medium--are not the cause of the observed abundance discrepancies. Instead, we conclude that the main cause of the discrepancy is enhanced ORL emission from cold ionized gas located in hydrogen-deficient clumps inside the main body of the nebulae. We have developed a new electron temperature diagnostic, based upon the relative intensities of the OII 4f-3d 4089A and 3p-3s 4649A recombination transitions. For six out of eight PNe for which both transitions are detected, we derive O2+ ORL electron temperatures of <300 K, very much less than the O2+ forbidden-line and Balmer jump temperatures derived for the same nebulae. These results provide direct observational evidence for the presence of H-deficient, cold plasma regions within the nebulae, consistent with gas cooled largely by infrared fine structure and recombination transitions.

We report the discovery of only the fourth massive WO star to be found in the Milky Way, and only the seventh identified within the Local Group. This has resulted from the first observations made in a programme of follow-up spectroscopy of candidate emission line stars from the AAO/UK Schmidt Southern Galactic Plane H-alpha Survey. The optical spectrum of this star, to become WR 93b in the Catalogue of Galactic Wolf-Rayet stars, is presented and described. WR 93b is classified as WO3 and is shown to be highly reddened (E(B-V) = 2.1 pm 0.1). A recombination line analysis of the emission lines yields the abundance ratios C/He = 0.95 and O/He = 0.13 (by number). Comparisons at near infrared wavelengths of reddening corrected photometry between WR 93b and both of Sand 2 (WO3, D = 49 kpc) and Sand 5 (WO2, D = 1.75 kpc) yields a consistent distance to WR 93b of 3.4 kpc. Positioned at Galactic co-ordinates l = 353.27, b = -0.85, the star is most likely located in the Scutum-Crux Arm of the inner Milky Way. We note that none of the four Galactic WO stars lies significantly beyond the Solar Circle (with two well inside). Estimation of the wind terminal velocity in WR 93b at 5750 km/s makes this star the current wind speed record holder among all non-degenerate stars.

Accumulating evidence points to a binary nature for the Wolf-Rayet ([WC]) central stars, a group that constitutes about 15% of all central stars of planetary nebula. From ISO observations, a dual dust chemistry (oxygen- and carbon-rich) has been shown to be almost exclusively associated with [WC] central stars, a fact that could be explained by O-rich dust residing in a disk, while the C-rich dust being more widely distributes. HST/STIS space resolved spectroscopy of the [WC10] central star CPD-568032, is interpreted as revealing a dust disk or torus around the central star. This, together with CPD-568032's variable lightcurve is taken as an indirect indication of binarity. Finally, we present here, for the first time, preliminary results from a radial velocity survey of central stars. Out of 18 stars with excellent data at least 8 are radial velocity variables. If these turn out to be binaries, it is likely that the central star binary fraction is as high as about 50%.

We have constructed a photoionization model, using the 3-D Monte Carlo code MOCASSIN for one of the hydrogen-deficient knots (J3) of the born-again planetary nebula Abell 30. The model consists of spherical knots, comprising a cold, dense, hydrogen-deficient core with very high metal abundances. The inner core is surrounded by a less dense hydrogen-deficient and metal-enriched gas envelope, with less extreme abundances. The envelope of the knot might have been formed by the mixing of the knot material with the surrounding nebular gas. This bi-chemistry, bi-density model did not produce enough heating to match the fluxes of the collisionally excited emission lines (CELs) and of the optical recombination lines (ORLs) observed in the spectrum of the knot. We therefore included heating by photoelectric emission from dust grains in the thermal equilibrium calculations, and found that dust-to-gas ratios of 0.077 and 0.107 by mass for the central core and the envelope of the knot, respectively, are sufficient to fit the spectrum. Surprisingly, photoelectric emission from grains is the dominant source of heating in the hot envelope of the knot, while heating by photoionization of helium and heavy elements dominates in the cold core. The chemical abundances inferred from our modelling indicate, in agreement with the empirical analysis of Wesson et al. (2003), that the C/O ratio in the two regions of the knot is less than unity, contrary to theoretical predictions for born-again nebulae.

We present spectrophotometry of 12 Galactic and 3 Magellanic Cloud planetary nebulae (PNe). Nine of the Galactic PNe were observed by scanning the slit across the PN. We use the fluxes of collisionally excited lines (CELs) to derive electron densities (D's) and temperatures (T's), and ionic abundances. We find that the D's derived from optical CEL ratios are systematically higher than those derived from the ratios of the IR fine-structure (FS) lines of [OIII], indicating the presence of significant density variations within the PNe. We also compare T's obtained from the ratio of optical nebular to auroral [OIII] lines with those obtained from the ratio of [OIII] optical to IR FS lines. We find that when the latter are derived using D's based on the [OIII] 52um/88um ratio, they yield values that are significantly higher than the optical [OIII] T's. Contrasting this, [OIII] optical/IR T's derived using the higher D's obtained from [ClIII] 5517A/5537A ratios show much closer agreement with optical [OIII] T's, implying that the observed [OIII] optical/IR ratios are significantly weighted by D's in excess of the critical densities of both [OIII] FS lines. Consistent with this, ionic abundances derived from [OIII] and [NIII] FS lines using D's from optical CELs show much better agreement with abundances derived for the same ions from optical and UV CELs than do abundances derived from the FS lines using the lower D's obtained from the 52um/88um ratios. The behaviour of T's obtained making use of the T-insensitive IR FS lines provides no support for significant T-fluctuations within the PNe that could be responsible for derived Balmer jump T's being lower than those obtained from the much more T-sensitive [OIII] optical lines.

We describe a multi-order spectrograph concept suitable for 8m-class telescopes, using the intrinsic spectral resolution of Superconducting Tunneling Junction detectors to sort the spectral orders. The spectrograph works at low orders, 1-5 or 1-6, and provides spectral coverage with a resolving power of R~8000 from the atmospheric cutoff at 320 nm to the long wavelength end of the infrared H or K band at 1800 nm or 2400 nm. We calculate that the spectrograph would provide substantial throughput and wavelength coverage, together with high time resolution and sufficient dynamic range. The concept uses currently available technology, or technologies with short development horizons, restricting the spatial sampling to two linear arrays; however an upgrade path to provide more spatial sampling is identified. All of the other challenging aspects of the concept - the cryogenics, thermal baffling and magnetic field biasing - are identified as being feasible.

The study of photoionized environments is fundamental to many astrophysical problems. Up to the present most photoionization codes have numerically solved the equations of radiative transfer by making the extreme simplifying assumption of spherical symmetry. Unfortunately very few real astronomical nebulae satisfy this requirement. To remedy these shortcomings, a self-consistent, three-dimensional radiative transfer code has been developed using Monte Carlo techniques. The code, Mocassin, is designed to build realistic models of photoionized nebulae having arbitraries geometry and density distributions with both the stellar and diffuse radiation fields treated self-consistently. In addition, the code is capable of tretating on or more exciting stars located at non-central locations. The gaseous region is approximated by a cuboidal Cartesian grid composed of numerous cells. The physical conditions within each grid cell are determined by solving the thermal equilibrium and ionization balance equations This requires a knowledge of the local primary and secondary radiation fields, which are calculated self-consistently by locally simulating the individual processes of ionization and recombination. The main structure and computational methods used in the Mocassin code are described in this paper. Mocassin has been benchmarked against established one-dimensional spherically symmetric codes for a number of standard cases, as defined by the Lexington/Meudon photoionization workshops (Pequignot et al., 1986; Ferland et al., 1995; Pequignot et al., 2001)\citep{pequignot86,ferland95, pequignot01}. The results obtained for the benchmark cases are satisfactory and are presented in this paper. A performance analysis has also been carried out and is discussed here.

We have analysed optical spectra of two of the hydrogen-deficient knots (J1 & J3) in the born-again planetary nebula Abell 30, together with UV spectra of knots J3 & J4. We determine electron temperatures in the knots based on several diagnostics. The [O III] nebular-to-auroral transition ratio yields temperatures of the order of 17 000 K. The weak temperature dependence of the ratios of helium lines at 4471, 5876 and 6678 Angstroms is used to derive a temperature of 8850 K for knot J3 and 4600 K for knot J1. Ratios of O II recombination lines, which directly measure the temperature in the coldest regions of the knots, are used to derive temperatures of 2500 K for knot J3, and just 500 K for knot J1. We calculate abundances both from collisionally excited lines and from the well-observed recombination spectra of C, N, O and Ne ions. The forbidden line abundances agree well with previous determinations, but the recombination line abundances are several hundred times higher. These results confirm the scenario proposed by Harrington and Feibelman (1984) in which the knots contain a cold core highly enriched in heavy elements. Forbidden lines are almost entirely emitted by the hot outer part of the knot, while recombination lines are emitted predominantly from the cold core. The C/O ratios we derive imply that the knots are oxygen-rich, contrary to theoretical predictions for born-again nebulae (Iben et al 1983).

The three-dimensional Monte Carlo photoionization code Mocassin has been applied to construct a realistic model of the planetary nebula NGC 3918. Three different geometric models were tried. The effects of the interaction of the diffuse fields coming from two adjacent regions of different densities were investigated. These are found to be non-negligible, even for the relatively uncomplicated case of a biconical geometry. We found that the ionization structure of low ionization species near the boundaries is particularly affected. It is found that all three models provided acceptable matches to the integrated nebular optical and ultraviolet spectrum. Large discrepancies were found between all of the model predictions of infrared fine-structure line fluxes and ISO SWS measurements. This was found to be largely due to an offset of ~14 arcsec from the centre of the nebula that affected all of the ISO observations of NGC 3918. For each model, we also produced projected emission-line maps and position-velocity diagrams from synthetic long-slit spectra, which could be compared to recent HST images and ground-based long-slit echelle spectra. This comparison showed that spindle-like model B provided the best match to the observations. We have therefore shown that although the integrated emission line spectrum of NGC 3918 can be reproduced by all three of the three-dimensional models investigated in this work, the capability of creating projected emission-line maps and position-velocity diagrams from synthetic long-slit spectra was crucial in allowing us to constrain the structure of this object.

We have obtained deep optical, long-slit spectrophotometry of the Galactic HII regions M 17, NGC 3576 and of the Magellanic Cloud HII regions 30 Doradus, LMC N11B and SMC N66, recording the optical recombination lines (ORLs) of CII, NII and OII. Temperature-insensitive ORL C2+/O2+ and N2+/O2 ratios are obtained for all nebulae except SMC N66. The ORL C2+/O2+ ratios show remarkable agreement within each galactic system, while also being in agreement with the corresponding CEL ratios. For all five nebulae, the O2+/H+ abundance derived from multiple OII ORLs is found to be higher than the corresponding value derived from the strong [OIII] 4959, 5007A CELs, by factors of 1.8--2.7 for four of the nebulae. The LMC N11B nebula exhibits a more extreme discrepancy factor for the O2+ ion, ~5. Thus these HII regions exhibit ORL/CEL abundance discrepancy factors that are similar to those previously encountered amongst planetary nebulae. Our optical CEL O2+/H+ abundances agree to within 20-30 per cent with published O2+/H+ abundances that were obtained from observations of infrared fine-structure lines. Since the low excitation energies of the latter make them insensitive to variations about typical nebular temperatures, fluctuations in temperature are ruled out as the cause of the observed ORL/CEL O2+ abundance discrepancies. We present evidence that the observed OII ORLs from these HII regions originate from gas of very similar density (<3500 cm-3) to that emitting the observed heavy-element optical and infrared CELs, ruling out models that employ high-density ionized inclusions in order to explain the abundance discrepancy. We consider a scenario whereby much of the heavy-element ORL emission originates from cold (<=500 K) metal-rich ionized regions.

We present Hubble Space Telescope ultraviolet and optical STIS spectroscopy of the [WCL] planetary nebula central star CPD-568032, obtained during its latest lightcurve minimum. The UV spectrum shows the central star's continuum light distribution to be split into two bright peaks separated by 0.10 arcsec. We interpret this finding as due to an edge-on disk or torus structure that obscures direct light from the star, which is seen primarily via its light scattered from the disk's rims or lobes. CPD-568032 is an archetype of dual dust chemistry [WCL] planetary nebulae, which exhibit strong infrared emission features from both carbon-rich and oxygen-rich materials, and for which the presence of a disk harboring the O-rich grains had been suggested. Our direct observation of an edge-on occulting dust structure around CPD-568032 provides strong support for such a model and for binary interactions being responsible for the correlation between the dual dust chemistry phenomenon in planetary nebulae and the presence of a hydrogen-deficient [WCL] Wolf-Rayet central star.

We present new HST/WFPC2 imagery and STIS long-slit spectroscopy of the planetary nebula NGC 7009. The primary goal was to obtain high spatial resolution of the intrinsic line ratio [O III] 4364/5008 and thereby evaluate the electron temperature (T_e) and the fractional mean-square T_e variation (t_A^2) across the nebula. The WFPC2 T_e map is rather uniform; almost all values are between 9000 - 11,000 K, with the higher T_e's closely coinciding with the inner He^++-zone. The results indicate very small values - <~ 0.01 - for t_A^2 throughout. Our STIS data allow an even more direct determination of T_e and t_A^2, albeit for a much smaller area than with WFPC2. We present results from binning the data along the slit into tiles that are 0.5'' square (matching the slit width). The average [O III] temperature using 45 tiles (excluding the central star and STIS fiducial bars) is 10,139 K; t_A^2 is 0.0035. The measurements of T_e reported here are an average along each line of sight. Therefore, despite finding remarkably low t_A^2, we cannot completely rule out temperature fluctuations along the line of sight as the cause of the large abundance discrepancy between heavy element abundances inferred from collisionally excited emission lines compared to those derived from recombination lines.

This is the first paper in a series of three where we present the first comprehensive inventory of solid state emission bands observed in a sample of 17 oxygen-rich circumstellar dust shells surrounding evolved stars. The data were taken with the Short and Long Wavelength Spectrographs on board of the Infrared Space Observatory (ISO) and cover the 2.4 to 195 micron wavelength range. The spectra show the presence of broad 10 and 18 micron bands that can be attributed to amorphous silicates. In addition, at least 49 narrow bands are found whose position and width indicate they can be attributed to crystalline silicates. Almost all of these bands were not known before ISO. We have measured the peak positions, widths and strengths of the individual, continuum subtracted bands. Based on these measurements, we were able to order the spectra in sequence of decreasing crystalline silicate band strength. We found that the strength of the emission bands correlates with the geometry of the circumstellar shell, as derived from direct imaging or inferred from the shape of the spectral energy distribution. This naturally divides the sample into objects that show a disk-like geometry (strong crystalline silicate bands), and objects whose dust shell is characteristic of an outflow (weak crystalline silicate bands). All stars with the 33.6 micron forsterite band stronger than 20 percent over continuum are disk sources. We define spectral regions (called complexes) where a concentration of emission bands is evident, at 10, 18, 23, 28, 33, 40 and 60 micron. We derive average shapes for these complexes and compare these to the individual band shapes of the programme stars.

We report the discovery of a double ring structure in the waist of the nebula surrounding eta Carinae. The rings are detected in the mid-IR dust continuum at wavelengths of 7.9, 11.9, 12.9 and 20 micronmeter. The dust in the rings has a temperature of about 300 K. The orientation of the rings is inclined with respect to the axis of the homunculus by either 37 or 58 degrees. The central star is not in the projected centre of the structure defined by the two rings. This geometry is reminiscent of that seen in SN1987A and some planetary nebulae. We discuss several possible origins for this remarkable geometry and its orientation.

The hydrogen-deficient carbon-rich [WCL] type central star HD167362 and its oxygen-rich planetary nebula (PN) SwSt~1 are investigated. The nebular chemistry might indicate a recent origin for the carbon-rich stellar spectrum. Its stellar and nebular properties might therefore provide further understanding of the origin of the [WCL] central star class. The UV-IR stellar spectra are modelled with state of the codes and show ~40kK central star with a wind and a C/O~3, indicative of efficient third dredge-up. The synthetic stellar flux distribution is used to model the high density, compact PN, which has a solar C/O ratio, is still enshrouded by 1200K and 230K dust shells and, reported here for the first time, in molecular hydrogen. Although it appears that the change in C/O ratio has been recent, the published spectroscopy since 1895 has been re-examined and no clear spectral change is seen. If an event occurred that has turned it into a hydrogen-deficient central star, it did not happen in the last 100 years.

We present the combined Infrared Space Observatory Short-Wavelength Spectrometer and Long-Wavelength Spectrometer 2.4--197 micron spectrum of the Planetary Nebula NGC 6302 which contains in addition to strong atomic lines, a series of emission features due to solid state components. The broad wavelength coverage enables us to more accurately identify and determine the properties of both oxygen- and carbon-rich circumstellar dust. A simple model fit was made to determine the abundance and typical temperature of the amorphous silicates, enstatite and forsterite. Forsterite and enstatite do have roughly the same abundance and temperature. The origin and location of the dust in a toroidal disk around the central star are discussed.

Far-infrared spectra of laboratory silicates are normally obtained at room temperature even though the grains responsible for astronomical silicate emission bands seen at wavelengths >20 micron are likely to be at temperatures below ~150 K. In order to investigate the effect of temperature on silicate spectra, we have obtained absorption spectra of powdered forsterite and olivine, along with two orthoenstatites and diopside clinopyroxene, at 3.5+-0.5 K and at room temperature (295+-2K). To determine the changes in the spectra the resolution must be increased from 1 to 0.25 cm^-1 at both temperatures since a reduction in temperature reduces the phonon density, thereby reducing the width of the infrared peaks. Several bands observed at 295 K split at 3.5 K. At 3.5 K the widths of isolated single bands in olivine, enstatites and diopside are ~ 90% of their 295 K-widths. However, in forsterite the 3.5-K-widths of the 31-, 49- and 69-micron bands are, respectively, 90%, 45% and 31% of their 295 K widths. Due to an increase in phonon energy as the lattice contracts, 3.5-K-singlet peaks occur at shorter wavelengths than do the corresponding 295-K peaks; the magnitude of the wavelength shift increases from \~ 0-0.2 micron at 25 micron to ~0.9 micron at 80 micron. Changes in the relative absorbances of spectral peaks are also observed. The temperature dependence of lambda_pk and bandwidth shows promise as a means to deduce characteristic temperatures of mineralogically distinct grain populations. In addition, the observed changes in band strength with temperature will affect estimates of grain masses and relative mineral abundances inferred using room-temperature laboratory data.

We report on advances in the study of the cores of NGC6302 and NGC6537 using infrared grating and echelle spectroscopy. In NGC6302, emission lines from species spanning a large range of ionization potential, and in particular [SiIX]3.934um, are interpreted using photoionization models (including CLOUDY), which allow us to reestimate the central star's temperature to be about 250000K. All of the detected lines are consistent with this value, except for [AlV] and [AlVI]. Aluminium is found to be depleted to one hundredth of the solar abundance, which provides further evidence for some dust being mixed with the highly ionized gas (with photons harder than 154eV). A similar depletion pattern is observed in NGC6537. Echelle spectroscopy of IR coronal ions in NGC6302 reveals a stratified structure in ionization potential, which confirms photoionization to be the dominant ionization mechanism. The lines are narrow (< 22km/s FWHM), with no evidence of the broad wings found in optical lines from species with similar ionization potentials, such as [NeV]3426A. We note the absence of a hot bubble, or a wind blown bipolar cavity filled with a hot plasma, at least on 1'' and 10km/s scales. We also provide accurate new wavelengths for several of the infrared coronal lines observed with the echelle.

Infrared spectra of a number of evolved O-rich stars have been obtained with the Short- and Long- Wavelength spectrometers on board the Infrared Space Observatory. The very broad wavelength coverage (2.4-197 micron) obtained by combining observations made with the two spectrometers includes practically all of the flux emitted by the sources, and allows us to determine the emission and absorption features of the dense circumstellar dust shells. Agreement between the fluxes obtained by the two instruments is generally very good; the largest discrepancies are probably due to source variability. Our sample of oxygen-rich AGB stars exhibits a wealth of spectral features due to crystalline silicates and crystalline water ice in emission and absorption. In this study a qualitative overview of all features due to crystalline silicates and water ice in these high mass loss rate objects is presented. It seems that there is a certain onset value for the mass loss rate above which these features appear in the spectrum. Moreover, crystalline silicate emission features have been detected for the first time at wavelengths where the amorphous silicates are still in absorption, implying different spatial distributions for the two materials. A spherically symmetric and an axi-symmetric geometry are proposed.

Pre-solar grains of silicon carbide found in meteorites and interpreted as having had an origin around carbon stars from their isotopic composition, have all been found to be of the beta-SiC polytype. Yet to date fits to the 11.3 microns SiC emission band of carbon stars had been obtained only for alpha-SiC grains. We present thin film infrared (IR) absorption spectra measured in a diamond anvil cell for both the alpha- and beta- polymorphs of synthetic SiC and compare the results with previously published spectra taken using the KBr matrix method. We find that our thin film spectra have positions nearly identical to those obtained previously from finely ground samples in KBr. Hence, we show that this discrepancy has arisen from inappropriate `KBr corrections' having been made to laboratory spectra of SiC particles dispersed in KBr matrices. We re-fit a sample of carbon star mid-IR spectra, using laboratory data with no KBr correction applied, and show that beta-SiC grains fit the observations, while alpha-SiC grains do not. The discrepancy between meteoritic and astronomical identifications of the SiC-type is therefore removed. This work shows that the diamond anvil cell thin film method can be used to produce mineral spectra applicable to cosmic environments without further manipulation.

Combined ISO SWS and LWS spectroscopy is presented of the late WC-type planetary nebula nucleus CPD-56 8032 and its carbon-rich nebula. The extremely broad coverage (2.4-197 microns) enables us to recognize the clear and simultaneous presence of emission features from both oxygen- and carbon- rich circumstellar materials. Removing a smooth continuum highlights bright emission bands characteristic of polycyclic aromatic hydrocarbons (hereafter PAHs) in the 3-14 micron region, bands from crystalline silicates longwards of 18 microns, and the 43- and 62-micron bands of crystalline water ice. We discuss the probable evolutionary state and history of this unusual object in terms of (a) a recent transition from an O-rich to a C-rich outflow following a helium shell flash; or (b) a carbon-rich nebular outflow encountering an O-rich comet cloud.

The morphology of planetary nebulae (PNe) provides an essential tool for understanding their origin and evolution, as it reflects both the dynamics of the gas ejected during the TP-AGB phase, and the central star energetics. Here we study the morphology of 27 Magellanic Cloud planetary nebulae (MCPNe) and present an analysis of their physical characteristics across morphological classes. Similar studies have been successfully carried out for galactic PNe, but were compromised by the uncertainty of individual PN distances. We present our own HST/FOC images of 15 Magellanic Cloud PNe (MCPNe) acquired through a narrow-band lambda 5007 [O III] filter. We use the Richardson-Lucy deconvolution technique on these pre-COSTAR images to achieve post-COSTAR quality. Three PNe imaged before and after COSTAR confirm the high reliability of our deconvolution procedure. We derive morphological classes, dimensions, and surface photometry for all these PNe. We have combined this sample with HST/PC1 images of 15 MCPNe, three of which are in common with the FOC set, acquired by Dopita et al. (1996), to obtain the largest MCPN sample ever examined from the morphological viewpoint. By using the whole database, supplemented with published data from the literature, we have analyzed the properties of the MCPNe and compared them to a typical, complete galactic sample. Morphology of the MCPNe is then correlated with PN density, chemistry, and evolution.

We present ultra-high resolution (0.32 km/s) spectra obtained with the 3.9m Anglo-Australian Telescope (AAT) and Ultra-High-Resolution Facility (UHRF), of interstellar NaI D1, D2, Ca II K, K I and CH absorption toward two high galactic latitude stars HD141569 and HD157841. We have compared our data with 21-cm observations obtained from the Leiden/Dwingeloo HI survey. We derive the velocity structure, column densities of the clouds represented by the various components and identify the clouds with ISM structures seen in the region at other wavelengths. We further derive abundances, linear depletions and H2 fractional abundances for these clouds, wherever possible. Toward HD141569, we detect two components in our UHRF spectra : a weak, broad component at - 15 km/s, seen only in CaII K absorption and another component at 0 km/s, seen in NaI D1, D2, Ca II K, KI and CH absorption. In the case of the HD157841 sightline, a total of 6 components are seen on our UHRF spectra in NaI D1, D2 Ca II K, K I and CH absorption. 2 of these 6 components are seen only in a single species.