Source author record

Karl D. Gordon

Karl D. Gordon 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

51works
7topics
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

51 published item(s)

preprint2024arXiv

Formation of the Methyl Cation by Photochemistry in a Protoplanetary Disk

Forty years ago it was proposed that gas phase organic chemistry in the interstellar medium was initiated by the methyl cation CH3+, but hitherto it has not been observed outside the Solar System. Alternative routes involving processes on grain surfaces have been invoked. Here we report JWST observations of CH3+ in a protoplanetary disk in the Orion star forming region. We find that gas-phase organic chemistry is activated by UV irradiation.

preprint2022arXiv

How do Spitzer IRAC Fluxes Compare to HST CALSPEC

An accurate tabulation of stellar brightness in physical units is essential for a multitude of scientific endeavors. The HST/CALSPEC database of flux standards contains many stars with spectral coverage in the 0.115--1 \micron\ range with some extensions to longer wavelengths of 1.7 or 2.5 \micron. Modeled flux distributions to 32 \micron\ for calibration of JWST complement the shorter wavelength HST measurements. Understanding the differences between IRAC observations and CALSPEC models is important for science that uses IR fluxes from multiple instruments, including JWST. The absolute flux of Spitzer IRAC photometry at 3.6--8 \micron\ agrees with CALSPEC synthetic photometry to 1\% for the three prime HST standards G191B2B, GD153, and GD71. For a set of 17--22 A-star standards, the average IRAC difference rises from agreement at 3.6 \micron\ to 3.4 $\pm$0.1\% brighter than CALSPEC at 8 \micron. For a smaller set of G-stars, the average of the IRAC photometry falls below CALSPEC by as much as 3.7 $\pm$0.3\% for IRAC1, while one G-star, P330E, is consistent with the A-star ensemble of IRAC/CALSPEC ratios.

preprint2022arXiv

METAL: The Metal Evolution, Transport, and Abundance in the Large Magellanic Cloud Hubble program. IV. Calibration of Dust Depletions vs Abundance Ratios in the Milky Way and Magellanic Clouds and Application to Damped Lyman-alpha Systems

The evolution of the metal content of the universe can be tracked through rest-frame UV spectroscopy of damped Ly-$α$ systems (DLAs). Gas-phase abundances in DLAs must be corrected for dust depletion effects, which can be accomplished by calibrating the relation between abundance ratios such as [Zn/Fe] and depletions (the fraction of metals in gas, as opposed to dust). Using samples of gas-phase abundances and depletions in the Milky Way (MW), LMC, and SMC, we demonstrate that the relation between [Zn/Fe] and other abundance ratios does not change significantly between these local galaxies and DLAs, indicating that [Zn/Fe] should trace depletions of heavy elements in those systems. The availability of photospheric abundances in young massive stars, a proxy for the total (gas+dust) metallicity of neutral gas, in the MW LMC, and SMC allows us to calibrate the relation between [Zn/Fe] and depletions in these nearby galaxies. We apply the local calibrations of depletions to DLA systems. We find that the fraction of metals in dust, the dust-to-gas-ratio, and total abundances are 2-5 times lower than inferred from previous depletion calibrations based on MW measurements and a different formalism. However, the trend of dust abundance vs. metallicity remains only slightly sub-linear for all existing depletion calibrations, contrary to what is inferred from FIR, 21 cm, and CO emission in nearby galaxies and predicted by chemical evolution models. Observational constraints on the FIR dust opacity and depletions at metallicities lower than 20\% solar will be needed to resolve this tension.

preprint2022arXiv

PDRs4All: A JWST Early Release Science Program on radiative feedback from massive stars

Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on 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. The dominant feedback processes can be probed by observations of the Photo-Dissociation Regions (PDRs) where the far-ultraviolet photons of massive stars create warm regions of gas and dust in the neutral atomic and molecular gas. PDR emission 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, proto-planetary disks and molecular cloud surfaces, globules, planetary nebulae, and star-forming regions. PDR emission dominates the infrared (IR) spectra of star-forming galaxies. Most of the Galactic and extragalactic observations obtained with the James Webb Space Telescope (JWST) will therefore arise in PDR emission. In this paper we present an Early Release Science program using the MIRI, NIRSpec, and NIRCam instruments dedicated to the observations of an emblematic and nearby PDR: the Orion Bar. These early JWST observations will provide template datasets designed to identify key PDR characteristics in JWST observations. These data will serve to benchmark PDR models and extend them into the JWST era. We also present the Science-Enabling products that we will provide to the community. These template datasets and Science-Enabling products will guide the preparation of future proposals on star-forming regions in our Galaxy and beyond and will facilitate data analysis and interpretation of forthcoming JWST observations.

preprint2022arXiv

Searching for TESS Photometric Variability of Possible JWST Spectrophotometric Standard Stars

We use data from the Transiting Exoplanet Survey Satellite (TESS) to search for, and set limits on, optical to near-infrared photometric variability of the well-vetted, candidate James Webb Space Telescope (JWST) spectrophotometric standards. Our search of 37 of these candidate standards has revealed measurable periodic variability in 15 stars. The majority of those show variability that is less than half a percent; however, four stars are observed to vary photometrically, from minimum to maximum flux, by more than 1% (the G dwarf HD 38949 and three fainter A dwarfs). Variability of this size would likely impact the error budget in the spectrophotometric calibration of the science instruments aboard JWST. For the 22 candidate standards with no detected variability, we report upper limits on the observed changes in flux. Despite some systematic noise, all stars brighter than 12 magnitude in the TESS band show a 3 sigma upper limit on the total change in brightness of less than half a percent on time scales between an hour and multiple weeks, empirically establishing their suitability as spectrophotometric standards. We further discuss the value and limits of high-cadence, high-precision photometric monitoring with TESS as a tool to vet the suitability of stars to act as spectrophotometric standards.

preprint2022arXiv

SpeX near-infrared spectroscopic extinction curves in the Milky Way

Interstellar dust extinction curves provide valuable information about dust properties, including the composition and size of the dust grains, and are essential to correct observations for the effects of interstellar dust. In this work, we measure a representative sample of near-infrared (NIR; 0.8-5.5 $μ$m) spectroscopic extinction curves for the first time, enabling us to investigate the extinction at wavelengths where it is usually only measured in broad photometric bands. We use IRTF/SpeX spectra of a sample of reddened and comparison stars to measure 15 extinction curves with the pair method. Our sample spans A(V) values from 0.78 to 5.65 and R(V) values from 2.43 to 5.33. We confirm that the NIR extinction curves are well fit by a power law, with indices and amplitudes differing from sight line to sight line. Our average diffuse NIR extinction curve can be represented by a single power law with index $α= 1.7$, but because of the sight line-to-sight line variations, the shape of any average curve will depend on the parental sample. We find that most of the variation in our sample can be linked to the ratio of total-to-selective extinction R(V), a rough measurement of the average dust grain size. Two sight lines in our sample clearly show the ice extinction feature at 3 $μ$m, which can be fitted by a modified Drude profile. We find tentative ice detections with slightly over 3$σ$ significance in two other sight lines. In our average diffuse extinction curve, we measure a 3$σ$ upper limit of A(ice)/A(V) = 0.0021 for this ice feature.

preprint2022arXiv

The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package

The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates for the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project.

preprint2022arXiv

The James Webb Space Telescope Absolute Flux Calibration. I. Program Design and Calibrator Stars

It is critical for James Webb Space Telescope (JWST) science that instrumental units are converted to physical units. We detail the design of the JWST absolute flux calibration program that has the core goal of ensuring a robust flux calibration internal to and between all the science instruments for both point and extended source science. This program will observe a sample of calibration stars that have been extensively vetted based mainly on Hubble Space Telescope, Spitzer Space Telescope, and Transiting Exoplanet Survey Satellite observations. The program uses multiple stars of three different, well understood types (hot stars, A dwarfs, and solar analogs) to allow for the statistical (within a type) and systematic (between types) uncertainties to be quantified. The program explicitly includes observations to calibrate every instrument mode, further vet the set of calibration stars, measure the instrumental repeatability, measure the relative calibration between subarrays and full frame, and check the relative calibration between faint and bright stars. For photometry, we have set up our calibration to directly support both the convention based on the band average flux density and the convention based on the flux density at a fixed wavelength.

preprint2022arXiv

The JWST Early Release Observations

