I will review recent developments in the modeling of high-redshift galaxy spectra, focusing in particular on the rest-frame ultraviolet and optical emission from young stellar populations and the interstellar medium.

In the quest to study early star-formation physics in the universe, one of the most sought after tracers is HeIIλ1640, with its presence in the lack of other metal emission/absorption lines generally being interpreted as evidence for metal-poor stellar populations. HeII ionizing photons are produced via sources of hard ionizing radiation and requires photons with energies ⩾ 54.4eV, however, traditional stellar population models lack sufficient ionising photons to match with current observations. Our analysis of z = 2 – 4 HeIIλ1640 emitters from deep 10-30h pointings from MUSE has shown that ISM properties inferred from multiple rest-UV diagnostics are not compatible with requirements necessary to reproduce HeIIλ1640 equivalent-widths. Thus, we have used latest generation of single, rotational, and binary stellar population models with realistic dust physics to explore rest-UV emission line diagnostics and link with H and He+ ionisation photon production efficiencies (ξion (H,He+)) in a variety of stellar/gas metallicities and star-formation histories. I will discus our latest results and show that including ‘exotic’ stellar phenomena such as extreme low-metallicity binary stars, X-ray binaries, and dust dissociation physics may be necessary to lessen the tension between models and observations.

Understanding the ionizing spectrum of low-metallicity galaxies is of great importance for modeling and interpreting emission line observations of early/distant galaxies.

Although a wide suite of stellar evolution, atmosphere, population synthesis, and photoionization models, taking many physical processes into account now exist, all models face a common problem: the inability to explain the presence of nebular Heii emission, which is observed in many low metallicity galaxies, both in UV and optical spectra. Several possible explanations have been proposed in the literature, including Wolf-Rayet (WR) stars, binaries, very massive stars, X-ray sources, or shocks. However, none has so far been able to explain the major observations.

We briefly discuss the He ii problem, available empirical data, and observed trends combining X-ray, optical and other studies. We present a simple and consistent physical model showing that X-ray binaries could explain the long-standing nebular Heii problem. Our model, described in Schaerer et al. (2019), successfully explains the observed trends and strength of nebular He ii emission in large samples of low metallicity galaxies and in individual galaxies, which have been studied in detail and with multi-wavelength observations. Our results have in particular important implications for the interpretation of galaxy spectra in the early Universe, which will be obtained with upcoming and future facilities.

Observations of star-forming galaxies in the distant Universe have confirmed the importance of massive stars in shaping galaxy emission and evolution. Distant stellar populations are unresolved, and the limited data available must be interpreted in the context of stellar population models. Understanding these populations, and their evolution with age and heavy element content is key to interpreting processes such as supernovae, cosmic reionization and the chemical enrichment of the Universe. With the upcoming launch of JWST and observations of galaxies within a billion years of the Big Bang, the uncertainties in modelling massive stars - particularly their interactions with binary companions - are becoming increasingly important to our interpretation of the high redshift Universe. In turn, observations of distant stellar populations provide ever stronger tests against which to gauge the success of, and flaws in, current massive star models. Here we briefly review the current status binary stellar population synthesis.

The birth of stars and the formation of galaxies are cornerstones of modern astrophysics. While much is known about how galaxies globally and their stars individually form and evolve, one fundamental property that affects both remains elusive. This is problematic because this key property, the stellar initial mass function (IMF), is a key tracer of the physics of star formation that underpins almost all of the unknowns in galaxy and stellar evolution. It is perhaps the greatest source of systematic uncertainty in star and galaxy evolution. The past decade has seen a growing number and variety of methods for measuring or inferring the shape of the IMF, along with progressively more detailed simulations, paralleled by refinements in the way the concept of the IMF is applied or conceptualised on different physical scales. This range of approaches and evolving definitions of the quantity being measured has in turn led to conflicting conclusions regarding whether or not the IMF is universal. Here I summarise the growing wealth of approaches to our understanding of this fundamental property that defines so much of astrophysics, and highlight the importance of considering potential IMF variations, reinforcing the need for measurements to quantify their scope and uncertainties carefully. I present a new framework to aid the discussion of the IMF and promote clarity in the further development of this fundamental field.

