We present a detailed analysis of deep VLT/X-Shooter observations of the planetary nebula Tc 1. We calculate gas temperature, density, extinction, and abundances for several species from the empirical analysis of the total line fluxes. In addition, a spatially resolved analysis of the most intense lines provides the distribution of such quantities across the nebula. The new data reveal that several lines exhibit a double peak spectral profile consistent with the blue- and red-shifted components of an expanding spherical shell. The study of such components allowed us to construct for the first time a three-dimensional morphological model, which reveals that Tc 1 is a slightly elongated spheroid with an equatorial density enhancement seen almost pole on. A few bright lines present extended wings (with velocities up to a few hundred km/s), but the mechanism producing them is not clear. We constructed photoionization models for the main shell of Tc 1. The models predict the central star temperature and luminosity, as well as the nebular density and abundances similar to previous studies. Our models indicate that Tc 1 is located at a distance of approximately 2 kpc. We report the first detection of the [Kr III] 6825 A emission line, from which we determine the Krypton abundance. Our model indicates that the main shell of Tc 1 is matter bounded; leaking H ionizing photons may explain the ionization of its faint AGB-remnant halo.

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.

Observations towards L1630 in the Orion B molecular cloud, comprising the iconic Horsehead Nebula, allow us to study the interplay between stellar radiation and a molecular cloud under relatively benign conditions, that is, intermediate densities and an intermediate UV radiation field. Contrary to the well-studied Orion Molecular Cloud 1 (OMC1), which hosts much harsher conditions, L1630 has little star formation. We aim to relate the [CII] fine-structure line emission to the physical conditions predominant in L1630 and compare it to studies of OMC1. The [CII] $158\,\mu\mathrm{m}$ emission from an area of $12' \times 17'$ in L1630 was observed using the upGREAT instrument onboard SOFIA. Of the [CII] emission from the mapped area 95%, $13\,L_{\odot}$, originates from the molecular cloud; the adjacent HII region contributes only 5%, that is, $1\,L_{\odot}$. From comparison with other data (CO (1-0)-line emission, far-infrared (FIR) continuum studies, emission from polycyclic aromatic hydrocarbons (PAHs)), we infer a gas density of the molecular cloud of $n_{\mathrm{H}}\sim 3\cdot 10^3\,\mathrm{cm^{-3}}$, with surface layers, including the Horsehead Nebula, having a density of up to $n_{\mathrm{H}}\sim 4\cdot 10^4\,\mathrm{cm^{-3}}$. The temperature of the surface gas is $T\sim 100\,\mathrm{K}$. The average [CII] cooling efficiency within the molecular cloud is $1.3\cdot 10^{-2}$. The fraction of the mass of the molecular cloud within the studied area that is traced by [CII] is only $8\%$. Our PDR models are able to reproduce the FIR-[CII] correlations and also the CO (1-0)-[CII] correlations. Finally, we compare our results on the heating efficiency of the gas with theoretical studies of photoelectric heating by PAHs, clusters of PAHs, and very small grains, and find the heating efficiency to be lower than theoretically predicted, a continuation of the trend set by other observations.

During its 2012 decline the R Coronae Borealis star (RCB) V854 Cen was spectroscopically monitored with X-shooter on the ESO Very Large Telescope. The obscured optical and near-infrared spectrum exhibits many narrow and several broad emission features, as previously observed. The envelope is spatially resolved along the slit and allows for a detailed study of the circumstellar material. In this Letter we report on the properties of a number of unidentified emission features (UFs), including the detection of a new one at 8692 {\AA}. These UFs have been observed in the Red Rectangle, but their chemical and physical nature is still a mystery. The previously known UFs behave similarly in the Red Rectangle and V854 Cen, but are not detected in six other observed RCBs. Possibly the presence of some hydrogen is required for the formation of their carrier(s). The $\lambda 8692$ UF is present in all RCBs. Its carrier is likely of a carbonaceous molecular nature, presumably different from that of the other UFs.

The formation process of massive stars is still poorly understood. Massive young stellar objects (mYSOs) are deeply embedded in their parental clouds, they are rare and thus typically distant, and their reddened spectra usually preclude the determination of their photospheric parameters. M17 is one of the best studied HII regions in the sky, is relatively nearby, and hosts a young stellar population. With X-shooter on the ESO Very Large Telescope we have obtained optical to near-infrared spectra of candidate mYSOs, identified by Hanson et al. (1997), and a few OB stars in this region. The large wavelength coverage enables a detailed spectroscopic analysis of their photospheres and circumstellar disks. We confirm the pre-main sequence (PMS) nature of six of the stars and characterise the O stars. The PMS stars have radii consistent with being contracting towards the main sequence and are surrounded by a remnant accretion disk. The observed infrared excess and the (double-peaked) emission lines provide the opportunity to measure structured velocity profiles in the disks. We compare the observed properties of this unique sample of young massive stars with evolutionary tracks of massive protostars by Hosokawa & Omukai (2009), and propose that these mYSOs near the western edge of the HII region are on their way to become main-sequence stars ($\sim 6 - 20$ $M_{\odot}$) after having undergone high mass-accretion rates (${\dot{M}_{\rm acc}} \sim 10^{-4} - 10^{-3}$ $M_{\odot}$ $\rm yr^{-1}$). Their spin distribution upon arrival at the zero age main sequence (ZAMS) is consistent with that observed for young B stars, assuming conservation of angular momentum and homologous contraction.

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.

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.

