Determining the efficiency with which gas is converted into stars in galaxies requires an accurate determination of the total reservoir of molecular gas mass. However, despite being the most abundant molecule in the Universe, H$_2$ is challenging to detect through direct observations and indirect methods have to be used to estimate the total molecular gas reservoir. These are often based on scaling relations from tracers such as CO or dust, and are generally calibrated in the Milky Way. Yet, evidence that these scaling relations are environmentally dependent is growing. In particular, the commonly used CO-to-H$_2$ conversion factor (X$_{\rm CO}$) is expected to be higher in metal-poor and/or strongly UV-irradiated environments. We use new SOFIA/FIFI-LS observations of far-infrared fine structure lines from the ionised and neutral gas and the Meudon photodissociation region model to constrain the physical properties and the structure of the gas in the massive star-forming region of 30 Doradus in the Large Magellanic Cloud, and determine the spatially resolved distribution of the total reservoir of molecular gas in the proximity of the young massive cluster R136. We compare this value with the molecular gas mass inferred from ground-based CO observations and dust-based estimates to quantify the impact of this extreme environment on commonly used tracers of the molecular gas. We find that the strong radiation field combined with the half-solar metallicity of the surrounding gas are responsible for a large reservoir of "CO-dark" molecular gas, leaving a large fraction of the total H$_2$ gas (> 75%) undetected when adopting a standard X$_{\rm CO}$ factor in this massive star-forming region.

The nearby ($\sim$500 kpc) metal-poor ([Fe/H] $\approx$ -1.2; $Z$ $\approx$ 30% $Z_{\odot}$) star-forming galaxy NGC 6822 has a metallicity similar to systems at the epoch of peak star formation. Through identification and study of dusty and dust-producing stars, it is therefore a useful laboratory to shed light on the dust life cycle in the early Universe. We present a catalog of sources combining near- and mid-IR photometry from the United Kingdom Infrared Telescope (UKIRT; $J$, $H$, and $K$) and the $Spitzer$ $Space$ $Telescope$ (IRAC 3.6, 4.5, 5.8, and 8.0 $\mu$m and MIPS 24 $\mu$m). This catalog is employed to identify dusty and evolved stars in NGC 6822 utilizing three color-magnitude diagrams (CMDs). With diagnostic CMDs covering a wavelength range spanning the near- and mid-IR, we develop color cuts using kernel density estimate (KDE) techniques to identify dust-producing evolved stars, including red supergiant (RSG) and thermally-pulsing asymptotic giant branch (TP-AGB) star candidates. In total, we report 1,292 RSG candidates, 1,050 oxygen-rich AGB star candidates, and 560 carbon-rich AGB star candidates with high confidence in NGC 6822. Our analysis of the AGB stars suggests a robust population inhabiting the central stellar bar of the galaxy, with a measured global stellar metallicity of [Fe/H] = -1.286 $\pm$ 0.095, consistent with previous studies. In addition, we identify 277 young stellar object (YSO) candidates. The detection of a large number of YSO candidates within a centrally-located, compact cluster reveals the existence of an embedded, high-mass star-formation region that has eluded previous detailed study. Spitzer I appears to be younger and more active than the other prominent star-forming regions in the galaxy.

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.

We have conducted ALMA CO isotopes and 1.3 mm continuum observations toward filamentary molecular clouds of the N159W-South region in the Large Magellanic Cloud with an angular resolution of $\sim$0"25 ($\sim$0.07 pc). Although the previous lower-resolution ($\sim$1") ALMA observations revealed that there is a high-mass protostellar object at an intersection of two line-shaped filaments in $^{13}$CO with the length scale of $\sim$10 pc, the spatially resolved observations, in particular, toward the highest column density part traced by the 1.3 mm continuum emission, the N159W-South clump, show complicated hub-filamentary structures. We also discovered that there are multiple protostellar sources with bipolar outflows along the massive filament. The redshifted/blueshifted components of the $^{13}$CO emission around the massive filaments/protostars have complementary distributions, which is considered to be a possible piece of evidence for a cloud-cloud collision. We propose a new scenario in which the supersonically colliding gas flow triggers the formation of both the massive filament and protostars. This is a modification of the earlier scenario of cloud-cloud collision, by Fukui et al., that postulated the two filamentary clouds occur prior to the high-mass star formation. A recent theoretical study of the shock compression in colliding molecular flows by Inoue et al. demonstrates that the formation of filaments with hub structure is a usual outcome of the collision, lending support for the present scenario. The theory argues that the filaments are formed as dense parts in a shock compressed sheet-like layer, which resembles $"$an umbrella with pokes.$"$

We present the ALMA observations of CO isotopes and 1.3 mm continuum emission toward the N159E-Papillon Nebula in the Large Magellanic Cloud (LMC). The spatial resolution is 0"25-0"28 (0.06-0.07 pc), which is a factor of 3 higher than the previous ALMA observations in this region. The high resolution allowed us to resolve highly filamentary CO distributions with typical widths of $\sim$0.1 pc (full width half maximum) and line masses of a few 100 $M_{\odot}$ pc$^{-1}$. The filaments (more than ten in number) show an outstanding hub-filament structure emanating from the nebular center toward the north. We identified for the first time two massive protostellar outflows of $\sim$10$^4$ yr dynamical age along one of the most massive filaments. The observations also revealed several pillar-like CO features around the Nebula. The H II region and the pillars have a complementary spatial distribution and the column density of the pillars is an order of magnitude higher than that of the pillars in the Eagle nebula (M16) in the Galaxy, suggesting an early stage of pillar formation with an age younger than $\sim$10$^5$ yr. We suggest that a cloud-cloud collision triggered the formation of the filaments and protostar within the last $\sim$2 Myr. It is possible that the collision is more recent, as part of the kpc-scale H I flows come from the tidal interaction resulting from the close encounter between the LMC and SMC $\sim$200 Myr ago as suggested for R136 by Fukui et al.

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.

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.

We present a new initiative by the Women in Astronomy Forum at Space Telescope Science Institute (STScI) to increase gender diversity and inclusion in STScI's scientific committees and the activities they generate. This initiative offers new and uniform guidelines on binary gender representation goals for each committee and recommendations on how to achieve them in a homogeneous way, as well as metrics and tools to track progress towards defined goals. While the new guidelines presented in the paper focus on binary gender representation, they can be adapted and implemented to support all minority groups. By creating diverse committees and making them aware of, and trained on implicit bias, we expect to create a diverse outcome in the activities they generate, which, in turn, will advance science further and faster.

Dust offers a unique probe of the interstellar medium (ISM) across multiple size, density, and temperature scales. Dust is detected in outflows of evolved stars, star-forming molecular clouds, planet-forming disks, and even in galaxies at the dawn of the Universe. These grains also have a profound effect on various astrophysical phenomena from thermal balance and extinction in galaxies to the building blocks for planets, and changes in dust grain properties will affect all of these phenomena. A full understanding of dust in all of its forms and stages requires a multi-disciplinary investigation of the dust life cycle. Such an investigation can be achieved with a statistical study of dust properties across stellar evolution, star and planet formation, and redshift. Current and future instrumentation will enable this investigation through fast and sensitive observations in dust continuum, polarization, and spectroscopy from near-infrared to millimeter wavelengths.

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.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

The aim of this study is to characterize the distribution and basic properties of the natal gas associated with high-mass young stellar objects (YSOs) in isolated environments in the Large Magellanic Cloud (LMC). High-mass stars usually form in Giant Molecular Clouds (GMCs) as part of a young stellar cluster, but some OB stars are observed far from GMCs. By examining the spatial coincidence between the high-mass YSOs and 12CO (J = 1-0) emission detected by NANTEN and Mopra observations, we selected ten high-mass YSOs that are located away from any of the NANTEN clouds but are detected by the Mopra pointed observations. The ALMA observations revealed that a compact molecular cloud whose mass is a few thousand solar masses or smaller is associated with the high-mass YSOs, which indicates that these compact clouds are the sites of high-mass star formation. The high-density and high-temperature throughout the clouds are explained by the severe photodissociation of CO due to the lower metallicity than in the Galaxy. The star formation efficiency ranges from several to as high as ~ 40%, indicating efficient star formation in these environments. The enhanced turbulence may be a cause of the efficient star formation therein, as judged from the gas velocity information and the association with the lower density gas.

Metal Evolution, Transport, and Abundance in the LMC (METAL) is a large Cycle 24 program on the Hubble Space Telescope aimed at measuring dust extinction properties and interstellar depletions in the Large Magellanic Cloud (LMC) at half-solar metallicity. The 101-orbit program is comprised of COS and STIS spectroscopy toward 33 LMC massive stars between 1150 A and 3180 A, and parallel WFC3 imaging in 7 NUV-NIR filters. The fraction of silicon in the gas-phase (depletion) obtained from the spectroscopy decreases with increasing hydrogen column density. Depletion patterns for silicon differ between the Milky Way, LMC, and SMC, with the silicon depletion level offsetting almost exactly the metallicity differences, leading to constant gas-phase abundances in those galaxies for a given hydrogen column density. The silicon depletion correlates linearly with the absolute-to-selective extinction, R$_V$, indicating a link between gas depletion and dust grain size. Extinction maps are derived from the resolved stellar photometry in the parallel imaging, which can be compared to FIR images from Herschel and Spitzer to estimate the emissivity of dust at LMC metallicity. The full METAL sample of depletions, UV extinction curves, and extinction maps will inform the abundance, size, composition, and optical properties of dust grains in the LMC, comprehensively improve our understanding of dust properties, and the accuracy with which dust-based gas masses, star formation rates and histories in nearby and high-redshift galaxies are estimated. This overview paper describes the goals, design, data reduction, and initial results of the METAL survey.

As nearby neighbors to the Milky Way, the Large and Small Magellanic Clouds (LMC and SMC) provide a unique opportunity to study star formation in the context of their galactic ecosystems. Thousands of young stellar objects (YSOs) have been characterized with large-scale Spitzer and Herschel surveys. In this paper, we present new near-IR spectroscopy of five high-mass YSOs in the LMC and one in the SMC. We detect multiple hydrogen recombination lines, as well as He I 2.058 $\mu$m, H$_2$, [Fe II], and [S III] in these highly excited sources. We estimate the internal extinction of each source and find that it is highest for sources with the youngest evolutionary classifications. Using line ratios, we assess the dominant excitation mechanism in the three sources where we detect both H$_2$ 2.12 $\mu$m and [Fe II] 1.64 $\mu$m. In each case, photoexcitation dominates over shock excitation. Finally, we detect CO bandhead absorption in one of our LMC sources. While this feature is often associated with evolved stars, this object is likely young with strong PAH and fine-structure emission lines tracing an H II region detected at longer wavelengths. Compared to high-mass YSOs in the Galaxy, our sources have higher bolometric and line luminosities, consistent with their selection as some of the brightest sources in the LMC and SMC.

The equatorial ring of Supernova (SN) 1987A has been exposed to forward shocks from the SN blast wave, and it has been suggested that these forward shocks have been causing on-going destruction of dust in the ring. We obtained SOFIA FORCAST 11.1, 19.7 and 31.5 micron photometry of SN\,1987A in 2016. Compared with Spitzer measurements 10 years earlier, the 31.5 micron flux has significantly increased. The excess at 31.5 micron appears to be related to the Herschel 70 micron excess, which was detected 5 years earlier. The dust mass needed to account for the the 31.5--70 micron excess is 3--7x10^-4 Msun, more than ten times larger than the ring dust mass (1x10^-5 Msun) estimate from the data 10-years earlier. We argue that dust grains are re-formed or grown in the post-shock regions in the ring after forward shocks have destroyed pre-existing dust grains in the ring and released refractory elements into gas. In the post-shock region, atoms can stick to surviving dust grains, and the dust mass may have increased (grain growth), or dust grains might have condensed directly from the gas. An alternative possibility is that the outer part of the expanding ejecta dust might have been heated by X-ray emission from the circumstellar ring. The future development of this excess could reveal whether grains are reformed in the post-shocked region of the ring or eject dust is heated by X-ray.

The Origins Space Telescope, one of four large Mission Concept studies sponsored by NASA for review in the 2020 US Astrophysics Decadal Survey, will open unprecedented discovery space in the infrared, unveiling our cosmic origins. We briefly describe in this article the key science themes and architecture for OST. With a sensitivity gain of up to a factor of 1,000 over any previous or planned mission, OST will open unprecedented discovery space, allow us to peer through an infrared window teeming with possibility. OST will fundamentally change our understanding of our cosmic origins - from the growth of galaxies and black holes, to uncovering the trail of water, to life signs in nearby Earth-size planets, and discoveries never imagined. Built to be highly adaptable, while addressing key science across many areas of astrophysics, OST will usher in a new era of infrared astronomy.

At a time when ALMA produces spectacular high resolution images of gas and dust in circumstellar disks, the next observational frontier in our understanding of planet formation and the chemistry of planet-forming material may be found in the mid- to far-infrared wavelength range. A large, actively cooled far-infrared telescope in space will offer enormous spectroscopic sensitivity improvements of 3-4 orders of magnitude, making it possible to uniquely survey certain fundamental properties of planet formation. Specifically, the Origins Space Telescope (OST), a NASA flagship concept to be submitted to the 2020 decadal survey, will provide a platform that allows complete surveys of warm and cold water around young stars of all masses and across all evolutionary stages, and to measure their total planet-forming gas mass using the ground-state line of HD. While this white paper is formulated in the context of the NASA Origins Space Telescope concept, it can be applied in general to inform any future space-based, cold far-infrared observatory.

The Origins Space Telescope (OST) is one of four mission concepts currently being studied by NASA in preparation for the Astrophysics 2020 Decadal Survey. With active cooling (~4 K), OST will be sensitive in mid- to far-IR wavelengths, using imaging and spectroscopy to probe the furthest reaches of our galaxies, trace the path of water through star and planet formation, and place thermochemical constraints on the atmospheres of exoplanets ranging in size from Jupiter to Earth. This contribution to the Exoplanet Science Strategy committee discusses the significant advancements that the OST Mid-Infrared Imager, Spectrometer, and Coronagraph (MISC) instrument can make in studying cool planetary atmospheres. We particularly focus on the atmospheres of transiting rocky planets in the habitable zones of mid-to-late M stars. We discuss how OST thermal infrared observations can significantly enhance our understanding of the temperature structure and molecular abundances of biologically interesting gases on these worlds, including O3, CH4, H2O, and CO2.

We report the first extragalactic detection of the complex organic molecules (COMs) dimethyl ether (CH$_3$OCH$_3$) and methyl formate (CH$_3$OCHO) with the Atacama Large Millimeter/submillimeter Array (ALMA). These COMs together with their parent species methanol (CH$_3$OH), were detected toward two 1.3 mm continuum sources in the N 113 star-forming region in the low-metallicity Large Magellanic Cloud (LMC). Rotational temperatures ($T_{\rm rot}\sim130$ K) and total column densities ($N_{\rm rot}\sim10^{16}$ cm$^{-2}$) have been calculated for each source based on multiple transitions of CH$_3$OH. We present the ALMA molecular emission maps for COMs and measured abundances for all detected species. The physical and chemical properties of two sources with COMs detection, and the association with H$_2$O and OH maser emission indicate that they are hot cores. The fractional abundances of COMs scaled by a factor of 2.5 to account for the lower metallicity in the LMC are comparable to those found at the lower end of the range in Galactic hot cores. Our results have important implications for studies of organic chemistry at higher redshift.

The Small Magellanic Cloud (SMC) provides the only laboratory to study the structure of molecular gas at high resolution and low metallicity. We present results from the Herschel Spectroscopic Survey of the SMC (HS$^{3}$), which mapped the key far-IR cooling lines [CII], [OI], [NII], and [OIII] in five star-forming regions, and new ALMA 7m-array maps of $^{12}$CO and $^{13}$CO $(2-1)$ with coverage overlapping four of the five HS$^{3}$ regions. We detect [CII] and [OI] throughout all of the regions mapped. The data allow us to compare the structure of the molecular clouds and surrounding photodissociation regions using $^{13}$CO, CO, [CII], and [OI] emission at $<10$" ($<3$ pc) scales. We estimate Av using far-IR thermal continuum emission from dust and find the CO/[CII] ratios reach the Milky Way value at high A$_{V}$ in the centers of the clouds and fall to $\sim{1/5-1/10}\times$ the Milky Way value in the outskirts, indicating the presence of translucent molecular gas not traced by bright CO emission. We estimate the amount of molecular gas traced by bright [CII] emission at low A$_{V}$ and bright CO emission at high A$_{V}$. We find that most of the molecular gas is at low A$_{V}$ and traced by bright [CII] emission, but that faint CO emission appears to extend to where we estimate the H$_{2}$-to-HI transition occurs. By converting our H$_{2}$ gas estimates to a CO-to-H$_{2}$ conversion factor ($X_{CO}$), we show that $X_{CO}$ is primarily a function of A$_{V}$, consistent with simulations and models of low metallicity molecular clouds.

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.

