Search

JWST/NIRCam obtained high angular-resolution (0.05-0.1''), deep near-infrared 1--5 micron imaging of Supernova (SN) 1987A taken 35 years after the explosion. In the NIRCam images, we identify: 1) faint H2 crescents, which are emissions located between the ejecta and the equatorial ring, 2) a bar, which is a substructure of the ejecta, and 3) the bright 3-5 micron continuum emission exterior to the equatorial ring. The emission of the remnant in the NIRCam 1-2.3 micron images is mostly due to line emission, which is mostly emitted in the ejecta and in the hot spots within the equatorial ring. In contrast, the NIRCam 3-5 micron images are dominated by continuum emission. In the ejecta, the continuum is due to dust, obscuring the centre of the ejecta. In contrast, in the ring and exterior to the ring, synchrotron emission contributes a substantial fraction to the continuum. Dust emission contributes to the continuum at outer spots and diffuse emission exterior to the ring, but little within the ring. This shows that dust cooling and destruction time scales are shorter than the synchrotron cooling time scale, and the time scale of hydrogen recombination in the ring is even longer than the synchrotron cooling time scale. With the advent of high sensitivity and high angular resolution images provided by JWST/NIRCam, our observations of SN 1987A demonstrate that NIRCam opens up a window to study particle-acceleration and shock physics in unprecedented details, probed by near-infrared synchrotron emission, building a precise picture of how a SN evolves., Comment: Accepted for publication in MNRAS; 18 pages

The nearby Supernova 1987A was accompanied by a burst of neutrino emission, which indicates that a compact object (a neutron star or black hole) was formed in the explosion. There has been no direct observation of this compact object. In this work, we observe the supernova remnant with JWST spectroscopy finding narrow infrared emission lines of argon and sulphur. The line emission is spatially unresolved and blueshifted in velocity relative to the supernova rest frame. We interpret the lines as gas illuminated by a source of ionizing photons located close to the center of the expanding ejecta. Photoionization models show that the line ratios are consistent with ionization by a cooling neutron star or pulsar wind nebula. The velocity shift could be evidence for a neutron star natal kick., Comment: Authors version of manuscript published in Science on 22 Feb 2024

There exist very few mid-infrared (IR) observations of supernovae (SNe) in general. Therefore, SN 1987A, the closest visible SN in 400 years, gives us the opportunity to explore the mid-IR properties of SNe, the dust in their ejecta and surrounding medium, and to witness the birth of a SN remnant (SNR). The James Webb Space Telescope (JWST), with its high spatial resolution and extreme sensitivity, gives a new view on these issues. We report on the first imaging observations obtained with the Mid-InfraRed Instrument (MIRI). We build temperature maps and discuss the morphology of the nascent SNR. Our results show that the temperatures in the equatorial ring (ER) are quite non-uniform. This could be due to dust destruction in some parts of the ring, as had been assumed in some previous works. We show that the IR emission extends beyond the ER, illustrating the fact that the shock wave has now passed through this ring to affect the circumstellar medium on a larger scale. Finally, while sub-mm Atacama Large Millimeter Array (ALMA) observations have hinted at the location of the compact remnant of SN 1987A, we note that our MIRI data have found no such evidence., Comment: 19 pages, 19 figures, 2 tables; Accepted for publication in the Astrophysical Journal (February 2, 2024)

PRIMA (The PRobe for-Infrared Mission for Astrophysics) is a concept for a far-infrared (IR) observatory. PRIMA features a cryogenically cooled 1.8 m diameter telescope and is designed to carry two science instruments enabling ultra-high sensitivity imaging and spectroscopic studies in the 24 to 235 microns wavelength range. The resulting observatory is a powerful survey and discovery machine, with mapping speeds better by 2 - 4 orders of magnitude with respect to its far-IR predecessors. The bulk of the observing time on PRIMA should be made available to the community through a General Observer (GO) program offering 75% of the mission time over 5 years. In March 2023, the international astronomy community was encouraged to prepare authored contributions articulating scientific cases that are enabled by the telescope massive sensitivity advance and broad spectral coverage, and that could be performed within the context of GO program. This document, the PRIMA General Observer Science Book, is the edited collection of the 76 received contributions., Comment: A. Moullet, T. Kataria, D. Lis, S. Unwin, Y. Hasegawa, E. Mills, C. Battersby, A. Roc, M. Meixner are the editors of the PRIMA General Observer Science Book. The book compiles 76 authored contributions. 399 pages

Supernova (SN) 1987A is the nearest supernova in $\sim$400 years. Using the {\em JWST} MIRI Medium Resolution Spectrograph, we spatially resolved the ejecta, equatorial ring (ER) and outer rings in the mid-infrared 12,927 days after the explosion. The spectra are rich in line and dust continuum emission, both in the ejecta and the ring. Broad emission lines (280-380~km~s$^{-1}$ FWHM) seen from all singly-ionized species originate from the expanding ER, with properties consistent with dense post-shock cooling gas. Narrower emission lines (100-170~km~s$^{-1}$ FWHM) are seen from species originating from a more extended lower-density component whose high ionization may have been produced by shocks progressing through the ER, or by the UV radiation pulse associated with the original supernova event. The asymmetric east-west dust emission in the ER has continued to fade, with constant temperature, signifying a reduction in dust mass. Small grains in the ER are preferentially destroyed, with larger grains from the progenitor surviving the transition from SN into SNR. The ER is fit with a single set of optical constants, eliminating the need for a secondary featureless hot dust component. We find several broad ejecta emission lines from [Ne~{\sc ii}], [Ar~{\sc ii}], [Fe~{\sc ii}], and [Ni~{\sc ii}]. With the exception of [Fe~{\sc ii}]~25.99$\mu$m, these all originate from the ejecta close to the ring and are likely being excited by X-rays from the interaction. The [Fe~{\sc ii}]~5.34$\mu$m to 25.99$\mu$m line ratio indicates a temperature of only a few hundred K in the inner core, consistent with being powered by ${}^{44}$Ti decay., Comment: 27 pages, 16 figures, 4 tables. Accepted ApJ

We present initial results from JWST NIRSpec integral field unit observations of the nearby Supernova (SN) 1987A. The observations provide the first spatially-resolved spectroscopy of the ejecta and equatorial ring (ER) over the 1-5 \mu m range. We construct 3D emissivity maps of the [Fe I] 1.443 \mu m line from the inner ejecta and the He I 1.083 \mu m line from the reverse shock (RS), where the former probes the explosion geometry and the latter traces the structure of the circumstellar medium. We also present a model for the integrated spectrum of the ejecta. The [Fe I] 3D map reveals a highly-asymmetric morphology resembling a broken dipole, dominated by two large clumps with velocities of ~2300 km/s. We also find evidence that the Fe-rich inner ejecta have started to interact with the RS. The RS surface traced by the He I line extends from just inside the ER to higher latitudes on both sides of the ER with a half-opening angle ~45 degrees, forming a bubble-like structure. The spectral model for the ejecta allows us to identify the many emission lines, including numerous H_2 lines. We find that the H_2 is most likely excited by far-UV emission, while the metal lines ratios are consistent with a combination of collisional excitation and recombination in the low-temperature ejecta. We also find several high-ionization coronal lines from the ER, requiring a temperature > 2 \times 10^6 K., Comment: Accepted for publication in ApJL

