Source author record

D. J. Stock

D. J. Stock 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

12works
3topics
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

12 published item(s)

preprint2016arXiv

Interpreting the subtle spectral variations of the 11.2 and 12.7 μm polycyclic aromatic hydrocarbon bands

We report new properties of the 11 and 12.7 μm emission complexes of polycyclic aromatic hydrocarbons (PAHs) by applying a Gaussian-based decomposition technique. Using high-resolution \textit{Spitzer} Space Telescope data, we study in detail the spectral and spatial characteristics of the 11 and 12.7 μm emission bands in maps of reflection nebulae NGC 7023 and NGC 2023 (North and South) and the star-forming region M17. Profile variations are observed in both the 11 and 12.7 μm emission bands. We identify a neutral contribution to the traditional 11.0 μm PAH band and a cationic contribution to the traditional 11.2 μm band, the latter of which affects the PAH class of the 11.2 μm emission in our sample. The peak variations of the 12.7 μm complex are explained by the competition between two underlying blended components. The spatial distributions of these components link them to cations and neutrals. We conclude that the 12.7 μm emission originates in both neutral and cationic PAHs, lending support to the use of the 12.7/11.2 intensity ratio as a charge proxy.

preprint2016arXiv

Polycyclic Aromatic Hydrocarbon emission in Spitzer/IRS maps I: Catalog and simple diagnostics

We present a sample of resolved galactic HII regions and photodissociation regions (PDRs) observed with the Spitzer infrared spectrograph (IRS) in spectral mapping mode between the wavelengths of 5--15 $μ$m. For each object we have spectral maps at a spatial resolution of $\sim$4" in which we have measured all of the mid-infrared emission and absorption features. These include the PAH emission bands, primarily at 6.2, 7.7, 8.6, 11.2 and 12.7 $μ$m, as well as the spectral emission lines of neon and sulfur and the absorption band caused by silicate dust at around 9.8 $μ$m. In this work we describe the data in detail, including the data reduction and measurement strategies, and subsequently present the PAH emission band intensity correlations for each of the objects and the sample as a whole. We find that there are distinct differences between the sources in the sample, with two main groups, the first comprising the HII regions and the second the reflection nebulae (RNe). Three sources, the reflection nebula NGC~7023, the Horsehead nebula PDR (an interface between the HII region IC~434 and the Orion B molecular cloud) and M 17, resist this categorization, with the Horsehead PDR points mimicking the RNe and the NGC~7023 fluxes displaying unique bifurcated appearance in our correlation plots. These discrepancies seem to be due to the very low radiation field experienced by the Horsehead PDR and the very clean separation between the PDR environment and a diffuse environment in the NGC~7023 observations.

preprint2016arXiv

The probability distribution functions of emission line flux measurements and their ratios

Many physical parameters in astrophysics are derived using the ratios of two observed quantities. If the relative uncertainties on measurements are small enough, uncertainties can be propagated analytically using simplifying assumptions, but for large normally distributed uncertainties, the probability distribution of the ratio become skewed, with a modal value offset from that expected in Gaussian uncertainty propagation. Furthermore, the most likely value of a ratio A/B is not equal to the reciprocal of the most likely value of B/A. The effect is most pronounced when the uncertainty on the denominator is larger than that on the numerator. We show that this effect is seen in an analysis of 12,126 spectra from the Sloan Digital Sky Survey. The intrinsically fixed ratio of the [O III] lines at 4959 and 5007 $Å$ is conventionally expressed as the ratio of the stronger line to the weaker line. Thus, the uncertainty on the denominator is larger, and non-Gaussian probability distributions result. By taking this effect into account, we derive an improved estimate of the intrinsic 5007/4959 ratio. We obtain a value of 3.012 $\pm$ 0.008, which is slightly but statistically significantly higher than the theoretical value of 2.98. We further investigate the suggestion that fluxes measured from emission lines at low signal to noise are strongly biased upwards. We were unable to detect this effect in the SDSS line flux measurements, and we could not reproduce the results of Rola and Pelat who first described this bias. We suggest that the magnitude of this effect may depend strongly on the specific fitting algorithm used.

