Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Superfluid Phonons in Neutron Star Core
Previous Article in Journal
Quantum String Cosmology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Intra-Day Variability Observations of Two Dozens of Blazars at 4.8 GHz

1
Xinjiang Astronomical Observatory, Chinese Academy of Sciences, 150 Science 1-Street, Urumqi 830011, China
2
Key Laboratory of Radio Astronomy, Chinese Academy of Sciences, Nanjing 210008, China
3
Physics and Astronomy Department, Qiannan Normal University for Nationalities, Duyun 558000, China
4
School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 100049, China
5
Max-Plank-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
*
Author to whom correspondence should be addressed.
Universe 2021, 7(1), 15; https://doi.org/10.3390/universe7010015
Submission received: 1 December 2020 / Revised: 8 January 2021 / Accepted: 10 January 2021 / Published: 14 January 2021
(This article belongs to the Section Galaxies and Clusters)

Abstract

:
Two dozens of radio loud active galactic nuclei (AGNs) have been observed with Urumqi 25 m radio telescope in order to search for intra-day variability (IDV). The target sources are blazars (namely flat spectrum radio quasars and BL Lac objects) which are mostly selected from the observing list of RadioAstron AGN monitoring campaigns. The observations were carried out at 4.8 GHz in two sessions of 8–12 February 2014 and 7–9 March respectively. We report the data reduction and the first results of observations. The results show that the majority of the blazars exhibit IDV in 99.9% confidence level, some of them show quite strong IDV. We find the strong IDV of blazar 1357 + 769 for the first time. The IDV at centimeter-wavelength is believed to be predominately caused by the scintillation of blazar emission through the local interstellar medium in a few hundreds parsecs away from Sun. No significant correlation between the IDV strength and either redshift or Galactic latitude is found in our sample. The IDV timescale along with source structure and brightness temperature analysis will be presented in a forthcoming paper.

1. Introduction

The intra-day variability (IDV) of active galactic nuclei (AGNs) at centimeter wavelengths was first discovered in 1980s [1,2]. Since then, more observations have been carried out, and it turns out that the IDV occurs only in flat spectrum radio sources (FSRSs) which are core-dominated AGNs, with a detection rate of 25–50% in the FSRSs [3,4] and with a higher rate of ∼60% in bright Fermi blazars [5].
For a long time, it is controversial if the origin of the radio IDV at centimeter wavelengths, is AGN-intrinsic variations or is due to scattering-induced variations of AGN signals by ionized interstellar medium in our Galaxy [6,7]. If the IDV is AGN-intrinsic, the timescale of a day or less will imply very high brightness temperatures of the jet-emitting region, several magnitudes (typically 10 3 5 times) greater than the inverse-Compton limit of T IC 10 11.5 K [8] or the equipartition temperature of T eq = 5 × 10 10 = 10 10.7 K [9]. On the other hand, the brightness temperatures ( T b ) measured from the RadioAstron (RA) space VLBI (SVLBI) observations, are 10–100 times larger than the IC limit, in a significant fraction of blazars [10,11,12,13]. The T b can be boosted by the Doppler beaming effect when the jet of blazar points closely to our line of sight, i.e., T b o b s = δ T b i n t , where T b o b s is the observed brightness temperature and T b i n t is brightness temperature in source frame. However, the Doppler factor δ is <60 for blazars observed from ground-based VLBI observations [14,15,16] and mostly δ 10 [11], which is not high enough to account for the T b excess in the space VLBI observations, implying that the IC limit might be violated in some of the blazars. Nevertheless, the apparent ‘brightness temperature’ implied from the IDV, as mentioned above, is much higher than the T b of blazars measured with the RA space baselines.
Today, it becomes clear that the IDV at centimeter wavelengths is not AGN-intrinsic, but caused by the interstellar scintillation (ISS) effect in our Galaxy [17]. The evidence is from the annual modulations of IDV timescales discovered in a couple of IDV sources that had been monitored for years, which can be well fitted by the Earth orbiting around the Sun and annually across the scintillation pattern (e.g., [18,19,20], see a summary on this in [21]).
Furthermore, the flux fluctuations caused by the ISS-induced IDV, may have contributed partly to the RA SVLBI observed T b excess (over the IC limit) of blazars [22]. So it is important to find and study the IDV from the blazars in the RA SVLBI observations [23]. In order to search for IDV in the RA sources in parallel with the RA SVLBI sessions, we carried out the IDV observations of the RA sources (plus a few blazars from other samples) with the Urumqi 25 m radio telescope in Xinjiang Astronomical Observatory. In this paper, we report the data reduction and the first results of the observations performed in search of IDV.

2. Observations and Data Reduction

