We study the properties of star-forming galaxies selected at 610 MHz with the GMRT in a survey covering $\sim$1.86 deg$^2$ down to a noise of $\sim$7.1\,$\mu$Jy / beam. These were identified by combining multiple classification diagnostics: optical, X-ray, infrared and radio data. Of the 1685 SFGs from the GMRT sample, 496 have spectroscopic redshifts whereas 1189 have photometric redshifts. We find that the IRRC of star-forming galaxies, quantified by the infrared-to-1.4 GHz radio luminosity ratio $\rm{q_{IR}}$, decreases with increasing redshift: $\rm{q_{IR}\,=\,2.86\pm0.04(1\,+\,z)^{-0.20\pm0.02}}$ out to $z \sim 1.8$. We use the $\rm{V/V_{max}}$ statistic to quantify the evolution of the co-moving space density of the SFG sample. Averaged over luminosity our results indicate $\rm{\langle V/V_{max} \rangle}$ to be $\rm{0.51\,\pm\, 0.06}$, which is consistent with no evolution in overall space density. However we find $\rm V/V_{max}$ to be a function of radio luminosity, indicating strong luminosity evolution with redshift. We explore the evolution of the SFGs radio luminosity function by separating the source into five redshift bins and comparing to theoretical model predictions. We find a strong redshift trend that can be fitted with a pure luminosity evolution of the form $\rm{L_{610\,MHz}\,\propto\,(\,1+\,z)^{(2.95\pm0.19)-(0.50\pm0.15)z}}$. We calculate the cosmic SFR density since $\rm{z \sim 1.5}$ by integrating the parametric fits of the evolved 610\,MHz luminosity function. Our sample reproduces the expected steep decline in the star formation rate density since $\rm{z\,\sim\,1}$.

[ABRIDGED] We aim to provide a holistic view on the typical size and kinematic evolution of massive early-type galaxies (ETGs), that encompasses their high-$z$ star-forming progenitors, their high-$z$ quiescent counterparts, and their configurations in the local Universe. Our investigation covers the main processes playing a relevant role in the cosmic evolution of ETGs. Specifically, their early fast evolution comprises: biased collapse of the low angular momentum gaseous baryons located in the inner regions of the host dark matter halo; cooling, fragmentation, and infall of the gas down to the radius set by the centrifugal barrier; further rapid compaction via clump/gas migration toward the galaxy center, where strong heavily dust-enshrouded star-formation takes place and most of the stellar mass is accumulated; ejection of substantial gas amount from the inner regions by feedback processes, which causes a dramatic puffing up of the stellar component. In the late slow evolution, passive aging of stellar populations and mass additions by dry merger events occur. We describe these processes relying on prescriptions inspired by basic physical arguments and by numerical simulations, to derive new analytical estimates of the relevant sizes, timescales, and kinematic properties for individual galaxies along their evolution. Then we obtain quantitative results as a function of galaxy mass and redshift, and compare them to recent observational constraints on half-light size $R_e$, on the ratio $v/\sigma$ between rotation velocity and velocity dispersion (for gas and stars) and on the specific angular momentum $j_\star$ of the stellar component; we find good consistency with the available multi-band data in average values and dispersion, both for local ETGs and for their $z\sim 1-2$ star-forming and quiescent progenitors.

Nuclear activity and star formation play relevant roles in the early stages of galaxy formation. We aim at identifying them in high redshift galaxies by exploiting high-resolution and sensitivity X-ray and mm data to confirm their presence and relative role in contributing to the galaxy SEDs and energy budget. We present the data, model and analysis in the X-ray and mm bands for two strongly lensed galaxies, SDP.9 and SDP.11, selected in the Herschel-ATLAS catalogues as having an excess emission in the mid-IR regime at z>1.5, suggesting nuclear activity in the early stages of galaxy formation. We observed both of them in X-ray with Chandra and analyzed the high-resolution mm data available in the ALMA Science Archive for SDP9, and, by combining the information available, we reconstructed the source morphology. Both the targets were detected in the X-ray, strongly indicating the presence of highly obscured nuclear activity. High resolution ALMA observations for SDP9 in continuum and CO(6-5) spectral line allowed us to estimate the lensed galaxy redshift to a better accuracy than pre-ALMA estimates and to model the emission of the optical, mm, and X-ray band emission for this galaxy. We demonstrated that the X-ray emission is generated in the nuclear environment and it strongly support the presence of nuclear activity in this object. Hence, we identified weak nuclear activity associated with high-z galaxies with large star formation rates, useful to extend the investigation of the relationship between star formation and nuclear activity to two intrinsically less luminous, high-z star forming galaxies than was possible so far. Given our results only for two objects, they solely cannot constrain the evolutionary models, but provide us with interesting hints and set an observational path towards addressing the role of star formation and nuclear activity in forming galaxies.

