Search NASA⌕ Search

SEARCH · Search NASA

Results for “BL Lacertae objects”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

First-principles Prediction of X-Ray Polarization from Magnetic Reconnection in High-frequency BL Lacertae Objects

Relativistic magnetic reconnection is a potential particle acceleration mechanism for high-frequency BL Lac objects (HBLs). The Imaging X-ray Polarimetry Explorer (IXPE) scheduled to launch in 2021 has the capability to probe the evolution of magnetic field in HBLs, examining the magnetic reconnection scenario for the HBL flares. In this paper, we make the first attempt to self-consistently predict HBL X-ray polarization signatures arising from relativistic magnetic reconnection via combined particle-in-cell and polarized radiation transfer simulations. Here, we find that although the intrinsic optical and X-ray polarization degrees are similar on average, the X-ray polarization is much more variable in both the polarization degree and angle (PD and PA). Given the sensitivity of the IXPE, it may obtain one to a few polarization data points for one flaring event of nearby bright HBLs Mrk 421 and 501. However, it may not fully resolve the highly variable X-ray polarization. Due to temporal depolarization, where the integration of photons with variable polarization states over a finite period of time can lower the detected PD, the measured X-ray PD can be considerably lower than the optical counterpart or even undetectable. The lower X-ray PD than the optical thus can be a characteristic signature of relativistic magnetic reconnection. For very bright flares where the X-ray polarization is well resolved, relativistic magnetic reconnection predicts smooth X-ray PA swings, which originate from large plasmoid mergers in the reconnection region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Synchrotron self-Compton radiation from magnetically dominated turbulent plasmas in relativistic jets

ABSTRACT Relativistic jets launched by rotating black holes are powerful emitters of non-thermal radiation. Extraction of the rotational energy via electromagnetic stresses produces magnetically dominated jets, which may become turbulent. Studies of magnetically dominated plasma turbulence from first principles show that most of the accelerated particles have small pitch angles, i.e. the particle velocity is nearly aligned with the local magnetic field. We examine synchrotron self-Compton radiation from anisotropic particles in the fast cooling regime. The small pitch angles reduce the synchrotron cooling rate and promote the role of inverse Compton (IC) cooling, which can occur in two different regimes. In the Thomson regime, both synchrotron and IC components have soft spectra, νFν ∝ ν1/2. In the Klein–Nishina regime, synchrotron radiation has a hard spectrum, typically νFν ∝ ν, over a broad range of frequencies. Our results have implications for the modelling of BL Lacertae objects (BL Lacs) and gamma-ray bursts (GRBs). BL Lacs produce soft synchrotron and IC spectra, as expected when Klein–Nishina effects are minor. The observed synchrotron and IC luminosities are typically comparable, which indicates a moderate anisotropy with pitch angles θ ≳ 0.1. Rare orphan gamma-ray flares may be produced when θ ≪ 0.1. The hard spectra of GRBs may be consistent with synchrotron radiation when the emitting particles are IC cooling in the Klein–Nishina regime, as expected for pitch angles θ ∼ 0.1. Blazar and GRB spectra can be explained by turbulent jets with a similar electron plasma magnetization parameter, σe ∼ 104, which for electron–proton plasmas corresponds to an overall magnetization σ = (me/mp)σe ∼ 10.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

K2 Optical Emission from OJ 287 and Other γ-Ray Blazars on Hours-to-weeks Timescales from 2014 to 2018

