Source author record

M. J. Barlow

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

110works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

110 published item(s)

preprint2026arXiv

WEAVE imaging spectroscopy of NGC 6720: an iron bar in the Ring

We present spatially resolved spectroscopic observations of the planetary nebula NGC 6720, the Ring Nebula, taken during the science verification phase of WEAVE, a new instrument mounted on the William Herschel Telescope on La Palma. We use the instrument's Large Integral Field Unit (LIFU) to obtain spectra of the Ring Nebula, covering its entire optically bright inner regions as well as parts of its much fainter outer molecular halo. We report the discovery of emission from [Fe~{\sc v}] and [Fe~{\sc vi}] confined to a narrow ``bar'' extending across the central regions of the nebula. No lines of other elements share this morphology or, at the spectral resolving power used ($R \sim 2500$), the same radial velocity. The extent to which iron in this bar is depleted is presently unclear; comparison with JWST-detected dust continuum emission suggests that some dust grain destruction may be occurring in the region, but there is currently no observational evidence for the $>$ 50~km\,s$^{-1}$ shock waves or $T > 10^6$~K X-ray emitting gas needed to enable this. Where the bar is located along the line of sight through the nebula, and how it was created, are new puzzles to be solved for this iconic planetary nebula.

preprint2023arXiv

First Measurement of Neutron Birefringence in Polarized $^{129}$Xe and $^{131}$Xe Nuclei

We present the first measurements of polarized neutron birefringence in transmission through nuclear-polarized $^{129}$Xe and $^{131}$Xe gas and determine the neutron incoherent scattering lengths $b_i(^{129}Xe)=0.186\pm(0.021)_{stat.}\pm(0.004)_{syst.}\space\text{ fm}$ and $b_i(^{131}Xe)=2.09\pm(0.29)_{stat.}\pm(0.12)_{syst.}\space\text{ fm}$ for the first time. These results determine the essential parameter needed for interpretation of spin-dependent neutron-scattering studies on polarized xenon ensembles, with possible future applications ranging from tests of time-reversal violation to mode-entangled neutron scattering experiments on nuclear-polarized systems.

preprint2023arXiv

Quantifying the dust in SN 2012aw and iPTF14hls with ORBYTS

Core-collapse supernovae (CCSNe) are potentially capable of producing large quantities of dust, with strong evidence that ejecta dust masses can grow significantly over extended periods of time. Red-blue asymmetries in the broad emission lines of CCSNe can be modelled using the Monte Carlo radiative transfer code DAMOCLES, to determine ejecta dust masses. To facilitate easier use of DAMOCLES, we present a Tkinter graphical user interface (GUI) running DAMOCLES. The GUI was tested by high school students as part of the Original Research By Young Twinkle Students (ORBYTS) programme, who used it to measure the dust masses formed at two epochs in two Type IIP CCSNe: SN 2012aw and iPTF14hls, demonstrating that a wide range of people can contribute significantly to scientific advancement. Bayesian methods were used to quantify uncertainties on our model parameters. From the presence of a red scattering wing in the day 1863 H$α$ profile of SN 2012aw, we were able to constrain the dust composition to large (radius $>0.1 μ$m) silicate grains, with a dust mass of $6.0^{+21.9}_{-3.6}\times10^{-4} M_\odot$. From the day 1158 H$α$ profile of SN 2012aw, we found a dust mass of $3.0^{+14}_{-2.5}\times10^{-4}$ M$_\odot$. For iPTF14hls, we found a day 1170 dust mass of 8.1 $^{+81}_{-7.6}\times10^{-5}$ M$_{\odot}$ for a dust composition consisting of 50% amorphous carbon and 50% astronomical silicate. At 1000 days post explosion, SN 2012aw and iPTF14hls have formed less dust than SN 1987A, suggesting that SN 1987A could form larger dust masses than other Type IIP's.

preprint2022arXiv

Spitzer and Herschel studies of dust in supernova remnants in the Small Magellanic Cloud

With the entire Small Magellanic Cloud (SMC) mapped by the Spitzer Space Telescope and Herschel Space Observatory, we were able to search 8-250 micron images in order to identify infrared (IR) emission associated with SMC supernova remnants (SNRs). A valid detection had to correspond with known X-ray, Halpha and radio emission from the SNRs. From the 24 known SNRs, we made 5 positive detections with another 5 possible detections. Two detections are associated with pulsars or pulsar wind nebula, and another three detections are part of the extended nebulous emission from the SNRs. We modelled dust emission where fast moving electrons are predicted to collide and heat dust grains which then radiate in IR. With known distance (62.44+-0.47kpc), measured SNR sizes, electron densities, temperatures from X-ray emission as well as hydrogen densities, the modelling of SMC SNRs is straightforward. If the higher range of hydrogen and electron densities were to be accepted, we would expect almost all SMC SNRs to be detected in the IR, at least at 24 micron, but the actual detection rate is only 25%. One possible and common explanation for this discrepancy is that small grains have been destroyed by the SNRs shockwave. However, within the uncertainties of hydrogen and electron densities, we find that infrared dust emission can be explained reasonably well, without invoking dust destruction. There is no conclusive evidence that SNRs destroy swept-up ISM dust.

