We present deep near-infrared $K_\textrm{s}$-band imaging for 35 of the 53 sources from the high-redshift ($z \gt 2$) radio galaxy candidate sample defined in Broderick et al. (2022, PASA, 39, e061). These images were obtained using the High-Acuity Widefield K-band Imager (HAWK-I) on the Very Large Telescope. Host galaxies are detected for 27 of the sources, with $K_\textrm{s} \approx 21.6$–23.0 mag (2$^{\prime\prime}$ diameter apertures; AB). The remaining eight targets are not detected to a median $3\unicode{x03C3}$ depth of $K_\textrm{s} \approx 23.3$ mag (2$^{\prime\prime}$ diameter apertures). We examine the radio and near-infrared flux densities of the 35 sources, comparing them to the known $z \gt 3$ powerful radio galaxies with 500-MHz radio luminosities $L_{500\,\textrm{MHz}} \gt 10^{27}$ W Hz$^{-1}$. By plotting 150-MHz flux density versus $K_\textrm{s}$-band flux density, we find that, similar to the sources from the literature, these new targets have large radio to near-infrared flux density ratios, but extending the distribution to fainter flux densities. Five of the eight HAWK-I deep non-detections have a median $3\unicode{x03C3}$ lower limit of $K_\textrm{s} \gtrsim 23.8$ mag (1$.\!^{\prime\prime}$5 diameter apertures); these five targets, along with a further source from Broderick et al. (2022, PASA, 39, e061) with a deep non-detection ($K_\textrm{s} \gtrsim 23.7$ mag; $3\unicode{x03C3}$; 2$^{\prime\prime}$ diameter aperture) in the Southern H-ATLAS Regions $K_\textrm{s}$-band Survey, are considered candidates to be ultra-high-redshift ($z \gt 5$) radio galaxies. The extreme radio to near-infrared flux density ratios ($\gt 10^5$) for these six sources are comparable to TN J0924$-$2201, GLEAM J0856$+$0223 and TGSS J1530$+$1049, the three known powerful radio galaxies at $z \gt 5$. For a selection of galaxy templates with different stellar masses, we show that $z \gtrsim 4.2$ is a plausible scenario for our ultra-high-redshift candidates if the stellar mass $M_\textrm{*} \gtrsim 10^{10.5}$ M$_\odot$. In general, the 35 targets studied have properties consistent with the previously known class of infrared-faint radio sources. We also discuss the prospects for finding more UHzRG candidates from wide and deep near-infrared surveys.

The baryon mass content (i.e. stellar and gas mass) of dark matter halos in the early Universe depends on both global factors – for example, ionising ultraviolet (UV) radiation background – and local factors – for example, star formation efficiency and assembly history. We use a lightweight semi-analytical model to investigate how both local and global factors impact the halo baryon mass content at redshifts of $z\geq 5$. Our model incorporates a time delay between when stars form and when they produce feedback of $0\leq t^d/\mathrm{Myr} \leq 30$, which can drive bursts of star formation, and a mass and redshift-dependent UV background, which captures the influence of cosmological reionisation on gas accretion onto halos. We use statistically representative halo assembly histories and assume that the cosmological gas accretion rate is proportional to the halo mass accretion rate. Delayed ($t^d$>0) feedback leads to oscillations in gas mass with cosmic time, behaviour that cannot be captured with instantaneous feedback ($t^d$=0). Highly efficient star formation drives stronger oscillations, while strong feedback impacts when oscillations occur; in contrast, inefficient star formation and weak feedback produce similar long-term behaviour to that observed in instantaneous feedback models. If the delayed feedback timescale is too long, a halo retains its gas reservoir but the feedback suppresses star formation. Our model predicts that lower mass systems (halo masses $m_\mathrm{h} \leq 10^7 \mathrm{M}_\odot$) at $z \leq 10$ should be strongly gas deficient ($m_\mathrm{g}\rightarrow 0$), whereas higher mass systems retain their gas reservoirs because they are sufficiently massive to continue accreting gas through cosmological reionisation. Interestingly, in higher mass halos, the median $m_\star/(m_\star+m_\mathrm{g}) \simeq 0.01-0.05$, but is a factor of 3–5 smaller when feedback is delayed. Our model does not include seed supermassive black hole feedback, which is necessary to explain massive quenched galaxies in the early Universe.