We present new Hubble Space Telescope WFC3/IR medium-band photometry of the compact elliptical galaxy M32, chemically resolving its thermally pulsating asymptotic giant branch stars. We find 2829 M-type stars and 57 C stars. The carbon stars are likely contaminants from M31. If carbon stars are present in M32 they are so in very low numbers. The uncorrected C/M ratio is 0.020 $\pm$ 0.003; this drops to less than 0.007 after taking into account contamination from M31. As the mean metallicity of M32 is just below solar, this low ratio of C to M stars is unlikely due to a metallicity ceiling for the formation of carbon stars. Instead, the age of the AGB population is likely to be the primary factor. The ratio of AGB to RGB stars in M32 is similar to that of the inner disc of M31 which contain stars that formed 1.5-4 Gyr ago. If the M32 population is at the older end of this age then its lack of C-stars may be consistent with a narrow mass range for carbon star formation predicted by some stellar evolution models. Applying our chemical classifications to the dusty variable stars identified with {\em Spitzer}, we find that the x-AGB candidates identified with Spitzer are predominately M-type stars. This substantially increases the lower limit to the cumulative dust-production rate in M32 to $>$ 1.20 $\times 10^{-5}$ ${\rm M}_{\odot} \, {\rm yr}^{-1}$., Comment: 10 pages, 7 figures, accepted MNRAS

We present medium-resolution (R $\sim$ 4000) YJ, H \& K band spectroscopy of candidate young stellar objects (YSOs) in NGC~346, the most active star-formation region in the metal-poor (Z = 1/5 Z$_{\sun}$) Small Magellanic Cloud. The spectra were obtained with the KMOS (K-Band Multi Object Spectrograph) integral field instrument on the Very Large Telescope. From our initial sample of 18 candidate high-mass YSOs previously identified from mid-IR photometry and radiative transfer model fits to their spectral energy distributions, approximately half were resolved into multiple components by our integral-field data. In total, we detect 30 continuum sources and extract reliable spectra for 12 of these objects. The spectra show various features including hydrogen recombination lines, and lines from H$_2$, He~{\sc i} and [Fe~{\sc ii}], which are indicative of accretion, discs and outflowing material in massive YSOs. We spectroscopically confirm the youthful nature of nine YSO candidates and identify two others as OB stars. All of the confirmed YSOs have Br$\gamma$ in emission, but no emission is seen from the CO bandhead, despite other disc tracers present in the spectra. He\,{\sc i}~1.083 $\mu$m emission is also detected at appreciably higher rates than for the Galaxy., Comment: 15 pages, 6 figures. Accepted for publication MNRAS

12CO and 13CO lines, as well as a mm-wave continuum, have been observed for a sample of 22 OH/IR stars in directions within 2 degrees of the Galactic Centre. Photometry data have been gathered from the literature to construct SEDs and to determine pulsational variability. Radiative transfer models have been used to interpret the data. All stars in the sample were detected in at least one CO line, and 8 objects were detected in 324 GHz continuum. Based on luminosity criteria, the sample is divided into 17 objects that most likely lie within the inner Galactic Bulge, and 5 objects that are most likely foreground objects. The median luminosity of the inner-Galactic-Bulge objects, 5600 Lsun, corresponds to an initial mass in the range 1.2-1.6 Msun, indicating that these OH/IR stars descend from solar-type stars. The objects in this sub-sample are further divided into two classes based on their SED characteristics: 11 objects have SEDs that are well matched by models invoking dust envelopes extending from a few stellar radii and outwards, while 6 objects are better modelled as having detached dust envelopes with inner radii in the range 200-600 au and warmer central stars. The former objects have periodic variability, while the latter objects are predominantly non-periodic. The median gas-mass-loss rate, gas terminal expansion velocity, gas-to-dust mass ratio, and circumstellar 12CO/13CO abundance ratio have been estimated to be 2x10{-5} Msun/yr, 18 km/s, 200 (excluding the sources with detached dust envelopes, which show markedly lower gas-to-dust ratios), and 5, respectively, for the inner-Galactic-Bulge objects. The inner-Galactic-Bulge OH/IR stars studied here constitute an excellent sample of equidistant objects for the purpose of understanding the evolution of the mass-loss-rate characteristics at the tip of the AGB., Comment: 32 pages, 17 figures

Variable stars in the compact elliptical galaxy M32 are identified, using three epochs of photometry from the Spitzer Space Telescope at 3.6 and 4.5 $\mu$m, separated by 32 to 381 days. We present a high-fidelity catalogue of sources detected in multiple epochs at both 3.6 and 4.5 $\mu$m, which we analysed for stellar variability using a joint probability error-weighted flux difference. Of these, 83 stars are identified as candidate large-amplitude, long-period variables, with 28 considered high-confidence variables. The majority of the variable stars are classified as asymptotic giant branch star candidates using colour-magnitude diagrams. We find no evidence supporting a younger, infrared-bright stellar population in our M32 field., Comment: 11 pages, 8 figures. Accepted for publication in MNRAS

Young Massive Clusters (YMCs) and Super Star Clusters (SSCs) represent an extreme mode of star formation. Far-infrared imaging of the Magellanic Clouds has identified one potential embedded SSC, HSO BMHERICC J72.971176-69.391112 (HH in short), in the south-west outskirts of the Large Magellanic Cloud. We present Gemini Flamingos 2 and GSAOI near-infrared imaging of a 3'x3' region around HH in order to characterize the stellar content of the cluster. The stellar content is probed down to 1.5 Msun. We find substantial dust extinction across the cluster region, extending up to A_K of 3. Deeply embedded stars are associated with ALMA-detected molecular gas suggesting that star formation is ongoing. The high spatial resolution of the GSAOI data allows identification of the central massive object associated with the 13CO ALMA observations and to detect fainter low-mass stars around the H30alpha ALMA source. The morphology of the molecular gas and the nebulosity from adjacent star formation suggest they have interacted covering a region of several pc. The total stellar content in the cluster is estimated from the intermediate- and high-mass stellar content to be at least 10000 Msun, less than R136 with up to 100 000 Msun within 4.7 pc radius, but places it in the regime of a super star cluster. Based on the extinction determination of individual stars we estimate a molecular gas mass in the vicinity of HH of 6600 Msun, suggesting more star formation can be expected., Comment: Acceepted, AJ

The Origins Space Telescope (Origins) is one of four science and technology definition studies selected by National Aeronautics and Space Administration (NASA) in preparation of the 2020 Astronomy and Astrophysics Decadal survey in the US. Origins will trace the history of our origins from the time dust and heavy elements permanently altered the cosmic landscape to present-day life. It is designed to answer three major science questions: How do galaxies form stars, make metals, and grow their central supermassive black holes from reionization? How do the conditions for habitability develop during the process of planet formation? Do planets orbiting M-dwarf stars support life? Origins operates at mid- to far-infrared wavelengths from ~2.8 to 588 {\mu}m, is more than 1000 times more sensitive than prior far-IR missions due to its cold (~4.5 K) aperture and state-of-the-art instruments., Comment: 24 pages, 12 figures, 5 tables. Submitted as White paper to Voyage 2050 Long-term planning of the ESA Science Programme. arXiv admin note: substantial text overlap with arXiv:1912.06213

The Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and large-area extragalactic fields. Origins operates in the wavelength range 2.8 to 588 microns and is 1000 times more sensitive than its predecessors due to its large, cold (4.5 K) telescope and advanced instruments. Origins was one of four large missions studied by the community with support from NASA and industry in preparation for the 2020 Decadal Survey in Astrophysics. This is the final study report., Comment: 376 pages

We present Herschel Space Observatory Photodetector Array Camera and Spectrometer (PACS) and Spectral and Photometric Imaging Receiver Fourier Transform Spectrometer (SPIRE FTS) spectroscopy of a sample of twenty massive Young Stellar Objects (YSOs) in the Large and Small Magellanic Clouds (LMC and SMC). We analyse the brightest far infrared (far-IR) emission lines, that diagnose the conditions of the heated gas in the YSO envelope and pinpoint their physical origin.We compare the properties of massive Magellanic and Galactic YSOs.We find that [OI] and [CII] emission, that originates from the photodissociation region associated with the YSOs, is enhanced with respect to the dust continuum in the Magellanic sample. Furthermore the photoelectric heating efficiency is systematically higher for Magellanic YSOs, consistent with reduced grain charge in low metallicity environments. The observed CO emission is likely due to multiple shock components. The gas temperatures, derived from the analysis of CO rotational diagrams, are similar to Galactic estimates. This suggests a common origin to the observed CO excitation, from low-luminosity to massive YSOs, both in the Galaxy and the Magellanic Clouds. Bright far-IR line emission provides a mechanism to cool the YSO environment. We find that, even though [OI], CO and [CII] are the main line coolants, there is an indication that CO becomes less important at low metallicity, especially for the SMC sources. This is consistent with a reduction in CO abundance in environments where the dust is warmer due to reduced ultraviolet-shielding. Weak H$_2$O and OH emission is detected, consistent with a modest role in the energy balance of wider massive YSO environments., Comment: 38 pages; accepted for publication by MNRAS; full integrated version; journal version will include appendices on-line only

We present a comprehensive study of massive young stellar objects (YSOs) in the metal-poor galaxy NGC 6822 using IRAC and MIPS data obtained from the {\em Spitzer Space Telescope}. We find over 500 new YSO candidates in seven massive star-formation regions; these sources were selected using six colour-magnitude cuts. Via spectral energy distribution fitting to the data with YSO radiative transfer models we refine this list, identifying 105 high-confidence and 88 medium-confidence YSO candidates. For these sources we constrain their evolutionary state and estimate their physical properties. The majority of our YSO candidates are massive protostars with an accreting envelope in the initial stages of formation. We fit the mass distribution of the Stage I YSOs with a Kroupa initial mass function and determine a global star-formation rate of 0.039 $M_{\odot} yr^{-1}$. This is higher than star-formation rate estimates based on integrated UV fluxes. The new YSO candidates are preferentially located in clusters which correspond to seven active high-mass star-formation regions which are strongly correlated with the 8 and 24 $\mu$m emission from PAHs and warm dust. This analysis reveals an embedded high-mass star-formation region, Spitzer I, which hosts the highest number of massive YSO candidates in NGC 6822. The properties of Spitzer I suggest it is younger and more active than the other prominent H\,{\sc ii} and star-formation regions in the galaxy., Comment: 16 pages, 10 figures, 9 tables; accepted for publication in MNRAS

Models of the chemical evolution of the interstellar medium, galaxies, and the Universe rely on our understanding of the amounts and chemical composition of the material returned by stars and supernovae. Stellar yields are obtained from stellar-evolution models, which currently lack predictive prescriptions of stellar mass loss, although it significantly affects stellar lifetimes, nucleosynthesis, and chemical ejecta. Galaxy properties are derived from observations of the integrated light of bright member stars. Stars in the late stages of their evolution are among the infrared-brightest objects in galaxies. An unrealistic treatment of the mass-loss process introduces significant uncertainties in galaxy properties derived from their integrated light. We describe current efforts and future needs and opportunities to characterize AGB outflows: driving mechanisms, outflow rates, underlying fundamental physical and chemical processes such as dust grain formation, and dependency of these on metallicity., Comment: Submitted to Astro2020 decadal survey; 8 pages, 2 figures

We present the molecular cloud properties of N55 in the Large Magellanic Cloud using $^{12}$CO(1-0) and $^{13}$CO(1-0) observations obtained with Atacama Large Millimeter Array. We have done a detailed study of molecular gas properties, to understand how the cloud properties of N55 differ from Galactic clouds. Most CO emission appears clumpy in N55, and molecular cores that have YSOs show larger linewidths and masses. The massive clumps are associated with high and intermediate mass YSOs. The clump masses are determined by local thermodynamic equilibrium and virial analysis of the $^{12}$CO and $^{13}$CO emissions. These mass estimates lead to the conclusion that, (a) the clumps are in self-gravitational virial equilibrium, and (b) the $^{12}$CO(1-0)-to-H$_2$ conversion factor, X$_{\rm CO}$, is 6.5$\times$10$^{20}$cm$^{-2}$(K km s$^{-1}$)$^{-1}$. This CO-to-H$_2$ conversion factor for N55 clumps is measured at a spatial scale of $\sim$0.67 pc, which is about two times higher than the X$_{\rm CO}$ value of Orion cloud at a similar spatial scale. The core mass function of N55 clearly show a turnover below 200M$_{\odot}$, separating the low-mass end from the high-mass end. The low-mass end of the $^{12}$CO mass spectrum is fitted with a power law of index 0.5$\pm$0.1, while for $^{13}$CO it is fitted with a power law index 0.6$\pm$0.2. In the high-mass end, the core mass spectrum is fitted with a power index of 2.0$\pm$0.3 for $^{12}$CO, and with 2.5$\pm$0.4 for $^{13}$CO. This power-law behavior of the core mass function in N55 is consistent with many Galactic clouds., Comment: Accepted to ApJ

The survey for DUST in Nearby Galaxies with Spitzer (DUSTiNGS) identified several candidate Asymptotic Giant Branch (AGB) stars in nearby dwarf galaxies and showed that dust can form even in very metal-poor systems (Z ~ 0.008 $Z_\odot$). Here, we present a follow-up survey with WFC3/IR on the Hubble Space Telescope (HST), using filters that are capable of distinguishing carbon-rich (C-type) stars from oxygen-rich (M-type) stars: F127M, F139M, and F153M. We include six star-forming DUSTiNGS galaxies (NGC 147, IC 10, Pegasus dIrr, Sextans B, Sextans A, and Sag DIG), all more metal-poor than the Magellanic Clouds and spanning 1 dex in metallicity. We double the number of dusty AGB stars known in these galaxies and find that most are carbon rich. We also find 26 dusty M-type stars, mostly in IC 10. Given the large dust excess and tight spatial distribution of these M-type stars, they are most likely on the upper end of the AGB mass range (stars undergoing Hot Bottom Burning). Theoretical models do not predict significant dust production in metal-poor M-type stars, but we see evidence for dust excess around M-type stars even in the most metal-poor galaxies in our sample (12+log(O/H) = 7.26-7.50). The low metallicities and inferred high stellar masses (up to ~10 $M_\odot$) suggest that AGB stars can produce dust very early in the evolution of galaxies (~30 Myr after they form), and may contribute significantly to the dust reservoirs seen in high-redshift galaxies., Comment: 18 pages, 9 figures, accepted for publication in ApJ

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