preprint2021arXiv

The efficiency of grain growth in the diffuse interstellar medium

Grain growth by accretion of gas-phase metals is a common assumption in models of dust evolution, but in dense gas, where the timescale is short enough for accretion to be effective, material is accreted in the form of ice mantles rather than adding to the refractory grain mass. It has been suggested that negatively-charged small grains in the diffuse interstellar medium (ISM) can accrete efficiently due to the Coulomb attraction of positively-charged ions, avoiding this issue. We show that this inevitably results in the growth of the small-grain radii until they become positively charged, at which point further growth is effectively halted. The resulting gas-phase depletions under diffuse ISM conditions are significantly overestimated when a constant grain size distribution is assumed. While observed depletions can be reproduced by changing the initial size distribution or assuming highly efficient grain shattering, both options result in unrealistic levels of far-ultraviolet extinction. We suggest that the observed elemental depletions in the diffuse ISM are better explained by higher initial depletions, combined with inefficient dust destruction by supernovae at moderate ($n_{\rm H} \sim 30 {\rm \, cm^{-3}}$) densities, rather than by higher accretion efficiences.

preprint2020arXiv

Constraining early-time dust formation in core-collapse supernovae

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.

preprint2020arXiv

Disentangling Dust Components in SN 2010jl: The First 1400 Days

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.

preprint2020arXiv