The continuity equation is developed for the stellar mass content of galaxies, and exploited to derive the stellar mass function of active and quiescent galaxies over the redshift range $z\sim 0-8$. The continuity equation requires two specific inputs gauged on observations: (i) the star formation rate functions determined on the basis of the latest UV+far-IR/sub-mm/radio measurements; (ii) average star-formation histories for individual galaxies, with different prescriptions for discs and spheroids. The continuity equation also includes a source term taking into account (dry) mergers, based on recent numerical simulations and consistent with observations. The stellar mass function derived from the continuity equation is coupled with the halo mass function and with the SFR functions to derive the star formation efficiency and the main sequence of star-forming galaxies via the abundance matching technique. A remarkable agreement of the resulting stellar mass function for active and quiescent galaxies, of the galaxy main sequence and of the star-formation efficiency with current observations is found; the comparison with data also allows to robustly constrain the characteristic timescales for star formation and quiescence of massive galaxies, the star formation history of their progenitors, and the amount of stellar mass added by in-situ star formation vs. that contributed by external merger events. The continuity equation is shown to yield quantitative outcomes that must be complied by detailed physical models, that can provide a basis to improve the (sub-grid) physical recipes implemented in theoretical approaches and numerical simulations, and that can offer a benchmark for forecasts on future observations with multi-band coverage, as it will become routinely achievable in the era of JWST.

We investigate the origin, the shape, the scatter, and the cosmic evolution in the observed relationship between specific angular momentum $j_\star$ and the stellar mass $M_\star$ in early-type (ETGs) and late-type galaxies (LTGs). Specifically, we exploit the observed star-formation efficiency and chemical abundance to infer the fraction $f_{\rm inf}$ of baryons that infall toward the central regions of galaxies where star formation can occur. We find $f_{\rm inf}\approx 1$ for LTGs and $\approx 0.4$ for ETGs with an uncertainty of about $0.25$ dex, consistent with a biased collapse. By comparing with the locally observed $j_\star$ vs. $M_\star$ relations for LTGs and ETGs we estimate the fraction $f_j$ of the initial specific angular momentum associated to the infalling gas that is retained in the stellar component: for LTGs we find $f_j\approx 1.11^{+0.75}_{-0.44}$, in line with the classic disc formation picture; for ETGs we infer $f_j\approx 0.64^{+0.20}_{-0.16}$, that can be traced back to a $z<1$ evolution via dry mergers. We also show that the observed scatter in the $j_{\star}$ vs. $M_{\star}$ relation for both galaxy types is mainly contributed by the intrinsic dispersion in the spin parameters of the host dark matter halo. The biased collapse plus mergers scenario implies that the specific angular momentum in the stellar components of ETG progenitors at $z\sim 2$ is already close to the local values, in pleasing agreement with observations. All in all, we argue such a behavior to be imprinted by nature and not nurtured substantially by the environment.

We investigate the astrophysics of radio-emitting star-forming galaxies and ac- tive galactic nuclei (AGNs), and elucidate their statistical properties in the radio band including luminosity functions, redshift distributions, and number counts at sub-mJy flux levels, that will be crucially probed by next-generation radio continuum surveys. Specifically, we exploit the model-independent approach by Mancuso et al. (2016a,b) to compute the star formation rate functions, the AGN duty cycles and the conditional probability of a star-forming galaxy to host an AGN with given bolometric luminosity. Coupling these ingredients with the radio emission properties associated to star formation and nuclear activity, we compute relevant statistics at different radio frequencies, and disentangle the relative con- tribution of star-forming galaxies and AGNs in different radio luminosity, radio flux, and redshift ranges. Finally, we highlight that radio-emitting star-forming galaxies and AGNs are expected to host supermassive black holes accreting with different Eddington ratio distributions, and to occupy different loci in the galaxy main sequence diagrams. These specific predictions are consistent with current datasets, but need to be tested with larger statistics via future radio data with multi-band coverage on wide areas, as it will become routinely achievable with the advent of the SKA and its precursors.

