Search

Using more than 100 galaxies in the MUSE Ultra Deep Field with spectroscopy from the Hubble Space Telescope's (HST) Wide Field Camera 3 and the Very Large Telescope's Multi Unit Spectroscopic Explorer, we extend the gas-phase mass-metallicity relation (MZR) at z approximate to 1-2 down to stellar masses of M-star approximate to 10(7.5)M(circle dot). The sample reaches 6 times lower in stellar mass and star formation rate (SFR) than previous HST studies at these redshifts, and we find that galaxy metallicities decrease to log(O/H) + 12 approximate to 7.8 +/- 0.1 (15% solar) at log(M-star/M-circle dot) approximate to 7.5, without evidence of a turnover in the shape of the MZR at low masses. We validate our strong-line metallicities using the direct method for sources with [O iii] lambda 4363 and [O iii] lambda 1666 detections, and find excellent agreement between the techniques. The [O iii] lambda 1666-based metallicities double existing measurements with a signal-to-noise ratio >= 5 for unlensed sources at z > 1, validating the strong-line calibrations up to z similar to 2.5. We confirm that the MZR resides similar to 0.3 dex lower in metallicity than local galaxies and is consistent with the fundamental metallicity relation if the low-mass slope varies with SFR. At lower redshifts (z similar to 0.5) our sample reaches similar to 0.5 dex lower in SFR than current calibrations and we find enhanced metallicities that are consistent with extrapolating the MZR to lower SFRs. Finally, we detect only an similar to 0.1 dex difference in the metallicities of galaxies in groups versus isolated environments. These results are based on robust calibrations and reach the lowest masses and SFRs that are accessible with HST, providing a critical foundation for studies with the Webb and Roman Space Telescopes.

One of the surprising early findings with JWST has been the discovery of a strong "roll-over" or a softening of the absorption edge of Ly$\alpha$ in a large number of galaxies at ($z\gtrsim 6$), in addition to systematic offsets from photometric redshift estimates and fundamental galaxy scaling relations. This has been interpreted as damped Ly$\alpha$ absorption (DLA) wings from high column densities of neutral atomic hydrogen (HI), signifying major gas accretion events in the formation of these galaxies. To explore this new phenomenon systematically, we assemble the JWST/NIRSpec PRImordial gas Mass AssembLy (PRIMAL) legacy survey of 494 galaxies at $z=5.5-13.4$. We characterize this benchmark sample in full and spectroscopically derive the galaxy redshifts, metallicities, star-formation rates, and ultraviolet slopes. We define a new diagnostic, the Ly$\alpha$ damping parameter $D_{\rm Ly\alpha}$ to measure and quantify the Ly$\alpha$ emission strength, HI fraction in the IGM, or local HI column density for each source. The JWST-PRIMAL survey is based on the spectroscopic DAWN JWST Archive (DJA-Spec). All the software, reduced spectra, and spectroscopically derived quantities and catalogs are made publicly available in dedicated repositories. The fraction of strong galaxy DLAs are found to be in the range $65-95\%$ at $z>5.5$. The fraction of strong Ly$\alpha$ emitters (LAEs) is found to increase with decreasing redshift, in qualitative agreement with previous observational results, and are predominantly associated with low-metallicity and UV faint galaxies. By contrast, strong DLAs are observed in galaxies with a variety of intrinsic physical properties. Our results indicate that strong DLAs likely reflect a particular early assembly phase of reionization-era galaxies, at which point they are largely dominated by pristine HI gas accretion. [abridged], Comment: 20 pages, 18 figures, 3 tables. Submitted to A&A. Comments welcome! All data and catalogs are available through the DAWN JWST Archive (DJA): https://dawn-cph.github.io/dja/ and https://github.com/keheintz/jwst-primal

Stellar populations and their distribution differ widely across the Galaxy, which is likely to affect planet demographics. Our local neighbourhood is dominated by young, metal-rich stars in the galactic thin disc, while the stellar halo and galactic bulge host a large fraction of older, metal-poor stars. We study the impact of these variations on planet populations in different regions of the Galaxy by combining a high-resolution galaxy formation simulation with state-of-the-art planet population synthesis models. We construct a population model to estimate occurrence rates of different planet types, based on the New Generation Planet Population Synthesis by Emsenhuber et al., 2021. We apply this to a simulated Milky Way--Analogue in the HESTIA galaxy formation simulation. We study the planet occurrence rate in the metal-rich regions of the inner Galaxy, i.e. in the galactic bulge and thin disc, and contrast them to the frequencies in the more distant, metal-poor region like the thick disc and stellar halo. We find that the planet demographics in the metal-poor regions of the Milky Way-Analogue, differ strongly from the planet populations in the more distant, metal-poor regions. The occurrence rate of giant planets ($>300 M_\bigoplus$) is 10 to 20 times larger in the thin disc compared to the thick disc, driven by the low amounts of solid material available for planet formation around metal-poor stars. Similarly, low-mass Earth-like planets around Sun-like stars are most abundant in the thick disc, being 1.5 times more frequent than in the thin disc. Moreover, low-mass planets are expected to be abundant throughout the galaxy, from the central regions to the outer halo, due to their formation processes being less dependent on stellar metallicity. The planet populations differ more strongly around Sun-like stars compared to dwarfs with masses 0.3 - 0.5 $M_\odot$., Comment: 19 pages, 13 figures, accepted by A&A. Code available at https://github.com/ChrisBoettner/skaro

The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational-wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational-wave observations by LISA to probe the universe.

The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA’s first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultra-compact stellar-mass binaries, massive black hole binaries, and extreme or interme-diate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.

We present a study of the metal-enriched halo gas, traced using Mg II and [O II] emission lines, in two large, blind galaxy surveys – the MUSE (Multi Unit Spectroscopic Explorer) Analysis of Gas around Galaxies (MAGG) and the MUSE Ultra Deep Field (MUDF). By stacking a sample of ≈600 galaxies (stellar masses M∗ ≈ 106–12 M☉), we characterize for the first time the average metal line emission from a general population of galaxy haloes at 0.7 ≤ z ≤ 1.5. The Mg II and [O II] line emission extends farther out than the stellar continuum emission, on average out to ≈25 and ≈45 kpc, respectively, at a surface brightness (SB) level of 10-20 erg s-1 cm-2 arcsec-2. The radial profile of the Mg II SB is shallower than that of the [O II], suggesting that the resonant Mg II emission is affected by dust and radiative transfer effects. The [O II] to Mg II SB ratio is ≈3 over ≈20–40 kpc, also indicating a significant in situ origin of the extended metal emission. The average SB profiles are intrinsically brighter by a factor of ≈2–3 and more radially extended by a factor of ≈1.3 at 1.0 < z ≤ 1.5 than at 0.7 ≤ z ≤ 1.0. The average extent of the metal emission also increases independently with increasing stellar mass and in overdense group environments. When considering individual detections, we find extended [O II] emission up to ≈50 kpc around ≈30–40 per cent of the group galaxies, and extended (≈30–40 kpc) Mg II emission around two z ≈ 1 quasars in groups, which could arise from outflows or environmental processes.

