Articles

CHARACTERISTICS OF GAMMA-RAY LOUD BLAZARS IN THE VLBA IMAGING AND POLARIMETRY SURVEY

, , , , , , , , and

Published 2010 December 9 © 2011. The American Astronomical Society. All rights reserved.
, , Citation J. D. Linford et al 2011 ApJ 726 16 DOI 10.1088/0004-637X/726/1/16

0004-637X/726/1/16

ABSTRACT

The radio properties of blazars detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope have been observed as part of the VLBA Imaging and Polarimetry Survey. This large, flux-limited sample of active galactic nuclei (AGNs) provides insights into the mechanism that produces strong γ-ray emission. At lower flux levels, radio flux density does not directly correlate with γ-ray flux. We find that the LAT-detected BL Lac objects tend to be similar to the non-LAT BL Lac objects, but that the LAT-detected FSRQs are often significantly different from the non-LAT FSRQs. The differences between the γ-ray loud and quiet FSRQs can be explained by Doppler boosting; these objects appear to require larger Doppler factors than those of the BL Lac objects. It is possible that the γ-ray loud FSRQs are fundamentally different from the γ-ray quiet FSRQs. Strong polarization at the base of the jet appears to be a signature for γ-ray loud AGNs.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

Several outstanding questions about γ-ray loud active galactic nuclei (AGNs) remain, including: "why are some sources γ-ray loud and others are γ-ray quiet?" The answer to this is likely related to Doppler boosting (Jorstad et al. 2001; Taylor et al. 2007; Kovalev et al. 2009; Lister et al. 2009a; Savolainen et al. 2010), though many details remain to be filled in. A related question is "where are the γ-rays produced: in or near the core or throughout the jets?" There are arguments for both locations (Sikora et al. 2009; Tavecchio et al. 2010; Marscher 2006; Marscher et al. 2008, 2010). In either case it still is unclear how the radio flaring, or the ejection of new components, is related to the gamma-ray activity, or whether there are multiple classes of γ-ray blazars (e.g., Sikora et al. 2001, 2002). In addition, Stawarz et al. (2008) argue that we should detect γ-rays from the compact (<kpc) lobes of young radio sources.

The Large Area Telescope (LAT; Atwood et al. 2009) on board the Fermi Gamma-ray Space Telescope is a wide-field telescope covering the energy range from about 20 MeV to more than 300 GeV. At high galactic latitudes 80% (106 of 132) of the γ-ray bright sources detected in the LAT Bright Source List derived from the first three months of Fermi observations (Abdo et al. 2009a) are associated with known AGNs. This trend continues in the LAT first-year catalog (1FGL; Abdo et al. 2010a), where 662 of 1043 high latitude sources are associated with AGNs with high confidence (P > 80%). The vast majority of the LAT γ-ray sources are blazars, with strong, compact radio emission. These blazars exhibit flat radio spectra, rapid variability, compact cores with one-sided parsec-scale jets, and superluminal motion in the jets (Marscher 2006). A handful of radio galaxies and Seyferts have also been identified in the First LAT AGN Catalog (1LAC; Abdo et al. 2010b). These radio galaxies (e.g., 3C84, M87, Cen A) have faint parsec-scale counterjets indicating that they are viewed at larger angles to the line of sight than the blazars. Also, there is a much larger fraction of unclassified sources in the one year catalog than was present in the three-month catalog.

The VLBA Imaging and Polarimetry Survey (VIPS; Helmboldt et al. 2007) was originally devised to obtain Very Long Baseline Array (VLBA) measurements of AGNs in the fifth data release of the Sloan Digital Sky Survey (DR5; Adelman-McCarthy et al. 2006). It consists of 5 GHz VLBA images and polarization data of 1127 AGNs located between declinations of 15° and 65°. The VIPS sources were drawn from the Cosmic Lens All-Sky Survey (CLASS; Myers et al. 2003). All VIPS sources were required to have CLASS flux densities of 85 mJy or more. The VLBA at 5 GHz gives us an angular resolution of just under 2 mas, which allows us to probe sub-parsec scale features of the objects. Using multi-epoch observations, we will be able to measure the apparent jet velocities.

Various programs, such as MOJAVE observing at 15 GHz (Lister et al. 2009a; Homan et al. 2009), the Boston University program observing at 22 and 43 GHz (Roca-Sogorb et al. 2010), and TANAMI (Ojha et al. 2010) observing at 8.4 and 22 GHz, are in place to monitor the radio jets from the brightest blazars such as 3C273, BL Lacertae, etc. The detection of Centaurus A, 3C84 and M87 by Fermi (Abdo et al. 2009a, 2009b, 2009c) and the detection of M87 at TeV energies (Aharonian et al. 2006) suggests that other classes of AGNs (e.g., radio galaxies) may form a significant population of γ-ray emitters. Many of the current programs are strongly biased toward the radio bright sample. VIPS has a flux limit roughly an order of magnitude below the MOJAVE survey and most other samples that have been used to study the properties of γ-ray blazars and so allows us to probe the extension of the radio core/γ-ray properties down to a fainter population. Many different radio-γ-ray correlations have been suggested (Taylor et al. 2007; Abdo et al. 2009a; Lister et al. 2009a; Kovalev et al. 2009; Giroletti et al. 2010). By examining a larger sample (1127 objects) we attempt to obtain more definitive insight into the properties of γ-ray loud blazars.

In Section 2, we define our sample. In Section 3, we present data on the γ-ray and radio properties of the LAT sources. In Section 4, we compare several parameters of the LAT and non-LAT sources. Notes on some individual objects of interest are presented in Section 5 and in Section 6 we discuss our results. Some conclusions are given in Section 7. Throughout this paper we assume H0 = 71 km s−1 Mpc−1 and ΛCDM cosmology (e.g., Hinshaw et al. 2009).

2. SAMPLE DEFINITION

We have compared the 11 month 1LAC to the VIPS (Helmboldt et al. 2007), and from the intersection compiled a sample of 102 LAT sources at |b|>10° associated with an AGN with a probability of more than 50% (although most of the probabilities are greater than 98%). All of these objects were observed with the VLBA from 1996 to 2006 (although 958 of the 1127 sources were observed during 2006) at 5 GHz. Each VIPS source was observed for approximately 500 s, with scans spread out over 11 hr. See Table 1 for a summary of our data. Of the LAT objects, 40 are BL Lac objects, 50 are Flat Spectrum Radio Quasars (FSRQs), and 12 are either radio galaxies or unclassified objects. The optical classifications were adopted from the Candidate Gamma-Ray Blazar Survey (CGRaBS; Healey et al. 2008), except that objects labeled in CGRaBS as "AGN" we call radio galaxies.

Table 1. LAT Sources in VIPS

