Search

We combine the deepest X-ray survey from the Chandra Deep Field-South (CDF-S) `7-Ms' survey with the deepest mid-infrared (5.6$ \mu m$) image from the JWST/MIRI Deep Imaging Survey (MIDIS) in the Hubble Ultra-Deep Field (HUDF) to study the infrared counterparts and point-source emission of 31 X-ray sources with a median, intrinsic, rest-frame X-ray luminosity of $\log_{10}(L_{\rm Xc}^{\rm 0.5-7keV})$=42.04$\pm$0.22 erg $\rm s^{-1}$. The sample includes 24 AGN with a redshift range, as set by the X-ray detectability, of $z \simeq 0.5-3$. Through a multi-wavelength morphological decomposition, employing three separate classifications (visual, parametric and non-parametric) we separate (where present) the luminosity of the point-like AGN component from the remainder of the host-galaxy emission. The unprecedented mid-infrared sensitivity and imaging resolution of MIRI allows, in many cases, the direct characterisation of point-like (i.e. unresolved) components in the galaxies' emission. We establish a broad agreement between the three morphological classifications. At least 70% of the X-ray sources, including some classified as galaxies, show unresolved emission in the MIRI images, with the unresolved-to-total flux fraction at rest-frame 2$\mu m$ ranging from $\sim$0.2 to $\sim$0.9. At high X-ray luminosities ($\log_{10}(L_{\rm Xc}$)>43 erg $\rm s^{-1}$) we derive a consistent rest-frame near-infrared 2$ \mu m$ point-source luminosity to that derived for local AGN, whilst at lower X-ray luminosity we identify an excess in the 2$ \mu m$ emission compared to pre-JWST studies. We speculate this offset may be driven by a combination of Compton-thick AGN components and nuclear starburst, merger driven activity. Our observations highlight the complex nature of X-ray sources in the distant Universe and demonstrate the power of JWST/MIRI in quantifying their nuclear infrared emission. (Abridged), Comment: 21 pages, 7 figures. Submitted to A&A

As part of the guaranteed time observations program Mid-Infrared Characterization Of Nearby Iconic galaxy Centers (MICONIC), we used the medium-resolution spectrometer (MRS) of the Mid-Infrared Instrument (MIRI) on board of the JWST to study the nearby merger NGC6240. We aim to characterise the dual active galactic nuclei (AGN), the ionised gas outflows and the main properties of the interstellar medium over a mapped area of 6.6"x7.7". We obtained integral field spectroscopic mid-infrared data of NGC6240, resolving both nuclei for the first time in the full 5-28{\mu}m spectral range. We modelled the emission lines through a kinematic decomposition, finding that the fine-structure lines in the southern (S) nucleus are broader than for the northern (N) nucleus (full width at half maximum of $\geq$1500 vs ~700 km s^{-1} on average). High excitation lines, such as [NeV], [NeVI], and [MgV], are clearly detected in the N nucleus. In the S nucleus, the same lines can be detected but only after a decomposition of the PAH features in the integrated spectrum, due to a combination of a strong mid-IR continuum, the broad emission lines, and the intense star formation (SF). The SF is distributed all over the mapped FoV of 3.5kpc x 4.1kpc (projected), with the maximum located around the S nucleus. Both nuclear regions appear to be connected by a bridge region detected with all the emission lines. Based on the observed MRS line ratios and the high velocity dispersion ({\sigma}~600 km s^{-1}), shocks are also dominating the emission in this system. We detected the presence of outflows as a bubble north-west from the N nucleus and at the S nucleus. We estimated a ionised mass outflow rate of 1.4$\pm$0.3 M yr^{-1} and 1.8$\pm$0.2 M yr^{-1}, respectively. Given the derived kinetic power of these outflows, both the AGN and the starburst could have triggered them. [Abridged], Comment: 23 pages, 10 figures. Accepted for publication in A&A on 26/11/2024

This paper presents a deep MIRI/JWST medium resolution spectroscopy (MRS) covering the rest-frame optical spectrum of the GN-z11 galaxy. The [OIII]5008 and H$\alpha$ emission lines are detected and spectroscopically resolved. The line profiles are well-modeled by a narrow Gaussian component with intrinsic FWHMs of 189$\pm$25 and 231$\pm$52 kms$^{-1}$, respectively. We do not find any evidence of a dominant broad H$\alpha$ emission line component tracing a Broad Line Region in a type 1 active galactic nuclei (AGN). However, a broad ($\sim$430-470 kms$^{-1}$) and weak ($<$ 20-30%) H$\alpha$ line component, tracing a minor AGN contribution in the optical, cannot be ruled out completely with the sensitivity of the present data. The physical and excitation properties of the ionized gas are consistent with a low-metallicity starburst forming stars at a rate of SFR(H$\alpha$)$=$24 $\pm$3$M_{\odot}$yr$^{-1}$. The electron temperature of the ionized gas is $T_{\mathrm{e}}$(O$^{++}$)$=$14000$\pm$2100K, while the direct-$T_{\mathrm{e}}$ gas-phase metallicity is 12+$\log$(O/H)$=$7.91$\pm$0.07 (Z=0.17$\pm$0.03Z$_{\odot}$). The optical line ratios locate GN-z11 in the starburst or AGN region but more consistent with those of local low-metallicity starbursts and high-$z$ luminous galaxies detected at redshifts similar to GN-z11. We conclude that the MRS optical spectrum of GN-z11 is consistent with that of a massive, compact, and low-metallicity starburst galaxy. Due to its high SFR and stellar mass surface densities, close to that of the densest stellar clusters, we speculate that GN-z11 could be undergoing a feedback-free, highly efficient starburst phase. Additional JWST data are needed to validate this scenario, and other recently proposed alternatives, to explain the existence of bright compact galaxies in the early Universe., Comment: accepted for publication in A&A (17 pages, 6 figures, 2 tables)

The Lyman-$\alpha$ (Lya) line of hydrogen is a well-known tracer of galaxies at high-z. However, the connection between Lya observables and galaxy properties has not fully been established, limiting the use of the line to probe the physics of galaxies. Here, we derive global neutral hydrogen gas (HI) properties of nearby Lya-emitting galaxies to assess the impact of HI on the Lya output of galaxies. We observed 21cm line emission using the VLA in D-array configuration (~55" resolution, ~38 kpc) for 37 star-forming galaxies with available Lya imaging from the Lyman Alpha Reference Samples (LARS and eLARS). We detect 21cm emission for 33/37 galaxies observed. We find no significant correlation of global HI properties with Lya luminosity, escape fraction or equivalent width derived with HST photometry. Additionally, both Lya-emitters and weak or non-emitters are distributed evenly along the HI parameter space of optically-selected z=0 galaxies. Around 74% of the sample is undergoing galaxy interaction, this fraction is higher for Lya-emitters (83% for galaxies with EW$\geq$20\r{A}) than for non or weak emitters (70%). Nevertheless, galaxies identified as interacting have Lya and HI properties statistically consistent with those of non-interacting galaxies. Our results show that global HI properties (on scales > 30kpc) have little direct impact on the Lya output from galaxies. Instead, HI likely regulates Lya emission on small scales: statistical comparisons of Lya and high angular resolution 21cm observations are required to fully assess the role of HI in Lya radiative transfer. While our study indicates that galaxy mergers could play a role in the emission of Lya photons in the local universe, especially for galaxies with high HI fractions, the line-of-sight through which a system is observed ultimately determines Lya observables., Comment: 15+35 pages, 5+33 figures, Accepted for publication in A&A

With a temperature of $\sim 285$ K WISE0855 is the coldest brown dwarf observed so far. Using the James Webb Space Telescope (JWST) we obtained observations that allow us to characterize WISE0855s atmosphere focusing on vertical variation in the water steam abundance, measuring trace gas abundances and receiving bulk parameters for this cold object. We observed the ultra cool dwarf WISE0855 using the Mid-Infrared Instrument Medium Resolution Spectrometer (MIRI/MRS) onboard JWST at a spectral resolution of up to 3750. We combined the observation with published data from the Near Infrared Spectrograph (NIRSpec) G395M and PRISM modes yielding a spectrum ranging from 0.8 to 22 um. We apply atmospheric retrievals using petitRADTRANS to measure atmospheric abundances, the pressure-temperature structure, radius and gravity of the brown dwarf. We also employ publicly available clear and cloudy self-consistent grid models to estimate bulk properties of the atmosphere such as the effective temperature, radius, gravity and metallicity. Atmospheric retrievals constrain a variable water abundance profile in the atmosphere, as predicted by equilibrium chemistry. We detect the 15NH3 isotopologue and infer a ratio of mass fraction of 14NH3/15NH3 = 332+63-43 for the clear retrieval. We measure the bolometric luminosity by integrating the presented spectrum and obtain a value of log(L/L$_{\odot}$) = -7.291+/-0.008. The detected water depletion indicates that water condenses out in the upper atmosphere due to the very low effective temperature of WISE0855. The height in the atmosphere where this occurs is covered by the MIRI/MRS data, and thus demonstrates the potential of MIRI to characterize cold gas giants atmospheres. Comparing the data to retrievals and self-consistent grid models, we do not detect signs for water ice clouds, although their spectral features have been predicted in previous studies., Comment: Submitted to A&A, 29 pages, 21 figures

We use the James Webb Space Telescope (JWST) and its Mid-Infrared Instrument (MIRI) (5-28 um), to study the embedded HH 211 flow. We map a 0.95'x0.22' region, covering the full extent of the blue-shifted lobe, the central protostellar region, and a small portion of the red-shifted lobe. The jet driving source is not detected even at the longest mid-IR wavelengths. The overall morphology of the flow consists of a highly collimated jet, mostly molecular (H2, HD) with an inner atomic ([FeI], [FeII], [SI], [NiII]) structure. The jet shocks the ambient medium, producing several large bow-shocks, rich in forbidden atomic and molecular lines, and is driving an H2 molecular outflow, mostly traced by low-J, v=0 transitions. Moreover, 0-0 S(1) uncollimated emission is also detected down to 2"-3" (~650-1000 au) from the source, tracing a cold (T=200-400 K), less dense and poorly collimated molecular wind. The atomic jet ([FeII] at 26 um) is detected down to ~130 au from source, whereas the lack of H2 emission close to the source is likely due to the large visual extinction. Dust continuum-emission is detected at the terminal bow-shocks, and in the blue- and red-shifted jet, being likely dust lifted from the disk. The jet shows an onion-like structure, with layers of different size, velocity, temperature, and chemical composition. Moreover, moving from the inner jet to the outer bow-shocks, different physical, kinematic and excitation conditions for both molecular and atomic gas are observed. The jet mass-flux rate, momentum, and momentum flux of the warm H2 component are up to one order of magnitude higher than those inferred from the atomic jet component. Our findings indicate that the warm H2 component is the primary mover of the outflow, namely it is the most significant dynamical component of the jet, in contrast to jets from more evolved YSOs, where the atomic component is dominant., Comment: Paper accepted in A&A for publication

We present JWST/MIRI MRS spatially resolved $\sim 5-28\,\mu$m observations of the central ~4-8kpc of the ultraluminous infrared galaxy and broad absorption line quasar Mrk231. These are part of the Mid-Infrared Characterization of Nearby Iconic galaxy Centers (MICONIC) program of the MIRI European Consortium guaranteed time observations. No high excitation lines (i.e., [MgV] at 5.61$\mu$m or [NeV] at 14.32$\mu$m) typically associated with the presence of an active galactic nucleus (AGN) are detected in the nuclear region of Mrk231. This is likely due to the intrinsically X-ray weak nature of its quasar. Some intermediate ionization potential lines, for instance, [ArIII] at 8.99$\mu$m and [SIV] at 10.51$\mu$m, are not detected either, even though they are clearly observed in a star-forming region ~920pc south-east of the AGN. Thus, the strong nuclear mid-infrared (mid-IR) continuum is also in part hampering the detection of faint lines in the nuclear region. The nuclear [NeIII]/[NeII]line ratio is consistent with values observed in star-forming galaxies. Moreover, we resolve for the first time the nuclear starburst in the mid-IR low-excitation line emission (size of ~400pc, FWHM). Several pieces of evidence also indicate that it is partly obscured even at these wavelengths. At the AGN position, the ionized and warm molecular gas emission lines have modest widths (W_80~300km/s). There are, however, weak blueshifted wings reaching velocities v_02~-400km/s in [NeII]. The nuclear starburst is at the center of a large (~8kpc), massive rotating disk with widely-spread, low velocity outflows. Given the high star formation rate of Mrk231, we speculate that part of the nuclear outflows and the large-scale non-circular motions observed in the mid-IR are driven by its powerful nuclear starburst., Comment: Accepted for publication in Astronomy & Astrophysics

Thanks to decades of observations using HST, the structure of galaxies at redshift $z>2$ has been widely studied in the rest-frame ultraviolet regime, which traces recent star formation from young stellar populations. But, we still have little information about the spatial distribution of the older, more evolved, stellar populations, constrained by the rest-frame infrared portion of galaxies' spectral energy distribution. We present the morphological characterization of a sample of 21 massive galaxies ($\log(M_{\star}/M_{\odot})>9.5$) at redshift $33$ massive galaxies show a smooth distribution of their rest-infrared light, strongly supporting the increasing number of regular disk galaxies already in place at early epochs. On the contrary, the ultraviolet structure obtained from HST observations is generally more irregular, catching the most recent episodes of star formation. Importantly, we find a segregation of morphologies across cosmic time, having massive galaxies at redshift $z>4$ later-type morphologies compared to $z\sim3$ galaxies. These findings suggest a transition phase in galaxy assembly and central mass build up already taking place at $z\sim3-4$. MIRI provides unique information about the structure of the mature stellar population of high-redshift galaxies, unveiling that massive galaxies beyond cosmic noon are prevalently compact disk galaxies with smooth mass distribution., Comment: Submitted to A&A. 7 pages, 4 figures. Comments welcome

Lyman break analogs in the local Universe serve as counterparts to Lyman break galaxies (LBGs) at high redshifts, which are widely regarded as major contributors to cosmic reionization in the early stages of the Universe. We studied XMM-Newton and Chandra observations of the nearby LBG analog Haro 11, which contains two X-ray-bright sources, X1 and X2. Both sources exhibit Lyman continuum (LyC) leakage, particularly X2. We analyzed the X-ray variability using principal component analysis (PCA) and performed spectral modeling of the X1 and X2 observations made with the Chandra ACIS-S instrument. The PCA component, which contributes to the X-ray variability, is apparently associated with variable emission features, likely from ionized superwinds. Our spectral analysis of the Chandra data indicates that the fainter X-ray source, X2 (X-ray luminosity $L_{\rm X} \sim 4 \times 10^{40} $ erg s$^{-1}$), the one with higher LyC leakage, has a much lower absorbing column ($N_{\rm H} \sim 1.2 \times 10^{21}$ cm$^{-2}$) than the heavily absorbed luminous source X1 ($L_{\rm X} \sim 9 \times 10^{40} $ erg s$^{-1}$ and $N_{\rm H} \sim 11.5 \times 10^{21}$ cm$^{-2}$). We conclude that X2 is likely less covered by absorbing material, which may be a result of powerful superwinds clearing galactic channels and facilitating the escape of LyC radiation. Much deeper X-ray observations are required to validate the presence of potential superwinds and determine their implications for the LyC escape., Comment: 15 pages, 8 figures, 4 tables, accepted by A&A

Very low-mass stars (those <0.3 solar masses) host orbiting terrestrial planets more frequently than other types of stars, but the compositions of those planets are largely unknown. We use mid-infrared spectroscopy with the James Webb Space Telescope to investigate the chemical composition of the planet-forming disk around ISO-ChaI 147, a 0.11 solar-mass star. The inner disk has a carbon-rich chemistry: we identify emission from 13 carbon-bearing molecules including ethane and benzene. We derive large column densities of hydrocarbons indicating that we probe deep into the disk. The high carbon to oxygen ratio we infer indicates radial transport of material within the disk, which we predict would affect the bulk composition of any planets forming in the disk., Comment: Published, 36 pages, 8 figures

Gas kinematics affect the radiative transfer and escape of hydrogen Lyman-$\alpha$ (Ly$\alpha$) emission from galaxies. We investigate this interplay empirically by relating the ionised gas kinematics of 42 galaxies in the extended Ly$\alpha$ Reference Sample with their Ly$\alpha$ escape fractions, $f_\rm{esc}$, Ly$\alpha$ equivalent widths, $\rm{EW}_\rm{Ly\alpha}$, and Ly$\alpha$ luminosities, $L_\rm{Ly\alpha}$. To this aim we use PMAS integral-field spectroscopic observations of the Balmer-$\alpha$ line. We calculate shearing velocities, $v_\rm{shear}$, and intrinsic velocity dispersions, $\sigma_0^\rm{obs}$ (empirically corrected for beam-smearing effects), as global kinematical measures for each galaxy. The sample is characterised by highly turbulent motions and more than half of the sample shows dispersion dominated kinematics. We uncover clear trends between Ly$\alpha$ observables and global kinematical statistics. We discuss statistically the importance of $v_\rm{shear}$, $\sigma_0^\rm{obs}$, and $v_\rm{shear}/\sigma_0^\rm{obs}$ for regulating the Ly$\alpha$ observables in comparison to other galaxy parameters. It emerges that $\sigma_0^\rm{obs}$ is the dominating parameter for $\rm{EW}_\rm{Ly\alpha}$ and that is as important as nebular extinction, gas covering fraction, and ionising photon production efficiency in regulating $f_\rm{esc}$. A simple scenario where the starburst age is simultaneously regulating turbulence, $\rm{EW}_\rm{Ly\alpha}$, and $f_\rm{esc}$ does not find support by our observations. However, we show that the small scale distribution of dust appears to be influenced by turbulence in some galaxies. In support of our observational result we discuss how turbulence is theoretically expected to play a significant role in modulating $f_\rm{esc}$. (abridged), Comment: Accepted version for publication in A&A

Context. Two protoplanets have recently been discovered within the PDS 70 protoplanetary disk. JWST/NIRCam offers a unique opportunity to characterize them and their birth environment at wavelengths difficult to access from the ground. Aims. We aim to image the circumstellar environment of PDS 70 at 1.87 $\mu$m and 4.83 $\mu$m, assess the presence of Pa-$\alpha$ emission due to accretion onto the protoplanets, and probe any IR excess indicative of heated circumplanetary material. Methods. We obtain non-coronagraphic JWST/NIRCam images of PDS 70 within the MINDS (MIRI mid-INfrared Disk Survey) program. We leverage the Vortex Image Processing (VIP) package for data reduction, and develop dedicated routines for optimal stellar PSF subtraction, unbiased imaging of the disk, and protoplanet flux measurement in this type of dataset. A radiative transfer model of the disk is used to disentangle the contributions from the disk and the protoplanets. Results. We re-detect both protoplanets and identify extended emission after subtracting a disk model, including a large-scale spiral-like feature. We interpret its signal in the direct vicinity of planet c as tracing the accretion stream feeding its circumplanetary disk, while the outer part of the feature may rather reflect asymmetric illumination of the outer disk. We also report a bright signal consistent with a previously proposed protoplanet candidate enshrouded in dust, near the 1:2:4 mean-motion resonance with planets b and c. The 1.87 $\mu$m flux of planet b is consistent with atmospheric model predictions, but not that of planet c. We discuss potential origins for this discrepancy, including significant Pa-$\alpha$ line emission. The 4.83 $\mu$m fluxes of planets b and c suggest enshrouding dust or heated CO emission from their circumplanetary environment., Comment: 6+11 pages, 3+10 figures (text+appendix). Accepted for publication in A&A v1: Version accepted for publication in A&A. v2: Tables 1 and C.1 were updated to reflect an erratum on the reported astrometry and photometry of protoplanet candidate d. The original version did not include the results from our best reduction, which corrects for a bias associated with self-subtraction

Despite the importance of the dusty star-forming galaxies (DSFGs) at $z$>2 for understanding the galaxy evolution in the early Universe, their stellar distributions traced by the near-IR emission were spatially unresolved until the arrival of the JWST. In this work we present, for the first time, a spatially-resolved morphological analysis of the rest-frame near-IR (~1.1-3.5$\mu$m) emission in DSFGs traced with the JWST/MIRI. In particular, we study the mature stellar component for the three DSFGs and a Lyman-break galaxy (LBG) present in an overdensity at $z$=4.05. Moreover, we use MIRI images along with UV to (sub)-mm ancillary photometric data to model their SEDs and extract their main physical properties. The sub-arcsec resolution MIRI images have revealed that the stellar component present a wide range of morphologies, from disc-like to compact and clump-dominated structures. These near-IR structures contrast with their UV emission, which is usually diffuse and off-centered. The SED fitting analysis shows that GN20 dominates the total SFR with a value ~2500 $M_\odot$yr$^{-1}$ while GN20.2b has the highest stellar mass in the sample ($M_*$~2$\times$10$^{11}$ $M_\odot$). The two DSFGs classified as LTGs (GN20 and GN20.2a) show high specific SFR (sSFR>30 Gyr$^{-1}$) placing them above the star-forming main sequence (SFMS) at z~4 by >0.5 dex while the ETG (i.e.,GN20.2b) is compatible with the high-mass end of the main sequence. When comparing with other DSFGs in overdensities at $z$~2-7 we observe that our objects present similar SFRs, depletion times and projected separations. Nevertheless, the effective radii computed for GN20 and GN20.2a are up to two times larger than those of isolated galaxies observed in CEERS and ALMA-HUDF at similar redshifts. We interpret this difference in size as an effect of rapid growth induced by the dense environment., Comment: 13 pages, 6 figures, 3 tables. Accepted for publication in A&A

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

We report observational evidence of highly turbulent ionized gas kinematics in a sample of 20 Lyman continuum (LyC) emitters (LCEs) at low redshift ($z\sim 0.3$). Detailed Gaussian modeling of optical emission line profiles in high-dispersion spectra consistently shows that both bright recombination and collisionally excited lines can be fitted as one or two narrow components with intrinsic velocity dispersion of $\sigma$ $\sim$ 40-100 km s$^{-1}$, in addition to a broader component with $\sigma \sim$ 100-300 km s$^{-1}$, which contributes up to $\sim$40% of the total flux and is preferentially blueshifted from the systemic velocity. We interpret the narrow emission as highly ionized gas close to the young massive star clusters and the broader emission as a signpost of unresolved ionized outflows, resulting from massive stars and supernova feedback. We find a significant correlation between the width of the broad emission and the LyC escape fraction, with strong LCEs exhibiting more complex and broader line profiles than galaxies with weaker or undetected LyC emission. We provide new observational evidence supporting predictions from models and simulations; our findings suggest that gas turbulence and outflows resulting from strong radiative and mechanical feedback play a key role in clearing channels through which LyC photons escape from galaxies. We propose that the detection of blueshifted broad emission in the nebular lines of compact extreme emission-line galaxies can provide a new indirect diagnostic of Lyman photon escape, which could be useful to identify potential LyC leakers in the epoch of reionization with the JWST., Comment: 11 Pages, 7 Figures. A&A in Press

Dusty star-forming galaxies emit most of their light at far-IR to mm wavelengths as their star formation is highly obscured. Far-IR and mm observations have revealed their dust, neutral and molecular gas properties. The sensitivity of JWST at rest-frame optical and near-infrared wavelengths now allows the study of the stellar and ionized gas content. We investigate the spatially resolved distribution and kinematics of the ionized gas in GN20, a dusty star forming galaxy at $z$=4.0548. We present deep MIRI/MRS integral field spectroscopy of the near-infrared rest-frame emission of GN20. We detect spatially resolved \paa, out to a radius of 6 kpc, distributed in a clumpy morphology. The star formation rate derived from \paa\ (144 $\pm$ 9 \msunperyear) is only 7.7 $\pm 0.5 $\% of the infrared star formation rate (1860 $\pm$ 90 \msunperyear). We attribute this to very high extinction (A$_V$ = 17.2 $\pm$ 0.4 mag, or A$_{V,mixed}$ = 44 $\pm$ 3 mag), especially in the nucleus of GN20, where only faint \paa\ is detected, suggesting a deeply buried starburst. We identify four, spatially unresolved, clumps in the \paa\ emission. Based on the double peaked \paa\ profile we find that each clump consist of at least two sub-clumps. We find mass upper limits consistent with them being formed in a gravitationally unstable gaseous disk. The UV bright region of GN20 does not have any detected \paa\ emission, suggesting an age of more than 10 Myrs for this region of the galaxy. From the rotation profile of \paa\ we conclude that the gas kinematics are rotationally dominated and the $v_{rot}/\sigma_{m} = 3.8 \pm 1.4$ is similar to low-redshift LIRGs. We speculate that the clumps seen in GN20 could contribute to building up the inner disk and bulge of GN20., Comment: 13 pages, 7 figures, accepted for publication in A&A

We study the resolved stellar populations and derive the star formation history of NGC 5474, a peculiar star-forming dwarf galaxy at a distance of $\sim 7$ Mpc, using Hubble Space Telescope Advanced Camera for Surveys data from the Legacy Extragalactic UV Survey (LEGUS) program. We apply an improved colour-magnitude diagram fitting technique based on the code SFERA and use the latest PARSEC-COLIBRI stellar models. Our results are the following. The off-centre bulge-like structure, suggested to constitute the bulge of the galaxy, is dominated by star formation (SF) activity initiated $14$ Gyr ago and lasted at least up to $1$ Gyr ago. Nevertheless, this component shows clear evidence of prolonged SF activity (lasting until $\sim 10$ Myr ago). We estimate the total stellar mass of the bulge-like structure to be $(5.0 \pm 0.3) \times 10^{8}$ \MSUN. Such a mass is consistent with published suggestions that this structure is in fact an independent system orbiting around and not within NGC 5474's disc. The stellar over-density located to the South-West of the bulge-like structure shows a significant SF event older than $1$ Gyr, while it is characterised by two recent peaks of SF, around $\sim10$ and $\sim100$ Myr ago. In the last Gyr, the behavior of the stellar disc is consistent with what is known in the literature as `gasping'. The synchronised burst at $10-35$ Myr in all components might hint to the recent gravitational interaction between the stellar bulge-like structure and the disc of NGC 5474., Comment: 18 pages, 12 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

