Search

We analysed ALMA FIR (1.3 mm) dust continuum and CO emission of 12 starburst galaxies at $z\sim 2.1-3.6$, selected for their extreme brightness in the rest-UV with $M_{\rm UV} = -23.4$ to $-24.7$. We also analysed VLT HAWK-I $H$- and $K_{\rm s}$-band images. The galaxies are characterised by negligible dust attenuations with blue UV spectral slopes ($-2.62$ to $-1.84$), very young stellar populations of $\sim 10$ Myr, and powerful starbursts with a high mean specific star formation rate of $\rm 112~Gyr^{-1}$, placing them $\sim 1.5$~dex above the main sequence at similar redshifts and stellar masses ($M_{\rm stars} \sim (1.5-4.6)\times 10^9~M_{\odot}$). The FIR dust continuum emission revealed in 9 galaxies yields IR luminosities of $(5.9-28.3)\times 10^{11}~L_{\odot}$ and large dust masses barely produced by SNe within the 10~Myr timescale. The CO emission detected in 8 galaxies evidence large molecular gas masses with a mean molecular gas fraction of 82%. The corresponding star formation efficiencies reach $\gtrsim 40$\%, with amazingly short molecular gas depletion timescales between <13 Myr and 71 Myr. These unique properties, never reported in previously studied galaxies, highlight that these galaxies are likely caught at the very beginning of their stellar mass build-up and undergo a very efficient and fast conversion of gas into stars that can only result from the gas collapse within very short free-fall times. We find that the feedback-free starburst model seems to be able to explain the formation of these galaxies. To reconcile the co-spatial FIR dust emission with the UV-bright unattenuated emission, we speculate about radiation-driven outflows which can temporarily remove dust at the location of the starburst and expel dust at large distances in line with the measured large FIR effective radii ($\rm 1.7~kpc - 5~kpc$) in comparison to very compact stellar radii., Comment: 20 pages, 10 figures, 3 tables. Submitted to Astronomy & Astrophysics

J1316+2614 at z=3.613 is the UV-brightest ($M_{UV}$ = -24.7) and strongest Lyman continuum (LyC, $f_{esc}^{LyC} \approx$ 90%) emitting star-forming galaxy known, showing also signatures of inflowing gas from its blue-dominated Ly$\alpha$ profile. Here, we present high-resolution imaging with the HST and VLT of the LyC, Ly$\alpha$, rest-UV, and optical emission of J1316+2614. Detailed analysis of the LyC and UV light distributions reveals compact yet resolved profiles, with LyC and UV morphologies showing identical half-light radii of $\simeq$ 220 pc. The continuum-subtracted Ly$\alpha$ emission reveals an extended filamentary structure of $\simeq$ 6.0 kpc oriented south-north with only weak/residual flux within the stellar core, suggesting a Ly$\alpha$ "hole". J1316+2614 presents remarkably high SFR and stellar mass surface densities of log($\Sigma_{SFR}$ [$M_{\odot}$/yr/kpc^2]) = 3.47$\pm$0.11 and log($\Sigma_{M}$ [$M_{\odot}$/pc^2]) = 4.20$\pm$0.06, respectively, which are among the highest observed in star-forming galaxies. Our findings indicate that J1316+2614 is a powerful, young, and compact starburst, leaking significant LyC photons due to the lack of gas and dust within the starburst. We explore the conditions for gas expulsion using a simple energetic balance and find that, given the strong binding force in J1316+2614, a high star formation efficiency ($\epsilon_{SF} \geq 0.7$) is necessary to remove the gas and explain its exposed nature. Our results thus suggest a close link between high $\epsilon_{SF}$ and high $f_{esc}^{LyC}$. This high efficiency can also naturally explain the remarkably high SFR, UV-luminosity, and efficient mass growth of J1316+2614, where at least 62% of its mass formed in the last 6 Myr. J1316+2614 may exemplify an intense, feedback-free starburst with a high $\epsilon_{SF}$, similar to those proposed for UV-bright galaxies at high redshifts. (ABRIDGED), Comment: 15 pages and 9 figures. Accepted for publication in A&A

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

We investigate the correlation between stellar mass (M*) and star formation rate (SFR) across the stellar mass range log10(M*/Msun)~6-11. We consider almost 50,000 star-forming galaxies at z~3-7, leveraging data from COSMOS/SMUVS, JADES/GOODS-SOUTH, and MIDIS/XDF. This is the first study spanning such a wide stellar mass range without relying on gravitational lensing effects. We locate our galaxies on the SFR-M* plane to assess how the location of galaxies in the star-formation main sequence (MS) and starburst (SB) region evolves with stellar mass and redshift. We find that the two star-forming modes tend to converge at log10(M*/Msun) < 7, with all galaxies found in the SB mode. By dissecting our galaxy sample in stellar mass and redshift, we show that the emergence of the star-formation MS is stellar-mass dependent: while in galaxies with log10(M*/Msun) > 9 the MS is already well in place at z = 5-7, for galaxies with log10(M*/Msun)~7-8 it only becomes significant at z<4. Overall, our results are in line with previous findings that the SB mode dominates amongst low stellar-mass galaxies. The earlier emergence of the MS for massive galaxies is consistent with galaxy downsizing., Comment: 15 pages, 4 figures. Submitted to ApJ

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

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

Astrometry is one of the main pillars of astronomy, and one of its oldest branches. Over the years, an increasing number of astrometric works by means of Hubble Space Telescope (HST) data have revolutionized our understanding of various phenomena. With the launch of JWST, it becomes almost instinctive to want to replicate or improve these results with data taken with the newest, state-of-the-art, space-based telescope. In this regard, the initial focus of the community has been on the Near-Infrared (NIR) detectors on board of JWST because of their high spatial resolution. This paper begins the effort to capture and apply what has been learned from HST to the Mid-InfraRed Instrument (MIRI) of JWST by developing the tools to obtain high-precision astrometry and photometry with its imager. We describe in detail how to create accurate effective point-spread-function (ePSF) models and geometric-distortion corrections, analyze their temporal stability, and test their quality to the extent of what is currently possible with the available data in the JWST MAST archive. We show that careful data reduction provides deep insight on the performance and intricacies of the MIRI imager, and of JWST in general. In an effort to help the community to devise new observing programs, we make our ePSF models and geometric-distortion corrections publicly available., Comment: 22 pages, 18 figures, 1 table. Accepted for publication in PASP. The ePSF models, geometric-distortion solutions and codes are available at the links provided in the manuscript

By using the ultra-deep \textit{JWST}/MIRI image at 5.6 $\mu m$ in the Hubble eXtreme Deep Field, we constrain the role of strong H$\alpha$-emitters (HAEs) during Cosmic Reionization at $z\simeq7-8$. Our sample of HAEs is comprised of young ($<35\;\rm Myr$) galaxies, except for one single galaxy ($\approx 300\;\rm Myr$), with low stellar masses ($\lesssim 10^{9}\;\rm M_{\odot}$). These HAEs show a wide range of UV-$\beta$ slopes, with a median value of $\beta = -2.15\pm0.21$ which broadly correlates with stellar mass. We estimate the ionizing photon production efficiency ($\xi_{ion,0}$) of these sources (assuming $f_{esc,LyC} = 0\%$), which yields a median value $\rm log_{10}(\xi_{ion,0}/(Hz\;erg^{-1})) = 25.50^{+0.10}_{-0.12}$. We show that $\xi_{ion,0}$ positively correlates with EW$_{0}$(H$\alpha$) and specific star formation rate (sSFR). Instead $\xi_{ion,0}$ weakly anti-correlates with stellar mass and $\beta$. Based on the $\beta$ values, we predict $f_{esc, LyC}=4\%^{+3}_{-2}$, which results in $\rm log_{10}(\xi_{ion}/(Hz\;erg^{-1})) = 25.55^{+0.11}_{-0.13}$. Considering this and related findings from the literature, we find a mild evolution of $\xi_{ion}$with redshift. Additionally, our results suggest that these HAEs require only modest escape fractions ($f_{esc, rel}$) of 6$-$15\% to reionize their surrounding intergalactic medium. By only considering the contribution of these HAEs, we estimated their total ionizing emissivity ($\dot{N}_{ion}$) as $\dot{N}_{ion} = 10^{50.53 \pm 0.45}; \text{s}^{-1}\text{Mpc}^{-3}$. When comparing their $\dot{N}_{ion}$ with "non-H$\alpha$ emitter" galaxies across the same redshift range, we find that that strong, young, and low-mass emitters may have played an important role during Cosmic Reionization., Comment: 21 pages, 11 Figures, 1 table. Accepted for publication in ApJ. Updated version 13/06/2024 -- found a mistake in the reported eq. 2 (now corrected)

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

The Medium-Resolution integral field Spectrometer (MRS) of MIRI on board JWST performs spectroscopy between 5 and 28~$\mu$m. The optics of the MRS introduce substantial distortion, and this needs to be rectified in order to reconstruct the observed astrophysical scene. We use data from the JWST/MIRI commissioning and cycle 1 calibration phase, to derive the MRS geometric distortion and astrometric solution, a critical step in the calibration of MRS data. These solutions come in the form of transform matrices that map the detector pixels to spatial coordinates of a local MRS coordinate system called $\alpha$/$\beta$, to the global JWST observatory coordinates V2/V3. For every MRS spectral band and each slice dispersed on the detector, the transform of detector pixels to $\alpha$/$\beta$ is fit by a two-dimensional polynomial, using a raster of point source observations. A polynomial transform is used to map the coordinates from $\alpha$/$\beta$ to V2/V3. We calibrated the distortion of all 198 discrete slices of the MIRI/MRS IFUs, and derived an updated Field of View (FoV) for each MRS spectral band. The precision of the distortion solution is estimated to be better than one tenth of a spatial resolution element, with a root mean square (rms) of 10 milli-arcsecond (mas) at 5 $\mu$m, to 23 mas at 27 $\mu$m. Finally we find that the wheel positioning repeatability causes an additional astrometric error of rms 30 mas. We have demonstrated the MRS astrometric calibration strategy and analysis enabling the calibration of MRS spectra, a critical step in the data pipeline especially for science with spatially resolved objects. The distortion calibration was folded into the JWST pipeline in Calibration Reference Data System (CRDS) context jwst\_1094.pmap. The distortion calibration precision meets the pre-launch requirement, and the estimated total astrometric uncertainty is 50 mas., Comment: 15 pages, 14 figures. Submitted to Astronomy & Astrophysics

