Empirical multi-wavelength prediction method for very high energy gamma-ray emitting BL Lacertae objects

Aim: We have collected the most complete multi-wavelength (6.0-6.0 x 10(18) cm) dataset of very high energy (VHE) gamma-ray emitting (TeV) BL Lacs, which are the most numerous extragalactic VHE sources. Using significant correlations between different bands, we aim to identify the best TeV BL Lac candidates that can be discovered by the current and next generation of imaging air Cherenkov telescopes.Methods: We formed five datasets from lower energy data, i.e. radio, mid-infrared, optical, X-rays, and GeV gamma-ray, and five VHE gamma-ray datasets to perform a correlation study between di ff erent bands and to construct the prediction method. The low energy datasets were averaged for individual sources, while the VHE gamma-ray data were divided into subsets according to the flux state of the source. We then looked for significant correlations and determined their best-fit parameters. Using the best-fit parameters we predicted the level of VHE gamma-ray flux for a sample of 182 BL Lacs, which have not been detected at TeV energies. We identified the most promising TeV BL Lac candidates based on the predicted VHE gamma-ray flux for each source.Results: We found 14 significant correlations between radio, mid-infrared, optical, gamma-ray, and VHE gamma-ray bands. The correlation between optical and VHE gamma-ray luminosity is established for the first time. We attribute this to the more complete sample and more accurate handling of host galaxy flux in our work. We found nine BL Lac candidates whose predicted VHE gamma-ray flux is high enough for detection in less than 25 h with current imaging air Cherenkov telescopes

Quasi-periodicities of BL Lacertae objects

We review the reports of possible year-long quasi-periodicities of BL Lac objects in the gamma-ray and optical bands, and present a homogeneous time analysis of the light curves of PKS2155-304, PG1553+113, and BL Lac. Based on results from a survey covering the entire Fermi gamma-ray sky we have estimated the fraction of possible quasi-periodic BL Lac objects. We compared the cyclical behaviour in BL Lac objects with that derived from the search of possible optical periodicities in quasars, and find that at z less than or similar to 1 the cosmic density of quasi-periodic BL Lac objects is larger than that of quasi-periodic quasars. If the BL Lac quasi-periodicities were due to a supermassive binary black hole (SBBH) scenario, there could be a tension with the upper limits on the gravitational wave background measured by the pulsar timing array. The argument clearly indicates the difficulties of generally associating quasi-periodicities of BL Lac objects with SBBHs

Multiwavelength Variability Power Spectrum Analysis of the Blazars 3C 279 and PKS 1510-089 on Multiple Timescales

We present the results of variability power spectral density (PSD) analysis using multiwavelength radio to GeV gamma-ray light curves covering timescales of decades/years to days/minutes for the blazars 3C 279 and PKS 1510-089. The PSDs are modeled as single power laws, and the best-fit spectral shape is derived using the "power spectral response" method. With more than 10 yr of data obtained with weekly/daily sampling intervals, most of the PSDs cover similar to 2-4 decades in temporal frequency; moreover, in the optical band, the PSDs cover similar to 6 decades for 3C 279 due to the availability of intranight light curves. Our main results are the following: (1) on timescales ranging from decades to days, the synchrotron and the inverse-Compton spectral components, in general, exhibit red-noise (slope similar to 2) and flicker-noise (slope similar to 1) type variability, respectively; (2) the slopes of gamma-ray variability PSDs obtained using a 3 hr integration bin and 3 weeks total duration exhibit a range between similar to 1.4 and similar to 2.0 (mean slope = 1.60 +/- 0.70), consistent within errors with the slope on longer timescales; (3) comparisons of fractional variability indicate more power on timescales <= 100 days at gamma-ray frequencies compared to longer wavelengths, in general (except between the gamma-ray and optical wavelengths for PKS 1510-089); (4) the normalization of intranight optical PSDs for 3C 279 appears to be a simple extrapolation from longer timescales, indicating a continuous (single) process driving the variability at optical wavelengths; and (5) the emission at optical/infrared wavelengths may involve a combination of disk and jet processes for PKS 1510-089

A Decade of Multiwavelength Observations of the TeV Blazar 1ES 1215+303: Extreme Shift of the Synchrotron Peak Frequency and Long-term Optical-Gamma-Ray Flux Increase

