Historically, GHz radio emission has been used extensively to characterize the star-formation activity in galaxies. In this work, we look for empirical relations amongst the radio luminosity, the infrared luminosity, and the CO-based molecular gas mass. We assemble a sample of 278 nearby galaxies with measurements of radio continuum and total infrared emission, and the $^{12}$CO (J = 1-0) emission line. We find a correlation between the radio continuum and the CO emission line (with a scatter of 0.36 dex), in a large sample of different kind of galaxies. Making use of this correlation, we explore the evolution of the molecular gas mass function and the cosmological molecular gas mass density in six redshift bins up to $z = 1.5$. These results agree with previous semi-analytic predictions and direct measurements: the cosmic molecular gas density increases up to $z=1.5$. In addition, we find a single plane across five orders of magnitude for the explored luminosities, with a scatter of 0.27 dex. These correlations are sufficiently robust to be used for samples where no CO measurements exist.

We present the highest resolution --- 15 pc (0.03'') --- ALMA $^{12}$CO(2-1) line emission and 1.3mm continuum maps, tracers of the molecular gas and dust, respectively, in the nearby merging galaxy system NGC 6240, that hosts two supermassive black holes growing simultaneously. These observations provide an excellent spatial match to existing Hubble optical and near-infrared observations of this system. A significant molecular gas mass, $\sim$9$\times$10$^9$M$_\odot$, is located in between the two nuclei, forming a clumpy stream kinematically dominated by turbulence, rather than a smooth rotating disk as previously assumed from lower resolution data. Evidence for rotation is seen in the gas surrounding the southern nucleus, but not in the northern one. Dynamical shells can be seen, likely associated with nuclear supernovae remnants. We further detect the presence of significant high velocity outflows, some of them reaching velocities $>$500 km/s, affecting a significant fraction, $\sim$11\% of the molecular gas in the nuclear region. Inside the spheres of influence of the northern and southern supermassive black holes we find molecular masses of 7.4$\times$10$^8$M$_\odot$ and 3.3$\times$10$^9$M$_\odot$, respectively. We are thus directly imaging the reservoir of gas that can accrete onto each supermassive black hole. These new ALMA maps highlight the critical need for high resolution observations of molecular gas in order to understand the feeding of supermassive black holes and its connection to galaxy evolution in the context of a major galaxy merger.

Using the IRAM NOrthern Extended Millimeter Array (NOEMA), we conducted a program to measure redshifts for 13 bright galaxies detected in the Herschel Astrophysical Large Area Survey (H-ATLAS) with $S_{500{\mu}\rm m}\ge$80 mJy. We report reliable spectroscopic redshifts for 12 individual sources, which are derived from scans of the 3 and 2 mm bands, covering up to 31 GHz in each band, and are based on the detection of at least two emission lines. The spectroscopic redshifts are in the range $2.08<z<4.05$ with a median value of $z=2.9\pm$0.6. The sources are unresolved or barely resolved on scales of 10 kpc. In one field, two galaxies with different redshifts were detected. In two cases the sources are found to be binary galaxies with projected distances of ~140 kpc. The linewidths of the sources are large, with a mean value for the full width at half maximum of 700$\pm$300 km/s and a median of 800 km/s. We analyse the nature of the sources with currently available ancillary data to determine if they are lensed or hyper-luminous $L_{\rm FIR} > 10^{13}\,L_\odot$ galaxies. We also present a reanalysis of the spectral energy distributions including the continuum flux densities measured at 3 and 2 mm to derive the overall properties of the sources. Future prospects based on these efficient measurements of redshifts of high-z galaxies using NOEMA are outlined, including a comprehensive survey of all the brightest Herschel galaxies.

The Hubble Frontier Fields offer an exceptionally deep window into the high-redshift universe, covering a substantially larger area than the Hubble Ultra-Deep field at low magnification and probing 1--2 mags deeper in exceptional high-magnification regions. We aim to leverage Atacama Large Millimetre Array (ALMA) band 6 ($\approx$263\,GHz) mosaics in the central portions of five Frontier Fields to characterize the infrared (IR) properties of $1582$ ultraviolet (UV)-selected Lyman-Break Galaxies (LBGs) at redshifts of $z {\sim}$2--8. We investigated individual and stacked fluxes and IR excess (IRX) values of the LBG sample as functions of stellar mass ($\mathrm{M}_{\bigstar}$), redshift, UV luminosity and slope $\beta$, and lensing magnification. Two (2) LBG candidates were individually detected above a significance of $4.1{-}\sigma$, while stacked samples of the remaining LBG candidates yielded no significant detections. We investigated our detections and upper limits in the context of the IRX-$\mathrm{M}_{\bigstar}$ and IRX-$\beta$ relations, probing at least one dex lower in stellar mass than past studies have done. Our upper limits exclude substantial portions of parameter space and they are sufficiently deep in a handful of cases to create mild tension with the typically assumed attenuation and consensus relations. We observe a clear and smooth trend between $\mathrm{M}_{\bigstar}$ and $\beta$, which extends to low masses and blue (low) $\beta$ values, consistent with expectations from previous works.

Dynamical black hole mass measurements in some gas-rich galaxy mergers indicate that they are overmassive relative to their host galaxy properties. Overmassive black holes in these systems present a conflict with the standard progression of galaxy merger - quasar evolution; an alternative explanation is that a nuclear concentration of molecular gas driven inward by the merger is affecting these dynamical black hole mass estimates. We test for the presence of such gas near the two black holes in NGC 6240 using long-baseline ALMA Band 6 observations (beam size 0"06 $\times$ 0"03 or 30 pc$\times$15 pc). We find (4.2-9.8) $\times10^{7}$ M$_{\odot}$ and (1.2-7.7) $\times10^{8}$ M$_{\odot}$ of molecular gas within the resolution limit of the original black hole mass measurements for the north and south black holes, respectively. In the south nucleus, this measurement implies that 6-89% of the original black hole mass measurement actually comes from molecular gas, resolving the tension in the original black hole scaling relations. For the north, only 5% to 11% is coming from molecular gas, suggesting the north black hole is actually overmassive. Our analysis provides the first measurement of significant molecular gas masses contaminating dynamical black hole mass measurements. These high central molecular gas densities further present a challenge to theoretical black hole accretion prescriptions, which often assume accretion proceeds rapidly through the central 10 pc.

In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the first image of the candidate super-massive black hole (SMBH) at the centre of the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image shows a ring of glowing plasma with a dark patch at the centre, which is interpreted as the shadow of the black hole. This breakthrough result, which represents a powerful confirmation of Einstein's theory of gravity, or general relativity, was made possible by assembling a global network of radio telescopes operating at millimetre wavelengths that for the first time included the Atacama Large Millimeter/ submillimeter Array (ALMA). The addition of ALMA as an anchor station has enabled a giant leap forward by increasing the sensitivity limits of the EHT by an order of magnitude, effectively turning it into an imaging array. The published image demonstrates that it is now possible to directly study the event horizon shadows of SMBHs via electromagnetic radiation, thereby transforming this elusive frontier from a mathematical concept into an astrophysical reality. The expansion of the array over the next few years will include new stations on different continents - and eventually satellites in space. This will provide progressively sharper and higher-fidelity images of SMBH candidates, and potentially even movies of the hot plasma orbiting around SMBHs. These improvements will shed light on the processes of black hole accretion and jet formation on event-horizon scales, thereby enabling more precise tests of general relativity in the truly strong field regime.

