The uncertainty region of the highly energetic neutrino IceCube200107A includes 3HSP J095507.9+355101 ($z$~=~0.557), an extreme blazar, which was detected in a high, very hard, and variable X-ray state shortly after the neutrino arrival. Following a detailed multi-wavelength investigation, we confirm that the source is a genuine BL Lac, contrary to TXS 0506+056, the first source so far associated with IceCube neutrinos, which is a "masquerading" BL Lac. As in the case of TXS0506+056, 3HSP J095507.9+355101 is also way off the so-called "blazar sequence". We consider 3HSP J095507.9+355101 a possible counterpart to the IceCube neutrino. Finally, we discuss some theoretical implications in terms of neutrino production.

We have analysed all the X-ray images centred on Gamma Ray Bursts generated by Swift over the last 15 years using automatic tools that do not require any expertise in X-ray astronomy, producing results in excellent agreement with previous findings. This work, besides presenting the largest medium-deep survey of the X-ray sky and a complete sample of blazars, wishes to be a step in the direction of achieving the ultimate goal of the Open Universe Initiative, that is to enable non expert people to fully benefit of space science data, possibly extending the potential for scientific discovery, currently confined within a small number of highly specialised teams, to a much larger population. We have used the Swift_deepsky Docker container encapsulated pipeline to build the largest existing flux-limited and unbiased sample of serendipitous X-ray sources. Swift_deepsky runs on any laptop or desktop computer with a modern operating system. The tool automatically downloads the data and the calibration files from the archives, runs the official Swift analysis software and produces a number of results including images, the list of detected sources, X-ray fluxes, SED data, and spectral slope estimations. We used our source list to build the LogN-LogS of extra-galactic sources, which perfectly matches that estimated by other satellites. Combining our survey with multi-frequency data we selected a complete radio flux-density limited sample of High Energy Peaked (HBL) blazars.

The association of two IceCube detections, the IceCube-170922A event and a neutrino flare, with the blazar TXS0506+056, has paved the way for the multimessenger quest for cosmic accelerators. IceCube has observed many other neutrinos but their origin remains unknown. To better understand the reason for the apparent lack of neutrino counterparts we have extended the comprehensive dissection of the sky area performed for the IceCube-170922A event to all the 70 public IceCube high-energy neutrinos that are well reconstructed and off the Galactic plane. Using the multi-frequency data available through the Open Universe platform, we have identified numerous candidate counterparts of IceCube events. We report here the classification of all the gamma-ray blazars found and the results of subsequent statistical tests. In addition, we have checked the 4LAC, 3FHL and 3HSP catalogues for potential counterparts. Following the dissection of all areas associated with IceCube neutrinos, we evaluate the data using a likelihood-ratio test and find a 3.23 sigma (post-trial) excess of HBLs and IBLs with a best-fit of 15 +/- 3.6 signal sources. This result, together with previous findings, consistently points to a growing evidence for a connection between IceCube neutrinos and blazars, the most energetic particle accelerators known in the Universe.

The Pacific Ocean Neutrino Experiment (P-ONE) is a new initiative with a vision towards constructing a multi-cubic kilometre neutrino telescope, to expand our observable window of the Universe to highest energies, installed within the deep Pacific Ocean underwater infrastructure of Ocean Networks Canada.

3HSP J095507.9+355101 is an extreme blazar to have been possibly associated with a high-energy neutrino, IceCube-200107A, which was detected one day before the blazar was observed and found to undergo a hard X-ray flare. Motivated by this observation, we perform a comprehensive study of the predicted multimessenger emission from 3HSP J095507.9+355101 during its recent X-ray flare, but also in the long term. We focus on single-zone leptohadronic models, where the broadband photon and high-energy neutrino emissions originate from same region of the jet, but we also explore alternative scenarios: (i) a blazar-core model, which considers neutrino production in the inner jet close to the accreting supermassive black hole, (ii) a hidden external-photon model, which considers neutrino production in the jet through interactions with photons from a possible weak broad line region, (iii) a one-zone proton synchrotron model, where high-energy protons produce $\gamma$-rays in the jet via synchrotron, and (iv) an intergalactic cascade scenario, where neutrinos are produced in the intergalactic medium by interactions of a high-energy cosmic-ray beam escaping the blazar jet. We find that the Poisson probability to detect one or more neutrinos in 8 years from 3HSP J095507.9+355101 is $\sim 2\%$ ($3\%$) when considering the most optimistic one-zone leptohadronic model (the multi-zone blazar-core model), while detection of one neutrino event during the X-ray flare is much less likely. If the association is real, then IceCube-Gen2 and other future detectors should be able to provide additional evidence for neutrino production in 3HSP J095507.9+355101 and other extreme blazars, assuming similar physical conditions.

The BL Lac object 4FGLJ0955.1+3551 has been suggested as a possible source of ultra energetic neutrinos detected by the Icecube observatory. The target was observed in January 2020 at the Large Binocular Telescope. Our spectroscopy (4100-8500 Ang) yields a firm redshift z = 0.557 as deduced by the absorption lines of the host galaxy. The upper limit of the minimum equivalent width on emission lines is 0.3 Ang. From the source image we are able to resolve the host galaxy for which we measure an absolute magnitude M(R) = -22.9 and Re = 8 kpc, that is values which are typical of the host galaxies of BL Lacs.

The most powerful persistent accelerators in the Universe are jetted active galaxies. Galaxies whose jets are directed towards Earth, so called blazars, dominate the extragalactic gamma-ray sky. Still, most of the highest-energy accelerators likely elude detection. These extreme blazars, whose radiated energy can peak beyond 10 TeV, are ideal targets to study particle acceleration and radiative processes, and may provide links to cosmic rays and astrophysical neutrinos. The growing number of extreme blazars observed at TeV energies has been critical for the emergence of gamma-ray cosmology, including measurements of the extragalactic background light, tight bounds on the intergalactic magnetic field, and constraints on exotic physics at energies inaccessible with human-made accelerators. Tremendous progress is expected in the decade to come, particularly with the deployment of the Cherenkov Telescope Array.

I summarise the concluding remarks I gave at the Multifrequency Behaviour of High Energy Cosmic Sources - XIII Workshop. That was not a summary talk and was meant to be provocative. I first give what I think the main message of the workshop was, then provide some (biased) highlights, touch upon the upcoming new facilities and the issues of "quantity vs. quality" and productivity in astronomy, and finally conclude with a look to the future. Astronomers who did not attend the workshop might still find the first two topics appealing and the last two thought-provoking.

High synchrotron peaked blazars (HSPs or HBLs) play a central role in very high energy (VHE) gamma-ray astronomy, and likely in neutrino astronomy. Currently, the largest compilation of HSP blazars, the 2WHSP sample, includes 1691 sources, but it is not complete neither in the radio nor in the X-ray band. In order to provide a more sizable and more accurate set of HSP blazars that is useful for future statistical studies and to plan for VHE/TeV observations, we present the largest sample of extreme and high synchrotron peaked (EHSP;HSP) blazars and blazar candidates: the 3HSP catalogue. The 3HSP catalogue includes 2013 sources, 88% of which with a redshift estimation, a much larger percentage than in any other list of HSP blazars. All new {\gamma}-ray detections are described in the First and Second Brazil ICRANet gamma-ray blazar catalogues (1BIGB & 2BIGB) also taking into account the 4FGL list of {\gamma}-ray sources published by the Fermi large area telescope (Fermi-LAT) team. Moreover, the cross-matching between the 2WHSP, 2FHL HSPs, and IceCube neutrinos position has suggested that HSPs are likely counterparts of neutrino events, which implies the 3HSP catalogue is useful also in that respect. The 3HSP catalogue features improved completeness compared to its predecessors, the 1WHSP and 2WHSP catalogues, and follows the track of their increasing relevance for VHE astronomy.

Motivated by the recently reported evidence of an association between a high-energy neutrino and a gamma-ray flare from the blazar TXS 0506+056, we calculate the expected high-energy neutrino signal from past, individual flares, from twelve blazars, selected in declinations favourable for detection with IceCube. To keep the number of free parameters to a minimum, we mainly focus on BL Lac objects and assume the synchrotron self-Compton mechanism produces the bulk of the high-energy emission. We consider a broad range of the allowed parameter space for the efficiency of proton acceleration, the proton content of BL Lac jets, and the presence of external photon fields. To model the expected neutrino fluence we use simultaneous multi-wavelength observations. We find that in the absence of external photon fields and with jet proton luminosity normalised to match the observed production rate of ultra-high-energy cosmic rays, individual flaring sources produce a modest neutrino flux in IceCube, $\lesssim10^{-3}$ muon neutrinos with energy exceeding 100 TeV, stacking ten years of flare periods selected in the >800 MeV Fermi energy range from each source. Under optimistic assumptions about the jet proton luminosity and in the presence of external photon fields, we find that the two most powerful sources in our sample, AO 0235+164, and OJ 287, would produce, in total, $\approx 3$ muon neutrinos during ten years of Fermi flaring periods, in future neutrino detectors with total instrumented volume $\sim$ten times larger than IceCube,or otherwise, constrain the proton luminosity of blazar jets.

The central nuclei of galaxies, where super-massive black holes (SMBHs) are thought to reside, can experience phases of activity when they become Active Galactic Nuclei (AGN). An AGN can eject winds, jets, and produce radiation across the entire electromagnetic spectrum. The fraction of the bolometric emission in the radio spans a factor of ~10^5 across the different AGN classes. The weakest radio sources, radio-quiet (RQ) AGN, are typically 1,000 times fainter than the radio-loud (RL) AGN, and represent the majority of the AGN population. In RL AGN, radio emission is essentially all produced by synchrotron emission from a relativistic jet. In contrast, in RQ AGN the absence of luminous jets allows us to probe radio emission from a wide range of possible mechanisms, from the host galaxy kpc scale down to the innermost region near the SMBHs: star formation, AGN driven wind, free-free emission from photo-ionized gas, low power jet, and the innermost accretion disc coronal activity. All these mechanisms can now be probed with unprecedented precision and spatial resolution, thanks to the current and forthcoming generation of highly sensitive radio arrays.

We present evidence that TXS 0506+056, the first plausible non-stellar neutrino source, despite appearances, is not a blazar of the BL Lac type but is instead a masquerading BL Lac, i.e., intrinsically a flat-spectrum radio quasar with hidden broad lines and a standard accretion disk. This re-classification is based on: (1) its radio and O II luminosities; (2) its emission line ratios; (3) its Eddington ratio. We also point out that the synchrotron peak frequency of TXS 0506+056 is more than two orders of magnitude larger than expected by the so-called "blazar sequence", a scenario which has been assumed by some theoretical models predicting neutrino (and cosmic-ray) emission from blazars. Finally, we comment on the theoretical implications this re-classification has on the location of the $\gamma$-ray emitting region and our understanding of neutrino emission in blazars.

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.

The Advanced Telescope for High Energy Astrophysics (Athena) is the X-ray observatory large mission selected by the European Space Agency (ESA), within its Cosmic Vision 2015-2025 programme, to address the "Hot and Energetic Universe" scientific theme (Nandra et al. 2013), and it is provisionally due for launch in the early 2030s. The Square Kilometer Array (SKA) is the next generation radio observatory and consists of two telescopes, one comprised of dishes operating at mid frequencies (SKA1-MID) and located in South Africa, and the other comprised of Log-Periodic antennas operating at low radio frequencies (SKA1-LOW), which will be located in Australia (Braun et al. 2017). The scientific commissioning of the radio telescope is planned to begin in 2021-2022. The SKA-Athena Synergy Team (SAST) has been tasked to single out the potential scientific synergies between Athena and SKA. The astrophysical community was involved in this exercise primarily through a dedicated SKA-Athena Synergy Workshop, which took place on April 24-25, 2017 at SKAO, Jodrell Bank, Manchester. The final result of the synergy exercise, this White Paper, describes in detail a number of scientific opportunities that will be opened up by the combination of Athena and SKA, these include: 1. the Cosmic Dawn; 2. the Evolution of black holes and galaxies; 3. Active galaxy feedback in galaxy clusters; 4. Non-thermal phenomena in galaxy clusters; 5. Detecting the cosmic web; 6. Black-hole accretion physics and astrophysical transients; 7. Galactic astronomy: stars, planets, pulsars and supernovae.

We present the dissection in space, time, and energy of the region around the IceCube-170922A neutrino alert. This study is motivated by: (1) the first association between a neutrino alert and a blazar in a flaring state, TXS 0506+056; (2) the evidence of a neutrino flaring activity during 2014 - 2015 from the same direction; (3) the lack of an accompanying simultaneous $\gamma$-ray enhancement from the same counterpart; (4) the contrasting flaring activity of a neighbouring bright $\gamma$-ray source, the blazar PKS 0502+049, during 2014 - 2015. Our study makes use of multi-wavelength archival data accessed through Open Universe tools and includes a new analysis of Fermi-LAT data. We find that PKS 0502+049 contaminates the $\gamma$-ray emission region at low energies but TXS 0506+056 dominates the sky above a few GeV. TXS 0506+056, which is a very strong (top percent) radio and $\gamma$-ray source, is in a high $\gamma$-ray state during the neutrino alert but in a low though hard $\gamma$-ray state in coincidence with the neutrino flare. Both states can be reconciled with the energy associated with the neutrino emission and, in particular during the low/hard state, there is evidence that TXS 0506+056 has undergone a hadronic flare with very important implications for blazar modelling. All multi-messenger diagnostics reported here support a single coherent picture in which TXS 0506+056, a very high energy $\gamma$-ray blazar, is the only counterpart of all the neutrino emissions in the region and therefore the most plausible first non-stellar neutrino and, hence, cosmic ray source.

We test the recently proposed idea that outflows associated with Active Galactic Nuclei (AGN) could be neutrino emitters in two complementary ways. First, we cross-correlate a list of 94 "bona fide" AGN outflows with the most complete and updated repository of IceCube neutrinos currently publicly available, assembled by us for this purpose. It turns out that AGN with outflows matched to an IceCube neutrino have outflow and kinetic energy rates, and bolometric powers larger than those of AGN with outflows not matched to neutrinos. Second, we carry out a statistical analysis on a catalogue of O III 5007 line profiles using a sample of 23,264 AGN at z < 0.4, a sub-sample of which includes mostly possible outflows sources. We find no significant evidence of an association between the AGN and the IceCube events, although we get the smallest p-values (~ 6 and 18 per cent respectively, pre-trial) for relatively high velocities and luminosities. Our results are consistent with a scenario where AGN outflows are neutrino emitters but at present do not provide a significant signal. This can be tested with better statistics and source stacking. A predominant role of AGN outflows in explaining the IceCube data appears in any case to be ruled out.

Active Galactic Nuclei (AGN) are traditionally divided empirically into two main classes: "radio-loud" and "radio-quiet" sources. These labels, which are more than fifty years old, are obsolete, misleading, and wrong. I argue that AGN should be classified based on a fundamentally physical rather than just an observational difference, namely the presence (or lack) of strong relativistic jets, and that we should use the terms "jetted" and "non-jetted" AGN instead.

