Blazar Jets Perturbed by Magneto-Gravitational Stresses in Supermassive Binaries

A. Cavaliere11affiliation: Astronomia, Accademia Nazionale dei Lincei, via della Lungara 10, I-00165 Roma, Italy 22affiliation: INAF/IAPS–Roma, via del Fosso del Cavaliere 100, I-00133 Roma, Italy , M. Tavani11affiliation: Astronomia, Accademia Nazionale dei Lincei, via della Lungara 10, I-00165 Roma, Italy 22affiliation: INAF/IAPS–Roma, via del Fosso del Cavaliere 100, I-00133 Roma, Italy 33affiliation: Universitá “Tor Vergata”, Dipartimento di Fisica, via della Ricerca Scientifica 1, I-00133 Roma, Italy 44affiliation: Gran Sasso Science Institute, viale Francesco Crispi 7, I-67100 L’Aquila, Italy , V. Vittorini22affiliation: INAF/IAPS–Roma, via del Fosso del Cavaliere 100, I-00133 Roma, Italy

We study particle acceleration and radiative processes in Blazar jets under recurring conditions set by gravitational perturbations in supermassive binary systems. We consider the action from a companion orbiting a primary black hole of , and perturbing its relativistic jet. We discuss how such conditions induce repetitive magneto-hydrodynamic stresses along the jet, and affect its inner electron acceleration and radiative processes. Specifically, we study how macroscopic perturbations related to increased jet ”magnetization” end up into higher radiative outputs in the optical, X-ray and gamma-ray bands. We find first an increase in magnetic field strength as gauged in the optical band from the Synchrotron emission of electrons accelerated in kinetic processes stimulated by reconnecting magnetic lines. The energetic electrons then proceed to up-scatter the Synchrotron photons to GeV energies after the canonical Synchrotron-Self Compton radiation process. Our model implies a specific, recurring pattern in the optical to gamma-ray emissions, made of high peaks and wide troughs. Progressing accelerations caused by spreading reconnections will produce an additional Synchrotron keV component. Such outbursts provide a diagnostics for enhanced acceleration of electrons which can up-scatter photons into the TeV range. We discuss how our model applies to the BL Lac object PG 1553+113, arguably the best candidate to now for high amplitude, recurring modulations in its gamma-ray emissions. We also consider other BL Lacs showing correlated keV - TeV radiations such as Mrk 421.

gamma rays: observations – BL Lac Objects; individual: PG 1553+113, Mrk 421.
slugcomment: Submitted to ApJ: July 27, 2016; accepted: Jan. 8, 2017

1 Introduction

Blazars are singled out among the Active Galactic Nuclei (AGN) by their relativistically collimated jets with bulk Lorentz factors (Urry & Padovani 1995). They are launched by a central super-massive black hole (SMBH) with mass , and are Doppler-boosted when a jet happens to be closely aligned with our line of sight.

Blazars are extreme in several respects. The observed outputs are very bright, up to erg/s in isotropic extrapolation, and strongly variable on diverse timescales from years to minutes depending on the observed bands. They show highly non-thermal spectra with energy distributions (SEDs) constituted by two humps: one peaking in the -UV bands of clear Synchrotron (S) origin; the other extending from hard X to gamma rays, of likely inverse Compton (IC) nature (for basics see Rybicki & Lightman 1979). Both are produced by highly relativistic electrons with random Lorentz factors up to that inhabit the jets.

Two main Blazar flavors are discerned (see Peterson 1997, Ghisellini 2016). The BL Lac-type sources feature two comparable spectral humps. The Flat Spectrum Radio Quasars (FSRQs), on the other hand, feature Compton-dominated spectra at gamma-ray energies, but also a conspicuous Big Blue Bump and the strong, broad optical emission lines common to many quasars; these thermal features yield evidence of a gas-rich environment surrounding the SMBH out to some pc.

Relativistic electrons within the jet emit by synchrotron (S) radiation an observed (and isotropically extrapolated) luminosity in terms of the density of relativistic electrons within the source size , and of the general magnetic field G threading the jet. On the other hand, the IC scattering operates on ”seed” photons on conserving their number while upgrading their energy by another factor , so as to yield photon energies and related luminosities .