1FGL Name VIPS Name R.A. Decl. Probability Type z F35 Δ F35 Core TB Opening Angle ΔP.A. S (8.5 GHz)
1FGL J0742.2+5443 J07426+5444 07:42:39.8 +54:44:24.679 1 FSRQ 0.723 2.2 0.3 5.26E+10 ... ... 142.8
1FGL J0746.6+2548 J07464+2549 07:46:25.9 +25:49:02.146 0.98 FSRQ 2.979 0.7 0.2 5.07E+11 ... ... 291.7
1FGL J0752.8+5353 J07530+5352 07:53:01.4 +53:52:59.636 1 BLL 0.2 1.2 0.3 6.79E+11 ... ... 1196.9
1FGL J0800.5+4407 J08011+4401 08:01:08.3 +44:01:10.148 0.92 ... ... 0.7 0.3 7.36E+10 ... ... 230.6
1FGL J0806.2+6148 J08053+6144 08:05:18.2 +61:44:23.704 0.94 FSRQ 3.033 0.6 0.2 1.30E+11 ... ... 725.7
1FGL J0809.4+3455 J08096+3455 08:09:38.9 +34:55:37.248 0.99 BLL 0.082 0.7 0.3 2.48E+10 ... ... 152.1
1FGL J0809.5+5219 J08098+5218 08:09:49.2 +52:18:58.252 0.99 BLL 0.138 2.4 0.3 4.83E+09 ... ... 154.2
1FGL J0815.0+6434 J08146+6431 08:14:39.2 +64:31:22.040 0.99 BLL ... 1.7 0.3 2.79E+10 ... ... 221.6
1FGL J0816.7+5739 J08163+5739 08:16:23.8 +57:39:09.509 1 BLL ... 0.7 0.3 4.16E+07 52 8 92
1FGL J0818.2+4222 J08182+4222 08:18:16.0 +42:22:45.408 1 BLL ... 8.7 0.6 3.76E+10 ... ... 1046.4
1FGL J0825.0+5555 J08247+5552 08:24:47.2 +55:52:42.662 0.99 FSRQ 1.417 0.9 0.3 7.86E+10 15 −7 1667.8
1FGL J0834.4+4221 J08338+4224 08:33:53.9 +42:24:01.859 0.98 FSRQ 0.249 0.9 0.2 8.66E+09 ... ... 560.8
1FGL J0856.6+2103 J08566+2057 08:56:39.7 +20:57:43.426 0.94 BLL ... 1 0.3 8.98E+09 ... ... 132.7
1FGL J0856.6+2103 J08569+2111 08:56:57.2 +21:11:43.640 0.96 FSRQ 2.098 1 0.3 1.70E+10 ... ... 332.8
1FGL J0910.7+3332 J09106+3329 09:10:37.0 +33:29:24.418 1 BLL 0.354 0.8 0.3 7.40E+10 ... ... 88.6
1FGL J0911.0+2247 J09107+2248 09:10:42.1 +22:48:35.565 0.98 FSRQ 2.661 2 0.3 8.82E+10 ... ... 178.8
1FGL J0912.3+4127 J09121+4126 09:12:11.6 +41:26:09.356 0.98 FSRQ 2.563 0.6 0.2 8.06E+09 ... ... 166.9
1FGL J0915.7+2931 J09158+2933 09:15:52.4 +29:33:23.982 1 BLL ... 2.1 0.3 1.87E+10 37 6 186.1
1FGL J0919.6+6216 J09216+6215 09:21:36.2 +62:15:52.185 0.82 FSRQ 1.446 1.1 0.3 2.61E+11 ... ... 1532
1FGL J0920.9+4441 J09209+4441 09:20:58.5 +44:41:53.988 1 FSRQ 2.19 14 0.7 8.58E+10 ... ... 1368.3
1FGL J0923.2+4121 J09235+4125 09:23:31.3 +41:25:27.429 0.97 RG 0.028 1.2 0.3 4.91E+09 27 −8 235.2
1FGL J0924.2+2812 J09238+2815 09:23:51.5 +28:15:24.966 0.97 FSRQ 0.744 1.3 0.3 9.87E+11 ... ... 218.4
1FGL J0929.4+5000 J09292+5013 09:29:15.4 +50:13:35.982 0.67 BLL ... 0.7 0.2 6.67E+10 ... ... 747.8
1FGL J0934.5+3929 J09341+3926 09:34:06.7 +39:26:32.125 0.93 BLL ... 0.6 0.2 1.51E+10 ... ... 188.3
1FGL J0937.7+5005 J09372+5008 09:37:12.3 +50:08:52.082 0.97 FSRQ 0.276 0.8 0.2 6.26E+10 ... ... 372.5
1FGL J0941.2+2722 J09418+2728 09:41:48.1 +27:28:38.818 0.97 FSRQ 1.306 0.9 0 1.94E+10 ... ... 268.4
1FGL J0949.8+1757 J09496+1752 09:49:39.8 +17:52:49.432 0.98 FSRQ 0.693 0.9 0 1.27E+11 ... ... 153
1FGL J0956.9+2513 J09568+2515 09:56:49.9 +25:15:16.047 0.99 FSRQ 0.712 0.7 0.2 1.55E+10 ... ... 1905
1FGL J0957.7+5523 J09576+5522 09:57:38.2 +55:22:57.740 1 FSRQ 0.896 10.5 0.6 1.20E+08 60 −33 1498.9
1FGL J1000.9+2915 J10011+2911 10:01:10.2 +29:11:37.543 0.99 BLL ... 1.4 0.3 3.59E+10 ... ... 262.1
1FGL J1015.1+4927 J10150+4926 10:15:04.1 +49:26:00.704 1 BLL 0.2 8.7 0.6 8.16E+09 16 −1 252.2
1FGL J1023.6+3937 J10231+3948 10:23:11.6 +39:48:15.378 0.89 FSRQ 1.254 0.6 0.2 9.53E+10 16 9 779.6
1FGL J1033.2+4116 J10330+4116 10:33:03.7 +41:16:06.234 0.99 FSRQ 1.117 1.1 0.3 4.84E+10 49 163 381.4
1FGL J1033.8+6048 J10338+6051 10:33:51.4 +60:51:07.342 1 FSRQ 1.401 2.2 0.3 1.16E+10 ... ... 430.4
1FGL J1043.1+2404 J10431+2408 10:43:09.0 +24:08:35.430 0.99 BLL 0.56 1.5 0.3 1.56E+10 ... ... 685.4
1FGL J1058.6+5628 J10586+5628 10:58:37.7 +56:28:11.180 1 BLL 0.143 5.7 0.5 9.49E+09 ... ... 189.8
1FGL J1104.4+3812 J11044+3812 11:04:27.3 +38:12:31.794 1 BLL 0.03 26.1 1 2.13E+11 15 −1 631.6
1FGL J1106.5+2809 J11061+2812 11:06:07.3 +28:12:47.045 0.97 ... ... 0.9 0.2 5.08E+10 ... ... 193.2
1FGL J1112.8+3444 J11126+3446 11:12:38.8 +34:46:39.124 0.99 FSRQ 1.949 0.9 0.3 6.95E+10 ... ... 160.5
1FGL J1123.9+2339 J11240+2336 11:24:02.7 +23:36:45.876 1 BLL ... 0.9 0 1.07E+10 ... ... 472.7
1FGL J1141.8+1549 J11421+1547 11:42:07.7 +15:47:54.202 0.99 ... ... 1 0.3 5.62E+10 20 −15 274
1FGL J1146.8+4004 J11469+3958 11:46:58.3 +39:58:34.307 0.98 FSRQ 1.089 1 0.3 4.80E+11 ... ... 571.5
1FGL J1150.2+2419 J11503+2417 11:50:19.2 +24:17:53.852 1 BLL 0.2 1 0.3 7.60E+10 17 7 745.7
1FGL J1151.6+5857 J11514+5859 11:51:24.7 +58:59:17.552 1 BLL ... 0.7 0.3 3.08E+09 35 10 102.2
1FGL J1152.1+6027 J11540+6022 11:54:04.5 +60:22:20.785 0.73 ... ... 0.6 0.3 1.71E+11 ... ... 151.6
1FGL J1202.9+6032 J12030+6031 12:03:03.5 +60:31:19.129 1 RG 0.065 1.1 0.3 7.22E+09 21 12 191.4
1FGL J1209.3+5444 J12089+5441 12:08:54.3 +54:41:58.190 0.97 FSRQ 1.344 0.6 0.2 1.17E+11 ... ... 262
1FGL J1209.4+4119 J12093+4119 12:09:22.8 +41:19:41.360 1 BLL ... 0.5 0.2 2.99E+10 ... ... 485.1
1FGL J1209.7+1806 J12098+1810 12:09:51.8 +18:10:06.796 0.97 FSRQ 0.845 0.8 0.2 4.41E+10 ... ... 136.3
1FGL J1217.7+3007 J12178+3007 12:17:52.1 +30:07:00.625 1 BLL 0.13 6.7 0.6 5.00E+10 19 0 334.9
1FGL J1220.2+3432 J12201+3431 12:20:08.3 +34:31:21.711 1 BLL ... 0.6 0.2 3.05E+10 ... ... 314.1
1FGL J1221.5+2814 J12215+2813 12:21:31.7 +28:13:58.497 1 BLL 0.102 6.9 0.5 2.30E+10 27 −12 893.1
1FGL J1225.8+4336 J12248+4335 12:24:51.5 +43:35:19.276 0.94 ... ... 0.8 0 1.80E+10 ... ... 220.8
1FGL J1225.8+4336 J12269+4340 12:26:57.9 +43:40:58.438 0.87 FSRQ 2.002 0.8 0 1.07E+10 ... ... 145.1
1FGL J1228.2+4855 J12288+4858 12:28:51.8 +48:58:01.293 0.91 FSRQ 1.722 1.3 0.3 1.89E+12 17 −6 265.5
1FGL J1230.4+2520 J12302+2518 12:30:14.1 +25:18:07.145 0.99 BLL 0.135 1.4 0.3 1.96E+10 30 −9 326.4
1FGL J1248.2+5820 J12483+5820 12:48:18.8 +58:20:28.725 1 BLL ... 4.5 0.4 1.54E+10 43 15 324.4
1FGL J1253.0+5301 J12531+5301 12:53:11.9 +53:01:11.741 1 BLL ... 3 0.4 3.05E+10 14 5 366
1FGL J1258.3+3227 J12579+3229 12:57:57.2 +32:29:29.321 0.98 FSRQ 0.806 0.8 0 2.26E+11 ... ... 477.1
1FGL J1303.0+2433 J13030+2433 13:03:03.2 +24:33:55.684 1 BLL ... 3.5 0.4 7.16E+10 ... ... 135.7
1FGL J1308.5+3550 J13083+3546 13:08:23.7 +35:46:37.160 0.99 FSRQ 1.055 2.2 0.3 2.15E+11 60 −4 527
1FGL J1312.4+4827 J13127+4828 13:12:43.4 +48:28:30.928 0.99 FSRQ 0.501 1.4 0.3 2.24E+10 ... ... 250.3
1FGL J1314.7+2346 J13147+2348 13:14:43.8 +23:48:26.701 1 BLL ... 1.4 0.3 6.62E+09 ... ... 158.2
1FGL J1317.8+3425 J13176+3425 13:17:36.5 +34:25:15.921 0.98 FSRQ 1.05 0.5 0.2 1.54E+10 ... ... 541.3
1FGL J1321.1+2214 J13211+2216 13:21:11.2 +22:16:12.098 0.99 FSRQ 0.943 1.1 0.3 1.83E+10 ... ... 323.6
1FGL J1331.0+5202 J13307+5202 13:30:42.6 +52:02:15.448 0.99 RG 0.688 0.5 0 5.43E+10 ... ... 171.9
1FGL J1332.9+4728 J13327+4722 13:32:45.2 +47:22:22.653 0.96 FSRQ 0.669 0.5 0.2 1.02E+11 ... ... 457.5
1FGL J1333.2+5056 J13338+5057 13:33:53.8 +50:57:35.914 0.93 ... ... 1.4 0.3 6.24E+09 ... ... 107.6
1FGL J1345.4+4453 J13455+4452 13:45:33.2 +44:52:59.581 0.99 FSRQ 2.534 1.6 0.3 2.01E+11 ... ... 197.5
1FGL J1351.0+3035 J13508+3034 13:50:52.7 +30:34:53.582 1 FSRQ 0.714 0.6 0.2 3.31E+10 ... ... 719.9
1FGL J1359.1+5539 J13590+5544 13:59:05.7 +55:44:29.362 0.98 FSRQ 1.014 0.9 0.2 9.56E+10 ... ... 143.8
1FGL J1421.0+5421 J14197+5423 14:19:46.6 +54:23:14.783 0.71 BLL 0.152 0.9 0.2 5.81E+10 20 4 2247.3
1FGL J1426.9+2347 J14270+2348 14:27:00.4 +23:48:00.045 1 BLL ... 10.2 0.6 8.28E+09 73 −84 241.6
1FGL J1433.9+4204 J14340+4203 14:34:05.7 +42:03:16.010 1 FSRQ 1.24 0.7 0.2 4.29E+10 ... ... 286.8
1FGL J1436.9+2314 J14366+2321 14:36:41.0 +23:21:03.297 0.96 FSRQ 1.545 0.5 0.2 1.40E+11 ... ... 650.7
1FGL J1438.7+3711 J14388+3710 14:38:53.6 +37:10:35.408 0.99 FSRQ 2.401 0.7 0.2 1.45E+09 35 22 315.4
1FGL J1451.0+5204 J14509+5201 14:51:00.0 +52:01:11.700 0.99 BLL ... 0.9 0.3 2.86E+09 ... ... 86.9
1FGL J1454.6+5125 J14544+5124 14:54:27.1 +51:24:33.734 1 BLL ... 1.1 0.3 6.16E+09 ... ... 96.3
1FGL J1503.3+4759 J15037+4759 15:03:48.0 +47:59:31.024 0.96 BLL ... 0.6 0.2 4.09E+10 ... ... 94.1
1FGL J1505.8+3725 J15061+3730 15:06:09.5 +37:30:51.128 0.98 FSRQ 0.674 0.7 0.2 1.40E+11 18 4 777.3
1FGL J1516.9+1928 J15169+1932 15:16:56.8 +19:32:13.010 1 BLL ... 0.6 0.3 3.94E+10 ... ... 609.3
1FGL J1522.1+3143 J15221+3144 15:22:10.0 +31:44:14.427 1 FSRQ 1.487 15.9 0.8 1.45E+11 ... ... 493.8
1FGL J1539.7+2747 J15396+2744 15:39:39.1 +27:44:38.288 1 FSRQ 2.19 1 0.2 1.29E+11 ... ... 291.5
1FGL J1542.9+6129 J15429+6129 15:42:56.9 +61:29:55.358 1 BLL ... 5.2 0.4 5.94E+10 ... ... 144
1FGL J1558.9+5627 J15588+5625 15:58:48.3 +56:25:14.108 1 BLL 0.3 1.8 0.3 1.30E+10 ... ... 169.7
1FGL J1604.3+5710 J16046+5714 16:04:37.4 +57:14:36.668 0.98 FSRQ 0.72 1.2 0.3 1.61E+10 ... ... 485.3
1FGL J1607.1+1552 J16071+1551 16:07:06.4 +15:51:34.495 1 RG 0.496 2.4 0.4 2.80E+10 13 −1 223.5
1FGL J1616.1+4637 J16160+4632 16:16:03.8 +46:32:25.239 0.96 FSRQ 0.95 0.6 0.2 2.77E+11 ... ... 125.5
1FGL J1635.0+3808 J16352+3808 16:35:15.5 +38:08:04.497 0.99 FSRQ 1.814 6.8 0.5 2.47E+10 ... ... 2423.3
1FGL J1637.9+4707 J16377+4717 16:37:45.1 +47:17:33.822 0.91 FSRQ 0.74 1 0 6.85E+10 ... ... 815.6
1FGL J1637.9+4707 J16399+4705 16:39:56.0 +47:05:23.575 0.51 FSRQ 0.86 1 0 2.28E+09 ... ... 108.7
1FGL J1647.4+4948 J16475+4950 16:47:34.9 +49:50:00.586 0.99 RG 0.047 1.1 0.3 2.47E+10 ... ... 189
1FGL J1653.9+3945 J16538+3945 16:53:51.8 +39:45:31.490 0.99 BLL 0.034 8.3 0.6 8.10E+10 17 12 1153.2
1FGL J1656.9+6017 J16568+6012 16:56:48.2 +60:12:16.455 0.96 FSRQ 0.623 0.8 0.2 2.51E+11 ... ... 226.5
1FGL J1709.6+4320 J17096+4318 17:09:41.1 +43:18:44.547 1 FSRQ 1.027 1.5 0.3 9.36E+10 40 10 197.2
1FGL J1724.0+4002 J17240+4004 17:24:05.4 +40:04:36.457 1 RG 1.049 2.9 0.4 6.08E+10 9 7 296.1
1FGL J1727.3+4525 J17274+4530 17:27:27.6 +45:30:39.743 0.97 FSRQ 0.714 1.2 0.3 1.95E+11 ... ... 1360.3
1FGL J1727.9+5010 J17283+5013 17:28:18.6 +50:13:10.480 0.98 BLL 0.055 0.9 0.2 4.82E+09 18 22 161.8
1FGL J1734.4+3859 J17343+3857 17:34:20.6 +38:57:51.446 1 FSRQ 0.976 6 0.5 1.52E+11 ... ... 1101.3
1FGL J1740.0+5209 J17406+5211 17:40:37.0 +52:11:43.413 0.78 FSRQ 1.379 3.4 0.4 5.75E+10 ... ... 1358.5
1FGL J1742.1+5947 J17425+5945 17:42:32.0 +59:45:06.729 0.97 BLL ... 0.6 0.2 3.86E+10 ... ... 116.2
1FGL J1749.0+4323 J17490+4321 17:49:00.4 +43:21:51.287 1 BLL ... 2.3 0.3 6.17E+10 21 −10 284.9