(abridged) We study and compare the warm (>100 K) and cold (<100 K) material toward the high-mass star-forming region IRAS 23385+6053 (IRAS 23385 hereafter) combining high angular resolution observations in the mid-infrared (MIR) with the JWST Observations of Young protoStars (JOYS) project and with the NOEMA at mm wavelengths at angular resolutions of 0.2"-1". The spatial morphology of atomic and molecular species is investigated by line integrated intensity maps. The temperature and column density of different gas components is estimated using H2 transitions (warm and hot component) and a series of CH3CN transitions as well as 3 mm continuum emission (cold component). Toward the central dense core in IRAS 23385 the material consists of relatively cold gas and dust (~50 K), while multiple outflows create heated and/or shocked H2 and show enhanced temperatures (~400 K) along the outflow structures. An energetic outflow with enhanced emission knots of [Fe II] and [Ni II] hints at J-type shocks, while two other outflows have enhanced emission of only H2 and [S I] caused by C-type shocks. The latter two outflows are also more prominent in molecular line emission at mm wavelengths (e.g., SiO, SO, H2CO, and CH3OH). Even higher angular resolution data are needed to unambiguously identify the outflow driving sources given the clustered nature of IRAS 23385. While most of the forbidden fine structure transitions are blueshifted, [Ne II] and [Ne III] peak at the source velocity toward the MIR source A/mmA2 suggesting that the emission is originating from closer to the protostar., Comment: 15 pages, 7 figures, accepted for publication in A&A

