On January 14, 2019, the Major Atmospheric Gamma Imaging Cherenkov telescopes detected GRB 190114C above 0.2 TeV, recording the most energetic photons ever observed from a gamma-ray burst. We use this unique observation to probe an energy dependence of the speed of light in vacuo for photons as predicted by several quantum gravity models. Based on a set of assumptions on the possible intrinsic spectral and temporal evolution, we obtain competitive lower limits on the quadratic leading order of speed of light modification.

Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the external medium generates external shock waves, responsible for the afterglow emission, which lasts from days to months, and occurs over a broad energy range, from the radio to the GeV bands. The afterglow emission is generally well explained as synchrotron radiation by electrons accelerated at the external shock. Recently, an intense, long-lasting emission between 0.2 and 1 TeV was observed from the GRB 190114C. Here we present the results of our multi-frequency observational campaign of GRB~190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from $5\times10^{-6}$ up to $10^{12}$\,eV. We find that the broadband spectral energy distribution is double-peaked, with the TeV emission constituting a distinct spectral component that has power comparable to the synchrotron component. This component is associated with the afterglow, and is satisfactorily explained by inverse Compton upscattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed TeV component are not atypical, supporting the possibility that inverse Compton emission is commonly produced in GRBs.

Gamma-ray bursts (GRBs) of the long-duration class are the most luminous sources of electromagnetic radiation known in the Universe. They are generated by outflows of plasma ejected at near the speed of light by newly formed neutron stars or black holes of stellar mass at cosmological distances. Prompt flashes of MeV gamma rays are followed by longer-lasting afterglow emission from radio waves to GeV gamma rays, due to synchrotron radiation by energetic electrons in accompanying shock waves. Although emission of gamma rays at even higher, TeV energies by other radiation mechanisms had been theoretically predicted, it had never been detected previously. Here we report the clear detection of GRB 190114C in the TeV band, achieved after many years of dedicated searches for TeV emission from GRBs. Gamma rays in the energy range 0.2--1 TeV are observed from about 1 minute after the burst (at more than 50 standard deviations in the first 20 minutes). This unambiguously reveals a new emission component in the afterglow of a GRB, whose power is comparable to that of the synchrotron component. The observed similarity in the radiated power and temporal behaviour of the TeV and X-ray bands points to processes such as inverse Compton radiation as the mechanism of the TeV emission, while processes such as synchrotron emission by ultrahigh-energy protons are disfavoured due to their low radiative efficiency.

Extreme High-frequency BL~Lacs (EHBL) feature their synchrotron peak of the broadband spectral energy distribution (SED) at $\nu_{\rm s} \geq $10$^{17}$\,Hz. The BL~Lac object 1ES~2344+514 was included in the EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In August 2016, 1ES~2344+514 was detected with the ground-based $\gamma$-ray telescope FACT during a high $\gamma$-ray state, triggering multi-wavelength (MWL) observations. We studied the MWL light curves of 1ES~2344+514 during the 2016 flaring state, using data from radio to VHE $\gamma$ rays taken with OVRO, KAIT, KVA, NOT, some telescopes of the GASP-WEBT collaboration at the Teide, Crimean, and St. Petersburg observatories, \textit{Swift}-UVOT, \textit{Swift}-XRT, \textit{Fermi}-LAT, FACT and MAGIC. With simultaneous observations of the flare, we built the broadband SED and studied it in the framework of a leptonic and an hadronic model. The VHE $\gamma$-ray observations show a flux level of 55\% of the Crab Nebula flux above 300\,GeV, similar to the historical maximum of 1995. The combination of MAGIC and \textit{Fermi}-LAT spectra provides an unprecedented characterization of the inverse-Compton peak for this object during a flaring episode. The $\Gamma$ index of the intrinsic spectrum in the VHE $\gamma$-ray band is $2.04\pm0.12_{\rm stat}\pm0.15_{\rm sys}$. We find the source in an extreme state with a shift of the position of the synchrotron peak to frequencies above or equal to $10^{18}$\,Hz

Context. It has become evident that one-zone synchrotron self-Compton models are not always adequate for very-high-energy (VHE) gamma-ray emitting blazars. While two-component models are performing better, they are difficult to constrain due to the large number of free parameters. Aims. In this work, we make a first attempt to take into account the observational constraints from Very Long Baseline Interferometry (VLBI) data, long-term light curves (radio, optical, and X-rays) and optical polarisation to limit the parameter space for a two-component model and test if it can still reproduce the observed spectral energy distribution (SED) of the blazars. Methods. We selected five TeV BL Lac objects based on the availability of VHE gamma-ray and optical polarisation data. We collected constraints for the jet parameters from VLBI observations. We evaluated the contributions of the two components to the optical flux by means of decomposition of long-term radio and optical light curves as well as modeling of the optical polarisation variability of the objects. We selected eight epochs for these five objects, based on the variability observed at VHE gamma rays, for which we constructed the SEDs that we then modeled with a two-component model. Results. We found parameter sets which can reproduce the broadband SED of the sources in the framework of two-component models considering all available observational constraints from VLBI observations. Moreover, the constraints obtained from the long-term behavior of the sources in the lower energy bands could be used to determine the region where the emission in each band originates. Finally, we attempted to use optical polarisation data to shed new light on the behavior of the two components in the optical band. Our observationally constrained two zone model allows explanation of the entire SED from radio to VHE with two co-located emission regions.

Aims: $\gamma$ rays can be used as a tracer in the search of sources of Galactic cosmic rays (CRs). We present deep observations of the Galactic Centre (GC) region with the MAGIC telescopes, which we use for inferring the underlying CR distribution. Methods: We observed the GC region for ${\approx}100$ hours with the MAGIC telescopes from 2012 to 2017, at high zenith angles (58-70~deg). This implies a larger energy threshold, but also an increased effective collection area compared to low zenith observations. Using new software, we derive instrument response and background models, enabling us to study the diffuse emission in the region. We use pre-existing data of the gas distribution in the GC region to derive the underlying distribution of CRs. Results: We obtain a significant detection for all four model components used to fit our data (Sgr~A*, ``Arc'', G0.9+0.1, and an extended component for the Galactic Ridge). We find that the diffuse component is best described as a power-law with index 2 and an exponential cut-off at around 20~TeV with the significance of the cut-off being only 2~$\sigma$. The derived cosmic-ray profile hints to a peak at the GC position, with a measured profile index of $1.2 \pm 0.3$, supporting the hypothesis of a CR accelerator at the GC. We argue that the measurements of this profile are presently limited by our knowledge of the gas distribution in the GC vicinity.

We present the first results from very-high-energy observations of the dwarf spheroidal satellite candidate Triangulum II with the MAGIC telescopes from 62.4 hours of good-quality data taken between August 2016 and August 2017. We find no gamma-ray excess in the direction of Triangulum II, and upper limits on both the differential and integral gamma-ray flux are presented. Currently, the kinematics of Triangulum II are affected by large uncertainties leading to a bias in the determination of the properties of its dark matter halo. Using a scaling relation between the annihilation J-factor and heliocentric distance of well-known dwarf spheroidal galaxies, we estimate an annihilation J-factor for Triangulum II for WIMP dark matter of $\log[J_{\text{ann}}({0.5^{\circ}})/$ GeV$^{2}$ cm$^{-5}] = 19.35 \pm 0.37$. We also derive a dark matter density profile for the object relying on results from resolved simulations of Milky Way sized dark matter halos. We obtain 95% confidence-level limits on the thermally averaged annihilation cross section for WIMP annihilation into various Standard Model channels. The most stringent limits are obtained in the $\tau^{+}\tau^{-}$ final state, where a cross section for annihilation down to $\langle \sigma_{\text{ann}} v \rangle = 3.05 \times 10^{-24}$ cm$^{3}$ s$^{-1}$ is excluded.

We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $\gamma$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $\gamma$-rays. Data are analyzed from a coordinated campaign across multiple instruments including VHE $\gamma$-ray (VERITAS, MAGIC), high-energy (HE) $\gamma$-ray (Fermi-LAT), X-ray (Swift}, RXTE, MAXI), optical (including the GASP-WEBT collaboration and polarization data) and radio (Mets\"ahovi, OVRO, UMRAO). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare `decline' epochs. The main flare statistics allow 2-minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of $\sim$25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor ($\delta \gtrsim 33$) and the size of the emission region ($ \delta^{-1}R_B \lesssim 3.8\times 10^{13}\,\,\mbox{cm}$) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10-minute-binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship: from linear to quadratic to lack of correlation to anti-correlation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain by the classic single-zone synchrotron self-Compton model.

We report on the first detection of very high-energy (VHE) gamma-ray emission from the Crab Nebula by a Cherenkov telescope in dual-mirror Schwarzschild-Couder (SC) configuration. The result has been achieved by means of the 4 m size ASTRI-Horn telescope, operated on Mt. Etna (Italy) and developed in the context of the Cherenkov Telescope Array Observatory preparatory phase. The dual-mirror SC design is aplanatic and characterized by a small plate scale, allowing us to implement large field of view cameras with small-size pixel sensors and a high compactness. The curved focal plane of the ASTRI camera is covered by silicon photo-multipliers (SiPMs), managed by an unconventional front-end electronics based on a customized peak-sensing detector mode. The system includes internal and external calibration systems, hardware and software for control and acquisition, and the complete data archiving and processing chain. The observations of the Crab Nebula were carried out in December 2018, during the telescope verification phase, for a total observation time (after data selection) of 24.4 h, equally divided into on- and off-axis source exposure. The camera system was still under commissioning and its functionality was not yet completely exploited. Furthermore, due to recent eruptions of the Etna Volcano, the mirror reflection efficiency was reduced. Nevertheless, the observations led to the detection of the source with a statistical significance of 5.4 sigma above an energy threshold of ~3 TeV. This result provides an important step towards the use of dual-mirror systems in Cherenkov gamma-ray astronomy. A pathfinder mini-array based on nine large field-of-view ASTRI-like telescopes is under implementation.

1ES 1959+650 is a bright TeV high-frequency-peaked BL Lac object exhibiting interesting features like "orphan" TeV flares and a broad emission in the high-energy regime, that are difficult to interpret using conventional one-zone Synchrotron Self-Compton (SSC) scenarios. We report the results from the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) observations in 2016 along with the multi-wavelength data from the Fermi Large Area Telescope (LAT) and Swift instruments. MAGIC observed 1ES 1959+650 with different emission levels in the very-high-energy (VHE, E >100 GeV) gamma-ray band during 2016. In the long-term data, the X-ray spectrum becomes harder with increasing flux and a hint of a similar trend is also visible in the VHE band. An exceptionally high VHE flux reaching ~ 3 times the Crab Nebula flux was measured by MAGIC on the 13th, 14th of June and 1st July 2016 (the highest flux observed since 2002). During these flares, the high-energy peak of the spectral energy distribution (SED) lies in the VHE domain and extends up to several TeV. The spectrum in the gamma-ray (both Fermi-LAT and VHE bands) and the X-ray bands are quite hard. On 13th June and 1st July 2016, the source showed rapid variations of the VHE flux within timescales of less than an hour. A simple one-zone SSC model can describe the data during the flares requiring moderate to high values of the Doppler factors (>=30-60). Alternatively, the high-energy peak of the SED can be explained by a purely hadronic model attributed to proton-synchrotron radiation with jet power L_{jet}~10^{46} erg/s and under high values of the magnetic field strength (~100 G) and maximum proton energy (~few EeV). Mixed lepto-hadronic models require super-Eddington values of the jet power. We conclude that it is difficult to get detectable neutrino emission from the source during the extreme VHE flaring period of 2016.

