Effects of Redshift on the Classifying Criteria of BL Lacertae Objects

We have collected a sample of 70 BL Lacs （33 radio-selected BL Lacs and 37 X-ray selected BL Lacs） with multi-waveband data for investigating the classifying criteria of BL Lacertae Objects. For each source, we estimate its luminosities in radio, optical and X-ray, the broad-band spectral index from radio to X-ray and the peak frequency of the synchrotron emission, and make a statistical analysis of the data obtained. Our main results are as follows： （1） The broad-band spectral index and the peak frequency have no correlation with the redshift, while they are inversely correlated with each other and they could be regarded as equivalent classifying criteria of BL Lac objects. （2） There are significant effects of the luminosity/redshift relation on the observed luminosity distribution in our sample, hence, if the radio luminosity is to be used as a classifying criterion of BL Lac objects, it should not be regarded as equivalent to the broad-band spectral index or the peak frequency. （3） Our resuits supply a specific piece of evidence for the suggestion that the use of luminosities always introduces a redshift bias to the data and show that the location of the peak frequency is not always linked to the luminosity of any wave band.

The Binary Black Hole Scenario for the BL Lacertae Object AO 0235+16

Recent analysis of the long term radio light curve of the extremely variableBL Lacertae object AO 0235+16 by Raiteri et al. have revealed the presence of recurrent outburstswith a period of ～ 5.7 +- 0.5 yr. Periodicity analysis of the optical light curve also showsevidence for a shorter period. Here we discuss whether such a behavior can be explained by a binaryblack hole model where the accretion disk of one of the supermassive black holes is precessing dueto the tidal effects of the companion. We estimate the mass of the accreting hole and analyzeconstraints on the secondary mass and the orbital parameters of the system. It is possible toprovide a viable interpretation of the available multiwavelength data.

On the weak Lyα emission line of BL Lacertae objects

Recent studies reveal that weak Lyα emission line may be a ubiquitous feature of nearby BL Laceratae objects(BL Lacs). We present a survey of the Lyα emission lines in BL Lacs, with a focus on data obtained by the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. Among the11 selected targets, seven show a clear detection of the intrinsic Lyα emission line, while one manifests a marginal detection. Two of these intrinsic Lyα lines are newly detected in this work(FBQS J1217+3007 and3 C 66 A). Most selected BL Lacs are either high or intermediate-frequency peaked BL Lacs, with the only exception being one low-frequency peaked BL Lac. A strong anti-correlation between the line equivalent width and the continuum luminosity was found. The possible sources of ionizing photons that are responsible for production of the Lyα emission line are studied based on a widely accepted detailed accretion-jet model. It is found that the clouds of the broad line region are unlikely to be located in the jet cone region with an inclination angle of less than 10 degrees. Contributions from the jet and disk to the emission lines are found to be comparable. Additionally, a possible way to constrain the accretion rate from the emission line is proposed.

BL Lacertae objects and the extragalactic γ-ray background

A tight correlation between γ-ray and radio emission is found for a sample of BL Lacertae （BL Lac） objects detected by the Fermi Gamma-ray Space Telescope （Fermi） and the Energetic Gamma-Ray Experiment Telescope （EGRET）. The γ-ray emission of BL Lac objects exhibits strong variability, and the detection rate of γ-ray BL Lac objects is low, which may be related to the γ-ray duty cycle of BL Lac objects. We estimate the γ-ray duty cycle, δγ ≌ 0.11, for BL Lac objects detected by EGRET and Fermi. Using the empirical relation of γ-ray emission with radio emission and the estimated γ-ray duty cycle δγ, we derive the γ-ray luminosity function （LF） of BL Lac objects from their radio LE Our derived γ-ray LF of BL Lac objects can almost reproduce that calculated with the recently released Fermi bright active galactic nuclei （AGN） sample. Comparison of the derived LF of the γ-ray BL Lac objects in this work with that derived by Abdo et al. （2009a） requires the γ-ray duty cycle of BL Lac objects to be almost luminosity-independent. We find that - 45% of the extragalactic diffuse γ-ray background （EGRB） is contributed by BL Lac objects. Combining the estimate of the quasar contribution to the EGRB in the previous work, we find that 77% of the EGRB is contributed by BL Lac objects and radio quasars.

A new three-stage evolution model for millimeter to centimeter wavelength outbursts in BL Lacertae