Silicate grain growth due to ion trapping in oxygen-rich supernova remnants like Cassiopeia A

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 $χ=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 ($χ\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 $χ=256$.

preprint2019arXiv

Dust masses and grain size distributions of a sample of Galactic pulsar wind nebulae

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

preprint2016arXiv

Atomic Gas in Debris Discs

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$_α$, 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 $β$ Pictoris system.

preprint2016arXiv

Herschel detects oxygen in the beta Pictoris debris disk

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.

preprint2016arXiv

Herschel observations of the nebula M1-67 around the Wolf-Rayet star WR 124

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.

preprint2016arXiv

Hi-GAL, the Herschel infrared Galactic Plane Survey: photometric maps and compact source catalogues. First data release for Inner Milky Way: +68°> l > -70°

(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° > l > -70° in a |b|<1° 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.

preprint2016arXiv

The extinction and dust-to-gas structure of the planetary nebula NGC 7009 observed with MUSE

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

preprint2016arXiv

The mass-loss rates of red supergiants at low metallicity: Detection of rotational CO emission from two red supergiants in the Large Magellanic Cloud

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.

preprint2015arXiv

Early Dust Formation and a Massive Progenitor for SN 2011ja?

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.

preprint2015arXiv

Herschel observations of extreme OH/IR stars - the isotopic ratios of oxygen as a sign-post for the stellar mass

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.

preprint2015arXiv

Modelling supernova line profile asymmetries to determine ejecta dust masses: SN 1987A from days 714 to 3604

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$α$ 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$μ$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.

preprint2015arXiv

The dust and gas content of the Crab Nebula

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.

preprint2015arXiv

TORUS-3DPDR: A self-consistent code treating three-dimensional photoionization and photodissociation regions

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.

preprint2014arXiv

A 3D extinction map of the Northern Galactic Plane based on IPHAS photometry

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 .

preprint2014arXiv

A Photodissociation Region study of NGC 4038

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χ\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μm and [OI] 146μ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μm/[OI] 146μ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.

preprint2014arXiv

A stubbornly large mass of cold dust in the ejecta of Supernova 1987A

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.

preprint2014arXiv

Discovery of Time Variation of the Intensity of Molecular Lines in IRC+10216 in The Submillimeter and Far Infrared Domains

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.

preprint2014arXiv

First release of the IPHAS Catalogue of New Extended Planetary Nebulae

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

preprint2014arXiv

Herschel imaging of the dust in the Helix Nebula (NGC 7293)

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.

preprint2014arXiv

Herschel spectral-mapping of the Helix Nebula (NGC 7293): Extended CO photodissociation and OH+ emission

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.

preprint2014arXiv

Herschel-PACS Measurements of Nitrogen Enrichment in Nebulae around Wolf-Rayet Stars

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.

preprint2014arXiv

New accurate measurement of 36ArH+ and 38ArH+ ro-vibrational transitions by high resolution IR absorption spectroscopy

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 $μ$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σ$) 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.

preprint2014arXiv

The Second Data Release of the INT Photometric H-Alpha Survey of the Northern Galactic Plane (IPHAS DR2)

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

preprint2014arXiv

The timing and location of dust formation in the remnant of SN 1987A

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\,$μ$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.

preprint2014arXiv

The VST Photometric Halpha Survey of the Southern Galactic Plane and Bulge (VPHAS+)

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.

preprint2014arXiv

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

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

preprint2013arXiv

A Herschel study of NGC 650

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.

preprint2013arXiv

A multiwavelength study of the Magellanic-type galaxy NGC 4449 - I. Modelling the spectral energy distribution, the ionization structure and the star formation history

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

preprint2013arXiv

Alignment in star-debris disc systems seen by Herschel

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.

preprint2013arXiv

An Overview of the Dwarf Galaxy Survey

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.

preprint2013arXiv

Carbon Monoxide in the Cold Debris of Supernova 1987A

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.

preprint2013arXiv

Detection of a Noble Gas Molecular Ion, {36}ArH^+, in the Crab Nebula

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.

preprint2013arXiv

Dust Production and Particle Acceleration in Supernova 1987A Revealed with ALMA

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.

preprint2013arXiv

First results of an H-alpha based search of classical Be stars in the Perseus Arm and beyond

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.

preprint2013arXiv

Herschel SPIRE and PACS observations of the red supergiant VY CMa: analysis of the molecular line spectra

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.

preprint2013arXiv

Mid-J CO Emission in Nearby Seyfert Galaxies

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.

preprint2013arXiv

OH/IR stars and their superwinds as observed by the Herschel Space Observatory

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.

preprint2013arXiv

Revealing the cold dust in low-metallicity environments: I - Photometry analysis of the Dwarf Galaxy Survey with Herschel

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.

preprint2013arXiv

The Lyman break analogue Haro 11: spatially resolved chemodynamics with VLT FLAMES

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.

preprint2012arXiv

3D-PDR: A new three-dimensional astrochemistry code for treating Photodissociation Regions

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.

preprint2012arXiv

A Cool Dust Factory in the Crab Nebula: A Herschel Study of the Filaments

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.

preprint2012arXiv

An independent distance estimate to CW Leo

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

preprint2012arXiv

Comet-like mineralogy of olivine crystals in an extrasolar proto-Kuiper belt

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.

preprint2012arXiv

Discovery of multiple dust shells beyond 1 arcmin in the circumstellar envelope of IRC+10216 using Herschel/PACS

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.

preprint2012arXiv

Dust in Historical Galactic Type Ia Supernova Remnants with Herschel

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.

preprint2012arXiv

Herschel images of Fomalhaut. An extrasolar Kuiper Belt at the height of its dynamical activity

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.

preprint2012arXiv

On the physical structure of IRC+10216. Ground-based and Herschel observations of CO and CCH

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.

preprint2012arXiv

Optical and Infrared Analysis of Type II SN 2006BC

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$α$ line profile to be red-shifted with respect to the host galaxy H$α$ 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$α$ 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.

preprint2012arXiv

The enigmatic nature of the circumstellar envelope and bow shock surrounding Betelgeuse as revealed by Herschel. I. Evidence of clumps, multiple arcs, and a linear bar-like structure

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]

preprint2012arXiv

The expanding dusty bipolar nebula around the nova V1280 Sco

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.

preprint2012arXiv

The merging dwarf galaxy UM 448: chemodynamics of the ionized gas from VLT integral field spectroscopy

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.

preprint2012arXiv

The nature of the interstellar medium of the starburst low-metallicity galaxy Haro11: a multi-phase model of the infrared emission

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

preprint2012arXiv

The Outer Wind of gamma Velorum

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.

preprint2011arXiv

A Multi-wavelength MOCASSIN Model of the Magellanic-type Galaxy NGC 4449

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.

preprint2011arXiv

Distinguishing between HII regions and planetary nebulae with Hi-GAL, WISE, MIPSGAL, and GLIMPSE

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.

preprint2011arXiv

Evidence for Pre-Existing Dust in the Bright Type IIn SN 2010jl

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 η Carinae.

preprint2011arXiv

Herschel Detects a Massive Dust Reservoir in Supernova 1987A

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.

preprint2011arXiv

Photometric and Spectroscopic Evolution of the Type IIP SN 2007it to Day 944

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.

preprint2011arXiv

The Circumstellar Environment of R Coronae Borealis: White Dwarf Merger or Final Helium Shell Flash?

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.

preprint2011arXiv

The effects of dust on the optical and infrared evolution of SN 2004et

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.

preprint2011arXiv

The spectroscopic properties of the nebulae around seven Galactic and LMC Wolf-Rayet stars

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.

preprint2011arXiv

Warm water vapour in the sooty outflow from a luminous carbon star

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.

preprint2010arXiv

870 um observations of evolved stars with LABOCA

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.

preprint2010arXiv