We investigate the degree of dust obscured star formation in 49 massive (${\rm log}_{10}(M_{\star}/{\rm M}_{\odot})>9$) Lyman-break galaxies (LBGs) at $z = 6.5$-$8$ observed as part of the ALMA Reionization Era Bright Emission Line Survey (REBELS) large program. By creating deep stacks of the photometric data and the REBELS ALMA measurements we determine the average rest-frame UV, optical and far-infrared (FIR) properties which reveal a significant fraction ($f_{\rm obs} = 0.4$-$0.7$) of obscured star formation, consistent with previous studies. From measurements of the rest-frame UV slope, we find that the brightest LBGs at these redshifts show bluer ($\beta \simeq -2.2$) colours than expected from an extrapolation of the colour-magnitude relation found at fainter magnitudes. Assuming a modified blackbody spectral-energy distribution (SED) in the FIR (with dust temperature of $T_{\rm d} = 46\,{\rm K}$ and $\beta_{\rm d} = 2.0$), we find that the REBELS sources are in agreement with the local ''Calzetti-like'' starburst Infrared-excess (IRX)-$\beta$ relation. By reanalysing the data available for 108 galaxies at $z \simeq 4$-$6$ from the ALPINE ALMA large program using a consistent methodology and assumed FIR SED, we show that from $z \simeq 4$-$8$, massive galaxies selected in the rest-frame UV have no appreciable evolution in their derived IRX-$\beta$ relation. When comparing the IRX-$M_{\star}$ relation derived from the combined ALPINE and REBELS sample to relations established at $z < 4$, we find a deficit in the IRX, indicating that at $z > 4$ the proportion of obscured star formation is lower by a factor of $\gtrsim 3$ at a given a $M_{\star}$. Our IRX-$\beta$ results are in good agreement with the high-redshift predictions of simulations and semi-analytic models for $z \simeq 7$ galaxies with similar stellar masses and SFRs., Comment: 18 pages, 10 figures, 2 tables (plus 1 figure and 2 tables in the appendix). Updated to match MNRAS accepted version after minor corrections

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.07\pm0.20$ 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.49^{+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\%^{+2}_{-1}$, which results in $\rm log_{10}(\xi_{ion}/(Hz\;erg^{-1})) = 25.52^{+0.09}_{-0.08}$. 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$-$14\% 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.59 \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. Revised version submitted in ApJ

A key to understanding the formation of the first galaxies is to quantify the content of the molecular gas as the fuel for star formation activity through the epoch of reionization. In this paper, we use the 158$\mu$m [CII] fine-structure emission line as a tracer of the molecular gas in the interstellar medium (ISM) in a sample of $z=6.5-7.5$ galaxies recently unveiled by the Reionization Era Bright Line Emission Survey, REBELS, with the Atacama Large Millimeter/submillimeter Array. We find substantial amounts of molecular gas ($\sim10^{10.5}\ M_\odot$) comparable to those found in lower redshift galaxies for similar stellar masses ($\sim10^{10}\ M_\odot$). The REBELS galaxies appear to follow the standard scaling relations of molecular gas to stellar mass ratio ($\mu_{\rm mol}$) and gas depletion timescale ($t_{\rm dep}$) with distance to the star-forming main-sequence expected from extrapolations of $z\sim1-4$ observations. We find median values at $z\sim7$ of $\mu_{\rm mol}=2.6_{-1.4}^{4.1}$ and $t_{\rm dep}=0.5_{-0.14}^{+0.26}$ Gyr, indicating that the baryonic content of these galaxies is gas-phase dominated and little evolution from $z\sim7$ to 4. Our measurements of the cosmic density of molecular gas, log$(\rho_{\rm mol}/(M_\odot {\rm Mpc}^{-3}))=6.34^{+0.34}_{-0.31}$, indicate a steady increase by an order of magnitude from $z\sim7$ to 4., Comment: Accepted for publication in Astronomy & Astrophysics. 11 pages, 7 figures. Typo fixed in an author's name

We present Atacama Large Millimeter/Submillimeter Array (ALMA) [CII] and $\sim158$ $\rm\mu m$ continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; $L_{\rm IR}=1.5^{+0.8}_{-0.5}\times10^{12}$ L$_\odot$) at $z=7.31$, spectroscopically confirmed by the Reionization Era Bright Emission Line Survey (REBELS) ALMA Large Programme. REBELS-25 has a significant stellar mass of $M_{*}=8^{+4}_{-2}\times10^{9}$ M$_\odot$. From dust-continuum and ultraviolet observations, we determine a total obscured + unobscured star formation rate of SFR $=199^{+101}_{-63}$ M$_\odot$ yr$^{-1}$. This is about four times the SFR estimated from an extrapolated main-sequence. We also infer a [CII]-based molecular gas mass of $M_{\rm H_2}=5.1^{+5.1}_{-2.6}\times10^{10}$ $M_\odot$, implying a molecular gas depletion time of $ t_{\rm depl, H_2}=0.3^{+0.3}_{-0.2}$ Gyr. We observe a [CII] velocity gradient consistent with disc rotation, but given the current resolution we cannot rule out a more complex velocity structure such as a merger. The spectrum exhibits excess [CII] emission at large positive velocities ($\sim500$ km s$^{-1}$), which we interpret as either a merging companion or an outflow. In the outflow scenario, we derive a lower limit of the mass outflow rate of 200 M$_\odot$ yr$^{-1}$, which is consistent with expectations for a star formation-driven outflow. Given its large stellar mass, SFR and molecular gas reservoir $\sim700$ Myr after the Big Bang, we explore the future evolution of REBELS-25. Considering a simple, conservative model assuming an exponentially declining star formation history, constant star formation efficiency, and no additional gas inflow, we find that REBELS-25 has the potential to evolve into a galaxy consistent with the properties of high-mass quiescent galaxies recently observed at $z\sim4$., Comment: Accepted for publication in MNRAS. 21 pages, 8 figures

We present the first observational infrared luminosity function (IRLF) measurement in the Epoch of Reionization (EoR) based on a UV-selected galaxy sample with ALMA spectroscopic observations. Our analysis is based on the ALMA large program Reionization Era Bright Emission Line Survey (REBELS), which targets 42 galaxies at $\mathrm{z=6.4-7.7}$ with [CII] 158$\micron$ line scans. 16 sources exhibit a dust detection, 15 of which are also spectroscopically confirmed through the [CII] line. The IR luminosities of the sample range from $\log L_{IR}/L_\odot=11.4$ to 12.2. Using the UVLF as a proxy to derive the effective volume for each of our target sources, we derive IRLF estimates, both for detections and for the full sample including IR luminosity upper limits. The resulting IRLFs are well reproduced by a Schechter function with the characteristic luminosity of $\log L_{*}/L_\odot=11.6^{+0.2}_{-0.1}$. Our observational results are in broad agreement with the average of predicted IRLFs from simulations at $z\sim7$. Conversely, our IRLFs lie significantly below lower redshift estimates, suggesting a rapid evolution from $z\sim4$ to $z\sim7$, into the reionization epoch. The inferred obscured contribution to the cosmic star-formation rate density at $z\sim7$ amounts to $\mathrm{log(SFRD/M_{\odot}/yr/Mpc^{3}) = -2.66^{+0.17}_{-0.14} }$ which is at least $\sim$10\% of UV-based estimates. We conclude that the presence of dust is already abundant in the EoR and discuss the possibility of unveiling larger samples of dusty galaxies with future ALMA and JWST observations., Comment: 9 pages, 5 figures

