We investigate the ultraviolet (UV) spectral properties of faint Lyman-$\alpha$ emitters (LAEs) in the redshift range 2.9<z<4.6 and provide material to prepare future observations of the faint Universe. We use data from the MUSE Hubble Ultra Deep Survey to construct mean rest-frame spectra of continuum-faint (median M$_{UV}$ of -18 and down to M$_{UV}$ of -16), low stellar mass (median value of $10^{8.4}$ and down to $10^{7}M_{\odot}$) LAEs at redshift z>3. We compute various averaged spectra of LAEs sub-sampled on the basis of their observational (e.g., Ly$\alpha$ strength, UV magnitude and spectral slope) and physical (e.g., stellar mass and star-formation rate) properties. We search for UV spectral features other than Ly$\alpha$, such as higher-ionization nebular emission lines and absorption features. We successfully observe the OIII]1666 and CIII]909 collisionally excited emission lines and the HeII1640 recombination feature, as well as the resonant CIV1550 doublet either in emission or P-Cygni. We compare the observed spectral properties of the different mean spectra and find the emission lines to vary with the observational and physical properties of the LAEs. In particular, the mean spectra of LAEs with larger Ly$\alpha$ equivalent widths, fainter UV magnitudes, bluer UV spectral slopes and lower stellar masses show the strongest nebular emission. The line ratios of these lines are similar to those measured in the spectra of local metal-poor galaxies, while their equivalent widths are weaker compared to the handful of extreme values detected in individual spectra of z>2 galaxies. This suggests that weak UV features are likely ubiquitous in high z, low-mass and faint LAEs. We publicly release the stacked spectra as they can serve as empirical templates for the design of future observations, such as those with the James Webb Space Telescope and the Extremely Large Telescope.

Obscured AGN represent a significant fraction of the entire AGN population, especially at high redshift (~70% at z=3--5). They are often characterized by the presence of large gas and dust reservoirs that are thought to sustain and possibly obscure vigorous star formation processes that make these objects shine at far-IR and sub-mm wavelengths. We exploit ALMA Cycle 4 observations of the continuum (~2.1mm) and high-J CO emission of a sample of six X-ray selected SMGs hosting an obscured AGN at z_spec>2.5 in the 7 Ms CDF-S. We measured the masses and sizes of the dust and molecular gas and we derived the gas density and column density on the basis of a uniform sphere geometry. Finally, we compared the measured column densities with those derived from the Chandra X-ray spectra. We detected both the continuum and line emission for three sources for which we measured both the flux density and size. For the undetected sources, we derived an upper limit on the flux density. We found that the detected galaxies are rich in gas and dust (molecular gas mass in the range <0.5 - 2.7 x 10^10 M_sun for {\alpha}_CO=0.8 and up to ~2 x 10^11~M_sun for {\alpha}_CO=6.5, and dust mass <0.9 - 4.9 x 10^8 M_sun) and compact (gas major axis 2.1-3.0 kpc, dust major axis 1.4-2.7 kpc). The column densities associated with the ISM are on the order of 10^(23-24) cm-2, which is comparable with those derived from the X-ray spectra. For the detected sources we also derived dynamical masses in the range 0.8 - 3.7 x 10^10 M_sun. We conclude that the ISM of high redshift galaxies can substantially contribute to nuclear obscuration up to the Compton-thick (>10^24 cm-2) regime. In addition, we found that all the detected sources show a velocity gradient reminding one rotating system, even though two of them show peculiar features in their morphology that can be associated with a chaotic, possibly merging, structure.

We explore the production and escape of ionizing photons in young galaxies by investigating the ultraviolet and optical emission-line properties of models of ionization-bounded and density-bounded HII regions, active-galactic-nucleus (AGN) narrow-line regions and radiative shocks computed all using the same physically-consistent description of element abundances and depletion on to dust grains down to very low metallicities. We compare these models with a reference sample of metal-poor star-forming galaxies and Lyman-continuum (LyC) leakers at various redshifts, which allows the simultaneous exploration of more spectral diagnostics than typically available at once for individual subsamples. We confirm that current single- and binary-star population synthesis models do not produce hard-enough radiation to account for the high-ionization emission of the most metal-poor galaxies. Introducing either an AGN or radiative-shock component brings models into agreement with observations. A published model including X-ray binaries is an attractive alternative to reproduce the observed rise in HeII4686/Hbeta ratio with decreasing oxygen abundance in metal-poor star-forming galaxies, but not the high observed HeII4686/Hbeta ratios of galaxies with large EW(Hbeta). A source of harder ionizing radiation appears to be required in these extreme objects, such as an AGN or radiative-shock component, perhaps linked to an initial-mass-function bias toward massive stars at low metallicity. This would also account for the surprisingly high [OI]/[OIII] ratios of confirmed LyC leakers relative to ionization-bounded models. We find no simple by-eye diagnostic of the nature of ionizing sources and the escape of LyC photon, which require proper simultaneous fits of several lines to be discriminated against.

The well-known quasar SDSS J095253.83+011421.9 (J0952+0114) at z=3.02 has one of the most peculiar spectra discovered so far, showing the presence of narrow Ly$\alpha$ and broad metal emission lines. Although recent studies have suggested that a Proximate Damped Ly$\alpha$ system (PDLA) causes this peculiar spectrum, the origin of the gas associated with the PDLA is unknown. Here we report the results of MUSE observations that reveal a new giant ($\approx$ 100 physical kpc) Lyman $\alpha$ nebula. The detailed analysis of the Ly$\alpha$ velocity, velocity dispersion, and surface brightness profiles suggests that the J0952+0114 Ly$\alpha$ nebula shares similar properties of other QSO nebulae previously detected with MUSE, implying that the PDLA in J0952+0144 is covering only a small fraction of the QSO emission solid angle. We also detected bright and spectrally narrow CIV$\lambda$1550 and HeII$\lambda$1640 extended emission around J0952+0114 with velocity centroids similar to the peak of the extended and central narrow Ly$\alpha$ emission. The presence of a peculiarly bright, unresolved, and relatively broad HeII$\lambda$1640 emission in the central region at exactly the same PDLA redshift hints at the possibility that the PDLA originates in a clumpy outflow with a bulk velocity of about 500 km/s. The smaller velocity dispersion of the large scale Ly$\alpha$ emission suggests that the high-speed outflow is confined to the central region. Lastly, the derived spatially resolved HeII/Ly$\alpha$ and CIV/Ly$\alpha$ maps show a positive gradient with the distance to the QSO hinting at a non-homogeneous ionization parameter distribution.

