Search

The Atacama Large Millimetre/submillimetre Array (ALMA) receivers and technical papers are cited fewer than once in every six publications. This citation shortage is impeding the development of future (sub)millimetre instruments. In an effort to facilitate the correct citations of ALMA receivers and technical papers, this memo provides a comprehensive list of papers for the scientific community. This list was produced in discussion with the scientific and instrumentalist community, based on a June 2024 survey at the European Southern Observatory workshop on the ALMA Wideband Sensitivity Upgrade, as well as with the ALMA technical staff. The authors now encourage the community to enhance their already-excellent ALMA science with the appropriate references to ensure future (sub)millimetre instrumentation can keep addressing the key questions about our Universe., Comment: 8 pages; 1 figure; 3 tables; to be updated as more papers become available. Suggestions are warmly welcomed. ALMA memos solely contains the opinion of the authors in discussion with the astronomer and instrument builder community, and does not reflect the official policy of the ALMA telescope

We report the highest-redshift detection of [O I] 63 $\mu$m from a luminous quasar, J2054-0005, at $z=6.04$ based on the Atacama Large Millimeter/sub-millimeter Array Band 9 observations. The [O I] 63 $\mu$m line luminosity is $(4.5\pm1.5) \times 10^{9}~L_{\rm \odot}$, corresponding to the [O I] 63 $\mu$m-to-far-infrared luminosity ratio of $\approx 6.7\times10^{-4}$, which is consistent with the value obtained in the local universe. Remarkably, [O I] 63 $\mu$m is as bright as [C II] 158 $\mu$m, resulting in the [O I]-to-[C II] line luminosity ratio of $1.3\pm0.5$. Based on a careful comparison of the luminosity ratios of [O I] 63 $\mu$m, [C II] 158 $\mu$m, and dust continuum emission to models of photo-dissociation regions, we find that J2054-0005 has a gas density log($n_{\rm H}$/cm$^{-3}$)$=3.7\pm0.3$ and an incident far-ultraviolet radiation field of log($G/G_{\rm 0}$)$= 3.0\pm0.1$, showing that [O I] 63 $\mu$m serves as an important coolant of the dense and warm gas in J2054-0005. A close examination of the [O I] and [C II] line profiles suggests that the [O I] line may be partially self-absorbed, however deeper observations are needed to verify this conclusion. Regardless, the gas density and incident radiation field are in a broad agreement with the values obtained in nearby star-forming galaxies and objects with [O I] 63 $\mu$m observations at $z=1-3$ with the Herschel Space Observatory. These results demonstrate the power of ALMA high-frequency observations targeting [O I] 63 $\mu$m to examine the properties of photo-dissociation regions in high-redshift galaxies., Comment: Submitted to PASJ. 5 figures and 3 tables. Comments are very welcome

We present new multi-frequency ALMA continuum observations of the massive [$\log_{10}(M_\star/M_\odot) = 10.3_{-0.2}^{+0.1}$], UV-luminous [$M_\mathrm{UV} = -21.7 \pm 0.2$] $z=7.31$ galaxy REBELS-25 in Bands 3, 4, 5, and 9. Combining the new observations with previously-taken data in Bands 6 and 8, we cover the dust continuum emission of the galaxy in six distinct bands -- spanning rest-frame $50-350\,\mu$m -- enabling simultaneous constraints on its dust mass ($M_\mathrm{dust}$), temperature ($T_\mathrm{dust}$) and emissivity index ($\beta_\mathrm{IR}$) via modified blackbody fitting. Given a fiducial model of optically thin emission, we infer a cold dust temperature of $T_\mathrm{dust} = 32_{-6}^{+9}\,$K and a high dust mass of $\log_{10}(M_\mathrm{dust}/M_\odot) = 8.2_{-0.4}^{+0.6}$, and moderately optically thick dust does not significantly alter these estimates. If we assume dust production is solely through supernovae (SNe), the inferred dust yield would be high, $y = 0.7_{-0.4}^{+2.3}\,M_\odot$ per SN. Consequently, we argue grain growth in the interstellar medium of REBELS-25 also contributes to its dust build-up. This is supported by the steep dust emissivity index $\beta_\mathrm{IR} = 2.5 \pm 0.4$ we measure for REBELS-25, as well as by its high stellar mass, dense interstellar medium, and metal-rich nature. Our results suggest that constraining the dust emissivity indices of high-redshift galaxies is important not only to mitigate systematic uncertainties in their dust masses and obscured star formation rates, but also to assess if dust properties evolve across cosmic time. We present an efficient observing setup to do so with ALMA, combining observations of the peak and Rayleigh-Jeans tail of the dust emission., Comment: 14 pages main text, 6 figures + appendices; re-submitted to MNRAS after a positive referee report

We introduce COBRA (Cosmology with Optimally factorized Bases of Radial Approximants), a novel framework for rapid computation of large-scale structure observables. COBRA separates scale dependence from cosmological parameters in the linear matter power spectrum while also minimising the number of necessary basis terms $N_b$, thus enabling direct and efficient computation of derived and nonlinear observables. Moreover, the dependence on cosmological parameters is efficiently approximated using radial basis function interpolation. We apply our framework to decompose the linear matter power spectrum in the standard $\Lambda$CDM scenario, as well as by adding curvature, dynamical dark energy and massive neutrinos, covering all redshifts relevant for Stage IV surveys. With only a dozen basis terms $N_b$, COBRA reproduces exact Boltzmann solver calculations to $\sim 0.1\%$ precision, which improves further to $0.02\%$ in the pure $\Lambda$CDM scenario. Using our decomposition, we recast the one-loop redshift space galaxy power spectrum in a separable minimal-basis form, enabling $\sim 4000$ model evaluations per second at $0.02\%$ precision on a single thread. This constitutes a considerable improvement over previously existing methods (e.g., FFTLog) opening a window for efficient computations of higher loop and higher order correlators involving multiple powers of the linear matter power spectra. The resulting factorisation can also be utilised in clustering, weak lensing and CMB analyses. Our implementation will be made public upon publication., Comment: 7+4 pages, to be submitted to PRL

