Showing 1–50 of 147 results for author: Dwek, E

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

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

JWST NIRCam Observations of SN 1987A: Spitzer Comparison and Spectral Decomposition

Authors: Richard G. Arendt, Martha L. Boyer, Eli Dwek, Mikako Matsuura, Aravind P. Ravi, Armin Rest, Roger Chevalier, Phil Cigan, Ilse De Looze, Guido De Marchi, Claes Fransson, Christa Gall, R. D. Gehrz, Haley L. Gomez, Tuomas Kangas, Florian Kirchschlager, Robert P. Kirshner, Josefin Larsson, Peter Lundqvist, Dan Milisavljevic, Sangwook Park, Nathan Smith, Jason Spyromilio, Tea Temim, Lifan Wang , et al. (2 additional authors not shown)

Abstract: JWST NIRCam observations at 1.5-4.5 $μ$m have provided broad and narrow band imaging of the evolving remnant of SN 1987A with unparalleled sensitivity and spatial resolution. Comparing with previous marginally spatially resolved Spitzer IRAC observations from 2004-2019 confirms that the emission arises from the circumstellar equatorial ring (ER), and the current brightness at 3.6 and 4.5 $μ$m was… ▽ More JWST NIRCam observations at 1.5-4.5 $μ$m have provided broad and narrow band imaging of the evolving remnant of SN 1987A with unparalleled sensitivity and spatial resolution. Comparing with previous marginally spatially resolved Spitzer IRAC observations from 2004-2019 confirms that the emission arises from the circumstellar equatorial ring (ER), and the current brightness at 3.6 and 4.5 $μ$m was accurately predicted by extrapolation of the declining brightness tracked by IRAC. Despite the regular light curve, the NIRCam observations clearly reveal that much of this emission is from a newly developing outer portion of the ER. Spots in the outer ER tend to lie at position angles in between the well-known ER hotspots. We show that the bulk of the emission in the field can be represented by 5 standard spectral energy distributions (SEDs), each with a distinct origin and spatial distribution. This spectral decomposition provides a powerful technique for distinguishing overlapping emission from the circumstellar medium (CSM) and the supernova (SN) ejecta, excited by the forward and reverse shocks respectively. △ Less

Submitted 22 September, 2023; originally announced September 2023.

Comments: Accepted for publication in ApJ. 16 pages, 12 figures. 2 animations not included here

JWST Discovery of Dust Reservoirs in Nearby Type IIP Supernovae 2004et and 2017eaw

Authors: Melissa Shahbandeh, Arkaprabha Sarangi, Tea Temim, Tamas Szalai, Ori D. Fox, Samaporn Tinyanont, Eli Dwek, Luc Dessart, Alexei V. Filippenko, Thomas G. Brink, Ryan J. Foley, Jacob Jencson, Justin Pierel, Szanna Zsiros, Armin Rest, WeiKang Zheng, Jennifer Andrews, Geoffrey C. Clayton, Kishalay De, Michael Engesser, Suvi Gezari, Sebastian Gomez, Shireen Gonzaga, Joel Johansson, Mansi Kasliwal , et al. (14 additional authors not shown)

Abstract: Supernova (SN) explosions have been sought for decades as a possible source of dust in the Universe, providing the seeds of galaxies, stars, and planetary systems. SN 1987A offers one of the most promising examples of significant SN dust formation, but until the James Webb Space Telescope (JWST), instruments have traditionally lacked the sensitivity at both late times (>1 yr post-explosion) and lo… ▽ More Supernova (SN) explosions have been sought for decades as a possible source of dust in the Universe, providing the seeds of galaxies, stars, and planetary systems. SN 1987A offers one of the most promising examples of significant SN dust formation, but until the James Webb Space Telescope (JWST), instruments have traditionally lacked the sensitivity at both late times (>1 yr post-explosion) and longer wavelengths (i.e., >10 um) to detect analogous dust reservoirs. Here we present JWST/MIRI observations of two historic Type IIP SNe, 2004et and SN 2017eaw, at nearly 18 and 5 yr post-explosion, respectively. We fit the spectral energy distributions as functions of dust mass and temperature, from which we are able to constrain the dust geometry, origin, and heating mechanism. We place a 90% confidence lower limit on the dust masses for SNe 2004et and 2017eaw of >0.014 and >4e-4 M_sun, respectively. More dust may exist at even colder temperatures or may be obscured by high optical depths. We conclude dust formation in the ejecta to be the most plausible and consistent scenario. The observed dust is radiatively heated to ~100-150 K by ongoing shock interaction with the circumstellar medium. Regardless of the best fit or heating mechanism adopted, the inferred dust mass for SN 2004et is the second highest (next to SN 1987A) inferred dust mass in extragalactic SNe thus far, promoting the prospect of SNe as potential significant sources of dust in the Universe. △ Less

Submitted 25 January, 2023; originally announced January 2023.

Comments: 12 pages, 7 figures, submitting to MNRAS

Dust contribution to the panchromatic galaxy emission

Authors: Dian P. Triani, Darren J. Croton, Manodeep Sinha, Edward N. Taylor, Camilla Pacifici, Eli Dwek

Abstract: We have developed a pipeline called \mentari to generate the far-ultraviolet to far-infrared spectral energy distribution (SED) of galaxies from the \dustysage semi-analytic galaxy formation model (SAM). \dustysage incorporates dust-related processes directly on top of the basic ingredients of galaxy formation like gas infall, cooling, star formation, feedback, and mergers. We derive a physically… ▽ More We have developed a pipeline called \mentari to generate the far-ultraviolet to far-infrared spectral energy distribution (SED) of galaxies from the \dustysage semi-analytic galaxy formation model (SAM). \dustysage incorporates dust-related processes directly on top of the basic ingredients of galaxy formation like gas infall, cooling, star formation, feedback, and mergers. We derive a physically motivated attenuation model from the computed dust properties in \dustysage, so each galaxy has a self-consistent set of attenuation parameters based on the complicated dust physics that occurred across the galaxy's assembly history. Then, we explore several dust emission templates to produce infrared spectra. Our results show that a physically-motivated attenuation model is better for obtaining a consistent multi-wavelength description of galaxy formation and evolution, compared to using a constant attenuation. We compare our predictions with a compilation of observations and find that the fiducial model is in reasonable agreement with: (i) the observed $z=0$ luminosity functions from the far-ultraviolet to far-infrared simultaneously, and hence (ii) the local cosmic SED in the same range, (iii) the rest-frame K-band luminosity function across $0 < z < 3$, and (iv) the rest-frame far-ultraviolet luminosity function across $0 < z < 1$. Our model underproduces the far-ultraviolet emission at $z=2$ and $z=3$, which can be improved by altering the AGN feedback and dust processes in \dustysage. However, this combination thus worses the agreement at $z=0$, which suggests that more detailed treatment of such processes is required. △ Less

Submitted 12 December, 2022; originally announced December 2022.

Comments: 20 pages, 14 figures, accepted for publication in MNRAS

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

SKYSURF: Constraints on Zodiacal Light and Extragalactic Background Light through Panchromatic HST All-Sky Surface-Brightness Measurements: II. First Limits on Diffuse Light at 1.25, 1.4, and 1.6 microns

Authors: Timothy Carleton, Rogier A. Windhorst, Rosalia O'Brien, Seth H. Cohen, Delondrae Carter, Rolf Jansen, Scott Tompkins, Richard G. Arendt, Sarah Caddy, Norman Grogin, Scott J. Kenyon, Anton Koekemoer, John MacKenty, Stefano Casertano, Luke J. M. Davies, Simon P. Driver, Eli Dwek, Alexander Kashlinsky, Nathan Miles, Rushabh Pawnikar, Nor Pirzkal, Aaron Robotham, Russell Ryan, Haley Abate, Hanga Andras-Letanovszky , et al. (6 additional authors not shown)

Abstract: We present the first results from the HST Archival Legacy project "SKYSURF." As described in Windhorst et al. 2022, SKYSURF utilizes the large HST archive to study the diffuse UV, optical, and near-IR backgrounds and foregrounds in detail. Here we utilize SKYSURF's first sky-surface brightness measurements to constrain the level of near-IR diffuse Extragalactic Background Light (EBL). Our sky-surf… ▽ More We present the first results from the HST Archival Legacy project "SKYSURF." As described in Windhorst et al. 2022, SKYSURF utilizes the large HST archive to study the diffuse UV, optical, and near-IR backgrounds and foregrounds in detail. Here we utilize SKYSURF's first sky-surface brightness measurements to constrain the level of near-IR diffuse Extragalactic Background Light (EBL). Our sky-surface brightness measurements have been verified to an accuracy of better than 1%, which when combined with systematic errors associated with HST, results in sky brightness uncertainties of $\sim$2-4% $\simeq$ 0.005 MJy/sr in each image. We put limits on the amount of diffuse EBL in three near-IR filters (F125W, F140W, and F160W) by comparing our preliminary sky measurements of $> 30,000$ images to Zodiacal light models, carefully selecting the darkest images to avoid contamination from stray light. In addition, we investigate the impact that instrumental thermal emission has on our measurements, finding that it has a limited impact on F125W and F140W measurements, whereas uncertainties in the exact thermal state of HST results in significant uncertainties in the level of astrophysical diffuse light in F160W images. When compared to the Kelsall et al. (1998) Zodiacal model, an isotropic diffuse background of $30$ nW m$^{-2}$ sr$^{-1}$ remains, whereas using the Wright (1998) Zodiacal model results in no discernible diffuse background. Based primarily on uncertainties in the foreground model subtraction, we present limits on the amount of diffuse EBL of 29 nW m$^{-2}$ sr$^{-1}$, 40 nW m$^{-2}$ sr$^{-1}$, and 29 nW m$^{-2}$ sr$^{-1}$ for F125W, F140W, and F160W respectively. While this light is generally isotropic, our modeling at this point does not distinguish between a cosmological origin or a Solar System origin (such as a dim, diffuse, spherical cloud of cometary dust). △ Less

Submitted 23 August, 2022; v1 submitted 12 May, 2022; originally announced May 2022.

Comments: Accepted with Windhorst et al. 2022 to AJ. Main figures are Fig. 10 and 11. Comments welcome!

SKYSURF: Constraints on Zodiacal Light and Extragalactic Background Light through Panchromatic HST All-Sky Surface-Brightness Measurements: I. Survey Overview and Methods

Authors: Rogier A. Windhorst, Timothy Carleton, Rosalia O'Brien, Seth H. Cohen, Delondrae Carter, Rolf Jansen, Scott Tompkins, Richard G. Arendt, Sarah Caddy, Norman Grogin, Anton Koekemoer, John MacKenty, Stefano Casertano, Luke J. M. Davies, Simon P. Driver, Eli Dwek, Alexander Kashlinsky, Scott J. Kenyon, Nathan Miles, Nor Pirzkal, Aaron Robotham, Russell Ryan, Haley Abate, Hanga Andras-Letanovszky, Jessica Berkheimer , et al. (13 additional authors not shown)