The IDV observations are carried out in 8 February 2014 05:12–12 February 2014 03:21 (94 h in total) and 7 March 2014 03:52–9 March 2014 17:53 (62 h in total) with the Urumqi 25 m radio telescope. Most of the blazars are selected from the 2014 February and March sessions of RadioAstron space VLBI program [23], with the declination > + 10 , plus several additional blazars selected from previous IDV observations. The observations were done with the 25 m telescope at the central frequency of 4.8 GHz with the bandwidth of 600 MHz, the lowest system temperature of 24 K, and the antenna sensitivity of ∼0.12 K/Jy.
The observations were made with the cross-scans mode consisting of 4 + 4 sub-scans in azimuth and elevation over the source position. After initial calibration of raw data, the intensity profile of each sub-scan is fitted with a baseline and a Gaussian function after subtracting the baseline, and then the data are averaged in 4 azimuth sub-scans and 4 elevation sub-scans respectively. The averaged azimuth amplitude is corrected with the pointing error in elevation direction, and vice verse. The antenna gain curve is applied for the data averaged from the azimuth and elevation scans, the gain curve (including the air-mass effect) is measured from the secondary calibrators in the same observation. A systematic time-dependent effect which is derived from the secondary calibrators is applied. Finally, the amplitude is converted to the absolute flux density with the average scale of the primary calibrators, 3C48, 3C286, and NGC7027 [24,25] which are included in the observation.
It takes about 5 min for each source with 4 + 4 observing sub-scans, and about 2 h for nearly 25 sources to complete one cycle of observation, if the sources are up (> 10 elevation) in the sky, and then we repeat the cycle of observation. So the sampling rate is about 0.5 per hour for each source.
To analyze the variability, we used the statistical quantities, such as the modulation index m (a measure of the relative strength of variations), the modified modulation index or fractional variability [26] m (a parameter that similar to m but take into account the error of flux density, when the part within the square root is less than 0, we set m = 0 , which imply that there is no variability), the variability amplitude Y which is assumed to be unbiased to the systematic noise, and the reduced chi-squared χ r 2 , to describe the variability as shown below, also see [27],
m [ % ] = 100 σ s < S > ,
m [ % ] = 100 σ s 2 < Δ S i 2 > < S > 2 ,
Y [ % ] = 3 ( m ) 2 ( m 0 ) 2 ,
χ r 2 = 1 N 1 ( S i < S > Δ S i ) 2 ,
where σ s is the standard deviation of flux densities from the mean flux density ( < S > ) in a light curve, m 0 is the mean modulation index of all calibrators in the same observation, N is the number of measurements in a light curve, S i and Δ S i denote individual flux density and its error.

3. Results and Discussion

3.1. Variability Characteristics

The results are listed in Table 1 and Table 2. In the February IDV session, 16 targets from the RA observing session are observed, plus 4 additional targets from other samples and 5 calibrators. In the March IDV session, 14 targets from the RA observing session are observed, plus 3 additional targets from other samples and 5 calibrators. There are more than 20 data points sampled for the majority of sources in about three-day observation of each session, some sources have fewer data points due to their low declinations or bad data being deleted by the calibration software (e.g., when Full Width of Half Power of source profile > 120% of antenna beam width, or antenna pointing error > 1/6 of beam width, or error of Gaussian fit > 10% of peak of the Gaussian, or amplitude of a sub-scan profile > 110% of average amplitude of 4 + 4 sub-scans; these can be caused by too large source, too weak source or interference, etc.). After careful calibrations as introduced above, we obtained flux densities and light curves for each source. Here, we only show a few sources for example, three sources with strong variability—1357 + 769, J1128 + 5925 and WISE J092915.43 + 501336.0 in Figure 1, two counterexamples—B2 1239 + 376 and 1633 + 382 with no sign of variability in Figure 2.
The same as our previous work [20], in order to ascertain whether or not an IDV is present, we applied a χ 2 test to each data set. We adopt the criterium that a data set with ‘a probability of ≤0.1% to be constant’ is considered to be variable (i.e., the ‘null hypothesis’ probability <1.0 × 10 3 in Table 1 or Table 2), which is equivalent to the confidence level of 99.9% to define the IDV. According to this definition, there are 14 out of 20 targets show IDV in session 1, all 17 targets show IDV in session 2.
In Table 1 and Table 2, we find there are a couple of prominent IDV sources from the two observing sessions, e.g., those with the modulation index m > 3 m 0 , the m 0 is the average modulation index of the calibrators, which is 0.64 for session 1, and 0.45 for session 2.
The modified modulation index m represents a variability strength, but not a real variability amplitude. The parameter Y is more like a variability amplitude, which ranges from a few per cent to ∼22% for the prominent IDV sources in Table 1 and Table 2.
We note that there are 11 common targets appeared in both session 1 and session 2, all of them have shown IDV in the two sessions except the quasar 1156 + 295 which only shows IDV in session 2, we give a short discussion for this in Section 3.3.
Some of the IDV sources show relatively fast variability, e.g., the source J1128 + 5925, WISE J092915.43 + 501336.0 exhibit a short timescale of a few hours and some of them e.g., the source 1357 + 769 shows a relatively longer timescale. A full analysis on the IDV timescales along with the source structure and brightness temperature derived from the RA space VLBI measurements will be presented in a forthcoming paper.

3.2. Comment on Individual Sources