The James Webb Space Telescope (JWST) Early Release Observations (EROs) is a set of public outreach products created to mark the end of commissioning and the beginning of science operations for JWST. Colloquially known as the "Webb First Images and Spectra", these products were intended to demonstrate to the worldwide public that JWST is ready for science, and is capable of producing spectacular results. The package was released on July 12, 2022, and included images and spectra of the galaxy cluster SMACS~J0723.3-7327 and distant lensed galaxies, the interacting galaxy group Stephan's Quintet, NGC 3324 in the Carina star-forming complex, the Southern Ring planetary nebula NGC 3132, and the transiting hot Jupiter WASP 96b. This paper describes the ERO technical design, observations, and scientific processing of data underlying the colorful outreach products.

preprint2021arXiv

Relations between mid-IR dust emission and UV extinction

We analyze low resolution Spitzer infrared (IR) 5-14 micron spectra of the diffuse emission toward a carefully selected sample of stars. The sample is composed of sight lines toward stars that have well determined ultraviolet (UV) extinction curves and which are shown to lie beyond effectively all of the extinguishing and emitting dust along their lines of sight. Our sample includes sight lines whose UV curve extinction curves exhibit a wide range of curve morphology and which sample a variety of interstellar environments. As a result, this unique sample enabled us to study the connection between the extinction and emission properties of the same grains, and to examine their response to different physical environments. We quantify the emission features in terms of the PAH model given by Draine & Li (2007) and a set of additional features, not known to be related to PAH emission. We compare the intensities of the different features in the Spitzer mid-IR spectra with the Fitzpatrick & Massa (2007) parameters which describe the shapes of UV to near-IR extinction curves. Our primary result is that there is a strong correlation between the area of the 2175 A UV bump in the extinction curves of the program stars and the strengths of the major PAH emission features in the mid-IR spectra for the same lines of sight.

preprint2019arXiv

An Analysis of the Shapes of Interstellar Extinction Curves. VII. Milky Way Spectrophotometric Optical-through-Ultraviolet Extinction and Its R-Dependence

We produce a set of 72 NIR through UV extinction curves by combining new HST/STIS optical spectrophotometry with existing IUE spectrophotometry (yielding gapless coverage from 1150 to 10000 Angstroms) and NIR photometry. These curves are used to determine a new, internally consistent, NIR through UV Milky Way Mean Curve and to characterize how the shapes of the extinction curves depend on R(V). We emphasize that while this dependence captures much of the curve variability, there remains considerable variation which is independent of R(V). We use the optical spectrophotometry to verify the presence of structure at intermediate wavelength scales in the curves. The fact that the optical through UV portions of the curves are sampled at relatively high resolution makes them very useful for determining how extinction affects different broad band systems, and we provide several examples. Finally, we compare our results to previous investigations.

preprint2019arXiv

Evidence of dust grain evolution from extinction mapping in the IC 63 photodissociation region

Photodissociation regions (PDRs) are parts of the ISM consisting of predominantly neutral gas, located at the interface between H II regions and molecular clouds. The physical conditions within these regions show variations on very short spatial scales, and therefore PDRs constitute ideal laboratories for investigating the properties and evolution of dust grains. We have mapped IC 63 at high resolution from the UV to the NIR (275 nm to 1.6 $μ$m), using the Hubble Space Telescope WFC3. Using a Bayesian SED fitting tool, we simultaneously derive a set of stellar ($T_\text{eff}$, $\log(g)$, distance) and extinction ($A_V$, $R_V$) parameters for 520 background stars. We present maps of $A_V$ and $R_V$ with a resolution of 25 arcsec based on these results. The extinction properties vary across the PDR, with values for $A_V$ between 0.5 and 1.4 mag, and a decreasing trend in $R_V$, going from 3.7 at the front of the nebula to values as low as 2.5 further in. This provides evidence for evolution of the dust optical properties. We fit two modified blackbodies to the MIR and FIR SED, obtained by combining the $A_V$ map with data from Spitzer and Herschel. We derive effective temperatures (30 K and 227 K) and the ratio of opacities at 160 $μ$m to V band $κ_{160} / κ_V$ ($7.0 \times 10^{-4}$ and $2.9 \times 10^{-9}$) for the two dust populations. Similar fits to individual pixels show spatial variations of $κ_{160} / κ_{V}$. The analysis of our HST data, combined with these Spitzer and Herschel data, provides the first panchromatic view of dust within a PDR.

preprint2016arXiv

Composite biasing in Monte Carlo radiative transfer

Biasing or importance sampling is a powerful technique in Monte Carlo radiative transfer, and can be applied in different forms to increase the accuracy and efficiency of simulations. One of the drawbacks of the use of biasing is the potential introduction of large weight factors. We discuss a general strategy, composite biasing, to suppress the appearance of large weight factors. We use this composite biasing approach for two different problems faced by current state-of-the-art Monte Carlo radiative transfer codes: the generation of photon packages from multiple components, and the penetration of radiation through high optical depth barriers. In both cases, the implementation of the relevant algorithms is trivial and does not interfere with any other optimisation techniques. Through simple test models, we demonstrate the general applicability, accuracy and efficiency of the composite biasing approach. In particular, for the penetration of high optical depths, the gain in efficiency is spectacular for the specific problems that we consider: in simulations with composite path length stretching, high accuracy results are obtained even for simulations with modest numbers of photon packages, while simulations without biasing cannot reach convergence, even with a huge number of photon packages.

preprint2016arXiv

The Panchromatic Hubble Andromeda Treasury XV. The BEAST: Bayesian Extinction and Stellar Tool

We present the Bayesian Extinction And Stellar Tool (BEAST), a probabilistic approach to modeling the dust extinguished photometric spectral energy distribution of an individual star while accounting for observational uncertainties common to large resolved star surveys. Given a set of photometric measurements and an observational uncertainty model, the BEAST infers the physical properties of the stellar source using stellar evolution and atmosphere models and constrains the line of sight extinction using a newly developed mixture model that encompasses the full range of dust extinction curves seen in the Local Group. The BEAST is specifically formulated for use with large multi-band surveys of resolved stellar populations. Our approach accounts for measurement uncertainties and any covariance between them due to stellar crowding (both systematic biases and uncertainties in the bias) and absolute flux calibration, thereby incorporating the full information content of the measurement. We illustrate the accuracy and precision possible with the BEAST using data from the Panchromatic Hubble Andromeda Treasury. While the BEAST has been developed for this survey, it can be easily applied to similar existing and planned resolved star surveys.

preprint2015arXiv

New Ultraviolet Extinction Curves for Interstellar Dust in M31

New low-resolution UV spectra of a sample of reddened OB stars in M31 were obtained with HST/STIS to study the wavelength dependence of interstellar extinction and the nature of the underlying dust grain populations. Extinction curves were constructed for four reddened sightlines in M31 paired with closely matching stellar atmosphere models. The new curves have a much higher S/N than previous studies. Direct measurements of N(H I) were made using the Ly$α$ absorption lines enabling gas-to-dust ratios to be calculated. The sightlines have a range in galactocentric distance of 5 to 14 kpc and represent dust from regions of different metallicities and gas-to-dust ratios. The metallicities sampled range from Solar to 1.5 Solar. The measured curves show similarity to those seen in the Milky Way and the Large Magellanic Cloud. The Maximum Entropy Method was used to investigate the dust composition and size distribution for the sightlines observed in this program finding that the extinction curves can be produced with the available carbon and silicon abundances if the metallicity is super-Solar.

preprint2015arXiv

The Mid-Infrared Instrument for the James Webb Space Telescope, X. Operations and Data Reduction

We describe the operations concept and data reduction plan for the Mid- Infrared Instrument (MIRI) for the James Webb Space Telescope (JWST). The overall JWST operations concept is to use Observation Templates (OTs) to provide a straightforward and intuitive way for users to specify observations. MIRI has four OTs that correspond to the four observing modes: 1.) Imaging, 2.) Coronagraphy, 3.) Low Resolution Spectroscopy, and 4.) Medium Resolution Spectroscopy. We outline the user choices and expansion of these choices into detailed instrument operations. The data reduction plans for MIRI are split into three stages, where the specificity of the reduction steps to the observation type increases with stage. The reduction starts with integration ramps: stage 1 yields uncalibrated slope images; stage 2 calibrates the slope images; and then stage 3 combines multiple calibrated slope images into high level data products (e.g. mosaics, spectral cubes, and extracted source information). Finally, we give examples of the data and data products that will be derived from each of the four different OTs.

