Source author record

Eli Dwek

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

34 published item(s)

preprint2021arXiv

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

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.

preprint2020arXiv

The origin of dust in galaxies across cosmic time

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.

preprint2019arXiv

Dust formation in AGN winds

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.

preprint2019arXiv

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

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.

preprint2016arXiv

Chandra Observes the End of an Era SN 1987A

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2015arXiv

Destruction of Interstellar Dust in Evolving Supernova Remnant Shock Waves

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

The Evolution of Dust Mass in the Ejecta of SN1987A

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.

preprint2015arXiv

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

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.

preprint2014arXiv

Dust Destruction Rates and Lifetimes in the Magellanic Clouds

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.

preprint2014arXiv

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

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.

preprint2014arXiv

Interstellar and Ejecta Dust in the Cas A Supernova Remnant

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}$.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

The GISMO 2-millimeter Deep Field in GOODS-N

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.

preprint2013arXiv

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

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.

preprint2012arXiv

Properties and Spatial Distribution of Dust Emission in the Crab Nebula

Recent infrared (IR) observations of freshly-formed dust in supernova remnants (SNRs) have yielded significantly lower dust masses than predicted by theoretical models and measured from high redshift observations. The Crab Nebula's pulsar wind is thought to be sweeping up freshly-formed supernova (SN) dust along with the ejected gas. The evidence for this dust was found in the form of an IR excess in the integrated spectrum of the Crab and in extinction against the synchrotron nebula that revealed the presence of dust in the filament cores. We present the first spatially resolved emission spectra of dust in the Crab Nebula acquired with the Infrared Spectrograph aboard the Spitzer Space Telescope. The IR spectra are dominated by synchrotron emission and show forbidden line emission from from S, Si, Ne, Ar, O, Fe, and Ni. We derived a synchrotron spectral map from the 3.6 and 4.5 microns images, and subtracted this contribution from our data to produce a map of the residual continuum emission from dust. The dust emission appears to be concentrated along the ejecta filaments and is well described by an amorphous carbon or silicate grain compositions. We find a dust temperature of 55+/- 4 K for silicates and 60 +/- 7 K for carbon grains. The total estimated dust mass is 0.0012-0.012 solar masses, well below the theoretical dust yield predicted for a core-collapse supernova. Our grain heating model implies that the dust grain radii are relatively small, unlike what is expected for dust grains formed in a Type IIP SN.

preprint2012arXiv

The Extragalactic Background Light and the Gamma-ray Opacity of the Universe

The extragalactic background light (EBL) is one of the fundamental observational quantities in cosmology. All energy releases from resolved and unresolved extragalactic sources, and the light from any truly diffuse background, excluding the cosmic microwave background (CMB), contribute to its intensity and spectral energy distribution. It therefore plays a crucial role in cosmological tests for the formation and evolution of stellar objects and galaxies, and for setting limits on exotic energy releases in the universe. The EBL also plays an important role in the propagation of very high energy gamma-rays which are attenuated en route to Earth by pair producing gamma-gamma interactions with the EBL and CMB. The EBL affects the spectrum of the sources, predominantly blazars, in the ~10 GeV to 10 TeV energy regime. Knowledge of the EBL intensity and spectrum will allow the determination of the intrinsic blazar spectrum in a crucial energy regime that can be used to test particle acceleration mechanisms and VHE gamma-ray production models. Conversely, knowledge of the intrinsic gamma-ray spectrum and the detection of blazars at increasingly higher redshifts will set strong limits on the EBL and its evolution. This paper reviews the latest developments in the determination of the EBL and its impact on the current understanding of the origin and production mechanisms of gamma-rays in blazars, and on energy releases in the universe. The review concludes with a summary and future directions in Cherenkov Telescope Array techniques and in infrared ground-based and space observatories that will greatly improve our knowledge of the EBL and the origin and production of very high energy gamma-rays.

preprint2011arXiv

Infrared and X-Ray Spectroscopy of the Kes 75 Supernova Remnant Shell: Characterizing the Dust and Gas Properties