We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta morphology is likely elongated. The difference of the CO and SiO line profiles is consistent with hydrodynamic simulations, which show that Rayleigh-Taylor instabilities cause mixing of gas, with heavier elements much more disturbed, making more elongated structure. We obtained isotopologue ratios of 28SiO/29SiO>13, 28SiO/30SiO>14, and 12CO/13CO>21, with the most likely limits of 28SiO/29SiO>128, 28SiO/30SiO>189. Low 29Si and 30Si abundances in SN 1987A are consistent with nucleosynthesis models that show inefficient formation of neutron-rich isotopes in a low metallicity environment, such as the Large Magellanic Cloud. The deduced large mass of HCO+ (~5x10^-6 Msun) and small SiS mass (<6x10^-5 Msun) might be explained by some mixing of elements immediately after the explosion. The mixing might have caused some hydrogen from the envelope to sink into carbon and oxygen-rich zones after the explosion, enabling the formation of a substantial mass of HCO+. Oxygen atoms may have penetrated into silicon and sulphur zones, suppressing formation of SiS. Our ALMA observations open up a new window to investigate chemistry, dynamics and explosive-nucleosynthesis in supernovae., Comment: Accepted by MNRAS, 16 pages

Context: Low- and intermediate-mass stars lose most of their stellar mass at the end of their lives on the asymptotic giant branch (AGB). Determining gas and dust mass-loss rates (MLRs) is important in quantifying the contribution of evolved stars to the enrichment of the interstellar medium. Aims: Attempt to, for the first time, spectrally resolve CO thermal line emission in a small sample of AGB stars in the Large Magellanic Cloud. Methods: ALMA was used to observe 2 OH/IR stars and 4 carbon stars in the LMC in the CO J= 2-1 line. Results: We present the first measurement of expansion velocities in extragalactic carbon stars. All four C-stars are detected and wind expansion velocities and stellar velocities are directly measured. Mass-loss rates are derived from modelling the spectral energy distribution and Spitzer/IRS spectrum with the DUSTY code. Gas-to-dust ratios are derived that make the predicted velocities agree with the observed ones. The expansion velocities and MLRs are compared to a Galactic sample of well-studied relatively low MLRs stars supplemented with "extreme" C-stars that have properties more similar to the LMC targets. Gas MLRs derived from a simple formula are significantly smaller than derived from the dust modelling, indicating an order of magnitude underestimate of the estimated CO abundance, time-variable mass loss, or that the CO intensities in LMC stars are lower than predicted by the formula derived for Galactic objects. This could be related to a stronger interstellar radiation field in the LMC. Conclusions: Although the LMC sample is small and the comparison to Galactic stars is non-trivial because of uncertainties in their distances it appears that for C stars the wind expansion velocities in the LMC are lower than in the solar neighbourhood, while the MLRs appear similar. This is in agreement with dynamical dust-driven wind models., Comment: Accepted for publication by A&A

Massive stars are key players in the evolution of galaxies, yet their formation pathway remains unclear. In this work, we use data from several galaxy-wide surveys to build an unbiased dataset of ~700 massive young stellar objects (MYSOs), ~200 giant molecular clouds (GMCs), and ~100 young (<10 Myr) optical stellar clusters (SCs) in the Large Magellanic Cloud. We employ this data to quantitatively study the location and clustering of massive star formation and its relation to the internal structure of GMCs. We reveal that massive stars do not typically form at the highest column densities nor centers of their parent GMCs at the ~6 pc resolution of our observations. Massive star formation clusters over multiple generations and on size scales much smaller than the size of the parent GMC. We find that massive star formation is significantly boosted in clouds near SCs. Yet, whether a cloud is associated with a SC does not depend on either the cloud's mass or global surface density. These results reveal a connection between different generations of massive stars on timescales up to 10 Myr. We compare our work with Galactic studies and discuss our findings in terms of GMC collapse, triggered star formation, and a potential dichotomy between low- and high-mass star formation., Comment: 13 pages, 7 figures, in press

This paper describes the beginning of the Far-Infrared Surveyor mission study for NASA's Astrophysics Decadal 2020. We describe the scope of the study, and the open process approach of the Science and Technology Definition Team. We are currently developing the science cases and provide some preliminary highlights here. We note key areas for technological innovation and improvements necessary to make a Far-Infrared Surveyor mission a reality., Comment: 8 pages, SPIE proceedings of the Space Telescopes and Instrumentation 2016: Optical, Infrared, and Millimeter Wave conference

We combine Spitzer and Herschel data of the star-forming region N11 in the Large Magellanic Cloud to produce detailed maps of the dust properties in the complex and study their variations with the ISM conditions. We also compare APEX/LABOCA 870um observations with our model predictions in order to decompose the 870um emission into dust and non-dust (free-free emission and CO(3-2) line) contributions. We find that in N11, the 870um can be fully accounted for by these 3 components. The dust surface density map of N11 is combined with HI and CO observations to study local variations in the gas-to-dust mass ratios. Our analysis leads to values lower than those expected from the LMC low-metallicity as well as to a decrease of the gas-to-dust mass ratio with the dust surface density. We explore potential hypotheses that could explain the low observed gas-to-dust mass ratios (variations in the XCO factor, presence of CO-dark gas or of optically thick HI or variations in the dust abundance in the dense regions). We finally decompose the local SEDs using a Principal Component Analysis (i.e. with no a priori assumption on the dust composition in the complex). Our results lead to a promising decomposition of the local SEDs in various dust components (hot, warm, cold) coherent with that expected for the region. Further analysis on a larger sample of galaxies will follow in order to understand how unique this decomposition is or how it evolves from one environment to another., Comment: 24 pages, 16 figures, accepted for publication in MNRAS

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.

The Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST) provides measurements over the wavelength range 5 to 28.5 microns. MIRI has, within a single 'package', four key scientific functions: photometric imaging, coronagraphy, single-source low-spectral resolving power (R ~ 100) spectroscopy, and medium-resolving power (R ~ 1500 to 3500) integral field spectroscopy. An associated cooler system maintains MIRI at its operating temperature of < 6.7 K. This paper describes the driving principles behind the design of MIRI, the primary design parameters, and their realization in terms of the 'as-built' instrument. It also describes the test program that led to delivery of the tested and calibrated Flight Model to NASA in 2012, and the confirmation after delivery of the key interface requirements.

