Showing 1–50 of 208 results for author: Barlow, M J

Deep JWST/NIRCam imaging of Supernova 1987A

Authors: Mikako Matsuura, M. Boyer, Richard G. Arendt, J. Larsson, C. Fransson, A. Rest, A. P. Ravi, S. Park, P. Cigan, T. Temim, E. Dwek, M. J. Barlow, P. Bouchet, G. Clayton, R. Chevalier, J. Danziger, J. De Buizer, I. De Looze, G. De Marchi, O. Fox, C. Gall, R. D. Gehrz, H. L. Gomez, R. Indebetouw, T. Kangas , et al. (24 additional authors not shown)

Abstract: JWST/NIRCam obtained high angular-resolution (0.05-0.1''), deep near-infrared 1--5 micron imaging of Supernova (SN) 1987A taken 35 years after the explosion. In the NIRCam images, we identify: 1) faint H2 crescents, which are emissions located between the ejecta and the equatorial ring, 2) a bar, which is a substructure of the ejecta, and 3) the bright 3-5 micron continuum emission exterior to the… ▽ More JWST/NIRCam obtained high angular-resolution (0.05-0.1''), deep near-infrared 1--5 micron imaging of Supernova (SN) 1987A taken 35 years after the explosion. In the NIRCam images, we identify: 1) faint H2 crescents, which are emissions located between the ejecta and the equatorial ring, 2) a bar, which is a substructure of the ejecta, and 3) the bright 3-5 micron continuum emission exterior to the equatorial ring. The emission of the remnant in the NIRCam 1-2.3 micron images is mostly due to line emission, which is mostly emitted in the ejecta and in the hot spots within the equatorial ring. In contrast, the NIRCam 3-5 micron images are dominated by continuum emission. In the ejecta, the continuum is due to dust, obscuring the centre of the ejecta. In contrast, in the ring and exterior to the ring, synchrotron emission contributes a substantial fraction to the continuum. Dust emission contributes to the continuum at outer spots and diffuse emission exterior to the ring, but little within the ring. This shows that dust cooling and destruction time scales are shorter than the synchrotron cooling time scale, and the time scale of hydrogen recombination in the ring is even longer than the synchrotron cooling time scale. With the advent of high sensitivity and high angular resolution images provided by JWST/NIRCam, our observations of SN 1987A demonstrate that NIRCam opens up a window to study particle-acceleration and shock physics in unprecedented details, probed by near-infrared synchrotron emission, building a precise picture of how a SN evolves. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: Accepted for publication in MNRAS; 18 pages

Emission lines due to ionizing radiation from a compact object in the remnant of Supernova 1987A

Authors: C. Fransson, M. J. Barlow, P. J. Kavanagh, J. Larsson, O. C. Jones, B. Sargent, M. Meixner, P. Bouchet, T. Temim, G. S. Wright, J. A. D. L. Blommaert, N. Habel, A. S. Hirschauer, J. Hjorth, L. Lenkić, T. Tikkanen, R. Wesson, A. Coulais, O. D. Fox, R. Gastaud, A. Glasse, J. Jaspers, O. Krause, R. M. Lau, O. Nayak , et al. (9 additional authors not shown)

Abstract: The nearby Supernova 1987A was accompanied by a burst of neutrino emission, which indicates that a compact object (a neutron star or black hole) was formed in the explosion. There has been no direct observation of this compact object. In this work, we observe the supernova remnant with JWST spectroscopy finding narrow infrared emission lines of argon and sulphur. The line emission is spatially unr… ▽ More The nearby Supernova 1987A was accompanied by a burst of neutrino emission, which indicates that a compact object (a neutron star or black hole) was formed in the explosion. There has been no direct observation of this compact object. In this work, we observe the supernova remnant with JWST spectroscopy finding narrow infrared emission lines of argon and sulphur. The line emission is spatially unresolved and blueshifted in velocity relative to the supernova rest frame. We interpret the lines as gas illuminated by a source of ionizing photons located close to the center of the expanding ejecta. Photoionization models show that the line ratios are consistent with ionization by a cooling neutron star or pulsar wind nebula. The velocity shift could be evidence for a neutron star natal kick. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: Authors version of manuscript published in Science on 22 Feb 2024

Journal ref: SCIENCE 22 Feb 2024 Vol 383, Issue 6685 pp. 898-903

JWST MIRI Imager Observations of Supernova SN 1987A

Authors: P. Bouchet, R. Gastaud, A. Coulais, M. J. Barlow, C. Fransson, P. J. Kavanagh, J. Larsson, T. Temim, O. C. Jones, A. S. Hirschauer, T. Tikkanen, J. A. D. L. Blommaert, O. D. Fox, A. Glasse, N. Habel, J. Hjorth, J. Jaspers, O. Krause, R. M. Lau, L. Lenkić, M. Meixner, O. Nayak, A. Rest, B. Sargent, R. Wesson , et al. (9 additional authors not shown)

Abstract: There exist very few mid-infrared (IR) observations of supernovae (SNe) in general. Therefore, SN 1987A, the closest visible SN in 400 years, gives us the opportunity to explore the mid-IR properties of SNe, the dust in their ejecta and surrounding medium, and to witness the birth of a SN remnant (SNR). The James Webb Space Telescope (JWST), with its high spatial resolution and extreme sensitivity… ▽ More There exist very few mid-infrared (IR) observations of supernovae (SNe) in general. Therefore, SN 1987A, the closest visible SN in 400 years, gives us the opportunity to explore the mid-IR properties of SNe, the dust in their ejecta and surrounding medium, and to witness the birth of a SN remnant (SNR). The James Webb Space Telescope (JWST), with its high spatial resolution and extreme sensitivity, gives a new view on these issues. We report on the first imaging observations obtained with the Mid-InfraRed Instrument (MIRI). We build temperature maps and discuss the morphology of the nascent SNR. Our results show that the temperatures in the equatorial ring (ER) are quite non-uniform. This could be due to dust destruction in some parts of the ring, as had been assumed in some previous works. We show that the IR emission extends beyond the ER, illustrating the fact that the shock wave has now passed through this ring to affect the circumstellar medium on a larger scale. Finally, while sub-mm Atacama Large Millimeter Array (ALMA) observations have hinted at the location of the compact remnant of SN 1987A, we note that our MIRI data have found no such evidence. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Comments: 19 pages, 19 figures, 2 tables; Accepted for publication in the Astrophysical Journal (February 2, 2024)

From total destruction to complete survival: Dust processing at different evolutionary stages in the supernova remnant Cassiopeia A

Authors: Florian Kirchschlager, Nina Sartorio, Ilse De Looze, M. J. Barlow, Franziska Schmidt, Felix Priestley

Abstract: The expanding ejecta of supernova remnants (SNRs) are believed to form dust in dense clumps of gas. Before the dust can be expelled into the interstellar medium and contribute to the interstellar dust budget, it has to survive the reverse shock that is generated through the interaction of the preceding supernova blast wave with the surrounding medium. The conditions under which the reverse shock h… ▽ More The expanding ejecta of supernova remnants (SNRs) are believed to form dust in dense clumps of gas. Before the dust can be expelled into the interstellar medium and contribute to the interstellar dust budget, it has to survive the reverse shock that is generated through the interaction of the preceding supernova blast wave with the surrounding medium. The conditions under which the reverse shock hits the clumps change with remnant age and define the dust survival rate. To study the dust destruction in the SNR Cassiopeia A, we conduct magnetohydrodynamical simulations of the evolution of a supernova blast wave and of the reverse shock. In a second step we use these evolving conditions to model clumps that are disrupted by the reverse shock at different remnant ages. Finally, we compute the amount of dust that is destroyed by the impact of the reverse shock. We find that most of the dust in the SNR is hit by the reverse shock within the first 350 yr after the SN explosion. While the dust destruction in the first 200 yr is almost complete, we expect greater dust survival rates at later times and almost total survival for clumps that are first impacted at ages beyond 1000 yr. Integrated over the entire evolution of the SNR, the dust mass shows the lowest survival fraction (17 per cent) for the smallest grains (1 nm) and the highest survival fraction (28 per cent) for the largest grains (1000 nm). △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: Accepted by MNRAS. Author accepted manuscript. Accepted on 31/01/2024. Deposited on 01/02/2024. 12 pages

A JWST Survey of the Supernova Remnant Cassiopeia A