We present an early analysis on the search for high redshift galaxies using the deepest public JWST imaging to date, the NGDEEP field. This data consists of 6-band NIRCam imaging on the Hubble Ultra Deep Field-Par2, covering a total area of 6.3 arcmin$^{2}$. Based on our initial reduction of the first half of this survey, we reach 5$\sigma$ depths up to mag = 29.5--29.9 between $1-5$ um. Such depths present an unprecedented opportunity to begin exploring the early Universe with JWST. As such, we find high redshift galaxies in this field by examining the spectral energy distribution of these systems and present 18 new $z > 8$ systems identified using two different photometric redshift codes: LePhare and EAZY, combined with other significance criteria. The highest redshift object in our sample is at $z=15.57^{+0.39}_{-0.38}$, which has a blue beta slope of $\beta=-3.25^{+0.41}_{-0.46}$ and a very low inferred stellar mass of $M_{*} = 10^{7.39}$ M_0,. We also discover a series of faint, low-mass dwarf galaxies with $M_{*} < 10^{8.5}$ M_0 at $z \sim 9$ that have blue colors and UV slopes. The structure of these galaxies is such that they all have very flat surface brightness profiles and small sizes $< 1 \ \mathrm{kpc}$. We also compare our results to theory, finding no significant disagreement with some CDM based models.The discovery of these objects, most of which are low luminosity and inferred stellar mass, demonstrates the power of probing continuously deeper into the Universe, pointing the way to deeper, or similar depth but wider area, surveys and demonstrate the critical need for JWST deep fields to explore this aspect of the early Universe., Comment: ApJL submitted. Added acknowledgements and slight change in title. Mosaics of our NGDEEP reductions are available here: tinyurl.com/yckcfnxu

We investigate the ionizing properties of the pair of bright Ly$\alpha$ emitting galaxies BDF521 and BDF2195 at z=7.012 in order to constrain their contribution to the formation of the BDF "reionized bubble". We obtain constraints on UV emission lines (CIV$\lambda 1548$ doublet, HeII$\lambda 1640$, OIII]$\lambda 1660$ doublet, and CIII]$\lambda 1909$ doublet) from deep VLT-XSHOOTER observations and compare them to those available for other high-redshift objects, and to models with mixed stellar and AGN emission. We use this spectroscopic information together with the photometry available in the field to constrain the physical properties of the two objects using the spectro-photometric fitting code BEAGLE. We do not detect any significant emission at the expected position of UV lines, with 3$\sigma$ upper limits of EW$\lesssim$2-7 AA rest-frame. We find that the two objects have lower CIII] emission than expected on the basis of the correlation between the Ly$\alpha$ and CIII] EWs. The EW limits on CIV and HeII emission exclude pure AGN templates at $\sim2-3\sigma$ significance, and only models with a $\lesssim$40% AGN contribution are compatible with the observations. The two objects are found to be relatively young ($\sim$20-30 Myrs) and metal-poor ($\lesssim 0.3 Z_{\odot}$) with stellar masses of a few $10^9M_{\odot}$. Their production rate of hydrogen ionizing photons per intrinsic UV luminosity is log($\xi_{ion}^*$/Hz erg$^{-1}$)=25.02-25.26, consistent with values typically found in high-redshift galaxies, but more than twice lower than values measured in $z>$7 galaxies with strong CIII] and/or optical line emission ($\simeq$25.6-25.7). The two BDF emitters have no evidence of higher than average ionizing capabilities and are not capable of reionizing their surroundings by their own means under realistic assumptions on the escape fraction of ionizing photons. (abridged), Comment: 9 pages, 4 figures, submitted to A&A

We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman Break Galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at $z\simeq 7$. We find that the galaxy dust mass, $M_d$, is regulated by the combination of SNII dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard timescale $\tau_0= 30$ Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of $\sim 0.07-0.1\%$ and a UV-to-total star formation rate relation such that $log (\psi_{\rm UV}) = -0.05 ~[log (\psi)]^{2} + 0.86 ~log(\psi) -0.05$ (implying that 55-80\% of the star formation is obscured) for REBELS galaxies with stellar mass $M_* = 10^{9-10} M_\odot$. This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at $z=7$. However, 2 out of the 13 systems show dust-to-stellar mass ratios ($\sim 0.94-1.1\%$) that are up to $18\times$ larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing $\tau_0$ to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor $ \sim 2$. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions., Comment: Accepted to MNRAS

We analyse FIR dust continuum measurements for 14 galaxies ($z\approx 7$) in the ALMA REBELS LP to derive their physical properties. Our model uses three input data: (a) the UV spectral slope, $\beta$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158\mu$m, $F_{158}$, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have (28-90.5)% of their star formation obscured; the total (UV+IR) star formation rates are in the range $31.5 < {\rm SFR}/ (M_\odot {\rm yr}^{-1}) < 129.5$. The sample-averaged dust mass and temperature are $(1.3\pm 1.1)\times 10^7 M_\odot$ and $52 \pm 11$ K, respectively. In some galaxies dust is abundant (REBELS-14, $M'_d \approx 3.4 \times 10^7 M_\odot$), or hot (REBELS-18, $T'_d \approx 67$ K). The dust distribution is compact ($<0.3$ kpc for 70% of the galaxies). The dust yield per supernova is $0.1 \le y_d/M_\odot \le 3.3$, with 70% of the galaxies requiring $y_d < 0.25 M_\odot$. Three galaxies (REBELS-12, 14, 39) require $y_d > 1 M_\odot$. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [CII] nicely follow the local $L_{\rm CII}-$SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is of $0.11\, y_P^{-2}$ Gyr, where $y_P$ is the ratio of the gas-to-stellar distribution radius. For some systems a solution simultaneously matching the observed ($\beta, F_{\rm UV}, F_{158}$) values cannot be found. This occurs when the index $I_m = (F_{158}/F_{\rm UV})/(\beta-\beta_{\rm int})$, where $\beta_{\rm int}$ is the intrinsic UV slope, exceeds $I_m^*\approx 1120$ for a MW curve. For these objects we argue that the FIR and UV emitting regions are not co-spatial, questioning the use of the IRX-$\beta$ relation., Comment: 22 pages, 7 figures, accepted by MNRAS. Comments welcome

ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature $T_d$ remains mostly unconstrained due to the few available FIR continuum data at redshift $z>5$. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star formation. Using a new method based on simultaneous [CII] 158$\mu$m line and underlying dust continuum measurements, we derive $T_ d$ in the continuum and [CII] detected $z\approx 7$ galaxies in the ALMA Large Project REBELS sample. We find $39\ \mathrm{K} < T_d < 58\ \mathrm{K}$, and dust masses in the narrow range $M_d = (0.9-3.6)\times 10^7 M_{\odot}$. These results allow us to extend for the first time the reported $T_d(z)$ relation into the Epoch of Reionization. We produce a new physical model that explains the increasing $T_ d(z)$ trend with the decrease of gas depletion time, $t_{dep}=M_g/\mathrm{SFR}$, induced by the higher cosmological accretion rate at early times; this hypothesis yields $T_d \propto (1+z)^{0.4}$. The model also explains the observed $T_d$ scatter at a fixed redshift. We find that dust is warmer in obscured sources, as a larger obscuration results in more efficient dust heating. For UV-transparent (obscured) galaxies, $T_d$ only depends on the gas column density (metallicity), $T_d \propto N_H^{1/6}$ ($T_d \propto Z^{-1/6}$). REBELS galaxies are on average relatively transparent, with effective gas column densities around $N_H \simeq (0.03-1)\times 10^{21} \mathrm{cm}^{-2}$. We predict that other high-$z$ galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated $T_d \gg 60$ K, are significantly obscured, low-metallicity systems. In fact $T_d$ is higher in metal-poor systems due to their smaller dust content, which for fixed $L_{ IR}$ results in warmer temperatures., Comment: Accepted for publication in MNRAS