Active Galactic Nuclei (AGN) are energetic astrophysical sources powered by accretion onto supermassive black holes in galaxies, and present unique observational signatures that cover the full electromagnetic spectrum over more than twenty orders of magnitude in frequency. The rich phenomenology of AGN has resulted in a large number of different "flavours" in the literature that now comprise a complex and confusing AGN "zoo". It is increasingly clear that these classifications are only partially related to intrinsic differences between AGN, and primarily reflect variations in a relatively small number of astrophysical parameters as well the method by which each class of AGN is selected. Taken together, observations in different electromagnetic bands as well as variations over time provide complementary windows on the physics of different sub-structures in the AGN. In this review, we present an overview of AGN multi-wavelength properties with the aim of painting their "big picture" through observations in each electromagnetic band from radio to gamma-rays as well as AGN variability. We address what we can learn from each observational method, the impact of selection effects, the physics behind the emission at each wavelength, and the potential for future studies. To conclude we use these observations to piece together the basic architecture of AGN, discuss our current understanding of unification models, and highlight some open questions that present opportunities for future observational and theoretical progress.

The Advanced Telescope for High ENergy Astrophysics (Athena) is the X-ray observatory mission selected by ESA within its Cosmic Vision 2015-2025 programme to address the Hot and Energetic Universe scientific theme. The ESO-Athena Synergy Team (EAST) has been tasked to single out the potential scientific synergies between Athena and optical/near-infrared (NIR) and sub/mm ground based facilities, in particular those of ESO (i.e., the VLT and ELT, ALMA and APEX), by producing a White Paper to identify and develop the: 1. needs to access ESO ground-based facilities to achieve the formulated Athena science objectives; 2. needs to access Athena to achieve the formulated science objectives of ESO facilities contemporary to Athena; 3. science areas where the synergetic use of Athena and ESO facilities in the late 2020s will result in scientific added value. Community input to the process happened primarily via a dedicated ESO - Athena Synergy Workshop that took place on Sept. 14 - 16, 2016 at ESO, Garching. This White Paper presents the results of the EAST's work, sorted by synergy area, and deals with the following topics: 1. the Hot Universe: Early groups and clusters and their evolution, Physics of the Intracluster medium, Missing baryons in cosmic filaments; 2. the Energetic Universe: Supermassive black hole (SMBH) history, SMBH accretion disks, Active Galactic Nuclei feedback - Molecular outflows, Ultra-fast outflows, Accretion Physics, Transient Science; 3. Observatory Science: Star Formation, Stars. It then discusses the optical-NIR-sub-mm perspective by providing details on VLT/MOONS, the E-ELT instruments, in particular the MOS, VISTA/4MOST, the ESO and ALMA archives, future ALMA and ESO developments, and finally the (likely) ESO - Athena astronomical scene in the 2020s. (abridged)

Similarly to the cosmic star formation history, the black hole accretion rate density of the Universe peaked at 1<z<3. This cosmic epoch is hence best suited for investigating the effects of radiative feedback from AGN. Observational efforts are underway to quantify the impact of AGN feedback, if any, on their host galaxies. Here we present a study of the molecular gas content of AGN hosts at z~1.5 using CO[2-1] line emission observed with ALMA for a sample of 10 AGNs. We compare this with a sample of galaxies without an AGN matched in redshift, stellar mass, and star formation rate. We detect CO in 3 AGNs with $\mathrm{L_{CO} \sim 6.3-25.1\times 10^{9} L_{\odot}}$ which translates to a molecular hydrogen gas mass of $\mathrm{2.5-10\times 10^{10} M_{\odot}}$ assuming conventional conversion factor of $\mathrm{\alpha_{CO}}\sim3.6$. Our results indicate a >99% probability of lower depletion time scales and lower molecular gas fractions in AGN hosts with respect to the non-AGN comparison sample. We discuss the implications of these observations on the impact that AGN feedback may have on star formation efficiency of z>1 galaxies.

We present a strong hint of a connection between high energy $\gamma$-ray emitting blazars, very high energy neutrinos, and ultra high energy cosmic rays. We first identify potential hadronic sources by filtering $\gamma$-ray emitters %from existing catalogs that are in spatial coincidence with the high energy neutrinos detected by IceCube. The neutrino filtered $\gamma$-ray emitters are then correlated with the ultra high energy cosmic rays from the Pierre Auger Observatory and the Telescope Array by scanning in $\gamma$-ray flux ($F_{\gamma}$) and angular separation ($\theta$) between sources and cosmic rays. A maximal excess of 80 cosmic rays (42.5 expected) is found at $\theta\leq10^{\circ}$ from the neutrino filtered $\gamma$-ray emitters selected from the second hard {\it Fermi}-LAT catalogue (2FHL) and for $F_\gamma\left(>50\:\mathrm{GeV}\right)\geq1.8\times10^{-11}\:\mathrm{ph}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}$. The probability for this to happen is $2.4 \times 10^{-5}$, which translates to $\sim 2.4 \times 10^{-3}$ after compensation for all the considered trials. No excess of cosmic rays is instead observed for the complement sample of $\gamma$-ray emitters (i.e. not in spatial connection with IceCube neutrinos). A likelihood ratio test comparing the connection between the neutrino filtered and the complement source samples with the cosmic rays favours a connection between neutrino filtered emitters and cosmic rays with a probability of $\sim1.8\times10^{-3}$ ($2.9\sigma)$ after compensation for all the considered trials. The neutrino filtered $\gamma$-ray sources that make up the cosmic rays excess are blazars of the high synchrotron peak type. More statistics is needed to further investigate these sources as candidate cosmic ray and neutrino emitters.

Aims. High Synchrotron Peaked blazars (HSPs) dominate the -ray sky at energies larger than a few GeV; however, only a few hundred blazars of this type have been catalogued so far. In this paper we present the 2WHSP sample, the largest and most complete list of HSP blazars available to date, which is an expansion of the 1WHSP catalog of gamma-ray source candidates off the Galactic plane. Methods. We cross-matched a number of multi-wavelength surveys (in the radio, infrared and X-ray bands) and applied selection criteria based on the radio to IR and IR to X-ray spectral slopes. To ensure the selection of genuine HSPs we examined the SED of each candidate and estimated the peak frequency of its synchrotron emission ($\nu_{\rm peak}$) using the ASDC SED tool, including only sources with Synchrotron peak $\nu_{\rm peak} > 10^{15}$ Hz (equivalent to $\nu_{\rm peak} > 4$ eV). Results. We have assembled the largest and most complete catalog of HSP blazars to date, which includes 1691 sources. A number of population properties, such as infrared colours, synchrotron peak, redshift distributions, and -ray spectral properties, have been used to characterise the sample and maximize completeness. We also derived the radio logN-logS distribution. This catalog has already been used to provide seeds to discover new very high energy objects within Fermi-LAT data and to look for the counterparts of neutrino and ultra high energy cosmic ray sources, showing its potential for the identification of promising high-energy gamma-ray sources and multi-messenger targets.

Radio astronomy has changed. For years it studied relatively rare sources, which emit mostly non-thermal radiation across the entire electromagnetic spectrum, i.e. radio quasars and radio galaxies. Now it is reaching such faint flux densities that it detects mainly star-forming galaxies and the more common radio-quiet active galactic nuclei. These sources make up the bulk of the extragalactic sky, which has been studied for decades in the infrared, optical, and X-ray bands. I follow the transformation of radio astronomy by reviewing the main components of the radio sky at the bright and faint ends, the issue of their proper classification, their number counts, luminosity functions, and evolution. The overall "big picture" astrophysical implications of these results, and their relevance for a number of hot topics in extragalactic astronomy, are also discussed. The future prospects of the faint radio sky are very bright, as we will soon be flooded with survey data. This review should be useful to all extragalactic astronomers, irrespective of their favourite electromagnetic band(s), and even stellar astronomers might find it somewhat gratifying.

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.

We explore the correlation of $\gamma$-ray emitting blazars with IceCube neutrinos by using three very recently completed, and independently built, catalogues and the latest neutrino lists. We introduce a new observable, namely the number of neutrino events with at least one $\gamma$-ray counterpart, $N_{\nu}$. In all three catalogues we consistently observe a positive fluctuation of $N_{\nu}$ with respect to the mean random expectation at a significance level of $0.4 - 1.3$ per cent. This applies only to extreme blazars, namely strong, very high energy $\gamma$-ray sources of the high energy peaked type, and implies a model-independent fraction of the current IceCube signal $\sim 10 - 20$ per cent. An investigation of the hybrid photon -- neutrino spectral energy distributions of the most likely candidates reveals a set of $\approx 5$ such sources, which could be linked to the corresponding IceCube neutrinos. Other types of blazars, when testable, give null correlation results. Although we could not perform a similar correlation study for Galactic sources, we have also identified two (further) strong Galactic $\gamma$-ray sources as most probable counterparts of IceCube neutrinos through their hybrid spectral energy distributions. We have reasons to believe that our blazar results are not constrained by the $\gamma$-ray samples but by the neutrino statistics, which means that the detection of more astrophysical neutrinos could turn this first hint into a discovery.

We present two ways of modeling the spectral energy distribution of blazars in the hadronic context and discuss the predictions of each "hadronic variant" on the spectral shape, the multi-wavelength variability, the cosmic-ray flux, and the high-energy neutrino emission. Focusing on the latter, we then present an application of the hadronic model to individual BL Lacs that were recently suggested to be the counterparts of some of the IceCube neutrinos, and conclude by discussing the contribution of the whole BL Lac class to the observed neutrino background.

We investigate the star formation properties of ~800 sources detected in one of the deepest radio surveys at 1.4 GHz. Our sample spans a wide redshift range (~0.1 - 4) and about four orders of magnitude in star formation rate (SFR). It includes both star forming galaxies (SFGs) and active galactic nuclei (AGNs), further divided into radio-quiet and radio-loud objects. We compare the SFR derived from the far infrared luminosity, as traced by Herschel, with the SFR computed from their radio emission. We find that the radio power is a good SFR tracer not only for pure SFGs but also in the host galaxies of RQ AGNs, with no significant deviation with redshift or specific SFR. Moreover, we quantify the contribution of the starburst activity in the SFGs population and the occurrence of AGNs in sources with different level of star formation. Finally we discuss the possibility of using deep radio survey as a tool to study the cosmic star formation history.

Blazars have been suggested as possible neutrino sources long before the recent IceCube discovery of high-energy neutrinos. We re-examine this possibility within a new framework built upon the blazar simplified view and a self-consistent modelling of neutrino emission from individual sources. The former is a recently proposed paradigm that explains the diverse statistical properties of blazars adopting minimal assumptions on blazars' physical and geometrical properties. This view, tested through detailed Monte Carlo simulations, reproduces the main features of radio, X-ray, and gamma-ray blazar surveys and also the extragalactic gamma-ray background at energies > 10 GeV. Here we add a hadronic component for neutrino production and estimate the neutrino emission from BL Lacs as a class, "calibrated" by fitting the spectral energy distributions of a preselected sample of BL Lac objects and their (putative) neutrino spectra. Unlike all previous papers on this topic, the neutrino background is then derived by summing up at a given energy the fluxes of each BL Lac in the simulation, all characterised by their own redshift, synchrotron peak energy, gamma-ray flux, etc. Our main result is that BL Lacs as a class can explain the neutrino background seen by IceCube above ~ 0.5 PeV while they only contribute ~ 10% at lower energies, leaving room to some other population(s)/physical mechanism. However, one cannot also exclude the possibility that individual BL Lacs still make a contribution at the ~ 20% level to the IceCube low-energy events. Our scenario makes specific predictions testable in the next few years.

We study the Extended Chandra Deep Field South (E-CDFS) Very Large Array sample, which reaches a flux density limit at 1.4 GHz of 32.5 microJy at the field centre and redshift ~ 4, and covers ~ 0.3 deg^2. Number counts are presented for the whole sample while the evolutionary properties and luminosity functions are derived for active galactic nuclei (AGN). The faint radio sky contains two totally distinct AGN populations, characterised by very different evolutions, luminosity functions, and Eddington ratios: radio-quiet (RQ)/radiative-mode, and radio-loud/jet-mode AGN. The radio power of RQ AGN evolves ~ (1+z)^2.5, similarly to star-forming galaxies, while the number density of radio-loud ones has a peak at ~ 0.5 and then declines at higher redshifts. The number density of radio-selected RQ AGN is consistent with that of X-ray selected AGN, which shows that we are sampling the same population. The unbiased fraction of radiative-mode RL AGN, derived from our own and previously published data, is a strong function of radio power, decreasing from ~ 0.5 at P_1.4GHz ~ 10^24 W/Hz to ~ 0.04$ at P_1.4GHz ~ 10^22 W/Hz. Thanks to our enlarged sample, which now includes ~ 700 radio sources, we also confirm and strengthen our previous results on the source population of the faint radio sky: star-forming galaxies start to dominate the radio sky only below ~ 0.1 mJy, which is also where radio-quiet AGN overtake radio-loud ones.

Blazars are the dominant type of extragalactic sources at microwave and at $\gamma$-ray energies. In the most energetic part of the electromagnetic spectrum (E>100GeV) a large fraction of high Galactic latitude sources are blazars of the High Synchrotron Peaked (HSP) type, that is BL Lac objects with synchrotron power peaking in the UV or in the X-ray band. HSP blazars are remarkably rare, with only a few hundreds of them expected to be above the sensitivity limits of currently available surveys. To find these very uncommon objects, we have devised a method that combines ALLWISE survey data with multi-frequency selection criteria. The sample was defined starting from a primary list of infrared colour-colour selected sources from the ALLWISE all sky survey database, and applying further restrictions on IR-radio and IR-X-ray flux ratios. Using a polynomial fit to the multi-frequency data (radio to X-ray) we estimated synchrotron peak frequencies and fluxes of each object. We assembled a sample including 992 sources, which is currently the largest existing list of confirmed and candidates HSP blazars. All objects are expected to radiate up to the highest $\gamma$-ray photon energies. In fact, 299 of these are confirmed emitters of GeV $\gamma$-ray photons (based on Fermi-LAT catalogues), and 36 have already been detected in the TeV band. The majority of sources in the sample are within reach of the upcoming Cherenkov Telescope Array (CTA), and many may be detectable even by the current generation of Cherenkov telescopes during flaring episodes. The sample includes 425 previously known blazars, 151 new identifications, and 416 HSP candidates (mostly faint sources) for which no optical spectra is available yet. The full 1WHSP catalogue is on-line at http://www.asdc.asi.it/1whsp/ providing a direct link to the SED building tool where multifrequency data can be easily visualised.