M87 is one of the closest (z=0.00436) extragalactic sources emitting at very-high-energies (VHE, E > 100 GeV). The aim of this work is to locate the region of the VHE gamma-ray emission and to describe the observed broadband spectral energy distribution (SED) during the low VHE gamma-ray state. The data from M87 collected between 2012 and 2015 as part of a MAGIC monitoring programme are analysed and combined with multi-wavelength data from Fermi-LAT, Chandra, HST, EVN, VLBA and the Liverpool Telescope. The averaged VHE gamma-ray spectrum can be fitted from 100GeV to 10TeV with a simple power law with a photon index of (-2.41 $\pm$ 0.07), while the integral flux above 300GeV is $(1.44 \pm 0.13) \times 10^{-12} cm^{-2} s^{-1}$. During the campaign between 2012 and 2015, M87 is generally found in a low emission state at all observed wavelengths. The VHE gamma-ray flux from the present 2012-2015 M87 campaign is consistent with a constant flux with some hint of variability ($\sim3\sigma$) on a daily timescale in 2013. The low-state gamma-ray emission likely originates from the same region as the flare-state emission. Given the broadband SED, both a leptonic synchrotron self Compton and a hybrid photo-hadronic model reproduce the available data well, even if the latter is preferred. We note, however, that the energy stored in the magnetic field in the leptonic scenario is very low suggesting a matter dominated emission region.

Aims: We aim to measure the Crab Nebula gamma-ray spectral energy distribution in the ~100 TeV energy domain and test the validity of existing leptonic emission models at these high energies. Methods: We use the novel very large zenith angle observations with the MAGIC telescope system to increase the collection area above 10 TeV. We also develop an auxiliary procedure of monitoring atmospheric transmission in order to assure proper calibration of the accumulated data. This employs recording of optical images of the stellar field next to the source position, which provides a better than 10% accuracy for the transmission measurements. Results: We demonstrate that MAGIC very large zenith angle observations yield a collection area larger than a square kilometer. In only ~56 hr of observations, we detect the gamma-ray emission from the Crab Nebula up to 100 TeV, thus providing the highest energy measurement of this source to date with Imaging Atmospheric Cherenkov Telescopes. Comparing accumulated and archival MAGIC and Fermi/LAT data with some of the existing emission models, we find that none of them provides an accurate description of the 1 GeV to 100 TeV gamma-ray signal.

We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

Mrk501 is a very high-energy (VHE) gamma-ray blazar located at z=0.034. During a period of two weeks in July 2014, the highest X-ray activity of Mrk501 was observed in ~14 years of operation of the Neil Gehrels Swift Gamma-ray Burst Observatory. We characterize the broadband variability of Mrk501 from radio to VHE gamma rays, and evaluate whether it can be interpreted within theoretical scenarios widely used to explain the broadband emission from blazars. The temporal evolution of the most prominent and variable segments of the SED is described with a one-zone synchrotron self-Compton model with variations in the break energy of the electron energy distribution (EED), and with some adjustments in the magnetic field strength and spectral shape of the EED. A narrow feature at ~3 TeV was observed in the VHE spectrum measured on 2014 July 19 (MJD 56857.98), which is the day with the highest X-ray flux ($>0.3$ keV) measured during the entire Swift mission. This feature is inconsistent with the classical analytic functions to describe the measured VHE spectra (power law, log-parabola, and log-parabola with exponential cutoff) at more than 3$\sigma$. A fit with a log-parabola plus a narrow component is preferred over the fit with a single log-parabola at more than 4$\sigma$, and a dedicated Monte Carlo simulation estimated the significance of this extra component to be larger than 3$\sigma$. Under the assumption that this VHE spectral feature is real, we show that it can be reproduced with three distinct theoretical scenarios: a) a pileup in the EED due to stochastic acceleration; b) a structured jet with two-SSC emitting regions, with one region dominated by an extremely narrow EED; and c) an emission from an IC pair cascade.

Extreme high-energy peaked BL Lac objects (EHBLs) are an emerging class of blazars. Their typical two-hump structured spectral energy distribution (SED) peaks at higher energies with respect to conventional blazars. Multi-wavelength (MWL) observations constrain their synchrotron peak in the medium to hard X-ray band. Their gamma-ray SED peaks above the GeV band, and in some objects it extends up to several TeV. Up to now, only a few EHBLs have been detected in the TeV gamma-ray range. In this paper, we report the detection of the EHBL 2WHSP J073326.7+515354, observed and detected during 2018 in TeV gamma rays with the MAGIC telescopes. The broad-band SED is studied within a MWL context, including an analysis of the Fermi-LAT data over ten years of observation and with simultaneous Swift-XRT, Swift-UVOT, and KVA data. Our analysis results in a set of spectral parameters that confirms the classification of the source as an EHBL. In order to investigate the physical nature of this extreme emission, different theoretical frameworks were tested to model the broad-band SED. The hard TeV spectrum of 2WHSP J073326.7+515354 sets the SED far from the energy equipartition regime in the standard one-zone leptonic scenario of blazar emission. Conversely, more complex models of the jet, represented by either a two-zone spine-layer model or a hadronic emission model, better represent the broad-band SED.

Extreme high-frequency peaked BL Lac objects (EHBLs) are blazars which exhibit extremely energetic synchrotron emission. They also feature non-thermal gamma-ray emission whose peak lies in the very high-energy (VHE, E > 100 GeV) range, and in some sources exceeds 1TeV: this is the case of hard-TeV EHBLs such as 1ES 0229+200. With the aim of increasing the EHBL population, ten targets were observed with the MAGIC telescopes from 2010 to 2017, for a total of 262 h of good quality data. The data were complemented by coordinated Swift observations. The X-ray data analysis confirms that all the sources but two are EHBLs. The sources show only a modest variability and a harder-when-brighter behavior, typical for this class of objects. At VHE gamma rays, three new sources were detected and a hint of signal was found for another new source. In each case the intrinsic spectrum is compatible with the hypothesis of a hard-TeV nature of these EHBLs. The broadband spectral energy distributions (SEDs) of all sources are built and modeled in the framework of a single-zone purely leptonic model. The VHE gamma-ray detected sources were also interpreted with a spine-layer and a proton synchrotron models. The three models provide a good description of the SEDs. However, the resulting parameters differ substantially in the three scenarios, in particular the magnetization parameter. This work presents a first mini-catalog of VHE gamma-ray and multi-wavelength observations of EHBLs.

Starburst galaxies and star-forming active galactic nuclei (AGN) are among the candidate sources thought to contribute appreciably to the extragalactic gamma-ray and neutrino backgrounds. NGC 1068 is the brightest of the star-forming galaxies found to emit gamma rays from 0.1 to 50 GeV. Precise measurements of the high-energy spectrum are crucial to study the particle accelerators and probe the dominant emission mechanisms. We have carried out 125 hours of observations of NGC 1068 with the MAGIC telescopes in order to search for gamma-ray emission in the very high energy band. We did not detect significant gamma-ray emission, and set upper limits at 95\% confidence level to the gamma-ray flux above 200 GeV f<5.1x10^{-13} cm^{-2} s ^{-1} . This limit improves previous constraints by about an order of magnitude and allows us to put tight constraints on the theoretical models for the gamma-ray emission. By combining the MAGIC observations with the Fermi-LAT spectrum we limit the parameter space (spectral slope, maximum energy) of the cosmic ray protons predicted by hadronuclear models for the gamma-ray emission, while we find that a model postulating leptonic emission from a semi-relativistic jet is fully consistent with the limits. We provide predictions for IceCube detection of the neutrino signal foreseen in the hadronic scenario. We predict a maximal IceCube neutrino event rate of 0.07 yr^{-1}.