Submillimeter and millimeter wavelengths provide a unique view of the Universe, from the gas and dust that fills and surrounds galaxies to the chromosphere of our own Sun. Current single-dish facilities have presented a tantalising view of the brightest (sub-)mm sources, and interferometers have provided the exquisite resolution necessary to analyse the details in small fields, but there are still many open questions that cannot be answered with current facilities. In this report we summarise the science that is guiding the design of the Atacama Large Aperture Submillimeter Telescope (AtLAST). We demonstrate how tranformational advances in topics including star formation in high redshift galaxies, the diffuse circumgalactic medium, Galactic ecology, cometary compositions and solar flares motivate the need for a 50m, single-dish telescope with a 1-2 degree field of view and a new generation of highly multiplexed continuum and spectral cameras. AtLAST will have the resolution to drastically lower the confusion limit compared to current single-dish facilities, whilst also being able to rapidly map large areas of the sky and detect extended, diffuse structures. Its high sensitivity and large field of view will open up the field of submillimeter transient science by increasing the probability of serendipitous detections. Finally, the science cases listed here motivate the need for a highly flexible operations model capable of short observations of individual targets, large surveys, monitoring programmes, target of opportunity observations and coordinated observations with other observatories. AtLAST aims to be a sustainable, upgradeable, multipurpose facility that will deliver orders of magnitude increases in sensitivity and mapping speeds over current and planned submillimeter observatories., Comment: 47 pages, 12 figures. For further details on AtLAST see https://atlast.uio.no

We present Atacama Large Millimetre/submillimetre Array (ALMA) observations of the [O$_{\rm III}$] 88 $\mu$m emission of a sample of thirteen galaxies at $z$ = 6 to 7.6 selected as [C$_{\rm II}$]-emitting companion sources of quasars. To disentangle the origins of the luminous Oxygen line in the $z$ > 6 Universe, we looked at emission-line galaxies that are selected through an excellent star-formation tracer [C$_{\rm II}$] with star-formation rates between 9 and 162 M$_{\odot}$/yr. Direct observations reveal [O$_{\rm III}$] emission in just a single galaxy (L$_{\rm [O_{\rm III}]}$/L$_{\rm [C_{\rm II}]}$ = 2.3), and a stacked image shows no [O$_{\rm III}$] detection, providing deep upper limits on the L$_{\rm [O_{\rm III}]}$/L$_{\rm [C_{\rm II}]}$ ratios in the $z > 6$ Universe (L$_{\rm [O_{\rm III}]}$/L$_{\rm [C_{\rm II}]}$ < 1.2 at 3${\sigma}$). While the fidelity of this sample is high, no obvious optical/near-infrared counterpart is seen in the JWST imaging available for four galaxies. Additionally accounting for low-redshift CO emitters, line stacking shows that our sample-wide result remains robust: The enhanced L$_{\rm [O_{\rm III}]}$/L$_{\rm [C_{\rm II}]}$ reported in the first billion years of the Universe is likely due to the selection towards bright, blue Lyman-break galaxies with high surface star-formation rates or young stellar populations. The deep upper limit on the rest-frame 90 $\mu$m continuum emission (< 141 $\mu$Jy at 3${\sigma}$), implies a low average dust temperature (T$_{\rm dust}$ < 30K) and high dust mass (M$_{\rm dust}$ ~ 10$^8$ M$_{\odot}$). As more normal galaxies are explored in the early Universe, synergy between JWST and ALMA is fundamental to further investigate the ISM properties of the a broad range of samples of high-$z$ galaxies., Comment: 20 pages; 13 figures; accepted for publication in MNRAS

Unveiling the emergence and prevalence of massive/bright galaxies during the epoch of reionization and beyond, within the first 600 million years of the Universe, stands as a pivotal pursuit in astronomy. Remarkable progress has been made by JWST in identifying an immense population of bright galaxies, which hints at exceptionally efficient galaxy assembly processes. However, the underlying physical mechanisms propelling their rapid growth remain unclear. With this in mind, millimeter and submillimeter-wave spectroscopic observations of redshifted far-infrared spectral lines, particularly the [O III] 88 micron and [C II] 158 micron lines, offers a crucial pathway to address this fundamental query. To this end, we develop a dual-polarization sideband-separating superconductor-insulator-superconductor (SIS) mixer receiver, FINER, for the Large Millimeter Telescope (LMT) situated in Mexico. Harnessing advancements from ALMA's wideband sensitivity upgrade (WSU) technology, FINER covers radio frequencies spanning 120-360 GHz, delivering an instantaneous intermediate frequency (IF) of 3-21 GHz per sideband per polarization, which is followed by a set of 10.24 GHz-wide digital spectrometers. At 40% of ALMA's light-collecting area, the LMT's similar atmospheric transmittance and FINER's 5 times wider bandwidth compared to ALMA culminate in an unparalleled spectral scanning capability in the northern hemisphere, paving the way for finer spectral-resolution detection of distant galaxies., Comment: 12 pages, 8 figures, and 3 tables. Proceedings paper presented in SPIE Astronomical Telescope and Instrumentation 2024

We report on deep SCUBA-2 observations at 850$\mu$m and NOEMA spectroscopic measurements at 2 mm of the environment surrounding the luminous, massive ($M_{*} \approx 2 \times 10^{11}$ M$_{\odot}$) Herschel-selected source HerBS-70. This source was revealed by previous NOEMA observations to be a binary system of dusty star-forming galaxies at $z= 2.3$, with the East component (HerBS-70E) hosting an Active Galactic Nucleus (AGN). The SCUBA-2 observations detected, in addition to the binary system, twenty-one sources at $> 3.5 \sigma$ over an area of $\sim 25$ square comoving Mpc with a sensitivity of $\sigma_{850} = 0.75$ mJy. The surface density of continuum sources around HerBS-70 is three times higher than for field galaxies. The NOEMA spectroscopic measurements confirm the protocluster membership of three of the nine brightest sources through their CO(4 - 3) line emission, yielding a volume density 36 times higher than for field galaxies. All five confirmed sub-mm galaxies in the HerBS-70 system have relatively short gas depletion times ($80 - 500$ Myr), indicating the onset of quenching for this protocluster core due to the depletion of gas. The dark matter halo mass of the HerBS-70 system is estimated around $5 \times{} 10^{13}$ M$_{\odot}$, with a projected current-day mass of $10^{15}$ M$_{\odot}$, similar to the local Virgo and Coma clusters. These observations support the claim that DSFGs, in particular the ones with observed multiplicity, can trace cosmic overdensities., Comment: 19 pages, 13 figures, accepted for publication in MNRAS