Notes. Column 1: 1FGL source name. Column 2: VIPS source name. Column 3: right ascension (J2000). Column 4: declination (J2000). Column 5: probability that the LAT source is associated with the VIPS source. Column 6: source classification. Column 7: redshift. Column 8: γ-ray flux in units of 10−9 photons cm−2 s−1 for 1–100 GeV. Note: this is not the full 100 MeV–100 GeV flux. Column 9: error in the γ-ray flux in the same units as the flux. Column 10: core brightness temperature as measured by automated program (Helmboldt et al. 2007). Column 11: opening angle in degrees. Column 12: change in jet position angle in degrees. Column 13: CLASS flux density at 8.5 GHz in mJy.

Download table as:  ASCIITypeset images: 1 2 3

To build a γ-ray quiet sample for comparison, we excluded all LAT sources from the VIPS sample, leaving 1018 objects in our non-LAT sample. Of these 1018 objects, 24 are BL Lac objects, 479 are FSRQs, and 515 are either radio galaxies or unclassified objects.

3. GAMMA-RAY FLUX AND RADIO FLUX DENSITY

The radio flux densities in this section, and those listed in Table 1, are measured at 8.5 GHz from CLASS (Myers et al. 2003). For this section we also include seven MOJAVE-detected γ-ray loud sources for which we had CLASS flux densities, increasing the size of our LAT sample to 109 sources.