We comment briefly on some prominent IDV sources which showed a modulation index m > 3 m 0 in one or two observing sessions, in the following. We note that all the target sources are blazars with core-dominated radio structures in milliarcsecond scales from the VLBI images (e.g., [28,29]).
PKS 0528 + 134: the source shows IDV in both sessions, see Table 1 and Table 2, the IDV strength is slightly increased from session 1 to session 2, while the mean flux density of 1.710 ± 0.034 Jy in session 1 decreases to 1.525 ± 0.034 Jy, at 4.8 GHz.
0716 + 714: this BL Lac object is a well known IDV source, although its IDV strength is moderate with relatively longer IDV timescales from previous observations and its IDV timescales can be fitted with the annual modulated ISS model [20]. It shows IDV in both sessions, and its IDV strength in session 1 is less than that in session 2. The total flux density of 1.255 ± 0.016 Jy in session 1 decreases to 1.193 ± 0.017 Jy in session 2.
0748 + 126: this blazar shows IDV in both sessions, its IDV strength in session 1 is higher than that in session 2. Its mean flux density is similar from session 1 to session 2.
OJ287: this is a BL Lac object, which is thought to be a supermassive binary black hole system from long-term optical periodic flux variations [30], but even the RA space VLBI still cannot resolve the binary. It is observed in session 2 only and shows an IDV.
1044 + 719: this blazar shows IDV in both sessions. The IDV strength increases from session 1 to session 2. Its mean flux density also increases from 2.581 ± 0.035 Jy in session 1 to 2.786 ± 0.049 Jy in session 2.
1357 + 769: this blazar shows strong IDVs in both sessions (upper panel of Figure 1), the m decreases from 5.28 in session 1 to 3.02 in session 2, while the mean flux density decreases from 0.461 ± 0.025 Jy in session 1 to 0.411 ± 0.013 Jy in session 2.
Table 1. The result of the IDV observation in 8–12 February 2014 with Urumqi 25 m telescope at 4.8 GHz. The columns are: (1) name; (2) source type (FSRQ: flat spectrum radio quasar, BL: BL Lac object, CSSQ: compact steep spectrum quasar, G: galaxy, PN: planetary nebula), most target sources’ type can be found in [31], except three targets with additional references in footnote; (3) number of flux density measurements; (4) mean flux density with the standard deviation; (5) modulation index m and modified modulation index m ; (6) relative variability amplitude Y; (7), (8) χ 2 and reduced χ 2 ; (9) the null probability (a probability that the light curve can be fitted with a constant value) and with detected IDV (y) or not (n) in the bracket; (10) redshift of source; (11) previously detected IDV (y) or not (n) with relative reference in the bracket, or no available observation/publication found (n.a.). The sources in the first part are from RadioAstron observing list, the sources in the second part are from other samples, the sources in the third part are calibrators which did not show IDV from previous observations.
Table 1. The result of the IDV observation in 8–12 February 2014 with Urumqi 25 m telescope at 4.8 GHz. The columns are: (1) name; (2) source type (FSRQ: flat spectrum radio quasar, BL: BL Lac object, CSSQ: compact steep spectrum quasar, G: galaxy, PN: planetary nebula), most target sources’ type can be found in [31], except three targets with additional references in footnote; (3) number of flux density measurements; (4) mean flux density with the standard deviation; (5) modulation index m and modified modulation index m ; (6) relative variability amplitude Y; (7), (8) χ 2 and reduced χ 2 ; (9) the null probability (a probability that the light curve can be fitted with a constant value) and with detected IDV (y) or not (n) in the bracket; (10) redshift of source; (11) previously detected IDV (y) or not (n) with relative reference in the bracket, or no available observation/publication found (n.a.). The sources in the first part are from RadioAstron observing list, the sources in the second part are from other samples, the sources in the third part are calibrators which did not show IDV from previous observations.
1234567891011
SourceTypeN < S > ± σ s m / m Y χ 2 χ r 2 P null (y/n)Redshifty/n (Ref.)
[Jy][%][%]
0202 + 319FSRQ271.675 ± 0.0201.19/0.953.0274.82.9 1.3 × 10 06 (y) *1.466n.a.
0300 + 470 BL a 212.548 ± 0.0291.14/0.912.8160.13.0 6.8 × 10 06 (y) *0.475nnnn (1)
B2 0340 + 362FSRQ100.452 ± 0.0051.11/0.002.709.51.1 3.9 × 10 01 (n)1.484yyyy (1)
PKS 0507 + 179FSRQ150.584 ± 0.0081.37/0.943.6330.62.2 6.4 × 10 03 (n)0.41yyyn (1)
PKS 0528 + 134FSRQ221.710 ± 0.0341.99/1.835.64149.87.1 1.9 × 10 21 (y)2.06nyny (1)
0642 + 449FSRQ232.553 ± 0.0271.06/0.792.5255.62.5 9.7 × 10 05 (y) *3.396nnnn (1)
0716 + 714BL431.255 ± 0.0161.27/0.993.30112.62.7 2.3 × 10 08 (y)0.3y (2)
0748 + 126FSRQ242.559 ± 0.0572.23/2.116.40237.910.3 1.3 × 10 37 (y)0.889nyyy (1)
0917 + 624FSRQ441.123 ± 0.0181.60/1.374.40186.04.3 1.1 × 10 19 (y)1.446nnnn (1)
0954 + 658FSRQ431.113 ± 0.0171.53/1.294.15164.33.9 2.2 × 10 16 (y)0.368nyyy (1)
1044 + 719FSRQ152.581 ± 0.0351.36/1.153.5857.54.1 3.3 × 10 07 (y)1.15nnyn (1)
1156 + 295FSRQ261.248 ± 0.0131.04/0.702.4545.41.8 7.5 × 10 03 (n)0.72469yyyy (1)
B2 1239 + 376FSRQ240.918 ± 0.0111.20/0.003.0347.82.1 1.8 × 10 03 (n)3.81908n.a
1357 + 769FSRQ510.461 ± 0.0255.42/5.2816.15972.419.4 3.6 × 10 17 (y) *1.585n.a.
PKS 1502 + 106FSRQ90.791 ± 0.0050.63/0.000.004.10.5 8.5 × 10 01 (n)1.83831n (3)
1637 + 574FSRQ391.308 ± 0.0161.22/0.933.1286.02.3 1.4 × 10 05 (y) *0.7506nnnn (1)
WISE J092915.43 + 501336.0 BL b 100.342 ± 0.0133.80/3.5011.2483.39.3 3.6 × 10 14 (y)0.37039yyyy (1)
BZQ J1128 + 5925FSRQ460.621 ± 0.0233.70/3.5410.94616.513.7 2.0 × 10 101 (y)1.795y (4)
1633 + 382FSRQ283.278 ± 0.0220.67/0.000.5629.11.1 3.6 × 10 01 (n)1.81309nyyn (1)
J2311 + 4543 FSRS c 240.662 ± 0.0162.42/2.226.99139.86.1 9.9 × 10 19 (y)1.447yyyy (1)
3C48CSSQ335.590 ± 0.0380.68/0.060.6535.31.1 3.1 × 10 02 0.367
0836 + 710FSRQ523.434 ± 0.0200.58/0.000.0034.70.7 9.6 × 10 01 2.172
0951 + 699 (M82)G493.557 ± 0.0210.59/0.000.0034.00.7 9.4 × 10 01 0.00068
3C286CSSQ287.511 ± 0.0460.61/0.000.0022.90.8 6.9 × 10 01 0.84993
NGC7027PN265.411 ± 0.0410.76/0.371.2034.51.4 9.7 × 10 02 0
Notes. References for sources’ type: a. [32]; b. [33]; c. [34]. References for previously detected IDV or not (Column 11). 1. [4], where the classification per epoch of variable (y) or nonvariable (n) in four epochs (with each of 3 days) observed with the VLA at 5 GHz around year 2002; 2. [5] 3. [3] 4. [35]. The * marked in Column 9 indicates new IDVs found from this observation.
Table 2. The result of the IDV observation in 7–9 March 2014 with Urumqi 25 m telescope at 4.8 GHz. The columns are same as that in Table 1. The * marked in Column 9 indicates new IDVs found from this observation.
Table 2. The result of the IDV observation in 7–9 March 2014 with Urumqi 25 m telescope at 4.8 GHz. The columns are same as that in Table 1. The * marked in Column 9 indicates new IDVs found from this observation.
SourceTypeN < S > ± σ s m / m Y χ 2 χ r 2 P null (y/n)Redshifty/n (Ref.)
[Jy][%][%]
PKS 0528 + 134FSRQ231.525 ± 0.0342.23/2.176.55405.818.4 2.6 × 10 72 (y)2.06nyny (1)
4C + 72.10 (0604 + 728)FSRQ430.622 ± 0.0091.45/1.234.12158.93.8 1.7 × 10 15 (y) *3.53n.a.
0716 + 714BL481.193 ± 0.0171.42/1.314.05322.06.9 1.2 × 10 42 (y)0.3y (2)
0748 + 126FSRQ192.525 ± 0.0431.70/1.634.93224.012.4 1.5 × 10 37 (y)0.889nyyy (1)
OJ287BL192.881 ± 0.0531.84/1.775.35274.015.2 1.1 × 10 47 (y)0.3056y (3)
0917 + 624FSRQ481.112 ± 0.0141.26/1.103.53206.44.4 7.4 × 10 22 (y)1.446nnnn (1)
0953 + 254FSRQ201.149 ± 0.0131.13/0.963.1170.63.7 7.4 × 10 08 (y)0.70721nyyn (1)
0954 + 658FSRQ481.038 ± 0.0131.25/1.093.50230.24.9 5.6 × 10 26 (y)0.368yyyy (1)
1044 + 719FSRQ482.786 ± 0.0491.76/1.685.10599.712.8 6.6 × 10 97 (y)1.15nnyn (1)
B2 1128 + 385FSRQ211.385 ± 0.0201.44/1.334.12144.47.2 7.3 × 10 21 (y)1.74046nyny (1)
1156 + 295FSRQ191.236 ± 0.0110.89/0.692.3048.92.7 1.1 × 10 04 (y)0.72469yyyy (1)
1357 + 769FSRQ380.411 ± 0.0133.16/3.029.39339.19.2 1.8 × 10 50 (y)1.585n.a.
1611 + 343FSRQ213.956 ± 0.0451.14/1.023.13113.25.7 5.1 × 10 15 (y)1.39915nynn (1)
3C345FSRQ216.366 ± 0.0570.90/0.752.3275.43.8 2.4 × 10 08 (y) *0.5928nnnn (1)
WISE J092915.43 + 501336.0BL180.347 ± 0.0257.20/7.1321.57879.251.7 5.9 × 10 176 (y)0.37039yyyy (1)
BZQ J1128 + 5925FSRQ340.630 ± 0.0335.24/5.1815.661721.552.2 0.0 × 10 00 (y)1.795y (4)
J2311 + 4543FSRQ210.671 ± 0.0162.38/2.277.02209.410.5 1.5 × 10 33 (y)1.447yyyy (1)
3C48 275.604 ± 0.0280.50/0.120.6428.51.1 3.3 × 10 01
0836 + 710 483.404 ± 0.0140.41/0.000.0032.90.7 9.4 × 10 01
0951 + 699 (M82) 483.548 ± 0.0140.39/0.000.0028.30.6 9.9 × 10 01
3C286 187.506 ± 0.0270.36/0.000.008.80.5 9.5 × 10 01
NGC7027 255.406 ± 0.0320.59/0.001.1537.21.6 4.2 × 10 02
WISE J092915.43 + 501336.0 (J0929 + 5013, B0925 + 504): the BL Lac object exhibits fast IDV in both sessions (bottom panel of Figure 1), m increases from 3.50 in session 1 to 7.13 in session 2. Its mean flux density is slightly increased from session 1 to session 2. The blazar has previously shown fast IDV with Urumqi 25 m and Efferlsberg 100 m telescope, its IDV timescales can be fitted with the annual modulated ISS model [36,37].
BZQ J1128 + 5925: the blazar shows a strong fast IDV in both sessions (middle panel of Figure 1), its m increases from 3.54 in session 1 to 5.18 in session 2. Its mean flux density increases slightly from 0.621 ± 0.023 Jy in session 1 to 0.630 ± 0.033 Jy in session 2. This fast IDV was first found with Urumqi 25 m telescope, and it has shown seasonal variations of IDV timescales and can be well fitted with the annual modulation model [35].
J2311 + 4543: the blazar shows a strong IDV in both sessions with both the m and mean flux density slightly increased from session 1 to session 2.

