In the past decade hundreds of galaxies have been identified in the Epoch of Reionisation, selected from their rest-frame UV light. Only a handful of these sources, however, have spectroscopic redshift determinations and we have limited understanding of their physical properties. ALMA is currently transforming this field by providing the first view of the dust obscured star-formation, the kinematics of these sources, the cool gas traced by [CII] and highly ionised gas traced by [OIII]. In this talk I will discuss new and recent results on the UV-bright galaxy population during the first billion years of cosmic time and what they imply for their observational and physical properties.

While the high-redshift component of the CANDELS survey was designed with the z ∼ 6–8 era in mind, these data do probe the far-UV of galaxies at even higher redshift. A few studies have ventured this far out, and have published conflicting results - some continue to find significant star-formation, while others conclude there is a steep decline in this quantity. Here I report on a new search for z = 9–10 galaxies, making significant use of the Spitzer/IRAC data in the CANDELS fields. We have discovered a larger number of galaxies in this epoch than previous works, implying the UV luminosity function, and thus the SFR density, may not evolve as steeply as previously thought. This implies that star-formation begins early in the universe. I will also report on a new study searching for the earliest quenched galaxies at 3 < z < 5, which are not predicted by models, yet may exist if galaxies form very early, and thus can approach their quenching phase quicker.

Over the last few years, great progress has been made in understanding the build-up of the first generations of galaxies based on deep optical and near-infrared imaging from the Hubble Space Telescope. However, HST only samples the rest-frame UV light of galaxies at z …4, providing only limited information on the dust obscuration and on stellar masses of these sources. Fortunately, several Spitzer/IRAC programs have complemented the extragalactic HST fields with ultra-deep imaging data, allowing for a rest-frame optical view on early galaxies. Together with first ALMA/ NOEMA (sub)mm observations on distant galaxies, we are starting to gain a more and more complete picture of galaxy star-formation and mass build-up in the early universe. In this talk, I will present an overview of our current understanding of normal star-forming galaxies at z > 3 based the combination of HST+Spitzer+ALMA/NOEMA data. In particular, I will show how HST as already pushed into JWST territory with the discovery and spectroscopic confirmation of a galaxy at z = 11.1 ± 0.1, only : 400 Myr after the Big Bang. I will also highlight some of the exciting possibilities that lie ahead with JWST to push the spectroscopic frontier to the cosmic dawn and to finally probe the physics of early galaxies.

Understanding properties of galaxies in the epoch of reionization (EoR) is a frontier in the modern astronomy. With the advent of ALMA, it has become possible to detect far-infrared fine structure lines (e.g. [CII] 158 μm and [OIII] 88 μm) and dust continuum emission in star-forming galaxies in the EoR. Among these lines, our team is focusing on [OIII] 88 μm observations in high-z galaxies. After the first detection of [OIII] in the epoch of reionization (EoR) in 2016 from our team at z = 7.21, there are now more than ten [OIII] detections at z > 6 up to z = 9.11. Interestingly, high-z galaxies typically have very high [OIII]-to-[CII] luminosity ratio ranging from 3 to 12 or higher, demonstrating [OIII] is a powerful tracer at high-z. The high luminosity ratios may imply that high-z galaxies have low gas-phase metallicity and/or high ionization states.

Gravitational lensing from galaxy clusters has great potential for deriving the prevalence and physical properties of ultra-faint galaxies at early times, with recent very impressive results from the Hubble Frontier Fields program. Important issues in deriving the most accurate results are accurate constraints on source sizes and a robust treatment of uncertainties in the magnification models. Using > 3300 z = 2 – 10 galaxies behind the 6 Hubble Frontier Fields clusters and a forwards modeling approach, I describe the efforts of my collaborators and me to map out the galaxy luminosity functions at ∼ − 13 mag from z ∼ 9 to z ∼ 2, i.e, a factor of 1000 below Lå and to the typical luminosity of galaxies suspected to drive cosmic reionization. Additionally, I discuss the constraints we can obtain on the properties of faint sources, in particular their stellar masses, mass-to-light ratios, colors, and stellar population ages. I conclude with a prospective on using cluster lenses to study the distant universe with the James Webb Space Telescope.

