We examine the 2008-2016 gamma-ray and optical light curves of a number of bright Fermi blazars. In a fraction of them, the periodograms show possible evidence of quasi-periodicities related in the two bands. This coincidence strengthens their physical meaning. Comparing with results from the periodicity search of quasars, the presence of quasi-periodicities in blazars suggests that the basic condition for its observability is related to the relativistic jet in the observer direction, but the overall picture remains uncertain.

Narrow-line Seyfert 1 (NLS1) galaxies provide us with unique insights into the drivers of AGN activity under extreme conditions. Given their low black hole (BH) masses and near-Eddington accretion rates, they represent a class of galaxies with rapidly growing supermassive BHs in the local universe. Here, we present the results from our multi-frequency radio monitoring of a sample of γ-ray loud NLS1 galaxies (γNLS1s), including systems discovered only recently, and featuring both the nearest and the most distant γNLS1s known to date. We also present high-resolution radio imaging of 1H 0323+342, which is remarkable for its spiral or ring-like host. Finally, we present new radio data of the candidate γ-emitting NLS1 galaxy RX J2314.9+2243, characterized by a very steep radio spectrum, unlike other γNLS1s.

With relatively small black hole masses and high accretion rates, narrow-line Seyfert 1 galaxies are thought to be young AGNs. About 7% of them are radio-loud narrow-line Seyfert 1 galaxies (RLNLS1s). RLNLS1s allow us to re-address some of the key questions regarding the physics of jet formation. As the first step of the systematic study on the jet properties of RLNLS1s, we present the radio structure of fourteen RLNLS1s from VLBA observations at 5 GHz in 2013. Although all these sources are very radio-loud with R > 100, their jet properties are diverse, in terms of their pc-scale morphology and overall radio spectral shape. The core brightness temperatures of our sources are significantly lower than those of blazars, therefore, the beaming effect is generally not significant, compared to blazars. This implies that the bulk jet speed may likely be low in our sources.

Low luminosity AGN (LLAGN) represent the vast majority of the AGN population in the near universe. However, they show intrinsic differences when compared to the Unified Model scenario, i.e. a lack of both the big blue bump and the IR bump in their spectral energy distribution (SED), which are the signatures of an accretion disk and a torus, respectively. In this work we present SED for a sample of six LLAGN in the nearby Universe obtained from subarcsec resolution observations. The nature of the IR continuum emission is discussed in the context of non-thermal radiation powered by compact jets.

We propose a thin accretion disk with magnetically driven winds model to investigate the high jet efficiency phenomenon in lobe-dominated luminous FRII radio galaxies, which was reported by Fernandes et al. (2011) and Punsly (2011). It was found that the temperature of a thin disk with winds was much lower comparing with that of a standard thin disk because the winds take the most of energy released in the disk away, resulting in a much smaller raditive efficiency. Therefore, the jet efficiency can be very high even for conventional jet power. Our model can explain the observations quite well.

We study the interaction of relativistic jets with their environment, using 3-dimen- sional relativistic particle-in-cell simulations for two cases of jet composition: (i) electron-proton (e− − p+) and (ii) electron-positron (e±) plasmas containing helical magnetic fields. We have performed simulations of “global” jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability and the Mushroom instability. We have found that these kinetic instabilities are suppressed and new types of instabilities can grow. For the e− − p+ jet, a recollimation-like instability occurs and jet electrons are strongly perturbed, whereas for the e± jet, a recollimation-like instability occurs at early times followed by kinetic instability and the general structure is similar to a simulation without a helical magnetic field. We plan to perform further simulations using much larger systems to confirm these new findings.

We present a study performed to understand the role of magnetic field during a historical outburst of IBL S5 0716+714 witnessed in early 2015. The two month long profile of the outburst reveals several episodes of sub-flares superimposed over the diminishing trend of the peak-flux. The broadband X-ray spectrum from Swift + Nustar exhibits a break at 4.9+0.4−0.5 keV which is consistent with the valley in the SED predicted by model. The spectral index variations, closely correlated to the flux variations during the first prominent sub-flare, indicate the alterations in the particle energy distribution coinciding the onset of the flare. The PA rotations approximately at the same epoch are used to constrain the field and its geometry.