The multi-frequency light curves of BL Lacertae during 1997.5 - 1999.5 have been modeled by four outbursts, each having a 3-stage evolution in the （Sm, vm） plane with distinct rising-plateau-decaying phases. It is shown that the observed light curves can be well fitted for the eight frequencies from 350 GHz to 4.8 GHz. The main characteristics of the model-fitting are; （1） the outbursts are found to have very flat spectra with an optically thin spectral index α （defined as Sv α u^-α） of about 0.15. This is consistent with the results previously obtained by Valtaoja et al. （1992）; （2） it is found that there is no spectral flattening between the rising-plateau phase and the decay phase. In other words, the optically thin spectral index does not change from the rising-plateau phase to the decay phase. These features are in contrast to the 3-stage shocked-in-jet model proposed by Marscher ＆ Gear （1985） for submm- IR-optical flares, in which a spectral flattening of △α = 0.5 is predicted when a transition occurs from the Compton/synchrotron phase （or rising-plateau phase） to the adiabatic phase （or decay phase） with α≥ 0.5 for the shock being non-radiative. We propose a new model to interpret the fitting results, suggesting that the 3-stage evolution of the mm-cm outbursts in BL Lacertae may be related to the process of shock formation and propagation in a highly collimated jet （for example, a ＇parabolic＇ jet）. In particular, during the rising phase, the thickness of the synchrotron-radiating region created by the shock may rapidly increase with time （relative to the jet width） due to the rapid injection of relativistic electrons and a magnetic field, and this leads to the observed behavior that the turnover flux density Sm rapidly increases while the turnover frequency um decreases. In the decay phase, the emitting plasma enters into a free expansion regime without further injection of relativistic electrons and a magnetic field （for example, when a transition from a collimated regime into a conical regime occurs）. The plateau phase is a short period between the two regimes with no distinct features determined.

Intra-day variability of BL Lacertae from 2016 to 2018

We monitored BL Lacertae in the B,V,R and I bands for 14 nights during the period of 2016-2018.The source showed significant intra-day variability on 12 nights.We performed color-magnitude analysis and found that the source exhibited bluer-when-brighter chromatism.This bluer-when-brighter behavior is at least partly caused by the larger variation amplitude at shorter wavelengths.The variations at different wavelengths are correlated well and show no inter-band time lag.

Is BZB J1450+5201 the most distant gamma-ray BL Lacertae object?

BL Lacertae （BL Lac） objects at high redshifts （z ≥ 2） are rarely detected. Through careful analysis of its Sloan Digital Sky Survey spectrum, BZB J1450＋5201 is confirmed to be a high-z BL Lac object with z ≥ 2.471 by identifying the Lyα 1216 and CIV 1548/1550 absorption lines. This indicates that BZB J1450＋5201 is the most distant BL Lac object discovered to date. Careful analysis of the five-year Fermi-LAT data of 2FGL J1451.0＋5159 shows that its γ-ray emission is robust with a confidence level of 6.2σ at 1-3 GeV and 6.7σ at 3-10 GeV. This analysis with the five-year data overcomes confusion with its bright neighbor, which is a problem when analyzing the two-year data. In addition, 2FGL J 1451.0＋5159 is confirmed to be associated with BZB J1450＋5201 using the five-year data. The analysis of multiwavelength data, from radio to γ-ray energies, indicates that BZB J1450＋5201 is an intermediate synchrotron peaked （ISP） source. Its multiwavelength properties are consistent with distributions of other ISP sources at lower redshifts in the second Fermi-LAT AGN catalog. The pure synchrotron self-Compton （SSC） model seems to be disfavored, but the scattering of weak external emission plus the SSC process can provide a satisfactory description of the broadband emission.

BL Lacertae: Hard Optical Spectrum and GeV γ-ray Emission

The spectral energy distribution (SED) of the γ-ray flare observed inJuly 1997 in BL Lacertae is re-considered. It is pointed out that the optical observa-tions made by Webb et al. showed the associated optical flare has a hard spectrum(the average spectral index αopt～ 0.48, F_v∝v^(-α)), and the ASCA observationsmade by Tanihata et al. showed very steep spectra in the soft X-ray band (0.7-1.5 keV) (α_x～3-4). We find that the flux densities and spectral indices in both theoptical and soft X-ray bands are closely consistent with a ‘canonical' synchrotronspectrum emitted by relativistic electrons of a power-law energy distribution witha high energy cutoff, and thus the peak of the SED of the synchrotron radiation (inrepresentation of vF_v) is located in the EUV -- soft X-ray bands. Therefore, theGeV γ-ray emission observed in the July 1997 outburst may be mainly due to thesynchrotron self-Compton (SSC) process, contrasting with the current explanationsin terms of external radiation Compton (ERC) process, in which the seed photonsare mostly taken to be the UV emission from the clouds of the broad emission lineregion. We argue that the hard optical spectra observed during the γ-ray outburstmay be an important signature for the acceleration of high energy electrons (γ_e～10~4)in the γ-ray emitting region.