Blazars are known for their variability on a wide range of timescales at all wavelengths. Most studies of TeV gamma-ray blazars focus on short timescales, especially during flares. With a decade of observations from the Fermi-LAT and VERITAS, we present an extensive study of the long-term multiwavelength radio-to-gamma-ray flux-density variability, with the addition of a couple of short-time radio-structure and optical polarization observations of the blazar 1ES 1215+303 (z = 0.130), with a focus on its gamma-ray emission from 100 MeV to 30 TeV. Multiple strong GeV gamma-ray flares, a long-term increase in the gamma-ray and optical flux baseline, and a linear correlation between these two bands are observed over the ten-year period. Typical HBL behaviors are identified in the radio morphology and broadband spectrum of the source. Three stationary features in the innermost jet are resolved by Very Long Baseline Array at 43.1, 22.2, and 15.3 GHz. We employ a two-component synchrotron self-Compton model to describe different flux states of the source, including the epoch during which an extreme shift in energy of the synchrotron peak frequency from infrared to soft X-rays is observed

Authenticating the Presence of a Relativistic Massive Black Hole Binary in OJ 287 Using Its General Relativity Centenary Flare: Improved Orbital Parameters

Results from regular monitoring of relativistic compact binaries like PSR 1913+16 are consistent with the dominant (quadrupole) order emission of gravitational waves (GWs). We show that observations associated with the binary black hole (BBH) central engine of blazar OJ 287 demand the inclusion of gravitational radiation reaction effects beyond the quadrupolar order. It turns out that even the effects of certain hereditary contributions to GW emission are required to predict impact flare timings of OJ 287. We develop an approach that incorporates this effect into the BBH model for OJ 287. This allows us to demonstrate an excellent agreement between the observed impact flare timings and those predicted from ten orbital cycles of the BBH central engine model. The deduced rate of orbital period decay is nine orders of magnitude higher than the observed rate in PSR 1913+16, demonstrating again the relativistic nature of OJ 287's central engine. Finally, we argue that precise timing of the predicted 2019 impact flare should allow a test of the celebrated black hole "no-hair theorem" at the 10% level

Stochastic Modeling of Multiwavelength Variability of the Classical BL Lac Object OJ287 on Timescales Ranging from Decades to Hours

We present the results of our power spectral density analysis for the BL Lac object OJ 287, utilizing the Fermi-LAT survey at high-energy.-rays, Swift-XRT in X-rays, several ground-based telescopes and the Kepler satellite in the optical, and radio telescopes at GHz frequencies. The light curves are modeled in terms of continuous-time autoregressive moving average (CARMA) processes. Owing to the inclusion of the Kepler data, we were able to construct for the first time the optical variability power spectrum of a blazar without any gaps across similar to 6 dex in temporal frequencies. Our analysis reveals that the radio power spectra are of a colored-noise type on timescales ranging from tens of years down to months, with no evidence for breaks or other spectral features. The overall optical power spectrum is also consistent with a colored noise on the variability timescales ranging from 117 years down to hours, with no hints of any quasi-periodic oscillations. The X-ray power spectrum resembles the radio and optical power spectra on the analogous timescales ranging from tens of years down to months. Finally, the.-ray power spectrum is noticeably different from the radio, optical, and X-ray power spectra of the source: we have detected a characteristic relaxation timescale in the Fermi-LAT data, corresponding to similar to 150 days, such that on timescales longer than this, the power spectrum is consistent with uncorrelated (white) noise, while on shorter variability timescales there is correlated (colored) noise

Observation of the black widow B1957+20 millisecond pulsar binary system with the MAGIC telescopes

