Next Article in Journal
The Giant Flares of the Microquasar Cygnus X-3: X-Rays States and Jets
Next Article in Special Issue
Multi-Band Intra-Night Optical Variability of BL Lacertae
Previous Article in Journal
The Jets of Microquasars during Giant Flares and Quiet State
Previous Article in Special Issue
Coordinated Micro-Variability CIRCE Polarimetry and SARA JKT Multi-Frequency Photometry Observations of the Blazar S5 0716+71
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Intra-Night Variability of OJ 287 with Long-Term Multiband Optical Monitoring

1
Department of Physics, Yunnan University, Kunming 650091, China
2
Key Laboratory of Astroparticle Physics of Yunnan Province, Yunnan University, Kunming 650091, China
3
Department of Physics and Technology, Kunming College, Kunming 650214, China
4
Department of Astronomy, Yunnan University, Kunming 650091, China
*
Author to whom correspondence should be addressed.
Galaxies 2017, 5(4), 85; https://doi.org/10.3390/galaxies5040085
Submission received: 28 October 2017 / Revised: 17 November 2017 / Accepted: 20 November 2017 / Published: 22 November 2017
(This article belongs to the Special Issue Microvariability of Blazars)

Abstract

:
We present long-term optical multi-band photometric monitoring of the blazar OJ 287 from 6 March 2010 to 3 April 2016, with high temporal resolution in the V R I -bands. The flux variations and colour-magnitude variations on long and short timescales were investigated to understand the emission mechanisms. In our observation, the major outbursts occurred in January 2016, as predicted by the binary pair of black holes model for OJ 287, with F v a r of 1.3∼2.1%, and variability amplitude (Amp) of 5.8∼9.0%. The intra-night variability (IDV) durations were from 18.5 to 51.3 min, and the minimal variability timescale was about 4.7 min. The colour-magnitude variation showed a weak positive correlation on the long timescale with Pearson’s r = 0 . 450 , while a negative correlation was found on intra-night timescales. We briefly discuss the possible physical mechanisms that are most likely to be responsible for the observed flux and colour-magnitude correlation variability.

1. Introduction

BL Lacertae (BL Lac) objects, which have either very weak or no emission lines [1], and flat-spectrum radio quasars (FSRQs) with strong emission lines [2,3] form a subclass of radio-loud active galactic nuclei (AGNs) known as blazars. Blazars are characterised by non-thermal emission, and strong and rapid flux variability across the entire electromagnetic spectrum, from radio- to gamma-rays. The emission is normally attributed to the relativistic jet oriented at a small angle from the line of sight [4]. Blazar flux variability timescales extending from a few minutes to years and even decades can be broadly divided into three classes: a large variation over hours to days is often known as intra-day or -night variability (IDV) or micro-variability [5,6]; variation on timescales of days to weeks, or even a few months, is considered short-term variability (STV); meanwhile, long-term variability (LTV) can have timescales from several months to years [7,8,9,10].
OJ 287 ( α = 08:54:48.9, δ = +20:06:30.6, J2000) is a blazar at redshift z = 0.306 [11]. Sillanpää et al. [12] pointed out for the first time that there is a double-peak structure in the cyclic optical outbursts of OJ 287 by using the optical V-band observations starting from 1890. The curve exhibited periodic outbursts at intervals of ∼12 years. A binary pair of super-massive black holes system was used to explain this quasi-periodic light curve. Sillanpää et al. [13] (1996a [13], 1996b [14]) reported the next double-peak outburst occurrence in 1994–1995, which occurred almost exactly at the predicted times. Then a double-peaked outburst was also seen during the next recurrence in 2005–2008 [15]. These predicted recurrences are usually interpreted as OJ 287 housing a binary pair of black holes system with a period of ∼12 years.
In the case of OJ 287, outbursts occurring roughly every 12 years are almost certainly produced, and it was predicted that OJ 287 should show a major outburst in 2015–2016 [16]. In our work, along with the confirmation of this predicted outburst during 2015–2016 reported by [10,17], we study the properties of variability and spectral variation in a short timescales.

2. Observations and Data Reduction

Our photometric observations of blazar OJ 287 were performed using the 1.02 m (YAO1.02) optical telescopes at the Yunnan Astronomical Observatory. Since 2009, this telescope has been equipped with an Andor DW436 CCD camera (2048 × 2048 pixels) at the Cassegrain focus (f = 13.3 m). The readout noise and gain are 6.33 electrons and 2.0 electrons/ADU in a 2 μ s readout time model, or 2.29 electrons and 1.4 electrons/ADU in a 16 μ s readout time model. The FOV of the CCD frame is 7.3 × 7.3 arcmin 2 with a pixel scale of 0.21 arcsec/pixel. We used standard Johnson broad-band filters with effective wavelength midpoints of V = 551 nm, R = 658 nm, and I = 806 nm.
We performed optical observations in the V-, R-, and I- bands in the cyclic mode. The exposure times from 50 to 300 s were chosen according to seeing, the filter and the telescope. The time resolutions (the time interval between two adjacent data points in the same band) were less the 15 min and were approximately 10 min in most cases. Therefore, these data in the cyclic mode were considered as quasi-simultaneous measurements, which were explored for analysing inter-band correlation and colour index. During our observation campaign (6 March 2010 to 3 April 2016), we observed for a total of 34 nights, to obtain 2255 CCD frames (shown in Table 1).
All images were reduced with standard image reduction and analysis facility (IRAF) 1 procedures after bias and flat-field corrections. Aperture photometry was performed on the source and comparison stars with apphot. Photometry of the source and comparison stars was performed with the same aperture, which was determined by the full width at half-maximum (FWHM) of the comparison stars and was the same for each observation. We compared the photometric results and found that an aperture radius of about 1.5 FWHM almost always provided the best S/N ratio, and thus we concentrated on this aperture for our analysis. Each night, we performed aperture photometry with different aperture radii.
For each CCD frame, the instrumental magnitudes of OJ 287 and the two comparison stars (listed in Table 2) were extracted. The brightness (magnitude) of OJ 287 was calculated as the average of that derived with respect to comparison stars 4 and 10, and the corresponding standard deviation was treated as the error of OJ 287 within each CCD frame. The deviation of the average of the differential instrumental magnitude of comparison stars 4 and 10, delta (star 4–star 10), was used to verify the stability of the comparison stars, which is also taken as the accuracy of the observations. Figure 1 shows the deviation of comparison stars 4 and 10 for the entire observational campaign. The deviations in bands V-, R- and I-bands are mainly distributed in ±0.025 mag.

3. Results

3.1. Variability

