• 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.
  • 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.
  • The Astropy Problem (1610.03159)

    Demitri Muna, Michael Alexander, Alice Allen, Richard Ashley, Daniel Asmus, Ruyman Azzollini, Michele Bannister, Rachael Beaton, Andrew Benson, G. Bruce Berriman, Maciej Bilicki, Peter Boyce, Joanna Bridge, Jan Cami, Eryn Cangi, Xian Chen, Nicholas Christiny, Christopher Clark, Michelle Collins, Johan Comparat, Neil Cook, Darren Croton, Isak Delberth Davids, Éric Depagne, John Donor, Leonardo A. dos Santos, Stephanie Douglas, Alan Du, Meredith Durbin, Dawn Erb, Daniel Faes, J. G. Fernández-Trincado, Anthony Foley, Sotiria Fotopoulou, Søren Frimann, Peter Frinchaboy, Rafael Garcia-Dias, Artur Gawryszczak, Elizabeth George, Sebastian Gonzalez, Karl Gordon, Nicholas Gorgone, Catherine Gosmeyer, Katie Grasha, Perry Greenfield, Rebekka Grellmann, James Guillochon, Mark Gurwell, Marcel Haas, Alex Hagen, Daryl Haggard, Tim Haines, Patrick Hall, Wojciech Hellwing, Edmund Christian Herenz, Samuel Hinton, Renee Hlozek, John Hoffman, Derek Holman, Benne Willem Holwerda, Anthony Horton, Cameron Hummels, Daniel Jacobs, Jens Juel Jensen, David Jones, Arna Karick, Luke Kelley, Matthew Kenworthy, Ben Kitchener, Dominik Klaes, Saul Kohn, Piotr Konorski, Coleman Krawczyk, Kyler Kuehn, Teet Kuutma, Michael T. Lam, Richard Lane, Jochen Liske, Diego Lopez-Camara, Katherine Mack, Sam Mangham, Qingqing Mao, David J. E. Marsh, Cecilia Mateu, Loïc Maurin, James McCormac, Ivelina Momcheva, Hektor Monteiro, Michael Mueller, Roberto Munoz, Rohan Naidu, Nicholas Nelson, Christian Nitschelm, Chris North, Juan Nunez-Iglesias, Sara Ogaz, Russell Owen, John Parejko, Vera Patrício, Joshua Pepper, Marshall Perrin, Timothy Pickering, Jennifer Piscionere, Richard Pogge, Radek Poleski, Alkistis Pourtsidou, Adrian M. Price-Whelan, Meredith L. Rawls, Shaun Read, Glen Rees, Hanno Rein, Thomas Rice, Signe Riemer-Sørensen, Naum Rusomarov, Sebastian F. Sanchez, Miguel Santander-García, Gal Sarid, William Schoenell, Aleks Scholz, Robert L. Schuhmann, William Schuster, Peter Scicluna, Marja Seidel, Lijing Shao, Pranav Sharma, Aleksandar Shulevski, David Shupe, Cristóbal Sifón, Brooke Simmons, Manodeep Sinha, Ian Skillen, Bjoern Soergel, Thomas Spriggs, Sundar Srinivasan, Abigail Stevens, Ole Streicher, Eric Suchyta, Joshua Tan, O. Grace Telford, Romain Thomas, Chiara Tonini, Grant Tremblay, Sarah Tuttle, Tanya Urrutia, Sam Vaughan, Miguel Verdugo, Alexander Wagner, Josh Walawender, Andrew Wetzel, Kyle Willett, Peter K. G. Williams, Guang Yang, Guangtun Zhu, Andrea Zonca
    Oct. 11, 2016 physics.soc-ph, astro-ph.IM
    The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.
  • The Magellanic Clouds provide the only laboratory to study the effect of metallicity and galaxy mass on molecular gas and star formation at high (~20 pc) resolution. We use the dust emission from HERITAGE Herschel data to map the molecular gas in the Magellanic Clouds, avoiding the known biases of CO emission as a tracer of H$_{2}$. Using our dust-based molecular gas estimates, we find molecular gas depletion times of ~0.4 Gyr in the LMC and ~0.6 SMC at 1 kpc scales. These depletion times fall within the range found for normal disk galaxies, but are shorter than the average value, which could be due to recent bursts in star formation. We find no evidence for a strong intrinsic dependence of the molecular gas depletion time on metallicity. We study the relationship between gas and star formation rate across a range in size scales from 20 pc to ~1 kpc, including how the scatter in molecular gas depletion time changes with size scale, and discuss the physical mechanisms driving the relationships. We compare the metallicity-dependent star formation models of Ostriker, McKee, and Leroy (2010) and Krumholz (2013) to our observations and find that they both predict the trend in the data, suggesting that the inclusion of a diffuse neutral medium is important at lower metallicity.
  • We define an appropriate problem for benchmarking dust emissivity calculations in the context of radiative transfer (RT) simulations, specifically including the emission from stochastically heated dust grains. Our aim is to provide a self-contained guide for implementors of such functionality, and to offer insights in the effects of the various approximations and heuristics implemented by the participating codes to accelerate the calculations. The benchmark problem definition includes the optical and calorimetric material properties, and the grain size distributions, for a typical astronomical dust mixture with silicate, graphite and PAH components; a series of analytically defined radiation fields to which the dust population is to be exposed; and instructions for the desired output. We process this problem using six RT codes participating in this benchmark effort, and compare the results to a reference solution computed with the publicly available dust emission code DustEM. The participating codes implement different heuristics to keep the calculation time at an acceptable level. We study the effects of these mechanisms on the calculated solutions, and report on the level of (dis)agreement between the participating codes. For all but the most extreme input fields, we find agreement within 10% across the important wavelength range from 3 to 1000 micron. We conclude that the relevant modules in RT codes can and do produce fairly consistent results for the emissivity spectra of stochastically heated dust grains.
  • We present a new method for deriving the stellar birth function (SBF) of resolved stellar populations. The SBF (stars born per unit mass, time, and metallicity) is the combination of the initial mass function (IMF), the star-formation history (SFH), and the metallicity distribution function (MDF). The framework of our analysis is that of Poisson Point Processes (PPPs), a class of statistical models suitable when dealing with points (stars) in a multidimensional space (the measurement space of multiple photometric bands). The theory of PPPs easily accommodates the modeling of measurement errors as well as that of incompleteness. Compared to most of the tools used to study resolved stellar populations, our method avoids binning stars in the color-magnitude diagram and uses the entirety of the information (i.e., the whole likelihood function) for each data point; the proper combination of the individual likelihoods allows the computation of the posterior probability for the global population parameters. This includes unknowns such as the IMF slope and combination of SFH and MDF, which are rarely solved for simultaneously in the literature, however entangled and correlated they might be. Our method also allows proper inclusion of nuisance parameters, such as distance and extinction distributions. The aim of this paper, is to assess the validity of this new approach under a range of assumptions, using only simulated data. Forthcoming work will show applications to real data. Although it has a broad scope of possible applications, we have developed this method to study multi-band HST observations of the Milky Way Bulge. Therefore we will focus on simulations with characteristics similar to those of the Galactic Bulge.
  • We present the analysis of supernova remnants (SNRs) in the Large Magellanic Cloud (LMC) and their influence on the environment at far-infrared (FIR) and submillimeter wavelengths. We use new observations obtained with the {\it Herschel} Space Observatory and archival data obtained with the {\it Spitzer} Space Telescope, to make the first FIR atlas of these objects. The SNRs are not clearly discernible at FIR wavelengths, however their influence becomes apparent in maps of dust mass and dust temperature, which we constructed by fitting a modified black-body to the observed spectral energy distribution in each sightline. Most of the dust that is seen is pre-existing interstellar dust in which SNRs leave imprints. The temperature maps clearly reveal SNRs heating surrounding dust, while the mass maps indicate the removal of 3.7$^{+7.5}_{-2.5}$ M$_{\odot}$ of dust per SNR. This agrees with the calculations by others that significant amounts of dust are sputtered by SNRs. Under the assumption that dust is sputtered and not merely pushed away, we estimate a dust destruction rate in the LMC of $0.037^{+0.075}_{-0.025}$ M$_\odot$ yr$^{-1}$ due to SNRs, yielding an average lifetime for interstellar dust of $2^{+4.0}_{-1.3}\times10^7$ yr. We conclude that sputtering of dust by SNRs may be an important ingredient in models of galactic evolution, that supernovae may destroy more dust than they produce, and that they therefore may not be net producers of long lived dust in galaxies.
  • The spatial variations of the gas-to-dust ratio (GDR) provide constraints on the chemical evolution and lifecycle of dust in galaxies. We examine the relation between dust and gas at 10-50 pc resolution in the Large and Small Magellanic Clouds (LMC and SMC) based on Herschel far-infrared (FIR), H I 21 cm, CO, and Halpha observations. In the diffuse atomic ISM, we derive the gas-to-dust ratio as the slope of the dust-gas relation and find gas-to-dust ratios of 380+250-130 in the LMC, and 1200+1600-420 in the SMC, not including helium. The atomic-to-molecular transition is located at dust surface densities of 0.05 Mo pc-2 in the LMC and 0.03 Mo pc-2 in the SMC, corresponding to AV ~ 0.4 and 0.2, respectively. We investigate the range of CO-to-H2 conversion factor to best account for all the molecular gas in the beam of the observations, and find upper limits on XCO to be 6x1020 cm-2 K-1 km-1 s in the LMC (Z=0.5Zo) at 15 pc resolution, and 4x 1021 cm-2 K-1 km-1 s in the SMC (Z=0.2Zo) at 45 pc resolution. In the LMC, the slope of the dust-gas relation in the dense ISM is lower than in the diffuse ISM by a factor ~2, even after accounting for the effects of CO-dark H2 in the translucent envelopes of molecular clouds. Coagulation of dust grains and the subsequent dust emissivity increase in molecular clouds, and/or accretion of gas-phase metals onto dust grains, and the subsequent dust abundance (dust-to-gas ratio) increase in molecular clouds could explain the observations. In the SMC, variations in the dust-gas slope caused by coagulation or accretion are degenerate with the effects of CO-dark H2. Within the expected 5--20 times Galactic XCO range, the dust-gas slope can be either constant or decrease by a factor of several across ISM phases. Further modeling and observations are required to break the degeneracy between dust grain coagulation, accretion, and CO-dark H2.
  • We have measured stellar photometry with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) and Advanced Camera for Surveys (ACS) in near ultraviolet (F275W, F336W), optical (F475W, F814W), and near infrared (F110W, F160W) bands for 117 million resolved stars in M31. As part of the Panchromatic Hubble Andromeda Treasury (PHAT) survey, we measured photometry with simultaneous point spread function fitting across all bands and at all source positions after precise astrometric image alignment (<5-10 milliarcsecond accuracy). In the outer disk, the photometry reaches a completeness-limited depth of F475W~28, while in the crowded, high surface brightness bulge, the photometry reaches F475W~25. We find that simultaneous photometry and optimized measurement parameters significantly increase the detection limit of the lowest resolution filters (WFC3/IR) providing color-magnitude diagrams that are up to 2.5 magnitudes deeper when compared with color-magnitude diagrams from WFC3/IR photometry alone. We present extensive analysis of the data quality including comparisons of luminosity functions and repeat measurements, and we use artificial star tests to quantify photometric completeness, uncertainties and biases. We find that largest sources of systematic error in the photometry are due to spatial variations in the point spread function models and charge transfer efficiency corrections. This stellar catalog is the largest ever produced for equidistant sources, and is publicly available for download by the community.
  • We search Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) and Wide Field Camera 3 (WFC3) broadband imaging data from the Panchromatic Hubble Andromeda Treasury (PHAT) survey to identify detections of cataloged planetary nebulae (PNe). Of the 711 PNe currently in the literature within the PHAT footprint, we find 467 detected in the broadband. For these 467 we are able to refine their astrometric accuracy from ~0."3 to 0."05. Using the resolution of HST, we are able to show that 152 objects currently in the catalogs are definitively not PNe, and we show that 32 objects thought to be extended in ground-based images are actually point-like and therefore good PN candidates. We also find one PN candidate that is marginally resolved. If this is a PN, it is up to 0.7 pc in diameter. With our new photometric data, we develop a method of measuring the level of excitation in individual PNe by comparing broadband and narrowband imaging and describe the effects of excitation on a PN's photometric signature. Using the photometric properties of the known PNe in the PHAT catalogs, we search for more PN, but do not find any new candidates, suggesting that ground-based emission-line surveys are complete in the PHAT footprint to F475W $\simeq$ 24.
  • In this paper we present a catalogue of Giant Molecular Clouds (GMCs) in the Andromeda (M31) galaxy extracted from the Hershel Exploitation of Local Galaxy Andromeda (HELGA) dataset. GMCs are identified from the Herschel maps using a hierarchical source extraction algorithm. We present the results of this new catalogue and characterise the spatial distribution and spectral energy properties of its clouds based on the radial dust/gas properties found by Smith et al (2012). 236 GMCs in the mass range 10^4-10^7 M_sol are identified, their cumulative mass distribution is found to be proportional to M^-1.45 in agreement with earlier studies. The GMCs appear to follow the same cloud mass to L_CO correlation observed in the Milky Way. However, comparison between this catalogue and interferometry studies also shows that the GMCs are substructured below the Herschel resolution limit suggesting that we are observing associations of GMCs. Following Gordon et al. (2006), we study the spatial structure of M31 by splitting the observed structure into a set of spiral arms and offset rings. We fit radii of 10.5 and 15.5 kpc to the two most prominent rings. We then fit a logarithmic spiral with a pitch angle of 8.9 deg to the GMCs not associated with either ring. Lastly, we comment upon the effects of deprojection on our results and investigate the effect different models for M31's inclination will have upon the projection of an unperturbed spiral arm system.
  • We investigate the far-infrared (IR) dust emission for 20 local star forming galaxies from the Key Insights on Nearby Galaxies: A Far-IR Survey with Herschel (KINGFISH) sample. We model the far-IR/submillimeter spectral energy distribution (SED) using images from Spitzer Space Telescope and Herschel Space Observatory. We calculate the cold dust temperature (T(cold)) and emissivity (beta) on a pixel by pixel basis (where each pixel ranges from 0.1-3 kpc^2) using a two temperature modified blackbody fitting routine. Our fitting method allows us to investigate the resolved nature of temperature and emissivity variations by modeling from the galaxy centers to the outskirts (physical scales of ~15-50 kpc, depending on the size of the galaxy). We fit each SED in two ways: (1) fit T(cold) and beta simultaneously, (2) hold beta constant and fit T(cold). We compare T(cold) and beta with star formation rates (calculated from L(Halpha) and L(24)), the luminosity of the old stellar population (traced through L(3.6), and the dust mass surface density (traced by 500 micron luminosity, L(500)). We find a significant trend between SFR/L(500) and T(cold), implying that the flux of hard UV photons relative to the amount of dust is significantly contributing to the heating of the cold, or diffuse, dust component. We also see a trend between L(3.6)/L(500) and beta, indicating that the old stellar population contributes to the heating at far-IR/submillimeter wavelengths. Finally, we find that when beta is held constant, T(cold) exhibits a strongly decreasing radial trend, illustrating that the shape of the far-IR SED is changing radially through a galaxy, thus confirming on a sample almost double in size the trends observed in Galametz et al. (2012).
  • We measure the extinction curve in the central 200 pc of M31 at mid-ultraviolet to near-infrared wavelengths (from 1928A to 1.5{\mu}m), using Swift/UVOT and HST WFC3/ACS observations in thirteen bands. Taking advantage of the high angular resolution of the HST WFC3 and ACS detectors, we develop a method to simultaneously determine the relative extinction and the fraction of obscured starlight for five dusty complexes located in the circumnuclear region. The extinction curves of these clumps (RV =2.4-2.5) are steeper than the average Galactic one (RV =3.1), but are similar to optical and near-infrared curves recently measured toward the Galactic Bulge (RV~2.5). This similarity suggests that steep extinction curves may be common in the inner bulge of galaxies. In the ultraviolet, the extinction curves of these clumps are also unusual. We find that one dusty clump (size <2 pc) exhibits a strong UV bump (extinction at 2175A), more than three standard deviation higher than that predicted by common models. Although the high stellar metallicity of the M31 bulge indicates that there are sufficient carbon and silicon to produce large dust grains, the grains may have been destroyed by supernova explosions or past activity of the central super-massive black hole, resulting in the observed steepened extinction curve.
  • 'Normal' galaxies observed at z>6, when the Universe was <1 billion years old, thus far show no evidence of the cold dust that accompanies star formation in the local Universe, where the dust-to-gas mass ratio is 1%. A prototypical example is 'Himiko' (z=6.6), which a mere 840 Myr after the Big Bang is forming stars at a rate of 30-100 Msun/yr, yielding a mass assembly time M^{star}/SFR 150x10^6 yr. Himiko is estimated to have a low fraction (2-3% of the Solar value) of elements heavier than helium (metallicity), and although its gas mass cannot be asserted at this time its dust-to-stellar mass ratio is constrained to be <0.05%. The local galaxy I Zw 18, with a metallicity 4% solar and forming stars less rapidly than Himiko but still vigorously for its mass (M^{star}/SFR 1.6x10^9 yr), is also very dust deficient and perhaps one of the best analogues of primitive galaxies accessible to detailed study. Here we report observations of dust emission from I Zw 18 from which we determine its dust mass to be 450-1800 Msun, yielding a dust-to-stellar mass ratio \approx 10^{-6}-10^{-5} and a dust-to-gas mass ratio 3.2-13x10^{-6}. If I Zw 18 is a reasonable analog of Himiko, then Himiko's dust mass is \approx 50,000 Msun, a factor of 100 below the current upper limit. These numbers are considerably uncertain, but if most high-z galaxies are more like Himiko than like the quasar host SDSS J114816.64+525150.3, then the prospects for detecting the gas and dust in them are much poorer than hitherto anticipated.
  • Motived by recent ISM studies, we present high quality MMT and Gemini spectroscopic observations of H II regions in the nearby spiral galaxies NGC 628 and NGC 2403 in order to measure their chemical abundance gradients. Using long-slit and multi-object mask optical spectroscopy, we obtained measurements of the temperature sensitive auroral lines [O III] {\lambda}4363 and/or [N II] {\lambda}5755 at a strength of 4{\sigma} or greater in 11 H II regions in NGC 628 and 7 regions in NGC 2403. These observations allow us, for the first time, to derive an oxygen abundance gradient in NGC 628 based solely on "direct" oxygen abundances of H II regions: 12 + log(O/H) = (8.43+/-0.03) + (-0.017+/-0.002) x Rg (dex/kpc), with a dispersion in log(O/H) of {\sigma} = 0.10 dex, from 14 regions with a radial coverage of ~2-19 kpc. This is a significantly shallower slope than found by previous "strong-line" abundance studies. In NGC 2403, we derive an oxygen abundance gradient of 12 + log(O/H) = (8.48+/-0.04) + (-0.032+/-0.007) x Rg (dex/kpc), with a dispersion in log(O/H) of {\sigma} = 0.07 dex, from 7 H II with a radial coverage of ~1-10 kpc. Additionally, we measure the N, S, Ne, and Ar abundances. We find the N/O ratio decreases with increasing radius for the inner disk, but reaches a plateau past R25 in NGC 628. NGC 2403 also has a negative N/O gradient with radius, but we do not sample the outer disk of the galaxy past R25 and so do not see evidence for a plateau. This bi-modal pattern measured for NGC 628 indicates dominant contributions from secondary nitrogen inside of the R25 transition and dominantly primary nitrogen farther out. As expected for {\alpha}-process elements, S/O, Ne/O, and Ar/O are consistent with constant values over a range in oxygen abundance.
  • Cosmic dust is present in many astrophysical objects, and recent observations across the electromagnetic spectrum show that the dust distribution is often strongly three-dimensional (3D). Dust grains are effective in absorbing and scattering ultraviolet (UV)/optical radiation, and they re-emit the absorbed energy at infrared wavelengths. Understanding the intrinsic properties of these objects, including the dust itself, therefore requires 3D dust radiative transfer (RT) calculations. Unfortunately, the 3D dust RT problem is nonlocal and nonlinear, which makes it one of the hardest challenges in computational astrophysics. Nevertheless, significant progress has been made in the past decade, with an increasing number of codes capable of dealing with the complete 3D dust RT problem. We discuss the complexity of this problem, the two most successful solution techniques [ray-tracing (RayT) and Monte Carlo (MC)], and the state of the art in modeling observational data using 3D dust RT codes. We end with an outlook on the bright future of this field.
  • We study the emission by dust and stars in the Large and Small Magellanic Clouds, a pair of low-metallicity nearby galaxies, as traced by their spatially resolved spectral energy distributions (SEDs). This project combines Herschel Space Observatory PACS and SPIRE far-infrared photometry with other data at infrared and optical wavelengths. We build maps of dust and stellar luminosity and mass of both Magellanic Clouds, and analyze the spatial distribution of dust/stellar luminosity and mass ratios. These ratios vary considerably throughout the galaxies, generally between the range $0.01\leq L_{\rm dust}/L_\ast\leq 0.6$ and $10^{-4}\leq M_{\rm dust}/M_\ast\leq 4\times10^{-3}$. We observe that the dust/stellar ratios depend on the interstellar medium (ISM) environment, such as the distance from currently or previously star-forming regions, and on the intensity of the interstellar radiation field (ISRF). In addition, we construct star formation rate (SFR) maps, and find that the SFR is correlated with the dust/stellar luminosity and dust temperature in both galaxies, demonstrating the relation between star formation, dust emission and heating, though these correlations exhibit substantial scatter.
  • We compare atomic gas, molecular gas, and the recent star formation rate (SFR) inferred from H-alpha in the Small Magellanic Cloud (SMC). By using infrared dust emission and local dust-to-gas ratios, we construct a map of molecular gas that is independent of CO emission. This allows us to disentangle conversion factor effects from the impact of metallicity on the formation and star formation efficiency of molecular gas. On scales of 200 pc to 1 kpc we find a characteristic molecular gas depletion time of ~1.6 Gyr, similar to that observed in the molecule-rich parts of large spiral galaxies on similar spatial scales. This depletion time shortens on much larger scales to ~0.6 Gyr because of the presence of a diffuse H-alpha component, and lengthens on much smaller scales to ~7.5 Gyr because the H-alpha and H2 distributions differ in detail. We estimate the systematic uncertainties in our measurement to be a factor of 2-3. We suggest that the impact of metallicity on the physics of star formation in molecular gas has at most this magnitude. The relation between SFR and neutral (H2+HI) gas surface density is steep, with a power-law index ~2.2+/-0.1, similar to that observed in the outer disks of large spiral galaxies. At a fixed total gas surface density the SMC has a 5-10 times lower molecular gas fraction (and star formation rate) than large spiral galaxies. We explore the ability of the recent models by Krumholz et al. (2009) and Ostriker et al. (2010) to reproduce our observations. We find that to explain our data at all spatial scales requires a low fraction of cold, gravitationally-bound gas in the SMC. We explore a combined model that incorporates both large scale thermal and dynamical equilibrium and cloud-scale photodissociation region structure and find that it reproduces our data well, as well as predicting a fraction of cold atomic gas very similar to that observed in the SMC.
  • We estimate the conversion factor relating CO emission to H2 mass, alpha_CO, in five Local Group galaxies that span approximately an order of magnitude in metallicity - M31, M 33, the Large Magellanic Cloud (LMC), NGC 6822, and the Small Magellanic Cloud (SMC). We model the dust mass along the line of sight from infrared (IR) emission and then solve for the alpha_CO that best allows a single gas-to-dust ratio (delta_GDR) to describe each system. This approach remains sensitive to CO-dark envelopes of H2 surrounding molecular clouds. In M 31, M 33, and the LMC we find alpha_CO \approx 3-9 M_sun pc^-2 (K km s^-1)^-1, consistent with the Milky Way value within the uncertainties. The two lowest metallicity galaxies in our sample, NGC 6822 and the SMC (12 + log(O/H) \approx 8.2 and 8.0), exhibit a much higher alpha_CO. Our best estimates are \alpha_NGC6822 \approx 30 M_sun/pc^-2 (K km s^-1)^-1 and \alpha_SMC \approx 70 M_sun/pc^-2 (K km s-1)-1. These results are consistent with the conversion factor becoming CO a strong function of metallicity around 12 + log(O/H) \sim 8.4 - 8.2. We favor an interpretation where decreased dust-shielding leads to the dominance of CO-free envelopes around molecular clouds below this metallicity.
  • Fundamental parameters and time-evolution of mass loss are investigated for post-main-sequence stars in the Galactic globular cluster 47 Tucanae (NGC 104). This is accomplished by fitting spectral energy distributions (SEDs) to existing optical and infrared photometry and spectroscopy, to produce a true Hertzsprung--Russell diagram. We confirm the cluster's distance as 4611 (+213, -200) pc and age as 12 +/- 1 Gyr. Horizontal branch models appear to confirm that no more RGB mass loss occurs in 47 Tuc than in the more-metal-poor omega Centauri, though difficulties arise due to inconsistencies between the models. Using our SEDs, we identify those stars which exhibit infrared excess, finding excess only among the brightest giants: dusty mass loss begins at a luminosity of ~ 1000 Lsun, becoming ubiquitous above 2000 Lsun. Recent claims of dust production around lower-luminosity giants cannot be reproduced, despite using the same archival Spitzer imagery.
  • The amount of molecular gas is a key for understanding the future star formation in a galaxy. Because H2 is difficult to observe directly in dense and cold clouds, tracers like CO are used. However, at low metallicities especially, CO only traces the shielded interiors of the clouds. mm dust emission can be used as a tracer to unveil the total dense gas masses. The comparison of masses deduced from the continuum SIMBA 1.2 mm emission and virial masses in a sample of giant molecular clouds (GMCs), in the SW region of the Small Magellanic Cloud (SMC), showed a discrepancy that is in need of an explanation. This study aims at better assessing possible uncertainties on the dust emission observed in the sample of GMCs from the SMC and focuses on the densest parts of the GMCs where CO is detected. New observations were obtained with the LABOCA camera on the APEX telescope. All GMCs previously observed in CO are detected and their emission at 870microns is compared to ancillary data. The different contributions to the sub-mmm emission are estimated, as well as dust properties, in order to deduce molecular cloud masses precisely. The (sub-)mm emission observed in the GMCs in the SW region of the SMC is dominated by dust emission and masses are deduced for the part of each cloud where CO is detected and compared to the virial masses. The mass discrepancy between both methods is confirmed at 870microns with the LABOCA observations: the virial masses are on average 4 times smaller than the masses of dense gas deduced from dust emission, contrary to what is observed for equivalent clouds in our Galaxy. The origin of this mass discrepancy in the SMC remains unkown. The direct interpretation of this effect is that the CO linewidth used to compute virial masses do not measure the full velocity distribution of the gas. Geometrical effects and uncertainties on the dust properties are also discussed.
  • Spitzer and IRAS observations of the LMC suggest an excess of FIR emission with respect to the gas surface density traced by 12CO and HI 21 cm emission lines. This "FIR excess" is noticeable near molecular clouds in the LMC, and has usually been interpreted as the presence of a self-shielded H2 component not traced by CO molecular clouds' envelopes. Based on Herschel observations, we examine the correlation between gas and dust at higher resolution than previously achieved. We consider three additional causes for the FIR excess: X factor, FIR dust emissivity, and gas-to-dust ratio variations between the diffuse and dense phases of the ISM. We examine the structure of NT80 and NT71, two molecular clouds detected in the NANTEN 12CO survey of the LMC. Dust surface density maps are derived from the HERITAGE data. The gas phase is traced by MAGMA 12CO and ATCA HI 21 cm observations of the LMC. The dust emissivity, gas-to-dust ratio, and X factor required to match the dust and gas surface densities are derived, and their correlations with the dust surface density are examined. The dust surface density is spatially correlated with the atomic and molecular gas phases. The dust temperature is consistently lower in the dense phase of the ISM than in the diffuse phase. We confirm variations in the ratio of FIR emission to gas surface density derived from HI and CO observations. There is an excess of FIR emission, spatially correlated with regions of intermediate HI and dust surface densities (Av = 1-2), and little or no CO. While there is no significant trend in the dust emissivity or gas-to-dust ratio with dust surface density, the X factor is enhanced at Av = 1-2. We conclude that H2 envelopes not traced by CO and X factor variations close to the CO boundary are more likely to cause these deviations between FIR emission and gas surface density than gas-to-dust ratio or emissivity variations.
  • Using the data obtained with the Spitzer Space telescope as part of the Surveying the Agents of a Galaxy's Evolution (SAGE) legacy survey, we have studied the variations of the dust composition and abundance across the Large Magellanic Cloud (LMC). Such variations are expected, as the explosive events which have lead to the formation of the many HI shells observed should have affected the dust properties. Using a model and comparing with a reference spectral energy distribution from our Galaxy, we deduce the relative abundance variations of small dust grains across the LMC. We examined the infrared color ratios as well as the relative abundances of very small grains (VSGs) and polycyclic aromatic hydrocarbons (PAHs) relative to the big grain (BG) abundance. Results show that each dust component could have different origins or evolution in the interstellar medium (ISM). The VSG abundance traces the star formation activity and could result from shattering of larger grains, whereas the PAH abundance increases around molecular clouds as well as in the stellar bar, where they could have been injected into the ISM during mass loss from old stars.
  • We explore the stellar population of M31 in a Spitzer Space Telescope survey utilizing IRAC and MIPS observations. Red supergiants are the brightest objects seen in the infrared; they are a prominent evolutionary phase. Due to their circumstellar envelopes, many of these radiate the bulk of their luminosity at IRAC wavelengths and do not stand out in the near infrared or optically. Going fainter, we see large numbers of luminous asymptotic giant branch (AGB) stars, many of which are known long period variables. Relative to M33, the AGB carbon star population of M31 appears sparse, but this needs to be spectroscopically confirmed.
  • Variable A in M33 is a member of a rare class of highly luminous, evolved stars near the upper luminosity boundary that show sudden and dramatic shifts in apparent temperature due to the formation of optically thick winds in high mass loss episodes. Recent optical and infrared spectroscopy and imaging reveal that its ``eruption'' begun in $\sim$1950 has ended, {\it lasting $\approx$ 45 yrs}. Our current observations show major changes in its wind from a cool, dense envelope to a much warmer state surrounded by low density gas with rare emission lines of Ca II, [Ca II] and K I. Its spectral energy distribution has unexpectedly changed, especially at the long wavelengths, with a significant decrease in its apparent flux, while the star remains optically obscured. We conclude that much of its radiation is now escaping out of our line of sight. We attribute this to the changing structure and distribution of its circumstellar ejecta corresponding to the altered state of its wind as the star recovers from a high mass loss event.