3.3. Discussion

We have tried to cross-correlate the parameters listed in Table 1 and Table 2. No correlation is found between the IDV strength ( m ) and either redshift or Galactic latitude, as shown in Figure 3a,b, with correlation parameters listed in Table 3. Three strong IDVs appear in area of relatively lower flux density < S > < 0.7 Jy as shown in Figure 4.
For 11 common target blazars in Table 1 and Table 2, there is no significant correlation between the changes of IDV strength m and the changes of average flux density < S > from session 1 (in February) to session 2 (in March), as the arrows shown in Figure 4 and the comments in Section 3.2.
The detection of an IDV or not from one epoch to another is dependent on many factors, e.g., the changes of physical properties of scattering plasma as well as the changes of radio core of AGN for the ISS model [4], and the telescope sensitivity/uncertainty, calibration procedure/precision, and also the definition of an IDV detection (slightly different criteria from different research groups). We can see that some of our detections (non-detections) of IDV have non-detections (detections) in previous observations e.g., in [4] as shown in Table 1 and Table 2, these occur mostly for weak or marginally detected IDVs. The non-detection of IDV in our session 1 but a weak IDV in session 2 for 1156 + 295 might be partly due to the higher residual uncertainty ( m 0 = 0.64 ) of observation in session 1 than ( m 0 = 0.45 ) in session 2. Except for the observation sensitivity (including weather and interferences) effects and the issue of different definition criteria (or confidence level) of IDV, it is believed that the ISS-induced IDVs are often intermittent in time, e.g., the strong IDV of 0917 + 624 has been quenched after year 2000 [38], although it shows weak IDV in this paper.
In the ISS models, the variability timescale of around a day, requires the angular size of radio-emitting component from sub-milliarcsecond to microarcsecond for the blazar emission to be scintillated by a nearby plasma so as to show an IDV (see [4,19,39]). Emission from the extended or long jet of a blazar is less possible to show an IDV. However, a detailed VLBI-structure and core fraction analysis for these blazars along with the brightness temperature of core components from the high-resolution RadioAstron data is planned in the forthcoming paper.