preprint2015arXiv

Herschel PACS and SPIRE spectroscopy of the Photodissociation Regions associated with S 106 and IRAS 23133+6050

Photodissociation regions (PDRs) contain a large fraction of all of the interstellar matter in galaxies. Classical examples include the boundaries between ionized regions and molecular clouds in regions of massive star formation, marking the point where all of the photons energetic enough to ionize hydrogen have been absorbed. In this paper we determine the physical properties of the PDRs associated with the star forming regions IRAS 23133+6050 and S 106 and present them in the context of other Galactic PDRs associated with massive star forming regions. We employ Herschel PACS and SPIRE spectroscopic observations to construct a full 55-650 μm spectrum of each object from which we measure the PDR cooling lines, other fine- structure lines, CO lines and the total far-infrared flux. These measurements are then compared to standard PDR models. Subsequently detailed numerical PDR models are compared to these predictions, yielding additional insights into the dominant thermal processes in the PDRs and their structures. We find that the PDRs of each object are very similar, and can be characterized by a two-phase PDR model with a very dense, highly UV irradiated phase (n $\sim$ 10^6 cm^(-3), G$_0$ $\sim$ 10^5) interspersed within a lower density, weaker radiation field phase (n $\sim$ 10^4 cm^(-3), G$_0$ $\sim$ 10^4). We employed two different numerical models to investigate the data, firstly we used RADEX models to fit the peak of the $^{12}$CO ladder, which in conjunction with the properties derived yielded a temperature of around 300 K. Subsequent numerical modeling with a full PDR model revealed that the dense phase has a filling factor of around 0.6 in both objects. The shape of the $^{12}$CO ladder was consistent with these components with heating dominated by grain photoelectric heating. An extra excitation component for the highest J lines (J > 20) is required for S 106.

preprint2015arXiv

Probing the Ionization States of Polycyclic Aromatic Hydrocarbons via the 15-20 μm Emission Bands

We report new correlations between ratios of band intensities of the 15-20 μm emission bands of polycyclic aromatic hydrocarbons (PAHs) in a sample of fifty-seven sources observed with Spitzer/IRS. This sample includes Large Magellanic Cloud point sources from the SAGE-Spec survey, nearby galaxies from the SINGS survey, two Galactic ISM cirrus sources and the spectral maps of the Galactic reflection nebulae NGC 2023 and NGC 7023. We find that the 16.4, 17.4 and 17.8 μm band intensities are inter-correlated in all environments. In NGC 2023 and NGC 7023 these bands also correlate with the 11.0 and 12.7 μm band intensities. The 15.8 μm band correlates only with the 15-20 μm plateau and the 11.2 μm emission. We examine the spatial morphology of these bands and introduce radial cuts. We find that these bands can be spatially organized into three sets: the 12.7, 16.4 and 17.8 μm bands; the 11.2, 15.8 μm bands and the 15-18 μm plateau; and the 11.0 and 17.4 μm bands. We also find that the spatial distribution of the 12.7, 16.4 and 17.8 μm bands can be reconstructed by averaging the spatial distributions of the cationic 11.0 μm and neutral 11.2 μm bands. We conclude that the 17.4 μm band is dominated by cations, the 15.8 μm band by neutral species, and the 12.7, 16.4 and 17.8 μm bands by a combination of the two. These results highlight the importance of PAH ionization for spatially differentiating sub-populations by their 15-20 μm emission variability.

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

The mid-infrared appearance of the Galactic Mini-Starburst W49A