We present MIRI/JWST medium-resolution spectroscopy (MRS) and imaging (MIRIM) of the lensed galaxy MACS1149-JD1 at a redshift of $z$=9.1092$\pm$0.0002 (Universe age about 530 Myr). We detect, for the first time, spatially resolved H$\alpha$ emission in a galaxy at a redshift above nine. The structure of the H$\alpha$ emitting gas consists of two clumps, S and N. The total H$\alpha$ luminosity implies an instantaneous star-formation of 5.3$\pm$0.4 $M_{\odot}$ yr$^{-1}$ for solar metallicities. The ionizing photon production efficiency, $\log(\zeta_\mathrm{ion})$, shows a spatially resolved structure with values of 25.55$\pm$0.03, 25.47$\pm$0.03, and 25.91$\pm$0.09 Hz erg$^{-1}$ for the integrated galaxy, and clumps S and N, respectively. The H$\alpha$ rest-frame equivalent width, EW$_{0}$(H$\alpha$), is 726$^{+660}_{-182}$ \'Angstrom for the integrated galaxy, but presents extreme values of 531$^{+300}_{-96}$ \'Angstrom and $\geq$1951 \'Angstrom for clumps S and N, respectively. The spatially resolved ionizing photon production efficiency is within the range of values measured in galaxies at redshift above six, and well above the canonical value (25.2$\pm$0.1 Hz erg$^{-1}$). The extreme difference of EW$_{0}$(H$\alpha$) for Clumps S and N indicates the presence of a recent (<5 Myrs) burst in clump N and a star formation over a larger period of time (e.g., $\sim$50 Myr) in clump S. Finally, clump S and N show very different H$\alpha$ kinematics with velocity dispersions of 56$\pm$4 km s$^{-1}$ and 113$\pm$33 km s$^{-1}$, likely indicating the presence of outflows or increased turbulence in the clump N. The dynamical mass, $M_\mathrm{dyn}$= (2.4$\pm$0.5)$\times$10$^{9}$ $M_{\odot}$, is within the range measured with the spatially resolved [OIII]88$\mu$m line., Comment: 11 pages, 6 figures, published version in A&A