We select a sample of 90 obscured (type2) AGN with 1.45<z<3.05 from the zCOSMOS-deep galaxy sample by 5 sigma-detection of the high-ionization CIV {\lambda}1549 narrow emission line. The presence of this feature in a galaxy spectrum is often associated with nuclear activity, and the selection effectiveness has been also confirmed by ultraviolet (UV) emission line ratio diagnostic diagrams. Applying the same selection technique, a sample of 102 unobscured (type 1) AGN was collected. Taking advantage of the large amount of multi-band data available in the COSMOS field, we investigate the properties of the CIV-selected type 2 AGN, focusing on their host galaxies, X-ray emission and UV emission lines. Finally, we investigate the physical properties of the ionized gas in the Narrow Line Region (NLR) of this type 2 AGN sample, combining the analysis of strong UV emission lines with predictions from photo-ionization models. We find that, in order to successfully reproduce the relative intensity of UV emission lines of the selected high-z type 2 AGN, two new ingredients in the photo-ionization models are fundamental,i.e. small inner radii of the NLR (~90pc for LAGN = 10^45erg/s) and the internal dissipative micro-turbulence of the gas emitting clouds (with vmicr~100km/s). With these modified models, we compute the gas-phase metallicity of the NLR, and our measurements indicate a statistically significant evolution of the metal content with redshift. Finally, we do not observe, in our CIV-selected type 2 AGN sample, a strong relationship between the NLR gas metallicity and the stellar mass of the host galaxy.

Deep optical spectroscopic surveys of galaxies provide us a unique opportunity to investigate rest-frame ultra-violet (UV) emission line properties of galaxies at ${z \sim 2-4.5}$. Here we combine VLT/MUSE Guaranteed Time Observations of the Hubble Deep Field South, Ultra Deep Field, COSMOS, and several quasar fields with other publicly available data from VLT/VIMOS and VLT/FORS2 to construct a catalogue of HeII${\lambda}$1640 emitters at ${z\sim2}$. The deepest areas of our MUSE pointings reach a ${3\sigma}$ line flux limit of 3.1${\times}$10-19 erg/ s/ cm$^2$. After discarding broad line active galactic nuclei we find 13 HeII${\lambda}$1640 detections from MUSE with a median MUV = $-20.1$ and 21 tentative HeII${\lambda}$1640 detections from other public surveys. Excluding Ly${\alpha}$, all except two galaxies in our sample show at least one other rest-UV emission line, with CIII]${\lambda}$1907,${\lambda}$1909 being the most prominent. We use multi-wavelength data available in the Hubble legacy fields to derive basic galaxy properties of our sample via spectral energy distribution fitting techniques. Taking advantage of the high quality spectra obtained by MUSE (${\sim10 - 30}$h of exposure time per pointing), we use photo-ionisation models to study the rest-UV emission line diagnostics of the HeII${\lambda}$1640 emitters. Line ratios of our sample can be reproduced by moderately sub-solar photo-ionisation models, however, we find that including effects of binary stars lead to degeneracies in most free parameters. Even after considering extra ionising photons produced by extreme sub-solar metallicity binary stellar models, photo-ionisation models are unable to reproduce rest-frame HeII${\lambda}$1640 equivalent widths (${\sim}$ 0.2 - 10 A), thus additional mechanisms are necessary in models to match the observed HeII${\lambda}$1640 properties.

We present a multi-wavelength study of seven AGN at spectroscopic redshift >2.5 in the 7 Ms Chandra Deep Field South, selected to have good FIR/sub-mm detections. Our aim is to investigate the possibility that the obscuration observed in the X-rays can be produced by the interstellar medium (ISM) of the host galaxy. Based on the 7 Ms Chandra spectra, we measured obscuring column densities N$_{H, X}$ in excess of 7x10$^{22}$ cm$^{-2}$ and intrinsic X-ray luminosities L$_{X}$>10$^{44}$ erg s$^{-1}$ for our targets, as well as equivalent widths for the Fe K$\alpha$ emission line EW>0.5-1 keV. We built the UV-to-FIR spectral energy distributions by using broad-band photometry from CANDELS and Herschel catalogs. By means of an SED decomposition technique, we derived stellar masses (M$_{*}$~10$^{11}$ Msun), IR luminosities (L$_{IR}$>10$^{12}$ Lsun), star formation rates (SFR~190-1680 Msun yr$^{-1}$) and AGN bolometric luminosities (L$_{bol}$~10$^{46}$ erg s$^{-1}$) for our sample. We used an empirically-calibrated relation between gas masses and FIR/sub-mm luminosities and derived M$_{gas}$~0.8-5.4x10$^{10}$ Msun. High-resolution (0.3-0.7'') ALMA data (when available, CANDELS data otherwise) were used to estimate the galaxy size and hence the volume enclosing most of the ISM under simple geometrical assumptions. These measurements were then combined to derive the column density associated with the ISM of the host, on the order of N$_{H, ISM}$~10$^{23-24}$ cm$^{-2}$. The comparison between the ISM column densities and those measured from the X-ray spectral analysis shows that they are similar. This suggests that, at least at high redshift, significant absorption on kpc scales by the dense ISM in the host likely adds to or substitutes that produced by circumnuclear gas on pc scales (i.e., the torus of unified models). The lack of unobscured AGN among our ISM-rich targets supports this scenario.