We present JWST NIRCam imaging of B14-65666 ("Big Three Dragons"), a bright Lyman-break galaxy system ($M_\text{UV}=-22.5$ mag) at $z=7.15$. The high angular resolution of NIRCam reveals the complex morphology of two galaxy components: galaxy E has a compact core (E-core), surrounded by diffuse, extended, rest-frame optical emission, which is likely to be tidal tails; and galaxy W has a clumpy and elongated morphology with a blue UV slope ($\beta_\text{UV}=-2.2\pm0.1$). The flux excess, F356W$-$F444W, peaks at the E-core ($1.05^{+0.08}_{-0.09}$ mag), tracing the presence of strong [OIII] 4960,5008 \r{A} emission. ALMA archival data show that the bluer galaxy W is brighter in dust continua than the redder galaxy E, while the tails are bright in [OIII] 88 $\mathrm{\mu m}$. The UV/optical and sub-mm SED fitting confirms that B14-65666 is a major merger in a starburst phase as derived from the stellar mass ratio (3:1 to 2:1) and the star-formation rate, $\simeq1$ dex higher than the star-formation main sequence at the same redshift. The galaxy E is a dusty ($A_\text{V}=1.2\pm0.1$ mag) starburst with a possible high dust temperature ($\ge63$-$68$ K). The galaxy W would have a low dust temperature ($\le27$-$33$ K) or patchy stellar-and-dust geometry, as suggested from the infrared excess (IRX) and $\beta_\text{UV}$ diagram. The high optical-to-FIR [OIII] line ratio of the E-core shows its lower gas-phase metallicity ($\simeq0.2$ Z$_{\odot}$) than the galaxy W. These results agree with a scenario where major mergers disturb morphology and induce nuclear dusty starbursts triggered by less-enriched inflows. B14-65666 shows a picture of complex stellar buildup processes during major mergers in the epoch of reionization., Comment: 18 pages, 6 figures, 4 tables. Submitted to ApJ

GHZ2/GLASS-z12 has been recently observed by JWST with both NIRSpec and MIRI spectrographs, making it the most distant galaxy ($z_{spec}=12.34$) with complete spectroscopic coverage from rest-frame UV to optical. It is identified as a strong CIV$_{1549}$ emitter with many detected emission lines (NIV], HeII, OIII], NIII], CIII], [OII], [NeIII], [OIII], and H$\alpha$), including a remarkable OIII$_{1333}$ Bowen fluorescence line. We analyze in this paper the joint NIRSpec+MIRI spectral data set. Combining six optical diagnostics (R2, R3, R23, O32, Ne3O2, and Ne3O2Hd), we find extreme ionization conditions, with O32 $=1.39 \pm 0.19$ and Ne3O2 $=0.37 \pm 0.18$ in stark excess compared to typical values in the ISM at lower redshifts. These line properties are compatible either with an AGN or with a compact, dense star-forming environment ($\Sigma_{\rm SFR}$ $\sim 10^2$-$10^3$ Msun/yr/kpc$^2$), with a high ionization parameter ($\log_{10}$(U) $=-1.75 \pm 0.16$), a high ionizing photon production efficiency $\log(\xi_{\rm ion}) = 25.7_{-0.1}^{+0.2}$, and a low, although not pristine, metal content between $5\%$ and $11\%$ Z$_\odot$ (confirmed by the T$_e$ method), indicating a rapid metal enrichment in the last few Myrs. These properties also suggest that a substantial amount of ionizing photons ($\sim 10\%$) are leaking outside. The general lessons learned from GHZ2 are the following: (i) the UV to optical combined nebular indicators are broadly in agreement with UV-only or optical-only indicators. (ii) UV+optical diagnostics fail to discriminate between an AGN and star-formation in a low metallicity, high density, and extreme ionization environment. (iii) comparing the nebular line ratios with local analogs may be approaching its limits at $z \gtrsim 10$, as this approach is potentially challenged by the unique conditions of star formation experienced by galaxies at these extreme redshifts., Comment: 14 pages, 3 figures, and 1 table. Resubmitted to ApJ. A section (3.4) has been added to the previous version

The James Webb Space Telescope (JWST) has discovered a surprising population of bright galaxies in the very early universe (< 500 Myrs after the Big Bang) that is hard to explain with conventional galaxy formation models and whose physical properties remain to be fully understood. Insight into the internal physics of galaxies is captured best via observations of excited-state atomic transitions of ionized gas, but beyond z~7-9, the brightest spectral signatures are redshifted into the mid-infrared regime, where observations are increasingly more difficult. Here, we present the first detection of a hydrogen recombination line (H{\alpha}) and doubly-ionized oxygen ([OIII]4959,5007{\AA}) at z>10 using the JWST Mid-Infrared Instrument, MIRI. These detections place the bright galaxy GHZ2/GLASS-z12 at z=12.33+/-0.02, making it the most distant astronomical object with direct spectroscopic detection of these lines and the brightest confirmed object at this epoch. These observations provide key insights into the conditions of this primeval galaxy, which shows hard ionizing conditions rarely seen in the local Universe and likely driven by compact, young (<30 Myr) star formation. Its oxygen-to-hydrogen abundance is close to a tenth of the solar value, indicating a rapid metal enrichment during the earliest phases of galaxy formation. This study confirms the unique conditions of the brightest and most distant galaxies recently discovered by JWST and the huge potential of mid-IR observations to characterize these systems, opening a range of new possibilities in the study of the very early Universe., Comment: 14 pages, 6 figures. Submitted for publication. Comments are most welcome