Several observations are revealing the widespread occurrence of mildly relativistic wide-angle AGN winds strongly interacting with the gas of their host galaxy. Such winds are potential cosmic-ray accelerators, as supported by gamma-ray observations of the nearby Seyfert galaxy NGC 1068 with the Fermi gamma-ray space telescope. The non-thermal emission produced by relativistic particles accelerated by the AGN-driven wind observed in the circum-nuclear molecular disk of such galaxy is invoked to produce the gamma-ray spectrum. The AGN wind model predicts a hard spectrum that extend in the very high energy band which differs significantly from those corresponding to other models discussed in the literature, like starburst or AGN jet. We present dedicated simulations of observations through the Cherenkov Telescope Array (CTA), the next-generation ground based gamma-ray observatory, of the very high energy spectrum of the Seyfert galaxy NGC 1068 assuming the AGN wind and the AGN jet models. We demonstrate that, considering 50 hours of observations, CTA can be effectively used to constrain the two different emission models, providing important insight into the physics governing the acceleration of particles in non-relativistic AGN-driven outflows. This analysis strongly motivates observations of Seyfert and starburst galaxies with CTA in order to test source population models of the extragalactic gamma-ray and neutrino backgrounds.

Recent advances in observations have provided a wealth of measurements of the expansions of outflows in galactic discs out to large radii in a variety of galactic hosts. To provide an updated baseline for the interpretation of such data, and to assess to what extent the present status of the modeling is consistent with the existing observations, we provide a compact two-dimensional description for the expansion of AGN-driven shocks in realistic galactic discs with exponential gas density profiles in a disc geometry. We derive solutions for the outflow expansion and the mass outflow rates in different directions with respect to the plane of the disc. These are expressed in terms of the global properties of the host galaxy and of the central AGN to allow for an easy and direct comparison with existing observations in a variety of galactic hosts with measured properties, and out to distances $\sim 10$ kpc from the centre. The results are compared with a state-of-the-art compilation of observed outflows in 19 galaxies with different measured gas and dynamical mass, allowing for a detailed, one-by-one comparison with the model predictions. The agreement we obtain for a wide range of host galaxy gas mass and AGN bolometric luminosity provides a quantitative systematic test for the modeling of AGN-driven outflows in galactic discs. We also consider a larger sample of galaxies with no reliable measurements of the gas and dynamical mass. In this case we perform a comparison of the model predictions for different bins of AGN luminosities assuming values for the gas mass and dynamical mass derived from scaling relations. The encouraging, quantitative agreement of the model predictions with a wide set of existing observations constitutes a baseline for the interpretation of forthcoming data, and for a more detailed treatment of AGN feedback in galaxy formation models.

We present a measurement of the extragalactic background light (EBL) based on a joint likelihood analysis of 32 gamma-ray spectra for 12 blazars in the redshift range z = 0.03 to 0.944, obtained by the MAGIC telescopes and Fermi-LAT. The EBL is the part of the diffuse extragalactic radiation spanning the ultraviolet, visible and infrared bands. Major contributors to the EBL are the light emitted by stars through the history of the universe, and the fraction of it which was absorbed by dust in galaxies and re-emitted at longer wavelengths. The EBL can be studied indirectly through its effect on very-high energy photons that are emitted by cosmic sources and absorbed via photon-photon interactions during their propagation across cosmological distances. We obtain estimates of the EBL density in good agreement with state-of-the-art models of the EBL production and evolution. The 1-sigma upper bounds, including systematic uncertainties, are between 13% and 23% above the nominal EBL density in the models. No anomaly in the expected transparency of the universe to gamma rays is observed in any range of optical depth.We also perform a wavelength-resolved EBL determination, which results in a hint of an excess of EBL in the 0.18 - 0.62 $\mu$m range relative to the studied models, yet compatible with them within systematics.

A population of globular clusters (GCs) has been recently established by the Fermi-LAT telescope as a new class of GeV $\gamma$-ray sources. Leptons accelerated to TeV energies, in the inner magnetospheres of MSPs or in their wind regions, should produce $\gamma$-rays through the inverse Compton scattering in the dense radiation field from the huge population of stars. We have conducted deep observations of the globular cluster M15 with the MAGIC telescopes and used 165 hrs in order to search for $\gamma$-ray emission. A strong upper limit on the TeV $\gamma$-ray flux $<3.2\times 10^{-13}\mathrm{cm^{-2}s^{-1}}$ above 300 GeV ($<0.26\%$ of the Crab nebula flux) has been obtained. We interpret this limit as a constraint on the efficiency of the acceleration of leptons in the magnetospheres of the MSPs. We constrain the injection rate of relativistic leptons, $\eta_{\rm e}$, from the MSPs magnetospheres and their surrounding. We conclude that $\eta_{\rm e}$ must be lower than expected from the modelling of high energy processes in MSP inner magnetospheres. For leptons accelerated with the power law spectrum in the MSP wind regions, $\eta_{\rm e}$ is constrained to be much lower than derived for the wind regions around classical pulsars. These constraints are valid for the expected range of magnetic field strengths within the GC and for the range of likely energies of leptons injected from the inner magnetospheres, provided that the leptons are not removed from the globular cluster very efficiently due to advection process. We discuss consequences of these constraints for the models of radiation processes around millisecond pulsars.

The mechanisms producing fast variability of the $\gamma$-ray emission in active galactic nuclei are under debate. The MAGIC telescopes detected a fast very high energy (VHE, E$>100$ GeV) $\gamma$-ray flare from BL Lacertae on 2015 June 15. The flare had a maximum flux of $(1.5\pm 0.3)\times 10^{-10}$ photons cm$^{-2}$ s$^{-1}$ and halving time of $26\pm8$ minutes. The MAGIC observations were triggered by a high state in the optical and high energy (HE, E$>100$ MeV) $\gamma$-ray bands. In this paper we present the MAGIC VHE $\gamma$-ray data together with multiwavelength data from radio, optical, X-rays, and HE $\gamma$ rays from 2015 May 1 to July 31. Well-sampled multiwavelength data allow us to study the variability in detail and compare it to the other epochs when fast VHE $\gamma$-ray flares have been detected from this source. Interestingly, we find that the behaviour in radio, optical, X-rays and HE $\gamma$-rays is very similar to two other observed VHE $\gamma$-ray flares. In particular, also during this flare there was an indication of rotation of the optical polarization angle and of activity at the 43\,GHz core. These repeating patterns indicate a connection between the three events. We also test modelling of the spectral energy distribution, based on constraints from the light curves and VLBA observations, with two different geometrical setups of two-zone inverse Compton models. In addition we model the $\gamma$-ray data with the star-jet interaction model. We find that all of the tested emission models are compatible with the fast VHE $\gamma$-ray flare, but all have some tension with the multiwavelength observations.

SNR G24.7+0.6 is a 9.5 kyrs radio and $\gamma$-ray supernova remnant evolving in a dense medium. In the GeV regime, SNR G24.7+0.6 (3FHL\,J1834.1--0706e/FGES\,J1834.1--0706) shows a hard spectral index ($\Gamma$$\sim$2) up to $200$\,GeV, which makes it a good candidate to be observed with Cherenkov telescopes such as MAGIC. We observed the field of view of \snr\ with the MAGIC telescopes for a total of 31 hours. We detect very high energy $\gamma$-ray emission from an extended source located 0.34\degr\ away from the center of the radio SNR. The new source, named \mgc\ is detected up to 5\,TeV, and its spectrum is well-represented by a power-law function with spectral index of $2.74 \pm 0.08$. The complexity of the region makes the identification of the origin of the very-high energy emission difficult, however the spectral agreement with the LAT source and overlapping position at less than 1.5$\sigma$ point to a common origin. We analysed 8 years of \fermi-LAT data to extend the spectrum of the source down to 60\,MeV. \fermi-LAT and MAGIC spectra overlap within errors and the global broad band spectrum is described by a power-law with exponential cutoff at $1.9\pm0.5$\,TeV. The detected $\gamma$-ray emission can be interpreted as the results of proton-proton interaction between the supernova and the CO-rich surrounding.

We report on observations of the pulsar / Be star binary system PSR J2032+4127 / MT91 213 in the energy range between 100 GeV and 20 TeV with the VERITAS and MAGIC imaging atmospheric Cherenkov telescope arrays. The binary orbit has a period of approximately 50 years, with the most recent periastron occurring on 2017 November 13. Our observations span from 18 months prior to periastron to one month after. A new, point-like, gamma-ray source is detected, coincident with the location of PSR J2032+4127 / MT91 213. The gamma-ray light curve and spectrum are well-characterized over the periastron passage. The flux is variable over at least an order of magnitude, peaking at periastron, thus providing a firm association of the TeV source with the pulsar / Be star system. Observations prior to periastron show a cutoff in the spectrum at an energy around 0.5 TeV. This result adds a new member to the small population of known TeV binaries, and it identifies only the second source of this class in which the nature and properties of the compact object are firmly established. We compare the gamma-ray results with the light curve measured with the X-ray Telescope (XRT) on board the Neil Gehrels \textit{Swift} Observatory and with the predictions of recent theoretical models of the system. We conclude that significant revision of the models is required to explain the details of the emission we have observed, and we discuss the relationship between the binary system and the overlapping steady extended source, TeV J2032+4130.

Theoretical models of galaxy formation suggest that the presence of an active galactic nucleus (AGN) is required to regulate the growth of its host galaxy through feedback mechanisms, produced by e.g. AGN-driven outflows. Although such outflows are common both at low and high redshift, a comprehensive picture is still missing. The peak epoch of galaxy assembly (1<z<3) has been poorly explored so far, and current observations in this redshift range are mostly limited to targets with high chances to be in an outflowing phase. This paper introduces SUPER (a SINFONI Survey for Unveiling the Physics and Effect of Radiative feedback), an ongoing ESO's VLT/SINFONI Large Programme. SUPER will perform the first systematic investigation of ionized outflows in a sizeable and blindly-selected sample of 39 X-ray AGN at z~2, which reaches high spatial resolutions (~2 kpc) thanks to the adaptive optics-assisted IFU observations. The outflow morphology and star formation in the host galaxy will be mapped through the broad component of [OIII] and the narrow component of Ha emission lines. The main aim of our survey is to infer the impact of outflows on the on-going star formation and to link the outflow properties to a number of AGN and host galaxy properties. We describe here the survey characteristics and goals, as well as the selection of the target sample. Moreover, we present a full characterization of its multi-wavelength properties: we measure, via spectral energy distribution fitting of UV-to-FIR photometry, stellar masses (4x10^9-2x10^11 Msun), star formation rates (25-680 Msun yr^-1) and AGN bolometric luminosities (2x10^44-8x10^47 erg s^-1), along with obscuring column densities (<2x10^24 cm^-2) and 2-10 keV luminosities (2x10^43-6x10^45 erg s^-1) derived through X-ray spectral analysis. Finally, we classify our AGN as jetted or non-jetted according to their radio and FIR emission.