The V, R and I light curves from 6 March 2010 to 3 April 2016 are shown in Figure 2. The average magnitudes in each band are V ¯ = 14 . 728 ± 0 . 352 , R ¯ = 14 . 351 ± 0 . 328 and I ¯ = 13 . 763 ± 0 . 309 . The amplitude variability (maximum minus minimum) in each filter are respectively Δ V = 1 . m 335 , Δ R = 1 . m 237 and Δ I = 1 . m 594 .
Presently, a number of statistical tests, such as the C-test, the F-test, the χ 2 -test and one-way analysis of variance (ANOVA), have been proposed to assess quantitatively whether there are IDVs [6,9,19,20,21,22,23]. Here, we apply the F-test and χ 2 -test to cross-check the intra-day light curves. We consider variability when the light curve satisfied both criteria described below.
The F-test is regarded as a powerful distribution statistic to check for the presence of variability, as introduced by [21]. When comparing two sample variances, the F-statistic values are calculated as:
F 1 = S B L - S t a r A 2 S S t a r A - S t a r B 2 , F 2 = S B L - S t a r B 2 S S t a r A - S t a r B 2
where S B L - S t a r A 2 , S B L - S t a r B 2 and S S t a r A - S t a r B 2 are the variances of the differential instrumental magnitudes of the blazar and comparison star A, the blazar and comparison star B, and comparison stars A and B, respectively. The F-statistic is compared with the F ν B L , ν c ( α ) critical value. The number of degrees of freedom for each sample, ν B L and ν c , are the number of measurements minus 1 ( N - 1 ), and α is the significance level set for the test. In this paper, the F-test was performed at two significance levels: 0 . 99 and 0 . 999 . If either F 1 or F 2 was larger than the critical value, the null hypothesis (no variability) was discarded.
The χ 2 -test was also used:
χ 2 = i = 1 N ( m i - m ¯ σ m i ) 2
where m i are the individual magnitudes, σ m i are their errors, and m ¯ denotes the mean value of magnitude for an observation night. This statistic is also compared against the critical value χ α , ν 2 obtained from the χ 2 probability function, where ν is the degrees of freedom and α is the significance level. χ 2 > χ α , ν 2 provides evidence of variability. OJ 287 was said to be variable (V) if the statistic from both tests was satisfied at the 0.999 level, while it was marked as “non-variable” (NV) if none of these criteria were met at the 0.99 level. For neither, it was marked as “probably variable” (PV). The results of the F- and χ 2 -test and the IDV observations in the V-, R- and I-bands are presented in Table 3. As shown in Table 3, 17 light curves were marked as variable among 100 light curves (18 light curves met the criterion with the F-test and 57 light curves with the χ 2 -test). As an illustration, the significant variability on two nights are presented, 6–8 January in 2016, which were marked as variable in the V-, R- and I-bands (shown in Figure 3). On January 6, 2016, the magnitudes showed a monotonically increasing trend. The magnitude changes of the V-, R- and I-bands were Δ V = 0 . m 074 , Δ R = 0 . m 072 and Δ I = 0 . m 063 in 252 min. On January 8, 2016, the same situation occurred; magnitudes continued to increase by Δ V = 0 . m 087 , Δ R = 0 . m 084 and Δ I = 0 . m 088 in 254 min.
For each IDV, the variability amplitude (Amp) could be calculated by [24]:
A m p = ( A m a x - A m i n ) 2 - 2 σ 2
where A m a x and A m i n are the maximum and minimum magnitude, respectively, of the light curve for the night being considered, and σ is the corresponding standard deviation.
Our observational data were obtained with different exposures and time intervals across 6 years. To avoid uncertainties caused by different observational set-ups, the fractional variability amplitude ( F v a r ) was calculated to show the intrinsic variability amplitude of the source by removing the effects of measurement noise; it is defined as follows [25]:
F v a r = 1 m ¯ S 2 - σ e r r 2 ¯
Here S 2 denotes the total variance of the light curve:
S 2 = 1 N - 1 i = 1 N ( m i - m ¯ ) 2
where σ e r r 2 ¯ is the mean error squared, and m ¯ is the mean value of the magnitude. The error on F v a r is
σ F v a r = 1 2 F v a r 1 N S 2 m ¯ 2
Our results of the IDV analysis are shown in Figure 4 and Table 3. Both the statistical distributions of two variability amplitudes show that the IDV amplitudes of OJ 287 during the whole observational campaign were small.
We extracted the variability for candidate light curves (marked as PV or V states) by an asymmetric flare template [26], which describes flares as exponentially rising and decaying by the function:
f i ( t ) = 2 F i e x p ( t i - t T r , i ) + e x p ( t - t i T d , i ) + F 0
where F 0 is the constant background flux level (could be a constant or linear model), F i is the flare normalisation, t i is the time of maximum flare, T r , i is the flux rising timescale, and T d , i is the flux decaying timescale.
Most of candidate light curves showed a monotonically increasing or decreasing trend, which meant the timescale of the variability was larger than the length of the data. Only two V light curves and four PV light curves showed complete flares. The results of the light curves are shown in Figure 5 and Table 4. On 17 March 2013, the durations of each flare were about 18.5, 21.1 and 23.1 min, respectively. The minimal variability timescale, defined as the minimum of the rising or decaying timescale, was about 5.3 min. On 20 March 2013, the duration was 24.3 min and the minimal variability timescale was about 2.7 min (with ± 4 . 688 min error), which is unreliable. On 5 January 2016, the durations were about 34.4, 45.6 and 51.3 min in the V-, R- and I-bands, respectively. The corresponding minimal variability timescales were about 14.6, 9.2 and 6.5 min. On 6 January 2016, the flux showed a monotonic increase in the V-, R- and I-bands, and in the I-band, there was a significant flare with the flux increasing. The duration was about 42.7 min and the minimal variability timescale was about 4.7 min. This was the identifiable minimum timescale in the whole observational campaign.

3.2. Colour Index

It is common that brightness variations are often associated with changes in spectral shapes. Thus we also investigated the correlations between the colour indices and magnitude variations. Here we consider whether variations in the V - R colour indices of OJ 287 changed with respect to variations in its brightness in the R-band on short- and long-term bases. These colour and magnitude plots of OJ 287 are displayed in Figure 6 and Figure 7. We calculated the best linear fit as shown by the straight lines in Figure 6 and Figure 7 for the colour index against magnitude, and the slope, intercept, linear Pearson correlation coefficients and corresponding null hypothesis probability values are also listed in Table 5. A positive Pearson’s r coefficient between the colour index and apparent magnitude of the blazar implies a positive correlation, which means the source tends to be bluer when it is brighter (BWB) or redder when fainter (RWF), while a negative Pearson’s r coefficient suggests the opposite correlation: redder when brighter (RWB) or bluer when fainter (BWF) behaviour. We found a weak positive correlation between the colour indices and the R-band magnitude on long timescales during these observations of OJ 287 (shown as Figure 6). Then we measured Pearson’s r coefficient between the colour indices and the R-band magnitude on intra-night timescales, as plotted in Figure 7. Ten nights showed a significant negative correlation between the colour indices and R-band magnitude on intra-night timescales, with probability values below 0.001; other nights showed nearly achromatic (ACH) behaviour. On 11–12 March 2010 and 12–13 April 2012, the OJ 287 colour-magnitude correlations tended to show BWF behaviour, while on 14–15 April 2012 they tended to show RWB behaviour.

4. Discussion and Conclusions