Abstract: We give an overview and describe the rationale, methods, and testing of the Hubble Space Telescope (HST) Archival Legacy project "SKYSURF." SKYSURF uses HST's unique capability as an absolute photometer to measure the ~0.2-1.7 $μ$m sky surface brightness (SB) from 249,861 WFPC2, ACS, and WFC3 exposures in ~1400 independent HST fields. SKYSURF's panchromatic dataset is designed to constrain the dis… ▽ More We give an overview and describe the rationale, methods, and testing of the Hubble Space Telescope (HST) Archival Legacy project "SKYSURF." SKYSURF uses HST's unique capability as an absolute photometer to measure the ~0.2-1.7 $μ$m sky surface brightness (SB) from 249,861 WFPC2, ACS, and WFC3 exposures in ~1400 independent HST fields. SKYSURF's panchromatic dataset is designed to constrain the discrete and diffuse UV to near-IR sky components: Zodiacal Light (ZL; inner Solar System), Kuiper Belt Objects (KBOs; outer Solar System), Diffuse Galactic Light (DGL), and the discrete plus diffuse Extragalactic Background Light (EBL). We outline SKYSURF's methods to: (1) measure sky-SB levels between its detected objects; (2) measure the integrated discrete EBL, most of which comes from AB$\simeq$17-22 mag galaxies; and (3) estimate how much diffuse light may exist in addition to the extrapolated discrete galaxy counts. Simulations of HST WFC3/IR images with known sky-values and gradients, realistic cosmic ray (CR) distributions, and star plus galaxy counts were processed with nine different algorithms to measure the "Lowest Estimated Sky-SB" (LES) in each image between the discrete objects. The best algorithms recover the inserted LES values within 0.2% when there are no image gradients, and within 0.2-0.4% when there are 5-10% gradients. SKYSURF requires non-standard re-processing of these HST images that includes restoring the lowest sky-level from each visit into each drizzled image. We provide a proof of concept of our methods from the WFC3/IR F125W images, where any residual diffuse light that HST sees in excess of the Kelsall et al. (1998) Zodiacal model prediction does not depend on the total object flux that each image contains. This enables us to present our first SKYSURF results on diffuse light in Carleton et al. (2022). △ Less

Submitted 25 August, 2022; v1 submitted 12 May, 2022; originally announced May 2022.

Comments: Accepted to AJ; see accompanying paper Carleton et al. 2022: arXiv:2205.06347. Comments welcome!

Spitzer's Last Look at Extragalactic Explosions: Long-Term Evolution of Interacting Supernovae

Authors: Tamás Szalai, Ori D. Fox, Richard G. Arendt, Eli Dwek, Jennifer E. Andrews, Geoffrey C. Clayton, Alexei V. Filippenko, Joel Johansson, Patrick L. Kelly, Kelsie Krafton, A. P. Marston, Jon C. Mauerhan, Schuyler D. Van Dyk

Abstract: Here we present new - and, nevertheless, last - mid-infrared (mid-IR) data for supernovae (SNe) based on measurements with the Spitzer Space Telescope. Comparing our recent 3.6 and 4.5 $μ$m photometry with previously published mid-IR and further multiwavelength datasets, we were able to draw some conclusions about the origin and heating mechanism of the dust in these SNe or in their environments,… ▽ More Here we present new - and, nevertheless, last - mid-infrared (mid-IR) data for supernovae (SNe) based on measurements with the Spitzer Space Telescope. Comparing our recent 3.6 and 4.5 $μ$m photometry with previously published mid-IR and further multiwavelength datasets, we were able to draw some conclusions about the origin and heating mechanism of the dust in these SNe or in their environments, as well as on possible connection with circumstellar matter (CSM) originating from pre-explosion mass-loss events in the progenitor stars. We also present new results regarding both certain SN classes and single objects. We highlight the mid-IR homogeneity of SNe Ia-CSM, which may be a hint of their common progenitor type and of their basically uniform circumstellar environments. Regarding single objects, it is worth highlighting the late-time interacting Type Ib SNe 2003gk and 2004dk, for which we present the first-ever mid-IR data, which seem to be consistent with clues of ongoing CSM interaction detected in other wavelength ranges. Our current study suggests that long-term mid-IR follow-up observations play a key role in a better understanding of both pre- and post-explosion processes in SNe and their environments. While Spitzer is not available any more, the expected unique data from the James Webb Space Telescope, as well as long-term near-IR follow-up observations of dusty SNe, can bring us closer to the hidden details of this topic. △ Less

Submitted 23 June, 2021; originally announced June 2021.

Comments: 25 pages, 10 figures, 4 tables; accepted for publication in ApJ

The infrared echo of SN2010jl and its implications for shock breakout characteristics

Authors: Eli Dwek, Arkaprabha Sarangi, Richard G. Arendt, Timothy Kallman, Demos Kazanas, Ori D. Fox

Abstract: SN 2010jl is a Type IIn core collapse supernova whose radiative output is powered by the interaction of the SN shock wave with its surrounding dense circumstellar medium (CSM). After day ~60, its light curve developed a NIR excess emission from dust. This excess could be a thermal IR echo from pre-existing CSM dust, or emission from newly-formed dust either in the cooling postshock region of the C… ▽ More SN 2010jl is a Type IIn core collapse supernova whose radiative output is powered by the interaction of the SN shock wave with its surrounding dense circumstellar medium (CSM). After day ~60, its light curve developed a NIR excess emission from dust. This excess could be a thermal IR echo from pre-existing CSM dust, or emission from newly-formed dust either in the cooling postshock region of the CSM, or in the cooling SN ejecta. Recent analysis has shown that dust formation in the CSM can commence only after day ~380, and has also ruled out newly-formed ejecta dust as the source of the NIR emission. The early (< 380 d) NIR emission can therefore only be attributed to an IR echo. The H-K color temperature of the echo is about 1250 K. The best fitting model requires the presence of about 1.6e-4 Msun of amorphous carbon dust at a distance of 2.2e16 cm from the explosion. The CSM-powered luminosity is preceded by an intense burst of hard radiation generated by the breakout of the SN shock through the stellar surface. The peak burst luminosity seen by the CSM dust is significantly reduced by Thomson scattering in the CSM, but still has the potential of evaporating the dust needed to produce the echo. We show that the survival of the echo-producing dust provides important constraints on the intensity, effective temperature, and duration of the burst. △ Less

Submitted 8 June, 2021; originally announced June 2021.

Comments: Accepted for publication in the ApJ, 20 pages, 1 table, and 17 figures

Exploring the relation between dust mass and galaxy properties using Dusty SAGE

Authors: Dian P. Triani, Manodeep Sinha, Darren J. Croton, Eli Dwek, Camilla Pacifici

Abstract: We explore the relation between dust and several fundamental properties of simulated galaxies using the Dusty SAGE semi-analytic model. In addition to tracing the standard galaxy properties, Dusty SAGE also tracks cold dust mass in the interstellar medium (ISM), hot dust mass in the halo and dust mass ejected by feedback activity. Based on their ISM dust content, we divide our galaxies into two ca… ▽ More We explore the relation between dust and several fundamental properties of simulated galaxies using the Dusty SAGE semi-analytic model. In addition to tracing the standard galaxy properties, Dusty SAGE also tracks cold dust mass in the interstellar medium (ISM), hot dust mass in the halo and dust mass ejected by feedback activity. Based on their ISM dust content, we divide our galaxies into two categories: ISM dust-poor and ISM dust-rich. We split the ISM dust-poor group into two subgroups: halo dust-rich and dust-poor (the latter contains galaxies that lack dust in both the ISM and halo). Halo dust-rich galaxies have high outflow rates of heated gas and dust and are more massive. We divide ISM dust-rich galaxies based on their specific star formation rate (sSFR) into star-forming and quenched subgroups. At redshift z=0, we find that ISM dust-rich galaxies have a relatively high sSFR, low bulge-to-total (BTT) mass ratio, and high gas metallicity. The high sSFR of ISM dust-rich galaxies allows them to produce dust in the stellar ejecta. Their metal-rich ISM enables dust growth via grain accretion. The opposite is seen in the ISM dust-poor group. Furthermore, ISM dust-rich galaxies are typically late-types, while ISM dust-poor galaxies resemble the early-type population, and we show how their ISM content evolves from being dust-rich to dust-poor. Finally, we investigate dust production from z=3 to z=0 and find that all groups evolve similarly, except for the quenched ISM dust-rich group. △ Less

Submitted 25 February, 2021; v1 submitted 24 February, 2021; originally announced February 2021.

Comments: 12 pages, 8 figure, accepted for publication in MNRAS

The Dynamics, Destruction, and Survival of Supernova-Formed Dust Grains

Authors: Jonathan D. Slavin, Eli Dwek, Mordecai-Mark Mac Low, Alex S. Hill

Abstract: Observations have demonstrated that supernovae efficiently produce dust. This is consistent with the hypothesis that supernovae and asymptotic giant branch stars are the primary producers of dust in the Universe. However, there has been a longstanding question of how much of the dust detected in the interiors of young supernova remnants can escape into the interstellar medium. We present new hydro… ▽ More Observations have demonstrated that supernovae efficiently produce dust. This is consistent with the hypothesis that supernovae and asymptotic giant branch stars are the primary producers of dust in the Universe. However, there has been a longstanding question of how much of the dust detected in the interiors of young supernova remnants can escape into the interstellar medium. We present new hydrodynamical calculations of the evolution of dust grains that were formed in dense ejecta clumps within a Cas A-like remnant. We follow the dynamics of the grains as they decouple from the gas after their clump is hit by the reverse shock. They are subsequently subject to destruction by thermal and kinetic sputtering as they traverse the remnant. Grains that are large enough ($\sim 0.25\,μ$m for silicates and $\sim 0.1\,μ$m for carbonaceous grains) escape into the interstellar medium while smaller grains get trapped and destroyed. However, grains that reach the interstellar medium still have high velocities, and are subject to further destruction as they are slowed down. We find that for initial grain size distributions that include large ($\sim 0.25 - 0.5\,μ$m) grains, 10--20\% of silicate grains can survive, while 30--50\% of carbonaceous grains survive even when the initial size distribution cuts off at smaller ($0.25\,μ$m) sizes. For a 19 M$_{\odot}$ star similar to the progenitor of Cas A, up to 0.1 M$_{\odot}$ of dust can survive if the dust grains formed are large. Thus we show that supernovae under the right conditions can be significant sources of interstellar dust. △ Less

Submitted 3 September, 2020; originally announced September 2020.

Comments: 22 pages, 8 figures, accepted for publication in ApJ

Dust masses, compositions, and luminosities in the nuclear disks and the diffuse circumnuclear medium of Arp220

Authors: Eli Dwek, Richard G. Arendt

