A programme of worldwide, multi-wavelength electromagnetic follow-up of sources detected by gravitational wave detectors is in place. Following the discovery of GW150914 and GW151226, wide field imaging of their sky localisations identified a number of candidate optical counterparts which were then spectrally classified. The majority of candidates were found to be supernovae at redshift ranges similar to the GW events and were thereby ruled out as a genuine counterpart. Other candidates ruled out include AGN and Solar System objects. Given the GW sources were black hole binary mergers, the lack of an identified electromagnetic counterpart is not surprising. However the observations show that it is possible to organise and execute a campaign that can eliminate the majority of potential counterparts. Finally we note the existence of a “classification gap” with a significant fraction of candidates going unclassified.

Basic questions about black holes, some of which are fairly old, include (1) What is a black hole? (2) Do black holes exist? And the answer to this depends a good deal on the answer to (1), (3) Where, when, why, and how have they formed? and (4) What are they good for? Here I attempt some elaboration of the questions and partial answers, noting that general relativity is required to described some of the phenomena, while dear old Isaac Newton is OK for others.

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.

The jets from active galactic nuclei exhibit stability which seems to be far superior compared to that of terrestrial and laboratory jets. They manage to propagate over distances up to a billion of initial jet radii. Yet this may not be an indication of some exotic physics but mainly a reflection of the specific environment these jets propagate through. The key property of this environment is a rapid decline of density and pressure along the jet, which promotes its rapid expansion. Such an expansion can suppress global instabilities, which require communication across the jet, and hence ensure its survival over huge distances. At kpc scales, some AGN jets do show signs of strong instabilities and even turn into plumes. This could be a result of the flattening of the external pressure distribution in their host galaxies or inside the radio lobes. In this regard, we discuss the possible connection between the stability issue and the Fanaroff-Riley classification of extragalactic radio sources. The observations of AGN jets on sub-kpc scale do not seem to support their supposed lack of causal connectivity. When interpreted using simple kinematic models, they reveal a rather perplexing picture with more questions than answers on the jets dynamics.

The first gravitational wave transient GW150914 was observed by Advanced LIGO on September 14th, 2015 at 09:50:45 Universal Time. In addition to follow-up electromagnetic observations, the detection of neutrinos will probe deeply and more on the nature of astrophysical sources, especially in the ultra-high energy regime. Neutrinos in the EeV energy range were searched in data collected at the surface detector of the Pierre Auger Observatory within ± 500 s and 1 day after the GW150914 event. No neutrino candidates were found. Based on this non-observation, we derive the first and only neutrino fluence upper limit at EeV energies for this event at 90% CL, and report constraints on existence of accretion disk around mergers.

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 (z l = 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 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.

Here we report ongoing efforts for an untriggered search of rapid optical transients of various astrophysical and non-astrophysical origins on time scales down to fractions of a second with Mini-MegaTORTORA. Mini-MegaTORTORA is a novel 9-channel wide-field optical monitoring system in operation since 2014 at Special Astrophysical Observatory on Russian Caucasus.

I describe recent advances in the formulation and modeling of relativistic stellar dynamics around a massive black hole, and the implications for the rates of infall (e.g. tidal disruption) and inspiral (e.g. gradual decay by gravitational wave emission) processes, and their branching ratios.

One might expect far away from physical black holes that quantum field quantisation performed in Minkowski space is a good approximation. Indeed, all experimental tests in particle colliders reveal no deviations so far. Nevertheless, the black holes should leave certain imprints of their presence in quantum processes. In this paper, we shall discuss several local imprints of small, primordial evaporating black holes in quantum electrodynamics in the weak gravity regime. Physically this can be interpreted as being macroscopic manifestations of vacuum fluctuations.

Tidal disruption events are characterized by the prompt release of very high energy electromagnetic radiation, which rapidly carries away a significant fraction of mass-energy. The time scale od prompt release is quite comparable to the time scale of gravitational collapse. We suggest that generation of very high energy radiation from pulsar nebulae may be an example of the relativistic coupling between gravity and electromagnetism. The main ingredient of our picture comes from the observation that the electron - ion energy exchange time scale is much longer than the electron ion energy exchange time scale, which leads to thermodynamic decoupling of electron and ion gases, and thus to their different temperatures. Under such conditions the dominating hotter component pushes the colder component further out and thus generates a global electric field which constrains electrons and ions by their mutually generated electric field, yet allows them to reach very high kinetic energies

The HAWC gamma-ray observatory is a wide field of view and high duty cycle γ-ray detector investigating the 0.1 - 100 TeV energy range. It has detected supermassive black holes in the near Universe, and is seeking to detect black hole related objects like gamma-ray bursts, Galactic binary systems, primordial black holes and gravitational wave mergers. Daily light curves of the BL Lac objects Mrk 421 and Mrk 501 are presented here, together with a compilation of studies of black hole related objects.

The BL Lac object H1722+119 was observed in the very high energy band (VHE, E > 100 GeV) by the MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov) telescopes (Aleksić et al. 2016a, b)) between 2013 May 17 and 22, following a state of high activity in the optical band measured by the KVA (Kungliga Vetenskapsakademien) telescope. Optical high states are often used to trigger MAGIC observations, which result in the VHE γ-ray signal detection (see e.g. Aleksić et al. 2015, Ahnen et al. 2016 and references therein).