The "blazar simplified view" is a new paradigm that explains well the diverse statistical properties of blazars observed over the entire electromagnetic spectrum on the basis of minimal assumptions on blazars' physical and geometrical properties. In this paper, the fourth in a series, we extend the predictions of this paradigm below the sensitivity of existing surveys and estimate the contribution of blazars to the X-ray and gamma-ray extragalactic backgrounds. We find that the integrated light from blazars can explain up to 100% of the cosmic background at energies larger than ~10 GeV, and contribute ~40 to 70% of the gamma-ray diffuse radiation between 100 MeV and 10 GeV. The contribution of blazars to the X-ray background, between 1 and 50 keV, is approximately constant and of the order of 4-5%. On the basis of an interpolation between the estimated flux at X-ray and gamma-ray energies we can expect that the contribution of blazars raises to ~10% at 100 keV, and continues to increase with energy until it becomes the dominant component at a few MeV. Finally, we show that a strong dependence of the synchrotron peak frequency on luminosity, as postulated by the "blazar sequence", is ruled out by the observational data as it predicts a gamma-ray background above a few GeV that is far in excess of the observed value.

We show that a blazar classification in BL Lacs and Flat Spectrum Radio Quasars may not be adequate when it relies solely on the equivalent widths (EWs) of optical lines. In fact, depending on redshift, some strong emission lines can fall in the infrared window and be missed. We selected a sample of BL Lacs with firm redshift identification and good visibility from Paranal. We targeted with the X-shooter spectrograph the five BL Lacs with z > 0.7, i.e., those for which the Halpha line, one of the strongest among blazars, falls outside the optical window and determined the EW of emission lines in both the infrared and optical bands. Two out of five sources show an observed Halpha EW > 5A (one has rest frame EW > 5A) and could be classified as FSRQs by one of the classification schemes used in the literature. A third object is border-line with an observed EW of 4.4 +/- 0.5A. In all these cases Halpha is the strongest emission line detected. The Halpha line of the other two blazars is not detected, but in one case it falls in a region strongly contaminated by sky lines and in the other one the spectrum is featureless. We conclude that a blazar classification based on EW width only can be inaccurate and may lead to an erroneous determination of blazar evolution. This effect is more severe for the BL Lac class, since FSRQs can be misclassified as BL Lacs especially at high redshifts (z > 0.7), where the latter are extremely rare.

The recent IceCube discovery of 0.1-1 PeV neutrinos of astrophysical origin opens up a new era for high-energy astrophysics. Although there are various astrophysical candidate sources, a firm association of the detected neutrinos with one (or more) of them is still lacking. A recent analysis of plausible astrophysical counterparts within the error circles of IceCube events showed that likely counterparts for nine of the IceCube neutrinos include mostly BL Lacs, among which Mrk 421. Motivated by this result and a previous independent analysis on the neutrino emission from Mrk 421, we test the BL Lac-neutrino connection in the context of a specific theoretical model for BL Lac emission. We model the spectral energy distribution (SED) of the BL Lacs selected as counterparts of the IceCube neutrinos using a one-zone leptohadronic model and mostly nearly simultaneous data. The neutrino flux for each BL Lac is self-consistently calculated, using photon and proton distributions specifically derived for every individual source. We find that the SEDs of the sample, although different in shape and flux, are all well fitted by the model using reasonable parameter values. Moreover, the model-predicted neutrino flux and energy for these sources are of the same order of magnitude as those of the IceCube neutrinos. In two cases, namely Mrk 421 and H 1914-194, we find a suggestively good agreement between the model prediction and the detected neutrino flux. Our predictions for all the BL Lacs of the sample are in the range to be confirmed or disputed by IceCube in the next few years of data sampling.

In this Chapter we present the motivation for undertaking both a wide and deep survey with the SKA in the context of studying AGN activity across cosmic time. With an rms down to 1 $\mu$Jy/beam at 1 GHz over 1,000 - 5,000 deg$^2$ in 1 year (wide tier band 1/2) and an rms down to 200 nJy/beam over 10 - 30 deg$^2$ in 2000 hours (deep tier band 1/2), these surveys will directly detect faint radio-loud and radio-quiet AGN (down to a 1 GHz radio luminosity of about $2\times10^{23}$ W/Hz at $z=6$). For the first time, this will enable us to conduct detailed studies of the cosmic evolution of radio AGN activity to the cosmic dawn ($z\gtrsim6$), covering all environmental densities.

Relativistic jets in active galactic nuclei (AGN) are among the most powerful astrophysical objects discovered to date. Indeed, jetted AGN studies have been considered a prominent science case for SKA, and were included in several different chapters of the previous SKA Science Book (Carilli & Rawlings 2004). Most of the fundamental questions about the physics of relativistic jets still remain unanswered, and await high-sensitivity radio instruments such as SKA to solve them. These questions will be addressed specially through analysis of the massive data sets arising from the deep, all-sky surveys (both total and polarimetric flux) from SKA1. Wide-field very-long-baseline-interferometric survey observations involving SKA1 will serve as a unique tool for distinguishing between extragalactic relativistic jets and star forming galaxies via brightness temperature measurements. Subsequent SKA1 studies of relativistic jets at different resolutions will allow for unprecedented cosmological studies of AGN jets up to the epoch of re-ionization, enabling detailed characterization of the jet composition, magnetic field, particle populations, and plasma properties on all scales. SKA will enable us to study the dependence of jet power and star formation on other properties of the AGN system. SKA1 will enable such studies for large samples of jets, while VLBI observations involving SKA1 will provide the sensitivity for pc-scale imaging, and SKA2 (with its extraordinary sensitivity and dynamic range) will allow us for the first time to resolve and model the weakest radio structures in the most powerful radio-loud AGN.

It has become apparent that active galactic nuclei (AGN) may have a significant impact on the growth and evolution of their host galaxies and vice versa but a detailed understanding of the interplay between these processes remains elusive. Deep radio surveys provide a powerful, obscuration-independent tool for measuring both star formation and AGN activity in high-redshift galaxies. Multiwavelength studies of deep radio fields show a composite population of star-forming galaxies and AGN, with the former dominating at the lowest flux densities (S$_{1.4\mathrm{GHz}}<$100~$\mu$Jy). The sensitivity and resolution of the SKA will allow us to identify, and separately trace, the total star formation in the bulges of individual high-redshift galaxies, the related nuclear activity and any star formation occurring on larger scales within a disc. We will therefore gain a detailed picture of the apparently simultaneous development of stellar populations and black holes in the redshift range where both star-formation and AGN activity peak (1$\leq$z$\leq$4). In this chapter we discuss the role of the SKA in studying the connection between AGN activity and galaxy evolution, and the most critical technical requirements for such of studies

We have recently proposed a simplified scenario for blazars in which these sources are classified as flat-spectrum radio quasars or BL Lacs according to the prescriptions of unified schemes, and to a varying combination of Doppler boosted radiation from the jet, emission from the accretion disk, the broad line region, and light from the host galaxy. This scenario has been thoroughly tested through detailed Monte Carlo simulations and reproduces all the main features of existing radio, X-ray, and gamma-ray surveys. In this paper we consider the case of very high energy emission (E > 100 GeV) extrapolating from the expectations for the GeV band, which are in full accordance with the Fermi-LAT survey results, and make detailed predictions for current and future Cherenkov facilities, including the Cherenkov Telescope Array. Our results imply that >~ 100 new blazars can be detected now at very high energy and up to z ~1, consistently with the very recent MAGIC detection of S4 0218+35 at z=0.944.

We present our latest results on the connection between accretion rate and relativistic jet power in AGN, by using a large sample which includes mostly blazars, but contains also some radio--galaxies. The jet power can be traced by $\gamma$--ray luminosity in the case of blazars, and radio luminosity for both classes. The accretion disc luminosity is instead traced by the broad emission lines. Among blazars, we find a correlation between broad line emission and the $\gamma$--ray or radio luminosities, suggesting a direct tight connection between jet power and accretion rate. We confirm that the observational differences between blazar subclasses reflect differences in the accretion regime, but with blazars only we cannot properly access the low--accretion regime. By introducing radio--galaxies, we succeed in observing the fingerprint of the transition between radiatively efficient and inefficient accretion discs in the jetted AGN family. The transition occurs at the standard critical value $L_{\rm d}/L_{\rm Edd}\sim10^{-2}$ and it appears smooth. Below this value, the ionizing luminosity emitted by the accretion structure drops significantly.

IceCube has recently reported the discovery of high-energy neutrinos of astrophysical origin, opening up the PeV (10^15 eV) sky. Because of their large positional uncertainties, these events have not yet been associated to any astrophysical source. We have found plausible astronomical counterparts in the GeV -- TeV bands by looking for sources in the available large area high-energy gamma-ray catalogues within the error circles of the IceCube events. We then built the spectral energy distribution of these sources and compared it with the energy and flux of the corresponding neutrino. Likely counterparts include mostly BL Lacs and two Galactic pulsar wind nebulae. On the one hand many objects, including the starburst galaxy NGC 253 and Centaurus A, despite being spatially coincident with neutrino events, are too weak to be reconciled with the neutrino flux. On the other hand, various GeV powerful objects cannot be assessed as possible counterparts due to their lack of TeV data. The definitive association between high-energy astrophysical neutrinos and our candidates will be significantly helped by new TeV observations but will be confirmed or disproved only by further IceCube data. Either way, this will have momentous implications for blazar jets, high-energy astrophysics, and cosmic-ray and neutrino astronomy.

We present our very recent results on the sub-mJy radio source populations at 1.4 GHz based on the Extended Chandra Deep Field South VLA survey, which reaches ~ 30 {\mu}Jy, with details on their number counts, evolution, and luminosity functions. The sub-mJy radio sky turns out to be a complex mix of star-forming galaxies and radio-quiet AGN evolving at a similar, strong rate and declining radio-loud AGN. While the well-known flattening of the radio number counts below 1 mJy is mostly due to star-forming galaxies, these sources and AGN make up an approximately equal fraction of the sub-mJy sky. Our results shed also light on a fifty-year-old issue, namely radio emission from radio-quiet AGN, and suggest that it is closely related to star formation, at least at z ~ 1.5 - 2. The implications of our findings for future, deeper radio surveys, including those with the Square Kilometre Array, are also discussed. One of the main messages, especially to non-radio astronomers, is that radio surveys are reaching such faint limits that, while previously they were mainly useful for radio quasars and radio galaxies, they are now detecting mostly star-forming galaxies and radio-quiet AGN, i.e., the bulk of the extragalactic sources studied in the infrared, optical, and X-ray bands.

The sub-mJy radio population is a mixture of active systems, that is star forming galaxies (SFGs) and active galactic nuclei (AGNs). We study a sample of 883 radio sources detected at 1.4 GHz in a deep Very Large Array survey of the Extended Chandra Deep Field South (E-CDFS) that reaches a best rms sensitivity of 6 microJy. We have used a simple scheme to disentangle SFGs, radio-quiet (RQ), and radio-loud (RL) AGNs based on the combination of radio data with Chandra X-ray data and mid-infrared observations from Spitzer. We find that at flux densities between about 30 and 100 microJy the radio population is dominated by SFGs (~60%) and that RQ AGNs become increasingly important over RL ones below 100 microJy. We also compare the host galaxy properties of the three classes in terms of morphology, optical colours and stellar masses. Our results show that both SFG and RQ AGN host galaxies have blue colours and late type morphology while RL AGNs tend to be hosted in massive red galaxies with early type morphology. This supports the hypothesis that radio emission in SFGs and RQ AGNs mainly comes from the same physical process: star formation in the host galaxy.

We have recently proposed a new scenario where blazars are classified as flat-spectrum radio quasars or BL Lacs according to the prescriptions of unified schemes, and to a varying combination of Doppler boosted radiation from the jet, emission from the accretion disk, the broad line region, and light from the host galaxy. This mix of different components leads to strong selection effects, which are properly taken into account in our scheme. We describe here the main features of our approach, which solves many long-standing issues of blazar research, give the most important results, and discuss its implications and testable predictions.

We have recently proposed a new simplified scenario where blazars are classified as flat-spectrum radio quasars (FSRQs) or BL Lacs according to the prescriptions of unified schemes, and to a varying combination of Doppler boosted radiation from the jet, emission from the accretion disk, the broad line region, and light from the host galaxy. Here we extend our approach, previously applied to radio and X-ray surveys, to the gamma-ray band and, through detailed Monte Carlo simulations, compare our predictions to Fermi-LAT survey data. Our simulations are in remarkable agreement with the overall observational results, including the percentages of BL Lacs and FSRQs, the fraction of redshift-less objects, and the redshift, synchrotron peak, and gamma-ray spectral index distributions. The strength and large scatter of the oft-debated observed Gamma-ray -- radio flux density correlation is also reproduced. In addition, we predict that almost 3/4 of Fermi-LAT BL Lacs, and basically all of those without redshift determination, are actually FSRQs with their emission lines swamped by the non-thermal continuum and as such should be considered. Finally, several of the currently unassociated high Galactic latitude Fermi sources are expected to be radio-faint blazars displaying a pure elliptical galaxy optical spectrum.

Deep radio observations at 1.4GHz for the Extended Chandra Deep Field South were performed in June through September of 2007 and presented in a first data release (Miller et al. 2008). The survey was made using six separate pointings of the Very Large Array (VLA) with over 40 hours of observation per pointing. In the current paper, we improve on the data reduction to produce a second data release (DR2) mosaic image. This DR2 image covers an area of about a third of a square degree and reaches a best rms sensitivity of 6 uJy and has a typical sensitivity of 7.4 uJy per 2.8" by 1.6" beam. We also present a more comprehensive catalog, including sources down to peak flux densities of five or more times the local rms noise along with information on source sizes and relevant pointing data. We discuss in some detail the consideration of whether sources are resolved under the complication of a radio image created as a mosaic of separate pointings each suffering some degree of bandwidth smearing, and the accurate evaluation of the flux densities of such sources. Finally, the radio morphologies and optical/near-IR counterpart identifications (Bonzini et al. 2012) are used to identify 17 likely multiple-component sources and arrive at a catalog of 883 radio sources, which is roughly double the number of sources contained in the first data release.

In the lead-up to the Square Kilometre Array (SKA) project, several next-generation radio telescopes and upgrades are already being built around the world. These include APERTIF (The Netherlands), ASKAP (Australia), eMERLIN (UK), VLA (USA), e-EVN (based in Europe), LOFAR (The Netherlands), Meerkat (South Africa), and the Murchison Widefield Array (MWA). Each of these new instruments has different strengths, and coordination of surveys between them can help maximise the science from each of them. A radio continuum survey is being planned on each of them with the primary science objective of understanding the formation and evolution of galaxies over cosmic time, and the cosmological parameters and large-scale structures which drive it. In pursuit of this objective, the different teams are developing a variety of new techniques, and refining existing ones. Here we describe these projects, their science goals, and the technical challenges which are being addressed to maximise the science return.

We study a sample of 883 sources detected in a deep Very Large Array survey at 1.4 GHz in the Extended Chandra Deep Field South. The paper focuses on the identification of their optical and infrared (IR) counterparts. We use a likelihood ratio technique that is particularly useful when dealing with deep optical images to minimize the number of spurious associations. We find a reliable counterpart for 95% of our radio sources. Most of the counterparts (74%) are detected at optical wavelengths, but there is a significant fraction (21%) only detectable in the IR. Combining newly acquired optical spectra with data from the literature we are able to assign a redshift to 81% of the identified radio sources (37% spectroscopic). We also investigate the X-ray properties of the radio sources using the Chandra 4 Ms and 250 ks observations. In particular, we use a stacking technique to derive the average properties of radio objects undetected in the Chandra images. The results of our analysis are collected in a new catalog containing the position of the optical/IR counterpart, the redshift information and the X-ray fluxes. It is the deepest multi-wavelength catalog of radio sources, which will be used for future study of this galaxy population.