4. Summary

We present the first results of the IDV observations of about two dozens of blazars with Urumqi 25 m radio telescope in 8–12 February 2014 and 7–9 March. Most of target sources are selected in parallel from the observing list of RadioAstron AGN monitoring campaign. The majority of the blazars exhibit IDV with a high confidence level, some of them show prominent IDVs in the two sessions. Seven new IDV sources are detected in our observations.
We find that the three strongest IDVs appear in relatively lower flux density < S > < 0.7 Jy in the sample, in which blazar 1357 + 769 is discovered as strong IDV for the first time. No significant correlation between the IDV strength and either redshift or Galactic latitude is found from this blazar sample.

Author Contributions

Conceptualization, writing—original draft, X.L.; analysis and plotting, X.W., N.C.; and review and editing, J.L., L.C., X.Y. and T.P.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Key R&D Program of China under grant No. 2018YFA0404602, and partly supported by Xiaofeng Yang’s Xinjiang Tianchi Bairen project and CAS Pioneer Hundred Talents Program, and the CAS ‘Light of West China’ Program under grant No. 2018-XBQNXZ-B-021 and the NSFC under grant No. U2031212 and 61931002. L.C. is also supported by the Youth Innovation Promotion Association of the CAS under grant No. 2017084.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank the referees for their helpful comments. This work is based on observations made with the Urumqi 25 m radio Telescope, which is operated by the Key Laboratory of Radio Astronomy, Chinese Academy of Sciences. The NASA/IPAC Extragalactic Database (NED) is funded by the National Aeronautics and Space Administration and operated by the California Institute of Technology.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Witzel, A.; Heeschen, D.S.; Schalinski, C.; Krichbaum, T. Kurzzeit-Variabilität extragalaktischer Radioquellen. Mitt. Astron. Ges. Hambg. 1986, 65, 239–241. [Google Scholar]
  2. Heeschen, D.S.; Krichbaum, T.P.; Schalinski, C.J.; Witzel, A. Rapid Variability of Extragalactic Radio Sources. Astron. J. 1987, 94, 1493–1507. [Google Scholar] [CrossRef]
  3. Quirrenbach, A.; Witzel, A.; Krichbaum, T.P.; Hummel, C.A.; Wegner, R.; Schalinski, C.J.; Ott, M.; Alberdi, A.; Rioja, M. Statistics of intraday variability in extragalactic radio sources. Astron. Astrophys. 1992, 258, 279–284. [Google Scholar]
  4. Lovell, J.E.J.; Rickett, B.J.; Macquart, J.-P.; Jauncey, D.L.; Bignall, H.E.; Kedziora-Chudczer, L.; Ojha, R.; Pursimo, T.; Dutka, M.; Senkbell, C.; et al. The Micro-Arcsecond Scintillation-Induced Variability (MASIV) Survey. II. The First Four Epochs. Astrophys. J. 2008, 689, 108–126. [Google Scholar] [CrossRef] [Green Version]
  5. Liu, X.; Song, H.-G.; Liu, J.; Ding, Z.; Marchili, N.; Krichbaum, T.P.; Fuhrmann, L.; Zensus, A.; An, T. Radio observations of the first three-months of Fermi AGN at 4.8 GHz. Res. Astron. Astrophys. 2012, 12, 147–157. [Google Scholar] [CrossRef] [Green Version]
  6. Wagner, S.J.; Witzel, A. Intraday Variability In Quasars and BL Lac Objects. Annu. Rev. Astron. Astrophys. 1995, 33, 163–198. [Google Scholar] [CrossRef]
  7. Jauncey, D.; Bignall, H.; Kedziora-Chudczer, L.; Koay, J.; Lovell, J.; Macquart, J.-P.; Ojha, R.; Pursimo, T.; Reynolds, C.; Rickett, B. Interstellar Scintillation and Scattering of Micro-arc-second AGN. Galaxies 2016, 4, 62. [Google Scholar] [CrossRef]
  8. Kellermann, K.I.; Pauliny-Toth, I.I.K. The Spectra of Opaque Radio Sources. Astrophys. J. 1969, 155, L71–L78. [Google Scholar] [CrossRef]
  9. Readhead, A.C.S. Equipartition Brightness Temperature and the Inverse Compton Catastrophe. Astrophys. J. 1994, 426, 51–59. [Google Scholar] [CrossRef]
  10. Lobanov, A. Brightness temperature constraints from interferometric visibilities. Astron. Astrophys. 2015, 574, A84. [Google Scholar] [CrossRef]
  11. Kovalev, Y.Y.; Kardashev, N.S.; Kellermann, K.I.; Lobanov, A.P.; Johnson, M.D.; Gurvits, L.I.; Voitsik, P.A.; Zensus, J.A.; Anderson, J.M.; Bach, U.; et al. RadioAstron Observations of the Quasar 3C273: A Challenge to the Brightness Temperature Limit. Astrophys. J. Lett. 2016, 820, L9. [Google Scholar] [CrossRef] [Green Version]
  12. Kravchenko, E.V.; Gómez, J.L.; Kovalev, Y.Y.; Lobanov, A.P.; Savolainen, T.; Bruni, G.; Fuentes, A.; Anderson, J.M.; Jorstad, S.G.; Marscher, A.P.; et al. Probing the Innermost Regions of AGN Jets and Their Magnetic Fields with RadioAstron. III. Blazar S5 0716+71 at Microarcsecond Resolution. Astrophys. J. 2020, 893, 68. [Google Scholar] [CrossRef]
  13. Zensus, J.A.; Vega-Garcia, L.; Ros, E.; Lobanov, A.P.; Perucho, M.; Bruni, G.; Kovalev, Y.Y. Imaging strong blazars with space VLBI. Adv. Space Res. 2020, 65, 725–730. [Google Scholar] [CrossRef]
  14. Cohen, M.H.; Lister, M.L.; Homan, D.C.; Kadler, M.; Kellermann, K.I.; Kovalev, Y.Y.; Vermeulen, R.C. Relativistic Beaming and the Intrinsic Properties of Extragalactic Radio Jets. Astrophys. J. 2007, 658, 232–244. [Google Scholar] [CrossRef] [Green Version]
  15. Liodakis, I.; Marchili, N.; Angelakis, E.; Fuhrmann, L.; Nestoras, I.; Myserlis, I.; Karamanavis, V.; Krichbaum, T.P.; Sievers, A.; Ungerechts, H.; et al. F-GAMMA: Variability Doppler factors of blazars from multiwavelength monitoring. Mon. Not. R. Astron. Soc. 2017, 466, 4625–4632. [Google Scholar] [CrossRef]
  16. Liodakis, I.; Hovatta, T.; Huppenkothen, D.; Kiehlmann, S.; Max-Moerbeck, W.; Readhead, A.C.S. Constraining the Limiting Brightness Temperature and Doppler Factors for the Largest Sample of Radio-bright Blazars. Astrophys. J. 2018, 866, 137. [Google Scholar] [CrossRef] [Green Version]
  17. Jauncey, D.L.; Koay, J.Y.; Bignall, H.; Macquart, J.-P.; Pursimo, T.; Giroletti, M.; Hovatta, T.; Kiehlmann, S.; Rickett, B.; Reahead, A.; et al. Interstellar scintillation, ISS, and intrinsic variability of radio AGN. Adv. Space Res. 2020, 65, 756–762. [Google Scholar] [CrossRef]
  18. Dennett-Thorpe, J.; de Bruyn, A.G. Annual modulation in the scattering of J1819+3845: Peculiar plasma velocity and anisotropy. Astron. Astrophys. 2003, 404, 113–132. [Google Scholar] [CrossRef] [Green Version]
  19. Bignall, H.E.; Macquart, J.-P.; Jauncey, D.L.; Lovell, J.E.J.; Tzioumis, A.K.; Kedziora-Chudczer, L. Rapid Interstellar Scintillation of PKS 1257-326: Two-Station Pattern Time Delays and Constraints on Scattering and Microarcsecond Source Structure. Astrophys. J. 2006, 652, 1050–1058. [Google Scholar] [CrossRef] [Green Version]
  20. Liu, X.; Song, H.-G.; Marchili, N.; Liu, B.-R.; Krichbaum, T.P.; Fuhrmann, L.; Zensus, J.A. Intra-day variability observations of S5 0716+714 over 4.5 years at 4.8 GHz. Astron. Astrophys. 2012, 543, A78. [Google Scholar] [CrossRef] [Green Version]
  21. Bignall, H.; Reynolds, C.; Stevens, J.; Bannister, K.; Johnston, S.; Tuntsov, A.V.; Walker, M.A.; Gulyaev, S.; Natusch, T.; Weston, S.; et al. Spica and the annual cycle of PKS B1322–110 scintillations. Mon. Not. R. Astron. Soc. 2019, 487, 4372–4381. [Google Scholar] [CrossRef] [Green Version]
  22. Liu, J.; Bignall, H.; Krichbaum, T.P.; Liu, X.; Kraus, A.; Kovalev, Y.; Sokolovsky, K.; Angelakis, E.; Zensus, J. Effelsberg Monitoring of a Sample of RadioAstron Blazars: Analysis of Intra-Day Variability. Galaxies 2018, 6, 49. [Google Scholar] [CrossRef] [Green Version]
  23. Kovalev, Y.Y.; Kardashev, N.S.; Sokolovsky, K.V.; Voitsik, P.A.; An, T.; Anderson, J.M.; Andrianov, A.S.; Avdeev, V.Y.; Bartel, N.; Bignall, H.E.; et al. Detection statistics of the RadioAstron AGN survey. Adv. Space Res. 2020, 65, 705–711. [Google Scholar] [CrossRef] [Green Version]
  24. Baars, J.W.M.; Genzel, R.; Pauliny-Toth, I.I.K.; Witzel, A. The Absolute Spectrum of Cas A; An Accurate Flux Density Scale and a Set of Secondary Calibrators. Astron. Astrophys. 1977, 61, 99–106. [Google Scholar]
  25. Ott, M.; Witzel, A.; Quirrenbach, A.; Krichbaum, T.P.; Standke, K.J.; Schalinski, C.J.; Hummel, C.A. An updated list of radio flux density calibrators. Astron. Astrophys. 1994, 284, 331–339. [Google Scholar]
  26. Schleicher, B.; Arbet-Engels, A.; Baack, D.; Balbo, M.; Biland, A.; Blank, M.; Bretz, T.; Bruegge, K.; Bulinski, M.; Buss, J.; et al. Fractional Variability—A Tool to Study Blazar Variability. Galaxies 2019, 7, 62. [Google Scholar] [CrossRef] [Green Version]
  27. Kraus, A.; Krichbaum, T.P.; Wegner, R.; Witzel, A.; Cimò, G.; Quirrenbach, A.; Britzen, S.; Fuhrmann, L.; Lobanov, A.P.; Naundorf, C.E.; et al. Intraday variability in compact extragalactic radio sources. II. Observations with the Effelsberg 100 m radio telescope. Astron. Astrophys. 2003, 401, 161–172. [Google Scholar] [CrossRef]
  28. Lister, M.L.; Aller, M.F.; Aller, H.D.; Homan, D.C.; Kellermann, K.I.; Kovalev, Y.Y.; Pushkarev, A.B.; Richards, J.L.; Ros, E.; Savolainen, T. MOJAVE: XIII. Parsec-scale AGN Jet Kinematics Analysis Based on 19 years of VLBA Observations at 15 GHz. Astron. J. 2016, 152, 12. [Google Scholar] [CrossRef] [Green Version]
  29. Gabányi, K.É.; Marchili, N.; Krichbaum, T.P.; Fuhrmann, L.; Müller, P.; Zensus, J.A.; Liu, X.; Song, H.G. Very Long Baseline Array observations of the intraday variable source J1128+592. Astron. Astrophys. 2009, 508, 161–171. [Google Scholar] [CrossRef]
  30. Valtonen, M.J. New Orbit Solutions for the Precessing Binary Black Hole Model of OJ 287. Astrophys. J. 2007, 659, 1074–1081. [Google Scholar] [CrossRef]
  31. Healey, S.E.; Romani, R.W.; Taylor, G.B.; Sadler, E.M.; Ricci, R.; Murphy, T.; Ulvestad, J.S.; Winn, J.N. CRATES: An All-Sky Survey of Flat-Spectrum Radio Sources. Astrophys. J. Suppl. 2007, 171, 61–71. [Google Scholar] [CrossRef] [Green Version]
  32. Véron-Cetty, M.-P.; Véron, P. A catalogue of quasars and active nuclei: 12th edition. Astron. Astrophys. 2006, 455, 773–777. [Google Scholar] [CrossRef] [Green Version]
  33. Plotkin, R.M.; Anderson, S.F.; Hall, P.B.; Margon, B.; Voges, W.; Schneider, D.P.; Stinson, G.; York, D.G. A Large Sample of BL Lac Objects from the SDSS and First. Astron. J. 2008, 135, 2453–2469. [Google Scholar] [CrossRef] [Green Version]
  34. Lazio, T.J.W.; Ojha, R.; Fey, A.L.; Kedziora-Chudczer, L.; Cordes, J.M.; Jauncey, D.L.; Lovell, J.E.J. Angular Broadening of Intraday Variable AGNs. II. Interstellar and Intergalactic Scattering. Astrophys. J. 2008, 672, 115–121. [Google Scholar] [CrossRef] [Green Version]
  35. Gabányi, K.É.; Marchili, N.; Krichbaum, T.P.; Britzen, S.; Fuhrmann, L.; Witzel, A.; Zensus, J.A.; Müller, P.; Liu, X.; Song, H.G.; et al. The IDV source J 1128+5925, a new candidate for annual modulation? Astron. Astrophys. 2007, 470, 83–95. [Google Scholar] [CrossRef] [Green Version]
  36. Liu, J.; Liu, X. Rapid variability of BL Lac 0925+504: Interstellar scintillation induced? Astrophys. Space Sci. 2015, 357, 165. [Google Scholar] [CrossRef] [Green Version]
  37. Liu, X.; Yang, P.P.; Liu, J.; Liu, B.R.; Hu, S.M.; kurtanidze, O.M.; Zola, S.; Kraus, A.; Krichbaum, T.P.; Su, R.Z.; et al. Radio and optical intra-day variability observations of five blazars. Mon. Not. R. Astron. Soc. 2017, 469, 2457–2463. [Google Scholar] [CrossRef]
  38. Liu, X.; Mi, L.-G.; Liu, J.; Cui, L.; Song, H.-G.; Krichbaum, T.P.; Kraus, A.; Fuhrmann, L.; Marchili, N.; Zensus, J.A. Intra-day variability observations and the VLBI structure analysis of quasar S4 0917+624. Astron. Astrophys. 2015, 578, 34. [Google Scholar] [CrossRef] [Green Version]
  39. Koay, J.Y.; Bignall, H.E.; Macquart, J.-P.; Jauncey, D.L.; Rickett, B.J.; Lovell, J.E.J. Detection of six rapidly scintillating active galactic nuclei and the diminished variability of J1819+3845. Astron. Astrophys. 2011, 534, L1. [Google Scholar] [CrossRef]