Fast radio bursts (FRBs) are bright flashes observed typically at GHz frequencies with millisecond duration, whose origin is likely extragalactic. Their nature remains mysterious, motivating searches for counterparts at other wavelengths. FRB 121102 is so far the only source known to repeatedly emit FRBs and is associated with a host galaxy at redshift $z \simeq 0.193$. We conducted simultaneous observations of FRB 121102 with the Arecibo and MAGIC telescopes during several epochs in 2016--2017. This allowed searches for millisecond-timescale burst emission in very-high-energy (VHE) gamma rays as well as the optical band. While a total of five FRBs were detected during these observations, no VHE emission was detected, neither of a persistent nature nor burst-like associated with the FRBs. The average integral flux upper limits above 100 GeV at 95% confidence level are $6.6 \times 10^{-12}~\mathrm{photons\ cm^{-2}\ s^{-1}}$ (corresponding to luminosity $L_{\rm VHE} \lesssim 10^{45}~\mathrm{erg\ s^{-1}}$) over the entire observation period, and $1.2 \times 10^{-7}~ \mathrm{photons\ cm^{-2}\ s^{-1}}$ ($L_{\rm VHE} \lesssim 10^{49}~\mathrm{erg\ s^{-1}}$) over the total duration of the five FRBs. We constrain the optical U-band flux to be below 8.6 mJy at 5-$\sigma$ level for 1-ms intervals around the FRB arrival times. A bright burst with U-band flux $29~\mathrm{mJy}$ and duration $\sim 12$ ms was detected 4.3 s before the arrival of one FRB. However, the probability of spuriously detecting such a signal within the sampled time space is 1.5% (2.2 $\sigma$, post-trial), i.e. consistent with the expected background. We discuss the implications of the obtained upper limits for constraining FRB models.

PKS 1510-089 is a flat spectrum radio quasar strongly variable in the optical and GeV range. We search for low-state VHE gamma-ray emission from PKS 1510-089. We aim to characterize and model the source in a broad-band context, which would provide a baseline over which high states and flares could be better understood. We use daily binned Fermi-LAT flux measurements of PKS 1510-089 to characterize the GeV emission and select the observation periods of MAGIC during low state of activity. For the selected times we compute the average radio, IR, optical, UV, X-ray and gamma-ray emission to construct a low-state spectral energy distribution of the source. The broadband emission is modelled within an External Compton scenario with a stationary emission region through which plasma and magnetic field are flowing. We perform also the emission-model-independent calculations of the maximum absorption in the broad line region (BLR) using two different models. Results. The MAGIC telescopes collected 75 hrs of data during times when the Fermi-LAT flux measured above 1 GeV was below 3x10-8cm-2s-1, which is the threshold adopted for the definition of a low gamma-ray activity state. The data show a strongly significant (9.5{\sigma}) VHE gamma-ray emission at the level of (4.27+-0.61stat)x10-12cm-2s-1 above 150GeV, a factor 80 smaller than the highest flare observed so far from this object. Despite the lower flux, the spectral shape is consistent with earlier detections in the VHE band. The broad-band emission is compatible with the EC scenario assuming a large emission region located beyond the BLR. For the first time the gamma-ray data allow us to place a limit on the location of the emission region during a low gamma-ray state of a FSRQ. For the used model of the BLR, the 95% C.L. on the location of the emission region allows us to place it at the distance >74% of the outer radius of the BLR.

Clusters of galaxies are the largest known gravitationally bound structures in the Universe, with masses around $10^{15}\ M_\odot$, most of it in the form of dark matter. The ground-based Imaging Atmospheric Cherenkov Telescope MAGIC made a deep survey of the Perseus cluster of galaxies using almost 400 h of data recorded between 2009 and 2017. This is the deepest observational campaign so far on a cluster of galaxies in the very high energy range. We search for gamma-ray signals from dark matter particles in the mass range between 200 GeV and 200 TeV decaying into standard model pairs. We apply an analysis optimized for the spectral and morphological features expected from dark matter decays and find no evidence of decaying dark matter. From this, we conclude that dark matter particles have a decay lifetime longer than $\sim10^{26}$~s in all considered channels. Our results improve previous lower limits found by MAGIC and represent the strongest limits on decaying dark matter particles from ground-based gamma-ray instruments.

A neutrino with energy of $\sim$290 TeV, IceCube-170922A, was detected in coincidence with the BL Lac object TXS~0506+056 during enhanced gamma-ray activity, with chance coincidence being rejected at $\sim 3\sigma$ level. We monitored the object in the very-high-energy (VHE) band with the MAGIC telescopes for $\sim$41 hours from 1.3 to 40.4 days after the neutrino detection. Day-timescale variability is clearly resolved. We interpret the quasi-simultaneous neutrino and broadband electromagnetic observations with a novel one-zone lepto-hadronic model, based on interactions of electrons and protons co-accelerated in the jet with external photons originating from a slow-moving plasma sheath surrounding the faster jet spine. We can reproduce the multiwavelength spectra of TXS 0506+056 with neutrino rate and energy compatible with IceCube-170922A, and with plausible values for the jet power of $\sim 10^{45} - 4 \times 10^{46} {\rm erg \ s^{-1}}$. The steep spectrum observed by MAGIC is concordant with internal $\gamma\gamma$ absorption above a few tens of GeV entailed by photohadronic production of a $\sim$290 TeV neutrino, corroborating a genuine connection between the multi-messenger signals. In contrast to previous predictions of predominantly hadronic emission from neutrino sources, the gamma-rays can be mostly ascribed to inverse Compton up-scattering of external photons by accelerated electrons. The X-ray and VHE bands provide crucial constraints on the emission from both accelerated electrons and protons. We infer that the maximum energy of protons in the jet co-moving frame can be in the range $\sim 10^{14}$ to $10^{18}$ eV.

The MAGIC telescopes observed S2 0109+22 in 2015 July during its flaring activity in high energy gamma rays observed by Fermi-LAT. We analyse the MAGIC data to characterise the very high energy (VHE) gamma-ray emission of S2 0109+22, which belongs to the subclass of intermediate synchrotron peak (ISP) BL Lac objects. We study the multi-frequency emission in order to investigate the source classification. Finally, we compare the source long-term behaviour to other VHE gamma-ray emitting (TeV) blazars. We performed a temporal and spectral analysis of the data centred around the MAGIC interval of observation (MJD 57225-57231). Long-term radio and optical data have also been investigated using the discrete correlation function. The redshift of the source is estimated through optical host-galaxy imaging and also using the amount of VHE gamma-ray absorption. The quasi-simultaneous multi-frequency spectral energy distribution (SED) is modelled with the conventional one-zone synchrotron self-Compton (SSC) model. MAGIC observations resulted in the detection of the source at a significance level of $5.3\,\sigma$. The VHE gamma-ray emission of S2 0109+22 is variable on a daily time scale. VHE gamma-ray luminosity of the source is lower than the average of TeV BL Lacs. The optical polarization, and long-term optical/radio behaviour of the source are different from the general population of TeV blazars. All these findings agree with the classification of the source as an ISP BL Lac object. We estimate the source redshift as $z = 0.36 \pm 0.07$. The SSC parameters describing the SED are rather typical for blazars.

We report on the detection of flaring activity from the Fanaroff-Riley I radio galaxy NGC 1275 in very-high-energy (VHE, E $>$ 100 GeV) gamma rays with the MAGIC telescopes. Observations were performed between 2016 September and 2017 February as part of a monitoring program. The brightest outburst with $\sim1.5$ times the Crab Nebula flux above 100 GeV (C.U.) was observed during the night between 2016 December 31 and 2017 January 1 (fifty times higher than the mean previously measured in two observational campaigns between 2009 and 2011). Significant variability of the day-by-day light curve was measured, the shortest flux-doubling time-scales was found to be of $(611\pm101)$ min. The combined spectrum of the MAGIC data during the strongest flare state and simultaneous data from the Fermi-LAT around 2017 January 1 follows a power-law with an exponential cutoff at the energy $(492\pm35)$ GeV. Simultaneous optical flux density measurements in the R-band obtained with the KVA telescope are also presented and the correlation between the optical and gamma-ray emission is investigated. Due to possible internal pair-production, the fast flux variability constrains the Doppler factor to values which are inconsistent with a large viewing angle as observed in the radio band. We investigate different scenarios for the explanation of fast gamma-ray variability, namely emission from: magnetospheric gaps, relativistic blobs propagating in the jet (mini-jets) or external cloud (or star) entering the jet. We find that the only plausible model to account for the luminosities here observed would be the production of gamma rays in a magnetospheric gap around the central black hole only in the eventuality of an enhancement of the magnetic field threading the hole from its equipartition value with the gas pressure in the accretion flow.

