On the High-Energy Neutrino Emission from Active Galactic Nuclei †

We review observational aspects of the active galactic nuclei and their jets in connection with the detection of high-energy neutrinos by the Antarctic IceCube Neutrino Observatory. We propose that a reoriented jet generated by the spin-flipping supermassive black hole in a binary merger is likely the source of such high-energy neutrinos. Hence they encode important information on the afterlife of coalescing supermassive black hole binaries. As the gravitational radiation emanating from them will be monitored by the future LISA space mission, high-energy neutrino detections could be considered a contributor to multi-messenger astronomy.


Introduction
Electromagnetic radiation has been providing a first observational window to the Universe for about 400 years, while a second window opened with the detection of the first flavour of neutrinos in 1956 (electron neutrino, [1]) and the ultra-high energy cosmic ray particles (UHECRs) in 1962 [2].The third observational window to the Universe opened up on 14th September 2015 [3], when the two Advanced LIGO observatories firstly detected a gravitational wave (GW) signal, emanated from the coalescence of two astrophysical black holes (BH).Up to date four GW detections together with a lower significance candidate from binary black hole sources, and a GW detection from colliding neutron stars were announced by the LIGO Scientific Collaboration and Virgo Collaboration [3][4][5][6][7][8].
High-frequency GWs (hfGW, of order of Hz to kHz) are emitted from the coalescence of astrophysical BHs or neutron stars, while low-frequency GWs (lfGW, of order of microHz to milliHz) are emitted from the merging of supermassive black holes (SMBHs).The LIGO and Virgo detectors are sensitive to hfGWs.The direct detection of lfGWs is the purpose of the future Laser Interferometer Space Antenna (LISA) space mission [9], an experiment facilitating astrophysics and cosmology research in the mHz regime.Providing priors and predictions to this experiment is of paramount importance.
The compact centre of some galaxies is radiating tremendous energy across the electromagnetic spectrum due to the accretion onto the central SMBH, believed to be the main engine of the activity in active galactic nuclei (AGN; for a review see e.g., [10]).Observations on binary black hole candidates in kpc-scale separated binaries are naturally feasible with spatially resolved AGN.However, one can reveal the presence of a sub-pc-scale separated binary in AGN only by indirect approaches.When an AGN is radio-loud, its periodic jet structure might indicate, that the jet emitter black hole interacts gravitationally with another black hole (e.g., [11][12][13][14][15][16][17]).In absence of strong relativistic jets, the compact centre of the galaxy is observed as a radio-quiet AGN.In case of these AGN their emission in optical wavelengths might reveal the binary nature of the centre of the galaxy (e.g., [18]).Evidence for a merger is detectable for a long time period, starting just before the actual merger, and persisting for quite a long period afterwards, estimated in thousands of years, based on the traversal time scale for the jet to penetrate initially the inner high density region of the host galaxy (e.g., [19]).

Active Galactic Nuclei and Their Jets
According to the accepted evolutionary models, galaxies gather their masses through a sequence of accretion and merger phases (e.g., [20][21][22]).As a result of the process, million-to-billion solar mass SMBHs formed at the centre of the massive galaxies [23], with some of them members of binary systems due to the long time-scale of the merger [24,25].An accretion disk may form around the SMBH due to the non-zero angular momentum of the infalling matter [26].As substantial energy is released in the accretion disk due to its efficient mass-energy conversion [27], the compact centre of the galaxy becomes an AGN (e.g., [28]).SMBHs are believed to be the main engine of the activity in AGNs [10].General relativistic modelling of magnetized jets from accreting BH showed, that provided the central spin is high enough [29], a pair of relativistic jets is launched from the immediate vicinity of the BH (e.g., [30]), composed mostly by electrons and positrons (in light jets), or electrons and protons (in heavy jets).These jets could be powered either by the rotational energy of the central SMBH (through the Blandford-Znajek process, [31]), or by the magnetized accretion disk wind accelerated by magneto-centrifugal forces (Blandford-Payne process, [32,33]).
Particles are accelerating on helical paths along the magnetic field lines, emitting synchrotron radiation (e.g., [34]).In the Fourier transform of the continuum spectrum of the synchrotron radiation the most powerful frequency is the critical frequency (e.g., [35]).The electron energies and magnetic field strengths typical to AGN render the critical frequency into the radio regime, so the collimated jets are observed as radio-loud AGN (approximately 10 per-cent of single AGN qualify [36,37]).
Apart from the synchrotron radiation in radio wavelengths, AGN and their jets emit energetic particles in the X-ray and γ-regime.According to the Pierre Auger and IceCube Collaborations the main sources of UHECRs and cosmic high-energy (HE) neutrinos are AGNs (e.g., [38][39][40][41][42]).Two interaction channels were proposed for neutrino production in AGN: one with protons interacting with ambient photons (from the accretion disk, synchrotron photons emitted from the jet, CMB photons strayed into the jet); and another with protons interacting with other protons within the jet or with protons of the external material trapped in the jet flow (e.g., [43][44][45]).