We analyse FIR dust continuum measurements for 14 galaxies (redshift z approximate to 7) in the ALMA Reionization Era Bright Emission Line Survey (REBELS) Large Program to derive their physical properties. Our model uses three input data, i.e. (a) the UV spectral slope, beta, (b) the observed UV continuum flux at 1500 angstrom, F-1500, (c) the observed continuum flux at approximate to 158 mu m, F-158, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have 28-90.5 per cent of their star formation obscured; the total (UV+IR) star formation rates are in the range 31.5 < SFR/(M-circle dot yr(-1)) < 129.5. The sample-averaged dust mass and temperature are (1.3 +/- 1.1) x 10(7) M-circle dot and 52 +/- 11 K, respectively. However, in some galaxies dust is particularly abundant (REBELS-14, M-d' approximate to 3.4 x 10(7) M-circle dot), or hot (REBELS-18, T-d' approximate to 67 K). The dust distribution is compact (<0.3 kpc for 70 per cent of the galaxies). The inferred dust yield per supernova is 0.1 <= y(d)/M-circle dot <= 3.3, with 70 per cent of the galaxies requiring y(d) < 0.25 M-circle dot. Three galaxies (REBELS-12, 14, 39) require y(d) > 1 M-circle dot, which is likely inconsistent with pure SN production, and might require dust growth via accretion of heavy elements from the interstellar medium. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [C II] nicely follow the local L-CII-SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is 0.11 y(p)(-2) Gyr, where y(p) is the ratio of the gas-to-stellar distribution radius. For some systems, a solution simultaneously matching the observed (beta, F-1500, F-158) values cannot be found. This occurs when the index I-m = (F-158/F-1500)/(beta - beta(int)), where beta(int) is the intrinsic UV slope, exceeds I-m* approximate to 1120 for an MW curve. F

We include a fully coupled treatment of metal and dust enrichment into the DELPHI semi-analytic model of galaxy formation to explain the dust content of 13 Lyman break galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at z similar or equal to 7. We find that the galaxy dust mass, M-d, is regulated by the combination of Type II supernova dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard time-scale tau(0) = 30 Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of similar to 0.07-0.1 per cent and a UV-to-total star formation rate relation such that log(psi(UV)) = -0.05 [log(psi)](2) + 0.86 log(psi) - 0.05 (implying that 55-80 per cent of the star formation is obscured) for REBELS galaxies with stellar mass M-* = 10(9)-10(10)M(circle dot). This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at z = 7. However, 2 out of the 13 systems show dust-to-stellar mass ratios (similar to 0.94-1.1 per cent) that are up to 18 times larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing to to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor of similar to 2. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions.

ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature T-d remains mostly unconstrained due to the few available FIR continuum data at redshift z > 5. This introduces large uncertainties in several properties of high-z galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star formation. Using a new method based on simultaneous [CII] 158-mu m line and underlying dust continuum measurements, we derive T-d in the continuum and [CII] detected z approximate to 7 galaxies in the ALMA Large Project REBELS sample. We find 39 < T-d < 58 K, and dust masses in the narrow range M-d = (0.9-3.6) x 10(7) M-circle dot. These results allow us to extend for the first time the reported T-d(z) relation into the Epoch of Reionization. We produce a new physical model that explains the increasing T-d(z) trend with the decrease of gas depletion time, t(d)(ep) = M-g/SFR, induced by the higher cosmological accretion rate at early times; this hypothesis yields T-d proportional to (1 + z)(0.4). The model also explains the observed T-d scatter at a fixed redshift. We find that dust is warmer in obscured sources, as a larger obscuration results in more efficient dust heating. For UV-transparent (obscured) galaxies, T-d only depends on the gas column density (metallicity), T-d proportional to N-H(1/6) (T-d proportional to Z(-1/6)). REBELS galaxies are on average relatively transparent, with effective gas column densities around N-H similar or equal to (0.03-1) x 10(21) cm(-2). We predict that other high-z galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated T-d >> 60 K, are significantly obscured, low-metallicity systems. In fact, T-d is higher in metal-poor systems due to their smaller dust content, which for fixed L-IR results in warmer temperatures.

We present the average [C ii] 158 mu m emission line sizes of UV-bright star-forming galaxies at z similar to 7. Our results are derived from a stacking analysis of [C ii] 158 mu m emission lines and dust continua observed by the Atacama Large Millimeter/submillimeter Array (ALMA), taking advantage of the large program Reionization Era Bright Emission Line Survey. We find that the average [C ii] emission at z similar to 7 has an effective radius r ( e ) of 2.2 +/- 0.2 kpc. It is greater than or similar to 2x larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for z less than or similar to 6 galaxies. Additionally, we compared the average [C ii] size with 4 < z < 6 galaxies observed by the ALMA Large Program to INvestigate [C ii] at Early times (ALPINE). By analyzing [C ii] sizes of 4 < z < 6 galaxies in two redshift bins, we find an average [C ii] size of r (e) = 2.2 +/- 0.2 kpc and r (e) = 2.5 +/- 0.2 kpc for z similar to 5.5 and z similar to 4.5 galaxies, respectively. These measurements show that star-forming galaxies, on average, show no evolution in the size of the [C ii] 158 mu m emitting regions at redshift between z similar to 7 and z similar to 4. This finding suggests that the star-forming galaxies could be morphologically dominated by gas over a wide redshift range.

Context. The Euclid mission is expected to discover thousands of z > 6 galaxies in three deep fields, which together will cover a ∼50 deg2 area. However, the limited number of Euclid bands (four) and the low availability of ancillary data could make the identification of z > 6 galaxies challenging. Aims. In this work we assess the degree of contamination by intermediate-redshift galaxies (z = 1–5.8) expected for z > 6 galaxies within the Euclid Deep Survey. Methods. This study is based on ∼176 000 real galaxies at z = 1–8 in a ∼0.7 deg2 area selected from the UltraVISTA ultra-deep survey and ∼96 000 mock galaxies with 25.3 ≤ H < 27.0, which altogether cover the range of magnitudes to be probed in the Euclid Deep Survey. We simulate Euclid and ancillary photometry from fiducial 28-band photometry and fit spectral energy distributions to various combinations of these simulated data. Results. We demonstrate that identifying z > 6 galaxies with Euclid data alone will be very effective, with a z > 6 recovery of 91% (88%) for bright (faint) galaxies. For the UltraVISTA-like bright sample, the percentage of z = 1–5.8 contaminants amongst apparent z > 6 galaxies as observed with Euclid alone is 18%, which is reduced to 4% (13%) by including ultra-deep Rubin (Spitzer) photometry. Conversely, for the faint mock sample, the contamination fraction with Euclid alone is considerably higher at 39%, and minimised to 7% when including ultra-deep Rubin data. For UltraVISTA-like bright galaxies, we find that Euclid (IE − YE) > 2.8 and (YE − JE) < 1.4 colour criteria can separate contaminants from true z > 6 galaxies, although these are applicable to only 54% of the co