We present second observations by K2 of OJ 287 and seven other γ-ray active galactic nuclei (AGNs) obtained in 2017–2018, second and third observations of the lobe-dominated, steep-spectrum quasar 3C 207, and observations of nine additional blazars not previously observed with K2. The AGNs were observed simultaneously with K2 and the Fermi Large Area Telescope for 51–81 days. Our full sample, observed in 2014–2018, contained 16 BL Lacertae objects (BL Lacs), nine flat-spectrum radio quasars (FSRQs), and four other γ-ray AGNs. Twelve BL Lacs and seven FSRQs exhibited fast, jagged light curves, while four BL Lacs and two FSRQs had slow, smooth light curves. Some objects changed their red-noise character significantly between repeated K2 observations. The optical characteristics of OJ 287 derived from the short-cadence K2 light curves changed between observations made before and after the predicted passage of the suspected secondary supermassive black hole through the accretion disk of the primary supermassive black hole. The average slopes of the periodogram power spectral densities of the BL Lacs' and FSRQs' light curves differed significantly, by ≈12%, with the BL Lac slopes being steeper, and a Kolmogorov–Smirnov test with a p-value of 0.039 indicates that these samples probably come from different populations; however, this result is not as strongly supported by power spectral response (PSRESP) analyses. Differences in the origin of the jets from the ergosphere or accretion disk in these two classes could produce such a disparity, as could different sizes or locations of emission regions within the jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Constraints on redshifts of blazars from extragalactic background light attenuation using Fermi - LAT data

ABSTRACT The extragalactic high-energy γ-ray sky is dominated by blazars, which are active galactic nuclei with their jets pointing towards us. Distance measurements are of fundamental importance yet for some of these sources are challenging because any spectral signature from the host galaxy may be outshone by the non-thermal emission from the jet. In this paper, we present a method to constrain redshifts for these sources that relies only on data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. This method takes advantage of the signatures that the pair-production interaction between photons with energies larger than approximately 10 GeV and the extragalactic background light leaves on γ-ray spectra. We find upper limits for the distances of 303 γ-ray blazars, classified as 157 BL Lacertae objects, 145 of uncertain class, and 1 flat-spectrum-radio quasar, whose redshifts are otherwise unknown. These derivations can be useful for planning observations with imaging atmospheric Cherenkov telescopes and also for testing theories of supermassive black hole evolution. Our results are applied to estimate the detectability of these blazars with the future Cherenkov Telescope Array, finding that at least 21 of them could be studied in a reasonable exposure of 20 h.

Astronomy & Astrophysics↗
DOE OSTI2021

Follow-up observations for IceCube-170922A: Detection of rapid near-infrared variability and intensive monitoring of TXS 0506+056

Abstract We present our follow-up observations to search for an electromagnetic counterpart of the IceCube high-energy neutrino IceCube-170922A. Monitoring observations of a likely counterpart, TXS 0506+056, are also described. First, we quickly took optical and near-infrared images of seven flat-spectrum radio sources within the IceCube error region right after the neutrino detection and found a rapid flux decline of TXS 0506+056 in Kanata/HONIR J-band data. Motivated by this discovery, intensive follow-up observations of TXS 0506+056 were continuously performed, including our monitoring imaging observations, spectroscopic observations, and polarimetric observations in optical and near-infrared wavelengths. TXS 0506+056 showed a large-amplitude (∼1.0 mag) variability in a time scale of several days or longer, although no significant variability was detected in a time scale of a day or shorter. TXS 0506+056 also showed a bluer-when-brighter trend in optical and near-infrared wavelengths. Structure functions of the variabilities were examined and indicate that TXS 0506+056 is not a special blazar in terms of optical variability. Polarization measurement results of TXS 0506+056 are also discussed.

(galaxies:) BL Lacertae objects: general↗
DOE OSTI2023

The spectra of IceCube Neutrino (SIN) candidate sources – III. Optical spectroscopy and source characterization of the full sample

A correlation between astrophysical high-energy neutrinos and blazars has been suggested by various authors. In particular, a likely association between IceCube events and intermediate- and high-energy peaked BL Lac objects has led to a sample of 47 objects having a high probability of being neutrino sources. In the first paper of this series we reported optical spectroscopy of 17 objects, which together with data taken from the literature covered 80 per cent of the sample. Here, we present spectroscopy obtained at large aperture telescopes of a further 17 objects (plus four additional targets), which completes the sample coverage. For 12 objects we are able to determine the redshift (0.07 < z < 1.6), while for the others we set a lower limit on it, based on either the robust detection of intervening absorption systems or on an estimation derived from the absence of spectral signatures of the host galaxy. With these new data we expand and reinforce the main results of our previous papers, namely the fact that in terms of their broad-band properties our sources appear to be indistinguishable from the rest of the blazar population and the relatively large (>34 per cent and possibly as high as 80 per cent) fraction of masquerading BL Lac objects, for which the low equivalent width of the emission lines is due to the brightness of the boosted continuum, rather than being an intrinsic property, in our sample.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Possible jet reorientation in low-frequency radio structures of blazars

