• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2023-05-29 10:27:20WeiZengTaoWenZhenLuGongShiChenFanWuHaoYangZhangandBenZhongDai

    Wei Zeng , Tao Wen, Zhen-Lu Gong, Shi Chen, Fan Wu, Hao-Yang Zhang, and Ben-Zhong Dai

    1 School of Physics & Astronomy, Yunnan University, Kunming 650091, China; zengwei@ynu.edu.cn, bzhdai@ynu.edu.cn

    2 Key Laboratory of Astroparticle Physics, Yunnan Province, Kunming 650091, China

    3 Astronomical Observatory, Yunnan University, Kunming 650500, China

    Abstract 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 V = 15.52 ±0.18 mag and R = 15.07 ±0.17 mag.The overall variability amplitudes were ΔV=0 74, Amp=70.27% and ΔR=0 72, 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 IDVs showed that the spectrums do not change with variabilities.

    Key words: (galaxies:) BL Lacertae objects: individual (3C 66A) – galaxies: photometry – galaxies: active

    1.Introduction

    3C 66A (or0219+428;α=02:22:39.6114,δ = +43:02:07.799, J2000) is one of the most luminous blazars at TeV γ-rays,which is a BL Lac type object classified as an intermediate synchrotron peaked blazar(Ackermann et al.2015) at redshift z=0.444 (Miller & McGimsey 1978) by detected Mg II emission line, or z=0.340 by detected host galaxy (Torres-Zafra et al.2016)).It was detected at VHE γrays (very high energy, VHE; E >100 GeV) by the Very Energetic Radiation Imaging Telescope Array System (VERITAS,Acciari et al.2009)and the Major Atmospheric Gamma Imaging Cherenkov Telescopes(MAGIC,Aleksi? et al.2011).

    In optical bands, the 1.56 m telescope of the Shanghai Astronomical Observatory (ShAO) had been monitored from 1996 December 11 to 2009 December 28(Fan et al.2018).The GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) carried out multiwavelength monitoring in 2003–2004 (B?ttcher et al.2005) and 2007–2008(B?ttcher et al.2009).The Steward Observatory of the University of Arizona in support of the Fermi γ-ray observations (Smith et al.2009) contributed optical data from years 2008 to 2018.Using an extremely wide field of view camera,the Zwicky Transient Facility (ZTF) (Bellm et al.2019; Masci et al.2019)scanning the entire Northern sky from 2018 March also carried out the optical monitoring of 3C 66A in ZTF_g,ZTF_r, and ZTF_i bands.

    Combining the above-mentioned sky survey data and additional observational data, 3C 66A intra-day variability(IDV) has been widely addressed by several authors with timescales ranging from 14.5 to 344.0 minutes (Gupta et al.2012; Kaur et al.2017; Fan et al.2018; Pasierb et al.2020).This source also exhibited quasi-periodic variability.Lainela et al.(1999) investigated a period of 65 days, similar to~50–57 days by B?ttcher et al.(2005), while different longterm analyses claimed periods of 156 days, ~2 yr (Fan et al.2018),~2.3 yr(Cheng et al.2022),~2.5 yr(Kaur et al.2017),and ~3 yr (Otero-Santos et al.2020; Cheng et al.2022).

    Blazars are characterized by non-thermal emission, strong and rapid variability, and corresponding timescales extending from a few minutes to years and even decades.So photometric monitoring is a powerful tool to study the structure and radiation mechanism of blazars by measuring their variability amplitude,periodic or aperiodic timescales,and color behavior.In this work, we performed optical monitoring of 3C 66A to study its flux and color variability on short-term timescales.In Section 2,we will describe the observations and data reduction procedures.Section 3 presents the results.Our discussions and conclusions are given in Section 4.

    2.Observations

    Figure 1.Contour map of YNUAO location.The observatory is located on the top of Yunshan Mountain,with an altitude of 2009.5 m.The upper right corner is the layout of the Chenggong Campus of Yunnan University.

    The one-meter ALT-AZ Ritchey Chret?ien Telescope at Yunnan University Astronomical Observatory (hereafter,YNUAO) was officially operational on 2021 December 27.It is operated by Department of Astronomy, Yunnan University.This is a Ritchey–Chretien telescope made by APM-Telescopes GmbH in Germany,equipped with Andor’s iKon-XL 231 CCD camera platform.

    2.1.Site Conditions

    The location of the YNUAO is 102° 51′ 02″ (E),2 4°4 9 ′1 2″(N),with an altitude of 2009.5 m (Figure 1).Figure 2 shows the data recorded from 2017 March 23 to April 13 to evaluate the site seeing conditions of YNUAO.

    2.2.1-meter Telescope and CCD Characteristics

    The one-meter telescope at YNUAO (YNUAO-1m) has a Ritchey ?

    Chretien design with an effective system focal ratio of f/8 by LOMO Optic Set.It has a hyperbolic primary mirror with an effective diameter of 1005 mm(1019.4 mm mechanical diameter), and a hyperbolic secondary mirror with an effective diameter of 369 mm (375 mm mechanical diameter).Two Nasmyth focal stages with image de-rotators are switched through the tunable tertiary mirror.The images have 80%encircled energy (EE80) of 0 5 maximum, over the entire corrected field of view.The telescope has been optimized for the wavelength range of 350–900 nm with a reflectance of at least 85%on each mirror.The baseline for the telescope mirror coatings (for primary, secondary, and tertiary mirrors) is aluminum with protection.The blind pointing accuracy is less than 10″rms (i.e., open-loop).Tracking accuracy and field derotation compensation will not degrade images at the periphery of the field of view by more than 1 0 rms over a minimum period of 1000 s.Figure 3 shows the distribution of displacement deviation of 1000 s of continuous exposure image over time (single 10 s exposure, 100 consecutive frames).The accuracy can reach 0 72 rms for open-loop tracking.

    The iKon-XL 231 CCD camera is used for photometry for YNUAO-1 m.It is a thermoelectrically(TE)cooled,very large area CCD camera platform with the e2v CCD231-84 Astro back-illuminated sensor, offering a very large 61.4×61.4 mm imaging area from a 4096×4108 array format and 15 μm pixel size.The minimum temperature achieved for the TE cooler is set to-100°C with water/liquid cooling.For a standard silicon sensor with mid-band AR coating, peak quantum efficiency could be larger than 95%.The high sensitivity model offers the absolute best CCD performance available, combining exceptionally low read noise of 2.2 electrons with a very large well depth of 350,000 electrons (in system readout rate 100KHz 18 bit).The minimum dark current achievable is 0.0002 electrons pixel-1second-1(@ sensor temperature of -100°C with 16°C water cooling).

    Figure 2.The curve of seeing at Chenggong Campus of Yunnan University from 2017 March 23 to April 12.

    Figure 3.The distribution of displacement deviation of 1000 s continuous exposure image(single 10 s exposure,100 consecutive frames).The left panel is open-loop tracking, and the right panel is closed-loop tracking.

    2.3.Data Reduction

    Our monitoring of 3C 66A with the YNUAO-1 m was from 2021 November 1 to 2022 February 27 in Johnson–Morgan system V and R bands in the alternate mode.The typical exposure times for V and R bands are 120 s and 80 s respectively.A total of 42 observation nights and 2533 CCD frames (including 1280 in the V band and 1253 in the R band) were obtained.All images were still reduced with the standard image reduction and analysis facility (IRAF).4IRAF is the Image Reduction and Analysis Facility, a general purpose software system for the reduction and analysis of astronomical data.IRAF was written by the National Optical Astronomy Observatories (NOAO) in Tucson,Arizona.However, development and maintenance of IRAF is discontinued since 2013 October.The iraf-community is still working on integrating the available patches into the IRAF source code and fixing bugs.Aperture photometry radii of the source and comparison stars are performed with the same value in each frame,which is determined by the FWHM of the source.As in our previous work, we concentrate on an aperture radius of about 1.5 FWHM for our analysis, which always provides the best signal-to-noise ratio(S/N)(Zeng et al.2017,2018,2019).We used the three comparison stars,star“2,”star“3”and star “5”listed in Table 1, to obtain magnitudes for 3C 66A, and then calculated the average and corresponding standard deviation as the magnitude of 3C 66A in each photometric observation.According to the deviation of the average of the differential instrumental magnitude of comparison stars “2”and “5,”Δ(star 2-star 5), we estimated the uncertainty in our observation campaign to be ±0.05 mag.Figure 4 shows the distribution of uncertainty in our observation in each band for the entire observational campaign.

    Figure 4.Distribution of uncertainty in our observation in each band for the entire observational campaign.The left panel is the distribution of the deviation of the differential instrumental magnitudes of comparison stars 2 and 5 in bands V and R.The right panel is the distribution of full width at half-maximum (FWHM) of 3C 66A for each observation.

    Table 1 The Magnitudes of the Comparison Stars

    3.Results

    3.1.Variability

    We used a variability parameter C, generalized by Romero et al.(1999), as

    Figure 5.Light curves from 2021 November 1 to 2022 February 27 in V and R bands.

    Table 2 Results of the Timescale Analysis are Shown in Figure 6

    where σSource-Star2, σSource-Star5, and σStar2-Star5are standard deviations of the differential instrumental magnitudes of 3C 66A and comparison star 2, 3C 66A and comparison star 5, and comparison stars 2 and 5, respectively.

    The χ2test is also used,

    To investigate thepresence ofvariability,theFtest(and χ2test) was comparedwith the Fα,ν(and)critical value,where ν is the degrees of freedom and α is the significance level.In this paper,both the F test and χ2test were performed at two significance levels: 0.99 and 0.999.As described by Romero et al.(1999), if C >2.576, then the source was classified as a variable at significance levels: 0.99.We considered an IDV to have occurred (variable, V), if the statistics from both tests were satisfied at the 0.999 level (or C >2.576).None of the criteria was met at the 0.99 level (or C <2.576): it was marked as non-variable (NV), nor was it marked as a probable variable (PV).The results of the F-test,C-test, χ2-test, and the IDV observations in V and R bands are presented in Table A1.

    Figure 6.Light curves of 3C 66A with flux in units of mJy and time in units of JD-JDstart.The red lines are the best fitting by an asymmetric flare template and the corresponding gray areas are 3σ confidence intervals.

    Table 3 Color Variability Epochs and Fit Parameters of Figures 11, 12, and 13

    During the 42 observation nights, two light curves were marked as “V”in the V band, while 12 were in the R band.Among the 12 nights IDV occurred,2021 December accounted for 6, and 2022 February accounted for 4.Some of the IDVs showed monotonically increasing or decreasing trends during observations.While for others shown flares, we fitted them with an asymmetric flare template (Danforth et al.2013):

    where Tr,imeans rise timescales and Td,imeans decay timescales.The results of the light curves are shown in Figure 6 and Table 2.The reliable minimal rise/decay timescale is about 6 minutes, while the durations of variations can last 12 minutes or more.

    3.2.Cross-Correlation Analysis

    Figure 7.The DCFs between the V and R bands in daily average magnitudes,the time lag was τcent=0.52±0.08 day.The dots represent the cross-correlation for the data,while the color lines show the distribution of random cross-correlations obtained by the Monte Carlo method with 1σ(red),2σ(orange),and 3σ(green).The left panel shows uncertainties with Model-dependent Monte Carlo Method and the right panel with Model-independent.

    In the variability analysis, we found that the light curves in the V and R bands showed consistency of variability, so we carried out the discrete correlation function (DCF) (Edelson &Krolik 1988) to analyze the correlation between the V and R bands with uncorrected light curves.While there is no general agreement on how to estimate the confidence levels (or significant levels) in DCFs, probably the method that has the widest acceptance is Monte Carlo (MC) simulations.For Model-dependent Monte Carlo Method, confidence is estimated by generating multiple randomized versions of time series with the same statistical properties (a simple power-law power spectral density, PSD ∝ν-α, α ≈2.0±0.5) and rerunning the DCF.All simulated time series are generated by a power law noise algorithm (Timmer & Koenig 1995), which randomizes both the amplitude and phase of the Fourier transform coefficients consistent with the statistical properties of the periodogram.Each simulated time series applied the same sampling,interpolated into the same time grid,and added corresponding errors.The distribution of DCFs of these simulated light curves is used to estimate the confidence levels of the data DCF (Max-Moerbeck et al.2014; Jackman et al.2018).

    The determination of a model responding to discrete time series is similarly fraught with uncertainty because the origin of irregular continuum variations is unknown, which leads to the uncertainty of the index α and simulated light curve model.Therefore,methods by which the uncertainties in DCFs can be estimated in a less model-dependent fashion that is simple to implement.

    In our photometric monitoring program, uncertainties in measured fluxes can be a significant source of error.The uncertainties in our observation campaign are normally distributed (as shown in Figure 4 left panel).We thus modify each flux by random Gaussian deviation based on the quoted error for each datum, which is regarded as “flux randomization”(FR).Cross-correlation results can be highly sensitive to the dispersion of data resulting from non-continuously sampled light curves and different time resolutions (exposure times).Thus, random subset selection (RSS) based on a randomly chosen subset of the original data points is implemented.Each MC realization is based on a number of data points that is typically 37% (1/e ≈0.37, Poisson probability) smaller than the real data set.The two elements can be combined in a single MC realization as Monte Carlo “FR/RSS”(Peterson et al.1998, 2004) method.

    We took 103MC simulations for both methods, shown in Figures 7, 8, and 9, the color lines show the 1σ (red), 2σ(orange), and 3σ (green) uncertainties.Then we fitted DCFs with a Gaussian function as: DCF (τ) =Ae-(τ-τcent)2(2σ2).Generally, the centroid of the DCF τcentis recommended as the lag rather than the location of the peak τpeak.During our observation campaign, variation in V and R bands showed a correlation, and time lag was τcent=0.52±0.08 day for daily average magnitudes and τcent=0.12±0.08 day for all observation frames.From this comparison, it can be found that in the optical V and R bands, the 3C 66A had an intra-day delay timescale.For 2021 December, the time lag was τcent=(5.35±1.02)×10-3day (7.70±1.47 minutes).While in 2022 February, it was τcent=(-3.60±0.43)×10-3day(-5.18±0.16 minutes).As time resolutions were ~3.90 minutes between V and R bands, these time lags were acceptable.

    3.3.The Spectral Properties

    Figure 8.The DCF between the V and R bands, the time lag was τcent=0.12±0.08 day.The dots represent the cross-correlation for all data, while the color lines show the distribution of random cross-correlations obtained by the Monte Carlo method with 1σ (red), 2σ (orange), and 3σ (green).

    Figure 9.The DCFs between the V and R bands.The left panel is from 2021 December,the time lag was τcent=(5.35±1.02)×10-3 day,and the right panel is from 2022 February, the time lag was τcent=(-3.60±0.43)×10-3 days.

    DCF analysis indicated that 3C 66A should have consistent spectral behavior during variability.Here we considered variations in the color indices V-R of 3C 66A with respect to variations in its brightness in the R band in 2021 December and 2022 February (as shown in Figure 10).Analysis of these two time periods showed that a weak positive correlation between the color index and brightness, with a slope of-0.40±0.05, r = -0.40 for 2021 December, and a slope of-0.16±0.05, r = -0.19 for 2022 February.The color bar represents the progression of time of color index variability with brightness, in units of Julian Day and from observation start time.Figures 11, 12, and 13 show the color–magnitude correlation for individual nights in consecutive months.Table 3 lists the slope and correlation coefficient for each color variability diagram.The intra-day analysis results showed that,in most cases, there was a negative correlation between the color index and the magnitude.However, the strength of this correlation had nothing to do with whether the variability occurred.Irregular minor variations in timescales of days superimposed on long-term trends led to a weak negative correlation between the color index and magnitude.

    Figure 10.Optical color index vs.magnitude plots of 3C 66A in 2021 December and 2022 February.The color bars indicate the progression of time in units of JD-JDstart.

    Figure 11.Optical color index vs.magnitude plots of 3C 66A in 2021 December.The color bars indicate the progression of time in units of JD-JDstart.

    Figure 12.Optical color index vs.magnitude plots of 3C 66A in January 2022.The color bars indicate the progression of time in units of JD-JDstart.

    4.Discussions and Conclusions

    3C 66A has been observed intensively in the optical band(Lainela et al.1999;B?ttcher et al.2005,2009;Tao et al.2011;Gupta et al.2012; Kaur et al.2017; Fan et al.2018; Otero-Santos et al.2020;Pasierb et al.2020 and references therein).It is also one of our first batches of photometric monitoring objects with the 1 m optical telescope at Yunnan University Astronomical Observatory (YNUAO-1m).Our 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 wereV= 15.52 ±0.18 mag andR= 15.07 ±0.17 mag.Corresponding variability amplitude in each band was respectively ΔV=0 74,Amp=70.27% and ΔR=0 72, Amp=68.56%, while distribution of Amp was around AmpV,median~10.32% and AmpR,median~9.63%(shown in Figure 14).During this period,3C 66A showed active states twice: 2021 December and 2022 February.The flux showed a continuous increasing or decreasing trend during the observation nights, which meant most intraday light curves showed monotonic variations(shown in Figure 15) with timescales greater than the observational duration.Combined with the minimal rise/decay timescale (shown in Figure 6 and Table 2), the durations of variations were ΔT ?11.99-179.67 minutes.These timescales were consistent with the results of ΔT=23.18-344.02 minutes(with Amp=13.2%–92.8%)by Fan et al.(2018)and ΔT=10–32 minutes (with ΔB=0 52–0 46) by Xie et al.(1992).

    Figure 13.Optical color index vs.magnitude plots of 3C 66A in 2022 February.The color bars indicate the progression of time in units of JD-JDstart.

    Figure 14.Distribution of the variability amplitude(Amp)in V and R bands in Table A1.

    Considering that the non-thermal radiation in the optical band is dominated by synchrotron radiation, we can estimate the strength of the magnetic field as Zeng et al.(2018, 2019),Yan et al.(2018).The radiation cooling timescale (in the observer frame) of typical particles emitting synchrotron radiation can be defined as (Chiaberge & Ghisellini 1999):

    where σTis the electron Thompson cross-section, and UB≡B2/8π is the magnetic energy density.The characteristic frequency can be expressed as:

    In the observer frame, our observed the minimum timescale must be longer than the radiative cooling timescale.Therefore,the strength of magnetic field should be satisfied with

    Figure 15.Intraday light curves in the R band with time in units of JD-JDstart.

    Figure 16.Light curves from 2021 July 17 to 2022 March 2.

    The minimal timescales in our observation campaign were 5.82±2.74 and 6.18±2.81 minutes in the V band (as shown in Table 2),which means the magnetic field strength of the jet should be larger than 7.28 G (δD~27.09).Assuming an equipartition parameter eB≡u(píng)B/ue,which refers to the energy density of the electrons in the co-moving frame, the magnetic field can be estimated asBeB=2.9eB27G for δ=15(B?ttcher et al.2005).Modeling of the overall SED with a Synchrotron Self-Compton(SSC)plus External Compton(EC)model with an external near-infrared radiation field, Abdo et al.(2011),Itoh et al.(2013) found comoving magnetic field B=0.21 ~0.35 G.Mohana et al.(2021) also determined the magnetic field as B=0.014 G(δD~27.09/21.95)through broad-band SED fitting of the two distinct states.Compared with the results by fitting the SED with the theoretical models,the magnetic field strength of the jet obtained by giving the minimum IDV timescale was much higher than expected.Therefore, the rapid variability with such a short timescale in our case may be caused by turbulence in the jet rather than by synchrotron radiation cooling.In the shock-in-jet model, a compact emission region (shock) propagates down the relativistic jet and interacts with a highly nonuniform portion in the jet flow.Any changes in the physical mechanisms of blazars, such as the accelerating particles or compressing magnetic fields, can trigger a shock that leads to optical photometric variabilities.Turbulence behind the shock-in-jet can be an excellent way to explain such minutes’timescale rapid variability of blazars (Wagner & Witzel 1995; B?ttcher& Dermer 2010; Wu et al.2012; B?ttcher 2021).

    In our observation campaign, most of the IDVs occurred in 2021 December and 2022 February.The correlation analysis with DCF between V and R bands revealed time lags of τcent=7.70±1.47 minutes in 2021 December, and τcent=-5.18±0.16 minutes in 2022 February.These inconsistent inter-band time lags were exactly as expected by the internal-shock model for blazars (B?ttcher & Dermer 2010).This semi-analytical model described the time-dependent radiative output from internal shocks arising from colliding relativistic shells in a blazar jet.Variations of specific parameters (such as the external radiation energy density, the electron equipartition fraction, the magnetic-field equipartition fraction, the electron acceleration efficiency parameter, the electron injection spectral index,and the relative Lorentz factor between the shells) may lead to changing signs of inter-band time lags,potentially explaining the lack of persistent trends of time lags in most blazars.

    The relationship between the color index and magnitude is used to ascertain whether variabilities are associated with spectral changes.We investigated the color–magnitude correlation for individual nights (shown in Figures 11, 12, and 13).The individual night results showed that, in most cases, there was a negative correlation between the color index and the magnitude.Irregular variations superimposed on long-term trends led to a weak negative correlation between the color index and magnitude.As shown in Figure 10, the analysis of two time periods with obvious IDVs showed a weak negative color–magnitude correlation (with a slope of -0.40±0.05,r = -0.40 for 2021 December, and a slope of -0.16±0.05,r = -0.19 for 2022 February), which suggested redder when brighter (RWB) or bluer when fainter (BWF) behavior.Gaur et al.(2019)found a weak positive correlation, with a slope of-0.021±0.019, r = 0.13 between B and R versus R magnitude.Such a weak correlation coefficient tended to consider that the spectrums do not change with variabilities.This conclusion was also consistent with the results of B?ttcher et al.(2005, 2009), no matter whether in a high optical state with R ~13.4 or a low brightness state with R ≥14.0,3C 66A did not display spectral variability.Such achromatic variation may be produced by geometric effects (lighthouse effect or microlensing effect) in the shock-in-jet model.When shocks moved toward an observer on helical trajectories in the jet, the flux variations could be produced by the sweeping beam,which usually produced achromatic variation (Camenzind & Krockenberger 1992;Wu et al.2007;Man et al.2016).Gravitational microlensing effects were from the external origin, which usually generated the light curves achromatically and symmetrically (Sun et al.2021).

    3C 66A, as a TeV blazar with optical quasi-periodic variation, was also the target of multiwavelength WEBT monitoring projects and Steward Observatory AGNs Optical Observational Program.Katzman Automatic Imaging Telescope(KAIT,Li et al.2003),the LX-200 0.4 m telescope from the St.Petersburg University, and the AZT-8 0.7 m telescope from the Crimean Astrophysical Observatory(SPB and CrAO)provided long-term optical monitoring data(Otero-Santos et al.2020).Now, as one of our first batches of photometric monitoring objects with YNUAO-1 m, it can well complement the vacant observation time of the ZTF monitoring project.Figure 16 shows the light curve from 2021 July 17 to 2022 March 2, transformations between ZTF magnitudes to Johnson–Morgan system R band following Lupton et al.(2005).5https://www.sdss.org/dr16/algorithms/sdssUBVRITransform/Our data filled in gaps in ZTF observations, and continuous observational data are critical for the accuracy of the quasi-periodic analysis and multiband correlation analysis.YUNAO-1 m will also conduct long-term monitoring with VHE sources, and provide high-quality optical photometric data in conjunction with various sky survey monitoring projects.

    Acknowledgments

    We acknowledge the support of the staff of the Lijiang 2.4 m telescope.This work is supported by the fund for the Youth Project of Basic Research Program of Yunnan Province(202001BB050012) and the Joint Foundation of Department of Science and Technology of Yunnan Province and Yunnan University (202201BF070001-020).This research uses data obtained through the one-meter optical telescopes at Yunnan University Astronomical Observatory (YNUAO-1m) which is funded by the “Yunnan University Development Plan for World-Class Astronomy Discipline.”

    Appendix

    Results of the IDV Observations of 3C 66A

    ORCID iDs

    Wei Zeng https://orcid.org/0000-0002-4144-5378

    Ben-Zhong Dai https://orcid.org/0000-0001-7908-4996

    叶爱在线成人免费视频播放| 听说在线观看完整版免费高清| 少妇的丰满在线观看| 老汉色av国产亚洲站长工具| 丁香欧美五月| 欧美av亚洲av综合av国产av| av国产免费在线观看| 无遮挡黄片免费观看| 日本 av在线| 1024香蕉在线观看| 久久久久久国产a免费观看| 亚洲欧美日韩卡通动漫| 国产精华一区二区三区| aaaaa片日本免费| 国产视频内射| 国产成人一区二区三区免费视频网站| 日韩有码中文字幕| 久久久久性生活片| 99久久成人亚洲精品观看| 欧美日本视频| 2021天堂中文幕一二区在线观| 中文字幕人妻丝袜一区二区| 桃红色精品国产亚洲av| 久久热在线av| 亚洲 国产 在线| 在线视频色国产色| 九九热线精品视视频播放| 亚洲五月婷婷丁香| 亚洲av电影不卡..在线观看| 亚洲成人久久爱视频| 白带黄色成豆腐渣| 日本一本二区三区精品| 精品久久久久久成人av| 伦理电影免费视频| 久久久久久久久免费视频了| 麻豆一二三区av精品| 特级一级黄色大片| 美女大奶头视频| 一个人免费在线观看的高清视频| 久久这里只有精品19| 91老司机精品| 1000部很黄的大片| 日韩三级视频一区二区三区| 国产精品美女特级片免费视频播放器 | 成在线人永久免费视频| 精品无人区乱码1区二区| 午夜精品一区二区三区免费看| 精品福利观看| 丰满的人妻完整版| 18禁裸乳无遮挡免费网站照片| 欧美日韩亚洲国产一区二区在线观看| 757午夜福利合集在线观看| 黄色 视频免费看| 18禁国产床啪视频网站| 青草久久国产| 国产在线精品亚洲第一网站| 不卡一级毛片| 亚洲熟妇中文字幕五十中出| 亚洲精华国产精华精| 婷婷丁香在线五月| 精品福利观看| 宅男免费午夜| 久久午夜亚洲精品久久| 脱女人内裤的视频| 大型黄色视频在线免费观看| 听说在线观看完整版免费高清| 狂野欧美白嫩少妇大欣赏| 日韩大尺度精品在线看网址| 日韩人妻高清精品专区| 亚洲狠狠婷婷综合久久图片| 日本撒尿小便嘘嘘汇集6| 久久久精品大字幕| 动漫黄色视频在线观看| 五月伊人婷婷丁香| 999久久久精品免费观看国产| 精品一区二区三区视频在线 | 国产精品免费一区二区三区在线| 国产野战对白在线观看| 熟女电影av网| 亚洲人成网站在线播放欧美日韩| 亚洲 欧美 日韩 在线 免费| 精品电影一区二区在线| 日韩 欧美 亚洲 中文字幕| 88av欧美| 久久久久久久精品吃奶| 99久久成人亚洲精品观看| 国产精品99久久99久久久不卡| av欧美777| 亚洲欧美日韩卡通动漫| 国产高清videossex| 国产精品一区二区三区四区免费观看 | 国内精品美女久久久久久| 美女黄网站色视频| 午夜免费激情av| 亚洲欧洲精品一区二区精品久久久| 亚洲精品一卡2卡三卡4卡5卡| 最新中文字幕久久久久 | 夜夜躁狠狠躁天天躁| 国产午夜福利久久久久久| 美女高潮喷水抽搐中文字幕| 亚洲欧美日韩高清专用| 国产 一区 欧美 日韩| 99久久精品一区二区三区| 国产乱人视频| 国产黄色小视频在线观看| 桃色一区二区三区在线观看| 啦啦啦观看免费观看视频高清| 午夜激情欧美在线| 色哟哟哟哟哟哟| 国产精品综合久久久久久久免费| 悠悠久久av| 搡老岳熟女国产| 亚洲18禁久久av| 欧美激情在线99| 久99久视频精品免费| 欧美极品一区二区三区四区| 日韩欧美在线乱码| 99re在线观看精品视频| 夜夜爽天天搞| 国产aⅴ精品一区二区三区波| 一二三四在线观看免费中文在| 国产精品久久久人人做人人爽| 精品久久久久久久末码| 99热精品在线国产| 午夜影院日韩av| 精品免费久久久久久久清纯| 亚洲男人的天堂狠狠| 网址你懂的国产日韩在线| 97人妻精品一区二区三区麻豆| 免费观看人在逋| 免费看日本二区| 日本熟妇午夜| 精品乱码久久久久久99久播| 在线免费观看不下载黄p国产 | 精品乱码久久久久久99久播| 成人av在线播放网站| 国产成人精品久久二区二区免费| 精品人妻1区二区| 亚洲一区二区三区不卡视频| 男女下面进入的视频免费午夜| 国产精品爽爽va在线观看网站| 男人和女人高潮做爰伦理| 久久这里只有精品19| 夜夜看夜夜爽夜夜摸| 国产精品亚洲一级av第二区| 色噜噜av男人的天堂激情| 国产精品久久久久久人妻精品电影| 97超级碰碰碰精品色视频在线观看| 国内精品久久久久精免费| 国产黄片美女视频| 亚洲国产精品999在线| 国产精品自产拍在线观看55亚洲| 国产淫片久久久久久久久 | 午夜精品久久久久久毛片777| 特大巨黑吊av在线直播| 日日摸夜夜添夜夜添小说| 日本成人三级电影网站| 神马国产精品三级电影在线观看| 亚洲国产精品成人综合色| 色尼玛亚洲综合影院| 国产精品亚洲一级av第二区| 琪琪午夜伦伦电影理论片6080| 国产69精品久久久久777片 | 麻豆国产av国片精品| 99热这里只有是精品50| 色在线成人网| 一本一本综合久久| 欧美乱码精品一区二区三区| 他把我摸到了高潮在线观看| 97超级碰碰碰精品色视频在线观看| 国内精品一区二区在线观看| 18禁国产床啪视频网站| 一本一本综合久久| 老司机深夜福利视频在线观看| 日本 av在线| 久久九九热精品免费| 国产激情欧美一区二区| 久久亚洲精品不卡| 狂野欧美白嫩少妇大欣赏| 搞女人的毛片| 免费av毛片视频| 欧美日韩一级在线毛片| 无限看片的www在线观看| 九九久久精品国产亚洲av麻豆 | 99国产精品一区二区三区| 国产精品久久久人人做人人爽| 国产欧美日韩一区二区精品| 婷婷亚洲欧美| 琪琪午夜伦伦电影理论片6080| www国产在线视频色| 小说图片视频综合网站| 精品国产乱码久久久久久男人| 欧美极品一区二区三区四区| 亚洲中文字幕一区二区三区有码在线看 | 性欧美人与动物交配| 婷婷精品国产亚洲av| 美女黄网站色视频| 午夜免费激情av| 熟女人妻精品中文字幕| 欧美在线黄色| 男女之事视频高清在线观看| 午夜福利成人在线免费观看| 国产精品99久久99久久久不卡| 日韩成人在线观看一区二区三区| 久久亚洲精品不卡| 午夜福利免费观看在线| 男插女下体视频免费在线播放| 国产成人aa在线观看| 天堂影院成人在线观看| 色老头精品视频在线观看| 村上凉子中文字幕在线| 国产精品爽爽va在线观看网站| 99久久精品国产亚洲精品| 成人精品一区二区免费| 不卡一级毛片| 精品久久久久久久久久久久久| 精品一区二区三区av网在线观看| 757午夜福利合集在线观看| 亚洲乱码一区二区免费版| 国产欧美日韩精品亚洲av| 久久久久九九精品影院| 午夜成年电影在线免费观看| 欧美激情久久久久久爽电影| 亚洲自拍偷在线| av在线天堂中文字幕| 桃红色精品国产亚洲av| 亚洲一区二区三区不卡视频| 亚洲熟妇熟女久久| 成年免费大片在线观看| 9191精品国产免费久久| 国模一区二区三区四区视频 | 精品99又大又爽又粗少妇毛片 | 成年女人看的毛片在线观看| 久久久久久久久久黄片| 精品久久久久久久末码| 久久精品亚洲精品国产色婷小说| 久久精品影院6| 亚洲第一欧美日韩一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 中亚洲国语对白在线视频| 小蜜桃在线观看免费完整版高清| 国产91精品成人一区二区三区| 国产美女午夜福利| 日韩中文字幕欧美一区二区| 在线观看美女被高潮喷水网站 | 高清毛片免费观看视频网站| 99国产精品一区二区三区| 男女床上黄色一级片免费看| 亚洲国产精品久久男人天堂| 村上凉子中文字幕在线| 国产成人aa在线观看| 精品福利观看| 少妇裸体淫交视频免费看高清| 欧美大码av| 日韩欧美精品v在线| 国产精华一区二区三区| 精品福利观看| 亚洲av电影在线进入| 伦理电影免费视频| 最好的美女福利视频网| 欧美日韩一级在线毛片| 99视频精品全部免费 在线 | 成年女人看的毛片在线观看| 亚洲国产中文字幕在线视频| 欧美+亚洲+日韩+国产| 最近最新中文字幕大全免费视频| 免费一级毛片在线播放高清视频| 男女做爰动态图高潮gif福利片| 久久性视频一级片| 怎么达到女性高潮| 白带黄色成豆腐渣| 在线播放国产精品三级| 特级一级黄色大片| 亚洲乱码一区二区免费版| 俄罗斯特黄特色一大片| 久久久久性生活片| 国产免费男女视频| 国产真人三级小视频在线观看| 色综合站精品国产| 香蕉av资源在线| 免费人成视频x8x8入口观看| 国内久久婷婷六月综合欲色啪| 国产亚洲av高清不卡| 18禁黄网站禁片免费观看直播| 啪啪无遮挡十八禁网站| 久久国产精品人妻蜜桃| 国产爱豆传媒在线观看| 国产精品自产拍在线观看55亚洲| e午夜精品久久久久久久| 亚洲无线观看免费| 日本 av在线| 久久久国产成人免费| 99久久精品一区二区三区| 一区二区三区激情视频| 国产精品一及| 亚洲五月婷婷丁香| 两个人看的免费小视频| 欧美日韩综合久久久久久 | 亚洲专区中文字幕在线| 亚洲中文av在线| 亚洲人成伊人成综合网2020| 一级毛片精品| 国内少妇人妻偷人精品xxx网站 | 精品一区二区三区av网在线观看| 99视频精品全部免费 在线 | 亚洲avbb在线观看| 手机成人av网站| 国产一区二区在线观看日韩 | 女人高潮潮喷娇喘18禁视频| 亚洲国产欧美人成| 国产麻豆成人av免费视频| 国产1区2区3区精品| 精品国产超薄肉色丝袜足j| 亚洲国产欧洲综合997久久,| 美女cb高潮喷水在线观看 | 男人舔女人下体高潮全视频| 日本免费a在线| 国产精品久久电影中文字幕| 久久久久国产精品人妻aⅴ院| 色在线成人网| 51午夜福利影视在线观看| 国产麻豆成人av免费视频| 午夜日韩欧美国产| 嫩草影视91久久| 久久久水蜜桃国产精品网| 亚洲色图 男人天堂 中文字幕| 国产aⅴ精品一区二区三区波| 女同久久另类99精品国产91| 国产一区二区在线观看日韩 | 国产高清激情床上av| 成在线人永久免费视频| 很黄的视频免费| 成年免费大片在线观看| 国产欧美日韩精品一区二区| 我的老师免费观看完整版| 国产欧美日韩一区二区三| 欧美黑人巨大hd| 亚洲国产中文字幕在线视频| 免费看光身美女| 国产精品av视频在线免费观看| 免费搜索国产男女视频| 校园春色视频在线观看| 亚洲最大成人中文| 国内揄拍国产精品人妻在线| 岛国在线免费视频观看| 精品一区二区三区四区五区乱码| 天堂影院成人在线观看| 亚洲自偷自拍图片 自拍| 18禁黄网站禁片免费观看直播| 国产一区二区在线av高清观看| 色播亚洲综合网| 国产精品久久电影中文字幕| 999精品在线视频| 91麻豆精品激情在线观看国产| 精品国产乱码久久久久久男人| 国产视频内射| 岛国在线观看网站| 国产野战对白在线观看| 欧美性猛交╳xxx乱大交人| 精品欧美国产一区二区三| 首页视频小说图片口味搜索| 精品国内亚洲2022精品成人| 亚洲av成人不卡在线观看播放网| 亚洲熟妇中文字幕五十中出| 最近在线观看免费完整版| 亚洲欧洲精品一区二区精品久久久| 日本熟妇午夜| 亚洲精品一区av在线观看| 亚洲五月天丁香| 国产1区2区3区精品| 巨乳人妻的诱惑在线观看| 免费观看的影片在线观看| 男女那种视频在线观看| 亚洲欧美精品综合一区二区三区| 日韩欧美国产一区二区入口| 五月玫瑰六月丁香| 十八禁人妻一区二区| 欧美中文日本在线观看视频| 天堂网av新在线| 亚洲人成网站在线播放欧美日韩| 亚洲国产看品久久| 日韩欧美免费精品| 国产美女午夜福利| 成人亚洲精品av一区二区| 黑人操中国人逼视频| 日韩免费av在线播放| 91九色精品人成在线观看| 少妇丰满av| 狂野欧美白嫩少妇大欣赏| 黑人欧美特级aaaaaa片| 伊人久久大香线蕉亚洲五| 欧美性猛交╳xxx乱大交人| 国产视频内射| 日日摸夜夜添夜夜添小说| 国产精品综合久久久久久久免费| 男女床上黄色一级片免费看| 又大又爽又粗| 国产精品久久电影中文字幕| 无限看片的www在线观看| 少妇人妻一区二区三区视频| 一本一本综合久久| 久久久久免费精品人妻一区二区| 美女大奶头视频| 色综合站精品国产| 好男人在线观看高清免费视频| 国内久久婷婷六月综合欲色啪| 国产综合懂色| 中文亚洲av片在线观看爽| 亚洲美女视频黄频| 久久久精品欧美日韩精品| 国产日本99.免费观看| 色综合婷婷激情| 女人高潮潮喷娇喘18禁视频| 亚洲中文av在线| 成人鲁丝片一二三区免费| 亚洲中文日韩欧美视频| 亚洲电影在线观看av| 99精品欧美一区二区三区四区| 十八禁人妻一区二区| www国产在线视频色| 99国产精品99久久久久| 1000部很黄的大片| 黄片小视频在线播放| 99国产综合亚洲精品| 97超级碰碰碰精品色视频在线观看| 99久久久亚洲精品蜜臀av| 亚洲国产欧美一区二区综合| 中文资源天堂在线| 久久99热这里只有精品18| 一本一本综合久久| 国产精品99久久久久久久久| 热99在线观看视频| 99热6这里只有精品| 成人国产综合亚洲| 90打野战视频偷拍视频| 国产久久久一区二区三区| av女优亚洲男人天堂 | 色在线成人网| 国产精华一区二区三区| 欧美精品啪啪一区二区三区| 久久精品国产清高在天天线| 很黄的视频免费| 黄色成人免费大全| or卡值多少钱| 真实男女啪啪啪动态图| www日本黄色视频网| 久99久视频精品免费| 亚洲色图av天堂| 90打野战视频偷拍视频| 九色国产91popny在线| 欧美在线黄色| 熟女人妻精品中文字幕| 国内精品一区二区在线观看| 久久精品综合一区二区三区| 日本在线视频免费播放| 亚洲精品一区av在线观看| 午夜福利在线观看吧| 美女扒开内裤让男人捅视频| 不卡av一区二区三区| 欧美+亚洲+日韩+国产| 国产午夜精品久久久久久| 欧美成狂野欧美在线观看| 天堂影院成人在线观看| 一二三四在线观看免费中文在| 色综合亚洲欧美另类图片| 在线观看午夜福利视频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲自偷自拍图片 自拍| 国产三级在线视频| 嫩草影院精品99| 亚洲一区高清亚洲精品| 搡老熟女国产l中国老女人| 日韩成人在线观看一区二区三区| 俄罗斯特黄特色一大片| 精品一区二区三区四区五区乱码| 欧美高清成人免费视频www| 亚洲精品在线观看二区| 成年女人永久免费观看视频| 很黄的视频免费| 亚洲国产日韩欧美精品在线观看 | 男女床上黄色一级片免费看| 一本精品99久久精品77| 久久精品夜夜夜夜夜久久蜜豆| 久久久久性生活片| 黄色 视频免费看| 午夜免费激情av| 欧美乱色亚洲激情| 亚洲熟妇中文字幕五十中出| 亚洲精品乱码久久久v下载方式 | 美女午夜性视频免费| 国产精品美女特级片免费视频播放器 | 免费看a级黄色片| 在线看三级毛片| 欧美又色又爽又黄视频| 别揉我奶头~嗯~啊~动态视频| av在线蜜桃| 精品欧美国产一区二区三| 午夜a级毛片| 熟妇人妻久久中文字幕3abv| 香蕉av资源在线| 精品久久久久久成人av| 日韩免费av在线播放| 国产高潮美女av| 久久性视频一级片| 亚洲精品一卡2卡三卡4卡5卡| а√天堂www在线а√下载| 一区二区三区国产精品乱码| 老汉色av国产亚洲站长工具| 亚洲欧洲精品一区二区精品久久久| 亚洲成人中文字幕在线播放| 欧美黄色淫秽网站| 香蕉av资源在线| 人妻久久中文字幕网| 男人舔奶头视频| 狂野欧美白嫩少妇大欣赏| 成年女人看的毛片在线观看| 久久天躁狠狠躁夜夜2o2o| 国产精品亚洲美女久久久| 久久中文字幕人妻熟女| 人人妻人人看人人澡| 色综合站精品国产| 国产精品野战在线观看| www.999成人在线观看| 在线播放国产精品三级| 巨乳人妻的诱惑在线观看| 女人被狂操c到高潮| 亚洲真实伦在线观看| 在线观看免费午夜福利视频| 国产极品精品免费视频能看的| 男人的好看免费观看在线视频| 日本成人三级电影网站| 精品免费久久久久久久清纯| 国产午夜精品论理片| 长腿黑丝高跟| 亚洲在线观看片| 美女高潮喷水抽搐中文字幕| 两个人看的免费小视频| 久久精品国产清高在天天线| 欧美大码av| 一进一出好大好爽视频| 国产精品 欧美亚洲| 精品国产三级普通话版| 国产激情久久老熟女| 国产精品香港三级国产av潘金莲| 成熟少妇高潮喷水视频| 91av网一区二区| 中出人妻视频一区二区| www日本黄色视频网| 国产精品香港三级国产av潘金莲| 51午夜福利影视在线观看| 免费在线观看视频国产中文字幕亚洲| 午夜福利免费观看在线| 日韩欧美在线乱码| 黄色视频,在线免费观看| 中文资源天堂在线| 搡老妇女老女人老熟妇| 中文资源天堂在线| 亚洲av第一区精品v没综合| 国产精品乱码一区二三区的特点| 男女之事视频高清在线观看| 99国产精品一区二区三区| 国产高清有码在线观看视频| 又黄又爽又免费观看的视频| 亚洲aⅴ乱码一区二区在线播放| 免费无遮挡裸体视频| 日本 欧美在线| 在线观看免费视频日本深夜| 午夜免费激情av| 999精品在线视频| 一级作爱视频免费观看| 91字幕亚洲| 国产成人av教育| 国产精品久久久久久亚洲av鲁大| 国产欧美日韩一区二区精品| 欧美一区二区精品小视频在线| 午夜精品在线福利| 999精品在线视频| 熟女人妻精品中文字幕| 成人av一区二区三区在线看| 国产欧美日韩一区二区精品| 国产精品综合久久久久久久免费| 黄色片一级片一级黄色片| 国内少妇人妻偷人精品xxx网站 | 成人永久免费在线观看视频| or卡值多少钱| 99国产精品一区二区三区| 黄频高清免费视频| 国语自产精品视频在线第100页| 国产成人影院久久av| 日本黄大片高清| 成年人黄色毛片网站| 国产一区二区三区视频了| 麻豆国产av国片精品| 亚洲精品在线观看二区| 成年女人毛片免费观看观看9| 热99在线观看视频| 男女那种视频在线观看| 亚洲国产欧美人成| 国产伦在线观看视频一区| 免费人成视频x8x8入口观看| 久久久久国产精品人妻aⅴ院| 99久久无色码亚洲精品果冻| 成年版毛片免费区| 久久天躁狠狠躁夜夜2o2o| 老汉色∧v一级毛片| 非洲黑人性xxxx精品又粗又长| 人人妻人人澡欧美一区二区| 一卡2卡三卡四卡精品乱码亚洲| 一进一出好大好爽视频| 高潮久久久久久久久久久不卡| 久久精品综合一区二区三区| 国产成+人综合+亚洲专区| 91久久精品国产一区二区成人 | 制服人妻中文乱码| 日韩三级视频一区二区三区|