We present the average [C ii] 158 mu m emission line sizes of UV-bright star-forming galaxies at z similar to 7. Our results are derived from a stacking analysis of [C ii] 158 mu m emission lines and dust continua observed by the Atacama Large Millimeter/submillimeter Array (ALMA), taking advantage of the large program Reionization Era Bright Emission Line Survey. We find that the average [C ii] emission at z similar to 7 has an effective radius r ( e ) of 2.2 +/- 0.2 kpc. It is greater than or similar to 2x larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for z less than or similar to 6 galaxies. Additionally, we compared the average [C ii] size with 4 < z < 6 galaxies observed by the ALMA Large Program to INvestigate [C ii] at Early times (ALPINE). By analyzing [C ii] sizes of 4 < z < 6 galaxies in two redshift bins, we find an average [C ii] size of r (e) = 2.2 +/- 0.2 kpc and r (e) = 2.5 +/- 0.2 kpc for z similar to 5.5 and z similar to 4.5 galaxies, respectively. These measurements show that star-forming galaxies, on average, show no evolution in the size of the [C ii] 158 mu m emitting regions at redshift between z similar to 7 and z similar to 4. This finding suggests that the star-forming galaxies could be morphologically dominated by gas over a wide redshift range.

The Reionization Era Bright Emission Line Survey (REBELS) is a cycle-7 ALMA Large Program (LP) that is identifying and performing a first characterization of many of the most luminous star-forming galaxies known in the z > 6.5 universe. REBELS is providing this probe by systematically scanning 40 of the brightest UV-selected galaxies identified over a 7 deg(2) area for bright [C ii](158 mu m) and [O iii](88 mu m) lines and dust-continuum emission. Selection of the 40 REBELS targets was done by combining our own and other photometric selections, each of which is subject to extensive vetting using three completely independent sets of photometry and template-fitting codes. Building on the observational strategy deployed in two pilot programs, we are increasing the number of massive interstellar medium (ISM) reservoirs known at z > 6.5 by similar to 4-5x to >30. In this manuscript, we motivate the observational strategy deployed in the REBELS program and present initial results. Based on the first-year observations, 18 highly significant >= 7 sigma [C ii](158 mu m) lines have already been discovered, the bulk of which (13/18) also show >= 3.3 sigma dust-continuum emission. These newly discovered lines more than triple the number of bright ISM-cooling lines known in the z > 6.5 universe, such that the number of ALMA-derived redshifts at z > 6.5 rival Ly alpha discoveries. An analysis of the completeness of our search results versus star formation rate (SFR) suggests an similar to 79% efficiency in scanning for [C ii](158 mu m) when the SFRUV+IR is >28 M (circle dot) yr(-1). These new LP results further demonstrate ALMA's efficiency as a "redshift machine," particularly in the Epoch of Reionization.

The neutral atomic gas content of individual galaxies at large cosmological distances has until recently been difficult to measure due to the weakness of the hyperfine H I 21 cm transition. Here we estimate the H I gas mass of a sample of main-sequence star-forming galaxies at z similar to 6.5-7.8 surveyed for [C II] 158 mu m emission as part of the Reionization Era Bright Emission Line Survey (REBELS), using a recent calibration of the [C II]-to-H I conversion factor. We find that the H I gas mass excess in galaxies increases as a function of redshift, with an average of M-HI /M-star approximate to 10, corresponding to H I gas mass fractions of f(HI) = M-HI /(M-star + M-HI) = 90%, at z approximate to 7. Based on the [C II] 158 mu m luminosity function (LF) derived from the same sample of galaxies, we further place constraints on the cosmic H I gas mass density in galaxies (rho(HI)) at this redshift, which we measure to be rho(HI) = 7.1(-3.0)(+6.4) x 10(6) M-circle dot Mpc(-3). This estimate is substantially lower by a factor of approximate to 10 than that inferred from an extrapolation of damped Ly alpha absorber (DLA) measurements and largely depends on the exact [C II] LF adopted. However, we find this decrease in rho(HI) to be consistent with recent simulations and argue that this apparent discrepancy is likely a consequence of the DLA sight lines predominantly probing the substantial fraction of H I gas in high-z galactic halos, whereas [C II] traces the H I in the ISM associated with star formation. We make predictions for this buildup of neutral gas in galaxies as a function of redshift, showing that at z greater than or similar to 5, only approximate to 10% of the cosmic H I gas content is confined in galaxies and associated with the star-forming ISM.

Star cluster formation in the early universe and their contribution to reionization remains to date largely unconstrained. Here we present JWST/NIRCam imaging of the most highly magnified galaxy known at z ~ 6, the Sunrise arc. We identify six young massive star clusters (YMCs) with measured radii spanning ~ 20 pc down to ~ 1 pc (corrected for lensing magnification), estimated stellar masses of ~ $10^{(6-7)}$ Msun, and with ages 1-30 Myr based on SED fitting to photometry measured in 8 filters extending to rest-frame 7000A. The resulting stellar mass surface densities are higher than 1000 Msun pc$^{-2}$ (up to a few $10^5$ Msun pc$^{-2}$) and their inferred dynamical ages qualify the majority of these systems as gravitationally-bound stellar clusters. The star cluster ages map the progression of star formation along the arc, with to evolved systems (>~ 10 Myr old) followed by very young clusters. The youngest stellar clusters (< 5 Myr) show evidence of prominent Hbeta + [OIII]4959,5007 emission, based on photometry, with equivalent widths larger than 1000 A rest-frame, and are hosted in a 200 pc sized star-forming complex. Such a region dominates the ionizing photon production, with a high efficiency log($\xi_{ion}$ [Hz erg$^{-1}$]) ~ 25.7. A significant fraction of the recently formed stellar mass of the galaxy (> 10-30 %) occurred in these YMCs. We speculate that such sources of ionizing radiation boost the ionizing photon production efficiency which eventually carve ionized channels that might favor the escape of Lyman continuum radiation. The survival of some of the clusters would make them the progenitors of massive and relatively metal-poor globular clusters in the local Universe., Comment: 18 pages, 8 figures, 3 tables. ApJ, Accepted

Recent ALMA large surveys unveiled the presence of significant dust continuum emission in star-forming galaxies at $z>4$. Unfortunately, such large programs -- i.e. ALPINE ($z\sim 5$) and REBELS ($z \sim 7$) -- only provide us with a single Far-Infrared (FIR) continuum data point for their individual targets. Therefore, high-$z$ galaxies FIR spectral energy densities (SEDs) remain mostly unconstrained, hinging on an assumption for their dust temperature ($T_{\rm d}$) in the SED fitting procedure. This introduces uncertainties in the inferred dust masses ($M_{\rm d }$), infrared luminosities ($L_{\rm IR}$), and obscured Star Formation Rate (SFR) fraction at $z > 4$. In this work we use a method that allows us to constrain $T_{\rm d}$ with a single band measurement by combining the $158\ \mathrm{\mu m}$ continuum information with the overlying [CII] emission line. We analyse the $21$ [CII] and FIR continuum detected $z\sim 5$ galaxies in ALPINE, finding a range of $T_{\rm d}=25-60\ \mathrm{K}$ and $M_{\rm d} = 0.6-25.1\ \times 10^{7}\ \mathrm{M_{\odot}}$. Given the measured stellar masses of ALPINE galaxies, the inferred dust yields are around $M_{\rm d}/M_{\star} = (0.2-8) \times 10^{-3}$, consistent with theoretical dust-production constraints. We find that $8$ out of $21$ ALPINE galaxies have $L_{\rm IR} \geq 10^{12}\ \mathrm{L_{\odot}}$, comparable to UltraLuminous IR Galaxies (ULIRGs). Relying on ultraviolet-to-optical SED fitting, the SFR was underestimated by up to $2$ orders of magnitude in $4$ of these $8$ ULIRGs-like galaxies. We conclude that these $4$ peculiar sources should be characterised by a two-phase interstellar medium structure with "spatially-segregated" FIR and ultraviolet emitting regions., Comment: Accepted for publication in MNRAS. 13 pages, 4 figures