Authors: Dan Milisavljevic, Tea Temim, Ilse De Looze, Danielle Dickinson, J. Martin Laming, Robert Fesen, John C. Raymond, Richard G. Arendt, Jacco Vink, Bettina Posselt, George G. Pavlov, Ori D. Fox, Ethan Pinarski, Bhagya Subrayan, Judy Schmidt, William P. Blair, Armin Rest, Daniel Patnaude, Bon-Chul Koo, Jeonghee Rho, Salvatore Orlando, Hans-Thomas Janka, Moira Andrews, Michael J. Barlow, Adam Burrows , et al. (21 additional authors not shown)

Abstract: We present initial results from a JWST survey of the youngest Galactic core-collapse supernova remnant Cassiopeia A (Cas A), made up of NIRCam and MIRI imaging mosaics that map emission from the main shell, interior, and surrounding circumstellar/interstellar material (CSM/ISM). We also present four exploratory positions of MIRI/MRS IFU spectroscopy that sample ejecta, CSM, and associated dust fro… ▽ More We present initial results from a JWST survey of the youngest Galactic core-collapse supernova remnant Cassiopeia A (Cas A), made up of NIRCam and MIRI imaging mosaics that map emission from the main shell, interior, and surrounding circumstellar/interstellar material (CSM/ISM). We also present four exploratory positions of MIRI/MRS IFU spectroscopy that sample ejecta, CSM, and associated dust from representative shocked and unshocked regions. Surprising discoveries include: 1) a web-like network of unshocked ejecta filaments resolved to 0.01 pc scales exhibiting an overall morphology consistent with turbulent mixing of cool, low-entropy matter from the progenitor's oxygen layer with hot, high-entropy matter heated by neutrino interactions and radioactivity, 2) a thick sheet of dust-dominated emission from shocked CSM seen in projection toward the remnant's interior pockmarked with small (approximately one arcsecond) round holes formed by knots of high-velocity ejecta that have pierced through the CSM and driven expanding tangential shocks, 3) dozens of light echoes with angular sizes between 0.1 arcsecond to 1 arcminute reflecting previously unseen fine-scale structure in the ISM. NIRCam observations place new upper limits on infrared emission from the neutron star in Cas A's center and tightly constrain scenarios involving a possible fallback disk. These JWST survey data and initial findings help address unresolved questions about massive star explosions that have broad implications for the formation and evolution of stellar populations, the metal and dust enrichment of galaxies, and the origin of compact remnant objects. △ Less

Submitted 10 June, 2024; v1 submitted 4 January, 2024; originally announced January 2024.

Comments: 27 pages, 10 figures, now published in ApJL

Journal ref: ApJL 965 (2024) L27 (21pp)

Oxygen abundance of gamma Vel from [O III] 88um Herschel/PACS spectroscopy

Authors: Paul A Crowther, M J Barlow, P Royer, D J Hillier, J M Bestenlehner, P W Morris, R Wesson

Abstract: We present Herschel PACS spectroscopy of the [O III] 88.4um fine-structure line in the nearby WC8+O binary system gamma Vel to determine its oxygen abundance. The critical density of this line corresponds to several 10^5 R* such that it is spatially extended in PACS observations at the 336 pc distance to gamma Vel. Two approaches are used, the first involving a detailed stellar atmosphere analysis… ▽ More We present Herschel PACS spectroscopy of the [O III] 88.4um fine-structure line in the nearby WC8+O binary system gamma Vel to determine its oxygen abundance. The critical density of this line corresponds to several 10^5 R* such that it is spatially extended in PACS observations at the 336 pc distance to gamma Vel. Two approaches are used, the first involving a detailed stellar atmosphere analysis of gamma Vel using CMFGEN, extending to Ne ~ 10^0 cm^-3 in order to fully sample the line formation region of [O III] 88.4um. The second approach involves the analytical model introduced by Barlow et al. and revised by Dessart et al, additionally exploiting ISO LWS spectroscopy of [O III] 51.8um. We obtain higher luminosities for the WR and O components of gamma Vel with respect to De Marco et al, log L/L_sun = 5.31 and 5.56, respectively, due to the revised (higher) interferometric distance. We obtain an oxygen mass fraction of X_O = 1.0+/- 0.3% for an outer wind volume filling factor of f = 0.5+/-0.25, favouring either standard or slightly reduced Kunz et al. rates for the ^12C(alpha, gamma)^16O reaction from comparison with BPASS binary population synthesis models. We also revisit neon and sulphur abundances in the outer wind of gamma Vel from ISO SWS spectroscopy of [S IV] 10.5um [Ne II] 12.8um and [Ne III] 15.5um. The inferred neon abundance X_Ne = 2.0-0.6+0.4%, is in excellent agreement with BPASS predictions, while the.sulphur abundance of X_S = 0.04 +/- 0.01% agrees with the solar abundance, as expected for unprocessed elements. △ Less

Submitted 10 January, 2024; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: 12 pages, 9 figures, plus Appendices, submitted to MNRAS, revised in response to referee report

Serendipitous detection of the dusty Type IIL SN 1980K with JWST/MIRI

Authors: Szanna Zsíros, Tamás Szalai, Ilse De Looze, Arkaprabha Sarangi, Melissa Shahbandeh, Ori D. Fox, Tea Temim, Dan Milisavljevic, Schuyler D. Van Dyk, Nathan Smith, Alexei V. Filippenko, Thomas G. Brink, WeiKang Zheng, Luc Dessart, Jacob Jencson, Joel Johansson, Justin Pierel, Armin Rest, Samaporn Tinyanont, Maria Niculescu-Duvaz, M. J. Barlow, Roger Wesson, Jennifer Andrews, Geoff Clayton, Kishalay De , et al. (17 additional authors not shown)

Abstract: We present mid-infrared (mid-IR) imaging of the Type IIL supernova (SN) 1980K with the James Webb Space Telescope (JWST) more than 40 yr post-explosion. SN 1980K, located in the nearby ($D\approx7$ Mpc) "SN factory" galaxy NGC 6946, was serendipitously captured in JWST/MIRI images taken of the field of SN 2004et in the same galaxy. SN 1980K serves as a promising candidate for studying the transiti… ▽ More We present mid-infrared (mid-IR) imaging of the Type IIL supernova (SN) 1980K with the James Webb Space Telescope (JWST) more than 40 yr post-explosion. SN 1980K, located in the nearby ($D\approx7$ Mpc) "SN factory" galaxy NGC 6946, was serendipitously captured in JWST/MIRI images taken of the field of SN 2004et in the same galaxy. SN 1980K serves as a promising candidate for studying the transitional phase between young SNe and older SN remnants and also provides a great opportunity to investigate its the close environment. SN 1980K can be identified as a clear and bright point source in all eight MIRI filters from F560W up to F2550W. We fit analytical dust models to the mid-IR spectral energy distribution that reveal a large amount ($M_d \approx 0.002 {M}_{\odot}$) of Si-dominated dust at $T_{dust}\approx 150$ K (accompanied by a hotter dust/gas component), and also computed numerical SED dust models. Radiative transfer modeling of a late-time optical spectrum obtained recently with Keck discloses that an even larger ($\sim 0.24-0.58~{M}_{\odot}$) amount of dust is needed in order for selective extinction to explain the asymmetric line profile shapes observed in SN 1980K. As a conclusion, with JWST, we may see i) pre-existing circumstellar dust heated collisionally (or, partly radiatively), analogous to the equatorial ring of SN 1987A, or ii) the mid-IR component of the presumed newly-formed dust, accompanied by much more colder dust present in the ejecta (as suggested by the late-time the optical spectra). △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 14 pages, 9 figures

Evidence for late-time dust formation in the ejecta of supernova SN~1995N from emission-line asymmetries

Authors: R. Wesson, A. M. Bevan, M. J. Barlow, I. De Looze, M. Matsuura, G. Clayton, J. Andrews

Abstract: We present a study of the dust associated with the core-collapse supernova SN~1995N. Infrared emission detected 14--15 years after the explosion was previously attributed to thermally echoing circumstellar material associated with the SN progenitor. We argue that this late-time emission is unlikely to be an echo, and is more plausibly explained by newly formed dust in the supernova ejecta, indirec… ▽ More We present a study of the dust associated with the core-collapse supernova SN~1995N. Infrared emission detected 14--15 years after the explosion was previously attributed to thermally echoing circumstellar material associated with the SN progenitor. We argue that this late-time emission is unlikely to be an echo, and is more plausibly explained by newly formed dust in the supernova ejecta, indirectly heated by the interaction between the ejecta and the CSM. Further evidence in support of this scenario comes from emission line profiles in spectra obtained 22 years after the explosion; these are asymmetric, showing greater attenuation on the red wing, consistent with absorption by dust within the expanding ejecta. The spectral energy distribution and emission line profiles at epochs later than $\sim$5000 days are both consistent with the presence of about 0.4~M$_\odot$ of amorphous carbon dust. The onset of dust formation is apparent in archival optical spectra, taken between 700 and 1700 days after the assumed explosion date. As this is considerably later than most other instances where the onset of dust formation has been detected, we argue that the explosion date must be later than previously assumed. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Comments: 14 pages, 16 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