With a projected size of about 450 kpc at z~2.3, the Slug Ly-alpha nebula is a rare laboratory to study, in emission, the properties of the intergalactic gas in the Cosmic Web. Since its discovery, the Slug has been the subject of several spectroscopic follow-ups to constrain the properties of the emitting gas. Here we report the results of a deep MUSE integral-field spectroscopic search for non-resonant, extended HeII1640 and metal emission. Extended HeII radiation is detected on scales of about 100 kpc, but only in some regions associated with the bright Ly-alpha emission and a continuum-detected source, implying large and abrupt variations in the line ratios across adjacent regions in projected space. The recent detection of associated H-alpha emission and similar abrupt variations in the Ly-alpha kinematics, strongly suggest that the HeII/Ly-alpha gradient is due to large variations in the physical distances between the associated quasar and these regions. This implies that the overall length of the emitting structure could extend to physical Mpc scales and be mostly oriented along our line of sight. At the same time, the relatively low HeII/Ly-alpha values suggest that the emitting gas has a broad density distribution that - if expressed in terms of a lognormal - implies dispersions as high as those expected in the interstellar medium of galaxies. These results strengthen the possibility that the density distribution of intergalactic gas at high-redshift is extremely clumpy and multiphase on scales below our current observational spatial resolution of a few physical kpc.

We compute synthetic optical and ultraviolet (UV) emission-line properties of galaxies in a full cosmological framework by coupling, in post-processing, new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies. Our self-consistent modelling accounts for nebular emission from young stars and accreting black holes (BHs). We investigate which optical- and UV-line diagnostic diagrams can best help to discern between the main ionizing sources, as traced by the ratio of BH accretion to star formation rates in model galaxies, over a wide range of redshifts. At low redshift, simulated star-forming galaxies, galaxies dominated by active galactic nuclei and composite galaxies are appropriately differentiated by standard selection criteria in the classical [OIII]$\lambda$5007/H$\beta$ versus [NII]$\lambda$6584/H$\alpha$ diagram. At redshifts $z \gt 1$, however, this optical diagram fails to discriminate between active and inactive galaxies at metallicities below $0.5\ Z_\odot$. To robustly classify the ionizing radiation of such metal-poor galaxies, which dominate in the early Universe, we confirm 3 previous, and propose 11 novel diagnostic diagrams based on equivalent widths and luminosity ratios of UV emission lines, such as EW(OIII]$\lambda$1663) versus OIII]$\lambda$1663/HeII$\lambda$1640, CIII]$\lambda$1908/HeII$\lambda$1640 versus OIII]$\lambda$1663/HeII$\lambda$1640, and CIV$\lambda$1550/CIII]$\lambda$1908 versus CIII]$\lambda$1908/CII$\lambda$2326. We formulate associated UV selection criteria and discuss some caveats of our results (e.g., uncertainties in the modelling of the HeII$\lambda$1640 line). These UV diagnostic diagrams are potentially important for the interpretation of high-quality spectra of very distant galaxies to be gathered by next-generation telescopes, such as the James Webb Space Telescope.

The physical origin of the near-ultraviolet MgII emission remains an under-explored domain, contrary to more typical emission lines detected in the spectra of star-forming galaxies. We explore the nebular and physical properties for a sample of 381 galaxies between 0.70 < z < 2.34 drawn from the MUSE Hubble Ultra Deep Survey. The spectra of these galaxies show a wide variety of profiles of the MgII 2796,2803 resonant doublet, from absorption to emission. We present a study on the main drivers for the detection of MgII emission in galaxy spectra. By exploiting photoionization models we verified that the emission-line ratios observed in galaxies with MgII in emission are consistent with nebular emission from HII regions. From a simultaneous analysis of MUSE spectra and ancillary HST information via spectral energy distribution (SED) fitting, we find that galaxies with MgII in emission have lower stellar masses, smaller sizes, bluer spectral slopes and lower optical depth than those with absorption. This leads us to suggest that MgII emission is a potential tracer of physical conditions not merely related to those of the ionized gas. We show that these differences in MgII emission/absorption can be explained in terms of a higher dust and neutral gas content in the interstellar medium (ISM) of galaxies showing MgII in absorption, confirming the extreme sensitivity of MgII to the presence of the neutral ISM. We conclude with an analogy between the MgII doublet and the Ly-alpha line, due to their resonant nature. Further investigations with current and future facilities, including JWST, are promising as the detection of MgII emission and its potential connection with Ly-alpha could provide new insights on the ISM content in the early Universe.

Measurements of the galaxy UV luminosity function at z>6 suggest that young stars hosted in low-mass star-forming galaxies produced the bulk of hydrogen-ionizing photons necessary to reionize the intergalactic medium (IGM) by redshift z~6. Whether star-forming galaxies dominated cosmic reionization, however, also depends on their stellar populations and interstellar medium properties, which set, among other things, the production rate of H-ionizing photons, $\xi_\text{ion}^\star$, and the fraction of these escaping into the IGM. Given the difficulty of constraining with existing observatories the physical properties of z>6 galaxies, in this work we focus on a sample of ten nearby objects showing UV spectral features comparable to those observed at z>6. We use the new-generation Beagle tool to model the UV-to-optical photometry and UV/optical emission lines of these Local 'analogues' of high-redshift galaxies, finding that our relatively simple, yet fully self-consistent, physical model can successfully reproduce the different observables considered. Our galaxies span a broad range of metallicities and are characterised by high ionization parameters, low dust attenuation, and very young stellar populations. Through our analysis, we derive a novel diagnostic of the production rate of H-ionizing photons per unit UV luminosity, $\xi_\text{ion}^\star$, based on the equivalent width of the bright $[\text{OIII}] \lambda 4959,5007$ doublet, which does not require measurements of H-recombination lines. This new diagnostic can be used to estimate $\xi_\text{ion}^\star$ from future direct measurements of the $[\text{OIII}] \lambda 4959,5007$ line using JWST/NIRSpec (out to z~9.5), and by exploiting the contamination by $\text{H}\beta + [\text{OIII}] \lambda 4959,5007$ of photometric observations of distant galaxies, for instance from existing Spitzer/IRAC data and from future ones with JWST/NIRCam.