We investigate galaxy formation in models with dark matter (DM) constituted by sterile neutrinos. Given their large parameter space, defined by the combinations of sterile neutrino mass $m_{\nu}$ and mixing parameter $\sin^2(2\theta)$ with active neutrinos, we focus on models with $m_{\nu}=7$ keV, consistent with the tentative 3.5 keV line detected in several X-ray spectra of clusters and galaxies. We consider i) two resonant production models with $\sin^2(2\theta)=5\,10^{-11}$ and $\sin^2(2\theta)=2\,10^{-10}$, to cover the range of mixing parameter consistent with the 3.5 keV line; ii) two scalar-decay models, representative of the two possible cases characterizing such a scenario: a freeze-in and a freeze-out case. We also consider thermal Warm Dark Matter with particle mass $m_X=3$ keV. Using a semi-analytic model, we compare the predictions for the different DM scenarios with a wide set of observables. We find that comparing the predicted evolution of the stellar mass function, the abundance of satellites of Milky Way-like galaxies, and the global star formation history of galaxies with observations does not allow to disentangle the effects of the baryonic physics from those related to the different DM models. On the other hand, the distribution of the stellar-to-halo mass ratios, the abundance of faint galaxies in the UV luminosity function at $z\gtrsim 6$, and the specific star formation and age distribution of local, low-mass galaxies constitute potential probes for the considered DM scenarios. We discuss how next observations with upcoming facilities will enable to rule out or to strongly support DM models based on sterile neutrinos.

We investigate the relation between the star formation rate (SFR) and the stellar mass (M), i.e. the Main Sequence (MS) relation of star-forming galaxies, at 1.3<=z<6 in the first four HST Frontier Fields, based on rest-frame UV observations. Gravitational lensing combined with deep HST observations allows us to extend the analysis of the MS down to stellar masses as low as logM/Msun~7.5 at z<~4 and logM/Msun~8 at higher redshifts, a factor of ~10 below most previous results. We perform an accurate simulation to take into account the effect of observational uncertainties on the MS and correct for the Eddington bias. This step allows us to reliably measure the MS and in particular its slope. While the normalization increases with redshift, we fit an unevolving and approximately linear slope. We nicely extend to lower masses the results of brighter surveys. Thanks to the large dynamic range in mass and by making use of the simulation, we analyzed any possible mass dependence of the dispersion around the MS. We find tentative evidence that the scatter decreases with increasing stellar masses, suggesting a larger variety of star formation histories in low mass galaxies. This trend agrees with theoretical predictions, and is explained as either a consequence of the smaller number of progenitors of low mass galaxies in a hierarchical scenario and/or of the efficient but intermittent stellar feedback processes in low mass halos. Finally, we observe an increase in the SFR per unit stellar mass with redshift milder than predicted by theoretical models, implying a still incomplete understanding of the processes responsible for galaxy growth.

We derive the contribution to the extragalactic gamma-ray background (EGB) from AGN winds and star-forming galaxies by including a physical model for the gamma-ray emission produced by relativistic protons accelerated by AGN-driven and supernova-driven shocks into a state-of-the-art semi-analytic model of galaxy formation. This is based on galaxy interactions as triggers of AGN accretion and starburst activity and on expanding blast wave as the mechanism to communicate outwards the energy injected into the interstellar medium by the active nucleus. We compare the model predictions with the latest measurement of the EGB spectrum performed by the Fermi-LAT in the range between 100 MeV and 820 GeV. We find that AGN winds can provide ~35$\pm$15% of the observed EGB in the energy interval E_{\gamma}=0.1-1 GeV, for ~73$\pm$15% at E_{\gamma}=1-10 GeV, and for ~60$\pm$20% at E_{\gamma}>10 GeV. The AGN wind contribution to the EGB is predicted to be larger by a factor of 3-5 than that provided by star-forming galaxies (quiescent plus starburst) in the hierarchical clustering scenario. The cumulative gamma-ray emission from AGN winds and blazars can account for the amplitude and spectral shape of the EGB, assuming the standard acceleration theory, and AGN wind parameters that agree with observations. We also compare the model prediction for the cumulative neutrino background from AGN winds with the most recent IceCube data. We find that for AGN winds with accelerated proton spectral index p=2.2-2.3, and taking into account internal absorption of gamma-rays, the Fermi-LAT and IceCube data could be reproduced simultaneously.

Feedback from accreting SMBHs is often identified as the main mechanism responsible for regulating star-formation in AGN host galaxies. However, the relationships between AGN activity, radiation, winds, and star-formation are complex and still far from being understood. We study scaling relations between AGN properties, host galaxy properties and AGN winds. We then evaluate the wind mean impact on the global star-formation history, taking into account the short AGN duty cycle with respect to that of star-formation. We first collect AGN wind observations for 94 AGN with detected massive winds at sub-pc to kpc spatial scales. We then fold AGN wind scaling relations with AGN luminosity functions, to evaluate the average AGN wind mass-loading factor as a function of cosmic time. We find strong correlations between the AGN molecular and ionised wind mass outflow rates and the AGN bolometric luminosity. The power law scaling is steeper for ionised winds (slope 1.29+/-0.38) than for molecular winds (0.76+/-0.06), meaning that the two rates converge at high bolometric luminosities. The molecular gas depletion timescale and the molecular gas fraction of galaxies hosting powerful AGN winds are 3-10 times shorter and smaller than those of main-sequence galaxies with similar SFR, stellar mass and redshift. These findings suggest that, at high AGN bolometric luminosity, the reduced molecular gas fraction may be due to the destruction of molecules by the wind, leading to a larger fraction of gas in the atomic ionised phase. The AGN wind mass-loading factor $\eta=\dot M_{OF}/SFR$ is systematically higher than that of starburst driven winds. Our analysis shows that AGN winds are, on average, powerful enough to clean galaxies from their molecular gas only in massive systems at z<=2, i.e. a strong form of co-evolution between SMBHs and galaxies appears to break down for the least massive galaxies.

We compute the non-thermal emissions produced by relativistic particles accelerated by the AGN-driven shocks in NGC 1068, and we compare the model predictions with the observed gamma-ray and radio spectra . The former is contributed by pion decay, inverse Compton scattering, and bremsstrahlung, while the latter is produced by synchrotron radiation. We derive the gamma-ray and radio emissions by assuming the standard acceleration theory, and we discuss how our results compare with those corresponding to other commonly assumed sources of gamma-ray and radio emissions, like Supernova remnants (SNR) or AGN jets. We find that the AGN-driven shocks observed in the circumnuclear molecular disk of such a galaxy provide a contribution to the gamma-ray emission comparable to that provided by the starburst activity when standard particle acceleration efficiencies are assumed, while they can yield the whole gamma-ray emission only when the parameters describing the acceleration efficiency and the proton coupling with the molecular gas are tuned to values larger than those assumed in standard, SNR-driven shocks. We discuss the range of acceleration efficiencies (for protons and electrons) and of proton calorimetric fractions required to account for the observed gamma-ray emission in the AGN outflow model. We further compare the neutrino flux expected in our model with constraints from current experiments, and we provide predictions for the detections by the upcoming KM3NeT neutrino telescope. This analysis strongly motivates observations of NGC 1068 at >TeV energies with current and future Cherenkov telescopes in order to gain insight into the nature of the gamma-rays source.

Recent observations indicate that the mass of Supermassive Black Holes (SMBHs) correlate differently with different galaxy stellar components. Comparing such observations with the results of "ab initio" galaxy formation models can provide insight on the mechanisms leading to the growth of SMBHs. Here we use a state-of-the-art semi-analytic model of galaxy formation to investigate the correlation of the different galaxy stellar components with the mass of the central SMBH. The stellar mass in the disc, in the bulge, and in the pseudo-bulge of galaxies is related to quiescent star formation, to galaxy interactions, and to the loss of angular momentum following disc instabilities, respectively. Consistently with recent findings, we find that while the predicted bulge masses are tightly correlated with the SMBH masses, the correlation between the latter and the galactic discs shows a much larger scatter, in particular when bulgeless galaxies are considered. In addition, we obtain that the predicted masses of pseudo-bulges shows little or no-correlation with the masses of SMBHs. We track the histories of merging, star formation, and SMBH accretion to investigate the physical processes at the origin of such findings within the context of cosmological models of galaxy formation. Finally, we discuss the effects of variations of our assumed fiducial model on the results.

AGN driven outflows are invoked in numerical simulations to reproduce several observed properties of local galaxies. The z > 1 epoch is of particular interest as it was during this time that the volume averaged star formation and the accretion rate of black holes were maximum. Radiatively driven outflows are therefore believed to be common during this epoch. We aim to trace and characterize outflows in AGN hosts with high mass accretion rates at z > 1 using integral field spectroscopy. We obtain spatially-resolved kinematics of the [OIII]5007 line in two targets which reveal the morphology and spatial extension of the outflows. We present J and H+K band SINFONI observations of 5 AGNs at 1.2 < z < 2.2. To maximize the chance of observing radiatively driven outflows, our sample was pre-selected based on peculiar values of the Eddington ratio and the hydrogen column density of the surrounding interstellar medium. We observe high velocity (~600-1900 km/s) and kiloparsec scale extended ionized outflows in at least 3 of our targets, using [OIII]5007 line kinematics tracing the AGN narrow line region. We estimate the total mass of the outflow, the mass outflow rate, and the kinetic power of the outflows based on theoretical models and report on the uncertainties associated with them. We find mass outflow rates of ~1-10 M_sun/yr for the sample presented in this paper. Based on the high star formation rates of the host galaxies, the observed outflow kinetic power and the expected power due to the AGN, we infer that both star formation and AGN radiation could be the dominant source for the outflows. The outflow models suffer from large uncertainties, hence we call for further detailed observations for an accurate determination of the outflow properties to confirm the exact source of these outflows.