Abstract: We present an analysis of the 4-2600 $μ$m spectral energy distributions (SEDs) of the west and east nuclei and the diffuse infrared (IR) region of the merger-driven starburst Arp 220. We examine several possible source morphologies and dust temperature distributions using a mixture of silicate and carbonaceous grains. From fits to the SEDs we derive dust masses, temperatures, luminosities, and dus… ▽ More We present an analysis of the 4-2600 $μ$m spectral energy distributions (SEDs) of the west and east nuclei and the diffuse infrared (IR) region of the merger-driven starburst Arp 220. We examine several possible source morphologies and dust temperature distributions using a mixture of silicate and carbonaceous grains. From fits to the SEDs we derive dust masses, temperatures, luminosities, and dust inferred gas masses. We show that the west and east nuclei are powered by central sources deeply enshrouded behind $\sim 10^{26}$ cm$^{-2}$ column densities of hydrogen with an exponential density distribution, and that the west and east nuclei are optically thick out to wavelengths of $\sim 1900$ and $\sim 770$ $μ$m, respectively. The nature of the central sources cannot be determined from our analysis. We derive star formation rates or black hole masses needed to power the IR emission, and show that the [C II] 158$μ$m line cannot be used as a tracer of the star formation rate in heavily obscured systems. Dust inferred gas masses are larger than those inferred from CO observations, suggesting either larger dust-to-H mass ratios, or the presence of hidden atomic H that cannot be inferred from CO observations. The luminosities per unit mass in the nuclei are $\sim 450$, in solar units, smaller that the Eddington limit of $\sim 1000 - 3000$ for an optically thick star forming region, suggesting that the observed gas outflows are primarily driven by stellar winds and supernova shock waves instead of radiation pressure on the dust. △ Less

Submitted 10 August, 2020; originally announced August 2020.

Comments: Accepted for publication in The Astrophysical Journal

The origin of dust in galaxies across cosmic time

Authors: Dian P. Triani, Manodeep Sinha, Darren J. Croton, Camilla Pacifici, Eli Dwek

Abstract: We study the dust evolution in galaxies by implementing a detailed dust prescription in the SAGE semi-analytical model for galaxy formation. The new model, called Dusty SAGE, follows the condensation of dust in the ejecta of type II supernovae and asymptotic giant branch (AGB) stars, grain growth in the dense molecular clouds, destruction by supernovae shocks, and the removal of dust from the ISM… ▽ More We study the dust evolution in galaxies by implementing a detailed dust prescription in the SAGE semi-analytical model for galaxy formation. The new model, called Dusty SAGE, follows the condensation of dust in the ejecta of type II supernovae and asymptotic giant branch (AGB) stars, grain growth in the dense molecular clouds, destruction by supernovae shocks, and the removal of dust from the ISM by star formation, reheating, inflows and outflows. Our model successfully reproduces the observed dust mass function at redshift z = 0 and the observed scaling relations for dust across a wide range of redshifts. We find that the dust mass content in the present Universe is mainly produced via grain growth in the interstellar medium (ISM). By contrast, in the early Universe, the primary production mechanism for dust is the condensation in stellar ejecta. The shift of the significant production channel for dust characterises the scaling relations of dust-to-gas (DTG) and dust-to-metal (DTM) ratios. In galaxies where the grain growth dominates, we find positive correlations for DTG and DTM ratios with both metallicity and stellar mass. On the other hand, in galaxies where dust is produced primarily via condensation, we find negative or no correlation for DTM and DTG ratios with either metallicity or stellar mass. In agreement with observation showing that the circumgalactic medium (CGM) contains more dust than the ISM, our model also shows the same trend for z < 4. Our semi-analytic model is publicly available at https: //github.com/dptriani/dusty-sage. △ Less

Submitted 12 February, 2020; originally announced February 2020.

Comments: 19 pages, 14 figures, accepted for publication in MNRAS

Journal ref: 2020MNRAS.tmp..415T

Final Spitzer IRAC Observations of the Rise and Fall of SN 1987A

Authors: Richard G. Arendt, Eli Dwek, Patrice Bouchet, I. John Danziger, Robert D. Gehrz, Sangwook Park, Charles E. Woodward

Abstract: Spitzer's final Infrared Array Camera (IRAC) observations of SN 1987A show the 3.6 and 4.5 $μ$m emission from the equatorial ring (ER) continues a period of steady decline. Deconvolution of the images reveals that the emission is dominated by the ring, not the ejecta, and is brightest on the west side. Decomposition of the marginally resolved emission also confirms this, and shows that the west si… ▽ More Spitzer's final Infrared Array Camera (IRAC) observations of SN 1987A show the 3.6 and 4.5 $μ$m emission from the equatorial ring (ER) continues a period of steady decline. Deconvolution of the images reveals that the emission is dominated by the ring, not the ejecta, and is brightest on the west side. Decomposition of the marginally resolved emission also confirms this, and shows that the west side of the ER has been brightening relative to the other portions of the ER. The infrared (IR) morphological changes resemble those seen in both the soft X-ray emission and the optical emission. The integrated ER light curves at 3.6 and 4.5 $μ$m are more similar to the optical light curves than the soft X-ray light curve, though differences would be expected if dust is responsible for this emission and its destruction is rapid. Future observations with the James Webb Space Telescope will continue to monitor the ER evolution, and will reveal the true spectrum and nature of the material responsible for the broadband emission at 3.6 and 4.5 $μ$m. The present observations also serendipitously reveal a nearby variable source, subsequently identified as a Be star, that has gone through a multi-year outburst during the course of these observations. △ Less

Submitted 27 December, 2019; originally announced December 2019.

Comments: 12 pages, 9 figures. Accepted for publication in the ApJ

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 formation in AGN winds

Authors: Arkaprabha Sarangi, Eli Dwek, Demos Kazanas

Abstract: Infrared observations of active galactic nucleus (AGN) reveal emission from the putative dusty circumnuclear 'torus' invoked by AGN unification, that is heated up by radiation from the central accreting black hole (BH). The strong 9.7 and 18 micron silicate features observed in the AGN spectra both in emission and absorption, further indicate the presence of such dusty environments. We present det… ▽ More Infrared observations of active galactic nucleus (AGN) reveal emission from the putative dusty circumnuclear 'torus' invoked by AGN unification, that is heated up by radiation from the central accreting black hole (BH). The strong 9.7 and 18 micron silicate features observed in the AGN spectra both in emission and absorption, further indicate the presence of such dusty environments. We present detailed calculations of the chemistry of silicate dust formation in AGN accretion disk winds. The winds considered herein are magnetohydrodynamic (MHD) winds driven off the entire accretion disk domain that extends from the BH vicinity to the radius of BH influence, of order of 1 to 100 pc depending on the mass of the resident BH. Our results indicate that these winds provide conditions conducive to the formation of significant amounts of dust, especially for objects accreting close to their Eddington limit, making AGN a significant source of dust in the universe, especially for luminous quasars. Our models justify the importance of a r to the power -1 density law in the winds for efficient formation and survival of dust grains. The dust production rate scales linearly with the mass of the central BH and varies as a power law of index between 2 to 2.5 with the dimensionless mass accretion rate. The resultant distribution of the dense dusty gas resembles a toroidal shape, with high column density and optical depths along the equatorial viewing angles, in agreement with the AGN unification picture. △ Less

Submitted 23 September, 2019; originally announced September 2019.

Comments: 22 pages, 17 figures (some have subfigures), Accepted for publication in The Astrophysical Journal

2 mm GISMO Observations of the Galactic Center. II. A Nonthermal Filament in the Radio Arc and Compact Sources

Authors: Johannes Staguhn, Richard G. Arendt, Eli Dwek, Mark R. Morris, Farhad Yusef-Zadeh, Dominic J. Benford, Attila Kovács, Junellie Gonzalez-Quiles

Abstract: We have used the Goddard IRAM 2-Millimeter Observer (GISMO) with the 30 m IRAM telescope to carry out a 2 mm survey of the Galaxy's central molecular zone (CMZ). These observations detect thermal emission from cold ISM dust, thermal free-free emission from ionized gas, and nonthermal synchrotron emission from relatively flat-spectrum sources. Archival data sets spanning $3.6 μ$m to 90 cm are used… ▽ More We have used the Goddard IRAM 2-Millimeter Observer (GISMO) with the 30 m IRAM telescope to carry out a 2 mm survey of the Galaxy's central molecular zone (CMZ). These observations detect thermal emission from cold ISM dust, thermal free-free emission from ionized gas, and nonthermal synchrotron emission from relatively flat-spectrum sources. Archival data sets spanning $3.6 μ$m to 90 cm are used to distinguish different emission mechanisms. After the thermal emission of dust is modeled and subtracted, the remaining 2 mm emission is dominated by free-free emission, with the exception of the brightest nonthermal filament (NTF) that runs though the middle of the bundle of filaments known as the Radio Arc. This is the shortest wavelength at which any NTF has been detected. The GISMO observations clearly trace this NTF over a length of ~0.2$^\circ$, with a mean 2 mm spectral index which is steeper than at longer wavelengths. The 2 mm to 6 cm (or 20 cm) spectral index steepens from $α\approx -0.2$ to $-0.7$ as a function distance from the Sickle H II region, suggesting that this region is directly related to the NTF. A number of unresolved (at $21''$) 2 mm sources are found nearby. One appears to be thermal dust emission from a molecular cloud that is associated with an enigmatic radio point source whose connection to the Radio Arc is still debated. The morphology and colors at shorter IR wavelengths indicate other 2 mm unresolved sources are likely to be compact H II regions. △ Less

Submitted 17 September, 2019; originally announced September 2019.

Comments: 17 pages, 11 figures. Accepted for publication in the ApJ

2 mm GISMO Observations of the Galactic Center. I. Dust Emission

Authors: Richard G. Arendt, Johannes Staguhn, Eli Dwek, Mark R. Morris, Farhad Yusef-Zadeh, Dominic J. Benford, Attila Kovács, Junellie Gonzalez-Quiles

Abstract: The Central Molecular Zone (CMZ), covering the inner ~1$^\circ$ of the Galactic plane has been mapped at 2 mm using the GISMO bolometric camera on the 30 m IRAM telescope. The $21''$ resolution maps show abundant emission from cold molecular clouds, from star forming regions, and from one of the Galactic center nonthermal filaments. In this work we use the Herschel Hi-GAL data to model the dust em… ▽ More The Central Molecular Zone (CMZ), covering the inner ~1$^\circ$ of the Galactic plane has been mapped at 2 mm using the GISMO bolometric camera on the 30 m IRAM telescope. The $21''$ resolution maps show abundant emission from cold molecular clouds, from star forming regions, and from one of the Galactic center nonthermal filaments. In this work we use the Herschel Hi-GAL data to model the dust emission across the Galactic center. We find that a single-temperature fit can describe the 160 -- 500 $μ$m emission for most lines of sight, if the long-wavelength dust emissivity scales as $λ^{-β}$ with $β\approx 2.25$. This dust model is extrapolated to predict the 2 mm dust emission. Subtraction of the model from the GISMO data provides a clearer look at the 2 mm emission of star-forming regions and the brightest nonthermal filament. △ Less