Feedback and outflows associated with a quasar phase are expected to be critical in quenching the most massive galaxies at high-z. Observations targeting the cool molecular and atomic phases, which dominate the mass and momentum budget of massive galaxy outflows and remove the direct fuel for star formation are, however, severely limited in high-z QSO hosts. We discuss two recent ALMA programs: one targeting molecular outflows in 3 z ∼ 6 QSO hosts using the OH 119 μm absorption line and another targeting the diffuse, predominantly atomic gas in the halos surrounding 5 QSO host between z ∼ 2 – 4 using the OH+(11 – 10) absorption line. Outflows are successfully detected in both samples and compared with outflows driven by high-z star-forming galaxies observed in the same lines. Both studies indicate that observing QSOs during the blow-out phase is crucial for studying the impact of the active nucleus on the ejection of gas from the host galaxy.

We report for the first time a relationship between galaxy kinematics and net Lyman-$\alpha$ equivalent width (net Ly$\alpha$ EW) in star-forming galaxies during the epoch of peak cosmic star formation. Building on the previously reported broadband imaging segregation of Ly$\alpha$-emitting and Ly$\alpha$-absorbing Lyman break galaxies (LBGs) at $z\sim2$ (Paper I in this series) and previously at $z\sim3$, we use the Ly$\alpha$ spectral type classification method to study the relationship between net Ly$\alpha$ EW and nebular emission-line kinematics in samples of $z\sim2$ and $z\sim3$ LBGs drawn from the literature for which matching rest-frame UV photometry, consistently measured net Ly$\alpha$ EWs, and kinematic classifications from integral field unit spectroscopy are available. We show that $z\sim2$ and $z\sim3$ LBGs segregate in colour-magnitude space according to their kinematic properties and Lyman-$\alpha$ spectral type and conclude that LBGs with Ly$\alpha$ dominant in absorption (aLBGs) are almost exclusively rotation-dominated (presumably disc-like) systems, and LBGs with Ly$\alpha$ dominant in emission (eLBGs) characteristically have dispersion-dominated kinematics. We quantify the relationship between the strength of rotational dynamic support (as measured using ${v}_{\mathrm{obs}}/2{\sigma }_{\mathrm{int}}$ and ${v}_{\mathrm{rot}}/{\sigma}_{\mathrm{0}}$) and net Ly$\alpha$ EW for subsets of our kinematic sample where these data are available, and demonstrate the consistency of our result with other properties that scale with net Ly$\alpha$ EW and kinematics. Based on these findings, we suggest a method by which large samples of rotation- and dispersion-dominated galaxies might be selected using broadband imaging in as few as three filters and/or net Ly$\alpha$ EW alone. If confirmed with larger samples, application of this method will enable an understanding of galaxy kinematic behaviour over large scales in datasets from current and future large-area and all-sky photometric surveys that will select hundreds of millions of LBGs in redshift ranges from $z\sim2-6$ across many hundreds to thousands of Mpc. Finally, we speculate that the combination of our result linking net Ly$\alpha$ EW and nebular emission-line kinematics with the known large-scale clustering behaviour of Ly$\alpha$-absorbing and Ly$\alpha$-emitting LBGs is evocative of an emergent bimodality of early galaxies that is consistent with a nascent morphology-density relation at $z\sim2-3$.

High-redshift Lyman break galaxies (LBGs) are efficiently selected in deep images using as few as three broadband filters, and have been shown to have multiple intrinsic and small- to large-scale environmental properties related to Lyman-$\alpha$. In this paper we demonstrate a statistical relationship between net Lyman-$\alpha$ equivalent width (net Ly$\alpha$ EW) and the optical broadband photometric properties of LBGs at $z\sim2$. We show that LBGs with the strongest net Ly$\alpha$ EW in absorption (aLBGs) and strongest net Ly$\alpha$ EW in emission (eLBGs) separate into overlapping but discrete distributions in $(U_n-\mathcal{R})$ colour and $\mathcal{R}$-band magnitude space, and use this segregation behaviour to determine photometric selection criteria by which sub-samples with a desired Ly$\alpha$ spectral type can be selected using data from as few as three broadband optical filters. We propose application of our result to current and future large-area and all-sky photometric surveys that will select hundreds of millions of LBGs across many hundreds to thousands of Mpc, and for which spectroscopic follow-up to obtain Ly$\alpha$ spectral information is prohibitive. To this end, we use spectrophotometry of composite spectra derived from a sample of 798 LBGs divided into quartiles on the basis of net Ly$\alpha$ EW to calculate selection criteria for the isolation of Ly$\alpha$-absorbing and Ly$\alpha$-emitting populations of $z\sim3$ LBGs using ugri broadband photometric data from the Vera Rubin Observatory Legacy Survey of Space and Time (LSST).