The present work proposes a method to calculate the AGN magnetic luminosity fraction to be converted into ultra high energy cosmic rays (UHECRs) luminosities for nine UHECR AGN Seyfert sources based on the respective observation of gamma ray upper limits. The motivation for such calculation comes attached to the fact that a fraction of the magnetic luminosity (LB) produced by Kerr black holes in some AGNs can produce the necessary energy to accelerate UHECRs beyond the GZK limit, observed, e.g., by the Auger experiment. Nevertheless, the direct detection of those UHECRs has a lack of information about the direction of the source from where those cosmic rays are coming, since charged particles are deflected by the intergalactic magnetic field. Such problem arises alternative methods to evaluate the luminosity of UHECRs from the observation of upper limits during the propagation.

We report a significant hardening of the Fermi-LAT gamma-ray spectrum from the core of Cen A at E > 2.4 GeV, suggesting there is a source of high energy particles in the core of Cen A which is in addition to the jet component. We show that the observed gamma-ray spectrum is compatible with either a spike in the dark matter halo profile or a population of millisecond pulsars. This work gives a strong indication of new gamma-ray production mechanisms in active galactic nuclei and could even provide evidence for the clustering of heavy dark matter particles around black holes.

In 2014, a 179-orbit allocation of Hubble Space Telescope time anchored a massive reverberation-mapping program on the well-studied Seyfert 1 galaxy NGC 5548. Supporting imaging and spectrophotometric observations were provided by Swift, Chandra, Spitzer, and a world-wide network of ground-based telescopes. Understanding the data remains a significant challenge, partly because the level of detail is far beyond what has been seen before and partly because the behavior of the AGN was not typical of its past behavior. Based on analysis to date, the following conclusions can be reached: (1) the AGN accretion disk has a temperature profile that is consistent with that predicted by the Shakura–Sunyaev model, but is about three times larger than expected; (2) at least part of the broad-line region appears to be a Keplerian disk seen at intermediate inclination, and (3) the broad-line emission response from the far side of the disk is weaker than expected.

About 100 AGNs have their black hole mass measured directly using the reverberation mapping technique over the past few decades. By now we have high enough numbers to explore unique subsamples within these objects and to study phenomena across variety of AGNs. I will review recent reverberation mapping studies which focus on high-redshift high-luminosity AGNs and on AGNs with super-Eddington accreting massive black holes. These studies enable to investigate the BLR size, mass, and luminosity relations in different subsamples of AGNs and to check whether there are differences in these relations in different types of AGNs. In particular I will discuss the following questions: Is the BLR size - luminosity relation the same over the whole AGNs luminosity range? Are there different relations for different types of AGNs? What are these studies teaching us about theory of accretion into black holes in AGNs?

In hierarchical structure formation scenarios, merging galaxies are expected to be seen in different phases of their coalescence. Simulations suggest that simultaneous activity of the supermassive black holes (SMBHs) in the centres of the merging galaxies may be expected at kpc-scale separations. Currently, there are no direct observational methods which allow the selection of a large number of such dual active galactic nuclei (AGN) candidates. SDSS J142507.32+323137.4 was reported as a promising candidate source based on its optical spectrum. Here we report on our sensitive e-MERLIN observations performed at 1.6 and at 5 GHz, which show that the optical spectrum of the source can be more straightforwardly explained with jet–cloud interactions instead of the dual AGN scenario.

We propose an observational test for gravitationally recoiling supermassive black holes in active galactic nuclei, based on a positive correlation between the velocities of black holes relative to their host galaxies, |Δv|, and their obscuring dust column densities, Σdust, both measured along the line of sight. Our findings using a set of toy models implemented to a Monte Carlo simulation imply that models of the galactic centre and of recoil dynamics can be tested by future observations of the potential Σdust–|Δv| correlation. We have also found that the fraction of obscured quasars decreases with |Δv|, for which the predicted trend can be compared to the observed fraction of type II quasars, and can further test combinations of models we may implement.