Galaxies are complicated physical systems which obey complex scaling relationships; as a result, properties measured from broadband photometry are often highly correlated, degenerate, or both. Therefore, the accuracy of basic properties like stellar masses and star formation rates (SFRs) depend on the accuracy of many second-order galaxy properties, including star formation histories (SFHs), stellar metallicities, dust properties, and many others. Here, we re-assess measurements of galaxy stellar masses and SFRs using a 14-parameter physical model built in the Prospector Bayesian inference framework. We find that galaxies are ∼0.2 dex more massive and have ∼0.2 dex lower star formation rates than classic measurements. These measurements lower the observed cosmic star formation rate density and increase the observed buildup of stellar mass, finally bringing these two metrics into agreement at the factor-of-two level at 0.5 < z < 2.5.

This article presents an up-dated analysis of synthetic optical and UV emission lines of simulated galaxies over cosmic time. The strong emission lines are derived from self-consistently coupling novel spectral models accounting for nebular emission from young stars, AGN and Post-AGB stars to cosmological zoom-in as well as large-scale simulations. Investigating the evolution of optical line-ratios in the BPT diagrams, the simulations can successfully reproduce the observed trend of [OIII]/Hβ ratio increasing from low to high redshifts, due to evolving star formation rate and gas metallicity. Standard selection criteria in the BPT diagrams can appropriately distinguish the main ionising source(s) of galaxies at low redshifts, but they are less reliable for metal-poor galaxies, dominating the early Universe. To robustly classify the ionising radiation of such metal-poor galaxies, diagnostic diagrams based on luminosity ratios of UV lines are discussed. The novel interface between simulations and observations is potentially important for the interpretation of high-quality spectra of very distant galaxies to be gathered by next-generation telescopes, such as the James Webb Space Telescope.

