Showing papers by "Karl D. Gordon published in 2017"

The SAGE-Spec Spitzer Legacy programme: the life-cycle of dust and gas in the Large Magellanic Cloud - Point source classification I

Abstract: We present the classification of 197 point sources observed with the Infrared Spectrograph in the SAGE-Spec Legacy programme on the Spitzer Space Telescope. We introduce a decision-tree method of object classification based on infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership and variability information, which is used to classify the SAGE-Spec sample of point sources. The decision tree has a broad application to mid-infrared spectroscopic surveys, where supporting photometry and variability information are available. We use these classifications to make deductions about the stellar populations of the Large Magellanic Cloud and the success of photometric classification methods. We find 90 asymptotic giant branch (AGB) stars, 29 young stellar objects, 23 post-AGB objects, 19 red supergiants, eight stellar photospheres, seven background galaxies, seven planetary nebulae, two H_(II) regions and 12 other objects, seven of which remain unclassified.

A New Stellar Atmosphere Grid and Comparisons with HST/STIS CALSPEC Flux Distributions

Abstract: The Space Telescope Imaging Spectrograph has measured the spectral energy distributions for several stars of types O, B, A, F, and G. These absolute fluxes from the CALSPEC database are fit with a new spectral grid computed from the ATLAS-APOGEE ATLAS9 model atmosphere database using a chi-square minimization technique in four parameters. The quality of the fits are compared for complete LTE grids by Castelli & Kurucz (CK04) and our new comprehensive LTE grid (BOSZ). For the cooler stars, the fits with the MARCS LTE grid are also evaluated, while the hottest stars are also fit with the NLTE Lanz & Hubeny OB star grids. Unfortunately, these NLTE models do not transition smoothly in the infrared to agree with our new BOSZ LTE grid at the NLTE lower limit of T eff = 15,000 K. The new BOSZ grid is available via the Space Telescope Institute MAST archive and has a much finer sampled IR wavelength scale than CK04, which will facilitate the modeling of stars observed by the James Webb Space Telescope. Our result for the angular diameter of Sirius agrees with the ground-based interferometric value.

Updated 34-band Photometry for the SINGS/KINGFISH Samples of Nearby Galaxies

Abstract: Alfred P. Sloan Foundation; Participating Institutions; NSF; U.S. Department of Energy, NASA; Japanese Monbukagakusho; Max Planck Society; Higher Education Funding Council for England; NASA

Updated 34-Band Photometry for the SINGS/KINGFISH Samples of Nearby Galaxies

Abstract: We present an update to the ultraviolet-to-radio database of global broadband photometry for the 79 nearby galaxies that comprise the union of the KINGFISH (Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel) and SINGS (Spitzer Infrared Nearby Galaxies Survey) samples. The 34-band dataset presented here includes contributions from observational work carried out with a variety of facilities including GALEX, SDSS, PS, NOAO, 2MASS, WISE, Spitzer, Herschel, Planck, JCMT, and the VLA. Improvements of note include recalibrations of previously-published SINGS BVRcIc and KINGFISH far-infrared/submillimeter photometry. Similar to previous results in the literature, an excess of submillimeter emission above model predictions is seen primarily for low-metallicity dwarf/irregular galaxies. This 34-band photometric dataset for the combined KINGFISH$+$SINGS sample serves as an important multi-wavelength reference for the variety of galaxies observed at low redshift. A thorough analysis of the observed spectral energy distributions is carried out in a companion paper.

Modeling dust emission in the Magellanic Clouds with Spitzer and Herschel

Abstract: Context. Dust modeling is crucial to infer dust properties and budget for galaxy studies. However, there are systematic disparities between dust grain models that result in corresponding systematic differences in the inferred dust properties of galaxies. Quantifying these systematics requires a consistent fitting analysis.Aims. We compare the output dust parameters and assess the differences between two dust grain models, the DustEM model and THEMIS. In this study, we use a single fitting method applied to all the models to extract a coherent and unique statistical analysis.Methods. We fit the models to the dust emission seen by Spitzer and Herschel in the Small and Large Magellanic Clouds (SMC and LMC). The observations cover the infrared (IR) spectrum from a few microns to the sub-millimeter range. For each fitted pixel, we calculate the full n-D likelihood based on a previously described method. The free parameters are both environmental (U , the interstellar radiation field strength; α ISRF , power-law coefficient for a multi-U environment; Ω∗ , the starlight strength) and intrinsic to the model (Y i : abundances of the grain species i ; α sCM20 , coefficient in the small carbon grain size distribution).Results. Fractional residuals of five different sets of parameters show that fitting THEMIS brings a more accurate reproduction of the observations than the DustEM model. However, independent variations of the dust species show strong model-dependencies. We find that the abundance of silicates can only be constrained to an upper-limit and that the silicate/carbon ratio is different than that seen in our Galaxy. In the LMC, our fits result in dust masses slightly lower than those found in the literature, by a factor lower than 2. In the SMC, we find dust masses in agreement with previous studies.