To explore the role environment plays in influencing galaxy evolution at high redshifts, we study $2.0\leq z<4.2$ environments using the FourStar Galaxy Evolution (ZFOURGE) survey. Using galaxies from the COSMOS legacy field with ${\rm log(M_{*}/M_{\odot})}\geq9.5$, we use a seventh nearest neighbour density estimator to quantify galaxy environment, dividing this into bins of low-, intermediate-, and high-density. We discover new high-density environment candidates across $2.0\leq z<2.4$ and $3.1\leq z<4.2$. We analyse the quiescent fraction, stellar mass and specific star formation rate (sSFR) of our galaxies to understand how these vary with redshift and environment. Our results reveal that, across $2.0\leq z<2.4$, the high-density environments are the most significant regions, which consist of elevated quiescent fractions, ${\rm log(M_{*}/M_{\odot})}\geq10.2$ massive galaxies and suppressed star formation activity. At $3.1\leq z<4.2$, we find that high-density regions consist of elevated stellar masses but require more complete samples of quiescent and sSFR data to study the effects of environment in more detail at these higher redshifts. Overall, our results suggest that well-evolved, passive galaxies are already in place in high-density environments at $z\sim2.4$, and that the Butcher–Oemler effect and SFR-density relation may not reverse towards higher redshifts as previously thought.

Active galactic nuclei (AGN) have been observed as far as redshift $z \sim 7$. They are crucial in investigating the early Universe as well as the growth of supermassive black holes at their centres. Radio-loud AGN with their jets seen at a small viewing angle are called blazars and show relativistic boosting of their emission. Thus, their apparently brighter jets are easier to detect in the high-redshift Universe. DES J014132.4–542749.9 is a radio-luminous but X-ray weak blazar candidate at $z = 5$. We conducted high-resolution radio interferometric observations of this source with the Australian Long Baseline Array at $1.7$ and $8.5$ GHz. A single, compact radio-emitting feature was detected at both frequencies with a flat radio spectrum. We derived the milliarcsecond-level accurate position of the object. The frequency dependence of its brightness temperature is similar to that of blazar sources observed at lower redshifts. Based on our observations, we can confirm its blazar nature. We compared its radio properties with those of two other similarly X-ray-weak and radio-bright AGN, and found that they show very different relativistic boosting characteristics.

Understanding properties of galaxies in the epoch of reionization (EoR) is a frontier in the modern astronomy. ALMA observations have demonstrated that i) some [O iii] 88 μm emitters have matured stellar populations at z>6, implying early star formation activity at z>10, and that ii) high-z star-forming galaxies typically have very high [O iii] 88 μm-to-[C ii] 158 μm luminosity ratios ranging from 3 to 12 or higher, indicating interstellar media of high-z galaxies could be highly ionized. We discuss initial results of a medium-sized JWST GO1 program that targets a sample of 12 z=6–8 ALMA [O iii] 88 μm emitters with NIRCam and NIRSPec IFU modes (GO-1840). Our JWST GO1 program, in conjunction with ALMA data, will characterize the stellar, nebular, and dust properties of these [O iii] 88 μm emitters and place this galaxies in the context of reionization.

We study the formation of the first massive galaxies and their observational properties for comparisons with JWST or ALMA data by performing cosmological radiative-hydrodynamics simulations. We find that galaxies in overdense regions have high star formation rates larger than ∼10 M⊙ yr-1 at z=10 and their stellar masses reach 109 M⊙. The star formation rates of our model galaxies at z≲12 nicely match recent JWST data. In addition, we show that the morphology of galaxies drastically changes with time via major mergers and stellar feedback.