ABSTRACT Blazars are active galactic nuclei (AGNs) observed at small viewing angles to the observer’s line of sight; hence most of them are compact sources in the radio regime. At low radio frequencies, we may observe additional extended structures, considered to be lobes propagating directly towards Earth. However, there are rare cases of blazars whose morphologies are interpreted in the framework of an episodic jet activity with a reorientation of the jet axis. In this paper, we aim to find other candidates for sources of this kind and present Low Frequency Array (LOFAR) observations of seven blazars with rare morphologies composed of kiloparsec-scale jets and elongated multihundred kiloparsec-scale features. Based on their radio images and physical parameters, we propose that jets could change the direction of propagation out of the sky plane during the lifetime of at least some of these sources and with such we may observe morphologically different AGN types simultaneously, similar to the previously discovered blazar SBS B1646+499. Several plausible origin scenarios are presented, including fast realignment and precession. We also derive the 144 MHz fluxes of the compact and extended structures, the total radio luminosities, the spectral indices for the diffuse haloes, and the core-dominance parameters. We also discuss the possibility that these objects are actually giant blazars with rare Fanaroff–Riley I (FR I) morphology. The previously discovered blazar SBS B1646+499 should be classified as a giant independently of its inclination with regard to the observer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Evidence of underdeveloped torus and broad-line region of weak emission line quasars based on their spectral energy distribution

To unravel the dominant cause of the weak emission line in a subset of optically selected radio-quiet ‘weak emission line quasars’ (WLQs), we have investigated the possibility of an underdeveloped broad line region (BLR). For this, we have modelled spectral energy distributions (SED) of 61 WLQs by using their optical and infrared (IR) photometric observations from Sloan Digital Sky Survey (SDSS) and WISE, respectively. SED fit consists of various emission components, including the luminosity from the dusty torus (L tor ). For comparison with the normal quasar, we have used a control sample of 55 Quasi-Stellar Objects (QSOs) for each WLQ matching in emission redshift and SDSS r-band. Based on our measurement of L tor , we found a decrement of 42 ± 2 per cent in IR luminosity in WLQs with respect to the control sample of normal QSOs. Using L tor /L bol as the measure of torus covering factor (CF tor ), we found a similar decrement in WLQ covering factor, with their CF tor distribution being significantly different with respect to the normal QSOs with a KS test P null of 4.27 × 10 −14 . As dusty torus and BLR covering factors are expected to be of a similar order in active galactic nucleus (AGN), our results suggest that the BLR in the WLQs is underdeveloped and could be a dominant cause of the weakness of their emission line. As a result, our analysis gives support to the models of WLQs based on the evolution scenario being in an early stage of AGNs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The spectra of IceCube neutrino (SIN) candidate sources – II. Source characterization

ABSTRACT Eight years after the first detection of high-energy astrophysical neutrinos by IceCube, we are still almost clueless as regards to their origin, although the case for blazars being neutrino sources is getting stronger. After the first significant association at the $3\!-\!3.5\, \sigma$ level in time and space with IceCube neutrinos, i.e. the blazar TXS 0506+056 at z = 0.3365, some of us have in fact selected a unique sample of 47 blazars, out of which ∼16 could be associated with individual neutrino track events detected by IceCube. Building upon our recent spectroscopy work on these objects, here we characterize them to determine their real nature and check if they are different from the rest of the blazar population. For the first time we also present a systematic study of the frequency of masquerading BL Lacs, i.e. flat-spectrum radio quasars with their broad lines swamped by non-thermal jet emission, in a γ-ray- and IceCube-selected sample, finding a fraction >24 per cent and possibly as high as 80 per cent. In terms of their broad-band properties, our sources appear to be indistinguishable from the rest of the blazar population. We also discuss two theoretical scenarios for neutrino emission, one in which neutrinos are produced in interactions of protons with jet photons and one in which the target photons are from the broad-line region. Both scenarios can equally account for the neutrino–blazar correlation observed by some of us. Future observations with neutrino telescopes and X-ray satellites will test them out.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Peculiar Variable X-Ray Spectrum of the Active Galactic Nucleus PKS 2005–489

