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.Comment: 11 pages, 5 figures, MNRAS accepte

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.Comment: Accepted in Ap

OH+^+ emission from cometary knots in planetary nebulae

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

Three-Dimensional Ionisation, Dust RT and Chemical Modelling of Planetary Nebulae

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

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\,$\mu$m are required to obtain a fit to the observed SED. This suggests that the dust mass increase at late times was caused by accretion onto and coagulation of the dust grains formed at earlier epochs. These findings provide further confirmation that core collapse supernovae can create large quantities of dust, and indicate that the reason for small dust masses being estimated in many cases is that the vast majority of the dust forms long after most supernovae have been detectable at mid-infrared wavelengths.Comment: 13 pages, 16 figures. Accepted for publication in MNRA

3D Photoionisation Modelling of NGC 6302

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

Modelling the ArH+^+ emission from the Crab Nebula

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