We present deep Chandra observations and Spitzer Space Telescope infrared (IR) spectroscopy of the shell in the composite supernova remnant (SNR) Kes 75 (G29.7-0.3). The remnant is composed of a central pulsar wind nebula and a bright partial shell in the south that is visible at radio, IR, and X-ray wavelengths. The X-ray emission can be modeled by either a single thermal component with a temperature of ~ 1.5 keV, or with two thermal components with temperatures of 1.5 and 0.2 keV. Previous studies suggest that the hot component may originate from reverse-shocked SN ejecta. However, our new analysis shows no definitive evidence for enhanced abundances of Si, S, Ar, Mg, and Fe, as expected from supernova (SN) ejecta, or for the IR spectral signatures characteristic of confirmed SN condensed dust, thus favoring a circumstellar or interstellar origin for the X-ray and IR emission. The X-ray and IR emission in the shell are spatially correlated, suggesting that the dust particles are collisionally heated by the X-ray emitting gas. The IR spectrum of the shell is dominated by continuum emission from dust with little, or no line emission. Modeling the IR spectrum shows that the dust is heated to a temperature of ~ 140 K by a relatively dense, hot plasma, that also gives rise to the hot X-ray emission component. The density inferred from the IR emission is significantly higher than the density inferred from the X-ray models, suggesting a low filling factor for this X-ray emitting gas. The total mass of the warm dust component is at least 0.013 solar masses, assuming no significant dust destruction has occurred in the shell. The IR data also reveal the presence of an additional plasma component with a cooler temperature, consistent with the 0.2 keV gas component. Our IR analysis therefore provides an independent verification of the cooler component of the X-ray emission.

preprint2011arXiv

Star and dust formation activities in AzTEC-3: A starburst galaxy at z = 5.3

Analyses of high-redshift ultraluminous infrared (IR) galaxies traditionally use the observed optical to submillimeter spectral energy distribution (SED) and estimates of the dynamical mass as observational constraints to derive the star formation rate (SFR), the stellar mass, and age of these objects. An important observational constraint neglected in the analysis is the mass of dust giving rise to the IR emission. In this paper we add this constraint to the analysis of AzTEC-3. Adopting an upper limit to the mass of stars and a bolometric luminosity for this object, we construct stellar and chemical evolutionary scenarios, constrained to produce the inferred dust mass and observed luminosity before the associated stellar mass exceeds the observational limit. We find that the model with a Top Heavy IMF provided the most plausible scenario consistent with the observational constraints. In this scenario the dust formed over a period of ~200 Myr, with a SFR of ~500 Msun/yr. These values for the age and SFR in AzTEC-3 are significantly higher and lower, respectively, from those derived without the dust mass constraint. However, this scenario is not unique, and others cannot be completely ruled out because of the prevailing uncertainties in the age of the galaxy, its bolometric luminosity, and its stellar and dust masses. A robust result of our models is that all scenarios require most of the radiating dust mass to have been accreted in molecular clouds. Our new procedure highlights the importance of a multiwavelength approach, and of the use of dust evolution models in constraining the age and the star formation activity and history in galaxies.

preprint2011arXiv

Strong New Constraints on the Extragalactic Background Light in the Near- to Mid-IR

Direct measurements of the extragalactic background light (EBL) in the near-IR to mid-IR waveband are extremely difficult due to an overwhelming foreground from the zodiacal light that outshines the faint cosmological diffuse radiation field by more than an order of magnitude. Indirect constraints on the EBL are provided by gamma-ray observations of AGN. Using the combination of the Fermi Gamma-Ray Space Telescope together with the current generation of ground-based air Cherenkov telescopes (H.E.S.S., MAGIC, and VERITAS) provides unprecedented sensitivity and spectral coverage for constraining the EBL in the near- to mid-IR. In this paper we present new limits on the EBL based on the analysis of the broad-band spectra of a select set of gamma-ray blazars covering 200 MeV to several TeV. The EBL intensity at 15 microns is constrained to be 1.36 +/- 0.58 nW m^-2 sr^-1. We find that the fast evolution and baseline EBL models of Stecker et al. (2006), as well as the model of Kneiske et al. (2004), predict significantly higher EBL intensities in the mid-IR (15 microns) than is allowed by the constraints derived here. In addition, the model of Franceschini et al. (2008) and the fiducial model of Dominguez et al. (2011) predict near- to mid-IR ratios smaller than that predicted by our analysis. Namely, their intensities in the near-IR are too low while their intensities in the mid-IR are marginally too high. All of the aforementioned models are inconsistent with our analysis at the >3 sigma level.