The simplest source structure envisages as seeds the very S photons that are radiated by a single population of relativistic electrons in one, homogeneous zone of the jet, so as to yield also a beamed flux of gamma rays . Such a synchrotron self-Compton (SSC) radiation process (e.g., Maraschi, Ghisellini & Celotti, 1992) is sufficient to account also for the gamma-ray outputs in tight correlation with O - X as featured by many BL Lacs in their quiescent states.

An interesting development is constituted by growing hints or signs of apparently recurrent emissions in a number of Blazars, BL Lacs in particular. On the origin and modeling of these phenomena we will focus our attention in the present paper.

2 A Binary Driver of Recurring Jet Instabilities

Interest in binary, massive BHs in AGNs started with Begelman et al. 1980, and is currently being spurred by the recent detections of bursts of gravitational waves constituting the events GW 150914 (Abbott et al. 2016a) and GW 151226 (Abbott et al. 2016b). These events have been shown to originate from binary systems of BHs with intermediate masses . The pairs have been caught in their late inspiral stages, as they were accelerating their orbital motions under increasing energy loss in gravitational waves to end up into catastrophic coalescence as predicted. In fact, extensive numerical work concerning binary BHs long pre-dated those events, aimed at including electromagnetic outputs from systems composed of two BHs and a warm disk or ring of gas (see Palenzuela et al. 2010, and the review by Baumgarte & Shapiro 2011). Similar events are expected also in binary systems of two SMBHs (see Colpi 2014, Volonteri et al. 2015).

Ever as such emissions are dramatically enhanced during the final coalescence, they are numerically found to set in gently at much lower rates during the early, long inspiral stage. Similar computations can be focused as well on investigating the initially slow evolution of nearly Keplerian orbits and their interactions with the accretion disk; we aim at understanding non-linear light curve shapes recurring on scales of a few years in the e.m. emissions, as suggested by current gamma-ray observations discussed in Sects. 3 and 4.

Figure 1: A schematic of our model for the SMBH binary system underlying a BL Lac Blazar. The primary SMBH produces a relativistic jet towering well above its formation region in the inner accretion disk (boundary marked by the thin curve). A SMBH companion orbits around the system barycenter close to the primary, following the red orbit. The gravitational force affects the jet base during its passes close to the periastron, modulating its emissions to a pattern with high peaks followed by long relaxations. The orbital size is of order for a total mass around , mass ratio , with an intrinsic period  years.

The emissions of all Blazars are governed by the hydrodynamic and magnetic properties of the jet’s relativistic plasma outflow (see Sect. 1). Jet formation and launching in many such sources appear to persist on gross average over timescales of many years. The considerable short term variations observed in their radiative outputs have often been considered to be fully erratic (see, e.g., Kelly et al. 2009). Jet formation and its short term instabilities are then attributed to random variability in the boundary conditions occurring at the base, such as accretion rate; alternatively, disk instabilities (see Nixon & King 2015, and refs. therein) have been computed and discussed. In the present paper we consider an additional possibility.

We will investigate the binary SMBH scenario depicted in Fig. 1, where the primary launches a main jet that dominates the observed fluxes, but is perturbed by the secondary to a considerable degree, This we shall compute from considering the binary dynamics of a system with a total mass of , mass ratio and a current orbital size around several milli-pc. Specifically, the bulk properties of the main jet are affected by episodes of shaking, squeezing and twisting caused by differential gravitational and electromagnetic stresses along the non-circular, or even not simply Keplerian orbits. The results include time-modulations of the average strength of the magnetic field, coordinated with changes in the local topology of B-lines by reconnecting and rearranging those that angularly diverge across narrow layers, as we shall discuss in the next Section.

Here we just note that passive geometrical effects of gentle beam redirection along steadily curved or helical jets as advocated, e.g., by Raiteri et al. (2015) and discussed by A15 may be re-interpreted in terms of initially mild and slow dynamical effects that take place during the quasi-period. Our discussion will focus instead onto the dynamical driving in the jet of sites suitable for magnetic reconnections to take place.

3 Unstable Jets

We consider how the progression from large-scale jet instabilities proceeds through an intermediate stage of tearing perturbations to end up into kinetic effects that spread out electron acceleration to high energies.