Submitted 17 September, 2019; originally announced September 2019.

Comments: 18 pages, 15 figures. Accepted for publication in the ApJ

Determination of the Cosmic Infrared Background from COBE/FIRAS and Planck HFI Observations

Authors: N. Odegard, J. L. Weiland, D. J. Fixsen, D. T. Chuss, E. Dwek, A. Kogut, E. R. Switzer

Abstract: New determinations are presented of the cosmic infrared background monopole brightness in the Planck HFI bands from 100 GHz to 857 GHz. Planck was not designed to measure the monopole component of sky brightness, so cross-correlation of the 2015 HFI maps with COBE/FIRAS data is used to recalibrate the zero level of the HFI maps. For the HFI 545 and 857 GHz maps, the brightness scale is also recali… ▽ More New determinations are presented of the cosmic infrared background monopole brightness in the Planck HFI bands from 100 GHz to 857 GHz. Planck was not designed to measure the monopole component of sky brightness, so cross-correlation of the 2015 HFI maps with COBE/FIRAS data is used to recalibrate the zero level of the HFI maps. For the HFI 545 and 857 GHz maps, the brightness scale is also recalibrated. Correlation of the recalibrated HFI maps with a linear combination of Galactic H I and H alpha data is used to separate the Galactic foreground emission and determine the cosmic infrared background brightness in each of the HFI bands. We obtain CIB values of 0.007 +- 0.014, 0.010 +- 0.019, 0.060 +- 0.023, 0.149 +- 0.017, 0.371 +- 0.018, and 0.576 +- 0.034 MJy/sr at 100, 143, 217, 353, 545, and 857 GHz, respectively. The estimated uncertainties for the 353 to 857 GHz bands are about 3 to 6 times smaller than those of previous direct CIB determinations at these frequencies. Our results are compared with integrated source brightness results from selected recent submillimeter and millimeter wavelength imaging surveys. △ Less

Submitted 23 May, 2019; v1 submitted 25 April, 2019; originally announced April 2019.

Comments: 29 pages, 18 figures, accepted for publication in ApJ; minor changes in v2 to match published version

Journal ref: Ap.J. 877, 40 (2019)

Taking Census of Massive, Star-Forming Galaxies formed <1 Gyr After the Big Bang

Authors: Caitlin M. Casey, Peter Capak, Johannes Staguhn, Lee Armus, Andrew Blain, Matthieu Bethermin, Jaclyn Champagne, Asantha Cooray, Kristen Coppin, Patrick Drew, Eli Dwek, Steven Finkelstein, Maximilien Franco, James Geach, Jacqueline Hodge, Maciej Koprowski, Claudia Lagos, Desika Narayanan, Alexandra Pope, David Sanders, Irene Shivaei, Sune Toft, Joaquin Vieira, Fabian Walter, Kate Whitaker , et al. (2 additional authors not shown)

Abstract: Two decades of effort have been poured into both single-dish and interferometric millimeter-wave surveys of the sky to infer the volume density of dusty star-forming galaxies (DSFGs, with SFR>100M$_\odot$ yr$^{-1}$) over cosmic time. Though obscured galaxies dominate cosmic star-formation near its peak at $z\sim2$, the contribution of such heavily obscured galaxies to cosmic star-formation is unkn… ▽ More Two decades of effort have been poured into both single-dish and interferometric millimeter-wave surveys of the sky to infer the volume density of dusty star-forming galaxies (DSFGs, with SFR>100M$_\odot$ yr$^{-1}$) over cosmic time. Though obscured galaxies dominate cosmic star-formation near its peak at $z\sim2$, the contribution of such heavily obscured galaxies to cosmic star-formation is unknown beyond $z\sim2.5$ in contrast to the well-studied population of Lyman-break galaxies (LBGs) studied through deep, space- and ground-based pencil beam surveys in the near-infrared. Unlocking the volume density of DSFGs beyond $z>3$, particularly within the first 1 Gyr after the Big Bang is critical to resolving key open questions about early Universe galaxy formation: (1) What is the integrated star-formation rate density of the Universe in the first few Gyr and how is it distributed among low-mass galaxies (e.g. Lyman-break galaxies) and high-mass galaxies (e.g. DSFGs and quasar host galaxies)? (2) How and where do the first massive galaxies assemble? (3) What can the most extreme DSFGs teach us about the mechanisms of dust production (e.g. supernovae, AGB stars, grain growth in the ISM) <1 Gyr after the Big Bang? We summarize the types of observations needed in the next decade to address these questions. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: Concept Paper submitted to Astro2020

The Evolution of Dust Opacity in Core Collapse Supernovae and the Rapid Formation of Dust in Their Ejecta

Authors: Eli Dwek, Arkaprabha Sarangi, Richard. G. Arendt

Abstract: Infrared (IR) observations of core collapse supernovae (CCSNe) have been used to infer the mass of dust that has formed in their ejecta. A plot of inferred dust masses versus supernova (SN) ages shows a trend of increasing dust mass with time, spanning a few decades of observations. This trend has been interpreted as evidence for the slow and gradual formation of dust in CCSNe. Observationally, th… ▽ More Infrared (IR) observations of core collapse supernovae (CCSNe) have been used to infer the mass of dust that has formed in their ejecta. A plot of inferred dust masses versus supernova (SN) ages shows a trend of increasing dust mass with time, spanning a few decades of observations. This trend has been interpreted as evidence for the slow and gradual formation of dust in CCSNe. Observationally, the trend exhibits a $ t^2$ behavior, exactly what is expected from an expanding optically thick ejecta. In this case, the observed dust resides in the infrared (IR)-thin "photosphere" of the ejecta, and constitutes only a fraction of the total dust mass. We therefore propose that dust formation proceeds very rapidly, condensing most available refractory elements within two years after the explosion. At early epochs, only a fraction of the dust emission escapes the ejecta accounting for the low observed dust mass. The ejecta's entire dust content is unveiled only a few decades after the explosion, with the gradual decrease in its IR opacity. Corroborating evidence for this picture includes the early depletions of refractory elements in the ejecta of SN1987A and the appearance of a silicate emission band around day 300 in SN~2004et. △ Less

Submitted 19 December, 2018; originally announced December 2018.

Comments: Accepted for publication in the ApJ Letters. The manuscript contains 7 pages, including 2 figures

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

Delayed Shock-induced Dust Formation in the Dense Circumstellar Shell Surrounding the Type IIn Supernova SN 2010jl

Authors: Arkaprabha Sarangi, Eli Dwek, Richard G Arendt

Abstract: The light curves of Type IIn supernovae are dominated by the radiative energy released through the interaction of the supernova shockwaves with their dense circumstellar medium (CSM). The ultraluminous Type IIn supernova SN 2010jl exhibits an infrared emission component that is in excess of the extrapolated UV-optical spectrum as early as a few weeks post-explosion. This emission has been attribut… ▽ More The light curves of Type IIn supernovae are dominated by the radiative energy released through the interaction of the supernova shockwaves with their dense circumstellar medium (CSM). The ultraluminous Type IIn supernova SN 2010jl exhibits an infrared emission component that is in excess of the extrapolated UV-optical spectrum as early as a few weeks post-explosion. This emission has been attributed by some as evidence for rapid formation of dust in the cooling postshock CSM. We investigate the physical processes that may inhibit or facilitate the formation of dust in the CSM. When only radiative cooling is considered, the temperature of the dense shocked gas rapidly drops below the dust condensation temperature. However, by accounting for the heating of the postshock gas by the downstream radiation from the shock, we show that dust formation is inhibited until the radiation from the shock weakens, as the shock propagates into the less dense outer regions of the CSM. In SN 2010jl dust formation can therefore only commence after day 380. Only the IR emission since that epoch can be attributed to the newly formed CSM dust. Observations on day 460 and later show that the IR luminosity exceeds the UV-optical luminosity. The post-shock dust cannot extinct the radiation emitted by the expanding SN shock. Therefore, its IR emission must be powered by an interior source, which we identify with the reverse shock propagating through the SN ejecta. IR emission before day 380 must therefore be an IR echo from preexisting CSM dust. △ Less

Submitted 18 April, 2018; originally announced April 2018.

Comments: Accepted for publication in The Astrophysical Journal, 26 Pages

On the Spatially Resolved Star Formation History in M51 I: Hybrid UV+IR Star Formation Laws and IR Emission from Dust Heated by Old Stars

Authors: R. T. Eufrasio, B. D. Lehmer, A. Zezas, E. Dwek, R. G. Arendt, A. Basu-Zych, T. Wiklind, M. Yukita, T. Fragos, A. E. Hornschemeier, L. Markwardt, A. Ptak, P. Tzanavaris

Abstract: We present Lightning, a new spectral energy distribution (SED) fitting procedure, capable of quickly and reliably recovering star formation history (SFH) and extinction parameters. The SFH is modeled as discrete steps in time. In this work, we assumed lookback times of 0-10 Myr, 10-100 Myr, 0.1-1 Gyr, 1-5 Gyr, and 5-13.6 Gyr. Lightning consists of a fully vectorized inversion algorithm to determin… ▽ More We present Lightning, a new spectral energy distribution (SED) fitting procedure, capable of quickly and reliably recovering star formation history (SFH) and extinction parameters. The SFH is modeled as discrete steps in time. In this work, we assumed lookback times of 0-10 Myr, 10-100 Myr, 0.1-1 Gyr, 1-5 Gyr, and 5-13.6 Gyr. Lightning consists of a fully vectorized inversion algorithm to determine SFH step intensities and combines this with a grid-based approach to determine three extinction parameters. We apply our procedure to the extensive FUV-to-FIR photometric data of M51, convolved to a common spatial resolution and pixel scale, and make the resulting maps publicly available. We recover, for M51a, a peak star formation rate (SFR) between 0.1 and 5 Gyr ago, with much lower star formation activity over the last 100 Myr. For M51b, we find a declining SFR toward the present day. In the outskirt regions of M51a, which includes regions between M51a and M51b, we recover a SFR peak between 0.1 and 1 Gyr ago, which corresponds to the effects of the interaction between M51a and M51b. We utilize our results to (1) illustrate how UV+IR hybrid SFR laws vary across M51, and (2) provide first-order estimates for how the IR luminosity per unit stellar mass varies as a function of the stellar age. From the latter result, we find that IR emission from dust heated by stars is not always associated with young stars, and that the IR emission from M51b is primarily powered by stars older than 5 Gyr. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Comments: Accepted for publication in ApJ

Constraints on the Progenitor of SN 2010jl and Pre-Existing Hot Dust in its Surrounding Medium

Authors: Eli Dwek, Richard G. Arendt, Ori D. Fox, Patrick L. Kelly, Nathan Smith, Schuyler D. Van Dyk, Alexei V. Filippenko, Jennifer Andrews, Isaac Shivvers