JWST observations of the Ring Nebula (NGC 6720): I. Imaging of the rings, globules, and arcs

Authors: R. Wesson, Mikako Matsuura, Albert A. Zijlstra, Kevin Volk, Patrick J. Kavanagh, Guillermo García-Segura, I. McDonald, Raghvendra Sahai, M. J. Barlow, Nick L. J. Cox, Jeronimo Bernard-Salas, Isabel Aleman, Jan Cami, Nicholas Clark, Harriet L. Dinerstein, K. Justtanont, Kyle F. Kaplan, A. Manchado, Els Peeters, Griet C. Van de Steene, Peter A. M. van Hoof

Abstract: We present JWST images of the well-known planetary nebula NGC 6720 (the Ring Nebula), covering wavelengths from 1.6$μ$m to 25 $μ$m. The bright shell is strongly fragmented with some 20 000 dense globules, bright in H$_2$, with a characteristic diameter of 0.2 arcsec and density $n_{\rm H} \sim 10^5$-$10^6$ cm$^{-3}$. The shell contains a thin ring of polycyclic aromatic hydrocarbon (PAH) emission.… ▽ More We present JWST images of the well-known planetary nebula NGC 6720 (the Ring Nebula), covering wavelengths from 1.6$μ$m to 25 $μ$m. The bright shell is strongly fragmented with some 20 000 dense globules, bright in H$_2$, with a characteristic diameter of 0.2 arcsec and density $n_{\rm H} \sim 10^5$-$10^6$ cm$^{-3}$. The shell contains a thin ring of polycyclic aromatic hydrocarbon (PAH) emission. H$_2$ is found throughout the shell and in the halo. H$_2$ in the halo may be located on the swept-up walls of a biconal polar flow. The central cavity is shown to be filled with high ionization gas and shows two linear structures. The central star is located 2 arcsec from the emission centroid of the cavity and shell. Linear features (`spikes') extend outward from the ring, pointing away from the central star. Hydrodynamical simulations are shown which reproduce the clumping and possibly the spikes. Around ten low-contrast, regularly spaced concentric arc-like features are present; they suggest orbital modulation by a low-mass companion with a period of about 280 yr. A previously known much wider companion is located at a projected separation of about 15 000 au; we show that it is an M2-M4 dwarf. The system is therefore a triple star. These features, including the multiplicity, are similar to those seen in the Southern Ring Nebula (NGC 3132) and may be a common aspect of such nebulae. △ Less

Submitted 21 August, 2023; v1 submitted 17 August, 2023; originally announced August 2023.

Comments: 25 pages, 23 figures. Submitted to Monthly Notices of the Royal Astronomical Society. Corrected typos in metadata

Ejecta, Rings, and Dust in SN 1987A with JWST MIRI/MRS

Authors: O. C. Jones, P. J. Kavanagh, M. J. Barlow, T. Temim, C. Fransson, J. Larsson, J. A. D. L. Blommaert, M. Meixner, R. M. Lau, B. Sargent, P. Bouchet, J. Hjorth, G. S. Wright, A. Coulais, O. D. Fox, R. Gastaud, A. Glasse, N. Habel, A. S. Hirschauer, J. Jaspers, O. Krause, Lenkić, O. Nayak, A. Rest, T. Tikkanen , et al. (9 additional authors not shown)

Abstract: Supernova (SN) 1987A is the nearest supernova in $\sim$400 years. Using the {\em JWST} MIRI Medium Resolution Spectrograph, we spatially resolved the ejecta, equatorial ring (ER) and outer rings in the mid-infrared 12,927 days after the explosion. The spectra are rich in line and dust continuum emission, both in the ejecta and the ring. Broad emission lines (280-380~km~s$^{-1}$ FWHM) seen from all… ▽ More Supernova (SN) 1987A is the nearest supernova in $\sim$400 years. Using the {\em JWST} MIRI Medium Resolution Spectrograph, we spatially resolved the ejecta, equatorial ring (ER) and outer rings in the mid-infrared 12,927 days after the explosion. The spectra are rich in line and dust continuum emission, both in the ejecta and the ring. Broad emission lines (280-380~km~s$^{-1}$ FWHM) seen from all singly-ionized species originate from the expanding ER, with properties consistent with dense post-shock cooling gas. Narrower emission lines (100-170~km~s$^{-1}$ FWHM) are seen from species originating from a more extended lower-density component whose high ionization may have been produced by shocks progressing through the ER, or by the UV radiation pulse associated with the original supernova event. The asymmetric east-west dust emission in the ER has continued to fade, with constant temperature, signifying a reduction in dust mass. Small grains in the ER are preferentially destroyed, with larger grains from the progenitor surviving the transition from SN into SNR. The ER is fit with a single set of optical constants, eliminating the need for a secondary featureless hot dust component. We find several broad ejecta emission lines from [Ne~{\sc ii}], [Ar~{\sc ii}], [Fe~{\sc ii}], and [Ni~{\sc ii}]. With the exception of [Fe~{\sc ii}]~25.99$μ$m, these all originate from the ejecta close to the ring and are likely being excited by X-rays from the interaction. The [Fe~{\sc ii}]~5.34$μ$m to 25.99$μ$m line ratio indicates a temperature of only a few hundred K in the inner core, consistent with being powered by ${}^{44}$Ti decay. △ Less

Submitted 29 February, 2024; v1 submitted 13 July, 2023; originally announced July 2023.

Comments: 27 pages, 16 figures, 4 tables. Accepted ApJ

Far-infrared Polarization of the Supernova Remnant Cassiopeia A with SOFIA HAWC+

Authors: Jeonghee Rho, Aravind P. Ravi, Le Ngoc Tram, Thiem Hoang, Jérémy Chastenet, Matthew Millard, Michael J. Barlow, Ilse De Looze, Haley L. Gomez, Florian Kirchschlager, Loretta Dunne

Abstract: We present polarization observations of the young supernova remnant (SNR) Cas A using the High-resolution Airborne Wideband Camera-Plus (HAWC+) instrument onboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The polarization map at 154 microns reveals dust grains with strong polarization fractions (5 - 30 percent), supporting previous measurements made over a smaller region of the… ▽ More We present polarization observations of the young supernova remnant (SNR) Cas A using the High-resolution Airborne Wideband Camera-Plus (HAWC+) instrument onboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The polarization map at 154 microns reveals dust grains with strong polarization fractions (5 - 30 percent), supporting previous measurements made over a smaller region of the remnant at 850 microns. The 154 microns emission and the polarization signal is coincident with a region of cold dust observed in the southeastern shell and in the unshocked central ejecta. The highly polarized far-IR emission implies the grains are large (greater than 0.14 microns) and silicate-dominated. The polarization level varies across the SNR, with an inverse correlation between the polarization degree and the intensity and smaller polarization angle dispersion for brighter SNR emission. Stronger polarization is detected between the bright structures. This may result from a higher collision rate between the gas and dust producing a lower grain alignment efficiency where the gas density is higher. We use the dust emission to provide an estimate of the magnetic field strength in Cas A using the Davis-Chandrasekhar-Fermi method. The high polarization level is direct evidence that grains are highly elongated and strongly aligned with the magnetic field of the SNR. The dust mass from the polarized region is 0.14+-0.04 Msun, a lower limit of the amount of dust present within the ejecta of Cas A. This result strengthens the hypothesis that core-collapse SNe are an important contributor to the dust mass in high redshift galaxies. △ Less

Submitted 13 April, 2023; originally announced April 2023.

Comments: MNRAS, accepted (18 pages with 14 figures)

JWST NIRSpec observations of Supernova 1987A -- from the inner ejecta to the reverse shock

Authors: J. Larsson, C. Fransson, B. Sargent, O. C. Jones, M. J. Barlow, P. Bouchet, M. Meixner, J. A. D. L. Blommaert, A. Coulais, O. D. Fox, R. Gastaud, A. Glasse, N. Habel, A. S. Hirschauer, J. Hjorth, J. Jaspers, P. J. Kavanagh, O. Krause, R. M. Lau, L. Lenkic, O. Nayak, A. Rest, T. Temim, T. Tikkanen, R. Wesson , et al. (1 additional authors not shown)

Abstract: We present initial results from JWST NIRSpec integral field unit observations of the nearby Supernova (SN) 1987A. The observations provide the first spatially-resolved spectroscopy of the ejecta and equatorial ring (ER) over the 1-5 μm range. We construct 3D emissivity maps of the [Fe I] 1.443 μm line from the inner ejecta and the He I 1.083 μm line from the reverse shock (RS), where the former pr… ▽ More We present initial results from JWST NIRSpec integral field unit observations of the nearby Supernova (SN) 1987A. The observations provide the first spatially-resolved spectroscopy of the ejecta and equatorial ring (ER) over the 1-5 μm range. We construct 3D emissivity maps of the [Fe I] 1.443 μm line from the inner ejecta and the He I 1.083 μm line from the reverse shock (RS), where the former probes the explosion geometry and the latter traces the structure of the circumstellar medium. We also present a model for the integrated spectrum of the ejecta. The [Fe I] 3D map reveals a highly-asymmetric morphology resembling a broken dipole, dominated by two large clumps with velocities of ~2300 km/s. We also find evidence that the Fe-rich inner ejecta have started to interact with the RS. The RS surface traced by the He I line extends from just inside the ER to higher latitudes on both sides of the ER with a half-opening angle ~45 degrees, forming a bubble-like structure. The spectral model for the ejecta allows us to identify the many emission lines, including numerous H_2 lines. We find that the H_2 is most likely excited by far-UV emission, while the metal lines ratios are consistent with a combination of collisional excitation and recombination in the low-temperature ejecta. We also find several high-ionization coronal lines from the ER, requiring a temperature > 2 \times 10^6 K. △ Less

Submitted 16 May, 2023; v1 submitted 7 February, 2023; originally announced February 2023.

Comments: Accepted for publication in ApJL

Mid-infrared imaging of Supernova 1987A

Authors: Mikako Matsuura, Roger Wesson, Richard G. Arendt, Eli Dwek, James M. De Buizer, John Danziger, Patrice Bouchet, M. J. Barlow, Phil Cigan, Haley L. Gomez, Jeonghee Rho, Margaret Meixner

Abstract: At a distance of 50 kpc, Supernova 1987A is an ideal target to study how a young supernova (SN) evolves in time. Its equatorial ring, filled with material expelled from the progenitor star about 20,000 years ago, has been engulfed with SN blast waves. Shocks heat dust grains in the ring, emitting their energy at mid-infrared (IR) wavelengths We present ground-based 10--18$μ$m monitoring of the rin… ▽ More At a distance of 50 kpc, Supernova 1987A is an ideal target to study how a young supernova (SN) evolves in time. Its equatorial ring, filled with material expelled from the progenitor star about 20,000 years ago, has been engulfed with SN blast waves. Shocks heat dust grains in the ring, emitting their energy at mid-infrared (IR) wavelengths We present ground-based 10--18$μ$m monitoring of the ring of SN 1987A from day 6067 to 12814 at a resolution of 0.5", together with SOFIA photometry at 10-30 $μ$m. The IR images in the 2000's (day 6067-7242) showed that the shocks first began brightening the east side of the ring. Later, our mid-IR images from 2017 to 2022 (day 10952-12714) show that dust emission is now fading in the east, while it has brightened on the west side of the ring. Because dust grains are heated in the shocked plasma, which can emit X-rays, the IR and X-ray brightness ratio represent shock diagnostics. Until 2007 the IR to X-ray brightness ratio remained constant over time, and during this time shocks seemed to be largely influencing the east side of the ring. However, since then, the IR to X-ray ratio has been declining, due to increased X-ray brightness. Whether the declining IR brightness is because of dust grains being destroyed or being cooled in the post-shock regions will require more detailed modelling. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: MNRAS accepted

Dust survival rates in clumps passing through the Cas A reverse shock -- II. The impact of magnetic fields

Authors: Florian Kirchschlager, Franziska D. Schmidt, M. J. Barlow, Ilse De Looze, Nina S. Sartorio

Abstract: Dust grains form in the clumpy ejecta of core-collapse supernovae where they are subject to the reverse shock, which is able to disrupt the clumps and destroy the grains. Important dust destruction processes include thermal and kinetic sputtering as well as fragmentation and grain vaporization. In the present study, we focus on the effect of magnetic fields on the destruction processes. We have pe… ▽ More Dust grains form in the clumpy ejecta of core-collapse supernovae where they are subject to the reverse shock, which is able to disrupt the clumps and destroy the grains. Important dust destruction processes include thermal and kinetic sputtering as well as fragmentation and grain vaporization. In the present study, we focus on the effect of magnetic fields on the destruction processes. We have performed magneto-hydrodynamical simulations using AstroBEAR to model a shock wave interacting with an ejecta clump. The dust transport and destruction fractions are computed using our post-processing code Paperboats in which the acceleration of grains due to the magnetic field and a procedure that allows partial grain vaporization have been newly implemented. For the oxygen-rich supernova remnant Cassiopeia A we found a significantly lower dust survival rate when magnetic fields are aligned perpendicular to the shock direction compared to the non-magnetic case. For a parallel field alignment, the destruction is also enhanced but at a lower level. The survival fractions depend sensitively on the gas density contrast between the clump and the ambient medium and on the grain sizes. For a low-density contrast of $100$, e.g., $5\,$nm silicate grains are completely destroyed while the survival fraction of $1\,μ$m grains is $86\,$per cent. For a high-density contrast of $1000$, $95\,$per cent of the $5\,$nm grains survive while the survival fraction of $1\,μ$m grains is $26\,$per cent. Alternative clump sizes or dust materials (carbon) have non-negligible effects on the survival rate but have a lower impact compared to density contrast, magnetic field strength, and grain size. △ Less

Submitted 16 February, 2023; v1 submitted 13 October, 2022; originally announced October 2022.

Comments: Accepted by MNRAS. Author accepted manuscript. Accepted on 23/01/2023. 24 pages, 21 Figures

Properties of shocked dust grains in supernova remnants

Authors: F. D. Priestley, H. Chawner, M. J. Barlow, I. De Looze, H. L. Gomez, M. Matsuura

Abstract: Shockwaves driven by supernovae both destroy dust and reprocess the surviving grains, greatly affecting the resulting dust properties of the interstellar medium (ISM). While these processes have been extensively studied theoretically, observational constraints are limited. We use physically-motivated models of dust emission to fit the infrared (IR) spectral energy distributions of seven Galactic s… ▽ More Shockwaves driven by supernovae both destroy dust and reprocess the surviving grains, greatly affecting the resulting dust properties of the interstellar medium (ISM). While these processes have been extensively studied theoretically, observational constraints are limited. We use physically-motivated models of dust emission to fit the infrared (IR) spectral energy distributions of seven Galactic supernova remnants, allowing us to determine the distribution of dust mass between diffuse and dense gas phases, and between large and small grain sizes. We find that the dense ($\sim 10^3 \,{\rm cm}^{-3}$), relatively cool ($\sim 10^3 \, {\rm K}$) gas phase contains $>90\%$ of the dust mass, making the warm dust located in the X-ray emitting plasma ($\sim 1 \,{\rm cm}^{-3}$/$10^6 \, {\rm K}$) a negligible fraction of the total, despite dominating the mid-IR emission. The ratio of small ($\lesssim 10 \, {\rm nm}$) to large ($\gtrsim 0.1 \, {\rm μm}$) grains in the cold component is consistent with that in the ISM, and possibly even higher, whereas the hot phase is almost entirely devoid of small grains. This suggests that grain shattering, which processes large grains into smaller ones, is ineffective in the low-density gas, contrary to model predictions. Single-phase models of dust destruction in the ISM, which do not account for the existence of the cold swept-up material containing most of the dust mass, are likely to greatly overestimate the rate of dust destruction by supernovae. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 13 pages, 12 figures. MNRAS accepted

SOFIA/HAWC+ observations of the Crab Nebula: dust properties from polarised emission

Authors: Jérémy Chastenet, Ilse De Looze, Brandon S. Hensley, Bert Vandenbroucke, Mike J. Barlow, Jeonghee Rho, Aravind P. Ravi, Haley L. Gomez, Florian Kirchschlager, Juan Macías-Pérez, Mikako Matsuura, Kate Pattle, Nicolas Ponthieu, Felix D. Priestley, Monica Relaño, Alessia Ritacco, Roger Wesson

Abstract: Supernova remnants (SNRs) are well-recognised dust producers, but their net dust production rate remains elusive due to uncertainties in grain properties that propagate into observed dust mass uncertainties, and determine how efficiently these grains are processed by reverse shocks. In this paper, we present a detection of polarised dust emission in the Crab pulsar wind nebula, the second SNR with… ▽ More Supernova remnants (SNRs) are well-recognised dust producers, but their net dust production rate remains elusive due to uncertainties in grain properties that propagate into observed dust mass uncertainties, and determine how efficiently these grains are processed by reverse shocks. In this paper, we present a detection of polarised dust emission in the Crab pulsar wind nebula, the second SNR with confirmed polarised dust emission after Cassiopeia A. We constrain the bulk composition of the dust with new SOFIA/HAWC+ polarimetric data in band C 89 um and band D 154 um. After correcting for synchrotron polarisation, we report dust polarisation fractions ranging between 3.7-9.6 per cent and 2.7-7.6 per cent in three individual dusty filaments at 89 and 154 um, respectively. The detected polarised signal suggests the presence of large (> 0.05-0.1 um) grains in the Crab Nebula. With the observed polarisation, and polarised and total fluxes, we constrain the temperatures and masses of carbonaceous and silicate grains. We find that the carbon-rich grain mass fraction varies between 12 and 70 per cent, demonstrating that carbonaceous and silicate grains co-exist in this SNR. Temperatures range from 40 K to 70 K and from 30 K to 50 K for carbonaceous and silicate grains, respectively. Dust masses range from 10^{-4} Msol to 10^{-2} Msol for carbonaceous grains and to 10^{-1} Msol for silicate grains, in three individual regions. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: Accepted for publication in MNRAS. 16 pages, 8 figures + 2 Appendices, 2 figures

Quantifying the dust in SN 2012aw and iPTF14hls with ORBYTS

Authors: Maria Niculescu-Duvaz, M. J. Barlow, W. Dunn, A. Bevan, Omar Ahmed, David Arkless, Jon Barker, Sidney Bartolotta, Liam Brockway, Daniel Browne, Ubaid Esmail, Max Garner, Wiktoria Guz, Scarlett King, Hayri Kose, Madeline Lampstaes-Capes, Joseph Magen, Nicole Morrison, Kyaw Oo, Balvinder Paik, Joanne Primrose, Danny Quick, Anais Radeka, Anthony Rodney, Eleanor Sandeman , et al. (10 additional authors not shown)

Abstract: 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 p… ▽ More 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. △ Less

Submitted 4 January, 2023; v1 submitted 1 June, 2022; originally announced June 2022.

Comments: Accepted by MNRAS on 28/11/2022, 10 pages, 6 figures. Author accepted manuscript

Dust masses for a large sample of core-collapse supernovae from optical emission line asymmetries: dust formation on 30-year timescales

Authors: Maria Niculescu-Duvaz, Michael J Barlow, Antonia Bevan, Roger Wesson, Danny Milisavljevic, Ilse De Looze, Geoff C. Clayton, Kelsie Krafton, Mikako Matsuura, Ryan Brady

Abstract: Modelling the red-blue asymmetries seen in the broad emission lines of core-collapse supernovae (CCSNe) is a powerful technique to quantify total dust mass formed in the ejecta at late times ($>5$ years after outburst) when ejecta dust temperatures become too low to be detected by mid-IR instruments. Following our success in using the Monte Carlo radiative transfer code DAMOCLES to measure the dus… ▽ More Modelling the red-blue asymmetries seen in the broad emission lines of core-collapse supernovae (CCSNe) is a powerful technique to quantify total dust mass formed in the ejecta at late times ($>5$ years after outburst) when ejecta dust temperatures become too low to be detected by mid-IR instruments. Following our success in using the Monte Carlo radiative transfer code DAMOCLES to measure the dust mass evolution in SN~1987A and other CCSNe, we present the most comprehensive sample of dust mass measurements yet made with DAMOCLES, for CCSNe aged between four and sixty years after outburst. Our sample comprises of multi-epoch late-time optical spectra taken with the Gemini GMOS and VLT X-Shooter spectrographs, supplemented by archival spectra. For the fourteen CCSNe that we have modelled, we confirm a dust mass growth with time that can be fit by a sigmoid curve which is found to saturate beyond an age of $\sim30$ years, at a mass of 0.23$^{+0.17}_{-0.12}$ M$_\odot$. An expanded sample including dust masses found in the literature for a further eleven CCSNe and six CCSN remnants, the dust mass at saturation is found to be 0.42$^{+0.09}_{-0.05}$~M$_\odot$. Uncertainty limits for our dust masses were determined from a Bayesian analysis using the affine invariant Markov Chain Monte Carlo ensemble sampler emcee with DAMOCLES. The best-fitting line profile models for our sample all required grain radii between 0.1 and 0.5 $μ$m. Our results are consistent with CCSNe forming enough dust in their ejecta to significantly contribute to the dust budget of the Universe. △ Less

Submitted 9 June, 2022; v1 submitted 29 April, 2022; originally announced April 2022.

Comments: Accepted by MNRAS, 32 pages plus Appendices. Author accepted manuscript. Accepted on 08/06/2022. Deposited on 09/06/2022

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

Authors: Mikako Matsuura, Victoria Ayley, Hannah Chawner, M. D. Filipovic, Warren Reid, F. D. Priestley, Andy Rigby, M. J. Barlow, Haley E. Gomez

Abstract: 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 wi… ▽ More 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. △ Less

Submitted 16 March, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: accepted for publication in MNRAS; 33 pages

Dust destruction and survival in the Cassiopeia A reverse shock

Authors: F. D. Priestley, M. Arias, M. J. Barlow, I. De Looze

Abstract: Core-collapse supernovae (CCSNe) produce large ($\gtrsim 0.1 \, {\rm M}_\odot$) masses of dust, and are potentially the primary source of dust in the Universe, but much of this dust may be destroyed before reaching the interstellar medium. Cassiopeia A (Cas A) is the only supernova remnant where an observational measurement of the dust destruction efficiency in the reverse shock is possible at pre… ▽ More Core-collapse supernovae (CCSNe) produce large ($\gtrsim 0.1 \, {\rm M}_\odot$) masses of dust, and are potentially the primary source of dust in the Universe, but much of this dust may be destroyed before reaching the interstellar medium. Cassiopeia A (Cas A) is the only supernova remnant where an observational measurement of the dust destruction efficiency in the reverse shock is possible at present. We determine the pre- and post-shock dust masses in Cas A using a substantially improved dust emission model. In our preferred models, the unshocked ejecta contains $0.6-0.8 \, {\rm M}_\odot$ of $0.1 \, {\rm μm}$ silicate grains, while the post-shock ejecta has $0.02-0.09 \, {\rm M}_\odot$ of $5-10 {\, {\rm nm}}$ grains in dense clumps, and $2 \times 10^{-3} \, {\rm M}_\odot$ of $0.1 \, {\rm μm}$ grains in the diffuse X-ray emitting shocked ejecta. The implied dust destruction efficiency is $74-94 \%$ in the clumps and $92-98 \%$ overall, giving Cas A a final dust yield of $0.05-0.30 \, {\rm M}_\odot$. If the unshocked ejecta grains are larger than $0.1 \, {\rm μm}$, the dust masses are higher, the destruction efficiencies are lower, and the final yield may exceed $0.5 \, {\rm M}_\odot$. As Cas A has a dense circumstellar environment and thus a much stronger reverse shock than is typical, the average dust destruction efficiency across all CCSNe is likely to be lower, and the average dust yield higher. This supports a mostly-stellar origin for the cosmic dust budget. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: 10 pages, 9 figures, MNRAS accepted

The impact of metallicity-dependent dust destruction on the dust-to-metals ratio in galaxies

Authors: F. D. Priestley, I. De Looze, M. J. Barlow

Abstract: The ratio of the mass of interstellar dust to the total mass of metals (the dust-to-metals/DTM ratio) tends to increase with metallicity. This can be explained by the increasing efficiency of grain growth in the interstellar medium (ISM) at higher metallicities, with a corollary being that the low DTM ratios seen at low metallicities are due to inefficient stellar dust production. This interpretat… ▽ More The ratio of the mass of interstellar dust to the total mass of metals (the dust-to-metals/DTM ratio) tends to increase with metallicity. This can be explained by the increasing efficiency of grain growth in the interstellar medium (ISM) at higher metallicities, with a corollary being that the low DTM ratios seen at low metallicities are due to inefficient stellar dust production. This interpretation assumes that the efficiency of dust destruction in the ISM is constant, whereas it might be expected to increase at low metallicity; the decreased cooling efficiency of low-metallicity gas should result in more post-shock dust destruction via thermal sputtering. We show that incorporating a sufficiently strong metallicity dependence into models of galaxy evolution removes the need for low stellar dust yields. The contribution of stellar sources to the overall dust budget may be significantly underestimated, and that of grain growth overestimated, by models assuming a constant destruction efficiency. △ Less

Submitted 13 October, 2021; originally announced October 2021.

Comments: 5 pages, 2 figures. MNRAS accepted

First deep images catalogue of extended IPHAS PNe

Authors: L. Sabin, M. A. Guerrero, G. Ramos-Larios, P. Boumis, A. A. Zijlstra, D. N. F. Awang Iskandar, M. J. Barlow, J. A. Toalá, Q. A. Parker, R. M. L. Corradi, R. A. H. Morris

Abstract: We present the first instalment of a deep imaging catalogue containing 58 True, Likely and Possible extended PNe detected with the Isaac Newton Telescope Photometric H$α$ Survey (IPHAS). The three narrow-band filters in the emission lines of H$α$, [N II] $λ$6584 Å and [O III] $λ$5007 Å used for this purpose allowed us to improve our description of the morphology and dimensions of the nebulae. In s… ▽ More We present the first instalment of a deep imaging catalogue containing 58 True, Likely and Possible extended PNe detected with the Isaac Newton Telescope Photometric H$α$ Survey (IPHAS). The three narrow-band filters in the emission lines of H$α$, [N II] $λ$6584 Å and [O III] $λ$5007 Å used for this purpose allowed us to improve our description of the morphology and dimensions of the nebulae. In some cases even the nature of the source has been reassessed. We were then able to unveil new macro-and micro-structures, which will without a doubt contribute to a more accurate analysis of these PNe. It has been also possible to perform a primary classification of the targets based on their ionization level. A Deep Learning classification tool has also been tested. We expect that all the PNe from the IPHAS catalogue of new extended planetary nebulae will ultimately be part of this deep H$α$, [N II] and [O III] imaging catalogue. △ Less

Submitted 31 August, 2021; originally announced August 2021.

Comments: 19 pages, 15 figures, 4 tables; accepted to MNRAS

High resolution H-alpha imaging of the Northern Galactic Plane, and the IGAPS images database

Authors: R. Greimel, J. E. Drew, M. Monguió, R. P. Ashley, G. Barentsen, J. Eislöffel, A. Mampaso, R. A. H. Morris, T. Naylor, C. Roe, L. Sabin, B. Stecklum, N. J. Wright, P. J. Groot, M. J. Irwin, M. J. Barlow, C. Fariña, A. Fernández-Martín, Q. A. Parker, S. Phillipps, S. Scaringi, A. A. Zijlstra

Abstract: The INT Galactic Plane Survey (IGAPS) is the merger of the optical photometric surveys, IPHAS and UVEX, based on data from the Isaac Newton Telescope (INT) obtained between 2003 and 2018. These capture the entire northern Galactic plane within the Galactic coordinate range, -5<b<+5 deg. and 30<l<215 deg. From the beginning, the incorporation of narrowband H-alpha imaging has been a unique and dist… ▽ More The INT Galactic Plane Survey (IGAPS) is the merger of the optical photometric surveys, IPHAS and UVEX, based on data from the Isaac Newton Telescope (INT) obtained between 2003 and 2018. These capture the entire northern Galactic plane within the Galactic coordinate range, -5<b<+5 deg. and 30<l<215 deg. From the beginning, the incorporation of narrowband H-alpha imaging has been a unique and distinctive feature of this effort. Alongside a focused discussion of the nature and application of the H-alpha data, we present the IGAPS world-accessible database of images for all 5 survey filters, i, r, g, U-RGO and narrowband H-alpha, observed on a pixel scale of 0.33 arcsec and at an effective (median) angular resolution of 1.1 to 1.3 arcsec. The background, noise, and sensitivity characteristics of the narrowband H-alpha filter images are outlined. Typical noise levels in this band correspond to a surface brightness at full one-arcsec resolution of around 2e-16 erg/cm2/s/arcsec2. Illustrative applications of the H-alpha data to planetary nebulae and Herbig-Haro objects are outlined and, as part of a discussion of mosaicking technique, we present a very large background-subtracted narrowband mosaic of the supernova remnant, Simeis 147. Finally we lay out a method that exploits the database via an automated selection of bright ionized diffuse interstellar emission targets for the coming generation of wide-field massive-multiplex spectrographs. Two examples of the diffuse H-alpha maps output from this selection process are presented and compared with previously published data. △ Less

Submitted 27 July, 2021; originally announced July 2021.

Comments: 21 pages, 14 main-text figures, 11 appendix figures. Images database and other supplementary items mentioned in the paper are available from http://www.igapsimages.org

Journal ref: A&A 655, A49 (2021)

The dust mass in Cassiopeia A from infrared and optical line flux differences

Authors: Maria Niculescu-Duvaz, Michael J. Barlow, Antonia Bevan, Danny Milisavljevic, Ilse De Looze

Abstract: The large quantities of dust that have been found in a number of high redshift galaxies have led to suggestions that core-collapse supernovae (CCSNe) are the main sources of their dust and have motivated the measurement of the dust masses formed by local CCSNe. For Cassiopeia~A, an oxygen-rich remnant of a Type~IIb CCSN, a dust mass of 0.6-1.1~M$_\odot$ has already been determined by two different… ▽ More The large quantities of dust that have been found in a number of high redshift galaxies have led to suggestions that core-collapse supernovae (CCSNe) are the main sources of their dust and have motivated the measurement of the dust masses formed by local CCSNe. For Cassiopeia~A, an oxygen-rich remnant of a Type~IIb CCSN, a dust mass of 0.6-1.1~M$_\odot$ has already been determined by two different methods, namely (a) from its far-infrared spectral energy distribution and (b) from analysis of the red-blue emission line asymmetries in its integrated optical spectrum. We present a third, independent, method for determining the mass of dust contained within Cas~A. This compares the relative fluxes measured in similar apertures from [O~{\sc iii}] far-infrared and visual-region emission lines, taking into account foreground dust extinction, in order to determine internal dust optical depths, from which corresponding dust masses can be obtained. Using this method we determine a dust mass within Cas~A of at least 0.99$^{+0.10}_{-0.09}$~M$_\odot$. △ Less

Submitted 1 June, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

Comments: Accepted by MNRAS, 17 pages, 12 figures. Author accepted manuscript. Accepted on 21/03/2021. Deposited on 22/03/2021

The efficiency of grain growth in the diffuse interstellar medium

Authors: F. D. Priestley, I. De Looze, M. J. Barlow

Abstract: 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 efficien… ▽ More 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. △ Less

Submitted 12 January, 2021; originally announced January 2021.

Comments: 9 pages, 11 figures. MNRAS accepted 11/01/21

H$α$ fluxes and extinction distances for planetary nebulae in the IPHAS survey of the Northern Galactic Plane

Authors: Thavisha E. Dharmawardena, M. J. Barlow, J. E. Drew, A. Seales, S. E. Sale, D. Jones, A. Mampaso, Q. A. Parker, L. Sabin, R. Wesson

Abstract: We report H$α$ filter photometry for 197 northern hemisphere planetary nebulae (PNe) obtained using imaging data from the IPHAS survey. H$α$+[N II] fluxes were measured for 46 confirmed or possible PNe discovered by the IPHAS survey and for 151 previously catalogued PNe that fell within the area of the northern Galactic Plane surveyed by IPHAS. After correcting for [N II] emission admitted by the… ▽ More We report H$α$ filter photometry for 197 northern hemisphere planetary nebulae (PNe) obtained using imaging data from the IPHAS survey. H$α$+[N II] fluxes were measured for 46 confirmed or possible PNe discovered by the IPHAS survey and for 151 previously catalogued PNe that fell within the area of the northern Galactic Plane surveyed by IPHAS. After correcting for [N II] emission admitted by the IPHAS H$α$ filter, the resulting H$α$ fluxes were combined with published radio free-free fluxes and H$β$ fluxes, in order to estimate mean optical extinctions to 143 PNe using ratios involving their integrated Balmer line fluxes and their extinction-free radio fluxes. Distances to the PNe were then estimated using three different 3D interstellar dust extinction mapping methods, including the IPHAS-based H-MEAD algorithm of Sale (2014). These methods were used to plot dust extinction versus distance relationships for the lines of sight to the PNe; the intercepts with the derived dust optical extinctions allowed distances to the PNe to be inferred. For 17 of the PNe in our sample reliable Gaia DR2 distances were available and these have been compared with the distances derived using three different extinction mapping algorithms as well as with distances from the nebular radius vs. H$α$ surface brightness relation of Frew et al. (2016). That relation and the H-MEAD extinction mapping algorithm yielded the closest agreement with the Gaia DR2 distances. △ Less

Submitted 4 December, 2020; originally announced December 2020.

Comments: Accepted by MNRAS, 13 pages, 3 figures. Author accepted manuscript. Accepted on 04/12/2020. Deposited on 04/12/2020

Revisiting the dust destruction efficiency of supernovae

Authors: F. D. Priestley, H. Chawner, M. Matsuura, I. De Looze, M. J. Barlow, H. L. Gomez

Abstract: Dust destruction by supernovae is one of the main processes removing dust from the interstellar medium (ISM). Estimates of the efficiency of this process, both theoretical and observational, typically assume a shock propagating into a homogeneous medium, whereas the ISM possesses significant substructure in reality. We self-consistently model the dust and gas properties of the shocked ISM in three… ▽ More Dust destruction by supernovae is one of the main processes removing dust from the interstellar medium (ISM). Estimates of the efficiency of this process, both theoretical and observational, typically assume a shock propagating into a homogeneous medium, whereas the ISM possesses significant substructure in reality. We self-consistently model the dust and gas properties of the shocked ISM in three supernova remnants (SNRs), using X-ray and infrared (IR) data combined with corresponding emission models. Collisional heating by gas with properties derived from X-ray observations produces dust temperatures too high to fit the far-IR fluxes from each SNR. An additional colder dust component is required, which has a minimum mass several orders of magnitude larger than that of the warm dust heated by the X-ray emitting gas. Dust-to-gas mass ratios indicate that the majority of the dust in the X-ray emitting material has been destroyed, while the fraction of surviving dust in the cold component is plausibly close to unity. As the cold component makes up virtually all the total dust mass, destruction timescales based on homogeneous models, which cannot account for multiple phases of shocked gas and dust, may be significantly overestimating actual dust destruction efficiencies, and subsequently underestimating grain lifetimes. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: 12 pages, 9 figures. MNRAS accepted 31/10/2020

Eta Carinae & the Homunculus: Far Infrared, Sub-millimeter Spectral Lines

Authors: T. R. Gull, P. W. Morris, J. H. Black, K. E. Nielsen, M. J. Barlow, P. Royer, B. M. Swinyard

Abstract: The evolved massive binary star Eta Carinae underwent eruptive mass loss events that formed the complex bi-polar Homunculus nebula harboring tens of solar masses of unusually nitrogen-rich gas and dust. Despite expectations for the presence of a significant molecular component to the gas, detections have been observationally challenged by limited access to the far-infrared and the intense thermal… ▽ More The evolved massive binary star Eta Carinae underwent eruptive mass loss events that formed the complex bi-polar Homunculus nebula harboring tens of solar masses of unusually nitrogen-rich gas and dust. Despite expectations for the presence of a significant molecular component to the gas, detections have been observationally challenged by limited access to the far-infrared and the intense thermal continuum. A spectral survey of the atomic and rotational molecular transitions was carried out with the Herschel Space Observatory, revealing a rich spectrum of broad emission lines originating in the ejecta. Velocity profiles of selected PACS lines correlate well with known substructures: H I in the central core; NH and weak [C II] within the Homunculus; and [N II] emissions in fast-moving structures external to the Homunculus. We have identified transitions from [O I], H I, and 18 separate light C- and O-bearing molecules including CO, CH, CH+, and OH, and a wide set of N-bearing molecules, NH, NH+, N2H+, NH2, NH3, HCN, HNC, CN, and N2H+. Half of these are new detections unprecedented for any early-type massive star environment. A very low ratio [12C/13C] LE 4 is estimated from five molecules and their isotopologues. We demonstrate that non-LTE effects due to the strong continuum are significant. Abundance patterns are consistent with line formation in regions of carbon and oxygen depletions with nitrogen enhancements, reflecting an evolved state of the erupting star with efficient transport of CNO-processed material to the outer layers. The results offer many opportunities for further observational and theoretical investigations of the molecular chemistry under extreme physical and chemical conditions around massive stars in their final stages of evolution. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: 37 pages, 16 figures, 5 tables, accepted MNRAS

A Galactic Dust Devil: far-infrared observations of the Tornado Supernova Remnant candidate

Authors: Hannah Chawner, Alex D. P. Howard, Haley L. Gomez, Mikako Matsuura, Felix Priestley, Mike J. Barlow, Ilse De Looze, Andreas Papageorgiou, Ken Marsh, Matt W. L. Smith, Alberto Noriega-Crespo, Jeonghee Rho, Loretta Dunne

Abstract: We present complicated dust structures within multiple regions of the candidate supernova remnant (SNR) the `Tornado' (G357.7-0.1) using observations with Spitzer and Herschel. We use Point Process Mapping, PPMAP, to investigate the distribution of dust in the Tornado at a resolution of 8", compared to the native telescope beams of 5-36". We find complex dust structures at multiple temperatures wi… ▽ More We present complicated dust structures within multiple regions of the candidate supernova remnant (SNR) the `Tornado' (G357.7-0.1) using observations with Spitzer and Herschel. We use Point Process Mapping, PPMAP, to investigate the distribution of dust in the Tornado at a resolution of 8", compared to the native telescope beams of 5-36". We find complex dust structures at multiple temperatures within both the head and the tail of the Tornado, ranging from 15 to 60K. Cool dust in the head forms a shell, with some overlap with the radio emission, which envelopes warm dust at the X-ray peak. Akin to the terrestrial sandy whirlwinds known as `Dust Devils', we find a large mass of dust contained within the Tornado. We derive a total dust mass for the Tornado head of 16.7 solar masses, assuming a dust absorption coefficient of kappa_300 =0.56m^2 kg^1, which can be explained by interstellar material swept up by a SNR expanding in a dense region. The X-ray, infra-red, and radio emission from the Tornado head indicate that this is a SNR. The origin of the tail is more unclear, although we propose that there is an X-ray binary embedded in the SNR, the outflow from which drives into the SNR shell. This interaction forms the helical tail structure in a similar manner to that of the SNR W50 and microquasar SS433. △ Less

Submitted 22 September, 2020; v1 submitted 17 September, 2020; originally announced September 2020.

Comments: 16 pages, 10 figures + 3 appendix figures. Accepted to be published in MNRAS

Constraining early-time dust formation in core-collapse supernovae

Authors: F. D. Priestley, A. Bevan, M. J. Barlow, I. De Looze

Abstract: 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… ▽ More 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. △ Less

Submitted 10 July, 2020; originally announced July 2020.

Comments: 14 pages, 13 figures. MNRAS accepted 10/07/20

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

Authors: A. M. Bevan, K. Krafton, R. Wesson, J. E. Andrews, E. Montiel, M. Niculescu-Duvaz, M. J. Barlow, I. De Looze, G. C. Clayton

Abstract: 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 evalu… ▽ More 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. △ Less

Submitted 3 April, 2020; originally announced April 2020.

Comments: Accepted by ApJ. Author accepted manuscript. Accepted on 26/03/2020. Deposited on 03/04/20. 21 pages, 10 figures

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

Authors: Florian Kirchschlager, M. J. Barlow, Franziska D. Schmidt

Abstract: 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 tra… ▽ More 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$. △ Less

Submitted 10 March, 2020; v1 submitted 6 March, 2020; originally announced March 2020.

Comments: Accepted by ApJ. Author accepted manuscript. Accepted on 06/03/2020. Deposited on 06/03/2020. 11 pages

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

Authors: F. D. Priestley, M. J. Barlow, I. De Looze, H. Chawner

Abstract: 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 s… ▽ More 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. △ Less

Submitted 4 December, 2019; originally announced December 2019.

Comments: 13 pages, 10 figures. MNRAS accepted 04/12/19

High angular resolution ALMA images of dust and molecules in the SN 1987A ejecta

Authors: Phil Cigan, Mikako Matsuura, Haley L. Gomez, Remy Indebetouw, Fran Abellán, Michael Gabler, Anita Richards, Dennis Alp, Tim Davis, Hans-Thomas Janka, Jason Spyromilio, M. J. Barlow, David Burrows, Eli Dwek, Claes Fransson, Bryan Gaensler, Josefin Larsson, P. Bouchet, Peter Lundqvist, J. M. Marcaide, C. -Y. Ng, Sangwook Park, Pat Roche, Jacco Th. van Loon, J. C. Wheeler , et al. (1 additional authors not shown)

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

Submitted 7 October, 2019; originally announced October 2019.

Comments: 32 pages, containing 19 figures and three appendices

Dust destruction by the reverse shock in the clumpy supernova remnant Cassiopeia A based on hydrodynamic simulations

Authors: Florian Kirchschlager, Franziska D. Schmidt, M. J. Barlow, Erica L. Fogerty, Antonia Bevan, Felix D. Priestley

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

Submitted 19 September, 2019; originally announced September 2019.

Comments: Conference proceeding

Dust survival rates in clumps passing through the Cas A reverse shock I: results for a range of clump densities

Authors: Florian Kirchschlager, Franziska D. Schmidt, M. J. Barlow, Erica L. Fogerty, Antonia Bevan, Felix D. Priestley

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

Submitted 28 August, 2019; originally announced August 2019.

Comments: Accepted by MNRAS. Author accepted manuscript. Accepted on 28/08/2019. Deposited on 28/08/19. 34 pages

Mid-infrared evolution of eta Car from 1968 to 2018

Authors: A. Mehner, W. -J. de Wit, D. Asmus, P. W. Morris, C. Agliozzo, M. J. Barlow, T. R. Gull, D. J. Hillier, G. Weigelt

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

Submitted 24 August, 2019; originally announced August 2019.

Comments: 5 pages, 4 Figures, accepted for publication in A&A

Journal ref: A&A 630, L6 (2019)

The dust content of the Crab Nebula

Authors: I. De Looze, M. J. Barlow, R. Bandiera, A. Bevan, M. F. Bietenholz, H. Chawner, H. L. Gomez, M. Matsuura, F. Priestley, R. Wesson

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

Submitted 5 June, 2019; originally announced June 2019.

Comments: 33 pages, 17 figures, Manuscript accepted for publication in MNRAS. Accepted on 31/05/2019. Deposited on 31/05/2019

The mass, location and heating of the dust in the Cassiopeia A supernova remnant

Authors: F. D. Priestley, M. J. Barlow, I. De Looze

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

Submitted 5 February, 2019; originally announced February 2019.

Comments: Accepted by MNRAS, 14 pages, 6 figures. Author accepted manuscript. Accepted on 05/02/2019. Deposited on 05/02/2019

An Imaging Spectroscopic Survey of the Planetary Nebula NGC 7009 with MUSE

Authors: J. R. Walsh, A. Monreal-Ibero, M. J. Barlow, T. Ueta, R. Wesson, A. A. Zijlstra, S. Kimeswenger, M. L. Leal-Ferreira, M. Otsuka

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

Submitted 9 October, 2018; originally announced October 2018.

Comments: Comments: 27 pages, 25 figures, 4 tables. Accepted for publication in A&A

Journal ref: A&A 620, A169 (2018)

SOFIA mid-infrared observations of Supernova 1987A in 2016 --- forward shocks and possible dust re-formation in the post-shocked region?

Authors: Mikako Matsuura, James M. De Buizer, Richard G. Arendt, Eli Dwek, M. J. Barlow, Antonia Bevan, Phil Cigan, Haley L. Gomez, Jeonghee Rho, Roger Wesson, Patrice Bouchet, John Danziger, Margaret Meixner

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

Submitted 8 October, 2018; originally announced October 2018.

Comments: 10 mages; MNRAS accented

A decade of ejecta dust formation in the Type IIn SN 2005ip

Authors: A. Bevan, R. Wesson, M. J. Barlow, I. De Looze, J. E. Andrews, G. C. Clayton, K. Krafton, M. Matsuura, D. Milisavljevic

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

Submitted 7 March, 2019; v1 submitted 24 September, 2018; originally announced September 2018.

Comments: Accepted by MNRAS, 17 pages, 11 figures. Author accepted manuscript. Accepted on 04/03/19. Deposited on 07/03/19

OH$^+$ emission from cometary knots in planetary nebulae

Authors: F. D. Priestley, M. J. Barlow

Abstract: 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 orde… ▽ More 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. △ Less

Submitted 27 April, 2018; originally announced April 2018.

Comments: Accepted by MNRAS, 11 pages, 15 figures. Author accepted manuscript. Accepted on 24/04/18. Deposited on 27/04/18

Modelling the ArH$^+$ emission from the Crab Nebula

Authors: F. D. Priestley, M. J. Barlow, S. Viti

Abstract: 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, $ζ_0$, is required to account for the lack of detected [C I] emission in published He… ▽ More 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, $ζ_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}$, $ζ= 10^7 ζ_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. △ Less

Submitted 5 September, 2017; originally announced September 2017.

Comments: Accepted by MNRAS. Author accepted manuscript. Accepted on 05/09/2017. Deposited on 05/09/17

$η$ Carinae's Dusty Homunculus Nebula from Near-Infrared to Submillimeter Wavelengths: Mass, Composition, and Evidence for Fading Opacity

Authors: Patrick W. Morris, Theodore R. Gull, D. John Hillier, M. J. Barlow, Pierre Royer, Krister Nielsen, John Black, Bruce Swinyard

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

Submitted 15 June, 2017; originally announced June 2017.

Comments: Accepted in ApJ

Journal ref: ApJ, 842, 79 (2017)

ALMA spectral survey of Supernova 1987A --- molecular inventory, chemistry, dynamics and explosive nucleosynthesis

Authors: M. Matsuura, R. Indebetouw, S. Woosley, V. Bujarrabal, F. J. Abellan, R. McCray, J. Kamenetzky, C. Fransson, M. J. Barlow, H. L. Gomez, P. Cigan, I De Looze, J. Spyromilio, L. Staveley-Smith, G. Zanardo, P. Roche, J. Larsson, S. Viti, J. Th. van Loon, J. C. Wheeler, M. Baes, R. Chevalier, P. Lundqvist, J. M. Marcaide, E. Dwek , et al. (4 additional authors not shown)

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

Submitted 7 April, 2017; originally announced April 2017.

Comments: Accepted by MNRAS, 16 pages

Atomic Gas in Debris Discs

Authors: Antonio S. Hales, M. J. Barlow, I. A. Crawford, S. Casassus

Abstract: 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 ava… ▽ More 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. △ Less

Submitted 16 December, 2016; originally announced December 2016.

Comments: Accepted 13/12/2016. Received 2/12/2016; Deposited on 22/11/2016. - 13 Pages, 9 Figures - MNRAS Advance Access published December 15, 2016

Dust masses for SN 1980K, SN1993J and Cassiopeia A from red-blue emission line asymmetries

Authors: Antonia Bevan, M. J. Barlow, D. Milisavljevic

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

Submitted 15 November, 2016; originally announced November 2016.

Comments: Author accepted manuscript. Accepted on 15/11/16. Deposited on 15/11/16

The dust mass in Cassiopeia A from a spatially resolved Herschel analysis

Authors: I. De Looze, M. J. Barlow, B. M. Swinyard, J. Rho, H. L. Gomez, M. Matsuura, R. Wesson

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

Submitted 8 January, 2017; v1 submitted 2 November, 2016; originally announced November 2016.

Comments: 37 pages, 24 figures, Manuscript published in MNRAS, including minor corrections. Accepted on 01/11/2016. Deposited on 08/01/2017

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

Authors: Mikako Matsuura, B. Sargent, Bruce Swinyard, Jeremy Yates, P. Royer, M. J. Barlow, Martha Boyer, L. Decin, Theo Khouri, Margaret Meixner, Jacco Th. van Loon, Paul M. Woods

Abstract: 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 i… ▽ More 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. △ Less

Submitted 4 August, 2016; originally announced August 2016.

Comments: 11 pages; accepted by MNRAS