We present a reduction and analysis of the \textit{James Webb Space Telescope} (JWST) SMACS~0723 field using new post-launch calibrations to conduct a search for ultra-high-redshift galaxies ($z > 9$) present within the Epoch of Reionisation. We conduct this search by modelling photometric redshifts in several ways for all sources and by applying conservative magnitude cuts ($m_{\rm F200W} < 28$) to identify strong Lyman breaks greater than 1 magnitude. We find four $z > 9$ candidate galaxies which have not previously been identified, with one object at $z = 11.5$, and another which is possibly a close pair of galaxies. We measure redshifts for candidate galaxies from other studies and find the recovery rate to be only 23 per cent, with many being assigned lower redshift, dusty solutions in our work. Most of our $z > 9$ sample show evidence for Balmer-breaks, or extreme emission lines from H$\beta$ and [OIII], demonstrating that the stellar populations could be advanced in age or very young depending on the cause of the F444W excess. We discuss the resolved structures of these early galaxies and find that the S\'{e}rsic indices reveal a mixture of light concentration levels, but that the sizes of all our systems are exceptionally small ($< 0.5$~kpc). These systems have stellar masses M$_{*} \sim 10^{9.0}$ M$_{\odot}$, with our $z \sim 11.5$ candidate a dwarf galaxy with a stellar mass M$_{*} \sim 10^{7.8}$ -- $10^{8.2}$ M$_{\odot}$. These candidate ultra high-redshift galaxies are excellent targets for future NIRSpec observations aimed to better understand their physical nature., Comment: Updated to accepted mnras manuscript, 12 pages, 8 figures, 6 tables, added section 3.3 and Table 5

The earliest JWST observations have revealed an unexpected abundance of super-early ($z>10$), massive ($M_*\approx 10^9\, M_\odot$) galaxies at the bright-end ($M_{\rm UV}\approx -21$) of the ultraviolet luminosity function (UV LF). We present a minimal physical model that explains the observed galaxy abundance at $z=10-14$. The model primarily combines (a) the halo mass function, with (b) an obscured star formation fraction prescription that is consistent with findings of the ALMA REBELS dusty galaxy survey. It has been successfully tested on well-known UV LFs up to $z=7$. The weak evolution from $z=7$ to $z\approx 14$ of the LF bright-end arises from a conspiracy between a decreasing dust attenuation, making galaxies brighter, that almost exactly compensates for the increasing shortage of their host halos. The model also predicts that galaxies at $z > 11$ should contain negligible amounts of dust. We speculate that dust could have been efficiently ejected during the very first phases of galaxy build-up., Comment: 6 pages, 4 figures. Analysis of early release JWST data. Comments welcome

The neutral atomic gas content of individual galaxies at large cosmological distances has until recently been difficult to measure due to the weakness of the hyperfine HI 21-cm transition. Here we estimate the HI gas mass of a sample of main-sequence star-forming galaxies at $z\sim 6.5 - 7.8$ surveyed for [CII]$-158\mu$m emission as part of the Reionization Era Bright Emission Line Survey (REBELS), using a recent calibration of the [CII]-to-HI conversion factor. We find that the HI gas mass excess in galaxies increases as a function of redshift, with an average of $M_{\rm HI} / M_\star \approx 10$, corresponding to HI gas mass fractions of $f_{\rm HI} = M_{\rm HI} / (M_\star + M_{\rm HI}) = 90\%$, at $z\approx 7$. Based on the [CII]-$158\mu$m luminosity function (LF) derived from the same sample of galaxies, we further place constraints on the cosmic HI gas mass density in galaxies ($\rho_{\rm HI}$) at this redshift, which we measure to be $\rho_{\rm HI} = 7.1^{+6.4}_{-3.0} \times 10^{6}\,M_{\odot}\,{\rm Mpc^{-3}}$. This estimate is substantially lower by a factor of $\approx 10$ than that inferred from an extrapolation of damped Lyman-$\alpha$ absorber (DLA) measurements, and largely depend on the exact [CII] LF adopted. However, we find this decrease in $\rho_{\rm HI}$ to be consistent with recent simulations and argue that this apparent discrepancy is likely a consequence of the DLA sightlines predominantly probing the substantial fraction of HI gas in high-$z$ galactic halos, whereas [CII] traces the HI in the ISM associated with star formation. We make predictions for this build-up of neutral gas in galaxies as a function of redshift, showing that at $z\gtrsim 5$ only $\approx 10\%$ of the cosmic HI gas content is confined in galaxies and associated with the star-forming ISM., Comment: The accepted version of the manuscript in press for ApJL

We present the average [CII] $158\,\rm{\mu m}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{\mu m}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective radius $r_e$ of $2.2\pm0.2\,\rm{kpc}$. It is $\gtrsim2\times$ larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for $z\lesssim6$ galaxies. Additionally, we compared the average [CII] size with $4

(Abridged) The Euclid mission is expected to discover thousands of z>6 galaxies in three Deep Fields, which together will cover a ~40 deg2 area. However, the limited number of Euclid bands and availability of ancillary data could make the identification of z>6 galaxies challenging. In this work, we assess the degree of contamination by intermediate-redshift galaxies (z=1-5.8) expected for z>6 galaxies within the Euclid Deep Survey. This study is based on ~176,000 real galaxies at z=1-8 in a ~0.7 deg2 area selected from the UltraVISTA ultra-deep survey, and ~96,000 mock galaxies with 25.3$\leq$H<27.0, which altogether cover the range of magnitudes to be probed in the Euclid Deep Survey. We simulate Euclid and ancillary photometry from the fiducial, 28-band photometry, and fit spectral energy distributions (SEDs) to various combinations of these simulated data. Our study demonstrates that identifying z>6 with Euclid data alone will be very effective, with a z>6 recovery of 91(88)% for bright (faint) galaxies. For the UltraVISTA-like bright sample, the percentage of z=1-5.8 contaminants amongst apparent z>6 galaxies as observed with Euclid alone is 18%, which is reduced to 4(13)% by including ultra-deep Rubin (Spitzer) photometry. Conversely, for the faint mock sample, the contamination fraction with Euclid alone is considerably higher at 39%, and minimized to 7% when including ultra-deep Rubin data. For UltraVISTA-like bright galaxies, we find that Euclid (I-Y)>2.8 and (Y-J)<1.4 colour criteria can separate contaminants from true z>6 galaxies, although these are applicable to only 54% of the contaminants, as many have unconstrained (I-Y) colours. In the most optimistic scenario, these cuts reduce the contamination fraction to 1% whilst preserving 81% of the fiducial z>6 sample. For the faint mock sample, colour cuts are infeasible., Comment: 27 pages, 14 figures, 9 tables; version accepted for publication by A&A