Abstract: A search for the progenitor of SN~2010jl, an unusually luminous core-collapse supernova of Type~IIn, using pre-explosion {\it Hubble}/WFPC2 and {\it Spitzer}/IRAC images of the region, yielded upper limits on the UV and near-infrared (IR) fluxes from any candidate star. These upper limits constrain the luminosity and effective temperature of the progenitor, the mass of any preexisting dust in its… ▽ More A search for the progenitor of SN~2010jl, an unusually luminous core-collapse supernova of Type~IIn, using pre-explosion {\it Hubble}/WFPC2 and {\it Spitzer}/IRAC images of the region, yielded upper limits on the UV and near-infrared (IR) fluxes from any candidate star. These upper limits constrain the luminosity and effective temperature of the progenitor, the mass of any preexisting dust in its surrounding circumstellar medium (CSM), and dust proximity to the star. A {\it lower} limit on the CSM dust mass is required to hide a luminous progenitor from detection by {\it Hubble}. {\it Upper} limits on the CSM dust mass and constraints on its proximity to the star are set by requiring that the absorbed and reradiated IR emission not exceed the IRAC upper limits. Using the combined extinction-IR emission constraints we present viable $M_d-R_1$ combinations, where $M_d$ and $R_1$ are the CSM dust mass and its inner radius. These depend on the CSM outer radius, dust composition and grain size, and the properties of the progenitor. The results constrain the pre-supernova evolution of the progenitor, and the nature and origin of the observed post-explosion IR emission from SN~2010jl. In particular, an $η$~Car-type progenitor will require at least 4~mag of visual extinction to avoid detection by the {\it Hubble}. This can be achieved with dust masses $\gtrsim 10^{-3}$~\msun\ (less than the estimated 0.2-0.5~\msun\ around $η$~Car) which must be located at distances of $\gtrsim 10^{16}$~cm from the star to avoid detection by {\it Spitzer}. △ Less

Submitted 16 August, 2017; originally announced August 2017.

Comments: Accepted for publication in the ApJ. 14 pages 10 figures. The complete figure set for Figure 10 (24 images) is available in the online journal

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

A Massive Shell of Supernova-formed Dust in SNR G54.1+0.3

Authors: Tea Temim, Eli Dwek, Richard G. Arendt, Kazimierz J. Borkowski, Stephen P. Reynolds, Patrick Slane, Joseph D. Gelfand, John C. Raymond

Abstract: While theoretical dust condensation models predict that most refractory elements produced in core-collapse supernovae (SNe) efficiently condense into dust, a large quantity of dust has so far only been observed in SN 1987A. We present the analysis of Spitzer Space Telescope, Herschel Space Observatory, Stratospheric Observatory for Infrared Astronomy (SOFIA), and AKARI observations of the infrared… ▽ More While theoretical dust condensation models predict that most refractory elements produced in core-collapse supernovae (SNe) efficiently condense into dust, a large quantity of dust has so far only been observed in SN 1987A. We present the analysis of Spitzer Space Telescope, Herschel Space Observatory, Stratospheric Observatory for Infrared Astronomy (SOFIA), and AKARI observations of the infrared (IR) shell surrounding the pulsar wind nebula in the supernova remnant G54.1+0.3. We attribute a distinctive spectral feature at 21 $μ$m to a magnesium silicate grain species that has been invoked in modeling the ejecta-condensed dust in Cas A, which exhibits the same spectral signature. If this species is responsible for producing the observed spectral feature and accounts for a significant fraction of the observed IR continuum, we find that it would be the dominant constituent of the dust in G54.1+0.3, with possible secondary contributions from other compositions, such as carbon, silicate, or alumina grains. The smallest mass of SN-formed dust required by our models is 1.1 $\pm$ 0.8 $\rm M_{\odot}$. We discuss how these results may be affected by varying dust grain properties and self-consistent grain heating models. The spatial distribution of the dust mass and temperature in G54.1+0.3 confirms the scenario in which the SN-formed dust has not yet been processed by the SN reverse shock and is being heated by stars belonging to a cluster in which the SN progenitor exploded. The dust mass and composition suggest a progenitor mass of 16$-$27 $\rm M_{\odot}$ and imply a high dust condensation efficiency, similar to that found for Cas A and SN 1987A. The study provides another example of significant dust formation in a Type IIP SN and sheds light on the properties of pristine SN-condensed dust. △ Less

Submitted 4 January, 2017; originally announced January 2017.

Comments: 16 pages, 3 tables, 11 figures, accepted for publication in ApJ

The Candidate Progenitor of the Type IIn SN 2010jl is Not an Optically Luminous Star

Authors: Ori D. Fox, Schuyler D. Van Dyk, Eli Dwek, Nathan Smith, Alexei V. Filippenko, Jennifer Andrews, Richard G. Arendt, Ryan J. Foley, Patrick L. Kelly, Adam A. Miller, Isaac Shivvers

Abstract: The nature of the progenitor star (or system) for the Type IIn supernova (SN) subclass remains uncertain. While there are direct imaging constraints on the progenitors of at least four Type IIn supernovae, one of them being SN 2010jl, ambiguities remain in the interpretation of the unstable progenitors and the explosive events themselves. A blue source in pre-explosion HST/WFPC2 images falls withi… ▽ More The nature of the progenitor star (or system) for the Type IIn supernova (SN) subclass remains uncertain. While there are direct imaging constraints on the progenitors of at least four Type IIn supernovae, one of them being SN 2010jl, ambiguities remain in the interpretation of the unstable progenitors and the explosive events themselves. A blue source in pre-explosion HST/WFPC2 images falls within the 5 sigma astrometric error circle derived from post-explosion ground-based imaging of SN 2010jl. At the time the ground-based astrometry was published, however, the SN had not faded sufficiently for post-explosion HST follow-up observations to determine a more precise astrometric solution and/or confirm if the pre-explosion source had disappeared, both of which are necessary to ultimately disentangle the possible progenitor scenarios. Here we present HST/WFC3 imaging of the SN 2010jl field obtained in 2014 and 2015, when the SN had faded sufficiently to allow for new constraints on the progenitor. The SN, which is still detected in the new images, is offset by 0."099 +/- 0."008 (24 +/- 2 pc) from the underlying and extended source of emission that contributes at least partially, if not entirely, to the blue source previously suggested as the candidate progenitor in the WFPC2 data. This point alone rules out the possibility that the blue source in the pre-explosion images is the exploding star, but may instead suggest an association with a young (<5-6 Myr) cluster and still argues for a massive (>30 solar masses) progenitor. We obtain new upper limits on the flux from a single star at the SN position in the pre-explosion WFPC2 and Spitzer/IRAC images that may ultimately be used to constrain the progenitor properties. △ Less

Submitted 1 November, 2016; originally announced November 2016.

Comments: Submitted to ApJ/AAS Journals; 6 Pages; 4 figures

Chandra Observes the End of an Era SN 1987A

Authors: Kari A. Frank, Svetozar A. Zhekov, Sangwook Park, Richard McCray, Eli Dwek, David N. Burrows

Abstract: Updated imaging and photometric results from Chandra observations of SN 1987A, covering the last 16 years, are presented. We find that the 0.5-2 keV light curve has remained constant at ~8x10^-12 erg s^-1 cm^-2 since 9500 days, with the 3-8 keV light curve continuing to increase until at least 10000 days. The expansion rate of the ring is found to be energy dependent, such that after day 6000 the… ▽ More Updated imaging and photometric results from Chandra observations of SN 1987A, covering the last 16 years, are presented. We find that the 0.5-2 keV light curve has remained constant at ~8x10^-12 erg s^-1 cm^-2 since 9500 days, with the 3-8 keV light curve continuing to increase until at least 10000 days. The expansion rate of the ring is found to be energy dependent, such that after day 6000 the ring expands faster in the 2-10 keV band than it does at energies <2 keV. Images show a reversal of the east-west asymmetry between 7000 and 8000 days after the explosion. The latest images suggest the southeastern side of the equatorial ring is beginning to fade. Consistent with the latest optical and infrared results, our Chandra analysis indicates the blast wave is now leaving the dense equatorial ring, which marks the beginning of a major change in the evolutionary phase of the supernova remnant 1987A. △ Less

Submitted 6 August, 2016; originally announced August 2016.

Iron: A Key Element for Understanding the Origin and Evolution of Interstellar Dust

Authors: Eli Dwek

Abstract: The origin and depletion of iron differ from all other abundant refractory elements that make up the composition of the interstellar dust. Iron is primarily synthesized in Type Ia supernovae (SNe Ia) and in core collapse supernovae (CCSN), and is present in the outflows from AGB stars. Only the latter two are observed to be sources of interstellar dust, since searches for dust in SN~Ia have provid… ▽ More The origin and depletion of iron differ from all other abundant refractory elements that make up the composition of the interstellar dust. Iron is primarily synthesized in Type Ia supernovae (SNe Ia) and in core collapse supernovae (CCSN), and is present in the outflows from AGB stars. Only the latter two are observed to be sources of interstellar dust, since searches for dust in SN~Ia have provided strong evidence for the absence of any significant mass of dust in their ejecta. Consequently, more than 65% of the iron is injected into the ISM in gaseous form. Yet, ultraviolet and X-ray observations along many lines of sight in the ISM show that iron is severely depleted in the gas phase compared to expected solar abundances. The missing iron, comprising about 90% of the total, is believed to be locked up in interstellar dust. This suggests that most of the missing iron must have precipitated from the ISM gas by cold accretion onto preexisting silicate, carbon, or composite grains. Iron is thus the only element that requires most of its growth to occur outside the traditional stellar condensation sources. This is a robust statement that does not depend on our evolving understanding of the dust destruction efficiency in the ISM. Reconciling the physical, optical, and chemical properties of such composite grains with their many observational manifestations is a major challenge for understanding the nature and origin of interstellar dust. △ Less

Submitted 6 May, 2016; originally announced May 2016.

Comments: Accepted for publication in The Astrophysical Journal (8 pages, 7 figures)

Dust Destruction by the Reverse Shock in the Cassiopeia A Supernova Remnant

Authors: Elisabetta R. Micelotta, Eli Dwek, Jonathan D. Slavin

Abstract: Core collapse supernovae (CCSNe) are important sources of interstellar dust, potentially capable of producing one solar mass of dust in their explosively expelled ejecta. However, unlike other dust sources, the dust has to survive the passage of the reverse shock, generated by the interaction of the supernova blast wave with its surrounding medium. Knowledge of the net amount of dust produced by C… ▽ More Core collapse supernovae (CCSNe) are important sources of interstellar dust, potentially capable of producing one solar mass of dust in their explosively expelled ejecta. However, unlike other dust sources, the dust has to survive the passage of the reverse shock, generated by the interaction of the supernova blast wave with its surrounding medium. Knowledge of the net amount of dust produced by CCSNe is crucial for understanding the origin and evolution of dust in the local and high-redshift universe. Our aim is to identify the dust destruction mechanisms in the ejecta, and derive the net amount of dust that survives the passage of the reverse shock. We use analytical models for the evolution of a supernova blast wave and of the reverse shock, with special application to the clumpy ejecta of the remnant of Cassiopeia A. We assume that the dust resides in cool oxygen-rich clumps that are uniformly distributed within the remnant and surrounded by a hot X-ray emitting plasma, and that the dust consists of silicates (MgSiO3) and amorphous carbon grains. The passage of the reverse shock through the clumps gives rise to a relative gas-grain motion and also destroys the clumps. Inside the ejecta clouds, dust is processed via kinetic sputtering, which is terminated either when the grains escape the clumps, or when the clumps are destroyed by the reverse shock. In either case, grain destruction proceeds thereafter by thermal sputtering in the hot ambient gas. We find that 11.8% and 15.9% of, respectively, the silicate and carbon dust survives the passage of the reverse shock by the time the shock has reached the center of the remnant. These fractions depend on the morphology of the ejecta and the medium into which the remnant is expanding, as well as the composition and size distribution of the grains that formed in the ejecta. Results will therefore differ for different types of supernovae. △ Less