preprint2015arXiv

The Panchromatic Hubble Andromeda Treasury VIII: A Wide-Area, High-Resolution Map of Dust Extinction in M31

We map the distribution of dust in M31 at 25pc resolution, using stellar photometry from the Panchromatic Hubble Andromeda Treasury. We develop a new mapping technique that models the NIR color-magnitude diagram (CMD) of red giant branch (RGB) stars. The model CMDs combine an unreddened foreground of RGB stars with a reddened background population viewed through a log-normal column density distribution of dust. Fits to the model constrain the median extinction, the width of the extinction distribution, and the fraction of reddened stars. The resulting extinction map has >4 times better resolution than maps of dust emission, while providing a more direct measurement of the dust column. There is superb morphological agreement between the new map and maps of the extinction inferred from dust emission by Draine et al. 2014. However, the widely-used Draine & Li (2007) dust models overpredict the observed extinction by a factor of ~2.5, suggesting that M31's true dust mass is lower and that dust grains are significantly more emissive than assumed in Draine et al. (2014). The discrepancy we identify is consistent with similar findings in the Milky Way by the Planck Collaboration (2015), but has a more complex dependence on parameters from the Draine & Li (2007) dust models. We also show that the discrepancy with the Draine et al. (2014) map is lowest where the interstellar radiation field has a harder spectrum than average. We discuss possible improvements to the CMD dust mapping technique, and explore further applications.

preprint2014arXiv

An Infrared Census of DUST in Nearby Galaxies with Spitzer (DUSTiNGS), I. Overview

Nearby resolved dwarf galaxies provide excellent opportunities for studying the dust-producing late stages of stellar evolution over a wide range of metallicity (-2.7 < [Fe/H] < -1.0). Here, we describe DUSTiNGS (DUST in Nearby Galaxies with Spitzer): a 3.6 and 4.5 micron post-cryogen Spitzer Space Telescope imaging survey of 50 dwarf galaxies within 1.5 Mpc that is designed to identify dust-producing Asymptotic Giant Branch (AGB) stars and massive stars. The survey includes 37 dwarf spheroidal, 8 dwarf irregular, and 5 transition-type galaxies. This near-complete sample allows for the building of statistics on these rare phases of stellar evolution over the full metallicity range. The photometry is >75% complete at the tip of the Red Giant Branch for all targeted galaxies, with the exception of the crowded inner regions of IC 10, NGC 185, and NGC 147. This photometric depth ensures that the majority of the dust-producing stars, including the thermally-pulsing AGB stars, are detected in each galaxy. The images map each galaxy to at least twice the half-light radius to ensure that the entire evolved star population is included and to facilitate the statistical subtraction of background and foreground contamination, which is severe at these wavelengths. In this overview, we describe the survey, the data products, and preliminary results. We show evidence for the presence of dust-producing AGB stars in 8 of the targeted galaxies, with metallicities as low as [Fe/H] = -1.9, suggesting that dust production occurs even at low metallicity.

preprint2014arXiv

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

The dust properties in the Large and Small Magellanic Clouds are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 micron. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a power- law emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two blackbodies with different temperatures, both modified by the same power-law emissivity (TTMBB). Using these models we investigate the origin of the submm excess; defined as the submillimeter (submm) emission above that expected from SMBB models fit to observations < 200 micron. We find that the BEMBB model produces the lowest fit residuals with pixel-averaged 500 micron submm excesses of 27% and 43% for the LMC and SMC, respectively. Adopting gas masses from previous works, the gas-to-dust ratios calculated from our the fitting results shows that the TTMBB fits require significantly more dust than are available even if all the metals present in the interstellar medium (ISM) were condensed into dust. This indicates that the submm excess is more likely to be due to emissivity variations than a second population of colder dust. We derive integrated dust masses of (7.3 +/- 1.7) x 10^5 and (8.3 +/- 2.1) times 10^4 M(sun) for the LMC and SMC, respectively. We find significant correlations between the submm excess and other dust properties; further work is needed to determine the relative contributions of fitting noise and ISM physics to the correlations.

preprint2014arXiv

Molecular hydrogen emission in the interstellar medium of the Large Magellanic Cloud

We present the detection and analysis of molecular hydrogen emission toward ten interstellar regions in the Large Magellanic Cloud. We examined low-resolution infrared spectral maps of twelve regions obtained with the Spitzer infrared spectrograph (IRS). The pure rotational 0--0 transitions of H$_2$ at 28.2 and 17.1${\,\rm μm}$ are detected in the IRS spectra for ten regions. The higher level transitions are mostly upper limit measurements except for three regions, where a 3$σ$ detection threshold is achieved for lines at 12.2 and 8.6${\,\rm μm}$. The excitation diagrams of the detected H$_2$ transitions are used to determine the warm H$_2$ gas column density and temperature. The single-temperature fits through the lower transition lines give temperatures in the range $86-137\,{\rm K}$. The bulk of the excited H$_2$ gas is found at these temperatures and contributes $\sim$5-17% to the total gas mass. We find a tight correlation of the H$_2$ surface brightness with polycyclic aromatic hydrocarbon and total infrared emission, which is a clear indication of photo-electric heating in photodissociation regions. We find the excitation of H$_2$ by this process is equally efficient in both atomic and molecular dominated regions. We also present the correlation of the warm H$_2$ physical conditions with dust properties. The warm H$_2$ mass fraction and excitation temperature show positive correlations with the average starlight intensity, again supporting H$_2$ excitation in photodissociation regions.

preprint2014arXiv

Spitzer View of Massive Star Formation in the Tidally Stripped Magellanic Bridge

The Magellanic Bridge is the nearest low-metallicity, tidally stripped environment, offering a unique high-resolution view of physical conditions in merging and forming galaxies. In this paper we present analysis of candidate massive young stellar objects (YSOs), i.e., {\it in situ, current} massive star formation (MSF) in the Bridge using {\it Spitzer} mid-IR and complementary optical and near-IR photometry. While we definitely find YSOs in the Bridge, the most massive are $\sim10 M_\odot$, $\ll45 M_\odot$ found in the Large Magellanic Cloud (LMC). The intensity of MSF in the Bridge also appears decreasing, as the most massive YSOs are less massive than those formed in the past. To investigate environmental effects on MSF, we have compared properties of massive YSOs in the Bridge to those in the LMC. First, YSOs in the Bridge are apparently less embedded than in the LMC: 81% of Bridge YSOs show optical counterparts, compared to only 56% of LMC sources with the same range of mass, circumstellar dust mass, and line-of-sight extinction. Circumstellar envelopes are evidently more porous or clumpy in the Bridge's low-metallicity environment. Second, we have used whole samples of YSOs in the LMC and the Bridge to estimate the probability of finding YSOs at a given \hi\ column density, N(HI). We found that the LMC has $\sim3\times$ higher probability than the Bridge for N(HI) $>10\times10^{20}$ cm$^{-2}$, but the trend reverses at lower N(HI). Investigating whether this lower efficiency relative to HI is due to less efficient molecular cloud formation, or less efficient cloud collapse, or both, will require sensitive molecular gas observations.

preprint2014arXiv

Techniques and Review of Absolute Flux Calibration from the Ultraviolet to the Mid-Infrared

The measurement of precise absolute fluxes for stellar sources has been pursued with increased vigor since the discovery of the dark energy and the realization that its detailed understanding requires accurate spectral energy distributions (SEDs) of redshifted Ia supernovae in the rest frame. The flux distributions of spectrophotometric standard stars were initially derived from the comparison of stars to laboratory sources of known flux but are now mostly based on calculated model atmospheres. For example, pure hydrogen white dwarf (WD) models provide the basis for the HST CALSPEC archive of flux standards. The basic equations for quantitative spectrophotometry and photometry are explained in detail. Several historical lab based flux calibrations are reviewed; and the SEDs of stars in the major on-line astronomical databases are compared to the CALSPEC reference standard spectrophotometry. There is good evidence that relative fluxes from the visible to the near-IR wavelength of ~2.5 micron are currently accurate to 1% for the primary reference standards; and new comparisons with lab flux standards show promise for improving that precision.

preprint2014arXiv

The Mass-Loss Return from Evolved Stars to the Large Magellanic Cloud II: Dust Properties for Oxygen-Rich Asymptotic Giant Branch Stars