The massive star forming region W49A represents one of the largest complexes of massive star formation present in the Milky Way and contains at least fifty young massive stars still enshrouded in their natal molecular cloud. We employ Spitzer/IRS spectral mapping observations of the northern part of W49A to investigate the mid-infrared (MIR) spatial appearance of the polycyclic aromatic hydrocarbon (PAH) bands, PAH plateau features, atomic lines and continuum emission. We examine the spatial variations of the MIR emission components in slices through two of the ultra compact-HII (UC-HII) regions. We find that the PAH bands reproduce known trends, with the caveat that the 6.2 $μ$m PAH band seems to decouple from the other ionized PAH bands in some of the UC-HII regions -- an effect previously observed only in one other object: the giant star forming region N66 in the LMC. Furthermore, we compare the nature of the emission surrounding W49A to that of `diffuse' sightlines. It is found that the surrounding emission can be explained by line of sight emission, and does not represent true `diffuse' material. Additionally, we examine the MIR appearance of star formation on various scales from UC-HII regions to starburst galaxies, including a discussion of the fraction of PAH emission in the 8 $μ$m IRAC filter. We find that the MIR appearance of W49A is that of a starburst on large scales yet its individual components are consistent with other galactic HII regions.

preprint2013arXiv

Extinction and PAH intensity variations across the HII region IRAS 12063-6259

The spatial variations in PAH band intensities are normally attributed to the physical conditions of the emitting PAHs, however in recent years it has been suggested that such variations are caused mainly by extinction. To resolve this question, we have obtained near-infrared (NIR), mid-infrared (MIR) and radio observations of the compact H ii region IRAS 12063-6259. We use these data to construct multiple independent extinction maps and also to measure the main PAH features (6.2, 7.7, 8.6 and 11.2 μm) in the MIR. Three extinction maps are derived: the first using the NIR hydrogen lines and case B recombination theory; the second combining the NIR data with radio data; and the third making use of the Spitzer/IRS MIR observations to measure the 9.8 μm silicate absorption feature using the Spoon method and PAHFIT (as the depth of this feature can be related to overall extinction). The silicate absorption over the bright, southern component of IRAS 12063-6259 is almost absent while the other methods find significant extinction. While such breakdowns of the relationship between the NIR extinction and the 9.8 μm absorption have been observed in molecular clouds, they have never been observed for HII regions. We then compare the PAH intensity variations in the Spitzer/IRS data after dereddening to those found in the original data. It was found that in most cases, the PAH band intensity variations persist even after dereddening, implying that extinction is not the main cause of the PAH band intensity variations.

preprint2012arXiv

Understanding and reducing statistical uncertainties in nebular abundance determinations

Whenever observations are compared to theories, an estimate of the uncertainties associated with the observations is vital if the comparison is to be meaningful. However, many determinations of temperatures, densities and abundances in photoionized nebulae do not quote the associated uncertainty. Those that do typically propagate the uncertainties using analytical techniques which rely on assumptions that generally do not hold. Motivated by this issue, we have developed NEAT (Nebular Empirical Analysis Tool), a new code for calculating chemical abundances in photoionized nebulae. The code carries out an analysis of lists of emission lines using long-established techniques to estimate the amount of interstellar extinction, calculate representative temperatures and densities, compute ionic abundances from both collisionally excited lines and recombination lines, and finally to estimate total elemental abundances using an ionization correction scheme. NEAT uses a Monte Carlo technique to robustly propagate uncertainties from line flux measurements through to the derived abundances. We show that for typical observational data, this approach is superior to analytic estimates of uncertainties. NEAT also accounts for the effect of upward biasing on measurements of lines with low signal to noise, allowing us to accurately quantify the effect of this bias on abundance determinations. We find not only that the effect can result in significant over-estimates of heavy element abundances derived from weak lines, but that taking it into account reduces the uncertainty of these abundance determinations. Finally, we investigate the effect of possible uncertainties in R, the ratio of selective to total extinction, on abundance determinations. We find that the uncertainty due to this parameter is negligible compared to the statistical uncertainties due to typical line flux measurement uncertainties.

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.

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.

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.