MIRI (the Mid-Infrared Instrument for the James Webb Space Telescope (JWST)) operates from 5 to 28.5 microns and combines over this range: 1.) unprecedented sensitivity levels; 2.) sub-arcsec angular resolution; 3.) freedom from atmospheric interference; 4.) the inherent stability of observing in space; and 5.) a suite of versatile capabilities including imaging, low and medium resolution spectroscopy (with an integral field unit), and coronagraphy. We illustrate the potential uses of this unique combination of capabilities with various science examples: 1.) imaging exoplanets; 2.) transit and eclipse spectroscopy of exoplanets; 3.) probing the first stages of star and planet formation, including identifying bioactive molecules; 4.) determining star formation rates and mass growth as galaxies are assembled; and 5.) characterizing the youngest massive galaxies. This paper is the introduction to a series of ten covering all aspects of the instrument.

We present the first results from the SAGE-Var program, a follow on to the Spitzer legacy program Surveying the Agents of Galaxy Evolution (SAGE; Meixner, et al. 2006). We obtained 4 epochs of photometry at 3.6 & 4.5 microns covering the bar of the Large Magellanic Cloud (LMC) and the central region of the Small Magellanic Cloud (SMC) in order to probe the variability of extremely red sources missed by variability surveys conducted at shorter wavelengths, and to provide additional epochs of observation for known variables. Our 6 total epochs of observations allow us to probe infrared variability on 15 different timescales ranging from ~20 days to ~5 years. Out of a full catalog of 1,717,554 (LMC) and 457,760 (SMC) objects, we find 10 (LMC) and 6 (SMC) large amplitude AGB variables without optically measured variability owing to circumstellar dust obscuration. The catalog also contains multiple observations of known AGB variables, type I and II Cepheids, eclipsing variables, R CrB stars and young stellar objects which will be discussed in following papers. Here we present infrared Period-Luminosity (PL) relations for classical Cepheids in the Magellanic Clouds, as well as improved PL relationships for AGB stars pulsating in the fundamental mode using mean magnitudes constructed from 6 epochs of observations.

The rate at which interstellar gas is converted into stars, and its dependence on environment, is one of the pillars on which our understanding of the visible Universe is build. We present a comparison of the surface density of young stars (Sigma_*) and dust surface density (Sigma_d) across NGC346 (N66) in 115 independent pixels of 6x6 pc^2. We find a correlation between Sigma_* and Sigma_d with a considerable scatter. A power law fit to the data yields a steep relation with an exponent of 2.6+-0.2. We convert Sigma_d to gas surface density (Sigma_g) and Sigma_* to star formation rate (SFR) surface densities (Sigma_SFR), using simple assumptions for the gas-to-dust mass ratio and the duration of star formation. The derived total SFR (4+-1 10^-3 M_sun/yr) is consistent with SFR estimated from the Ha emission integrated over the Ha nebula. On small scales the Sigma_SFR derived using Ha systematically underestimates the count-based Sigma_SFR, by up to a factor of 10. This is due to ionizing photons escaping the area, where the stars are counted. We find that individual 36 pc^2 pixels fall systematically above integrated disc-galaxies in the Schmidt-Kennicutt diagram by on average a factor of ~7. The NGC346 average SFR over a larger area (90 pc radius) lies closer to the relation but remains high by a factor of ~3. The fraction of the total mass (gas plus young stars) locked in young stars is systematically high (~10 per cent) within the central 15 pc and systematically lower outside (2 per cent), which we interpret as variations in star formation efficiency. The inner 15 pc is dominated by young stars belonging to a centrally condensed cluster, while the outer parts are dominated by a dispersed population. Therefore, the observed trend could reflect a change of star formation efficiency between clustered and non-clustered star-formation., Comment: Accepted for publication in MNRAS

The DUSTiNGS survey (DUST in Nearby Galaxies with Spitzer) is a 3.6 and 4.5 micron imaging survey of 50 nearby dwarf galaxies designed to identify dust-producing Asymptotic Giant Branch (AGB) stars and massive stars. Using two epochs, spaced approximately six months apart, we identify a total of 526 dusty variable AGB stars (sometimes called "extreme" or x-AGB stars; [3.6]-[4.5]>0.1 mag). Of these, 111 are in galaxies with [Fe/H] < -1.5 and 12 are in galaxies with [Fe/H] < -2.0, making them the most metal-poor dust-producing AGB stars known. We compare these identifications to those in the literature and find that most are newly discovered large-amplitude variables, with the exception of approximately 30 stars in NGC 185 and NGC 147, one star in IC 1613, and one star in Phoenix. The chemical abundances of the x-AGB variables are unknown, but the low metallicities suggest that they are more likely to be carbon-rich than oxygen-rich and comparisons with existing optical and near-IR photometry confirms that 70 of the x-AGB variables are confirmed or likely carbon stars. We see an increase in the pulsation amplitude with increased dust production, supporting previous studies suggesting that dust production and pulsation are linked. We find no strong evidence linking dust production with metallicity, indicating that dust can form in very metal-poor environments., Comment: 18 pages, 16 figures, accepted for publication in ApJ

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

We present an analysis of mid-and far-infrared (IR) spectrum and spectral energy distribution (SED) of the LBV R71 in the LMC.This work aims to understand the overall contribution of high-mass LBVs to the total dust-mass budget of the interstellar medium (ISM) of the LMC and compare this with the contribution from low-mass asymptotic giant branch (AGB) stars. As a case study, we analyze the SED of R71. We compiled all the available photometric and spectroscopic observational fluxes from various telescopes for a wide wavelength range (0.36 -- 250\,$\mu$m). We determined the dust composition from the spectroscopic data, and derived the ejected dust mass, dust mass-loss rate, and other dust shell properties by modeling the SED of R71. We noted nine spectral features in the dust shell of R71 by analyzing Spitzer spectroscopic data. Among these, we identified three new crystalline silicate features. We computed our model spectrum by using 3D radiative transfer code MCMax. Our model calculation shows that dust is dominated by amorphous silicates, with some crystalline silicates, metallic iron, and a very tiny amount of polycyclic aromatic hydrocarbon (PAH) molecules. The presence of both silicates and PAHs indicates that the dust has a mixed chemistry. We derived a dust mass of 0.01 M$_\odot$, from which we arrive at a total ejected mass of $\approx$ 5 M$_\odot$. This implies a time-averaged dust mass-loss rate of 2.5$\times$10$^{-6}$ M$_\odot$\,yr$^{-1}$ with an explosion about 4000 years ago. We assume that the other five confirmed dusty LBVs in the LMC loose mass at a similar rate, and estimate the total contribution to the mass budget of the LMC to be $\approx$ 10$^{-5}$ M$_\odot$\,yr$^{-1}$, which is comparable to the contribution by all the AGB stars in the LMC. Based on our analysis on R71, we speculate that LBVs as a class may be an important dust source in the ISM of the LMC., Comment: 10 pages, 6 figures, 2 tables