We present JHK$_{s}$ observations of the metal-poor ([Fe/H] $<$ -1.40) Dwarf-irregular galaxies, Leo A and Sextans A obtained with the WIYN High-Resolution Infrared Camera at Kitt Peak. Their near-IR stellar populations are characterized by using a combination of colour-magnitude diagrams and by identifying long-period variable stars. We detected red giant and asymptotic giant branch stars, consistent with membership of the galaxy's intermediate-age populations (2-8 Gyr old). Matching our data to broadband optical and mid-IR photometry we determine luminosities, temperatures and dust-production rates (DPR) for each star. We identify 32 stars in Leo A and 101 stars in Sextans A with a DPR $>10^{-11}$ $M_\odot \,{\rm yr}^{-1}$, confirming that metal-poor stars can form substantial amounts of dust. We also find tentative evidence for oxygen-rich dust formation at low metallicity, contradicting previous models that suggest oxygen-rich dust production is inhibited in metal-poor environments. The total rates of dust injection into the interstellar medium of Leo A and Sextans A are (8.2 $\pm$ 1.8) $\times 10^{-9}$ $M_\odot \,{\rm yr}^{-1}$ and (6.2 $\pm$ 0.2) $\times 10^{-7}$ $M_\odot \,{\rm yr}^{-1}$, respectively. The majority of this dust is produced by a few very dusty evolved stars, and does not vary strongly with metallicity.

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.