I. A 3D externally illuminated slab benchmark for dust radiative transfer

Abstract: Context. The radiative transport of photons through arbitrary three-dimensional (3D) structures of dust is a challenging problem due to the anisotropic scattering of dust grains and strong coupling between different spatial regions. The radiative transfer problem in 3D is solved using Monte Carlo or Ray Tracing techniques as no full analytic solution exists for the true 3D structures. Aims. We provide the first 3D dust radiative transfer benchmark composed of a slab of dust with uniform density externally illuminated by a star. This simple 3D benchmark is explicitly formulated to provide tests of the different components of the radiative transfer problem including dust absorption, scattering, and emission. Methods. The details of the external star, the slab itself, and the dust properties are provided. This benchmark includes models with a range of dust optical depths fully probing cases that are optically thin at all wavelengths to optically thick at most wavelengths. The dust properties adopted are characteristic of the diffuse Milky Way interstellar medium. This benchmark includes solutions for the full dust emission including single photon (stochastic) heating as well as two simplifying approximations: One where all grains are considered in equilibrium with the radiation field and one where the emission is from a single effective grain with size-distribution-averaged properties. A total of six Monte Carlo codes and one Ray Tracing code provide solutions to this benchmark. Results. The solution to this benchmark is given as global spectral energy distributions (SEDs) and images at select diagnostic wavelengths from the ultraviolet through the infrared. Comparison of the results revealed that the global SEDs are consistent on average to a few percent for all but the scattered stellar flux at very high optical depths. The image results are consistent within 10%, again except for the stellar scattered flux at very high optical depths. The lack of agreement between different codes of the scattered flux at high optical depths is quantified for the first time. Convergence tests using one of the Monte Carlo codes illustrate the sensitivity of the solutions to various model parameters. Conclusions. We provide the first 3D dust radiative transfer benchmark and validate the accuracy of this benchmark through comparisons between multiple independent codes and detailed convergence tests.

Dust Abundance Variations in the Magellanic Clouds: Probing the Life-cycle of Metals with All-sky Surveys

Abstract: Observations and modeling suggest that the dust abundance (gas-to-dust ratio, G/D) depends on (surface) density. The variations of the G/D provide constraints on the timescales for the different processes involved in the lifecycle of metals in galaxies. Recent G/D measurements based on Herschel data suggest a factor 5---10 decrease in the dust abundance between the dense and diffuse interstellar medium (ISM) in the Magellanic Clouds. However, the relative nature of the Herschel measurements precludes definitive conclusions on the magnitude of those variations. We investigate the variations of the dust abundance in the LMC and SMC using all-sky far-infrared surveys, which do not suffer from the limitations of Herschel on their zero-point calibration. We stack the dust spectral energy distribution (SED) at 100, 350, 550, and 850 microns from IRAS and Planck in intervals of gas surface density, model the stacked SEDs to derive the dust surface density, and constrain the relation between G/D and gas surface density in the range 10---100 \\Msu pc$^{-2}$ on $\\sim$ 80 pc scales. We find that G/D decreases by factors of 3 (from 1500 to 500) in the LMC and 7 (from 1.5$\\times 10^4$ to 2000) in the SMC between the diffuse and dense ISM. The surface density dependence of G/D is consistent with elemental depletions and with simple modeling of the accretion of gas-phase metals onto dust grains. This result has important implications for the sub-grid modeling of galaxy evolution, and for the calibration of dust-based gas mass estimates, both locally and at high-redshift.

A New Stellar Atmosphere Grid and Comparisons with HST/STIS Calspec Flux Distributions