The Orion star-forming region is the nearest active high-mass star-forming region and has created a large superbubble, the Orion-Eridanus superbubble. Recent work by Ochsendorf et al. has extended the accepted boundary of the superbubble. We fit Kompaneets models of superbubbles expanding in exponential atmospheres to the new, larger shape of the Orion-Eridanus superbubble. We find that this larger morphology of the superbubble is consistent with the evolution of the superbubble being primarily controlled by expansion into the exponential Galactic disk ISM if the superbubble is oriented with the Eridanus side farther from the Sun than the Orion side. Unlike previous Kompaneets model fits that required abnormally small scale heights for the Galactic disk (<40 pc), we find morphologically consistent models with scale heights of 80 pc, similar to that expected for the Galactic disk.

We performed a detailed spectroscopic analysis of the fullerene C60-containing planetary nebula (PN) Lin49 in the Small Magellanic Cloud using XSHOOTER at the ESO VLT and the Spitzer/IRS instruments. We derived nebular abundances for nine elements. We used TLUSTY to derive photospheric parameters for the central star. Lin49 is C-rich and metal-deficient PN (Z~0.0006). The nebular abundances are in good agreement with Asymptotic Giant Branch nucleosynthesis models for stars with initial mass 1.25 Msun and metallicity Z = 0.001. Using the TLUSTY synthetic spectrum of the central star to define the heating and ionising source, we constructed the photoionisation model with CLOUDY that matches the observed spectral energy distribution (SED) and the line fluxes in the UV to far-IR wavelength ranges simultaneously. We could not fit the ~1-5 um SED using a model with 0.005-0.1 um-sized graphite grains and a constant hydrogen density shell owing to the prominent near-IR excess, while at other wavelengths the model fits the observed values reasonably well. We argue that the near-IR excess might indicate either (1) the presence of very small particles in the form of small carbon clusters, small graphite sheets, or fullerene precursors, or (2) the presence of a high-density structure surrounding the central star. We found that SMC C60 PNe show a near-IR excess component to lesser or greater degree. This suggests that these C60 PNe might maintain a structure nearby their central star.

The Orion-Eridanus superbubble is the prototypical superbubble due to its proximity and evolutionary state. Here, we provide a synthesis of recent observational data from WISE and Planck with archival data, allowing to draw a new and more complete picture on the history and evolution of the Orion-Eridanus region. We discuss the general morphological structures and observational characteristics of the superbubble, and derive quantitative properties of the gas- and dust inside Barnard's Loop. We reveal that Barnard's Loop is a complete bubble structure which, together with the lambda Ori region and other smaller-scale bubbles, expands within the Orion-Eridanus superbubble. We argue that the Orion-Eridanus superbubble is larger and more complex than previously thought, and that it can be viewed as a series of nested shells, superimposed along the line of sight. During the lifetime of the superbubble, HII region champagne flows and thermal evaporation of embedded clouds continuously mass-load the superbubble interior, while winds or supernovae from the Orion OB association rejuvenate the superbubble by sweeping up the material from the interior cavities in an episodic fashion, possibly triggering the formation of new stars that form shells of their own. The steady supply of material into the superbubble cavity implies that dust processing from interior supernova remnants is more efficient than previously thought. The cycle of mass-loading, interior cleansing, and star formation repeats until the molecular reservoir is depleted or the clouds have been disrupted. While the nested shells come and go, the superbubble remains for tens of millions of years.

Recent studies of dust in the interstellar medium have challenged the capabilities and validity of current dust models, indicating that the properties of dust evolve as it transits between different phases of the interstellar medium. We conduct a multi-wavelength study of the dust emission from the ionized gas of the IC 434 emission nebula, and combine this with modeling, from large scales that provide insight into the history of the IC 434/L1630 region, to small scales that allow us to infer quantitative properties of the dust content inside the H II region. The dust enters the H II region through momentum transfer with a champagne flow of ionized gas, set up by a chance encounter between the L1630 molecular cloud and the star cluster of $\sigma$ Ori. We observe two clearly separated dust populations inside the ionized gas, that show different observational properties, as well as contrasting optical properties. Population A is colder ($\sim$ 25 K) than predicted by widely-used dust models, its temperature is insensitive to an increase of the impinging radiation field, is momentum-coupled to the gas, and efficiently absorbs radiation pressure to form a dust wave at 1.0 pc ahead of $\sigma$ Ori AB. Population B is characterized by a constant [20/30] flux ratio throughout the HII region, heats up to $\sim$ 75 K close to the star, and is less efficient in absorbing radiation pressure, forming a dust wave at 0.1 pc from the star. We conclude that the dust inside IC 434 is bimodal. The characteristics of population A are remarkable and can not be explained by current dust models. Population B are grains that match the classical description of spherical, compact dust. Our results confirm recent work that stress the importance of variations in the dust properties between different regions of the interstellar medium.