Non-resonant FeII* 2365, 2396, 2612, 2626 emission can potentially trace galactic winds in emission and provide useful constraints to wind models. From the 3'x3' mosaic of the Hubble Ultra Deep Field (UDF) obtained with the VLT/MUSE integral field spectrograph, we identify a statistical sample of 40 FeII* emitters and 50 MgII 2796, 2803 emitters from a sample of 271 [OII] 3726, 3729 emitters with reliable redshifts from z = 0.85 - 1.5 down to 2E-18 (3 sigma) ergs/s/cm^2 (for [OII]), covering the stellar mass range 10^8 - 10^11 Msun. The FeII* and MgII emitters follow the galaxy main sequence, but with a clear dichotomy. Galaxies with masses below 10^9 Msun and star formation rates (SFRs) of <1 Msun/year have MgII emission without accompanying FeII* emission, whereas galaxies with masses above 10^10 Msun and SFRs >10 Msun/year have FeII* emission without accompanying MgII emission. Between these two regimes, galaxies have both MgII and FeII* emission, typically with MgII P-Cygni profiles. Indeed, the MgII profile shows a progression along the main sequence from pure emission to P-Cygni profiles to strong absorption, due to resonant trapping. Combining the deep MUSE data with HST ancillary information, we find that galaxies with pure MgII emission profiles have lower star formation rate surface densities than those with either MgII P-Cygni profiles or FeII* emission. These spectral signatures produced through continuum scattering and fluorescence, MgII P-Cygni profiles and FeII* emission, are better candidates for tracing galactic outflows than pure MgII emission, which may originate from HII regions. We compare the absorption and emission rest-frame equivalent widths for pairs of FeII transitions to predictions from outflow models and find that the observations consistently have less total re-emission than absorption, suggesting either dust extinction or non-isotropic outflow geometries.

The C III] 1907,1909 emission doublet has been proposed as an alternative to Lyman-alpha in redshift confirmations of galaxies at z > 6 since it is not attenuated by the largely neutral intergalactic medium at these redshifts and is believed to be strong in the young, vigorously star-forming galaxies present at these early cosmic times. We present a statistical sample of 17 C III]-emitting galaxies beyond z~1.5 using 30 hour deep VLT/MUSE integral field spectroscopy covering 2 square arcminutes in the Hubble Deep Field South (HDFS) and Ultra Deep Field (UDF), achieving C III] sensitivities of ~2e-17 erg/s/cm^2 in the HDFS and ~7e-18 erg/s/cm^2 in the UDF. The rest-frame equivalent widths range from 2 to 19 Angstroms. These 17 galaxies represent ~3% of the total sample of galaxies found between 1.5 < z < 4. They also show elevated star formation rates, lower dust attenuation, and younger mass-weighted ages than the general population of galaxies at the same redshifts. Combined with deep slitless grism spectroscopy from the HST/WFC3 in the UDF, we can tie the rest-frame ultraviolet C III] emission to rest-frame optical emission lines, namely [O III] 5007, finding a strong correlation between the two. Down to the flux limits that we observe (~1e-18 erg/s/cm^2 with the grism data in the UDF), all objects with a rest-frame [O III] 4959,5007 equivalent width in excess of 250 Angstroms, the so-called Extreme Emission Line Galaxies, have detections of C III] in our MUSE data. More detailed studies of the C III]-emitting population at these intermediate redshifts will be crucial to understand the physical conditions in galaxies at early cosmic times and to determine the utility of C III] as a redshift tracer.

Galaxies occupy different regions of the [OIII]$\lambda5007$/H$\beta$-versus-[NII]$\lambda6584$/H$\alpha$ emission-line ratio diagram in the distant and local Universe. We investigate the origin of this intriguing result by modelling self-consistently, for the first time, nebular emission from young stars, accreting black holes (BHs) and older, post-asymptotic-giant-branch (post-AGB) stellar populations in galaxy formation simulations in a full cosmological context. In post-processing, we couple new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies to explore which galaxy physical properties drive the cosmic evolution of the optical-line ratios [OIII]$\lambda5007$/H$\beta$, [NII]$\lambda6584$/H$\alpha$, [SII]$\lambda\lambda6717,6731$/H$\alpha$ and [OI]$\lambda6300$/H$\alpha$. The line ratios of simulated galaxies agree well with observations of both star-forming and active local SDSS galaxies. Towards higher redshifts, at fixed galaxy stellar mass, the average [OIII]/H$\beta$ increases and [NII]/H$\alpha$, [SII]/H$\alpha$ and [OI]/H$\alpha$ decrease -- widely consistent with observations. At fixed stellar mass, we identify star formation history, which controls nebular emission from young stars via the ionization parameter, as the primary driver of the cosmic evolution of [OIII]/H$\beta$ and [NII]/H$\alpha$. For [SII]/H$\alpha$ and [OI]/H$\alpha$, this applies only to redshifts above $z=1.5$, the evolution at lower redshift being driven in roughly equal parts by nebular emission from AGN and post-AGB stars. Instead, changes in the hardness of ionizing radiation, ionized-gas density, the prevalence of BH accretion relative to star formation and the dust-to-metal mass ratio (whose impact on the gas-phase N/O ratio we model at fixed O/H) play at most a minor role in the cosmic evolution of simulated galaxy line ratios.

Nearby dwarf galaxies provide a unique laboratory in which to test stellar population models below $Z_\odot/2$. Such tests are particularly important for interpreting the surprising high-ionization UV line emission detected at $z>6$ in recent years. We present HST/COS ultraviolet spectra of ten nearby metal-poor star-forming galaxies selected to show He II emission in SDSS optical spectra. The targets span nearly a dex in gas-phase oxygen abundance ($7.8<12+\log\mathrm{O/H}<8.5$) and present uniformly large specific star formation rates (sSFR $\sim 10^2$ $\mathrm{Gyr}^{-1}$). The UV spectra confirm that metal-poor stellar populations can power extreme nebular emission in high-ionization UV lines, reaching C III] equivalent widths comparable to those seen in systems at $z\sim 6-7$. Our data reveal a marked transition in UV spectral properties with decreasing metallicity, with systems below $12+\log\mathrm{O/H}\lesssim 8.0$ ($Z/Z_\odot \lesssim 1/5$) presenting minimal stellar wind features and prominent nebular emission in He II and C IV. This is consistent with nearly an order of magnitude increase in ionizing photon production beyond the $\mathrm{He^+}$-ionizing edge relative to H-ionizing flux as metallicity decreases below a fifth solar, well in excess of standard stellar population synthesis predictions. Our results suggest that often neglected sources of energetic radiation such as stripped binary products and very massive O-stars produce a sharper change in the ionizing spectrum with decreasing metallicity than expected. Consequently, nebular emission in C IV and He II powered by these stars may provide useful metallicity constraints in the reionization era.