Abstract: The Space Telescope Imaging Spectrograph (STIS) has measured the spectral energy distributions (SEDs) for several stars of types O, B, A, F, and G. These absolute fluxes from the CALSPEC database are fit with a new spectral grid computed from the ATLAS-APOGEE ATLAS9 model atmosphere database using a chi-square minimization technique in four parameters. The quality of the fits are compared for complete LTE grids by Castelli & Kurucz (CK04) and our new comprehensive LTE grid (BOSZ). For the cooler stars, the fits with the MARCS LTE grid are also evaluated, while the hottest stars are also fit with the NLTE Lanz & Hubeny OB star grids. Unfortunately, these NLTE models do not transition smoothly in the infrared to agree with our new BOSZ LTE grid at the NLTE lower limit of Teff =15,000K. The new BOSZ grid is available via the Space Telescope Institute MAST archive and has a much finer sampled IR wavelength scale than CK04, which will facilitate the modeling of stars observed by the James Webb Space Telescope (JWST). Our result for the angular diameter of Sirius agrees with the ground- based interferometric value.

Dust Abundance Variations in the Magellanic Clouds: Probing the Lifecycle of Metals with All-Sky Surveys

Abstract: Observations and modeling suggest that the dust abundance (gas-to-dust ratio, G/D) depends on (surface) density. The variations of the G/D provide constraints on the timescales for the different processes involved in the lifecycle of metals in galaxies. Recent G/D measurements based on Herschel data suggest a factor 5---10 decrease in the dust abundance between the dense and diffuse interstellar medium (ISM) in the Magellanic Clouds. However, the relative nature of the Herschel measurements precludes definitive conclusions on the magnitude of those variations. We investigate the variations of the dust abundance in the LMC and SMC using all-sky far-infrared surveys, which do not suffer from the limitations of Herschel on their zero-point calibration. We stack the dust spectral energy distribution (SED) at 100, 350, 550, and 850 microns from IRAS and Planck in intervals of gas surface density, model the stacked SEDs to derive the dust surface density, and constrain the relation between G/D and gas surface density in the range 10---100 \Msu pc$^{-2}$ on $\sim$ 80 pc scales. We find that G/D decreases by factors of 3 (from 1500 to 500) in the LMC and 7 (from 1.5$\times 10^4$ to 2000) in the SMC between the diffuse and dense ISM. The surface density dependence of G/D is consistent with elemental depletions and with simple modeling of the accretion of gas-phase metals onto dust grains. This result has important implications for the sub-grid modeling of galaxy evolution, and for the calibration of dust-based gas mass estimates, both locally and at high-redshift.

TRUST I: A 3D externally illuminated slab benchmark for dust radiative transfer

Abstract: The radiative transport of photons through arbitrary three-dimensional (3D) structures of dust is a challenging problem due to the anisotropic scattering of dust grains and strong coupling between different spatial regions. The radiative transfer problem in 3D is solved using Monte Carlo or Ray Tracing techniques as no full analytic solution exists for the true 3D structures. We provide the first 3D dust radiative transfer benchmark composed of a slab of dust with uniform density externally illuminated by a star. This simple 3D benchmark is explicitly formulated to provide tests of the different components of the radiative transfer problem including dust absorption, scattering, and emission. This benchmark includes models with a range of dust optical depths fully probing cases that are optically thin at all wavelengths to optically thick at most wavelengths. This benchmark includes solutions for the full dust emission including single photon (stochastic) heating as well as two simplifying approximations: One where all grains are considered in equilibrium with the radiation field and one where the emission is from a single effective grain with size-distribution-averaged properties. A total of six Monte Carlo codes and one Ray Tracing code provide solutions to this benchmark. Comparison of the results revealed that the global SEDs are consistent on average to a few percent for all but the scattered stellar flux at very high optical depths. The image results are consistent within 10%, again except for the stellar scattered flux at very high optical depths. The lack of agreement between different codes of the scattered flux at high optical depths is quantified for the first time. We provide the first 3D dust radiative transfer benchmark and validate the accuracy of this benchmark through comparisons between multiple independent codes and detailed convergence tests.

The Small Magellanic Cloud Investigation of Dust and Gas Evolution (SMIDGE): The Dust Extinction Curve from Red Clump Stars