preprint2010arXiv

Disentangling the Origin and Heating Mechanism of Supernova Dust: Late-Time Spitzer Spectroscopy of the Type IIn SN 2005ip

This paper presents late-time near-infrared and {\it Spitzer} mid-infrared photometric and spectroscopic observations of warm dust in the Type IIn SN 2005ip in NGC 2906. The spectra show evidence for two dust components with different temperatures. Spanning the peak of the thermal emission, these observations provide strong constraints on the dust mass, temperature, and luminosity, which serve as critical diagnostics for disentangling the origin and heating mechanism of each component. The results suggest the warmer dust has a mass of $\sim 5 \times 10^{-4}~$\msolar, originates from newly formed dust in the ejecta, or possibly the cool, dense shell, and is continuously heated by the circumstellar interaction. By contrast, the cooler component likely originates from a circumstellar shock echo that forms from the heating of a large, pre-existing dust shell $\sim 0.01 - 0.05$~\msolar~by the late-time circumstellar interaction. The progenitor wind velocity derived from the blue edge of the He 1 1.083 \micron~P Cygni profile indicates a progenitor eruption likely formed this dust shell $\sim$100 years prior to the supernova explosion, which is consistent with a Luminous Blue Variable (LBV) progenitor star.

preprint2010arXiv

Five Years of Mid-Infrared Evolution of the Remnant of SN 1987A: The Encounter Between the Blast Wave and the Dusty Equatorial Ring

We have used the Spitzer satellite to monitor the mid-IR evolution of SN 1987A over a 5 year period spanning the epochs between days 6000 and 8000 since the explosion. The supernova (SN) has evolved into a supernova remnant (SNR) and its radiative output is dominated by the interaction of the SN blast wave with the pre-existing equatorial ring (ER). The mid-IR spectrum is dominated by emission from ~180 K silicate dust, collisionally-heated by the hot X-ray emitting gas with a temperature and density of ~5x10^6 K and 3x10^4 cm-3, respectively. The mass of the radiating dust is ~1.2x10^(-6) Msun on day 7554, and scales linearly with IR flux. The infrared to soft-X-ray flux ratio is roughly constant with a value of 2.5. Gas-grain collisions therefore dominate the cooling of the shocked gas. The constancy of of this ratio suggests that very little grain processing or gas cooling have occurred throughout this epoch. The shape of the dust spectrum remained unchanged during the observations while the total flux increased with a time dependence of t^(0.87), t being the time since the first encounter between the blast wave and the ER. These observations are consistent with the transitioning of the blast wave from free expansion to a Sedov phase as it propagates into the main body of the ER.

preprint2010arXiv

Observing Supernova 1987A with the Refurbished Hubble Space Telescope

Observations with the Hubble Space Telescope (HST), conducted since 1990, now offer an unprecedented glimpse into fast astrophysical shocks in the young remnant of supernova 1987A. Comparing observations taken in 2010 using the refurbished instruments on HST with data taken in 2004, just before the Space Telescope Imaging Spectrograph failed, we find that the Ly-a and H-a lines from shock emission continue to brighten, while their maximum velocities continue to decrease. We observe broad blueshifted Ly-a, which we attribute to resonant scattering of photons emitted from hotspots on the equatorial ring. We also detect NV~λλ1239,1243 A line emission, but only to the red of Ly-A. The profiles of the NV lines differ markedly from that of H-a, suggesting that the N^{4+} ions are scattered and accelerated by turbulent electromagnetic fields that isotropize the ions in the collisionless shock.

preprint2010arXiv

The chemistry of population III supernova ejecta: II - The nucleation of molecular clusters as a diagnostic for dust in the early universe