PKS 2005-489 is a well-known, bright southern BL Lac object that has been detected up to TeV energies. In a low-flux state it exhibits the expected multiwavelength double-peaked spectrum in the radio- γ -ray band. The high-flux state shows extreme flux variations in the X-ray band with a hardening as well as a peculiar curved feature in the spectrum. Thus far, PKS 2005-489 is the only source to exhibit such a feature. To study the X-ray variability further, we obtained the first hard X-ray spectrum of the source with the Nuclear Spectroscopic Telescope Array. We compare quasi-simultaneous radio, optical, UV, soft and hard X-ray, and γ -ray data of PKS 2005-489 to archival data in order to study its broadband behavior. We find a very consistent quiet state in the spectral energy distribution, with little variation in spectral shape or flux between the 2012 and 2020 data. A possible explanation for the peculiar X-ray spectrum in the flaring state is an additional component in the jet, possibly accelerated via magnetic reconnection, that is not co-spatial to the low-flux state emission region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Fourth Catalog of Active Galactic Nuclei Detected by the Fermi Large Area Telescope

We report the fourth catalog of active galactic nuclei (AGNs) detected by the Fermi Gamma-ray Space Telescope Large Area Telescope (4LAC) between 2008 August 4 and 2016 August 2 contains objects located at high Galactic latitudes ( ). It includes 85% more sources than the previous 3LAC catalog based on 4 yr of data. AGNs represent at least 79% of the high-latitude sources in the fourth Fermi -Large Area Telescope Source Catalog (4FGL), which covers the energy range from 50 MeV to 1 TeV. In addition, gamma-ray AGNs are found at low Galactic latitudes. Most of the 4LAC AGNs are blazars (98%), while the remainder are other types of AGNs. The blazar population consists of 24% Flat Spectrum Radio Quasars (FSRQs), 38% BL Lac-type objects, and 38% blazar candidates of unknown types (BCUs). On average, FSRQs display softer spectra and stronger variability in the gamma-ray band than BL Lacs do, confirming previous findings. All AGNs detected by ground-based atmospheric Cerenkov telescopes are also found in the 4LAC.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Multiwavelength Observations of the Blazar VER J0521+211 during an Elevated TeV Gamma-Ray State

We report on a long-lasting, elevated gamma-ray flux state from VER J0521+211 observed by VERITAS, MAGIC, and Fermi-LAT in 2013 and 2014. The peak integral flux above 200 GeV measured with the nightly binned light curve is (8.8 ± 0.4) × 10 –7 photons m –2 s –1 , or ~37% of the Crab Nebula flux. Multiwavelength observations from X-ray, UV, and optical instruments are also presented. A moderate correlation between the X-ray and TeV gamma-ray fluxes was observed, and the X-ray spectrum appeared harder when the flux was higher. Using the gamma-ray spectrum and four models of the extragalactic background light (EBL), a conservative 95% confidence upper limit on the redshift of the source was found to be z ≤ 0.31. Unlike the gamma-ray and X-ray bands, the optical flux did not increase significantly during the studied period compared to the archival low-state flux. The spectral variability from optical to X-ray bands suggests that the synchrotron peak of the spectral energy distribution (SED) may become broader during flaring states, which can be adequately described with a one-zone synchrotron self-Compton model varying the high-energy end of the underlying particle spectrum. The synchrotron peak frequency of the SED and the radio morphology of the jet from the MOJAVE program are consistent with the source being an intermediate-frequency-peaked BL Lac object.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Simultaneous observations of the blazar PKS 2155-304 from ultra-violet to TeV energies