High-Energy Neutrino Emission in AGN
Various scenarios for HE neutrino production in AGN jets have been suggested (e.g., [46][47][48][49]); it is plausible to assume that the observed cosmic HE neutrino background originates from some combination of these.The common property of the models is that there has to be a process inducing shock-waves in the accretion disk or the jet, that accelerates protons to energies exceeding the threshold of pion-production, that decay further to HE neutrinos, among other particles.It means the incoming proton has to have at least a relativistic mass of 1.3 times its rest mass, and thus a kinetic energy around 290 MeV, or moving with at least 0.65c (e.g., [50]).This condition is naturally satisfied in AGN jets, since the Lorentz factor of the turbulent flow is usually pretty high in such systems.
In [51] three of the present authors proposed a scenario of binary black hole evolution leading to the emission of HE neutrinos.The scenario relies on the final coalescence of two SMBHs, followed by the emission of lfGWs, an enormous shock-wave and the formation of a new energetic jet.The freshly started jet has high Lorentz factor, and is able to accelerate the particles leading to the emission of HE neutrinos.By cross-correlating both the Parkes Catalogue and the Second Planck Catalogue of Compact Sources with the arrival direction of the track-type neutrinos detected by the IceCube Neutrino Observatory, we identified four track events, all with flat spectrum radio sources, and of those again two with flat spectrum radio sources with their spectrum extending to near THz.The brightest of these two, the flat-spectrum blazar PKS 0723-008 emerged as a good candidate for the high-energy neutrino event 5 [51].
The following quantitative picture illustrates what may be common.We observe UHECR particles, and two scenarios are viable to explain them, with probably both realized in nature [52]: one is based on the relativistic jets and their shock waves in radio galaxies (e.g., [46]), and the other on Gamma Ray Bursts (GRBs) and their relativistic jets.Both sites of origin may involve mergers of black holes or neutron stars [51,53], or the ab initio formation of black holes.In a series of papers by Mannheim, Falcke and Markoff (e.g., [43,54,55]) it has been shown that at a distance of order 10 3.5 to 10 5 Schwarzschild radii from a central black hole a reconfinement shock system occurs, starting the routine jet flow (e.g., [56,57]).So when a magnetic jet system is newly established in a new direction [58] during a merger event, the jet flow has to be re-established (the characteristic velocity for this may be the relativistic speed of sound c/ √ 3), particles have be newly injected, and built up as a viable population.The observation of a time lag between the GW event and the HE photons [59] of 1.74 s is compatible with this characteristic distance for a central black hole of a mass of order 1 to 2 Solar masses: the characteristic distance introduced above yields a range of expected time lags of 10 −0.45 to 10 +1.05 s for a neutron star merger resulting in a black hole, fully compatible with the observations; the additional time lags due to the build-up of a viable population of energetic particles may be as fast or even faster.This time scales with the mass of the central black hole, and so the merger of two super-massive black holes of order 10 8 Solar masses the corresponding time scale is then about 3 years.There is an uncertainty in this, as the final jet direction may already be established before the actual merger occurs [58], but the actual merger will surely re-energize the jet, since the jet gets a new foot structure, switching from the dominant compact object of a pair to a single merged compact object.Another uncertainty is the relativistic geometry effects, like boosting, ignored here; relativistic effects could shorten the lag time-scale.The interaction of the newly energized population of hadronic particles then produces both HE photons and neutrinos.
Apart from AGN, for example the interaction of a freely expanding SNIIn-ejecta with the circumstellar medium gives rise to a shock-wave propagating in the dense SN environment, which may accelerate protons to PeV energies.The cumulative neutrino emission from such process is found to be ∼10 per-cent of the observed IceCube neutrino flux above 60 TeV [60].
The merger of a neutron star binary may result in the formation of a millisecond magnetar remnant surrounded by a low-mass ejecta shell.A portion of the magnetar's rotational energy is deposited behind the ejecta in a pulsar wind nebula.Ions in this pulsar wind may also be accelerated to ultra-high energies, providing HE cosmic rays and neutrinos.The predicted neutrino spectrum from such sources peaks just below the 90% confidence level sensitivity curve of IceCube Observatory [61].In optimistic scenarios the cumulative neutrino background resulting from these events may be observed by the IceCube in the near future, in less optimistic scenarios the diffuse flux is potentially within the reach of next-generation neutrino telescopes [61].