We put forward a scenario where blazars are classified as flat-spectrum radio quasars, BL Lacs, low synchrotron, or high synchrotron peaked objects according to a varying combination of Doppler boosted radiation from the jet, emission from the accretion disk, the broad line region, and light from the host galaxy. We thoroughly test this new approach, which builds upon unified schemes, using Monte Carlo simulations and show that it can provide simple answers to a number of long-standing open issues. We also demonstrate that selection effects play a very important role in the diversity observed in radio and X-ray samples and in the correlation between luminosity and peak frequency of the synchrotron power (the so-called "blazar sequence"). It turns out that sources so far classified as BL Lacs on the basis of their observed weak, or undetectable, emission lines are of two physically different classes: intrinsically weak-lined objects, more common in X-ray selected samples, and heavily diluted broad-lined sources, more frequent in radio selected samples, which explains some of the confusion in the literature.

We study the quasi-simultaneous near-IR, optical, UV, and X-ray photometry of eleven gamma-ray selected blazars for which redshift estimates larger than 1.2 have been recently provided. Four of these objects turn out to be high-power blazars with the peak of their synchrotron emission between ~ 3 x 10^15 and ~ 10^16 Hz, and therefore of a kind predicted to exist but never seen before. This discovery has important implications for our understanding of physical processes in blazars, including the so-called "blazar sequence", and might also help constraining the extragalactic background light through gamma-ray absorption since two sources are strongly detected even in the 10 - 100 GeV Fermi-LAT band. Based on our previous work and their high powers, these sources are very likely high-redshift flat-spectrum radio quasars with their emission lines swamped by the non-thermal continuum.

In order to trace the instantaneous star formation rate at high redshift, and hence help understanding the relation between the different emission mechanisms related to star formation, we combine the recent 4 Ms Chandra X-ray data and the deep VLA radio data in the Extended Chandra Deep Field South region. We find 268 sources detected both in the X-ray and radio band. The availability of redshifts for $\sim 95$ of the sources in our sample allows us to derive reliable luminosity estimates and the intrinsic properties from X-ray analysis for the majority of the objects. With the aim of selecting sources powered by star formation in both bands, we adopt classification criteria based on X-ray and radio data, exploiting the X-ray spectral features and time variability, taking advantage of observations scattered across more than ten years. We identify 43 objects consistent with being powered by star formation. We also add another 111 and 70 star forming candidates detected only in the radio or X-ray band, respectively. We find a clear linear correlation between radio and X-ray luminosity in star forming galaxies over three orders of magnitude and up to $z \sim 1.5$. We also measure a significant scatter of the order of 0.4 dex, higher than that observed at low redshift, implying an intrinsic scatter component. The correlation is consistent with that measured locally, and no evolution with redshift is observed. Using a locally calibrated relation between the SFR and the radio luminosity, we investigate the L_X(2-10keV)-SFR relation at high redshift. The comparison of the star formation rate measured in our sample with some theoretical models for the Milky Way and M31, two typical spiral galaxies, indicates that, with current data, we can trace typical spirals only at z<0.2, and strong starburst galaxies with star-formation rates as high as $\sim 100 M_\odot yr^{-1}$, up to $z\sim 1.5$.

We propose a scenario where blazars are classified as flat-spectrum radio quasars (FSRQs), BL Lacs, low synchrotron, or high synchrotron peaked objects according to a varying mix of the Doppler boosted radiation from the jet, the emission from the accretion disk, the broad line region, and the light from the host galaxy. In this framework the peak energy of the synchrotron power (nu_peak) in blazars is independent of source type and of radio luminosity. We test this new approach, which builds upon unified schemes, using extensive Monte Carlo simulations and show that it can provide simple answers to a number of long-standing issues including, amongst others, the different cosmological evolution of BL Lacs selected in the radio and X-ray bands, the larger nu_peak values observed in BL Lacs, the fact that high synchrotron peaked blazars are always of the BL Lac type, and the existence of FSRQ/BL Lac transition objects. Objects so far classified as BL Lacs on the basis of their observed weak, or undetectable, emission lines are of two physically different classes: intrinsically weak lined objects, more common in X-ray selected samples, and heavily diluted broad lined sources, more frequent in radio selected samples, which explains some of the confusion in the literature. We also show that strong selection effects are the main cause of the diversity observed in radio and X-ray samples, and that the correlation between luminosity and nu_peak, that led to the proposal of the "blazar sequence", is also a selection effect arising from the comparison of shallow radio and X-ray surveys, and to the fact that high nu_peak-high radio power objects have never been considered because their redshift is not measurable.

We present the evolutionary properties and luminosity functions of the radio sources belonging to the Chandra Deep Field South VLA survey, which reaches a flux density limit at 1.4 GHz of 43 microJy at the field center and redshift ~5, and which includes the first radio-selected complete sample of radio-quiet active galactic nuclei (AGN). We use a new, comprehensive classification scheme based on radio, far- and near-IR, optical, and X-ray data to disentangle star-forming galaxies from AGN and radio-quiet from radio-loud AGN. We confirm our previous result that star-forming galaxies become dominant only below 0.1 mJy. The sub-mJy radio sky turns out to be a complex mix of star-forming galaxies and radio-quiet AGN evolving at a similar, strong rate; non-evolving low-luminosity radio galaxies; and declining radio powerful (P > 3 10^24 W/Hz) AGN. Our results suggest that radio emission from radio-quiet AGN is closely related to star formation. The detection of compact, high brightness temperature cores in several nearby radio-quiet AGN can be explained by the co-existence of two components, one non-evolving and AGN-related and one evolving and star-formation-related. Radio-quiet AGN are an important class of sub-mJy sources, accounting for ~30% of the sample and ~60% of all AGN, and outnumbering radio-loud AGN at < 0.1 mJy. This implies that future, large area sub-mJy surveys, given the appropriate ancillary multi-wavelength data, have the potential of being able to assemble vast samples of radio-quiet AGN by-passing the problems of obscuration, which plague the optical and soft X-ray bands.

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

I present simple but robust estimates of the types of sources making up the faint, sub-microJy radio sky. These include, not surprisingly, star-forming galaxies and radio-quiet AGN but also two "new" populations, that is low radio power ellipticals and dwarf galaxies, the latter likely constituting the most numerous component of the radio sky. I then estimate for the first time the X-ray, optical, and mid-infrared fluxes these objects are likely to have, which are very important for source identification and the synergy between the upcoming SKA and its various pathfinders with future missions in other bands. On large areas of the sky the SKA, and any other radio telescope producing surveys down to at least the microJy level, will go deeper than all currently planned (and past) sky surveys, with the possible exception of the optical ones from PAN-STARRS and the LSST. SPICA, JWST, and in particular the Extremely Large Telescopes (ELTs) will be a match to the next generation radio telescopes but only on small areas and above ~ 0.1 - 1 microJy (at 1.4 GHz), while even IXO will only be able to detect a small (tiny) fraction of the microJy (nanoJy) population in the X-rays. On the other hand, most sources from currently planned all-sky surveys, with the likely exception of the optical ones, will have a radio counterpart within the reach of the SKA. JWST and the ELTs might turn out to be the main, or perhaps even the only, facilities capable of securing optical counterparts and especially redshifts of microJy radio sources. Because of their sensitivity, the SKA and its pathfinders will have a huge impact on a number of topics in extragalactic astronomy including star-formation in galaxies and its co-evolution with supermassive black holes, radio-loudness and radio-quietness in AGN, dwarf galaxies, and the main contributors to the radio background.[ABRIDGED]

I highlight the synergies of the Wide Field X-ray Telescope (WFXT) with the next generation radio surveys, including those to be obtained with the Australian Square Kilometre Array Pathfinder and the Square Kilometre Array, and discuss the overlap between the X-ray and radio source populations. WFXT will benefit greatly from the availability of deep radio catalogues with very high astrometric precision, while on the other hand WFXT data will be vital for the identification of faint radio sources down to ~ 50 microJy.

(Abridged) We carried out an extensive search to identify the counterparts of all the sources listed in the WMAP 3-yr catalogue using literature and archival data. Our work led to the identification of 309 WMAP sources, 98% of which are blazars, radio quasars or radio galaxies. At present, 15 objects still remain without identification due to the lack of optical spectroscopic data or a clear radio counterpart. Our results allow us to define a flux limited sample of 203 high Galactic latitude microwave sources ($f_{41GHz} \ge 1$ Jy, $|b_{\rm II}| > 15^\circ$) which is virtually completely identified (99%). The microwave band is ideally suited for blazar statistical studies since this is the part of the em spectrum that is least affected by the superposition of spectral components of different origin. Using this data-set we derived number counts, luminosity functions and cosmological evolution of blazars and radio galaxies at microwave frequencies. Our results are in good agreement with those found at radio frequencies. The 5 GHz bivariate blazar luminosity functions are similar to those derived from the DXRBS survey, which shows that this sample is representative of the blazar population at 41 GHz. Microwave selected broad- lined quasars are about 6 times more abundant than BL Lacs, a ratio that is similar to, or larger than, that seen at radio and gamma-ray frequencies, once spectral selection effects are taken into account. This strongly suggests that the mechanism responsible for the generation of gamma-rays is, at first order, the same in all blazar types. Our results confirm the findings of Giommi & Colafrancesco (2004, 2006) that blazars and radio galaxies are the largest contaminants of the CMB anisotropy maps. We predict that these sources are also bright gamma-ray sources, most of which will be detected by AGILE and FERMI.

We present the multiwavelength properties of 266 cataloged radio sources identified with 20 and 6 cm VLA deep observations of the CDFS at a flux density limit of 42 \mu Jy at the field centre at 1.4 GHz. These new observations probe the faint end of both the star formation and radio galaxy/AGN population. X-ray data, including upper limits, turn out to be a key factor in establishing the nature of faint radio sources. We find that, while the well-known flattening of the radio number counts below 1 mJy is mostly due to star forming galaxies, these sources and AGN make up an approximately equal fraction of the sub--millijansky sky, contrary to some previous results. We have also uncovered a population of distant AGN systematically missing from many previous studies of sub-millijansky radio source identifications. The AGN include radio galaxies, mostly of the low-power, Fanaroff-Riley I type, and a significant radio-quiet component, which amounts to approximately one fifth of the total sample. We also find that radio detected, X-ray AGN are not more heavily obscured than the X-ray detected AGN. This argues against the use of radio surveys as an efficient way to search for the missing population of strongly absorbed AGN.

We discuss the X-ray properties of the radio sources detected in a deep 1.4 and 5 GHz VLA Radio survey of the Extended Chandra Deep Field South (E-CDFS). Among the 266 radio sources detected, we find 89 sources (1/3 of the total) with X-ray counterparts in the catalog of the 1Ms exposure of the central 0.08 deg^2 (Giacconi et al. 2002; Alexander et al. 2003) or in the catalog of the 250 ks exposure of the 0.3 deg^2 E-CDFS field (Lehmer et al. 2005). For 76 (85%) of these sources we have spectroscopic or photometric redshifts, and therefore we are able to derive their intrinsic properties from X-ray spectral analysis, namely intrinsic absorption and total X-ray luminosities. We find that the population of submillijansky radio sources with X-ray counterparts is composed of a mix of roughly 1/3 star forming galaxies and 2/3 AGN. The distribution of intrinsic absorption among X-ray detected radio sources is different from that of the X-ray selected sample. Namely, the fraction of low absorption sources is at least two times larger than that of X-ray selected sources in the CDFS. This is mostly due to the larger fraction of star forming galaxies present among the X-ray detected radio sources. If we investigate the distribution of intrinsic absorption among sources with L_X>10^42 erg s^-1 in the hard 2-10 keV band (therefore in the AGN luminosity regime), we find agreement between the X-ray population with and without radio emission. In general, radio detected X-ray AGN are not more heavily obscured than the non radio detected AGN. This argues against the use of radio surveys as an efficient way to search for the missing population of strongly absorbed AGN.

We present a detailed analysis of 256 radio sources from our deep (flux density limit of 42 microJy at the field centre at 1.4 GHz) Chandra Deep Field South 1.4 and 5 GHz VLA survey. The radio population is studied by using a wealth of multi-wavelength information in the radio, optical, and X-ray bands. The availability of redshifts for ~ 80% of the sources in our complete sample allows us to derive reliable luminosity estimates for the majority of the objects. X-ray data, including upper limits, for all our sources turn out to be a key factor in establishing the nature of faint radio sources. Due to the faint optical levels probed by this study, we have uncovered a population of distant Active Galactic Nuclei (AGN) systematically missing from many previous studies of sub-millijansky radio source identifications. We find that, while the well-known flattening of the radio number counts below 1 mJy is mostly due to star forming galaxies, these sources and AGN make up an approximately equal fraction of the sub-millijansky sky, contrary to some previous results. The AGN include radio galaxies, mostly of the low-power, Fanaroff-Riley I type, and a significant radio-quiet component, which amounts to approximately one fifth of the total sample. The ratio of radio to optical luminosity depends more on radio luminosity, rather than being due to optical absorption.

We report 20 and 6 cm VLA deep observations of the CDF-S including the Extended CDF-S. We discuss the radio properties of 266 cataloged radio sources, of which 198 are above a 20 cm completeness level reaching down to 43 microJy at the center of the field. Survey observations made at 6 cm over a more limited region covers the original CDF-S to a comparable level of sensitivity as the 20 cm observations. Of 266 cataloged radio sources, 52 have X-ray counterparts in the CDF-S and a further 37 in the E-CDF-S area not covered by the 1 Megasecond exposure. Using a wide range of material, we have found optical or infrared counterparts for 254 radio sources, of which 186 have either spectroscopic or photometric redshifts (Paper II). Three radio sources have no apparent counterpart at any other wavelength. Measurements of the 20 cm radio flux density at the position of each CDF-S X-ray source detected a further 30 radio sources above a conservative 3-sigma detection limit. X-ray and sub-mm observations have been traditionally used as a measure of AGN and star formation activity, respectively. These new observations probe the faint end of both the star formation and radio galaxy/AGN population, as well as the connection between the formation and evolution of stars and SMBHs. Both of the corresponding gravitational and nuclear fusion driven energy sources can lead to radio synchrotron emission. AGN and radio galaxies dominate at high flux densities. Although emission from star formation becomes more prominent at the microjansky levels reached by deep radio surveys, even for the weakest sources, we still find an apparent significant contribution from low luminosity AGN as well as from star formation.

We present the optical and infrared identifications of the 266 radio sources detected at 20 cm with the Very Large Array in the Chandra Deep Field South (Kellermann et al. 2008). Using deep i-band Advanced Camera for Surveys, R-band Wide Field Imager, K-band SOFI/NTT, K-band ISAAC/VLT and Spitzer imaging data, we are able to find reliable counterparts for 254 (~95%) VLA sources. Twelve radio sources remain unidentified and three of them are ``empty fields''. Using literature and our own data we are able to assign redshifts to 186 (~70%) radio sources: 108 are spectroscopic redshifts and 78 reliable photometric redshifts. Based on the rest frame colors and morphological distributions of the host galaxies we find evidences for a change in the submillijansky radio source population: a) above ~ 0.08 mJy early-type galaxies are dominating; b) at flux densities below ~0.08 mJy, starburst galaxies become dominant.