We model multi-wavelength broadband UBVIJHKs and Spitzer IRAC and MIPS photometry and IRS spectra from the SAGE and SAGE-Spec observing programs of two oxygen-rich asymptotic giant branch (O-rich AGB) stars in the Large Magellanic Cloud (LMC) using radiative transfer models of dust shells around stars. We chose a star from each of the bright and faint O-rich AGB populations found by earlier studies of the SAGE sample in order to derive a baseline set of dust properties to be used in the construction of an extensive grid of radiative transfer models of the O-rich AGB stars found in the SAGE surveys. From the bright O-rich AGB population we chose HV 5715, and from the faint O-rich AGB population we chose SSTISAGE1C J052206.92-715017.6 (SSTSAGE052206). We found the complex indices of refraction of oxygen-deficient silicates from Ossenkopf et al (1992) and a "KMH"-like grain size distribution with gamma of -3.5, a_min of 0.01 microns, and a_0 of 0.1 microns to be reasonable dust properties for these models. There is a slight indication that the dust around the faint O-rich AGB may be more silica-rich than that around the bright O-rich AGB. Simple models of gas emission suggest a relatively extended gas envelope for the faint O-rich AGB star modeled, consistent with the relatively large dust shell inner radius for the same model. The models compute the dust mass-loss rates for the two stars to be 2.0x10^-9 M_Sun/yr and 2.3x10^-9 M_Sun/yr, respectively. When a dust-to-gas mass ratio of 0.002 is assumed for SSTSAGE052206 and HV 5715, the dust mass-loss rates imply total mass-loss rates of 1.0x10^-6 M_Sun/yr and 1.2x10^-6 M_Sun/yr, respectively. These properties of the dust shells and stars, as inferred from our models of the two stars, are found to be consistent with properties observed or assumed by detailed studies of other O-rich AGB stars in the LMC and elsewhere.

preprint2013arXiv

Herschel Exploitation of Local Galaxy Andromeda (HELGA) III: The Star Formation Law in M31

We present a detailed study of how the Star Formation Rate (SFR) relates to the interstellar medium (ISM) of M31 at ~140pc scales. The SFR is calculated using the far-ultraviolet and 24um emission, corrected for the old stellar population in M31. We find a global value for the SFR of 0.25+/-0.05Msun/yr and compare this with the SFR found using the total far-infrared (FIR) luminosity. There is general agreement in regions where young stars dominate the dust heating. Atomic hydrogen (HI) and molecular gas (traced by carbon monoxide, CO) or the dust mass is used to trace the total gas in the ISM. We show that the global surface densities of SFR and gas mass place M31 amongst a set of low-SFR galaxies in the plot of Kennicutt (1998b). The relationship between SFR and gas surface density is tested in six radial annuli across M31, assuming a power law relationship with index, N. The star formation law using total gas traced by HI and CO gives a global index of N=2.03+/-0.04, with a significant variation with radius; the highest values are observed in the 10kpc ring. We suggest that this slope is due to HI turning molecular at ~10Msun/pc2. When looking at H2 regions, we measure a higher mean SFR suggesting a better spatial correlation between H2 and SF. We find N~0.6 with consistent results throughout the disk - this is at the low end of values found in previous work and argues against a superlinear SF law on small scales.

preprint2012arXiv

The Panchromatic Hubble Andromeda Treasury IV. A Probabilistic Approach to Inferring the High Mass Stellar Initial Mass Function and Other Power-law Functions

We present a probabilistic approach for inferring the parameters of the present day power-law stellar mass function (MF) of a resolved young star cluster. This technique (a) fully exploits the information content of a given dataset; (b) accounts for observational uncertainties in a straightforward way; (c) assigns meaningful uncertainties to the inferred parameters; (d) avoids the pitfalls associated with binning data; and (e) is applicable to virtually any resolved young cluster, laying the groundwork for a systematic study of the high mass stellar MF (M > 1 Msun). Using simulated clusters and Markov chain Monte Carlo sampling of the probability distribution functions, we show that estimates of the MF slope, α, are unbiased and that the uncertainty, Δα, depends primarily on the number of observed stars and stellar mass range they span, assuming that the uncertainties on individual masses and the completeness are well-characterized. Using idealized mock data, we compute the lower limit precision on α and provide an analytic approximation for Δα as a function of the observed number of stars and mass range. We find that ~ 3/4 of quoted literature uncertainties are smaller than the theoretical lower limit. By correcting these uncertainties to the theoretical lower limits, we find the literature studies yield <α>=2.46 with a 1-σ dispersion of 0.35 dex. We verify that it is impossible for a power-law MF to obtain meaningful constraints on the upper mass limit of the IMF. We show that avoiding substantial biases in the MF slope requires: (1) including the MF as a prior when deriving individual stellar mass estimates; (2) modeling the uncertainties in the individual stellar masses; and (3) fully characterizing and then explicitly modeling the completeness for stars of a given mass. (abridged)

preprint2011arXiv

A Population of Accreted SMC Stars in the LMC

We present an analysis of the stellar kinematics of the Large Magellanic Cloud based on ~5900 new and existing velocities of massive red supergiants, oxygen-rich and carbon-rich AGB stars, and other giants. After correcting the line-of-sight velocities for the LMC's space motion and accounting for asymmetric drift in the AGB population, we derive a rotation curve that is consistent with all of the tracers used, as well as that of published HI data. The amplitude of the rotation curve is v_0=87+/-5 km s^-1 beyond a radius R_0=2.4+/-0.1 kpc, and has a position angle of the kinematic line of nodes of theta=142 degrees +/-5 degrees. By examining the outliers from our fits, we identify a population of 376 stars, or >~5% of our sample, that have line-of-sight velocities that apparently oppose the sense of rotation of the LMC disk. We find that these kinematically distinct stars are either counter-rotating in a plane closely aligned with the LMC disk, or rotating in the same sense as the LMC disk, but in a plane that is inclined by 54 degrees +/- 2 degrees to the LMC. Their kinematics clearly link them to two known HI arms, which have previously been interpreted as being pulled out from the LMC. We measure metallicities from the Ca triplet lines of ~1000 LMC field stars and 30 stars in the kinematically distinct population. For the LMC field, we find a median [Fe/H]=-0.56 +/- 0.02 with dispersion of 0.5 dex, while for the kinematically distinct stars the median [Fe/H] is -1.25 +/- 0.13 with a dispersion of 0.7 dex. The metallicity differences provide strong evidence that the kinematically distinct population originated in the SMC. This interpretation has the consequence that the HI arms kinematically associated with the stars are likely falling into the LMC, instead of being pulled out.

preprint2011arXiv

Optimal Cosmic-Ray Detection for Nondestructive Read Ramps

Cosmic rays are a known problem in astronomy, causing both loss of data and data inaccuracy. The problem becomes even more extreme when considering data from a high-radiation environment, such as in orbit around Earth or outside the Earth's magnetic field altogether, unprotected, as will be the case for the James Webb Space Telescope (JWST). For JWST, all the instruments employ nondestructive readout schemes. The most common of these will be "up the ramp" sampling, where the detector is read out regularly during the ramp. We study three methods to correct for cosmic rays in these ramps: a two-point difference method, a deviation from the fit method, and a y-intercept method. We apply these methods to simulated nondestructive read ramps with single-sample groups and varying combinations of flux, number of samples, number of cosmic rays, cosmic-ray location in the exposure, and cosmic-ray strength. We show that the y-intercept method is the optimal detection method in the read-noise-dominated regime, while both the y-intercept method and the two-point difference method are best in the photon-noise-dominated regime, with the latter requiring fewer computations.

preprint2011arXiv

Surveying the Agents of Galaxy Evolution in the Tidally-Stripped, Low Metallicity Small Magellanic Cloud (SAGE-SMC) II. Cool Evolved Stars

We investigate the infrared (IR) properties of cool, evolved stars in the Small Magellanic Cloud (SMC), including the red giant branch (RGB) stars and the dust-producing red supergiant (RSG) and asymptotic giant branch (AGB) stars using observations from the Spitzer Space Telescope Legacy program entitled: "Surveying the Agents of Galaxy Evolution in the Tidally-stripped, Low Metallicity SMC", or SAGE-SMC. The survey includes, for the first time, full spatial coverage of the SMC bar, wing, and tail regions at infrared (IR) wavelengths (3.6 - 160 microns). We identify evolved stars using a combination of near-IR and mid-IR photometry and point out a new feature in the mid-IR color-magnitude diagram that may be due to particularly dusty O-rich AGB stars. We find that the RSG and AGB stars each contribute ~20% of the global SMC flux (extended + point-source) at 3.6 microns, which emphasizes the importance of both stellar types to the integrated flux of distant metal-poor galaxies. The equivalent SAGE survey of the higher-metallicity Large Magellanic Cloud (SAGE-LMC) allows us to explore the influence of metallicity on dust production. We find that the SMC RSG stars are less likely to produce a large amount of dust (as indicated by the [3.6]-[8] color). There is a higher fraction of carbon-rich stars in the SMC, and these stars appear to able to reach colors as red as their LMC counterparts, indicating that C-rich dust forms efficiently in both galaxies. A preliminary estimate of the dust production in AGB and RSG stars reveals that the extreme C-rich AGB stars dominate the dust input in both galaxies, and that the O-rich stars may play a larger role in the LMC than in the SMC.

preprint2011arXiv

Surveying the Agents of Galaxy Evolution in the Tidally-Stripped, Low Metallicity Small Magellanic Cloud (SAGE-SMC). I. Overview

The Small Magellanic Cloud (SMC) provides a unique laboratory for the study of the lifecycle of dust given its low metallicity (~1/5 solar) and relative proximity (~60 kpc). This motivated the SAGE-SMC (Surveying the Agents of Galaxy Evolution in the Tidally-Stripped, Low Metallicity Small Magellanic Cloud) Spitzer Legacy program with the specific goals of studying the amount and type of dust in the present interstellar medium, the sources of dust in the winds of evolved stars, and how much dust is consumed in star formation. This program mapped the full SMC (30 sq. deg.) including the Body, Wing, and Tail in 7 bands from 3.6 to 160 micron using the IRAC and MIPS instruments on the Spitzer Space Telescope. The data were reduced, mosaicked, and the point sources measured using customized routines specific for large surveys. We have made the resulting mosaics and point source catalogs available to the community. The infrared colors of the SMC are compared to those of other nearby galaxies and the 8 micron/24 micron ratio is somewhat lower and the 70 micron/160 micron ratio is somewhat higher than the average. The global infrared spectral energy distribution shows that the SMC has ~3X lower aromatic emission/PAH (polycyclic aromatic hydrocarbon) abundances compared to most nearby galaxies. Infrared color-magnitude diagrams are given illustrating the distribution of different asymptotic giant branch stars and the locations of young stellar objects. Finally, the average spectral energy distribution (SED) of HII/star formation regions is compared to the equivalent Large Magellanic Cloud average HII/star formation region SED. These preliminary results are expanded in detail in companion papers.

preprint2011arXiv

The Emission by Dust and Stars of Nearby Galaxies in the Herschel KINGFISH Survey

Using new far-infrared imaging from the Herschel Space Observatory with ancillary data from ultraviolet to submillimeter wavelengths, we estimate the total emission from dust and stars of 62 nearby galaxies in the KINGFISH survey in a way that is as empirical and model-independent as possible. We collect and exploit these data in order to measure from the spectral energy distributions (SEDs) precisely how much stellar radiation is intercepted and re-radiated by dust, and how this quantity varies with galaxy properties. By including SPIRE data, we are more sensitive to emission from cold dust grains than previous analyses at shorter wavelengths, allowing for more accurate estimates of dust temperatures and masses. The dust/stellar flux ratio, which we measure by integrating the SEDs, has a range of nearly three decades. The inclusion of SPIRE data shows that estimates based on data not reaching these far-IR wavelengths are biased low. We find that the dust/stellar flux ratio varies with morphology and total IR luminosity. We also find that dust/stellar flux ratios are related to gas-phase metallicity, while the dust/stellar mass ratios are less so. The substantial scatter between dust/stellar flux and dust/stellar mass indicates that the former is a poor proxy of the latter. Comparing the dust/stellar flux ratios and dust temperatures, we show that early-types tend to have slightly warmer temperatures than spiral galaxies, which may be due to more intense interstellar radiation fields, or to different dust grain compositions. Finally, we show that early-types and early-type spirals have a strong correlation between the dust/stellar flux ratio and specific star formation rate, which suggests that the relatively bright far-IR emission of some of these galaxies is due to ongoing star formation and the radiation field from older stars.

preprint2011arXiv

Young, UV-bright Stars Dominate Dust Heating in Star Forming Galaxies

In star forming galaxies, dust plays a significant role in shaping the ultraviolet (UV) through infrared (IR) spectrum. Dust attenuates the radiation from stars, and re-radiates the energy through equilibrium and non-equilibrium emission. Polycyclic aromatic hydrocarbons (PAH), graphite, and silicates contribute to different features in the spectral energy distribution; however, they are all highly opaque in the same spectral region - the UV. Compared to old stellar populations, young populations release a higher fraction of their total luminosity in the UV, making them a good source of the energetic UV photons that can power dust emission. However, given their relative abundance, the question of whether young or old stellar populations provide most of these photons that power the infrared emission is an interesting question. Using three samples of galaxies observed with the Spitzer Space Telescope and our dusty radiative transfer model, we find that young stellar populations (on the order of 100 million years old) dominate the dust heating in star forming galaxies, and old stellar populations (13 billion years old) generally contribute less than 20% of the far-IR luminosity.

preprint2010arXiv

A Spitzer Space Telescope far-infrared spectral atlas of compact sources in the Magellanic Clouds. II. The Small Magellanic Cloud

We present 52-93 micron spectra, obtained with the Spitzer Space Telescope, of luminous compact far-IR sources in the SMC. These comprise 9 Young Stellar Objects (YSOs), the compact HII region N81 and a similar object within N84, and two red supergiants (RSGs). The spectra of the sources in N81 (of which we also show the ISO-LWS spectrum between 50-170 micron) and N84 both display strong [OI] 63-micron and [OIII] 88-micron fine-structure line emission. We attribute these lines to strong shocks and photo-ionized gas, respectively, in a ``champagne flow'' scenario. The nitrogen content of these two HII regions is very low, definitely N/O<0.04 but possibly as low as N/O<0.01. Overall, the oxygen lines and dust continuum are weaker in star-forming objects in the SMC than in the LMC. We attribute this to the lower metallicity of the SMC compared to that of the LMC. Whilst the dust mass differs in proportion to metallicity, the oxygen mass differs less; both observations can be reconciled with higher densities inside star-forming cloud cores in the SMC than in the LMC. The dust in the YSOs in the SMC is warmer (37-51 K) than in comparable objects in the LMC (32-44 K). We attribute this to the reduced shielding and reduced cooling at the low metallicity of the SMC. On the other hand, the efficiency of the photo-electric effect to heat the gas is found to be indistinguishable to that measured in the same manner in the LMC, 0.1-0.3%. This may result from higher cloud-core densities, or smaller grains, in the SMC. The dust associated with the two RSGs in our SMC sample is cool, and we argue that it is swept-up interstellar dust, or formed (or grew) within the bow-shock, rather than dust produced in these metal-poor RSGs themselves. Strong emission from crystalline water ice is detected in at least one YSO. (abridged)

preprint2010arXiv

A survey of diffuse interstellar bands in the Andromeda galaxy: optical spectroscopy of M31 OB stars

We present the largest sample to-date of intermediate-resolution blue-to-red optical spectra of B-type supergiants in M31 and undertake the first survey of diffuse interstellar bands (DIBs) in this galaxy. Spectral classifications, radial velocities and interstellar reddenings are presented for 34 stars in three regions of M31. Radial velocities and equivalent widths are given for the 5780 and 6283 DIBs towards 11 stars. Equivalent widths are also presented for the following DIBs detected in three sightlines in M31: 4428, 5705, 5780, 5797, 6203, 6269, 6283, 6379, 6613, 6660, and 6993. All of these M31 DIB carriers reside in clouds at radial velocities matching those of interstellar Na I and/or H I. The relationships between DIB equivalent widths and reddening (E(B-V)) are consistent with those observed in the local ISM of the Milky Way. Many of the observed sightlines show DIB strengths (per unit reddening) which lie at the upper end of the range of Galactic values. DIB strengths per unit reddening are found (with 68% confidence), to correlate with the interstellar UV radiation field strength. The strongest DIBs are observed where the interstellar UV flux is lowest. The mean Spitzer 8/24 micron emission ratio in our three fields is slightly lower than that measured in the Milky Way, but we identify no correlation between this ratio and the DIB strengths in M31. Interstellar oxygen abundances derived from the spectra of three M31 H II regions in one of the fields indicate that the average metallicity of the ISM in that region is 12 + log[O/H] = 8.54 +- 0.18, which is approximately equal to the value in the solar neighbourhood.