Within the early Universe, `extreme' star formation may have been the norm rather than the exception. Super Star Clusters (SSCs; $M_\star$ $\gtrsim$ 10$^5$ M$_\odot$) are thought to be the modern-day analogs of globular clusters, relics of a cosmic time ($z$ $\gtrsim$ 2) when the Universe was filled with vigorously star-forming systems. The giant HII region 30 Doradus in the Large Magellanic Cloud (LMC) is often regarded as a benchmark for studies of extreme star formation. Here, we report the discovery of a massive embedded star forming complex spanning $\sim$ 500 pc in the unexplored southwest region of the LMC, which manifests itself as a younger, embedded twin of 30 Doradus. Previously known as N79, this region has a star formation efficiency exceeding that of 30 Doradus by a factor of $\sim$ 2 as measured over the past $\lesssim$ 0.5 Myr. Moreover, at the heart of N79 lies the most luminous infrared (IR) compact source discovered with large-scale IR surveys of the LMC and Milky Way, possibly a precursor to the central SSC of 30 Doradus, R136. The discovery of a nearby candidate SSC may provide invaluable information to understand how extreme star formation proceeds in the current and high-redshift Universe.

We have performed Atacama Large Millimeter/submillimeter Array (ALMA) observations in $^{12}$CO($J=2-1$), $^{13}$CO($J=2-1$), C$^{18}$O($J=2-1$), $^{12}$CO($J=3-2$), $^{13}$CO($J=3-2$), and CS($J=7-6$) lines toward the active star-forming region N83C in the Small Magellanic Cloud (SMC), whose metallicity is $\sim$ 1/5 of the Milky Way (MW). The ALMA observations first reveal sub-pc scale molecular structures in $^{12}$CO($J=2-1$) and $^{13}$CO($J=2-1$) emission. We found strong CO peaks associated with young stellar objects (YSOs) identified by the $Spitzer$ Space Telescope, and also found that overall molecular gas is distributed along the edge of the neighboring HII region. We derived a gas density of $\sim 10^4$ cm$^{-3}$ in molecular clouds associated with YSOs based on the virial mass estimated from $^{12}$CO($J=2-1$) emission. This high gas density is presumably due to the effect of the HII region under the low-metallicity (accordingly small-dust content) environment in the SMC; far-UV radiation from the HII region can easily penetrate and photo-dissociate the outer layer of $^{12}$CO molecules in the molecular clouds, and thus only the innermost parts of the molecular clouds are observed even in $^{12}$CO emission. We obtained the CO-to-H$_2$ conversion factor $X_{\rm CO}$ of $7.5 \times 10^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$ in N83C based on virial masses and CO luminosities, which is four times larger than that in the MW, 2 $\times 10^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$. We also discuss the difference in the nature between two high-mass YSOs, each of which is associated with a molecular clump with a mass of about a few $\times 10^3 M_{\odot}$.

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.

The Infrared Spectrograph (IRS) on the {\em Spitzer Space Telescope} observed nearly 800 point sources in the Large Magellanic Cloud (LMC), taking over 1,000 spectra. 197 of these targets were observed as part of the Sage-Spec Spitzer Legacy program; the remainder are from a variety of different calibration, guaranteed time and open time projects. We classify these point sources into types according to their infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership, and variability information, using a decision-tree classification method. We then refine the classification using supplementary information from the astrophysical literature. We find that our IRS sample is comprised substantially of YSO and H\,{\sc ii} regions, post-Main Sequence low-mass stars: (post-)AGB stars and planetary nebulae and massive stars including several rare evolutionary types. Two supernova remnants, a nova and several background galaxies were also observed. We use these classifications to improve our understanding of the stellar populations in the Large Magellanic Cloud, study the composition and characteristics of dust species in a variety of LMC objects, and to verify the photometric classification methods used by mid-IR surveys. We discover that some widely-used catalogues of objects contain considerable contamination and others are missing sources in our sample.

Galactic star formation scaling relations show increased scatter from kpc to sub-kpc scales. Investigating this scatter may hold important clues to how the star formation process evolves in time and space. Here, we combine different molecular gas tracers, different star formation indicators probing distinct populations of massive stars, and knowledge on the evolutionary state of each star forming region to derive star formation properties of $\sim$150 star forming complexes over the face of the Large Magellanic Cloud. We find that the rate of massive star formation ramps up when stellar clusters emerge and boost the formation of subsequent generations of massive stars. In addition, we reveal that the star formation efficiency of individual GMCs declines with increasing cloud gas mass ($M_\mathrm{cloud}$). This trend persists in Galactic star forming regions, and implies higher molecular gas depletion times for larger GMCs. We compare the star formation efficiency per freefall time ($\epsilon_\mathrm{ff}$) with predictions from various widely-used analytical star formation models. We show that while these models can produce large dispersions in $\epsilon_\mathrm{ff}$ similar to observations, the origin of the model-predicted scatter is inconsistent with observations. Moreover, all models fail to reproduce the observed decline of $\epsilon_\mathrm{ff}$ with increasing $M_\mathrm{cloud}$ in the LMC and the Milky Way. We conclude that analytical star formation models idealizing global turbulence levels, cloud densities, and assuming a stationary SFR are inconsistent with observations from modern datasets tracing massive star formation on individual cloud scales. Instead, we reiterate the importance of local stellar feedback in shaping the properties of GMCs and setting their massive star formation rate.

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.

The Mid-Infrared Instrument (MIRI) for the {\em James Webb Space Telescope} (JWST) will revolutionize our understanding of infrared stellar populations in the Local Volume. Using the rich {\em Spitzer}-IRS spectroscopic data-set and spectral classifications from the Surveying the Agents of Galaxy Evolution (SAGE)-Spectroscopic survey of over a thousand objects in the Magellanic Clouds, the Grid of Red supergiant and Asymptotic giant branch star ModelS ({\sc grams}), and the grid of YSO models by Robitaille et al. (2006), we calculate the expected flux-densities and colors in the MIRI broadband filters for prominent infrared stellar populations. We use these fluxes to explore the {\em JWST}/MIRI colours and magnitudes for composite stellar population studies of Local Volume galaxies. MIRI colour classification schemes are presented; these diagrams provide a powerful means of identifying young stellar objects, evolved stars and extragalactic background galaxies in Local Volume galaxies with a high degree of confidence. Finally, we examine which filter combinations are best for selecting populations of sources based on their JWST colours.

We present an updated general purpose nuclear equation of state (EoS) for use in simulations of core-collapse supernovae, neutron star mergers and black hole collapse. This EoS is formulated in the context of Density Functional Theory (DFT) and is generalized to include all DFT EoSs consistent with known nuclear and astrophysical constraints. This EoS also allows for the possibility of the formation of material with a net proton excess ($Y_p > 0.5$) and has an improved treatment of the nuclear statistical equilibrium and the transition to heavy nuclei as the density approaches nuclear matter density. We include the effects of pions in the regime above nuclear matter density and incorporate all of the known mesonic and baryonic states at high temperature. We analyze how a 3-body nuclear force term in the DFT at high densities stiffens the EoS to satisfy the maximum neutron star constraint, however the density dependence of the symmetry anergy and the formation of pions at high temperatures allows for a softening of the central core in supernova collapse calculations leading to a robust explosion. We also add the possibility of a transition to a QCD chiral-symmetry-restoration and deconfinement phase at densities above nuclear matter density. This paper details the physics, and constraints on, this new EoS and presents an illustration of its implementation in both neutron stars and core-collapse supernova simulations. We present the first results from core-collapse supernova simulations with this EoS.

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.

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.

We present $\mathrm{^{12}CO}$ and $\mathrm{^{13}CO}$ molecular gas data observed by ALMA, massive early stage young stellar objects identified by applying color-magnitude cuts to \textit{Spitzer} and \textit{Herschel} photometry, and low-mass late stage young stellar objects identified via H$\mathrm{\alpha}$ excess. Using dendrograms, we derive properties for the molecular cloud structures. This is the first time a dendrogram analysis has been applied to extragalactic clouds. The majority of clumps have a virial parameter equal to unity or less. The size-linewidth relations of $\mathrm{^{12}CO}$ and $\mathrm{^{13}CO}$ show the clumps in this study have a larger linewidth for a given size (by factor of 3.8 and 2.5, respectively) in comparison to several, but not all, previous studies. The larger linewidths in 30 Doradus compared to typical Milky Way quiescent clumps are probably due to the highly energetic environmental conditions of 30 Doradus. The slope of the size-linewidth relations of $\mathrm{^{12}CO}$, 0.65 $\pm$ 0.04, and $\mathrm{^{13}CO}$, 0.97 $\pm$ 0.12, are on the higher end but consistent within 3$\mathrm{\sigma}$ of previous studies. Massive star formation occurs in clumps with high masses ($> 1.83 \times 10^{2}\;\mathrm{M_{\odot}}$), high linewidths (v $> 1.18\;\mathrm{km/s}$), and high mass densities ($> 6.67 \times 10^{2}\;\mathrm{M_{\odot}\;pc^{-2}}$). The majority of embedded, massive young stellar objects are associated with a clump. However the majority of more evolved, low-mass young stellar objects are not associated with a clump.

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.

Using the PACS and SPIRE spectrometers on-board the Herschel Space Observatory, we obtained spectra of two red supergiants (RSGs) in the Large Magellanic Cloud (LMC). Multiple rotational CO emission lines (J=6-5 to 15-14) and 15 H2O lines were detected from IRAS 05280-6910, and one CO line was detected from WOH G64. This is the first time CO rotational lines have been detected from evolved stars in the LMC. Their CO line intensities are as strong as those of the Galactic RSG, VY CMa. Modelling the CO lines and the spectral energy distribution results in an estimated mass-loss rate for IRAS 05280-6910 of 3x10^-4 Msun per yr. The model assumes a gas-to-dust ratio and a CO-to-H2 abundance ratio is estimated from the Galactic values scaled by the LMC metallicity ([Fe/H]~-0.3), i.e., that the CO-to-dust ratio is constant for Galactic and LMC metallicities within the uncertainties of the model. The key factor determining the CO line intensities and the mass-loss rate found to be the stellar luminosity.

We summarize the impact of sterile neutrino dark matter on core-collapse supernova explosions. We explore various oscillations between electron neutrinos or mixed $\mu-\tau$ neutrinos and right-handed sterile neutrinos that may occur within a core-collapse supernova. In particular, we consider sterile neutrino masses and mixing angles that are consistent with sterile neutrino dark matter candidates as indicated by recent X-ray flux measurements. We find that the interpretation of the observed 3.5 keV X-ray excess as due to a decaying 7 keV sterile neutrino that comprises 100\% of the dark matter would have almost no observable effect on supernova explosions. However, in the more realistic case in which the decaying sterile neutrino comprises only a small fraction of the total dark matter density due to the presence of other sterile neutrino flavors, WIMPs, etc., a larger mixing angle is allowed. In this case a 7 keV sterile neutrino could have a significant impact on core-collapse supernovae. We also consider mixing between $\mu-\tau$ neutrinos and sterile neutrinos. We find, however, that this mixing does not significantly alter the explosion and has no observable effect on the neutrino luminosities at early times.

We present Herschel SPIRE Fourier Transform Spectrometer (FTS) observations of N159W, an active star-forming region in the Large Magellanic Cloud (LMC). In our observations, a number of far-infrared cooling lines including CO(4-3) to CO(12-11), [CI] 609 and 370 micron, and [NII] 205 micron are clearly detected. With an aim of investigating the physical conditions and excitation processes of molecular gas, we first construct CO spectral line energy distributions (SLEDs) on 10 pc scales by combining the FTS CO transitions with ground-based low-J CO data and analyze the observed CO SLEDs using non-LTE radiative transfer models. We find that the CO-traced molecular gas in N159W is warm (kinetic temperature of 153-754 K) and moderately dense (H2 number density of (1.1-4.5)e3 cm-3). To assess the impact of the energetic processes in the interstellar medium on the physical conditions of the CO-emitting gas, we then compare the observed CO line intensities with the models of photodissociation regions (PDRs) and shocks. We first constrain the properties of PDRs by modelling Herschel observations of [OI] 145, [CII] 158, and [CI] 370 micron fine-structure lines and find that the constrained PDR components emit very weak CO emission. X-rays and cosmic-rays are also found to provide a negligible contribution to the CO emission, essentially ruling out ionizing sources (ultraviolet photons, X-rays, and cosmic-rays) as the dominant heating source for CO in N159W. On the other hand, mechanical heating by low-velocity C-type shocks with ~10 km/s appears sufficient enough to reproduce the observed warm CO.

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 present the ALMA Band 3 and Band 6 results of 12CO(2-1), 13$CO(2-1), H30alpha recombination line, free-free emission around 98 GHz, and the dust thermal emission around 230 GHz toward the N159 East Giant Molecular Cloud (N159E) in the Large Magellanic Cloud (LMC). LMC is the nearest active high-mass star forming face-on galaxy at a distance of 50 kpc and is the best target for studing high-mass star formation. ALMA observations show that N159E is the complex of filamentary clouds with the width and length of ~1 pc and 5 pc - 10 pc, respectively. The total molecular mass is 0.92 x 10^5 Msun from the 13CO(2-1) intensity. N159E harbors the well-known Papillon Nebula, a compact high-excitation HII region. We found that a YSO associated with the Papillon Nebula has the mass of 35 Msun and is located at the intersection of three filamentary clouds. It indicates that the formation of the high-mass YSO was induced by the collision of filamentary clouds. Fukui et al. 2015 reported a similar kinematic structure toward a YSO in the N159 West region which is another YSO that has the mass larger than 35 Msun in these two regions. This suggests that the collision of filamentary clouds is a primary mechanism of high-mass star formation. We found a small molecular hole around the YSO in Papillon Nebula with sub-pc scale. It is filled by free-free and H30alpha emission. Temperature of the molecular gas around the hole reaches ~ 80 K. It indicates that this YSO has just started the distruction of parental molecular cloud.

The lifecycle of dust in the interstellar medium (ISM) is heavily influenced by outflows from asymptotic giant branch (AGB) and red supergiant (RSG) stars, a large fraction of which is contributed by a few very dusty sources. We compute the dust input to the Small Magellanic Cloud (SMC) by fitting the multi-epoch mid-infrared spectral energy distributions (SEDs) of AGB/RSG candidates with models from the {\em G}rid of {\em R}SG and {\em A}GB {\em M}odel{\em S} (GRAMS) grid, allowing us to estimate the luminosities and dust-production rates (DPRs) of the entire population. By removing contaminants, we guarantee a high-quality dataset with reliable DPRs and a complete inventory of the dustiest sources. We find a global AGB/RSG dust-injection rate of $(1.3\pm 0.1)\times 10^{-6}$ \msunperyr, in agreement with estimates derived from mid-infrared colours and excess fluxes. As in the LMC, a majority (66\%) of the dust arises from the extreme AGB stars, which comprise only $\approx$7\% of our sample. A handful of far-infrared sources, whose 24 \mic\ fluxes exceed their 8 \mic\ fluxes, dominate the dust input. Their inclusion boosts the global DPR by $\approx$1.5$\times$, making it necessary to determine whether they are AGB stars. Model assumptions, rather than missing data, are the major sources of uncertainty; depending on the choice of dust shell expansion speed and dust optical constants, the global DPR can be up to $\approx$10 times higher. Our results suggest a non-stellar origin for the SMC dust, barring as yet undiscovered evolved stars with very high DPRs.

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.

Using observations from the {\em Herschel} Inventory of The Agents of Galaxy Evolution (HERITAGE) survey of the Magellanic Clouds, we have found thirty five evolved stars and stellar end products that are bright in the far-infrared. These twenty eight (LMC) and seven (SMC) sources were selected from the 529 evolved star candidates in the HERITAGE far-infrared point source catalogs. Our source identification method is based on spectral confirmation, spectral energy distribution characteristics, careful examination of the multiwavelength images and includes constraints on the luminosity, resulting in a thoroughly vetted list of evolved stars. These sources span a wide range in luminosity and hence initial mass. We found thirteen low- to intermediate mass evolved stars, including asymptotic giant branch (AGB) stars, post-AGB stars, planetary nebulae and a symbiotic star. We also identify ten high mass stars, including four of the fifteen known B[e] stars in the Magellanic Clouds, three extreme red supergiants which are highly enshrouded by dust, a Luminous Blue Variable, a Wolf-Rayet star and two supernova remnants. Further, we report the detection of nine probable evolved objects which were previously undescribed in the literature. These sources are likely to be among the dustiest evolved objects in the Magellanic Clouds. The {\em Herschel} emission may either be due to dust produced by the evolved star or it may arise from swept-up ISM material.

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.

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.

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.

We have carried out 13CO(J=2-1) observations of the active star-forming region N159 West in the LMC with ALMA. We have found that the CO distribution at a sub-pc scale is highly elongated with a small width. These elongated clouds called "filaments" show straight or curved distributions with a typical width of 0.5-1.0pc and a length of 5-10pc. All the known infrared YSOs are located toward the filaments. We have found broad CO wings of two molecular outflows toward young high-mass stars in N159W-N and N159W-S, whose dynamical timescale is ~10^4 yrs. This is the first discovery of protostellar outflow in external galaxies. For N159W-S which is located toward an intersection of two filaments we set up a hypothesis that the two filaments collided with each other ~10^5 yrs ago and triggered formation of the high-mass star having ~37 Mo. The colliding clouds show significant enhancement in linewidth in the intersection, suggesting excitation of turbulence in the shocked interface layer between them as is consistent with the magneto-hydro-dynamical numerical simulations (Inoue & Fukui 2013). This turbulence increases the mass accretion rate to ~6x10^-4 Mo yr^-1, which is required to overcome the stellar feedback to form the high-mass star.

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.

We present a study of the composition of gas and dust in the Large and Small Magellanic Clouds (LMC and SMC, together -- the MCs) as measured by UV absorption spectroscopy. We have measured P II and Fe II along 85 sightlines toward the MCs using archival FUSE observations. For 16 of those sightlines, we have measured Si II, Cr II, and Zn II from new HST COS observations. We have combined these measurements with H I and H$_2$ column densities and reference stellar abundances from the literature to derive gas-phase abundances, depletions, and gas-to-dust ratios (GDRs). 80 of our 84 P measurements and 13 of our 16 Zn measurements are depleted by more than 0.1 decades, suggesting that P and Zn abundances are not accurate metallicity indicators at and above the metallicity of the SMC. The maximum P and Zn depletions are the same in the MW, LMC, and SMC. Si, Cr, and Fe are systematically less depleted in the SMC than in the MW or LMC. The minimum Si depletion in the SMC is consistent with zero. Our depletion-derived GDRs broadly agree with GDRs from the literature. The GDR varies from location to location within a galaxy by a factor of up to 2 in the LMC and up to 5 in the SMC. This variation is evidence of dust destruction and/or growth in the diffuse neutral phase of the interstellar medium.

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

We present our study on the infrared variability of point sources in the Small Magellanic Cloud (SMC). We use the data from the Spitzer Space Telescope Legacy Program "Surveying the Agents of Galaxy Evolution in the Tidally Stripped, Low Metallicity Small Magellanic Cloud" (SAGE-SMC) and the "Spitzer Survey of the Small Magellanic Cloud" (S$^{3}$MC) survey, over three different epochs, separated by several months to three years. Variability in the thermal infrared is identified using a combination of Spitzer's IRAC 3.6, 4.5, 5.8, and 8.0 $\mu$m bands, and the MIPS 24 $\mu$m band. An error-weighted flux difference between each pair of three epochs ("variability index") is used to assess the variability of each source. A visual source inspection is used to validate the photometry and image quality. Out of $\sim$2 million sources in the SAGE-SMC catalog, 814 meet our variability criteria. We matched the list of variable star candidates to the catalogs of SMC sources classified with other methods, available in the literature. Carbon-rich Asymptotic Giant Branch (AGB) stars make up the majority (61%) of our variable sources, with about a third of all of our sources being classified as extreme AGB stars. We find a small, but significant population of oxygen-rich AGB (8.6%), Red Supergiant (2.8%), and Red Giant Branch (<1%) stars. Other matches to the literature include Cepheid variable stars (8.6%), early-type stars (2.8%), young-stellar objects (5.8%), and background galaxies (1.2%). We found a candidate OH maser star, SSTISAGE1C J005212.88-730852.8, which is a variable O-rich AGB star, and would be the first OH/IR star in the SMC, if confirmed. We measured the infrared variability of a rare RV Tau variable (a post-AGB star) that has recently left the AGB phase. Fifty nine variable stars from our list remain unclassified.

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.

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.

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.

The nature, composition, abundance, and size distribution of dust in galaxies is determined by the rate at which it is created in the different stellar sources and destroyed by interstellar shocks. Because of their extensive wavelength coverage, proximity, and nearly face-on geometry, the Magellanic Clouds (MCs) provide a unique opportunity to study these processes in great detail. In this paper we use the complete sample of supernova remnants (SNRs) in the MCs to calculate the lifetime and destruction efficiencies of silicate and carbon dust in these galaxies. We find dust lifetimes of 22 +- 13 Myr (30 +- 17 Myr) for silicate (carbon) grains in the LMC, and 54 +- 32 Myr (72 +- 43 Myr) for silicate (carbon) grains in the SMC. The significantly shorter lifetimes in the MCs, as compared to the Milky Way, are explained as the combined effect of their lower total dust mass, and the fact that the dust-destroying isolated SNe in the MCs seem to be preferentially occurring in regions with higher than average dust-to-gas (D2G) mass ratios. We also calculate the supernova rate and the current star formation rate in the MCs, and use them to derive maximum dust injection rates by asymptotic giant branch stars and core collapse supernovae. We find that the injection rates are an order of magnitude lower than the dust destruction rates by the SNRs. This supports the conclusion that, unless the dust destruction rates have been considerably overestimated, most of the dust must be reconstituted from surviving grains in dense molecular clouds. More generally, we also discuss the dependence of the dust destruction rate on the local D2G mass ratio, the ambient gas density and metallicity, as well as the application of our results to other galaxies and dust evolution models.

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

We present a comprehensive spectral and morphological analysis of the remnant of Supernova (SN) 1987A with the Australia Telescope Compact Array (ATCA) and the Atacama Large Millimeter/submillimeter Array (ALMA). The non-thermal and thermal components of the radio emission are investigated in images from 94 to 672 GHz ($\lambda$ 3.2 mm to 450 $\mu$m), with the assistance of a high-resolution 44 GHz synchrotron template from the ATCA, and a dust template from ALMA observations at 672 GHz. An analysis of the emission distribution over the equatorial ring in images from 44 to 345 GHz highlights a gradual decrease of the east-to-west asymmetry ratio with frequency. We attribute this to the shorter synchrotron lifetime at high frequencies. Across the transition from radio to far infrared, both the synchrotron/dust-subtracted images and the spectral energy distribution (SED) suggest additional emission beside the main synchrotron component ($S_{\nu}\propto\nu^{-0.73}$) and the thermal component originating from dust grains at $T\sim22$ K. This excess could be due to free-free flux or emission from grains of colder dust. However, a second flat-spectrum synchrotron component appears to better fit the SED, implying that the emission could be attributed to a pulsar wind nebula (PWN). The residual emission is mainly localised west of the SN site, as the spectral analysis yields $-0.4\lesssim\alpha\lesssim-0.1$ across the western regions, with $\alpha\sim0$ around the central region. If there is a PWN in the remnant interior, these data suggest that the pulsar may be offset westward from the SN position.

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

We have made core-collapse supernova simulations that allow oscillations between electron neutrinos (or their anti particles) with right-handed sterile neutrinos. We have considered a range of mixing angles and sterile neutrino masses including those consistent with sterile neutrinos as a dark matter candidate. We examine whether such oscillations can impact the core bounce and shock reheating in supernovae. We identify the optimum ranges of mixing angles and masses that can dramatically enhance the supernova explosion by efficiently transporting electron anti-neutrinos from the core to behind the shock where they provide additional heating leading to much larger explosion kinetic energies. We show that this effect can cause stars to explode that otherwise would have collapsed. We find that an interesting periodicity in the neutrino luminosity develops due to a cycle of depletion of the neutrino density by conversion to sterile neutrinos that shuts off the conversion, followed by a replenished neutrino density as neutrinos transport through the core.

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.

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.

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.

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.

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

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

This paper reports variations of polycyclic aromatic hydrocarbons (PAHs) features that were found in Spitzer Space Telescope spectra of carbon-rich post-asymptotic giant branch (post-AGB) stars in the Large Magellanic Cloud (LMC). The paper consists of two parts. The first part describes our Spitzer spectral observing programme of 24 stars including post-AGB candidates. The latter half of this paper presents the analysis of PAH features in 20 carbon-rich post-AGB stars in the LMC, assembled from the Spitzer archive as well as from our own programme. We found that five post-AGB stars showed a broad feature with a peak at 7.7 micron, that had not been classified before. Further, the 10--13 micron PAH spectra were classified into four classes, one of which has three broad peaks at 11.3, 12.3 and 13.3 micron rather than two distinct sharp peaks at 11.3 and 12.7 micron, as commonly found in HII regions. Our studies suggest that PAHs are gradually processed while the central stars evolve from post-AGB phase to PNe, changing their composition before PAHs are incorporated into the interstellar medium. Although some metallicity dependence of PAH spectra exists, the evolutionary state of an object is more significant than its metallicity in determining the spectral characteristics of PAHs for LMC and Galactic post-AGB stars.

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.

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.

We present ALMA observations of 30 Doradus -- the highest resolution view of molecular gas in an extragalactic star formation region to date (~0.4pc x 0.6pc). The 30Dor-10 cloud north of R136 was mapped in 12CO 2-1, 13CO 2-1, C18O 2-1, 1.3mm continuum, the H30alpha recombination line, and two H2CO 3-2 transitions. Most 12CO emission is associated with small filaments and clumps (<1pc, ~1000 Msun at the current resolution). Some clumps are associated with protostars, including "pillars of creation" photoablated by intense radiation from R136. Emission from molecular clouds is often analyzed by decomposition into approximately beam-sized clumps. Such clumps in 30 Doradus follow similar trends in size, linewidth, and surface density to Milky Way clumps. The 30 Doradus clumps have somewhat larger linewidths for a given size than predicted by Larson's scaling relation, consistent with pressure confinement. They extend to higher surface density at a given size and linewidth compared to clouds studied at 10pc resolution. These trends are also true of clumps in Galactic infrared-dark clouds; higher resolution observations of both environments are required. Consistency of clump masses calculated from dust continuum, CO, and the virial theorem reveals that the CO abundance in 30 Doradus clumps is not significantly different from the LMC mean, but the dust abundance may be reduced by ~2. There are no strong trends in clump properties with distance from R136; dense clumps are not strongly affected by the external radiation field, but there is a modest trend towards lower dense clump filling fraction deeper in the cloud.

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.

We present an updated and improved equation of state (which we call the NDL EoS) for use in neutron-star structure and supernova simulations. This EoS is based upon a framework originally developed by Bowers & Wilson, but there are numerous changes. Among them are: (1) a reformulation in the context of density functional theory; (2) the possibility of the formation of material with a net proton excess (Ye > 0.5); (3) an improved treatment of the nuclear statistical equilibrium and the transition to heavy nuclei as the density approaches nuclear matter density; (4) an improved treatment of the effects of pions in the regime above nuclear matter density including the incorporation of all the known mesonic and baryonic states at high temperature; (5) the effects of 3-body nuclear forces at high densities; and (6) the possibility of a first-order or crossover transition to a QCD chiral symmetry restoration and deconfinement phase at densities above nuclear matter density. This paper details the physics of, and constraints on, this new EoS and describes its implementation in numerical simulations. We show comparisons of this EoS with other equations of state commonly used in supernova collapse simulations.

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.

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.

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.

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 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)

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.

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.

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.

We report optical and mid-infrared photometry of SN 1980K between 2004 and 2010, which show slow monotonic fading consistent with previous spectroscopic and photometric observations made 8 to 17 years after outburst. The slow rate-of-change over two decades suggests that this evolution may result from scattered and thermal light echoes off of extended circumstellar material. We present a semi- analytic dust radiative-transfer model that uses an empirically corrected effective optical depth to provide a fast and robust alternative to full Monte-Carlo radiative transfer modeling for homogenous dust at low to intermediate optical depths. We find that unresolved echoes from a thin circumstellar shell 14-15 lt-yr from the progenitor, and containing about 0.02 Msun of carbon-rich dust, can explain the broadband spectral and temporal evolution. The size, mass and dust composition are in good agreement with the contact discontinuity observed in scattered echoes around SN 1987A. The origin of slowly-changing high-velocity [O I] and Halpha lines is also considered. We propose an origin in shocked high-velocity metal-rich clumps of ejecta, rather than arising in the impact of ejecta on slowly-moving circumstellar material, as is the case with hot spots in SN 1987A.

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.

We present BVRIJHK band photometry of 6 core-collapse supernovae, SNe 1999bw, 2002hh, 2003gd, 2004et, 2005cs, and 2006bc measured at late epochs (>2 yrs) based on Hubble Space Telescope (HST), Gemini north, and WIYN telescopes. We also show the JHK lightcurves of a supernova impostor SN 2008S up to day 575. Of our 43 HST observations in total, 36 observations are successful in detecting the light from the SNe alone and measuring magnitudes of all the targets. HST observations show a resolved scattered light echo around SN 2003gd at day 1520 and around SN 2002hh at day 1717. Our Gemini and WIYN observations detected SNe 2002hh and 2004et, as well. Combining our data with previously published data, we show VRIJHK-band lightcurves and estimate decline magnitude rates at each band in 4 different phases. Our prior work on these lightcurves and other data indicate that dust is forming in our targets from day ~300-400, supporting SN dust formation theory. In this paper we focus on other physical properties derived from the late time light curves. We estimate 56Ni masses for our targets (0.5-14 x 10^{-2} Msun) from the bolometric lightcurve of each for days ~150-300 using SN 1987A as a standard (7.5 x 10^{-2} Msun). The flattening or sometimes increasing fluxes in the late time light curves of SNe 2002hh, 2003gd, 2004et and 2006bc indicate the presence of light echos. We estimate the circumstellar hydrogen density of the material causing the light echo and find that SN 2002hh is surrounded by relatively dense materials (n(H) >400 cm^{-3}) and SNe 2003gd and 2004et have densities more typical of the interstellar medium (~1 cm^{-3}). The 56Ni mass appears well correlated with progenitor mass with a slope of 0.31 x 10^{-2}, supporting the previous work by Maeda et al. (2010), who focus on more massive Type II SNe. The dust mass does not appear to be correlated with progenitor mass.

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.

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.

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.

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

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.

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

The total dust return rate from AGB and RSG star outflows is an important constraint to galactic chemical evolution models. However, this requires detailed radiative transfer (RT) modeling of individual stars, which becomes impractical for large data sets. Another approach is to select the best-fit spectral energy distribution (SED) from a grid of dust shell models, allowing for a faster determination of the luminosities and mass-loss rates for entire samples. We have developed the Grid of RSG and AGB ModelS (GRAMS) to measure the mass-loss return from evolved stars. The models span the range of stellar, dust shell and grain properties relevant to evolved stars. In this paper we present the carbon-star grid and compare our results with data of Large Magellanic Cloud (LMC) carbon stars from the SAGE and SAGE-Spec programs. We generate spherically symmetric dust shell models using the 2Dust code, with hydrostatic models for the central stars. We explore five values of the inner radius R_in of the dust shell (1.5, 3, 4.5, 7 and 12 R_star). We use amorphous carbon dust mixed with 10% silicon carbide by mass. The grain sizes follows a KMH distribution. The models span 26 values of 11.3 um optical depth, ranging from 0.001 to 4. For each model, 2Dust calculates the output SED from 0.2 to 200 um. Over 12,000 models have dust temperatures below 1800 K. The GRAMS synthetic photometry is in good agreement with SAGE photometry for LMC carbon-rich and extreme AGB star candidates, as well as spectroscopically confirmed carbon stars from the SAGE-Spec study. Our models reproduce the IRAC colors of most of the extreme AGB star candidates, consistent with the expectation that a majority of these enshrouded stars have carbon-rich dust. Finally, we fit the SEDs of some well-studied carbon stars and compare the resulting luminosities and mass-loss rates with those from previous studies.

We present new Spitzer/IRS spectra of two hot R Coronae Borealis (RCB) stars, one in the Galaxy,V348 Sgr, and one lying in the LMC, HV 2671. These two objects may constitute a link between the RCB stars and the late Wolf-Rayet ([WCL]) class of central stars of planetary nebula (CSPNe) such as CPD -56 8032 that has little or no hydrogen in their atmospheres. HV 2671 and V348 Sgr are members of a rare subclass that has significantly higher effective temperatures than most RCB stars, but sharing the traits of hydrogen deficiency and dust formation that define the cooler RCB stars. The [WC] CSPNe star, CPD -56 8032, displays evidence for dual-dust chemistry showing both PAHs and crystalline silicates in its mid-IR spectrum. HV 2671 shows strong PAH emission but shows no sign of having crystalline silicates. The spectrum of V348 Sgr is very different from those of CPD -56 8032 and HV 2671. The PAH emission seen strongly in the other two stars is not present. Instead, the spectrum is dominated by a broad emission centered at about 8.2 micron. The mid-IR spectrum of CPD -56 8032 shows emission features that may be associated with C60. The other two stars do not show evidence for C60. HV 2671 has also been detected by Herschel/PACS and SPIRE. V348 Sgr and CPD -56 8032 have been detected by AKARI/FIS. These data were combined with Spitzer, IRAS, 2MASS and other photometry to produce their spectral energy distributions from the visible to the far-IR. Monte Carlo radiative transfer modeling was used to study the circumstellar dust around these stars. HV 2671 and CPD -56 8032 require both a flared inner disk with warm dust and an extended diffuse envelope with cold dust to to fit their SEDs. The SED of V348 Sgr can be fit with a much smaller disk and envelope.

SN 2010jl was an extremely bright, Type IIn SNe which showed a significant IR excess no later than 90 days after explosion. We have obtained Spitzer 3.6 and 4.5 \mum and JHK observations of SN 2010jl \sim90 days post explosion. Little to no reddening in the host galaxy indicated that the circumstellar material lost from the progenitor must lie in a torus inclined out of the plane of the sky. The likely cause of the high mid-IR flux is the reprocessing of the initial flash of the SN by pre-existing circumstellar dust. Using a 3D Monte Carlo Radiative Transfer code, we have estimated that between 0.03-0.35 Msun of dust exists in a circumstellar torus around the SN located 6 \times 10 ^17 cm away from the SN and inclined between 60-80\cdot to the plane of the sky. On day 90, we are only seeing the illumination of approximately 5% of this torus, and expect to see an elevated IR flux from this material up until day \sim 450. It is likely this dust was created in an LBV-like mass loss event of more than 3 Msun, which is large but consistent with other LBV progenitors such as {\eta} Carinae.

We present a galaxy (SAGE1CJ053634.78-722658.5) at a redshift of 0.14 of which the IR is entirely dominated by emission associated with the AGN. We present the 5-37 um Spitzer/IRS spectrum and broad wavelength SED of SAGE1CJ053634, an IR point-source detected by Spitzer/SAGE (Meixner et al 2006). The source was observed in the SAGE-Spec program (Kemper et al., 2010) and was included to determine the nature of sources with deviant IR colours. The spectrum shows a redshifted (z=0.14+-0.005) silicate emission feature with an exceptionally high feature-to-continuum ratio and weak polycyclic aromatic hydrocarbon (PAH) bands. We compare the source with models of emission from dusty tori around AGNs from Nenkova et al. (2008). We present a diagnostic diagram that will help to identify similar sources based on Spitzer/MIPS and Herschel/PACS photometry. The SED of SAGE1CJ053634 is peculiar because it lacks far-IR emission and a clear stellar counterpart. We find that the SED and the IR spectrum can be understood as emission originating from the inner ~10 pc around an accreting black hole. There is no need to invoke emission from the host galaxy, either from the stars or from the interstellar medium, although a possible early-type host galaxy cannot be excluded based on the SED analysis. The hot dust around the accretion disk gives rise to a continuum, which peaks at 4 um, whereas the strong silicate features may arise from optically thin emission of dusty clouds within ~10 pc around the black hole. The weak PAH emission does not appear to be linked to star formation, as star formation templates strongly over-predict the measured far-IR flux levels. The SED of SAGE1CJ053634 is rare in the local universe but may be more common in the more distant universe. The conspicuous absence of host-galaxy IR emission places limits on the far-IR emission arising from the dusty torus alone.

The absolute flux calibration of the James Webb Space Telescope will be based on a set of stars observed by the Hubble and Spitzer Space Telescopes. In order to cross-calibrate the two facilities, several A, G, and white dwarf (WD) stars are observed with both Spitzer and Hubble and are the prototypes for a set of JWST calibration standards. The flux calibration constants for the four Spitzer IRAC bands 1-4 are derived from these stars and are 2.3, 1.9, 2.0, and 0.5% lower than the official cold-mission IRAC calibration of Reach et al. (2005), i.e. in agreement within their estimated errors of ~2%. The causes of these differences lie primarily in the IRAC data reduction and secondarily in the SEDs of our standard stars. The independent IRAC 8 micron band-4 fluxes of Rieke et al. (2008) are about 1.5 +/- 2% higher than those of Reach et al. and are also in agreement with our 8 micron result.

We present a study of the dust properties of 12 classical and superbubble HII regions in the Large Magellanic Cloud. We use infrared photometry from Spitzer (8, 24, 70, and 160 \mum bands), obtained as part of the Surveying the Agents of a Galaxy's Evolution (SAGE) program, along with archival spectroscopic classifications of the ionizing stars to examine the role of stellar sources on dust heating and processing. Our infrared observations show surprisingly little correlation between the emission properties of the dust and the effective temperatures or bolometric magnitudes of stars in the HII regions, suggesting that the HII region evolutionary timescale is not on the order of the dust processing timescale. We find that the infrared emission of superbubbles and classical HII regions shows little differentiation between the two classes, despite the significant differences in age and morphology. We do detect a correlation of the 24 \mum emission from hot dust with the ratio of 70 to 160 \mum flux. This correlation can be modeled as a trend in the temperature of a minority hot dust component, while a majority of the dust remains significantly cooler.

SN 2007it is a bright, Type IIP supernova which shows indications of both pre-existing and newly formed dust. The visible photometry shows a bright late-time luminosity, powered by the 0.09 M$_{\sun}$ of $^{56}$Ni present in the ejecta. There is also a sudden drop in optical brightness after day 339, and a corresponding brightening in the IR due to new dust forming in the ejecta. CO and SiO emission, generally thought to be precursors to dust formation, may have been detected in the mid-IR photometry of SN 2007it. The optical spectra show stronger than average [O I] emission lines and weaker than average [Ca II] lines, which may indicate a 16 - 27 M$_{\sun}$ progenitor, on the higher end of expected Type IIP masses. Multi-component [O I] lines are also seen in the optical spectra, most likely caused by an asymmetric blob or a torus of oxygen core material being ejected during the SN explosion. Interaction with circumstellar material prior to day 540 may have created a cool dense shell between the forward and reverse shocks where new dust is condensing. At late times there is also a flattening of the visible lightcurve as the ejecta luminosity fades and a surrounding light echo becomes visible. Radiative transfer models of SN 2007it SEDs indicate that up to 10$^{-4}$ M$_{\sun}$ of new dust has formed in the ejecta, which is consistent with the amount of dust formed in other core collapse supernovae.

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.

We report the discovery of a supernova (SN) with the highest apparent energy output to date and conclude that it represents an extreme example of the Type IIn subclass. The SN, which was discovered behind the Large Magellanic Cloud at z = 0.289 by the SuperMACHO microlensing survey, peaked at M_R = -21.5 mag and only declined by 2.9 mag over 4.7 years after the peak. Over this period, SN 2003ma had an integrated bolometric luminosity of 4 x 10^51 ergs, more than any other SN to date. The radiated energy is close to the limit allowed by conventional core-collapse explosions. Optical spectra reveal that SN 2003ma has persistent single-peaked intermediate-width hydrogen lines, a signature of interaction between the SN and a dense circumstellar medium. The light curves show further evidence for circumstellar interaction, including a long plateau with a shape very similar to the classic SN IIn 1988Z -- however, SN 2003ma is ten times more luminous at all epochs. The fast velocity measured for the intermediate-width H_alpha component (~6000 km/s) points towards an extremely energetic explosion (> 10^52 ergs), which imparts a faster blast-wave speed to the post-shock material and a higher luminosity from the interaction than is observed in typical SNe IIn. Mid-infrared observations of SN 2003ma suggest an infrared light echo is produced by normal interstellar dust at a distance ~0.5 pc from the SN.

We present the photometric catalogs for the star-forming cluster NGC 602 in the wing of the Small Magellanic Cloud covering a range of wavelengths from optical HST/ACS (F555W, F814W) and SMARTS/ANDICAM (V, I) to infrared (Spitzer/IRAC 3.6, 4.5, 5.8, and 8 micron and MIPS 24 micron). Combining this with IRSF (InfraRed Survey Facility) near-infrared photometry (J, H, Ks), we compare the young main sequence (MS) and pre-main sequence (PMS) populations prominent in the optical with the current young stellar object (YSO) populations revealed by the infrared (IR). We analyze the MS and PMS population with isochrones in color-magnitude diagrams to derive ages and masses. The optical data reveal ~565 PMS candidates, low mass Stage III YSOs. We characterize ~40 YSOs by fitting their spectral energy distributions (SEDs) to a grid of models (Robitaille et al. 2007) to derive luminosities, masses and evolutionary phase (Stage I-III). The higher resolution HST images reveal that ~70% of the YSO candidates are either multiples or protoclusters. For YSOs and PMS sources found in common, we find a consistency in the masses derived. We use the YSO mass function to derive a present-day star-formation rate of ~0.2-1.0 Msun/yr/kpc^2, similar to the rate derived from the optical star formation history suggesting a constant star formation rate for this region. We demonstrate a progression of star formation from the optical star cluster center to the edge of the star forming dust cloud. We derive lifetimes of a few 10^5 years for the YSO Stages I and II.

We have estimated elemental abundances of the planetary nebula (PN) Hen2-436 in the Sagittarius (Sgr) spheroidal dwarf galaxy using ESO/VLT FORS2, Magellan/MMIRS, and Spitzer/IRS spectra. We have detected candidates of [F II] 4790A, [Kr III] 6826A, and [P II] 7875A lines and successfully estimated the abundances of these elements ([F/H]=+1.23, [Kr/H]=+0.26, [P/H]=+0.26) for the first time. We present a relation between C, F, P, and Kr abundances among PNe and C-rich stars. The detections of F and Kr support the idea that F and Kr together with C are synthesized in the same layer and brought to the surface by the third dredge-up. We have estimated the N^2+ and O^2+ abundances using optical recombination lines (ORLs) and collisionally excited lines (CELs). The discrepancy between the abundance derived from the O ORL and that derived from the O CEL is >1 dex. To investigate the status of the central star of the PN, nebula condition, and dust properties, we construct a theoretical SED model with CLOUDY. By comparing the derived luminosity and temperature of the central star with theoretical evolutionary tracks, we conclude that the initial mass of the progenitor is likely to be ~1.5-2.0 Msun and the age is ~3000 yr after the AGB phase. The observed elemental abundances can be explained by a theoretical nucleosynthesis model with a star of initial mass 2.25 Msun, Z=0.008 and LMC compositions. We have estimated the dust mass to be 2.9x10^-4 Msun (amorphous carbon only) or 4.0x10^-4 Msun (amorphous carbon and PAH). Based on the assumption that most of the observed dust is formed during the last two thermal pulses and the dust-to-gas mass ratio is 5.58x10^-3, the dust mass-loss rate and the total mass-loss rate are <3.1x10^-8 Msun / yr and <5.5x10^-6 Msun / yr, respectively. Our estimated dust mass-loss rate is comparable to a Sgr dwarf galaxy AGB star with similar metallicity and luminosity.

We present the comparison of the three most important ice constituents (water, CO and CO2) in the envelopes of massive Young Stellar Objects (YSOs), in environments of different metallicities: the Galaxy, the Large Magellanic Cloud (LMC) and, for the first time, the Small Magellanic Cloud (SMC). We present observations of water, CO and CO2 ice in 4 SMC and 3 LMC YSOs (obtained with Spitzer-IRS and VLT/ISAAC). While water and CO2 ice are detected in all Magellanic YSOs, CO ice is not detected in the SMC objects. Both CO and CO2 ice abundances are enhanced in the LMC when compared to high-luminosity Galactic YSOs. Based on the fact that both species appear to be enhanced in a consistent way, this effect is unlikely to be the result of enhanced CO2 production in hotter YSO envelopes as previously thought. Instead we propose that this results from a reduced water column density in the envelopes of LMC YSOs, a direct consequence of both the stronger UV radiation field and the reduced dust-to-gas ratio at lower metallicity. In the SMC the environmental conditions are harsher, and we observe a reduction in CO2 column density. Furthermore, the low gas-phase CO density and higher dust temperature in YSO envelopes in the SMC seem to inhibit CO freeze-out. The scenario we propose can be tested with further observations.

The objective of this paper is to construct a catalog providing the dust properties and the star formation efficiency (SFE) of the molecular clouds in the Large Magellanic Cloud (LMC). We use the infrared (IR) data obtained with the Spitzer telescope as part of the ``Surveying the Agents of a Galaxy's Evolution'' (SAGE) Legacy survey as well as the IRAS data. We also work with extinction (Av) maps of the LMC. A total of 272 molecular clouds have been detected in the LMC in a previous molecular survey, accounting for 230 giant molecular clouds and 42 smaller clouds. We perform correlations between the IR emission/extinction, and atomic and molecular gas tracers. We compare the atomic gas that surrounds the molecular cloud with the molecular gas in the cloud. Using a dust emission model, we derive the physical properties of dust in and outside the clouds: equilibrium temperature, emissivity and extinction. We also determine the luminosity of the interstellar radiation field intercepted by the cloud, and the total IR luminosity. Statistically, we do not find any significant difference in the dust properties between the atomic and the molecular phases. In particular we do not find evidence for a systematic decrease of the dust temperature in the molecular phase, with respect to the surrounding, presumably atomic gas. This is probably because giant molecular clouds are the sites of star formation, which heat the dust, while the smallest clouds are unresolved. The ratio between the IR luminosity and the cloud mass (LDust/Mgas) does not seem to correlate with Mgas. The highest value of the ratio we derived is 18.1 Lsol/Msol in the 30 Doradus region, which is known to be the most prominent star formation region of the LMC, while the most likely value is 0.5 and is representative of quiescent clouds. We provide a prescription to associate the various stages of star formation with its LDust/Mgas.

We present a study of the infrared properties of 4922 spectroscopically confirmed massive stars in the Large and Small Magellanic Clouds, focusing on the active OB star population. Besides OB stars, our sample includes yellow and red supergiants, Wolf-Rayet stars, Luminous Blue Variables (LBVs) and supergiant B[e] stars. We detect a distinct Be star sequence, displaced to the red, and find a higher fraction of Oe and Be stars among O and early-B stars in the SMC, respectively, when compared to the LMC, and that the SMC Be stars occur at higher luminosities. We also find photometric variability among the active OB population and evidence for transitions of Be stars to B stars and vice versa. We furthermore confirm the presence of dust around all the supergiant B[e] stars in our sample, finding the shape of their spectral energy distributions (SEDs) to be very similar, in contrast to the variety of SED shapes among the spectrally variable LBVs.

We present results of our study of the infrared properties of massive stars in the Large and Small Magellanic Clouds, which are based on the Spitzer SAGE surveys of these galaxies. We have compiled catalogs of spectroscopically confirmed massive stars in each galaxy, as well as photometric catalogs for a subset of these stars that have infrared counterparts in the SAGE database, with uniform photometry from 0.3 to 24 microns in the UBVIJHKs+IRAC+MIPS24 bands. These catalogs enable a comparative study of infrared excesses of OB stars, classical Be stars, yellow and red supergiants, Wolf-Rayet stars, Luminous Blue Variables and supergiant B[e] stars, as a function of metallicity, and provide the first roadmaps for interpreting luminous, massive, resolved stellar populations in nearby galaxies at infrared wavelengths.

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

We present a 2Dust model for the dust shell around a LMC long-period variable (LPV) previously studied as part of the OGLE survey. OGLE LMC LPV 28579 (SAGE J051306.40-690946.3) is a carbon-rich asymptotic giant branch (AGB) star for which we have photometry and spectra from the Spitzer SAGE and SAGE-Spec programs along with UBVIJHK_s photometry. By modeling this source, we obtain a baseline set of dust properties to be used in the construction of a grid of models for carbon stars. We reproduce its spectral energy distribution using a mixture of AmC and SiC (15% by mass). The grain sizes are distributed according to the KMH model. The best-fit model has an optical depth of 0.28 for the shell at the peak of the SiC feature, with R_in~1430 R_sun or 4.4 R_star. The temperature at this inner radius is 1310 K. Assuming an expansion velocity of 10 km s^-1, we obtain a dust mass-loss rate of 2.5x10^-9 M_sun yr-1. We calculate a 15% variation in this rate by testing the fit sensitivity against variation in input parameters. We also present a simple model for the molecular gas in the extended atmosphere that could give rise to the 13.7 \mu m feature seen in the spectrum. We find that a combination of CO and C_2H_2 gas at an excitation temperature of about 1000 K and column densities of 3x10^21 cm^-2 and 10^19 cm^-2 respectively are able to reproduce the observations. Given that the excitation temperature is close to T_dust(R_in), most of the molecular contribution probably arises from the inner shell region. The luminosity corresponding to the first epoch of SAGE observations is 6580 L_sun. For an effective temperature of about 3000 K, this implies a stellar mass of 1.5-2 M_sun and an age of 1-2.5 Gyr. For a gas:dust ratio of 200, we obtain a gas mass-loss rate of 5.0x10^-7 M_sun yr^-1, consistent with the gas mass-loss rates estimated from the period, color and 8 \mu m flux of the source.

The LMC star, SSTISAGE1C J050756.44-703453.9, was first noticed during a survey of EROS-2 lightcurves for stars with large irregular brightness variations typical of the R Coronae Borealis (RCB) class. However, the visible spectrum showing emission lines including the Balmer and Paschen series as well as many Fe II lines is emphatically not that of an RCB star. This star has all of the characteristics of a typical UX Ori star. It has a spectral type of approximately A2 and has excited an H II region in its vicinity. However, if it is an LMC member, then it is very luminous for a Herbig Ae/Be star. It shows irregular drops in brightness of up to 2 mag, and displays the reddening and "blueing" typical of this class of stars. Its spectrum, showing a combination of emission and absorption lines, is typical of a UX Ori star that is in a decline caused by obscuration from the circumstellar dust. SSTISAGE1C J050756.44-703453.9 has a strong IR excess and significant emission is present out to 500 micron. Monte Carlo radiative transfer modeling of the SED requires that SSTISAGE1C J050756.44-703453.9 has both a dusty disk as well as a large extended diffuse envelope to fit both the mid- and far-IR dust emission. This star is a new member of the UX Ori subclass of the Herbig Ae/Be stars and only the second such star to be discovered in the LMC.

The HII complex N159 in the Large Magellanic Cloud (LMC) is used to study massive star formation in different environments, as it contains three giant molecular clouds (GMCs) that have similar sizes and masses but exhibit different intensities of star formation. We identify candidate massive young stellar objects (YSOs) using infrared photometry, and model their SEDs to constrain mass and evolutionary state. Good fits are obtained for less evolved Type I, I/II, and II sources. Our analysis suggests that there are massive embedded YSOs in N159B, a maser source, and several ultracompact HII regions. Massive O-type YSOs are found in GMCs N159-E and N159-W, which are associated with ionized gas, i.e., where massive stars formed a few Myr ago. The third GMC, N159-S, has neither O-type YSOs nor evidence of previous massive star formation. This correlation between current and antecedent formation of massive stars suggests that energy feedback is relevant. We present evidence that N159-W is forming YSOs spontaneously, while collapse in N159-E may be triggered. Finally, we compare star formation rates determined from YSO counts with those from integrated H-alpha and 24 micron luminosities and expected from gas surface densities. Detailed dissection of extragalactic GMCs like the one presented here is key to revealing the physics underlying commonly used star formation scaling laws.

We combine variability information from the MAssive Compact Halo Objects (MACHO) survey of the Large Magellanic Cloud (LMC) with infrared photometry from the Spitzer Space Telescope Surveying the Agents of a Galaxy's Evolution (SAGE) survey to create a dataset of ~30 000 variable red sources. We photometrically classify these sources as being on the first ascent of the Red Giant Branch (RGB), or as being in one of three stages along the Asymptotic Giant Branch (AGB): oxygen-rich, carbon-rich, or highly reddened with indeterminate chemistry ("extreme" AGB candidates). We present linear period-luminosity relationships for these sources using 8 separate infrared bands (J, H, K, 3.6, 4.5, 5.8, 8.0, and 24 micron) as proxies for the luminosity. We find that the wavelength dependence of the slope of the period-luminosity relationship is different for different photometrically determined classes of AGB stars. Stars photometrically classified as O-rich show the least variation of slope with wavelength, while dust enshrouded extreme AGB stars show a pronounced trend toward steeper slopes with increasing wavelength. We find that O-rich AGB stars pulsating in the fundamental mode obey a period-magnitude relation with a slope of -3.41 +/- 0.04 when magnitude is measured in the 3.6 micron band, in contrast to C-rich AGB stars, which obey a relation of slope -3.77 +/- 0.05.

We present a catalog of 5324 massive stars in the Small Magellanic Cloud (SMC), with accurate spectral types compiled from the literature, and a photometric catalog for a subset of 3654 of these stars, with the goal of exploring their infrared properties. The photometric catalog consists of stars with infrared counterparts in the Spitzer, SAGE-SMC survey database, for which we present uniform photometry from 0.3-24 um in the UBVIJHKs+IRAC+MIPS24 bands. We compare the color magnitude diagrams and color-color diagrams to those of the Large Magellanic Cloud (LMC), finding that the brightest infrared sources in the SMC are also the red supergiants, supergiant B[e] (sgB[e]) stars, luminous blue variables, and Wolf-Rayet stars, with the latter exhibiting less infrared excess, the red supergiants being less dusty and the sgB[e] stars being on average less luminous. Among the objects detected at 24 um are a few very luminous hypergiants, 4 B-type stars with peculiar, flat spectral energy distributions, and all 3 known luminous blue variables. We detect a distinct Be star sequence, displaced to the red, and suggest a novel method of confirming Be star candidates photometrically. We find a higher fraction of Oe and Be stars among O and early-B stars in the SMC, respectively, when compared to the LMC, and that the SMC Be stars occur at higher luminosities. We estimate mass-loss rates for the red supergiants, confirming the correlation with luminosity even at the metallicity of the SMC. Finally, we confirm the new class of stars displaying composite A & F type spectra, the sgB[e] nature of 2dFS1804 and find the F0 supergiant 2dFS3528 to be a candidate luminous blue variable with cold dust.

Recently, we have discovered an error in our Monte-Carlo spectral fitting routine, more specifically where the errors on the fluxes were rescaled to get a reduced chi2 of 1. The rescaled errors were too big, resulting in too wide a range of good fits in our 100 step Monte-Carlo routine. This problem affects Figs. 7-9 and Tables A.1, A.2 in Gielen et al. (2008), Table 3 in Gielen et al. (2009a), and Table 4 in Gielen et al. (2009b). We corrected for this error and present the new values and errors in the tables below. The new values and errors nearly all fall within the old error range. Our best chi2 values and overall former scientific results are not affected. With these new errors some possible new trends in the dust parameters might be observed. These will be discussed in an upcoming paper where we extend the sample presented in Gielen et al. (2008) with newly obtained SPITZER-IRS data.

The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) of the Magellanic Clouds will use dust emission to investigate the life cycle of matter in both the Large and Small Magellanic Clouds (LMC and SMC). Using the Herschel Space Observatory's PACS and SPIRE photometry cameras, we imaged a 2x8 square degree strip through the LMC, at a position angle of ~22.5 degrees as part of the science demonstration phase of the Herschel mission. We present the data in all 5 Herschel bands: PACS 100 and 160 {\mu}m and SPIRE 250, 350 and 500 {\mu}m. We present two dust models that both adequately fit the spectral energy distribution for the entire strip and both reveal that the SPIRE 500 {\mu}m emission is in excess of the models by 6 to 17%. The SPIRE emission follows the distribution of the dust mass, which is derived from the model. The PAH-to-dust mass (f_PAH) image of the strip reveals a possible enhancement in the LMC bar in agreement with previous work. We compare the gas mass distribution derived from the HI 21 cm and CO J=1-0 line emission maps to the dust mass map from the models and derive gas-to-dust mass ratios (GDRs). The dust model, which uses the standard graphite and silicate optical properties for Galactic dust, has a very low GDR = 65(+15,-18) making it an unrealistic dust model for the LMC. Our second dust model, which uses amorphous carbon instead of graphite, has a flatter emissivity index in the submillimeter and results in a GDR = 287(+25,-42) that is more consistent with a GDR inferred from extinction.

Massive evolved stars can produce large amounts of dust, and far-infrared (IR) data are essential for determining the contribution of cold dust to the total dust mass. Using Herschel, we search for cold dust in three very dusty massive evolved stars in the Large Magellanic Cloud: R71 is a Luminous Blue Variable, HD36402 is a Wolf-Rayet triple system, and IRAS05280-6910 is a red supergiant. We model the spectral energy distributions using radiative transfer codes and find that these three stars have mass-loss rates up to 10^-3 solar masses/year, suggesting that high-mass stars are important contributors to the life-cycle of dust. We found far-IR excesses in two objects, but these excesses appear to be associated with ISM and star-forming regions. Cold dust (T < 100 K) may thus not be an important contributor to the dust masses of evolved stars.

We present clumps of dust emission from Herschel observations of the Large Magellanic Cloud (LMC) and their physical and statistical properties. We catalog cloud features seen in the dust emission from Herschel observations of the LMC, the Magellanic type irregular galaxy closest to the Milky Way, and compare these features with HI catalogs from the ATCA+Parkes HI survey. Using an automated cloud-finding algorithm, we identify clouds and clumps of dust emission and examine the cumulative mass distribution of the detected dust clouds. The mass of cold dust is determined from physical parameters that we derive by performing spectral energy distribution fits to 250, 350, and 500 micronm emission from SPIRE observations using DUSTY and GRASIL radiative transfer calculation with dust grain size distributions for graphite/silicate in low-metallicity extragalactic environments. The dust cloud mass spectrum follows a power law distribution with an exponent of gamma=-1.8 for clumps larger than 400 solar mass and is similar to the HI mass distribution. This is expected from the theory of ISM structure in the vicinity of star formation.

We investigate the dust associated with the supernova remnant (SNR) N49 in the Large Magellanic Cloud (LMC) as observed with the Herschel Space Observatory. N49 is unusually bright because of an interaction with a molecular cloud along its eastern edge. We have used PACS and SPIRE to measure the far IR flux densities of the entire SNR and of a bright region on the eastern edge of the SNR where the SNR shock is encountering the molecular cloud. Using these fluxes supplemented with archival data at shorter wavelengths, we estimate the dust mass associated with N49 to be about 10 Msun. The bulk of the dust in our simple two-component model has a temperature of 20-30 K, similar to that of nearby molecular clouds. Unfortunately, as a result of the limited angular resolution of Herschel at the wavelengths sampled with SPIRE, the uncertainties are fairly large. Assuming this estimate of the dust mass associated with the SNR is approximately correct, it is probable that most of the dust in the SNR arises from regions where the shock speed is too low to produce significant X-ray emission. The total amount of warm 50-60 K dust is ~0.1 or 0.4 Msun, depending on whether the dust is modeled in terms of carbonaceous or silicate grains. This provides a firm lower limit to the amount of shock heated dust in N49.

The properties of the dust grains (e.g., temperature and mass) can be derived from fitting far-IR SEDs (>100 micron). Only with SPIRE on Herschel has it been possible to get high spatial resolution at 200 to 500 micron that is beyond the peak (~160 micron) of dust emission in most galaxies. We investigate the differences in the fitted dust temperatures and masses determined using only <200 micron data and then also including >200 micron data (new SPIRE observations) to determine how important having >200 micron data is for deriving these dust properties. We fit the 100 to 350 micron observations of the Large Magellanic Cloud (LMC) point-by-point with a model that consists of a single temperature and fixed emissivity law. The data used are existing observations at 100 and 160 micron (from IRAS and Spitzer) and new SPIRE observations of 1/4 of the LMC observed for the HERITAGE Key Project as part of the Herschel Science Demonstration phase. The dust temperatures and masses computed using only 100 and 160 micron data can differ by up to 10% and 36%, respectively, from those that also include the SPIRE 250 & 350 micron data. We find that an emissivity law proportional to lambda^-1.5 minimizes the 100-350 micron fractional residuals. We find that the emission at 500 micron is ~10% higher than expected from extrapolating the fits made at shorter wavelengths. We find the fractional 500 micron excess is weakly anti-correlated with MIPS 24 micron flux and the total gas surface density. This argues against a flux calibration error as the origin of the 500 micron excess. Our results do not allow us to distinguish between a systematic variation in the wavelength dependent emissivity law or a population of very cold dust only detectable at lambda > 500 micron for the origin of the 500 micron excess.

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.

We demonstrate the unique capabilities of Herschel to study very young luminous extragalactic young stellar objects (YSOs) by analyzing a central strip of the Large Magellanic Cloud obtained through the HERITAGE Science Demonstration Program. We combine PACS 100 and 160, and SPIRE 250, 350, and 500 microns photometry with 2MASS (1.25-2.17 microns) and Spitzer IRAC and MIPS (3.6-70 microns) to construct complete spectral energy distributions (SEDs) of compact sources. From these, we identify 207 candidate embedded YSOs in the observed region, ~40% never-before identified. We discuss their position in far-infrared color-magnitude space, comparing with previously studied, spectroscopically confirmed YSOs and maser emission. All have red colors indicating massive cool envelopes and great youth. We analyze four example YSOs, determining their physical properties by fitting their SEDs with radiative transfer models. Fitting full SEDs including the Herschel data requires us to increase the size and mass of envelopes included in the models. This implies higher accretion rates (greater than or equal to 0.0001 M_sun/yr), in agreement with previous outflow studies of high-mass protostars. Our results show that Herschel provides reliable longwave SEDs of large samples of high-mass YSOs; discovers the youngest YSOs whose SEDs peak in Herschel bands; and constrains the physical properties and evolutionary stages of YSOs more precisely than was previously possible.

We study the structure of the medium surrounding sites of high-mass star formation to determine the interrelation between the HII regions and the environment from which they were formed. The density distribution of the surroundings is key in determining how the radiation of the newly formed stars interacts with the surrounds in a way that allows it to be used as a star formation tracer. We present new Herschel/SPIRE 250, 350 and 500 mum data of LHA 120-N44 and LHA 120-N63 in the LMC. We construct average spectral energy distributions (SEDs) for annuli centered on the IR bright part of the star formation sites. The annuli cover ~10-~100 pc. We use a phenomenological dust model to fit these SEDs to derive the dust column densities, characterise the incident radiation field and the abundance of polycyclic aromatic hydrocarbon molecules. We see a factor 5 decrease in the radiation field energy density as a function of radial distance around N63. N44 does not show a systematic trend. We construct a simple geometrical model to derive the 3-D density profile of the surroundings of these two regions. Herschel/SPIRE data have proven very efficient in deriving the dust mass distribution. We find that the radiation field in the two sources behaves very differently. N63 is more or less spherically symmetric and the average radiation field drops with distance. N44 shows no systematic decrease of the radiation intensity which is probably due to the inhomogeneity of the surrounding molecular material and to the complex distribution of several star forming clusters in the region.

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

The Magellanic Mopra Assessment (MAGMA) is a high angular resolution CO mapping survey of giant molecular clouds (GMCs) in the Large and Small Magellanic Clouds using the Mopra Telescope. Here we report on the basic physical properties of 125 GMCs in the LMC that have been surveyed to date. The observed clouds exhibit scaling relations that are similar to those determined for Galactic GMCs, although LMC clouds have narrower linewidths and lower CO luminosities than Galactic clouds of a similar size. The average mass surface density of the LMC clouds is 50 Msol/pc2, approximately half that of GMCs in the inner Milky Way. We compare the properties of GMCs with and without signs of massive star formation, finding that non-star-forming GMCs have lower peak CO brightness than star-forming GMCs. We compare the properties of GMCs with estimates for local interstellar conditions: specifically, we investigate the HI column density, radiation field, stellar mass surface density and the external pressure. Very few cloud properties demonstrate a clear dependence on the environment; the exceptions are significant positive correlations between i) the HI column density and the GMC velocity dispersion, ii) the stellar mass surface density and the average peak CO brightness, and iii) the stellar mass surface density and the CO surface brightness. The molecular mass surface density of GMCs without signs of massive star formation shows no dependence on the local radiation field, which is inconsistent with the photoionization-regulated star formation theory proposed by McKee (1989). We find some evidence that the mass surface density of the MAGMA clouds increases with the interstellar pressure, as proposed by Elmegreen (1989), but the detailed predictions of this model are not fulfilled once estimates for the local radiation field, metallicity and GMC envelope mass are taken into account.

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

SN 2007od exhibits characteristics that have rarely been seen in a Type IIP supernova (SN). Optical V band photometry reveals a very steep brightness decline between the plateau and nebular phases of ~4.5 mag, likely due to SN 2007od containing a low mass of 56Ni. The optical spectra show an evolution from normal Type IIP with broad Halpha emission, to a complex, four component Halpha emission profile exhibiting asymmetries caused by dust extinction after day 232. This is similar to the spectral evolution of the Type IIn SN 1998S, although no early-time narrow (~200 km s-1) Halpha component was present in SN 2007od. In both SNe, the intermediate-width Halpha emission components are thought to arise in the interaction between the ejecta and its circumstellar medium (CSM). SN 2007od also shows a mid-IR excess due to new dust. The evolution of the Halpha profile and the presence of the mid-IR excess provide strong evidence that SN 2007od formed new dust before day 232. Late-time observations reveal a flattening of the visible lightcurve. This flattening is a strong indication of the presence of a light echo, which likely accounts for much of the broad, underlying Halpha component seen at late-times. We believe the multi-peaked Halpha emission is consistent with the interaction of the ejecta with a circumstellar ring or torus (for the inner components at \pm1500 km s-1), and a single blob or cloud of circumstellar material out of the plane of the CSM ring (for the outer component at -5000 km s-1). The most probable location for the formation of new dust is in the cool dense shell created by the interaction between the expanding ejecta and its CSM. Monte Carlo radiative transfer modeling of the dust emission from SN 2007od implies that up to 4x 10-4Msun of new dust has formed. This is similar to the amounts of dust formed in other CCSNe such as SNe 1999em, 2004et, and 2006jc.

Using Spitzer IRAC observations from the SAGE-SMC Legacy program and archived Spitzer IRAC data, we investigate dust production in 47 Tuc, a nearby massive Galactic globular cluster. A previous study detected infrared excess, indicative of circumstellar dust, in a large population of stars in 47 Tuc, spanning the entire Red Giant Branch (RGB). We show that those results suffered from effects caused by stellar blending and imaging artifacts and that it is likely that no stars below about 1 mag from the tip of the RGB are producing dust. The only stars that appear to harbor dust are variable stars, which are also the coolest and most luminous stars in the cluster.

We investigate the dust size and dust shell structure of the bipolar proto-planetary nebula M 1--92 by means of radiative transfer modeling. Our models consists of a disk and bipolar lobes that are surrounded by an AGB shell, each component having different dust characteristics. The upper limit of the grain size $a_\mathrm{max}$ in the lobes is estimated to be $0.5 \mu$m from the polarization value in the bipolar lobe. The $a_\mathrm{max}$ value of the disk is constrained with the disk mass (0.2 $M_{\sun}$), which was estimated from a previous CO emission line observation. We find a good model with $a_\mathrm{max}=1000.0 \mu$m, which provides an approximated disk mass of 0.15 $M_{\sun}$. Even taking into account uncertainties such as the gas-to-dust mass ratio, a significantly larger dust of $a_\mathrm{max}>100.0 \mu$m, comparing to the dust in the lobe, is expected. We also estimated the disk inner radius, the disk outer radius, and the envelope mass to be 30 $R_\star$(=9 AU), 4500 AU, and 4 $M_{\sun}$, respectively, where $v_\mathrm{exp}$ is the expansion velocity. If the dust existing in the lobes in large separations from the central star undergoes little dust processing, the dust sizes preserves the ones in the dust formation. Submicron-sized grains are found in many objects besides M 1--92, suggesting that the size does not depend much on the object properties, such as initial mass of the central star and chemical composition of the stellar system. On the other hand, the grain sizes in the disk do. Evidence of large grains has been reported in many bipolar PPNs, including M 1--92. This result suggests that disks play an important role in grain growth.

SN 2008S erupted in early 2008 in the grand design spiral galaxy NGC 6946. The progenitor was detected by Prieto et al. in Spitzer Space Telescope images taken over the four years prior to the explosion, but was not detected in deep optical images, from which they inferred a self-obscured object with a mass of about 10 Msun. We obtained Spitzer observations of SN 2008S five days after its discovery, as well as coordinated Gemini and Spitzer optical and infrared observations six months after its outburst. We have constructed radiative transfer dust models for the object before and after the outburst, using the same r^-2 density distribution of pre-existing amorphous carbon grains for all epochs and taking light-travel time effects into account for the early post-outburst epoch. We rule out silicate grains as a significant component of the dust around SN 2008S. The inner radius of the dust shell moved outwards from its pre-outburst value of 85 AU to a post-outburst value of 1250 AU, attributable to grain vaporisation by the light flash from SN 2008S. Although this caused the circumstellar extinction to decrease from Av = 15 before the outburst to 0.8 after the outburst, we estimate that less than 2% of the overall circumstellar dust mass was destroyed. The total mass-loss rate from the progenitor star is estimated to have been (0.5-1.0)x10^-4 Msun yr^-1. The derived dust mass-loss rate of 5x10^-7 Msun yr^-1 implies a total dust injection into the ISM of up to 0.01 Msun over the suggested duration of the self-obscured phase. We consider the potential contribution of objects like SN 2008S to the dust enrichment of galaxies.

We present spectroscopic observations of a sample of 15 embedded young stellar objects (YSOs) in the Large Magellanic Cloud (LMC). These observations were obtained with the Spitzer Infrared Spectrograph (IRS) as part of the SAGE-Spec Legacy program. We analyze the two prominent ice bands in the IRS spectral range: the bending mode of CO_2 ice at 15.2 micron and the ice band between 5 and 7 micron that includes contributions from the bending mode of water ice at 6 micron amongst other ice species. The 5-7 micron band is difficult to identify in our LMC sample due to the conspicuous presence of PAH emission superimposed onto the ice spectra. We identify water ice in the spectra of two sources; the spectrum of one of those sources also exhibits the 6.8 micron ice feature attributed to ammonium and methanol. We model the CO_2 band in detail, using the combination of laboratory ice profiles available in the literature. We find that a significant fraction (> 50%) of CO_2 ice is locked in a water-rich component, consistent with what is observed for Galactic sources. The majority of the sources in the LMC also require a pure-CO_2 contribution to the ice profile, evidence of thermal processing. There is a suggestion that CO_2 production might be enhanced in the LMC, but the size of the available sample precludes firmer conclusions. We place our results in the context of the star formation environment in the LMC.

Aims: By studying the photospheric abundances of 4 RV Tauri stars in the LMC, we test whether the depletion pattern of refractory elements, seen in similar Galactic sources, is also common for extragalactic sources. Since this depletion process probably only occurs through interaction with a stable disc, we investigate the circumstellar environment of these sources. Methods: A detailed photospheric abundance study was performed using high-resolution UVES optical spectra. To study the circumstellar environment we use photometric data to construct the spectral energy distributions of the stars, and determine the geometry of the circumstellar environment, whereas low-resolution Spitzer-IRS infrared spectra are used to trace its mineralogy. Results: Our results show that, also in the LMC, the photospheres of RV Tauri stars are commonly affected by the depletion process, although it can differ significantly in strength from source to source. From our detailed disc modelling and mineralogy study, we find that this process, as in the Galaxy, appears closely related to the presence of a stable Keplerian disc. The newly studied extragalactic objects have similar observational characteristics as Galactic post-AGB binaries surrounded by a dusty disc, and are therefore also believed to be part of a binary system. One source shows a very small infrared excess, atypical for a disc source, but still has evidence for depletion. We speculate this could point to the presence of a very evolved disc, similar to debris discs seen around young stellar objects.

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

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

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

We present a catalog of 1750 massive stars in the Large Magellanic Cloud, with accurate spectral types compiled from the literature, and a photometric catalog for a subset of 1268 of these stars, with the goal of exploring their infrared properties. The photometric catalog consists of stars with infrared counterparts in the Spitzer SAGE survey database, for which we present uniform photometry from 0.3-24 microns in the UBVIJHKs+IRAC+MIPS24 bands. The resulting infrared color-magnitude diagrams illustrate that the supergiant B[e], red supergiant and luminous blue variable (LBV) stars are among the brightest infrared point sources in the Large Magellanic Cloud, due to their intrinsic brightness, and at longer wavelengths, due to dust. We detect infrared excesses due to free-free emission among ~900 OB stars, which correlate with luminosity class. We confirm the presence of dust around 10 supergiant B[e] stars, finding the shape of their spectral energy distributions (SEDs) to be very similar, in contrast to the variety of SED shapes among the spectrally variable LBVs. The similar luminosities of B[e] supergiants (log L/Lo>=4) and the rare, dusty progenitors of the new class of optical transients (e.g. SN 2008S and NGC 300 OT), plus the fact that dust is present in both types of objects, suggests a common origin for them. We find the infrared colors for Wolf-Rayet stars to be independent of spectral type and their SEDs to be flatter than what models predict. The results of this study provide the first comprehensive roadmap for interpreting luminous, massive, resolved stellar populations in nearby galaxies at infrared wavelengths.

We report near-infrared photometric measurements of 35 Herbig Ae/Be candidate stars obtained with direct imaging and aperture photometry. Observations were made through the broadband J, H, and K' filters, with each source imaged in at least one of the wavebands. We achieved subarcsecond angular resolution for all observations, providing us with the opportunity to search for close binary candidates and extended structure. The imaging revealed five newly identified binary candidates and one previously resolved T Tauri binary among the target sources with separations of <~2.5". Separate photometry is provided for each of the binary candidate stars. We detect one extended source that has been identified as a protoplanetary nebula. Comparing our magnitudes to past measurements yields significant differences for some sources, possibly indicating photometric variability. H-band finding charts for all of our sources are provided to aid follow-up high-resolution imaging.

An updated Science Vision for the SOFIA project is presented, including an overview of the characteristics and capabilities of the observatory and first generation instruments. A primary focus is placed on four science themes: 'The Formation of Stars and Planets', 'The Interstellar Medium of the Milky Way', 'Galaxies and the Galactic Center' and 'Planetary Science'.

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 report on an analysis of the gas and dust budget in the the interstellar medium (ISM) of the Large Magellanic Cloud (LMC). Recent observations from the Spitzer Space Telescope enable us to study the mid-infrared dust excess of asymptotic giant branch (AGB) stars in the LMC. This is the first time we can quantitatively assess the gas and dust input from AGB stars over a complete galaxy, fully based on observations. The integrated mass-loss rate over all intermediate and high mass-loss rate carbon-rich AGB candidates in the LMC is 8.5x10^-3 solar mass per year, up to 2.1x10^-2 solar mass per year. This number could be increased up to 2.7x10^-2 solar mass per year, if oxygen-rich stars are included. This is overall consistent with theoretical expectations, considering the star formation rate when these low- and intermediate-mass stars where formed, and the initial mass functions. AGB stars are one of the most important gas sources in the LMC, with supernovae (SNe), which produces about 2-4x10^-2 solar mass per year. At the moment, the star formation rate exceeds the gas feedback from AGB stars and SNe in the LMC, and the current star formation depends on gas already present in the ISM. This suggests that as the gas in the ISM is exhausted, the star formation rate will eventually decline in the LMC, unless gas is supplied externally. Our estimates suggest `a missing dust-mass problem' in the LMC, which is similarly found in high-z galaxies: the accumulated dust mass from AGB stars and possibly SNe over the dust life time (400--800 Myrs) is significant less than the dust mass in the ISM. Another dust source is required, possibly related to star-forming regions.