To measure the mass loss from dusty oxygen-rich (O-rich) evolved stars in the Large Magellanic Cloud (LMC), we have constructed a grid of models of spherically-symmetric dust shells around stars with constant mass-loss rates using 2Dust. These models will constitute the O-rich model part of the "Grid of Red supergiant and Asymptotic giant branch star ModelS" (GRAMS). This model grid explores 4 parameters - stellar effective temperature from 2100 K - 4700 K; luminosity from 10^3-10^6 L_Sun; dust shell inner radii of 3, 7, 11, and 15 R_Star; and 10.0 micron optical depth from 10^-4 to 26. From an initial grid of ~1200 2Dust models, we create a larger grid of ~69,000 models by scaling to cover the luminosity range required by the data. These models are offered to the public on a website. The matching in color-magnitude diagrams and color-color diagrams to observed O-rich asymptotic giant branch (AGB) and red supergiant (RSG) candidate stars from the SAGE and SAGE-Spec LMC samples and a small sample of OH/IR stars is generally very good. The extreme AGB star candidates from SAGE are more consistent with carbon-rich (C-rich) than O-rich dust composition. Our model grid suggests lower limits to the mid-infrared colors of the dustiest AGB stars for which the chemistry could be O-rich. Finally, the fitting of GRAMS models to SEDs of sources fit by other studies provides additional verification of our grid and anticipates future, more expansive efforts., Comment: 49 pages, 15 figures, appeared in the 20 February 2011 issue of Astrophysical Journal

We observe a sample of 8 evolved stars in the Galactic Bulge in the CO J = 2 - 1 line using the Submillimeter Array (SMA) with angular resolution of 1 - 4 arcseconds. These stars have been detected previously at infrared wavelengths, and several of them have OH maser emission. We detect CO J = 2 - 1 emission from three of the sources in the sample: OH 359.943 +0.260, [SLO2003] A12, and [SLO2003] A51. We do not detect the remaining 5 stars in the sample because of heavy contamination from the galactic foreground CO emission. Combining CO data with observations at infrared wavelengths constraining dust mass loss from these stars, we determine the gas-to-dust ratios of the Galactic Bulge stars for which CO emission is detected. For OH 359.943 +0.260, we determine a gas mass-loss rate of 7.9 (+/- 2.2) x 10^-5 M_Sun/year and a gas-to-dust ratio of 310 (+/- 89). For [SLO2003] A12, we find a gas mass-loss rate of 5.4 (+/- 2.8) x 10^-5 M_Sun/year and a gas-to-dust ratio of 220 (+/- 110). For [SLO2003] A51, we find a gas mass-loss rate of 3.4 (+/- 3.0) x 10^-5 M_Sun/year and a gas-to-dust ratio of 160 (+/- 140), reflecting the low quality of our tentative detection of the CO J = 2 - 1 emission from A51. We find the CO J = 2 - 1 detections of OH/IR stars in the Galactic Bulge require lower average CO J = 2 - 1 backgrounds., Comment: 40 pages, 16 figures, appeared in the 1 March 2013 issue of the Astrophysical Journal

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., Comment: 37 pages, 3 figures, appeared in the 10 June 2010 issue of the Astrophysical Journal

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

We report spectroscopic and imaging observations of rotational transitions of cold CO and SiO in the ejecta of SN1987A, the first such emission detected in a supernova remnant. In addition to line luminosities for the CO J=1-0, 2-1, 6-5, and 7-6 transitions, we present upper limits for all other transitions up to J=13-12, collectively measured from the Atacama Large Millimeter Array (ALMA), the Atacama Pathfinder EXperiment (APEX), and the Herschel Spectral and Photometric Imaging REceiver (SPIRE). Simple models show the lines are emitted from at least 0.01 solar masses of CO at a temperature > 14 K, confined within at most 35% of a spherical volume expanding at ~ 2000 km/s. Moreover, we locate the emission within 1'' of the central debris. These observations, along with a partial observation of SiO, confirm the presence of cold molecular gas within supernova remnants and provide insight into the physical conditions and chemical processes in the ejecta. Furthermore, we demonstrate the powerful new window into supernova ejecta offered by submillimeter observations., Comment: Accepted to the Astrophysical Journal Letters, 6 pages, 3 figures

The Magellanic Quasars Survey (MQS) has now increased the number of quasars known behind the Magellanic Clouds by almost an order of magnitude. All survey fields in the Large Magellanic Cloud (LMC) and 70% of those in the Small Magellanic Cloud (SMC) have been observed. The targets were selected from the third phase of the Optical Gravitational Lensing Experiment (OGLE-III) based on their optical variability, mid-IR and/or X-ray properties. We spectroscopically confirmed 758 (565 LMC and 193 SMC) quasars behind the Clouds, of which 94% (527 LMC and 186 SMC) are newly identified. The MQS quasars have long-term (12 years and growing for OGLE), high-cadence light curves, enabling unprecedented variability studies of quasars. The MQS quasars also provide a dense reference grid for measuring both the internal and bulk proper motions of the Clouds, and 50 quasars are bright enough (I<18 mag) for absorption studies of the interstellar/galactic (ISM/IGM) medium of the Clouds., Comment: 32 pages that include 10 figures and 4 tables (full Tables 2 and 3 included in the source file); submitted to ApJ

We summarize several new developments in the nuclear equation of state for supernova simulations and neutron stars. We discuss an updated and improved Notre-Dame-Livermore Equation of State (NDL EoS) for use in supernovae simulations. This Eos contains many updates. Among them are the effects of 3- body nuclear forces at high densities and the possible transition to a QCD chiral and/or super-conducting color phase at densities. We also consider the neutron star equation of state and neutrino transport in the presence of strong magnetic fields. We study a new quantum hadrodynamic (QHD) equation of state for neutron stars (with and without hyperons) in the presence of strong magnetic fields. The parameters are constrained by deduced masses and radii. The calculated adiabatic index for these magnetized neutron stars exhibit rapid changes with density. This may provide a mechanism for star-quakes and flares in magnetars. We also investigate the strong magnetic field effects on the moments of inertia and spin down of neutron stars. The change of the moment of inertia associated with emitted magnetic flares is shown to match well with observed glitches in some magnetars. We also discuss a perturbative calculation of neutrino scattering and absorption in hot and dense hyperonic neutron-star matter in the presence of a strong magnetic field. The absorption cross-sections show a remarkable angular dependence in that the neutrino absorption strength is reduced in a direction parallel to the magnetic field and enhanced in the opposite direction. The pulsar kick velocities associated with this asymmetry comparable to observed pulsar velocities and may affect the early spin down rate of proto-neutron star magnetars with a toroidal field configuration., Comment: 6 pages, revised 1 figure, International Symposium "Exotic Nuclear Structure From Nucleons" (ENSFN 2012)

We present a study of the infrared/submm emission of the LMC star forming complex N158-N159-N160. Combining observations from the Spitzer Space Telescope (3.6-70um), the Herschel Space Observatory (100-500um) and LABOCA (870um) allows us to work at the best angular resolution available now for an extragalactic source. We observe a remarkably good correlation between SPIRE and LABOCA emission and resolve the low surface brightnesses emission. We use the Spitzer and Herschel data to perform a resolved Spectral Energy Distribution (SED) modelling of the complex. Using MBB, we derive a global emissivity index beta_c of 1.47. If beta cold is fixed to 1.5, we find an average temperature of 27K. We also apply the Galliano et al. (2011) modelling technique (and amorphous carbon to model carbon dust) to derive maps of the star formation rate, the mean starlight intensity, the fraction of PAHs or the dust mass surface density of the region. We observe that the PAH fraction strongly decreases in the HII regions. This decrease coincides with peaks in the mean radiation field intensity map. The dust surface densities follow the FIR distribution, with a total dust mass of 2.1x10^4 Msolar (2.8 times less than when using graphite grains) in the resolved elements we model. We find a non-negligible amount of dust in the molecular cloud N159 South (showing no massive SF). We also investigate the drivers of the Herschel/PACS and SPIRE submm colours as well as the variations in the gas-to-dust mass ratio (G/D) and the XCO conversion factor in the region N159. We finally model individual regions to analyse variations in the SED shape across the complex and the 870um emission in more details. No measurable submm excess emission at 870um seems to be detected in these regions., Comment: Accepted for publication in MNRAS; 23 pages, 11 figures, 2 tables