Dusty star-forming galaxies (DSFGs) detected at $z > 4$ provide important examples of the first generations of massive galaxies. However, few examples with spectroscopic confirmation are currently known, with Hershel struggling to detect significant numbers of $z > 6$ DSFGs. NGP6_D1 is a bright 850 $\mu m$ source (12.3 $\pm$ 2.5 mJy) with no counterparts at shorter wavelengths (a SPIRE dropout). Interferometric observations confirm it is a single source, with no evidence for any optical or NIR emission, or nearby likely foreground lensing sources. No $>3\sigma$ detected lines are seen in both LMT RSR and IRAM 30m EMIR spectra of NGP6_D1 across 32 $GHz$ of bandwidth despite reaching detection limits of $\sim 1 mJy/500 km~s^{-1}$, so the redshift remains unknown. Template fitting suggests that NGP6_D1 is most likely between $z = 5.8$ and 8.3. SED analysis finds that NGP6_D1 is a ULIRG, with a dust mass $\sim 10^8$ - $10^9$ $M_{\odot}$ and a SFR of $\sim$ 500 $M_{\odot}~yr^{-1}$. We place upper limits on the gas mass of NGP6_D1 of $M_{H2}$ $ < (1.1~\pm~3.5) \times 10^{11}$ $M_{\odot}$, consistent with a gas-to-dust ratio of $\sim$ 100 - 1000. We discuss the nature of NGP6_D1 in the context of the broader submm population, and find that comparable SPIRE dropouts account for $\sim$ 20% of all SCUBA-2 detected sources, but with a similar flux density distribution to the general population.

Follow-up observations of (sub-)mm-selected gravitationally-lensed systems have allowed a more detailed study of the dust-enshrouded phase of star-formation up to very early cosmic times. Here, the case of the gravitationally lensed merger in HATLAS J142935.3-002836 (also known as H1429-0028; z_lens=0.218, z_bkg=1.027) is revisited following recent developments in the literature and new APEX observations targeting two carbon monoxide (CO) rotational transitions J_up=3 and 6. We show that the line-profiles comprise three distinct velocity components, where the fainter high-velocity one is less magnified and more compact. The modelling of the observed spectral line energy distribution of CO J_up=2 to 6 and [CI]3P_1-3P_0 assumes a large velocity gradient scenario, where the analysis is based on four statistical approaches. Since the detected gas and dust emission comes exclusively from only one of the two merging components (the one oriented North-South, NS), we are only able to determine upper-limits for the companion. The molecular gas in the NS component in H1429-0028 is found to have a temperature of ~70K, a volume density of log(n[/cm3])~3.7, to be expanding at ~10km/s/pc, and amounts to M_H2=4(-2,+3)*1e9 Msun. The CO to H2 conversion factor is estimated to be alpha_CO=0.4(-0.2,+0.3) Msun/(K.km/s.pc2). The NS galaxy is expected to have a factor of >10x more gas than its companion (M_H2<3e8 Msun). Nevertheless, the total amount of molecular gas in the system comprises only up to 15 per cent (1sigma upper-limit) of the total (dynamical) mass.

We present a detailed description of the special procedures for calibration and quality assurance of Atacama Large Millimeter/submillimeter Array (ALMA) observations in Very Long Baseline Interferometry (VLBI) mode. These procedures are required to turn the phased ALMA array into a fully calibrated VLBI station. As an illustration of these methodologies, we present full-polarization observations carried out with ALMA as a phased array at 3mm (Band 3) and 1.3mm (Band 6) as part of Cycle-4. These are the first VLBI science observations conducted with ALMA and were obtained during a 2017 VLBI campaign in concert with other telescopes worldwide as part of the Global mm-VLBI Array (GMVA, April 1-3) and the Event Horizon Telescope (EHT, April 5-11) in ALMA Bands 3 and 6, respectively.

We present an exploration of the expected detection of the earliest Active Galactic Nuclei (AGN) in the Universe from state-of-art galaxy formation and evolution semi-analytic models and hydro-dynamical simulations. We estimate the number and radiative characteristics of Super Massive Black Holes (SMBHs) at $z\geq 6$, a redshift range that will be intensively explored by the next generation of telescopes, in particular in the radio through the Square Kilometre Array (SKA) and at high energies with ESA's Athena X-ray Observatory. We find that Athena will be able to observe over 5000 AGN/$\rm deg^2$ at the Epoch of Re-ionization (EoR), $6\leq z \leq 10$. Similarly, for the same redshift range the models/simulations suggest that SKA will detect at least 400 AGN/$\rm deg^2$. Additionally, we stress the importance of the volume of the simulation box as well as the initial physical conditions of the models/simulations on their effect on the luminosity functions (LFs) and the creation of the most massive SMBHs that we currently observe at the EoR. Furthermore, following the evolution of the accretion mode of the SMBHs in each model/simulation, we show that, while the quasar dominates over the radio mode at the EoR, detection at radio wavelengths still reaches significant numbers even at the highest redshifts. Finally, we present the effect that the radiative efficiency has on the LFs by comparing results produced with a constant value for the radiative efficiency and more complex calculations based on the spin of each SMBH.