We present empirical relations describing excess emission from evolved stars in the Large Magellanic Cloud (LMC) using data from the SAGE (Surveying the Agents of a Galaxy's Evolution) survey which includes the IRAC 3.6, 4.5, 5.8 and 8.0 \mu m and MIPS 24, 70 and 160 \mu m bands. We combine the SAGE data with the Two Micron All Sky Survey (2MASS; J, H and Ks) and the optical Magellanic Cloud Photometric Survey (MCPS; U, B, V and I) point source catalogs to create complete spectral energy distributions (SEDs) of the asymptotic giant branch (AGB) star candidates in the LMC. AGB star outflows are among the main producers of dust in a galaxy, and this mass loss results in an excess in the fluxes observed in the 8 and 24 \mic m bands. We identify oxygen-rich, carbon-rich and extreme AGB star populations in our sample based on their 2MASS and IRAC colors. We calculate excesses fluxes in the mid-IR bands by comparison of the SEDs with model photospheres. We find about 16,000 O-rich, 6300 C-rich and 1000 extreme sources with reliable 8 \mu m excesses, and about 4500 O-rich, 5300 C-rich and 960 extreme sources with reliable 24 \mic m excesses. The excesses are in the range 0.1 mJy--5 Jy. The 8 and 24 \mic m excesses for all three types of AGB candidates show a general increasing trend with luminosity. The dust color temperature derived from the ratio of the 8 and 24 \mic m excesses decreases with an increase in excess, while the 24 \mic m optical depth increases with excess. The extreme AGB candidates are the major contributors to the mass loss, and we estimate the total AGB mass-loss return to the LMC to be (5.9--13)\times 10^{-3} Msun yr^{-1}.

We present a mid-infrared spectroscopic data cube of the central part of 30 Doradus, observed with Spitzer's IRS and MIPS/SED mode. Aromatic dust emission features and emission lines from molecular and atomic hydrogen are detected but not particularly strong. The dominant spectral features are emission lines from moderately ionized species of argon, neon, and sulphur, which are used to determine the physical conditions in the ionized gas. The ionized gas excitation shows strong variations on parsec scales, some of which can plausibly be associated with individual hot stars. We fit the ionic line strengths with photoionization and shock models, and find that photoionization dominates in the region. The ionization parameter U traces the rim of the central bubble, as well as highlighting isolated sources of ionization, and at least one quiescent clump. The hardness of the ionizing radiation field T_rad reveals several "hot spots" that are either the result of individual very hot stars or trace the propagation of the diffuse ionizing field through the surrounding neutral cloud. Consistent with other measurements of giant molecular hydrogen regions, log(U) ranges between -3 and -0.75, and T_rad between 30000 and 85000K.

In late 2003, the young eruptive variable star V1647 Orionis optically brightened by over 5 magnitudes, stayed bright for around 26 months, and then decline to its pre-outburst level. In August 2008 the star was reported to have unexpectedly brightened yet again and we herein present the first detailed observations of this new outburst. Photometrically, the star is now as bright as it ever was following the 2003 eruption. Spectroscopically, a pronounced P Cygni profile is again seen in Halpha with an absorption trough extending to -700 km/s. In the near-infrared, the spectrum now possesses very weak CO overtone bandhead absorption in contrast to the strong bandhead emission seen soon after the 2003 event. Water vapor absorption is also much stronger than previously seen. We discuss the current outburst below and relate it to the earlier event.

Deep HST/ACS photometry of the young cluster NGC 602, located in the remote low density "wing" of the Small Magellanic Cloud, reveals numerous pre-main sequence stars as well as young stars on the main sequence. The resolved stellar content thus provides a basis for studying the star formation history into recent times and constraining several stellar population properties, such as the present day mass function, the initial mass function and the binary fraction. To better characterize the pre-main sequence population, we present a new set of model stellar evolutionary tracks for this evolutionary phase with metallicity appropriate for the Small Magellanic Cloud (Z = 0.004). We use a stellar population synthesis code, which takes into account a full range of stellar evolution phases to derive our best estimate for the star formation history in the region by comparing observed and synthetic color-magnitude diagrams. The derived present day mass function for NGC 602 is consistent with that resulting from the synthetic diagrams. The star formation rate in the region has increased with time on a scale of tens of Myr, reaching $0.3-0.7 \times 10^{-3} M_\odot yr^{-1}$ in the last 2.5 Myr, comparable to what is found in Galactic OB associations. Star formation is most complete in the main cluster but continues at moderate levels in the gas-rich periphery of the nebula.

(abridged) We study photoelectric heating throughout the Large Magellanic Cloud. We quantify the importance of the [CII] cooling line and the photoelectric heating process of various environments in the LMC and investigate which parameters control the extent of photoelectric heating. We use the BICE [CII] map and the Spitzer/SAGE infrared maps. We examine the spatial variations in the efficiency of photoelectric heating: photoelectric heating rate over power absorbed by grains. We correlate the photoelectric heating efficiency and the emission from various dust constituents and study the variations as a function of H\alpha emission, dust temperatures, and the total infrared luminosity. From this we estimate radiation field, gas temperature, and electron density. We find systematic variations in photoelectric efficiency. The highest efficiencies are found in the diffuse medium, while the lowest coincide with bright star-forming regions (~1.4 times lower). The [CII] line emission constitutes 1.32% of the far infrared luminosity across the whole of the LMC. We find correlations between the [CII] emission and ratios of the mid infrared and far infrared bands, which comprise various dust constituents. The correlations are interpreted in light of the spatial variations of the dust abundance and by the local environmental conditions that affect the dust emission properties. As a function of the total infrared surface brightness, S_{TIR}, the [CII] surface brightness can be described as: S_{[CII]}=1.25 S_{TIR}^{0.69} [10^{-3} erg s^{-1} cm^{-2} sr^{-1}]. The [CII] emission is well-correlation with the 8 micrometer emission, suggesting that the polycyclic aromatic hydrocarbons play a dominant role in the photoelectric heating process.

The Tail region of the Small Magellanic Cloud (SMC) was imaged using the MIPS instrument on the Spitzer Space Telescope as part of the SAGE-SMC Spitzer Legacy. Diffuse infrared emission from dust was detected in all the MIPS bands. The Tail gas-to-dust ratio was measured to be 1200 +/- 350 using the MIPS observations combined with existing IRAS and HI observations. This gas-to-dust ratio is higher than the expected 500-800 from the known Tail metallicity indicating possible destruction of dust grains. Two cluster regions in the Tail were resolved into multiple sources in the MIPS observations and local gas-to-dust ratios were measured to be ~440 and ~250 suggests dust formation and/or significant amounts of ionized gas in these regions. These results support the interpretation that the SMC Tail is a tidal tail recently stripped from the SMC that includes gas, dust, and young stars.

We present initial results and source lists of variable sources in the Large Magellanic Cloud (LMC) for which we detect thermal infrared variability from the SAGE (Surveying the Agents of a Galaxy's Evolution) survey, which had 2 epochs of photometry separated by three months. The SAGE survey mapped a 7 degree by 7 degree region of the LMC using the IRAC and the MIPS instruments on board Spitzer. Variable sources are identified using a combination of the IRAC 3.6, 4.5, 5.8, 8.0 \micron bands and the MIPS 24 \micron bands. An error-weighted flux difference between the two epochs is used to assess the variability. Of the ~ 3 million sources detected at both epochs we find ~ 2,000 variable sources for which we provide electronic catalogs. Most of the variable sources can be classified as asymptotic giant branch (AGB) stars. A large fraction (> 66%) of the extreme AGB stars are variable and only smaller fractions of carbon-rich (6.1%) and oxygen-rich (2.0%) stars are detected as variable. We also detect a population of variable young stellar object candidates.

We introduce the SMC in space and time, a large coordinated space and ground-based program to study star formation processes and history, as well as variable stars, structure, kinematics and chemical evolution of the whole SMC. Here, we present the Colour-Magnitude Diagrams(CMDs) resulting from HST/ACS photometry, aimed at deriving the star formation history (SFH) in six fields of the SMC. The fields are located in the central regions, in the stellar halo, and in the wing toward the LMC. The CMDs are very deep, well beyond the oldest Main Sequence Turn-Off, and will allow us to derive the SFH over the entire Hubble time.

The results from a Hubble Space Telescope (HST) snapshot survey of post-AGB objects are shown. The aim of the survey is to complement existing HST images of PPN and to connect various types of nebulosities with physical and chemical properties of their central stars. Nebulosities are detected in 15 of 33 sources. Images and photometric and geometric measurements are presented. For sources with nebulosities we see a morphological bifurcation into two groups, DUPLEX and SOLE, as previous studies have found. We find further support to the previous results suggesting that this dichotomy is caused by a difference in optical thickness of the dust shell. The remaining 18 sources are classified as stellar post-AGB objects, because our observations indicate a lack of nebulosity. We show that some stellar sources may in fact be DUPLEX or SOLE based on their infrared colors. The cause of the differences among the groups are investigated. We discuss some evidence suggesting that high progenitor-mass AGB stars tend to become DUPLEX post-AGB objects. Intermediate progenitor-mass AGB stars tend to be SOLE post-AGB objects. Most of the stellar sources probably have low mass progenitors and do not seem to develop nebulosities during the post-AGB phase and therefore do not become planetary nebulae.

We present IRAC and MIPS images and photometry of a sample of previously known planetary nebulae (PNe) from the SAGE survey of the Large Magellanic Cloud (LMC) performed with the Spitzer Space Telescope. Of the 233 known PNe in the survey field, 185 objects were detected in at least two of the IRAC bands, and 161 detected in the MIPS 24 micron images. Color-color and color-magnitude diagrams are presented using several combinations of IRAC, MIPS, and 2MASS magnitudes. The location of an individual PN in the color-color diagrams is seen to depend on the relative contributions of the spectral components which include molecular hydrogen, polycyclic aromatic hydrocarbons (PAHs), infrared forbidden line emission from the ionized gas, warm dust continuum, and emission directly from the central star. The sample of LMC PNe is compared to a number of Galactic PNe and found to not significantly differ in their position in color-color space. We also explore the potential value of IR PNe luminosity functions (LFs) in the LMC. IRAC LFs appear to follow the same functional form as the well-established [O III] LFs although there are several PNe with observed IR magnitudes brighter than the cut-offs in these LFs.

Deep F555W and F814W Hubble Space Telescope ACS images are the basis for a study of the present day mass function (PDMF) of NGC346, the largest active star forming region in the Small Magellanic Cloud (SMC). We find a PDMF slope of Gamma=-1.43+/-0.18 in the mass range 0.8-60 Mo, in excellent agreement with the Salpeter Initial Mass Function (IMF) in the solar neighborhood. Caveats on the conversion of the PDMF to the IMF are discussed. The PDMF slope changes, as a function of the radial distance from the center of the NGC 346 star cluster, indicating a segregation of the most massive stars. This segregation is likely primordial considering the young age (~3 Myr) of NGC346, and its clumpy structure which suggests that the cluster has likely not had sufficient time to relax. Comparing our results for NGC346 with those derived for other star clusters in the SMC and the Milky Way (MW), we conclude that, while the star formation process might depend on the local cloud conditions, the IMF does not seem to be affected by general environmental effects such as galaxy type, metallicity, and dust content.

NGC 602 is a young stellar cluster located in a peripheral region of the Small Magellanic Cloud known as the wing. Far from the main body of the galaxy and abutting the Magellanic Bridge, the SMC's wing is characterized by low gas and stellar content. With deep optical imaging from the Advanced Camera for Surveys (ACS) aboard the Hubble Space Telescope (HST), we have discovered an extensive pre-Main Sequence (PMS) population, with stellar masses in the range 0.6-3 Solar masses. These low mass PMS stars formed coevally with the central cluster about 4 Myr ago. Spitzer Space Telescope (Spitzer) images of the same region from the Infrared Array Camera (IRAC) also reveal a population of Young Stellar Objects (YSOs), some of which are still embedded in nebular material and most of which likely formed even more recently than the young stars detected with HST/ACS imaging. We infer that star formation started in this region about 4 Myr ago with the formation of the central cluster and gradually propagated towards the outskirts where star formation is presently ongoing.

We present new, very late-time optical photometry and spectroscopy of the interesting Type II-P supernova, SN 2002hh, in NGC 6946. Gemini/GMOS-N has been used to acquire visible spectra at six epochs between 2004 August and 2006 July, following the evolution of the SN from age 661 to 1358 days. Few optical spectra of Type II supernovae with ages greater than one year exist. In addition, g'r'i' images were acquired at all six epochs. The spectral and photometric evolution of SN 2002hh has been very unusual. Measures of the brightness of this SN, both in the R and I bands as well as in the H-alpha emission flux, show no significant fading over an interval of nearly two years. The most straightforward explanation for this behavior is that the light being measured comes not only from the SN itself but also from an echo off of nearby dust. Echoes have been detected previously around several SNe but these echoes, at their brightest, were ~8 mag below the maximum brightness of the SN. At V~21 mag, the putative echo dominates the light of SN 2002hh and is only ~4 mag below the outburst's peak brightness. There is an estimated 6 magnitudes of total extinction in V towards SN 2002hh. The proposed explanation of a differential echo/SN absorption is inconsistent with the observed BVRI colors.

The structure of the dusty circumstellar envelopes (CSEs) of proto-planetary nebulae (PPNs) reveals the mass-loss history of these sources and how such histories may differ for elliptical (SOLE) and bipolar (DUPLEX) PPNs. To study the PPN structures via dust-scattered linearly polarized starlight, we have compiled the imaging-polarimetric data for all 18 evolved stars that have been obtained to date with NICMOS on-board the Hubble Space Telescope (HST). This alternative imaging technique provides a unique way to probe the distribution of dust grains that scatter light around evolved stars. The new perspective gained from the imaging-polarimetric data has revealed several new aspects to the structures of PPNs. Point-symmetry is a prevalent imaging-polarimetric characteristic resulting from the azimuthal density gradient in the CSEs. Among these point-symmetric nebulae, three detailed morphological types can be differentiated by their polarized intensity, I_pol, and polarization strength, P. While the azimuthal density gradient is reversed above and below the equatorial plane in optically thicker bipolar nebulae, there is no gradient reversal in optically thinner elliptical nebulae. The equatorial plane of the system defined by the integrated angle of polarization is not necessarily orthogonal to the axis of the apparent bipolar structure in the total intensity data.

We present the results of a coronagraphic scattered-light imaging survey of six young disk candidate stars using the Hubble Space Telescope Advanced Camera for Surveys. The observations made use of the 1.8" occulting spot through the F606W (broad V) filter. Circumstellar material was imaged around HD 97048, a Herbig Ae/Be star located in the Chamaeleon I dark cloud at a distance of 180 pc. The material is seen between ~2" (360 AU) and ~4" (720 AU) from the star in all directions. A V-band azimuthally-averaged radial surface brightness profile peaks at r = 2" with a value of 19.6 +/- 0.2 mag arcsec^-2 and smoothly decreases with projected distance from the star as r^(-3.3 +/- 0.5). An integrated flux of 16.8 +/- 0.1 mag is measured between 2" and 4", corresponding to a scattered-light fractional luminosity lower limit of 8.4 x 10^-4. Filamentary structure resembling spiral arms similar to that seen in Herbig Ae/Be disks is observed. Such structure has been attributed to the influence of orbiting planets or stellar encounters. Average surface brightness upper limits are determined for the five non-detections: HD 34282, HD 139450, HD 158643, HD 159492, and HD 195627. Possible reasons for the non-detections are disks that are too faint or disks hidden by the occulter.

In the quest of understanding how star formation occurs and propagates in the low metallicity environment of the Small Magellanic Cloud (SMC), we acquired deep F555W (~V), and F814W (~I) HST/ACS images of the young and massive star forming region NGC 346. These images and their photometric analysis provide us with a snapshot of the star formation history of the region. We find evidence for star formation extending from ~10 Gyr in the past until ~150 Myr in the field of the SMC. The youngest stellar population (~3 +/- 1 Myr) is associated with the NGC 346 cluster. It includes a rich component of low mass pre-main sequence stars mainly concentrated in a number of sub-clusters, spatially co- located with CO clumps previously detected by Rubio et al. (2000). Within our analysis uncertainties, these sub-clusters appear coeval with each other. The most massive stars appear concentrated in the central sub-clusters, indicating possible mass segregation. A number of embedded clusters are also observed. This finding, combined with the overall wealth of dust and gas, could imply that star formation is still active. An intermediate age star cluster, BS90, formed ~4.3 +/-0.1 Gyr ago, is also present in the region. Thus, this region of the SMC has supported star formation with varying levels of intensity over much of the cosmic time.

Color-magnitude diagrams (CMDs) are presented for the Spitzer SAGE (Surveying the Agents of a Galaxy's Evolution) survey of the Large Magellanic Cloud (LMC). IRAC and MIPS 24 um epoch one data are presented. These data represent the deepest, widest mid-infrared CMDs of their kind ever produced in the LMC. Combined with the 2MASS survey, the diagrams are used to delineate the evolved stellar populations in the Large Magellanic Cloud as well as Galactic foreground and extragalactic background populations. Some 32000 evolved stars brighter than the tip of the red giant branch are identified. Of these, approximately 17500 are classified as oxygen-rich, 7000 carbon-rich, and another 1200 as ``extreme'' asymptotic giant branch (AGB) stars. Brighter members of the latter group have been called ``obscured'' AGB stars in the literature owing to their dusty circumstellar envelopes. A large number (1200) of luminous oxygen--rich AGB stars/M supergiants are also identified. Finally, there is strong evidence from the 24 um MIPS channel that previously unexplored, lower luminosity oxygen-rich AGB stars contribute significantly to the mass loss budget of the LMC (1200 such sources are identified).

New Spitzer imaging observations have revealed the structure around the Mira variable star R Hya to be a one-sided parabolic arc 100 arcsec to the West stretching from North to South. We successfully model R Hya and its surroundings in terms of an interaction of the stellar wind from an asymptotic giant branch (AGB) star with the interstellar medium (ISM) the star moves through. Our three-dimensional hydrodynamic simulation reproduces the structure as a bow shock into the oncoming ISM. We propose this as another explanation of detached shells around such stars which should be considered alongside current theories of internal origin. The simulation predicts the existence of a tail of ram-pressure-stripped AGB material stretching downstream. Indications for such a tail behind R Hya are seen in IRAS maps.

We present the first results of the MIRIAD (MIPS [Multiband Imaging Photometer for Spitzer] Infra-Red Imaging of AGB [asymptotic giant branch] Dustshells) project using the Spitzer Space Telescope. The primary aim of the project is to probe the material distribution in the extended circumstellar envelopes (CSE) of evolved stars and recover the fossil record of their mass loss history. Hence, we must map the whole of the CSEs plus the surrounding sky for background subtraction, while avoiding the central star that is brighter than the detector saturation limit. With our unique mapping strategy, we have achieved better than one MJy/sr sensitivity in three hours of integration and successfully detected a faint (< 5 MJy/sr), extended (~400 arcsec) far-infrared nebula around the AGB star R Hya. Based on the parabolic structure of the nebula, the direction of the space motion of the star with respect to the nebula shape, and the presence of extended H alpha emission co-spatial to the nebula, we suggest that the detected far-IR nebula is due to a bow shock at the interface of the interstellar medium and the AGB wind of this moving star. This is the first detection of the stellar-wind bow-shock interaction for an AGB star and exemplifies the potential of Spitzer as a tool to examine the detailed structure of extended far-IR nebulae around bright central sources. \

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

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

With central densities way above the density of atomic nuclei, neutron stars contain matter in one of the densest forms found in the universe. Depending of the density reached in the cores of neutron stars, they may contain stable phases of exotic matter found nowhere else in space. This article gives a brief overview of the phases of ultra-dense matter predicted to exist deep inside neutron stars and discusses the equation of state associated with such matter.

We report the discovery of a rich population of low mass stars in the young, massive star forming region N66/NGC346 in the Small Magellanic Cloud, from deep V, I and H alpha images taken with the HST/ACS. These stars have likely formed together with the NGC346 cluster, ~3-5 Myr ago. Their magnitude and colors are those of pre-main sequence stars in the mass range 0.6-3 Mo, mostly concentrated in the main cluster, but with secondary subclusters spread over a region across ~45 pc. These subclusters appear to be spatially coincident with previously known knots of molecular gas identified in ground based and ISO observations. We show that N66/NGC346 is a complex region, being shaped by its massive stars, and the observations presented here represent a key step towards the understanding of how star formation occurred and has progressed in this low metallicity environment.

We present the results of proper-motion measurements of the dust shell structure in the Egg Nebula (AFGL 2688, CRL 2688, V1610 Cyg), based on the archived two-epoch imaging-polarimetry data in the 2 micron band taken with NICMOS on-board the Hubble Space Telescope. We measured the amount of motions of local structures in the nebula by determining their relative shifts over an interval of 5.5 years. The dynamical age of the nebula is found to be roughly 350 years based on the overall motion of the nebula that exhibits a Hubble-law-esque linear relation between the measured proper motion and the projected radial distance from the origin of the expansion. By adopting the de-projected velocity of 45 km/s at the tips of the bipolar lobes, our proper-motion measurements indicate that the distance to the Egg Nebula is about 420 pc and that the lobes are inclined at 7.7 degrees with respect to the plane of the sky. The refined distance estimate yields the luminosity of the central star of 3.3 x 10^3 L_sun, the total shell mass of 1.2 M_sun, and the mass loss rate (the upper limit) of 3.6 x 10^(-3) M_sun/yr. Assuming 0.6 M_sun central post-AGB stellar mass, the initial mass of the Egg is 1.8 M_sun. Upon analysis, we also discovered that the central star of the Egg Nebula has proper motion of its own at the rate of (14, -10) mas/yr.

We present HST/NICMOS imaging of the H_2 2.12 \mu m emission in 5 fields in the Helix Nebula ranging in radial distance from 250-450" from the central star. The images reveal arcuate structures with their apexes pointing towards the central star. Comparison of these images with comparable resolution ground based images reveals that the molecular gas is more highly clumped than the ionized gas line tracers. From our images, we determine an average number density of knots in the molecular gas ranging from 162 knots/arcmin^2 in the denser regions to 18 knots/arcmin^2 in the lower density outer regions. Using this new number density, we estimate that the total number of knots in the Helix to be ~23,000 which is a factor of 6.5 larger than previous estimates. The total neutral gas mass in the Helix is 0.35 M_\odot assuming a mass of \~1.5x10^{-5} M_\odot for the individual knots. The H_2 intensity, 5-9x10^{-5} erg s^{-1} cm^{-2} sr^{-1}, remains relatively constant with projected distance from the central star suggesting a heating mechanism for the molecular gas that is distributed almost uniformly in the knots throughout the nebula. The temperature and H_2 2.12 \mu m intensity of the knots can be approximately explained by photodissociation regions (PDRs) in the individual knots; however, theoretical PDR models of PN under-predict the intensities of some knots by a factor of 10.

Dust emission from the Type II supernova SN 2002hh in NGC 6946 has been detected at mid-infrared wavelengths by the Spitzer Space Telescope from 590 to 758 days after outburst and confirmed by higher angular resolution Gemini-N mid-IR observations. The day-600 5.8-24 micron emission can be fit by a 290-K blackbody having a luminosity of 1.6E+7 L_Sol. The minimum emitting radius of 1.1E+17 cm is too large for the emitting dust to have been formed in the supernova ejecta. Using radiative transfer models and realistic dust grain parameters, fits to the observed flux distribution could be obtained with an optically thick dust shell having a mass of 0.10-0.15 M_Sol, corresponding to a total dust+gas mass in excess of 10 M_Sol, suggesting a massive M supergiant or luminous blue variable precursor to this self-obscured object.

The Eskimo Nebula (NGC 2392) is a double-shell planetary nebula (PN) known for the exceptionally large expansion velocity of its inner shell, ~90 km/s, and the existence of a fast bipolar outflow with a line-of-sight expansion velocity approaching 200 km/s. We have obtained XMM-Newton observations of the Eskimo and detected diffuse X-ray emission within its inner shell. The X-ray spectra suggest thin plasma emission with a temperature of ~2x10^6 K and an X-ray luminosity of L_X = (2.6+/-1.0)x10^31 (d/1150 pc)^2 ergs/s, where d is the distance in parsecs. The diffuse X-ray emission shows noticeably different spatial distributions between the 0.2-0.65 keV and 0.65-2.0 keV bands. High-resolution X-ray images of the Eskimo are needed to determine whether its diffuse X-ray emission originates from shocked fast wind or bipolar outflows.

Using NICMOS on HST, we have performed imaging polarimetry of proto-planetary nebulae. Our objective is to study the structure of optically thin circumstellar shells of post-asymptotic giant branch stars by separating dust-scattered, linearly polarized star light from unpolarized direct star light. This unique technique allows us to probe faint reflection nebulae around the bright central star, which can be buried under the point-spread-function of the central star in conventional imaging. Our observations and archival search have yielded polarimetric images for five sources: IRAS 07134+1005 (HD 56126), IRAS 06530-0213, IRAS 04296+3429, IRAS (Z)02229+6208, and IRAS 16594-4656. These images have revealed the circumstellar dust distribution in an unprecedented detail via polarized intensity maps, providing a basis to understand the 3-D structure of these dust shells. We have observationally confirmed the presence of the inner cavity caused by the cessation of AGB mass loss and the internal shell structures which is strongly tied to the progenitor star's mass loss history on the AGB. We have also found that equatorial enhancement in these circumstellar shells comes with various degrees of contrast, suggesting a range of optical depths in these optically thin shells. Our data support the interpretation that the dichotomy of PPN morphologies is due primarily to differences in optical depth and secondary to the inclination effect. The polarization maps reveal a range of inclination angles for these optically thin reflection nebulae, dispelling the notion that elliptical nebulae are pole-on bipolar nebulae.

We present millimeter interferometry images of the CO J=1-0 line emission arising in the circumstellar envelope of HD 56126 (a.k.a. IRAS 07134+1005), which is one of the best studied 21-micron proto-planetary nebulae (PPNs). The CO emission extends from 1.2" to 7" in radius from the central star and appears consistent with a simple expanding envelope, as expected for a post-AGB star. We quantitatively model the molecular envelope using a radiative transfer code that we have modified for detached shells. Our best fit model reveals that two sequential winds created the circumstellar envelope of HD 56126: an AGB wind that lasted 6500 years with a mass-loss rate of 5.1x10^{-6} M_odot yr^{-1} and a more intense superwind that lasted 840 years with a mass-loss rate of 3x10^{-5} M_odot yr^{-1} and that ended the star's life on the AGB 1240 years ago. The total mass of this envelope is 0.059 M_odot which indicates a lower limit progenitor mass for the system of 0.66 M_odot, quite reasonable for this low-metallicity star which probably resides in the thick disk of the Galaxy. Comparison with images of the dust emission reveal a similar structure with the gas in the inner regions. Using 2-DUST, we model the dust emission of this source so that the model is consistent with the CO emission model and find an average gas-to-dust mass ratio of 75.

Through HST imaging of the inner part of the main-ring of the Helix Nebula together with CTIO 4-m images of the fainter outer parts, we have an unprecedented-quality view of the nearest bright planetary nebula. These images have allowed determination that the main-ring of the nebula is composed of an inner-disk of about 499\arcsec diameter (0.52 pc) surrounded by an outer-ring (in reality a torus) of 742\arcsec diameter (0.77 pc) whose plane is highly inclined to the plane of the disk. This outer-ring is surrounded by an outermost-ring of 1500\arcsec (1.76 pc) diameter which is flattened on the side colliding with the ambient interstellar medium. The inner-disk has an extended distribution of low density gas along its rotational axis of symmetry and the disk is optically thick to ionizing radiation, as is the outer-ring. Published radial velocities of the knots provides support for the two-component structure of the main-ring of the nebula and to the idea that the knots found there are expanding along with the nebular material from which it recently originated. There is a change in the morphology of the knots as a function of the distance from the local ionization front. This supports a scenario in which the knots are formed in or near the ionization front and are then sculpted by the stellar radiation from the central star as the ionization front advances beyond them.

We present the preliminary results from our HST/NICMOS imaging polarimetry on optically thin PPNs. The data show the dust distribution in an unprecedented detail. The structure of these faint PPNs is revealed via polarization properties of the dust-scattered light without much interference from the bright central star, proving the unique advantages of the technique over simple imaging.

We present JHK' photometric measurements of 78 objects mostly consisting of proto-planetary nebula candidates. Photometric magnitudes are determined by means of imaging and aperture photometry. Unlike the observations with a photometer with a fixed-sized beam, the method of imaging photometry permits accurate derivation of photometric values because the target sources can be correctly identified and confusion with neighboring sources can be easily avoided. Of the 78 sources observed, we report 10 cases in which the source seems to have been misidentified or confused by nearby bright sources. We also present nearly two dozen cases in which the source seems to have indicated a variability which prompts a follow-up monitoring. There are also a few sources that show previously unreported extendedness. In addition, we present H band finding charts of the target sources.

We have developed a general purpose dust radiative transfer code for an axisymmetric system, 2-Dust, motivated by the recent increasing availability of high-resolution images of circumstellar dust shells at various wavelengths. This code solves the equation of radiative transfer following the principle of long characteristic in a 2-D polar grid while considering a 3-D radiation field at each grid point. A solution is sought through an iterative scheme in which self-consistency of the solution is achieved by requiring a global luminosity constancy throughout the shell. The dust opacities are calculated through Mie theory from the given size distribution and optical properties of the dust grains. The main focus of the code is to obtain insights on (1) the global energetics of dust grains in the shell (2) the 2-D projected morphologies that are strongly dependent on the mixed effects of the axisymmetric dust distribution and inclination angle of the shell. Here, test models are presented with discussion of the results. The code can be supplied with a user-defined density distribution function, and thus, is applicable to a variety of dusty astronomical objects possessing the axisymmetric geometry.

We use detailed radiative transfer models to investigate the differences between the Star-Obvious Low-level-Elongated proto-planetary nebulae (SOLE PPNs) and DUst-Prominent Longitudinally-EXtended proto-planetary nebulae (DUPLEX PPNs) which are two subclasses of PPNs suggested by Ueta, Meixner, & Bobrowsky (2000). We select one SOLE PPN, HD 161796, and one DUPLEX PPN, IRAS 17150$-$3224, both of which are well studied and representative of their PPN classes. Using an axisymmetric dust shell radiative transfer code, we model these two sources in detail and constrain their mass-loss histories, inclination angles and dust composition. The physical parameters derived for HD 161796 and IRAS 17150$-$3224 demonstrate that they are physically quite different and that their observed differences cannot be attributed to inclination angle effects. If these calculations reflect a more general truth about SOLE vs. DUPLEX PPNs, then these two subclasses of PPNs are physically distinct with the SOLE PPNs derived from low mass progenitors and DUPLEX PPNs derived from high mass progenitors.

We present high-sensitivity near-IR images of a carbon-rich proto-planetary nebula, AFGL 618, obtained with the Subaru Telescope. These images have revealed ``bullets'' and ``horns'' extending farther out from the edges of the previously known bipolar lobes. The spatial coincidence between these near-IR micro-structures and the optical collimated outflow structure, together with the detection of shock-excited, forbidden IR lines of atomic species, strongly suggests that these bullets and horns represent the locations from which [\ion{Fe}{2}] IR lines arise. We have also discovered CO clumps moving at $> 200$ km s$^{-1}$ at the positions of the near-IR bullets by re-analyzing the existing $^{12}$CO $J=1-0$ interferometry data. These findings indicate that the near-IR micro-structures represent the positions of shocked surfaces at which fast-moving molecular clumps interface with the ambient circumstellar shell.

We report sub-arcsecond imaging of extended mid-infrared emission from a proto-planetary nebula (PPN), \iras 22272+5435, performed at the MMT observatory with its newly upgraded 6.5 m aperture telescope and at the Keck observatory. The mid-infrared emission structure is resolved into two emission peaks separated by $0\arcsec.5 - 0\arcsec.6$ in the MMT 11.7 $\um$ image and in the Keck 7.9, 9.7, and 12.5 $\um$ images, corroborating the predictions based on previous multi-wavelength morphological studies and radiative transfer calculations. The resolved images show that the PPN dust shell has a toroidal structure with the $0\arcsec.5$ inner radius. In addition, an unresolved mid-IR excess appears at the nebula center. Radiative transfer model calculations suggest that the highly equatorially-enhanced ($\rho_{\rm eq}/\rho_{\rm pole} = 9$) structure of the PPN shell was generated by an axisymmetric superwind with ${\dot M}_{\rm sw} = 4 \times 10^{-6} M_{\odot}$ yr$^{-1}$, which was preceded by a spherical asymptotic giant branch (AGB) wind with ${\dot M}_{\rm AGB} = 8 \times 10^{-7} M_{\odot}$ yr$^{-1}$. These model calculations also indicate that the superwind shell contains larger dust grains than the AGB wind shell. The unresolved mid-infrared excess may have been produced by a post-AGB mass loss at a rate of $2 - 6 \times 10^{-7} M_{\odot}$ yr$^{-1}$. While the central star left the AGB about 380 years ago after the termination of the superwind, the star seems to have been experiencing an ambient post-AGB mass loss with a sudden, increased mass ejection about 10 years ago.

We will review our mid-infrared and HST imaging surveys of the circumstellar dust shells of proto-planetary nebulae. While optical imaging indirectly probes the dust distribution via dust-scattered starlight, mid-IR imaging directly maps the distribution of warm dust grains. Both imaging surveys revealed preferencially axisymmetric nature of PPN dust shells, suggesting that axisymmetry in planetary nebulae sets in by the end of the asymptotic giant branch phase, most likely by axisymmetric superwind mass loss. Moreover, both surveys yielded two morphological classes which have one-to-one correspondence between the two surveys, indicating that the optical depth of circumstellar dust shells plays an equally important role as the inclination angle in determining the morphology of the PPN shells.

We investigate spatial distributions of specific dust components in the circumstellar shells of a proto-planetary nebula candidate, HD 179821, and a planetary nebula, BD$+30^{\circ}$3639, by means of spectral imaging. With high-resolution ground-based images and ISO spectra in the mid-infrared, we can derive ``dust feature only'' maps by subtracting synthesized continuum maps from the observed images at the feature wavelength. Such spatially detailed information will help to develop models for the evolution of dust grains around evolved stars.

We report detection of an extended mid-infrared emission from \iras 18576+0341 (AFGL 2298). The object shows a dusty circumstellar shell that has diameter of $\age 7\arcsec$ at 10.3 and 18.0 $\um$. The dust nebula shows two emission peaks concentrically elongated and symmetrically oriented on the opposite sides of the third, central peak, which appears to be the central star of the system. The observed mid-infrared morphology indicates that the circumstellar dust shell has an equatorially-enhanced material distribution, which is a common signature of stellar objects that have experienced mass loss. Radiative transfer model calculations suggest that the central star is an extremely bright ($\lstar = 10^{6.4} \lsun$) star at a distance of about 10 kpc: this object is best described as a new luminous blue variable candidate. The circumstellar dust shell seems to have been generated by an equatorially-enhanced mass loss process with ${\dot M} \ge 6.8 \times 10^{-6} \msun$ yr$^{-1}$ and ${\dot M}_{\rm pole}/{\dot M}_{\rm eq} \sim 0.5$.

This review describes some general properites of proto-planetary nebulae with particular emphasis on the recent work of morpholgical studies. The weight of observational evidence shows that proto-planetary nebulae (PPNe) are most certainly axisymmetric like planetary nebulae. Recent work suggests two subclasses of PPNe optical morphology, DUst-Prominent Longitudinally-EXtended (DUPLEX) and Star-Obvious Low-level Elongated (SOLE). Radiative transfer models of an example DUPLEX PPN and SOLE PPN, presented here, support the interpretation that DUPLEX and SOLE are two physically distinct types of PPNe. The DUPLEX PPNe and SOLE PPNe may well be the precursors to bipolar and elliptical PNe, respectively.

We summarize results of our {\it HST} imaging survey of proto-planetary nebula (PPN) candidates, in which we discovered two types of axisymmetric reflection nebulosities. The Star-Obvious Low-level-Elongated (SOLE) nebulae show a bright central star embedded in a smooth, faint extended nebulosity, whereas the DUst-Prominent Longitudinally-EXtended (DUPLEX) nebulae have remarkable bipolar structure with a completely or partially obscured central star. The intrinsic axisymmetry of the PPN reflection nebulosities demonstrates that the axisymmetry often found in planetary nebulae predates the PPN phase. We suggest that the major cause for the apparent dual morphology is the optical depth difference in the circumstellar dust shell rather than the inclination angle effect.

I-band images of the OH maser source OH44.8-2.3 are presented. Having a well-determined distance, this object is shown to be very young and physically small, with an asymmetric compact bipolar nebula already present. In the morphological terminology of Meixner et al. (2000) and Ueta et al. (2000), this object would be classified as a DUPLEX nebula because of its high optical depth and pole-to-equator density contrast.

We report the results from an optical imaging survey of proto-planetary nebula candidates using the HST. We exploited the high resolving power and wide dynamic range of HST and detected nebulosities in 21 of 27 sources. All detected reflection nebulosities show elongation, and the nebula morphology bifurcates depending on the degree of the central star obscuration. The Star-Obvious Low-level-Elongated (SOLE) nebulae show a bright central star embedded in a faint, extended nebulosity, whereas the DUst-Prominent Longitudinally-EXtended (DUPLEX) nebulae have remarkable bipolar structure with a completely or partially obscured central star. The intrinsic axisymmetry of these proto-planetary nebula reflection nebulosities demonstrates that the axisymmetry frequently found in planetary nebulae predates the proto-planetary nebula phase, confirming previous independent results. We suggest that axisymmetry in proto-planetary nebulae is created by an equatorially enhanced superwind at the end of the asymptotic giant branch phase. We discuss that the apparent morphological dichotomy is caused by a difference in the optical thickness of the circumstellar dust/gas shell with a differing equator-to-pole density contrast. Moreover, we show that SOLE and DUPLEX nebulae are physically distinct types of proto-planetary nebulae, with a suggestion that higher mass progenitor AGB stars are more likely to become DUPLEX proto-planetary nebulae.

We present the data from a mid-infrared imaging survey of 66 proto-planetary nebula candidates using two mid-IR cameras (MIRAC2 and Berkcam) at the NASA Infrared Telescope Facility and the United Kingdom Infrared Telescope. The goal of this survey is to determine the size, flux, and morphology of the mid-IR emission regions, which sample the inner regions of the circumstellar dust shells of proto-planetary nebulae. We imaged these proto-planetary nebulae with narrow-band filters ($\Delta\lambda / \lambda \sim 10%$) at wavelengths of notable dust features. With our typical angular resolution of 1\arcsec, we resolve 17 sources, find 48 objects unresolved, and do not detect 1 source. For several sources, we checked optical and infrared associations and positions of the sources. In table format, we list the size and flux measurements for all the detected objects and show figures of all the resolved sources. Images for all the detected objects are available on line in FITS format from the Astronomy Digital Image Library at the National Center for Supercomputing Application. The proto-planetary nebula candidate sample includes, in addition to the predominant proto-planetary nebulae, extreme asymptotic giant branch stars, young planetary nebulae, a supergiant, and a luminous blue variable. We find that dust shells which are cooler ($\rm T \sim 150$ K) and brighter in the infrared are more easily resolved. Eleven of the seventeen resolved sources are extended and fall into one of two types of mid-IR morphological classes: core/elliptical or toroidal. Core/elliptical structures show unresolved cores with lower surface brightness elliptical nebulae. Toroidal structures show limb-brightened peaks suggesting equatorial density enhancements. We argue that core/ellipticals have denser dust shells than toroidals.

We present arcsecond resolution mid-infrared (8--13 $\mu$m) images and photometry of four young stellar objects (YSOs)~--- L1551-IRS5, HL~Tau, AS~205, and AS~209 (V1121~Oph)~--- taken with the Berkeley Mid-Infrared Camera. For AS~205, a known T Tauri binary, we also present near-infrared JHK images and HKL$^{\prime}$ speckle imaging data. All three single stars are unresolved in our mid-IR images. Our data is the first to resolve in the mid-IR both components of the close binary AS~205 (projected separation $\sim$1.3$^{\prime \prime}$ (210~A.U.)). AS~205~North is the IR brighter star in our data while published observations find it to be the optically fainter star. Pre-main sequence evolutionary models suggest the AS~205 system is non-coeval; we discuss possible explanations for this result and comment on the circumstellar material and evolutionary status of this young binary. Nearly all of our objects exhibit changes in their mid-IR flux in measurements separated by intervals of days up to many years; the variations range from 30--300\%. We speculate that the cause of the variability lies in the accretion disks of these objects; the data suggest disk accretion rate fluctuations of nearly an order of magnitude. The existence of large mid-IR variability argues that simultaneous multiwavelength observations are needed for a proper analysis of YSO spectral energy distributions.