We have observed the Extended Chandra Deep Field South (E-CDF-S) using a mosaic of six deep Very Large Array (VLA) pointings at 1.4GHz. In this paper, we present the survey strategy, description of the observations, and the first data release. The observations were performed during June through September of 2007 and included from 15 to 17 ``classic'' VLA antennas and 6 to 11 that had been retrofitted for the Expanded VLA (EVLA). The first data release consists of a 34.1' x 34.1' image and the attendant source catalog. The image achieves an rms sensitivity of 6.4 uJy per 2.8" x 1.6" beam in its deepest regions, with a typical sensitivity of 8 uJy. The catalog is conservative in that it only lists sources with peak flux densities greater than seven times the local rms noise, yet it still contains 464 sources. Nineteen of these are complex sources consisting of multiple components. Cross matching of the catalog to prior surveys of the E-CDF-S confirms the linearity of the flux density calibration, albeit with a slight possible offset (a few percent) in scale. Improvements to the data reduction and source catalog are ongoing, and we intend to produce a second data release in January 2009.

This paper presents XMM-Newton and Chandra X-ray spectroscopy of ten flat-spectrum radio quasars (FSRQ) which are candidates to have an X-ray spectrum dominated by jet synchrotron emission. In all these FSRQ, which are less strongly relativistically beamed than blazars, a considerable contribution from a power-law component similar to that present in radio-quiet quasars is required to explain the observed X-ray fluxes and X-ray spectral slopes. And as in radio-quiet quasars, their relatively high optical/UV fluxes can be accounted for by a significant contribution from thermal accretion disk emission. The lack of success in finding radio quasars with synchrotron X-rays is attributed to the adopted selection criteria, which were based on the multiwavelength flux ratios of BL Lacertae (BL Lac) objects. A refined selection technique, which additionally involves radio imaging, is proposed to search for these important candidates for the Gamma Ray Large Area Space Telescope (GLAST). On the other hand, the discovered FSRQ with their strong accretion disk signatures are expected to be important probes for studies of the poorly known accretion disk - jet connection.

We computed the chemical evolution of spiral bulges hosting Seyfert nuclei, based on updated chemical and spectro-photometrical evolution models for the bulge of our Galaxy, made predictions about other quantities measured in Seyferts, and modeled the photometry of local bulges. The chemical evolution model contains detailed calculations of the Galactic potential and of the feedback from the central supermassive black hole, and the spectro-photometric model covers a wide range of stellar ages and metallicities. We followed the evolution of bulges in the mass range 10^9 - 10^{11} Msun by scaling the star formation efficiency and the bulge scalelength as in the inverse-wind scenario for elliptical galaxies, and considering an Eddington limited accretion onto the central supermassive black hole. We successfully reproduced the observed black hole-host bulge mass relation. The observed nuclear bolometric luminosity is reproduced only at high redshift or for the most massive bulges; in the other cases, at z = 0 a rejuvenation mechanism is necessary. The black hole feedback is in most cases not significant in triggering the galactic wind. The observed high star formation rates and metal overabundances are easily achieved, as well as the constancy of chemical abundances with redshift and the bulge present-day colours. Those results are not affected if we vary the index of the stellar IMF from x=0.95 to x=1.35; a steeper IMF is instead required in order to reproduce the colour-magnitude relation and the present K-band luminosity of the bulge.

The multi-frequency Sedentary Survey is a flux limited, statistically well-defined sample of highly X-ray dominated BL Lacertae objects (HBLs) which includes 150 sources. In this paper, the third of the series, we report the results of a dedicated optical spectroscopy campaign that, together with results from other independent optical follow up programs, led to the spectroscopic identification of all sources in the sample. We carried out a systematic spectroscopic campaign for the observation of all unidentified objects of the sample using the ESO 3.6m, the KPNO 4m, and the TNG optical telescopes. We present new identifications and optical spectra for 76 sources, 50 of which are new BL Lac objects, 18 are sources previously referred as BL Lacs but for which no redshift information was available, and 8 are broad emission lines AGNs. We find that the multi-frequency selection technique used to build the survey is highly efficient (about 90%) in selecting BL Lacs objects. We present positional and spectroscopic information for all confirmed BL Lac objects. Our data allowed us to determined 36 redshifts out of the 50 new BL Lacs and 5 new redshifts for the previously known objects. The redshift distribution of the complete sample is presented and compared with that of other BL Lacs samples. For 26 sources without recognizable absorption features, we calculated lower limits to the redshift using a method based on simulated optical spectra with different ratios between jet and galaxy emission. For a subsample of 38 object with high-quality spectra, we find a correlation between the optical spectral slope, the 1.4 GHz radio luminosity, and the Ca H&K break value, indicating that for powerful/beamed sources the optical light is dominated by the non-thermal emission from the jet.

We present the preliminary results from two new surveys of blazars that have direct implications on the GLAST detection of extragalactic sources from two different perspectives: microwave selection and a combined deep X-ray/radio selection. The first one is a 41 GHz flux-limited sample extracted from the WMAP 3-yr catalog of microwave point sources. This is a statistically well defined sample of about 200 blazars and radio galaxies, most of which are expected to be detected by GLAST. The second one is a new deep survey of Blazars selected among the radio sources that are spatially coincident with serendipitous sources detected in deep X-ray images (0.3-10 keV) centered on the Gamma Ray Bursts (GRB) discovered by the Swift satellite. This sample is particularly interesting from a statistical viewpoint since a) it is unbiased as GRBs explode at random positions in the sky, b) it is very deep in the X-ray band (\fx \simgt $10^{-15}$ \ergs) with a position accuracy of a few arc-seconds, c) it will cover a fairly large (20-30 square deg.) area of sky, d) it includes all blazars with radio flux (1.4 GHz) larger than 10 mJy, making it approximately two orders of magnitude deeper than the WMAP sample and about one order of magnitude deeper than the deepest existing complete samples of radio selected blazars, and e) it can be used to estimate the amount of unresolved GLAST high latitude gamma-ray background and its anisotropy spectrum.

I describe the different classes of Active Galactic Nuclei (AGN) and the basic tenets of unified schemes. I then review the properties of the extragalactic sources detected in the GeV and TeV bands, showing that the vast majority of them belong to the very rare blazar class. I further discuss the kind of AGN GLAST is likely to detect, making some predictions going from the obvious to the likely, all the way to the less probable.

We report the discovery of ROXA J081009.9+384757.0 = SDSS J081009.9+384757.0, a z=3.95 blazar with a highly unusual Spectral Energy Distribution (SED). This object was first noticed as a probable high $f_x/f_r$, high-luminosity blazar within the error region of a $\approx 10^{-12}$ erg/s cm$^2$ ROSAT source which, however, also included a much brighter late-type star. We describe the results of a recent Swift observation that establishes beyond doubt that the correct counterpart of the X-ray source is the flat spectrum radio quasar. With a luminosity well in excess of $10^{47}$ erg/s, ROXA J081009.9+384757.0 is therefore one of the most luminous blazars known. We consider various possibilities for the nature of the electromagnetic emission from this source. In particular, we show that the SED is consistent with that of a blazar with synchrotron power peaking in the hard X-ray band. If this is indeed the case, the combination of high-luminosity and synchrotron peak in the hard-X-ray band contradicts the claimed anti-correlation between luminosity and position of the synchrotron peak usually referred to as the "blazar sequence". An alternative possibility is that the X-rays are not due to synchrotron emission, in this case the very peculiar SED of ROXA J081009.9+384757.0 would make it the first example of a new class of radio loud AGN.

Our knowledge of the blazar surface densities and luminosity functions, which are fundamental parameters, relies still on samples at relatively high flux limits. As a result, our understanding of this rare class of active galactic nuclei is mostly based on relatively bright and intrinsically luminous sources. We present the radio number counts, evolutionary properties, and luminosity functions of the faintest blazar sample with basically complete (~ 95%) identifications. Based on the Deep X-ray Radio Blazar Survey (DXRBS), it includes 129 flat-spectrum radio quasars (FSRQ) and 24 BL Lacs down to a 5 GHz flux and power ~ 50 mJy and ~ 10^{24} W/Hz, respectively, an order of magnitude improvement as compared to previously published (radio-selected) blazar samples. DXRBS FSRQ are seen to evolve strongly, up to redshift ~ 1.5, above which high-power sources show a decline in their comoving space density. DXRBS BL Lacs, on the other hand, do not evolve. High-energy (HBL) and low-energy (LBL) peaked BL Lacs share the same lack of cosmological evolution, which is at variance with some previous results. The observed luminosity functions are in good agreement with the predictions of unified schemes, with FSRQ getting close to their expected minimum power. Despite the fact that the large majority of our blazars are FSRQ, BL Lacs are intrinsically ~ 50 times more numerous. Finally, the relative numbers of HBL and LBL in the radio and X-ray bands are different from those predicted by the so-called "blazar sequence" and support a scenario in which HBL represent a small minority (~ 10%) of all BL Lacs.

We analyse the 267 radio sources from our deep (flux limit of 42 microJy at the field center at 1.4 GHz) Chandra Deep Field South 1.4 and 5 GHz VLA survey. The radio population is studied by using a wealth of multi-wavelength information, including morphology and spectral types, in the radio, optical, and X-ray bands. The availability of redshifts for ~ 70% of our sources allows us to derive reliable luminosity estimates for the majority of the objects. Contrary to some previous results, we find that star-forming galaxies make up only a minority (~ 1/3) of sub-mJy sources, the bulk of which are faint radio galaxies, mostly of the Fanaroff-Riley I type.

The "blazar sequence" posits that the most powerful BL Lacertae objects and flat-spectrum radio quasars should have relatively small synchrotron peak frequencies, nu_peak, and that the least powerful such objects should have the highest nu_peak values. This would have strong implications for our understanding of jet formation and physics and the possible detection of powerful, moderately high-redshift TeV blazars. I review the validity of the blazar sequence by using the results of very recent surveys and compare its detailed predictions against observational data. I find that the blazar sequence in its simplest form is ruled out. However, powerful flat-spectrum radio quasars appear not to reach the nu_peak typical of BL Lacs. This could indeed be related to some sort of sequence, although it cannot be excluded that it is instead due to a selection effect.

The large majority of BL Lacertae objects belonging to the 1 Jy sample, the class prototype for radio-selected sources, are thought to emit most of their synchrotron power in the far IR band. Ironically, this spectral region is very sparsely sampled, with only a minority of the objects having IRAS data (most of them being upper limits or low-quality detections). We aim at filling this IR gap by presenting new, simultaneous ISOCAM and ISOPHOT observations over the 7 - 200 micron range for half the sample. A measurement of the position of the synchrotron peak frequency, nu_peak, can provide information about particle acceleration mechanisms and constrain the inverse Compton radiation that will be detected by up-coming new gamma-ray missions. We have observed 17 1 Jy BL Lacertae objects with the camera and the photometer on board the Infrared Space Observatory (ISO) satellite. The ISOPHOT data reduction was done employing a novel correction, which mitigates the effect of chopping for faint sources. Using our new ISO data, complemented by nearly-simultaneous radio and optical observations for 10 and 4 objects respectively, and other multi-frequency data, we have built the spectral energy distributions of our sources (plus a previously published one) and derived the rest-frame nu_peak. Its distribution is centered at 10^13 Hz (30 micron) and is very narrow, with 60% of the BL Lacs in the 1 - 3 10^13 Hz range. Given our set of simultaneous infrared data, these represent the best determinations available of the synchrotron peak frequencies for low-energy peaked BL Lacs. A comparison with previous such estimates, based on non-simultaneous optical and near IR data, may indicate strong nu_peak variations in a number of sources, possibly associated with large flares as observed in the high-energy peaked BL Lac MKN 501. (abridged)

The Astrophysical Virtual Observatory (AVO) initiative, jointly funded by the European Commission and six European organisations, had the task of creating the foundations of a regional scale infrastructure by conducting a research and demonstration programme on the VO scientific requirements and necessary technologies. The AVO project is now formally concluded. I highlight AVO's main achievements and then describe its successor, the EURO-VO project. With its three new interlinked structures, the Data Centre Alliance, the Facility Centre, and the Technology Centre, the EURO-VO is the logical next step for the deployment of an operational VO in Europe.

We present the first deep VLA radio images of flat-spectrum radio quasars (FSRQ) with multiwavelength emission properties similar to those of BL Lacs with synchrotron X-rays. Our observations of twenty-five of these sources show that their radio morphologies are similar to those of other radio quasars. However, their range of extended powers is more similar to that of BL Lacertae objects (BL Lacs) and extends down to the low values typical of FR I radio galaxies. Five out of our nine lobe-dominated sources have extended radio powers in the range typical of both FR I and FR II radio galaxies, but their extended radio structure is clearly FR II-like. Therefore, we have not yet found a large population of radio quasars hosted by FR Is. Two thirds of our sources have a core-dominated radio morpholgy and thus X-rays likely dominated by the jet. We find that their ratios of radio core to total X-ray luminosity are low and in the regime indicative of synchrotron X-rays. This result shows that also blazars with strong emission lines can produce jets of high-energy synchrotron emission and undermines at least in part the ``blazar sequence'' scenario which advocates that particle Compton cooling by an external radiation field governs the frequency of the synchrotron emission peak.

The Virtual Observatory (VO) will revolutionise the way we do Astronomy by allowing easy access to all astronomical data and by making the handling and analysis of datasets at various locations across the globe much simpler and faster. I report here on the need for the VO and its status in Europe, concentrating on the recently started EURO-VO project, and then give two specific applications of VO tools to high-energy astrophysics.

An analysis of SNIa events in early type galaxies from the Cappellaro et al (1999) database provides conclusive evidence that the rate of type Ia Supernovae (SNe) in radio-loud galaxies is about 4 times higher than the rate measured in radio-quiet galaxies, i.e. SNIa-rate$(radio-loud galaxies) = 0.43^{+0.19}_{-0.14}h^2_{75}$ SNu as compared to SNIa-rate$(radio-quiet galaxies) = 0.11^{+0.06}_{-0.03}h^2_{75}$ SNu. The actual value of the enhancement is likely to be in the range $\sim 2-7$ (P$\sim 10^{-4}$). This finding puts on robust empirical grounds the results obtained by Della Valle & Panagia (2003) on the basis of a smaller sample of SNe. We analyse the possible causes of this result and conclude that the enhancement of SNIa explosion rate in radio-loud galaxies has the same origin as their being strong radio sources, but there is no causality link between the two phenomena. We argue that repeated episodes of interaction and/or mergers of early type galaxies with dwarf companions, on times-scale of about 1 Gyr, are responsible for inducing both strong radio activity observed in $\sim$14% of early type galaxies and to supply an adequate number of SNIa progenitors to the stellar population of ellipticals.