Blazars, being dominated by a relativistic jet oriented at a small angle from the line of sight, display variability ranging from a few hundredths of a magnitude to more than several magnitudes. Their variability timescales range from minutes to years. Either factors intrinsic to the jet or due to changes in overall jet power could explain such variability. The shock-in-jet model can be used to explain most of variabilities in blazars [10,27,28,29,30]. Any changes in the physical mechanisms of blazars, such as to the magnetic field, electron density or velocity, can trigger a shock that leads to flares when propagating along the relativistic jet. Turbulence behind the shock-in-jet can be a good way to explain the fast variability in the IDV of blazars [10,31,32,33]. However, when blazars are in a low state, the emission from the accretion disk dominates the jet emission. Instability in the accretion disk, such as at “hot spots”, can also produce IDVs [9,10,34,35,36]. Perturbations on the accretion disks transfer into the relativistic jet and, being Doppler amplified, might also produce fluctuations [9,10,37,38].
Previous colour index studies suggest apparent dichotomy responses: BL Lacs become BWB and FSRQs become RWB [9,10,39,40,41,42,43,44]. The BWB behavior in BL Lacs can be interpreted as electrons being accelerated to preferentially higher energies before radiatively cooling, and the RWB behavior in FSRQs is explained by the addition of redder, non-thermal jet emission to an already bluer, thermal disk component. These results are often preferentially observed during flaring, which causes the selection effect that the two extreme ends of the colour-magnitude correlation are most regularly reported. Ruan et al. [45] point out that spectral variability could be the result of hot spots in the accretion disk emission by using a sample of 604 variable quasars. Isler et al. [46] suggest colour-magnitude variability should be a continuum rather than a dichotomy. They use empirical data from 3C 279 to explain how blazars can smoothly evolve from a jet-quiescent, disk-dominated colour profile to an actively jet-dominated state and associated colour profile and back to a jet-quiescent state. Our results of colour-magnitude correlations on long and short timescales tend to be interpreted in the framework of a relatively bluer accretion disk and redder jet change in intensity. In addition, the jet intensity could change in a short timescales, while the disk intensity changed more slowly. Thus the intra-night colour behavior was alternatively BWF, RWB or ACH, which can be described by the variation in contributions of the non-thermal relativistic jet in OJ 287.
We carried out multi-band optical photometric monitoring of the blazar OJ 287 over 34 observation nights between 6 March 2010 and 3 April 2016. We searched for flux and colour-magnitude variation on IDV timescales. The IDV amplitudes of OJ 287 were small for the observed flux variabilities. The blazar was in a highly active and variable state during January 2016, showing relatively large amplitude flux variations: F v a r of 1.3−2.1%, and Amp of 5.8−9.0%. OJ 287 showed a weak positive correlation on the long timescale, while negative correlations were found on IDV timescales. Our present analysis cannot distinguish between the possible physical mechanisms in detail. Very dense and highly precise simultaneous multi-band observations are necessary. Thus OJ 287 should continue to be monitored whenever possible.

Acknowledgments

This work is partially supported by the National Science Foundation of China (Grants U1531131, 11063003 and 11433004), and the Top Talents Program of Yunnan Province, Yunnan University Research Innovation Fund for Graduate Students and National Undergraduate Innovation Training Program 201610673009. We acknowledge the support of the staff of the 1.02 m optical telescopes at the Yunnan Astronomical Observatory.

Author Contributions