We explore the relationship between active galactic nuclei and star formation in a sample of 513 optically luminous type 1 quasars up to redshifts of $\sim$4 hosting extremely high star formation rates (SFRs). The quasars are selected to be individually detected by the \textit{Herschel} SPIRE instrument at $> $3$\sigma$ at 250 $\mu$m, leading to typical SFRs of order of 1000 M$_{\odot}$yr$^{-1}$. We find the average SFRs to increase by almost a factor 10 from $z\sim0.5$ to $z\sim3$, mirroring the rise in the comoving SFR density over the same epoch. However, we find that the SFRs remain approximately constant with increasing accretion luminosity for accretion luminosities above 10$^{12}$ L$_{\odot}$. We also find that the SFRs do not correlate with black hole mass. Both of these results are most plausibly explained by the existence of a self-regulation process by the starburst at high SFRs, which controls SFRs on time-scales comparable to or shorter than the AGN or starburst duty cycles. We additionally find that SFRs do not depend on Eddington ratio at any redshift, consistent with no relation between SFR and black hole growth rate per unit black hole mass. Finally, we find that high-ionisation broad absorption line (HiBAL) quasars have indistinguishable far-infrared properties to those of classical quasars, consistent with HiBAL quasars being normal quasars observed along a particular line of sight, with the outflows in HiBAL quasars not having any measurable effect on the star formation in their hosts.

We investigate the relation between star formation rates ($\dot{M}_{s}$) and AGN properties in optically selected type 1 quasars at $2<z<3$ using data from Herschel and the SDSS. We find that $\dot{\rm{M}}_s$ remains approximately constant with redshift, at $300\pm100~\rm{M}_{\odot}$yr$^{-1}$. Conversely, $\dot{\rm{M}}_s$ increases with AGN luminosity, up to a maximum of $\sim600~\rm{M}_{\odot}$yr$^{-1}$, and with CIV FWHM. In context with previous results, this is consistent with a relation between $\dot{\rm{M}}_s$ and black hole accretion rate ($\dot{\rm{M}}_{bh}$) existing in only parts of the $z-\dot{\rm{M}}_{s}-\dot{\rm{M}}_{bh}$ plane, dependent on the free gas fraction, the trigger for activity, and the processes that may quench star formation. The relations between $\dot{\rm{M}}_s$ and both AGN luminosity and CIV FWHM are consistent with star formation rates in quasars scaling with black hole mass, though we cannot rule out a separate relation with black hole accretion rate. Star formation rates are observed to decline with increasing CIV equivalent width. This decline can be partially explained via the Baldwin effect, but may have an additional contribution from one or more of three factors; $M_i$ is not a linear tracer of L$_{2500}$, the Baldwin effect changes form at high AGN luminosities, and high CIV EW values signpost a change in the relation between $\dot{\rm{M}}_s$ and $\dot{\rm{M}}_{bh}$. Finally, there is no strong relation between $\dot{\rm{M}}_s$ and Eddington ratio, or the asymmetry of the CIV line. The former suggests that star formation rates do not scale with how efficiently the black hole is accreting, while the latter is consistent with CIV asymmetries arising from orientation effects.

In the context of observations of the rest-frame ultraviolet and optical emission from distant galaxies, we explore the emission-line properties of photoionization models of active and inactive galaxies. Our aim is to identify new line-ratio diagnostics to discriminate between gas photoionization by active galactic nuclei (AGN) and star formation. We use a standard photoionization code to compute the emission from AGN narrow-line regions and compare this with calculations of the nebular emission from star-forming galaxies achieved using the same code. We confirm the appropriateness of widely used optical spectral diagnostics of nuclear activity versus star formation and explore new diagnostics at ultraviolet wavelengths. We find that combinations of a collisionally excited metal line or line multiplet, such as CIV 1548,1551, OIII]1661,1666, NIII]1750, [SiIII]1883+[SiIII]1892 and [CIII]1907+CIII]1909, with the HeII 1640 recombination line are individually good discriminants of the nature of the ionizing source. Diagrams involving at least 3 of these lines allow an even more stringent distinction between active and inactive galaxies, as well as valuable constraints on interstellar gas parameters and the shape of the ionizing radiation. Several line ratios involving Ne-based emission lines, such as [NeIV]2424, [NeIII]3343 and [NeV]3426, are also good diagnostics of nuclear activity. Our results provide a comprehensive framework to identify the sources of photoionization and physical conditions of the ionized gas from the ultraviolet and optical nebular emission from galaxies. This will be particularly useful to interpret observations of high-redshift galaxies with future facilities, such as the James Webb Space Telescope and extremely large ground-based telescopes.