The Virtual Observatory is now mature enough to produce cutting-edge science results. The exploitation of astronomical data beyond classical identification limits with interoperable tools for statistical identification of sources has become a reality. I present the discovery of 68 optically faint, obscured (i.e., type 2) active galactic nuclei (AGN) candidates in the two GOODS fields using the Astrophysical Virtual Observatory (AVO) prototype. Thirty-one of these sources have high estimated X-ray powers (>10^44 erg/s) and therefore qualify as optically obscured quasars, the so-called QSO 2. The number of these objects in the GOODS fields is now 40, an improvement of a factor > 4 when compared to the only 9 such sources previously known. By going ~ 3 magnitudes fainter than previously known type 2 AGN in the GOODS fields the AVO is sampling a region of redshift -- power space much harder to reach with classical methods. I also discuss the AVO move to our next phase, the EURO-VO, and our short-term plans to continue doing science with the Virtual Observatory.

The multi-frequency `Sedentary Survey' is a deep, statistically complete, radio flux limited sample comprising 150 BL Lacertae objects distinguished by their extremely high X-ray to radio flux ratio, ranging from five hundred to over five thousand times that of typical BL Lacs discovered in radio surveys. This paper presents the final, 100% identified, catalog together with the optical, X-ray and broad-band SEDs constructed combining literature multi-frequency data with non-simultaneous optical observations and BeppoSAX X-ray data, when available. The SEDs confirm that the peak of the synchrotron power in these objects is located at very high energies. BeppoSAX wide band X-ray observations show that, in most cases, the X-ray spectra are convex and well described by a logarithmic parabola model peaking (in a E f(E) vs E representation) between 0.02 to several keV. Owing to the high synchrotron energies involved most of the sources in the catalog are likely to be TeV emitters, with the closest and brightest ones probably detectable by the present generation of Cherenkov telescopes. The optical spectrum of about one fourth of the sources is totally featureless. Because this implies that the non-thermal emission must be well above that of the host galaxy, these objects are likely to be the most powerful sources in the survey and therefore be examples of the yet unreported high radio luminosity-high energy peaked BL Lacs.

The Virtual Observatory (VO) will revolutionise the way we do Astronomy, by allowing easy access to all astronomical data and by making the handling and analysis of datasets at various locations across the globe much simpler and faster. I report here on the need for the VO and its status in Europe, including the first ever VO-based astronomical paper, and then give two specific applications of VO tools to open problems of A-stars research.

We use Virtual Observatory (VO) tools to identify optically faint, obscured (i.e., type 2) active galactic nuclei (AGN) in the two Great Observatories Origins Deep Survey (GOODS) fields. By employing publicly available X-ray and optical data and catalogues we discover 68 type 2 AGN candidates. The X-ray powers of these sources are estimated by using a previously known correlation between X-ray luminosity and X-ray-to-optical flux ratio. Thirty-one of our candidates have high estimated powers (Lx > 10^44 erg/s) and therefore qualify as optically obscured quasars, the so-called ``QSO 2''. Based on the derived X-ray powers, our candidates are likely to be at relatively high redshifts, z ~ 3, with the QSO 2 at z ~ 4. By going ~ 3 magnitudes fainter than previously known type 2 AGN in the two GOODS fields we are sampling a region of redshift - power space which was previously unreachable with classical methods. Our method brings to 40 the number of QSO 2 in the GOODS fields, an improvement of a factor ~ 4 when compared to the only 9 such sources previously known. We derive a QSO 2 surface density down to 10^-15 erg/cm^2/s in the 0.5 - 8 keV band of >~ 330/deg^2, ~ 30% of which is made up of previously known sources. This is larger than current estimates and some predictions and suggests that the surface density of QSO 2 at faint flux limits has been underestimated. This work demonstrates that VO tools are mature enough to produce cutting-edge science results by exploiting astronomical data beyond ``classical'' identification limits (R <~ 25) with interoperable tools for statistical identification of sources using multiwavelength information.

We report results from {\it XMM-Newton} observations of thirteen X-ray bright BL Lacertae objects, selected from the {\it Einstein} Slew Survey sample. The spectra are generally well fit by power-law models, with four objects having hard ($\alpha<1; F_\nu \propto \nu^{-\alpha}$) spectra that indicates synchrotron peaks at $>5$ keV. None of our spectra show line features, indicating that soft X-ray absorption ``notches'' must be rare amongst BL Lacs, rather than common or ubiquitous as had previously been asserted. We find significant curvature in most of the spectra. This curvature is almost certainly intrinsic, as it appears nearly constant from 0.5 to 6 keV, an observation which is inconsistent with the small columns seen in these sources.

The Great Observatories Origins Deep Survey provides significant constraints on the space density of less luminous QSOs at high redshift, which is particularly important to understand the interplay between the formation of galaxies and super-massive black holes and to measure the QSO contribution to the UV ionizing background. We present the results of a search for high-z QSOs, identified in the two GOODS fields on the basis of deep imaging in the optical (with HST) and X-ray (Chandra), and discuss the allowed space density of QSOs in the early universe.

Blazars are currently separated into BL Lacertae objects (BL Lacs) and flat spectrum radio quasars (FSRQ) based on the strength of their emission lines. This is done rather arbitrarily by defining a diagonal line in the Ca H&K break value -- equivalent width plane, following Marcha et al. We readdress this problem and put the classification scheme for blazars on firm physical grounds. We study ~100 blazars and radio galaxies from the Deep X-ray Radio Blazar Survey (DXRBS) and 2 Jy radio survey and find a significant bimodality for the narrow emission line [OIII] 5007. This suggests the presence of two physically distinct classes of radio-loud AGN. We show that all radio-loud AGN, blazars and radio galaxies, can be effectively separated into weak- and strong-lined sources using the [OIII] 5007 -- [OII] 3727 equivalent width plane. This plane allows one to disentangle orientation effects from intrinsic variations in radio-loud AGN. Based on DXRBS, the strongly beamed sources of the new class of weak-lined radio-loud AGN are made up of BL Lacs at the ~75 per cent level, whereas those of the strong-lined radio-loud AGN include mostly (~97 per cent) quasars.

We present new BeppoSAX LECS and MECS observations, covering the energy range 0.1 - 10 keV (observer's frame), of four BL Lacertae objects selected from the 1 Jy sample. All sources display a flat (alpha_x ~ 0.7) X-ray spectrum, which we interpret as inverse Compton emission. One object shows evidence for a low-energy steepening (Delta alpha_x ~ 0.9) which is likely due to the synchrotron component merging into the inverse Compton one around ~ 2 keV. A variable synchrotron tail would explain why the ROSAT spectra of our sources are typically steeper than the BeppoSAX ones (Delta alpha_x} ~ 0.7). The broad-band spectral energy distributions fully confirm this picture and model fits using a synchrotron inverse Compton model allow us to derive the physical parameters (intrinsic power, magnetic field, etc.) of our sources. By combining the results of this paper with those previously obtained on other sources we present a detailed study of the BeppoSAX properties of a well-defined sub-sample of 14 X-ray bright (f_x (0.1 - 10 keV) > 3 x 10^{-12} erg/cm^2/s) 1-Jy BL Lacs. We find a very tight proportionality between nearly simultaneous radio and X-ray powers for the 1-Jy sources in which the X-ray band is dominated by inverse Compton emission, which points to a strong link between X-ray and radio emission components in these objects.

We present a sample of 17 high-redshift (3.5<z<5.2) QSO candidates in the 320 sq.arcmin area of the Great Observatories Origins Deep Survey, selected in the magnitude range 22.45<z_{850}<25.25 using deep imaging with the Advanced Camera for Surveys onboard the Hubble Space Telescope and the Advanced CCD Imaging Spectrometer onboard the Chandra X-ray Observatory. On the basis of seven spectroscopic and ten photometric redshifts we estimate that the final sample will contain between two and four QSOs with 4<z<5.2. A dearth of high-redshift, moderate-luminosity (M_{145}=~-23) QSOs is observed with respect to predictions based on a) the extrapolation of the z~2.7 luminosity function (LF), according to a pure luminosity evolution calibrated by the results of the Sloan Digital Sky Survey; and b) a constant universal efficiency in the formation of super-massive black holes (SMBHs) in dark-matter halos. Evidence is gathered in favor of a density evolution of the LF at high redshift and of a suppression of the formation or feeding of SMBHs in low-mass halos.

We describe a possible new class of X-ray sources that have robust detections in ultra-deep Chandra data, yet have no detections at all in our deep multi-band GOODS Hubble Space Telescope (HST) ACS images, which represent the highest quality optical imaging obtained to date on these fields. These extreme X-ray / Optical ratio sources ("EXO"s) have values of Fx/Fopt at least an order of magnitude above those generally found for other AGN, even those that are harbored by reddened hosts. We thus infer two possible scenarios: (1) if these sources lie at redshifts z < 6, then their hosts need to be exceedingly underluminous, or more reddened, compared with other known sources; (2) if these sources lie above z ~ 6-7, such that even their Lyman-alpha emission is redshifted out of the bandpass of our ACS z(850) filter, then their optical and X-ray fluxes can be accounted for in terms of relatively normal L* hosts and moderate-luminosity AGN.

We use statistical results from a large sample of about 500 blazars, based on two surveys, the Deep X-ray Radio Blazar Survey (DXRBS), nearly complete, and the RASS-Green Bank survey (RGB), to provide new constraints on the spectral energy distribution of blazars, particularly flat-spectrum radio quasars (FSRQ). This reassessment is prompted by the discovery of a population of FSRQ with spectral energy distribution similar to that of high-energy peaked BL Lacs. The fraction of these sources is sample dependent, being ~ 10% in DXRBS and ~ 30% in RGB (and reaching ~ 80% for the Einstein Medium Sensitivity Survey). We show that these ``X-ray strong'' radio quasars, which had gone undetected or unnoticed in previous surveys, indeed are the strong-lined counterparts of high-energy peaked BL Lacs and have synchrotron peak frequencies, nu_peak, much higher than ``classical'' FSRQ, typically in the UV band for DXRBS. Some of these objects may be 100 GeV - TeV emitters, as are several known BL Lacs with similar broadband spectra. Our large, deep, and homogeneous DXRBS sample does not show anti-correlations between nu_peak and radio, broad line region, or jet power, as expected in the so-called ``blazar sequence'' scenario. However, the fact that FSRQ do not reach X-ray-to-radio flux ratios and nu_peak values as extreme as BL Lacs and the elusiveness of high nu_peak - high-power blazars suggest that there might be an intrinsic, physical limit to the synchrotron peak frequency that can be reached by strong-lined, powerful blazars. Our findings have important implications for the study of jet formation and physics and its relationship to other properties of active galactic nuclei.

The relation between accretion disc (thermal emission) and jet (non-thermal emission) in blazars is still a mystery as, typically, the beamed jet emission swamps the disc even in the ultraviolet band where disc emission peaks. In this paper we estimate the accretion disc component for 136 flat-spectrum radio quasars selected from the Deep X-ray Radio Blazar Survey. We do this by deriving the accretion disc spectrum from the mass and accretion rate onto the central black hole for each object, estimated using the emission line widths and the power emitted from the broad line region. We find that non-thermal emission dominates the optical/UV band of our sources. The thermal component, in fact, is, on average, ~ 15 per cent of the total and > 90 per cent of the objects in the sample have a thermal component < 0.5 of the total luminosity. We then estimate the integrated disc and kinetic jet powers and find that, on average, the disc luminosity is ~ 1 to 20 times the jet power (depending on the uncertainties in the estimation of the latter quantity). A comparison with previous, independent results favours a scenario in which jet and disk powers are of the same order of magnitude. Extraction of energy from a rotating black hole via the ``Blandford-Znajek'' mechanism fails to explain the estimated jet power in the majority of our sources. Finally, we find that the typical masses for our sources are ~ 5 10^8 solar masses and that, contrary to previous claims, about one quarter of our radio quasars have relatively small (< 3 10^8 solar masses) black hole mass.

Using three complete, radio flux limited, blazar samples we compare the LogN-LogS and the preliminary radio luminosity function of the general population of BL Lacs to those of the subclass of high energy synchrotron peaked (HBL) BL Lacs. We also examine recent results on the cosmological evolution in different samples of BL Lacs and we investigate the controversial issue of the correlation between the synchrotron peak frequency and radio luminosity in BL Lacertae objects. We find that the fraction of HBL objects is approximately the same at all observed radio fluxes and luminosities implying that there cannot be any strong correlation between the position of the synchrotron peak and radio luminosity. The amount of cosmological evolution in BL Lacs is confirmed to be low and negative at low radio fluxes, although the large number of objects without redshift prevents a precise estimation. At high radio fluxes the amount of cosmological evolution is zero or slightly positive but this could be induced by a possible contamination with Flat Spectrum Radio Quasars.

We present new BeppoSAX LECS, MECS, and PDS observations of four flat-spectrum radio quasars (FSRQ) having effective spectral indices alpha_ro and alpha_ox typical of high-energy peaked BL Lacs. Our sources have X-ray-to-radio flux ratios on average ~ 70 times larger than ``classical'' FSRQ and lie at the extreme end of the FSRQ X-ray-to-radio flux ratio distribution. The collected data cover the energy range 0.1 - 10 keV (observer's frame), reaching ~ 100 keV for one object. The BeppoSAX band in one of our sources, RGB J1629+4008, is dominated by synchrotron emission peaking at ~ 2 x 10^16 Hz, as also shown by its steep (energy index alpha_x ~ 1.5) spectrum. This makes this object the FIRST known FSRQ whose X-ray emission is not due to inverse Compton radiation. Two other sources display a flat BeppoSAX spectrum (alpha_x ~ 0.7), with weak indications of steepening at low X-ray energies. The combination of BeppoSAX and ROSAT observations, (non-simultaneous) multifrequency data, and a synchrotron inverse Compton model suggest synchrotron peak frequencies ~ 10^15 Hz, although a better coverage of their spectral energy distributions is needed to provide firmer values. If confirmed, these values would be typical of ``intermediate'' BL Lacs for which the synchrotron and inverse Compton components overlap in the BeppoSAX band. Our sources, although firmly in the radio-loud regime, have powers more typical of high-energy peaked BL Lacs than of FSRQ, and indeed their radio powers put them near the low-luminosity end of the FSRQ luminosity function. We discuss this in terms of an anti-correlation between synchrotron peak frequency and total power, based on physical arguments, and also as possibly due to a selection effect.

We investigate why BL Lacertae objects (BL Lacs) have values of the Ca H&K break (a stellar absorption feature) lower than low-power radio galaxies and if its use is justified to separate the two classes. For this purpose we relate this parameter to the radio and optical core emissions, as well as to the X-ray powers, for a sample of ~90 radio sources. We find that the Ca H&K break value decreases with increasing jet powers, and that it also anti-correlates with the radio core dominance parameter but not with extended radio emission. Based on this we conclude that the Ca H&K break value of BL Lacs and radio galaxies is a suitable indicator of orientation. From the luminosity ratios between objects with low and high Ca H&K break values we constrain the average Lorentz factors for BL Lacs and low-power radio galaxies in the radio and X-ray band to Gamma ~ 2 -- 4 and derive average viewing angles for the galaxies. Our values are in agreement with results from independent methods. We find that the correlations between Ca H&K break and radio core and X-ray luminosity hold independently for low- (LBL) and high-energy peaked BL Lacs (HBL). We derive average viewing angles for their parent populations, which turn out to be similar to the ones for our entire sample, and compare for the first time the luminosities of LBL and HBL at different orientations.

The high redshift galaxy MS1512-cB58 (z=2.724, m_V=20.64) has been observed with the very efficient high resolution echelle spectrograph VLT/UVES. Although this is a very challenging observational program for a Southern hemisphere telescope (the galaxy is located at +36 deg declination), high resolution spectra (FWHM ~ 26 km/s) have revealed, with unprecedented detail along a galaxy sight line, the Lyman-alpha forest due to intervening clouds in the intergalactic medium (IGM). The mean depression D_A due to IGM absorption blueward of the galaxy Ly-alpha wavelength and the number density dn/dz of Ly-alpha clouds have been compared with equivalent results obtained for QSO sight lines at similar redshifts. Our results indicate a possible excess of absorption close to the galaxy. The mean depression at ~ 150 h_65^-1 Mpc comoving (Omega_m=0.3, Omega_Lambda=0.7) from the galaxy is D_A=0.36+/-0.03, to be compared with 0.22+/-0.04, expected from a best fit to QSO sight lines. In the same region (z=2.610), the number density of lines with HI column density in excess of 10^14 atoms/cm2 is also about 3 sigma larger than expected. This high density region is at least 60 h_65^-1 Mpc comoving wide, but the large Ly-alpha absorption of the galaxy itself prevents us from detecting a possible structure extending down to the galaxy. This excess of Ly-alpha clouds is suggestive of two possible scenarios. One is the presence of a super cluster of Ly-alpha clouds not associated with cB58. The other is a high density of gas associated with the environment of cB58. Indeed, a hint of the complexity of cB58 and possibly its environment is given by the huge velocity range (almost 1000 km/s) between the optical emission of star forming regions and UV absorption of its interstellar medium.

We present new BeppoSAX observations of seven BL Lacertae objects selected from the 1 Jy sample plus one additional source. The collected data cover the energy range 0.1 - 10 keV (observer's frame), reaching ~ 50 keV for one source (BL Lac). All sources characterized by a peak in their multifrequency spectra at infrared/optical energies (i.e., of the LBL type) display a relatively flat (alpha_x ~ 0.9) X-ray spectrum, which we interpret as inverse Compton emission. Four objects (2/3 of the LBL) show some evidence for a low-energy steepening which is likely due to the synchrotron tail merging into the inverse Compton component around ~ 1 - 3 keV. If this were generally the case with LBL, it would explain why the 0.1 - 2.4 keV ROSAT spectra of our sources are systematically steeper than the BeppoSAX ones (by ~ 0.5 in alpha_x). The broad-band spectral energy distributions fully confirm this picture and a synchrotron inverse Compton model allows us to derive the physical parameters (intrinsic power, magnetic field, etc.) of our sources. Combining our results with those obtained by BeppoSAX on BL Lacs covering a wide range of synchrotron peak frequencies, we confirm and clarify the dependence of the X-ray spectral index on synchrotron peak frequency originally found in ROSAT data.

As an overview of blazar research, we briefly discuss some of the themes of this conference, including the various types of blazars, blazar samples and selection effects, unified schemes, evolution, and jet physics.

I address the need for deep blazar surveys by showing that our current understanding of blazars is based on a relatively small number of intrinsically luminous sources. I then review the on-going deeper surveys, addressing in particular their limits and limitations. Finally, I present some preliminary results on the evolutionary properties of faint blazars as derived from the Deep X-ray Radio Blazar Survey (DXRBS).

We present the results of our methods to "mine" the blazar sky, i.e., select blazar candidates with very high efficiency. These are based on the cross-correlation between public radio and X-ray catalogs and have resulted in two surveys, the Deep X-ray Radio Blazar Survey (DXRBS) and the "Sedentary" BL Lac survey. We show that data mining is vital to select sizeable, deep samples of these rare active galactic nuclei and we touch upon the identification problems which deeper surveys will face.

We have searched the archived, pointed ROSAT Position Sensitive Proportional Counter data for blazars by correlating the WGACAT X-ray database with several publicly available radio catalogs, restricting our candidate list to serendipitous X-ray sources with a flat radio spectrum (alpha_r <= 0.70). This makes up the Deep X-ray Radio Blazar Survey (DXRBS). Here we present new identifications and spectra for 106 sources, including 86 radio-loud quasars, 11 BL Lacertae objects, and 9 narrow-line radio galaxies. Together with our previously published objects and already known sources, our sample now contains 298 identified objects: 234 radio-loud quasars (181 flat-spectrum quasars: FSRQ [alpha_r <= 0.50] and 53 steep-spectrum quasars: SSRQ), 36 BL Lacs, and 28 narrow-line radio galaxies. Redshift information is available for 96% of these. Thus our selection technique is ~ 90% efficient at finding radio-loud quasars and BL Lacs. Reaching 5 GHz radio fluxes ~ 50 mJy and 0.1-2.0 keV X-ray fluxes a few x 10^-14 erg/cm^2/s, DXRBS is the faintest and largest flat-spectrum radio sample with nearly complete (~ 85%) identification. We review the properties of the DXRBS blazar sample, including redshift distribution and coverage of the X-ray-radio power plane for quasars and BL Lacs. Additionally, we touch upon the expanded multiwavelength view of blazars provided by DXRBS. By sampling for the first time the faint end of the radio and X-ray luminosity functions, this sample will allow us to investigate the blazar phenomenon and the validity of unified schemes down to relatively low powers.

The rareness of blazars, combined with the previous history of relatively shallow, single-band surveys, has dramatically colored our perception of these objects. Despite a quarter-century of research, it is not at all clear whether current samples can be combined to give us a relatively unbiased view of blazar properties, or whether they present a view so heavily affected by biases inherent in single-band surveys that a synthesis is impossible. We will use the coverage of X-ray/radio flux space for existing surveys to assess their biases. Only new, deeper blazar surveys approach the level needed in depth and coverage of parameter space to give us a less biased view of blazars. These surveys have drastically increased our knowledge of blazars' properties. We will specifically review the discovery of ``blue'' blazars, objects with broad emission lines but broadband spectral characteristics similar to HBL BL Lac objects.

We review the main results from several radio, X-ray and multi-frequency surveys on the topic of cosmological evolution of BL Lacertae objects. Updated findings on BL Lac evolution following the recent identification of many sources in the ``Sedentary Multi-Frequency survey'' are also discussed. By means of extensive Monte Carlo simulations we test some possible explanations for the peculiar cosmological evolution of BL Lacs. We find that a dependence of the relativistic Doppler factor on radio luminosity (as expected within the beaming scenario) may induce low values of V/V_max and that both edge effects at the low luminosity end of the BL Lacs radio luminosity function, and incompleteness at faint optical magnitudes may be the cause of the low V/V_max found for extreme HBL sources in X-ray selected samples.

Blazars can be classified on the basis of their overall Spectral Energy Distribution (SED). BL Lac objects are usually divided in LBL or HBL (Low or High energy peaked BL Lacs), according to the peak frequency of the synchrotron emission, if in the optical or UV-soft-X band respectively. FSRQs instead are characterized by synchrotron peaks mainly at IR-optical frequencies, similarly to LBLs. Here we report on recent BeppoSAX observations which are unveiling the high energy branch of the range of synchrotron peak frequencies. Four new "extreme" HBLs have been discovered, one of which (1ES 1426+428) peaks near or above 100 keV, in a quiescent state.

The HST NICMOS Camera 2 was used for H-band imaging of 12 BL Lacertae objects taken from the larger sample observed with the WFPC2 in the R band (Urry et al. 2000; Scarpa et al. 2000). Ten of the 12 BL Lacs are clearly resolved, and the detected host galaxies are large, bright ellipticals with average H-band absolute magnitude M=-26.2+-0.45 mag and effective radius 10+-5 kpc. The rest-frame integrated color of the host galaxies is on average R-H=2.3+-0.3, consistent with the value for both radio galaxies and normal, non-active elliptical galaxies, and indicating the dominant stellar population is old. The host galaxies tend to be bluer in their outer regions than in their cores, with average color gradient Delta(R-H)/Delta(log r)=-0.2 mag, again consistent with results for normal non-active elliptical galaxies. The infrared Kormendy relation, derived for the first time for BL Lac host galaxies, is m(e) = 3.8*log(R)+14.8 (where m(e) is the surface brightness at the effective radius R), fully in agreement with the relation for normal ellipticals. The close similarity between BL Lac host galaxies and normal ellipticals suggests the active nucleus has surprisingly little effect on the host galaxy. This supports a picture in which all elliptical galaxies harbor black holes which can be actively accreting for some fraction of their lifetime.

BL Lacertae objects are extreme extragalactic sources characterized by the emission of strong and rapidly variable nonthermal radiation over the entire electromagnetic spectrum. Synchrotron emission followed by inverse Compton scattering in a relativistic beaming scenario is generally thought to be the mechanism powering these objects. ...

We report intensity variations of different spectral components in the BL Lac object S5 0716+714 detected during coordinated BeppoSAX and optical observations in 1996 and 1998. The transition between synchrotron and inverse Compton emission has been clearly detected as sharp X-ray spectral breaks at around 2-3 keV on both occasions. Correlated optical and soft X-ray variability was found during the second BeppoSAX pointing when intensive optical monitoring could be arranged. The hard (Compton) component changed by a factor of 2 between the two observations, but remained stable within each exposure. During events of rapid variability S5 0716+714 showed spectral steepening with intensity, a behaviour rarely observed in BL Lacs. We interpret these findings as the probable consequence of a shift of the synchrotron peak emission from the IR/optical band to higher energies, causing the synchrotron tail to push into the soft X-ray band more and more as the source brightens.

We present the results of six BeppoSAX observations of 1ES 2344+514, five of which were taken within a week. 1ES 2344+514, one of the few known TeV BL Lac objects, was detected by the BeppoSAX Narrow Field Instruments between 0.1 to \~50 keV. During the first five closely spaced observations 1ES 2344+514 showed large amplitude luminosity variability, associated with spectacular spectral changes, particularly when compared to the last observation when the source was found to be several times fainter, with a much steeper X-ray spectrum. The energy dependent shape of the lightcurve and the spectral changes imply a large shift (factor of 30 or more in frequency) of the peak of the synchrotron emission. At maximum flux the peak was located at or above 10 keV, making 1ES 2344+514 the second blazar (after MKN501) with the synchrotron peak in the hard X-ray band. The shift, and the corresponding increase in luminosity, might be due to the onset of a second synchrotron component extending from the soft to the hard X-ray band where most of the power is emitted. Rapid variability on a timescale of approximately 5000 seconds has also been detected when the source was brightest.

We have assembled a multi-frequency database by cross-correlating the NVSS catalog of radio sources with the RASSBSC list of soft X-ray sources, obtaining optical magnitude estimates from the Palomar and UK Schmidt surveys as provided by the APM and COSMOS on-line services. By exploiting the nearly unique broad-band properties of High-Energy Peaked (HBL) BL Lacs we have statistically identified a sample of 218 objects that is expected to include about 85% of BL Lacs and that is therefore several times larger than all other published samples of HBLs. Using a subset (155 objects) that is radio flux limited and statistically well-defined we have derived the \vovm distribution and the LogN-LogS of extreme HBLs (fx/fr >= 3E-10 erg/cm2/s/Jansky) down to 3.5 mJy. We find that the LogN-LogS flattens around 20 mJy and that <V/Vm> = 0.42 +/- 0.02. This extends to the radio band earlier results, based on much smaller X-ray selected samples, about the anomalous cosmological observational properties of HBL BL Lacs. A comparison with the expected radio LogN-LogS of all BL Lacs (based on a beaming model) shows that extreme HBLs make up roughly 2% of the BL Lac population, independently of radio flux. This result, together with the flatness of the radio logN-logS at low fluxes, is in contrast with the predictions of a recent model which assumes an anti-correlation between peak frequency and bolometric luminosity. The extreme fx/fr flux ratios and high X-ray fluxes of these BL Lacs makes them good candidate TeV sources, some of the brighter (and closer) ones possibly detectable with the current generation of Cerenkov telescopes.

BeppoSAX observed 8 times the gamma-ray bright blazar PKS 0528+134 in Feb. and Mar. 1997, during a multiwavelength campaign involving EGRET, onboard CGRO, and ground based telescopes. The source was in its faintest X-ray state observed so far, with an unabsorbed [2-10] keV flux of 2.5E-12 erg/cm2/s. The spectrum can be fitted with a power law, with energy index alpha_X=0.49+-0.07 between 0.1 and 10 keV. The high energy instrument PDS detected a flux in the 15-100 keV band which appears disconnected from the lower X-ray energies, most likely due to a contaminating source. We consider our findings in the context of the overall spectral energy distribution and discuss their implications for synchrotron and inverse Compton models.

Active Galactic Nuclei (AGN) are now known to be strong gamma-ray emitters. After briefly describing AGN classification and the main ideas behind unified schemes, I summarize the main properties of blazars (that is BL Lacs and flat-spectrum radio quasars) and their connection with relativistic beaming. Finally, I address the question of why blazars, despite being extreme and very rare objects, are the only AGN detected at high (E > 100 MeV) energies, and touch upon the relevance of TeV astronomy for AGN research.

I briefly summarize the main tenets of unified schemes of BL Lacs and low-luminosity radio galaxies, discussing in particular the evolution of this field after the Como 1988 meeting. I also examine some of the open problems and complications of the simplest scheme. Finally, the question of the existence of two classes of BL Lacs and our related change of perspective in the past few years are also addressed.

We present new BeppoSAX LECS, MECS and PDS observations of five lobe-dominated, broad-line active galactic nuclei selected from the 2-Jy sample of southern radio sources. These include three radio quasars and two broad-line radio galaxies. ROSAT PSPC data, available for all the objects, are also used to better constrain the spectral shape in the soft X-ray band. The collected data cover the energy range 0.1 - 10 keV, reaching ~ 50 keV for one source (Pictor A). The main result from the spectral fits is that all sources have a hard X-ray spectrum with energy index alpha_x ~ 0.75 in the 2 - 10 keV range. This is at variance with the situation at lower energies where these sources exhibit steeper spectra. Spectral breaks ~ 0.5 at 1 - 2 keV characterize in fact the overall X-ray spectra of our objects. The flat, high-energy slope is very similar to that displayed by flat-spectrum/core-dominated quasars, which suggests that the same emission mechanism (most likely inverse Compton) produces the hard X-ray spectra in both classes. Finally, a (weak) thermal component is also present at low energies in the two broad-line radio galaxies included in our study.

We describe a new WWW interface that allows the cross-correlation of an Active Galactic Nuclei catalog with various archives (HST, IUE, EUVE) available at the Multimission Archive at the Space Telescope Science Institute. Details of the catalog and the interface are provided.

We present X-ray spectra obtained with BeppoSAX (Satellite per Astronomia X) of 10 BL Lac objects, selected from the Einstein Medium Sensitivity and Slew Surveys. We find that in about half of the objects a fit in the 0.1-10 keV range with a single power law and free absorption yields values of NH larger than the Galactic ones. In most of these cases, however, broken power law fits with NH fixed at the Galactic values yield an alternative, better description of the data and indicate a steepening of the spectrum with increasing energy. One object (1ES1101-232) is detected up to ~100 keV. Its spectral energy distribution (SED) peaks in the medium energy X-ray band. For each object we compute the peak frequency of the SED from multifrequency data. The spectral indices alpha_x in the 2-10 keV band (Fnu prop to nu^{-alpha_x}) are smaller (i.e. flatter spectrum) for objects with higher peak frequencies. We therefore confirm and extend to higher energies the behavior already known for X-ray selected BL Lac objects in the ROSAT band. We do not find spectral indices smaller than 1; however, the flat distribution of alpha_x and the correlation between alpha_x and peak frequency found from our data suggest that a number of objects may exist, which in the quiescent status have flatter spectrum and peak frequency in the hard X-ray range.

We have undertaken a survey for blazars by correlating the ROSAT WGACAT database with publicly available radio catalogs, restricting our candidate list to serendipitous flat-spectrum sources (alpha_r <= 0.7). We discuss here our survey methods, identification procedure and first results. Our survey is found to be ~ 95% efficient at finding blazars, a figure which is comparable to or greater than that achieved by other radio and X-ray survey techniques. DXRBS provides a much more uniform coverage of the parameter space occupied by blazars than any previous survey. Particularly important is the identification of a large population of flat-spectrum radio quasars with ratios of X-ray to radio luminosity >~ 10^-6 (alpha_rx <~ 0.78) and of many low-luminosity flat-spectrum radio quasars. Moreover, DXRBS fills in the region of parameter space between X-ray selected and radio-selected samples of BL Lacs.

The results of a series of five \sax observations of the TeV BL Lac object 1ES2344+514 are briefly presented. Large amplitude luminosity variability, associated to impressive spectral changes in the hard X-rays, have been found. The shape of the lightcurve depends on energy, with the flare starting and ending in the hard band, but with maximum intensity possibly reached earlier in the soft X-rays. The luminosity and spectral changes may be due to a shift of the peak of the synchrotron emission from the soft X-rays to the hard X-ray band similar to that detected during \sax observations of MKN 501.

We present new BeppoSAX LECS, MECS, and PDS observations of five lobe-dominated, broad-line active galactic nuclei selected from the 2 Jy sample of southern radio sources. These include three radio quasars and two broad-line radio galaxies. ROSAT PSPC data, available for all the objects, are also used to better constrain the spectral shape in the soft X-ray band. The collected data cover the 0.1 - 10 keV energy range, reaching 40 keV for one source. Detailed spectral fitting shows that all sources have a flat hard X-ray spectrum with energy index alpha_x ~ 0.75 in the 2 - 10 keV energy range. This is a new result, which is at variance with the situation at lower energies where these sources exhibit steeper spectra. Spectral breaks ~0.5 at 1 - 2 keV characterize the overall X-ray spectra of our objects. The flat, high-energy slope is very similar to that displayed by flat-spectrum/core-dominated quasars, which suggests that the same emission mechanism (most likely inverse Compton) produces the hard X-ray spectra in both classes. Contrary to the optical evidence for some of our sources, no absorption above the Galactic value is found in our sample. Finally, a (weak) thermal component is also present at low energies in the two broad-line radio galaxies included in our study.

We have undertaken a survey of archived, pointed ROSAT PSPC data for blazars by correlating the ROSAT WGACAT database with several publicly available radio catalogs, restricting our candidate list to serendipitous flat radio spectrum sources (alpha_r <= 0.70). Here we discuss our survey methods, identification procedure and first results. Our survey is found to be ~ 95% efficient at finding flat-spectrum radio-loud quasars (FSRQs, 59 of our first 85 IDs) and BL Lacertae objects (22 of our first 85 IDs), a figure which is comparable to or greater than that achieved by other radio and X-ray survey techniques. The identifications presented here show that all previous samples of blazars (even when taken together) did not representatively survey the blazar population, missing critical regions of (L_X,L_R) parameter space within which large fractions of the blazar population lie. Particularly important is the identification of a large population of FSRQs (>~ 25% of DXRBS FSRQs) with ratios of X-ray to radio luminosity >~ 10^-6 (alpha_rx <~ 0.78). In addition, due to our greater sensitivity, DXRBS has already more than doubled the number of FSRQs in complete samples with 5 GHz (radio) luminosities between 10^31.5 and 10^33.5 erg/s/Hz and fills in the region of parameter space between X-ray selected and radio-selected samples of BL Lacs. DXRBS is the very first sample to contain statistically significant numbers of blazars at low luminosities, approaching what should be the lower end of the FSRQ luminosity function.

We present the first results of BeppoSAX observations of a small sample of low brightness FRI radio galaxies. The flux of all the targets is consistent with a thermal spectrum, as due to the presence of hot intracluster gas or galactic corona. Moreover in three sources a non thermal absorbed spectrum can be present in the MECS spectrum at energies larger than 7 keV, while for a fourth object a high energy flux has been detected in the PDS instrument at energies larger than 15 keV. This component could be related to the inner AGN surrounded by an obscuring torus.

We report on the BeppoSAX observations of the gamma-ray blazar PKS 0528+134 performed in Feb and Mar 1997, during a multiwavelength campaign involving EGRET and ground based telescopes. The source was in a faint and hard state, with energy spectral index alpha=0.48+-0.04 between 0.1 and 10 keV, and [2-10] keV flux of 2.7E-12 erg/cm2/s. No significant variability was observed. The source was detected in the 20-120 keV band by the PDS, with a flux lying slightly above the extrapolation from lower X-ray energies. Comparing this low state with previous higher states of the source, there is an indication that the X-ray spectrum hardens and the gamma-ray spectrum steepens when the source is fainter.

We have used the WGA catalog of ROSAT PSPC X-ray sources to study the X-ray spectrum of about 500 quasars in the redshift interval 0.1--4.1, detected with a signal to noise better than 7. We have parameterized the PSPC spectrum in terms of two `effective energy spectral indices', alpha_Soft (0.1-0.8 keV), and alpha_Hard (0.4-2.4 keV), which allows for the different Galactic N_H along the quasars line of sight. We have used these data to explore the questions raised by the initial PSPC high redshift quasar studies, and in particular the occurrence of low X-ray energy cut-offs in high redshift radio-loud quasars. We have also studied the emission spectra of a large sample of radio-loud and radio-quiet quasars and studied their differences. We find that low energy X-ray cut-offs are more commonly (and perhaps exclusively) found in radio-loud quasars. Therefore the low energy X-ray cut-offs are physically associated with the quasars, and not with intervening systems, since those would affect radio-quiet and radio-loud equally. We suggest that photoelectric absorption is a likely origin of the these cut-offs.