preprint2010arXiv

Spitzer Analysis of HII Region Complexes in the Magellanic Clouds: Determining a Suitable Monochromatic Obscured Star Formation Indicator

HII regions are the birth places of stars, and as such they provide the best measure of current star formation rates (SFRs) in galaxies. The close proximity of the Magellanic Clouds allows us to probe the nature of these star forming regions at small spatial scales. We aim to determine the monochromatic IR band that most accurately traces the bolometric IR flux (TIR), which can then be used to estimate an obscured SFR. We present the spatial analysis, via aperture/annulus photometry, of 16 LMC and 16 SMC HII region complexes using the Spitzer IRAC and MIPS bands. UV rocket data and SHASSA H-alpha data are also included. We find that nearly all of the LMC and SMC HII region SEDs peak around 70um, from ~10 to ~400 pc from the central sources. As a result, the sizes of HII regions as probed by 70um is approximately equal to the sizes as probed by TIR (about 70 pc in radius); the radial profile of the 70um flux, normalized by TIR, is constant at all radii (70um ~ 0.45 TIR); the 1-sigma standard deviation of the 70um fluxes, normalized by TIR, is a lower fraction of the mean (0.05 to 0.12 out to ~220 pc) than the normalized 8, 24, and 160um normalized fluxes (0.12 to 0.52); and these results are invariant between the LMC and SMC. From these results, we argue that 70um is the most suitable IR band to use as a monochromatic obscured star formation indicator because it most accurately reproduces the TIR of HII regions in the LMC and SMC and over large spatial scales. We also explore the general trends of the 8, 24, 70, and 160um bands in the LMC and SMC HII region SEDs, radial surface brightness profiles, sizes, and normalized (by TIR) radial flux profiles. We derive an obscured SFR equation that is modified from the literature to use 70um luminosity, SFR(Mo/yr) = 9.7(0.7)x10^{-44} L(70)(ergs/s), which is applicable from 10 to 300 pc distance from the center of an HII region.

preprint2010arXiv

The SAGE-Spec Spitzer Legacy program: The life-cycle of dust and gas in the Large Magellanic Cloud

The SAGE-Spec Spitzer Legacy program is a spectroscopic follow-up to the SAGE-LMC photometric survey of the Large Magellanic Cloud carried out with the Spitzer Space Telescope. We present an overview of SAGE-Spec and some of its first results. The SAGE-Spec program aims to study the life cycle of gas and dust in the Large Magellanic Cloud, and to provide information essential to the classification of the point sources observed in the earlier SAGE-LMC photometric survey. We acquired 224.6 hours of observations using the InfraRed Spectrograph and the SED mode of the Multiband Imaging Photometer for Spitzer. The SAGE-Spec data, along with archival Spitzer spectroscopy of objects in the Large Magellanic Cloud, are reduced and delivered to the community. We discuss the observing strategy, the specific data reduction pipelines applied and the dissemination of data products to the scientific community. Initial science results include the first detection of an extragalactic "21 um" feature towards an evolved star and elucidation of the nature of disks around RV Tauri stars in the Large Magellanic Cloud. Towards some young stars, ice features are observed in absorption. We also serendipitously observed a background quasar, at a redshift of z~0.14, which appears to be host-less.

preprint2010arXiv

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

We present the classification of 197 point sources observed with the Infrared Spectrograph in the SAGE-Spec Legacy program 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 HII regions and 12 other objects, seven of which remain unclassified.

preprint2009arXiv

A Spitzer Space Telescope far-infrared spectral atlas of compact sources in the Magellanic Clouds. I. The Large Magellanic Cloud

[abridged] We present 52-93 micron spectra obtained with Spitzer in the MIPS-SED mode, of a representative sample of luminous compact far-IR sources in the LMC. These include carbon stars, OH/IR AGB stars, post-AGB objects and PNe, RCrB-type star HV2671, OH/IR red supergiants WOHG064 and IRAS05280-6910, B[e] stars IRAS04530-6916, R66 and R126, Wolf-Rayet star Brey3a, Luminous Blue Variable R71, supernova remnant N49, a large number of young stellar objects, compact HII regions and molecular cores, and a background galaxy (z~0.175). We use the spectra to constrain the presence and temperature of cold dust and the excitation conditions and shocks within the neutral and ionized gas, in the circumstellar environments and interfaces with the surrounding ISM. Evolved stars, including LBV R71, lack cold dust except in some cases where we argue that this is swept-up ISM. This leads to an estimate of the duration of the prolific dust-producing phase ("superwind") of several thousand years for both RSGs and massive AGB stars, with a similar fractional mass loss experienced despite the different masses. We tentatively detect line emission from neutral oxygen in the extreme RSG WOHG064, with implications for the wind driving. In N49, the shock between the supernova ejecta and ISM is revealed by its strong [OI] 63-micron emission and possibly water vapour; we estimate that 0.2 Msun of ISM dust was swept up. Some of the compact HII regions display pronounced [OIII] 88-micron emission. The efficiency of photo-electric heating in the interfaces of ionized gas and molecular clouds is estimated at 0.1-0.3%. We confirm earlier indications of a low nitrogen content in the LMC. Evidence for solid state emission features is found in both young and evolved object; some of the YSOs are found to contain crystalline water ice.

preprint2009arXiv

Dust Production and Mass Loss in the Galactic Globular Cluster NGC 362

We investigate dust production and stellar mass loss in the Galactic globular cluster NGC 362. Due to its close proximity to the Small Magellanic Cloud (SMC), NGC 362 was imaged with the IRAC and MIPS cameras onboard the Spitzer Space Telescope as part of the Surveying the Agents of Galaxy Evolution (SAGE-SMC) Spitzer Legacy program. We detect several cluster members near the tip of the Red Giant Branch that exhibit infrared excesses indicative of circumstellar dust and find that dust is not present in measurable quantities in stars below the tip of the Red Giant Branch. We modeled the spectral energy distribution (SED) of the stars with the strongest IR excess and find a total cluster dust mass-loss rate of 3.0(+2.0/-1.2) x 10^-9 solar masses per year, corresponding to a gas mass-loss rate of 8.6(+5.6/-3.4) x 10^-6 solar masses per year, assuming [Fe/H] = -1.16. This mass loss is in addition to any dust-less mass loss that is certainly occurring within the cluster. The two most extreme stars, variables V2 and V16, contribute up to 45% of the total cluster dust-traced mass loss. The SEDs of the more moderate stars indicate the presence of silicate dust, as expected for low-mass, low-metallicity stars. Surprisingly, the SED shapes of the stars with the strongest mass-loss rates appear to require the presence of amorphous carbon dust, possibly in combination with silicate dust, despite their oxygen-rich nature. These results corroborate our previous findings in omega Centauri.

preprint2009arXiv

FUSE Measurements of Far-Ultraviolet Extinction. III. The Dependence on R(V) and Discrete Feature Limits from 75 Galactic Sightlines

We present a sample of 75 extinction curves derived from FUSE far-ultraviolet spectra supplemented by existing IUE spectra. The extinction curves were created using the standard pair method based on a new set of dereddened FUSE+IUE comparison stars. Molecular hydrogen absorption features were removed using individualized H_2 models for each sightline. The general shape of the FUSE extinction (8.4 micron^-1 < lambda^-1 < 11 micron^-1) was found to be broadly consistent with extrapolations from the IUE extinction (3.3 micron-1 < lambda^-1 < 8.6 micron^-1) curve. Significant differences were seen in the strength of the far-UV rise and the width of the 2175 A bump. All the FUSE+IUE extinction curves had positive far-UV slopes giving no indication that the far-UV rise was turning over at the shortest wavelengths. The dependence of A(lambda)/A(V) versus R(V)^-1 in the far-UV using the sightlines in our sample was found to be stronger than tentatively indicated by previous work. We present an updated R(V) dependent relationship for the full UV wavelength range (3.3 micron^-1 <= lambda^-1 <= 11 micron^-1). Finally, we searched for discrete absorption features in the far-ultraviolet. We found a 3 sigma upper limit of ~0.12 A(V) on features with a resolution of 250 (~4 A width) and 3 sigma upper limits of ~0.15 A(V) for lambda^-1 < 9.6 micron^-1 and ~0.68 A(V) for lambda^-1 > 9.6 micron^-1 on features with a resolution of 10^4 (~0.1 A width).