The very brightest z > 6 galaxies are ideal laboratories for studying the physical properties of star-forming objects into the epoch of reionization. Selected from degree-scale, ground-based fields, these rare objects provide a key insight into early dust production and may harbour faint AGN. Targeted follow-up of small samples have unexpectedly shown both Lyman-alpha emission and other rest-frame UV lines (e.g CIV and HeII), suggesting unique star-formation conditions (or AGN) at early times. Furthermore, ALMA observations have revealed that 75% of the star-formation in these galaxies may be obscured. I will talk about HST/ALMA follow-up of bright z ∼ 7 LBGs in COSMOS and present new results from even brighter samples from z = 6 – 9 selected over ∼ 5 deg2. The power of both ALMA and JWST, coupled with the intrinsic luminosity of these sources, will provide a unique insight into the formation and evolution of vigorously star-forming galaxies in the first billion years.

Distant luminous Lyman-α emitters are excellent targets for detailed observations of galaxies in the epoch of reionisation. Spatially resolved observations of these galaxies allow us to simultaneously probe the emission from young stars, partially ionised gas in the interstellar medium and to constrain the properties of the surrounding hydrogen in the circumgalactic medium. We review recent results from (spectroscopic) follow-up studies of the rest-frame UV, Lyman-α and [CII] emission in luminous galaxies observed ∼500 Myr after the Big Bang with ALMA, HST/WFC3 and VLT/X-SHOOTER. These galaxies likely reside in early ionised bubbles and are complex systems, consisting of multiple well separated and resolved components where traces of metals are already present.

Ultra-deep observations of blank fields with the Hubble Space Telescope have made important inroads in characterizing galaxy populations at redshift z = 6 – 10. Gravitational lensing by massive galaxy clusters offers a new route to identify the faintest sources at the epoch of reionization. In particular, thanks to the Hubble Frontier Fields program, we robustly pushed the detection limit down to MAB = − 15 mag at z ∼ 6. I will present the latest results based on the complete dataset of the HFF clusters and parallel fields, and their implications on the ability of galaxies to reionize the Universe. I will also discuss the results of a comprehensive end-to-end modeling effort towards constraining the systematic uncertainties of the lens models, which are currently the last hurdle before extending the UV LF to fainter luminosities. Finally, I will discuss the great discoveries awaiting combination of such cosmic lenses with the upcoming James Webb Space Telescope and the exciting opportunity to probe the turnover of the UV LF, hence the limit of the star formation process at those early epochs.

Characterising primeval galaxies entails the challenging goal of observing galaxies with modest star formation rates (SFR < 100 Mȯyr−1) and approaching the beginning of the reionisation epoch (z > 6). To date a large number of primeval galaxies have been identified thanks to deep near-infrared surveys. However, to further our understanding on the formation and evolution of such primeval objects, we must investigate their nature and physical properties through multi-band spectroscopic observations. Information on dust content, metallicity, interactions with the surrounding environment, and outflows can be obtained with ALMA observations of far-infrared (FIR) lines such as the [Cii] at 158 μm and [Oiii] at 88 μm. Here, we, thus, discuss the recent results unveiled by ALMA observations and present new [Cii] observations of BDF-3299, a star-forming galaxy at z = 7.1 showing a spatial and spectral offset between the rest-frame UV and the FIR lines emission.

The opportunities offered by JWST and the ELT for the detection and study of forming/just formed globular clusters at high redshifts are illustrated, also alluding at the unique insight we may get on the very early stages of galaxy formation.

The Atacama Large (Sub-)millimeter Array (ALMA) has provided glimpse of the interstellar medium (ISM) properties of galaxies at the Epoch of Reionization (EoR); however, detailed understanding of their internal structure is still lacking. We present properties of molecular cloud complexes (MCCs) in a prototypical galaxy at this epoch studied in cosmological zoom-in simulations (Leung et al. 2019c). Typical MCC mass and size are comparable to nearby spirals and starburst galaxies (Mgas∼106.5Mȯ and R≃45–100 pc). MCCs are highly supersonic, with velocity dispersion of σgas≃20–100 km s−1 and pressure of P/kB ≃107.6Kcm−3, which are comparable to gas-rich starburst galaxies. In addition, we perform stability analysis to understand the origin and dynamical properties of MCCs. We find that MCCs are globally stable in the main disk of Althæa. Densest regions where star formation is expected to take place in clumps and cores on even smaller scales instead have lower virial parameter and Toomre-Q values. Detailed studies of the star-forming gas dynamics at the EoR thus require a spatial resolution of < 40 pc ( ≃ 0.01″), which is within reach of ALMA, to complement studies of stellar populations at EoR using the James Webb Space Telescope (JWST).