A synchrotron self-Compton scenario for the very high energy γ-ray emission of the intermediate BL Lacertae object W Comae

W Comae has significant variability in multi-wavelengthes, from radio to gamma-ray bands. A bright outburst in optical and X-ray bands was observed in 1998, and most recently, a strong TeV flare was detected by VERITAS in 2008. It is the first TeV intermediate-frequency-peaked BL Lacertae source. I find that both the broadband spectral energy distributions （SEDs） which were quasi-simultaneously obtained during the TeV flare and during the optical/X-ray outburst are well fit by using a single-zone synchrotron ＋ synchrotron-self-Compton model. The satisfactory fitting requires a large beaming factor, i.e., δ- 25 and δ- 20 for the TeV flare and the optical/X-ray outburst, respectively, suggesting that both the optical/X-ray outburst and the TeV flare are from a relativistic jet. The size of the emission region of the TeV flare is three times larger than that of the optical/X-ray outburst, and the strength of the magnetic field for the TeV flare is - 14 times smaller than that of the X-ray/optical outburst, likely indicating that the region of the TeV flare is more distant from the core than that of the X-ray/optical outburst. The inverse Compton component of the TeV flare peaks around 1.3 GeV, but it is around 20 MeV for the X-ray/optical outburst, lower than that for the TeV flare by two orders of magnitude. The model predicts that the optical/X-ray outburst might be accompanied by a strong MeV/GeV emission, but the TeV flare may be not associated with the X-ray/optical outburst. The GeV emission is critical for characterizing the SEDs of the optical/X-ray outburst and the TeV flare. The predicted GeV flux is above the sensitivity of Fermi/LAT, and it could be verified with the observations by Fermi/LAT in the near future.

A Quasi-periodic Oscillation of ~4.6 yr in the Radio Light Curves of Blazar PKS 0607-157

We present periodicity search analyses on long-term radio light curves at 4.8,8,and 14.5 GHz of blazar PKS 0607–157 observed by the University of Michigan Radio Astronomical Observatory telescope.The highly variable radio emissions are approximately distributed as a log-normal probability distribution function.The Power Spectral Density for the radio light curves can be well characterized by a power-law model.Using the Weighted Wavelet Z-transform and Lomb-Scargle periodogram methods,significant Quasi-periodic Oscillation(QPO)of∼4.6 yr in the radio light curve has been observed above the 3σconfidence level,which presents an interesting case among blazar QPO phenomena.We explore three plausible physical models to explain the observed QPOs:a supermassive binary black hole system,Lense-Thirring precession of the disk,and helical motion of plasma blobs within the jet.

Comprehensive Study of the Blazars from Fermi-LAT LCR: The Log-Normal Flux Distribution and Linear rms–Flux Relation

Fermi-LAT LCR provides continuous and regularly sampled gamma-ray light curves, spanning about 14 yr, for a large sample of blazars. The log-normal flux distribution and linear rms–flux relation of the light curves for a few Fermi blazars have been examined in previous studies. However, the probability that blazars exhibit the log-normal flux distribution and linear rms–flux relation in their gamma-ray light curves has not been systematically explored.In this study, we comprehensively research the distribution of γ-ray flux and the statistical characteristics on a large sample of 1414 variable blazars from the Fermi-LAT LCR catalog, including 572 FSRQs, 477 BL Lacs, and 365BCUs, and statistically compare their flux distributions with normal and log-normal distributions. The results indicate that the probability of not rejecting log-normal is 42.05% for the large sample, and there is still a 2.05%probability of not rejecting normality, based on the joint of Kolmogorov–Smirnov, Shapiro–Wilk, and Normality tests. We further find that the probability that BL Lacs conform to the log-normal distribution is higher than that of FSRQs. Besides, after removing sources with less than 200 data points from this large sample, a sample of 549blazars, which is still a large sample compared to the previous studies, was obtained. Based on dividing the light curves into segments every 20 points(or 40 points, or one year), we fitted the linear rms–flux relation of these three different sets and found that the Pearson correlation coefficients are all close to 1 for most blazars. This result indicates a strong linear correlation between the rms and the flux of these 549 blazars. The log-normal distribution and linear rms–flux relation indicate that the variability of the γ-ray flux for most blazars is a non-linear and multiplicative process.