The freezeout of gas-phase species onto cold dust grains can drastically alter the chemistry and the heating-cooling balance of protostellar material. In contrast to well-known species such as carbon monoxide (CO), the freezeout of various carriers of elements with abundances $<10^{-5}$ has not yet been well studied. Our aim here is to study the depletion of chlorine in the protostellar core, OMC-2 FIR 4. We observed transitions of HCl and H2Cl+ towards OMC-2 FIR 4 using the Herschel Space Observatory and Caltech Submillimeter Observatory facilities. Our analysis makes use of state of the art chlorine gas-grain chemical models and newly calculated HCl-H$_{2}$ hyperfine collisional excitation rate coefficients. A narrow emission component in the HCl lines traces the extended envelope, and a broad one traces a more compact central region. The gas-phase HCl abundance in FIR 4 is 9e-11, a factor of only 0.001 that of volatile elemental chlorine. The H2Cl+ lines are detected in absorption and trace a tenuous foreground cloud, where we find no depletion of volatile chlorine. Gas-phase HCl is the tip of the chlorine iceberg in protostellar cores. Using a gas-grain chemical model, we show that the hydrogenation of atomic chlorine on grain surfaces in the dark cloud stage sequesters at least 90% of the volatile chlorine into HCl ice, where it remains in the protostellar stage. About 10% of chlorine is in gaseous atomic form. Gas-phase HCl is a minor, but diagnostically key reservoir, with an abundance of <1e-10 in most of the protostellar core. We find the 35Cl/37Cl ratio in OMC-2 FIR 4 to be 3.2\pm0.1, consistent with the solar system value.

This paper presents a spectral survey of diffuse interstellar bands (DIBs) in the NIR range, from 0.9 to 2.5 micron. The observations were designed to detect new DIBs, confirm previously proposed NIR DIBs, and characterise their behaviour with respect to known line-of-sight properties (including the optical DIBs present in our spectra). X-shooter at the VLT was used to obtained medium-resolution spectra of eight known DIB targets and one telluric reference star. In addition to the known 9577, 9632, 10780, 11797, and 13175 Angstroms NIR DIBs, we confirm 9 out of the 13 NIR DIBs that were presented by Geballe and co-workers in 2011. Furthermore, we report 12 new NIR DIB candidates. The strengths of the strongest NIR DIBs show a general correlation with reddening, E(B-V), but with a large scatter. Several NIR DIBs are more strongly correlated with the 5780 Angstroms DIB strength than with E(B-V); this is especially the case for the 15268 Angstroms DIB. The NIR DIBs are strong: the summed equivalent widths of the five strongest NIR DIBs represent a few percent of the total equivalent width of the entire average DIB spectrum (per unit reddening). The NIR DIBs towards the translucent cloud HD 147889 are all weak with respect to the general trend. No direct match was found between observed NIR DIBs and laboratory matrix-isolation spectroscopic data of PAHs. The strong correlation between the 5780-15268 DIB pair implies that (Nf)_5780 / (Nf)_15268 = 14. However, the reduced strength of the 15268 Angstroms DIB in HD 147889 rules out a common carrier for these two DIBs. Since the ionisation fraction for small PAHs in this translucent cloud is known to be low compared to diffuse clouds, the weakness of the 15268 Angstrom DIB suggests that an ionised species could be the carrier of this NIR DIB. (abridged).

Observations with the Spitzer Space Telescope and the WISE satellite have revealed a prominent arc-like structure at 50" ($\simeq0.1$ pc) from the O9.5V/B0.5V system $\sigma$ Ori AB. We attribute this dust structure to the interaction of radiation pressure from the star with dust carried along by the IC 434 photo-evaporative flow of ionized gas from the dark cloud L1630. We have developed a quantitative model for the interaction of a dusty ionized flow with nearby (massive) stars where radiation pressure stalls dust, piling it up at an appreciable distance (> 0.1 pc), and force it to flow around the star. The model demonstrates that for the conditions in IC 434, the gas will decouple from the dust and will keep its original flow lines. We argue that this dust structure is the first example of a dust wave created by a massive star moving through the interstellar medium. Dust waves (and bow waves) stratify dust grains according to their radiation pressure opacity, which reflects the size distribution and composition of the grain material. Comparison of our model with observations implies that dust-gas coupling through Coulomb interaction is less important than previously thought, challenging our understanding of grain dynamics in hot, ionized regions of space. We describe the difference between dust (and bow) waves and classical bow shocks. We conclude that dust waves and bow waves should be common around stars showing the weak-wind phenomenon and that these structures are best observed at mid-IR to FIR wavelengths. In particular, dust waves and bow waves are most efficiently formed around weak-wind stars moving through a high density medium. Moreover, they provide a unique opportunity to study the direct interaction between a (massive) star and its immediate surroundings.

Massive stars drive the evolution of the interstellar medium through their radiative and mechanical energy input. After their birth, they form bubbles of hot gas surrounded by a dense shell. Traditionally, the formation of bubbles is explained through the input of a powerful stellar wind, even though direct evidence supporting this scenario is lacking. Here we explore the possibility that interstellar bubbles seen by the Spitzer- and Herschel space telescopes, blown by stars with log(L/L_sun) < 5.2, form and expand due to the thermal pressure accompanying ionization of the surrounding gas. We show that density gradients in the natal cloud or a puncture in the swept up shell lead to an ionized gas flow through the bubble into the general interstellar medium, which is traced by a dust wave near the star, demonstrating the importance of radiation pressure during this phase. Dust waves provide a natural explanation for the presence of dust inside H II bubbles, offer a novel method to study dust in H II regions and provide direct evidence that bubbles are relieving their pressure into the ISM through a champagne flow, acting as a probe of the radiative interaction of a massive star with its surroundings. We create a parameter space connecting the ambient density, the ionizing source luminosity, and the position of the dust wave, while using the well-studied H II bubbles RCW 120 and RCW 82 as benchmarks of our model. Finally, we briefly examine the implications of our study for the environments of super star clusters formed in UltraLuminous InfraRed Galaxies (ULIRGs), merging galaxies and the early Universe, which occur in very luminous and dense environments and where radiation pressure is expected to dominate the dynamical evolution.

The optical to near-infrared (300 - 2500 nm) spectrum of the candidate massive Young Stellar Object (YSO) B275, embedded in the star forming region M17, has been observed with X-shooter on the ESO Very Large Telescope. The spectrum includes both photospheric absorption lines and emission features (H and Ca triplet emission lines, 1st and 2nd overtone CO bandhead emission), as well as an infrared excess indicating the presence of a (flaring) circumstellar disk. The strongest emission lines are double-peaked with a peak separation ranging between 70 and 105 km s^-1 and provide information on the physical structure of the disk. The underlying photospheric spectrum is classified as B6-B7, significantly cooler than a previous estimate based on modeling of the spectral energy distribution. This discrepancy is solved by allowing for a larger stellar radius (i.e. a bloated star), and thus positioning the star above the main sequence. This constitutes the first firm spectral classification of an early-B pre-main-sequence (PMS) star. We discuss the position of B275 in the Hertzsprung-Russell diagram in terms of PMS evolution. Although the position is consistent with PMS tracks of heavily accreting protostars (Macc > 10^-5 Msun yr^-1), the fact that the photosphere of the object is detectable suggests that the current mass-accretion rate is not very high.

Based on the SED, Herbig stars have been categorized into two observational groups, reflecting their overall disk structure: group I members have disks with a higher degree of flaring than their group II counterparts. We investigate the 5-35 um Spitzer IRS spectra of a sample of 13 group I sources and 20 group II sources. We focus on the continuum emission to study the underlying disk geometry. We have determined the [30/13.5] and [13.5/7] continuum flux ratios. The 7-um flux excess with respect to the stellar photosphere is measured, as a marker for the strength of the near-IR emission produced by the inner disk. We have compared our data to self-consistent passive-disk model spectra, for which the same quantities were derived. We confirm the literature result that the difference in continuum emission between group I and II sources can largely be explained by a different amount of small dust grains. However, we report a strong correlation between the [30/13.5] and [13.5/7] flux ratios for Meeus group II sources. Moreover, the [30/13.5] flux ratio decreases with increasing 7-um excess for all targets in the sample. To explain these correlations with the models, we need to introduce an artificial scaling factor for the inner disk height. In roughly 50% of the Herbig Ae/Be stars in our sample, the inner disk must be inflated by a factor 2 to 3 beyond what hydrostatic calculations predict. The total disk mass in small dust grains determines the degree of flaring. We conclude, however, that for any given disk mass in small dust grains, the shadowing of the outer (tens of AU) disk is determined by the scale height of the inner disk (1 AU). The inner disk partially obscures the outer disk, reducing the disk surface temperature. Here, for the first time, we prove these effects observationally.