Gamma-Ray Bursts (GRBs) are the most violent explosions in the Universe, releasing a huge amount of energy in few seconds. While our understanding of the prompt and the afterglow phases has increased with Swift and Fermi, we have very few information about their High Energy (HE, E ≲ 100) emission components. This requires a ground-based experiment able to perform fast follow-up with enough sensitivity above ~ 50 GeV. The MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov) telescopes have been designed to perform fast follow-up on GRBs thanks to fast slewing movement and low energy threshold (~ 50 GeV). Since the beginning of the operations, MAGIC followed-up 89 GRBs in good observational conditions. In this contribution the MAGIC GRBs follow-up campaign and the results which could be obtained by detecting HE and Very High Energy (VHE, E ≳ 100 GeV) γ-rays from GRBs will be reviewed.

Since the first investigations into accretion onto black holes, astrophysicists have proposed effective Newtonian-like potentials to mimic the strong-field behavior of matter near a Schwarzschild or Kerr black hole. On the other hand, the fields of neutron stars or black holes in many of the alternative gravity theories differ from the idealized Schwarzschild or Kerr field which would require a number of new potentials. To resolve this, we give a Newtonian-like Hamiltonian which almost perfectly mimics the behavior of test particles in any given stationary space-time. The properties of the Hamiltonian are excellent in static space-times such as the Schwarzschild black hole, but become worse for space-times with gravito-magnetic or dragging effects such as near the Kerr black hole.

Recently, the Laser Interferometer Gravitational-Wave Observatory (LIGO) has detected the gravitational wave (GW) event, GW150914, as a result of the merger of a ~30 M ⊙ black hole (BH) binary. So far, the merger of stellar-mass BHs has been thought to result from the evolution of binary stars. Here, we propose a novel path of the merger stemming from non-binary isolated stars. In our previous studies, we have found that multiple non-binary stellar-mass BHs whose separations are larger than 1000 AU can merge with each other under the gas-rich environments through the gas dynamical friction and three-body interaction. In this case, a considerable amount of gas can accrete onto BHs before the merger, that is, the initial mass of BHs can be lower than 30 M ⊙. To explore this possibility, we perform post-Newtonian N-body simulations on mergers of accreting stellar-mass BHs. Based on our simulations, we find that the BH merger in GW150914 from smaller seed BHs is likely to occur in galactic nuclear regions or dense interstellar cloud cores. Furthermore, we roughly estimate event rates to be ~0.4 yr−1 in galactic nuclear regions and ~8 yr−1 in dense interstellar cloud cores.

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.

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.

Our understanding of gamma-ray bursts (GRBs) has come a long way in the past fifty years since their first detection. We now know that GRBs arise in distant galaxies and that there are at least two distinct sub-classes, the long-duration class being produced by some rare massive star core collapse and the short-duration class likely by compact binary mergers involved neutron stars. In both cases, the final remnant will be a stellar-mass black-hole or a massive neutron star. The bursts themselves are associated with ultra-relativistic jetted outflows created by these events, and their afterglows by the impact of these outflows on the surrounding circumburst material. Increasingly GRBs are also being used as probes of the universe, both for understanding galaxy evolution back to the era of reionization, and for the physics of gravitational wave sources. However, many aspects of GRBs remain poorly understood, some pointers to which are given here.

We analyze the fluctuations in the X-ray flux of 20 AGN (mainly Seyfert 1 galaxies) monitored by RXTE and XMM-Newton with a sampling frequency ranging from hours to years, using structure function (SF) analysis. We derive SFs over four orders of magnitude in the time domain (0.03-300 days). Most objects show a characteristic time scale, where the SF flattens or changes slope. For 10 objects with published power-spectral density (PSD) the break time scales in the SF and PSD are similar and show a good correlation. We also find a significant correlation between the SF timescale and the mass of the central black hole, determined for most objects by reverberation mapping.