We study the formation of molecular precursors to dust in the ejecta of Population III supernovae using a chemical kinetic approach. Our work focuses on zero-metallicity 20 Msun and 170 Msun progenitors, and we consider fully-macroscopically mixed and unmixed ejecta. The nucleation stage for small silica, metal oxides and sulphides, pure metal, and carbon clusters is described with a new chemical reaction network. We consider the effect of the pressure dependence of critical nucleation rates, and test the impact of microscopically-mixed He+ on carbon dust formation. The unmixed ejecta of a 170 Msun progenitor supernova synthesizes ~ 5.6 Msun of small clusters, while its 20 Msun counterpart produces ~ 0.103 Msun. Our results point to smaller amounts of dust formed in the ejecta of Pop. III supernovae by a factor ~ 5 compared to values derived by previous studies, and to different dust chemical composition. Such deviations result from some erroneous assumptions made, the inappropriate use of classical nucleation theory to model dust formation, and the omission of the synthethis of molecules in supernova ejecta. Unmixed ejecta of massive Pop. III supernovae chiefly form silica and/or silicates, and pure silicon grains whereas their lower mass counterparts form a dust mixture dominated by silica and/or silicates, pure silicon and iron sulphides. Amorphous carbon can only condense in ejecta where the carbon-rich zone is deprived of He+. The first dust enrichment to the primordial gas in the early universe from Pop. III massive supernova comprises primarily pure silicon, silica and silicates. If carbon dust is present at redshift z> 6, alternative dust sources must be considered.

preprint2010arXiv

The Origin of Dust in the Early Universe: Probing the Star Formation History of Galaxies by their Dust Content

Two distinct scenarios for the origin of the ~ 4 \times 10^8 M\odot of dust observed in the high-redshift (z = 6.4) quasar J1148+5251 have been proposed. The first assumes that this galaxy is much younger than the age of the universe at that epoch so that only supernovae (SNe) could have produced this dust. The second scenario assumes a significantly older galactic age, so that the dust could have formed in lower-mass asymptotic giant branch (AGB) stars. Presenting new integral solutions for the chemical evolution of metals and dust in galaxies, we offer a critical evaluation of these two scenarios, and observational consequences that can discriminate between the two. We show that AGB stars can produce the inferred mass of dust in this object, however, the final mass of surviving dust depends on the galaxy's star formation history (SFH). In general supernovae cannot produce the observed amount of dust unless the average SN event creates over ~ 1 M\odot of dust in its ejecta. However, special SFHs can be constructed in which SNe can produce the inferred dust mass with a reasonable average dust yield of ~ 0.15 M\odot. The two scenarios propose different origins for the galaxy's spectral energy distribution, different star formation efficiencies and stellar masses, and consequently different comoving number densities of J1148+5251-type hyperluminous infrared (IR) objects. The detection of diagnostic mid-IR fine structure lines, and more complete surveys determining the comoving number density of these objects can discriminate between the two scenarios.

preprint2007arXiv

Evolution of Dust in Primordial Supernova Remnants: Can Dust Grains Formed in the Ejecta Survive and be Injected into the Early Interstellar Medium?

We investigate the evolution of dust that formed at Population III supernova (SN) explosions and its processing through the collisions with the reverse shocks resulting from the interaction of the SN ejecta with the ambient medium. In particular, we investigate the transport of the shocked dust within the SN remnant (SNR), and its effect on the chemical composition, the size distribution, and the total mass of dust surviving in SNR. We find that the evolution of the reverse shock, and hence its effect on the processing of the dust depends on the thickness of the envelope retained by the progenitor star. Furthermore, the transport and survival of the dust grains depend on their initial radius, a_{ini}, and composition: For Type II SNRs expanding into the interstellar medium (ISM) with a density of n_{H,0}=1 cm^{-3}, small grains with a_{ini} < ~ 0.05 micron are completely destroyed by sputtering in the postshock flow, while grains with a_{ini}= 0.05--0.2 micron are trapped into the dense shell behind the forward shock. Very large grains of a_{ini} > ~ 0.2 micron are ejected into the ISM without decreasing their sizes significantly. We find that the total mass fraction of dust that is destroyed by the reverse shock ranges from 0.2 to 1.0, depending on the energy of the explosion and the density of the ambient ISM. The results of our calculations have significant impact on the abundance pattern of subsequent generation of stars that form in the dense shell of primordial SNRs.