Feedback released during the growth of supermassive black holes is expected to play a key role in shaping black hole-host galaxy co-evolution. Powerful, accretion disc driven winds have been invoked to explain both observed scaling relations (e.g., M − σ) and large-scale outflows with mass outflow rates of ~ 100 − 1000 M⊙ yr−1 and momentum rates of up to ~ 30 LAGN/c. Critically, how these winds couple to the host galaxy depends on if they are momentum or energy conserving. I outline observational signatures that could distinguish between these regimes and discuss their roles in establishing galaxy properties. Furthermore, I discuss high-resolution simulations exploring feedback in a multi-phase medium, highlighting how structural properties of galaxies can impact feedback efficiency. Finally, feedback, in the form of collimated jets, is expected to regulate cooling in galaxy clusters. I discuss new simulations of jet feedback using the moving-mesh code AREPO and outline the scope of our new study.

The blazar QSO B0218+357 is the first gravitationally lensed blazar detected in the very high energy (VHE, E > 100 GeV) gamma-ray spectral range (Ahnen et al. 2016). It is gravitationally lensed by the intervening galaxy B0218+357G (zl = 0.68466 ± 0.00004, Carilli et al. 1993), which splits the blazar emission into two components, spatially indistinguishable by gamma-ray instruments, but separated by a 10-12 days delay. In July 2014 a flare from QSO B0218+357 was observed by the Fermi-LAT (Large Area Telescope, Atwood et al. 2009, Ackermann et al. 2012), and followed-up by the MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes, a stereoscopic system of two 17m Imaging Atmospheric Cherenkov Telescopes located on La Palma, Canary Islands (Aleksić et al. 2016a, 2016b), during the expected time of arrival of the delayed component of the emission. MAGIC could not observe the leading image due to the Full Moon. The MAGIC and Fermi-LAT observations were accompanied by optical data from KVA telescope at La Palma, and X-ray observations by Swift-XRT (Fig. 1 left). Variability in gamma-rays was of the order of one day, while no variability correlated with gamma-rays was observed at lower energies. The flux ratio of the leading to trailing image in HE gamma-rays was larger than in the flare of QSO B0218+357 observed by Fermi-LAT in 2012 (Cheung et al. 2014). Changes in the observed flux ratio can be caused by gravitational microlensing on individual stars in the host galaxy (Neronov et al. 2015), or by other compact objects like for ex. clumps in giant molecular clouds (Sitarek & Bednarek 2016).

We studied the imprints that feedback from Active Galactic Nuclei (AGN) leaves on the intracluster plasma during the assembly history of galaxy clusters. To this purpose we used state-of-the-art cosmological hydrodynamical simulations based on an updated version of the Tree-PM SPH GADGET-3 code, comparing three sets of simulations with different prescriptions for the physics of baryons (including AGN and/or stellar feedback). We explore the effect of these different physics, in particular AGN feedback, on IntraCluster medium (ICM) properties observed via Sunyaev-Zel’dovich (SZ) effect using an extended set of galaxy clusters (~100 clusters with M500 masses above 5 × 1013M⊙/h). Some of the main findings are that the scaling relation between the integrated SZ flux and the galaxy cluster total mass is in good accordance with several observed samples, especially for massive clusters, and does not show any clear redshift evolution, with the slope of the relation close to the theoretical one in the AGN feedback case. As for the scatter of this relation, we obtain a mild dependence on the cluster dynamical state.

We present the most up-to-date and complete multi-wavelength correlation analysis on luminosity properties of TeV BL Lacs. Correlation function (power law or linear) parameters are calculated based on linear regression method. Using the lower energy luminosities of a sample of 182 non-TeV BL Lacs and the generated functions, minimum level of VHE gamma-ray emission was calculated for each non-TeV BL Lacs. This multi wavelength prediction method gives us a list of best candidates to be observed with current generation of Imaging Air Cherenkov Telescopes.

The blazar, OJ287, underwent a predicted period of flaring in late 2015, triggering multiple observing campaigns to conduct intensive monitoring. We present the polarisation data taken during this event with the Ringo3 polarimeter on the Liverpool Telescope and discuss two possible causes of the flares.

We analyzed the light curves (LCs) of several radio-quiet and radio-loud quasars belonging to the same parameter space volume in the 4D Eigenvector 1 (4DE1) quasar sequence, using data from the Catalina Real-time Transient Survey (CRTS). We report preliminary results on detected variability pattern, and discuss possible cases of periodic variability.

Highly accreting quasars are possible cosmological probes, as their Eddington ratio is expected to saturate toward values of order unity. We present preliminary estimates of redshift- independent source luminosities and the Hubble diagram for quasars in the redshift range 0.1 ≲ z ≲ 2.6.