3.1. Redshift Selection Effect

Throughout the rest of this paper, we will use the nonparametric Spearman test (e.g., Press et al. 1986) to look for correlations between the LAT-detected and non-LAT-detected objects. The Spearman correlation coefficient (ρs) has a range of 0 < |ρs| < 1. A high value of ρs indicates a significant correlation. This is a powerful test for statistical correlation, but it does not test an actual physical correlation. In order to be certain our correlations are physically significant, we make sure that a redshift selection effect is not adding a bias to our data.

Imagine a population of sources covering a wide range of redshifts in which the radio and γ-ray emission is not physically correlated, but in which there is a correlation of radio flux, and gamma-ray flux, with redshift. This is actually what one would naively expect because the more distant objects will be fainter. Such a population will show a strong correlation of radio flux with gamma-ray flux via the Spearman test, but this does not indicate a significant physical correlation between these two observables.

However, by investigating the relationship between radio flux density and redshift for the LAT-detected sources (see Figure 1) we can rule out a correlation between the two. We calculated the Spearman ρs values for the S8.5z relationship for the BL Lac objects and FSRQs separately. For the BL Lac objects, the ρs value is 0.0737, with a 79% probability that random sampling would produce this same ρs value. For the FSRQs, the ρs value is 0.0450, with a 74% probability of getting the same value by random sampling. Therefore, there is no significant correlation between the radio flux density and redshift for our objects.

Figure 1.

Figure 1. CLASS 8.5 GHz radio flux density (in mJy) vs. redshift of LAT-detected sources. BL Lac objects are the top plot, FSRQs are the bottom plot.

Standard image High-resolution image

For completeness, we also tested the correlation between γ-ray flux and redshift. The ρs values are −0.319 for BL Lac objects and 0.104 for FSRQs, with the probability of getting the same values from random sampling of 21% and 45%, respectively. So, there is no significant correlation between γ-ray flux and redshift, either. See Figure 2 for a plot of γ-ray flux versus redshift.

Figure 2.

Figure 2. LAT γ-ray flux (100 MeV–100 GeV) vs. redshift. The γ-ray fluxes are in units of 10−9 photons cm−2 s−1. BL Lac objects are the top plot, FSRQs are the bottom plot.

Standard image High-resolution image

3.2. Gamma-ray Flux versus Radio Flux Density

In Figure 3, we plot the LAT flux versus the flux density at 8.5 GHz from the CLASS. The LAT fluxes are broadband fluxes from 100 MeV to 100 GeV. The BL Lac objects have a large range in their flux density and γ-ray flux, but one can see a large fraction of them have relatively low radio flux density (S8.5< 400 mJy) but high γ-ray flux. The FSRQ objects also have a very large range in both their radio flux density and their γ-ray flux, but the higher flux density FSRQs also tend to have higher γ-ray flux. Overall, there is a lack of strong γ-ray emitters among the weak radio sources, but the weak γ-ray emitters can be bright or faint at radio wavelengths. The median γ-ray flux is 3.5 × 10−8 photons cm−2 s−1 for BL Lac objects, 3.7 × 10−8 photons cm−2 s−1 for FSRQs, and 3.8 × 10−8 photons cm−2 s−1 for the radio galaxies and unclassified objects.

Figure 3.

Figure 3. LAT γ-ray flux (100 MeV–100 GeV) vs. CLASS radio flux density at 8.5 GHz. The γ-ray fluxes are in units of 10−9 photons cm−2 s−1. The black squares are BL Lac objects, the blue triangles are FSRQs, and the red circles are radio galaxies and unclassified objects.

Standard image High-resolution image