We describe an algorithm for application of the classic `drizzle' technique to produce 3d spectral cubes using data obtained from the slicer-type integral field unit (IFU) spectrometers on board the James Webb Space Telescope. This algorithm relies upon the computation of overlapping volume elements (composed of two spatial dimensions and one spectral dimension) between the 2d detector pixels and the 3d data cube voxels, and is greatly simplified by treating the spatial and spectral overlaps separately at the cost of just 0.03% in spectrophotometric fidelity. We provide a matrix-based formalism for the computation of spectral radiance, variance, and covariance from arbitrarily dithered data and comment on the performance of this algorithm for the Mid-Infrared Instrument's Medium Resolution IFU Spectrometer (MIRI MRS). We derive a series of simplified scaling relations to account for covariance between cube spaxels in spectra extracted from such cubes, finding multiplicative factors ranging from 1.5 to 3 depending on the wavelength range and kind of data cubes produced. Finally, we discuss how undersampling produces periodic amplitude modulations in the extracted spectra in addition to those naturally produced by fringing within the instrument; reducing such undersampling artifacts below 1% requires a 4-point dithering strategy and spectral extraction radii of 1.5 times the PSF FWHM or greater., Comment: 16 pages, 12 figures. Revised version resubmitted to AJ

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

During the commissioning of {\em JWST}, the Medium-Resolution Spectrometer (MRS) on the Mid-Infrared Instrument (MIRI) observed the planetary nebula SMP LMC 058 in the Large Magellanic Cloud. The MRS was designed to provide medium resolution (R = $\lambda$/$\Delta\lambda$) 3D spectroscopy in the whole MIRI range. SMP LMC 058 is the only source observed in {\em JWST} commissioning that is both spatially and spectrally unresolved by the MRS and is a good test of {\em JWST's} capabilities. The new MRS spectra reveal a wealth of emission lines not previously detected in this planetary nebula. From these lines, the spectral resolving power ($\lambda$/$\Delta\lambda$) of the MRS is confirmed to be in the range R $=$ 4000 to 1500, depending on the MRS spectral sub-band. In addition, the spectra confirm that the carbon-rich dust emission is from SiC grains and that there is little to no time evolution of the SiC dust and emission line strengths over a 17-year epoch. These commissioning data reveal the great potential of the MIRI MRS for the study of circumstellar and interstellar material., Comment: 12 pages, 7 figures, 3 tables, accepted MNRAS

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

Using MIRI on-board JWST we present mid-infrared sub-arcsec imaging (MIRIM) and spectroscopy (MRS) of the hyperluminous infrared system SPT0311-58 at z=6.9. MIRI observations are compared with existing ALMA far-infrared continuum and [CII]158$\mu$m imaging. Even though the ALMA observations suggests very high star formation rates (SFR) in the eastern (E) and western (W) galaxies of the system, the H$\alpha$ line is not detected. This, together with the detection of the Pa$\alpha$ line, implies very high optical nebular extinction with lower limits of 4.2 (E) and 3.9 mag (W), and even larger 5.6 (E) and 10.0 (W) for SED derived values. The extinction-corrected Pa$\alpha$ SFRs are 383 and 230M$_{\odot}$yr$^{-1}$ for the E and W galaxies, respectively. This represents 50% of the SFRs derived from the [CII]158$\mu$m line and infrared light for the E galaxy and as low as 6% for the W galaxy. The MIRI observations reveal a clumpy stellar structure, with each clump having 3 to 5 $\times$10$^{9}$M$_\mathrm{\odot}$, leading to a total stellar mass of 2.0 and 1.5$\times$10$^{10}$M$_\mathrm{\odot}$ for the E and W galaxies, respectively. The specific SFR in the stellar clumps ranges from 25 to 59Gyr$^{-1}$, which are 3 to 10 times larger than the values measured in galaxies of similar mass at redshifts 6 to 8. The overall gas mass fraction is $M_\mathrm{gas}$/$M_*\sim3$, similar to that of z=4.5-6 star-forming galaxies. The observed properties of SPT0311-58 such as the clumpy distribution at sub(kpc) scales and the very high average extinction are similar to those observed in low- and intermediate-z LIRGs and ULIRGs, even though SPT0311-58 is observed only 800 Myr after the Big Bang. Massive, heavily obscured, clumpy starburst systems like SPT0311-58 likely represent the early phases in the formation of massive high-z bulge/spheroids and luminous quasars., Comment: 18 pages, 8 figures, final published version in A&A

The MRS is one of the four observing modes of JWST/MIRI. Using JWST in-flight data of unresolved (point) sources, we can derive the MRS absolute spectral response function (ASRF) starting from raw data. Spectral fringing plays a critical role in the derivation and interpretation of the MRS ASRF. In this paper, we present an alternative way to calibrate the data. Firstly, we aim to derive a fringe correction that accounts for the dependence of the fringe properties on the MIRI pupil illumination and detector pixel sampling of the point spread function. Secondly, we aim to derive the MRS ASRF using an absolute flux calibrator observed across the full 5 to 28 $\mu$m wavelength range of the MRS. Thirdly, we aim to apply the new ASRF to the spectrum of a G dwarf and compare with the output of the JWST/MIRI default data reduction pipeline. Finally, we examine the impact of the different fringe corrections on the detectability of molecular features in the G dwarf and K giant. The absolute flux calibrator HD 163466 (A-star) is used to derive tailored point source fringe flats at each of the default dither locations of the MRS. The fringe-corrected point source integrated spectrum of HD 163466 is used to derive the MRS ASRF using a theoretical model for the stellar continuum. A cross-correlation is run to quantify the uncertainty on the detection of CO, SiO, and OH in the K giant and CO in the G dwarf for different fringe corrections. The point-source-tailored fringe correction and ASRF are found to perform at the same level as the current corrections, beating down the fringe contrast to the sub-percent level, whilst mitigating the alteration of real molecular features. The same tailored solutions can be applied to other MRS unresolved targets. A pointing repeatability issue in the MRS limits the effectiveness of the tailored fringe flats is at short wavelengths., Comment: 17 pages, 15 figures, updated to accepted version (results did not change)

We present the discovery and analysis of J1316+2614 at z=3.6130, a UV-bright star-forming galaxy ($M_{\rm UV} \simeq -24.7$) with large escape of Lyman continuum (LyC) radiation. J1316+2614 is a young ($\simeq 10$ Myr) star-forming galaxy with $SFR \simeq 500 M_{\odot}$ yr$^{-1}$ and a starburst mass of log($M_{\star}/M_{\odot}) \simeq 9.7$. It shows a very steep UV continuum, $\beta_{\rm UV} \simeq -2.59 \pm 0.05$, consistent with residual dust obscuration, $E(B-V)\simeq 0$. LyC emission is detected with high significance ($\simeq 17 \sigma$) down to $830$\r{A}, for which a very high relative (absolute) LyC escape fraction $f_{\rm esc} \rm (LyC) \simeq 0.92$ ($\simeq 0.87$) is inferred. The contribution of a foreground or AGN contamination to the LyC signal is discussed but is unlikely. J1316$+$2614 is the most powerful ionizing source known among the star-forming galaxy population, both in terms of production ($Q_{\rm H} \approx 10^{56}$ s$^{-1}$) and escape of ionizing photons ($f_{\rm esc} \rm (LyC) \approx 0.9$). Nebular emission in Ly$\alpha$, H$\beta$, and other rest-frame optical lines are detected, but these are weak ($EW_{0} \rm [H\beta] \simeq 35$\r{A}), with their strengths reduced roughly by $\simeq 90\%$. J1316+2614 is the first case known where the effect of large escape of ionizing photons on the strength of nebular lines and continuum emission is clearly observed. Gas inflows are detected in J1316+2614 from the blue-dominated peak Ly$\alpha$ emission (with a blue-to-red peak line ratio $I_{\rm blue}/I_{\rm red} \simeq 3.7$) and redshifted ISM absorption ($\simeq 100$ km s$^{-1}$). Our results suggest that J1316+2614 is undergoing a gas compaction event, possibly representing a short-lived phase in the evolution of massive and compact galaxies, where strong gas inflows have triggered an extreme star formation episode and nearly $100\%$ LyC photons are escaping., Comment: Accepted for publication in MNRAS. 18 pages, 10 figures, 5 tables

During the commissioning of the James Webb Space Telescope (JWST), the mid-infrared instrument (MIRI) observed NGC6552 with the MIRI Imager and the medium-resolution spectrograph (MRS). NGC6552 is an active galactic nucleus (AGN) at redshift 0.0266 classified as a Seyfert 2 nucleus in the optical, and Compton-thick AGN in X-rays. This work exemplifies and demonstrates the MRS capabilities to study the mid-infrared (mid-IR) spectra and characterize the physical conditions and kinematics of the ionized and molecular gas in the nuclear regions of nearby galaxies. We obtained the nuclear, circumnuclear, and central mid-IR spectra of NGC6552. They provide the first clear observational evidence for a nuclear outflow in NGC6552. The outflow contributes to 67$\pm$7% of the total line flux independent of the ionization potential (27 to 187 eV) and critical densities (10$^4$ to 4$\times$10$^{6}$ cm$^{-3}$), showing an average blue-shifted peak velocity of -127$\pm$45 kms$^{-1}$ and an outflow maximal velocity of 698$\pm$80 kms$^{-1}$. Since the mid-IR photons penetrate dusty regions as efficiently as X-ray keV photons, we interpret these results as the evidence for a highly ionized, non-stratified, AGN-powered, and fast outflowing gas in a low density environment (few 10$^{3}$ cm$^{-3}$) located very close (<0.2kpc) to the Compton-thick AGN. Nine pure rotational molecular Hydrogen lines are detected and spectrally resolved, and exhibit symmetric Gaussian profiles, consistent with the galactic rotation, and with no evidence of outflowing H$_{2}$ material. We detect a warm H$_{2}$ mass of $1.9\pm1.1\times10^7 M_{\odot}$ in the central region (1.8 kpc in diameter) of the galaxy, with almost 30% of that mass in the circum-nuclear region. Line ratios confirm that NGC6552 has a Seyfert nucleus with a black hole mass estimated in the range of 0.6 to 6 million solar masses., Comment: 13 pages, 5 figures, 5 tables, accepted in A&A

We present JWST/MIRI MRS spectroscopy of NGC7319, the largest galaxy in the Stephan's Quintet, observed as part of the Early Release Observations (ERO). NGC7319 hosts a type 2 active galactic nucleus (AGN) and a low-power radio jet (L_1.4GHz=3.3x10^22 W Hz^-1) with two asymmetric radio hotspots at 430 pc (N2) and 1.5 kpc (S2) projected distances from the unresolved radio core. The MRS data suggest that the molecular material in the disk of the galaxy decelerates the jet and causes this length asymmetry. We find enhanced emission from warm and hot H_2 (T_w=330+-40 K, T_h=900+-60 K) and ionized gas at the intersection between the jet axis and dust lanes in the disk. This emission is coincident with the radio hotspot N2, the hotspot closer to the core, suggesting that the jet-interstellar medium (ISM) interaction decelerates the jet. Conversely, the mid-infrared emission at the more distant hotspot is fainter, more highly ionized, and with lower H_2 excitation, suggesting a more diffuse atomic environment where the jet can progress to farther distances. At the N2 radio hotspot, the ionized gas mass (M_ion=(2.4-12)x10^5 Msun) is comparable to that of the warm H_2, but the former is more turbulent (sigma_ion~300 vs. sigma_H2~150 km/s), so the mechanical energy of the ionized gas is ~1.3-10 times higher. From these estimates, we find that only <1% of the jet energy remains as mechanical energy in these two ISM phases at N2. We also find extended (r>0.3-1.5 kpc) high-ionization emission ([MgV], [NeVI], and [NeV]) close to the radio hotspots. This initial analysis of NGC7319 shows the potential of MIRI/MRS to investigate the AGN feedback mechanisms due to radio jets and their radiation field in the, often heavily dust-enshrouded, central regions of galaxies. Understanding these mechanisms is an essential ingredient in the development of cosmological simulations of galaxy evolution., Comment: 10 pages, 8 figures. Small changes to match the accepted for publication version in A&A

The Mid-Infrared Instrument (MIRI) on-board JWST will provide imaging, coronagraphy, low-resolution spectroscopy and medium-resolution spectroscopy at unprecedented sensitivity levels in the mid-infrared wavelength range. The Medium-Resolution Spectrometer (MRS) of MIRI is an integral field spectrograph that provides diffraction-limited spectroscopy between 4.9 and 28.3 um, within a FOV varying from 13 to 56" square. From ground testing, we calculate the physical parameters essential to general observers and calibrating the wavelength solution and resolving power of the MRS is critical for maximising the scientific performance of the instrument. We have used ground-based observations of discrete spectral features in combination with Fabry-Perot etalon spectra to characterize the wavelength solution and spectral resolving power of the MRS. We present the methodology used to derive the MRS spectral characterisation, which includes the precise wavelength coverage of each MRS sub-band, computation of the resolving power as a function of wavelength, and measuring slice-dependent spectral distortions. The resolving power varies from R3500 in channel 1 to R1500 in channel 4. Based on the ground test data, the wavelength calibration accuracy is estimated to be below one tenth of a pixel, with small systematic shifts due to the target position within a slice for unresolved sources, that have a maximum amplitude of about 0.25 spectral resolution elements. Based on ground test data, the MRS complies with the spectral requirements for both the R and wavelength accuracy for which it was designed. We also present the commissioning strategies and targets that will be followed to update the spectral characterisation of the MRS., Comment: 13 pages. Accepted for publication in A&A

BOSS-EUVLG1 is the most ultraviolet (UV) and Ly$\alpha$ luminous galaxy detected so far in the Universe, going through a very active starburst phase, and forming stars at a rate (SFR) of 955 $\pm$ 118 M$_{\odot}$ yr$^{-1}$. We report the detection of a broad H$\alpha$ component carrying 25% of the total H$\alpha$ flux. The broad H$\alpha$ line traces a fast and massive ionized gas outflow characterized by a total mass, $\log(M_{out}[M_{\odot}]),$ of 7.94 $\pm$ 0.15, an outflowing velocity (V$_{out}$) of 573 $\pm$ 151 km s$^{-1}$, and an outflowing mass rate ($\dot{M}_{out}$) of 44 $\pm$ 20 M$_{\odot}$ yr$^{-1}$. The presence of the outflow in BOSS-EUVLG1 is also supported by the identification of blueshifted UV absorption lines in low and high ionization states. The energy involved in the H$\alpha$ outflow can be explained by the ongoing star formation without the need for an Active Galactic Nucleus. The derived low mass loading factor ($\eta$= 0.05 $\pm$ 0.03) indicates that although massive, this phase of the outflow can not be relevant for the quenching of the star formation. In addition, only a small fraction ($\leq$ 15%) of the ionized outflowing material with velocities above 372 km s$^{-1}$ could escape the gravitational potential, and enrich the surrounding circum-galactic medium at distances above tens of kpc. The ionized phase of the outflow does not carry the mass and energy to play a relevant role neither in the evolution of the host galaxy nor in the enrichment of the intergalactic medium. Other phases of the outflow could be carrying most of the outflow energy and mass in the form of hot X-ray emitting gas as predicted by some recent simulations. The expected emission of the extended X-ray emitting halo associated with the outflow in BOSS-EUVLG1 and similar galaxies could be detected with the future X-ray observatory, {\it ATHENA} but could not be resolved spatially., Comment: 8 pages, 3 figures, accepted in A&A

We report the discovery of BOSS-EUVLG1 at z=2.469, by far the most luminous, almost un-obscured star-forming galaxy known at any redshift. First classified as a QSO within the Baryon Oscillation Spectroscopic Survey, follow-up observations with the Gran Telescopio Canarias reveal that its large luminosity, MUV = -24.40 and log(L_Lya/erg s-1) = 44.0, is due to an intense burst of star-formation, and not to an AGN or gravitational lensing. BOSS-EUVLG1 is a compact (reff = 1.2 kpc), young (4-5 Myr) starburst with a stellar mass log(M*/Msun) = 10.0 +/- 0.1 and a prodigious star formation rate of ~1000 Msun yr-1. However, it is metal- and dust-poor (12+log(O/H) = 8.13 +/- 0.19, E(B-V) = 0.07, log(LIR/LUV) < -1.2), indicating that we are witnessing the very early phase of an intense starburst that has had no time to enrich the ISM. BOSS-EUVLG1 might represent a short-lived (<100 Myrs), yet important phase of star-forming galaxies at high redshift that has been missed in previous surveys. Within a galaxy evolutionary scheme, BOSS-EUVLG1 could likely represent the very initial phases in the evolution of massive quiescent galaxies, even before the dusty star-forming phase., Comment: 6 pages, 3 figures, 1 table. Accepted for publication in MNRAS Letters

We present deep rest-frame UV spectroscopic observations using the Gran Telescopio Canarias of six gravitationally lensed Lya emitters (LAEs) at $2.36

This work explores from a statistical point of view the rest-frame Far-ultraviolet (FUV) to Far-infrared (FIR) emission of a population of Lyman break galaxies (LBGs) at $z\sim3$ that cannot be individually detected from current FIR observations. We perform a stacking analysis over a sample of $\sim$17000 LBGs at redshift $2.5 10.25$)., Comment: 20 pages, 12 figures, Accepted for publication in A&A

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

$\mathit{Herschel}$ extragalactic surveys offer a unique opportunity to efficiently select a significant number of rare and massive dusty objects, and thus gain insight into the prodigious star-forming activity that takes place in the very distant Universe. To search for $z\geq4$ dusty star-forming galaxies, in this work we consider red SPIRE objects with fluxes rising from 250 $\mu$m to $500\:\mu$m (so-called "500 $\mu$m-risers"). We aim to implement a novel method to obtain a statistical sample of "500 $\mu$m-risers" and fully evaluate our selection inspecting different models of galaxy evolution. We consider one of the largest and deepest ${\it Herschel}$ surveys, the Herschel Virgo Cluster Survey. We develop a novel selection algorithm which links the source extraction and spectral energy distribution fitting. We select 133 "500 $\mu$m-risers" over 55 deg$^{2}$, imposing the criteria: $S_{500}>S_{350}>S_{250}$, $S_{250}>13.2$ mJy and $S_{500}>$30 mJy. Differential number counts are in a fairly good agreement with models, displaying better match than other existing samples. In order to interpret the statistical properties of selected sources, which has been proven as a very challenging task due the complexity of observed artefacts, we make end-to-end simulations including physical clustering and lensing. The estimated fraction of strongly lensed sources is $24^{+6}_{-5}\%$ based on models. We present the faintest known statistical sample of "500 $\mu$m-risers" and show that noise and strong lensing have crucial impact on measured counts and redshift distribution of selected sources. We estimate the flux-corrected star formation rate density at $44$ sources., Comment: 28 pages, 13 figures, version accepted for publication in A&A

In this work we investigate the far-UV to NIR shape of the dust attenuation curve of a sample of IR selected dust obscured (U)LIRGs at z$\sim$2. The spectral energy distributions (SEDs) are fitted with CIGALE, a physically-motivated spectral synthesis model based on energy balance. Its flexibility allows us to test a wide range of different analytical prescriptions for the dust attenuation curve, including the well-known Calzetti and Charlot & Fall curves, and modified versions of them. The attenuation curves computed under the assumption of our reference double power-law model are in very good agreement with those derived, in previous works, with radiative transfer (RT) SED fitting. We investigate the position of our galaxies in the IRX-$\beta$ diagram and find this to be consistent with grayer slopes, on average, in the UV. We also find evidence for a flattening of the attenuation curve in the NIR with respect to more classical Calzetti-like recipes. This larger NIR attenuation yields larger derived stellar masses from SED fitting, by a median factor of $\sim$ 1.4 and up to a factor $\sim$10 for the most extreme cases. The star formation rate appears instead to be more dependent on the total amount of attenuation in the galaxy. Our analysis highlights the need for a flexible attenuation curve when reproducing the physical properties of a large variety of objects., Comment: 20 pages, 14 figures

We present \ci\,(2--1) and multi-transition $^{12}$CO observations of a dusty star-forming galaxy, ACT\,J2029+0120, which we spectroscopically confirm to lie at $z$\,=\,2.64. We detect CO(3--2), CO(5--4), CO(7--6), CO(8--7), and \ci\,(2--1) at high significance, tentatively detect HCO$^{+}$(4--3), and place strong upper limits on the integrated strength of dense gas tracers (HCN(4--3) and CS(7--6)). Multi-transition CO observations and dense gas tracers can provide valuable constraints on the molecular gas content and excitation conditions in high-redshift galaxies. We therefore use this unique data set to construct a CO spectral line energy distribution (SLED) of the source, which is most consistent with that of a ULIRG/Seyfert or QSO host object in the taxonomy of the \textit{Herschel} Comprehensive ULIRG Emission Survey. We employ RADEX models to fit the peak of the CO SLED, inferring a temperature of T$\sim$117 K and $n_{\text{H}_2}\sim10^5$ cm$^{-3}$, most consistent with a ULIRG/QSO object and the presence of high density tracers. We also find that the velocity width of the \ci\ line is potentially larger than seen in all CO transitions for this object, and that the $L'_{\rm C\,I(2-1)}/L'_{\rm CO(3-2)}$ ratio is also larger than seen in other lensed and unlensed submillimeter galaxies and QSO hosts; if confirmed, this anomaly could be an effect of differential lensing of a shocked molecular outflow., Comment: Accepted for publication in ApJ

We explore from a statistical point of view the far-infrared (far-IR) and sub-millimeter (sub-mm) properties of a large sample of LBGs (22,000) at z~3 in the COSMOS field. The large number of galaxies allows us to split it in several bins as a function of UV luminosity, UV slope, and stellar mass to better sample their variety. We perform stacking analysis in PACS (100 and 160 um), SPIRE (250, 350 and 500 um) and AzTEC (1.1 mm) images. Our stacking procedure corrects the biases induced by galaxy clustering and incompleteness of our input catalogue in dense regions. We obtain the full IR spectral energy distributions (SED) of subsamples of LBGs and derive the mean IR luminosity as a function of UV luminosity, UV slope, and stellar mass. The average IRX is roughly constant over the UV luminosity range, with a mean of 7.9 (1.8 mag). However, it is correlated with UV slope, and stellar mass. We investigate using a statistically-controlled stacking analysis as a function of (stellar mass, UV slope) the dispersion of the IRX-UVslope and IRX-M* plane. Our results enable us to study the average relation between star-formation rate (SFR) and stellar mass, and we show that our LBG sample lies on the main sequence of star formation at z~3., Comment: Accepted to A&A, 17 Pages, 14 Figures, 2 Tables

We report contributions to cosmic infrared background (CIB) intensities originating from known galaxies and their faint companions at submillimeter wavelengths. Using the publicly-available UltraVISTA catalog, and maps at 250, 350, and 500 {\mu}m from the \emph{Herschel} Multi-tiered Extragalactic Survey (HerMES), we perform a novel measurement that exploits the fact that uncatalogued sources may bias stacked flux densities --- particularly if the resolution of the image is poor --- and intentionally smooth the images before stacking and summing intensities. By smoothing the maps we are capturing the contribution of faint (undetected in K_S ~ 23.4) sources that are physically associated, or correlated, with the detected sources. We find that the cumulative CIB increases with increased smoothing, reaching 9.82 +- 0.78, 5.77 +- 0.43, and 2.32 +- 0.19$\, \rm nW m^{-2} sr^{-1}$ at 250, 350, and 500 {\mu}m at 300 arcsec FWHM. This corresponds to a fraction of the fiducial CIB of 0.94 +- 0.23, 1.07 +- 0.31, and 0.97 +- 0.26 at 250, 350, and 500 {\mu}m, where the uncertainties are dominated by those of the absolute CIB. We then propose, with a simple model combining parametric descriptions for stacked flux densities and stellar mass functions, that emission from galaxies with log(M/Msun) > 8.5 can account for the most of the measured total intensities, and argue against contributions from extended, diffuse emission. Finally, we discuss prospects for future survey instruments to improve the estimates of the absolute CIB levels, and observe any potentially remaining emission at z > 4., Comment: Accepted to ApJL. 6 Pages, 3 figures