Submitted 3 May, 2016; v1 submitted 8 February, 2016; originally announced February 2016.

Comments: A&A, in press, 20 pages (including two appendices), 15+1 figures, one additional reference in the Introduction, two paragraphs added to the Discussion section

Journal ref: A&A 590, A65 (2016)

Infrared Continuum and Line Evolution of the Equatorial Ring around SN 1987A

Authors: Richard G. Arendt, Eli Dwek, Patrice Bouchet, I. John Danziger, Kari A. Frank, Robert D. Gehrz, Sangwook Park, Charles E. Woodward

Abstract: Spitzer observations of SN 1987A have now spanned more than a decade. Since day ~4,000, mid-infrared (mid-IR) emission has been dominated by that from shock-heated dust in the equatorial ring (ER). From 6,000 to 8,000 days after the explosion, Spitzer observations included broadband photometry at 3.6 - 24 micron, and low and moderate resolution spectroscopy at 5 - 35 micron. Here we present later… ▽ More Spitzer observations of SN 1987A have now spanned more than a decade. Since day ~4,000, mid-infrared (mid-IR) emission has been dominated by that from shock-heated dust in the equatorial ring (ER). From 6,000 to 8,000 days after the explosion, Spitzer observations included broadband photometry at 3.6 - 24 micron, and low and moderate resolution spectroscopy at 5 - 35 micron. Here we present later Spitzer observations, through day 10,377, which include only the broadband measurements at 3.6 and 4.5 micron. These data show that the 3.6 and 4.5 micron brightness has clearly begun to fade after day ~8,500, and no longer tracks the X-ray emission as well as it did at earlier epochs. This can be explained by the destruction of the dust in the ER on time scales shorter than the cooling time for the shocked gas. We find that the evolution of the late time IR emission is also similar to the now fading optical emission. We provide the complete record of the IR emission lines, as seen by Spitzer prior to day 8,000. The past evolution of the gas as seen by the IR emission lines seems largely consistent with the optical emission, although the IR [Fe II] and [Si II] lines show different, peculiar velocity structures. △ Less

Submitted 11 January, 2016; originally announced January 2016.

Comments: Accepted for publication in AJ. 22 pages, 18 figures

Submillimeter Observations of CLASH 2882 and the Evolution of Dust in this Galaxy

Authors: Eli Dwek, Johannes Staguhn, Richard G. Arendt, Attila Kovács, Roberto Decarli, Eiichi Egami, Michał J. Michałowski, Timothy D. Rawle, Sune Toft, Fabian Walter

Abstract: Two millimeter observations of the MACS J1149.6+2223 cluster have detected a source that was consistent with the location of the lensed MACS1149-JD galaxy at z=9.6. A positive identification would have rendered this galaxy as the youngest dust forming galaxy in the universe. Follow up observation with the AzTEC 1.1 mm camera and the IRAM NOrthern Extended Millimeter Array (NOEMA) at 1.3 mm have no… ▽ More Two millimeter observations of the MACS J1149.6+2223 cluster have detected a source that was consistent with the location of the lensed MACS1149-JD galaxy at z=9.6. A positive identification would have rendered this galaxy as the youngest dust forming galaxy in the universe. Follow up observation with the AzTEC 1.1 mm camera and the IRAM NOrthern Extended Millimeter Array (NOEMA) at 1.3 mm have not confirmed this association. In this paper we show that the NOEMA observations associate the 2 mm source with [PCB2012] 2882 ([PCB2012] 2882 is the NED-searchable name for this source.), source number 2882 in the Hubble Space Telescope ( HST) Cluster Lensing and Supernova (CLASH) catalog of MACS J1149.6+2223. This source, hereafter referred to as CLASH 2882, is a gravitationally lensed spiral galaxy at z=0.99. We combine the GISMO 2 mm and NOEMA 1.3 mm fluxes with other (rest frame) UV to far-IR observations to construct the full spectral energy distribution (SED) of this galaxy, and derive its star formation history, and stellar and interstellar dust content. The current star formation rate of the galaxy is 54/mu Msun yr-1, and its dust mass is about 5 10^7/mu Msun, where mu is the lensing magnification factor for this source, which has a mean value of 2.7. The inferred dust mass is higher than the maximum dust mass that can be produced by core collapse supernovae (CCSN) and evolved AGB stars. As with many other star forming galaxies, most of the dust mass in CLASH 2882 must have been accreted in the dense phases of the ISM. △ Less

Submitted 30 September, 2015; originally announced September 2015.

Comments: Accepted for publication in The Astrophysical Journal

What Powers the 3000-Day Light Curve of SN 2006gy?

Authors: Ori D. Fox, Nathan Smith, S. Mark Ammons, Jennifer Andrews, K. Azalee Bostroem, S. Bradley Cenko, Geoffrey C. Clayton, Eli Dwek, Alexei V. Filippenko, Joseph S. Gallagher, Patrick L. Kelly, Jon C. Mauerhan, Adam M. Miller, Schuyler D. Van Dyk

Abstract: SN 2006gy was the most luminous SN ever observed at the time of its discovery and the first of the newly defined class of superluminous supernovae (SLSNe). The extraordinary energetics of SN 2006gy and all SLSNe (>10^51 erg) require either atypically large explosion energies (e.g., pair-instability explosion) or the efficient conversion of kinetic into radiative energy (e.g., shock interaction). T… ▽ More SN 2006gy was the most luminous SN ever observed at the time of its discovery and the first of the newly defined class of superluminous supernovae (SLSNe). The extraordinary energetics of SN 2006gy and all SLSNe (>10^51 erg) require either atypically large explosion energies (e.g., pair-instability explosion) or the efficient conversion of kinetic into radiative energy (e.g., shock interaction). The mass-loss characteristics can therefore offer important clues regarding the progenitor system. For the case of SN 2006gy, both a scattered and thermal light echo from circumstellar material (CSM) have been reported at later epochs (day ~800), ruling out the likelihood of a pair-instability event and leading to constraints on the characteristics of the CSM. Owing to the proximity of the SN to the bright host-galaxy nucleus, continued monitoring of the light echo has not been trivial, requiring the high resolution offered by the Hubble Space Telescope (HST) or ground-based adaptive optics (AO). Here we report detections of SN 2006gy using HST and Keck AO at ~3000 days post-explosion and consider the emission mechanism for the very late-time light curve. While the optical light curve and optical spectral energy distribution are consistent with a continued scattered-light echo, a thermal echo is insufficient to power the K'-band emission by day 3000. Instead, we present evidence for late-time infrared emission from dust that is radiatively heated by CSM interaction within an extremely dense dust shell, and we consider the implications on the CSM characteristics and progenitor system. △ Less

Submitted 21 September, 2015; originally announced September 2015.

Comments: Submitted to MNRAS, 14 pages

The Evolution of Dust Mass in the Ejecta of SN1987A

Authors: Eli Dwek, Richard G. Arendt

Abstract: We present a new analysis of the infrared (IR) emission from the ejecta of SN1987A covering days 615, 775, 1144, 8515, and 9090 after the explosion. We show that the observations are consistent with the rapid formation of about 0.4 Msun of dust, consisting of mostly silicates, near day 615, and evolving to about 0.45 Msun of composite dust grains consisting of ~0.4 Msun of silicates and ~ 0.05 Msu… ▽ More We present a new analysis of the infrared (IR) emission from the ejecta of SN1987A covering days 615, 775, 1144, 8515, and 9090 after the explosion. We show that the observations are consistent with the rapid formation of about 0.4 Msun of dust, consisting of mostly silicates, near day 615, and evolving to about 0.45 Msun of composite dust grains consisting of ~0.4 Msun of silicates and ~ 0.05 Msun of amorphous carbon after day ~8500. The proposed scenario challenges previous claims that dust in SN ejecta is predominantly carbon, and that it grew from an initial mass of ~1e-3 Msun, to over 0.5 Msun by cold accretion. It alleviates several problems with previous interpretations of the data: (1) it reconciles the abundances of silicon, magnesium, and carbon with the upper limits imposed by nucleosynthesis calculations; (2) it eliminates the requirement that most of the dust observed around day 9000 grew by cold accretion onto the1e-3 Msun of dust previously inferred for days 615 and 775 after the explosion; and (3) establishes the dominance of silicate over carbon dust in the SN ejecta. At early epochs, the IR luminosity of the dust is powered by the radioactive decay of 56Co, and at late times by at least (1.3-1.6)e-4 Msun 44Ti. Even if only a fraction greater than ~10% of the silicate dust survives the injection into the ISM, the observations firmly establish the role of core collapse SNe as the major source of thermally condensed silicate dust in the universe. △ Less

Submitted 24 July, 2015; originally announced July 2015.

Comments: Accepted for publication in The Astrophysical Journal. The manuscript has 10 pages, 6 figures, and 4 tables

Spectral Confusion for Cosmological Surveys of Redshifted CII Observations

Authors: A. Kogut, E. Dwek, S. H. Moseley

Abstract: Far infrared cooling lines are ubiquitous features in the spectra of star forming galaxies. Surveys of redshifted fine-structure lines provide a promising new tool to study structure formation and galactic evolution at redshifts including the epoch of reionization as well as the peak of star formation. Unlike neutral hydrogen surveys, where the 21 cm line is the only bright line, surveys of red-sh… ▽ More Far infrared cooling lines are ubiquitous features in the spectra of star forming galaxies. Surveys of redshifted fine-structure lines provide a promising new tool to study structure formation and galactic evolution at redshifts including the epoch of reionization as well as the peak of star formation. Unlike neutral hydrogen surveys, where the 21 cm line is the only bright line, surveys of red-shifted fine-structure lines suffer from confusion generated by line broadening, spectral overlap of different lines, and the crowding of sources with redshift. We use simulations to investigate the resulting spectral confusion and derive observing parameters to minimize these effects in pencil-beam surveys of red-shifted far-IR line emission. We generate simulated spectra of the 17 brightest far-IR lines in galaxies, covering the 150 to 1300 micron wavelength region corresponding to redshifts 0 < z < 7, and develop a simple iterative algorithm that successfully identifies the 158 micron [CII] line and other lines. Although the [CII] line is a principal coolant for the interstellar medium, the assumption that the brightest observed lines in a given line of sight are always [CII] lines is a poor approximation to the simulated spectra once other lines are included. Blind line identification requires detection of fainter companion lines from the same host galaxies, driving survey sensitivity requirements. The observations require moderate spectral resolution 700 < R < 4000 with angular resolution between 20 arcsec and 10 armin, sufficiently narrow to minimize confusion yet sufficiently large to include a statistically meaningful number of sources. △ Less