Here we report the results of the first ever contemporaneous multi-wavelength observation campaign on the BL Lac object PKS 2155−304 involving Swift , NuSTAR , Fermi -LAT, and H.E.S.S. The use of these instruments allows us to cover a broad energy range, which is important for disentangling the different radiative mechanisms. The source, observed from June 2013 to October 2013, was found in a low flux state with respect to previous observations but exhibited highly significant flux variability in the X-rays. The high-energy end of the synchrotron spectrum can be traced up to 40 keV without significant contamination by high-energy emission. A one-zone synchrotron self-Compton model was used to reproduce the broadband flux of the source for all the observations presented here but failed for previous observations made in April 2013. A lepto-hadronic solution was then explored to explain these earlier observational results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Long-term multi-wavelength study of 1ES 0647+250

The BL Lac object 1ES 0647+250 is one of the few distant γ-ray emitting blazars detected at very high energies (VHEs; &100 GeV) during a non-flaring state. It was detected with the MAGIC telescopes during a period of low activity in the years 2009-2011 as well as during three flaring activities in the years 2014, 2019, and 2020, with the highest VHE flux in the last epoch. An extensive multi-instrument data set was collected as part of several coordinated observing campaigns over these years. We aim to characterise the long-term multi-band flux variability of 1ES 0647+250, as well as its broadband spectral energy distribution (SED) during four distinct activity states selected in four different epochs, in order to constrain the physical parameters of the blazar emission region under certain assumptions. We evaluated the variability and correlation of the emission in the different energy bands with the fractional variability and the Z transformed discrete correlation function, as well as its spectral evolution in X-rays and γ rays. Owing to the controversy in the redshift measurements of 1ES 0647+250 reported in the literature, we also estimated its distance in an indirect manner through a comparison of the GeV and TeV spectra from simultaneous observations with Fermi-LAT and MAGIC during the strongest flaring activity detected to date. Moreover, we interpret the SEDs from the four distinct activity states within the framework of one-component and two-component leptonic models, proposing specific scenarios that are able to reproduce the available multi-instrument data. We find significant long-term variability, especially in X-rays and VHE γ rays. Furthermore, significant (3-4σ) correlations were found between the radio, optical, and high-energy (HE) γ-ray fluxes, with the radio emission delayed by about ~400 days with respect to the optical and γ-ray bands. The spectral analysis reveals a harder-when-brighter trend during the non-flaring state in the X-ray domain. However, no clear patterns were observed for either the enhanced states or the HE (30 MeV < E < 100 GeV) and VHE γ-ray emission of the source. The indirect estimation of the redshift yielded a value of z = 0.45 ± 0.05, which is compatible with some of the values reported in the literature. The SEDs related to the low-activity state and the three flaring states of 1ES 0647+250 can be described reasonably well with the both one-component and two-component leptonic scenarios. However, the long-term correlations indicate the need for an additional radio-producing region located about 3.6 pc downstream from the gamma-ray producing region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Investigation of the correlation patterns and the Compton dominance variability of Mrk 421 in 2017

We present a detailed characterisation and theoretical interpretation of the broadband emission of the paradigmatic TeV blazar Mrk 421, with a special focus on the multi-band flux correlations. The dataset has been collected through an extensive multi-wavelength campaign organised between 2016 December and 2017 June. The instruments involved are MAGIC, FACT, Fermi-LAT, Swift, GASP-WEBT, OVRO, Medicina, and Metsähovi. Additionally, four deep exposures (several hours long) with simultaneous MAGIC and NuSTAR observations allowed a precise measurement of the falling segments of the two spectral components. The very-high-energy (VHE; E> 100 GeV) gamma rays and X-rays are positively correlated at zero time lag, but the strength and characteristics of the correlation change substantially across the various energy bands probed. The VHE versus X-ray fluxes follow different patterns, partly due to substantial changes in the Compton dominance for a few days without a simultaneous increase in the X-ray flux (i.e., orphan gamma-ray activity). Studying the broadband spectral energy distribution (SED) during the days including NuSTAR observations, we show that these changes can be explained within a one-zone leptonic model with a blob that increases its size over time. The peak frequency of the synchrotron bump varies by two orders of magnitude throughout the campaign. Our multi-band correlation study also hints at an anti-correlation between UV-optical and X-ray at a significance higher than 3σ. A VHE flare observed on MJD 57788 (2017 February 4) shows gamma-ray variability on multi-hour timescales, with a factor ten increase in the TeV flux but only a moderate increase in the keV flux. The related broadband SED is better described by a two-zone leptonic scenario rather than by a one-zone scenario. We find that the flare can be produced by the appearance of a compact second blob populated by high energetic electrons spanning a narrow range of Lorentz factors, from γ' min =2x10 4 to γ' max =6x10 5 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Probing circular polarization and magnetic field structure in active galactic nuclei

The composition and magnetic field morphology of relativistic jets can be studied using circular polarization (CP) measurement. Recent three-dimensional relativistic magnetohydrodynamic (3D RMHD) simulations coupled with radiative transfer (RT) calculations make strong predictions about the level (and morphology) of the jet’s CP emission. These simulations show that the sign of CP and the electric vector position angle (EVPA) are both sensitive to the jet’s magnetic field morphology within the radio core. We probe this theory by exploring whether the jet’s radio core EVPA orientation is consistent with the observed sign of the core CP in deep full-track polarimetric observations. Based on a selection of sources from earlier MOJAVE observations, we aim to probe the nature of linear polarization (LP) and CP in the innermost regions of jets from a small sample of nine blazars. This sample includes sources that have exhibited: (i) positive CP; (ii) negative CP; or (iii) positive and negative CP simultaneously in the radio core region. By coupling deep polarimetric observations of a carefully selected sample of blazars with state-of-the-art RMHD and RT calculations, we hope to gain a deeper understanding of the physics of blazar jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The relation between optical and γ-ray emission in BL Lac sources

ABSTRACT The relativistic jets produced by some Active Galactic Nuclei (AGNs) are among the most efficient persistent sources of non-thermal radiation and represent an ideal laboratory for studying high-energy interactions. In particular, when the relativistic jet propagates along the observer’s line of sight, the beaming effect produces dominant signatures in the observed spectral energy distribution (SED), from the radio domain up to the highest energies, with the further possibility of resulting in radiation-particle multimessenger associations. In this work, we investigate the relationships between the emission of γ rays and the optical spectra of a sample of AGN, selected from BL Lac sources detected by the Fermi Large Area Telescope (Fermi-LAT). We find that there is a close relationship between the optical and γ-ray spectral indices. Despite all the limitations due to the non-simultaneity of the data, this observation strongly supports a substantial role of Synchrotron-Self Compton (SSC) radiation in a single zone leptonic scenario for most sources. This result simplifies the application of theoretical models to explore the physical parameters of the jets in this type of sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Periodic Gamma-Ray Modulation of the Blazar PG 1553+113 Confirmed by Fermi-LAT and Multiwavelength Observations

A 2.1 yr periodic oscillation of the gamma-ray flux from the blazar PG 1553+113 has previously been tentatively identified in ∼7 yr of data from the Fermi Large Area Telescope. After 15 yr of Fermi sky-survey observations, doubling the total time range, we report >7 cycle gamma-ray modulation with an estimated significance of 4σ against stochastic red noise. Independent determinations of oscillation period and phase in the earlier and the new data are in close agreement (chance probability <0.01). Pulse timing over the full light curve is also consistent with a coherent periodicity. Multiwavelength new data from Swift X-Ray Telescope, Burst Alert Telescope, and UVOT, and from KAIT, Catalina Sky Survey, All-Sky Automated Survey for Supernovae, and Owens Valley Radio Observatory ground-based observatories as well as archival Rossi X-Ray Timing Explorer satellite-All Sky Monitor data, published optical data of Tuorla, and optical historical Harvard plates data are included in our work. Optical and radio light curves show clear correlations with the gamma-ray modulation, possibly with a nonconstant time lag for the radio flux. We interpret the gamma-ray periodicity as possibly arising from a pulsational accretion flow in a sub-parsec binary supermassive black hole system of elevated mass ratio, with orbital modulation of the supplied material and energy in the jet. Other astrophysical scenarios introduced include instabilities, disk and jet precession, rotation or nutation, and perturbations by massive stars or intermediate-mass black holes in polar orbit.

79 ASTRONOMY AND ASTROPHYSICS↗