The triggering mechanisms for Active Galactic Nuclei (AGN) are still debated. Some of the most popular ones include galaxy interactions (IT) and disk instabilities (DI). Using an advanced semi analytic model (SAM) of galaxy formation, coupled to accurate halo occupation distribution modeling, we investigate the imprint left by each separate triggering process on the clustering strength of AGN at small and large scales. Our main results are as follows: i) DIs, irrespective of their exact implementation in the SAM, tend to fall short in triggering AGN activity in galaxies at the center of halos with $M_h>10^{13.5} h^{-1}M_{\odot}$. On the contrary, the IT scenario predicts abundance of active, central galaxies that generally agrees well with observations at every halo mass. ii) The relative number of satellite AGN in DIs at intermediate-to-low luminosities is always significantly higher than in IT models, especially in groups and clusters. The low AGN satellite fraction predicted for the IT scenario might suggest that different feeding modes could simultaneously contribute to the triggering of satellite AGN. iii) Both scenarios are quite degenerate in matching large-scale clustering measurements, suggesting that the sole average bias might not be an effective observational constraint. iv) Our analysis suggests the presence of both a mild luminosity and a more consistent redshift dependence in the AGN clustering, with AGN inhabiting progressively less massive dark matter halos as the redshift increases. We also discuss the impact of different observational selection cuts in measuring AGN clustering, including possible discrepancies between optical and X-ray surveys.

Using an advanced semi analytic model (SAM) for galaxy formation, we have investigated the statistical effects of assuming two different mechanisms for triggering AGN activity on the properties of AGN host galaxies. We have considered a first accretion mode where AGN activity is triggered by disk instabilities (DI) in isolated galaxies, and a second feeding mode where such an activity is triggered by galaxy mergers and fly-by events (interactions, IT). We obtained the following results:i) for hosts with $M_* \lesssim 10^{11} M_{\bigodot}$, both DI and IT modes are able to account for the observed AGN hosts stellar mass function; for more massive hosts, the DI scenario predicts a lower space density than the IT model, lying below the observational estimates for z>0.8.ii) The analysis of the color-magnitude diagram (CMD) of AGN hosts for redshift z < 1.5 can provide a good observational test to effectively discriminate between the DI and IT mode, since DIs are expected to yield AGN host galaxy colors skewed towards bluer colors, while in the IT scenario the majority of hosts are expected to reside in the red sequence.iii) While both IT and DI scenarios can account for AGN triggered in main sequence or starburst galaxies, DIs fail in triggering AGN activity in passive galaxies.iv) The two modes are characterized by a different duration of the AGN phase, with DIs lasting even on time scales $\sim $ Gyr, much longer with respect to the IT scenario.v) The scatter of the $SFR-L_{bol}$ relation could represent another crucial diagnostics to discriminate between the two triggering modes, since the DI scenario predicts an appreciably lower scatter of the relation than the IT scenario. vi) Disk instabilities are not able to account for the observed fraction of AGN in groups for z < 1 and clusters for z < 0.7, while the IT scenario provides a good match to observational data.

The process that quenched star formation in galaxies at intermediate and high redshift is still the subject of considerable debate. One way to investigate this puzzling issue is to study the number density of quiescent galaxies at z~2, and its dependence on mass. Here we present the results of a new study based on very deep Ks-band imaging (with the HAWK-I instrument on the VLT) of two HST CANDELS fields (the UKIDSS Ultra-deep survey (UDS) field and GOODS-South). The new HAWK-I data (taken as part of the HUGS VLT Large Program) reach detection limits of Ks>26 (AB mag). We select a sample of passively-evolving galaxies in the redshift range 1.4<z<2.5. Thanks to the depth and large area coverage of our imaging, we have been able to extend the selection of quiescent galaxies a magnitude fainter than previous analyses. Through extensive simulations we demonstrate, for the first time, that the observed turn-over in the number of quiescent galaxies at K>22 is real. This has enabled us to establish unambiguously that the number counts of quiescent galaxies at z~2 flatten and slightly decline at magnitudes fainter than Ks~22(AB mag.). We show that this trend corresponds to a stellar mass threshold $M_*10^{10.8}\,{\rm M_{\odot}}$ below which the mechanism that halts the star formation in high-redshift galaxies seems to be inefficient. Finally we compare the observed pBzK number counts with those of quiescent galaxies extracted from four different semi-analytic models. We find that none of the models provides a statistically acceptable description of the number density of quiescent galaxies at these redshifts. We conclude that the mass function of quiescent galaxies as a function of redshift continues to present a key and demanding challenge for proposed models of galaxy formation and evolution.

Using a semi analytic model for galaxy formation we investigate the effects of Black Hole accretion triggered by disk instabilities (DI) in isolated galaxies on the evolution of AGN. Specifically, we took on, developed and expanded the Hopkins & Quataert (2011) model for the mass inflow following disk perturbations, and compare the corresponding evolution of the AGN population with that arising in a scenario where galaxy interactions trigger AGN (IT mode). We extended and developed the DI model by including different disk surface density profiles, to study the maximal contribution of DI to the evolution of the AGN population. We obtained the following results: i) for luminosities corresponding to $M_{1450}\gtrsim -26$ the DI mode can provide the BH accretion needed to match the observed AGN luminosity functions up to $z \approx 4.5$; in such a luminosity range and redshift, it can compete with the IT scenario as the main driver of cosmological evolution of AGN; ii) The DI scenario cannot provide the observed abundance of high-luminosity QSO with $M_{1450}\lesssim -26$ AGN, as well as the abundance of high-redhshift $z \approx 4.5$ QSOs with $M_{1450}\lesssim -24$, while the IT scenario provides an acceptable match up to $z \approx 6$, as found in our earliest works; iii) The dispersion of the distributions of Eddington ratio for low- and intermediate-luminosity AGN (bolometric $L_{AGN}$ = $10^{43}$ - $10^{45}$ erg/s) is predicted to be much smaller in the DI scenario compared to the IT mode; iv) The above conclusions are robust with respect to the explored variants of the Hopkins & Quataert (2011) model. We discuss the physical origin of our findings, and how it is possible to pin down the dominant fueling mechanism in the low-intermediate luminosity range $M_{1450}\gtrsim -26$ where both the DI and the IT modes are viable candidates as drivers for the AGN evolution.

We present near-infrared spectra, obtained with SINFONI and XShooter observations at ESO VLT, of nine dusty, red QSOs at 1.2<z<2.6. The sources are hard X-ray detected, characterized by cold absorption (N_H>10^{21} - 10^{22} cm^{-2}) and show a broad Ha component in the NIR spectra. We complement this sample with twelve additional sources taken from the literature with similar properties resulting in a total sample of 21 X-ray obscured, intermediate type (1.8-1.9), dusty reddened QSOs. From the broad Ha line we have computed the BH masses through the virial formula and derived Eddington ratios. Moreover, from optical/IR multi-component SED fitting we have derived the stellar mass of their host galaxies and their SFRs. We find that most of the sources in our sample are hosted in starburst and main sequence star-forming galaxies with Eddington ratios lambda>0.1. We find a strong trend with the BH mass i.e. less massive objects are scattered below and above the local relation while the most massive ones are mainly located above it. We also studied the evolution of these sources on the M_BH-M_star plane compared to a sample of optically blue type--1 QSOs and we find that obscured red QSOs show a ratio of M_BH to M_star that increases with redshift which is consistent with or slightly lower than what has been found for blue QSOs. These sources may represent the blow-out phase at the end of the rapid BH growth and immediately preceding the classical blue QSOs typically sampled in optical surveys. They in fact show evidence of outflows in the ionized gas component, but their BH has already fully formed.

We investigate various galaxy population properties of the massive X-ray luminous galaxy cluster XDCP J0044.0-2033 at z=1.58, which constitutes the most extreme matter density peak at this redshift currently known. We analyze deep VLT/HAWK-I NIR data in the J- and Ks-bands, complemented by Subaru imaging in i and V, Spitzer observations at 4.5 micron, and new spectroscopic observations with VLT/FORS2. We detect a cluster-associated excess population of about 90 galaxies, which follows a centrally peaked, compact NFW galaxy surface density profile with a concentration of c200~10. Based on the Spitzer 4.5 micron imaging data, we measure a stellar mass fraction of fstar,500=(3.3+-1.4)% consistent with local values. The total J- and Ks-band galaxy luminosity functions of the core region yield characteristic magnitudes J* and Ks* consistent with expectations from simple z_f=3 burst models. However, a detailed look at the morphologies and color distributions of the spectroscopically confirmed members reveals that the most massive galaxies are undergoing a very active mass assembly epoch through merging processes. Consequently, the bright end of the cluster red-sequence is not in place, while at intermediate magnitudes [Ks*,Ks*+1.6] a red-locus population is present, which is then sharply truncated at magnitudes fainter than Ks*+1.6. The dominant cluster core population comprises post-quenched galaxies transitioning towards the red-sequence at intermediate magnitudes, while additionally a significant blue cloud population of faint star-forming galaxies is present even in the densest central regions. Our observations lend support to the scenario in which the dominant effect of the dense z~1.6 cluster environment is an accelerated mass assembly timescale through merging activity that is responsible for driving core galaxies across the mass quenching threshold of log(Mstar/Msun)~10.4.

We use deep Herschel PACS and SPIRE observations in GOODSS, GOODSN and COSMOS to estimate the average dust mass (Mdust) of galaxies on a redshift-stellar mass (Mstar)-SFR grid. We study the scaling relations between Mdust, Mstar and SFR at z<=2.5. No clear evolution of Mdust is observed at fixed SFR and Mstar. We find a tight correlation between SFR and Mdust, likely a consequence of the Schmidt-Kennicutt (S-K) law. The Mstar-Mdust correlation observed by previous works flattens or sometimes disappears when fixing the SFR. Most of it likely derives from the combination of the Mdust-SFR and Mstar-SFR correlations. We then investigate the gas content as inferred by converting Mdust by assuming that the dust/gas ratio scales linearly with the gas metallicity. All galaxies in the sample follow, within uncertainties, the same SFR-Mgas relation (integrated S-K law), which broadly agrees with CO-based results for the bulk of the population, despite the completely different approaches. The majority of galaxies at z~2 form stars with an efficiency (SFE=SFR/Mgas) ~5 times higher than at z~0. It is not clear what fraction of such variation is an intrinsic redshift evolution and what fraction arises from selection effects. The gas fraction (fgas) decreases with Mstar and increases with SFR, and does not evolve with z at fixed Mstar and SFR. We explain these trends by introducing a universal relation between fgas, Mstar and SFR, non-evolving out to z~2.5. Galaxies move across this relation as their gas content evolves in time. We use the 3D fundamental fgas-Mstar-SFR relation and the redshift evolution of the Main Sequence to estimate the evolution of fgas in the average population of galaxies as a function of z and Mstar, and we find evidence a downsizing scenario.