The Spearman ρs values we found were 0.217 for the BL Lac objects, 0.318 for the FSRQs, and −0.413 for the radio galaxies and unclassified objects. For the BL Lac objects and radio galaxies/unclassified objects, the chances that random samplings would result in the same ρs values were greater than 16%, indicating a lack of correlation between radio flux density and γ-ray flux. For the FSRQs, the chance of a random sampling having the same ρs was only 1.7%, which indicates that there is a significant correlation between 8.5 GHz flux density and γ-ray flux for the FSRQs. However, when we exclude the brightest 10% of sources in the radio, we find that the remaining FSRQs no longer correlate. The weaker 90% of sources have ρs's of 0.135 for BL Lac objects and 0.159 for FSRQs, with the chances of random samples having the same ρs values of 18% for BL Lac objects and 33% for FSRQs.

The lack of correlation for BL Lac objects is counter to what has been reported in other studies. Kovalev et al. (2009) find a definite correlation for their MOJAVE sources at 15 GHz. Their radio data were taken closer to the LAT measurements, and the γ-ray blazars are suspected to have higher radio flux densities during or shortly after γ-ray emission (Pushkarev et al. 2010). However, their sample is only about 1/3 the size of ours and it is biased to sources with high radio flux density (>1.5 Jy). Very strong flux–flux correlations are found by Ghirlanda et al. (2010) at 20 GHz. They have a larger sample than ours with 230 LAT sources, about half of which are FSRQs. However, their sample is also mostly comprised of bright sources. More than half of their sources have S20 > 500 mJy, compared to our sample where about 1/3 have S8.5 > 500 mJy.

To see if we could reproduce the correlation seen by others, we applied the Spearman test to the 36 sources with 8.5 GHz flux densities above 500 mJy. This subsample is roughly the brightest third of our overall sample and contains 11 BL Lac objects and 25 FSRQs. At first, we did not break the subsample into optical types. The result is a ρs of 0.390. The chance of a random sampling having the same ρs value as our strong source subsample was only 1.9%, indicating that this is a significant correlation. When the subsample was broken up by optical type, the ρs's were 0.214 for BL Lac objects and 0.422 for FSRQs. The chances of random samples having the same ρs's were 69% and 3.7%, respectively. So we find that FSRQs exhibit a significant correlation between radio flux density and γ-ray flux but BL Lac objects do not, even when limited to the strongest radio sources. Any sample that contains mostly high radio flux density objects will automatically be dominated by FSRQs and will likely show a correlation between radio flux density and γ-ray flux.

3.3. Radio Flux Density

The median flux density at 8.5 GHz is 247 mJy for the LAT BL Lac objects and 276 mJy for the non-LAT BL Lac objects. The LAT FSRQs appear to have significantly higher flux densities than the non-LAT FSRQs, with a median of 377 mJy for LAT FSRQs compared to a median of 238 mJy for non-LAT FSRQs. The LAT radio galaxies and unclassified objects also have higher 8.5 GHz flux densities than the non-LAT sources. The median flux density for LAT radio galaxies/unclassified objects is 207 mJy and the median flux density for non-LAT radio galaxies/unclassified objects is 119 mJy. The Kolomogorov–Smirnov (K-S) test probability that the LAT and non-LAT BL Lac objects belong to the same parent population is about 81%. The K-S test probability that the LAT and non-LAT FSRQs are related is only 0.3%. The radio galaxy/unclassified objects K-S test probability is 0.01%.

4. MORPHOLOGICAL COMPARISON

4.1. Source Classes

VIPS sources are classified as point sources (PSs), short jets (SJETs), long jets (LJETs), complex (CPLX), or compact symmetric object candidates (CSO) by the automatic classification scheme of Helmboldt et al. (2007). See Table 2 for the source classifications for both LAT and non-LAT sources in VIPS. From the classifications, it appears that SJETs and PS BL Lac objects are less likely to produce γ-rays. About 75% of the LAT BL Lac objects are classified as LJET, compared to only 45% for the non-LAT BL Lac objects. The major difference between LAT and non-LAT FSRQs is the lack of CSOs in the LAT sample. There is also a lack of CSOs among the LAT radio galaxies and objects without an optical classification.

Table 2. Source Classifications

LAT/Non-LAT Optical Type LJET SJET PS CPLX CSO Unidentified
LAT-detected              
  BL Lac objects 30 (75%) 4 (10%) 4 (10%) 2 (5%) ... ...
  FSRQs 22 (44%) 11 (22%) 16 (32%) 1 (2%) ... ...
  RG/Other 8 (67%) 1 (8%) 3 (25%) ... ... ...
Non-LAT-detected              
  BL Lac objects 11 (46%) 7 (29%) 6 (25%) ... ... ...
  FSRQs 188 (39%) 121 (25%) 136 (28%) 2 (∼1%) 30 (6%) 2 (∼1%) 
  RG/Other 214 (42%) 98 (19%) 111 (21%) 11 (2%) 71 (14%) 10 (2%)

Note. LJET, long jet; SJET, short jet; PS, point source; CPLX, complex; CSO, compact symmetric object candidate.

Download table as:  ASCIITypeset image

4.2. Polarization

Polarization is detected from the cores of nearly half (49 of 102) of the LAT sources in our sample. We define the core as the bright compact component at the base of the jet. We detected polarization in the core of 17 of the 40 BL Lac objects, 24 of the 50 FSRQs, and 8 of the 12 radio galaxies and unclassified objects. For the non-LAT sources, core polarization is detected in only 270 of the total 1018 sources. There are 10 of 24 non-LAT BL Lac objects, 158 of 479 FSRQs, and 102 of 515 radio galaxies or unclassified objects detected. The LAT FSRQs and radio galaxies/unclassified objects tend to have polarization in their cores more than their non-LAT counterparts. In contrast, the LAT BL Lac objects do not appear to be any more likely to have polarized cores than the non-LAT BL Lac objects. The distributions of fractional polarization in the cores for LAT and non-LAT FSRQs, and radio galaxies/unclassified objects are shown in Figure 4. The K-S test probabilities that these distributions were drawn from the same parent populations are 84.6% for the BL Lac objects, 39.4% for the FSRQs, and 76.9% for the radio galaxies/unclassified objects.

Figure 4.

Figure 4. Distributions of fractional polarization in the cores of LAT-detected FSRQs and radio galaxies/unclassified objects (top) and non-LAT FSRQs and radio galaxies/unclassified objects (bottom).

Standard image High-resolution image

Our polarization results are different than those reported by Hovatta et al. (2010) with their MOJAVE data at 15 GHz. They found that the median fractional polarization for the cores of LAT-detected sources was higher than for non-LAT sources. Their median values for their 2008–2009 observations are 2.5% for LAT-detected sources and 1.86% for non-LAT sources. Our overall median values are 3.5% for LAT sources and 4.4% for non-LAT sources. However, our result is not of high significance as the standard deviations are 2.5% for LAT sources and 3.4% for non-LAT sources. It should also be noted that our VIPS data are not contemporaneous with LAT observations. It is possible that polarization increases with γ-ray emission. In fact, the MOJAVE data indicate that the polarization levels were stronger during the LAT detections than in previous years (2002–2008), but only by a factor of 1.2.

4.3. Core Brightness Temperature

We have obtained the brightness temperatures from automatic modelfits to source components as described in Helmboldt et al. (2007). For sources with two or fewer components we have further refined our modelfitting procedure by fitting to the visibility data directly (Taylor et al. 2007). Only image plane modelfitting was carried out in Helmboldt et al. (2007) due to the tendency for the automatic visibility modelfitting to go awry for complicated sources.

