Research Article | Open Access

# Long-Term Optical and Spectral Variability of FSRQ 3C454.3

**Academic Editor:**Geza Kovacs

#### Abstract

3C454.3 is a widely studied and monitored flat spectrum radio** quasar**. We use the observations from the Small and Moderate Aperture Research Telescope System** (SMARTS)** to analyze the long-term optical and spectral variabilities of 3C454.3 at B, V, R, J, and K (BVRJK) bands. Based on the relation: , we calculate the multiband spectral indices ()** and analyze** the relations between and , where is the flux density at band (=B, V, R, J, K). Main results are as follows. (1) The largest variations at BVRJK bands are **B = 3.37****0.08 mag, with the timescale ** 1015 days,** at B band**; **V = 3.31****0.07 mag, with the timescale ** 1014 days,** at V band**; **R = 3.62****0.12 mag, with the timescale ** 679 days,** at R band**; **J = 4.08****0.01 mag, with the timescale ** 763 days,** at J band**; and **K = 5.03****0.03 mag, with the timescale ** 2715 days,** at K band**. (2) We analyze the long-term BVRJK lightcurves and obtain the quasiperiodicities: days, days, and years. (3) Multiband lightcurves** show time** delays: = days, = days, = days, and = days. (4) The relations between and show strong correlations, which are typical RWB behaviors; when the source turns to be brighter, the spectral indices turn to be redder.

#### 1. Introduction

Blazars show some extreme properties, such as violently optical variability, core dominance, and superluminal motion [1, 2]. Blazars can be divided into two subclasses: BL Lacs and FSRQs (flat spectrum radio quasars). BL Lacs are characterized by featureless optical spectra or weak emission line [3], and FSRQs are composed of the flat-spectrum radio spectrum and show** typically** broad emission lines [1].

Blazars have relativistic** jets. Their jets** can produce the nonthermal emission, which can dominate the total observed emission.** Many works analyzed the relations between spectral indices (or spectrum) and flux densities (or brightness) (e.g.,** [4ā10]**). Those studies can explain the reasons behind the characteristic of blazar variability and can also help us to constrain the emitting region.** BL Lacs and FSRQs show different correlations. Generally, BL Lacs show that when the sources become brighter, the spectrums become harder and when the sources become fainter, the spectrums become softer. FSRQs show that when the sources become brighter, the spectrums become softer (redder). The redder-when-brighter** (RWB)** behavior of FSRQs can be explained as the contribution from less variable, bluer accretion disk to the variable, redder jet emission.

Kirk and Mastichiadis (1999) [11] put forward that, for blazars, the spectral shape can be influenced by the synchrotron particle acceleration and synchrotron cooling. The spectral shape** can** reflect the intrinsic variability of jet [12ā15]. Fiorucci, Ciprini, and Tosti (2004) [16] pointed out that the optical spectrum of QSOs (subclass of Active Galactic Nuclei) consisted of two components, the first one was variable (, with a flatter slope), which came from synchrotron emission, and the other part was stable (), which might come from the thermal emission. Fan et al. (2014) [17] pointed out that the correlation between brightness and spectral index can be influenced by the brightness intensity.

**For blazars, the optical variable timescales can divide** into three types: (1) the intraday variability (IDV), with the timescales from minutes to hours; (2) short-term variability, with the timescales of days; and (3) long-term variability, which are generally periodic, with the timescales from months to years [18]. The variable timescales can help us to study the physical processes of blazars. The IDV timescales can constrain the size of the emitting region. The long-term periodicity can give us the information about the central black hole mass of the blazars. For example, based on a slim disk, the periodicity can be determined as [18], where is the viscosity coefficient and P is the periodicity, in unit of years.

FSRQ 3C454.3 (z=0.859, [19]) displays strong variabilities over the entire electromagnetic bands (e.g., [20ā25]). Based on the B-band** lightcurve**, Webb et al. (1988) [26] obtained three periodicities, 0.8, 3.0, and 6.4 years. Based on the University of Michigan Radio Astronomy Observatory (UMRAO) database, Ciaramella et al. (2004) [27] obtained a periodicity of years and Fan et al. (2007) [28] obtained two periodicities, years and years. Based on UMRAO and MRO (Metshovi Radio Observatory database), Kudryavtseva and Pyatunina (2006) [29] obtained two periodicities, years and years. Woo et al. (2002) [30] availed the optical luminosity to calculate the central black hole mass** being**āā, where is the solar mass.

This paper is arranged as follows: in Section 2, we introduce observations and data reductions; in Section 3, results are proposed; Section 4 includes discussion and conclusion.

#### 2. Observations and Data Reductions

The data used in this work are collected from** the Small and Moderate Aperture Research Telescope System (SMARTS), which is** located at the Cerro Tololo Inter-American Observatory (CTIO) in Chile. The data reduction and processing** had** been described by Brown et al. (1989) [31] and Bonning et al. (2012) [32].

For 3C454.3, the monitoring time is in the range from June 23, 2008 to November 30, 2015, and the observations cover five optical and near-IR bands, B, V, R, J, and K bands.** The BVRJK lightcurves have been shown in** Figure 1, which are built by(1)753 observations at B band, within days, the largest variations B = 3.370.08 mag (from mag to mag);(2)739 observations at V band, within days, the largest variations V = 3.310.07 mag (from mag to mag);(3)745 observations at R band, within days, the largest variations R = 3.620.12 mag (from mag to mag);(4)744 observations at J band, within days, the largest variations J = 4.080.01 mag (from mag to mag);(5)680 observations at K band, within days, the largest variations K = 5.030.03 mag (from mag to mag).

#### 3. Results

##### 3.1. Periodicity Analysis

It is very important to choose suitable methods to analyze the long-term optical variability. Considering the uneven** lightcurves**, we avail the power spectrum to deal with this question and choose the common part as the quasiperiodicity.

The mostly common tool for periodicity analysis of both evenly and unevenly sampled signals is the Periodogram method, which is an estimator of the signal energy in the frequency domain by Deeming (1975) [33]. Lomb (1976) [34] introduced a modified form of this method, and additionally it was elaborated by Scargle (1982) [35], which can be described in the following. Considering a series with N points, ** is** the frequency and is a variable timescale. Their mean and deviation are given by and . The normalized Lombās , i.e., the power spectrum as a function of the angular frequency is defined asand is defined by the equation:The errors of the period are estimated by the half width at half-maximum (HWHM) of the minimum standard deviation (). We use the Periodogram method to analyze the BVRJK lightcurves and obtain the periodicities days, days, years. The periodic signals are shown in Figure 2. Red noise is a random signal, which has been filtered in order to generate a lot of energy at low frequencies. In order to check the strength of the periodic signal, we compare the periodic signal with the red noise. The red noise [36] with the noise levels 80%, 90%, 95%, and 99% is shown in Figure 2. Based on the results, we can find that at the whole five bands, is higher than 99% noise level, at J and K bands, is higher than 90% noise level, and only K band are higher than 99%.

##### 3.2. The Spectral Indices

We use the following method to obtain the spectral indices. Firstly, we make Galactic Extinction correction, using = , = , = , = , and = , which are collected from NED (http://ned.ipac.caltech.edu/); then we convert the magnitudes () into flux densities (), where is frequency; and at last, based on the relation , the spectral indices () can be calculated by the linear relation: const.

After calculation, there are 752 spectral indices (), which are in the range from to , with the averaged value .

##### 3.3. Relations between Flux Densities and Spectral Indices

**Because the linear relations can clearly show the interdependence between the two parameters, so many works used the linear correlations to analyze the relations between spectral indices** (**) and flux densities** ().** To compare with the others, we use linear correlations to analyze**āā, with the correlation coefficient r and the chance probability p, the slope k, and the intercept b. When the absolute value of r is higher than 0.5 and p is lower than 0.05, the correlations show strong correlation.

**The results are as follows:**āat B band, , with r=0.68, ;āat V band, , with r=0.70, ;āat R band, , with r=0.72, ;āat J band, , with r=0.75, ;āat K band, , with r=0.91, .

The fitting results are shown in Figure 3, with the red lines standing for the fitting lines.

#### 4. Discussion

##### 4.1. Optical Variability

**The optically variable timescale** is an important physical quantity** and** is often used to probe the physics process of blazars. Based on the optical-infrared lightcurves, we obtain the spectral indices () and analyze the long-term variability; see Figure 4 (the upper subpicture). The** spectral variability** shows** quasiperiodic** properties, years, years, and years,** among which**āā** shows** the strongest signal and is consistent with the result calculated from the lightcurves; see Figure 4 (the lower subpicture). The upper results show** that** the lightcurves and** spectral variability** have the same long-term variation** tendencies**.

If the long-term periodicity** ( P)** is caused by a slim disk, the periodicity can be expressed as , where is the viscosity coefficient [18], is periodicity (in unit of years), and is the central black hole mass.

**The mass**calculated from this method is about ,

**which is smaller than**the others, such as calculated by Woo et al. (2002) [30] and calculated by Sbarrato et al. (2012) [37].

Long-term periodicity might come from the influence of binary black hole. OJ287 is considered as a binary black system [38]. Some authors [39, 40] consider PKS 1510-089 as a binary black hole.

3C454.3 might be a binary black hole system.** Based on the variability timescales**, Li et al. (2007) [40] gave a method about how to calculate the masses of the primary black hole (, PBH) and secondary black hole (, SBH). The methods are as follows.(1)The distance (a) between the two black holes can be calculated from the following relation:āwhere is variable** timescale**, is the** quasiperiodicity**, and is the ratio between the jet diameter and radius of secondary black hole.(2)The radio (q) between the mass from the primary black hole and the one from the secondary black hole is(3)For the binary black hole system, the following relation can be deduced

= 0.685 days is collected from [41].** Because the quasiperiodicity 1.24 years (****) are shown both in the lightcurves and in the spectral variability, so we choose the periodicity being 1.24 years to calculate the binary black hole masses. Based on the** expression (3), (4), (5), we can obtainSo, the masses** should be**

If the total mass is [30], we can obtain . When Sillanpaa et al. (1988) [38] analyzed the binary black hole system of OJ287, q = 200 was used. If we take q=200,** the binary black holes masses** are , .

##### 4.2. Time Delay among Different Bands

The** analysis** about** the multiband time delays** can help us to study the emission properties** and** can reflect the terms of electron cooling timescales [42, 43].** For S5 0716+714,** Gupta et al. (2012) [44]** calculated the** time delays** between V and R bands and obtained **** days and -1.4 days. For S5 0716+714,** Wu et al. (2012) [45]** obtained the time delays**, = minutes and = minutes. For 1ES 1959+650, Yuan et al. (2015) [9]** obtained**āā days, days, and days.

We avail the DCF method [46, 47] to explore time delays among different bands. In order to get the delay time more accurately, we use** three methods** to fit the DCF results: local polynomial regression (LOESS) [48], local Kernel regression (LOCFIT) [49], and Gaussian regression. The fitting curves have been added in Figure 5, in which red lines stand for Gaussian regression, blue lines stand for Locfit regression, and green lines stand for Loess regression, based on** which** we can obtain days, days, days, and days. Our results show that the B-band** lightcurves** are ahead of the other bands.

##### 4.3. Relation between Flux Density and Spectral Index

In this work, we analyze the relations between** flux densities and spectral indices**. At the whole five bands (B,V,R,J,K), and show strong correlations. When the source turn to be brighter, the spectrum turn to be redder, which is a typical RWB behavior and consistent with the other FSRQs [32, 50, 51].

We check those distributions and find that there** lie break points; see** Figure 6. We use a broken power law to fit them and obtain the break points: mJy, mJy, mJy, mJy, mJy,** which are noted in** Figure 6. When the source becomes fainter, the relations show correlations and when the source becomes brighter,** their relations turn to be anticorrelations**.

In this paper, we use the observations of 3C454.3 from SMARTs to analyze the lightcurves and calculate the optical spectral indices. Our results show that except K band, the largest variances at B, V, R, and J bands are similar (about 3.5 mag) within about 3 years. The long-term bright variability and spectral variability have the same quasiperiodicities. There lie time delays among different optical and near-IR bands (B, V, R, J, K).

#### Data Availability

The text formatting data used to support the findings of this study have been deposited in the āhttp://www.astro.yale.edu/smarts/glast/home.phpā repository. The data on the upper website are public and may be used by other investigators. If others intend to make use of the data, please inform the SMARTs group by email to glast@elilists.yale.edu and send copies of any resulting publications, including telegrams. Please acknowledge them in your paper by including a citation to the Bonning et al. (2012), mentioned above, and with the following: āThis paper has made use of up-to-date SMARTS optical/near-infrared lightcurves that are available at http://www.astro.yale.edu/smarts/glast/home.php.ā Please also include āSMARTSā as a facility keyword.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

The work is partially supported by the National Natural Science Foundation of China (NSFC 11403006, NSFC 11733001, NSFC U1831119, NSFC U1531245 and NSFC 11503004), Science and Technology Program of Guangzhou (201707010401), Innovative Training Program of Guangzhou University, Innovation Team in Guangdong Province (2014KCXD014), and the financial support for the key subject of Guangzhou City. This paper has made use of up-to-date SMARTS optical/near-infrared lightcurves that are available at http://www.astro.yale.edu/smarts/glast/home.php.

#### References

- C. M. Urry and P. Padovani, āUnified schemes for radio-loud active galactic nuclei,ā
*Publications of the Astronomical Society of the Pacific*, vol. 107, pp. 803ā882, 1995. View at: Publisher Site | Google Scholar - M. H. Ulrich, L. Maraschi, and C. M. Urry, āVariability of active galactic nuclei,ā
*Annual Review of Astronomy and Astrophysics*, vol. 35, pp. 445ā502, 1997. View at: Google Scholar - M. Stickel, J. W. Fried, H. Kuehr, P. Padovani, and C. M. Urry, āThe complete sample of 1 Jansky BL Lacertae objects. I - Summary properties,ā
*The Astrophysical Journal*, vol. 374, p. 431, 1991. View at: Publisher Site | Google Scholar - R. A. Edelson, J. H. Krolik, and G. F. Pike, āBroad-band properties of the CfA Seyfert galaxies. III. Ultraviolet variability,ā
*The Astrophysical Journal*, vol. 359, no. 1, pp. 86ā97, 1990. View at: Publisher Site | Google Scholar - D. Trevese and F. Vagnetti, āQuasar spectral slope variability in the optical band,ā
*The Astrophysical Journal*, vol. 564, pp. 624ā630, 2002. View at: Publisher Site | Google Scholar - F. Vagnetti, D. Trevese, and R. Nesci, āSpectral slope variability of BL Lacertae objects in the optical band,ā
*The Astrophysical Journal*, vol. 590, no. 1 I, pp. 123ā127, 2003. View at: Publisher Site | Google Scholar - B. Z. Dai, X. H. Li, Z. M. Liu et al., āThe long-term multiband optical observations and colour index for the quasar 3C 273,ā
*Monthly Notices of the Royal Astronomical Society*, vol. 392, no. 3, pp. 1181ā1192, 2009. View at: Publisher Site | Google Scholar - K. Poon, J. H. Fan, and J. N. Fu, āThe optical microvariability and spectral changes of the BL Lacertae object S5 0716+714,ā
*The Astrophysical Journal Supplement Series*, vol. 185, pp. 511ā525, 2009. View at: Google Scholar - Y. H. Yuan, J. H. Fan, and H. J. Pan, āOptical photometry of the BL Lac object 1ES 1959+650,ā
*The Astronomical Journal*, vol. 150, no. 3, p. 67, 2015. View at: Publisher Site | Google Scholar - Y. H. Yuan, J. H. Fan, H. J. Pan, Z. Y. Peng, Y. Liu, and H. B. Xiao, āOptical monitoring and spectral analysis of the Seyfert galaxy NGC1275,ā
*Astrophysics and Space Science*, vol. 360, p. 9, 2015. View at: Google Scholar - J. G. Kirk and A. Mastichiadis, āVariability patterns of synchrotron and inverse Compton emission in blazars,ā
*Astroparticle Physics*, vol. 11, no. 1-2, pp. 45ā48, 1999. View at: Publisher Site | Google Scholar - A. Mastichiadis and J. G. Kirk, āVariability in the synchrotron self-Compton model of blazar emission,ā
*Astronomy & Astrophysics*, vol. 320, pp. 19ā25, 1997. View at: Google Scholar - J. G. Kirk, F. M. Rieger, and A. Mastichiadis, āParticle acceleration and synchrotron emission in blazar jets,ā
*Astronomy & Astrophysics*, vol. 333, no. 2, pp. 452ā458, 1998. View at: Google Scholar - W. Cui, āXāRay flaring activity of markarian 421,ā
*The Astrophysical Journal*, vol. 605, no. 2, pp. 662ā669, 2004. View at: Publisher Site | Google Scholar - R. Sato, J. Kataoka, T. Takahashi, G. M. Madejski, S. RĆ¼gamer, and S. J. Wagner,
*The Astrophysical Journal*, vol. 680, no. 1, pp. L9āL12, 2008. View at: Publisher Site - M. Fiorucci, S. Ciprini, and G. Tosti, āThe continuum spectral features of blazars in the optical region,ā
*Astronomy & Astrophysics*, vol. 419, no. 1, pp. 25ā34, 2004. View at: Publisher Site | Google Scholar - J. H. Fan, O. Kurtanidze, Y. Liu, G. M. Richter, R. Chanishvili, and Y. H. Yuan, āOptical monitoring of two brightest nearby quasars, PHL 1811 and 3C 273,ā
*The Astrophysical Journal Supplement Series*, vol. 213, no. 2, 2014. View at: Google Scholar - J. Fan, āOptical variability of blazars,ā
*Chinese Journal of Astronomy and Astrophysics*, vol. 5, no. S1, pp. 213ā223, 2005. View at: Publisher Site | Google Scholar - N. Jackson and W. A. Browne, āOptical properties of quasars. I - Observations. II - Emission-line geometry and radio properties,ā
*Monthly Notices of the Royal Astronomical Society*, vol. 250, pp. 414ā421, 1991. View at: Publisher Site | Google Scholar - S. G. Jorstad, A. P. Marscher, V. M. Larionov et al., āFlaring behavior of the quasar 3C 454.3 across the electromagnetic spectrum,ā
*The Astrophysical Journal*, vol. 715, no. 1, pp. 362ā384, 2010. View at: Publisher Site | Google Scholar - Y. Zhou, D.-H. Yan, and B.-Z. Dai, āThe optical variability properties of flat spectrum radio quasar 3C 454.3,ā
*New Astronomy*, vol. 36, pp. 19ā25, 2015. View at: Publisher Site | Google Scholar - C. M. Raiteri, āWEBT and XMM-Newton observations of 3C 454.3 during the post-outburst phase. Detection of the little and big blue bumps,ā
*Astronomy & Astrophysics*, vol. 473, pp. 819ā827, 2007. View at: Publisher Site | Google Scholar - C. M. Raiteri, āRadio-to-UV monitoring of AO 0235+164 by the WEBT and Swift during the 2006-2007 outburst,ā
*Astronomy & Astrophysics*, vol. 480, pp. 339ā347, 2008. View at: Google Scholar - E. W. Bonning, C. Bailyn, C. M. Urry et al., āCorrelated variability in the blazar 3C 454.3,ā
*The Astrophysical Journal*, vol. 697, no. 2, pp. L81āL85, 2009. View at: Publisher Site | Google Scholar - C. M. Raiteri, M. Villata, M. F. Aller, and M. A. Gurwell, āThe long-lasting activity of 3C 454.3. GASP-WEBT and satellite observations in 2008-2010,ā
*Astronomy & Astrophysics*, vol. 534, p. 87, 2011. View at: Google Scholar - J. R. Webb, A. G. Smith, R. J. Leacock, G. L. Fitzgibbons, P. P. Gombola, and D. W. Shepherd, āOptical observations of 22 violently variable extragalactic sources - 1968-1986,ā
*The Astronomical Journal*, vol. 95, p. 374, 1988. View at: Publisher Site | Google Scholar - A. Ciaramella, āA multifrequency analysis of radio variability of blazars,ā
*Astronomy & Astrophysics*, vol. 419, pp. 485ā500, 2004. View at: Google Scholar - J. H. Fan, Y. Liu, and Y. H. Yuan, āRadio variability properties for radio sources,ā
*Astronomy & Astrophysics*, vol. 462, pp. 547ā552, 2007. View at: Google Scholar - N. A. Kudryavtseva and T. B. Pyatunina, āA search for periodicity in the light curves of selected blazars,ā
*Astronomy Reports*, vol. 50, pp. 1ā11, 2006. View at: Google Scholar - J. H. Woo and C. M. Urry, āActive galactic nucleus black hole masses and bolometric luminosities,ā
*The Astrophysical Journal*, vol. 579, p. 530, 2002. View at: Publisher Site | Google Scholar - L. M. Brown, E. I. Robson, W. K. Gear, and M. G. Smith, āMultifrequency observations of blazars. IV - The variability of the radio to ultraviolet continuum,ā
*The Astrophysical Journal*, vol. 340, p. 150, 1989. View at: Publisher Site | Google Scholar - E. Bonning, C. M. Urry, C. Bailyn et al., āSmarts optical and infrared monitoring of 12 gamma-ray bright blazars,ā
*The Astrophysical Journal*, vol. 756, no. 1, p. 13, 2012. View at: Publisher Site | Google Scholar - T. J. Deeming, āFourier analysis with unequally-spaced data,ā
*Astrophysics and Space Science*, vol. 36, pp. 137ā156, 1975. View at: Google Scholar - N. R. Lomb, āLeast-squares frequency analysis of unequally spaced data,ā
*Astrophysics and Space Science*, vol. 39, no. 2, pp. 447ā462, 1976. View at: Publisher Site | Google Scholar - J. D. Scargle, āStudies in astronomical time series analysis. IIāstatistical aspects of spectral analysis of unevenly spaced data,ā
*The Astrophysical Journal*, vol. 263, pp. 835ā853, 1982. View at: Publisher Site | Google Scholar - M. Schulz and M. Mudelsee, āREDFIT: Estimating red-noise spectra directly from unevenly spaced paleoclimatic time series,ā
*Computers & Geosciences*, vol. 28, no. 3, pp. 421ā426, 2002. View at: Publisher Site | Google Scholar - T. Sbarrato, G. Ghisellini, L. Maraschi, and M. Colpi, āThe relation between broad lines and
*Ī³*-ray luminosities in Fermi blazars,ā*Monthly Notices of the Royal Astronomical Society*, vol. 421, no. 2, pp. 1764ā1778, 2012. View at: Publisher Site | Google Scholar - A. Sillanpaa, S. Haarala, M. J. Valtonen, B. Sundelius, and G. G. Byrd, āOJ 287 - Binary pair of supermassive black holes,ā
*The Astrophysical Journal*, vol. 325, p. 628, 1988. View at: Publisher Site | Google Scholar - M. F. Gu, X. Cao, and D. R. Jiang, āRelation between blue/ultraviolet continuum shape and the ratio of radio-to-optical emission for B3-VLA quasars,ā
*Monthly Notices of the Royal Astronomical Society*, vol. 321, p. 369, 2001. View at: Publisher Site | Google Scholar - J. Li, J. H. Fan, and Y. H. Yuan, āThe parameters of binary black hole systemin PKS 1510-089,ā
*ChPhy*, vol. 16, p. 877, 2007. View at: Google Scholar - X. Fan, S. Li, N. Liao et al., āOptical and Gamma-Ray Variability Behaviors of 3C 454.3 from 2006 to 2011,ā
*The Astrophysical Journal*, vol. 856, no. 1, p. 80, 2018. View at: Publisher Site | Google Scholar - J. M. Bai and M. G. Lee, āRadio/X-ray offsets of large-scale jets caused by synchrotron time lags,ā
*The Astrophysical Journal*, vol. 585, pp. 113ā116, 2003. View at: Publisher Site | Google Scholar - B. K. Zhang, S. Wang, X. Y. Zhao, B. Z. Dai, and M. Zha, āLong-term optical and infrared variability of the BL Lac object PKS 0537 - 441,ā
*Monthly Notices of the Royal Astronomical Society*, vol. 428, pp. 3630ā3635, 2013. View at: Google Scholar - A. C. Gupta, T. P. Krichbaum, and P. J. Wiita, āMultiwavelength intraday variability of the BL Lacertae S5 0716+714,ā
*Monthly Notices of the Royal Astronomical Society*, vol. 425, pp. 1357ā1370, 2012. View at: Google Scholar - J. Wu, M. BĆ¶ttcher, X. Zhou, X. He, J. Ma, and Z. Jiang, āSimultaneous Bā²Vā²Rā² monitoring of bl lacertae object S5 0716+714 and detection of inter-band time delay,ā
*The Astronomical Journal*, vol. 143, no. 5, 2012. View at: Google Scholar - R. A. Edelson and J. H. Krolik, āThe discrete correlation function - A new method for analyzing unevenly sampled variability data,ā
*The Astrophysical Journal*, vol. 333, pp. 646ā659, 1988. View at: Publisher Site | Google Scholar - B. R. Hufnagel and J. N. Bregman, āOptical and radio variability in blazars,ā
*The Astrophysical Journal*, vol. 386, pp. 473ā484, 1992. View at: Publisher Site | Google Scholar - W. S. Cleveland,
*Visualizing Data*, Hobart Press, 1993. - C. Loader,
*Local regression and likelihood*, Statistics and Computing, Springer-Verlag, New York, 1999. View at: MathSciNet - A. RamĆrez, J. A. De Diego, D. Dultzin-Hacyan, and J. N. GonzĆ”lez-PĆ©rez, āOptical variability of PKS 0736+017,ā
*Astronomy & Astrophysics*, vol. 421, no. 1, pp. 83ā89, 2004. View at: Publisher Site | Google Scholar - M. F. Gu, C.-U. Lee, S. Pak et al., āMulti-colour optical monitoring of eight red blazars,ā
*Astronomy & Astrophysics*, vol. 450, pp. 39ā51, 2006. View at: Google Scholar

#### Copyright

Copyright © 2018 Yu-Hai Yuan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.