Abstract: We use Hubble Space Telescope (HST) observations of red clump stars taken as part of the Small Magellanic Cloud Investigation of Dust and Gas Evolution (SMIDGE) program to measure the average dust extinction curve in a ~ 200 pc x 100 pc region in the southwest bar of the Small Magellanic Cloud (SMC). The rich information provided by our 8-band ultra-violet through near-infrared photometry allows us to model the color-magnitude diagram of the red clump accounting for the extinction curve shape, a log-normal distribution of $A_{V}$, and the depth of the stellar distribution along the line of sight. We measure an extinction curve with $R_{475} = A_{475}/(A_{475}-A_{814})$ = 2.65 $\pm$ 0.11. This measurement is significantly larger than the equivalent values of published Milky Way $R_{V}$ = 3.1 ($R_{475} = 1.83$) and SMC Bar $R_{V}$ = 2.74 ($R_{475} = 1.86$) extinction curves. Similar extinction curve offsets in the Large Magellanic Cloud (LMC) have been interpreted as the effect of large dust grains. We demonstrate that the line-of-sight depth of the SMC (and LMC) introduces an apparent "gray" contribution to the extinction curve inferred from the morphology of the red clump. We show that no gray dust component is needed to explain extinction curve measurements when a full-width half-max depth of 10 $\pm$ 2 kpc in the stellar distribution of the SMC (5 $\pm$ 1 kpc for the LMC) is considered, which agrees with recent studies of Magellanic Cloud stellar structure. The results of our work demonstrate the power of broad-band HST imaging for simultaneously constraining dust and galactic structure outside the Milky Way.

Erratum: “Dust and Gas in the Magellanic Clouds from the HERITAGE Herschel Key Project. I. Dust Properties and Insights into the Origin of the Submm Excess Emission” ( 2014, ApJ, 797, 85 )

Abstract: Marc Sauvage, Jonathan Seale, Marta Sewiło, Kirill Tchernyshyov, and Ramin Skibba 1 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 2 Sterrenkundig Observatorium, Universiteit Gent, Gent, Belgium 3 Observatoire astronomique de Strasbourg, Université de Strasbourg, CNRS, UMR 7550, 11 rue de lUniversité, F-67000 Strasbourg, France 4 Department of Astronomy, 475 North Charter Street, University of Wisconsin, Madison, WI 53706, USA 5 CESR, Université de Toulouse, UPS, 9 Avenue du Colonel Roche, F-31028 Toulouse, Cedex 4, France 6 Université de Toulouse, UPS-OMP, IRAP, 31028 Toulouse Cedex 4, France 7 Department of Astronomy, Lab for Millimeter-wave Astronomy, University of Maryland, College Park, MD 20742, USA 8 Observational Cosmology Lab, Code 665, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 9 Oak Ridge Associated Universities (ORAU), Oak Ridge, TN 37831, USA 10 Louisiana State University, Department of Physics & Astronomy, 233-A Nicholson Hall, Tower Dr., Baton Rouge, LA 70803, USA 11 Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA 12 Raytheon Company, 1151 East Hermans Road, Tucson, AZ 85756, USA 13 Department of Physics, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan 14 European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching-bei-Mnchen, Germany 15 CEA, Laboratoire AIM, Irfu/SAp, Orme des Merisiers, F-91191 Gif-sur-Yvette, France 16 Max-Planck-Institut für Astronomie, Königstuhl 17, D-69117 Heidelberg, Germany 17 Department of Astronomy, University of Virginia, and National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903, USA 18 Sterrewacht Leiden, Leiden University, P.O. Box 9513, NL-2300 RA Leiden, The Netherlands 19 National Astronomical Observatory of Japan, Osawa, Mitaka, Tokyo, 181-8588, Japan 20 314 Physics Building, Department of Physics and Astronomy, University of Missouri, Columbia, MO 65211, USA 21 Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK 22 Department of Astrophysics, Graduate School of Science, Osaka Prefecture University, Sakai, Osaka 599-8531, Japan 23 CNRS, Observatoire de Paris—Lab. GEPI, Bat. 11, 5, place Jules Janssen, 92195 Meudon CEDEX, France 24 Departamento de Astronomía, Universidad de Chile, Casilla 36-D, Santiago, Chile 25 The Johns Hopkins University, Department of Physics and Astronomy, 366 Bloomberg Center, 3400 N. Charles Street, Baltimore, MD 21218, USA 26 Center for Astrophysics and Space Sciences, Department of Physics, University of California, 9500 Gilman Dr, La Jolla, San Diego, CA 92093, USA Received 2017 February 10; published 2017 March 7

Dust Emission at 8 and 24 μm as Diagnostics of H II Region Radiative Transfer

Abstract: National Science Foundation AST-1210285 CONICYT (Chile) through FONDECYT 1140839 project BASAL PFB-06