There is strong evidence that the brightness temperature for the cores of the LAT-detected sources may be higher than the γ-ray quiet population (see the FSRQ and radio galaxy/unclassified distributions in Figure 5). The BL Lac objects are the exception to this trend. The median core brightness temperature for BL Lac objects detected by LAT is about 2.6 × 1010 K, while the median for remaining non-LAT BL Lac objects is 2.7 × 1010 K. The K-S test probability for these two distributions coming from the same parent population is 93.6%. The median core brightness temperature for FSRQ objects in the LAT sample is 8.2 × 1010 K, while the non-LAT FSRQ median is only about 2.5 × 1010 K. The K-S test probability that the two FSRQ distributions come from the same parent population is very low (0.001%). The median core brightness temperature for LAT radio galaxies and unclassified objects is 3.9 × 1010 K, and for non-LAT objects the median is 1.0 × 1010 K. The K-S test probability that these distributions are drawn from the same parent sample is only 4.4%.

Figure 5.

Figure 5. Distributions of core brightness temperatures for LAT-detected FSRQs and radio galaxies/unclassified objects (top) compared to the distributions of core brightness temperature for non-LAT FSRQs and radio galaxies/unclassified objects(bottom).

Standard image High-resolution image

Interestingly, γ-ray flux does not appear to have any relation to core brightness temperature. We applied the Spearman test to the γ-ray versus core brightness temperature distributions and found correlation coefficients of 0.093 for BL Lac objects, 0.111 for FSRQs, and −0.181 for radio galaxies and unclassified objects.

4.4. Jet Opening Angle

We have measured a mean opening half-angle by the following procedure (Taylor et al. 2007): We measure the separation of each jet component from the core along the jet axis (taken to be a linear fit to the component positions) and the distance of each component from the jet axis, i.e., x' and y' positions in a rotated coordinate system with the jet axis along the x'-axis. For each component, we measure its extent from its center along a line perpendicular to the jet axis using the parameters of its elliptical fit, and then deconvolve this using the extent of the Gaussian restoring beam along the same line. The opening half-angle measured from each component is then taken to be

where dr is the deconvolved Gaussian size perpendicular to the jet axis. After measuring this for each jet component, we average them to get a single value. This is only done for sources with more than two total components, (i.e., at least two jet components). Unfortunately, we can only make quantitative estimates for 30 of the LAT sources; 16 BL Lac, 9 FSRQ, and 5 radio galaxy or unclassified objects. The distributions of opening angles are shown in Figure 6. Although this is based on very small statistics there is evidence that the LAT source jets have unusually large opening angles. This is especially true for the FSRQ detected by LAT. The MOJAVE (Pushkarev et al. 2009) and TANAMI (Ojha et al. 2010) groups both reported larger opening angles for γ-ray loud blazars. Taylor et al. (2007) also found evidence for larger opening angles in EGRET sources.

Figure 6.

Figure 6. Distribution of opening angles for LAT-detected BL Lac objects, FSRQs, and radio galaxies/unclassified objects (top) compared to the distribution of opening angles for non-LAT objects (bottom).

Standard image High-resolution image

We found 10 LAT-detected objects with opening angles greater than 30°. The large opening angle objects are BL Lac objects and FSRQs, five of each. See Figure 7 for contour maps of these 10 sources.

Figure 7.

Figure 7. Contour maps of the 10 LAT sources with opening angles larger than 30°. All data were taken with the VLBA at 5 GHz. Bottom contours are at 0.5 mJy beam−1, except for the first object which has a bottom contour level is 1 mJy beam−1. The restoring beam is shown as an ellipse in the lower left corner of each panel.

Standard image High-resolution image

4.5. Core Fraction

We define the core fraction as the ratio of flux density in the core compared to the total flux density at 5 GHz. A high core fraction can be an indication of a high Doppler factor. Taylor et al. (2007) reported that EGRET sources tended to have greater core fractions than sources not detected by EGRET. With the LAT sources, we find that this is not the case. The median core fractions are 0.868 for LAT BL Lac objects and 0.896 for non-LAT BL Lac objects, 0.939 for LAT FSRQs and 0.929 for non-LAT FSRQs, and 0.889 for LAT radio galaxies/unclassified objects and 0.886 for non-LAT radio galaxies/unclassified objects. The K-S test probabilities that the distributions were drawn from the same parent distribution are 26.4% for BL Lac objects, 16.1% for FSRQs, and 23.5% for radio galaxies/unclassified objects, respectively.

4.6. Jet Bending, Length, and Brightness Temperature

We can also compare the amount of bending in the jet on the parsec scale for the LAT sources compared to non-LAT sources. Unfortunately, as with the jet opening angle measurements, we are limited to only 30 LAT sources and the results are inconclusive.

We compared the distribution of jet length in LAT sources with non-LAT sources. There does not appear to be an appreciable difference between the two populations of FSRQ an radio galaxies/unclassified objects. The K-S test results are 63.3% for the FSRQ objects, and 64.6% for the radio galaxy and unclassified objects. Interestingly, the BL Lac objects in the LAT sample appear to have longer jets than the non-LAT BL Lac objects. The distributions for the BL Lac objects are shown in Figure 8. The K-S test probability that the two are drawn from the same distribution is 2.1%. The FSRQ and AGNs/other objects have essentially the same distributions.

Figure 8.

Figure 8. Distribution of jet lengths for LAT-detected BL Lac object (top) compared to the distribution of jet lengths for non-LAT BL Lac objects (bottom).

Standard image High-resolution image

Unlike the core brightness temperatures, the jet brightness temperatures (formally, the brightness temperature of the brightest jet component) for the LAT and non-LAT distributions look fairly similar. The K-S test probabilities for the BL Lac object, FSRQ, and radio galaxy/other pairs of distributions are 34.2%, 92.1%, and 83.5%, respectively.

5. NOTES ON INDIVIDUAL SOURCES

5.1. Low Core Brightness Temperature

J09576+5522. With a core brightness temperature of only 9.7 × 108 K, this FSRQ has a γ-ray flux of (106.5 ± 7.4) × 10−9 photons cm−2 s−1, nearly three times the median for FSRQs. Also known as [HB89] 0954+556, it has a redshift of 0.896 (Abazajian et al. 2004) and a 8.5 GHz flux density is 1.5 Jy. This FSRQ has been studied in the UV (Pian et al. 2005), radio, X-ray (Tavecchio et al. 2007), and was detected in the γ-rays by EGRET (Hartman et al. 1999). Tavecchio et al. (2007) noted that X-ray flux was only found at the terminal portion of the jet. The automated program initially reported the core brightness temperature for this object to be 1.2 × 108 K. We suspected program did not generate an accurate model for this source. When we did the model fitting by hand in DIFMAP, we calculated a core brightness temperature of 9.7 × 108 K. This is the lowest core temperature of all objects in our LAT sample. This object is featured in Figure 7 because its jet opening angle is about 60°.

J08163+5739. The automated program found that this BL Lac object had a core brightness temperature of only 4.16 × 107 K. Again, when we did the model fitting by hand we calculated a core brightness temperature of 1.6 × 109 K. Even with this higher core brightness temperature, this object is still the coolest BL Lac object in our LAT sample. This object has a γ-ray flux of (15.2 ± 0.3) × 10−9 photons cm−2 s−1. It is a compact and weak radio source with a radio flux density of 92 mJy at 8.5 GHz. There is no measured redshift for this object. It is also known as SBS 0812+578 and 1RXS J081624.6+573910. It has a small jet extending to the east. This object is featured in Figure 7 because its jet opening angle is approximately 52°.

5.2. High Redshift