We used wide area surveys over 39 deg$^2$ by the HerMES collaboration, performed with the Herschel Observatory SPIRE multi-wavelength camera, to estimate the low-redshift, $0.02<z<0.5$, monochromatic luminosity functions (LFs) of galaxies at 250, 350 and 500$\,\mu$m. SPIRE flux densities were also combined with Spitzer photometry and multi-wavelength archival data to perform a complete SED fitting analysis of SPIRE detected sources to calculate precise k-corrections, as well as the bolometric infrared (8-1000$\,\mu$m) luminosity functions and their low-$z$ evolution from a combination of statistical estimators. Integration of the latter prompted us to also compute the local luminosity density (LLD) and the comoving star formation rate density (SFRD) for our sources, and to compare them with theoretical predictions of galaxy formation models. The luminosity functions show significant and rapid luminosity evolution already at low redshifts, $0.02<z<0.2$, with L$_{IR}^* \propto (1+z)^{6.0\pm0.4}$ and $\Phi_{IR}^* \propto (1+z)^{-2.1\pm0.4}$, L$_{250}^* \propto (1+z)^{5.3\pm0.2}$ and $\Phi_{250}^* \propto (1+z)^{-0.6\pm0.4}$ estimated using the IR bolometric and the 250$\,\mu$m LFs respectively. Converting our IR LD estimate into an SFRD assuming a standard Salpeter IMF and including the unobscured contribution based on the UV dust-uncorrected emission from local galaxies, we estimate a SFRD scaling of SFRD$_0+0.08 z$, where SFRD$_0\simeq (1.9\pm 0.03)\times 10^{-2} [\mathrm{M}_\odot\,\mathrm{Mpc}^{-3}]$ is our total SFRD estimate at $z\sim0.02$.

We present Keck/MOSFIRE observations of UV metal lines in four bright gravitationally-lensed z~6-8 galaxies behind the cluster Abell 1703. The spectrum of A1703-zd6, a highly-magnified star forming galaxy with a Lyman-alpha redshift of z=7.045, reveals a confident detection of the nebular CIV emission line (unresolved with FWHM < 125 km/s). UV metal lines are not detected in the three other galaxies. At z~2-3, nebular CIV emission is observed in just 1% of UV-selected galaxies. The presence of strong CIV emission in one of the small sample of galaxies targeted in this paper may indicate hard ionizing spectra are more common at z~7. The total estimated equivalent width of the CIV doublet (38 A) and CIV/Lyman-alpha flux ratio (0.3) are comparable to measurements of narrow-lined AGNs. Photoionization models show that the nebular CIV line can also be reproduced by a young stellar population, with very hot metal poor stars dominating the photon flux responsible for triply ionizing carbon. Regardless of the origin of the CIV, we show that the ionizing spectrum of A1703-zd6 is different from that of typical galaxies at z~2, producing more H ionizing photons per unit 1500A luminosity and a larger flux density at 30-50 eV. If such extreme radiation fields are typical in UV-selected systems at z>7, it would indicate that reionization-era galaxies are more efficient ionizing agents than previously thought. Alternatively, we suggest that the small sample of Lyman-alpha emitters at z>7 may trace a rare population with intense radiation fields capable of ionizing their surrounding hydrogen distribution. Additional constraints on high ionization emission lines in galaxies with and without Lyman-alpha detections will help clarify whether hard ionizing spectra are common in the reionization era.

We present a comprehensive study of the silicate features at 9.7 and 18 micron of a sample of almost 800 active galactic nuclei (AGN) with available spectra from the Spitzer InfraRed Spectrograph (IRS). We measure the strength of the silicate feature at 9.7 micron, S9.7, before and after subtracting the host galaxy emission from the IRS spectra. The numbers of type 1 and 2 AGN with the feature in emission increase by 20 and 50%, respectively, once the host galaxy is removed, while 35% of objects with the feature originally in absorption exhibit it in even deeper absorption. The peak of S9.7, lambda_peak, has a bimodal distribution when the feature is in emission, with about 65% of the cases showing lambda_peak > 10.2 micron. Silicates can appear in emission in objects with mid-infrared (MIR) luminosity spanning over six orders of magnitude. The derived distributions of the strength of the silicate features at 9.7 and 18 micron provide a solid test bed for modeling the dust distribution in AGN. Clumpiness is needed in order to produce absorption features in unobscured AGN and can also cause the silicates to be in absorption at 9.7 micron and in emission at 18 micron in type 1 sources. We find the `cosmic' silicates of Ossenkopf et al. to be more consistent with the observations than Draine's `astronomical' silicates. Finally, we discuss the possibility of a foreground absorber to explain the deep silicate absorption features in the MIR spectra of some type 2 AGN.

We study a sample of Herschel-PACS selected galaxies within the GOODS-South and the COSMOS fields in the framework of the PACS Evolutionary Probe (PEP) project. Starting from the rich multi-wavelength photometric data-sets available in both fields, we perform a broad-band Spectral Energy Distribution (SED) decomposition to disentangle the possible active galactic nucleus (AGN) contribution from that related to the host galaxy. We find that 37 per cent of the Herschel-selected sample shows signatures of nuclear activity at the 99 per cent confidence level. The probability to reveal AGN activity increases for bright ($L_{\rm 1-1000} > 10^{11} \rm L_{\odot}$) star-forming galaxies at $z>0.3$, becoming about 80 per cent for the brightest ($L_{\rm 1-1000} > 10^{12} \rm L_{\odot}$) infrared (IR) galaxies at $z \geq 1$. Finally, we reconstruct the AGN bolometric luminosity function and the super-massive black hole growth rate across cosmic time up to $z \sim 3$ from a Far-Infrared (FIR) perspective. This work shows general agreement with most of the panchromatic estimates from the literature, with the global black hole growth peaking at $z \sim 2$ and reproducing the observed local black hole mass density with consistent values of the radiative efficiency $\epsilon_{\rm rad}$ ($\sim$0.07).