Observation of Possible Sources of Low-Frequency GWs
At present, radio interferometry, notably the technique of very long baseline interferometry (VLBI) provides the finest angular resolution images of AGN jets.The periodicities observed in the parsec-scale jets are consistent with SMBH binaries in the inspiral evolution stage (e.g., [16,17,62]).These sources often exhibit significant periodic variability in their optical brightness or periodic modulation of the shape of their emission lines due to the orbital motion (e.g., [63][64][65][66]).In this stage gravitational radiation plays the leading role in the dissipation of the total energy and the reduction of the separation of the binary, surpassing the effect of dynamical friction [67].The gravitational lifetime of such systems is million to billion years and the GW signal of their merger is considerable only after a long time-scale, hence yet undetectable by LISA.This begs the question of how can even tighter binaries with decade-long gravitational lifetime be identified as sources for LISA.Observation of such sources in the electromagnetic spectrum or via HE neutrinos or UHECRs is important for breaking degeneracies in their parameter estimation.
LISA will be sensitive in a frequency range from 0.1 mHz to 1 Hz, being able to detect GWs emmited by massive binaries with total mass m(= m 1 + m 2 ) ∈ [10 4 ÷ 10 7 ]M Sun , mass ratio ν(= m 2 /m 1 ) ∈ [1/10 ÷ 1] [9].The equation of motion of the binary is (e.g., [68]): where r is the binary separation (r = |r|), ε = Gm/c 2 /r is the post-Newtonian (PN) parameter, employed as expansion paramter, with G the gravitational constant, c the speed of the light and O(ε n ) represents the nth PN order.Since the gravitational radiation tends to circularize the orbits [69], we discuss circular orbits.In that case the orbital period is given by With the gravitational radiation the dominant dissipative effect the binary radiates its orbital angular momentum L at rate L, from which the inverse of its gravitational lifetime emerges as [58]: where η = m 1 m 2 /m 2 = ν/(1 + ν) 2 is the symmetric mass ratio.The frequency of the GW to leading order is given by: which is twice the orbital frequency f T = 1/T (T is the orbital period).We plot in Figure 1  The PN-parameter ε was expressed with the total mass and the orbital period using Equation ( 2), and the resulted function ε(m, T) was substituted to Equation (3).Sub-parsec separated binary systems coalescing presumably on time-scale of up to several decades are the potential sources of the lfGWs: these are either binaries at the high-mass end of the LISA range with up to several days orbital period, or its low-mass end with several hour orbital period (see Figure 1).Such close SMBH systems are impossible to resolve due to their closeness, even with the astronomical techniques giving the finest angular resolutions.Indirect methods, based on the studies of X-shaped radio galaxies, accretion disks with central cavities, ripples in FeKα X-ray lines, HE neutrino emission are rather explored.
A theoretically promising technique relies on the investigation of the FeKα line.The inner edge of the accretion disk is hot enough to excite this emission line (E = 6.4 keV, λ = 0.19 nm).The orbital motion of the secondary BH opens a gap in the disk, and it affects the line profile.The authors of [70] have simulated the FeKα line of the accretion disk of a SMBH in a 10, 20, 50, 90 Schwarzschild radii separated binary (with 10 per-cent error on the separation), and their simulations resulted ripples on the theoretical line-shape induced by such orbital motion, which do not appear on the unperturbed iron line.Considering the closeness of these systems they could be possible sources of strong GW bursts over human lifetime.However, the X-ray observations are far more difficult in technique, than e.g., optical observations, and the S/N of the present X-ray spectra does not allow to fit such rippled FeKα models.Additionally, the inner edge of the accretion disk may depend not only on the SMBH mass and spin, but also on the size of the jet-launching region or on the perturbations induced by an accompanying SMBH.Hence the plasma temperature at the modified inner edge cannot reach the level necessary to excite the FeKα emission line.
An important question arises then, on which other ways could we identify such merging binaries.For this purpose particle physical processes following the final coalescence of SMBH binary systems can prove useful.

Detection of Neutrinos by the IceCube Neutrino Observatory
IceCube, located at the Amundsen-Scott South Pole Station, uses one cubic kilometre of the Antarctic ice as detector matter.The detector strings penetrate up to 2500 m deep underground, and are sensitive to the Cherenkov radiation of electrons and muons emerging from the interaction of electron and muon neutrinos with the ice.Muons leave track-type paths, such that their origin can be located in the sky within about 1 degree circle.Electrons leave shower-type paths, the directional uncertainty of their origin being several ten degrees.Background events (99.9% of the events) producing first light in the veto region of IceCube are discarded as entering tracks (usually from cosmic ray muons entering the detector).The deposited energy is also a veto, only the HE (TeV, PeV) neutrinos are considered to be extragalactic.Up to date the IceCube Collaboration has published the parameters of 55 cosmic HE neutrino detections [71][72][73], from which 15 events were track-type.An excellent review of the high-energy neutrino astrophysics with IceCube can be found in [74].
Identification of the origin of these HE neutrinos [71][72][73] is an actively pursued topic in astroparticle physics consisting of constructing models for the neutrino emission, then comparing their predictions with actual detections.In statistical sense there are two fundamental ways to adjust a model to data.First, determine the best-fitting parameters.Second, answer how likely is the model with a set of parameters, given the data.This method relies on the Fisher concept of likelihood of a set of parameters having an assigned set of values being proportional to the probability, that with them the totality of observations emerge.
The energy spectrum of the UHECRs is measured to be power-law (e.g., [75]).It is assumed that the neutrino spectrum essentially traces the proton spectrum (e.g., [49]), and the source of neutrinos is assumed to emit them according to a power-law energy spectrum.The question is then what is the event rate for a neutrino model predicting a flux given the energy, solid angle and observing frequency.Then a hypothesis test follows, where the null-hypothesis (H 0 ) is that the data consists solely of background atmospheric neutrino events, and the alternative hypothesis (H S ) states that the data consists of atmospheric neutrino events as well as astrophysical neutrino events produced by a source with some strength and energy spectrum [76].These latter are outputs of a quantitative model leading to the neutrino emission.The likelihood of obtaining the data given each hypothesis is calculable; the usually employed test statistic is the log of the likelihood ratio λ (e.g., [76]): λ = −2 log P(Data|H 0 ) P(Data|H S ) .( 5) Large (small) values of λ indicate that data is less (rather) compatible with the background hypothesis H 0 .The probability density functions (PDF) P(Data|H 0 ) and P(Data|H S ) are calculated using knowledge of the spatial and energy distribution of background and astrophysical neutrino events [76].The source PDF is a Gaussian, and it depends on the space angle difference between the direction of the reconstructed event ( x i ) and the source ( x s ), the reconstruction error estimate (σ), and the probability distribution function of the neutrino spectrum given its spectral index γ (P(E i |γ)): where there are tabulated values for P(E i |γ) (e.g., [76]).The background PDF depends on the probability of the atmospheric neutrino flux (φ atm ) given the energy (E i ), over the solid angle (Ω): where P(E i |φ atm ) is precalculated assuming φ atm [77].Then the source and background PDFs are combined, and the likelihood L( x s , n s , γ) is evaluated over all events in the band for a source supposed to emit neutrinos: Then the test statistics λ written in terms of L becomes [76] λ = −2sign( ns ) log L( x s , 0) L( x s , ns , γ) The fraction of signal events (n s ) as well as the source spectral index (γ) are unknown and must be determined by maximizing the likelihood L. This is done by minimizing the quantity − log L with respect to the unknown quantities n s and γ, obtaining the best value of each of these parameters ns and γ.