Radiative Feedback from Massive Stars as Traced by Multiband Imaging and Spectroscopic Mosaics

Abstract: Massive stars disrupt their natal molecular cloud material by dissociating molecules, ionizing atoms and molecules, and heating the gas and dust. These processes drive the evolution of interstellar matter in our Galaxy and throughout the Universe from the era of vigorous star formation at redshifts of 1-3, to the present day. Much of this interaction occurs in Photo-Dissociation Regions (PDRs) where far-ultraviolet photons of these stars create a largely neutral, but warm region of gas and dust. PDR emission dominates the IR spectra of star-forming galaxies and also provides a unique tool to study in detail the physical and chemical processes that are relevant for most of the mass in inter- and circumstellar media including diffuse clouds, protoplanetary disk- and molecular cloud surfaces, globules, planetary nebulae, and starburst galaxies. We propose to provide template datasets designed to identify key PDR characteristics in JWST spectra in order to guide the preparation of Cycle 2 proposals on star-forming regions in our Galaxy and beyond. We plan to obtain the first spatially resolved, high spectral resolution IR observations of a PDR using NIRCam, NIRSpec, and MIRI. These data will test widely used theoretical models and extend them into the JWST era. We have engaged the broader community as exemplified by the supporting large international team of 138 scientists. We will assist the community interested in JWST observations of PDRs through science-enabling products that will guide observational planning and allow fast data analysis. We will train the community through telecons and dedicated workshops.

The Small Magellanic Cloud Investigation of Dust and Gas Evolution (SMIDGE): The Dust Extinction Curve from Red Clump Stars

Abstract: We use Hubble Space Telescope (HST) observations of red clump stars taken as part of the Small Magellanic Cloud Investigation of Dust and Gas Evolution (SMIDGE) program to measure the average dust extinction curve in a ~ 200 pc x 100 pc region in the southwest bar of the Small Magellanic Cloud (SMC). The rich information provided by our 8-band ultra-violet through near-infrared photometry allows us to model the color-magnitude diagram of the red clump accounting for the extinction curve shape, a log-normal distribution of $A_{V}$, and the depth of the stellar distribution along the line of sight. We measure an extinction curve with $R_{475} = A_{475}/(A_{475}-A_{814})$ = 2.65 $\pm$ 0.11. This measurement is significantly larger than the equivalent values of published Milky Way $R_{V}$ = 3.1 ($R_{475} = 1.83$) and SMC Bar $R_{V}$ = 2.74 ($R_{475} = 1.86$) extinction curves. Similar extinction curve offsets in the Large Magellanic Cloud (LMC) have been interpreted as the effect of large dust grains. We demonstrate that the line-of-sight depth of the SMC (and LMC) introduces an apparent "gray" contribution to the extinction curve inferred from the morphology of the red clump. We show that no gray dust component is needed to explain extinction curve measurements when a full-width half-max depth of 10 $\pm$ 2 kpc in the stellar distribution of the SMC (5 $\pm$ 1 kpc for the LMC) is considered, which agrees with recent studies of Magellanic Cloud stellar structure. The results of our work demonstrate the power of broad-band HST imaging for simultaneously constraining dust and galactic structure outside the Milky Way.

Dust emission at 8-mic and 24-mic as Diagnostics of HII Region Radiative Transfer

Abstract: We use the Spitzer SAGE survey of the Magellanic Clouds to evaluate the relationship between the 8-mic PAH emission, 24-mic hot dust emission, and HII region radiative transfer. We confirm that in the higher-metallicity Large Magellanic Cloud, PAH destruction is sensitive to optically thin conditions in the nebular Lyman continuum: objects identified as optically thin candidates based on nebular ionization structure show 6 times lower median 8-mic surface brightness (0.18 mJy arcsec^-2) than their optically thick counterparts (1.2 mJy arcsec^-2). The 24-mic surface brightness also shows a factor of 3 offset between the two classes of objects (0.13 vs 0.44 mJy arcsec^-2, respectively), which is driven by the association between the very small dust grains and higher density gas found at higher nebular optical depths. In contrast, PAH and dust formation in the low-metallicity Small Magellanic Cloud is strongly inhibited such that we find no variation in either 8-mic or 24-mic emission between our optically thick and thin samples. This is attributable to extremely low PAH and dust production together with high, corrosive UV photon fluxes in this low-metallicity environment. The dust mass surface densities and gas-to-dust ratios determined from dust maps using Herschel HERITAGE survey data support this interpretation.