W.Z. and B.-Z.D. are responsible for the writing and the motivation of the paper; X.-F.G., S-B.Y., Z.L., D.-D.W. and Z.-H.K are responsible for observation and reduction; Z.-J.J. and Q.-J.Z. are responsible for data analysis and discussion.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stocke, J.T.; Morris, S.L.; Gioia, I.M.; Maccacaro, T.; Schild, R.; Wolter, A.; Fleming, T.A.; Henry, J.P. The Einstein Observatory Extended Medium-Sensitivity Survey. II—The optical identifications. Astrophys. J. Suppl. Ser. 1991, 76, 813–874. [Google Scholar] [CrossRef]
  2. Blandford, R.D.; Rees, M.J. Extended and compact extragalactic radio sources—Interpretation and theory. Phys. Scr. 1978, 17, 265–274. [Google Scholar] [CrossRef]
  3. Ghisellini, G.; Villata, M.; Raiteri, C.M.; Bosio, S.; de Francesco, G.; Latini, G.; Maesano, M.; Massaro, E.; Montagni, F.; Nesci, R.A. Optical-IUE observations of the gamma-ray loud BL Lacertae object S5 0716+714: Data and interpretation. Astron. Astrophys. 1997, 327, 61–71. [Google Scholar]
  4. Urry, C.M.; Padovani, P. Unified Schemes for Radio-Loud Active Galactic Nuclei. Publ. Astron. Soc. Pac. 1995, 107, 803–845. [Google Scholar] [CrossRef]
  5. Wagner, S.J.; Witzel, A. Intraday Variability in Quasars and BL Lac Objects. Ann. Rev. Astron. Astrophys. 1995, 33, 163–198. [Google Scholar] [CrossRef]
  6. Goyal, A.; Wiita, P.J.; Anupama, G.C.; Sahu, D.K.; Sagar, R.; Joshi, S. Intra-night optical variability of core dominated radio quasars: The role of optical polarization. Astron. Astrophys. 2012, 544, 37–64. [Google Scholar] [CrossRef]
  7. Gupta, A.C.; Banerjee, D.P.K.; Ashok, N.M.; Joshi, U.C. Near infrared intraday variability of Mrk 421. Astron. Astrophys. 2004, 422, 505–508. [Google Scholar] [CrossRef]
  8. Fan, J.H.; Zhang, Y.W.; Qian, B.C.; Tao, J.; Liu, Y.; Hua, T.X. Photometric Monitoring of OJ 287 from 2002 to 2007. Astrophys. J. Suppl. Ser. 2009, 181, 466–472. [Google Scholar] [CrossRef]
  9. Dai, B.Z.; Zeng, W.; Jiang, Z.J.; Fan, Z.H.; Hu, W.; Zhang, P.F.; Yang, Q.Y.; Yan, D.H.; Wang, D.; Zhang, L. Long-term Multi-band Photometric Monitoring of Blazar S5 0716+714. Astrophys. J. Suppl. Ser. 2015, 218, 18–39. [Google Scholar] [CrossRef]
  10. Gupta, A.C.; Agarwal, A.; Mishra, A.; Gaur, H.; Wiita, P.J.; Gu, M.F.; Kurtanidze, O.M.; Damljanovic, G.; Uemura, M.; Semkov, E.; et al. Multiband optical variability of the blazar OJ 287 during its outbursts in 2015–2016. Mon. Not. R. Astron. Soc. 2017, 465, 4423–4433. [Google Scholar] [CrossRef]
  11. Sitko, M.L.; Junkkarinen, V.T. Continuum and line fluxes of OJ287 at minimum light. Publ. Astron. Soc. Pac. 1985, 97, 1158–1162. [Google Scholar] [CrossRef]
  12. Sillanpää, A.; Haarala, S.; Valtonen, M.J.; Sundelius, B.; Byrd, G.G. OJ 287—Binary pair of supermassive black holes. Astrophys. J. 1988, 325, 628–634. [Google Scholar] [CrossRef]
  13. Sillanpää, A.; Takalo, L.O.; Pursimo, T.; Lehto, H.J.; Nilsson, K.; Teerikorpi, P.; Heinaemaeki, P.; Kidger, M.; De Diego, J.A.; Gonzalez-Perez, J.N.; et al. Confirmation of the 12-year optical outburst cycle in blazar OJ 287. Astron. Astrophys. 1996, 305, 17–20. [Google Scholar]
  14. Sillanpää, A.; Takalo, L.O.; Pursimo, T.; Nilsson, K.; Heinamaki, P.; Katajainen, S.; Pietila, H.; Hanski, M.; Rekola, R.; Kidger, M.; et al. Double-peak structure in the cyclic optical outbursts of blazar OJ 287. Astron. Astrophys. 1996, 315, 13–16. [Google Scholar]
  15. Valtonen, M.J.; Nilsson, K.; Villforth, C.; Lehto, H.J.; Takalo, L.O.; Lindfors, E.; Sillanpää, A.; Hentunen, V.P.; Mikkola, S.; Zola, S.; et al. Tidally Induced Outbursts in OJ 287 during 2005–2008. Astrophys. J. 2009, 698, 781–785. [Google Scholar] [CrossRef]
  16. Valtonen, M.J.; Mikkola, S.; Lehto, H.J.; Gopakumar, A.; Hudec, R.; Polednikova, J. Testing the Black Hole No-hair Theorem with OJ287. Astrophys. J. 2011, 742, 22–33. [Google Scholar] [CrossRef]
  17. Valtonen, M.J.; Zola, S.; Ciprini, S.; Gopakumar, A.; Matsumoto, K.; Sadakane, K.; Kidger, M.; Gazeas, K.; Nilsson, K.; Berdyugin, A.; et al. Primary Black Hole Spin in OJ 287 as Determined by the General Relativity Centenary Flare. Astrophys. J. Lett. 2016, 819, 37–42. [Google Scholar] [CrossRef]
  18. Fiorucci, M.; Tosti, G. VRI photometry of stars in the fields of 12 BL Lacertae objects. Astron. Astrophys. Suppl. Ser. 1996, 116, 403–407. [Google Scholar] [CrossRef]
  19. Agarwal, A.; Gupta, A.C. Multiband optical variability studies of BL Lacertae. Mon. Not. R. Astron. Soc. 2015, 450, 541–551. [Google Scholar] [CrossRef]
  20. Dai, B.Z.; Li, X.H.; Liu, Z.M.; Zhang, B.K.; Na, W.W.; Wu, Y.F.; Hao, J.M.; Xiang, Y.; Jiang, Z.J.; Zhang, L. The long-term multiband optical observations and colour index for the quasar 3C 273. Mon. Not. R. Astron. Soc. 2009, 392, 1181–1192. [Google Scholar] [CrossRef]
  21. De Diego, J.A. Testing Tests on Active Galactic Nucleus Microvariability. Astron. J. 2010, 139, 1269–1282. [Google Scholar] [CrossRef]
  22. Goyal, A.; Gopal-Krishna; Wiita, P.J.; Stalin, C.S.; Sagar, R. Improved characterization of intranight optical variability of prominent AGN classes. Mon. Not. R. Astron. Soc. 2013, 435, 1300–1312. [Google Scholar] [CrossRef]
  23. Xiong, D.R.; Bai, J.M.; Zhang, H.J.; Fan, J.H.; Gu, M.F.; Yi, T.F.; Zhang, X. Multicolor Optical Monitoring of the Quasar 3C 273 from 2005 to 2016. Astrophys. J. Suppl. Ser. 2017, 229, 21–38. [Google Scholar] [CrossRef]
  24. Heidt, J.; Wagner, S.J. Statistics of optical intraday variability in a complete sample of radio-selected BL Lacertae objects. Astron. Astrophys. 1996, 305, 42–52. [Google Scholar]
  25. Edelson, R.; Turner, T.J.; Pounds, K.; Vaughan, S.; Markowitz, A.; Marshall, H.; Dobbie, P.; Warwick, R. X-ray Spectral Variability and Rapid Variability of the Soft X-ray Spectrum Seyfert 1 Galaxies Arakelian 564 and Ton S180. Astrophys. J. 2002, 568, 610–626. [Google Scholar] [CrossRef]
  26. Danforth, C.W.; Nalewajko, K.; France, K.; Keeney, B.A. A Fast Flare and Direct Redshift Constraint in Far-ultraviolet Spectra of the Blazar S5 0716+714. Astrophys. J. 2013, 764, 57–63. [Google Scholar] [CrossRef]
  27. Marscher, A.P.; Gear, W.K. Models for high-frequency radio outbursts in extragalactic sources, with application to the early 1983 millimeter-to-infrared flare of 3C 273. Astrophys. J. 1985, 298, 114–127. [Google Scholar] [CrossRef]
  28. Spada, M.; Ghisellini, G.; Lazzati, D.; Celotti, A. Internal shocks in the jets of radio-loud quasars. Mon. Not. R. Astron. Soc. 2001, 325, 1559–1570. [Google Scholar] [CrossRef]
  29. Graff, P.B.; Georganopoulos, M.; Perlman, E.S.; Kazanas, D. A Multizone Model for Simulating the High-Energy Variability of TeV Blazars. Astrophys. J. 2008, 689, 68–78. [Google Scholar] [CrossRef]
  30. Joshi, M.; Böttcher, M. Time-dependent Radiation Transfer in the Internal Shock Model Scenario for Blazar Jets. Astrophys. J. 2011, 727, 21–40. [Google Scholar] [CrossRef]
  31. Marscher, A.P. Turbulent, Extreme Multi-zone Model for Simulating Flux and Polarization Variability in Blazars. Astrophys. J. 2014, 780, 87–96. [Google Scholar] [CrossRef]
  32. Calafut, V.; Wiita, P.J. Modeling the Emission from Turbulent Relativistic Jets in Active Galactic Nuclei. J. Astrophys. Astron. 2015, 36, 255–268. [Google Scholar] [CrossRef]
  33. Pollack, M.; Pauls, D.; Wiita, P.J. Variability in Active Galactic Nuclei from Propagating Turbulent Relativistic Jets. Astrophys. J. 2016, 820, 12–23. [Google Scholar] [CrossRef]
  34. Chakrabarti, S.K.; Wiita, P.J. Spiral shocks in accretion disks as a contributor to variability in active galactic nuclei. Astrophys. J. 1993, 411, 602–609. [Google Scholar] [CrossRef]
  35. Mangalam, A.V.; Wiita, P.J. Accretion disk models for optical and ultraviolet microvariability in active galactic nuclei. Astrophys. J. 1993, 406, 420–429. [Google Scholar] [CrossRef]
  36. Fan, J.H.; Rieger, F.M.; Hua, T.X.; Joshi, U.C.; Li, J.; Wang, Y.X.; Zhou, J.L.; Yuan, Y.H.; Su, J.B.; Zhang, Y.W. A possible disk mechanism for the 23-day QPO in Mkn 501. Astropart. Phys. 2008, 28, 508–515. [Google Scholar] [CrossRef]
  37. Webb, W.; Malkan, M. Rapid Optical Variability in Active Galactic Nuclei and Quasars. Astrophys. J. 2000, 540, 652–677. [Google Scholar] [CrossRef]
  38. Wiita, P.J. Accretion Disks, Jets and Blazar Variability. Astron. Soc. Pac. 2006, 350, 183–190. [Google Scholar]
  39. Massaro, E.; Nesci, R.; Maesano, M.; Montagni, F.; D’Alessio, F. Fast variability of BL Lacertae at 1mum. Mon. Not. R. Astron. Soc. 1998, 299, 47–50. [Google Scholar] [CrossRef]
  40. Ghosh, K.K.; Ramsey, B.D.; Sadun, A.C.; Soundararajaperumal, S. Optical Variability of Blazars. Astrophys. J. Suppl. Ser. 2000, 127, 11–26. [Google Scholar] [CrossRef]
  41. Trévese, D.; Vagnetti, F. Quasar Spectral Slope Variability in the Optical Band. Astrophys. J. 2002, 564, 624–630. [Google Scholar] [CrossRef]
  42. Villata, M.; Raiteri, C.M.; Kurtanidze, O.M.; Nikolashvili, M.G.; Ibrahimov, M.A.; Papadakis, I.E.; Tsinganos, K.; Sadakane, K.; Okada, N.; Takalo, L.O.; et al. The WEBT BL Lacertae Campaign 2000. Astron. Astrophys. 2002, 390, 407–421. [Google Scholar] [CrossRef]
  43. Ramírez, A.; de Diego, J.A.; Dultzin-Hacyan, D.; González-Pérez, J.N. Optical variability of PKS 0736+017. Astron. Astrophys. 2004, 421, 83–89. [Google Scholar] [CrossRef]
  44. Gu, M.F.; Lee, C.-U.; Pak, S.; Yim, H.S.; Fletcher, A.B. Multi-colour optical monitoring of eight red blazars. Astron. Astrophys. 2006, 450, 39–51. [Google Scholar] [CrossRef]
  45. Ruan, J.J.; Anderson, S.F.; Dexter, J.; Agol, E. Evidence for Large Temperature Fluctuations in Quasar Accretion Disks From Spectral Variability. Astrophys. J. 2014, 783, 105–115. [Google Scholar] [CrossRef]
  46. Isler, J.C.; Urry, C.M.; Coppi, P.; Bailyn, C.; Brady, M.; MacPherson, E.; Buxton, M.; Hasan, I. A Consolidated Framework of the Color Variability in Blazars: Long-term Optical/Near-infrared Observations of 3C 279. Astrophys. J. 2017, 844, 107–114. [Google Scholar] [CrossRef]