Submitted 1 May, 2015; originally announced May 2015.

Comments: 7 pages including 7 figures. Accepted for publication in The Astrophysical Journal

Elemental depletions in the Magellanic Clouds and the evolution of depletions with metallicity

Authors: Kirill Tchernyshyov, Margaret Meixner, Jonathan Seale, Andrew Fox, Scott D. Friedman, Eli Dwek, Frédéric Galliano, Kenneth Sembach

Abstract: We present a study of the composition of gas and dust in the Large and Small Magellanic Clouds (LMC and SMC, together -- the MCs) as measured by UV absorption spectroscopy. We have measured P II and Fe II along 85 sightlines toward the MCs using archival FUSE observations. For 16 of those sightlines, we have measured Si II, Cr II, and Zn II from new HST COS observations. We have combined these mea… ▽ More We present a study of the composition of gas and dust in the Large and Small Magellanic Clouds (LMC and SMC, together -- the MCs) as measured by UV absorption spectroscopy. We have measured P II and Fe II along 85 sightlines toward the MCs using archival FUSE observations. For 16 of those sightlines, we have measured Si II, Cr II, and Zn II from new HST COS observations. We have combined these measurements with H I and H$_2$ column densities and reference stellar abundances from the literature to derive gas-phase abundances, depletions, and gas-to-dust ratios (GDRs). 80 of our 84 P measurements and 13 of our 16 Zn measurements are depleted by more than 0.1 decades, suggesting that P and Zn abundances are not accurate metallicity indicators at and above the metallicity of the SMC. The maximum P and Zn depletions are the same in the MW, LMC, and SMC. Si, Cr, and Fe are systematically less depleted in the SMC than in the MW or LMC. The minimum Si depletion in the SMC is consistent with zero. Our depletion-derived GDRs broadly agree with GDRs from the literature. The GDR varies from location to location within a galaxy by a factor of up to 2 in the LMC and up to 5 in the SMC. This variation is evidence of dust destruction and/or growth in the diffuse neutral phase of the interstellar medium. △ Less

Submitted 30 March, 2015; originally announced March 2015.

Comments: Submitted to ApJ

Destruction of Interstellar Dust in Evolving Supernova Remnant Shock Waves

Authors: Jonathan D. Slavin, Eli Dwek, Anthony P. Jones

Abstract: Supernova generated shock waves are responsible for most of the destruction of dust grains in the interstellar medium (ISM). Calculations of the dust destruction timescale have so far been carried out using plane parallel steady shocks, however that approximation breaks down when the destruction timescale becomes longer than that for the evolution of the supernova remnant (SNR) shock. In this pape… ▽ More Supernova generated shock waves are responsible for most of the destruction of dust grains in the interstellar medium (ISM). Calculations of the dust destruction timescale have so far been carried out using plane parallel steady shocks, however that approximation breaks down when the destruction timescale becomes longer than that for the evolution of the supernova remnant (SNR) shock. In this paper we present new calculations of grain destruction in evolving, radiative SNRs. To facilitate comparison with the previous study by Jones et al. (1996), we adopt the same dust properties as in that paper. We find that the efficiencies of grain destruction are most divergent from those for a steady shock when the thermal history of a shocked gas parcel in the SNR differs significantly from that behind a steady shock. This occurs in shocks with velocities >~ 200 km/s for which the remnant is just beginning to go radiative. Assuming SNRs evolve in a warm phase dominated ISM, we find dust destruction timescales are increased by a factor of ~2 compared to those of Jones et al. (1996), who assumed a hot gas dominated ISM. Recent estimates of supernova rates and ISM mass lead to another factor of ~3 increase in the destruction timescales, resulting in a silicate grain destruction timescale of ~2-3 Gyr. These increases, while not able resolve the problem of the discrepant timescales for silicate grain destruction and creation, are an important step towards understanding the origin, and evolution of dust in the ISM. △ Less

Submitted 3 February, 2015; originally announced February 2015.

Comments: 30 pages, 8 figures, accepted for publication in the Astrophysical Journal

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

Authors: M. Matsuura, E. Dwek, M. J. Barlow, B. Babler, M. Baes, M. Meixner, Jose Cernicharo, Geoff C. Clayton, L. Dunne, C. Fransson, Jacopo Fritz, Walter Gear, H. L. Gomez, M. A. T. Groenewegen, R. Indebetouw, R. J. Ivison, A. Jerkstrand, V. Lebouteiller, T. L. Lim, P. Lundqvist, C. P. Pearson, J Roman-Duval, P. Royer, Lister Staveley-Smith, B. M. Swinyard , et al. (11 additional authors not shown)

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

Submitted 26 November, 2014; originally announced November 2014.

Comments: ApJ accepted, 8 pages

Dust Destruction Rates and Lifetimes in the Magellanic Clouds

Authors: Tea Temim, Eli Dwek, Kirill Tchernyshyov, Martha L. Boyer, Margaret Meixner, Christa Gall, Julia Roman-Duval

Abstract: The nature, composition, abundance, and size distribution of dust in galaxies is determined by the rate at which it is created in the different stellar sources and destroyed by interstellar shocks. Because of their extensive wavelength coverage, proximity, and nearly face-on geometry, the Magellanic Clouds (MCs) provide a unique opportunity to study these processes in great detail. In this paper w… ▽ More The nature, composition, abundance, and size distribution of dust in galaxies is determined by the rate at which it is created in the different stellar sources and destroyed by interstellar shocks. Because of their extensive wavelength coverage, proximity, and nearly face-on geometry, the Magellanic Clouds (MCs) provide a unique opportunity to study these processes in great detail. In this paper we use the complete sample of supernova remnants (SNRs) in the MCs to calculate the lifetime and destruction efficiencies of silicate and carbon dust in these galaxies. We find dust lifetimes of 22 +- 13 Myr (30 +- 17 Myr) for silicate (carbon) grains in the LMC, and 54 +- 32 Myr (72 +- 43 Myr) for silicate (carbon) grains in the SMC. The significantly shorter lifetimes in the MCs, as compared to the Milky Way, are explained as the combined effect of their lower total dust mass, and the fact that the dust-destroying isolated SNe in the MCs seem to be preferentially occurring in regions with higher than average dust-to-gas (D2G) mass ratios. We also calculate the supernova rate and the current star formation rate in the MCs, and use them to derive maximum dust injection rates by asymptotic giant branch stars and core collapse supernovae. We find that the injection rates are an order of magnitude lower than the dust destruction rates by the SNRs. This supports the conclusion that, unless the dust destruction rates have been considerably overestimated, most of the dust must be reconstituted from surviving grains in dense molecular clouds. More generally, we also discuss the dependence of the dust destruction rate on the local D2G mass ratio, the ambient gas density and metallicity, as well as the application of our results to other galaxies and dust evolution models. △ Less

Submitted 17 November, 2014; originally announced November 2014.

Comments: 15 pages, 8 figures, 5 tables, accepted to ApJ

Star Formation Histories Across the Interacting Galaxy NGC 6872, the Largest-Known Spiral

Authors: Rafael T. Eufrasio, Eli Dwek, Richard G. Arendt, Duilia F. de Mello, Dimitri Gadotti, Fernanda Urrutia-Viscarra, Claudia Mendes de Oliveira, Dominic Benford

Abstract: NGC 6872, hereafter the Condor, is a large spiral galaxy that is interacting with its closest companion, the S0 galaxy IC 4970. The extent of the Condor provides an opportunity for detailed investigation of the impact of the interaction on the current star formation rate and its history across the galaxy, on the age and spatial distribution of its stellar population, and on the mechanism that driv… ▽ More NGC 6872, hereafter the Condor, is a large spiral galaxy that is interacting with its closest companion, the S0 galaxy IC 4970. The extent of the Condor provides an opportunity for detailed investigation of the impact of the interaction on the current star formation rate and its history across the galaxy, on the age and spatial distribution of its stellar population, and on the mechanism that drive the star formation activity. To address these issues we analyzed the far-ultraviolet (FUV) to near-infrared (near-IR) spectral energy distribution (SED) of 17, 10 kpc diameter, regions across the galaxy, and derived their star formation history, current star formation rate, and stellar population and mass. We find that most of the star formation takes place in the extended arms, with very little star formation in the central 5 kpc of the galaxy, in contrast to what was predicted from previous numerical simulations. There is a trend of increasing star formation activity with distance from the nucleus of the galaxy, and no evidence for a recent increase in the current star formation rate (SFR) due to the interaction. The nucleus itself shows no significant current star formation activity. The extent of the Condor also provides an opportunity to test the applicability of a single standard prescription for conversion of the FUV + IR (22 micron) intensities to a star formation rate for all regions. We find that the conversion factor differs from region to region, arising from regional differences in the stellar populations. △ Less

Submitted 10 September, 2014; originally announced September 2014.

Comments: 14 pages. 10 figures. Accepted for publication in ApJ

Rapid formation of large dust grains in the luminous supernova SN 2010jl

Authors: Christa Gall, Jens Hjorth, Darach Watson, Eli Dwek, Justyn R. Maund, Ori Fox, Giorgos Leloudas, Daniele Malesani, Avril C. Day-Jones

Abstract: The origin of dust in galaxies is still a mystery. The majority of the refractory elements are produced in supernova explosions but it is unclear how and where dust grains condense and grow, and how they avoid destruction in the harsh environments of star-forming galaxies. The recent detection of 0.1-0.5 solar masses of dust in nearby supernova remnants suggests in situ dust formation, while other… ▽ More The origin of dust in galaxies is still a mystery. The majority of the refractory elements are produced in supernova explosions but it is unclear how and where dust grains condense and grow, and how they avoid destruction in the harsh environments of star-forming galaxies. The recent detection of 0.1-0.5 solar masses of dust in nearby supernova remnants suggests in situ dust formation, while other observations reveal very little dust in supernovae the first few years after explosion. Observations of the bright SN 2010jl have been interpreted as pre-existing dust, dust formation or no dust at all. Here we report the rapid (40-240 days) formation of dust in its dense circumstellar medium. The wavelength dependent extinction of this dust reveals the presence of very large (> 1 micron) grains, which are resistant to destructive processes. At later times (500-900 days), the near-IR thermal emission shows an accelerated growth in dust mass, marking the transition of the supernova from a circumstellar- to an ejecta-dominated source of dust. This provides the link between the early and late dust mass evolution in supernovae with dense circumstellar media. △ Less

Submitted 16 July, 2014; originally announced July 2014.