We investigate the ionizing properties of the pair of bright Ly$\alpha$ emitting galaxies BDF521 and BDF2195 at z=7.012 in order to constrain their contribution to the formation of the BDF "reionized bubble". We obtain constraints on UV emission lines (CIV$\lambda 1548$ doublet, HeII$\lambda 1640$, OIII]$\lambda 1660$ doublet, and CIII]$\lambda 1909$ doublet) from deep VLT-XSHOOTER observations and compare them to those available for other high-redshift objects, and to models with mixed stellar and AGN emission. We use this spectroscopic information together with the photometry available in the field to constrain the physical properties of the two objects using the spectro-photometric fitting code BEAGLE. We do not detect any significant emission at the expected position of UV lines, with 3$\sigma$ upper limits of EW$\lesssim$2-7 AA rest-frame. We find that the two objects have lower CIII] emission than expected on the basis of the correlation between the Ly$\alpha$ and CIII] EWs. The EW limits on CIV and HeII emission exclude pure AGN templates at $\sim2-3\sigma$ significance, and only models with a $\lesssim$40% AGN contribution are compatible with the observations. The two objects are found to be relatively young ($\sim$20-30 Myrs) and metal-poor ($\lesssim 0.3 Z_{\odot}$) with stellar masses of a few $10^9M_{\odot}$. Their production rate of hydrogen ionizing photons per intrinsic UV luminosity is log($\xi_{ion}^*$/Hz erg$^{-1}$)=25.02-25.26, consistent with values typically found in high-redshift galaxies, but more than twice lower than values measured in $z>$7 galaxies with strong CIII] and/or optical line emission ($\simeq$25.6-25.7). The two BDF emitters have no evidence of higher than average ionizing capabilities and are not capable of reionizing their surroundings by their own means under realistic assumptions on the escape fraction of ionizing photons. (abridged), Comment: 9 pages, 4 figures, accepted for publication in A&A

We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman Break Galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at $z\simeq 7$. We find that the galaxy dust mass, $M_d$, is regulated by the combination of SNII dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard timescale $\tau_0= 30$ Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of $\sim 0.07-0.1\%$ and a UV-to-total star formation rate relation such that $log (\psi_{\rm UV}) = -0.05 ~[log (\psi)]^{2} + 0.86 ~log(\psi) -0.05$ (implying that 55-80\% of the star formation is obscured) for REBELS galaxies with stellar mass $M_* = 10^{9-10} M_\odot$. This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at $z=7$. However, 2 out of the 13 systems show dust-to-stellar mass ratios ($\sim 0.94-1.1\%$) that are up to $18\times$ larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing $\tau_0$ to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor $ \sim 2$. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions., Comment: Accepted to MNRAS

ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature $T_d$ remains mostly unconstrained due to the few available FIR continuum data at redshift $z>5$. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star formation. Using a new method based on simultaneous [CII] 158$\mu$m line and underlying dust continuum measurements, we derive $T_ d$ in the continuum and [CII] detected $z\approx 7$ galaxies in the ALMA Large Project REBELS sample. We find $39\ \mathrm{K} < T_d < 58\ \mathrm{K}$, and dust masses in the narrow range $M_d = (0.9-3.6)\times 10^7 M_{\odot}$. These results allow us to extend for the first time the reported $T_d(z)$ relation into the Epoch of Reionization. We produce a new physical model that explains the increasing $T_ d(z)$ trend with the decrease of gas depletion time, $t_{dep}=M_g/\mathrm{SFR}$, induced by the higher cosmological accretion rate at early times; this hypothesis yields $T_d \propto (1+z)^{0.4}$. The model also explains the observed $T_d$ scatter at a fixed redshift. We find that dust is warmer in obscured sources, as a larger obscuration results in more efficient dust heating. For UV-transparent (obscured) galaxies, $T_d$ only depends on the gas column density (metallicity), $T_d \propto N_H^{1/6}$ ($T_d \propto Z^{-1/6}$). REBELS galaxies are on average relatively transparent, with effective gas column densities around $N_H \simeq (0.03-1)\times 10^{21} \mathrm{cm}^{-2}$. We predict that other high-$z$ galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated $T_d \gg 60$ K, are significantly obscured, low-metallicity systems. In fact $T_d$ is higher in metal-poor systems due to their smaller dust content, which for fixed $L_{ IR}$ results in warmer temperatures., Comment: Accepted for publication in MNRAS

We analyse FIR dust continuum measurements for 14 galaxies ($z\approx 7$) in the ALMA REBELS LP to derive their physical properties. Our model uses three input data: (a) the UV spectral slope, $\beta$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158\mu$m, $F_{158}$, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have (28-90.5)% of their star formation obscured; the total (UV+IR) star formation rates are in the range $31.5 < {\rm SFR}/ (M_\odot {\rm yr}^{-1}) < 129.5$. The sample-averaged dust mass and temperature are $(1.3\pm 1.1)\times 10^7 M_\odot$ and $52 \pm 11$ K, respectively. In some galaxies dust is abundant (REBELS-14, $M'_d \approx 3.4 \times 10^7 M_\odot$), or hot (REBELS-18, $T'_d \approx 67$ K). The dust distribution is compact ($<0.3$ kpc for 70% of the galaxies). The dust yield per supernova is $0.1 \le y_d/M_\odot \le 3.3$, with 70% of the galaxies requiring $y_d < 0.25 M_\odot$. Three galaxies (REBELS-12, 14, 39) require $y_d > 1 M_\odot$. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [CII] nicely follow the local $L_{\rm CII}-$SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is of $0.11\, y_P^{-2}$ Gyr, where $y_P$ is the ratio of the gas-to-stellar distribution radius. For some systems a solution simultaneously matching the observed ($\beta, F_{\rm UV}, F_{158}$) values cannot be found. This occurs when the index $I_m = (F_{158}/F_{\rm UV})/(\beta-\beta_{\rm int})$, where $\beta_{\rm int}$ is the intrinsic UV slope, exceeds $I_m^*\approx 1120$ for a MW curve. For these objects we argue that the FIR and UV emitting regions are not co-spatial, questioning the use of the IRX-$\beta$ relation., Comment: 22 pages, 7 figures, accepted by MNRAS. Comments welcome

The Reionization Era Bright Emission Line Survey (REBELS) is a cycle-7 ALMA Large Program (LP) that is identifying and performing a first characterization of many of the most luminous star-forming galaxies known in the z > 6.5 universe. REBELS is providing this probe by systematically scanning 40 of the brightest UV-selected galaxies identified over a 7 deg(2) area for bright [C ii](158 mu m) and [O iii](88 mu m) lines and dust-continuum emission. Selection of the 40 REBELS targets was done by combining our own and other photometric selections, each of which is subject to extensive vetting using three completely independent sets of photometry and template-fitting codes. Building on the observational strategy deployed in two pilot programs, we are increasing the number of massive interstellar medium (ISM) reservoirs known at z > 6.5 by similar to 4-5x to >30. In this manuscript, we motivate the observational strategy deployed in the REBELS program and present initial results. Based on the first-year observations, 18 highly significant >= 7 sigma [C ii](158 mu m) lines have already been discovered, the bulk of which (13/18) also show >= 3.3 sigma dust-continuum emission. These newly discovered lines more than triple the number of bright ISM-cooling lines known in the z > 6.5 universe, such that the number of ALMA-derived redshifts at z > 6.5 rival Ly alpha discoveries. An analysis of the completeness of our search results versus star formation rate (SFR) suggests an similar to 79% efficiency in scanning for [C ii](158 mu m) when the SFRUV+IR is >28 M (circle dot) yr(-1). These new LP results further demonstrate ALMA's efficiency as a "redshift machine," particularly in the Epoch of Reionization.