Despite significant progress over the past decades, all state-of-the-art population synthesis (PS) codes suffer from deficiencies limiting their potential of gaining sharp insights into the star formation history (SFH) and Chemical Enrichment History (CEH) of star-forming galaxies, i.e. the neglect of nebular continuum and, the lack of a mechanism to ensure consistency between the best-fitting SFH and the observed nebular characteristics (ONC; Balmer-lines, Balmer/Paschen jumps). These introduce biases in their recovered physical properties (stellar mass M* and sSFR). FADO is a novel self-consistent PS code employing genetic optimization, publicly available (http://www.spectralsynthesis.org), capable of identifying the SFH & CEH that reproduce the ONC of a galaxy, alleviating degeneracies in the spectral fits. The current version of FADO (v1.b) uses standard BPT emission-line ratios for the classification of low redshift (z) galaxies. Whereas this permits a reliable distinction between star-forming, Composite, Seyfert and LINERs, it is inapplicable to many intermediate-z galaxies. We present an adaptation of FADO (version v1.c) to classify higher z galaxies employing the “Blue Diagram” (e.g., Lamareille 2010) for which the most prominent blue emission-lines (< [OIII]5007Å) are observable while the Hα and [NII] are inaccessible. FADO v1.c was applied to synthetic spectra simulating the evolution of galaxies formed at higher-z with different SFHs. FADO can recover the physical and evolutionary properties of galaxies, such as M* and mean age/metallicity, with an accuracy significantly better than purely-stellar codes. An outline of FADO v1.c and applications to local and intermediate-z galaxies will be presented.

The mass-SFR relation of galaxies encodes information of present and historical star formation in the galaxy population. We expect the intrinsic scatter in the relation to increase to low mass where SFR becomes more stochastic. Measurements at z ‰ 4 from the Hubble Frontier fields have hinted at this (Santini et al., 2017), however, with the added uncertainty of lensing magnification we await JWST to provide robust measurements. Even with data-sets provided by JWST, uncertainties on mass and SFR estimates are often large, potentially covariant and dependent on assumptions used. I will present our method of Bayesian hierarchical modelling of the mass-SFR relation that self-consistently propagates uncertainties on mass and SFR estimates to uncertainties on the mass-SFR relation parameters. I will expose the biases imposed by standard SED-modelling practices, and address to what significance we can measure an increase in intrinsic scatter to low masses with JWST.

Current observations suggest an accelerated evolution of the cosmic star formation rate density for 8 < z < 10, indicating that galaxy assembly experienced an extremely intense phase during the first ∼600Myr years of cosmic time. We performed a systematic search of ultrabright star-forming galaxies at z ∼ 8 over the COSMOS/UltraVISTA survey, identifying 16 candidate Lyman-break galaxies. The still large uncertainties on the associated volume density do not yet allow us to ascertain whether a different star-formation efficiency (SFE) existed at early cosmic epochs. Leveraging the deepest Spitzer/IRAC data available from the GREATS program over the CANDELS/GOODS fields, we also constructed stacked SEDs of sub- L* LBGs at z ∼ 8. We find extreme nebular line emission (EW0 (Hα) ∼ 1000Å), high specific star-formation rates (∼10/Gyr) and indication of an inverse Balmer break. These results point toward very young ages (<100 Myr), and, combined with measurements at lower redshifts, that the SFE evolved only marginally during the first ∼1.5Gyr of cosmic history.

Spectral analysis is nowadays a widely used tool to investigate the evolution of galaxies. Assessing the reliability of this approach is crucia, motivating a through analysis. In this poster, a comparative study between two widely tools, FADO and STECKMAP, is performed, focusing on the discrepancies between the different approaches. Both codes use different methods to extract the best fit, allowing the possibility to disentangle possible biases introduced in the analysis. Our analysis showed that where nebular emission is not negligible, the results obtained with methods taking into account such a component are more reliable, and this can be very important when moving at higher redshift, where stellar populations are younger. In particular, this is true for starburst systems, where a huge amount of stars are forming almost at the same epoch. This is an important aspect to take into account the future survey, as JWST for example, which will provide the community with medium resolution spectra of galaxies at redshift 3-4 and even higher.

To understand how the nature of the ionizing sources and the leakage of ionizing photons in high-redshift galaxies can be constrained from their emission-line spectra, we compare emission-line models of star-forming galaxies including leakage of ionizing radiation, active galactic nuclei (AGN) and radiative shocks, with observations of galaxies at various redshifts with properties expected to approach those of primeval galaxies.

The number of quasars known within the first billion years of the universe (z > 6) has increased significantly over the last five years. Many of these recently discovered quasars are ideal targets for observatories in the southern hemisphere such as ALMA. I will review the current status of the highest-redshift quasars and their environments, highlighting main achievements and limitations. I will then discuss how synergistic JWST/ALMA observations will shed light onto the properties and formation of some of the most extreme environments in the early universe.

I will present results from our on-going large area survey of high-redshift quasars, which has discovered more than 20 new quasars at z > 6.5, at the epoch of reionization, forming the first large statistical sample of EoR quasars. I will discuss the rapid evolution of quasar density at that epoch, which suggests that we are witnessing the emergence of the first supermassive black hole population. I will also present multiwavelength followup observation results, especially from ALMA and Chandra, which reveals a diverse environment of quasar activities and yields new insights into the supermassive black hole/massive galaxy co-evolution.

Luminous quasars are powered by accretion onto supermassive black holes. Such luminous quasars have been discovered up to the highest redshifts, z > 7. Here we discuss recent observations of the host galaxies of luminous quasars at z ≳ 6. We do not find a correlation between ongoing black hole growth and star-formation rate in the high redshift quasars, possibly indicating that black holes and their hosts do not co-evolve. We further show that even with high spatial resolution observations of the gas kinematics, dynamical mass estimates remain highly uncertain and should be used with caution.

We leverage new ultra-deep, high resolution, multi-frequency radio imaging at 6 and 3 GHz with the unique datasets available in the GOODS-S/HUDF region in order to assess the AGN fraction in a faint radio-selected sample. For AGN identification, we adopt a multi-wavelength approach, combining X-ray and (mid-)infrared (IR) selections with radio identification such as X-ray to radio excess, flat radio spectral slopes, and the radio-IR correlation. We identify AGN in 43% of our radio sample, yielding an AGN source density of ∼ 1 arcmin−2. This AGN fraction is likely underestimated, as 1) our shallower 3 GHz data is biased against flat radio spectrum sources and 2) all of our selections may be biased against the most heavily obscured AGN. The James Webb Space Telescope’s Mid-Infrared Instrument (MIRI) will address the latter issue and we briefly outline our Cycle 1 Guaranteed Time Observation (GTO) program to search for heavily obscured AGN.