Galaxy formation in the first billion years mark a time of great upheaval in the history of the Universe: the first galaxies started both the ‘metal age’ as well as the era of cosmic reionization. I will start by reviewing the dust production mechanisms and dust masses for high-redshift galaxies which will be revolutionized in the ALMA era. I will then show how the JWST will be an invaluable experiment to shed light on the impact of reionization feedback on early galaxy formation. As we look forward towards the era of 21cm cosmology, I will highlight the crucial and urgent synergies required between 21cm facilities (such as the SKA) and galaxy experiments (JWST, E-ELT and Subaru to name a few) to understand the physics of the epoch of reionization that remains a crucial frontier in the field of astrophysics and physical cosmology. Time permitting, I will try to give a flavour of how the assembly of early galaxies, accessible with the forthcoming JWST, can provide a powerful testbed for Dark Matter models beyond ‘Cold Dark Matter’.

In this invited talk, we discuss the physics of the lifecycle of dust in the context of galaxy formation simulations. After outlining the basic physical processes, we apply algorithms for the formation, growth, and destruction of dust in the ISM to a state-of-the-art cosmological simulation to develop a model for the evolution of the dust to gas and dust to metals ratios in galaxies. We show that while modern simulations are able to match the observed dust mass function at redshift z = 0, most models underpredict the observed mass function at high-redshift (z = 2). We then show the power of these techniques by expanding our model to include a spectrum of dust sizes, and make initial predictions for extinction laws in local galaxies.

Recent observations have successfully detected UV or infrared flux from galaxies at the epoch of reionization. However, the origin of their radiative properties has not been fully understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present theoretical predictions of multi-wavelength radiative properties of the first galaxies at z = 6–15. We find that most of the gas and dust are ejected from star-forming regions due to supernova (SN) feedback, which allows UV photons to escape. We show that the peak of SED rapidly shifts between UV and infrared wavelengths on a timescale of 100 Myr due to intermittent star formation and feedback. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. In addition, we find that the escape fraction of ionizing photons also changes between 1–40% at z > 10. The mass fraction of H ii region changes with star formation history, resulting in fluctuations of metal lines and Lyman-α line luminosities. In the starbursting phase of galaxies with a halo mass ∼1011Mȯ (1012Mȯ), the simulated galaxy has L[OIII] ∼ 1042 (1043) erg s−1, which is consistent with the observed star-forming galaxies at z > 7. Our simulations suggest that deep [Cii] observation with ALMA can trace the distribution of neutral gas extending over ∼20 physical kpc. We also find that the luminosity ratio L[OIII]/L[CII] decreases with bolometric luminosity due to metal enrichment. Our simulations show that the combination of multi-wavelength observations by ALMA and JWST will be able to reveal the multi-phase ISM structure and the transition from starbursting to outflowing phases of high-z galaxies.

Cosmological hydrodynamical simulations have become an important theoretical tool for understanding the formation and evolution of the first galaxies during cosmic dawn, between redshifts 5 and 15. I will introduce the FirstLight database of about 300 zoom-in simulations with a resolution of 10 parsecs. This database agrees well with observed UV luminosity functions and stellar mass functions. I will discuss the origin and evolution of the star-forming main sequence of galaxies and the main drivers of the star formation histories at these early epochs. I will show simulated SEDs from UV to IR, including stellar and nebular emission. The rest-frame UV spectra show steep slopes and a high production efficiency of Lyman continuum photons. These properties are consistent with young stellar populations with low metallicities. Simulated recombination lines allow us to link the physical conditions of the gas around these stellar populations with observables, like equivalent widths in OIII or Hα or BPT diagrams at high-z. These simulations are making predictions that will be tested for the first time in future deep fields with the James Webb Space Telescope (JWST). I will finally discuss preliminary results involving JWST mock fields and predictions for ALMA observations by post-processing FirstLight snapshots with Powderday radiative transfer code.