Discussion: Reorienting Jets, as Sources of High-Energy Neutrinos
If the spins are misaligned with the orbital angular momentum, a spin-orbit type precession will cause the spin vectors to slowly rotate about the orbital angular momentum vector in the inspiral phase [78,79].For the typical mass ratio of the merging SMBHs ν = m 2 /m 1 ∈ [1/30 ÷ 1/3] (ν ≤ 1) the spin of the primary SMBH is dominant [58], and the spin of the secondary SMBH can be neglected.
In binary SMBHs with the typical mass ratio the dominant spin-flips, and the jet is reoriented along the spin-flipped dominant spin.This process is the most probable driving mechanism for the formation of the X-shaped radio galaxies [80].In them the remnant of the old jet appears together with the freshly made jet, forming an X-shape.The old jet-pair has steep continuum spectrum emanated from slowly moving charged particles.The jet-pair formed after the spin-flip is bright and has flat continuum spectrum due to the energetic and fast, charged particles.This morphology can reveal either a recent merger of two SMBHs or the presence of a second active SMBH in the nucleus.
In [51] the blazar PKS 0723-008 was identified as an excellent candidate source of a track-type HE neutrino event, as its spectrum is flat up to high Planck frequencies, and it produced a fivefold-increased radio flux density through the last decade, supporting the scenario of the track-type HE neutrino event related to a binary SMBH merger, more specifically to the plunge of a newly formed jet into the environment of the binary SMBH.The observable consequences are: (i) the high Lorentz factor of the freshly made jet due to the gravitational shock-wave, (ii) enhanced radiation in all EM frequencies, (iii) the spectrum of the AGN is flat up to THz frequencies due to the energetic synchrotron-radiating electrons, (iv) its radio flux density is increasing as the electrons are speeding, leading to enhanced synchrotron radiation, (v) neutrino emission.
These were all verified to hold true.The sources of neutrinos might be energetic proton-proton collisions, in which the proton energy is larger than the limit to create pions in the process (larger then 290 MeV, assuming two protons moving with equal and opposite velocities).The pion decays into a muon and a muon neutrino, while the muon decays into an electron, an anti-muon neutrino and an electron neutrino.If the kinetical energy of the colliding protons cannot supply for the rest-mass of the pion, they fuse into a diproton, which β + decays into a deuterium, while emitting a positron, an electron-neutrino, and a γ-ray photon.Therefore the detection of the muon neutrinos has to be explained with HE proton collisions.
Finally we pinpoint here the usefulness of constructing a quantitative model leading to the emission of HE neutrinos, which would give the probability of their detections from a given source.By fitting related parameters of the model (neutrino event fraction rate, spectral index) to already existing neutrino and AGN data, such as the arrival direction and energy distribution of HE neutrinos, position and energy spectrum of AGNs, one could predict how likely an AGN, given its spectral properties, will emit a neutrino at a particular time due to a post-merger SMBH binary at its centre.