We spectroscopically confirm the $M_{\rm UV} = -20.5$ mag galaxy GHZ2/GLASS-z12 to be at redshift $z=12.34$. The source was selected via NIRCam photometry in GLASS-JWST ERS data, providing the first evidence of a surprising abundance of bright galaxies at $z \gtrsim 10$. The NIRSpec PRISM spectrum shows detections of N IV, C IV, He II, O III, C III, O II, and Ne III lines, and the first detection at high-redshift of the O III Bowen fluorescence line at 3133 {\AA} rest-frame. The prominent C IV line with rest-frame equivalent width (EW) $\approx 46$ {\AA} puts GHZ2 in the category of extreme C IV emitters. GHZ2 displays UV lines with EWs that are only found in active galactic nuclei (AGNs) or composite objects at low/intermediate redshifts. The UV line-intensity ratios are compatible both with AGNs and star formation in a low-metallicity environment, with the low limit on the [Ne IV]/[N IV] ratio favoring a stellar origin of the ionizing photons. We discuss a possible scenario in which the high ionizing output is due to low metallicity stars forming in a dense environment. We estimate a metallicity $\lesssim 0.1 Z/{\rm Z}_{\odot}$, a high ionization parameter logU $> -2$, a N/O abundance 4--5 times the solar value, and a subsolar C/O ratio similar to the recently discovered class of nitrogen-enhanced objects. Considering its abundance patterns and the high stellar mass density ($10^4$~M$_{\odot}$~pc$^{-2}$), GHZ2 is an ideal formation site for the progenitors of today's globular clusters. The remarkable brightness of GHZ2 makes it a ``Rosetta stone'' for understanding the physics of galaxy formation within just 360 Myr after the Big Bang., Comment: Accepted for publication in The Astrophysical Journal. 20 pages, 7 figures, 2 tables

During the most active period of star formation in galaxies, which occurs in the redshift range 1

We have carried out HST snapshot observations at 1.1 $\mu$m of 281 candidate strongly lensed galaxies identified in the wide-area extragalactic surveys conducted with the Herschel space observatory. Our candidates comprise systems with flux densities at $500\,\mu$m$ S_{500}\geq 80$ mJy. We model and subtract the surface brightness distribution for 130 systems, where we identify a candidate for the foreground lens candidate. After combining visual inspection, archival high-resolution observations, and lens subtraction, we divide the systems into different classes according to their lensing likelihood. We confirm 65 systems to be lensed. Of these, 30 are new discoveries. We successfully perform lens modelling and source reconstruction on 23 systems, where the foreground lenses are isolated galaxies and the background sources are detected in the HST images. All the systems are successfully modelled as a singular isothermal ellipsoid. The Einstein radii of the lenses and the magnifications of the background sources are consistent with previous studies. However, the background source circularised radii (between 0.34 kpc and 1.30 kpc) are $\sim$3 times smaller than the ones measured in the sub-mm/mm for a similarly selected and partially overlapping sample. We compare our lenses with those in the SLACS survey, confirming that our lens-independent selection is more effective at picking up fainter and diffuse galaxies and group lenses. This sample represents the first step towards characterising the near-IR properties and stellar masses of the gravitationally lensed dusty star-forming galaxies., Comment: 57 pages, 18 figures, 11 tables. Accepted for publication in MNRAS

Intrinsic alignment (IA) of source galaxies represents an important contaminant for upcoming cosmic shear surveys. In particular, it is expected on general grounds that IA contains a B-mode while the weak lensing signal does not. Thus, a detection of B-modes offers the possibility to study directly the IA signal of the sources. Galaxy clusters exhibit strong IA and are therefore a natural candidate to look for a B-mode signal. We forecast the signal-to-noise ratio (SNR) for B-modes from IA of galaxy clusters in the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). We use a perturbative model for the IA multipoles based on the Effective Field Theory of Intrinsic Alignments (EFT of IA), which has recently been validated against N-body simulations. We forecast SNR $\approx 12$ and find that this detectability is not significantly impacted by different analysis choices. Lastly, we also apply our forecast to clusters in the redMaPPer SDSS and DESY1 samples. We find SNR $\approx 5$ and SNR $\approx 3$, respectively, suggesting a detection is within reach, provided accurate redshift information is available., Comment: 10 pages, 5 figures, accepted for publication in the OJA

The James Webb Space Telescope (JWST) has revealed extremely distant galaxies at unprecedentedly early cosmic epochs from its deep imaging using the technique of photometric redshift estimation, with its subsequent spectroscopy confirming their redshifts unambiguously, demonstrating the ability of JWST to probe the earliest galaxies, one of its major scientific goals. However, as larger samples continue to be followed up spectroscopically, it has become apparent that nearly all photometric redshifts at these epochs are biased high with confidence >>99%, for as yet unclear reasons. Here we show that this is the same statistical effect that was predicted in different contexts by Sir Arthur Eddington in 1913, in that there exist more lower redshift galaxies to be scattered upwards than the reverse. The bias depends on the shape of the intrinsic redshift distribution, but as an approximate heuristic, all ultra-high photometric redshift estimates must be corrected downwards by up to one standard deviation., Comment: Nature Astronomy, published 03 October 2023. Authors' revised version. 4 pages, one figure. Published version at https://www.nature.com/articles/s41550-023-02093-8

The mean free path of ionizing photons, $\lambda_{\rm mfp}$, is a critical parameter for modeling the intergalactic medium (IGM) both during and after reionization. We present direct measurements of $\lambda_{\rm mfp}$ from QSO spectra over the redshift range $5

The z-GAL survey observed 137 bright Herschel-selected targets with the IRAM NOrthern Extended Millimeter Array, with the aim to measure their redshift and study their properties. Several of them have been resolved into multiple sources. Consequently, robust spectroscopic redshifts have been measured for 165 individual galaxies in the range 0.8

(Abridged) We present the dust properties of 125 bright Herschel galaxies selected from the z-GAL survey. The large instantaneous bandwidth of NOEMA provides an exquisite sampling of the underlying dust continuum emission at 2 and 3 mm in the observed frame, with flux densities in at least four side bands for each source. Together with the available Herschel 250, 350, and 500 micron and SCUBA-2 850 micron flux densities, the spectral energy distribution of each source can be analyzed from the far-infrared to the millimeter, with a fine sampling of the Rayleigh-Jeans tail. This wealth of data provides a solid basis to derive robust dust properties, in particular the dust emissivity index, beta, and the dust temperature, T(dust). In order to demonstrate our ability to constrain the dust properties, we used a flux-generated mock catalog and analyzed the results under the assumption of an optically thin and optically thick modified black body emission. For the z-GAL sources, we report a range of dust emissivities with beta ~ 1.5 - 3 estimated up to high precision with relative uncertainties that vary in the range 7% - 15%, and an average of 2.2 +/- 0.3. We find dust temperatures varying from 20 to 50 K with an average of T(dust) ~ 30 K for the optically thin case and ~38 K in the optically thick case. For all the sources, we estimate the dust masses and apparent infrared luminosities (based on the optically thin approach). An inverse correlation is found between T(dust) and beta, which is similar to what is seen in the local Universe. Finally, we report an increasing trend in the dust temperature as a function of redshift at a rate of 6.5 +/- 0.5 K/z for this 500 micron-selected sample. Based on this study, future prospects are outlined to further explore the evolution of dust temperature across cosmic time., Comment: Accepted for publication on A&A; 32 pages; 26 figures

