Observations suggest that there is a significant fraction of O-stars in the field of the Milky Way that appear to have formed in isolation or in low mass clusters ($<$100 $M_\odot$). The existence of these high-mass stars that apparently formed in the field challenges the generally accepted paradigm, which requires star formation to occur in clustered environments. In order to understand the physical conditions for the formation of these stars, it is necessary to observe isolated high-mass stars while they are still forming. With the $Hubble$ $Space$ $Telescope$, we observe the seven most isolated massive ($>$8 $M_\odot$) young stellar objects (MYSOs) in the Large Magellanic Cloud (LMC). The observations show that while these MYSOs are remote from other MYSOs, OB associations, and even from known giant molecular clouds, they are actually not isolated at all. Imaging reveals $\sim$100 to several hundred pre--main-sequence (PMS) stars in the vicinity of each MYSO. These previously undetected PMS stars form prominent compact clusters around the MYSOs, and in most cases they are also distributed sparsely across the observed regions. Contrary to what previous high-mass field star studies show, these observations suggest that high-mass stars may not be able to form in clusters with masses less than 100 $M_\odot$. If these MYSOs are indeed the best candidates for isolated high-mass star formation, then the lack of isolation is at odds with random sampling of the IMF. Moreover, while isolated MYSOs may not exist, we find evidence that isolated clusters containing O-stars can exist, which in itself is rare.

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.

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 spatial variations of the gas-to-dust ratio (GDR) provide constraints on the chemical evolution and lifecycle of dust in galaxies. We examine the relation between dust and gas at 10-50 pc resolution in the Large and Small Magellanic Clouds (LMC and SMC) based on Herschel far-infrared (FIR), H I 21 cm, CO, and Halpha observations. In the diffuse atomic ISM, we derive the gas-to-dust ratio as the slope of the dust-gas relation and find gas-to-dust ratios of 380+250-130 in the LMC, and 1200+1600-420 in the SMC, not including helium. The atomic-to-molecular transition is located at dust surface densities of 0.05 Mo pc-2 in the LMC and 0.03 Mo pc-2 in the SMC, corresponding to AV ~ 0.4 and 0.2, respectively. We investigate the range of CO-to-H2 conversion factor to best account for all the molecular gas in the beam of the observations, and find upper limits on XCO to be 6x1020 cm-2 K-1 km-1 s in the LMC (Z=0.5Zo) at 15 pc resolution, and 4x 1021 cm-2 K-1 km-1 s in the SMC (Z=0.2Zo) at 45 pc resolution. In the LMC, the slope of the dust-gas relation in the dense ISM is lower than in the diffuse ISM by a factor ~2, even after accounting for the effects of CO-dark H2 in the translucent envelopes of molecular clouds. Coagulation of dust grains and the subsequent dust emissivity increase in molecular clouds, and/or accretion of gas-phase metals onto dust grains, and the subsequent dust abundance (dust-to-gas ratio) increase in molecular clouds could explain the observations. In the SMC, variations in the dust-gas slope caused by coagulation or accretion are degenerate with the effects of CO-dark H2. Within the expected 5--20 times Galactic XCO range, the dust-gas slope can be either constant or decrease by a factor of several across ISM phases. Further modeling and observations are required to break the degeneracy between dust grain coagulation, accretion, and CO-dark H2.

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

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.

We present parsec scale interferometric maps of HCN(1-0) and HCO$^{+}$(1-0) emission from dense gas in the star-forming region 30 Doradus, obtained using the Australia Telescope Compact Array. This extreme star-forming region, located in the Large Magellanic Cloud (LMC), is characterized by a very intense ultraviolet ionizing radiation field and sub-solar metallicity, both of which are expected to impact molecular cloud structure. We detect 13 bright, dense clumps within the 30 Doradus-10 giant molecular cloud. Some of the clumps are aligned along a filamentary structure with a characteristic spacing that is consistent with formation via the varicose fluid instability. Our analysis shows that the filament is gravitationally unstable and collapsing to form stars. There is a good correlation between HCO$^{+}$ emission in the filament and signatures of recent star formation activity including H$_{2}$O masers and young stellar objects (YSOs). YSOs seem to continue along the same direction of the filament toward the massive compact star cluster R136 in the southwest. We present detailed comparisons of clump properities (masses, linewidths, sizes) in 30Dor-10 to those in other star forming regions of the LMC (N159, N113, N105, N44). Our analysis shows that the 30Dor-10 clumps have similar mass but wider linewidths and similar HCN/HCO$^{+}$ (1-0) line ratios as clumps detected in other LMC star-forming regions. Our results suggest that the dense molecular gas clumps in the interior of 30Dor-10 are well-shielded against the intense ionizing field that is present in the 30Doradus region.

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.

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 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 report the results of a high spatial (parsec) resolution HCO+ (J = 1-0) and HCN (J = 1-0) emission survey toward the giant molecular clouds of the star formation regions N105, N113, N159, and N44 in the Large Magellanic Cloud. The HCO+ and HCN observations at 89.2 and 88.6 GHz, respectively, were conducted in the compact configuration of the Australia Telescope Compact Array. The emission is imaged into individual clumps with masses between 10^2 and 10^4 solar masses and radii of <1 pc to ~2 pc. Many of the clumps are coincident with indicators of current massive star formation, indicating that many of the clumps are associated with deeply-embedded forming stars and star clusters. We find that massive YSO-bearing clumps tend to be larger (>1 pc), more massive (M > 10^3 solar masses), and have higher surface densities (~1 g cm^-2), while clumps without signs of star formation are smaller (<1 pc), less massive (M < 10^3 solar masses), and have lower surface densities (~0.1 g cm^-2). The dearth of massive (M >10^3 solar masses) clumps not bearing massive YSOs suggests the onset of star formation occurs rapidly once the clump has attained physical properties favorable to massive star formation. Using a large sample of LMC massive YSO mid-IR spectra, we estimate that ~2/3 of the massive YSOs for which there are Spitzer mid-IR spectra are no longer located in molecular clumps; we estimate that these young stars/clusters have destroyed their natal clumps on a time scale of at least 3 x 10^{5}$ yrs.