The assessment of the relationship between radio continuum luminosity and star formation rate (SFR) is of crucial importance to make reliable predictions for the forthcoming ultra-deep radio surveys and to allow a full exploitation of their results to measure the cosmic star formation history. We have addressed this issue by matching recent accurate determinations of the SFR function up to high redshifts with literature estimates of the 1.4 GHz luminosity functions of star forming galaxies (SFGs). This was done considering two options, proposed in the literature, for the relationship between the synchrotron emission ($L_{\rm synch}$), that dominates at 1.4 GHz, and the SFR: a linear relation with a decline of the $L_{\rm synch}$/SFR ratio at low luminosities or a mildly non-linear relation at all luminosities. In both cases we get good agreement with the observed radio luminosity functions but, in the non-linear case, the deviation from linearity must be small. The luminosity function data are consistent with a moderate increase of the $L_{\rm synch}$/SFR ratio with increasing redshift, indicated by other data sets, although a constant ratio cannot be ruled out. A stronger indication of such increase is provided by recent deep 1.4 GHz counts, down to $\mu$Jy levels. This is in contradiction with models predicting a decrease of that ratio due to inverse Compton cooling of relativistic electrons at high redshifts. Synchrotron losses appear to dominate up to $z\simeq 5$. We have also updated the Massardi et al. (2010) evolutionary model for radio loud AGNs.

By means of the updated PARSEC database of evolutionary tracks of massive stars, we compute the integrated stellar light, the ionizing photon budget and the supernova rates of young simple stellar populations (SSPs), for different metallicities and IMF upper mass limits. Using CLOUDY we compute and include in the SSP spectra the neb- ular emission contribution. We also revisit the thermal and non-thermal radio emission contribution from young stars. Using GRASIL we can thus predict the panchromatic spectrum and the main recombination lines of any type of star-forming galaxy, including the effects of dust absorption and re-emission. We check the new models against the spectral energy distributions (SEDs) of selected well-observed nearby galaxies. From the best-fit models we obtain a consistent set of star formation rate (SFR) calibrations at wavelengths ranging from ultraviolet (UV) to radio. We also provide analytical calibrations that take into account the dependence on metallcity and IMF upper mass limit of the SSPs. We show that the latter limit can be well constrained by combining information from the observed far infrared, 24 {\mu}m, 33 GHz and H{\alpha} luminosities. Another interesting property derived from the fits is that, while in a normal galaxy the attenuation in the lines is significantly higher than that in the nearby continuum, in individual star bursting regions they are similar, supporting the notion that this effect is due to an age selective extinction. Since in these conditions the Balmer decrement method may not be accurate, we provide relations to estimate the attenuation from the observed 24 {\mu}m or 33 GHz fluxes. These relations can be useful for the analysis of young high redshift galaxies.

We provide an holistic view of galaxy evolution at high redshift z>4, that incorporates the constraints from various astrophysical/cosmological probes, including the estimate of the cosmic SFR density from UV/IR surveys and long GRB rates, the cosmic reionization history after the latest Planck measurements, and the missing satellites issue. We achieve this goal in a model-independent way by exploiting the SFR functions derived by Mancuso et al. (2016) on the basis of an educated extrapolation of the latest UV/far-IR data from HST/Herschel, and already tested against a number of independent observables. Our SFR functions integrated down to an UV magnitude limit M_UV<-13 (or SFR limit around 10^-2 M_sun/yr) produces a cosmic SFR density in excellent agreement with recent determinations from IR surveys and, taking into account a metallicity ceiling Z<Z_sun/2, with the estimates from long GRB rates. They also yield a cosmic reionization history consistent with that implied by the recent measurements of the Planck mission on the electron scattering optical depth tau_es~0.058; remarkably, this result is obtained under a conceivable assumption regarding the average value f_esc~0.1 of the escape fraction for ionizing photons. We demonstrate via the abundance matching technique that the above constraints concurrently imply galaxy formation to become inefficient within dark matter halos of mass below a few 10^8 M_sun; pleasingly, such a limit is also required not to run into the missing satellite issue. Finally, we predict a downturn of the galaxy luminosity function faintward of M_UV<-12, and stress that its detailed shape, as plausibly probed in the next future by the JWST, will be extremely informative on the astrophysics of galaxy formation in small halos, or even on the microscopic nature of the dark matter.