Figure 1. Three strong IDV sources with m 3 in both observations. (Left upper panel) the light curve of 1357 + 769 at 4.8 GHz in 8–12 February 2014, the flux density I (left y-axis) and the ratio of flux density over mean flux density < I > (right y-axis) versus time (MJD) in this and following figures; (right upper panel) the light curve of 1357 + 769 at 4.8 GHz in 7–9 March 2014. (Left middle panel) the light curve of J1128 + 5925 at 4.8 GHz in 8–12 February 2014; (right middle panel) the light curve of J1128 + 5925 at 4.8 GHz in 7–9 March 2014. (Left bottom panel) the light curve of J0929 + 5013 at 4.8 GHz in 8–12 February 2014; (right bottom panel) the light curve of J0929 + 5013 at 4.8 GHz in 7–9 March 2014.
Figure 1. Three strong IDV sources with m 3 in both observations. (Left upper panel) the light curve of 1357 + 769 at 4.8 GHz in 8–12 February 2014, the flux density I (left y-axis) and the ratio of flux density over mean flux density < I > (right y-axis) versus time (MJD) in this and following figures; (right upper panel) the light curve of 1357 + 769 at 4.8 GHz in 7–9 March 2014. (Left middle panel) the light curve of J1128 + 5925 at 4.8 GHz in 8–12 February 2014; (right middle panel) the light curve of J1128 + 5925 at 4.8 GHz in 7–9 March 2014. (Left bottom panel) the light curve of J0929 + 5013 at 4.8 GHz in 8–12 February 2014; (right bottom panel) the light curve of J0929 + 5013 at 4.8 GHz in 7–9 March 2014.
Universe 07 00015 g001
Figure 2. Two example sources with no variability. (Left panel) the light curve of B2 1239 + 376 at 4.8 GHz in 8–12 February 2014; (right panel) the light curve of 1633 + 382 at 4.8 GHz in 8–12 February 2014.
Figure 2. Two example sources with no variability. (Left panel) the light curve of B2 1239 + 376 at 4.8 GHz in 8–12 February 2014; (right panel) the light curve of 1633 + 382 at 4.8 GHz in 8–12 February 2014.
Universe 07 00015 g002
Figure 3. (a) The modified modulation index ( m in per cent) versus redshift (z) of the target sources observed in session 1 (blue circle) and session 2 (red square); (b) the modified modulation index ( m in per cent) versus Galactic latitude of the target sources observed in session 1 (blue circle) and session 2 (red square).
Figure 3. (a) The modified modulation index ( m in per cent) versus redshift (z) of the target sources observed in session 1 (blue circle) and session 2 (red square); (b) the modified modulation index ( m in per cent) versus Galactic latitude of the target sources observed in session 1 (blue circle) and session 2 (red square).
Universe 07 00015 g003
Figure 4. The modified modulation index ( m in per cent) versus average flux density ( < S > , in Jy) of the target sources observed in session 1 (blue circle) and session 2 (red square). Arrows indicate from session 1 to session 2 for 11 common target sources.
Figure 4. The modified modulation index ( m in per cent) versus average flux density ( < S > , in Jy) of the target sources observed in session 1 (blue circle) and session 2 (red square). Arrows indicate from session 1 to session 2 for 11 common target sources.
Universe 07 00015 g004
Table 3. Results of correlation analysis for modified modulation index ( m ), redshift (z) and Galactic latitude. Column 1: Sample (Full—all target sources in Table 1 and Table 2; S1—target sources in Table 1; S2—target sources in Table 2). Column 2 and 3: Pearson correlation coefficient and the probability of no correlation. Column 4 and 5: Spearman correlation coefficient and the probability of no correlation.
Table 3. Results of correlation analysis for modified modulation index ( m ), redshift (z) and Galactic latitude. Column 1: Sample (Full—all target sources in Table 1 and Table 2; S1—target sources in Table 1; S2—target sources in Table 2). Column 2 and 3: Pearson correlation coefficient and the probability of no correlation. Column 4 and 5: Spearman correlation coefficient and the probability of no correlation.
SamplePearson p null Spearman p null
12345
Modified modulation index ( m ) versus redshift (z)
Full−0.1210.475−0.0050.978
S1−0.1730.465−0.2420.305
S2−0.0320.9030.2790.277
Modified modulation index ( m ) versus Galactic latitude
Full0.0160.926−0.1280.449
S1−0.0450.850−0.1220.607
S2−0.0380.886−0.2720.290
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, X.; Wang, X.; Chang, N.; Liu, J.; Cui, L.; Yang, X.; Krichbaum, T.P. Intra-Day Variability Observations of Two Dozens of Blazars at 4.8 GHz. Universe 2021, 7, 15. https://doi.org/10.3390/universe7010015

AMA Style

Liu X, Wang X, Chang N, Liu J, Cui L, Yang X, Krichbaum TP. Intra-Day Variability Observations of Two Dozens of Blazars at 4.8 GHz. Universe. 2021; 7(1):15. https://doi.org/10.3390/universe7010015

Chicago/Turabian Style

Liu, Xiang, Xin Wang, Ning Chang, Jun Liu, Lang Cui, Xiaofeng Yang, and Thomas P. Krichbaum. 2021. "Intra-Day Variability Observations of Two Dozens of Blazars at 4.8 GHz" Universe 7, no. 1: 15. https://doi.org/10.3390/universe7010015

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