J07464+2549. With a redshift of 2.979 (Abazajian et al. 2004), this FSRQ has a γ-ray flux (42.3 ± 7.9) × 10−9 photons cm−2 s−1. Its core brightness temperature is 5.07 × 1011 K, which is well above the median for the LAT-detected FSRQs. It is also known as B2 0743+25 and RX J0746.4+2549. It is a relatively compact and bright radio source with a flux density of 731 mJy at 8.5 GHz (see Figure 9). It has a very small jet extending to the north. In 2005, the Burst Alert Telescope on board the Swift satellite detected a flare of 15–195 keV emission from this object (Sambruna et al. 2006).

Figure 9.

Figure 9. Contour maps of the objects in Sections 5.2 and 5.3. All data were taken with the VLBA at 5 GHz. Bottom contours are at 0.5 mJy beam−1. The restoring beam is shown as an ellipse in the lower left corner of each panel.

Standard image High-resolution image

J08053+6144. With a redshift of 3.033 (Sowards-Emmerd et al. 2005), this FSRQ has the highest redshift among the LAT-detected objects. It has a γ-ray flux 36.0 ± 8.9 × 10−9 photons cm−2 s−1. Its core brightness temperature is 1.30 × 1011 K, which is slightly higher than the median for LAT-detected FSRQs. This object is fairly bright in the radio with a flux density of 725 mJy at 8.5 GHz. It has a short jet extending to the southeast (see Figure 9).

5.3. Radio Galaxies

Our sample contains six sources identified as radio galaxies, although in some cases there is controversy in this identification as described below.

J09235+4125. This object is also known as B3 0920+416 and GB6 J0923+4125. It has a γ-ray flux of (40.6 ± 10.0) × 10−9 photons cm−2 s−1 and a redshift of 0.028 (Falco et al. 1998). Its peak radio flux density is 120 mJy at 8.5 GHz. It has a short jet extending to the east (see Figure 9). The NASA/IPAC Extragalactic Database (NED) classifies this object as an FSRQ. However, this is probably an error. Falco et al. (1998) call it a "late-type galaxy."

J12030+6031. This object is also known as SBS 1200+608 and RX J1203.0+6031. It has a γ-ray flux of (44.8 ± 0.3) × 10−9 photons cm−2 s−1 and a redshift of 0.065 (Falco et al. 1998). Its peak radio flux density is 105 mJy at 8.5 GHz. It has a short jet extending to the south (see Figure 9). Véron-Cetty & Véron (2006) classify this object as a low ionization nuclear emission region (LINER) galaxy, sometimes referred to as a Seyfert 3.

J13307+5202. This radio galaxy is also known as CGRaBS J1330+5202. It has an upper limit to its γ-ray flux of 41.5 × 10−9 photons cm−2 s−1 and a redshift of 0.688 (Healey et al. 2008). It has a peak radio flux density of 127 mJy at 8.5 GHz and a jet extending to the southwest (see Figure 9). Sowards-Emmerd et al. (2005) list it as an unknown optical type.

J16071+1551. This object is also known as [HB89] 1604+159. It is an EGRET-detected source known as 3EG J1605+1553. It has a LAT γ-ray flux of (33.9 ± 7.7) × 10−9 photons cm−2 s−1 and a redshift of 0.496 (Adelman-McCarthy et al. 2008). It has a peak radio flux density of 214 mJy at 8.5 GHz and a jet extending to the east (see Figure 9). Véron-Cetty & Véron (2006) classify this object as a BL Lac object.

J16475+4950. This object is also known as SBS 1646+499 and 1RXS J164735.4+495001. It has a γ-ray flux of 33.3 ± 8.5 × 10−9 photons cm−2 s−1 and a redshift of 0.047 (Falco et al. 1998). It has a peak radio flux density of 131 mJy at 8.5 GHz. It has a very short jet extending to the southeast (see Figure 9). The spectrum in Marcha et al. (1996) indicates that this object is a Seyfert 1 blazar.

J17240+4004. This object is also known as B2 1722+40 and RX J1724.0+4004. It has a γ-ray flux of (47.1 ± 9.0) × 10−9 photons cm−2 s−1 and a redshift of 1.049 (Vermeulen et al. 1996). This is the brightest radio galaxy in our sample with a peak flux density of 323 mJy at 8.5 GHz. It has a long jet extending to the northwest (see Figure 9). Véron-Cetty & Véron (2006) tentatively classify this object as a BL Lac object. However, the spectrum in Vermeulen et al. (1996) does not support this classification. The spectrum has only narrow lines, so it cannot be an FSRQ, and the equivalent widths are greater than 5 Å, so it cannot be a BL Lac object.

6. DISCUSSION

One of the primary explanations for some AGNs being γ-ray loud is Doppler Boosting. The orientation and material speed of the AGN jets have a large effect on their observed fluxes due to the Doppler effect. The kinematic Doppler factor for a moving source is defined as

where β is the bulk velocity v/c, Γ is the Lorentz factor, Γ = (1 − β2)−1/2, and θ is the angle between the velocity vector and the line of sight. Emissions from an approaching source are blueshifted by

where ν' is the frequency in the source's rest frame. The intensity enhancement, commonly known as the "Doppler boosting," is given in Rybicki & Lightman (1979) as

Because LAT measures the flux across a large bandwidth (1–100 GeV), we have to integrate over frequency:

Therefore, a small difference in material velocity or orientation angle leads to a large difference in intensity. See Urry & Padovani (1995) for a more thorough discussion of Doppler enhancement.

6.1. BL Lac Objects

Almost two-thirds of the BL Lac objects in VIPS are detected by LAT, and the only significant difference between the LAT and non-LAT populations is in their jet lengths. We have shown that the LAT BL Lac objects have longer jets than the non-LAT BL Lac objects (see Section 4.6). There are no differences in the fraction of polarized BL Lac objects or the distributions of the polarization. The LAT and non-LAT BL Lac objects have nearly identical core brightness temperature distributions. Why are the non-LAT BL Lac objects quiet in the γ-ray? It is possible that the non-LAT BL Lac objects are fundamentally different from the LAT BL Lac objects, but it must be in some subtle way. It seems more likely that all the BL Lac objects are emitting in the γ-ray, but some are emitting below the threshold for the LAT to detect.

Initially, we suspected that the non-LAT BL Lac objects might have been too far away to detect their γ-rays. The more distant BL Lac objects definitely have lower γ-ray flux, but it is clear that the non-LAT BL Lac objects are not at higher redshifts than the LAT BL Lac objects. While we only have redshift measurements for 16 of the LAT BL Lac objects and 15 of the non-LAT BL Lac objects, that is still enough to rule out redshift as the cause of LAT non-detection. The maximum measured redshift for the LAT BL Lac objects is 0.56. Nearly half (11/24) of the non-LAT BL Lac objects have redshifts lower than 0.56. Also, there is no strong correlation between redshift and γ-ray flux for the BL Lac objects (see Section 3.1).

The prevalence of LJET BL Lac objects is an interesting clue to this puzzle. The K-S test result of 2.1% for the jet length distributions (see Section 4.6) is only of modest statistical significance, but the fact that a larger percentage of the LAT BL Lac objects are classified as LJET than the non-LAT BL Lac objects (see Table 2) supports it. One would expect the jets to appear shorter if they are pointed more directly at us. However, if the γ-ray loud BL Lac objects have higher Doppler factors than their γ-ray quiet counterparts, that would indicate that the LAT BL Lac objects have more powerful jets than the non-LAT BL Lac objects. A more powerful jet will appear longer even when viewed at a small angle.