A Herschel PACS and SPIRE study of the dust content of the Cassiopeia A supernova remnant

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.

preprint2010arXiv

A high resolution line survey of IRC+10216 with Herschel. First results: Detection of warm silicon dicarbide SiC2

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.

preprint2010arXiv

A search for Ejecta Nebulae around Wolf-Rayet Stars using the SHS H$α$ survey

Recent large scale Galactic Plane H$α$ 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$α$ 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.

preprint2010arXiv

Chemical and radiative transfer modeling of the ISO-LWS Fabry-Perot spectra of Orion-KL water lines

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.

preprint2010arXiv

Detection of anhydrous hydrochloric acid, HCl, in IRC+10216 with the Herschel SPIRE and PACS spectrometers

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.

preprint2010arXiv

Herschel and SCUBA-2 imaging and spectroscopy of a bright, lensed submillimetre galaxy at z = 2.3

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.

preprint2010arXiv

Herschel images of NGC 6720: H2 formation on dust grains

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.

preprint2010arXiv

Herschel PACS and SPIRE imaging of CW Leo

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.

preprint2010arXiv

Herschel photometric observations of the low metallicity dwarf galaxy NGC 1705

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.

preprint2010arXiv

Herschel photometric observations of the nearby low metallicity irregular galaxy NGC 6822

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.

preprint2010arXiv

Herschel-SPIRE FTS spectroscopy of the carbon-rich objects AFGL 2688, AFGL 618 and NGC 7027

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.

preprint2010arXiv

In-flight calibration of the Herschel-SPIRE instrument

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 μm, and an imaging Fourier transform spectrometer (FTS) covering 194-671 μ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.

preprint2010arXiv

Mapping the interstellar medium in galaxies with Herschel/SPIRE

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.

preprint2010arXiv

On the origin of M81 group extended dust emission

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 &μ&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.

preprint2010arXiv

PACS and SPIRE Spectroscopy of the Red Supergiant VY CMa

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.

preprint2010arXiv

Probing the molecular interstellar medium of M82 with Herschel-SPIRE spectroscopy

We present the observations of the starburst galaxy M82 taken with the Herschel SPIRE Fourier Transform Spectrometer. The spectrum (194-671 μ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.

preprint2010arXiv

Silicon in the dust formation zone of IRC +10216 as observed with PACS and SPIRE on board Herschel

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 μm) and SPIRE (194-672 μ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.

preprint2010arXiv

SN 2007od: A Type IIP SN with Circumstellar Interaction

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.

preprint2010arXiv

SPIRE imaging of M82: cool dust in the wind and tidal streams

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.

preprint2010arXiv

Testing PDR models against ISO fine structure line data for extragalactic sources

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.

preprint2010arXiv

The central region of spiral galaxies as seen by Herschel. M81, M99 and M100

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.

preprint2010arXiv

The dust morphology of the elliptical Galaxy M86 with SPIRE

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.

preprint2010arXiv

The Herschel Space Observatory view of dust in M81

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.

preprint2010arXiv

The Herschel-SPIRE instrument and its in-flight performance

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.

preprint2010arXiv

The orbital period of V458 Vulpeculae, a post double common-envelope nova

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 λ5412 and HeII λ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.

preprint2010arXiv

The Vega Debris Disc: A view from Herschel

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.

preprint2010arXiv

The β Pictoris disk imaged by Herschel PACS and SPIRE

We obtained Herschel PACS and SPIRE images of the thermal emission of the debris disk around the A5V star β Pic. The disk is well resolved in the PACS filters at 70, 100, and 160 μm. The surface brightness profiles between 70 and 160 μm show no significant asymmetries along the disk, and are compatible with 90% of the emission between 70 and 160 μ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 μ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 (β = 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.

preprint2009arXiv

A VLT VIMOS integral field spectroscopic study of perturbed blue compact galaxies: UM 420 and UM 462

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.

preprint2009arXiv

A VLT VIMOS study of the anomalous BCD Mrk 996: mapping the ionised gas kinematics and abundances

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α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]

preprint2009arXiv

Galactic AGB stars from the IPHAS survey

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.

preprint2009arXiv

IPHAS and the symbiotic stars. II. New discoveries and a sample of the most common mimics

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.

preprint2009arXiv

Near-IR spectra of IPHAS extremely red Galactic AGB stars

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.

preprint2009arXiv

The destruction and survival of dust in the shell around SN 2008S

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.

preprint2009arXiv

The global gas and dust budget of the Large Magellanic Cloud: AGB stars and supernovae, and the impact on the ISM evolution

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.

preprint2006arXiv

Massive-Star Supernovae as Major Dust Factories

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.

preprint2006arXiv

The Water Vapor Abundance in Orion KL Outflows

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.

preprint2004arXiv

Discovery of a WO star in the Scutum-Crux arm of the inner Galaxy

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.