Comments: 62 pages, 13 figures, 1 table. Author version of the Letter to Nature, published online July 9 2014 (Nature, 511, 7509, pp. 326-329 (2014)), prior to the final editorial changes to conform to Journal style; includes Methods and Extended Data Figures and the Supplementary Information. See published version http://www.nature.com/nature/journal/v511/n7509/full/nature13558.html

Journal ref: Nature, Volume 511, Issue 7509, pp. 326-329 (17 July 2014)

Dust formation, evolution, and obscuration effects in the very high-redshift universe

Authors: Eli Dwek, Johannes Staguhn, Richard G. Arendt, Attila Kovacs, Ting Su, Dominic J. Benford

Abstract: The evolution of dust at redshifts z>9, and consequently the dust properties, differs greatly from that in the local universe. In contrast to the local universe, core collapse supernovae (CCSNe) are the only source of thermally-condensed dust. Because of the low initial dust-to-gas mass ratio, grain destruction rates are low, so that CCSNe are net producers of interstellar dust. Galaxies with larg… ▽ More The evolution of dust at redshifts z>9, and consequently the dust properties, differs greatly from that in the local universe. In contrast to the local universe, core collapse supernovae (CCSNe) are the only source of thermally-condensed dust. Because of the low initial dust-to-gas mass ratio, grain destruction rates are low, so that CCSNe are net producers of interstellar dust. Galaxies with large initial gas mass or high mass infall rate will therefore have a more rapid net rate of dust production comported to galaxies with lower gas mass, even at the same star formation rate. The dust composition is dominated by silicates, which exhibit a strong rise in the UV opacity near the Lyman break. This "silicate-UV break" may be confused with the Lyman break, resulting in a misidentification of a galaxies' photometric redshift. In this paper we demonstrate these effects by analyzing the spectral energy distribution (SED) of MACS1149-JD, a lensed galaxy at z=9.6. A potential 2mm counterpart of MACS1149-JD has been identified with GISMO. While additional observations are required to corroborate this identification, we use this possible association to illustrate the physical processes and the observational effects of dust in the very high redshift universe. △ Less

Submitted 5 May, 2014; originally announced May 2014.

Comments: Accepted for publication in ApJ Letters

Interstellar and Ejecta Dust in the Cas A Supernova Remnant

Authors: Richard G. Arendt, Eli Dwek, Gladys Kober, Jeonghee Rho, Una Hwang

Abstract: Infrared continuum observations provide a means of investigating the physical composition of the dust in the ejecta and swept up medium of the Cas A supernova remnant. Using low resolution Spitzer IRS spectra (5-35 $μ$m), and broad-band Herschel PACS imaging (70, 100, and 160 $μ$m), we identify characteristic dust spectra, associated with ejecta layers that underwent distinct nuclear burning histo… ▽ More Infrared continuum observations provide a means of investigating the physical composition of the dust in the ejecta and swept up medium of the Cas A supernova remnant. Using low resolution Spitzer IRS spectra (5-35 $μ$m), and broad-band Herschel PACS imaging (70, 100, and 160 $μ$m), we identify characteristic dust spectra, associated with ejecta layers that underwent distinct nuclear burning histories. The most luminous spectrum exhibits strong emission features at $\sim9$ and 21 $μ$m and is closely associated with ejecta knots with strong Ar emission lines. The dust features can be reproduced by magnesium silicate grains with relatively low Mg to Si ratios. Another dust spectrum is associated with ejecta having strong Ne emission lines. It has no indication of any silicate features, and is best fit by Al$_2$O$_3$ dust. A third characteristic dust spectrum shows features that are best matched by magnesium silicates with a relatively high Mg to Si ratio. This dust is primarily associated with the X-ray emitting shocked ejecta, but it is also evident in regions where shocked interstellar or circumstellar material is expected. However, the identification of dust composition is not unique, and each spectrum includes an additional featureless dust component of unknown composition. Colder dust of indeterminate composition is associated with emission from the interior of the SNR, where the reverse shock has not yet swept up and heated the ejecta. Most of the dust mass in Cas A is associated with this unidentified cold component, which is $\lesssim0.1$ $M_{\odot}$. The mass of warmer dust is only $\sim 0.04$ $M_{\odot}$. △ Less

Submitted 12 March, 2014; originally announced March 2014.

Comments: 45 pages. 21 Figures. Accepted for publication in ApJ

Journal ref: 2014 ApJ, 786, 55

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

Authors: R. Indebetouw, M. Matsuura, E. Dwek, G. Zanardo, M. J. Barlow, M. Baes, P. Bouchet, D. N. Burrows, R. Chevalier, G. C. Clayton, C. Fransson, B. Gaensler, R. Kirshner, M. Lakicevic, K. S. Long, P. Lundqvist, I. Marti-Vidal, J. Marcaide, R. McCray, M. Meixner, C. Y. Ng, S. Park, G. Sonneborn, L. Staveley-Smith, C. Vlahakis , et al. (1 additional authors not shown)

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

Submitted 14 December, 2013; originally announced December 2013.

Comments: ApJL accepted

The GISMO 2-millimeter Deep Field in GOODS-N

Authors: Johannes G. Staguhn, Attila Kovacs, Richard G. Arendt, Dominic J. Benford, Roberto Decarli, Eli Dwek, Dale J. Fixsen, Gene C. Hilton, Kent D. Irwin, Christine A. Jhabvala, Alexander Karim, Samuel Leclercq, Stephen F. Maher, Timothy M. Miller, S. Harvey Moseley, Elmer H. Sharp, Fabian Walter, Edward J. Wollack

Abstract: We present deep continuum observations using the GISMO camera at a wavelength of 2 mm centered on the Hubble Deep Field (HDF) in the GOODS-N field. These are the first deep field observations ever obtained at this wavelength. The 1 sigma sensitivity in the innermost approx. 4 arcminutes of the 7 utes map is approx. 135 uJy/beam, a factor of three higher in flux/beam sensitivity than the deepest av… ▽ More We present deep continuum observations using the GISMO camera at a wavelength of 2 mm centered on the Hubble Deep Field (HDF) in the GOODS-N field. These are the first deep field observations ever obtained at this wavelength. The 1 sigma sensitivity in the innermost approx. 4 arcminutes of the 7 utes map is approx. 135 uJy/beam, a factor of three higher in flux/beam sensitivity than the deepest available SCUBA 850 um observations, and almost a factor of four higher in flux / beam sensitivity than the combined MAMBO/AzTEC 1.2 mm observations of this region. Our source extraction algorithm identifies 12 sources directly, and another 3 through correlation with known sources at 1.2 mm and 850 um. Five of the directly detected GISMO sources have counterparts in the MAMBO/AzTEC catalog, and four of those also have SCUBA counterparts. HDF850.1, one of the first blank-field detected submillimeter galaxies, is now detected at 2 mm. The median redshift of all sources with counterparts of known redshifts is med(z) = 2.91 +/- 0.94. Statistically, the detections are most likely real for 5 of the seven 2 mm sources without shorter wavelength counterparts, while the probability for none of them being real is negligible. △ Less

Submitted 7 June, 2014; v1 submitted 6 November, 2013; originally announced November 2013.

Comments: Accepted for publication in ApJ

The suppression of star formation by powerful active galactic nuclei

Authors: M. J. Page, M. Symeonidis, J. D. Vieira, B. Altieri, A. Amblard, V. Arumugam, H. Aussel, T. Babbedge, A. Blain, J. Bock, A. Boselli, V. Buat, N. Castro-Rodr'iguez, A. Cava, P. Chanial, D. L. Clements, A. Conley, L. Conversi, A. Cooray, C. D. Dowell, E. N. Dubois, J. S. Dunlop, E. Dwek, S. Dye, S. Eales , et al. (53 additional authors not shown)

Abstract: The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results… ▽ More The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time. △ Less

Submitted 15 October, 2013; originally announced October 2013.

Comments: Published in Nature, May 2012. Includes supplementary information

The Importance of Physical Models for Deriving Dust Masses and Grain Size Distributions in Supernova Ejecta I: Radiatively Heated Dust in the Crab Nebula

Authors: Tea Temim, Eli Dwek

Abstract: Recent far-infrared (IR) observations of supernova remnants (SNRs) have revealed significantly large amounts of newly-condensed dust in their ejecta, comparable to the total mass of available refractory elements. The dust masses derived from these observations assume that all the grains of a given species radiate at the same temperature, regardless of the dust heating mechanism or grain radius. In… ▽ More Recent far-infrared (IR) observations of supernova remnants (SNRs) have revealed significantly large amounts of newly-condensed dust in their ejecta, comparable to the total mass of available refractory elements. The dust masses derived from these observations assume that all the grains of a given species radiate at the same temperature, regardless of the dust heating mechanism or grain radius. In this paper, we derive the dust mass in the ejecta of the Crab Nebula, using a physical model for the heating and radiation from the dust. We adopt a power-law distribution of grain sizes and two different dust compositions (silicates and amorphous carbon), and calculate the heating rate of each dust grain by the radiation from the pulsar wind nebula (PWN). We find that the grains attain a continuous range of temperatures, depending on their size and composition. The total mass derived from the best-fit models to the observed IR spectrum is 0.019-0.13 solar masses, depending on the assumed grain composition. We find that the power-law size distribution of dust grains is characterized by a power-law index of 3.5-4.0 and a maximum grain size larger than 0.1 microns. The grain sizes and composition are consistent with what is expected for dust grains formed in a Type IIP SN. Our derived dust mass is at least a factor of two less than the mass reported in previous studies of the Crab Nebula that assumed more simplified two-temperature models. The results of this study show that a physical model resulting in a realistic distribution of dust temperatures can constrain the dust properties and affect the derived dust masses. Our study may also have important implications for deriving grain properties and mass estimates in other SNRs and for the ultimate question of whether SNe are major sources of dust in the Galactic interstellar medium (ISM) and in external galaxies. △ Less

Submitted 3 July, 2013; v1 submitted 21 February, 2013; originally announced February 2013.

Comments: 9 pages, 2 tables, 8 figures, Accepted to The Astrophysical Journal

The Herschel PEP/HerMES Luminosity Function. I: Probing the Evolution of PACS selected Galaxies to z~4

Authors: C. Gruppioni, F. Pozzi, G. Rodighiero, I. Delvecchio, S. Berta, L. Pozzetti, G. Zamorani, P. Andreani, A. Cimatti, O. Ilbert, E. Le Floch, D. Lutz, B. Magnelli, L. Marchetti, P. Monaco, R. Nordon, S. Oliver, P. Popesso, L. Riguccini, I. Roseboom, D. J. Rosario, M. Sargent, M. Vaccari, B. Altieri, H. Aussel , et al. (78 additional authors not shown)

Abstract: We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity… ▽ More We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies. △ Less

Submitted 25 July, 2013; v1 submitted 21 February, 2013; originally announced February 2013.

Comments: 32 pages, 20 figures, 9 tables. Published in MNRAS. Replaced Fig. 8 (small scaling bug in the previous version)