Long-Term Optical and Spectral Variability of FSRQ 3C454.3
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.370.08 mag, with the timescale 1015 days, at B band; V = 3.310.07 mag, with the timescale 1014 days, at V band; R = 3.620.12 mag, with the timescale 679 days, at R band; J = 4.080.01 mag, with the timescale 763 days, at J band; and K = 5.030.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.
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 , and FSRQs are composed of the flat-spectrum radio spectrum and show typically broad emission lines .
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)  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)  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)  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 . 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 , where is the viscosity coefficient and P is the periodicity, in unit of years.
FSRQ 3C454.3 (z=0.859, ) displays strong variabilities over the entire electromagnetic bands (e.g., [20–25]). Based on the B-band lightcurve, Webb et al. (1988)  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)  obtained a periodicity of years and Fan et al. (2007)  obtained two periodicities, years and years. Based on UMRAO and MRO (Metshovi Radio Observatory database), Kudryavtseva and Pyatunina (2006)  obtained two periodicities, years and years. Woo et al. (2002)  availed the optical luminosity to calculate the central black hole mass being , where is the solar mass.
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)  and Bonning et al. (2012) .
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.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) . Lomb (1976)  introduced a modified form of this method, and additionally it was elaborated by Scargle (1982) , 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  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.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 , 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)  and calculated by Sbarrato et al. (2012) .
3C454.3 might be a binary black hole system. Based on the variability timescales, Li et al. (2007)  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 . 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
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)  calculated the time delays between V and R bands and obtained days and -1.4 days. For S5 0716+714, Wu et al. (2012)  obtained the time delays, = minutes and = minutes. For 1ES 1959+650, Yuan et al. (2015)  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) , local Kernel regression (LOCFIT) , 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).
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 email@example.com 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.
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.
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
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, 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
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
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
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
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
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
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
T. J. Deeming, “Fourier analysis with unequally-spaced data,” Astrophysics and Space Science, vol. 36, pp. 137–156, 1975.View at: 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
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
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
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