We base our discussion on the canonical SSC radiation process recalled in Sect. 1. As to the primary S emission, we adopt standard values of the average magnetic field G in the source frame, at the height of some cm above the primary SMBH. These are modulated in time by factors of about three on timescales of a few years by compressing and bending B as discussed below.

The other basic component for S emission is constituted by highly relativistic electrons. Their acceleration occurs in regions within the jet that are affected by magnetic field lines packed, sheared and reconnected so as to induce macroscopic and kinetic effects. For a review covering the magnetic reconnection (MR) theory and the ongoing numerical simulations see Kagan et al. (2015), also Melzani et al. (2014). For recent theoretical developments on reconnecting structures in collisionless plasmas, see Coppi (2016); for detailed observational results concerning reconnections in a magnetospheric environment, see Burch et al. (2016).

The intermediate output of MRs and associated tearing instabilities is constituted by strings of ”magnetic islands” or plasmoids, many of which then merge into a few giant ones. Meanwhile, on the kinetic side sharply sheared, or even annihilating B-lines induce strong if confined E-lines; these repeatedly and efficiently accelerate electrons in situ to high values of , particularly around the giant, coalescing plasmoids and in the intervening gaps. Next, we focus on the key parameters governing these two, related sides of the processes occurring in magnetized jets.

3.1 The macroscopic side

The relevant parameter for stability on large scales is provided by the bulk magnetization of the jet


This depends on the average field that appears in the magnetic stress , and on the kinetic stress dominated by the jet protons (see Celotti & Ghisellini 2008).

The condition yields stability of a jet on large scales, for example, on the scales cm numerically investigated by, e.g., Mignone et al. (2013); on the other hand, values are found by the same authors to promote jet instabilities. Thus the magnetization constitutes the main parameter governing the overall stability of a jet configuration on large scales.

Using values of G, related to the accretion rates on the primary SMBH, and values , as derived, e.g., by Paggi et al. (2009), Tavani et al. (2015), and observed by Hovatta et al. (2009), we find to hold in FSRQs, while prevails in BL Lacs. Thus the latter appear to sit on the brink of instability, and may be driven unstable by minor dynamical perturbations such as may be induced by a binary BH companion. On evaluating at (details are given in the Appendix) the period fraction with strong primary acceleration while revolving around the barycenter, it is seen that the jet cannot remain connected over scales cm against dynamical perturbations that shake its base on an effective timescale of a few months. This also sets an upper limit to any spread along the axis of the sources of correlated emissions.

3.2 The kinetic side

At the opposite extreme constituted by the local kinetic level, acceleration of electrons to high energies is governed by repeated MR events. The acceleration process has been parametrized and numerically computed in terms of the electron local magnetization (see the review of Kagan et al. 2015)


which depends on the local magnetic field and electron density , given their mass .

Many numerical calculations have been carried out in the range of values from tens to up to several hundreds (see Sironi & Spitkovsky 2014, Melzani et al. 2014). Remarkably, for increasing the resulting electron energy distributions approach flat power-laws up to, or asymptotically exceeding limited by computing times (e.g., Melzani et al. 2014). We stress that in the simplest case the relation


holds, so that the electron energies can grow to values and possibly beyond as soon as bulk instability conditions are established by values111Note that Eq. 3 as it stands applies for values of . For larger values, the electron contribution to the total kinetic energy should not be neglected; is correspondingly decreased, and so is the effective . This feature may imply self-regulation to limiting energies . of . In other words, conditions conducive to bulk instability act as a trigger for strong particle accelerations to occur. Note, however, that the electron magnetization in Eq. 2 affects the particle energy distribution function only in the presence of topological configurations of the local magnetic field that include sheared layers or close reversals conducive to tearing and reconnections in the jet’s collisionless plasma.

Extensive and detailed computations have recently been presented by Yuan et al. 2016 concerning the link between macroscopic field configurations with higher magnetization and kinetic conditions conducive to acceleration. We conclude that a slow time shift of macroscopic conditions toward increased and sheared B topology can end up into triggering and enhancing local accelerations governed by high values of .

3.3 The binary drive

Thus current knowledge backs our expectation that in BL Lac jets – by themselves in near-critical conditions – even modest perturbations that increase B can easily drive the structure beyond the threshold of macroscopic instabilities. Soon after, the conditions will progress toward strong kinetic instabilities if the macroscopic modes induce local torsion and shear of B. The ensuing scale for field shearing will be close to the longitudinal mode wavelengths experienced by the jet, that are of order of cm as numerically computed, e.g., by Mignone et al. 2013 (see also Striani et al. 2016).

In the case of a binary SMBH system, we envision as driving agent the dynamical perturbations of the magnetic jet configuration by the companion orbiting around the accretion disk, from which the main jet is launched. When the perturbations occur, the field B not only is locally compressed to larger strength a few times , but also is locally sheared so that regions of sharply different magnetic field lines are drawn together.

In a nutshell, binary perturbations drive macroscopic jet instabilities that offer many suitable locations for MR to occur and trigger effective electron acceleration; the ensuing emissions are detailed next.

4 Emissions and Their Interpretations

The basic feature to be expected from the binary drive of jet instabilities as discussed in Sect. 3 concerns a recurring, non-linear spike–trough pattern in the light curves, reflecting the stroke–relaxation nature of the dynamical drive along a Keplerian elliptical orbit.

4.1 Light curves

To be specific, we shall refer to the BL Lac object PG 1553+113 as an interesting candidate for recurrent non-linear gamma-ray emission as observed and discussed by Ackermann et al. (2015, henceforth A15). The source features in its continuously monitored gamma-ray emission high amplitude modulations that recur with an apparent period yr; at the source redshift (see Danforth et al. 2010) this corresponds to an intrinsic value yr. Three full (in fact, 3.5), consecutive periods are continuously covered by the existing gamma-ray observations of PG 1553+113 presented in A15; these data extend the dynamic range (number of cycles and amplitudes) of previous blazar periodicity analyses such as that of Vaughan et al. 2016. From the A15 analysis of PG 1553+113, a next gamma-ray peak is expected in early 2017.

Figure 2: The gamma-ray light curve of PG 1553+113 (green data points and error bars from A15) is compared with the predictions from our model, in which the gamma-ray flux is proportional to the gravitational force after (cf. Eq. A1 in the Appendix) along the elliptical orbit at a distance from the focus. We used the apparent period yr (see Sect 4.2), and three values for the eccentricity: (blue curve), (black curve) and (red curve). For the orbit would be circular and the light curve just flat, apart from minor random fluctuations.

In the Appendix we derive in detail a simple model for the gamma-ray light curve based on the gravitational modulation in a binary system, as impressed by the companion on the primary BH and its associated jet. Our result is presented here in Fig. 2, along with the 7-year data concerning PG 1553+113 from A15. It is seen that our model yields a non-linear pattern of the light curve with prominent peaks and wide troughs, that for an eccentricity is quantitatively very close to the overlaid data points.

4.2 Spectral distributions

The basic spectral change expected from compressed and sheared magnetic fields in the jet is twofold. To begin with, we expect an increased overall SSC emissivity, led by the optical S emission and followed after some weeks by IC radiation at GeV energies. Such a direct correlation is well compatible with the data as discussed by A15, but is also expected to be appreciably smeared out along the jet axis by outflow and long cooling times for both component sources (see Table 1 and its caption).

In detail, the implied enhancements of S, optical radiation and of its IC counterpart at GeV energies are shown by the spectra marked by red color in Fig. 3 with overlaid data. Such flux variations in both the R-band and the GeV component range by a factor around 3 from minima to maxima. In conclusion, the emissions in the optical, the soft X and the GeV bands appear to generally correlate along the spectra in shape and timing in the way described by the canonical SSC radiation model.

On the other hand, the behavior at keV energies of PG 1553+113 appears to stray outside the above framework, and has been reported as ”enigmatic” (see, e.g., Raiteri et al. 2915). For one, the increases during episodes of optical rise differ in softer X-ray ranges: they rise from a factor of 3 in the range 0.3 - 2 keV, but are up to factors 5 - 10 in the 2-10 keV range as reported by Aleksic et al. (2015) for the episode around MJD 56000. For two, other BL Lacs such as Mkn 421 appear to share with PG 1553+113 such keV features as missing correlation with neighboring X rays, and positive correlations with TeV emissions (Balokovic et al., 2016).

4.3 An additional spectral component in X rays

To interpret such a behavior, we recall the contents of Sect. 3. In particular, in Sect. 3.1 we noted that a BL Lac jet is normally metastable, so that a limited increase of can drive it first toward macroscopic, and then into kinetic instability conditions.

In a first stage the near-critical bulk magnetization is increased by a factor of a few, as deduced from the optical light curves enhanced during their 100-day rise. This implies a drift from marginal values of into the fully critical range for macroscopic instabilities. If these produce twisted and inverted lines a second stage arises, driven by increasing local magnetizations after Eq. 3. These enhance MRs and cause extensive particle acceleration. Based on current numerical simulations that are exploring different regimes of ion-electron plasmas (see Melzani et al. 2014, Kogan et al. 2015), we derive that the maximal electron energies can approach values for values of exceeding . Such enhanced accelerations will add a stronger and harder spectral component in the X-ray range around and above keV, as indeed detected in at least two episodes within the currently available data concerning PG 1553+113 (A15). The X-ray flux increase amounts to factors of 5-10, i.e., substantially larger than the related rise of optical and GeV emissions. Clearly, here one-zone, single population modelling of the X-ray emission together with the optical and GeV emissions would not be adequate. Our picture, on the other hand, implies additional particle acceleration triggered in nearby zones with increasing values of . The corresponding SEDs for one of such episodes are presented in Fig. 3.

Table 1: parameters for modelling the spectral states of PG 1553+113.

State Comp. B(G) R(cm)
low gray 15 0.25 0.6 2 3.9

enhanced soft
red 15 0.29 0.6 2 3.9
enhanced hard blue dashed 15 0.29 0.1 2 4.3

The Table summarizes the parameters used in our SSC spectral modelling of PG 1553+113 (see Fig. 3.) The Table provides values for the bulk Lorentz factor , the average magnetic field inducing Synchrotron emission, the source size , the electron energy distribution with normalization K, upper break at and lower bound , and power-law indices below and above the break. The associated cooling times for S and for IC (see Rybicki & Lightman 1979) read d .

Figure 3: Spectral energy distribution (SED) of the emissions by electrons accelerated in the jet of the BL Lac PG 1553+113 . Spectral data are from Ackermann et al. (2015) and Aleksic et al. (2015). The gray line shows the spectral state during normal conditions. The blue, dashed curve represents the harder, increased Synchrotron and IC components resulting from the extra electron acceleration discussed in Sect. 4.3. The total emission in the enhanced state is represented by the blue solid curve. EBL de-absorption in the TeV data after Aleksic et al. 2012.

5 Discussion and Conclusions

We have argued in Sect. 2 that dynamic events in Blazar jets - and specifically in BL Lac Objects - originate from an underlying binary system of two Super Massive Black Holes with total mass . Binary dynamics constitutes, as we stated in Sect. 2 and explained in Sect. 3, the prime mover for a chain of events in the jet: it recurrently drives bulk plasma instabilities in the metastable magnetized jet launched by the primary; in turn, these trigger local tearing instabilities and in situ electron acceleration to high energies corresponding to, and at times exceeding .

Such local instabilities and accelerations are known to arise in narrow layers of a collisionless plasma, where reconnections occur among adjacent -lines that are strongly sheared or even reversed, as we detailed in Sects. 3.1 and 3.2. Such processes are bound to be triggered during dynamic episodes of large-scale squeezing and shearing of a magnetized jet. A conducive context for such episodes to occur is provided by binary SMBH dynamics (as discussed in Sect. 3.3), that also provides the simplest interpretation for the recurrent, non-linear pattern of peaks and trough in the light curves of the BL Lac source PG 1553+113 as discussed in Sect. 4.1.

We have analyzed in Sect. 4.2 the different spectral states that arise in this source and exhibit repetition. They start with a recurrent increase of the optical Synchrotron and of the inverse Compton GeV radiations produced by the same population of energetic electrons. We have stressed in Sect. 4.3 the additional spectral component that has been observed in the keV range as the optical light curve approaches a peak. We traced it back to spreading accelerations triggered by additional reconnections stimulated by binary dynamics. We expect this component to correlate with TeV enhancement, as occurring in other BL Lacs.

So we conclude that at times two spectral states overlap in the emissions from PG 1553+113 . One is well represented as quasi-periodic emission, that envisages the optical - soft X-ray components and the GeV photons as linked components of the same SSC radiation process; the latter therefore should correlate with the former, after a lag of some weeks. The other, more energetically demanding state requires additional and higher electron accelerations that produce a harder X-ray peak, as a component (see Section 4.3) observed in two out of three monitored instances. this we expect to correlate with TeV enhancements, not necessarily with the optical – GeV radiations.

Next we summarize our overall conclusions, and specifically two points. First, we have described how a BL Lac jet, normally sitting on the brink of instability, can undergo recurrent stresses from the binary dynamics sufficient to stimulate cyclic emissions. In fact, increasing bulk magnetization leads to large scale instabilities that cascade down to small scale MRs and ultimately percolate to the kinetic level; there they accelerate the electrons that emit from the optical to the gamma-ray band, and sometimes beyond.

Second, the binary Keplerian motion drives strokes and allows relaxations corresponding to orbital arcs close to, and far away from the periastron, respectively. So it imprints a specific non-linearpattern of sharp peaks and wide troughs in the light curves at different frequencies; these are particularly sharp in gamma rays that have been continuously monitored from space. Such a pattern can provide diagnostics for orbital parameters, in particular for the orbital eccentricity as detailed in the Appendix, and used in Fig. 2 to derive for PG 1553+113 a value .

As to predictions and prospects, in the framework of the next gamma-ray peak expected for PG 1553 in early 2017 by A15, we stress the following points. (1) The new peak should appear in reasonable phase agreement (within the time binning) with the previous three occurrences. (2) Thus the optical emission should be monitored as continuously as possible across the peak phase range, with the help of telescopes strategically coordinated. The optical emission may be quite structured as shown by Figs. 1 and 2 of A15; in our picture we expect the main optical increase to begin several days to some weeks before the smoother and better sampled gamma-ray peak, which can sound a last call (as it were) for aimed X-ray observations. (3) In fact, the X-ray extra component in the keV range is expected to be triggered by expanding reconnections that occur just at, or on the approach to the periastron. Current data have the peak of additional X-ray emission corresponding in one case to the gamma-ray peak near MJD 57000 (to within the time binning); in another case, preceding he gamma-ray peak by 100 days. Our model contemplates enhanced keV emission somewhere between the rise of the optical and of the gamma-ray light curves, on account of the spread along the jet of the S and the IC main sources implied by their long cooling times (see Sect. 4.3, and Table 1).

It is interesting to recall that a similar radiative behavior to PG 1553+113 is shared by other BL Lacs, in particular by Mrk 421 with its long puzzling, X-ray outbursts correlated with TeV activity but not with optical and GeV emissions (se Raiteri et al. 2015, Balokovic et al. 2015). Some of these sources (listed, e.g., by Sandrinelli et al. 2014) have been proposed also as candidates for recurring behavior in the optical and gamma-ray bands.

We note that the dominance of possible BL Lacs’ periodicities over FSRQs’ in the set of current candidates may just reflect selection favoring the former on account of their generally smaller redshifts and lack of disturbing background from BLs and BBB pointed out in Sect. 1. On the other hand, intrinsic association of BL Lacs with binary sources is suggested by a number of circumstances. First, as a consequence of hierarchical galaxy formation, the BL Lacs as aged AGNs (Cavaliere & D’Elia 2002, Böttcher & Dermer 2002) must have witnessed several merging events involving their host nuclei. Thus, expanding on the pioneering suggestion by Begelman et al. (1980) such nuclei would be more massive and likely to have accrued several SMBHs to eventually form a dominant binary pair. In addition, the intrinsic meta-stability of BL Lac jets (see Sect. 3) makes the latter susceptible even to minor dynamic disturbances in forming a recurrent binary-driven source. Furthermore, the absence in BL Lac spectra of thermal features as recalled above points toward gas-poor conditions in the vicinity of the central SMBH binary; these prevent local energy dissipations from occurring, and allow for longer binary lifetimes implying higher intrinsic statistics (see Colpi 2914 and refs. therein), limited mainly by e.m. losses, and eventually cut off by a gravitational catastrophe (see Sect. 2).

Finally, we consider the A15 prediction concerning a next maximum of PG 1553+113 in early 2017, with a spread of some 20-40 days due to time binning; this would mark an obvious plus for the periodic view. On the quantitative side, extending the sequence of fully observed cycles from 3 to 4 would considerably strengthen the confidence level recalled in Sect. 4. Thus by early 2017 one had better be prepared with astrophysical analysis of the currently known cycles as undertaken here, and with coordinated observations of keV – TeV emissions. We add that evidence of any secondary jet would provide a tell-tale signature of binarity that deserves a close search, despite a number of difficulties: the expected weak power ratio even in the favorable case of both jets fed by accretion at a common fraction of their Eddington rates; the uncertainty about the jet relative angle; and the range in the time phase caused by precessional shifts.

These tasks are demanding, but the reward of confirmed evidence for binarity would be worth of keen observational efforts. It is genarally agree that understanding even a single SMBH binary would open a new window on the vexing issue of co-evolution of galaxies and their central BHs (see Lapi et al. 2014, Volonteri et al. 2015). In addition, a few confirmed SMBH binaries with assessed orbital parameters such as proposed here may be of keen value for planning searches of giant, slow bursts of gravitational waves with space interferometers such as eLISA (see eLISA Consortium, 2013). In fact, the path from host galaxy mergers to a sub-pc binary doomed to eventually coalesce, is complex (see Colpi 2015) and currently still fraught with a number of uncertain steps (”stalling” in particular), that not even massive numerical computations or simulations have been able to fix to now.

We are indebted to Bruno Coppi for enlightening discussions on magnetic field reconnecting structures in collisionless plasmas. We thank our referee for recommending a generally conservative attitude towards the quasi-periodicity of PG 1553+113), and for pressing us to specify the role of light curve profiles, as we did in the revised MS. Work partially supported by the ASI grant no. I/028/12/2.


Abbott et al., 2016a, Phys. Rev. Letters, 116, 061102

Abbott et al., 2016b, Phys. Rev. Letters, 116, 241103

Ackermann, M., Ajello, M., Albert, A., et al., 2015, ApJL, 813, 41 (A15)

Aleksic J., Alvarez, E.A., Antonelli, L.A. et al., 2012, ApJ, 748, 46

Aleksic J., Ansoldi, S., Antonelli, L.A. et al., 2015, MNRAS, 450, 4399

Balokovic, M., Paneque, D., Madejski, G. et al., 2016, 819, 156

Baumgarte, T.W. & Shapiro, S.L., 2011, Physics Today, 64, 32

Begelman, M., Blandford, R.D., & Rees M.J., 1980, Nature, 207, 307

Boettcher, M. & Dermer, C.D., 2002 ApJ, 564, 86

Burch, J.L., Torbert, R.B., Phan, T.D. et al., 2016, Science, 352, 1189

Cavaliere, A. & D’Elia, V., 2002, ApJ, 571, 226

Celotti, A & Ghisellini, G., 2008, MNRAS, 385, 283

Cohen, D.P, Romani, R.W., Filippenko, A.V., et al, 2014, ApJ, 797, 137

Colpi, M., 2014, Space Science Reviews, 183, 189

Coppi, B., 2016, Plasma Phys. Rep., 42, 383

Danforth, C.W., Keeney, B.A., Stocke, J.T., Shull, J.M., Yao, Y., 2010, ApJ, 720, 976

eLISA Consortium, 2013, The Gravitational Universe, arXiv:1305.5720v1

Ghisellini, G., 2016, arXiv:1609.08606

Hovatta, T., Valtaoja, E., Tornikoski, M., Lahteenmäki, 2009, A., 2009, A&A 494, 527

Kagan, D., Sironi, L., Cerutti, B., Giannios, D., 2015, SSRv, 191, 545

Kelly. B. C., Bechtold, J. & Siemiginowska, A., 2009, ApJ, 698, 895

Krawczynski, H., Hughes, S.B., Horan, D., et al., 2004, ApJ, 601, 151

Landau, L. & Lifshitz, E.M., 1976, Mechanics, Oxford, Pergamon Press

Lapi, A., Raimundo, S., Aversa, R., et al., 2014, ApJ, 782: 69

Maraschi, L., Ghisellini, G. & Celotti, A., 1992, ApJ, 397, L5

Melzani, M., Walder, R., Folini, D., Winisdoerffer, C., Favre, J.M. et al., 2014, A&A, 570A, 111

Mignone, A., Striani, E., Tavani, M., and Ferrari, A., 2013, MNRAS, 436, 1102

Nixon, C., & King, A., 2015, arXiv:1505.07827

Palenzuela, C., Garrett, T., Lehner, L., Liebling, S.L., 2010, PhRvD, 82, 404

Peterson, B.M., 1997, An Introduction to Active Galactic Nuclei, Cambridge Univ. Press

Peterson, B.M., 2006, The Broad-Line Region in Active Galactic Nuclei, Lect. Notes Phys. 693, 77 (Springer)

Raiteri, C. M., Stamerra, A., Villata, M., et al., 2015, MNRAS, 454, 353

Rybicki, G.H., Lightman, A.P., 1979, Radiative Processes in Astrophysics (Wiley)

Sandrinelli, A., Covino, S., Treves, A. 2014, A &A, 562, 79

Sironi,L. & Spitkovsky, A., 2014, ApJL, 783, 21

Striani, E., Mignone, A., Vaidya, B., Bodo, G., Ferrari, A., 2016, MNRAS, 462, 2970

Tavani, M., Vittorini, V., Cavaliere, A., 2015, ApJ, 814, 51

Vaughan, S., Uttley, P., Markowitz, A.G., et al., 2016, MNRAS, 461, 3145

Volonteri, M., Bogdanovic, T., Dotti, M., Colpi, M., 2015, Proc. XXIXth IAU General Assembly, arXiv:1509.09027v1

Yuan, Y, Nalewaiko, K., Zrake, J., East, W. E., Blandford, R., 2016, ApJ vol. 828; also preprint arXiv:1604.03179.

Appendix A Appendix: Deriving the Light Curves and the Corresponding Orbital Eccentricities

In Sect. 3 we consider a gravitationally bound binary system with separation of several milli-pc. Drawing on time-honored dynamical relations holding for Keplerian, elliptical orbits of the secondary (see, e.g., Landau & Lifshitz 1976, and their Figure 11) compared with the light curves observed by A15 (cf, their Figs. 1,   2), we detail how we deduce our expected light curve and the corresponding orbital eccentricity .

We shall follow two ways. For semplicity, in both we neglect the primary motion around the system’s center of mass, and in particular the size of its obit that scales as . We base on the gravitational stress from the secondary being proportional to the force between the two BHs, with enhanced emission also proportional to it.

1) The ratio of the force at periastron to that at the opposite vertex reads


In our approach, this also yields the ratio of gamma-ray fluxes at maxima and minima.

From the gamma-ray light-curve in A15 (their Fig. 1) we read out the observed value for the flux ratio between peaks and troughs. After Eq. 4 this corresponds to an orbital eccentricity .

In addition, in Fig. 4 we plot the full dependence of flux that we expect along the orbit at a distance from the focus, parameterized with three eccentricity values: . The data are overplotted, and agreement is found for

2) A second approach considers the time fraction along the orbit most affected by such an effect. For an estimate, we concentrate on the points on the elliptical orbit marking the intersections with the so-called latus rectum , i.e., the orbital distance corresponding to the ”true anomaly” angle , with at periastron. The time spent on the orbital arc circling the focus from through the periastron out to the opposite phase, is given by


Fig. 4 shows the ratio as a function of orbital eccentricities. From the observations recalled above we see that the strongly enhanced emission lasts for days, corresponding to a fraction . Interpreting this feature as caused by the gravitational perturbation of a secondary BH in an elliptical orbit, we deduce from Fig. 4 a value of the eccentricity .

Thus we conclude that on interpreting the repetitive pattern of enhanced high-energy emissions from PG 1553+113 in terms of stresses induced by binary dynamics, both the above approaches lead to evaluate for the orbital eccentricity a value close to .

Note that for a circular orbit the expected light curve is just flat; in fact, from Eq. 4 one obtains . Concerning Fig. 2 in the main text, we stress that the dynamically stimulated (as it were) emission maxima largely exceed the spontaneous random fluctuations best noticed in the troughs.

Figure 4: The ratio plotted as a function of the eccentricity for a Keplerian elliptical orbit. This quantity marks the strongly active fraction of the orbital period, that we take to be the interval between the intersections of the ”latus rectum”, each side of the periastron.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description