We have selected quasars with X-ray colors suggestive of a low energy cut-off, from the ROSAT PSPC pointed archive. We examine the radio and optical properties of these 13 quasars. Five out of the seven quasars with good optical spectra show associated optical absorption lines, with two having high delta-v candidate systems. Two other cut-off quasars show reddening associated with the quasar. We conclude that absorption is highly likely to be the cause of the X-ray cut-offs, and that the absorbing material associated with the quasars, not intervening along the line-of-sight. The suggestion that Gigahertz Peaked Sources are associated with X-ray cut-offs remains unclear with this expanded sample.

Active Galactic Nuclei (AGN) are now known to be strong gamma-ray emitters. After briefly describing the different classes of AGN and the basic tenets of unified schemes, I discuss the role of blazars (that is BL Lacs and flat-spectrum radio quasars) as gamma-ray sources. The main properties of blazars and their connection with relativistic beaming are then summarized. Finally, I address the question of why blazars, despite being extreme and very rare objects, are the only AGN detected at very high (E > 100 MeV) energies, and touch upon the relevance of TeV astronomy for AGN research.

After briefly summarizing the main tenets of unified schemes of Active Galactic Nuclei, I review some recent results in the field of unification of radio-loud sources, both for the low-luminosity (BL Lacs and Fanaroff-Riley type I radio galaxies) and high-luminosity (radio quasars and Fanaroff-Riley type II radio galaxies) populations.