We provide a novel, unifying physical interpretation on the origin, the average shape, the scatter, and the cosmic evolution for the main sequences of starforming galaxies and active galactic nuclei at high redshift z $\gtrsim$ 1. We achieve this goal in a model-independent way by exploiting: (i) the redshift-dependent SFR functions based on the latest UV/far-IR data from HST/Herschel, and re- lated statistics of strong gravitationally lensed sources; (ii) deterministic evolutionary tracks for the history of star formation and black hole accretion, gauged on a wealth of multiwavelength observations including the observed Eddington ratio distribution. We further validate these ingredients by showing their consistency with the observed galaxy stellar mass functions and AGN bolometric luminosity functions at different redshifts via the continuity equation approach. Our analysis of the main sequence for high-redshift galaxies and AGNs highlights that the present data are consistently interpreted in terms of an in situ coevolution scenario for star formation and black hole accretion, envisaging these as local, time coordinated processes.

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.

We have worked out predictions for the radio counts of star-forming galaxies down to nJy levels, along with redshift distributions down to the detection limits of the phase 1 Square Kilometer Array MID telescope (SKA1-MID) and of its precursors. Such predictions were obtained by coupling epoch dependent star formation rate (SFR) functions with relations between SFR and radio (synchrotron and free-free) emission. The SFR functions were derived taking into account both the dust obscured and the unobscured star-formation, by combining far-infrared (FIR), ultra-violet (UV) and H_alpha luminosity functions up to high redshifts. We have also revisited the South Pole Telescope (SPT) counts of dusty galaxies at 95\,GHz performing a detailed analysis of the Spectral Energy Distributions (SEDs). Our results show that the deepest SKA1-MID surveys will detect high-z galaxies with SFRs two orders of magnitude lower compared to Herschel surveys. The highest redshift tails of the distributions at the detection limits of planned SKA1-MID surveys comprise a substantial fraction of strongly lensed galaxies. We predict that a survey down to 0.25 microJy at 1.4 GHz will detect about 1200 strongly lensed galaxies per square degree, at redshifts of up to 10. For about 30% of them the SKA1-MID will detect at least 2 images. The SKA1-MID will thus provide a comprehensive view of the star formation history throughout the re-ionization epoch, unaffected by dust extinction. We have also provided specific predictions for the EMU/ASKAP and MIGHTEE/MeerKAT surveys.

We have combined determinations of the epoch-dependent star formation rate (SFR) function with relationships between SFR and radio (synchrotron and free-free) emission to work out detailed predictions for the counts and the redshift distributions of star-forming galaxies detected by planned Square Kilometer Array (SKA) surveys. The evolving SFR function comes from recent models fitting the far-infrared (FIR) to millimeter-wave luminosity functions and the ultraviolet (UV) luminosity functions up to z=10, extended to take into account additional UV survey data. We used very deep 1.4 GHz number counts from the literature to check the relationship between SFR and synchrotron emission, and the 95 GHz South Pole Telescope (SPT) counts of dusty galaxies to test the relationship between SFR and free-free emission. We show that the SKA will allow us to investigate the SFRs of galaxies down to few Msun/yr up to z=10, thus extending by more than two orders of magnitude the high-z SFR functions derived from Herschel surveys. SKA1-MID surveys, down to microJy levels, will detect hundreds of strongly lensed galaxies per square degree; a substantial fraction of them will show at least two images above the detection limits.