While unobscured and radio-quiet active galactic nuclei are regularly being found at redshifts $z > 6$ , their obscured and radio-loud counterparts remain elusive. We build upon our successful pilot study, presenting a new sample of low-frequency-selected candidate high-redshift radio galaxies (HzRGs) over a sky area 20 times larger. We have refined our selection technique, in which we select sources with curved radio spectra between 72–231 MHz from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. In combination with the requirements that our GLEAM-selected HzRG candidates have compact radio morphologies and be undetected in near-infrared $K_{\rm s}$ -band imaging from the Visible and Infrared Survey Telescope for Astronomy Kilo-degree Infrared Galaxy (VIKING) survey, we find 51 new candidate HzRGs over a sky area of approximately $1200\ \mathrm{deg}^2$ . Our sample also includes two sources from the pilot study: the second-most distant radio galaxy currently known, at $z=5.55$ , with another source potentially at $z \sim 8$ . We present our refined selection technique and analyse the properties of the sample. We model the broadband radio spectra between 74 MHz and 9 GHz by supplementing the GLEAM data with both publicly available data and new observations from the Australia Telescope Compact Array at 5.5 and 9 GHz. In addition, deep $K_{\rm s}$ -band imaging from the High-Acuity Widefield K-band Imager (HAWK-I) on the Very Large Telescope and from the Southern Herschel Astrophysical Terahertz Large Area Survey Regions $K_{\rm s}$ -band Survey (SHARKS) is presented for five sources. We discuss the prospects of finding very distant radio galaxies in our sample, potentially within the epoch of reionisation at $z \gtrsim 6.5$ .

Long gamma-ray bursts (LGRBs) and superluminous supernovae (SLSNe) are expected to result from massive star deaths. However to date, there has been no direct observational measurement of their cloud collapse timescales nor progenitor lifetimes to help constrain their mass. Our analyses of z 2 LGRB afterglow spectra and Hubble Space Telescope images find a higher fraction of host galaxies that are interacting, have a close companion, and/or may have experienced a recent galaxy ‘fly by’ as compared to the general z 2 galaxy population. A smaller set of z 2 SLSNe suggests a similar result. Under the hypothesis that galaxy interactions induce cloud collapse and star formation near their closest approach, we explore measurements of the host and companion galaxy velocities and separations at the time of the LRGB/SLSN event as a direct physical means to measure the timescale of cloud collapse plus progenitor star lifetime.

The FAST Ultra-Deep Survey (FUDS) is a blind survey that aims for the direct detection of H i in galaxies at redshifts $z<0.42$ . The survey uses the multibeam receiver on the Five-hundred-metre Aperture Spherical Telescope (FAST) to map six regions, each of size $0.72\ \textrm{deg}^2$ at high sensitivity ( ${\sim}50\,\mu \textrm{Jy}$ ) and high-frequency resolution (23 kHz). The survey will enable studies of the evolution of galaxies and their H i content with an eventual sample size of ${\sim}1\,000$ . We present the science goals, observing strategy, the effects of radio frequency interference at the FAST site, our mitigation strategies and the methods for calibration, data reduction and imaging as applied to initial data. The observations and reductions for the first field, FUDS0, are completed, with around 128 H i galaxies detected in a preliminary analysis. Example spectra are given in this paper, including a comparison with data from the overlapping GAL2577 field of Arecibo Ultra-Deep Survey.

We present an overview of the project “The Physics of Galaxy Assembly: IFS observations of high-z galaxies”, a Guaranteed Time Observations (GTO) programme of the James Webb Space Telescope (JWST). It an ambitious project aimed at investigating the internal structure of distant galaxies with the NIRSpec integral field spectrograph (IFS), having allocated 273 hours of JWST prime time. The NIRSpec capability will provide us with spatially resolved spectroscopy in the 1-5 μm range of a sample of over forty galaxies and Active Galactic Nuclei in the redshift range 3 < z 9. IFS observations of individual galaxies will enable us to investigate in detail the most important physical processes driving galaxy evolution across the cosmic epoch. More in detail, the main specific objectives are: to trace the distribution of star formation, to map the resolved properties of the stellar populations, to trace the gas kinematics (i.e. velocity fields, velocity dispersion) and, hence, determine dynamical masses and also identify non-virial motions (outflow and inflows), and to map metallicity gradients and dust attenuation.

The Atacama Large Millimetre/Sub-millimetre Array (ALMA) is obtaining the deepest observations of early galaxies ever achieved at (sub-)millimetre wavelengths, and detecting the dust emission of young galaxies in the first billion years of cosmic history, well in the epoch of reionization. Here I review some of the latest results from these observations, with special focus on the REBELS large programme, which targets a sample of 40 star-forming galaxies at z ⋍ 7. ALMA detects significant amounts of dust in very young galaxies, and this dust might have different properties to dust in lower-redshift galaxies. I describe the evidence for this, and discuss theoretical/modelling efforts to explain the dust properties of these young galaxies. Finally, I describe two additional surprising results to come out of the REBELS survey: (i) a new population of completely dust-obscured galaxies at z ⋍ 7, and (ii) the prevalence of spatial offsets between the ultraviolet and infrared emission of UV-bright, high-redshift star-forming galaxies.