We present evidence in favour of a link between the luminosity radiatively dissipated in the central engine of radio-loud Active Galactic Nuclei and the kinetic power in their jets. This piece of evidence is based on the relation we find between the luminosity in broad emission lines and the kinetic power in pc-scale radio jets, for a sample of radio-loud quasars for which suitable data are available in the literature. We find that the ionizing luminosity and the kinetic one are of the same order of magnitude, suggesting that the processes responsible for them are somehow related. A strong magnetic field in equipartition with the radiation field could be responsible for regulating both processes. BL Lac objects seem to follow a similar behaviour, but with comparatively fainter broad line emission.

Active Galactic Nuclei (AGN) are ideal sources for multi-wavelength studies as their emission can cover almost 20 orders of magnitude in frequency from the radio to the gamma-ray band. After reviewing their basic properties, I will assess how well we know the multifrequency spectra of AGN as a class. I will then briefly illustrate how currently available and forthcoming sky surveys will help in addressing some of the open questions of AGN studies. Finally, an analysis of the problem of the missing Type 2 QSO will exemplify the dangers of monochromatic sky surveys for AGN.

Combining information on absorbing material in AGN from X-ray and the UV creates a powerful investigative tool. Here we give examples from both low and high redshift.

We study the X-ray spectra of 114 flat-spectrum radio quasars (FSRQ) using the hardness ratios as given in the WGA catalogue of ROSAT sources. This sample includes all WGA FSRQ with high-quality data and comprises about 20 per cent of presently known such objects, which makes this the largest FSRQ sample ever studied in the X-ray band. We find that FSRQ have a distribution of energy spectral indices ranging between 0 and 3 with a mean value alpha_x ~ 1. This is consistent with that of low-energy cutoff BL Lacs (LBL; alpha_x ~ 1.1), generally found in radio surveys, but significantly different from that of high-energy cutoff BL Lacs (HBL), normally selected in the X-ray band, which display steeper X-ray spectra (alpha_x ~ 1.5). The shape of the optical-to-X-ray continuum is concave (that is alpha_x < alpha_ox) for the majority of FSRQ, as found for LBL, supporting a dominance of inverse Compton emission in the X-ray band in most objects. Our results are at odds with previous studies of the X-ray spectra of FSRQ, which were however plagued by low spectral resolution and/or small number statistics and selection effects, and have important implications for the proposed connections between FSRQ and BL Lacs.

We study the X-ray spectra of 85 BL Lacertae objects using the hardness ratios as given in the WGA catalogue of {\it ROSAT} sources. This sample includes all WGA BL Lacs with high-quality data and comprises about 50 per cent of presently known BL Lacs. We find that BL Lacs have energy power-law spectral indices between 0 and 3 with a mean value $\alpha_{\rm x}\sim 1.4$. Significant differences, however, are present between high-energy cutoff BL Lacs (HBLs), normally selected in the X-ray band, and low-energy cutoff BL Lacs (LBLs), generally found in radio surveys. HBLs have steeper X-ray spectral slopes ($\alpha_{\rm x}\sim 1.5$) well correlated with $\alpha_{\rm ox}$, and anti-correlated with the X-ray-to-radio flux ratio and cutoff frequency, with convex overall broad-band spectra. LBLs, on the other hand, have flatter X-ray spectra ($\alpha_{\rm x}\sim 1.1$) and concave optical-X-ray continuum. We interpret these results in terms of different mechanisms being responsible for the X-ray emission in the two classes, namely synchrotron and inverse Compton for HBLs and LBLs respectively. The observed differences are consistent with the hypothesis that HBLs and LBLs are powered by essentially the same non-thermal engines differing mainly in their synchrotron cutoff energy.

We present a catalogue of 233 BL Lacertae objects compiled through an extensive bibliographic search updated to mid-1995. A large fraction of the sources listed in the catalogue belongs to well-defined samples and can be used for statistical purposes. A smaller fraction consists of miscellaneous (but confirmed) BL Lacs and of objects classified as BL Lac candidates. We discuss the selection criteria of the different samples, report the discovery of two previously unnoticed BL Lacs in the Palomar-Green survey, and comment on the possible association of some of the still unidentified high galactic latitude gamma-ray (EGRET) sources with BL Lacs. Some statistical properties of the catalogue are also briefly discussed.

The appearance of active galactic nuclei (AGN) depends so strongly on orientation that our current classification schemes are dominated by random pointing directions instead of more interesting physical properties. Light from the centers of many AGN is obscured by optically thick circumnuclear matter and in radio-loud AGN, bipolar jets emanating from the nucleus emit light that is relativistically beamed along the jet axes. Understanding the origin and magnitude of radiation anisotropies in AGN allows us to unify different classes of AGN; that is, to identify each single, underlying AGN type that gives rise to different classes through different orientations. This review describes the unification of radio-loud AGN, which include radio galaxies, quasars, and blazars. We describe the classification and properties of AGN and summarize the evidence for anisotropic emission. We outline the two most plausible unified schemes for radio-loud AGN, one linking quasars and luminous radio galaxies and another linking BL~Lac objects and less luminous radio galaxies. Using the formalism appropriate to samples biased by relativistic beaming, we show the population statistics for two schemes are in accordance with available data. We analyze the possible connections between low- and high-luminosity radio-loud AGN. We review potential difficulties with unification and conclude that none currently constitutes a serious problem. We discuss likely complications to unified schemes that are suggested by realistic physical considerations; these will be important to consider when more comprehensive data for larger complete samples become available. We conclude with a list of the ten questions we believe are the most pressing in this field.

We explain the properties of X-ray selected BL Lacertae objects, under the assumption that they constitute the small minority of the BL Lac population with energy cutoff located in the UV/X-ray band, as suggested by their multifrequency spectra. In particular, we derive their X-ray luminosity function, log N-log S, and radio flux distribution starting from radio selected BL Lacs using a method that does not require any free parameters. Our hypothesis is at variance with the most common explanation for the existence of two classes of BL Lacs, namely that X-ray selected objects are viewed significantly off the beaming axis while radio-selected ones have their jets more aligned with the line of sight. In this case, in fact, X-ray selected sources should be intrinsically more numerous than objects selected in the radio band at a fixed value of some isotropic luminosity such as that of the diffuse radio emission. While also accounting for most of the properties of BL Lacs detected in current X-ray surveys, this latter scenario would predict their radio emission to be dominated by the extended (unbeamed) component, in contrast with observations. We also address the question of the negative evolution apparently exhibited by X-ray selected BL Lacs, showing that neither recognition problems nor a break in the X-ray spectrum are likely to be responsible for it. Our hypothesis explains in a simple way the bimodal distribution of BL Lacs in the alpha_ro - alpha_ ox diagram. It also makes specific predictions for the X-ray, radio, and optical properties of BL Lacs which will be tested by future deep X-ray surveys.