preprint2009arXiv

Spitzer View of Young Massive Stars in the LMC HII Complex N44

The HII complex N44 in the Large Magellanic Cloud (LMC) provides an excellent site to perform a detailed study of star formation in a mild starburst, as it hosts three regions of star formation at different evolutionary stages and it is not as complicated and confusing as the 30 Doradus giant HII region. We have obtained Spitzer Space Telescope observations and complementary ground-based 4m uBVIJK observations of N44 to identify candidate massive young stellar objects (YSOs). We further classify the YSOs into Types I, II, and III, according to their spectral energy distributions (SEDs). In our sample of 60 YSO candidates, ~65% of them are resolved into multiple components or extended sources in high-resolution ground-based images. We have modeled the SEDs of 36 YSOs that appear single or dominant within a group. We find good fits for Types I and I/II YSOs,but Types II and II/III YSOs show deviations between their observed SEDs and models that do not include PAH emission. We have also found that some Type III YSOs have central holes in their disk components. YSO counterparts are found in four ultracompact HII regions and their stellar masses determined from SED model fits agree well with those estimated from the ionization requirements of the HII regions. The distribution of YSOs is compared with those of the underlying stellar population and interstellar gas conditions to illustrate a correlation between the current formation of O-type stars and previous formation of massive stars. Evidence of triggered star formation is also presented.

preprint2009arXiv

Warm and Dense Molecular Gas in the N159 Region: 12CO J=4-3 and 13CO J=3-2 Observations with NANTEN2 and ASTE

New 12CO J=4-3 and 13CO J=3-2 observations of the N159 region in the Large Magellanic Cloud have been made. The 12CO J=4-3 distribution is separated into three clumps. These new measurements toward the three clumps are used in coupled calculations of molecular rotational excitation and line radiation transfer, along with other transitions of the 12CO as well as the isotope transitions of 13CO. The temperatures and densities are determined to be ~70-80K and ~3x10^3 cm-3 in N159W and N159E and ~30K and ~1.6x10^3 cm-3 in N159S. These results are compared with the star formation activity. The N159E clump is associated with embedded cluster(s) as observed at 24 micron and the derived high temperature is explained as due to the heating by these sources. The N159E clump is likely responsible for a dark lane in a large HII region by the dust extinction. The N159W clump is associated with embedded clusters mainly toward the eastern edge of the clump only. These clusters show offsets of 20"-40" from the 12CO J=4-3 peak and are probably responsible for heating indicated by the derived high temperature. The N159W clump exhibits no sign of star formation toward the 12CO J=4-3 peak position and its western region. We suggest that the N159W peak represents a pre-star-cluster core of ~105M_sol which deserves further detailed studies. Note that recent star formation took place between N159W and N159E as indicated by several star clusters and HII regions, while the natal molecular gas toward the stars have already been dissipated by the ionization and stellar winds of the OB stars. The N159S clump shows little sign of star formation as is consistent with the lower temperature and somewhat lower density. The N159S clump is also a candidate for future star formation.

preprint2007arXiv

Absolute Calibration and Characterization of the Multiband Imaging Photometer for Spitzer. II. 70 micron Imaging

The absolute calibration and characterization of the Multiband Imaging Photometer for Spitzer (MIPS) 70 micron coarse- and fine-scale imaging modes are presented based on over 2.5 years of observations. Accurate photometry (especially for faint sources) requires two simple processing steps beyond the standard data reduction to remove long-term detector transients. Point spread function (PSF) fitting photometry is found to give more accurate flux densities than aperture photometry. Based on the PSF fitting photometry, the calibration factor shows no strong trend with flux density, background, spectral type, exposure time, or time since anneals. The coarse-scale calibration sample includes observations of stars with flux densities from 22 mJy to 17 Jy, on backgrounds from 4 to 26 MJy sr^-1, and with spectral types from B to M. The coarse-scale calibration is 702 +/- 35 MJy sr^-1 MIPS70^-1 (5% uncertainty) and is based on measurements of 66 stars. The instrumental units of the MIPS 70 micron coarse- and fine-scale imaging modes are called MIPS70 and MIPS70F, respectively. The photometric repeatability is calculated to be 4.5% from two stars measured during every MIPS campaign and includes variations on all time scales probed. The preliminary fine-scale calibration factor is 2894 +/- 294 MJy sr^-1 MIPS70F^-1 (10% uncertainty) based on 10 stars. The uncertainty in the coarse- and fine-scale calibration factors are dominated by the 4.5% photometric repeatability and the small sample size, respectively. The 5-sigma, 500 s sensitivity of the coarse-scale observations is 6-8 mJy. This work shows that the MIPS 70 micron array produces accurate, well calibrated photometry and validates the MIPS 70 micron operating strategy, especially the use of frequent stimulator flashes to track the changing responsivities of the Ge:Ga detectors.

preprint2007arXiv

Star Formation in NGC 5194 (M51a). II. The Spatially-Resolved Star Formation Law

We have studied the relationship between the star formation rate (SFR) surface density and gas surface density in the spiral galaxy M51a (NGC 5194), using multi-wavelength data obtained as part of the Spitzer Infrared Nearby Galaxies Survey (SINGS). We introduce a new SFR index based on a linear combination of H-alpha emission-line and 24 micron continuum luminosities, that provides reliable extinction-corrected ionizing fluxes and SFR densities over a wide range of dust attenuations. The combination of these extinction-corrected SFR densities with aperture synthesis HI and CO maps has allowed us to probe the form of the spatially-resolved star formation law on scales of 0.5 to 2 kpc. We find that the resolved SFR vs gas surface density relation is well represented by a Schmidt power law, which is similar in form and dispersion to the disk-averaged Schmidt law. We observe a comparably strong correlation of the SFR surface density with the molecular gas surface density, but no significant correlation with the surface density of atomic gas. The best-fitting slope of the Schmidt law varies from N = 1.37 to 1.56, with zeropoint and slope that change systematically with the spatial sampling scale. We tentatively attribute these variations to the effects of areal sampling and averaging of a nonlinear intrinsic star formation law. Our data can also be fitted by an alternative parametrization of the SFR surface density in terms of the ratio of gas surface density to local dynamical time, but with a considerable dispersion.

preprint2006arXiv

Massive-Star Supernovae as Major Dust Factories

We present late-time optical and mid-infrared observations of the Type-II supernova 2003gd in NGC 628. Mid-infrared excesses consistent with cooling dust in the ejecta are observed 499-678 days after outburst, and are accompanied by increasing optical extinction and growing asymmetries in the emission-line profiles. Radiative-transfer models show that up to 0.02 solar masses of dust has formed within the ejecta, beginning as early as 250 days after outburst. These observations show that dust formation in supernova ejecta can be efficient and that massive-star supernovae can be major dust producers throughout the history of the Universe.

preprint2006arXiv

Spitzer Survey of the Large Magellanic Cloud, Surveying the Agents of a Galaxy's Evolution (SAGE) I: Overview and Initial Results

We are performing a uniform and unbiased, ~7x7 degrees imaging survey of the Large Magellanic Cloud (LMC), using the IRAC and MIPS instruments on board the Spitzer Space Telescope in order to survey the agents of a galaxy's evolution (SAGE), the interstellar medium (ISM) and stars in the LMC. The detection of diffuse ISM with column densities >1.2x10^21 H cm^-2 permits detailed studies of dust processes in the ISM. SAGE's point source sensitivity enables a complete census of newly formed stars with masses >3 solar masses that will determine the current star formation rate in the LMC. SAGE's detection of evolved stars with mass loss rates >1x10^-8 solar masses per year will quantify the rate at which evolved stars inject mass into the ISM of the LMC. The observing strategy includes two epochs in 2005, separated by three months, that both mitigate instrumental artifacts and constrain source variability. The SAGE data are non-proprietary. The data processing includes IRAC and MIPS pipelines and a database for mining the point source catalogs, which will be released to the community in support of Spitzer proposal cycles 4 and 5. We present initial results on the epoch 1 data with a special focus on the N79 and N83 region. The SAGE epoch 1 point source catalog has ~4 million sources. The point source counts are highest for the IRAC 3.6 microns band and decrease dramatically towards longer wavelengths consistent with the fact that stars dominate the point source catalogs and that the dusty objects, e.g. young stellar objects and dusty evolved stars that detected at the longer wavelengths, are rare in comparison. We outline a strategy for identifying foreground MW stars, that may comprise as much as 18% of the source list, and background galaxies, that may comprise ~12% of the source list.