Because we find that the LAT BL Lac objects have longer jets than non-LAT BL Lac objects, the orientation angle for LAT BL Lac objects cannot be smaller than non-LAT BL Lac objects. The non-LAT BL Lac objects could have slightly lower jet speeds. This can explain why the LAT BL Lac objects have longer jets and why they are γ-ray loud. Unfortunately, this theory is not supported by recent work by Lister et al. (2009b), where they find that two non-LAT BL Lac objects have higher maximum jet material velocities than LAT BL Lac objects. However, they only had 10 BL Lac objects in their LAT-detected sample, so this is still an open question.

BL Lac objects are known to have high variability in their flux densities. The non-LAT BL Lac objects could currently be in low γ-ray luminosity states, but may enter states of enhanced γ-ray activity at a later time. It is possible that the jet material speed is related to this flux variability. Based on the number of BL Lac objects in our sample, the duty cycle would need to be about 2/3 over the course of one year.

6.2. FSRQs

The LAT FSRQs seem to be substantially different from the non-LAT FSRQs. Only about 9% (50/529) of the FSRQs in VIPS are γ-ray loud. This indicates that the LAT FSRQs are special in some way. Recall from Section 3.2 that γ-ray flux correlates with radio flux density for FSRQs, especially for high flux density sources. Also, in general the LAT FSRQs have significantly higher core brightness temperatures, with the notable exception of J09576+5522, and they appear to have larger opening angles. So, brighter FSRQs are more likely to produce γ-rays. We showed in Section 3.1 that γ-ray flux does not correlate with redshift, so we do not expect the LAT FSRQs to have lower redshift than their non-LAT counterparts. We explored the possibility that LAT FSRQs were at low redshift, but found that they cover a range similar to the non-LAT FSRQs.

The γ-ray flux from FSRQs can be explained with Doppler boosting, but requires a more dramatic difference than with the BL Lac objects. Lister et al. (2009b) reported that LAT-detected FSRQs had higher median jet speeds than LAT-detected BL Lac objects by more than a factor of two. It is also possible that there are two subclasses for FSRQs: γ-ray loud and γ-ray quiet. The difference in γ-ray detection rates for FSRQs could also be due to variations in the γ-ray emission, however, the duty cycle required (about 10%) is much more extreme than for the BL Lac objects.

A larger fraction of the LAT FSRQs have polarization in their cores compared to the non-LAT FSRQs. We detected polarization in the cores of 48% of the LAT FSRQs compared to 33% of the non-LAT FSRQs, a difference that is significant at about the 1σ level. This points to an intrinsic difference in γ-ray quiet versus γ-ray loud FSRQs. The higher polarization could indicate more organized magnetic fields, possibly pointing to a greater spin of the supermassive black hole (Meier 2001). A higher black hole spin would also lead to higher material velocities in the jets (Meier 1999).

6.3. Radio Galaxies

There is a definite lack of CSOs among the LAT-detected objects. This is not surprising from a Doppler boosting point of view. For an object to appear symmetric, we need to see it close to "edge on" so we can see both jets. With the jets pointed nearly perpendicular to our line of sight, we would not expect to see much Doppler boosting in the radiation. However, Stawarz et al. (2008) suggested that CSOs should account for a "relatively numerous class of extragalactic sources" for Fermi due to ultrarelativistic electrons in the lobes of these compact objects. We have yet to detect any γ-rays from a CSO.

The LAT radio galaxies and unclassified objects are a puzzle. They do not appear to differ significantly from their non-LAT counterparts in any appreciable way. Their opening angles appear to be very similar to the non-LAT objects and they are not particularly bright in the radio. The LAT radio galaxies/unclassified objects may have higher core temperatures than the non-LAT objects. The one difference that stands out is that we detected polarization in the cores of 67% of the LAT radio galaxies/unclassified objects compared to only 20% for their non-LAT counterparts. Unfortunately, our current sample is too small to make any definite claims about these objects. There is also controversy about the optical types of several of the objects we label as radio galaxies.

6.4. Constraints on the γ-ray Emitting Region

The major differences between LAT and non-LAT sources are confined to the cores of the sources. The LAT BL Lac objects have stronger jets, suggesting higher ejection velocities from the cores. The LAT FSRQs have higher core temperatures and more frequently exhibit polarized cores. The LAT radio galaxies/unclassified objects are also more polarized than their non-LAT counterparts. We suspect that the γ-rays are being produced very near to the cores of these objects. To put a limit on the size of the γ-ray emitting region, we measured the sizes of the cores for several objects. Using all of the objects for which we have measured redshifts, the median upper bound on the core size is 4.54 pc with a standard deviation of 10.98 pc. However, a better estimate for the upper limit on the size of the emitting region should come from measurements of the objects nearest to us. Limiting ourselves to objects with z ⩽ 0.1, we have four BL Lac objects and three radio galaxies (a LINER, a Seyfert 1, and an object that may actually be an FSRQ). Using this subset of sources, we find a median upper bound on the core size of 0.58 pc with a standard deviation of 0.34 pc. Therefore, the γ-ray emission region must be within about 0.9 pc of the central engine.

7. CONCLUSIONS

We compared the distributions of core polarization, core brightness temperature, jet opening angles, core fractions, jet bending, jet length, and 8.5 GHz radio flux densities for LAT and non-LAT sources. The only significant difference between the LAT and non-LAT BL Lac objects is that the γ-ray loud BL Lac objects appear to have longer jets. Also, we find no correlation between BL Lac radio flux density and γ-ray flux. We suspect that all BL Lac objects are producing γ-rays, but some are below LAT's threshold due to jet material speed, or possibly variability in the γ-ray luminosity. In contrast, γ-ray bright FSRQs are considerably different from the γ-ray quiet population. The LAT FSRQs appear to be extreme sources with higher core brightness temperatures, larger opening angles, and greater core polarization than their non-LAT counterparts. There is also a significant correlation between the FSRQ radio flux density and γ-ray flux, especially for brighter sources. The γ-ray loud FSRQs can be explained by Doppler boosting, but the orientation and/or jet material speeds must be significantly different than the γ-ray quiet FSRQs. One possible way to explain this difference is that the γ-ray loud FSRQs may have greater angular momentum in their central supermassive black holes, leading to higher jet velocities and greater polarization. Due to the small number of radio galaxies in our sample, combined with the controversy over some of their optical types, there is little we can say about γ-ray bright radio galaxies. However, we do find two intriguing clues about them. First, there are no γ-ray loud CSOs in our sample. Second, the radio galaxies that we do detect in the γ-rays tend to have strong polarizations in their cores.

Because the differences we find between the γ-ray loud and γ-ray quiet objects are related to the core, we suspect that the γ-ray radiation originates within the core (i.e., at the base of the jet). By calculating the core sizes of the LAT objects with measured redshifts, we put an upper limit of 15 pc on the size of this emitting region. A better limit of 0.9 pc results from using the core sizes of the lowest redshift objects in our LAT sample.

As more data are collected more will become clear about γ-ray loud blazars. We, and several other groups, are currently conducting further observations. We will soon have more information about the velocities of the jet material and more associations of AGNs with γ-ray sources.

We thank Matt Lister and Ann Wehrle for useful discussions. We also thank Teddy Cheung, Marcello Giroletti, Yuri Kovalev, Frank Schinzel, David Thompson, and Steve Tremblay for their helpful comments. Also, we express our appreciation to the anonymous referee for their insightful and helpful criticism. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. We thank NASA for support under FERMI grant GSFC 21078/FERMI08-0051.

Please wait… references are loading.
10.1088/0004-637X/726/1/16