We present a suite of high-resolution cosmological zoom-in simulations of galaxies at z⩾ 5using the state-of-the-art models for the multi-phase ISM, star formation, and stellar feedback from the FIRE project. We present a series of key results from these simulations, including the stellar mass–halo mass relation, the ultraviolet luminosity functions, dust attenuation and dust temperatures, the ubiquitous formation of bound star clusters, morphology and clumpiness, and the escape fractions of ionizing photons from high-redshift galaxies. We discuss how different simulations in the literature agree and disagree and what observations are most useful for testing the models in the era of ALMA and JWST.

Reionization represents an important epoch in the history in the Universe, when the first stars and galaxies gradually ionize the neutral hydrogen in the intergalactic medium (IGM). Understanding the nature of the ionizing sources, the associated ionization of the IGM, and its impact on subsequent structure formation and galaxy evolution by means of radiative feedback effects, represent key outstanding questions in current astrophysics. High-redshift galaxy observations and simulations have significantly extended our knowledge on the nature of high-redshift galaxies. However, essential properties such as the escape fraction of ionizing photons from galaxies into the IGM and their dependency on galactic properties remain essentially unknown, but determine significantly the distribution and time evolution of the ionized regions during reionization. Analyzing this ionization topology by means of the neutral hydrogen sensitive 21cm signal with radio interferometers such as SKA offers a complementary and unique opportunity to determine the nature of these first galaxies. I will show results from a self-consistent semi-numerical model of galaxy evolution and reionization, and discuss the potential of inferring galactic properties with the 21cm signal as well as the impact of reionization on the high-redshift galaxy population and its evolution.

We present an empirical model built on a high-resolution N-body dark matter simulation. We assume a redshift-independent star-formation efficiency for each halo to convert the accretion rate into a star-formation rate. Our model is calibrated using the z = 4 UV luminosity function (UVLF) and successfully predicts the observed UVLF at z = 5 – 10. We present predictions at z = 5 – 10 for UV luminosity and stellar mass functions, JWST number counts, the stellar-to-halo mass relation and star-formation histories. We combine this model with bleeding-edge reionization constraints (from z > 7 quasars, z ∼ 7 Ly α line-profiles, the updated Planck τ) to find new perspectives on the Epoch of Reionization (EoR). We find MUV < − 13.5 galaxies need an average fesc = 0.22 ± 0.05 to drive reionization and a highly compressed timeline: the IGM neutral fraction is [0.9, 0.5, 0.1] at z = [8.4 ± 0.2, 7.0 ± 0.2, 6.3 ± 0.2]. Inspired by the newly assembled sample of Lyman Continuum leakers that unanimously displays higher-than-average star-formation surface density (sigma), we fit a model tying fesc to sigma. Since sigma grows by > 2.5 dex over z = 0 – 8, our model explains the humble values of fesc at low-z. We find, strikingly, that < 5% of galaxies with MUV < − 18 account for > 80% of the reionization budget. We predict leakers like COLA1 (z = 6.6, MUV = − 21.5) become common towards the EoR and that the protagonists of reionization are not hiding across the faint-end of the luminosity function but are already known to us.

Accurate predictions of the physics of interstellar medium (ISM) are vital for understanding galaxy formation and evolution. Modelling photoionized regions with complex geometry produces realistic ionization structures within the nebulae, providing the necessary physical predictions to interpret observational data. 3D photoionization codes built with Monte Carlo techniques provide powerful tools to produce the ionizing radiation field with fractal geometry. We present a high-resolution Monte Carlo modelling of a nebula with fractal geometry, and will further show how nebular geometry influences the emission-line behaviours. Our research has important implications for studies of emission-line ratios in high redshift galaxies.

We update the Paris-Durham shock model, a state-of-the-art magnetohydrodynamic (MHD) shock code developed with a focus on molecular chemistry, in order to account for the self-generated UV field produced in shocks at velocities in the range 25-50 km/s. In these shocks there is significant excitation of atomic Hydrogen, with a large flux of Lyα photons escaping ahead of the shock to heat, ionize and drive molecular chemistry in a large slab of preshock gas.