B1957+20 is a millisecond pulsar located in a black-widow-type compact binary system with a low-mass stellar companion. The interaction of the pulsar wind with the companion star wind and/or the interstellar plasma is expected to create plausible conditions for acceleration of electrons to TeV energies and subsequent production of very high-energy gamma-rays in the inverse Compton process. We performed extensive observations with the Major Atmospheric Gamma Imaging Cherenkov Telescopes (MAGIC) telescopes of B1957+20. We interpret results in the framework of a few different models, namely emission from the vicinity of the millisecond pulsar, the interaction of the pulsar and stellar companion wind region or bow shock nebula. No significant steady very high-energy gamma-ray emission was found. We derived a 95 per cent confidence level upper limit of 3.0 x 10(-12) cm(-2) s(-1) on the average gamma-ray emission from the binary system above 200 GeV. The upper limits obtained with the MAGIC constrain, for the first time, different models of the high-energy emission in B1957+20. In particular, in the inner mixed wind nebula model with mono-energetic injection of electrons, the acceleration efficiency of electrons is constrained to be below similar to 2-10 per cent of the pulsar spin-down power. For the pulsar emission, the obtained upper limits for each emission peak are well above the exponential cut-off fits to the Fermi-LAT data, extrapolated to energies above 50 GeV. The MAGIC upper limits can rule out a simple power-law tail extension through the sub-TeV energy range for the main peak seen at radio frequencies

Constraining Lorentz Invariance Violation Using the Crab Pulsar Emission Observed up to TeV Energies by MAGIC

Spontaneous breaking of Lorentz symmetry at energies on the order of the Planck energy or lower is predicted by many quantum gravity theories, implying non-trivial dispersion relations for the photon in vacuum. Consequently, gamma-rays of different energies, emitted simultaneously from astrophysical sources, could accumulate measurable differences in their time of flight until they reach the Earth. Such tests have been carried out in the past using fast variations of gamma-ray flux from pulsars, and more recently from active galactic nuclei and gamma-ray bursts. We present new constraints studying the gamma-ray emission of the galactic Crab Pulsar, recently observed up to TeV energies by the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) collaboration. A profile likelihood analysis of pulsar events reconstructed for energies above 400. GeV finds no significant variation in arrival time as their energy increases.. Ninety-five percent CL limits are obtained on the effective Lorentz invariance violating energy scale at the level of E-QG1 > 5.5 x 10(17) GeV (4.5 x 10(17) GeV) for a linear, and E-QG2 > 5.9 x 10(10) GeV (5.3 x 10(10) GeV) for a quadratic scenario, for the subluminal and the superluminal cases, respectively. A substantial part of this study is dedicated to calibration of the test statistic, with respect to bias and coverage properties. Moreover, the limits take into account systematic uncertainties, which are found to worsen the statistical limits by about 36%-42%. Our constraints would have been much more stringent if the intrinsic pulse shape of the pulsar between 200 GeV and 400 GeV was understood in sufficient detail and allowed inclusion of events well below 400 GeV

Search for very high-energy gamma-ray emission from the microquasar Cygnus X-1 with the MAGIC telescopes

The microquasar Cygnus X-1 displays the two typical soft and hard X-ray states of a black hole transient. During the latter, Cygnus X-1 shows a one-sided relativistic radio-jet. Recent detection of the system in the high energy (HE; E greater than or similar to 60 MeV) gamma-ray range with FermiLAT associates this emission with the outflow. Former MAGIC observations revealed a hint of flaring activity in the very high-energy (VHE; E greater than or similar to 100 GeV) regime during this X-ray state. We analyse similar to 97 h of Cygnus X-1 data taken with the MAGIC telescopes between July 2007 and October 2014. To shed light on the correlation between hard X-ray and VHE gamma rays as previously suggested, we study each main X-ray state separately. We perform an orbital phase-folded analysis to look for variability in the VHE band. Additionally, to place this variability behaviour in a multiwavelength context, we compare our results with Fermi-LAT, AGILE, Swift-BAT, MAXI, RXTE-ASM, AMI and RATAN-600 data. We do not detect Cygnus X-1 in the VHE regime. We establish upper limits for each X-ray state, assuming a power-law distribution with photon index Gamma = 3.2. For steady emission in the hard and soft X-ray states, we set integral upper limits at 95 per cent confidence level for energies above 200 GeV at 2.6 x 10(-12) photons cm(-2) s(-1) and 1.0 x 10(-11) photons cm(-2) s(-1), respectively. We rule out steady VHE gamma-ray emission above this energy range, at the level of theMAGIC sensitivity, originating in the interaction between the relativistic jet and the surrounding medium, while the emission above this flux level produced inside the binary still remains a valid possibility

Incremental Fermi Large Area Telescope Fourth Source Catalog

We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi Large Area Telescope (LAT) catalog of gamma-ray sources. Based on the first 12 years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss the degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources