This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end.

We selected a sample of 28 Type 1 AGNs for which a black hole mass has been inferred using the reverberation mapping technique and single epoch scaling relations. All 28 sources show clear evidence of the "Big Blue Bump" in the optical-UV band whose emission is produced by an accretion disk (AD) around a supermassive black hole. We fitted the spectrum of these sources with the relativistic thin AD model KERRBB in order to infer the black hole masses and compared them with those from Reverberation mapping and Single epoch methods, discussing the possible uncertainties linked to such a model by quantifying their weight on our results. We find that for the majority of the sources, KERRBB is a good description of the AD emission for a wide wavelength range. The overall uncertainty on the black hole mass estimated through the disk fitting procedure is $\sim 0.45$ dex (which includes the uncertainty on fitting parameters such as e.g. spin and viewing angle), comparable to the systematic uncertainty of reverberation mapping and single epoch methods; however, such an uncertainty can be $\leq 0.3$ dex if one of the parameters of the fit is well constrained. Although all of the estimates are affected by large uncertainties, the masses inferred using the three methods are compatible if the dimensionless scale factor $f$ (linked to the unknown kinematics and geometry of the Broad Line Region) is assumed to be larger than one. For the majority of the sources, the comparison between the results coming from the three methods favors small spin values. To check the goodness of the KERRBB results, we compared them with those inferred with other models, such as AGNSED, a model that also accounts for the emission originating from an X-ray corona: using two sources with a good data coverage in the X band, we find that the masses estimated with the two models differ at most by a factor of $\sim 0.2$ dex.

We present the discovery of PSO J030947.49+271757.31, the radio brightest (23.7 mJy at 1.4 GHz) active galactic nucleus (AGN) at z>6.0. It was selected by cross-matching the NRAO VLA Sky Survey and the Panoramic Survey Telescope and Rapid Response System PS1 databases and its high-z nature was confirmed by a dedicated spectroscopic observation at the Large Binocular Telescope. A pointed Neil Gehrels $Swift$ Observatory XRT observation allowed us to measure a flux of $\sim$3.4$\times$10$^{-14}$ erg s$^{-1}$ cm$^{-2}$ in the [0.5-10] keV energy band, which also makes this object the X-ray brightest AGN ever observed at z>6.0. Its flat radio spectrum ($\alpha_{\nu r}$<0.5), very high radio loudness (R>10$^3$), and strong X-ray emission, compared to the optical, support the hypothesis of the blazar nature of this source. %i.e. a radio-loud (RL) AGN with the relativistic jet pointed toward us. Assuming that this is the only blazar at this redshift in the surveyed area of sky, we derive a space density of blazars at z$\sim$6 and with M$_{1450 \mbox{\AA}}$ < -25.1 of 5.5$^{+11.2}_{-4.6}$$\times$10$^{-3}$ Gpc$^{-3}$. From this number, and assuming a reasonable value of the bulk velocity of the jet ($\Gamma$=10), we can also infer a space density of the entire radio-loud AGN population at z$\sim$6 with the same optical/UV absolute magnitude of 1.10$^{+2.53}_{-0.91}$ Gpc$^{-3}$. Larger samples of blazars will be necessary to better constrain these estimates.

Accretion onto black holes at rates above the Eddington limit has long been discussed in the context of supermassive black hole (SMBH) formation and evolution, providing a possible explanation for the presence of massive quasars at high redshifts (z$\gtrsim$7), as well as having implications for SMBH growth at later epochs. However, it is currently unclear whether such `super-Eddington' accretion occurs in SMBHs at all, how common it is, or whether every SMBH may experience it. In this work, we investigate the observational consequences of a simplistic model for super-Eddington accretion flows -- an optically thick, geometrically thin accretion disc (AD) where the inner-most parts experience severe photon-trapping, which is enhanced with increased accretion rate. The resulting spectral energy distributions (SEDs) show a dramatic lack of rest-frame UV, or even optical, photons. Using a grid of model SEDs spanning a wide range in parameter space (including SMBH mass and accretion rate), we find that large optical quasar surveys (such as SDSS) may be missing most of these luminous systems. We then propose a set of colour selection criteria across optical and infra-red colour spaces designed to select super-Eddington SEDs in both wide-field surveys (e.g., using SDSS, 2MASS and WISE) and deep & narrow-field surveys (e.g., COSMOS). The proposed selection criteria are a necessary first step in establishing the relevance of advection-affected super-Eddington accretion onto SMBHs at early cosmic epochs.

While massive black holes (MBHs) are known to inhabit all massive galaxies, their ubiquitous presence in dwarf galaxies has not been confirmed yet, with only a limited number of sources detected so far. Recently, some studies proposed infrared emission as an alternative way to identify MBHs in dwarfs, based on a similar approach usually applied to quasars. In this study, by accurately combining optical and infrared data taking into account resolution effects and source overlapping, we investigate in detail the possible limitations of this approach with current ground-based facilities, finding a quite low ($\sim$0.4 per cent) fraction of active MBH in dwarfs that are luminous in mid-infrared, consistent with several previous results. Our results suggest that the infrared selection is strongly affected by several limitations that make the identification of MBHs in dwarf galaxies currently prohibitive, especially because of the very poor resolution compared to optical surveys, and the likely contamination by nearby sources, although we find a few good candidates worth further follow-ups. Optical, X-ray and radio observations, therefore, still represent the most secure way to search for MBH in dwarfs.

We present the discovery and properties of DESJ014132.4-542749.9 (DES0141-54), a new powerful radio-loud active galactic nucleus (AGN) in the early Universe (z=5.0). It was discovered by cross-matching the first data release of the Dark Energy Survey (DES DR1) with the Sidney University Molonglo Survey (SUMSS) radio catalog at 0.843 GHz. This object is the first radio-loud AGN at high redshift discovered in the DES. The radio properties of DES0141-54, namely its very large radio-loudness (R>10$^{4}$), the high radio luminosity (L$_{0.8 GHz}$=1.73$\times$10$^{28}$ W Hz$^{-1}$), and the flatness of the radio spectrum ($\alpha$=0.35) up to very high frequencies (120 GHz in the source's rest frame), classify this object as a blazar, meaning, a radio-loud AGN observed along the relativistic jet axis. However, the X--ray luminosity of DESJ0141-54 is much lower compared to those of the high redshift (z$\geq$4.5) blazars discovered so far. Moreover its X-ray-to-radio luminosity ratio (log($\frac{L_{[0.5-10]keV}}{L_{1.4GHz}}$)=9.96$\pm$0.30 Hz) is small also when compared to lower redshift blazars: only 2\% of the low-z population has a similar ratio. By modeling the spectral energy distribution we found that this peculiar X--ray weakness and the powerful radio emission could be related to a particularly high value of the magnetic field. Finally, the mass of the central black hole is relatively small (M$_{BH}$ = 3-8 $\times$10$^8$ M$_{\odot}$) compared to other confirmed blazars at similar redshift, making DES0141-54 the radio-loud AGN that host the smallest supermassive black hole ever discovered at z$\geq$5.

We observed with the NuSTAR satellite 3 blazars at z>2, detected in gamma-rays by Fermi/LAT and in the soft X-rays, but not yet observed above 10 keV. The flux and slope of their X-ray continuum, together with Fermi/LAT data allows us to estimate their total electromagnetic output and peak frequency. For some of them we can study the source in different states, and investigate the main cause of the observed different spectral energy distribution. We then collected all blazars at redshift greater than 2 observed by NuSTAR, and confirm that these hard and luminous X-ray blazars are among the most powerful persistent sources in the Universe. We confirm the relation between the jet power and the disk luminosity, extending it at the high energy end.

We investigate the properties of the most distant quasars ULASJ134208.10+092838.61 ($z = 7.54$), ULASJ112001.48+064124.3 ($z = 7.08$) and DELSJ003836.10-152723.6 ($z = 7.02$) studying their Optical-UV emission that shows clear evidence of the presence of an accretion disk. We model such emission applying the relativistic disk models KERRBB and SLIMBH for which we have derived some analytical approximations to describe the observed emission as a function of the black hole mass, accretion rate, spin and the viewing angle. We found that: 1] our black hole mass estimates are compatible with the ones found using the virial argument but with a smaller uncertainty; 2] assuming that the virial argument is a reliable method to have a black hole mass measurement (with no systematic uncertainties involved), we found an upper limit for the black hole spin of the three sources: very high spin values are ruled out; 3] our Eddington ratio estimates are smaller than those found in previous studies by a factor $\sim 2$: all sources are found to be sub-Eddington. Using our results, we explore the parameter space (efficiency, accretion rate) to describe the possible evolution of the black hole assuming a $\sim 10^{2-4} M_{\odot}$ seed: if the black hole in these sources formed at redshift $z = 10 - 20$, we found that the accretion has to proceed at the Eddington rate with a radiative efficiency $\eta \sim 0.1$ in order to reach the observed masses in less than $\sim 0.7$ Gyr.

Volonteri et al. (2011) found that the number of radio-loud quasars above redshift 4 calculated from the luminosity function (based upon Swift/BAT observations) is much smaller than the number estimated from the known high-redshift beamed sources, blazars, assuming that for every beamed source with a Lorentz factor of $\Gamma$, statistically $2 \Gamma^2$ non-beamed sources should exist. To explain the missing misaligned (non-beamed) population of high-redshift sources, they proposed various explanations, involving heavy optical obscuration and significantly different Lorentz factors at early cosmological epochs. Our EVN observations targeting high-redshift ($z>4$) blazar candidates revealed 3 sources not showing relativistic beaming, but rather kpc-scale double structures. These three sources have significant radio emission resolved out with the EVN, while they are compact on $\sim 5-10$ arcsec scale. Our dual-frequency ($1.5$ and $5$ GHz) e-MERLIN observations of these three sources revealed a rich morphology, bending jets, and hot spots with possible sites of interaction between the jets and the surrounding medium at intermediate scales.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

We study a sample of Narrow-line Seyfert 1 galaxies (NLS1) in their optical and radio features, in order to understand the differences between their radio silent, radio-loud and radio-quiet subclasses. We first show that the different redshift and mass distributions of radio-loud and -quiet NLS1s could be ascribed to observational biases. We then present a geometrical model according to which most of the different observational features of radio-loud and radio-quiet NLS1s are ascribed to the orientation of an intrinsically structured jet. We estimate the fraction of intrinsically jetted sources among NLS1s that justifies the observed radio-detected population. Noticeably, under the assumptions of the geometrical model, we derive a fraction of jetted sources significantly larger than in standard AGN.

Virial-based methods for estimating active supermassive black hole masses are now commonly used on extremely large spectroscopic quasar catalogues. Most spectral analyses, though, do not pay enough attention to a detailed continuum decomposition. To understand how this affects virial mass estimate results, we test the influence of host galaxy light on them, along with Balmer continuum component. A detailed fit with the new spectroscopic analysis software QSFit demonstrated that the presence or absence of continuum components do not affect significantly the virial-based results for our sample. Taking or not in consideration a host galaxy component, instead, affects the emission line fitting in a more pronounced way at lower redshifts, where in fact we observe dimmer quasars and more visible host galaxies.

We investigate the global properties of the radiation emitted by the accretion disk around Kerr black holes. Using the Kerr blackbody (KERRBB) numerical model, we build an analytic approximation of the disk emission features focusing on the pattern of the produced radiation as a function of the black hole spin, mass, accretion rate and viewing angle. The assumption of having a geometrically thin disk limits our analysis to systems emitting below ~0.3 of the Eddington luminosity. We apply this analytical model to four blazars (whose jets are pointing at us) at high redshift that show clear signatures of disk emission. For them, we derive the black hole masses as a function of spin. If these jetted sources are powered by the black hole rotation, they must have high spin values, further constraining their masses.

We present QSFit, a new software package to automatically perform the analysis of Active Galactic Nuclei (AGN) optical spectra. The software provides luminosity estimates for the AGN continuum, the Balmer continuum, both optical and UV iron blended complex, host galaxy and emission lines, as well as width, velocity offset and equivalent width of 20 emission lines. Improving on a number of previous studies on AGN spectral analysis, QSFit fits all the components simultaneously, using an AGN continuum model which extends over the entire available spectrum, and is thus a probe of the actual AGN continuum whose estimates are scarcely influenced by localized features (e.g. emission lines) in the spectrum. We used QSFit to analyze 71,251 optical spectra of Type 1 AGN at z < 2 (obtained by the Sloan Digital Sky Survey, SDSS) and to produce a publicly available catalog of AGN spectral properties. Such catalog allowed us (for the first time) to estimate the AGN continuum slope and the Balmer continuum luminosity on a very large sample, and to show that there is no evident correlation between these quantities the redshift. All data in the catalog, the plots with best fitting model and residuals, and the IDL code we used to perform the analysis, are available on a dedicated website. The whole fitting process is customizable for specific needs, and can be extended to analyze spectra from other data sources. The ultimate purpose of QSFit is to allow astronomers to run standardized recipes to analyze the AGN data, in a simple, replicable and shareable way.

A brief Chandra observation of the ultraluminous quasar, SDSS J010013.02+280225.8 at redshift 6.326, showed it to be a relatively bright, soft X-ray source with a count rate of about 1 ct/ks. In this paper we present results for the quasar from a 65ks XMM-Newton observation, which well constrains its spectral shape. The quasar is clearly detected with a total of $\sim$ 460 net counts in the 0.2-10 keV band. The spectrum is characterised by a simple power-law model with photon index of $\Gamma = 2.30^{+0.10}_{-0.10}$, and the intrinsic 2-10 keV luminosity is $3.14\times10^{45}$ erg $\text{s}^{-1}$. The 1 $\sigma$ upper limit to any intrinsic absorption column density is $N_{H} = 6.07\times 10^{22} {\text{cm}}^{-2}$. No significant iron emission lines were detected. We derive the X-ray-to-optical flux ratio $\alpha_{\text{ox}}$ of $-1.74\pm$0.01, consistent with the values found in other quasars of comparable ultraviolet luminosity. We did not detect significant flux variations either in the XMM-Newton exposure or between XMM-Newton and XMM-Newton observations, which are separated by $\sim$ 8 months. The X-ray observation enables the bolometric luminosity to be calculated after modelling the spectral energy distribution: the accretion rate is found to be sub-Eddington.

In an effort to understand the cause of the apparent depletion in the number density of radio-loud AGNs at $z>3$, this work investigates the viability of the so-called Cosmic Microwave Background (CMB) quenching mechanism of intrinsically jetted, high-z AGNs, whereby Inverse Compton scattering of CMB photons off electrons within the extended lobes results in a substantial dimming of the lobe synchrotron emission at GHz frequencies, while simultaneously boosting their diffuse X-ray signal. We focus on five $z>3.5$ radio galaxies that have sufficiently deep Chandra exposure (> 50 ks) to warrant a meaningful investigation of any extended X-ray emission. For those objects with evidence for statistically significant extended X-ray lobes (4C 41.17 and 4C 03.24), we combine the Chandra measurements with literature data at lower frequencies to assemble the systems' Spectral Energy Distributions (SEDs), and utilize state-of-the-art SED modelling (Ghisellini et al. 2015) -- including emission from the disk, torus, jet, hotspots, and lobes -- to infer their physical parameters. For both radio galaxies, the magnetic energy density in the hotspots is found to exceed the energy density in CMB photons, wheres the opposite is true for the lobes. This implies that any extended synchrotron emission likely originates from the hotspots themselves, rather than the lobes. Conversely, Inverse Compton scattering of CMB photons dominates the extended X-ray emission from the lobes, which are effectively "radio-quenched". As a result, CMB quenching is effective in these systems in spite of the fact that the observed X-ray to radio luminosity ratio does not bear the signature $(1+z)^4$ dependence of the CMB energy density.

Observational searches for faint active nuclei at $z > 6$ have been extremely elusive, with a few candidates whose high-$z$ nature is still to be confirmed. Interpreting this lack of detections is crucial to improve our understanding of high-$z$ supermassive black holes (SMBHs) formation and growth. In this work, we present a model for the emission of accreting BHs in the X-ray band, taking into account super-Eddington accretion, which can be very common in gas-rich systems at high-$z$. We compute the spectral energy distribution for a sample of active galaxies simulated in a cosmological context, which represent the progenitors of a $z \sim 6$ SMBH with $M_{\rm BH} \sim 10^9 \, M_\odot$. We find an average Compton thick fraction of $\sim 45\%$ and large typical column densities ($N_H \gtrsim 10^{23} \rm \, cm^2$). However, faint progenitors are still luminous enough to be detected in the X-ray band of current surveys. Even accounting for a maximum obscuration effect, the number of detectable BHs is reduced at most by a factor 2. In our simulated sample, observations of faint quasars are mainly limited by their very low active fraction ($f_{\rm act} \sim 1 \%$), which is the result of short, super-critical growth episodes. We suggest that to detect high-$z$ SMBHs progenitors, large area surveys with shallower sensitivities, such as Cosmos Legacy and XMM-LSS+XXL, are to be preferred with respect to deep surveys probing smaller fields, such as CDF-S.

Hard X-ray observations are crucial to study the non-thermal jet emission from high-redshift, powerful blazars. We observed two bright z>2 flat spectrum radio quasars (FSRQs) in hard X-rays to explore the details of their relativistic jets and their possible variability. S5 0014+81 (at z=3.366) and B0222+185 (at z=2.690) have been observed twice by the Nuclear Spectroscopic Telescope Array (NuSTAR) simultaneously with Swift/XRT, showing different variability behaviours. We found that NuSTAR is instrumental to explore the variability of powerful high-redshift blazars, even when no gamma-ray emission is detected. The two sources have proven to have respectively the most luminous accretion disk and the most powerful jet among known blazars. Thanks to these properties, they are located at the extreme end of the jet-accretion disk relation previously found for gamma-ray detected blazars, to which they are consistent.

At redshift larger than 3 there is a disagreement between the number of blazars (whose jet is pointing at us) and the number of expected parents (whose jet is pointing elsewhere). Now we strengthen this claim because (i) the number of blazars identified within the SDSS+FIRST survey footprint increased, demanding a more numerous parent population, and (ii) the detected blazars have a radio flux large enough to be above the FIRST flux limit even if the jet is slightly misaligned. The foreseen number of these slightly misaligned jets, in principle detectable, is much larger than the radio-detected sources in the FIRST+SDSS survey (at redshift larger than 4). This argument is independent of the presence of an isotropic radio component, such as the hot spot or the radio lobe, and does not depend on the bulk Lorentz factor Gamma. We propose a scenario that ascribes the lack of slightly misaligned sources to an over-obscuration of the nucleus by a "bubble" of dust, possibly typical of the first high-redshift quasars.

The very existence of more than a dozen of high-redshift (z>4) blazars indicates that a much larger population of misaligned powerful jetted AGN was already in place when the Universe was <1.5 Gyr old. Such parent population proved to be very elusive, and escaped direct detection in radio surveys so far. High redshift blazars themselves seem to be failing in producing extended radio-lobes, raising questions about the connection between such class and the vaster population of radio-galaxies. We show that the interaction of the jet electrons with the intense cosmic microwave background (CMB) radiation explains the lack of extended radio emission in high redshift blazars and in their parent population, helping to explain the apparently missing misaligned counterparts of high redshift blazars. On the other hand, the emission from the more compact and more magnetised hot spots are less affected by the enhanced CMB energy density. By modelling the spectral energy distribution of blazar lobes and hot spots we find that most of them should be detectable by low frequency deep radio observations, e.g., by LOw-Frequency ARray for radio astronomy (LOFAR) and by relatively deep X-ray observations with good angular resolution, e.g., by the Chandra satellite. At high redshifts, the emission of a misaligned relativistic jet, being de-beamed, is missed by current large sky area surveys. The isotropic flux produced in the hot spots can be below ~1 mJy and the isotropic lobe radio emission is quenched by the CMB cooling. Consequently, even sources with very powerful jets can go undetected in current radio surveys, and misclassified as radio-quiet AGNs.

The most powerful blazars are the flat spectrum radio quasars whose emission is dominated by a Compton component peaking between a few hundred keV and a few hundred MeV. We selected two bright blazars, PKS 2149-306 at redshift z=2.345 and S5 0836+710 at z=2.172, in order to observe them in the hard X-ray band with the NuSTAR satellite. In this band the Compton component is rapidly rising almost up to the peak of the emission. Simultaneous soft-X-rays and UV-optical observations were performed with the Swift satellite, while near-infrared (NIR) data were obtained with the REM telescope. To study their variability, we repeated these observations for both sources on a timescale of a few months. While no fast variability was detected during a single observation, both sources were found to be variable in the X-ray band, up to 50%, between the two observations, with larger variability at higher energies. No variability was detected in the optical/NIR band. These data together with Fermi-LAT, WISE and other literature data are then used to study the overall spectral energy distributions (SEDs) of these blazars. Although the jet non-thermal emission dominates the SED, it leaves the UV band unhidden, allowing us to detect the thermal emission of the disc and to estimate the mass of the black hole. The non-thermal emission is well reproduced by a one-zone leptonic model. The non-thermal radiative processes are synchrotron, self-Compton and external Compton using seed photons from both the broad-line region (BLR) and the torus. We find that our data are better reproduced if we assume that the location of the dissipation region of the jet, R_diss, is in-between the torus, (at R_torus), and the BLR (R_torus>R_diss>R_BLR). The observed variability is explained by changing a minimum number of model parameters by a very small amount.

High energy observations of extreme BL Lac objects, such as 1ES 0229+200 or 1ES 0347-121, recently focused interest both for blazar and jet physics and for the implication on the extragalactic background light and intergalactic magnetic field estimate. However, the number of these extreme highly peaked BL Lac objects (EHBL) is still rather small. Aiming at increase their number, we selected a group of EHBL candidates starting from the BL Lac sample of Plotkin et al. (2011), considering those undetected (or only barely detected) by the Large Area Telescope onboard Fermi and characterized by a high X-ray vs. radio flux ratio. We assembled the multi-wavelength spectral energy distribution of the resulting 9 sources, profiting of publicly available archival observations performed by the Swift, Galex and Fermi satellites, confirming their nature. Through a simple one-zone synchrotron self-Compton model we estimate the expected VHE flux, finding that in the majority of cases it is within the reach of present generation of Cherenkov arrays or of the forthcoming Cherenkov Telescope Array (CTA).

The radio-loud quasar SDSS J013127.34-032100.1at a redshift z=5.18 is one of the most distant radio-loud objects. The radio to optical flux ratio (i.e. the radio-loudness) of the source is large, making it a promising blazar candidate. Its overall spectral energy distribution, completed by the X-ray flux and spectral slope derived through Target of Opportunity Swift/XRT observations, is interpreted by a non-thermal jet plus an accretion disc and molecular torus model. We estimate that its black hole mass is (1.1+-0.2)1e10 Msun. for an accretion efficiency eta=0.08, scaling roughly linearly with eta. Although there is a factor ~2 of systematic uncertainty, this black hole mass is the largest found at these redshifts in a radio loud object. We derive a viewing angle between 3 and 5 degrees. This implies that there must be other (hundreds) sources with the same black hole mass of SDSS J013127.34-032100.1, but whose jets are pointing away from Earth. We discuss the problems posed by the existence of such large black hole masses at such redshifts, especially in jetted quasars. In fact, if they are associated to rapidly spinning black holes, the accretion efficiency is high, implying a slower pace of black hole growth with respect to radio-quiet quasars.

We have selected SDSS J222032.50+002537.5 and SDSS J142048.01+120545.9 as best blazar candidates out of a complete sample of extremely radio-loud quasars at z>4, with highly massive black holes. We observed them and a third serendipitous candidate with similar features (PMN J2134-0419) in the X-rays with the Swift/XRT telescope, to confirm their blazar nature. We observed strong and hard X-ray fluxes (i.e. $\alpha_X<0.6$, where $F(\nu)\propto\nu^{-\alpha_X}$ in the 0.3-10keV observed energy range, ~1-40keV rest frame) in all three cases. This allowed us to classify our candidates as real blazars, being characterized by large Lorentz factors (~13) and very small viewing angles (~3deg). All three sources have black hole masses exceeding 10^9Msun and their classification provides intriguing constraints on supermassive black hole formation and evolution models. We confirm our earlier suggestion that there are different formation epochs of extremely massive black holes hosted in jetted (z~4) and non-jetted systems (z~2.5).

Theoretical models for the production of relativistic jets from active galactic nuclei predict that jet power arises from the spin and mass of the central black hole, as well as the magnetic field near the event horizon. The physical mechanism mechanism underlying the contribution from the magnetic field is the torque exerted on the rotating black hole by the field amplified by the accreting material. If the squared magnetic field is proportional to the accretion rate, then there will be a correlation between jet power and accretion luminosity. There is evidence for such a correlation, but inadequate knowledge of the accretion luminosity of the limited and inhomogeneous used samples prevented a firm conclusion. Here we report an analysis of archival observations of a sample of blazars (quasars whose jets point towards Earth) that overcomes previous limitations. We find a clear correlation between jet power as measured through the gamma-ray luminosity, and accretion luminosity as measured by the broad emission lines, with the jet power dominating over the disk luminosity, in agreement with numerical simulations. This implies that the magnetic field threading the black hole horizon reaches the maximum value sustainable by the accreting matter.

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.

The radio-loud quasar SDSS J114657.79+403708.6 at a redshift z=5.0 is one of the most distant radio-loud objects. The IR-optical luminosity and spectrum suggest that its black hole has a very large mass: M=(5+-1)x 1e9 Msun. The radio-loudness (ratio of the radio to optical flux) of the source is large (around 100), suggesting that the source is viewed at small angles from the jet axis, and could be a blazar. The X-ray observations fully confirm this hypothesis, due to the high level and hardness of the flux. This makes SDSS J114657.79+403708.6 the third most distant blazar known, after Q0906+693 (z=5.47) and B2 1023+25 (z=5.3). Among those, SDSS J114657.79+403708.6 has the largest black hole mass, setting interesting constraints on the mass function of heavy (larger than one billion solar masses) black holes at high redshifts.

We discuss how the interaction between the electrons in a relativistic jet and the Cosmic Microwave Background (CMB) affects the observable properties of radio-loud AGN at early epochs. At high z the magnetic energy density in the radio lobes of powerful radio-loud quasars can be exceeded by the energy density of the CMB (because of its (1+z)^4 dependance). In this case, relativistic electrons cool preferentially by scattering off CMB photons, rather than by synchrotron. Thus, sources sharing the same intrinsic properties have different extended radio and X-ray luminosities when located at different z: more distant sources are less luminous in radio and more luminous in X-rays than their closer counterparts. Instead, in compact regions where the local magnetic field still exceeds the CMB in terms of energy density, synchrotron radiation would be unaffected by the presence of the CMB. Such regions include the compact inner jet and the so-called hot spots in the radio lobes. The decrease in radio luminosity is larger in misaligned sources, whose radio flux is dominated by the extended isotropic component. These sources can fail detection in current flux limited radio surveys, and therefore they are possibly under-represented in the associated samples. As the cooling time is longer for lower energy electrons, the radio luminosity deficit due to the CMB photons is less important at low radio frequencies. Therefore objects not detected so far in current surveys at a few GHz could be picked up by low frequency deep surveys, such as LOFAR and SKA. Until then, we can estimate the number of high redshift radio-loud AGNs through the census of their aligned proxies, i.e., blazars. Indeed, their observed radio emission arises in the inner and strongly magnetized compact core of the relativistic jet, and not affected by inverse Compton scattering off CMB photons.

Mkn 421 is a nearby active galactic nucleus dominated at all wavelengths by a very broad non-thermal continuum thought to arise from a relativistic jet seen at a small angle to the line of sight. Its spectral energy distribution peaks in the X-ray and TeV gamma-ray bands, where the energy output is dominated by cooling of high-energy electrons in the jet. In order to study the electron distribution and its evolution, we carried out a dedicated multi-wavelength campaign, including extensive observations by the recently launched highly sensitive hard X-ray telescope NuSTAR, between December 2012 and May 2013. Here we present some initial results based on NuSTAR data from January through March 2013, as well as calibration observations conducted in 2012. Although the observations cover some of the faintest hard X-ray flux states ever observed for Mkn 421, the sensitivity is high enough to resolve intra-day spectral variability. We find that in this low state the dominant flux variations are smooth on timescales of hours, with typical intra-hour variations of less than 5%. We do not find evidence for either a cutoff in the hard X-ray spectrum, or a rise towards a high-energy component, but rather that at low flux the spectrum assumes a power law shape with a photon index of approximately 3. The spectrum is found to harden with increasing brightness.

B2 1023+25 is an extremely radio-loud quasar at z=5.3 which was first identified as a likely high-redshift blazar candidate in the SDSS+FIRST quasar catalog. Here we use the Nuclear Spectroscopic Telescope Array (NuSTAR) to investigate its non-thermal jet emission, whose high-energy component we detected in the hard X-ray energy band. The X-ray flux is ~5.5x10^(-14) erg cm^(-2)s^(-1) (5-10keV) and the photon spectral index is Gamma_X=1.3-1.6. Modeling the full spectral energy distribution, we find that the jet is oriented close to the line of sight, with a viewing angle of ~3deg, and has significant Doppler boosting, with a large bulk Lorentz factor ~13, which confirms the identification of B2 1023+25 as a blazar. B2 1023+25 is the first object at redshift larger than 5 detected by NuSTAR, demonstrating the ability of NuSTAR to investigate the early X-ray Universe and to study extremely active supermassive black holes located at very high redshift.

The search for extremely massive high redshift blazars is essential to put strong constraints on the supermassive black holes formation theories. Up to now, the few blazars known to have a redshift larger than 4 have been discovered serendipitously. We try a more systematic approach. Assuming radio-loudness as a proxy for the jet orientation, we select a sample of extremely radio-loud quasars. We measure their black hole masses with a method based on fitting the thermal emission from the accretion disc. We achieve a precision of a factor of two for our measures, thanks to the observations performed with the Gamma-Ray Burst Optical Near-Infrared Detector (GROND). The infrared to optical GROND data allow us to observe directly the peak of the disc emission, thus constraining the overall disc luminosity. We obtain a small range of masses, that peaks at 10^{9.3}Msun. If some of our candidates will be confirmed as blazars, these results would introduce interesting constraints on the mass function of extremely massive black holes at very high redshift. Moreover, all our blazar candidates have high accretion rates. This result, along with the high masses, opens an interesting view on the need of a fast growth of the heaviest black holes at very high redshift.

In powerful radio-quiet active galactic nuclei (AGN), black holes heavier than one billion solar masses form at a redshift ~1.5-2. Supermassive black holes in jetted radio-loud AGN seems to form earlier, at a redshift close to 4. The ratio of active radio-loud to radio-quiet AGN hosting heavy black holes is therefore a rather a strong function of redshift. We report on some recent evidence supporting this conclusion, gathered from the Burst Alert Telescope (BAT, onboard Swift) and by the Large Area Telescope (LAT, onboard Fermi). We suggest that the more frequent occurrence of relativistic jets in the most massive black holes at high redshifts, compared to later times, could be due to the average black hole spin being greater in the distant past, or else to the jet helping a fast accretion rate (or some combination of the two scenarios). We emphasize that the large total accretion efficiency of rapidly spinning black holes inhibits a fast growth, unless a large fraction of the available gravitational energy of the accreted mass is not converted into radiation, but used to form and maintain a powerful jet.

We observed 5 gamma-ray loud blazars at redshift greater than 2 with the X-Ray Telescope (XRT) and the UltraViolet and Optical Telescope (UVOT) onboard the Swift satellite, and the Gamma-Ray burst Optical Near-Infrared Detector (GROND) instrument. These observations were quasi simultaneous, usually within a few hours. For 4 of these blazars the near-IR to UV data show the presence of an accretion disc, and we could reliably estimate its accretion rate and black hole mass. One of them, PKS 1348+007, was found in an extraordinarily high IR-optical state, almost two orders of magnitude brighter than at the epoch of the Sloan Digital Sky Survey observations. For all the 5 quasars the physical parameters of the jet emitting zone, derived by applying a one-zone emission model, are similar to that found for the bulk of other gamma-ray loud quasars. With our observations we have X-ray data for the full sample of blazars at z>2 present in the Fermi 2-yrs (2LAC) catalog. This allows to have a rather complete view of the spectral energy distribution of all high-redshift Fermi blazars, and to draw some conclusions about their properties, and especially about the relation between the accretion rate and the jet power.

The radio-loud quasar SDSS J102623.61+254259.5, at a redshift z=5.3, is one of the most distant radio-loud objects. Since its radio flux exceeds 100 mJy at a few GHz, it is also one of the most powerful radio-loud sources. We propose that this source is a blazar, i.e. we are seeing its jet at a small viewing angle. This claim is based on the spectral energy distribution of this source, and especially on its strong and hard X-ray spectrum, as seen by Swift, very typical of powerful blazars. Observations by the Gamma-Ray Burst Optical/Near-Infrared Detector (GROND) and by theWide-field Infrared Survey Explorer (WISE) allow to establish the thermal nature of the emission in the near IR-optical band. Assuming that this is produced by a standard accretion disk, we derive that it emits a luminosity of L_d \simeq 9 \times 10^46 erg s^{-1} and that the black hole has a mass between 2 and 5 billion solar masses. This poses interesting constraints on the mass function of heavy (> 10^9 M_sun) black holes at high redshifts.

Many blazars detected by the Fermi satellite, observed spectroscopically in the optical, are line-less, and have been classified as BL Lac objects. Optical-UV photometry of nearly one hundred of them allowed to determine the redshift for a handful of objects and redshift upper limits for the great majority. A few of these are candidates to be "blue quasars", namely flat spectrum radio quasars whose broad emission lines are hidden by an overwhelming synchrotron emission peaking in the UV. This implies that the emitting electrons have high energies. In turn, this requires relatively weak radiative cooling, a condition that can be met if the main radiative dissipation of the jet power occurs outside the broad line region. We confirm this hypothesis by studying and modelling the spectral energy distributions of the 4 "blue quasars" recently discovered. Furthermore, we discuss the distribution of Fermi blazars in the gamma-ray spectral index -- gamma-ray luminosity plane, and argue that "blue quasars" objects are a minority within the blazar populations.

We selected all radio-quiet AGN in the latest release of Sloan digital sky survey quasar catalog, with redshift in the range 0.56-0.73. About 4000 (~80%) of these have been detected in all four IR-bands of WISE (Wide-field Infrared Survey Explorer). This is the largest sample suitable to study the disc-torus connection. We find that the torus reprocesses on average ~1/3-1/2 of the accretion disc luminosity.

We study the relation between the mass accretion rate, the jet power, and the black hole mass of blazars. To this aim, we make use of the Sloan Digital Sky Survey (SDSS) and the 11 months catalog of blazars detected at energies larger than 100 MeV by the Large Area Telescope (LAT) onboard the Fermi satellite. This allows to construct a relatively large sample of blazars with information about the luminosity (or upper limits) of their emission lines used as a proxy for the strength of the disc luminosity and on the luminosity of the high energy emission, used as a proxy for the jet power. We find a good correlation between the luminosity of the broad lines and the gamma-ray luminosity as detected by Fermi, both using absolute values of the luminosities and normalising them to the Eddington value. The data we have analyzed confirm that the division of blazars into BL Lacs and Flat Spectrum Radio Quasars (FSRQs) is controlled by the line luminosity in Eddington units. For small values of this ratio the object is a BL Lac, while it is a FSRQs for large values. The transition appears smooth, but a much larger number of objects is needed to confirm this point.

The gamma-ray emission of blazar jets shows a pronounced variability and this feature provides limits to the size and to the speed of the emitting region. We study the gamma-ray variability of bright blazars using data from the first 18 months of activity of the Large Area Telescope on the Fermi Gamma-Ray Space Telescope. From the daily light-curves of the blazars characterized by a remarkable activity, we firstly determine the minimum variability time-scale, giving an upper limit for the size of the emitting region of the sources, assumed to be spheroidal blobs in relativistic motion. These regions must be smaller than ~10^-3 parsec. Another interesting time-scale is the duration of the outbursts. We conclude that they cannot correspond to radiation produced by a single blob moving relativistically along the jet, but they are either the signature of emission from a standing shock extracting energy from a modulated jet, or the superposition of a number of flares occurring on a shorter time-scale. We also derive lower limits on the bulk Lorentz factor needed to make the emitting region transparent for gamma-rays interacting through photon-photon collisions.