We report ALMA Cycle 0 observations at 1.3mm of LESS J033229.4-275619 (XID403), an Ultraluminous Infrared Galaxy at $z=4.75$ in the Chandra Deep Field South hosting a Compton-thick QSO. The source is not resolved in our data at a resolution of $\sim$0.75 arcsec, placing an upper-limit of 2.5 kpc to the half-light radius of the continuum emission from heated-dust. After deconvolving for the beam size, however, we found a $\sim3\sigma$ indication of an intrinsic source size of $0.27\pm0.08$ arcsec (Gaussian FWHM), which would correspond to $r_{half}\sim0.9\pm0.3$ kpc. We build the far-IR SED of XID403 by combining datapoints from both ALMA and Herschel and fit it with a modified blackbody spectrum. For the first time, we measure the dust temperature $T_d=58.5\pm5.3$ K in this system, which is comparable to what has been observed in other high-z submillimeter galaxies. The measured star formation rate is SFR=$1020\pm150$ $M_{\odot}$ yr$^{-1}$, in agreement with previous estimates at lower S/N. Based on the measured SFR and source size, we constrain the SFR surface density to be $\Sigma_{SFR}>26\;M_{\odot}$yr$^{-1}$kpc$^{-2}$ ($\sim200\;M_{\odot}$yr$^{-1}$kpc$^{-2}$ for $r_{half}\sim0.9$ kpc). The compactness of this starburst is comparable to what has been observed in other local and high-z starburst galaxies. If the gas mass measured from previous [CII] and CO(2-1) observations at low resolution is confined within the same dust region, assuming $r_{half}\sim0.9\pm0.3$ kpc, this would produce a column density of $N_H\sim0.3-1.1\times10^{24}$cm$^{-2}$ towards the central SMBH, similar to the column density of $\approx1.4\times10^{24}$cm$^{-2}$ measured from the X-rays. Then, in principle, if both gas and dust were confined on sub-kpc scales, this would be sufficient to produce the observed X-ray column density without any need of a pc-scale absorber [abridged].

We study the mid- to far-IR properties of a 24um-selected flux-limited sample (S24 > 5mJy) of 154 intermediate redshift (<z>~0.15), infrared luminous galaxies, drawn from the 5MUSES survey. By combining existing mid-IR spectroscopy and new Herschel SPIRE submm photometry from the HerMES program, we derived robust total infrared luminosity (LIR) and dust mass (Md) estimates and infered the relative contribution of the AGN to the infrared energy budget of the sources. We found that the total infrared emission of galaxies with weak 6.2um PAH emission (EW<0.2um) is dominated by AGN activity, while for galaxies with EW>0.2um more than 50% of the LIR arises from star formation. We also found that for galaxies detected in the 250-500um Herschel bands an AGN has a statistically insignificant effect on the temperature of the cold dust and the far-IR colours of the host galaxy, which are primarily shaped by star formation activity. For star-forming galaxies we reveal an anti-correlation between the LIR-to-rest-frame 8um luminosity ratio, IR8 = LIR\L8, and the strength of PAH features. We found that this anti-correlation is primarily driven by variations in the PAHs emission, and not by variations in the 5-15um mid-IR continuum emission. Using the [NeIII]/[NeII] line ratio as a tracer of the hardness of the radiation field, we confirm that galaxies with harder radiation fields tend to exhibit weaker PAH features, and found that they have higher IR8 values and higher dust-mass-weighted luminosities (LIR/Md), the latter being a proxy for the dust temperature (Td). We argue that these trends originate either from variations in the environment of the star-forming regions or are caused by variations in the age of the starburst. Finally, we provide scaling relations that will allow estimating LIR, based on single-band observations with the mid-infrared instrument, on board the upcoming JWST.

In this work we explore the impact of the presence of an active galactic nucleus (AGN) on the mid- and far-infrared (IR) properties of galaxies as well as the effects of simultaneous AGN and starburst activity in these same galaxies. To do this we apply a multi-component, multi-band spectral synthesis technique to a sample of 250 micron selected galaxies of the Herschel Multi-tiered Extragalactic Survey, with IRS spectra available for all galaxies. Our results confirm that the inclusion of the IRS spectra plays a crucial role in the spectral analysis of galaxies with an AGN component improving the selection of the best-fit hot dust model (torus). We find a correlation between the obscured star formation rate (SFR) derived from the IR luminosity of the starburst component, SFR_IR and SFR_PAH, derived from the luminosity of the PAH features, L_PAH, with SFR_FIR taking higher values than SFR_PAH. The correlation is different for AGN- and starburst-dominated objects. The ratio of L_PAH to that of the starburst component, L_PAH/L_SB, is almost constant for AGN-dominated objects but decreases with increasing L_SB for starburst-dominated objects. SFR_FIR increases with the accretion luminosity, L_acc, with the increase less prominent for the very brightest, unobscured AGN-dominated sources. We find no correlation between the masses of the hot and cold dust components. We interpret this as a non-constant fraction of gas driven by the gravitational effects to the AGN while the starburst is ongoing. We also find no evidence of the AGN affecting the temperature of the cold dust component, though this conclusion is mostly based on objects with a non-dominant AGN component. We conclude that our findings do not provide evidence that the presence of AGN affects the star formation process in the host galaxy, but rather that the two phenomena occur simultaneously over a wide range of luminosities.

One of the remaining open issues in the context of the analysis of Active Galactic Nuclei (AGN) is the evidence that nuclear gravitational accretion is often accompanied by a concurrent starburst (SB) activity. What is, in this picture, the role played by the obscuring dust around the nucleus and what do the state of the art AGN torus models have to say? Can the IR data provided by Spitzer and Herschel help us in extensively investigate both phenomena and, if so, how and with what limitations? In this paper we present our contribution to the efforts of answering these questions. We show some of the main results coming from a comparative study of various AGN SED modeling approaches, focusing mostly on the much-debated issue about the morphology of the dust distribution in the toroidal structure surrounding the AGN. We found that the properties of dust in AGN as measured by matching observations (be it broad band IR photometry or IR spectra) with models, strongly depend on the choiceof the dust distribution. Then, we present the spectral energy distribution (SED) fitting procedure we developed, making make the best use of Spitzer and Herschel SPIRE mid- and far-IR observations, to dig into the role played by the possible presence of an AGN on the host galaxy's properties.