We report on further observations of the field of the quasar Q1218+0832. Geier et al. 2019 presented the discovery of the quasar resulting from a search for quasars reddened and dimmed by dust in foreground damped Lyman-alpha absorbers (DLAs). The DLA is remarkable by having a very large HI column density close to 10^22 cm^-2 . Its dust extinction curve shows the 2175 AA bump known from the Local Group. It also shows absorption from cold gas exemplified by CI and CO molecules. For this paper, we present narrow-band observations of the field of Q1218+0832 and also use an archival Hubble Space Telescope (HST) image to search for the galaxy counterpart of the DLA. No emission from the DLA galaxy is found in either the narrow-band imaging or in the HST image. In the HST image, we could probe down to an impact parameter of 0.3 arcsec and a 3-sigma detection limit of 26.8 mag per arcsec^2. In the narrow-band image, we probed down to a 0 arcsec impact parameter and detected nothing down to a 3-sigma detection limit of about 3x10-17 erg s^-1 cm^-2 . We did detect a bright Lyman-alpha emitter 59 arcsec south of Q1218+0832 with a flux of 3x10^-16 erg s^-1 cm^-2 . We conclude that the DLA galaxy must be located at a very small impact parameter (<0.3 arcsec, 2.5 kpc) or it is optically dark. Also, the DLA galaxy most likely is part of a galaxy group., Comment: 9 pages, 10 figures. Accepted for publication in A&A. Corrections from the proofs have been added

Terrestrial and sub-Neptune planets are expected to form in the inner ($<10~$AU) regions of protoplanetary disks. Water plays a key role in their formation, although it is yet unclear whether water molecules are formed in-situ or transported from the outer disk. So far Spitzer Space Telescope observations have only provided water luminosity upper limits for dust-depleted inner disks, similar to PDS 70, the first system with direct confirmation of protoplanet presence. Here we report JWST observations of PDS 70, a benchmark target to search for water in a disk hosting a large ($\sim54~$AU) planet-carved gap separating an inner and outer disk. Our findings show water in the inner disk of PDS 70. This implies that potential terrestrial planets forming therein have access to a water reservoir. The column densities of water vapour suggest in-situ formation via a reaction sequence involving O, H$_2$, and/or OH, and survival through water self-shielding. This is also supported by the presence of CO$_2$ emission, another molecule sensitive to UV photodissociation. Dust shielding, and replenishment of both gas and small dust from the outer disk, may also play a role in sustaining the water reservoir. Our observations also reveal a strong variability of the mid-infrared spectral energy distribution, pointing to a change of inner disk geometry., Comment: To appear in Nature on 24 July 2023. 21 pages, 10 figures; includes extended data. Part of the JWST MINDS Guaranteed Time Observations program's science enabling products. Spectra downloadable on Zenodo at https://zenodo.org/record/7991022

Luminous infrared galaxies at high redshifts ($z$>4) include extreme starbursts that build their stellar mass over short periods of time (>100 Myr). These galaxies are considered to be the progenitors of massive quiescent galaxies at intermediate redshifts ($z\sim$2) but their stellar structure and buildup is unknown. Here, we present the first spatially resolved near-infrared imaging of GN20, one of the most luminous dusty star-forming galaxies known to date, observed at an epoch when the Universe was only 1.5 Gyr old. The 5.6$\mu$m image taken with the JWST Mid-Infrared Instrument (MIRI/JWST) shows that GN20 is a very luminous galaxy (M$_\mathrm{1.1\mu m,AB}$=$-$25.01), with a stellar structure composed of a conspicuous central source and an extended envelope. The central source is an unresolved nucleus that carries 9% of the total flux. The nucleus is co-aligned with the peak of the cold dust emission, and offset by 3.9 kpc from the ultraviolet stellar emission. The diffuse stellar envelope is similar in size to the clumpy CO molecular gas distribution. The centroid of the stellar envelope is offset by 1 kpc from the unresolved nucleus, suggesting GN20 is involved in an interaction or merger event supported by its location as the brightest galaxy in a proto-cluster. The stellar size of GN20 is larger by a factor of about 3-5 than known spheroids, disks, and irregulars at $z\sim$4, while its size and low S\'ersic index are similar to those measured in dusty, infrared luminous galaxies at $z\sim$2 of the same mass. GN20 has all the ingredients necessary for evolving into a massive spheroidal quiescent galaxy at intermediate $z$: it is a large, luminous galaxy at $z$=4.05 involved in a short and massive starburst centred in the stellar nucleus and extended over the entire galaxy, out to radii of 4 kpc, and likely induced by the interaction or merger with a member of the proto-cluster., Comment: 7 pages, 4 figures

Carbon is an essential element for life but how much can be delivered to young planets is still an open question. The chemical characterization of planet-forming disks is a crucial step in our understanding of the diversity and habitability of exoplanets. Very low-mass stars ($<0.2~M_{\odot}$) are interesting targets because they host a rich population of terrestrial planets. Here we present the JWST detection of abundant hydrocarbons in the disk of a very low-mass star obtained as part of the MIRI mid-INfrared Disk Survey (MINDS). In addition to very strong and broad emission from C$_2$H$_2$ and its $^{13}$C$^{12}$CH$_2$ isotopologue, C$_4$H$_2$, benzene, and possibly CH$_4$ are identified, but water, PAH and silicate features are weak or absent. The lack of small silicate grains implies that we can look deep down into this disk. These detections testify to an active warm hydrocarbon chemistry with a high C/O ratio in the inner 0.1 au of this disk, perhaps due to destruction of carbonaceous grains. The exceptionally high C$_2$H$_2$/CO$_2$ and C$_2$H$_2$/H$_2$O column density ratios suggest that oxygen is locked up in icy pebbles and planetesimals outside the water iceline. This, in turn, will have significant consequences for the composition of forming exoplanets., Comment: version submitted to Nature Astronomy

Aims: The JWST program JOYS (JWST Observations of Young protoStars) aims at characterizing the physical and chemical properties of young high- and low-mass star-forming regions, in particular the unique mid-infrared diagnostics of the warmer gas and solid-state components. We present early results from the high-mass star formation region IRAS23385+6053. Methods: The JOYS program uses the MIRI MRS with its IFU to investigate a sample of high- and low-mass star-forming protostellar systems. Results: The 5 to 28mum MIRI spectrum of IRAS23385+6053 shows a plethora of features. While the general spectrum is typical for an embedded protostar, we see many atomic and molecular gas lines boosted by the higher spectral resolution and sensitivity compared to previous space missions. Furthermore, ice and dust absorption features are also present. Here, we focus on the continuum emission, outflow tracers like the H2, [FeII] and [NeII] lines as well as the potential accretion tracer Humphreys alpha HI(7--6). The short-wavelength MIRI data resolve two continuum sources A and B, where mid-infrared source A is associated with the main mm continuum peak. The combination of mid-infrared and mm data reveals a young cluster in its making. Combining the mid-infrared outflow tracer H2, [FeII] and [NeII] with mm SiO data shows a complex interplay of at least three molecular outflows driven by protostars in the forming cluster. Furthermore, the Humphreys alpha line is detected at a 3-4sigma level towards the mid-infrared sources A and B. Following Rigliaco et al. (2015), one can roughly estimate accretion luminosities and corresponding accretion rates between ~2.6x10^-6 and ~0.9x10^-4 M_sun/yr. This is discussed in the context of the observed outflow rates. Conclusions: The analysis of the MIRI MRS observations for this young high-mass star-forming region reveals connected outflow and accretion signatures., Comment: 12 pages, 9 figures, accepted for Astronomy & Astrophysics, the paper is also available at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

We make use of \textit{JWST} medium and broad-band NIRCam imaging, along with ultra-deep MIRI $5.6 \rm \mu m$ imaging, in the Hubble eXtreme Deep Field (XDF) to identify prominent line emitters at $z\simeq 7-8$. Out of a total of 58 galaxies at $z\simeq 7-8$, we find 18 robust candidates ($\simeq$31\%) for (H$\beta$ + [OIII]) emitters, based on their enhanced fluxes in the F430M and F444W filters, with EW$_{0}$(H$\beta$ +[OIII]) $\simeq 87 - 2100$ {\AA}. Among these emitters, 16 lie in the MIRI coverage area and 12 exhibit a clear flux excess at $5.6 \, \rm \mu m$, indicating the simultaneous presence of a prominent H$\alpha$ emission line with EW$_{0}$(H$\alpha$) $\simeq 200-3000$ {\AA}. This is the first time that H$\alpha$ emission can be detected in individual galaxies at $z>7$. The H$\alpha$ line, when present, allows us to separate the contributions of H$\beta$ and [OIII] to the (H$\beta$ +[OIII]) complex, and derive H$\alpha$-based star formation rates (SFRs). We find that in most cases [OIII]/H$\beta > 1$. Instead, two galaxies have [OIII]/H$\beta < 1$, indicating that the NIRCam flux excess is mainly driven by H$\beta$. This could potentially imply extremely low metallicities. Most prominent line emitters are very young starbursts or galaxies on their way to/from the starburst cloud. They make for a cosmic SFR density $\rm log_{10}(\rho_{SFR_{H\alpha}}) \simeq -2.35$, which is about a quarter of the total value ($\rm log_{10}(\rho_{SFR_{tot}}) \simeq -1.76$) at $z\simeq 7-8$. Therefore, the strong H$\alpha$ emitters likely had a significant role in reionization., Comment: 18 pages, 11 figures. The paper has been accepted for publication in ApJ

We study the young star cluster populations in 23 dwarf and irregular galaxies observed by the HST Legacy ExtraGalactic Ultraviolet Survey (LEGUS), and examine relationships between the ensemble properties of the cluster populations and those of their host galaxies: star formation rate (SFR) density ($\Sigma_{\rm SFR}$). A strength of this analysis is the availability of SFRs measured from temporally resolved star formation histories which provide the means to match cluster and host-galaxy properties on several timescales (1-10, 1-100, and 10-100~Myr). Nevertheless, studies of this kind are challenging for dwarf galaxies due to the small numbers of clusters in each system. We mitigate these issues by combining the clusters across different galaxies with similar $\Sigma_{\rm SFR}$ properties. We find good agreement with a well-established relationship ($M_{V}^{brightest}$-SFR), but find no significant correlations between $\Sigma_{\rm SFR}$ and the slopes of the cluster luminosity function, mass function, nor the age distribution. We also find no significant trend between the the fraction of stars in bound clusters at different age ranges ($\Gamma_{1-10}$, $\Gamma_{10-100}$, and $\Gamma_{1-100}$) and $\Sigma_{\rm SFR}$ of the host galaxy. Our data show a decrease in $\Gamma$ over time (from 1-10 to 10-100~Myr) suggesting early cluster dissolution, though the presence of unbound clusters in the youngest time bin makes it difficult to quantify the degree of dissolution. While our data do not exhibit strong correlations between $\Sigma_{\rm SFR}$ and ensemble cluster properties, we cannot rule out that a weak trend might exist given the relatively large uncertainties due to low number statistics and the limited $\Sigma_{\rm SFR}$ range probed., Comment: 30 pages, 15 figures, 7 tables, Accepted for publication in MNRAS

Context: Outflows from low-mass star-forming galaxies are a fundamental ingredient for models of galaxy evolution and cosmology. Aims: The onset of kpc-scale ionised filaments in the halo of the metal-poor compact dwarf SBS 0335-052E was previously not linked to an outflow. We here we investigate whether these filaments provide evidence for an outflow. Methods: We obtained new VLT/MUSE WFM and deep NRAO/VLA B-configuration 21cm data of the galaxy. The MUSE data provide morphology, kinematics, and emission line ratios H$\beta$/H$\alpha$ and [\ion{O}{iii}]$\lambda5007$/H$\alpha$ of the low surface-brightness filaments, while the VLA data deliver morphology and kinematics of the neutral gas in and around the system. Both datasets are used in concert for comparisons between the ionised and the neutral phase. Results: We report the prolongation of a lacy filamentary ionised structure up to a projected distance of 16 kpc at $\mathrm{SB}_\mathrm{H\alpha} = 1.5\times10^{-18}$erg s$^{-1}$ cm$^{-2}$arcsec$^{-2}$. The filaments exhibit unusual low H$\alpha$/H$\beta \approx 2.4$ and low [\ion{O}{iii}]/H$\alpha \sim 0.4 - 0.6$ typical of diffuse ionised gas. They are spectrally narrow ($\sim 20$ km s$^{-1}$) and exhibit no velocity sub-structure. The filaments extend outwards of the elongated \ion{H}{I} halo. On small scales the $N_\mathrm{HI}$ peak is offset from the main star-forming sites. Morphology and kinematics of \ion{H}{I} and \ion{H}{II} reveal how star-formation driven feedback interacts differently with the ionised and the neutral phase. Conclusions: We reason that the filaments are a large scale manifestation of star-formation driven feedback, namely limb-brightened edges of a giant outflow cone that protrudes through the halo of this gas-rich system. A simple toy model of such a conical-structure is found to be commensurable with the observations., Comment: Accepted version in A&A after language editing. 22 pages, 24 figures

The CLusters in the Uv as EngineS (CLUES) survey is a Cosmic Origins Spectrograph (COS) campaign aimed at acquiring the 1130 to 1770 {\AA}, restframe spectroscopy of very young (<20 Myr) and massive (>10^4 solar masses) star clusters in galaxies that are part of the Hubble treasury program Legacy ExtraGalactic Uv Survey (LEGUS). In this first paper of a series, we describe the CLUES sample consisting of 20 young star clusters and report their physical properties as derived by both multi-wavelength photometry and far-UV (FUV) spectroscopy with Hubble Space Telescope (HST). Thanks to the synergy of the two different datasets we build a coherent picture of the diverse stellar populations found in each region (with sizes of 40 to 160 pc). We associate the FUV-brightest stellar population to the central targeted star cluster and the other modeled population to the diffuse stars that are included in the COS aperture. We observe better agreement between photometric and spectroscopic ages for star clusters younger than 5 Myr. For clusters older than 5 Myr, photometry and spectroscopy measurements deviate, with the latter producing older ages, due to the degeneracy of photometric models. FUV spectroscopy enables us to better constrain the stellar metallicities, a parameter that optical colors are insensitive to. Finally, the derived E(B-V) are quite similar, with a tendency for FUV spectroscopy to favor solutions with higher extinctions. The recovered masses are in agreement within a factor of 2 for all the clusters., Comment: Accepted for publication on the Astronomical Journal on 16th September 2022. 26 pages, 12 figures