We performed multiwavelength observations of the young planetary nebula (PN) M1-11 and obtained its elemental abundances, dust mass, and the evolutionary status of the central star. The AKARI/IRC, VLT/VISIR, and Spitzer/IRS spectra show features due to carbon-rich dust, such as the 3.3, 8.6, and 11.3 um features due to polycyclic aromatic hydrocarbons (PAHs), a smooth continuum attributable to amorphous carbon, and the broad 11.5 and 30 um features often ascribed to SiC and MgS, respectively. We also report the presence of an unidentified broad feature at 16-22 um, similar to the feature found in Magellanic Cloud PNe with either C-rich or O-rich gas-phase compositions. We identify for the first time in M1-11 spectral lines at 8.5 (blended with PAH), 17.3, and 18.9 um that we attribute to the C60 fullerene. This identification is strengthened by the fact that other Galactic PNe in which fullerenes are detected, have similar central stars, similar gas-phase abundances, and a similar dust composition to M1-11. The weak radiation field due to the relatively cool central stars in these PNe may provide favorable conditions for fullerenes to survive in the circumstellar medium. Using the photo-ionization code CLOUDY, combined with a modified blackbody, we have fitted the ~0.1-90 um spectral energy distribution and determined the dust mass in the nebula to be ~3.5x10^{-4} Msun$. Our chemical abundance analysis and SED model suggest that M1-11 is perhaps a C-rich PN with C/O ratio in the gas-phase of +0.19 dex, and that it evolved from a 1-1.5 Msun star., Comment: Accepted for publication by the Astrophysical Journal

We present new observations of 34 YSO candidates in the SMC. The anchor of the analysis is a set of Spitzer-IRS spectra, supplemented by groundbased 3-5 micron spectra, Spitzer and NIR photometry, optical spectroscopy and radio data. The sources' SEDs and spectral indices are consistent with embedded YSOs; prominent silicate absorption is observed in the spectra of at least ten sources, silicate emission is observed towards four sources. PAH emission is detected towards all but two sources. Based on band ratios (in particular the strength of the 11.3 micron and the weakness of the 8.6 micron bands) PAH emission towards SMC YSOs is dominated by predominantly small neutral grains. Ice absorption is observed towards fourteen sources in the SMC. The comparison of H2O and CO2 ice column densities for SMC, LMC and Galactic samples suggests that there is a significant H2O column density threshold for the detection of CO2 ice. This supports the scenario proposed by Oliveira et al. (2011), where the reduced shielding in metal-poor environments depletes the H2O column density in the outer regions of the YSO envelopes. No CO ice is detected towards the SMC sources. Emission due to pure-rotational 0-0 transitions of H2 is detected towards the majority of SMC sources, allowing us to estimate rotational temperatures and column densities. All but one source are spectroscopically confirmed as SMC YSOs. Of the 33 YSOs identified in the SMC, 30 sources populate different stages of massive stellar evolution. The remaining three sources are classified as intermediate-mass YSOs with a thick dusty disc and a tenuous envelope still present. We propose one of the sources is a D-type symbiotic system, based on the presence of Raman, H and He emission lines in the optical spectrum, and silicate emission in the IRS-spectrum. This would be the first dust-rich symbiotic system identified in the SMC. (abridged), Comment: 34 pages, 13 figures, 4 tables, accepted for publication in MNRAS main journal

We present nebular phase optical imaging and spectroscopy and near/mid-IR imaging of the Type II SN 2006bc. Observations reveal the central wavelength of the symmetric H$\alpha$ line profile to be red-shifted with respect to the host galaxy H$\alpha$ emission by day 325. Such an phenomenon has been argued to result from an asymmetric explosion in the iron-peak elements resulting in a larger mass of $^{56}$Ni and higher excitation of hydrogen on the far side of the SN explosion. We also observe a gradual blue-shifting of this H$\alpha$ peak which is indicative of dust formation in the ejecta. Although showing a normal peak brightness, V $\sim$ -17.2, for a core-collapse SN, 2006bc fades by $\sim$6 mag during the first 400 days suggesting either a relatively low $^{56}$Ni yield, an increase in extinction due to new dust, or both. A short duration flattening of the light curve is observed from day 416 to day 541 suggesting an optical light echo. Based on the narrow time window of this echo, we discuss implications on the location and geometry of the reflecting ISM. With our radiative transfer models, we find an upper limit of 2 x 10$^{-3}$ M$_{\odot}$ of dust around SN 2006bc. In the event that all of this dust were formed during the SN explosion, this quantity of dust is still several orders of magnitude lower than that needed to explain the large quantities of dust observed in the early universe., Comment: 6 pages, 10 figures, accepted for publication in ApJ

We present near infra-red light curves of supernova (SN) 2011fe in M101, including 34 epochs in H band starting fourteen days before maximum brightness in the B-band. The light curve data were obtained with the WIYN High-Resolution Infrared Camera (WHIRC). When the data are calibrated using templates of other Type Ia SNe, we derive an apparent H-band magnitude at the epoch of B-band maximum of 10.85 \pm 0.04. This implies a distance modulus for M101 that ranges from 28.86 to 29.17 mag, depending on which absolute calibration for Type Ia SNe is used., Comment: 9 pages, 3 figures, emulateapj style, accepted for publication in The Astrophysical Journal