We present the properties of an extensive sample of molecular clouds in the Large Magellanic Cloud (LMC) mapped at 11 pc resolution in the CO(1-0) line. We identify clouds as regions of connected CO emission, and find that the distributions of cloud sizes, fluxes and masses are sensitive to the choice of decomposition parameters. In all cases, however, the luminosity function of CO clouds is steeper than dN/dL \propto L^{-2}, suggesting that a substantial fraction of mass is in low-mass clouds. A correlation between size and linewidth, while apparent for the largest emission structures, breaks down when those structures are decomposed into smaller structures. We argue that the correlation between virial mass and CO luminosity is the result of comparing two covariant quantities, with the correlation appearing tighter on larger scales where a size-linewidth relation holds. The virial parameter (the ratio of a cloud's kinetic to self-gravitational energy) shows a wide range of values and exhibits no clear trends with the CO luminosity or the likelihood of hosting young stellar object (YSO) candidates, casting further doubt on the assumption of virialization for molecular clouds in the LMC. Higher CO luminosity increases the likelihood of a cloud harboring a YSO candidate, and more luminous YSOs are more likely to be coincident with detectable CO emission, confirming the close link between giant molecular clouds and massive star formation.

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 present and categorize Spitzer IRS spectra of 294 objects in the Large Magellanic Cloud (LMC) to create the largest and most complete catalog of massive young stellar object (YSO) spectra in the LMC. Target sources were identified from infrared photometry and multi-wavelength images indicative of young, massive stars highly enshrouded in their natal gas and dust clouds. Our sample primarily consists of 277 objects we identify as having spectral features indicative of embedded YSOs. The remaining sources are comprised of 7 C-rich evolved sources, 8 sources dominated by broad silicate emission, and 1 source with multiple broad emission features. Those with YSO-like spectra show a range of spectral features including polycyclic aromatic hydrocarbon emission, deep silicate absorption, fine-structure lines, and ice absorption features. Based upon the relative strengths of these features, we have classified the YSO candidates into several distinct categories using the widely-used statistical procedure known as principal component analysis. We propose that these categories represent a spectrum of evolutionary stages during massive YSO formation. We conclude that massive pre-main sequence stars spend a majority of their massive, embedded lives emitting in the UV. Half of the sources in our study have features typical of compact HII regions, suggesting that massive YSOs can create a detectable compact HII region half-way through the formation time present in our sample. This study also provides a check on commonly used source-selection procedures including the use of photometry to identify YSOs. We determine a high success rate (>95%) of identifying objects with YSO-like spectra can be achieved through careful use of infrared CMDs, SEDs, and image inspections.

Using Spitzer IRAC and MIPS observations of the Large Magellanic Cloud, we have identified 13 objects that have extremely red mid-IR colors. Follow-up Spitzer IRS observations of seven of these sources reveal varying amounts of SiC and C2H2 absorption as well as the presence of a broad MgS feature in at least two cases, indicating that these are extreme carbon stars. Preliminary estimates find these objects have luminosities of 4-11x10^3 Lsol and preliminary model fitting gives mass-loss rates between 4x10^-5 and 2x10^-4 Msol/yr, higher than any known carbon-rich AGB star in the LMC. These spectral and physical properties require careful reconsideration of dust condensation and mass-loss processes for carbon stars in low metallicity environments.

We present Spitzer IRAC images that indicate the presence of cavities cut into the dense outer envelope surrounding very young pre-main sequence stars. These young stellar objects (YSOs) characterized by an outflow represent the earliest stages of star formation. Mid-infrared photons thermally created by the central protostar/disk are scattered by dust particles within the outflow cavity itself into the line of sight. We observed this scattered light from 27 nearby, cavity-resolved YSOs, and quantified the shape of the outflow cavities. Using the grid models of Robitaille et al. (2006), we matched model spectral energy distributions (SEDs) to the observed SEDs of the 27 cataloged YSOs using photometry from IRAC, MIPS, and IRAS. This allows for the estimation of geometric and physical properties such as inclination angle, cavity density, and accretion rate. By using the relative parameter estimates determined by the models, we are able to deduce an evolutionary picture for outflows. Our work supports the concept that cavities widen with time, beginning as a thin jet-like outflow that widens to reveal the central protostar and disk until the protostellar envelope is completely dispersed by outflow and accretion.