The Ly$\alpha$ emission line is one of the main observables of galaxies at high redshift, but its output depends strongly on the neutral gas distribution and kinematics around the star-forming regions where UV photons are produced. We present observations of Ly$\alpha$ and 21-cm HI emission at comparable scales with the goal to qualitatively investigate how the neutral interstellar medium (ISM) properties impact Ly$\alpha$ transfer in galaxies. We have observed 21-cm HI at the highest angular resolution possible (~ 3" beam) with the VLA in two local galaxies from the Lyman Alpha Reference Sample. We contrast this data with HST Ly$\alpha$ imaging and spectroscopy, and MUSE and PMAS ionized gas observations. In LARS08, high intensity Ly$\alpha$ emission is co-spatial with high column density HI where dust content is the lowest. The Ly$\alpha$ line is strongly redshifted, consistent with velocity redistribution which allows Ly$\alpha$ escape from high column density neutral medium with low dust content. In eLARS01, high intensity Ly$\alpha$ emission is located in regions of low column density HI, below the HI data sensitivity limit ($<2\times10^{20}\,$cm$^{-2}$). The perturbed ISM distribution with low column density gas in front of the Ly$\alpha$ emission region plays an important role in the escape. In both galaxies, the faint Ly$\alpha$ emission ($\sim 1\times10^{-16}$erg.s$^{-1}$cm$^{-2}$arcsec$^{-2}$) traces intermediate H$\alpha$ emission regions where HI is found, regardless of the dust content. Dust seems to modulate, but not prevent, the formation of a faint Ly$\alpha$ halo. This study suggests the existence of scaling relations between dust, H$\alpha$, HI, and Ly$\alpha$ emission in galaxies., Comment: 28 pages, 14 figures, accepted for publication in ApJ

The properties that govern the production and escape of hydrogen ionizing photons (Lyman continuum, LyC; with energies >13.6 eV) in star-forming galaxies are still poorly understood, but they are key to identifying and characterizing the sources that reionized the Universe. Here we empirically explore the relationship between the hardness of ionizing radiation and the LyC leakage in a large sample of low-$z$ star-forming galaxies from the recent Hubble Space Telescope Low-$z$ Lyman Continuum Survey. Using Sloan Digital Sky Survey stacks and deep XShooter observations, we investigate the hardness of the ionizing spectra ($Q_{\rm He^+}/Q_{\rm H}$) between 54.4 eV (He$^{+}$) and 13.6 eV (H) from the optical recombination lines HeII 4686A and H$\beta$ 4861A for galaxies with LyC escape fractions spanning a wide range, $f_{\rm esc} \rm (LyC) \simeq 0 - 90\%$. We find that the observed intensity of HeII/H$\beta$ is primarily driven by variations in the metallicity, but is not correlated with LyC leakage. Both very strong ($ \simeq 0.5$) and nonleakers ($ < f_{\rm esc} \rm (LyC) > \simeq 0$) present similar observed intensities of HeII/H$\beta$ at comparable metallicity, between $\simeq 0.01$ and $\simeq 0.02$ for $12 + \log({\rm O/H}) > 8.0$ and $<8.0$, respectively. Our results demonstrate that $Q_{\rm He^+}/Q_{\rm H}$ does not correlate with $f_{\rm esc} \rm (LyC)$, which implies that strong LyC emitters do not show harder ionizing spectra than nonleakers at similar metallicity., Comment: 6 pages, 2 figures, 1 tables. Accepted for publication in A&A Letters. v2: updated to match published version, including new title

Star formation is a clustered process that regulates the structure and evolution of galaxies. We investigate this process in the dwarf galaxy Haro 11, forming stars in three knots (A, B, C). The exquisite resolution of HST imaging allows us to resolve the starburst into tens of bright star clusters. We derive masses between $10^5$ and $10^7\,\rm M_{\odot}$ and ages younger than 20 Myr, using photometric modeling. We observe that the clustered star formation has propagated from knot C (the oldest) through knot A (in between) towards knot B (the youngest). We use aperture-matched ultraviolet and optical spectroscopy (HST + MUSE) to independently study the stellar populations of Haro 11 and determine the physical properties of the stellar populations and their feedback in 1 kpc diameter regions. We discuss these results in light of the properties of the ionised gas within the knots. We interpret the broad blue-shifted components of the optical emission lines as outflowing gas ($v_{max} \sim 400$ km/s). The strongest outflow is detected in knot A with a mass-rate of $\dot{M}_{out}\sim 10\,\rm M_{\odot}/yr$, ten times higher than the star-formation in the same region. Knot B hosts a young and not fully developed outflow, whereas knot C has likely been already evacuated. Because Haro 11 has properties similar to high-redshift unresolved galaxies, our work can additionally aid the understanding of star formation at high redshift, a window that will be opened by upcoming facilities., Comment: Accepted for publication in MNRAS on 22nd Dec 2022: 18 pages, 9 figures

We present the Hubble imaging Probe of Extreme Environments and Clusters (HiPEEC) survey. We fit HST NUV to NIR broadband and H$\alpha$ fluxes, to derive star cluster ages, masses, extinctions and determine the star formation rate (SFR) of 6 merging galaxies. These systems are excellent laboratories to trace cluster formation under extreme gas physical conditions, rare in the local universe, but typical for star-forming galaxies at cosmic noon. We detect clusters with ages of 1-500 Myr and masses that exceed $10^7$ M$_\odot$. The recent cluster formation history and their distribution within the host galaxies suggest that systems like NGC34, NGC1614, NGC4194 are close to their final coalescing phase, while NGC3256, NGC3690, NGC6052 are at an earlier/intermediate stage. A Bayesian analysis of the cluster mass function in the age interval 1-100 Myr provides strong evidence in 4 of the 6 galaxies that an exponentially truncated power law better describes the observed mass distributions. For two galaxies, the fits are inconclusive due to low number statistics. We determine power-law slopes $\beta \sim-1.5$ to $-2.0$, and truncation masses, M$_c$, between $10^6$ and a few times $10^7$ M$_\odot$, among the highest values reported in the literature. Advanced mergers have higher M$_c$ than early/intermediate merger stage galaxies, suggesting rapid changes in the dense gas conditions during the merger. We compare the total stellar mass in clusters to the SFR of the galaxy, finding that these systems are among the most efficient environments to form star clusters in the local universe., Comment: Accepted for publication in MNRAS; data and catalogues will be released at this URL http://dx.doi.org/10.17909/t9-cn0b-ht83

We study the global dust and (molecular) gas content in the Lyman Alpha Reference Sample (LARS), i.e. 14 local star-forming galaxies. We characterize their interstellar medium and relate newly derived properties to quantities relevant for Ly$\alpha$ escape. We observed LARS galaxies with Herschel/PACS, SOFIA/FIFI-LS, the IRAM 30m telescope and APEX, targeting far-infrared (FIR) continuum and emission lines of [C II]158$\mu$m, [O I]63$\mu$m, [O III]88$\mu$m and low-J CO lines. Using Bayesian methods we derive dust model parameters and estimate total gas masses for all LARS galaxies, taking into account a metallicity-dependent gas-to-dust ratio. Star formation rates were estimated from FIR, [C II]158$\mu$m and [O I]63$\mu$m luminosities. LARS covers a wide dynamic range in the derived properties, with FIR-based star formation rates from $\sim$0.5-100 $M_{\odot}\ yr^{-1}$, gas fractions between $\sim$15-80% and gas depletion times ranging from a few hundred Myr up to more than 10 Gyr. The distribution of LARS galaxies in the $\Sigma_{gas}$ vs. $\Sigma_{SFR}$ (Kennicutt-Schmidt plane) is thus quite heterogeneous. However, we find that LARS galaxies with the longest gas depletion times, i.e. relatively high gas surface densities ($\Sigma_{gas}$) and low star formation rate densities ($\Sigma_{SFR}$), have by far the highest Ly$\alpha$ escape fraction. A strong $\sim$linear relation is found between Ly$\alpha$ escape fraction and the total gas (HI+H$_2$) depletion time. We argue that the Ly$\alpha$ escape in those galaxies is driven by turbulence in the star-forming gas that shifts the Ly$\alpha$ photons out of resonance close to the places where they originate. We further report on an extreme [C II]158$\mu$m excess in LARS 5, corresponding to $\sim$14$\pm$3% of the FIR luminosity, i.e. the most extreme [C II]-to-FIR ratio observed in a non-AGN galaxy to date.

Ly$\alpha$ emission is a standard tracer of starburst galaxies at high redshift. However, a number of local Ly$\alpha$ emitters (LAEs) are X-ray sources, suggesting a possible origin of Ly$\alpha$ photons other than young, hot stars, and which may be active at much later ages relative to the parent starburst. Resolved, nearby LAEs offer the opportunity to discriminate between diffuse X-ray emission arising from supernova-heated gas, high-mass X-ray binaries (HMXBs), or low-luminosity active galactic nuclei (LLAGN). We examine archival X-ray imaging from Chandra and XMM-Newton for 11 galaxies with spatially resolved Ly$\alpha$ imaging to determine the luminosity, morphology, and spectral hardness of the X-ray sources. The data are consistent with 9 of the 12, bright Ly$\alpha$ sources being driven by luminous, $10^{40}$ erg s$^{-1}$ X-ray sources. Half of the 8 Chandra sources are unresolved. The data suggest that nuclear activity, whether from LLAGN or nuclear starbursts, may play an important role in Ly$\alpha$ emission. Our results also suggest a significant link between Ly$\alpha$ emission and HMXBs, ULXs, and/or LLAGN, which would imply that Ly$\alpha$ may be generated over timescales 1 - 2 orders of magnitude longer than produced by photoionization from OB stars. This highlights a critical need to quantify the relative contributions of different sources across cosmic time, to interpret Ly$\alpha$ observations and the resulting properties of distant galaxies., Comment: 10 pages, 3 figures, 2 tables. Accepted by ApJ Letters

The James Webb Space Telescope will provide deep imaging and spectroscopy for sources at redshifts above 6, covering the Epoch of Reionization (EoR, 6 < z < 10). The Mid-IR instrument (MIRI) integral field spectrograph (MRS) will be the only instrument on board JWST able to observe the brightest optical emission lines H$\alpha$ and [OIII]0.5007$\mu$m at redshifts above 7 and 9, respectively. This paper presents a study of the H$\alpha$ fluxes predicted by FIRSTLIGHT cosmological simulations for galaxies at redshifts of 6.5 to 10.5, and its detectability with MIRI. Deep (40 ks) spectroscopic integrations with MRS will be able to detect (S/N > 5) EoR sources at redshifts above 7 with intrinsic star formation rates of more than 2 M$_{\odot}$ yr$^{-1}$, and stellar masses above 4-9 $\times$ 10$^7$ M$_{\odot}$. In addition, the paper presents realistic MRS simulated observations of the expected (rest-frame) optical and near-infrared spectra for some spectroscopically confirmed EoR sources detected by ALMA as [OIII]88$\mu$m emitters. The MRS simulated spectra cover a wide range of low metallicities from about 0.2 to 0.02Z$_{\odot}$, and different [OIII]88$\mu$m/[OIII]0.5007$\mu$m line ratios. The simulated 10ks MRS spectra show S/N in the range of 5 to 90 for H$\beta$, [OIII]0.4959,0.5007$\mu$m, H$\alpha$ and HeI1.083$\mu$m emission lines of MACS1149-JD1 at z = 9.11, independent of metallicity. In addition, deep 40 ks simulated spectra of the luminous merger candidate B14-65666 at z=7.15 shows the MRS capabilities of detecting, or putting strong upper limits, on the [NII]0.6584$\mu$m, [SII]0.6717,0.6731$\mu$m, and [SIII]0.9069,0.9532$\mu$m emission lines. In summary, MRS will enable the detailed study of key physical properties like internal extinction, instantaneous star formation, hardness of the ionising continuum, and metallicity, in bright (intrinsic or lensed) EoR sources., Comment: 11 pages, 6 figures, Accepted for publication in A&A

We have used the capability of the MUSE instrument to explore the impact of stellar feedback at large scales in Haro 11, a galaxy under extreme starburst condition and one of the first galaxies where Lyman continuum (LyC) has been detected. Using Ha, [OIII] and [OI] emission lines from deep MUSE observations, we have constructed a sequence of velocity-dependent maps of the Ha emission, the state of the ionised gas and a tracer of fast shocks. These allowed us to investigate the ionisation structure of the galaxy in 50 kms^2 bins over a velocity range of -400 to 350 kms. The ionised gas in Haro 11 is assembled by a rich arrangement of structures, such as superbubbles, filaments, arcs and galactic ionised channels, whose appearances change drastically with velocity. The central star forming knots and the star forming dusty arm are the main engines that power the strong mechanical feedback in this galaxy, although with different impact on the ionisation structure. Haro 11 appears to leak LyC radiation in many directions. We found evidence of a kpc-scale fragmented superbubble, that may have cleared galactic-scale channels in the ISM. Additionally, the southwestern hemisphere is highly ionised in all velocities, hinting at a density bound scenario. A compact kpc-scale structure of lowly ionised gas coincides with the diffuse Lya emission and the presence of fast shocks. Finally, we find evidence that a significant fraction of the ionised gas mass may escape the gravitational potential of the galaxy., Comment: Pubisched version. Ionisation values were corrected after discovering a bug in a code used to generate the maps

Stellar feedback strongly affects the interstellar medium (ISM) of galaxies. Stellar feedback in the first galaxies likely plays a major role in enabling the escape of LyC photons, which contribute to the re-ionization of the Universe. Nearby starburst galaxies serve as local analogues allowing for a spatially resolved assessment of the feedback processes in these galaxies. We characterize the feedback effects from the star clusters in the local high-redshift analogue ESO 338-IG04 on the ISM and compare the results with the properties of the most massive clusters. We use high quality VLT/MUSE optical integral field data to derive the physical properties of the ISM such as ionization, density, shocks, and perform new fitting of the spectral energy distributions of the brightest clusters in ESO 338-IG04 from HST imaging. ESO 338-IG04 has a large ionized halo which we detect to a distance of 9 kpc. We identify 4 Wolf-Rayet (WR) clusters based on the blue and red WR bump. We follow previously identified ionization cones and find that the ionization of the halo increases with distance. Analysis of the galaxy kinematics shows two complex outflows driven by the numerous young clusters in the galaxy. We find a ring of shocked emission traced by an enhanced [OI]/H$\alpha$ ratio surrounding the starburst and at the end of the outflow. Finally we detect nitrogen enriched gas associated with the outflow, likely caused by the WR stars in the massive star clusters. Photo-ionization dominates the central starburst and sets the ionization structure of the entire halo, resulting in a density bounded halo, facilitating the escape of LyC photons. Outside the central starburst, shocks triggered by an expanding super bubble become important. The shocks at the end of the outflow suggest interaction between the hot outflowing material and the more quiescent halo gas., Comment: Accepted for publication in Astronomy and Astrophysics, 22 pages, 15 figures

Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5x10^52 erg in the infrared and radio, but was not luminous at optical or X-ray wavelengths. We interpret this as a TDE with much of its emission re-radiated at infrared wavelengths by dust. Efficient reprocessing by dense gas and dust may explain the difference between theoretical predictions and observed luminosities of TDEs. The radio observations resolve an expanding and decelerating jet, probing the jet formation and evolution around a SMBH., Comment: Published online by Science on 14 Jun 2018

It has recently been established that the properties of young star clusters (YSCs) can vary as a function of the galactic environment in which they are found. We use the cluster catalogue produced by the Legacy Extragalactic UV Survey (LEGUS) collaboration to investigate cluster properties in the spiral galaxy M51. We analyse the cluster population as a function of galactocentric distance and in arm and inter-arm regions. The cluster mass function exhibits a similar shape at all radial bins, described by a power law with a slope close to $-2$ and an exponential truncation around $10^5\ \rm{M}_{\odot}$ . While the mass functions of the YSCs in the spiral arm and inter-arm regions have similar truncation masses, the inter-arm region mass function has a significantly steeper slope than the one in the arm region; a trend that is also observed in the giant molecular cloud mass function and predicted by simulations. The age distribution of clusters is dependent on the region considered, and is consistent with rapid disruption only in dense regions, while little disruption is observed at large galactocentric distances and in the inter-arm region. The fraction of stars forming in clusters does not show radial variations, despite the drop in the $H_2$ surface density measured as function of galactocentric distance. We suggest that the higher disruption rate observed in the inner part of the galaxy is likely at the origin of the observed flat cluster formation efficiency radial profile., Comment: 22 pages, 18 figures, 9 tables (+appendix: 6 pages, 3 figures, 4 tables). Accepted for publication in MNRAS

We investigate the relation between gas and star formation in sub-galactic regions, ~360 pc to ~1.5 kpc in size, within the nearby starburst dwarf NGC4449, in order to separate the underlying relation from the effects of sampling at varying spatial scales. Dust and gas mass surface densities are derived by combining new observations at 1.1 mm, obtained with the AzTEC instrument on the Large Millimeter Telescope, with archival infrared images in the range 8-500 micron from the Spitzer Space Telescope and the Herschel Space Observatory. We extend the dynamic range of our mm (and dust) maps at the faint end, using a correlation between the far-infrared/millimeter colors F(70)/F(1100) [and F(160)/F(1100)] and the mid-infrared color F(8)/F(24) that we establish for the first time for this and other galaxies. Supplementing our data with maps of the extinction-corrected star formation rate (SFR) surface density, we measure both the SFR-molecular gas and the SFR-total gas relations in NGC4449. We find that the SFR-molecular gas relation is described by a power law with exponent that decreases from ~1.5 to ~1.2 for increasing region size, while the exponent of the SFR-total gas relation remains constant with value ~1.5 independent of region size. We attribute the molecular law behavior to the increasingly better sampling of the molecular cloud mass function at larger region sizes; conversely, the total gas law behavior likely results from the balance between the atomic and molecular gas phases achieved in regions of active star formation. Our results indicate a non-linear relation between SFR and gas surface density in NGC4449, similar to what is observed for galaxy samples., Comment: 63 pages, 4 tables, 17 figures. Accepted for publication on the Astrophysical Journal

Recently acquired WFC3 UV (F275W and F336W) imaging mosaics under the Legacy Extragalactic UV Survey (LEGUS) combined with archival ACS data of M51 are used to study the young star cluster (YSC) population of this interacting system. Our newly extracted source catalogue contains 2834 cluster candidates, morphologically classified to be compact and uniform in colour, for which ages, masses and extinction are derived. In this first work we study the main properties of the YSC population of the whole galaxy, considering a mass-limited sample. Both luminosity and mass functions follow a power law shape with slope -2, but at high luminosities and masses a dearth of sources is observed. The analysis of the mass function suggests that it is best fitted by a Schechter function with slope -2 and a truncation mass at $1.00\pm0.12\times10^5\ M_\odot$. Through Monte Carlo simulations we confirm this result and link the shape of the luminosity function to the presence of a truncation in the mass function. A mass limited age function analysis, between 10 and 200 Myr, suggests that the cluster population is undergoing only moderate disruption. We observe little variation in the shape of the mass function at masses above $1\times10^4\ M_\odot$, over this age range. The fraction of star formation happening in the form of bound clusters in M51 is $\sim20\%$ in the age range 10 to 100 Myr and little variation is observed over the whole range from 1 to 200 Myr., Comment: 25 pages, 19 figures, 10 tables. Accepted for publication in MNRAS (14 Sep. 2017)

We report on the discovery of ionised gas filaments in the circum-galactic halo of the extremely metal-poor compact starburst SBS 0335-052E in a 1.5h integration with the MUSE integral-field spectrograph. We detect these features in H${\alpha}$ and [OIII] emission down to surface-brightness levels of $5 \times 10^{-19}$erg s$^{-1}$cm$^{-2}$arcsec$^{-2}$. The filaments have projected diameters of 2.1 kpc and extend more than 9 kpc to the north and north-west from the main stellar body. We also detect extended nebular HeII $\lambda$4686 emission that brightens towards the north-west at the rim of a star-burst driven super-shell, suggestive of a locally enhanced UV radiation field due to shocks. We also present a velocity field of the ionised gas. The filaments appear to connect seamlessly in velocity space to the kinematical disturbances caused by the shell. Similar to high-$z$ star-forming galaxies, the ionised gas in this galaxy is dispersion dominated. We argue that the filaments were created via feedback from the starburst and that these ionised structures in the halo may act as escape channels for Lyman continuum radiation in this gas-rich system., Comment: Revised version after peer review. Accepted for publication in A&A letters

We test the predictions of spectral synthesis models based on seven different massive-star prescriptions against Legacy ExtraGalactic UV Survey (LEGUS) observations of eight young massive clusters in two local galaxies, NGC 1566 and NGC 5253, chosen because predictions of all seven models are available at the published galactic metallicities. The high angular resolution, extensive cluster inventory and full near-ultraviolet to near-infrared photometric coverage make the LEGUS dataset excellent for this study. We account for both stellar and nebular emission in the models and try two different prescriptions for attenuation by dust. From Bayesian fits of model libraries to the observations, we find remarkably low dispersion in the median E(B-V) (~0.03 mag), stellar masses (~10^4 M_\odot) and ages (~1 Myr) derived for individual clusters using different models, although maximum discrepancies in these quantities can reach 0.09 mag and factors of 2.8 and 2.5, respectively. This is for ranges in median properties of 0.05-0.54 mag, 1.8-10x10^4 M_\odot and 1.6-40 Myr spanned by the clusters in our sample. In terms of best fit, the observations are slightly better reproduced by models with interacting binaries and least well reproduced by models with single rotating stars. Our study provides a first quantitative estimate of the accuracies and uncertainties of the most recent spectral synthesis models of young stellar populations, demonstrates the good progress of models in fitting high-quality observations, and highlights the needs for a larger cluster sample and more extensive tests of the model parameter space., Comment: Accepted for publication in MNRAS (14 Jan. 2016). 30 pages, 16 figures, 9 tables