We analyse FIR dust continuum measurements for 14 galaxies (redshift z approximate to 7) in the ALMA Reionization Era Bright Emission Line Survey (REBELS) Large Program to derive their physical properties. Our model uses three input data, i.e. (a) the UV spectral slope, beta, (b) the observed UV continuum flux at 1500 angstrom, F-1500, (c) the observed continuum flux at approximate to 158 mu m, F-158, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have 28-90.5 per cent of their star formation obscured; the total (UV+IR) star formation rates are in the range 31.5 < SFR/(M-circle dot yr(-1)) < 129.5. The sample-averaged dust mass and temperature are (1.3 +/- 1.1) x 10(7) M-circle dot and 52 +/- 11 K, respectively. However, in some galaxies dust is particularly abundant (REBELS-14, M-d' approximate to 3.4 x 10(7) M-circle dot), or hot (REBELS-18, T-d' approximate to 67 K). The dust distribution is compact ( 1 M-circle dot, which is likely inconsistent with pure SN production, and might require dust growth via accretion of heavy elements from the interstellar medium. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [C II] nicely follow the local L-CII-SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is 0.11 y(p)(-2) Gyr, where y(p) is the ratio of the gas-to-stellar distribution radius. For some systems, a solution simultaneously matching the observed (beta, F-1500, F-158) values cannot be found. This occurs when the index I-m = (F-158/F-1500)/(beta - beta(int)), where beta(int) is the intrinsic UV slope, exceeds I-m* approximate to 1120 for an MW curve. For these objects, we argue that the FIR and UV emitting regions are not co-spatial, questioning the use of the IRX-beta relation.

We include a fully coupled treatment of metal and dust enrichment into the DELPHI semi-analytic model of galaxy formation to explain the dust content of 13 Lyman break galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at z similar or equal to 7. We find that the galaxy dust mass, M-d, is regulated by the combination of Type II supernova dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard time-scale tau(0) = 30 Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of similar to 0.07-0.1 per cent and a UV-to-total star formation rate relation such that log(psi(UV)) = -0.05 [log(psi)](2) + 0.86 log(psi) - 0.05 (implying that 55-80 per cent of the star formation is obscured) for REBELS galaxies with stellar mass M-* = 10(9)-10(10)M(circle dot). This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at z = 7. However, 2 out of the 13 systems show dust-to-stellar mass ratios (similar to 0.94-1.1 per cent) that are up to 18 times larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing to to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor of similar to 2. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions.

We analyse FIR dust continuum measurements for 14 galaxies (redshift z approximate to 7) in the ALMA Reionization Era Bright Emission Line Survey (REBELS) Large Program to derive their physical properties. Our model uses three input data, i.e. (a) the UV spectral slope, beta, (b) the observed UV continuum flux at 1500 angstrom, F-1500, (c) the observed continuum flux at approximate to 158 mu m, F-158, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have 28-90.5 per cent of their star formation obscured; the total (UV+IR) star formation rates are in the range 31.5 < SFR/(M-circle dot yr(-1)) < 129.5. The sample-averaged dust mass and temperature are (1.3 +/- 1.1) x 10(7) M-circle dot and 52 +/- 11 K, respectively. However, in some galaxies dust is particularly abundant (REBELS-14, M-d' approximate to 3.4 x 10(7) M-circle dot), or hot (REBELS-18, T-d' approximate to 67 K). The dust distribution is compact ( 1 M-circle dot, which is likely inconsistent with pure SN production, and might require dust growth via accretion of heavy elements from the interstellar medium. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [C II] nicely follow the local L-CII-SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is 0.11 y(p)(-2) Gyr, where y(p) is the ratio of the gas-to-stellar distribution radius. For some systems, a solution simultaneously matching the observed (beta, F-1500, F-158) values cannot be found. This occurs when the index I-m = (F-158/F-1500)/(beta - beta(int)), where beta(int) is the intrinsic UV slope, exceeds I-m* approximate to 1120 for an MW curve. For these objects, we argue that the FIR and UV emitting regions are not co-spatial, questioning the use of the IRX-beta relation.

We include a fully coupled treatment of metal and dust enrichment into the DELPHI semi-analytic model of galaxy formation to explain the dust content of 13 Lyman break galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at z similar or equal to 7. We find that the galaxy dust mass, M-d, is regulated by the combination of Type II supernova dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard time-scale tau(0) = 30 Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of similar to 0.07-0.1 per cent and a UV-to-total star formation rate relation such that log(psi(UV)) = -0.05 [log(psi)](2) + 0.86 log(psi) - 0.05 (implying that 55-80 per cent of the star formation is obscured) for REBELS galaxies with stellar mass M-* = 10(9)-10(10)M(circle dot). This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at z = 7. However, 2 out of the 13 systems show dust-to-stellar mass ratios (similar to 0.94-1.1 per cent) that are up to 18 times larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing to to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor of similar to 2. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions.

Over the past decades, rest-frame ultraviolet (UV) observations have provided large samples of UV luminous galaxies at redshift (z) greater than 6, during the so-called epoch of reionization. While a few of these UV identified galaxies revealed significant dust reservoirs, very heavily dust-obscured sources at these early times have remained elusive. They are limited to a rare population of extreme starburst galaxies, and companions of rare quasars. These studies conclude that the contribution of dust-obscured galaxies to the cosmic star formation rate density at $z>6$ is sub-dominant. Recent ALMA and Spitzer observations have identified a more abundant, less extreme population of obscured galaxies at $z=3-6$. However, this population has not been confirmed in the reionization epoch so far. Here, we report the discovery of two dust-obscured star forming galaxies at $z=6.6813\pm0.0005$ and $z=7.3521\pm0.0005$. These objects are not detected in existing rest-frame UV data, and were only discovered through their far-infrared [CII] lines and dust continuum emission as companions to typical UV-luminous galaxies at the same redshift. The two galaxies exhibit lower infrared luminosities and star-formation rates than extreme starbursts, in line with typical star-forming galaxies at $z\sim7$. This population of heavily dust-obscured galaxies appears to contribute 10-25 per cent to the $z>6$ cosmic star formation rate density., Comment: 16 pages, 7 figures, 4 tables, published in Nature

The locations of Ly-$\alpha$ emitting galaxies (LAEs) at the end of the Epoch of Reionisation (EoR) are expected to correlate with regions of ionised hydrogen, traced by the redshifted 21~cm hyperfine line. Mapping the neutral hydrogen around regions with detected and localised LAEs offers an avenue to constrain the brightness temperature of the Universe within the EoR by providing an expectation for the spatial distribution of the gas, thereby providing prior information unavailable to power spectrum measurements. We use a test set of 12 hours of observations from the Murchison Widefield Array (MWA) in extended array configuration, to constrain the neutral hydrogen signature of 58 LAEs, detected with the Subaru Hypersuprime Cam in the \textit{Silverrush} survey, centred on $z$=6.58. We assume that detectable emitters reside in the centre of ionised HII bubbles during the end of reionization, and predict the redshifted neutral hydrogen signal corresponding to the remaining neutral regions using a set of different ionised bubble radii. A prewhitening matched filter detector is introduced to assess detectability. We demonstrate the ability to detect, or place limits upon, the amplitude of brightness temperature fluctuations, and the characteristic HII bubble size. With our limited data, we constrain the brightness temperature of neutral hydrogen to $\Delta{\rm T}_B<$30 mK ($<$200 mK) at 95% (99%) confidence for lognormally-distributed bubbles of radii, $R_B =$ 15$\pm$2$h^{-1}$cMpc., Comment: 9 pages, 6 figures, accepted for publication in MNRAS