preprint2006arXiv

The Spectral Energy Distribution of Dust Emission in the Edge-on Spiral Galaxy NGC 4631 as seen with Spitzer and JCMT

We explore the nature of variations in dust emission within an individual galaxy using 3.6 - 160 micron Spitzer Space Telescope observations and 450 and 850 micron James Clerk Maxwell Telescope observations of the edge-on Sd spiral galaxy NGC 4631 with the goals of understanding the relation between polycyclic aromatic hydrocarbons (PAH) and dust emission, studying the variations in the colors of the dust emission, and searching for possible excess submillimeter emission compared to what is anticipated based on dust models applied to the mid- and far-infrared data. PAH emission at 8 micron is found to correlate best with hot dust emission at 24 micron on kiloparsec scales, although the relation breaks down on scales equal to hundreds of parsecs, possibly because of differences in the mean free paths between the photons that excite the PAHs and heat the dust and possibly because the PAHs are destroyed by the hard radiation fields in the centers of some star formation regions. The ratio of 8 micron PAH emission to 160 micron cool dust emission appears to be a function of radius. The 70/160 and 160/450 micron flux density ratios are remarkably constant even though the surface brightness varies by a factor of 25 in each wave band, which suggests that the emission is from dust heated by a nearly-uniform radiation field. Globally, we find an excess of 850 micron emission relative to what would be predicted by dust models. The 850 micron excess is highest in regions with low 160 micron surface brightness, although the strength and statistical significance of this result depends on the model fit to the data. We rule out variable emissivity functions or ~4 K dust as the possible origins of this 850 micron emission, but we do discuss the other possible mechanisms that could produce the emission.

preprint2006arXiv

Warm Dust and Spatially Variable PAH Emission in the Dwarf Starburst Galaxy NGC 1705

We present Spitzer observations of the dwarf starburst galaxy NGC 1705 obtained as part of SINGS. The galaxy morphology is very different shortward and longward of ~5 microns: short-wavelength imaging shows an underlying red stellar population, with the central super star cluster (SSC) dominating the luminosity; longer-wavelength data reveals warm dust emission arising from two off-nuclear regions offset by ~250 pc from the SSC. These regions show little extinction at optical wavelengths. The galaxy has a relatively low global dust mass (~2E5 solar masses, implying a global dust-to-gas mass ratio ~2--4 times lower than the Milky Way average). The off-nuclear dust emission appears to be powered by photons from the same stellar population responsible for the excitation of the observed H Alpha emission; these photons are unassociated with the SSC (though a contribution from embedded sources to the IR luminosity of the off-nuclear regions cannot be ruled out). Low-resolution IRS spectroscopy shows moderate-strength PAH emission in the 11.3 micron band in the eastern peak; no PAH emission is detected in the SSC or the western dust emission complex. There is significant diffuse 8 micron emission after scaling and subtracting shorter wavelength data; the spatially variable PAH emission strengths revealed by the IRS data suggest caution in the interpretation of diffuse 8 micron emission as arising from PAH carriers alone. The metallicity of NGC 1705 falls at the transition level of 35% solar found by Engelbracht and collaborators; the fact that a system at this metallicity shows spatially variable PAH emission demonstrates the complexity of interpreting diffuse 8 micron emission. A radio continuum non-detection, NGC 1705 deviates significantly from the canonical far-IR vs. radio correlation. (Abridged)

preprint2005arXiv

An Extremely Large Excess of 18O in the Hydrogen-Deficient Carbon Star, HD 137613

We report the discovery of a uniquely large excess of 18O in the hydrogen-deficient carbon (HdC) star, HD 137613, based on a spectrum of the first overtone bands of CO at 2.3-2.4 micron in which three strong absorption bands of 12C18O are clearly present. Bands of 12C16O also are present but no bands of 13C16O or 12C17O are seen. We estimate an isotopic ratio 16O/18O ~1. The Solar value of this ratio is ~500. Neither He-core burning nor He-shell flash burning can produce the isotopic ratios of oxygen and carbon observed in HD 137613. However, a remarkable similarity exists between the observed abundances and those found in the outer layers of the broad He-shell of early-AGB stars, soon after the end of He-core burning. It is not known how the outer envelope down to the He-shell could be lost but some mechanism of enhanced mass loss must be involved. HD 137613 may be a post-early-AGB star with the outer layers of the former He-burning shell as its photosphere. The unusual elemental abundances of the HdC stars resemble those of the R Coronae Borealis (RCB) stars, but HdC stars do not produce clouds of dust that produce declines in brightness. None of the other RCB or HdC stars observed shows significant 18O.

preprint2005arXiv

Discovery of Five New R Coronae Borealis Stars in the MACHO Galactic Bulge Database

We have identified five new R Coronae Borealis (RCB) stars in the Galactic bulge using the MACHO Project photometry database, raising the total number of known Galactic RCB stars to about 40. We have obtained spectra to confirm the identifications. The fact that four out of the five newly identified RCB stars are ``cool'' (T(eff) < 6000 K) rather than ``warm'' (T(eff) > 6000 K) suggests that the preponderance of warm RCB stars among the existing sample is a selection bias. These cool RCB stars are redder and fainter than their warm counterparts and may have been missed in surveys done with blue plates. Based on the number of new RCB stars discovered in the MACHO bulge fields, there may be ~250 RCB stars in the reddened "exclusion" zone toward the bulge.

preprint2005arXiv

Spitzer Observations Of The Supergiant Shell Region In IC 2574

We present spatially resolved Spitzer imaging of the supergiant shell region of the M81 group dwarf galaxy IC 2574 obtained as part of the Spitzer Infrared Nearby Galaxies Survey. This region harbors one of the best nearby examples of a kinematically distinct HI shell, with an associated remnant stellar cluster; the shell is initiating sequential star formation as it interacts with the surrounding interstellar medium. This region dominates the infrared luminosity of IC 2574 and is spatially resolved in all Spitzer imaging bands. We study the differences in dust temperature as a function of local environment and compare local star formation rates as inferred from H Alpha and total infrared luminosities. We find that the strong H Alpha sources are associated with regions of warm dust; however, the most luminous infrared and H Alpha sources are not necessarilyco-spatial. The coolest dust is found in the regions farthest from the rim of the shell; these regions show the best agreement between star formation rates derived from H Alpha and from total infrared luminosities (although discrepancies at the factor of 3-4 level still exist). There is considerable variation in the radio-far infrared correlation in different regions surrounding the shell. The low dust content of the region may influence the scatter seen in these relations; these data demonstrate that the expanding shell is dramatically affecting its surroundings by triggering star formation and altering the dust temperature.

preprint2005arXiv

Understanding Radio-Selected Thermal Sources in M 33: Ultraviolet, Optical, Near-Infrared, Spitzer Mid-Infrared, and Radio Observations

We present ultraviolet, optical, near-infrared, Spitzer mid-infrared, and radio images of 14 radio-selected objects in M 33. These objects are thought to represent the youngest phase of star cluster formation. We have detected the majority of cluster candidates in M 33 at all wavelengths. From the near-IR images, we derived ages 2-10 Myr, K_S-band extinctions (A_K_S) of 0-1 mag, and stellar masses of 10^3-10^4 M_solar. We have generated spectral energy distributions (SEDs) of each cluster from 0.1 micron to 160 microns. From these SEDs, we have modeled the dust emission around these star clusters to determine the dust masses (1-10^3 M_solar) and temperatures (40-90 K) of the clusters' local interstellar medium. Extinctions derived from the JHK_S, Halpha, and UV images are similar to within a factor of 2 or 3. These results suggest that eleven of the fourteen radio-selected objects are optically-visible young star clusters with a surrounding H II region, that two are background objects, possibly AGN, and that one is a Wolf-Rayet star with a surrounding H II region.