We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA Band-7, as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE) in COSMOS and ECDFS. These sources have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z=6. We have looked for counterparts in all the available multi-wavelength and photometric redshift catalogues, and in deeper near- and mid-IR source lists and maps, to identify optically dark sources with no matches in the public catalogues. Our ALMA blind survey allows us to push further the study of the nature and evolution of dusty galaxies at high-z, identifying luminous and massive sources to redshifts and faint luminosities never probed before by any far-infrared surveys. The ALPINE data are the first ones to sample the faint-end of the infrared LF, showing little evolution from z=2.5 to z=6, and a flat slope up to the highest redshifts. The SFRD obtained by integrating the luminosity function remains almost constant between z=2 and 6, and significantly higher than the optical/UV derivations, showing an important contribution of dusty galaxies and obscured star formation up to high-z. About 16. per cent of the ALPINE serendipitous continuum sources are optically+near-IR dark (7 show a counterpart only in the mid-IR and no HST or near-IR identification, while 2 are detected as [CII] emitters at z=5). The 7 HST and near-IR dark galaxies with mid-IR counterpart contribute for about 15 per cent of the total SFRD at z=5 and dominate the high-mass end of the stellar mass function at z>3.

We present a new modeling of the X-ray luminosity function (XLF) of Active Galactic Nuclei (AGN) out to z$\sim$3, dissecting the contribution of main-sequence (MS) and starburst (SB) galaxies. For each galaxy population, we convolved the observed galaxy stellar mass (M$_{\star}$) function with a grid of M$_{\star}$-independent Eddington ratio ($\lambda_{\rm EDD}$) distributions, normalised via empirical black hole accretion rate (BHAR) to star formation rate (SFR) relations. Our simple approach yields an excellent agreement with the observed XLF since z$\sim$3. We find that the redshift evolution of the observed XLF can only be reproduced through an intrinsic flattening of the $\lambda_{\rm EDD}$ distribution, and with a positive shift of the break $\lambda^{*}$, consistent with an anti-hierarchical behavior. The AGN accretion history is predominantly made by massive (10$^{10}<$M$_{\star}<$10$^{11}$ M$_{\odot}$) MS galaxies, while SB-driven BH accretion, possibly associated with galaxy mergers, becomes dominant only in bright quasars, at $\log$(L$_{\rm X}$/erg s$^{-1}$)$>$44.36 + 1.28$\cdot$(1+z). We infer that the probability of finding highly-accreting ($\lambda_{\rm EDD}>$ 10%) AGN significantly increases with redshift, from 0.4% (3.0%) at z=0.5 to 6.5% (15.3%) at z=3 for MS (SB) galaxies, implying a longer AGN duty cycle in the early Universe. Our results strongly favor a M$_{\star}$-dependent ratio between BHAR and SFR, as BHAR/SFR $\propto$ M$_{\star}^{0.73[+0.22,-0.29]}$, supporting a non-linear BH buildup relative to the host. Finally, this framework opens potential questions on super-Eddington BH accretion and different $\lambda_{\rm EDD}$ prescriptions for understanding the cosmic BH mass assembly.

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

We study the efficiency of galactic feedback in the early Universe by stacking the [C II] 158 um emission in a large sample of normal star-forming galaxies at 4 < z < 6 from the ALMA Large Program to INvestigate [C II] at Early times (ALPINE) survey. Searching for typical signatures of outflows in the high-velocity tails of the stacked [C II] profile, we observe (i) deviations from a single-component Gaussian model in the combined residuals and (ii) broad emission in the stacked [C II] spectrum, with velocities of |v|<~ 500 km/s. The significance of these features increases when stacking the subset of galaxies with star formation rates (SFRs) higher than the median (SFRmed = 25 Msun/yr), thus confirming their star-formation-driven nature. The estimated mass outflow rates are comparable to the SFRs, yielding mass-loading factors of the order of unity (similarly to local star-forming galaxies), suggesting that star-formation-driven feedback may play a lesser role in quenching galaxies at z > 4. From the stacking analysis of the datacubes, we find that the combined [C II] core emission (|v|< 200 km/s) of the higher-SFR galaxies is extended on physical sizes of ~ 30 kpc (diameter scale), well beyond the analogous [C II] core emission of lower-SFR galaxies and the stacked far-infrared continuum. The detection of such extended metal-enriched gas, likely tracing circumgalactic gas enriched by past outflows, corroborates previous similar studies, confirming that baryon cycle and gas exchanges with the circumgalactic medium are at work in normal star-forming galaxies already at early epochs.

We investigate a sample of 152 dusty sources at 1:5 < z < 2:5 to understand the connection of enhanced Star-Formation-Rate (SFR) and Black-Hole-Accretion-Rate (BHAR). The sources are Herschel-selected, having stellar masses M*>10^10Msun and SFR (100-1000Msun/yr) elevated(>4?) above the star-forming "main sequence", classifying them as Starbursts (SB). Through a multiwavelength fitting approach (including a dusty torus component), we divided the sample into active SBs (dominated by an AGN emission, SBs-AGN, ? 23% of the sample) and purely star-forming SBs (SBs-SFR). We visually inspected their HST/UV-restframe maps: SBs-SFR are generally irregular and composite systems; ? 50% of SBs-AGN are instead dominated by regular compact morphologies. We then found archival ALMA continuum counterparts for 33 galaxies (12 SBs-AGN and 21 SBs-SFR). For these sources we computed dust masses, and, with standard assumptions, we also guessed total molecular gas-masses. SBs turn to be gas rich systems (fgas = Mgas=Mgas/(Mgas+M*) \sim 20%-70%), and the gas fractions of the two SB classes are very similar (fgas = 43 +/-4% and fgas = 42+/- 2%). Our results show that SBs are consistent with a mixture of: 1) highly star-forming merging systems (dominating the SBs-SFR), and 2) primordial galaxies, rapidly growing their M* together with their Black Hole (mainly the more compact SBs-AGN). Anyway, feedback effects have not reduced their fgas yet. Indeed, SBs at z = 2, with relatively low bolometric AGN luminosities in the range 10^44 < Lbol(AGN) < 10^46 erg/s (compared to bright optical and X-ray quasars), are still relatively far from the epoch when the AGN feedback will quench the SFR in the host and will substantially depress the gas fractions.