We estimate the fraction of AGNs hosted in starburst galaxies (f_bursty) as a function of the AGN luminosity predicted under the assumption that starburst events and AGN activity are triggered by galaxy interactions during their merging histories. The latter are described through Monte Carlo realizations, and are connected to star formation and BH accretion using a semi-analytic model of galaxy formation in a cosmological framework. The predicted fraction f_bursty increases steeply with AGN luminosity from <0.2 at L_X < 10^44 erg/s to >0.9 at L_X > 10^45 erg/s over a wide redshift interval from z=0 to z=6. We compare the model predictions with new measurements of f_bursty from a sample of X-ray selected AGNs in the XMM-COSMOS field at 0.3< z< 2, and from a sample of QSOs (L_X > 10^45 erg/s) in the redshift range 2< z< 6.5. We find preliminary indications that under conservative assumptions half of the QSO hosts are starburst galaxies. This result provide motivation for future systematic studies of the stellar properties of high luminosity AGN hosts in order to constrain AGN triggering mechanisms.

Recent measurements of the abundance of AGN with low-luminosities (L_X< 10^44 erg/s in the 2-10 keV energy band) at high redshifts z>4 provide a serious challenge for Cold Dark Matter (CDM) models based on interaction-driven fueling of AGN. Using a semi-analytic model of galaxy formation we investigate how such observations fit in a Warm Dark Matter (WDM) scenario of galaxy formation, and compare the results with those obtained in the standard CDM scenario with different efficiencies for the stellar feedback. Taking on our previous exploration of galaxy formation in WDM cosmology, we assume as a reference case a spectrum which is suppressed - compared to the standard CDM case - below a cut-off scale ~ 0.2$ Mpc corresponding (for thermal relic WDM particles) to a mass m_X=0.75 keV. We run our fiducial semi-analytic model with such a WDM spectrum to derive AGN luminosity functions from z~6 to the present over a wide range of luminosities (10^43< L_X/erg/s < 10^46 in the 2-10 keV X-ray band), to compare with recent observations and with the results in the CDM case. When compared with the standard CDM case, the luminosity distributions we obtain assuming a WDM spectrum are characterized by a similar behaviour at low redshift, and by a flatter slope at faint magnitudes for z>3, which provide an excellent fit to present observations. We discuss how such a result compares with CDM models with maximized feedback efficiency, and how future deep AGN surveys will allow for a better discrimination between feedback and cosmological effects on the evolution of AGN in interaction-driven models for AGN fueling.

An increasing amount of observational evidence supports the notion that there are two modes of star formation: a quiescent mode in disk-like galaxies, and a starburst mode, which is generally interpreted as driven by merging. Using a semi-analytic model of galaxy formation, we derive the relative contribution to the cosmic star formation rate density of quiescently starforming and starburst galaxies, predicted under the assumption that starburst events are triggered by galaxy encounters (merging and fly-by kind) during their merging histories. We show that, within this framework, quiescently starforming galaxies dominate the cosmic star formation rate density at all redshifts. The contribution of the burst-dominated starforming galaxies increases with redshift, rising from <5% at low redshift (z<0.1) to ~20% at z>5. We estimated that the fraction of the final (z=0) galaxy stellar mass which is formed through the burst component of star formation is ~10% for 10^10 M_\odot<M_*<10^11.5 M_\odot. Starburst galaxies, selected according to their distance from the galaxy main sequence, account for ~10% of the star formation rate density in the redshift interval 1.5<z<2.5, i.e. at the cosmic peak of the star formation activity.

We investigate for the first time the effects of a Warm Dark Matter (WDM) power spectrum on the statistical properties of galaxies using a semi-analytic model of galaxy formation. The WDM spectrum we adopt as a reference case is suppressed - compared to the standard Cold Dark Matter (CDM) case - below a cut-off scale ~ 1 Mpc corresponding (for thermal relic WDM particles) to a mass m_X=0.75 keV. This ensures consistency with present bounds provided by the microwave background WMAP data and by the comparison of hydrodynamical N-body simulations with observed Lyman-{\alpha} forest. We run our fiducial semi-analytic model with such a WDM spectrum to derive galaxy luminosity functions (in B, UV, and K bands) and the stellar mass distributions over a wide range of cosmic epochs, to compare with recent observations and with the results in the CDM case. The predicted color distribution of galaxies in the WDM model is also checked against the data. When compared with the standard CDM case, the luminosity and stellar mass distributions we obtain assuming a WDM spectrum are characterized by: i) a flattening of the faint end slope and ii) a sharpening of the cutoff at the bright end for z \lesssim 0.8. We discuss how the former result is directly related to the smaller number of low-mass haloes collapsing in the WDM scenario, while the latter is related to the smaller number of satellite galaxies accumulating in massive haloes at low redshift, thus suppressing the accretion of small lumps on the central, massive galaxies. These results shows how a adopting a WDM power spectrum may contribute to solve two major problems of CDM galaxy formation scenarios, namely, the excess of predicted faint (low mass) galaxies at low and - most of all - high redshifts, and the excess of bright (massive) galaxies at low redshifts.

We compute the number density of massive Black Holes (BHs) at the centre of galaxies at z=6 in different Dynamical Dark Energy (DDE) cosmologies, and compare it with existing observational lower limits, to derive constraints on the evolution of the Dark Energy equation of state parameter w. Our approach only assumes the canonical scenario for structure formation from the collapse of overdense regions of the Dark Matter dominated primordial density field on progressively larger scales; the Black Hole accretion and merging rate have been maximized in the computation so as to obtain robust constraints on w and on its look-back time derivative w_a. Our results provide independent constraints complementary to those obtained by combining Supernovae, Cosmic Microwave Background and Baryonic Acoustic Oscillations; while the latter concern combinations of w_0 and w_a leaving the time evolution of the state parameter w_a highly unconstrained, the BH abundance mainly provide upper limits on w_a, only weakly depending on w_0. Combined with the existing constraints, our results significantly restrict the allowed region in DDE parameter space, ruling out DDE models not providing cosmic time and fast growth factor large enough to allow for the building up of the observed abundance of BHs; in particular, models with -1.2 \leq w_0 \leq -1 and positive redshift evolution w_a > 0.8 - completely consistent with previous constraints - are strongly disfavoured by our independent constraints from BH abundance. Such range of parameters corresponds to "Quintom" DDE models, with w crossing -1 starting from larger values.

We present detection and analysis of faint X-ray sources in the Chandra deep field south (CDFS) using the 4 Msec Chandra observation and adopting a new detection algorithm, based on a targeted search at the position of known high-z galaxies. This optimized technique results in the identification of 54 z>3 AGNs, 29 of which are new detections. Applying stringent completeness criteria, we derive AGN luminosity functions in the redshift bins 3-4, 4-5 and >5.8 and for 42.75<log L(2-10 keV)<44.5. We join this data with the luminous AGN luminosity functions from optical surveys and find that the evolution of the high-z, wide luminosity range luminosity function can be best modeled by pure luminosity evolution with L* decreasing from 6x10^44 ergs/s at z=3 to L*=2x10^44 ergs/s at z=6. We compare the high-z luminosity function with the prediction of theoretical models using galaxy interactions as AGN triggering mechanism. We find that these models are broadly able to reproduce the high-z AGN luminosity functions. A better agreement is found assuming a minimum Dark Matter halo mass for black hole formation and growth. We compare our AGN luminosity functions with galaxy mass functions to derive high-z AGN duty cycle using observed Eddington ratio distributions to derive black hole masses. We find that the duty cycle increases with galaxy stellar mass and redshift by a factor 10-30 from z=0.25 to z=4-5. We also report on the detection of a large fraction of highly obscured, Compton thick AGN at z>3 (18+17-10%). Their optical counterparts are not strongly reddened and we thus conclude that the size of the X-ray absorber is likely smaller than the dust sublimation radius. We finally report the discovery of a highly star-forming galaxy at z=3.47. If confirmed, this would be one of the farthest objects in which stellar sources are detected in X-rays.

In the standard scenario for galaxy evolution young star-forming galaxies transform into red bulge-dominated spheroids, where star formation has been quenched. To explain such a transformation, a strong negative feedback generated by accretion onto a super-massive black hole is often invoked. The depletion of gas resulting from quasar-driven outflows should eventually stop star-formation in the host galaxy and lead the black hole to "suicide" for starvation. Direct observational evidence for a major quasar feedback onto the host galaxy is still missing, since outflows previously observed in quasars are generally associated with the ionized component of the gas, which only accounts for a minor fraction of the total gas content, and typically occurring in the central regions. We used the IRAM PdB Interferometer to observe CO(1-0) in Mrk 231, the closest QSO known. Thanks to the wide band we detect broad wings of the CO line, with velocities up to 750 km/s and spatially resolved on the kpc scale. Such broad CO wings trace a giant molecular outflow of about 700 MSun/yr, far larger than the ongoing star-formation rate (~200 MSun/yr) observed in the host galaxy. This wind will totally expel the cold gas reservoir in Mrk 231 in about 1e7 yrs, therefore halting the star-formation activity on the same timescale. The inferred kinetic energy in the molecular outflow is ~1.2e44 erg/s, corresponding to a few percent of the AGN bolometric luminosity, very close to the fraction expected by models ascribing quasar feedback to highly supersonic shocks generated by radiatively accelerated nuclear winds. The direct observational evidence for quasar feedback reported here provides solid support to the scenarios ascribing the observed properties of local massive galaxies to quasar-induced large scale winds.