We present results from the first application of the Grid of Red Supergiant and Asymptotic Giant Branch ModelS (GRAMS) model grid to the entire evolved stellar population of the Large Magellanic Cloud (LMC). GRAMS is a pre-computed grid of 80,843 radiative transfer (RT) models of evolved stars and circumstellar dust shells composed of either silicate or carbonaceous dust. We fit GRAMS models to ~30,000 Asymptotic Giant Branch (AGB) and Red Supergiant (RSG) stars in the LMC, using 12 bands of photometry from the optical to the mid-infrared. Our published dataset consists of thousands of evolved stars with individually determined evolutionary parameters such as luminosity and mass-loss rate. The GRAMS grid has a greater than 80% accuracy rate discriminating between Oxygen- and Carbon-rich chemistry. The global dust injection rate to the interstellar medium (ISM) of the LMC from RSGs and AGB stars is on the order of 1.5x10^(-5) solar masses/yr, equivalent to a total mass injection rate (including the gas) into the ISM of ~5x10^(-3) solar masses/yr. Carbon stars inject two and a half times as much dust into the ISM as do O-rich AGB stars, but the same amount of mass. We determine a bolometric correction factor for C-rich AGB stars in the K band as a function of J - K color, BC(K) = -0.40(J-K)^2 + 1.83(J-K) + 1.29. We determine several IR color proxies for the dust mass-loss rate (MLR) from C-rich AGB stars, such as log (MLR) = (-18.90)/((K-[8.0])+3.37)-5.93. We find that a larger fraction of AGB stars exhibiting the `long-secondary period' phenomenon are O-rich than stars dominated by radial pulsations, and AGB stars without detectable mass-loss do not appear on either the first-overtone or fundamental-mode pulsation sequences., Comment: 19 pages, 19 figures

We estimate the total dust input from the cool evolved stars in the Small Magellanic Cloud (SMC), using the 8 micron excess emission as a proxy for the dust-production rate. We find that Asymptotic Giant Branch (AGB) and red supergiant (RSG) stars produce (8.6-9.5) x 10^7 solar masses per year of dust, depending on the fraction of far-infrared sources that belong to the evolved star population (with 10%-50% uncertainty in individual dust-production rates). RSGs contribute the least (<4%), while carbon-rich AGB stars (especially the so-called "extreme" AGB stars) account for 87%-89% of the total dust input from cool evolved stars. We also estimate the dust input from hot stars and supernovae (SNe), and find that if SNe produce 10^-3 solar masses of dust each, then the total SN dust input and AGB input are roughly equivalent. We consider several scenarios of SNe dust production and destruction and find that the interstellar medium (ISM) dust can be accounted for solely by stellar sources if all SNe produce dust in the quantities seen around the dustiest examples and if most SNe explode in dense regions where much of the ISM dust is shielded from the shocks. We find that AGB stars contribute only 2.1% of the ISM dust. Without a net positive contribution from SNe to the dust budget, this suggests that dust must grow in the ISM or be formed by another unknown mechanism., Comment: 10 pages, 5 figures. Accepted for publication in ApJ

We present an analysis of multi-epoch observations of the Type II-P supernova SN 2004et. New and archival optical spectra of SN 2004et are used to study the evolution of the Halpha and [O I] 6300A line profiles between days 259 and 646. Mid-infrared imaging was carried out between 2004 to 2010. We include Spitzer `warm' mission photometry at 3.6 and 4.5um obtained on days 1779, 1931 and 2151, along with ground-based and HST optical and near-infrared observations obtained between days 79 and 1803. Multi-wavelength light curves are presented, as well as optical-infrared spectral energy distributions (SEDs) for multiple epochs. Starting from about day 300, the optical light curves provide evidence for an increasing amount of circumstellar extinction attributable to newly formed dust, with the additional extinction reaching 0.8-1.5 magnitudes in the V-band by day 690. The overall SEDs were fitted with multiple blackbody components, in order to investigate the luminosity evolution of the supernova, and then with Monte Carlo radiative transfer models using smooth or clumpy dust distributions, in order to estimate how much new dust condensed in the ejecta. The luminosity evolution was consistent with the decay of 56Co in the ejecta up until about day 690, after which an additional emission source is required, in agreement with the findings of Kotak et al. (2009). Clumped dust density distributions consisting of 20% amorphous carbons and 80% silicates by mass were able to match the observed optical and infrared SEDs, with dust masses that increased from 8x10^{-5} Msun on day 300 to 1.5x10^{-3} Msun on day 690, still significantly lower than the values needed for core collapse supernovae to make a significant contribution to the dust enrichment of galaxies., Comment: 24 pages, 12 figures, 9 tables, published in MNRAS

Aims: In this paper, we perform detailed modelling of the Spitzer and Herschel observations of the LMC, in order to: (i) systematically study the uncertainties and biases affecting dust mass estimates; and to (ii) explore the peculiar ISM properties of the LMC. Methods: To achieve these goals, we have modelled the spatially resolved SEDs with two alternate grain compositions, to study the impact of different submillimetre opacities on the dust mass. We have rigorously propagated the observational errors (noise and calibration) through the entire fitting process, in order to derive consistent parameter uncertainties. Results: First, we show that using the integrated SED leads to underestimating the dust mass by ~50 % compared to the value obtained with sufficient spatial resolution, for the region we studied. This might be the case, in general, for unresolved galaxies. Second, we show that Milky Way type grains produce higher gas-to-dust mass ratios than what seems possible according to the element abundances in the LMC. A spatial analysis shows that this dilemma is the result of an exceptional property: the grains of the LMC have on average a larger intrinsic submm opacity (emissivity index beta~1.7 and opacity kappa_abs(160 microns)=1.6 m2/kg) than those of the Galaxy. By studying the spatial distribution of the gas-to-dust mass ratio, we are able to constrain the fraction of unseen gas mass between ~10, and ~100 % and show that it is not sufficient to explain the gas-to-dust mass ratio obtained with Milky Way type grains. Finally, we confirm the detection of a 500 microns extended emission excess with an average relative amplitude of ~15 %, varying up to 40 %. This excess anticorrelates well with the dust mass surface density. Although we do not know the origin of this excess, we show that it is unlikely the result of very cold dust, or CMB fluctuations., Comment: 33 pages, 36 figures, accepted for publication in Astronomy and Astrophysics

Aims. In this paper we study the Spitzer and TIMMI2 infrared spectra of post-AGB disc sources, both in the Galaxy and the LMC. Using the observed infrared spectra we determine the mineralogy and dust parameters of the discs, and look for possible differences between the Galactic and extragalactic sources. Methods. Modelling the full spectral range observed allows us to determine the dust species present in the disc and different physical parameters such as grain sizes, dust abundance ratios, and the dust and continuum temperatures. Results. We find that all the discs are dominated by emission features of crystalline and amorphous silicate dust. Only a few sample sources show features due to CO2 gas or carbonaceous molecules such as PAHs and C60 fullerenes. Our analysis shows that dust grain processing in these discs is strong, resulting in large average grain sizes and a very high crystallinity fraction. However, we do not find any correlations between the derived dust parameters and properties of the central source. There also does not seem to be a noticeable difference between the mineralogy of the Galactic and LMC sources. Even though the observed spectra are very similar to those of protoplanetary discs around young stars, showing similar mineralogy and strong grain processing, we do find evidence for differences in the physical and chemical processes of the dust processing., Comment: Accepted for publication in Astronomy and Astrophysics, 16 pages, 12 pages online material, 23 figures

We report far-infrared and submillimeter observations of Supernova 1987A, the star that exploded on February 23, 1987 in the Large Magellanic Cloud, a galaxy located 160,000 lightyears away. The observations reveal the presence of a population of cold dust grains radiating with a temperature of ~17-23 K at a rate of about 220 solar luminosity. The intensity and spectral energy distribution of the emission suggests a dust mass of ~0.4-0.7 solar mass. The radiation must originate from the SN ejecta and requires the efficient precipitation of all refractory material into dust. Our observations imply that supernovae can produce the large dust masses detected in young galaxies at very high redshifts., Comment: This is the author's version of the work. It is posted here by permission of the AAAS for personal use, not for redistribution. To appear in Science; available in Science Express on 7th July, DOI:10.1126/science.1205983 at http://www.sciencemag.org/lookup/doi/10.1126/science.1205983 ; 14 pages, 2 figures, 4 tables. High resolution images are available from the authors