preprint2007arXiv

The origin and evolution of dust in high-redshift galaxies

Dusty hyperluminous galaxies in the early universe provide unique environments for studying the role of massive stars in the formation and destruction of dust. At redshifts above 6, when the universe was less than 1 Gyr old, dust could have only condensed in the explosive ejecta of Type II supernovae (SNe), since most of the progenitors of the AGB stars, the major alternative source of interstellar dust, did not have time to evolve off the main sequence. We present analytical models for the evolution of the gas, dust, and metals in high redshift galaxies, with a special application to SDSS J1148+5251, a hyperluminous quasar at z = 6.4. We show that an average SN must condense at least 1 Msun of dust to account for the mass of dust in this object, when grain destruction by supernova remnants (SNRs) is taken into account. This required yield is in excess of ~0.05 Msun, the largest mass of dust inferred from infrared observations of Cas A. If the yield of Cas A is typical, then other processes, such as accretion onto preexisting grains in molecular clouds is needed to produce the mass of dust in J1148+5251. For such process to be effective, SNR must significantly increase, presumably by non-evaporative grain-grain collisions during the late stages of their evolution, the number of nucleation centers onto which refractory elements can condense in molecular clouds.

preprint2003arXiv

New Interstellar Dust Models Consistent with Extinction, Emission, and Abundance Constraints

We present new interstellar dust models which have been derived by simultaneously fitting the far-ultraviolet to near-infrared extinction, the diffuse infrared (IR) emission and, unlike previous models, the elemental abundance constraints on the dust for different interstellar medium abundances, including solar, F and G star, and B star abundances. The fitting problem is a typical ill-posed inversion problem, in which the grain size distribution is the unknown, which we solve by using the method of regularization. The dust model contains various components: PAHs, bare silicate, graphite, and amorphous carbon particles, as well as composite particles containing silicate, organic refractory material, water ice, and voids. The optical properties of these components were calculated using physical optical constants. As a special case, we reproduce the Li & Draine (2001) results, however their model requires an excessive amount of silicon, magnesium, and iron to be locked up in dust: about 50 ppm (atoms per million of H atoms), significantly more than the upper limit imposed by solar abundances of these elements, about 34, 35, and 28 ppm, respectively. A major conclusion of this paper is that there is no unique interstellar dust model that simultaneously fits the observed extinction, diffuse IR emission, and abundances constraints.

preprint1998arXiv

The Submillimeter Frontier: A Space Science Imperative

A major goal of modern astrophysics is to understand the processes by which the universe evolved from its initial simplicity, as seen in measurements of the Cosmic Microwave Background, to the universe we see today, with complexity on all scales. While the diffuse background measurements of COBE reveal the importance of the far infrared and submillimeter in early galaxy and star formation, the understanding of the development of complex structure requires high resolution imaging and spectroscopy. We present a concept for a space mission called SPECS, the Submillimeter Probe of the Evolution of Cosmic Structure, which provides these capabilities. SPECS is a cold, spatial and spectral Michelson interferometer with adjustable baselines ranging up to 1 km. It has Hubble sensitivity and angular resolution in the far IR and submillimeter, spectral resolution up to 10,000, and a 14' field of view. SPECS will be able to image thermal dust continuum and infrared cooling and diagnostic line emission over a wide range of redshifts, providing extinction-free astrophysical probes of young galaxies and early cosmic structures, and measures of the luminosity and heavy element formation history of the universe. SPECS would also have the potential to improve vastly our knowledge of protostars, protoplanetary systems, Active Galactic Nuclei and other objects in the local universe. We recommend that a concerted effort be made during the next decade to develop and test certain critical technologies (photon-counting far IR detectors, formation flying spacecraft, cold, lightweight mirrors, and active coolers), so that SPECS can be deployed in about 2015.