We present Atacama Large Millimeter/submillimeter Array (ALMA) resolved observations of molecular gas in galaxies up to $z=0.35$ to characterise the role of global galactic dynamics on the global interstellar medium (ISM) properties. These observations consist of a sub-sample of 39 galaxies taken from the Valpara\'iso ALMA Line Emission Survey (VALES). From the CO($J=1-0)$ emission line, we quantify the kinematic parameters by modelling the velocity fields. We find that the IR luminosity increases with the rotational to dispersion velocity ratio ($V_{\rm rot}/\sigma_v$, corrected for inclination). We find a dependence between $V_{\rm rot}/\sigma_v$ and the [CII]/IR ratio, suggesting that the so-called `[CII] deficit' is related to the dynamical state of the galaxies. We find that global pressure support is needed to reconcile the dynamical mass estimates with the stellar masses in our systems with low $V_{\rm rot}/\sigma_v$ values. The star formation rate (SFR) is weakly correlated with the molecular gas fraction ($f_{\rm H_2}$) in our sample, suggesting that the release of gravitational energy from cold gas may not be the main energy source of the turbulent motions seen in the VALES galaxies. By defining a proxy of the `star formation efficiency' parameter as the SFR divided by the CO luminosity (SFE$'\equiv$ SFR/L$'_{\rm CO}$), we find a constant SFE$'$ per crossing time ($t_{\rm cross}$). We suggest that $t_{\rm cross}$ may be the controlling timescale in which the star formation occurs in dusty $z\sim0.03-0.35$ galaxies.

We present a 69 arcmin$^2$ ALMA survey at 1.1mm, GOODS-ALMA, matching the deepest HST-WFC3 H-band part of the GOODS-South field. We taper the 0"24 original image with a homogeneous and circular synthesized beam of 0"60 to reduce the number of independent beams - thus reducing the number of purely statistical spurious detections - and optimize the sensitivity to point sources. We extract a catalogue of galaxies purely selected by ALMA and identify sources with and without HST counterparts down to a 5$\sigma$ limiting depth of H=28.2 AB (HST/WFC3 F160W). ALMA detects 20 sources brighter than 0.7 mJy in the 0"60 tapered mosaic (rms sensitivity =0.18 mJy/beam) with a purity greater than 80%. Among these detections, we identify three sources with no HST nor Spitzer-IRAC counterpart, consistent with the expected number of spurious galaxies from the analysis of the inverted image; their definitive status will require additional investigation. An additional three sources with HST counterparts are detected either at high significance in the higher resolution map, or with different detection-algorithm parameters ensuring a purity greater than 80%. Hence we identify in total 20 robust detections. Our wide contiguous survey allows us to push further in redshift the blind detection of massive galaxies with ALMA with a median redshift of $z$=2.92 and a median stellar mass of M$_{\star}$ = 1.1 $\times 10^{11}$M$_\odot$. Our sample includes 20% HST-dark galaxies (4 out of 20), all detected in the mid-infrared with IRAC. The near-infrared based photometric redshifts of two of them $z\sim$4.3 and 4.8) suggest that these sources have redshifts $z$>4. At least 40% of the ALMA sources host an X-ray AGN, compared to 14% for other galaxies of similar mass and redshift. The wide area of our ALMA survey provides lower values at the bright end of number counts than single-dish telescopes

We present the detection at 89 $\mu$m (observed frame) of the {\it Herschel}-selected gravitationally lensed starburst galaxy HATLASJ1429-0028 (also known as G15v2.19) in 15 minutes with the High-resolution Airborne Wideband Camera-plus (HAWC+) onboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The spectacular lensing system consists of an edge-on foreground disk galaxy at $z$ = 0.22 and a nearly complete Einstein ring of an intrinsic ultra-luminous infrared galaxy at $z$ = 1.03. Is this high IR luminosity powered by pure star formation (SF) or also an active galactic nucleus (AGN)? Previous nebular line diagnostics indicate that it is star-formation dominated. We perform a 27-band multi-wavelength spectral energy distribution modeling (SED) including the new SOFIA/HAWC+ data to constrain the fractional AGN contribution to the total IR luminosity. The AGN fraction in the IR turns out to be negligible. In addition, J1429-0028 serves as a testbed for comparing SED results from different models/templates and SED codes (MAGPHYS, SED3FIT, and CIGALE). We stress that star formation history is the dominant source of uncertainty in the derived stellar mass (as high as a factor of $\sim$ 10) even in the case of extensive photometric coverage. Furthermore, the detection of a source at $z$ $\sim$ 1 with SOFIA/HAWC+ demonstrates the potential of utilizing this facility for distant galaxy studies including the decomposition of SF/AGN components, which cannot be accomplished with other current facilities.

Hyperluminous infrared galaxies (HyLIRGs) provide crucial "long lever arm" constraints on galaxy evolution. H-ATLAS $J084933.4+021443$, a $z=2.41$ binary HyLIRG with at least two additional luminous companion galaxies, is thus an optimal test-ground for studies of star formation and galaxy evolution during "cosmic noon". We have used ALMA to obtain resolved imaging and kinematics of atomic and molecular emission lines, and rest-frame $340$ to $1160$GHz continuum emission, for the known luminous component galaxies in H-ATLAS $J084933.4+021443$: W, T, M, C. All four component galaxies are spatially ($\sim 0 .\!\!^{''} 3$ or $2.5$kpc) resolved in CO J:7-6, [C$_\mathrm{I}$] 2-1, H$_2$O and the millimetre (mm) to sub-mm continuum. Rotation-dominated gas kinematics is confirmed in W and T. The significant extension to component T, in gas and continuum, along its kinematic minor axis, is attributable to its lensing magnification. Spatially resolved sub-mm spectral energy distributions reveal that component W is well fit with greybody emission from dust at a single temperature over the full extent of the galaxy, despite it containing a powerful AGN, while component T requires an additional component of hotter nuclear dust and additional sources of emission in the mm. We confirm that [C$_\mathrm{I}$] 2-1 can be used as a rough tracer of warm/dense molecular gas in extreme systems, though the [C$_\mathrm{I}$] 2-1/CO luminosity ratio increases sub-linearly. We obtain an exquisite and unprecedented resolved ($2.5$-kpc-scale) "warm/dense molecular gas" Schmidt-Kennicutt (SK) relationship for components W and T. Gas exhaustion times for all apertures in W (T) are $1-4$Gyr ($0.5-2$Gyr). Both W and T follow a resolved "warm/dense gas" SK relationship with power law $n\sim1.7$, significantly steeper than the $n\sim1$ found previously via "cold" molecular gas in nearby "normal" star-forming galaxies.

We present optical and near-IR Integral Field Unit (IFU) and ALMA band 6 observations of the nearby dual Active Galactic Nuclei (AGN) Mrk 463. At a distance of 210 Mpc, and a nuclear separation of $\sim$4 kpc, Mrk 463 is an excellent laboratory to study the gas dynamics, star formation processes and supermassive black hole (SMBH) accretion in a late-stage gas-rich major galaxy merger. The IFU observations reveal a complex morphology, including tidal tails, star-forming clumps, and emission line regions. The optical data, which map the full extent of the merger, show evidence for a biconical outflow and material outflowing at $>$600 km s$^{-1}$, both associated with the Mrk 463E nucleus, together with large scale gradients likely related to the ongoing galaxy merger. We further find an emission line region $\sim$11 kpc south of Mrk 463E that is consistent with being photoionized by an AGN. Compared to the current AGN luminosity, the energy budget of the cloud implies a luminosity drop in Mrk 463E by a factor 3-20 over the last 40,000 years. The ALMA observations of $^{12}$CO(2-1) and adjacent 1mm continuum reveal the presence of $\sim$10$^{9}$M$_\odot$ in molecular gas in the system. The molecular gas shows velocity gradients of $\sim$800 km/s and $\sim$400 km/s around the Mrk 463E and 463W nuclei, respectively. We conclude that in this system the infall of $\sim$100s $M_\odot$/yr of molecular gas is in rough balance with the removal of ionized gas by a biconical outflow being fueled by a relatively small, $<$0.01% of accretion onto each SMBH.

The Atacama Millimeter/submillimeter Array (ALMA) Phasing Project (APP) has developed and deployed the hardware and software necessary to coherently sum the signals of individual ALMA antennas and record the aggregate sum in Very Long Baseline Interferometry (VLBI) Data Exchange Format. These beamforming capabilities allow the ALMA array to collectively function as the equivalent of a single large aperture and participate in global VLBI arrays. The inclusion of phased ALMA in current VLBI networks operating at (sub)millimeter wavelengths provides an order of magnitude improvement in sensitivity, as well as enhancements in u-v coverage and north-south angular resolution. The availability of a phased ALMA enables a wide range of new ultra-high angular resolution science applications, including the resolution of supermassive black holes on event horizon scales and studies of the launch and collimation of astrophysical jets. It also provides a high-sensitivity aperture that may be used for investigations such as pulsar searches at high frequencies. This paper provides an overview of the ALMA Phasing System design, implementation, and performance characteristics.

[abridged] The Hubble and Spitzer Space Telescope surveys of the Frontier Fields (FF) provide extremely deep images around six massive, strong-lensing clusters of galaxies. The ALMA FF survey aims to cover the same fields at 1.1mm, with maps reaching (unlensed) sensitivities of $<$70$\mu$Jy, in order to explore the properties of background dusty star-forming galaxies. We report on the multi-wavelength photometric analysis of all 12 significantly detected ($>$5$\sigma$) sources in the first three FF clusters observed by ALMA, based on data from Hubble and Spitzer, the VLT and Herschel. We measure the total photometry in all available bands and determine the photometric redshifts and the physical properties of the counterparts via SED-fitting. In particular, we carefully estimate the FIR photometry using 1.1mm priors to limit the misidentification of blended FIR counterparts, which strongly affect some flux estimates in previous FIR catalogs. We identify robust near-infrared (NIR) counterparts for all 11 sources with K$_s$ detection, the majority of which are quite red, with eight having $F814W-K_s\gtrsim 4$ and five having $F160W-[4.5]\gtrsim3$. From the FIR point of view, all our objects have $z_{phot}$$\sim$1--3, whereas based on the optical SED one object prefers a high-$z$ solution ($z\geq\ $7). Five objects among our sample have spectroscopic redshifts from the GLASS survey for which we can reproduce their SEDs with existing templates. This verification confirms the validity of our photometric redshift methodology. The mean redshift of our sample is $z_{phot}$=1.99$\pm$0.27. All 1.1mm selected objects are massive (10.0$<\log[M_{\star}(M_{\odot})]<$ 11.5), with high star formation rates ($<\log[SFR(M_{\odot}/yr)]> \approx$1.6) and high dust contents (8.1 $<\log[M_{dust} (M_{\odot})]<$8.8), consistent with previous ALMA surveys.

We aim to characterize the relationship between dust properties. We also aim to provide equations to estimate accurate dust properties from limited observational datasets. We assemble a sample of 1,630 nearby (z<0.1) galaxies-over a large range of Mstar, SFR - with multi-wavelength observations available from wise, iras, planck and/or SCUBA. The characterization of dust emission comes from SED fitting using Draine & Li dust models, which we parametrize using two components (warm and cold ). The subsample of these galaxies with global measurements of CO and/or HI are used to explore the molecular and/or atomic gas content of the galaxies. The total Lir, Mdust and dust temperature of the cold component (Tc) form a plane that we refer to as the dust plane. A galaxy's sSFR drives its position on the dust plane: starburst galaxies show higher Lir, Mdust and Tc compared to Main Sequence and passive galaxies. Starburst galaxies also show higher specific Mdust (Mdust/Mstar) and specific Mgas (Mgas/Mstar). The Mdust is more closely correlated with the total Mgas (atomic plus molecular) than with the individual components. Our multi wavelength data allows us to define several equations to estimate Lir, Mdust and Tc from one or two monochromatic luminosities in the infrared and/or sub-millimeter. We estimate the dust mass and infrared luminosity from a single monochromatic luminosity within the R-J tail of the dust emission, with errors of 0.12 and 0.20dex, respectively. These errors are reduced to 0.05 and 0.10 dex, respectively, if the Tc is used. The Mdust is correlated with the total Mism (Mism \propto Mdust^0.7). For galaxies with Mstar 8.5<log(Mstar/Msun) < 11.9, the conversion factor \alpha_850mum shows a large scatter (rms=0.29dex). The SF mode of a galaxy shows a correlation with both the Mgass and Mdust: high Mdust/Mstar galaxies are gas-rich and show the highest SFRs.

Dusty star-forming galaxies are among the most prodigious systems at high redshift (z>1), characterized by high star formation rates and huge dust reservoirs. The bright end of this population has been well characterized in recent years, but considerable uncertainties remain for fainter dusty star-forming galaxies, which are responsible for the bulk of star formation at high redshift and thus play a key role in galaxy growth and evolution. In this first paper of our series, we describe our methods for finding high redshift faint dusty galaxies using millimeter observations with ALMA. We obtained ALMA 1.1 mm mosaic images for three strong-lensing galaxy clusters from the Frontier Fields survey. The 2'x2' mosaics overlap with the deep HST WFC3/IR footprints and encompass the high magnification regions of each cluster. The combination of extremely high ALMA sensitivity and the magnification power of these clusters allows us to systematically probe the sub-mJy population of dusty star-forming galaxies over a large surveyed area. We present a description of the reduction and analysis of the ALMA continuum observations for the galaxy clusters Abell 2744 (z=0.308), MACSJ0416.1-2403 (z=0.396) and MACSJ1149.5+2223 (z=0.543), for which we reach observed rms sensitivities of 55, 59 and 71 $\mu$Jy beam$^{-1}$ respectively. We detect 12 dusty star-forming galaxies at S/N>5.0 across the three clusters, all of them presenting coincidence with NIR detected counterparts in the HST images. None of the sources fall close to the lensing caustics, thus they are not strongly lensed. The observed 1.1 mm flux densities for the total sample of galaxies range from 0.41 to 2.82 mJy, with observed effective radii spanning <0.05" to 0.37"$\pm$0.21". The lensing-corrected sizes of the detected sources appear to be in the same range as those measured in brighter samples, albeit with possibly larger dispersion.

We present a source-plane reconstruction of a ${\it Herschel}$ and ${\it Planck}$-detected gravitationally-lensed dusty star-forming galaxy (DSFG) at $z=1.68$ using {\it Hubble}, Sub-millimeter Array (SMA), and Keck observations. The background sub-millimeter galaxy (SMG) is strongly lensed by a foreground galaxy cluster at $z=0.997$ and appears as an arc of length $\sim 15^{\prime \prime}$ in the optical images. The continuum dust emission, as seen by SMA, is limited to a single knot within this arc. We present a lens model with source plane reconstructions at several wavelengths to show the difference in magnification between the stars and dust, and highlight the importance of a multi-wavelength lens models for studies involving lensed DSFGs. We estimate the physical properties of the galaxy by fitting the flux densities to model SEDs leading to a magnification-corrected star formation rate of $390 \pm 60$ M$_{\odot}$ yr$^{-1}$ and a stellar mass of $1.1 \pm 0.4\times 10^{11}$ M$_{\odot}$. These values are consistent with high-redshift massive galaxies that have formed most of their stars already. The estimated gas-to-baryon fraction, molecular gas surface density, and SFR surface density have values of $0.43 \pm 0.13$, $350 \pm 200$ M$_{\odot}$ pc$^{-2}$, and $\sim 12 \pm 7~$M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$, respectively. The ratio of star formation rate surface density to molecular gas surface density puts this among the most star-forming systems, similar to other measured SMGs and local ULIRGS.

We have recently suggested that gas accretion can be studied using host galaxies of gamma-ray bursts (GRBs). We obtained the first ever far-infrared (FIR) line observations of a GRB host, namely Herschel/PACS resolved [CII] 158 um and [OI] 63 um spectroscopy, as well as APEX CO(2-1) and ALMA CO(1-0) observations of the GRB 980425 host. It has elevated [CII]/FIR and [OI]/FIR ratios and higher values of star formation rate (SFR) derived from line ([CII], [OI], Ha) than from continuum (UV, IR, radio) indicators. [CII] emission exhibits a normal morphology, peaking at the galaxy center, whereas [OI] is concentrated close to the GRB position and the nearby Wolf-Rayet region. The high [OI] flux indicates high radiation field and gas density. The [CII]/CO luminosity ratio of the GRB 980425 host is close to the highest values found for local star-forming galaxies. Its CO-derived molecular gas mass is low given its SFR and metallicity, but the [CII]-derived molecular gas mass is close to the expected value. The [OI] and HI concentrations, and the high radiation field and density are consistent with the hypothesis of a very recent (at most a few tens of Myr ago) inflow of atomic gas triggering star formation. Dust has not had time to build up (explaining high line-to-continuum ratios). Such a recent enhancement of star-formation would indeed manifest itself in high SFR_line/SFR_continuum ratios, because the line indicators are sensitive only to recent (<10 Myr) activity, whereas the continuum indicators measure the SFR averaged over much longer periods (~100 Myr). Other GRB hosts exhibit a mean SFR_line/SFR_continuum of 1.74+-0.32. This is consistent with a very recent enhancement of star formation being common among GRB hosts, so galaxies which have recently experienced inflow of gas may preferentially host stars exploding as GRBs. Hence GRB hosts may be used to investigate recent gas accretion.

Exploiting the sensitivity and spatial resolution of the Atacama Large Millimeter/submillimeter Array (ALMA), we have studied the morphology and the physical scale of the interstellar medium - both gas and dust - in SGP38326, an unlensed pair of interacting starbursts at $z= 4.425$. SGP38326 is the most luminous star bursting system known at $z > 4$ with an IR-derived ${\rm SFR \sim 4300 \,} M_\odot \, {\rm yr}^{-1}$. SGP38326 also contains a molecular gas reservoir among the most massive ever found in the early Universe, and it is the likely progenitor of a massive, red-and-dead elliptical galaxy at $z \sim 3$. Probing scales of $\sim 0.1"$ or $\sim 800 \, {\rm pc}$ we find that the smooth distribution of the continuum emission from cool dust grains contrasts with the more irregular morphology of the gas, as traced by the [CII] fine structure emission. The gas is also extended over larger physical scales than the dust. The velocity information provided by the resolved [CII] emission reveals that the dynamics of the two components of SGP38326 are compatible with disk-like, ordered rotation, but also reveals an ISM which is turbulent and unstable. Our observations support a scenario where at least a subset of the most distant extreme starbursts are highly dissipative mergers of gas-rich galaxies.

In recent years there has been a paradigm shift from centralised to geographically distributed resources. Individual entities are no longer able to host or afford the necessary expertise in-house, and, as a consequence, society increasingly relies on widespread collaborations. Although such collaborations are now the norm for scientific projects, more technical structures providing support to a distributed scientific community without direct financial or other material benefits are scarce. The network of European ALMA Regional Centre (ARC) nodes is an example of such an internationally distributed user support network. It is an organised effort to provide the European ALMA user community with uniform expert support to enable optimal usage and scientific output of the ALMA facility. The network model for the European ARC nodes is described in terms of its organisation, communication strategies and user support.

We present Hubble Space Telescope (HST) WFC3 imaging and grism spectroscopy observations of the {\it Herschel}-selected gravitationally-lensed starburst galaxy HATLASJ1429-0028. The lensing system consists of an edge-on foreground disk galaxy at $z=0.218$ with a nearly complete Einstein ring of the infrared luminous galaxy at $z=1.027$. The WFC3 spectroscopy with G102 and G141 grisms, covering the wavelength range of 0.8 to 1.7 $\mu$m, resulted in detections of H$\alpha$+[NII], H$\beta$, [SII], and [OIII] for the background galaxy from which we measure line fluxes and ratios. The Balmer line ratio H$\alpha$/H$\beta$ of $7.5 \pm 4.4$, when corrected for [NII], results in an extinction for the starburst galaxy of $E(B-V)=0.8 \pm 0.5$. The $H\alpha$ based star-formation rate, when corrected for extinction, is $60 \pm 50$ M$_{\odot}$ yr$^{-1}$, lower than the instantaneous star-formation rate of 390 $\pm$ 90 M$_{\odot}$ yr$^{-1}$ from the total IR luminosity. We also compare the nebular line ratios of HATLASJ1429-0028 with other star-forming and sub-mm bright galaxies. The nebular line ratios are consistent with an intrinsic ultra-luminous infrared galaxy with no evidence for excitation by an active galactic nucleus (AGN). We estimate the metallicity, $12 + log(O/H)$, of HATLASJ1429-0028 to be 8.49 $\pm$ 0.16. Such a low value is below the average relations for stellar mass vs. metallicity of galaxies at $z \sim 1$ for a galaxy with stellar mass of $\sim 2 \times 10^{11}$ M$_{\odot}$. The combination of high stellar mass, lack of AGN indicators, low metallicity, and the high star-formation rate of HATLASJ1429-0028 suggest that this galaxy is currently undergoing a rapid formation.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

[Abridged] Aims: This work focuses on one lensed system, HATLAS J142935.3-002836 (H1429-0028), selected in the Herschel-ATLAS field. Gathering a rich, multi-wavelength dataset, we aim to confirm the lensing hypothesis and model the background source's morphology and dynamics, as well as to provide a full physical characterisation. Methods: Multi-wavelength high-resolution data is utilised to assess the nature of the system. A lensing-analysis algorithm which simultaneously fits different wavebands is adopted to characterise the lens. The background galaxy dynamical information is studied by reconstructing the 3-D source-plane of the ALMA CO(J:4-3) transition. Near-IR imaging from HST and Keck-AO allows to constrain rest-frame optical photometry independently for the foreground and background systems. Physical parameters (such as stellar and dust masses) are estimated via modelling of the spectral energy distribution taking into account source blending, foreground obscuration, and differential magnification. Results: The system comprises a foreground edge-on disk galaxy (at z_sp=0.218) with an almost complete Einstein ring around it. The background source (at z_sp=1.027) is magnified by a factor of ~8-10 depending on wavelength. It is comprised of two components and a tens of kpc long tidal tail resembling the Antennae merger. As a whole, the system is a massive stellar system (1.32[-0.41,+0.63] x1E11 Mo) forming stars at a rate of 394+-90 Mo/yr, and has a significant gas reservoir M_ISM = 4.6+-1.7 x1E10 Mo. Its depletion time due to star formation alone is thus expected to be tau_SF=M_ISM/SFR=117+-51 Myr. The dynamical mass of one of the components is estimated to be 5.8+-1.7 x1E10 Mo, and, together with the photometric total mass estimate, it implies that H1429-0028 is a major merger system (1:2.8[-1.5,+1.8]).

I make a short revision of Cosmology questions which ALMA was built to address. Without diving into much detail, I point out the ALMA specifications and strategies which are expected to provide a better handle of: the temperature evolution of the Cosmic Microwave Background (CMB) and the properties of its secondary anisotropies (such as the thermal and kinetic Sunyaev-Zel'dovich and the Ostriker-Vishniac effects); variability of dimensionless fundamental constants; Ho and galaxy initial mass function by means of strong gravitational lensing; black hole science with the greatly expected Event Horizon Telescope.

Chandra data in the COSMOS, AEGIS-XD and 4Ms CDFS are combined with optical/near-IR photometry to determine the rest-frame U-V vs V-J colours of X-ray AGN hosts at mean redshifts 0.40 and 0.85. This combination of colours (UVJ) provides an efficient means of separating quiescent from star-forming, including dust reddened, galaxies. Morphological information emphasises differences between AGN split by their UVJ colours. AGN in quiescent galaxies are dominated by spheroids, while star-forming hosts are split between bulges and disks. The UVJ diagram of AGN hosts is then used to set limits on the accretion density associated with evolved and star-forming systems. Most of the black hole growth since z~1 is associated with star-forming hosts. Nevertheless, ~15-20% of the X-ray luminosity density since z~1, is taking place in the quiescent region of the UVJ diagram. For the z~0.40 subsample, there is tentative evidence (2sigma significance), that AGN split by their UVJ colours differ in Eddington ratio. AGN in star-forming hosts dominate at high Eddington ratios, while AGN in quiescent hosts become increasingly important as a fraction of the total population toward low Eddington ratios. At higher redshift, z~0.8, such differences are significant at the 2sigma level only at Eddington ratios >1e-3. These findings are consistent with scenarios in which diverse accretion modes are responsible for the build-up of SMBHs at the centres of galaxies. We compare our results with the GALFORM semi-analytic model, which postulates two black hole fuelling modes, the first linked to star-formation and the second occuring in passive galaxies. GALFORM predicts a larger fraction of black hole growth in quiescent galaxies at z<1, compared to the data. Relaxing the strong assumption of the model that passive AGN hosts have zero star-formation rate could reconcile this disagreement.

We present a catalogue of bulgeless galaxies, which includes 19225 objects selected in four of the deepest, largest multi-wavelength datasets available -- COSMOS, AEGIS, GEMS and GOODS -- at intermediate redshift ($0.4 \leq z \leq 1.0$). The morphological classification was provided by the Advanced Camera for Surveys General Catalogue (ACS-GC), which used publicly available data obtained with the ACS instrument on the Hubble Space Telescope. Rest-frame photometric quantities were derived using kcorrect. We analyse the properties of the sample and the evolution of pure-disc systems with redshift. Very massive [$\log (M_\star/M_{\odot}) > 10.5$] bulgeless galaxies contribute to ~30% of the total galaxy population number density at $z \geq 0.7$, but their number density drops substantially with decreasing redshift. We show that only a negligible fraction of pure discs appear to be quiescent systems, and red sequence bulgeless galaxies show indications of dust-obscured star formation. X-ray catalogues were used to search for X-ray emission within our sample. After visual inspection and detailed parametric morphological fitting we identify 30 AGN that reside in galaxies without a classical bulge. The finding of such peculiar objects at intermediate redshift shows that while AGN growth in merger-free systems is a rare event (0.2% AGN hosts in this sample of bulgeless galaxies), it can indeed happen relatively early in the Universe history.

AIMS. This work addresses the AGN IR-selection dependency on intrinsic source luminosity and obscuration, in order to identify and charaterise biases which could affect conclusions in studies. METHODS. We study IR-selected AGN in the Chandra Deep Field South (CDFS) survey and in the Cosmological Survey (COSMOS). The AGN sample is divided into low and high X-ray luminosity classes, and into unobscured (type-1) and obscured (type-2) classes by means of X-ray and optical spectroscopy data. Specifically in the X-rays regime, we adopt the intrinsic luminosity taking into account the estimated column density (N_H). We also take the chance to highlight important differences resulting from adopting different methods to assess AGN obscuration. RESULTS. In agreement with previous studies, we also find that AGN IR-selection efficiency shows a decrease with decreasing source AGN X-ray luminosity. For the intermediate-luminosity AGN population (43.3<~log(L_X[erg/s])<~44$), the efficiency also worsens with increasing obscuration. The same sample also shows an evolution with cosmic time of the obscured fraction at the highest X-ray luminosities, independently of adopted type-1/type-2 classification method. CONCLUSIONS. We confirm that AGN IR-selection is genuinely biased toward unobscured AGN, but only at intermediate luminosities. At the highest luminosities, where AGN IR-selection is more efficient, there is no obscuration bias. We show that type-1 AGN are intrinsically more luminous than type-2 AGN only at z<~1.6, thus resulting in more type-1 AGN being selected the shallower the IR survey is. [ABRIDGED: Conclusion on survey strategy. See manuscript.] Finally, when the James Webb Space Telescope comes online, the broad-band filters 2.0um, 4.4um, 7.7um, and 18um will be key to disentangle AGN from non-AGN dominated SEDs at depths where spectroscopy becomes impractical.

[Abridged] We study the contribution of hot-dust to the luminosity density of galaxies and its evolution with cosmic time. Using the Spitzer-IRAC data in the COSMOS field, we estimate the contribution from hot-dust at rest-frame 4.2um (from ~0 < z < ~0.2 up to ~0.5 < z < ~0.9). This wavelength corresponds to black-body temperature of ~690K. The contribution due to stellar emission is estimated from the rest-frame 1.6um luminosity (assumed to result from stellar emission alone) and subtracted from the mid-infrared luminosity of galaxies to measure hot-dust emission. In order to attempt the study of the 3.3um-PAH feature, we use the rest-frame 4.2um to infer the hot-dust flux at 3.3um. This study is performed for different spectral types of galaxies: early-type, late-type, starburst, and IR-selected AGN. We find that: (a) the decrease of the hot-dust luminosity density since ~0.5 < z < ~1 is steeper (by at least ~0.5dex) compared to that of the cold-dust, giving support to the scenario where galaxy obscuration increases with redshift, as already proposed in the literature; (b) hot-dust and PAH emission evolution seems to be correlated with stellar mass, where rest-frame 1.6um luminous non-AGN galaxies (i.e., massive systems) show a stronger decrement (with decreasing redshift) in hot-dust and PAH emission than the less luminous (less massive) non-AGN galaxies; (c) despite comprising < ~3% of the total sample, AGN contribute as much as a third to the hot-dust luminosity density at z < 1 and clearly dominate the bright-end of the total hot-dust Luminosity Density Function at ~0.5 < z < ~0.9; (d) the average dust-to-total luminosity ratio increases with redshift, while PAH-to-total luminosity ratio remains fairly constant; (e) at M_1.6 > -25, the dust-to-total and PAH-to-total luminosity ratios increase with decreasing luminosity, but deeper data is required to confirm this result.

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.

Using a sample of BzK-selected galaxies at z~2 identified from the CFHT/WIRCAM near-infrared survey of GOODS-North, we discuss the relation between star formation rate (SFR), specific star formation rate (SSFR), and stellar mass (M_{*}), and the clustering of galaxies as a function of these parameters. For star-forming galaxies (sBzKs), the UV-based SFR, corrected for extinction, scales with the stellar mass as SFR ~ M_{*}^{alpha} with alpha = 0.74+/-0.20 down to M_{*} ~ 10^{9} M_{solar}, indicating a weak dependence on the stellar mass of the SSFR. We also measure the angular correlation function and hence infer the correlation length for sBzK galaxies as a function of M_{*}, SFR, and SSFR, as well as K-band apparent magnitude. We show that passive galaxies (pBzKs) are more strongly clustered than sBzK galaxies at a given stellar mass, mirroring the color-density relation seen at lower redshifts. We also find that the correlation length of sBzK galaxies ranges from 4 to 20 h^{-1}Mpc, being a strong function of M_{K}, M_{*}, and SFR. On the other hand, the clustering dependence on SSFR changes abruptly at 2x10^{-9} yr^{-1}, which is the typical value for "main sequence" star-forming galaxies at z~2. We show that the correlation length reaches a minimum at this characteristic value, and is larger for galaxies with both smaller and larger SSFRs; a dichotomy that is only marginally implied from the predictions of the semi-analytical models. Our results suggest that there are two types of environmental effects at work at z~2. Stronger clustering for relatively quiescent galaxies implies that the environment has started to play a role in quenching star formation. At the same time, stronger clustering for galaxies with elevated SSFRs ("starbursts") might be attributed to an increased efficiency for galaxy interactions and mergers in dense environments.

[Abridged] It is widely accepted that observations at mid-infrared (mid-IR) wavelengths enable the selection of galaxies with nuclear activity, which may not be revealed even in the deepest X-ray surveys. In this work new near- and mid-IR color diagnostics are explored, aiming for improved efficiency - better completeness and less contamination - in selecting AGN out to very high redshifts. We restrict our study to the James Webb Space Telescope wavelength range (0.6-27um). The criteria are created based on the predictions by state-of-the-art galaxy and AGN templates covering a wide variety of galaxy properties, and tested against control samples with deep multi-wavelength coverage (ranging from the X-rays to radio frequencies). We show that the colors Ks-[4.5], [4.5]-[8.0], and [8.0]-[24] are ideal as AGN/non-AGN diagnostics at, respectively, z<~1, 1<~z<~2.5, and z>~2.5-3. However, when the source redshift is unknown, these colors should be combined. We thus develop an improved IR criterion (using Ks and IRAC bands, KI) as a new alternative at z<~2.5. KI does not show improved completeness (50-60% overall) in comparison to commonly used IRAC-based AGN criteria, but is less affected by non-AGN contamination (revealing a >50-90% level of successful AGN selection). We also propose KIM (using Ks, IRAC, and MIPS-24um bands, KIM), which aims to select AGN hosts from local distances to as far back as the end of reionization (0<z<~7) with reduced non-AGN contamination. However, the necessary testing-constraints and the small control-sample sizes prevent the confirmation of its improved efficiency at z<~2.5. Overall, KIM shows a ~30-40% completeness and a >70-90% level of successful AGN selection. KI and KIM are built to be reliable against a ~10-20% error in flux, are based on existing filters, and are suitable for immediate use.

We investigate how star-forming galaxies typically assemble their masses at high redshift. Using the deep multi-wavelength coverage of the GOODS dataset, we measure stellar mass of a large sample of star-forming galaxies at z~4 and 5, and make a robust determination of stellar mass function (SMF). We report a broad correlation between stellar mass and UV luminosity, such that more UV-luminous galaxies are, on average, more massive. However, the correlation has a substantial intrinsic scatter evidenced by a non-negligible number of UV-faint but massive galaxies. Furthermore, the low-mass end of the SMF does not rise as steeply as the UV luminosity function (alpha_UVLF} -(1.7-1.8), alpha_SMF -(1.3-1.4)) of the same galaxies. In a smooth formation scenario where star formation (SF) is sustained at the observed rates for a long time, these galaxies would have accumulated more mass (by a factor of ~3) than observed and therefore the SMF would mirror more closely that of the UVLF. The relatively shallow slope of the SMF is due to the fact that many of the UV-selected galaxies are not massive enough, and therefore are too faint in their rest-frame optical bands, to be detected in the current observations. Our results favor a episodic formation history in which SF in low-mass galaxies vary significantly over time, a scenario favored by galaxy clustering. Our findings for the UV-faint galaxies are in contrast with those found for more UV-luminous galaxies, which exhibit tighter SFR-M_star correlations. The discrepancy may suggest that galaxies at different luminosities may have different evolutionary paths.

A new set of color selection criteria (VJL) analogous with the BzK method is designed to select both star-forming galaxies (SFGs) and passively-evolving galaxies (PEGs) at 2.3<z<3.5 by using rest-frame UV--optical (V-J vs. J-L) colors. The criteria are thoroughly tested with theoretical stellar population synthesis models and real galaxies with spectroscopic redshifts to evaluate their efficiency and contamination. We apply the well-tested VJL criteria to the HST/WFC3 Early Release Science field and study the physical properties of selected galaxies. The redshift distribution of selected SFGs peaks at z~2.7, slightly lower than that of Lyman Break Galaxies at z~3. Comparing the observed mid-infrared fluxes of selected galaxies with the prediction of pure stellar emission, we find that our VJL method is effective at selecting massive dusty SFGs that are missed by the Lyman Break Technique. About half of the star formation in massive (M_{star}>10^{10}M_{Sun}) galaxies at 2.3<z<3.5 is contributed by dusty (extinction E(B-V)>0.4) SFGs, which however, only account for ~20% of the number density of massive SFGs. We also use the mid-infrared fluxes to clean our PEG sample, and find that galaxy size can be used as a secondary criterion to effectively eliminate the contamination of dusty SFGs. The redshift distribution of the cleaned PEG sample peaks at z~2.5. We find 6 PEG candidates at z>3 and discuss possible methods to distinguish them from dusty contamination. We conclude that at least part of our candidates are real PEGs at z~3, implying that this type of galaxies began to form their stars at z>5. We measure the integrated stellar mass density of PEGs at z~2.5 and set constraints on it at z>3. We find that the integrated stellar mass density grows by at least about factor of 10 in 1 Gyr at 3<z<5 and by another factor of 10 in next 3.5 Gyr (1<z<3).

The main focus of this thesis is the IR spectral regime, which since the 70's and 80's has revolutionised our understanding of the Universe. A multi-wavelength analysis on Extremely Red Galaxy populations is first presented in one of the most intensively observed patch of the sky, the Chandra Deep Field South. By adopting a purely statistical methodology, we consider all the photometric and spectroscopic information available on large samples of Extremely Red Objects (EROs, 553 sources), IRAC EROs (IEROs, 259 sources), and Distant Red Galaxies (DRGs, 289 sources). We derive general properties: redshift distributions, AGN host fraction, star-formation rate densities, dust content, morphology, mass functions and mass densities. The results point to the fact that EROs, IEROs, and DRGs all belong to the same population, yet seen at different phases of galaxy evolution. The second part of this thesis is dedicated to the AGN selection in the IR, with particular relevance to the James Webb Space Telescope, to be launched in 2018. We develop an improved IR criterion (using K and IRAC bands) as an alternative to existing IR AGN criteria for the z<2.5 regime, and develop another IR criterion which reliably selects AGN hosts at 0<z<7 (using K, Spitzer-IRAC, and Spitzer-MIPS24um bands, KIM). The ability to track AGN activity since the end of reionization holds great advantages for the study of galaxy evolution. The thesis then focus on the importance of dust. Based on deep IR data on the Cosmological Survey, we derive rest-frame 1.6, 3.3, and 6.2um luminosity functions and their dependency on redshift. We estimate the dust contribution to those wavelengths and show that the hot dust luminosity density evolves since z=1-2 with a much steeper drop than the star-formation history of the Universe. Future prospects are finally discussed in the last chapter.

Ultra Steep Spectrum (USS) radio sources have been successfully used to select powerful radio sources at high redshifts (z>~2). Typically restricted to large-sky surveys and relatively bright radio flux densities, it has gradually become possible to extend the USS search to sub-mJy levels, thanks to the recent appearance of sensitive low-frequency radio facilities. Here a first detailed analysis of the nature of the faintest USS sources is presented. By using Giant Metrewave Radio Telescope and Very Large Array radio observations of the Lockman Hole at 610 MHz and 1.4 GHz, a sample of 58 USS sources, with 610 MHz integrated fluxes above 100 microJy, is assembled. Deep infrared data at 3.6 and 4.5 micron from the Spitzer Extragalactic Representative Volume Survey (SERVS) is used to reliably identify counterparts for 48 (83%) of these sources, showing an average total magnitude of [3.6](AB)=19.8 mag. Spectroscopic redshifts for 14 USS sources, together with photometric redshift estimates, improved by the use of the deep SERVS data, for a further 19 objects, show redshifts ranging from z=0.1 to z=2.8, peaking at z~0.6 and tailing off at high redshifts. The remaining 25 USS sources, with no redshift estimate, include the faintest [3.6] magnitudes, with 10 sources undetected at 3.6 and 4.5 micron (typically [3.6]>22-23 mag, from local measurements), which suggests the likely existence of higher redshifts among the sub-mJy USS population. The comparison with the Square Kilometre Array Design Studies Simulated Skies models indicate that Fanaroff-Riley type I radio sources and radio-quiet Active Galactic Nuclei may constitute the bulk of the faintest USS population, and raises the possibility that the high efficiency of the USS technique for the selection of high redshift sources remains even at the sub-mJy level.

With this paper, we release accurate photometric redshifts for 1692 counterparts to Chandra sources in the central square degree of the COSMOS field. The availability of a large training set of spectroscopic redshifts that extends to faint magnitudes enabled photometric redshifts comparable to the highest quality results presently available for normal galaxies. We demonstrate that morphologically extended, faint X-ray sources without optical variability are more accurately described by a library of normal galaxies (corrected for emission lines) than by AGN-dominated templates, even if these sources have AGN-like X-ray luminosities. Preselecting the library on the bases of the source properties allowed us to reach an accuracy sigma_(Delta z/(1+z_spec)) \sim0.015 with a fraction of outliers of 5.8% for the entire Chandra-COSMOS sample. In addition, we release revised photometric redshifts for the 1735 optical counterparts of the XMM-detected sources over the entire 2 sq. deg.of COSMOS. For 248 sources, our updated photometric redshift differs from the previous release by Delta z>0.2. These changes are predominantly due to the inclusion of newly available deep H-band photometry H_AB=24 mag. We illustrate once again the importance of a spectroscopic training sample and how an assumption about the nature of a source together with the number and the depth of the available bands influence the accuracy of the photometric redshifts determined for AGN. These considerations should be kept in mind when defining the observational strategies of upcoming large surveys targeting AGN, such as eROSITA at X-ray energies and ASKAP/EMU in the radio band.

EMU is a wide-field radio continuum survey planned for the new Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The primary goal of EMU is to make a deep (rms ~ 10 microJy/beam) radio continuum survey of the entire Southern Sky at 1.3 GHz, extending as far North as +30 degrees declination, with a resolution of 10 arcsec. EMU is expected to detect and catalogue about 70 million galaxies, including typical star-forming galaxies up to z~1, powerful starbursts to even greater redshifts, and AGNs to the edge of the visible Universe. It will undoubtedly discover new classes of object. This paper defines the science goals and parameters of the survey, and describes the development of techniques necessary to maximise the science return from EMU.

We present a multi-wavelength analysis of the properties of Extremely Red Galaxy (ERG) populations, selected in the GOODS-South/Chandra Deep Field South field. By using all the photometric and spectroscopic information available on large deep samples of EROs (645 sources), IEROs (294 sources), and DRGs (350 sources), we derive redshift distributions, identify AGN powered and Star-formation powered galaxies, and, using the radio observations of this field, estimate robust (AGN- and dust-unbiased) Star Formation Rate Densities (SFRD) for these populations. We also investigate the properties of "pure" (galaxies that conform to only one of the three ERG criteria considered) and "combined" (galaxies that verify all three criteria) sub-populations. Overall, a large number of AGN are identified (up to ~30%, based on X-rays, and mid-infrared criteria), the majority of which are type-2 (obscured) objects. Among ERGs with no evidence for AGN activity, we identify sub-populations covering a wide range of average star-formation rates, from below 10 Mo/yr to as high as 200 Mo/yr. Applying a redshift separation (1<z<2 and 2<z<3) we find significant evolution (an increase of a factor of 2 or higher) of SFRD for EROs and DRGs, while none is observed for IEROs. The former populations can contribute more than 20\% to the global SFRD at 2<z<3. The emission from AGN activity is typically not strong in the ERG population, with AGN increasing the average radio luminosity of ERG sub-populations by, nominally, less than 20%. AGN are common, however, and, if no discrimination is attempted, this could significantly increase the SFRD estimate (by over 100% in some cases). Thus, and while the contribution of star forming processes to the radio luminosity in galaxies with AGN remains uncertain, a comprehensive identification of AGN in these populations is necessary to obtain meaningful results.