Based on optical/near-IR Magellan FIRE spectra of 25 starburst galaxies at 0.5 < z < 0.9, Calabr\`o et al.(2018) showed that their attenuation properties can be explained by a single-parameter sequence of total obscurations ranging from A(V)=2 to A(V)=30 towards the starburst core centers in a mixed stars and dust configuration. We investigate here the origin of this sequence for the same sample. We show that total attenuations anti-correlate with the starburst sizes in radio (3 GHz) with a significance larger than 5sigma and a scatter of 0.26 dex. More obscured and compact starbursts also show enhanced N2 (=[NII]/Halpha) ratios and larger line velocity widths that we attribute to an increasing shock contribution toward later merger phases, driven by deeper gravitational potential wells at the coalescence. Additionally, the attenuation is also linked to the equivalent width (EW) of hydrogen recombination lines, which is sensitive to the luminosity weighted age of the relatively unobscured stellar populations. Overall, the correlations among A(V), radio size, line width, N2 and EW of Balmer/Paschen lines converge towards suggesting an evolutionary sequence of merger stages: all of these quantities are likely to be good time-tracers of the merger phenomenon, and their large spanned range appears to be characteristic of the different merger phases. Half of our sample at higher obscurations have radio sizes approximately 3 times smaller than early type galaxies at the same redshift, suggesting that, in analogy with local Ultraluminous Infrared galaxies (ULIRGs), these cores cannot be directly forming elliptical galaxies. Finally, we detect mid-IR AGN torus for half of our sample and additional X-ray emission for 6 starbursts; intriguingly, the latter have systematically more compact sizes, suggestive of emerging AGNs towards later merger stages, possibly precursors of a later QSO phase.

We obtained optical/near-IR rest-frame Magellan FIRE spectra (including Pa$\beta$ and Pa$\gamma$) of 25 starburst galaxies at 0.5<z<0.9, with average star formation rates (SFR) x7 above the Main Sequence (MS). We find that Paschen-to-Balmer line ratios saturate around a constant value corresponding to $A_{\rm V}\sim$2-3 mag, while line to IR luminosity ratios suggest a large range of more extreme obscurations and appear to be uncorrelated to the former. This behavior is not consistent with standard attenuation laws derived for local and distant galaxies, while being remarkably consistent with observations of starburst cores in which young stars and dust are homogeneously mixed. This model implies $A_{\rm V}=$2-30 mag attenuation to the center of starburst cores, with a median of ~9 mag (a factor of 4000). X-ray hardness ratios for 6 AGNs in our sample and column densities derived from observed dust masses and radio sizes independently confirm this level of attenuation. In these conditions observed optical/near-IR emission comes from surface regions, while inner starburst cores are invisible. We thus attribute the high [NII]/H$\alpha$ ratios to widespread shocks from accretion, turbulence and dynamic disturbances rather than to AGNs. The large range of optical depths demonstrates that substantial diversity is present within the starburst population, possibly connected to different merger phases or progenitor properties. The majority of our targets are, in fact, morphologically classified as mergers. We argue that the extreme obscuration provides in itself smoking gun evidence of their merger origin, and a powerful tool for identifying mergers at even higher redshifts.

For a sample of star forming galaxies in the redshift interval 0.15$<$z$<$0.3, we study how both the relative strength of the AGN infra-red emission, compared to that due to the star formation (SF), and the numerical fraction of AGNs, change as a function of the total stellar mass of the hosting galaxy group (M$^{*}_{\mathrm{group}}$), between $10^{10.25}$ and $10^{11.9}$M$_{\odot}$. Using a multi-component SED fitting analysis, we separate the contribution of stars, AGN torus and star formation to the total emission at different wavelengths. This technique is applied to a new multi-wavelength data-set in the SIMES field (23 not redundant photometric bands), spanning the wavelength range from the UV (GALEX) to the far-IR (Herschel) and including crucial AKARI and WISE mid-IR observations (4.5 \mu m$<\lambda<$24 \mu m), where the BH thermal emission is stronger. This new photometric catalog, that includes our best photo-z estimates, is released through the NASA/IPAC Infrared Science Archive (IRSA). Groups are identified through a friends of friends algorithm ($\sim$62% purity, $\sim$51% completeness). We identified a total of 45 galaxies requiring an AGN emission component, 35 of which in groups and 10 in the field. We find BHAR$\propto ($M$^{*}_{\mathrm{group}})^{1.21\pm0.27}$ and (BHAR/SFR)$\propto ($M$^{*}_{\mathrm{group}})^{1.04\pm0.24}$ while, in the same range of M$^{*}_{\mathrm{group}}$, we do not observe any sensible change in the numerical fraction of AGNs. Our results indicate that the nuclear activity (i.e. the BHAR and the BHAR/SFR ratio) is enhanced when galaxies are located in more massive and richer groups.