Remarkable progress has been made in the last few years in understanding the global properties of galaxies and how they evolve through cosmic time. Major focus has been given to studies of how the availability of molecular gas regulates star-forming activity and galaxy growth, the eventual quenching of star formation, and how these mechanisms evolve through cosmic time. Most of these advances have been made thanks to ALMA and the upgraded capabilities of NOEMA. In this contribution, I briey review the latest constraints on the molecular gas content based on dierent tracers of the interstellar medium (ISM; dust continuum and CO, [CI] and [CII] line emission), including recent determinations of the molecular gas fraction, gas depletion timescales, and molecular gas cosmic density provided by the recent ALMA programs out to z ∼ 7. Finally, I concentrate on recent and ongoing studies aiming to spatially and kinematically resolve the cold ISM and star formation activity down to kpc scales in galaxies out to z ∼ 6 – 7, which represent an unprecedented view of the galaxy assembly and feedback processes in the early universe.

We present an overview of the project “The Physics of Galaxy Assembly: IFS observations of high-z galaxies”, a Guaranteed Time Observations (GTO) programme of the James Webb Space Telescope (JWST). It an ambitious project aimed at investigating the internal structure of distant galaxies with the NIRSpec integral field spectrograph (IFS), having allocated 273 hours of JWST prime time. The NIRSpec capability will provide us with spatially resolved spectroscopy in the 1-5 μm range of a sample of over forty galaxies and Active Galactic Nuclei in the redshift range 3 < z < 9. IFS observations of individual galaxies will enable us to investigate in detail the most important physical processes driving galaxy evolution across the cosmic epoch. More in detail, the main specific objectives are: to trace the distribution of star formation, to map the resolved properties of the stellar populations, to trace the gas kinematics (i.e. velocity fields, velocity dispersion) and, hence, determine dynamical masses and also identify non-virial motions (outflow and inflows), and to map metallicity gradients and dust attenuation.

We processed the catalogue data for all snapshots of the Illustris TNG100 cosmological simulation and collected every calculated property of the galaxies formed at different redshifts. With this dataset we can statistically analyze parameters for galaxy samples at given redshifts, as well as trace sample parameters over the entire time range of the simulation. Focusing first on star formation rate (SFR) and metallicity, we see the cosmic star formation history with the mean maximum at around z ≈ 1.6 and the reionization bump at around z ≈ 5, while metallicity increases. For a sample of strongly star-forming galaxies with SFR > 10 M⊙ yr−1 we found different characteristics compared to the whole sample. The mean metallicity of highly star-forming galaxies is higher and changes less, and the mean SFR has its maximum at around the reionization bump.

Feedback effects by supernovae (SNe) and active galactic nuclei (AGNs) are believed to be essential for galaxy evolution and shaping present-day galaxies, but their exact mechanisms on galactic scales and their impact on CGM/IGM are not well understood yet. In galaxy formation simulations, it is still challenging to resolve sub-parsec scales, and we need to implement subgrid models to account for the physics on small scales. In this article, we summarize some of the efforts to build more physically based feedback models, discuss about pushing the resolution to its limits in galaxy simulations, testing galaxy formation codes under the AGORA code comparison project, and how to probe the impact of feedback using cosmological hydrodynamic simulations via Lyα absorption and CGM/IGM tomography technique. We also discuss our future directions of research in this field and how we make progress by comparing our simulations with observations.

Existence of the cold-mode gas accretion along with the hot-mode accretion can explain the diversity in the galactic star formation history across galaxy mass. We examine the role of various physical processes in producing the observed diversity.

Understanding the nucleosynthetic origin of nitrogen and the evolution of the N/O ratio in the interstellar medium is crucial for a comprehensive picture of galaxy chemical evolution at high-redshift because most observational metallicity (O/H) estimates are implicitly dependent on the N/O ratio. The observed N/O at high-redshift shows an overall constancy with O/H, albeit with a large scatter. We show that these heretofore unexplained features can be explained by the pre-supernova wind yields from rotating massive stars (M≳10M⊙,ν/νcrit≳0.4). Our models naturally produce the observed N/O plateau, as well as the scatter at low O/H. We find the scatter to arise from varying star formation efficiency. However, the models that have supernovae dominated yields produce a poor fit to the observed N/O at low O/H. This peculiar abundance pattern at low O/H suggests that dwarf galaxies are most likely to be devoid of SNe yields and are primarily enriched by pre-supernova wind abundances.