We combine far-infrared photometry from Herschel (PEP/HERMES) with deep mid-infrared spectroscopy from Spitzer to investigate the nature and the mass assembly history of a sample of 31 Luminous and Ultraluminous Infrared Galaxies at z~1 and 2 selected in GOODS-S with 24 $\mu$m fluxes between 0.2 and 0.5 mJy. We model the data with a self-consistent physical model (GRASIL) which includes a state-of-the-art treatment of dust extinction and reprocessing. We find that all of our galaxies appear to require massive populations of old (>1 Gyr) stars and, at the same time, to host a moderate ongoing activity of SF (SFR < 100 M$_{\odot}$/yr). The bulk of the stars appear to have been formed a few Gyr before the observation in essentially all cases. Only five galaxies of the sample require a recent starburst superimposed on a quiescent star formation history (SFH). We also find discrepancies between our results and those based on optical-only SED fitting for the same objects; by fitting their observed Spectral Energy Distributions with our physical model we find higher extinctions (by $\Delta$A_{V} ~ 0.81 and 1.14) and higher stellar masses (by $\Delta$Log(M*) ~ 0.16 and 0.36 dex) for z~1 and z~2 (U)LIRGs, respectively. The stellar mass difference is larger for the most dust obscured objects. We also find lower SFRs than those computed from L_{IR} using the Kennicutt relation due to the significant contribution to the dust heating by intermediate-age stellar populations through 'cirrus' emission (~73% and ~66% of total L_{IR} for z~1 and z~2 (U)LIRGs, respectively).

The geometry of the dust distribution within the inner regions of Active Galactic Nuclei (AGN) is still a debated issue and relates directly with the AGN unified scheme. Traditionally, models discussed in the literature assume one of two distinct dust distributions in what is believed to be a toroidal region around the Supermassive Black Holes: a continuous distribution, customarily referred to as smooth, and a concentration of dust in clumps or clouds, referred to as clumpy. In this paper we perform a thorough comparison between two of the most popular models in the literature, namely the smooth models by Fritz. et al. 2006 and the clumpy models by Nenkova et al. 2008a, in their common parameters space. Particular attention is paid to the silicate features at ~9.7 and ~18 micron, the width of the infrared bump, the near-infrared index and the luminosity at 12.3 micron, all previously reported as possible diagnostic tools to distinguish between the two dust distributions. We find that, due to the different dust chemical compositions used in the two models, the behaviour of the silicate features at 9.7 and 18 micron is quite distinct between the two models. The width of the infrared bump and the peak of the infrared emission can take comparable values, their distributions do, however, vary. The near-infrared index is also quite different, due partly to the primary sources adopted by the two models. Models with matched parameters do not produce similar SEDs and virtually no random parameter combinations can result in seemingly identical SEDs.

We use Herschel-PACS far-infrared data, combined with previous multi-band information and mid-IR spectra, to properly account for the presence of an active nucleus and constrain its energetic contribution in luminous infrared (IR) sources at z\sim2. The sample is composed of 24 sources in the GOODS-South field, with typical IR luminosity of 10^{12} Lo. Data from the 4 Ms Chandra X-ray imaging in this field are also used to identify and characterize AGN emission. We reproduce the observed spectral energy distribution (SED), decomposed into a host-galaxy and an AGN component. A smooth-torus model for circum-nuclear dust is used to account for the direct and re-processed contribution from the AGN. We confirm that galaxies with typical L_{8-1000um}\sim10^{12}Lo at z\sim2 are powered predominantly by star-formation. An AGN component is present in nine objects (\sim35% of the sample) at the 3sigma confidence level, but its contribution to the 8-1000 um emission accounts for only \sim5% of the energy budget. The AGN contribution rises to \sim23% over the 5-30 um range (in agreement with Spitzer IRS results) and to \sim60% over the narrow 2-6 um range. The presence of an AGN is confirmed by X-ray data for 3 (out of nine) sources, with X-ray spectral analysis indicating the presence of significant absorption, i.e. NH\sim10^{23}-10^{24} cm^{-2}. An additional source shows indications of obscured AGN emission from X-ray data. The comparison between the mid-IR--derived X-ray luminosities and those obtained from X-ray data suggests that obscuration is likely present also in the remaining six sources that harbour an AGN according to the SED-fitting analysis.

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

We present a long (~76 ks) Chandra observation of IRAS 09104+4109, a hyper-luminous galaxy, optically classified as a Type 2 AGN hosted in a cD galaxy in a cluster at z=0.442. We also report on the results obtained by fitting its broad-band spectral energy distribution. The Compton-thick nature of this source (which has been often referred to as an "archetype" of Compton-thick Type 2 quasars) was formerly claimed on the basis of its marginal detection in the PDS instrument onboard BeppoSAX, being then disputed using XMM-Newton data. Both Chandra analysis and optical/mid-IR spectral fitting are consistent with the presence of heavy (~1-5 10^{23} cm^{-2}) but not extreme (Compton-thick) obscuration. However, using the mid-IR and the [OIII] emission as proxies of the nuclear hard X-ray luminosity suggests the presence of heavier obscuration. The 54-month Swift BAT map shows excess hard X-ray emission likely related to a nearby (z=0.009) Type 2 AGN, close enough to IRAS 09104+4109 to significantly enhance and contaminate its emission in the early BeppoSAX PDS data.

Two main modes of star formation are know to control the growth of galaxies: a relatively steady one in disk-like galaxies, defining a tight star formation rate (SFR)-stellar mass sequence, and a starburst mode in outliers to such a sequence which is generally interpreted as driven by merging. Such starburst galaxies are rare but have much higher SFRs, and it is of interest to establish the relative importance of these two modes. PACS/Herschel observations over the whole COSMOS and GOODS-South fields, in conjunction with previous optical/near-IR data, have allowed us to accurately quantify for the first time the relative contribution of the two modes to the global SFR density in the redshift interval 1.5<z<2.5, i.e., at the cosmic peak of the star formation activity. The logarithmic distributions of galaxy SFRs at fixed stellar mass are well described by Gaussians, with starburst galaxies representing only a relatively minor deviation that becomes apparent for SFRs more than 4 times higher than on the main sequence. Such starburst galaxies represent only 2% of mass-selected star forming galaxies and account for only 10% of the cosmic SFR density at z~2. Only when limited to SFR>1000M(sun)/yr, off-sequence sources significantly contribute to the SFR density (46+/-20%). We conclude that merger-driven starbursts play a relatively minor role for the formation of stars in galaxies, whereas they may represent a critical phase towards the quenching of star formation and morphological transformation in galaxies.