(Abridged) Using the IRAM NOEMA interferometer, we measures the redshifts of 126 bright galaxies detected in the Herschel H-ATLAS, HeLMS, and HerS surveys. We report reliable spectroscopic redshifts for a total of 124 of the Herschel-selected galaxies. The redshifts are estimated from scans of the 3 and 2-mm bands (and, in one case, the 1-mm band) and are based on the detection of at least two emission lines. Together with the Pilot Programme (Neri et al. 2020), including spectroscopic redshifts of 11 sources, our survey has derived precise redshifts for 135 bright Herschel-selected galaxies, making it the largest sample of high-z galaxies with robust redshifts to date. Most emission lines detected are from 12CO (mainly from J=2-1 to 5-4), with some sources seen in [CI] and H2O emission lines. The spectroscopic redshifts are in the range 0.8

The intrinsic alignment of galaxies is an important ingredient for modelling weak-lensing measurements, and a potentially valuable cosmological and astrophysical signal. In this paper, we present HYMALAIA: a new model to predict the intrinsic alignments of biased tracers. HYMALAIA is based on a perturbative expansion of the statistics of the Lagrangian shapes of objects, which is then advected to Eulerian space using the fully non-linear displacement field obtained from $N$-body simulations. We demonstrate that HYMALAIA is capable of consistently describing monopole and quadrupole of halo shape-shape and matter-shape correlators, and that, without increasing the number of free parameters, it does so more accurately than other perturbatively inspired models such as the non-linear alignment (NLA) model and the tidal-alignment-tidal-torquing (TATT) model., Comment: 17 pages, 9 figures. Accepted version published in MNRAS

Molecular outflows are expected to play a key role in galaxy evolution at high redshift. To study the impact of outflows on star formation at the epoch of reionization, we performed sensitive ALMA observations of OH 119 $\mu$m toward J2054-0005, a luminous quasar at $z=6.04$. The OH line is detected and exhibits a P-Cygni profile that can be fitted with a broad blue-shifted absorption component, providing unambiguous evidence of an outflow, and an emission component at near-systemic velocity. The mean and terminal outflow velocities are estimated to be $v_\mathrm{out}\approx670~\mathrm{km~s}^{-1}$ and $1500~\mathrm{km~s}^{-1}$, respectively, making the molecular outflow in this quasar one of the fastest at the epoch of reionization. The OH line is marginally spatially resolved for the first time in a quasar at $z>6$, revealing that the outflow extends over the central 2 kpc region. The mass outflow rate is comparable to the star formation rate ($\dot{M}_\mathrm{out}/\mathrm{SFR}\sim2$), indicating rapid ($\sim10^7~\mathrm{yr}$) quenching of star formation. The mass outflow rate in a sample star-forming galaxies and quasars at $4

The protoclusters in the epoch of reionization, traced by galaxies overdensity regions, are ideal laboratories for studying the process of stellar assembly and cosmic reionization. We present the spectroscopic confirmation of the core of the most distant protocluster at $z = 7.88$, A2744-z7p9OD, with the James Webb Space Telescope NIRSpec integral field unit spectroscopy. The core region includes as many as 4 galaxies detected in [OIII] 4960 \AA\ and 5008 \AA\ in a small area of $\sim 3\arcsec \times 3\arcsec$, corresponding to $\sim$ 11 kpc $\times$ 11 kpc, after the lensing magnification correction. Three member galaxies are also tentatively detected in dust continuum in Atacama Large Millimeter/submillimeter Array Band 6, which is consistent with their red ultraviolet continuum slopes, $\beta \sim -1.3$. The member galaxies have stellar masses in the range of log($M_{*}/M_{\rm \odot}$) $\sim 7.6-9.2$ and star formation rates of $\sim 3-50$ $M_{\rm \odot}$ yr$^{-1}$, showing a diversity in their properties. FirstLight cosmological simulations reproduce the physical properties of the member galaxies including the stellar mass, [OIII] luminosity, and dust-to-stellar mass ratio, and predict that the member galaxies are on the verge of merging in a few to several tens Myr to become a large galaxy with $M_{\rm *}\sim 6\times10^{9} M_{\rm \odot}$. The presence of a multiple merger and evolved galaxies in the core region of A2744-z7p9OD indicates that environmental effects are already at work 650 Myr after the Big Bang., Comment: 12 pages, 4 figures, 1 table, published in The Astrophysical Journal Letters

We test the regime of validity of the effective field theory (EFT) of intrinsic alignments (IA) at the one-loop level by comparing with 3D halo shape statistics in N-body simulations. This model is based on the effective field theory of large-scale structure (EFT of LSS) and thus a theoretically well-motivated extension of the familiar non-linear alignment (NLA) model and the tidal-alignment-tidal-torquing (TATT) model. It contains a total of $8$ free bias parameters. Specifically, we measure the dark matter halo shape-shape multipoles $P_{EE}^{(0)}(k), P_{EE}^{(2)}(k), P_{BB}^{(0)}(k), P_{BB}^{(2)}(k)$ as well as the matter-shape multipoles $P_{\delta E}^{(0)}(k), P_{\delta E}^{(2)}(k)$ from the simulations and perform a joint fit to determine the largest wavenumber $k_{\text{max}}$ up to which the theory predictions from the EFT of IA are consistent with the measurements. We find that the EFT of IA is able to describe intrinsic alignments of dark matter halos up to $k_\text{max}=0.30\,h/\text{Mpc}$ at $z=0$. This demonstrates a clear improvement over other existing alignment models like NLA and TATT, which are only accurate up to $k_\text{max}=0.05\,h/\text{Mpc}$ . We examine the posterior distributions of the higher-order bias parameters, and show that their inclusion is necessary to describe intrinsic alignments in the quasi-linear regime. Further, the EFT of IA is able to accurately describe the auto-spectrum of intrinsic alignment B-modes, in contrast to the other alignment models considered., Comment: 35 pages, 2 tables, 6 figures

We present the results of 300 pc resolution ALMA imaging of the [OIII] 88 $\mu$m line and dust continuum emission from a $z = 8.312$ Lyman break galaxy MACS0416_Y1. The velocity-integrated [OIII] emission has three peaks which are likely associated with three young stellar clumps of MACS0416_Y1, while the channel map shows a complicated velocity structure with little indication of a global velocity gradient unlike what was found in [CII] 158 $\mu$m at a larger scale, suggesting random bulk motion of ionized gas clouds inside the galaxy. In contrast, dust emission appears as two individual clumps apparently separating or bridging the [OIII]/stellar clumps. The cross correlation coefficient between dust and ultraviolet-related emission (i.e., [OIII] and ultraviolet continuum) is unity on a galactic scale, while it drops at < 1 kpc, suggesting well mixed geometry of multi-phase interstellar media on sub-kpc scales. If the cutoff scale characterizes different stages of star formation, the cutoff scale can be explained by gravitational instability of turbulent gas. We also report on a kpc-scale off-center cavity embedded in the dust continuum image. This could be a superbubble producing galactic-scale outflows, since the energy injection from the 4 Myr starburst suggested by a spectral energy distribution analysis is large enough to push the surrounding media creating a kpc-scale cavity., Comment: 18 pages, 6 figures, 1 table, published in ApJ

We present integral field, far-infrared (FIR) spectroscopy of Mrk 54, a local Lyman Continuum Emitter (LCE), obtained with FIFI-LS on the Stratospheric Observatory for Infrared Astronomy. This is only the second time, after Haro 11, that [C II] 158 $\mu$m and [O III] 88 $\mu$m spectroscopy of the known LCEs have been obtained. We find that Mrk 54 has a strong [C II] emission that accounts for $\sim1$% of the total FIR luminosity, whereas it has only moderate [O III] emission, resulting in the low [O III]/[C II] luminosity ratio of $0.22\pm0.06$. In order to investigate whether [O III]/[C II] is a useful tracer of $f_{\rm esc}$ (LyC escape fraction), we examine the correlations of [O III]/[C II] and (i) the optical line ratio of $\rm O_{32} \equiv$ [O III] 5007 \AA/[O II] 3727 \AA, (ii) specific star formation rate, (iii) [O III] 88 $\mu$m/[O I] 63 $\mu$m ratio, (iv) gas phase metallicity, and (v) dust temperature based on a combined sample of Mrk 54 and the literature data from the Herschel Dwarf Galaxy Survey and the LITTLE THINGS Survey. We find that galaxies with high [O III]/[C II] luminosity ratios could be the result of high ionization (traced by $\rm O_{32}$), bursty star formation, high ionized-to-neutral gas volume filling factors (traced by [O III] 88 $\mu$m/[O I] 63 $\mu$m), and low gas-phase metallicities, which is in agreement with theoretical predictions. We present an empirical relation between the [O III]/[C II] ratio and $f_{\rm esc}$ based on the combination of the [O III]/[C II] and $\rm O_{32}$ correlation, and the known relation between $\rm O_{32}$ and $f_{\rm esc}$. The relation implies that high-redshift galaxies with high [O III]/[C II] ratios revealed by ALMA may have $f_{\rm esc}\gtrsim0.1$, significantly contributing to the cosmic reionization., Comment: 14 pages, 5 figures, Accepted for publication in ApJ

We analyse the molecular and atomic emission lines of 71 bright Herschel-selected galaxies between redshifts 1.4 to 4.6 detected by the Atacama Large Millimetre/submillimetre Array. These lines include a total of 156 CO, [C I], and H2O emission lines. For 46 galaxies, we detect two transitions of CO lines, and for these galaxies we find gas properties similar to those of other dusty star-forming galaxy (DSFG) samples. A comparison to photo-dissociation models suggests that most of Herschel-selected galaxies have similar interstellar medium conditions as local infrared-luminous galaxies and high-redshift DSFGs, although with denser gas and more intense far-ultraviolet radiation fields than normal star-forming galaxies. The line luminosities agree with the luminosity scaling relations across five orders of magnitude, although the star-formation and gas surface density distributions (i.e., Schmidt-Kennicutt relation) suggest a different star-formation phase in our galaxies (and other DSFGs) compared to local and low-redshift gas-rich, normal star-forming systems. The gas-to-dust ratios of these galaxies are similar to Milky Way values, with no apparent redshift evolution. Four of 46 sources appear to have CO line ratios in excess of the expected maximum (thermalized) profile, suggesting a rare phase in the evolution of DSFGs. Finally, we create a deep stacked spectrum over a wide rest-frame frequency (220-890 GHz) that reveals faint transitions from HCN and CH, in line with previous stacking experiments., Comment: 30 pages, 17 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal. Comments are warmly welcomed

We present the results of a survey of CO(1-0) emission in 14 infrared luminous dusty star forming galaxies (DSFGs) at 2 < z < 4 with the NSF's Karl G. Jansky Very Large Array. All sources are detected in CO(1-0), with an ~1arcsec angular resolution. Seven sources show extended and complex structure. We measure CO luminosities of $({\mu})L'_{CO(1-0)}=0.4-2.9x10^{11}$ K km s$^{-1}$ pc$^2$, and molecular gas masses of (${\mu}$)M$_{H2}$ = 1.3 - 8.6 x 10$^{11}$ Mo, where ({\mu}) is the magnification factor. The derived molecular gas depletion times of t$_{\rm dep}$ = 40 - 460 Myr, cover the expected range of both normal star forming galaxies and starbursts. Comparing to the higher-J CO transitions previously observed for the same sources, we find CO temperature brightness ratios of r$_{32/10}$ = 0.4 - 1.4, r$_{43/10}$ = 0.4 - 1.7, and r$_{54/10}$ = 0.3 - 1.3. We find a wide range of CO spectral line energy distributions (SLEDs), in agreement with other high-z DSFGs, with the exception of three sources that are most comparable to the Cloverleaf and APM08279+5255. Based on radiative transfer modelling of the CO SLEDs we determine densities of n$_{H2}$ = 0.3 - 8.5 x 10$^3$ cm$^{-3}$ and temperatures of T$_K$ = 100 - 200 K. Lastly, four sources are detected in the continuum, three have radio emission consistent with their infrared derived star formation rates, while HerBS-70E requires an additional synchrotron radiation component from an active galactic nucleus. Overall, we find that even though the sample is similarly luminous in the infrared, by tracing the CO(1-0) emission a diversity of galaxy and excitation properties are revealed, demonstrating the importance of CO(1-0) observations in combination to higher-J transitions., Comment: 25pages , 13 figures, Accepted for publication at ApJ

We present 101 and 151 GHz ALMA continuum images for 85 fields selected from Herschel observations that have 500 micron flux densities >80 mJy and 250-500 micron colours consistent with z > 2, most of which are expected to be gravitationally lensed or hyperluminous infrared galaxies. Approximately half of the Herschel 500 micron sources were resolved into multiple ALMA sources, but 11 of the 15 brightest 500 micron Herschel sources correspond to individual ALMA sources. For the 37 fields containing either a single source with a spectroscopic redshift or two sources with the same spectroscopic redshift, we examined the colour temperatures and dust emissivity indices. The colour temperatures only vary weakly with redshift and are statistically consistent with no redshift-dependent temperature variations, which generally corresponds to results from other samples selected in far-infrared, submillimetre, or millimetre bands but not to results from samples selected in optical or near-infrared bands. The dust emissivity indices, with very few exceptions, are largely consistent with a value of 2. We also compared spectroscopic redshifts to photometric redshifts based on spectral energy distribution templates designed for infrared-bright high-redshift galaxies. While the templates systematically underestimate the redshifts by ~15%, the inclusion of ALMA data decreases the scatter in the predicted redshifts by a factor of ~2, illustrating the potential usefulness of these millimetre data for estimating photometric redshifts., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Infrared (IR), sub-millimetre (sub-mm) and millimetre (mm) databases contain a huge quantity of high quality data. However, a large part of these data are photometric, and are thought not to be useful to derive a quantitative information on the nebular emission of galaxies. The aim of this project is first to identify galaxies at z > 4-6, and in the epoch of reionization from their sub-mm colours. We also aim at showing that the colours can be used to try and derive physical constraints from photometric bands, when accounting for the contribution from the IR fine structure lines to these photometric bands. We model the flux of IR fine structure lines with CLOUDY, and add them to the dust continuum emission with CIGALE. Including or not emission lines in the simulated spectral energy distribution (SED) modifies the broad band emission and colours. The introduction of the lines allows to identify strong star forming galaxies at z > 4 - 6 from the log10 (PSW_250um/PMW_350um) versus log10 (LABOCA_870um/PLW_500um) colour-colour diagramme. By comparing the relevant models to each observed galaxy colour, we are able to roughly estimate the fluxes of the lines, and the associated nebular parameters. This method allows to identify a double sequence in a plot built from the ionization parameter and the gas metallicity. The HII and photodissociation region (PDR) fine structure lines are an essential part of the SEDs. It is important to add them when modelling the spectra, especially at z > 4 - 6 where their equivalent widths can be large. Conversely, we show that we can extract some information on strong IR fine structure lines and on the physical parameters related to the nebular emission from IR colour-colour diagrams., Comment: Paper accepted in Astronomy and Astrophysics on 10 November 2022

The James Webb Space Telescope (JWST) has discovered a surprising abundance of bright galaxy candidates in the very early Universe ($< 500$ Myrs after the Big Bang), calling into question current galaxy formation models. Spectroscopy is needed to confirm the primeval nature of these candidates, as well as to understand how the first galaxies form stars and grow. Here we present deep spectroscopic and continuum ALMA observations towards GHZ2/GLASS-z12, one of the brightest and most robust candidates at $z > 10$ identified in the GLASS-JWST Early Release Science Program. We detect a $5.8 \sigma$ line, offset 0.5" from the JWST position of GHZ2/GLASS-z12 that, associating it with the [OIII] 88 micron transition, implies a spectroscopic redshift of $z = 12.117 \pm 0.001$. We verify the detection using extensive statistical tests. The oxygen line luminosity places GHZ2/GLASS-z12 above the [OIII]-SFR relation for metal-poor galaxies, implying an enhancement of [OIII] emission in this system while the JWST-observed emission is likely a lower-metallicity region. The lack of dust emission seen by these observations is consistent with the blue UV slope observed by JWST, which suggest little dust attenuation in galaxies at this early epoch. Further observations will unambiguously confirm the redshift and shed light on the origins of the wide and offset line and physical properties of this early galaxy. This work illustrates the synergy between JWST and ALMA and paves the way for future spectroscopic surveys of $z > 10$ galaxy candidates., Comment: Accepted to MNRAS. Comments welcome

We present spatially-resolved morphological properties of [CII] 158 $\mu$m, [OIII] 88 $\mu$m, dust, and rest-frame ultraviolet (UV) continuum emission for A1689-zD1, a strongly lensed, sub-L* galaxy at $z=7.13$, by utilizing deep Atacama Large Millimeter/submillimeter Array (ALMA) and Hubble Space Telescope (HST) observations. While the [OIII] line and UV continuum are compact, the [CII] line is extended up to a radius of $r \sim 12$ kpc. Using multi-band rest-frame far-infrared (FIR) continuum data ranging from 52-400 $\mu$m, we find an average dust temperature and emissivity index of $T_{\rm dust} = 41^{+17}_{-14}$ K and $\beta = 1.7^{+1.1}_{-0.7}$, respectively, across the galaxy. We find slight differences in the dust continuum profiles at different wavelengths, which may indicate that the dust temperature decreases with distance. We map the star-formation rate (SFR) via IR and UV luminosities and determine a total SFR of $37\pm 1~M_\odot~{\rm yr}^{-1}$ with an obscured fraction of $87\%$. While the [OIII] line is a good tracer of the SFR, the [CII] line shows deviation from the local $L_{\rm [CII]}$-SFR relations in the outskirts of the galaxy. Finally, we observe a clear difference in the line profile between [CII] and [OIII], with significant residuals ($\sim 5\sigma$) in the [OIII] line spectrum after subtracting a single Gaussian model. This suggests a possible origin of the extended [CII] structure from the cooling of hot ionized outflows. The extended [CII] and high-velocity [OIII] emission may both contribute in part to the high $L_{\rm [OIII]}$/$L_{\rm [CII]}$ ratios recently reported in $z>6$ galaxies., Comment: 22 pages, 14 figures, accepted to ApJ, presented at AAS240

We present methods to (i) graphically identify robust redshifts using emission lines in the (sub)mm regime, (ii) evaluate the capabilities of different (sub)mm practices for measuring spectroscopic redshifts, and (iii) optimise future (sub)mm observations towards increasing the fraction of robust redshifts. Using this publicly-available code (https://github.com/tjlcbakx/redshift-search-graphs), we discuss scenarios where robust redshifts can be identified using both single- and multiple-line detections, as well as scenarios where the redshift remains ambiguous, even after the detection of multiple lines. Using the redshift distribution of (sub)mm samples, we quantify the efficiencies of various practices for measuring spectroscopic redshifts, including interferometers, as well as existing and future instruments specifically designed for redshift searches. Finally, we provide a method to optimise the observation strategy for future (sub)mm spectroscopic redshift searches with the Atacama Large Millimetre/submillimetre Array, where 2 mm proves indispensable for robust redshifts in the z = 2 - 4 region., Comment: 9 pages, 8 figures, resubmitted and implemented the comments of the referee

We present spectroscopic measurements for 71 galaxies associated with 62 of the brightest high-redshift submillimeter sources from the Southern fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), while targeting 85 sources which resolved into 142. We have obtained robust redshift measurements for all sources using the 12-m Array and an efficient tuning of ALMA to optimise its use as a redshift hunter, with 73 per cent of the sources having a robust redshift identification. Nine of these redshift identifications also rely on observations from the Atacama Compact Array. The spectroscopic redshifts span a range $1.41

We have carried out the first spatially-resolved investigation of the multi-phase interstellar medium (ISM) at high redshift, using the z=4.24 strongly-lensed sub-millimetre galaxy H-ATLASJ142413.9+022303 (ID141). We present high-resolution (down to ~350 pc) ALMA observations in dust continuum emission and in the CO(7-6), H_2O (2_{1,1} - 2_{0,2}), CI(1-0) and CI(2-1) lines, the latter two allowing us to spatially resolve the cool phase of the ISM for the first time. Our modelling of the kinematics reveals that the system appears to be dominated by a rotationally-supported gas disk with evidence of a nearby perturber. We find that the CI(1-0) line has a very different distribution to the other lines, showing the existence of a reservoir of cool gas that might have been missed in studies of other galaxies. We have estimated the mass of the ISM using four different tracers, always obtaining an estimate in the range (3.2-3.8) x 10^{11} M_sol, significantly higher than our dynamical mass estimate of (0.8-1.3) x 10^{11} M_sol. We suggest that this conflict and other similar conflicts reported in the literature is because the gas-to-tracer ratios are ~4 times lower than the Galactic values used to calibrate the ISM in high-redshift galaxies. We demonstrate that this could result from a top-heavy initial mass function and strong chemical evolution. Using a variety of quantitative indicators, we show that, extreme though it is at z=4.24, ID141 will likely join the population of quiescent galaxies that appears in the Universe at z~3., Comment: This article has been accepted for publication in MNRAS, published by Oxford University Press on behalf of the Royal Astronomical Society. 24 pages

We present a feasibility study for the high-redshift galaxy part of the Science Verification Campaign with the 220-440 GHz DESHIMA 2.0 integrated superconducting spectrometer on the ASTE telescope. The first version of the DESHIMA 2.0 chip has been recently manufactured and tested in the lab. Based on these realistic performance measurements, we evaluate potential target samples and prospects for detecting the [CII] and CO emission lines. The planned observations comprise two distinct, but complementary objectives: (1) acquiring spectroscopic redshifts for dusty galaxies selected in far-infrared/mm-wave continuum surveys; (2) multi-line observations to infer physical conditions in dusty galaxies., Comment: Expanded version of the article for the conference proceedings for the 9th International Workshop on Low Temperature Detectors (submitted to JLTP, 1st November 2021). Added minor corrections and an extra panel for Fig. 2 (2021/11/2). Updated Figure 3 (2022/1/13)

Integrated superconducting spectrometer (ISS) technology will enable ultra-wideband, integral-field spectroscopy for (sub)millimeter-wave astronomy, in particular, for uncovering the dust-obscured cosmic star formation and galaxy evolution over cosmic time. Here we present the development of DESHIMA 2.0, an ISS for ultra-wideband spectroscopy toward high-redshift galaxies. DESHIMA 2.0 is designed to observe the 220-440 GHz band in a single shot, corresponding to a redshift range of $z$=3.3-7.6 for the ionized carbon emission ([C II] 158 $\mu$m). The first-light experiment of DESHIMA 1.0, using the 332-377 GHz band, has shown an excellent agreement among the on-sky measurements, the lab measurements, and the design. As a successor to DESHIMA 1.0, we plan the commissioning and the scientific observation campaign of DESHIMA 2.0 on the ASTE 10-m telescope in 2023. Ongoing upgrades for the full octave-bandwidth system include the wideband 347-channel chip design and the wideband quasi-optical system. For efficient measurements, we also develop the observation strategy using the mechanical fast sky-position chopper and the sky-noise removal technique based on a novel data-scientific approach. In the paper, we show the recent status of the upgrades and the plans for the scientific observation campaign., Comment: 9 pages, 4 figures, 1 table, accepted for publication in the Journal of Low Temperature Physics for the LTD19 special issue

We report ALMA Band 9 continuum observations of the normal, dusty star-forming galaxy A1689-zD1 at $z = 7.13$, resulting in a $\sim$4.6$\sigma$ detection at $702$ GHz. For the first time these observations probe the far infrared (FIR) spectrum shortward of the emission peak of a galaxy in the Epoch of Reionization (EoR). Together with ancillary data from earlier works, we derive the dust temperature, $T_{\rm d}$, and mass, $M_{\rm d}$, of A1689-zD1 using both traditional modified blackbody spectral energy density fitting, and a new method that relies only on the [CII] $158\ \mathrm{\mu m}$ line and underlying continuum data. The two methods give $T_{\rm d} = (42^{+13}_{-7}, 40^{+13}_{-7}$) K, and $M_{\rm d} = (1.7^{+1.3}_{-0.7}, 2.0^{+1.8}_{-1.0})\,\times{}\,10^{7} \,M_{\odot}$. Band 9 observations improve the accuracy of the dust temperature (mass) estimate by $\sim 50$% (6 times). The derived temperatures confirm the reported increasing $T_{\rm d}$-redshift trend between $z=0$ and $8$; the dust mass is consistent with a supernova origin. Although A1689-zD1 is a normal UV-selected galaxy, our results, implying that $\sim$85% of its star formation rate is obscured, underline the non-negligible effects of dust in EoR galaxies., Comment: Published in MNRAS