We derive the growth of SMBHs relative to the stellar content of their host galaxy predicted under the assumption of BH accretion triggered by galaxy encounters occurring during their merging histories. We show that, within this framework, the ratio Gamma=(M_BH/M_*)(z)/(M_BH/M_*)(z=0) between the Black Hole mass and the galactic stellar mass (normalized to the local value) depends on both BH mass and redshift. While the average value and the spread of Gamma(z) increase with z, such an effect is larger for massive BHs, reaching values Gamma=5 for massive Black Holes (M>10^9 M_{\odot}) at z>4, in agreement with recent observations of high-redshift QSOs; this is due to the effectiveness of interactions in triggering BH accretion in high-density environments at high redshifts. To test such a model against observations, we worked out specific predictions for sub-samples of the simulated galaxies corresponding to the different observational samples for which measurements of Gamma have been obtained. We found that for Broad Line AGNs at 1<z<2 values of Gamma=2 are expected, with a mild trend toward larger value for increasing BH mass. Instead, when we select from our Monte Carlo simulations only extremely gas rich, rapidly star forming galaxies at 2<z<3, we find low values 0.3<Gamma<1.5, consistent with recent observational findings on samples of sub-mm galaxies; in the framework of our model, these objects end up at z=0 in low-to-intermediate mass BHs (M<10^9 M_{\odot}), and they do not represent typical paths leading to local massive galaxies. The latter have formed preferentially through paths passing above the local M_*-M_BH relation. We discuss how the global picture emerging from the model is consistent with a downsizing scenario, where massive BHs accrete a larger fraction of their final mass at high redshifts z>4.

There have been recent claims that a significant fraction of type 2 AGN accrete close or even above the Eddington limit. In type 2 AGN the bolometric luminosity (L_b) is generally inferred from the [OIII] emission line luminosity (L_OIII). The key issue, in order to estimate the bolometric luminosity in these AGN, is therefore to know the bolometric correction to be applied to L_OIII. A complication arises from the fact that the observed L_OIII is affected by extinction, likely due to dust within the narrow line region. The extinction-corrected [OIII] luminosity (L^c_OIII) is a better estimator of the nuclear luminosity than L_OIII. However, so far only the bolometric correction to be applied to the uncorrected L_OIII has been evaluated. This paper is devoted to estimate the bolometric correction C_OIII=L_b/L^c_OIII in order to derive the Eddington ratios for the type 2 AGN in a sample of SDSS objects. We have collected from the literature 61 sources with reliable estimate of both L^c_OIII and X-ray luminosities (L_X). To estimate C_OIII, we combined the observed correlation between L^c_OIII and L_X with the X-ray bolometric correction. We found, contrary to previous studies, a linear correlation between L^c_OIII and L_X. We estimated C_OIII using the luminosity-dependent X-ray bolometric correction of Marconi et al. (2004), and we found a mean value of C_OIII in the luminosity ranges log L_OIII=38-40, 40-42, and 42-44 of 87, 142 and 454 respectively. We used it to calculate the Eddington ratio distribution of type 2 SDSS AGN at 0.3<z<0.4 and we found that these sources are not accreting near their Eddington limit, contrary to previous claims.

The evidence of a decrease with increasing luminosity of the fraction f_{abs} of absorbed and Compton-thin among X-ray (2-10 keV) selected AGN is observationally rather well supported, while that of an increase of f_{abs} with redshift is rather controversial. In Lamastra, Perola & Matt (2006) the gravitational effect of the SMBH on the molecular interstellar gas, in the central region of the host galaxy, was shown to predict an anti-correlation between f_{abs} and M_{BH}. The most recent findings on the distribution of the Eddington ratio \lambda=L_b/L_E as a function of M_{BH} and z are used to convert that relationship into one between f_{abs} and both bolometric (L_b) and X-ray (L_X) luminosities at various values of z. The findings on \lambda(M_{BH},z) are properly treated in order to ensure completeness in the prediction of f_{abs} above a certain luminosity, at values of z=0.1, 0.35, 0.7 and >1. To verify the consequence of these findings alone, we adopted in a first istance a distribution of gas surface density \Sigma, observed in a sample of local spiral galaxies, irrespective of the galaxy morphological type and z. Assuming in the \lambda(M_{BH},z) distribution the Eddington limit, \lambda=1, as a ``natural'' cut-off, the predictions are consistent with the existence of an anti-correlation between f_{abs} and L_X, but fail to reproduce an increase of f_{abs} with z. Because the early type galaxies on average are much poorer in molecular gas than late type ones, a quantitative agreement with the local value of f_{abs} requires the existence of a correlation between \Sigma and the central activity. An increase of typical values of \Sigma with z, correlated with the activity, might explain an increase of f_{abs} with z. However, at the highest luminosities f_{abs} could hardly exceed about 0.3.

Context: The counts of galaxies and AGN in the mid infra-red (MIR) bands are important instruments for studying their cosmological evolution. However, the classic spectral line ratios techniques can become misleading when trying to properly separate AGN from starbursts or even from apparently normal galaxies. Aims: We use X-ray band observations to discriminate AGN activity in previously classified MIR-selected starburst galaxies and to derive updated AGN1 and (Compton thin) AGN2 counts at 15 um. Methods: XMM observations of the ELAIS-S1 15um sample down to flux limits ~2x10^-15 erg cm^-2 s^-1 (2-10 keV band) were used. We classified as AGN all those MIR sources with a unabsorbed 2-10 keV X-ray luminosity higher that ~10^42 erg/s. Results: We find that at least about 13(+/-6) per cent of the previously classified starburst galaxies harbor an AGN. According to these figures, we provide an updated estimate of the counts of AGN1 and (Compton thin) AGN2 at 15 um. It turns out that at least 24% of the extragalactic sources brighter than 0.6 my at 15 um are AGN (~13% contribution to the extragalactic background produced at fluxes brighter than 0.6 mJy).

The BH mass (and the related Eddington ratio) in broad line AGN is usually evaluated by combining estimates (often indirect) of the BLR radius and of the FWHM of the broad lines, under the assumption that the BLR clouds are in Keplerian motion around the BH. Such an evaluation depends on the geometry of the BLR. There are two major options for the BLR configuration: spherically symmetric or ``flattened''. In the latter case the inclination to the line of sight becomes a relevant parameter. This paper is devoted to evaluate the bias on the estimate of the Eddington ratio when a spherical geometry is assumed (more generally when inclination effects are ignored), while the actual configuration is ``flattened'', as some evidence suggests. This is done as a function of luminosity and redshift, on the basis of recent results which show the existence of a correlation between the fraction of obscured AGN and these two parameters up to at least z=2.5. The assumed BLR velocity field is akin to the ``generalized thick disk'' proposed by Collin et al. (2006). Assuming an isotropic orientation in the sky, the mean value of the bias is calculated as a function of luminosity and redshift. It is demonstrated that, on average, the Eddington ratio obtained assuming a spherical geometry is underestimated for high luminosities, and overestimated for low luminosities. This bias converges for all luminosities at z about 2.7, while nothing can be said on this bias at larger redshifts due to the lack of data. The effects of the bias, averaged over the luminosity function of broad line AGN, have been calculated. The results imply that the bias associated with the a-sphericity of the BLR make even worse the discrepancy between the observations and the predictions of evolutionary models.

The fraction of AGN with photoelectric absorption in the X-rays ranging from NH of 10^{22} up to about 10^{24} cm^{-2} (Compton-thin) appears observationally to be anticorrelated to their luminosity Lx. This recently found evidence is used to investigate the location of the absorbing gas. The molecular torus invoked in the unified picture of AGN, while it can be regarded as confirmed on several grounds to explain the Compton-thick objects, do not conform to this new constraint, at least in its physical models as developed so far. In the frame of observationally based evidence that in Compton-thin sources the absorbing gas might be located far away from the X-ray source, it is shown that the gravitational effects of the black hole (BH) on the molecular gas in a disk, within 25-450 pc (depending on the BH mass, from 10^6 to 10^9 M_solar, leads naturally to the observed anticorrelation, under the assumption of a statistical correlation between the BH mass and Lx. Its normalization is also reproduced provided that the surface density, Sigma, of this gas is larger than about 150-200 M_solar pc^{-2}, and assuming that the bolometric luminosity is one tenth of the Eddington limit. Interestingly, the required values are consistent with the value of the 300 pc molecular disk in our own galaxy, namely 500 M_solar pc^{-2}. In a sample of nearby galaxies from the BIMA SONG survey, it is found that half of the objects have central Sigma larger than 150 M_solar pc${-2}. Given the simplicity of the proposed model, this finding is very encouraging, waiting for future higher resolution surveys in CO on more distant galaxies.

We have determined the cosmological evolution of the density of AGN and of their Nh distribution as a function of the un-absorbed 2-10 keV luminosity up to redshift 4. We used the HELLAS2XMM sample combined with other published catalogs, yielding a total of 508 AGN. Our best fit is obtained with a luminosity-dependent density evolution (LDDE) model where low luminosity (Lx~10^43 erg s^-1) AGN peak at z~0.7, while high luminosity AGN (Lx>10^45 erg s^-1) peak at z~2.0. A pure luminosity evolution model (PLE) can instead be rejected. There is evidence that the fraction of absorbed (Nh>10^22 cm^-2) AGN decreases with the intrinsic X-ray luminosity, and increases with the redshift. Our best fit solution provides a good fit to the observed counts, the cosmic X-ray background, and to the observed fraction of absorbed AGN as a function of the flux in the 10^-15<S(2-10)<10^-10 erg s^-1 cm^-2 range. We find that the absorbed, high luminosity (Lx>10^44 erg s^-1) AGN have a density of 267 deg^-2 at fluxes S(2-10)>10^-15 erg s^-1 cm^-2. Using these results, we estimate a density of supermassive black holes in the local Universe of 3.2 (h^2_70) x 10^5 M_sol Mpc^-3, which is consistent with the recent measurements of the black hole mass function in the local galaxies.