We exploit the continuity equation approach and the `main sequence' star-formation timescales to show that the observed high abundance of galaxies with stellar masses > a few 10^10 M_sun at redshift z>4 implies the existence of a galaxy population featuring large star formation rates (SFRs) > 10^2 M_sun/yr in heavily dust-obscured conditions. These galaxies constitute the high-redshift counterparts of the dusty star-forming population already surveyed for z<3 in the far-IR band by the Herschel space observatory. We work out specific predictions for the evolution of the corresponding stellar mass and SFR functions out to z~10, elucidating that the number density at z<8 for SFRs >30 M_sun/yr cannot be estimated relying on the UV luminosity function alone, even when standard corrections for dust extinction based on the UV slope are applied. We compute the number counts and redshift distributions (including galaxy-scale gravitational lensing) of this galaxy population, and show that current data from AzTEC-LABOCA, SCUBA-2 and ALMA-SPT surveys are already digging into it. We substantiate how an observational strategy based on a color preselection in the far-IR or (sub-)mm band with Herschel and SCUBA-2, supplemented by photometric data via on-source observations with ALMA, can allow to reconstruct the bright end of the SFR functions out to z~8. In parallel, such a challenging task can be managed by exploiting current UV surveys in combination with (sub-)mm observations by ALMA and NIKA2 and/or radio observations by SKA and its precursors.

[abridged] We investigate the coevolution of galaxies and hosted supermassive black holes throughout the history of the Universe by a statistical approach based on the continuity equation and the abundance matching technique. Specifically, we present analytical solutions of the continuity equation without source term to reconstruct the supermassive black hole (BH) mass function from the AGN luminosity functions. Such an approach includes physically-motivated AGN lightcurves tested on independent datasets, which describe the evolution of the Eddington ratio and radiative efficiency from slim- to thin-disc conditions. We nicely reproduce the local estimates of the BH mass function, the AGN duty cycle as a function of mass and redshift, along with the Eddington ratio function and the fraction of galaxies with given stellar mass hosting an AGN with given Eddington ratio. We exploit the same approach to reconstruct the observed stellar mass function at different redshift from the UV and far-IR luminosity functions associated to star formation in galaxies. These results imply that the buildup of stars and BHs in galaxies occurs via in-situ processes, with dry mergers playing a marginal role at least for stellar masses < 3 10^11 M_sun and BH masses < 10^9 M_sun, where the statistical data are more secure and less biased by systematic errors. In addition, we develop an improved abundance matching technique to link the stellar and BH content of galaxies to the gravitationally dominant dark matter component. The resulting relationships constitute a testbed for galaxy evolution models, highlighting the complementary role of stellar and AGN feedback in the star formation process. Finally, the clustering properties of BHs and galaxies are found to be in full agreement with current observations, so further validating our results from the continuity equation.

We exploit the recent, wide samples of far-infrared (FIR) selected galaxies followed-up in X rays and of X-ray/optically selected active galactic nuclei (AGNs) followed-up in the FIR band, along with the classic data on AGN and stellar luminosity functions at high redshift z>1.5, to probe different stages in the coevolution of supermassive black holes (BHs) and host galaxies. The results of our analysis indicate the following scenario: (i) the star formation in the host galaxy proceeds within a heavily dust-enshrouded medium at an almost constant rate over a timescale ~0.5-1 Gyr, and then abruptly declines due to quasar feedback; over the same timescale, (ii) part of the interstellar medium loses angular momentum, reaches the circum-nuclear regions at a rate proportional to the star formation and is temporarily stored into a massive reservoir/proto-torus wherefrom it can be promptly accreted; (iii) the BH grows by accretion in a self-regulated regime with radiative power that can slightly exceed the Eddington limit L/L_Edd< 4, particularly at the highest redshifts; (iv) for massive BHs the ensuing energy feedback at its maximum exceeds the stellar one and removes the interstellar gas, thus stopping the star formation and the fueling of the reservoir; (v) afterwards, if the latter has retained enough gas, a phase of supply-limited accretion follows exponentially declining with a timescale of about 2 e-folding times. We show that the ratio of the FIR luminosity of the host galaxy to the bolometric luminosity of the AGN maps the various stages of the above sequence. Finally, we discuss how the detailed properties and the specific evolution of the reservoir can be investigated via coordinated, high-resolution observations of starforming, strongly-lensed galaxies in the (sub-)mm band with ALMA and in the X-ray band with Chandra and the next generation X-ray instruments.

We have used new, deep, visible and near infrared observations of the compact starburst cluster in the giant HII region NGC 3603 and its surroundings with the WFC3 on HST and HAWK-I on the VLT to study in detail the physical properties of its intermediate mass (~ 1 - 3 M_sun) stellar population. We show that after correction for differential extinction and actively accreting stars, and the study of field star contamination, strong evidence remains for a continuous spread in the ages of pre-main sequence stars in the range ~ 2 to ~ 30 Myr within the temporal resolution available. Existing differences among presently available theoretical models account for the largest possible variation in shape of the measured age histograms within these limits. We also find that this isochronal age spread in the near infrared and visible Colour-Magnitude Diagrams cannot be reproduced by any other presently known source of astrophysical or instrumental scatter that could mimic the luminosity spread seen in our observations except, possibly, episodic accretion. The measured age spread and the stellar spatial distribution in the cluster are consistent with the hypothesis that star formation started at least 20-30 Myrs ago progressing slowly but continuously up to at least a few million years ago. All the stars in the considered mass range are distributed in a flattened oblate spheroidal pattern with the major axis oriented in an approximate South-East - North-West direction, and with the length of the equatorial axis decreasing with increasing age. This asymmetry is most likely due to the fact that star formation occurred along a filament of gas and dust in the natal molecular cloud oriented locally in this direction.