Photometric Monitoring of Blazar 3C 66A with the Yunnan University Astronomical Observatory 1 m Telescope

3C 66A is one of our first batches of photometric monitoring objects with the 1 m optical telescope at Yunnan University Astronomical Observatory.In the present work,the observational campaign was performed from 2021 November 1 to 2022 February 27 in the Johnson-Morgan system V and R bands.The average magnitudes in each band were■=15.52±0.18 mag and■=15.07±0.17 mag.The overall variability amplitudes wereΔV=■,Amp=70.27%andΔR=■,Amp=68.56%,respectively.Most of the intraday variabilities(IDVs)occurred in 2021 December and 2022 February.The minimal rise/decay timescale was about 6 minutes(5.82±2.74 minutes and 6.18±2.81 minutes on 2022 February 11,6.99±3.70 minutes and 6.17±2.91 minutes on 2022 February 12).Durations of these rapid variabilities were from 11.99 to 179.67 minutes.The discrete correlation function analyses between V and R bands showed significantly correlated variability.Color index analysis of ID Vs showed that the spectrums do not change with variabilities.

Fermi Blazars in the Zwicky Transient Facility Survey:A Correlation Study

We report on our correlation analysis for theγ-ray blazars detected with the Large Area Telescope(LAT)onboard the Fermi Gamma-ray Space Telescope(Fermi),for which we use the Fermi-LAT data and the optical zr and zg band data from the Zwicky Transient Facility survey in the time period approximately from 2018 July to 2022February.We start from the full sample of the Fermi-LAT blazars and find~38%of them as significant variables in bothγ-ray and optical bands.Among the variables,which are 705 blazars,95 are selected as our targets based on the requirement for theγ-ray brightnesses.Then 15 out of the 95 blazars are found to show theγ-ray-optical correlations in the~3.7 yr time period.Most of the 15 blazars are classified as low synchrotron peaked(LSP)ones with the time lags between theγ-ray and optical bands in a range from-19 to+25 days.The results generally fit the often-considered emission scenario for the LSP blazars.Two cases of showing~-100 days time lags are discussed,which require farther studies for confirmation.Our analysis results suggest that among the detected significant variables,~23%LSPs can showγ-ray and optical correlations,while their averageγ-ray fluxes are above~4×10^(-8)photon cm^(-2)s^(-1).

Optical quasi-periodic oscillation and color behavior of blazar PKS 2155–304

PKS 2155-304 is a well studied BL Lac object in the southern sky. The historical optical data during different periods have been collected and compiled. Light curves spanning 35 yr have been constructed. The R-band light curve has been ana- lyzed by means of three methods： the epoch folding method, the Jurkevich method and the discrete correlation function method. It is derived that there is an evident periodic component of 317 d （i.e. 0.87 yr） superposed on a long-term trend with large- amplitude variation in the light curve. The variability of this source is accompanied by a slight color variation, and the brightness and color index are correlated with each other. On a long-term time scale, PKS 2155-304 exhibits a tendency of bluer-when- brighter, which means the spectrum becomes flatter when the source brightens.

The estimations of four basic parameters for gamma-ray loud blazars

The method used in our previous papers is adopted to estimate four basic parameters （the central black hole mass （M）, the boosting factor （or Doppler factor） （6）, the propagation angle （Φ） and the distance along the axis to the site of the γ-ray production （d）） for 59 γ-ray loud blazars （20 BL Lacertae objects and 39 flat spectrum radio quasars）. The central black hole masses estimated for this sample are in a range of from 107 Me to 109 MG. In the case of black hole mass, there is no clear difference between BL Lacertae objects and flat spectrum radio quasars, which is consistent with the previous results suggesting that the central black hole masses do not play an important role in the evolutionary sequence of blazars.

Statistics of Superluminal Motion in Active Galactic Nuclei

We have collected an up-to-date sample of 123 superluminal sources （84 quasars, 27 BL Lac objects and 12 galaxies） and calculated the apparent velocities （βapp） for 224 components in the sources with the A-CDM model. We checked the relationships between their proper motions, redshifts,βapp and 5 GHz flux densities. Our analysis shows that the radio emission is strongly boosted by the Doppler effect. The superluminal motion and the relativistic beaming boosting effect are, to some extent, the same in active galactic nuclei.