1
IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.
Figure 1. Distribution of the deviation of the average of the differential instrumental magnitudes of comparison stars 4 and 10 in bands V, R and I.
Figure 1. Distribution of the deviation of the average of the differential instrumental magnitudes of comparison stars 4 and 10 in bands V, R and I.
Galaxies 05 00085 g001
Figure 2. Light curves from 6 March 2010 to 3 April 2016 in V-, R- and I-bands. The error bars of some data points are smaller than the symbol.
Figure 2. Light curves from 6 March 2010 to 3 April 2016 in V-, R- and I-bands. The error bars of some data points are smaller than the symbol.
Galaxies 05 00085 g002
Figure 3. Light curves of two significant micro-variabilities on 6–8 January in 2016, which were marked as variable in V-, R- and I-bands.
Figure 3. Light curves of two significant micro-variabilities on 6–8 January in 2016, which were marked as variable in V-, R- and I-bands.
Galaxies 05 00085 g003
Figure 4. Distribution of the fractional variability amplitude ( F v a r ) and variability amplitude (Amp) in V, R and I in Table 3.
Figure 4. Distribution of the fractional variability amplitude ( F v a r ) and variability amplitude (Amp) in V, R and I in Table 3.
Galaxies 05 00085 g004
Figure 5. Light curves of OJ 287 on 17 March 2013 for R-band; 20 March 2013 for V-band; 5 January 2016 for V-band; 5 January 2016 for R-band; 5 January 2016 I-band; and 6 January 2016 for I-band. The lines are the best fitting by an asymmetric flare template.
Figure 5. Light curves of OJ 287 on 17 March 2013 for R-band; 20 March 2013 for V-band; 5 January 2016 for V-band; 5 January 2016 for R-band; 5 January 2016 I-band; and 6 January 2016 for I-band. The lines are the best fitting by an asymmetric flare template.
Galaxies 05 00085 g005
Figure 6. Optical colour index vs magnitude plots of OJ 287 during our monitoring. Fitted value for slope is 0.041 ± 0.003, intercept is −0.212 ± 0.039 and Pearson’s r coefficient is 0.450 with Pro 6 . 080 × 10 - 46 .
Figure 6. Optical colour index vs magnitude plots of OJ 287 during our monitoring. Fitted value for slope is 0.041 ± 0.003, intercept is −0.212 ± 0.039 and Pearson’s r coefficient is 0.450 with Pro 6 . 080 × 10 - 46 .
Galaxies 05 00085 g006
Figure 7. Optical colour index vs magnitude plots of OJ 287. The colour bars indicate the progression of time in units of JD start . Detail about the colour-magnitude correlation variability can be found in Table 5.
Figure 7. Optical colour index vs magnitude plots of OJ 287. The colour bars indicate the progression of time in units of JD start . Detail about the colour-magnitude correlation variability can be found in Table 5.
Galaxies 05 00085 g007
Table 1. A short summary of the observations in each band performed from 2010 and 2016. In total we observed for 34 nights and obtained 2255 CCD frames.
Table 1. A short summary of the observations in each band performed from 2010 and 2016. In total we observed for 34 nights and obtained 2255 CCD frames.
YearNightsV(N)R(N)I(N)
201010183196197
20124133136136
20137142153156
20145424655
20152192539
20166199201197
Table 2. The magnitudes of the comparison stars. Comparison stars (VRI) from Fiorucci & Tosti [18].
Table 2. The magnitudes of the comparison stars. Comparison stars (VRI) from Fiorucci & Tosti [18].
StarVRI
ID(mag)(mag)(mag)
414.18 (0.04)13.74 (0.04)13.28 (0.04)
1014.60 (0.05)14.34 (0.05)14.03 (0.05)
Table 3. Results of the intra-day/night variability (IDV) observations of OJ 287.
Table 3. Results of the intra-day/night variability (IDV) observations of OJ 287.
Date Band N Magnitude (mag) F-Test χ 2 -TestAmp %Flux mJy F var %State
F 1 , F 2 , F 0.99 , F 0.999 χ 2 , χ 0.99 2 , χ 0.999 2
6 March 2010V3314.468 ± 0.0112.996, 1.858, 2.318, 3.092219.915, 53.486, 62.4874.4806.172 ± 0.0630.908 ± 0.001PV
R3314.136 ± 0.0091.790, 2.214, 2.318, 3.09242.286, 53.486, 62.4873.2056.520 ± 0.0540.000 ± 0.000PV
I3313.543 ± 0.0102.039, 0.685, 2.318, 3.09279.848, 53.486, 62.4873.28410.079 ± 0.0880.444 ± 0.002PV
7 March 2010V2514.485 ± 0.0091.783, 1.285, 2.659, 3.735144.765, 42.980, 51.1793.3726.075 ± 0.0500.653 ± 0.001PV
R2514.152 ± 0.0070.798, 1.501, 2.659, 3.73521.153, 42.980, 51.1792.7486.425 ± 0.0390.000 ± 0.000NV
I2513.559 ± 0.0050.709, 0.365, 2.659, 3.73534.276, 42.980, 51.1792.2739.932 ± 0.0500.000 ± 0.000NV
10 March 2010 R1214.301 ± 0.0140.936, 0.502, 4.462, 7.7613.314, 24.725, 31.2644.0015.600 ± 0.0710.000 ± 0.000NV
I1213.713 ± 0.0091.427, 0.706, 4.462, 7.76146.593, 24.725, 31.2642.9008.625 ± 0.0690.562 ± 0.002PV
11 March 2010V2014.643 ± 0.0051.146, 0.579, 3.027, 4.47423.267, 36.191, 43.8201.5685.254 ± 0.0270.000 ± 0.000NV
R1914.268 ± 0.0207.582, 8.658, 3.128, 4.683309.459, 34.805, 42.3126.0295.773 ± 0.1051.740 ± 0.003V
I2013.712 ± 0.0071.382, 1.689, 3.027, 4.47417.720, 36.191, 43.8202.5978.634 ± 0.0590.000 ± 0.000NV
12 March 2010V2614.636 ± 0.0092.454, 2.456, 2.604, 3.62933.089, 44.314, 52.6203.8795.287 ± 0.0460.376 ± 0.003NV
R2614.263 ± 0.0161.599, 3.275, 2.604, 3.629128.718, 44.314, 52.6205.7275.799 ± 0.0841.275 ± 0.002PV
I2613.704 ± 0.0061.287, 1.618, 2.604, 3.62939.641, 44.314, 52.6201.7098.690 ± 0.0470.265 ± 0.002NV
13 March 2010V1114.679 ± 0.0301.028, 2.299, 4.849, 8.75429.427, 23.209, 29.58812.0605.083 ± 0.1442.381 ± 0.007PV
R1314.308 ± 0.0423.205, 6.351, 4.155, 7.005121.479, 26.217, 32.90915.8975.566 ± 0.2153.605 ± 0.008PV
I1313.744 ± 0.0191.566, 3.960, 4.155, 7.00563.668, 26.217, 32.9097.5978.382 ± 0.1471.523 ± 0.004PV
8 April 2010V1614.602 ± 0.0040.781, 1.767, 3.522, 5.5352.780, 30.578, 37.6971.3955.454 ± 0.0200.000 ± 0.000NV
R1614.261 ± 0.0031.034, 0.769, 3.522, 5.5353.770, 30.578, 37.6971.0155.807 ± 0.0160.000 ± 0.000NV
I1613.717 ± 0.0040.662, 1.544, 3.522, 5.53530.781, 30.578, 37.6971.1998.587 ± 0.0280.197 ± 0.001PV
9 April 2010V1614.664 ± 0.0071.056, 1.612, 3.522, 5.53517.115, 30.578, 37.6972.2165.154 ± 0.0310.000 ± 0.000NV
R1614.317 ± 0.0051.194, 1.316, 3.522, 5.53511.264, 30.578, 37.6971.6625.517 ± 0.0250.000 ± 0.000NV
I1613.772 ± 0.0061.201, 3.118, 3.522, 5.53539.401, 30.578, 37.6972.3448.165 ± 0.0460.416 ± 0.001PV
10 April 2010V2114.680 ± 0.0080.991, 1.034, 2.938, 4.29031.456, 37.566, 45.3152.5045.079 ± 0.0380.000 ± 0.000NV
R2114.333 ± 0.0060.987, 1.140, 2.938, 4.29019.336, 37.566, 45.3152.1625.434 ± 0.0280.000 ± 0.000NV
I2113.786 ± 0.0071.137, 1.657, 2.938, 4.29030.316, 37.566, 45.3152.1858.060 ± 0.0520.000 ± 0.000NV
11 April 2010V1514.644 ± 0.0070.816, 1.661, 3.698, 5.93013.450, 29.141, 36.1231.7765.247 ± 0.0310.000 ± 0.000NV
R1514.297 ± 0.0040.923, 0.620, 3.698, 5.9308.547, 29.141, 36.1231.6515.620 ± 0.0210.000 ± 0.000NV
I1513.749 ± 0.0051.502, 0.635, 3.698, 5.93018.913, 29.141, 36.1231.6538.341 ± 0.0390.113 ± 0.004NV
12 April 2012V3514.301 ± 0.0110.867, 1.752, 2.258, 2.98346.245, 56.061, 65.2475.2087.197 ± 0.0750.000 ± 0.000NV
R3613.936 ± 0.0100.781, 1.054, 2.231, 2.93449.662, 57.342, 66.6193.4407.834 ± 0.0700.000 ± 0.000NV
I3713.358 ± 0.0121.144, 0.759, 2.205, 2.88864.849, 58.619, 67.9855.87611.959 ± 0.1260.627 ± 0.002PV
13 April 2012V2414.322 ± 0.0100.830, 1.150, 2.719, 3.85326.900, 41.638, 49.7284.3407.060 ± 0.0630.000 ± 0.000NV
R2513.964 ± 0.0101.058, 2.296, 2.659, 3.73538.305, 42.980, 51.1794.0257.636 ± 0.0680.234 ± 0.005PV
I2513.389 ± 0.0060.494, 0.577, 2.659, 3.73528.780, 42.980, 51.1792.77511.616 ± 0.0640.000 ± 0.000NV
14 April 2012V3814.268 ± 0.0100.925, 0.927, 2.181, 2.84435.374, 59.893, 69.3464.2207.417 ± 0.0680.000 ± 0.000PV
R3913.906 ± 0.0101.039, 0.903, 2.157, 2.80358.151, 61.162, 70.7034.4638.054 ± 0.0770.000 ± 0.000NV
I3913.333 ± 0.0100.938, 1.674, 2.157, 2.80396.397, 61.162, 70.7033.84612.235 ± 0.1130.658 ± 0.001PV
15 April 2012V3614.157 ± 0.0160.763, 1.056, 2.231, 2.934114.645, 57.342, 66.6196.6468.219 ± 0.1170.650 ± 0.004PV
R3613.795 ± 0.0161.698, 1.443, 2.231, 2.934123.812, 57.342, 66.6195.5638.922 ± 0.1301.025 ± 0.002PV
I3513.224 ± 0.0151.085, 1.146, 2.258, 2.983150.305, 56.061, 65.2475.93213.528 ± 0.1861.048 ± 0.002PV
16 March 2013V914.484 ± 0.0382.330, 0.912, 6.029, 12.04629.022, 20.090, 26.1249.1376.085 ± 0.2132.830 ± 0.010PV
R914.119 ± 0.0100.116, 0.919, 6.029, 12.0463.322, 20.090, 26.1242.8376.620 ± 0.0640.000 ± 0.000NV
I913.546 ± 0.0180.488, 1.091, 6.029, 12.04627.032, 20.090, 26.1245.12010.056 ± 0.1650.877 ± 0.007PV
17 March 2013V2714.509 ± 0.0180.584, 1.063, 2.554, 3.53234.531, 45.642, 54.0527.4485.942 ± 0.0960.000 ± 0.000NV
R2914.149 ± 0.0191.534, 2.097, 2.464, 3.36263.400, 48.278, 56.8927.5996.444 ± 0.1120.882 ± 0.005PV
I2913.569 ± 0.0180.474, 1.714, 2.464, 3.36271.269, 48.278, 56.8927.7639.843 ± 0.1601.176 ± 0.003PV
19 March 2013V1414.771 ± 0.0270.877, 1.384, 3.905, 6.40915.011, 27.688, 34.5289.8504.669 ± 0.1141.090 ± 0.010NV
R1714.395 ± 0.0392.594, 4.396, 3.372, 5.20574.936, 32.000, 39.25215.1795.138 ± 0.1893.168 ± 0.007PV
I1713.837 ± 0.0170.762, 0.926, 3.372, 5.20531.938, 32.000, 39.2526.7107.691 ± 0.1240.801 ± 0.006NV
20 March 2013V2214.820 ± 0.0262.317, 2.289, 2.857, 4.12757.574, 38.932, 46.7979.7604.464 ± 0.1081.803 ± 0.005PV
R2314.454 ± 0.0244.800, 3.911, 2.785, 3.98357.654, 40.289, 48.2689.4234.863 ± 0.1061.564 ± 0.004V
I2413.871 ± 0.0282.890, 3.196, 2.719, 3.853189.788, 41.638, 49.72812.2917.458 ± 0.1912.310 ± 0.004PV
21 March 2013V2415.074 ± 0.0311.165, 0.624, 2.719, 3.85338.513, 41.638, 49.72811.0543.533 ± 0.1011.050 ± 0.011NV
R2614.680 ± 0.0381.070, 0.460, 2.604, 3.62964.444, 44.314, 52.62017.8093.951 ± 0.1400.000 ± 0.000PV
I2614.103 ± 0.0281.161, 0.612, 2.604, 3.62932.312, 44.314, 52.62012.6616.022 ± 0.1480.958 ± 0.009NV
23 March 2013V1815.102 ± 0.0864.072, 4.066, 3.242, 4.924155.018, 33.409, 40.79026.4383.451 ± 0.2717.236 ± 0.014V
R2114.716 ± 0.0774.759, 4.914, 2.938, 4.290197.476, 37.566, 45.31524.1183.831 ± 0.2656.347 ± 0.012V
I2314.153 ± 0.0481.671, 1.774, 2.785, 3.983112.549, 40.289, 48.26819.8155.755 ± 0.2503.583 ± 0.008PV
24 March 2013V2814.981 ± 0.0453.927, 3.201, 2.507, 3.44397.643, 46.963, 55.47622.6813.851 ± 0.1533.383 ± 0.006V
R2814.598 ± 0.0480.224, 1.097, 2.507, 3.443116.079, 46.963, 55.47622.8194.262 ± 0.1852.327 ± 0.011PV
I2814.012 ± 0.0221.513, 2.013, 2.507, 3.443103.345, 46.963, 55.4767.2336.547 ± 0.1361.623 ± 0.004PV
31 March 2014V714.729 ± 0.1994.712, 8.024, 8.466, 20.03081.483, 16.812, 22.45860.8834.936 ± 0.89117.063 ± 0.051PV
R714.430 ± 0.20956.965, 52.442, 8.466, 20.03072.512, 16.812, 22.45862.6915.053 ± 0.82415.795 ± 0.045V
I813.786 ± 0.1081.572, 0.982, 6.993, 15.019297.725, 18.475, 24.32236.0978.102 ± 0.7857.760 ± 0.030PV
1 April 2014V1215.164 ± 0.2420.029, 1.201, 4.462, 7.76134.159, 24.725, 31.26498.8433.347 ± 0.9740.000 ± 0.000PV
R1114.725 ± 0.1440.091, 1.260, 4.849, 8.75487.322, 23.209, 29.58849.3523.824 ± 0.5212.436 ± 0.163PV
I1314.105 ± 0.2310.717, 2.427, 4.155, 7.005166.964, 26.217, 32.909104.2316.133 ± 1.15913.975 ± 0.050PV
2 April 2014V1514.869 ± 0.07411.715, 10.409, 3.698, 5.930324.871, 29.141, 36.12330.3304.277 ± 0.3106.913 ± 0.014V
R1514.458 ± 0.04510.458, 10.827, 3.698, 5.930204.227, 29.141, 36.12316.5054.849 ± 0.2053.999 ± 0.008V
I1513.826 ± 0.0333.458, 2.949, 3.698, 5.930291.375, 29.141, 36.1239.4007.772 ± 0.2382.869 ± 0.006PV
3 April 2014V414.941 ± 0.0330.625, 1.733, 29.457, 141.1085.923, 11.345, 16.2667.1403.994 ± 0.1211.760 ± 0.018NV
R614.481 ± 0.26559.894, 73.206, 10.967, 29.7521001.118, 15.086, 20.51580.6094.884 ± 1.17423.894 ± 0.070V
I813.750 ± 0.1218.107, 3.941, 6.993, 15.019443.834, 18.475, 24.32237.4648.385 ± 1.00611.503 ± 0.031PV
4 April 2014V414.958 ± 0.1246.919, 2.946, 29.457, 141.10850.771, 11.345, 16.26625.3013.954 ± 0.44310.089 ± 0.044PV
R714.478 ± 0.0842.329, 4.181, 8.466, 20.03049.708, 16.812, 22.45823.6294.772 ± 0.3506.598 ± 0.022PV
I1113.701 ± 0.155101.854, 98.563, 4.849, 8.7541163.937, 23.209, 29.58842.1378.815 ± 1.37915.563 ± 0.034V
27 March 2015V1714.650 ± 0.0530.129, 0.562, 3.372, 5.20518.000, 32.000, 39.25221.6365.228 ± 0.2580.000 ± 0.000NV
R1714.270 ± 0.0390.431, 1.543, 3.372, 5.20535.530, 32.000, 39.25213.7755.766 ± 0.2101.923 ± 0.012PV
I2113.699 ± 0.0300.401, 0.751, 2.938, 4.29023.917, 37.566, 45.31512.0738.736 ± 0.2420.000 ± 0.000NV
28 March 2015V214.715 ± 0.0689.5464.928 ± 0.3063.685 ± 0.052
R814.384 ± 0.0340.439, 0.354, 6.993, 15.01910.458, 18.475, 24.3227.1355.187 ± 0.1620.000 ± 0.000NV
I1813.802 ± 0.0530.800, 0.458, 3.242, 4.92495.727, 33.409, 40.79017.6947.954 ± 0.3840.970 ± 0.040PV
5 January 2016V4915.233 ± 0.0093.775, 5.041, 1.977, 2.48939.766, 73.683, 84.0374.1243.052 ± 0.0240.000 ± 0.000PV
R4914.830 ± 0.0093.676, 3.784, 1.977, 2.489104.805, 73.683, 84.0373.3683.441 ± 0.0290.527 ± 0.001V
I4914.201 ± 0.0071.467, 1.474, 1.977, 2.489109.280, 73.683, 84.0373.2645.500 ± 0.0340.413 ± 0.001PV
6 January 2016V5015.260 ± 0.01811.933, 15.218, 1.963, 2.465165.498, 74.919, 85.3516.8722.975 ± 0.0501.346 ± 0.002V
R5014.852 ± 0.02017.948, 18.593, 1.963, 2.465395.483, 74.919, 85.3518.9723.371 ± 0.0611.635 ± 0.002V
I5014.218 ± 0.0188.024, 10.653, 1.963, 2.465771.064, 74.919, 85.3515.7895.417 ± 0.0871.547 ± 0.002V
8 January 2016V5215.167 ± 0.02437.909, 37.202, 1.936, 2.419251.989, 77.386, 87.9687.9983.243 ± 0.0731.956 ± 0.003V
R5114.764 ± 0.02225.587, 28.567, 1.949, 2.441340.977, 76.154, 86.6617.7883.655 ± 0.0761.870 ± 0.002V
I5114.134 ± 0.0239.004, 10.514, 1.949, 2.4411089.284, 76.154, 86.6618.1085.852 ± 0.1262.084 ± 0.002V
1 April 2016V1214.305 ± 0.0071.185, 0.758, 4.462, 7.7617.032, 24.725, 31.2642.4397.171 ± 0.0480.000 ± 0.000NV
R1213.912 ± 0.0071.625, 1.915, 4.462, 7.76111.356, 24.725, 31.2642.1908.010 ± 0.0510.000 ± 0.000NV
I1113.297 ± 0.0073.768, 4.057, 4.849, 8.754231.306, 23.209, 29.5882.14812.651 ± 0.0780.577 ± 0.001PV
2 April 2016V1014.364 ± 0.0241.266, 2.796, 5.351, 10.10793.338, 21.666, 27.8776.4326.793 ± 0.1511.516 ± 0.007PV
R1013.969 ± 0.0211.404, 0.485, 5.351, 10.107192.535, 21.666, 27.8774.7747.604 ± 0.1450.602 ± 0.013PV
I713.353 ± 0.0200.863, 2.103, 8.466, 20.03084.205, 16.812, 22.4585.46512.017 ± 0.2171.561 ± 0.006PV
3 April 2016V2614.616 ± 0.0090.253, 1.588, 2.604, 3.62932.455, 44.314, 52.6203.7595.385 ± 0.0430.000 ± 0.000NV
R2914.226 ± 0.0080.412, 1.525, 2.464, 3.36242.732, 48.278, 56.8923.5565.997 ± 0.0470.246 ± 0.003NV
I2913.609 ± 0.0060.517, 0.920, 2.464, 3.36280.325, 48.278, 56.8922.3709.491 ± 0.0490.239 ± 0.001PV
Table 4. The fitting parameters and their errors of Figure 5.
Table 4. The fitting parameters and their errors of Figure 5.
Date Band ID F i t i T r , i T d , i
+2,456,000minmin
17 March 2013R10.192 ± 0.052369.03710.612 ± 4.8617.916 ± 3.970
20.291 ± 0.048369.10415.788 ± 3.8655.326 ± 1.983
30.342 ± 0.049369.16413.585 ± 3.2289.547 ± 3.588
20 March 2013V10.159 ± 0.061372.136 ± 0.00221.615 ± 12.0642.712 ± 4.668
5 January 2016V10.063 ± 0.0091393.282 ± 0.00714.594 ± 6.45119.763 ± 7.906
R10.052 ± 0.0141393.275 ± 0.0069.160 ± 6.09036.480 ± 16.328
I10.056 ± 0.0151393.274 ± 0.0046.536 ± 4.59644.752 ± 17.941
6 January 2016I10.079 ± 0.0161394.282 ± 0.00337.986 ± 14.5154.672 ± 2.890
Table 5. Color-magnitude dependencies and colour-magnitude correlation coefficients on short timescales of Figure 7.
Table 5. Color-magnitude dependencies and colour-magnitude correlation coefficients on short timescales of Figure 7.
DateData PointsSlopeInterceptPearson’s CoefficientPro
6 March 201049−0.460 ± 0.1886.833 ± 2.661−0.3360.018
7 March 201022−0.145 ± 0.2352.382 ± 3.324−0.1370.545
11 March 201022−0.954 ± 0.06313.989 ± 0.895−0.9590.000
12 March 201035−0.935 ± 0.10813.702 ± 1.539−0.8330.000
8 April 201016−0.311 ± 0.2954.774 ± 4.210−0.2710.310
9 April 201016−0.671 ± 0.3469.950 ± 4.951−0.4600.073
10 April 201040−1.026 ± 0.24515.047 ± 3.506−0.5620.000
11 April 201029−0.783 ± 0.29211.547 ± 4.171−0.4590.012
12 April 201269−0.918 ± 0.14713.157 ± 2.053−0.6060.000
13 April 201224−1.151 ± 0.20916.434 ± 2.917−0.7620.000
14 April 201275−0.623 ± 0.1069.033 ± 1.473−0.5670.000
15 April 201271−0.541 ± 0.1047.827 ± 1.440−0.5290.000
17 March 201327−0.730 ± 0.17910.689 ± 2.528−0.6330.000
20 March 201321−0.419 ± 0.1826.421 ± 2.626−0.4680.033
21 March 201324−0.891 ± 0.17913.471 ± 2.628−0.7280.000
27 March 201514−1.225 ± 0.29317.866 ± 4.182−0.7700.001
5 January 201649−0.296 ± 0.0954.795 ± 1.411−0.4130.003
6 January 201650−0.157 ± 0.0562.742 ± 0.839−0.3730.008
8 January 2016510.011 ± 0.0410.234 ± 0.6000.0400.778
3 April 201625−0.161 ± 0.1172.679 ± 1.659−0.2760.181

Share and Cite

MDPI and ACS Style

Zeng, W.; Zhao, Q.-J.; Jiang, Z.-J.; Kong, Z.-H.; Liu, Z.; Wang, D.-D.; Geng, X.-F.; Yang, S.-B.; Dai, B.-Z. Intra-Night Variability of OJ 287 with Long-Term Multiband Optical Monitoring. Galaxies 2017, 5, 85. https://doi.org/10.3390/galaxies5040085

AMA Style

Zeng W, Zhao Q-J, Jiang Z-J, Kong Z-H, Liu Z, Wang D-D, Geng X-F, Yang S-B, Dai B-Z. Intra-Night Variability of OJ 287 with Long-Term Multiband Optical Monitoring. Galaxies. 2017; 5(4):85. https://doi.org/10.3390/galaxies5040085

Chicago/Turabian Style

Zeng, Wei, Qing-Jiang Zhao, Ze-Jun Jiang, Zhi-Hui Kong, Zhen Liu, Dong-Dong Wang, Xiong-Fei Geng, Shen-Bang Yang, and Ben-Zhong Dai. 2017. "Intra-Night Variability of OJ 287 with Long-Term Multiband Optical Monitoring" Galaxies 5, no. 4: 85. https://doi.org/10.3390/galaxies5040085

APA Style

Zeng, W., Zhao, Q. -J., Jiang, Z. -J., Kong, Z. -H., Liu, Z., Wang, D. -D., Geng, X. -F., Yang, S. -B., & Dai, B. -Z. (2017). Intra-Night Variability of OJ 287 with Long-Term Multiband Optical Monitoring. Galaxies, 5(4), 85. https://doi.org/10.3390/galaxies5040085

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop