Search

The evolution of protoplanetary disks in regions with massive OB stars is influenced by externally driven winds that deplete the outer parts of disks. These winds have previously been studied via forbidden oxygen emission lines, which also arise in isolated disks in low-mass star forming-regions (SFRs) with weak external UV fields in photoevaporative or magnetic (internal) disk winds. It is crucial to determine how to disentangle external winds from internal ones. Here, we report a proxy for unambiguously identifying externally driven winds with a forbidden line of neutral atomic carbon, [C i] 8727 A. We compare for the first time the spatial location of the emission in the [O i] 5577 A, [O i] 6300 A, and [C i] 8727 A lines traced by VLT/MUSE-NFM, with the ALMA Band 7 continuum disk emission in a sample of 12 proplyds in the Orion Nebula Cluster (ONC). We confirm that the [O i] 5577 A emission is co-spatial with the disk emission, whereas the [O i] 6300 A is emitted both on the disk surface and on the ionization front of the proplyds. We show for the first time that the [C i] 8727 A line is also co-spatial with the disk surface in proplyds, as seen in the MUSE and ALMA data comparison. To verify whether the [C i] 8727 A line is detected in regions where external photoevaporation is not expected, we examine VLT/X-Shooter spectra for young stars in low-mass SFRs. Although the [O i] lines are well detected in all these targets, there is <<10% detection rate in the case of the [C i] 8727 A line. This number increases substantially to a ~40% detection rate in sigma-Orionis, a region with intermediate UV radiation. The spatial location of the [C i] 8727 A line emission and the lack of its detection in isolated disks in low-mass SFRs strongly suggest that this line is a tell-trace tracer of externally driven photoevaporative winds, which agrees with recent excitation models., Comment: Accepted to A&A on October 25, 2024

We present new JWST observations of the nearby, prototypical edge-on, spiral galaxy NGC 891. The northern half of the disk was observed with NIRCam in its F150W and F277W filters. Absorption is clearly visible in the mid-plane of the F150W image, along with vertical dusty plumes that closely resemble the ones seen in the optical. A $\sim 10 \times 3~{\rm kpc}^2$ area of the lower circumgalactic medium (CGM) was mapped with MIRI F770W at 12 pc scales. Thanks to the sensitivity and resolution of JWST, we detect dust emission out to $\sim 4$ kpc from the disk, in the form of filaments, arcs, and super-bubbles. Some of these filaments can be traced back to regions with recent star formation activity, suggesting that feedback-driven galactic winds play an important role in regulating baryonic cycling. The presence of dust at these altitudes raises questions about the transport mechanisms at play and suggests that small dust grains are able to survive for several tens of million years after having been ejected by galactic winds in the disk-halo interface. We lay out several scenarios that could explain this emission: dust grains may be shielded in the outer layers of cool dense clouds expelled from the galaxy disk, and/or the emission comes from the mixing layers around these cool clumps where material from the hot gas is able to cool down and mix with these cool cloudlets. This first set of data and upcoming spectroscopy will be very helpful to understand the survival of dust grains in energetic environments, and their contribution to recycling baryonic material in the mid-plane of galaxies., Comment: Accepted for publication in Astronomy & Astrophysics; 16 pages, 8 figures

BlueMUSE is a blue-optimised, medium spectral resolution, panoramic integral field spectrograph under development for the Very Large Telescope (VLT). With an optimised transmission down to 350 nm, spectral resolution of R$\sim$3500 on average across the wavelength range, and a large FoV (1 arcmin$^2$), BlueMUSE will open up a new range of galactic and extragalactic science cases facilitated by its specific capabilities. The BlueMUSE consortium includes 9 institutes located in 7 countries and is led by the Centre de Recherche Astrophysique de Lyon (CRAL). The BlueMUSE project development is currently in Phase A, with an expected first light at the VLT in 2031. We introduce here the Top Level Requirements (TLRs) derived from the main science cases, and then present an overview of the BlueMUSE system and its subsystems fulfilling these TLRs. We specifically emphasize the tradeoffs that are made and the key distinctions compared to the MUSE instrument, upon which the system architecture is built., Comment: 20 pages, 10 figures, proceedings of the SPIE astronomical telescopes and instrumentation conference, Yokohama, 16-21 June

We present the results of direct-method metallicity measurements in the disk and outflow of the low-metallicity starburst galaxy NGC 1569. We use Keck Cosmic Web Imager observations to map the galaxy across 54$\arcsec$ (800 pc) along the major axis and 48$\arcsec$ (700 pc) along the minor axis with a spatial resolution of 1$\arcsec$ ($\sim$15 pc). We detect common strong emission lines ([\ion{O}{III}] $\lambda$5007, H$\beta$, [\ion{O}{II}] $\lambda$3727) and the fainter [\ion{O}{III}] $\lambda$4363 auroral line, which allows us to measure electron temperature ($T_e$) and metallicity. Theory suggests that outflows drive metals out of the disk driving observed trends between stellar mass and gas-phase metallicity. Our main result is that the metallicity in the outflow is similar to that of the disk, $Z_{\rm out} / Z_{\rm ISM} \approx 1$. This is consistent with previous absorption line studies in higher mass galaxies. Assumption of a mass-loading factor of $\dot{M}_{\rm out}/{\rm SFR}\sim3$ makes the metal-loading of NGC 1569 consistent with expectations derived from the mass-metallicity relationship. Our high spatial resolution metallicity maps reveal a region around a supermassive star cluster (SSC-B) with distinctly higher metallicity and higher electron density, compared to the disk. Given the known properties of SSC-B the higher metallicity and density of this region are likely the result of star formation-driven feedback acting on the local scale. Overall, our results are consistent with the picture in which metal-enriched winds pollute the circumgalactic medium surrounding galaxies, and thus connect the small-scale feedback processes to large-scale properties of galaxy halos., Comment: 15 pages, 11 figures, accepted by MNRAS

In the Orion Nebula Cluster (ONC), protoplanetary disks exhibit ionized gas clouds in the form of a striking teardrop shape as massive stars irradiate the disk material. We present the first spatially and spectrally resolved observations of 12 proplyds, using Integral Field Spectroscopy observations performed with the MUSE instrument in Narrow Field Mode (NFM) on the VLT. We present the morphology of the proplyds in seven emission lines and measure the radius of the ionization front (I-front) of the targets in four tracers, covering transitions of different ionization states for the same element. We also derive stellar masses for the targets. The measurements follow a consistent trend of increasing I-front radius for a decreasing strength of the far-UV radiation as expected from photoevaporation models. By analyzing the ratios of the I-front radii as measured in the emission lines of Ha, [OI] 6300, [OII] 7330, and [OIII] 5007, we observe the ionization stratification, that is, the most ionized part of the flow being the furthest from the disk (and closest to the UV source). The ratios of I-front radii scale in the same way for all proplyds in our sample regardless of the incident radiation. We show that the stratification can help constrain the densities near the I-front by using a 1D photoionization model. We derive the upper limits of photoevaporative mass-loss rates by assuming ionization equilibrium, and estimate values decreasing towards lower impinging radiation. We do not find a correlation between Mloss and stellar mass. The highest mass-loss rate is for the proplyd 244-440. These values of Mloss, combined with estimates of the disk mass with ALMA, confirm previous estimates of the short lifetime of these proplyds. This work demonstrates the potential of this dataset and offers a new set of observables to be used to test current and future models of external photoevaporation., Comment: Submitted to A&A on December 18, 2023. Accepted on April 2, 2024

The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integral field spectrograph (IFS). In scientific capability these requirements place WST far ahead of existing and planned facilities. Given the current investment in deep imaging surveys and noting the diagnostic power of spectroscopy, WST will fill a crucial gap in astronomical capability and work synergistically with future ground and space-based facilities. This white paper shows that WST can address outstanding scientific questions in the areas of cosmology; galaxy assembly, evolution, and enrichment, including our own Milky Way; origin of stars and planets; time domain and multi-messenger astrophysics. WST's uniquely rich dataset will deliver unforeseen discoveries in many of these areas. The WST Science Team (already including more than 500 scientists worldwide) is open to the all astronomical community. To register in the WST Science Team please visit https://www.wstelescope.com/for-scientists/participate, Comment: 194 pages, 66 figures. Comments are welcome (wstelescope@gmail.com)

We measure resolved (kiloparsec-scale) outflow properties in a sample of 10 starburst galaxies from the DUVET sample, using Keck/KCWI observations of H$\beta$ and [OIII]~$\lambda$5007. We measure $\sim450$ lines-of-sight that contain outflows, and use these to study scaling relationships of outflow velocity ($v_{\rm out}$), mass-loading factor ($\eta$; mass outflow rate per SFR) and mass flux ($\dot{\Sigma}_{\rm out}$; mass outflow rate per area) with co-located SFR surface density ($\Sigma_{\rm SFR}$) and stellar mass surface density ($\Sigma_{\ast}$). We find strong, positive correlations of $\dot{\Sigma}_{\rm out} \propto \Sigma_{\rm SFR}^{1.2}$ and $\dot{\Sigma}_{\rm out} \propto \Sigma_{\ast}^{1.7}$. We also find shallow correlations between $v_{\rm out}$ and both $\Sigma_{\rm SFR}$ and $\Sigma_{\ast}$. Our resolved observations do not suggest a threshold in outflows with $\Sigma_{\rm SFR}$, but rather we find that the local specific SFR ($\Sigma_{\rm SFR}/\Sigma_\ast$) is a better predictor of where outflows are detected. We find that outflows are very common above $\Sigma_{\rm SFR}/\Sigma_\ast\gtrsim 0.1$~Gyr$^{-1}$ and rare below this value. We argue that our results are consistent with a picture in which outflows are driven by supernovae, and require more significant injected energy in higher mass surface density environments to overcome local gravity. The correlations we present here provide a statistically robust, direct comparison for simulations and higher redshift results from JWST., Comment: 14 pages, 7 figures, plus 4 figures in appendix, submitted to MNRAS

We present the first Chandra X-ray observations of H72.97-69.39, a highly-embedded, potential super-star cluster (SSC) in its infancy located in the star-forming complex N79 of the Large Magellanic Cloud. We detect particularly hard, diffuse X-ray emission that is coincident with the young stellar object (YSO) clusters identified with JWST, and the hot gas fills cavities in the dense gas mapped by ALMA. The X-ray spectra are best fit with either a thermal plasma or power-law model, and assuming the former, we show that the X-ray luminosity of L_X = (1.5 +- 0.3)e34 erg/s is a factor of ~20 below the expectation for a fully-confined wind bubble. Our results suggest that stellar wind feedback produces diffuse hot gas in the earliest stages of massive star cluster formation and that wind energy can be lost quickly via either turbulent mixing followed by radiative cooling or by physical leakage., Comment: 10 pages, 6 figures, 2 tables, submitted to ApJ

Stellar feedback in dwarf galaxies remains, to date, poorly explored, yet is crucial to understanding galaxy evolution in the early Universe. In particular, pre-supernova feedback has recently been found to play a significant role in regulating and disrupting star formation in larger spiral galaxies, but it remains uncertain if it also plays this role in dwarfs. We study the ionised gas properties and stellar content of individual star-forming regions across three nearby, low-metallicity, dwarf starburst galaxies (J0921, KKH046, and Leo P) to investigate how massive stars influence their surroundings and how this influence changes as a function of environment. To achieve this, we extracted integrated spectra of 30 HII regions from archival VLT/MUSE integral field spectroscopic observations of these three dwarf starburst galaxies. We fitted the HII regions' main emission lines with Gaussian profiles to derive their oxygen abundances, electron densities, and luminosities, and we used the Stochastically Lighting Up Galaxies (SLUG) code to derive the stellar mass, age, and bolometric luminosity of the stellar populations driving the HII regions. We then quantified two pre-supernova stellar feedback mechanisms, namely the direct radiation pressure and photoionisation feedback, and explored how feedback strength varies with HII region properties. Our findings suggest that stellar feedback has less of an impact on evolved regions, with both the pressure of the ionised gas and the direct radiation pressure decreasing as a function of HII region size. We also find that these stellar feedback mechanisms are dependent on the metallicity of the HII regions. These findings extend results from stellar feedback studies of more massive star-forming galaxies to the low-mass, low-metallicity regime., Comment: The paper has been published in A&A

The Large Magellanic Cloud (LMC) is home to many HII regions, which may lead to significant outflows. We examine the LMC's multiphase gas ($T\sim10^{4-5}$ K) in HI, SII, SiIV, and CIV using 110 stellar sight lines from the HST's Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) program. We develop a continuum fitting algorithm based on the concept of Gaussian Process regression and identify reliable LMC interstellar absorption over $v_{\rm helio}=175-375$ km s$^{-1}$. Our analyses show disk-wide ionized outflows in SiIV and CIV across the LMC with bulk velocities of $|v_{\rm out, bulk}|\sim20-60$ km s$^{-1}$, which indicates that most of the outflowing mass is gravitationally bound. The outflows' column densities correlate with the LMC's star formation rate surface densities ($\Sigma_{\rm SFR}$), and the outflows with higher $\Sigma_{\rm SFR}$ tend to be more ionized. Considering outflows from both sides of the LMC as traced by CIV, we conservatively estimate a total outflow rate of $\dot{M}_{\rm out}\gtrsim 0.03~M_\odot {\rm yr}^{-1}$ and a mass loading factor of $\eta\gtrsim 0.15$. We compare the LMC's outflows with those detected in starburst galaxies and simulation predictions, and find a universal scaling relation of $|v_{\rm out, bulk}|\propto \Sigma_{\rm SFR}^{0.23}$ over a wide range of star-forming conditions ($\Sigma_{\rm SFR}\sim10^{-4.5}-10^{2}~M_\odot {\rm yr}^{-1} {\rm kpc}^{-2}$). Lastly, we find that the outflows are co-rotating with the LMC's young stellar disk and the velocity field does not seem to be significantly impacted by external forces; we thus speculate on the existence of a bow shock leading the LMC, which may have shielded the outflows from ram pressure as the LMC orbits the Milky Way., Comment: Accepted for publication in ApJ. Key findings can be found in Figures 9-12. Main updates from the previous version include new estimates on star formation rate surface densities. Normalized SII, SiIV, and CIV line spectra derived for this work are published as a High Level Science Product called LMC-FLOWS (doi: 10.17909/hz0m-np43), available on website: https://archive.stsci.edu/hlsp/lmc-flows

Protostellar jets and outflows are signposts of active star formation. In H II regions, molecular tracers like CO only reveal embedded portions of the outflow. Outside the natal cloud, outflows are dissociated, ionized, and eventually completely ablated, leaving behind only the high-density jet core. Before this process is complete, there should be a phase where the outflow is partially molecular and partially ionized. In this paper, we capture the HH 900 outflow while this process is in action. New observations from the ERIS/SPIFFIER near-IR integral field unit (IFU) spectrograph using the K-middle filter ($\lambda$=2.06-2.34 $\mu$m) reveal H$_2$ emission from the dissociating outflow and Br-$\gamma$ tracing its ionized skin. Both lines trace the wide-angle outflow morphology but H$_2$ only extends $\sim$5000 au into the H II region while Br-$\gamma$ extends the full length of the outflow ($\sim$12,650 au), indicating rapid dissociation of the molecules. H$_2$ has higher velocities further from the driving source, consistent with a jet-driven outflow. Diagnostic line ratios indicate that photoexcitation, not just shocks, contributes to the excitation in the outflow. We argue that HH 900 is the first clear example of an evaporating molecular outflow and predict that a large column of neutral material that may be detectable with ALMA accompanies the dissociating molecules. Results from this study will help guide the interpretation of near-IR images of externally irradiated jets and outflows such as those obtained with the James Webb Space Telescope (JWST) in high-mass star-forming regions where these conditions may be common., Comment: MNRAS, accepted

We present the first results of the eXtreme UV Environments (XUE) James Webb Space Telescope (JWST) program, that focuses on the characterization of planet forming disks in massive star forming regions. These regions are likely representative of the environment in which most planetary systems formed. Understanding the impact of environment on planet formation is critical in order to gain insights into the diversity of the observed exoplanet populations. XUE targets 15 disks in three areas of NGC 6357, which hosts numerous massive OB stars, among which some of the most massive stars in our Galaxy. Thanks to JWST we can, for the first time, study the effect of external irradiation on the inner ($< 10$ au), terrestrial-planet forming regions of proto-planetary disks. In this study, we report on the detection of abundant water, CO, CO$_2$, HCN and C$_2$H$_2$ in the inner few au of XUE 1, a highly irradiated disk in NGC 6357. In addition, small, partially crystalline silicate dust is present at the disk surface. The derived column densities, the oxygen-dominated gas-phase chemistry, and the presence of silicate dust are surprisingly similar to those found in inner disks located in nearby, relatively isolated low-mass star-forming regions. Our findings imply that the inner regions of highly irradiated disks can retain similar physical and chemical conditions as disks in low-mass star-forming regions, thus broadening the range of environments with similar conditions for inner disk rocky planet formation to the most extreme star-forming regions in our Galaxy., Comment: Accepted for publication in ApJ Letters. 20 pages, 7 figures

Massive star-forming regions are thought to be the most common birth environments in the Galaxy and the only birth places of very massive stars. Their presence in the stellar cluster alters the conditions within the cluster impacting at the same time the evolution of other cluster members. In principle, copious amounts of ultraviolet radiation produced by massive stars can remove material from outer parts of the protoplanetary disks around low- and intermediate-mass stars in the process of external photoevaporation, effectively reducing the planet-formation capabilities of those disks. Here, we present deep VLT/MUSE observations of low-mass stars in Trumpler 14, one of the most massive, young, and compact clusters in the Carina Nebula Complex. We provide spectral and stellar properties of 717 sources and based on the distribution of stellar ages derive the cluster age of $\sim$1~Myr. The majority of the stars in our sample have masses $\leqslant$1~$M_\odot$, what makes our spectroscopic catalogue the most deep to date in term of masses, and proves that detailed investigations of low-mass stars are possible in the massive but distant regions. Spectroscopic studies of low-mass members of the whole Carina Nebula Complex are missing. Our work provides an important step forward towards filling this gap and set the stage for follow-up investigation of accretion properties in Trumpler 14., Comment: Accepted for publication in A&A, 27 pages, 28 figures

We present new observations of the Mystic Mountains cloud complex in the Carina Nebula using the ALMA Atacama Compact Array (ACA) to quantify the impact of strong UV radiation on the structure and kinematics of the gas. Our Band~6 observations target CO, $^{13}$CO, and C$^{18}$O; we also detect DCN J=3-2 and $^{13}$CS J=5-4. A dendrogram analysis reveals that the Mystic Mountains are a coherent structure, with continuous emission over $-$10.5 km s$^{-1}$ $<$ v < $-$2 km s$^{-1}$. We perform multiple analyses to isolate non-thermal motions in the Mystic Mountains including computing the turbulent driving parameter, $b$, which indicates whether compressive or solenoidal modes dominate. Each analysis yields values similar to other pillars in Carina that have been observed in a similar way but are subject to an order of magnitude less intense ionizing radiation. We find no clear correlation between the velocity or turbulent structure of the gas and the incident radiation, in contrast to other studies targeting different regions of Carina. This may reflect differences in the initial densities of regions that go on to collapse into pillars and those that still look like clouds or walls in the present day. Pre-existing over-densities that enable pillar formation may also explain why star formation in the pillars appears more evolved (from the presence of jets) than in other heavily-irradiated but non-pillar-like regions. High resolution observations of regions subject to an array of incident radiation are required to test this hypothesis., Comment: MNRAS, accepted

The ionising radiation of young and massive stars is a crucial form of stellar feedback. Most ionising (Lyman-continuum; LyC, $\lambda < 912A$) photons are absorbed close to the stars that produce them, forming compact HII regions, but some escape into the wider galaxy. Quantifying the fraction of LyC photons that escape is an open problem. In this work, we present a semi-novel method to estimate the escape fraction by combining broadband photometry of star clusters from the Legacy ExtraGalactic UV Survey (LEGUS) with HII regions observed by the Star formation, Ionized gas, and Nebular Abundances Legacy Survey (SIGNALS) in the nearby spiral galaxy NGC 628. We first assess the completeness of the combined catalogue, and find that 49\% of HII regions lack corresponding star clusters as a result of a difference in the sensitivities of the LEGUS and SIGNALS surveys. For HII regions that do have matching clusters, we infer the escape fraction from the difference between the ionising power required to produce the observed HII luminosity and the predicted ionising photon output of their host star clusters; the latter is computed using a combination of LEGUS photometric observations and a stochastic stellar population synthesis code SLUG (Stochastically Lighting Up Galaxies). Overall, we find an escape fraction of $f_{esc} = 0.09^{+0.06}_{-0.06}$ across our sample of 42 HII regions; in particular, we find HII regions with high $f_{esc}$ are predominantly regions with low H$\alpha$-luminosity. We also report possible correlation between $f_{esc}$ and the emission lines [O ii]/[N ii] and [O ii]/H$\beta$., Comment: Accepted for publication at MNRAS

Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity (Z). The Hubble Space Telescope has devoted 500 orbits to observe 250 massive stars at low Z in the ultraviolet (UV) with the COS and STIS spectrographs under the ULLYSES program. The complementary ``X-Shooting ULLYSES'' (XShootU) project provides enhanced legacy value with high-quality optical and near-infrared spectra obtained with the wide-wavelength coverage X-shooter spectrograph at ESO's Very Large Telescope. We present an overview of the XShootU project, showing that combining ULLYSES UV and XShootU optical spectra is critical for the uniform determination of stellar parameters such as effective temperature, surface gravity, luminosity, and abundances, as well as wind properties such as mass-loss rates in function of Z. As uncertainties in stellar and wind parameters percolate into many adjacent areas of Astrophysics, the data and modelling of the XShootU project is expected to be a game-changer for our physical understanding of massive stars at low Z. To be able to confidently interpret James Webb Space Telescope (JWST) spectra of the first stellar generations, the individual spectra of low Z stars need to be understood, which is exactly where XShootU can deliver., Comment: Accepted in A&A - 35 Pages, 12 Figures, 4 Tables, 2 Large Tables

In this work we present the highest spatial and spectral resolution integral field observations to date of the bipolar jet from the Orion proplyd 244-440 using MUSE NFM) observations on the VLT. We observed a previously unreported chain of six distinct knots in a roughly S-shaped pattern, and by comparing them with HST images we estimated proper motions in the redshifted knots of 9.5 mas yr$^{-1}$ with an inclination angle of $73^{\circ}$, though these quantities could not be measured for the blueshifted lobe. Analysis of the [FeII] and [NiII] lines suggests jet densities on the order of $\sim 10^5$ cm$^{-3}$. We propose that the observed S-shaped morphology originates from a jet launched by a smaller source with $M_\star < 0.2$ M$_{\odot}$ in orbital motion around a larger companion of $M_\star \simeq 0.5$ M$_{\odot}$ at a separation of 30-40 au. The measured luminosities of the knots using the [OI]$\lambda6300$ and [SII]$\lambda6731$ lines were used to estimate a lower limit to the mass-loss rate in the jet of $1.3 \times 10^{-11}$ M$_{\odot}$ yr$^{-1}$ and an upper limit of $10^{-9}$ M$_{\odot}$ yr$^{-1}$, which is typical for low-mass driving sources. While the brightness asymmetry between the redshifted and blueshifted lobes is consistent with external irradiation, further analysis of the [NiII] and [FeII] lines suggests that photoionization of the jet is not likely to be a dominant factor, and that the emission is dominated by collisional excitation. The dynamical age of the jet compared to the anticipated survival time of the proplyd demonstrates that photoevaporation of the proplyd occurred prior to jet launching, and that this is still an active source. These two points suggest that the envelope of the proplyd may shield the jet from the majority of external radiation, and that photoionization of the proplyd does not appear to impact the ability of a star to launch a jet., Comment: 30 pages, 20 figures

Stars strongly impact their environment, and shape structures on all scales throughout the universe, in a process known as ``feedback''. Due to the complexity of both stellar evolution and the physics of larger astrophysical structures, there remain many unanswered questions about how feedback operates, and what we can learn about stars by studying their imprint on the wider universe. In this white paper, we summarize discussions from the Lorentz Center meeting `Bringing Stellar Evolution and Feedback Together' in April 2022, and identify key areas where further dialogue can bring about radical changes in how we view the relationship between stars and the universe they live in., Comment: Accepted to the Publications of the Astronomical Society of the Pacific

Binary neutron star mergers are one of the ultimate events of massive binary star evolution, and our understanding of their parent system is still in its infancy. Upcoming gravitational wave detections, coupled with multi-wavelength follow-up observations, will allow us to study an increasing number of these events by characterising their neighbouring stellar populations and searching for their progenitors. Stellar evolution simulations are essential to this work but they are also based on numerous assumptions. Additionally, the models used to study the host galaxies differ from those used to characterise the progenitors and are typically based on single star populations. Here we introduce a framework to perform an end-to-end analysis and deploy it to the first binary neutron star merger - GW170817. With the Binary Population And Spectral Synthesis (BPASS) codes we are able to retrieve the physical properties of the host galaxy NGC 4993 as well as infer progenitor candidates. In our simulations there is a >98% chance that GW170817 originated from a stellar population with Z=0.010 born between 5 and 12.5 Gyrs ago. By carefully weighing the stellar genealogies we find that GW170817 most likely came from a binary system born with a 13-24 Msol primary and 10-12 Msol secondary which underwent two or three common envelope events over their lifetime., Comment: Main text (3 Figures, 10 pages) and Supplementary Information (23 pages, 14 Figures). Published in Nature Astronomy

We characterize massive stars (M>8 M_sun) in the nearby (D~0.8 Mpc) extremely metal-poor (Z~5% Z_sun) galaxy Leo A using Hubble Space Telescope ultra-violet (UV), optical, and near-infrared (NIR) imaging along with Keck/LRIS and MMT/Binospec optical spectroscopy for 18 main sequence OB stars. We find that: (a) 12 of our 18 stars show emission lines, despite not being associated with an H II region, suggestive of stellar activity (e.g., mass loss, accretion, binary star interaction), which is consistent with previous predictions of enhanced activity at low metallicity; (b) 6 are Be stars, which are the first to be spectroscopically studied at such low metallicity -- these Be stars have unusual panchromatic SEDs; (c) for stars well-fit by the TLUSTY non-local thermodynamic equilibrium (non-LTE) models, the photometric and spectroscopic values of T_eff and log(g) agree to within ~0.01 dex and ~0.18 dex, respectively, indicating that NUV/optical/NIR imaging can be used to reliably characterize massive (M ~ 8-30 M_sun) main sequence star properties relative to optical spectroscopy; (d) the properties of the most massive stars in H II regions are consistent with constraints from previous nebular emission line studies; and (e) 13 stars with M>8 M_sun are >40 pc from a known star cluster or H II region. Our sample comprises ~50% of all known massive stars at Z < 10% Z_sun with derived stellar parameters, high-quality optical spectra, and panchromatic photometry., Comment: 35 pages, 18 figures

Baryonic feedback is expected to play a key role in regulating the star formation of low-mass galaxies by producing galaxy-scale winds associated with mass-loading factors $\beta\!\sim\!1\!-\!50$. We have tested this prediction using a sample of 19 nearby systems with stellar masses $10^7\!<\!M_\star/{\rm M}_{\odot}\!<\!10^{10}$, mostly lying above the main sequence of star-forming galaxies. We used MUSE@VLT optical integral field spectroscopy to study the warm ionised gas kinematics of these galaxies via a detailed modelling of their H$\alpha$ emission line. The ionised gas is characterised by irregular velocity fields, indicating the presence of non-circular motions of a few tens of km/s within galaxy discs, but with intrinsic velocity dispersion of $40$-$60$ km/s that are only marginally larger than those measured in main-sequence galaxies. Galactic winds, defined as gas at velocities larger than the galaxy escape speed, encompass only a few percent of the observed fluxes. Mass outflow rates and loading factors are strongly dependent on $M_\star$, star formation rate (SFR), SFR surface density and specific SFR. For $M_\star$ of $10^8$ M$_\odot$ we find $\beta\simeq0.02$, which is more than two orders of magnitude smaller than the values predicted by theoretical models of galaxy evolution. In our galaxy sample, baryonic feedback stimulates a gentle gas cycle rather than causing a large-scale blow out., Comment: 26 pages, 12 figures, accepted by A&A

We study pre-supernova feedback in a sample of $\sim$ 4700 HII regions in the nearby spiral galaxy M83, identified on their H$\alpha$ emission. We pectroscopically identify Wolf-Rayet (WR) stars populating the star-forming regions. For each HII region, we compute the pressure of ionised gas ($P_{\rm ion}$) and the direct radiation pressure ($P_{\rm dir}$) acting in the region, and investigate how they vary with galactocentric distance, with the physical properties of the region, and with the pressure of the galactic environment ($P_\mathrm{DE}$). For a subset of $\sim$ 500 regions, we also investigate the link between the pressure terms and the properties of the cluster population (age, mass, and LyC flux). We find that $P_{\rm ion}$ dominates over $P_{\rm dir}$ by at least a factor of 10 on average over the disk. Both pressure terms are strongly enhanced and become almost comparable in the central starburst region. In the disk ($R \geq 0.15~R_e$), we observe that $P_{\rm dir}$ stays approximately constant with galactocentric distance. $P_{\rm dir}$ is positively correlated with an increase in radiation field strength (linked to the negative metallicity gradient in the galaxy), while it decreases in low extinction regions. $P_{\rm ion}$ decreases constantly for increasing galactocentric distances. In general, we observe that HII regions near the center are underpressured with respect to their surroundings, whereas regions in the disk are overpressured and hence expanding. We find that regions hosting younger clusters or having more mass in young star clusters have a higher internal pressure, indicating that clustered star formation is likely playing a dominant role in setting the pressure. Finally, we estimate that only 13 % of HII regions hosting young clusters and WR stars have $f_{\rm esc} \geq 0$.[Abridged], Comment: Accepted for publication in A&A

One of the fundamental factors regulating the evolution of galaxies is stellar feedback. However, we still do not have strong observational constraints on the relative importance of the different feedback mechanisms (e.g. radiation, ionised gas pressure, stellar winds) in driving HII region evolution and molecular cloud disruption. In this letter, we constrain the relative importance of the various feedback mechanisms from young massive star populations by resolving HII regions across the disk of the nearby star-forming galaxy NGC 1672. We combine measurements of ionised gas nebular lines obtained by PHANGS-MUSE, with high-resolution imaging from the HST in both the narrow-band H{\alpha} and broad-band filters. We identify a sample of 40 isolated, compact HII regions in the HST H{\alpha} image, for which we measure the sizes that were previously unresolved in seeing-limited ground-based observations. Additionally, we identify the ionisation source(s) for each HII region from catalogues produced as part of the PHANGS-HST survey. We find that the HII regions investigated are mildly dominated by thermal or wind pressure, yet their elevation above the radiation pressure is within the expected uncertainty range. We see that radiation pressure provides a substantially higher contribution to the total pressure than previously found in the literature over similar size scales. In general, we find higher pressures within more compact HII regions, which is driven by the inherent size scaling relations of each pressure term, albeit with significant scatter introduced by the variation in the stellar population properties (e.g. luminosity, mass, age, metallicity). For nearby galaxies, here we provide a promising approach that could yield the statistics required to map out how the importance of different stellar feedback mechanisms evolve over the lifetime of an HII region., Comment: 5(+4) pages, 6 figures, 1 table (available online), accepted for publication in A&A Letters

Giant molecular clouds (GMCs) are the sites of star formation and stellar feedback in galaxies. Their properties set the initial conditions for star formation and their lifecycles determine how feedback regulates galaxy evolution. In recent years, the advent of high-resolution telescopes has enabled systematic GMC-scale studies of the molecular interstellar medium in nearby galaxies, now covering a wide range of physical conditions and allowing for the first studies of how GMC properties depend on galactic environment. These observational developments have been accompanied by numerical simulations of improving resolution that are increasingly accurately accounting for the effects of the galactic-scale environment on GMCs, while simultaneously improving the treatment of the small-scale processes of star-formation and stellar feedback within them. The combination of these recent developments has greatly improved our understanding of the formation, evolution, and destruction of GMCs. We review the current state of the field, highlight current open questions, and discuss promising avenues for future studies., Comment: 37 pages, 5 figures; To appear in Protostars and Planets VII; Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

Here we describe a new study of the SNRs and SNR candidates in nearby face-on spiral galaxy M83, based primarily on MUSE integral field spectroscopy. Our revised catalog of SNR candidates in M83 has 366 objects, 81 of which are reported here for the first time. Of these, 229 lie within the MUSE observation region, 160 of which have spectra with [SII]:Halpha ratios exceeding 0.4, the value generally accepted as confirmation that an emission nebula is shock-heated. Combined with 51 SNR candidates outside the MUSE region with high [SII]:Halpha ratios, there are 211 spectroscopically-confirmed SNRs in M83, the largest number of confirmed SNRs in any external galaxy. MUSE's combination of relatively high spectral resolution and broad wavelength coverage has allowed us to explore two other properties of SNRs that could serve as the basis of future SNR searches. Specifically, most of the objects identified as SNRs on the basis of [SII]:Halpha ratios exhibit more velocity broadening and lower ratios of [SIII]:[SII] emission than HII regions. A search for nebulae with the very broad emission lines expected from young, rapidly expanding remnants revealed none, except for the previously identified B12-174a. The SNRs identified in M83 are, with few exceptions, middle-aged ISM-dominated ones. Smaller diameter candidates show a larger range of velocity broadening and a larger range of gas densities than the larger diameter objects, as expected if the SNRs expanding into denser gas brighten and then fade from view at smaller diameters than those expanding into a more tenuous ISM, Comment: 28 page, including 14 figures and 4 tables; accepted for publication in ApJ

We present a large VLT/MUSE mosaic (3.8 x 3.8 kpc) of the nearby spiral galaxy M83, with a spatial resolution ~20 pc. We obtained the kinematics of the stars and ionised gas, and compared them with molecular gas kinematics from ALMA CO(2-1). We separated the ionised gas into HII regions and diffuse ionised gas (DIG) and determined the fraction of Ha luminosity originating from the DIG (f_DIG). We observe that both stars and gas trace the galactic disk rotation, as well as a fast-rotating nuclear component, likely connected to secular processes driven by the galactic bar. In the gas kinematics, we observe a stream east of the nucleus, redshifted with respect to the disk. The stream is surrounded by an extended ionised gas region with enhanced velocity dispersion and a high ionisation state, which is largely consistent with being ionised by slow shocks. We interpret this feature as either the superposition of the disk and an extraplanar layer of DIG, or as a bar-driven inflow of shocked gas. A double Gaussian component fit to the Ha line also reveals the presence of a nuclear biconic structure whose axis of symmetry is perpendicular to the bar. The two cones appear blue- and redshifted along the line of sight and stand out for having an Ha emission separated by up to 200 km s-1 from that of the disk, and a high velocity dispersion ~80-200 km s-1. At the far end of the cones, we observe that the gas is consistent with being ionised by shocks. These features had never been observed before in M83; we postulate that they are tracing a starburst-driven outflow shocking into the surrounding ISM. Finally, we obtain f_DIG ~ 13% in our field of view. We inspect the emission of the HII regions and DIG in `BPT' diagrams, finding that in HII regions photoionisation accounts for 99.8% of the Ha flux, whereas the DIG has a mixed contribution from photoionisation (94.9%) and shocks (5.1%). [abridged], Comment: Accepted for publication in A&A

The radio continuum emission is a dust-unbiased tracer of both the thermal and non-thermal processes in the interstellar medium. We present new maps of the free-free and synchrotron emission in the Magellanic Clouds (MCs) at 0.166, 1.4, and 4.8GHz with no prior assumption about the radio non-thermal spectrum. The maps were constructed using a de-reddened H$\alpha$ map as a template for the thermal radio emission, which we subtract from maps of the total radio continuum emission. To de-redden the H$\alpha$ emission, it is important to know the fraction of dust surface density that attenuates the H$\alpha$ emission along the line-of-sight, $f_\text{d}$. This fraction is obtained by comparing the dust opacity obtained through the infrared emission spectrum and the Balmer decrement method. In star-forming regions, the median $f_\text{d}$ is about 0.1 which is by a factor of 3 lower than that in diffuse regions. We obtain a global thermal radio fraction, $f_\text{th}$, of 30 per cent (35 per cent) in the LMC (SMC) at 1.4GHz. Furthermore, we present maps of the equipartition magnetic field strength with average values of $\simeq10.1\mu$G in the LMC and $\simeq\,5.5\,\mu$G in the SMC. The magnetic field is proportional to the star formation rate to a power of 0.24 and 0.20 for the LMC and SMC, respectively. This study shows that the non-thermal processes control the interstellar medium in the MCs., Comment: Accepted for publication in MNRAS

Integral field units enable resolved studies of a large number of star-forming regions across entire nearby galaxies, providing insight on the conversion of gas into stars and the feedback from the emerging stellar populations over unprecedented dynamic ranges in terms of spatial scale, star-forming region properties, and environments. We use the VLT/MUSE legacy data set covering the central $35$ arcmin$^{2}$ (${\sim}12$ kpc$^{2}$) of the nearby galaxy NGC 300 to quantify the effect of stellar feedback as a function of the local galactic environment. We extract spectra from emission line regions identified within dendrograms, combine emission line ratios and line widths to distinguish between HII regions, planetary nebulae, and supernova remnants, and compute their ionised gas properties, gas-phase oxygen abundances, and feedback-related pressure terms. For the HII regions, we find that the direct radiation pressure ($P_\mathrm{dir}$) and the pressure of the ionised gas ($P_{HII}$) weakly increase towards larger galactocentric radii, i.e. along the galaxy's (negative) abundance and (positive) extinction gradients. While the increase of $P_{HII}$ with galactocentric radius is likely due to higher photon fluxes from lower-metallicity stellar populations, we find that the increase of $P_\mathrm{dir}$ is likely driven by the combination of higher photon fluxes and enhanced dust content at larger galactocentric radii. In light of the above, we investigate the effect of increased pre-supernova feedback at larger galactocentric distances (lower metallicities and increased dust mass surface density) on the ISM, finding that supernovae at lower metallicities expand into lower-density environments, thereby enhancing the impact of supernova feedback., Comment: Accepted for publication in MNRAS

Analyzing the dynamical state of nearby young massive star clusters is essential understanding star cluster formation and evolution during their earliest stages. In this work we analyze the stellar and gas kinematics of the young massive star cluster Westerlund 2 (Wd2) using data from the integral field unit MUSE and complement them with proper motions from the Gaia DR2. The mean gas radial velocity of $15.9\,{\rm km}\,{\rm s}^{-1}$ agrees with the assumption that Wd2 is the result of a cloud-cloud collision. The gas motions show the expansion of the HII region, driven by the radiation from the many OB stars in the cluster center. The velocity profile of the cluster member stars reveal an increasing velocity dispersion with decreasing stellar mass and that the low-mass stars show five distinct velocity groups. Based on their spatial correlation with the cluster's two clumps, we concluded that this is the imprint of the initial cloud collapse that formed Wd2. A thorough analysis of the dynamical state of Wd2, which determines a dynamical mass range of $M_{\rm dyn,Wd2} = (7.5 \pm 1.9)\cdot10^4 - (4.4 \pm 1.1)\cdot10^5\,{\rm M}_\odot$ and exceeds the photometric mass by at least a factor of two leads to the conclusion that Wd2 is not massive enough to remain gravitationally bound. Additionally we also identify 22 runaway candidates with peculiar velocities between 30 and $546\,{\rm km}\,{\rm s}^{-1}$., Comment: 23 pages, 14 figures, 4 tables, accepted to AJ

We combine MUSE and ALMA observations with theoretical models to evaluate how a tadpole-shaped globule located in the Carina Nebula has been influenced by its environment. This globule is now relatively small (radius ~2500 au), hosts a protostellar jet+outflow (HH 900) and, with a blue-shifted velocity of ~10 km/s, is travelling faster than it should be if its kinematics were set by the turbulent velocity dispersion of the precursor cloud. Its outer layers are currently still subject to heating, but comparing the internal and external pressures implies that the globule is in a post-collapse phase. Intriguingly the outflow is bent, implying that the YSO responsible for launching it is comoving with the globule, which requires that the star formed after the globule was up to speed since otherwise it would have been left behind. We conclude that the most likely scenario is one in which the cloud was much larger before being subject to radiatively-driven implosion, which accelerated the globule to the high observed speeds under the photoevaporative rocket effect and triggered the formation of the star responsible for the outflow. The globule may now be in a quasi-steady state following collapse. Finally, the HH 900 YSO is likely $\gtrsim$1 M$_{\odot}$ and may be the only star forming in the globule. It may be that this process of triggered star formation has prevented the globule from fragmenting to form multiple stars (e.g., due to heating) and has produced a single higher mass star., Comment: MNRAS, accepted

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the tadpole, a small globule in the Carina Nebula that hosts the HH 900 jet+outflow system. Our data include $^{12}$CO, $^{13}$CO, C$^{18}$O J=2-1, $^{13}$CO, C$^{18}$O J=3-2, and serendipitous detections of DCN J=3-2 and CS J=7-6. With angular resolution comparable to the Hubble Space Telescope (HST), our data reveal for the first time the bipolar molecular outflow in CO, seen only inside the globule, that is launched from the previously unseen jet-driving protostar (the HH 900 YSO). The biconical morphology joins smoothly with the externally irradiated outflow seen in ionized gas tracers outside the globule, tracing the overall morphology of a jet-driven molecular outflow. Continuum emission at the location of the HH 900 YSO appears to be slightly flattened perpendicular to outflow axis. Model fits to the continuum have a best-fit spectral index of $\sim 2$, suggesting cold dust and the onset of grain growth. In position-velocity space, $^{13}$CO and C$^{18}$O gas kinematics trace a C-shaped morphology, similar to infall profiles seen in other sources, although the global dynamical behaviour of the gas remains unclear. Line profiles of the CO isotopologues display features consistent with externally heated gas. We estimate a globule mass of $\sim 1.9$ M$_{\odot}$, indicating a remaining lifetime of $\sim 4$ Myr, assuming a constant photoevaporation rate. This long globule lifetime will shield the disk from external irradiation perhaps prolonging its life and enabling planet formation in regions where disks are typically rapidly destroyed., Comment: MNRAS, accepted

We present the very first wide-field, 11 by 11 arcmin, optical spectral mapping of M16, one of the most famous star-forming regions in the Galaxy. The data were acquired with the new imaging Fourier transform spectrograph SITELLE mounted on the Canada-France-Hawaii Telescope (CFHT). We obtained three spectral cubes with R=10'000 (SN1 filter), 1500 (SN2 filter) and 600 (SN3 filter), centered on the Pillars of Creation and the HH216 flow, covering the main optical nebular emission lines: [O II]3726,29 (SN1), Hb, [O III]4959,5007 (SN2), [N II]6548,84, Ha, and [S II]6717,31 (SN3). We validate the performance, calibration, and data reduction of SITELLE, and analyze the structures in the large field-of-view in terms of their kinematics and nebular emission. We compared the SITELLE data to MUSE integral field observations and other spectroscopic and narrow-band imaging data to validate the performance of SITELLE. We computed gas-phase metallicities via the strong-line method, performed a pixel-by-pixel fit to the main emission lines to derive kinematics of the ionized gas, computed the mass-loss rate of the Eastern pillar (the Spire), and combined the SITELLE data with near-infrared narrow-band imaging to characterize the HH216 flow. The comparison with previously published fluxes demonstrates very good agreement. We disentangle the dependence of the gas-phase metallicities (derived via abundance-tracing line ratios) on the degree of ionization and obtain metallicities that are in excellent agreement with the literature. We confirm the bipolar structure of HH216, find evidence for episodic accretion from the source of the flow, and identify its likely driving source. We compute the mass-loss rate of the Spire pillar on the East side of the H II region and find excellent agreement with the correlation between the mass-loss rate and the ionizing photon flux from the nearby cluster NGC6611., Comment: 21 pages, accepted for publication in A&A

The INT Galactic Plane Survey (IGAPS) is the merger of the optical photometric surveys, IPHAS and UVEX, based on data from the Isaac Newton Telescope (INT) obtained between 2003 and 2018. Here, we present the IGAPS point source catalogue. It contains 295.4 million rows providing photometry in the filters, i, r, narrow-band Halpha, g and U_RGO. The IGAPS footprint fills the Galactic coordinate range, |b| < 5deg and 30deg < l < 215deg. A uniform calibration, referred to the Pan-STARRS system, is applied to g, r and i, while the Halpha calibration is linked to r and then is reconciled via field overlaps. The astrometry in all 5 bands has been recalculated on the Gaia DR2 frame. Down to i ~ 20 mag (Vega system), most stars are also detected in g, r and Halpha. As exposures in the r band were obtained within the IPHAS and UVEX surveys a few years apart, typically, the catalogue includes two distinct r measures, r_I and r_U. The r 10sigma limiting magnitude is ~21, with median seeing 1.1 arcsec. Between ~13th and ~19th magnitudes in all bands, the photometry is internally reproducible to within 0.02 magnitudes. Stars brighter than r=19.5 have been tested for narrow-band Halpha excess signalling line emission, and for variation exceeding |r_I-r_U| = 0.2 mag. We find and flag 8292 candidate emission line stars and over 53000 variables (both at >5sigma confidence). The 174-column catalogue will be available via CDS Strasbourg., Comment: 28 pages, 22 figures

We present MUSE Integral Field Unit (IFU) observations of five individual HII regions in two giant (>100 pc in radius) star-forming complexes in the low-metallicity ($Z$~0.33 $Z_{\odot}$) nearby (D ~ 2 Mpc) dwarf spiral galaxy NGC 300. We combine the IFU data with high spatial resolution HST photometry to demonstrate the extraction of stellar spectra and the classification of individual stars from ground-based data at the distance of 2 Mpc. For the two star-forming complexes, in which no O-type stars had previously been identified, we find a total of 13 newly identified O-type stars in the mass range 15-50 M$_{\odot}$, as well as 4 Wolf-Rayet stars. We use the derived massive stellar content to analyze the impact of stellar feedback on the HII regions. As already found for HII regions in the Magellanic Clouds, the dynamics of the analyzed NGC 300 HII regions are dominated by a combination of the pressure of the ionized gas and stellar winds. By comparing the derived ionized gas mass loading factors to the total gas mass loading factor across the NGC 300 disk, we find that the latter is an order of magnitude higher, either indicating very early evolutionary stages for these HII regions, or being a direct result of the multi-phase nature of feedback-driven bubbles. Moreover, we analyze the relation between the star formation rate and the pressure of the ionized gas as derived from small (<100 pc) scales, as both quantities are systematically overestimated when derived on galactic scales. With the wealth of ongoing and upcoming IFU instruments and programs, this study serves as a pathfinder for the systematic investigation of resolved stellar feedback in nearby galaxies, and it delivers the necessary analysis tools to enable massive stellar content and feedback studies sampling an unprecedented range of HII region properties across entire galaxies in the nearby Universe., Comment: Accepted for publication in ApJ

Forming high-mass stars have a significant effect on their natal environment. Their feedback pathways, including winds, outflows, and ionising radiation, shape the evolution of their surroundings which impacts the formation of the next generation of stars. They create or reveal dense pillars of gas and dust towards the edges of the cavities they clear. They are modelled in feedback simulations, and the sizes and shapes of the pillars produced are consistent with those observed. However, these models predict measurably different kinematics which provides testable discriminants. Here we present the first ALMA Compact Array (ACA) survey of 13 pillars in Carina, observed in $^{12}$CO, $^{13}$CO and C$^{18}$O J=2-1, and the 230 GHz continuum. The pillars in this survey were chosen to cover a wide range in properties relating to the amount and direction of incident radiation, proximity to nearby irradiating clusters and cloud rims, and whether they are detached from the cloud. With these data, we are able to discriminate between models. We generally find pillar velocity dispersions of $<$ 1 km s$^{-1}$ and that the outer few layers of molecular emission in these pillars show no significant offsets from each other, suggesting little bulk internal motions within the pillars. There are instances where the pillars are offset in velocity from their parental cloud rim, and some with no offset, hinting at a stochastic development of these motions., Comment: 24 Pages, 19 figures. Accepted to MNRAS

We present new MUSE/VLT observations of a small globule in the Carina H II region that hosts the HH 900 jet+outflow system. Data were obtained with the GALACSI ground-layer adaptive optics system in wide-field mode, providing spatially-resolved maps of diagnostic emission lines. These allow us to measure the variation of the physical properties in the globule and jet+outflow system. We find high temperatures ($T_e \approx 10^4$ K), modest extinction ($A_V \approx 2.5$ mag), and modest electron densities ($n_e \approx 200$ cm$^{-3}$) in the ionized gas. Higher excitation lines trace the ionized outflow; both the excitation and ionization in the outflow increase with distance from the opaque globule. In contrast, lower excitation lines that are collisionally de-excited at densities $\gtrsim 10^4$ cm$^{-3}$ trace the highly collimated protostellar jet. Assuming the globule is an isothermal sphere confined by the pressure of the ionization front, we compute a Bonnor-Ebert mass of $\sim 3.7$ M$_{\odot}$. This is two orders of magnitude higher than previous mass estimates, calling into question whether small globules like the Tadpole contribute to the bottom of the IMF. Derived globule properties are consistent with a cloud that has been and/or will be compressed by the ionization front on its surface. At the estimated globule photoevaporation rate of $\sim 5 \times 10^{-7}$ M$_{\odot}$ yr$^{-1}$, the globule will be completely ablated in $\sim 7$ Myr. Stars that form in globules like the Tadpole will emerge into the H II later and may help resolve some of the temporal tension between disk survival and enrichment., Comment: MNRAS, accepted

We mapped the Galactic young massive star cluster Westerlund 2 (Wd2) with the integral field spectrograph MUSE (spatial resolution: 0.2arcsec/px, spectral resolution: $\Delta \lambda$ = 1.25A, wavelength range 4600-9350A) mounted on the VLT, as part of an on-going study to measure the stellar and gas kinematics of the cluster region. In this paper we present the fully reduced dataset and introduce our new Python package "MUSEpack", which we developed to measure stellar radial velocities with an absolute precision of 1-2km/s without the necessity of a spectral template library. This novel method uses the two-dimensional spectra and an atomic transition line library to create templates around strong absorption lines for each individual star. The code runs fully automatically on multi-core machines, which makes it possible to efficiently determine stellar radial velocities of a large number of stars with the necessary precision to measure the velocity dispersion of young star clusters. MUSEpack also provides an enhanced method for removing telluric lines in crowded fields without sky exposures and a Python wrapper for ESO's data reduction pipeline. We observed Wd2 with a total of 11 short and 5 long exposures to cover the bright nebular emission and OB stars, as well as the fainter pre-main sequence stars down to ~1M$_\odot$. The survey covers an area of ~11arcmin$^2$ (15.8pc$^2$). In total, we extracted 1,725 stellar spectra with a mean S/N>5 per pixel. A typical radial velocity (RV) uncertainty of 4.78km/s, 2.92km/s, and 1.1km/s is reached for stars with a mean S/N>10, S/N>20, S/N>50 per pixel, respectively. Depending on the number of spectral lines used to measure the RVs, it is possible to reach RV accuracies of 0.9km/s, 1.3km/s, and 2.2km/s with $\geq5$, 3-4, and 1-2 spectral lines, respectively. The combined statistical uncertainty on the radial velocity measurements is 1.10km/s., Comment: accepted to AJ; 19 pages, 10 figures, 4 tables; MUSEpack is available from https://github.com/pzeidler89/MUSEpack; The documentation can be found on: https://musepack.readthedocs.io/en/latest/index.html

SIGNALS, the Star formation, Ionized Gas, and Nebular Abundances Legacy Survey, is a large observing program designed to investigate massive star formation and HII regions in a sample of local extended galaxies. The program will use the imaging Fourier transform spectrograph SITELLE at the Canada-France-Hawaii Telescope. Over 355 hours (54.7 nights) have been allocated beginning in fall 2018 for eight consecutive semesters. Once completed, SIGNALS will provide a statistically reliable laboratory to investigate massive star formation, including over 50 000 resolved HII regions : the largest, most complete, and homogeneous database of spectroscopically and spatially resolved extragalactic HII regions ever assembled. For each field observed, three datacubes covering the spectral bands of the filters SN1 (363 -386 nm), SN2 (482 - 513 nm), and SN3 (647 - 685 nm) are gathered. The spectral resolution selected for each spectral band is 1000, 1000, and 5000, respectively. As defined, the project sample will facilitate the study of small-scale nebular physics and many other phenomena linked to star formation at a mean spatial resolution of 20 pc. This survey also has considerable legacy value for additional topics including planetary nebulae, diffuse ionized gas, andsupernova remnants. The purpose of this paper is to present a general outlook of the survey, notably the observing strategy, galaxy sample, and science requirements., Comment: 19 pages, 14 figures, submitted to MNRAS

The physics of star formation and the deposition of mass, momentum, and energy into the interstellar medium by massive stars (`feedback') are the main uncertainties in modern cosmological simulations of galaxy formation and evolution. These processes determine the properties of galaxies, but are poorly understood on the $\lesssim$100 pc scale of individual giant molecular clouds (GMCs) resolved in modern galaxy formation simulations. The key question is why the timescale for depleting molecular gas through star formation in galaxies ($t_{\rm dep}\approx 2$ Gyr) exceeds the dynamical timescale of GMCs by two orders of magnitude. Either most of a GMC's mass is converted into stars over many dynamical times, or only a small fraction turns into stars before the GMC is dispersed on a dynamical timescale. Here we report our observation that molecular gas and star formation are spatially de-correlated on GMC scales in the nearby flocculent spiral galaxy NGC300, contrary to their tight correlation on galactic scales. We demonstrate that this de-correlation implies rapid evolutionary cycling between GMCs, star formation, and feedback. We apply a novel statistical method to quantify the evolutionary timeline and find that star formation is regulated by efficient stellar feedback, driving GMC dispersal on short timescales (~1.5 Myr) due to radiation and stellar winds, prior to supernova explosions. This feedback limits GMC lifetimes to about one dynamical timescale (~10 Myr), with integrated star formation efficiencies of only 2-3%. Our findings reveal that galaxies consist of building blocks undergoing vigorous, feedback-driven lifecycles, that vary with the galactic environment and collectively define how galaxies form stars. Systematic applications of this multi-scale analysis to large galaxy samples will provide key input for a predictive, bottom-up theory of galaxy formation and evolution., Comment: Published in the May 23 issue of Nature; this is the authors' version before final edits; low-resolution versions of the Supplementary Videos are available as ancillary files in the arXiv source (see the link in the top right of this page), whereas the full-resolution versions are available with the Nature publication

We recently reported the discovery of a candidate jet-driving microquasar (S10) in the nearby spiral galaxy NGC 300. However, in the absence of kinematic information, we could not reliably determine the jet power or the dynamical age of the jet cavity. Here, we present optical MUSE integral field unit (IFU) observations of S10, which reveal a bipolar line-emitting jet structure surrounding a continuum-emitting central source. The optical jet lobes of S10 have a total extent of $\sim$ 40 pc and a shock velocity of $\sim$ 150 km s$^{-1}$. Together with the jet kinematics, we exploit the MUSE coverage of the Balmer H$\beta$ line to estimate the density of the surrounding matter and therefore compute the jet power to be $P_{jet}\approx$ 6.3 $\times$ 10$^{38}$ erg s$^{-1}$. An optical analysis of a microquasar jet bubble and a consequent robust derivation of the jet power have been possible only in a handful of similar sources. This study therefore adds valuable insight into microquasar jets, and demonstrates the power of optical integral field spectroscopy in identifying and analysing these objects., Comment: Accepted for publication in MNRAS

We present the discovery of a powerful candidate microquasar in NGC 300, associated with the S 10 optical nebula (previously classified as a supernova remnant). Chandra images show four discrete X-ray knots aligned in the plane of the sky over a length of $\approx$150 pc. The X-ray emission from the knots is well fitted with a thermal plasma model at a temperature of $\approx$0.6 keV and a combined 0.3-8 keV luminosity of $\approx$10$^{37}$ erg s$^{-1}$. The X-ray core, if present at all, does not stand out above the thermal emission of the knots: this suggests that the accreting compact object is either currently in a dim state or occulted from our view. We interpret the emission from the knots as the result of shocks from the interaction of a jet with the interstellar medium (possibly over different epochs of enhanced activity). Cooler shock-heated gas is likely the origin also of the optical bubble and lobes near the X-ray structure, detected in images from the Hubble Space Telescope and the Very Large Telescope. In the radio bands, we observed the region with the Australia Telescope Compact Array, and discovered an elongated radio nebula (about 170 $\times$ 55 pc in size) with its major axis aligned with the chain of Chandra sources. The radio nebula has an integrated 5.5 GHz radio luminosity of $\approx$10$^{34}$ erg s$^{-1}$ for a distance of 1.88 Mpc. The morphology, size and luminosity of the extended X-ray, optical and radio structure suggest that NGC 300-S 10 belongs to the same class of powerful ($P_{\rm jet} > 10^{39}$ erg s$^{-1}$) microquasars as SS 433, Ho II X-1 and NGC 7793-S 26., Comment: 19 pages, 8 figures, 8 tables. Accepted for publication in MNRAS

We present MUSE integral field data of two HII region complexes in the Large Magellanic Cloud (LMC), N44 and N180. Both regions consist of a main superbubble and a number of smaller, more compact HII regions that formed on the edge of the superbubble. For a total of 11 HII regions, we systematically analyse the radiative and mechanical feedback from the massive O-type stars on the surrounding gas. We exploit the integral field property of the data and the coverage of the HeII$\lambda$5412 line to identify and classify the feedback-driving massive stars, and from the estimated spectral types and luminosity classes we determine the stellar radiative output in terms of the ionising photon flux $Q_{0}$. We characterise the HII regions in terms of their sizes, morphologies, ionisation structure, luminosity and kinematics, and derive oxygen abundances via emission line ratios. We analyse the role of different stellar feedback mechanisms for each region by measuring the direct radiation pressure, the pressure of the ionised gas, and the pressure of the shock-heated winds. We find that stellar winds and ionised gas are the main drivers of HII region expansion in our sample, while the direct radiation pressure is up to three orders of magnitude lower than the other terms. We relate the total pressure to the star formation rate per unit area, $\Sigma_{SFR}$, for each region and find that stellar feedback has a negative effect on star formation, and sets an upper limit to $\Sigma_{SFR}$ as a function of increasing pressure., Comment: Accepted for publication in MNRAS, 27 pages, 21 figures

Highly collimated parsec-scale jets, generally linked to the presence of an accretion disk, are a commonly observed phenomenon from revealed low-mass young stellar objects. In the past two decades, only a very few of these objects have been directly (or indirectly) observed towards high-mass (M > 8 M$_{\odot}$) young stellar objects, adding to the growing evidence that disk-mediated accretion is a phenomenon also occurring in high-mass stars, the formation mechanism of which is still poorly understood. Of the observed jets from massive young stars, none is in the optical regime (due to these being typically highly obscured by their native material), and none are found outside of the Milky Way. Here, we report the detection of HH 1177, the first extragalactic optical ionized jet originating from a massive young stellar object located in the Large Magellanic Cloud. The jet is highly collimated over the entire measured extent of at least 10 pc, and has a bipolar geometry. The presence of a jet indicates ongoing, disk-mediated accretion, and together with the high degree of collimation, this system is therefore likely to be an up-scaled version of low-mass star formation. We conclude that the physics governing jet launching and collimation is independent of stellar mass., Comment: 9 pages, 5 figures, 2 tables

We present an analysis of the effect of feedback from O- and B-type stars with data from the integral field spectrograph MUSE mounted on the Very Large Telescope of pillar-like structures in the Carina Nebular Complex, one of the most massive star-forming regions in the Galaxy. For the observed pillars, we compute gas electron densities and temperatures maps, produce integrated line and velocity maps of the ionised gas, study the ionisation fronts at the pillar tips, analyse the properties of the single regions, and detect two ionised jets originating from two distinct pillar tips. For each pillar tip we determine the incident ionising photon flux $Q_\mathrm{0,pil}$ originating from the nearby massive O- and B-type stars and compute the mass-loss rate $\dot{M}$ of the pillar tips due to photo-evaporation caused by the incident ionising radiation. We combine the results of the Carina data set with archival MUSE data of a pillar in NGC 3603 and with previously published MUSE data of the Pillars of Creation in M16, and with a total of 10 analysed pillars, find tight correlations between the ionising photon flux and the electron density, the electron density and the distance from the ionising sources, and the ionising photon flux and the mass-loss rate. The combined MUSE data sets of pillars in regions with different physical conditions and stellar content therefore yield an empirical quantification of the feedback effects of ionising radiation. In agreement with models, we find that $\dot{M}\propto Q_\mathrm{0,pil}^{1/2}$., Comment: Accepted for publication in MNRAS, 36 pages, 30 figures