Accretion and jet power in active galactic nuclei

The classical diagrams of radio loudness and jet power as a function of mass and accretion rate of the central spacetime singularity in active galactic nuclei are reanalyzed by including the data of the recently discovered powerful relativistic jets in Narrow-Line Seyfert 1 Galaxies. The results are studied in the light of the known theories of relativistic jets, indicating that, although the Blandford-Znajek mechanism is sufficient to explain the power radiated by BL Lac Objects, it falls to completely account for the power from quasars and Narrow-Line Seyfert 1 Galaxies. This favors the scenario outlined by Cavaliere ＆ D＇Elia of a composite jet, with a magnetospheric core plus a hydromagnetic component emerging when the accretion power increases and the disk becomes radiation-pressure dominated. A comparison with Galactic compact objects is also made, finding some striking similarities, indicating that since neutron stars are low-mass jet systems analogous to black holes, Narrow-Line Seyfert 1 Galaxies are low-mass counterparts of blazars.

Spectral energy distributions for TeV blazars

In this work, we collected a sample of 69 TeV blazars from TeVCat, obtained their multi-wavelength observations, and fitted their spectral energy distributions by using a second degree poly-nomial function. The structure parameters of synchrotron bumps for 68 blazars and those of inverse- Compton bumps for 56 blazars were derived. Then, we conducted statistical analysis on the parame-ters （curvature, peak frequency, peak luminosity, bolometric luminosity and X/γ-ray spectral indexes）. From our analysis and discussions, we can conclude the following： （1） There is a clear positive correla-tion between the synchrotron peak frequency, log （Vp）^s, and the inverse-Compton peak frequency, log （Vp）^IC, and between the synchrotron peak luminosity, log （Vp）^sL（Vp）^s, and the inverse-Compton peak luminosity,log（Vp）^IC L（Vp）^IC. （2） The correlation between the peak frequency and the curvature of synchrotron bump is clearly different from that of the inverse-Compton bump, which further indicates that there are differ-ent emission mechanisms between them. （3） There is a correlation between log （Vp）^IC and γ-ray spectral index, αγ, for the TeV blazars： log （Vp）^IC = -（4.59 ± 0.30）αγ ＋ （32.67 ±0.59）, which is consistent with previous work. （4） An ＂L-shape＂ relation is found between log（Vp）^s Up and αx for both TeV blazars and Fermi blazars. A significant correlation between log（Vp）^s and X-ray photon index （αx） is found for the TeV blazars with high synchrotron peak frequency： log（Vp）^s = -（3.20 ± 0.34）αx ＋ （24.33 ± 0.79）, while the correlation is positive for low synchrotron peaked TeV sources. （5） In the αx -αγ diagram, there is also an ＂L-shape.＂ The anti-correlation is consistent with results available in the literature, and we also find a positive correlation between them. （6） Inverse-Compton dominant sources have luminous bolometric luminosities.

A Possible 456-Day Optical Period of OQ 530

The post-1994 observations of the blazar OQ 530 in optical BVRI bands, and radio observations at 22 GHz and 37 GHz were collected. The date compensated discrete Fourier transform （DCDFr） and CLEANest methods were used to search for possible periodicities. Two possible periods of 456±100 days and 258±29 days were found in the optical bands. The existence of possible correlations between the optical and radio emissions was investigated by means of discrete correlation function （DCF） analysis, and no significant correlation was found. Some possible mechanisms for the periodic variability are discussed.

Optical Photometric Observations of 7-Ray Loud Blazars

We report results of our optical photometric observations of ten gamma-ray loud blazers, namely: 0219+428 (3C66A), PKS 0420-014 (OA 129), S5 0716+714, 0754+100 (OI 090.4), 0827+243 (OJ248), 1652+398 (Mrk 501), 2200+420 (BL Lacertae), 2230+114 (CTA 102), 2251+158 (3C 454.3) and 2344+514. The observations were carried out in September-October, 2000 using the 70 cm optical telescope at Abstumani Observatory, Georgia. We found intra-day variations in 0420-014, S5 0716+714, BL Lacertae and CTA 102. A variation of 0.3 magnitude over a time scale of about 3 hours was observed in the R passband in BL Lacertae on JD 2451827. We did not detect any variation in 3C 66A, Mrk 501, or 3C 454.3 during our observations. Nor did we detect any clear evidence of variation in 1ES 2344+514 during our two weeks' observing run of the TeV gamma-ray source.