I will present evidence of QSO-driven outflows in the early Universe, resulting from the stacking analysis of ALMA observations of the [CII] emission line for a sample of 50 QSOs at z ∼ 5 – 7. The high sensitivity reached by our analysis allowed us to assess that very broad wings are on average present in the [CII] profile, and extend beyond velocities of 1000 km/s in systems with low and high SFR. Such wings are therefore tracing QSO-driven [CII] outflows, with associated mass outflow rates of 100 – 200Mȯ / yr. I will discuss how these outflows relate to those observed in lower-z AGNs and give an estimate of their spatial extent. Furthermore, I will focus on the high-resolution ALMA observation of a hyper luminous QSO at z ∼ 4.5, revealing an exceptional overdensity with multiple companions as close as 2 kpc. These crowded surroundings, and the QSO host galaxy itself, are forming stars at a very high rate (hundreds of Mȯ / yr), suggesting that a significant fraction of the stellar mass assembly at early epochs might have taken place in the companions. I will discuss how the BH and host-galaxy masses are growing in this multi- source system, which likely represents the cradle of what would be a giant galaxy at z = 0.

We present ALMA [CII] line and far-infrared (FIR) continuum observations of seven z > 6 low-luminosity quasars (M1450 > −25 mag) discovered by our on-going Subaru Hyper Suprime-Cam survey. The [CII] line was detected in all targets with luminosities of ∼(2−10) × 108 L⊙, about one order of magnitude smaller than optically luminous quasars. Also found was a wide scatter of FIR continuum luminosity, ranging from LFIR < 1011L⊙ to ∼2 × 1012L⊙. With the [CII]-based dynamical mass, we suggest that a significant fraction of low-luminosity quasars are located on or even below the local Magorrian relation, particularly at the massive end of the galaxy mass distribution. This is a clear contrast to the previous finding that luminous quasars tend to have overmassive black holes relative to the relation. Our result is expected to show a less-biased nature of the early co-evolution of black holes and their host galaxies.

Several scenarios have been proposed to describe the physical connection between galaxies and their central active galactic nuclei (AGN). This connection could act on a range of spatial scales and vary across cosmic time. In these proceedings, we consider black hole and galaxy growth and whether that growth is affected by AGN feedback both based on statistical approaches – which reveal general population trends – and based on an individual case study – which gives us a more detailed insight on the physical processes at play. For the statistical approach, we showcase a low-redshift (0.04 < z < 0.2) SDSS sample with AGN classification based on a combination of emission-line diagnostic diagrams, and for which we account for sample selection by using a V/Vmax approach. The trends on the star formation rate - stellar mass (SFR – M*) plane suggest that the most likely connection is a common gas reservoir for star formation and AGN, and that they both decline as the gas reservoir is consumed. The trends established at low-redshift could act as a local benchmark against which to compare higher redshift studies. As a complementary approach, we use a detailed case study of a nearby AGN host with integral field spectroscopy from the VLT/MUSE instrument in order to spatially resolve the interplay between AGN feedback and the host galaxy. We find that the galaxy substructure likely plays a role by collimating and/or obscuring the outflows and radiation from the central engine. Ongoing and future work with 3D spectroscopy will enable us to learn more about galaxy and black hole coevolution. Lastly, we briefly discuss lessons learnt from both approaches.

Galaxies and nuclei in dense environment at high redshift provide a good laboratory to investigate accelerated, most extreme evolution of galaxies. The SSA22 proto-cluster at z = 3.1 is known to have a three-dimensional 50 (comoving) Mpc-scale filamentary structure, traced by Lyα emitters, which makes the field a suitable target in this regard. To identify dust-obscured star-formation, a contiguous 20 arcmin2 region at the node of the cosmic structure was observed in ALMA band 6. In total 57 ALMA sources have been identified above 5σ, which makes the field one of the richest field in ALMA-identified (sub)millimeter galaxies. The follow-up spectroscopy confirmed about 20 sources as exact proto-cluster members so far. Together with high X-ray AGN fraction, our results suggest that the vigorous star formation activity and the growth of super massive black holes occurred simultaneously in the densest regions at z ∼ 3.