• Combined with very-long-baseline interferometry measurements, the observations of fast TeV gamma-ray flares probe the structure and emission mechanism of blazar jets. However, only a handful of such flares have been detected to date, and only within the last few years have these flares been observed from lower-frequency-peaked BL~Lac objects and flat-spectrum radio quasars. We report on a fast TeV gamma-ray flare from the blazar BL~Lacertae observed by VERITAS, with a rise time of $\sim$2.3~hr and a decay time of $\sim$36~min. The peak flux above 200 GeV is $(4.2 \pm 0.6) \times 10^{-6} \;\text{photon} \;\text{m}^{-2}\; \text{s}^{-1}$ measured with a 4-minute-binned light curve, corresponding to $\sim$180\% of the flux which is observed from the Crab Nebula above the same energy threshold. Variability contemporaneous with the TeV gamma-ray flare was observed in GeV gamma-ray, X-ray, and optical flux, as well as in optical and radio polarization. Additionally, a possible moving emission feature with superluminal apparent velocity was identified in VLBA observations at 43 GHz, potentially passing the radio core of the jet around the time of the gamma-ray flare. We discuss the constraints on the size, Lorentz factor, and location of the emitting region of the flare, and the interpretations with several theoretical models which invoke relativistic plasma passing stationary shocks.
  • We aim to better understand the dynamics within relativistic magneto-hydrodynamical flows in the extreme environment and close vicinity of supermassive black holes. To do so, we analyze the peculiar radio galaxy 3C 111, for which long-term polarimetric observations are available. We make use of the high spatial resolution of the VLBA network and the MOJAVE monitoring program, which provides high data quality also for single sources and allows us to study jet dynamics on parsec scales in full polarization with an evenly sampled time-domain. We additionally consider data from the IRAM 30-m Telescope as well as the SMA. Jet properties such as the electric vectors, the (polarized) flux density, feature size, and brightness temperature, describe a complex evolution of the polarized jet. The electric vector position angles (EVPAs) of features traveling down the jet perform a large and smooth rotation of $\gtrsim 180^{\circ}$ across a distance of about 20 pc. In contrast, the EVPAs are strongly variable within the first parsecs of the jet. We find a tendency towards transverse EVPAs across the jet with a local anomaly of aligned vectors in between. The transverse extent of the flow decreases coincident with a jump in brightness temperature around where we observe the EVPAs to turn into alignment with the jet flow. Also the gradients of the feature size and particle density with distance steepen in that region. We interpret the propagating polarized features with shocks and the observed local anomalies with the interaction of these shocks with a recollimation shock of the underlying flow. Together with a sheared magnetic field, this shock-shock interaction can explain the large rotation of the EVPA. The superimposed variability of the EVPAs close to the core is likely related to a clumpy Faraday screen, which also contributes significantly to the observed EVPA rotation in that region.
  • Observations of fast TeV $\gamma$-ray flares from blazars reveal the extreme compactness of emitting regions in blazar jets. Combined with very-long-baseline radio interferometry measurements, they probe the structure and emission mechanism of the jet. We report on a fast TeV $\gamma$-ray flare from BL Lacertae observed by VERITAS, with a rise time of about 2.3 hours and a decay time of about 36 minutes. The peak flux at $>$200 GeV measured with the 4-minute binned light curve is $(4.2 \pm 0.6) \times 10^{-6} \;\text{photons} \;\text{m}^{-2}\, \text{s}^{-1}$, or $\sim$180% the Crab Nebula flux. Variability in GeV $\gamma$-ray, X-ray, and optical flux, as well as in optical and radio polarization was observed around the time of the TeV $\gamma$-ray flare. A possible superluminal knot was identified in the VLBA observations at 43 GHz. The flare constrains the size of the emitting region, and is consistent with several theoretical models with stationary shocks.
  • We have measured the angular sizes of radio cores of active galactic nuclei (AGN) and analyzed their sky distributions and frequency dependencies to study synchrotron opacity in AGN jets and the strength of angular broadening in the interstellar medium. We have used archival very long baseline interferometry (VLBI) data of more than 3000 compact extragalactic radio sources observed at frequencies, $\nu$, from 2 to 43 GHz to measure the observed angular size of VLBI cores. We have found a significant increase in the angular sizes of the extragalactic sources seen through the Galactic plane ($|b|\lesssim10^\circ$) at 2, 5 and 8 GHz, about 1/3 of which show significant scattering. These sources are mainly detected in directions to the Galactic bar, the Cygnus region, and a region with galactic longitudes $220^\circ\lesssim l\lesssim260^\circ$ (the Fitzgerald window). The strength of interstellar scattering of the AGNs is found to correlate with the Galactic H$\alpha$ intensity, free-electron density, and Galactic rotation measure. The dependence of scattering strengths on source redshift is insignificant, suggesting that the dominant scattering screens are located in our Galaxy. The observed angular size of Sgr A$^\ast$ is found to be the largest among thousands of AGN observed over the sky; we discuss possible reasons of this strange result. Excluding extragalactic radio sources with significant scattering, we find that angular size of opaque cores in AGN scales typically as $\nu^{-1}$ confirming predictions of a conical synchrotron jet model with equipartition.
  • We present 15 GHz stacked VLBA images of 373 jets associated with active galactic nuclei (AGN) having at least five observing epochs within a 20 yr time interval 1994-2015 from the MOJAVE programme and/or its precursor, the 2 cm VLBA Survey. These data are supplemented by 1.4 GHz single-epoch VLBA observations of 135 MOJAVE AGNs to probe larger scale jet structures. The typical jet geometry is found to be close to conical on scales from hundreds to thousands of parsecs, while a number of galaxies show quasi-parabolic streamlines on smaller scales. A true jet geometry in a considerable fraction of AGNs appears only after stacking epochs over several years. The jets with significant radial accelerated motion undergo more active collimation. We have analysed total intensity jet profiles transverse to the local jet ridgeline and derived both apparent and intrinsic opening angles of the flows, with medians of $21.5\deg$ and $1.3\deg$, respectively. The Fermi LAT-detected gamma-ray AGNs in our sample have, on average, wider apparent and narrower intrinsic opening angle, and smaller viewing angle than non LAT-detected AGNs. We have established a highly significant correlation between the apparent opening angle and gamma-ray luminosity, driven by Doppler beaming and projection effects.
  • J. Aleksić, S. Ansoldi, L. A. Antonelli, P. Antoranz, C. Arcaro, A. Babic, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, E. Carmona, A. Carosi, P. Colin, E. Colombo, J. L. Contreras, J. Cortina, S. Covino, P. Da Vela, F. Dazzi, A. De Angelis, G. De Caneva, B. De Lotto, E. de O. na Wilhelmi, C. Delgado Mendez, F. Di Pierro, D. Dominis Prester, D. Dorner, M. Doro, S. Einecke, D. Eisenacher, D. Elsaesser, A. Fernández-Barral, D. Fidalgo, M. V. Fonseca, L. Font, K. Frantzen, C. Fruck, D. Galindo, R. J. García López, M. Garczarczyk, D. Garrido Terrats, M. Gaug, N. Godinović, A. González Mu. noz, S. R. Gozzini, D. Hadasch, Y. Hanabata, M. Hayashida, J. Herrera, J. Hose, D. Hrupec, W. Idec, V. Kadenius, H. Kellermann, M. L. Knoetig, K. Kodani, Y. Konno, J. Krause, H. Kubo, J. Kushida, A. La Barbera, D. Lelas, N. Lewandowska, E. Lindfors, S. Lombardi, F. Longo, M. López, R. López-Coto, A. López-Oramas, E. Lorenz, I. Lozano, M. Makariev, K. Mallot, G. Maneva, K. Mannheim, L. Maraschi, B. Marcote, M. Mariotti, M. Martínez, D. Mazin, U. Menzel, J. M. Miranda, R. Mirzoyan, A. Moralejo, P. Munar-Adrover, D. Nakajima, V. Neustroev, A. Niedzwiecki, M. Nievas Rosillo, K. Nilsson, K. Nishijima, K. Noda, R. Orito, A. Overkemping, S. Paiano, M. Palatiello, D. Paneque, R. Paoletti, J. M. Paredes, X. Paredes-Fortuny, M. Persic, J. Poutanen, P. G. Prada Moroni, E. Prandini, I. Puljak, R. Reinthal, W. Rhode, M. Ribó, J. Rico, J. Rodriguez Garcia, T. Saito, K. Saito, K. Satalecka, V. Scalzotto, V. Scapin, T. Schweizer, S. N. Shore, A. Sillanpää, J. Sitarek, I. Snidaric, D. Sobczynska, A. Stamerra, T. Steinbring, M. Strzys, L. Takalo, H. Takami, F. Tavecchio, P. Temnikov, T. Terzić, D. Tescaro, M. Teshima, J. Thaele, D. F. Torres, T. Toyama, A. Treves, P. Vogler, M. Will, R. Zanin, S. Buson, F. D'Ammando, A. Lähteenmäki, T. Hovatta, Y. Y. Kovalev, M. L. Lister, W. Max-Moerbeck, C. Mundell, A. B. Pushkarev, E. Rastorgueva-Foi, A. C. S. Readhead, J. L. Richards, J. Tammi, D. A. Sanchez, M. Tornikoski, T. Savolainen, I. Steele
    June 6, 2016 astro-ph.HE
    1ES 1011+496 $(z=0.212)$ was discovered in very high energy (VHE, E >100 GeV) $\gamma$-rays with MAGIC in 2007. The absence of simultaneous data at lower energies led to a rather incomplete characterization of the broadband spectral energy distribution (SED). We study the source properties and the emission mechanisms, probing whether a simple one-zone synchrotron-self-Compton (SSC) scenario is able to explain the observed broadband spectrum. We analyzed VHE to radio data from 2011 and 2012 collected by MAGIC, $Fermi$-LAT, $Swift$, KVA, OVRO, and Mets\"ahovi in addition to optical polarimetry data and radio maps from the Liverpool Telescope and MOJAVE. The VHE spectrum was fit with a simple power law with a photon index of $3.69\pm0.22$ and a flux above 150 GeV of $(1.46\pm0.16)\times10^{-11}$ ph cm$^{-2}$ s$^{-1}$. 1ES 1011+496 was found to be in a generally quiescent state at all observed wavelengths, showing only moderate variability from radio to X-rays. A low degree of polarization of less than 10% was measured in optical, while some bright features polarized up to 60% were observed in the radio jet. A similar trend in the rotation of the electric vector position angle was found in optical and radio. The radio maps indicated a superluminal motion of $1.8\pm0.4\,c$, which is the highest speed statistically significantly measured so far in a high-frequency-peaked BL Lac. For the first time, the high-energy bump in the broadband SED of 1ES 1011+496 could be fully characterized from 0.1 GeV to 1 TeV which permitted a more reliable interpretation within the one-zone SSC scenario. The polarimetry data suggest that at least part of the optical emission has its origin in some of the bright radio features, while the low polarization in optical might be due to the contribution of parts of the radio jet with different orientations of the magnetic field to the optical emission.
  • We report on the acceleration properties of 329 features in 95 blazar jets from the MOJAVE VLBA program. Nearly half the features and three-quarters of the jets show significant changes in speed and/or direction. In general, apparent speed changes are distinctly larger than changes in direction, indicating that changes in the Lorentz factors of jet features dominate the observed speed changes rather than bends along the line of sight. Observed accelerations tend to increase the speed of features near the jet base, $\lesssim 10-20$ parsecs projected, and decrease their speed at longer distances. The range of apparent speeds at fixed distance in an individual jet can span a factor of a few, indicating that shock properties and geometry may influence the apparent motions; however, we suggest that the broad trend of jet features increasing their speed near the origin is due to an overall acceleration of the jet flow out to de-projected distances of order $10^2$ parsecs, beyond which the flow begins to decelerate or remains nearly constant in speed. We estimate intrinsic rates of change of the Lorentz factors in the galaxy frame of order $\dot{\Gamma}/\Gamma \simeq 10^{-3}$ to $10^{-2}$ per year which can lead to total Lorentz factor changes of a factor of a few on the length scales observed here. Finally, we also find evidence for jet collimation at projected distances of $\lesssim 10$ parsecs in the form of the non-radial motion and bending accelerations that tend to better align features with the inner jet.
  • Parsec-scale VLBA images of BL Lac at 15 GHz show that the jet contains a permanent quasi-stationary emission feature 0.26 mas (0.34 pc projected) from the core, along with numerous moving features. In projection, the tracks of the moving features cluster around an axis at position angle -166.6 deg that connects the core with the standing feature. The moving features appear to emanate from the standing feature in a manner strikingly similar to the results of numerical 2-D relativistic magneto-hydrodynamic (RMHD) simulations in which moving shocks are generated at a recollimation shock. Because of this, and the close analogy to the jet feature HST-1 in M87, we identify the standing feature in BL Lac as a recollimation shock. We assume that the magnetic field dominates the dynamics in the jet, and that the field is predominantly toroidal. From this we suggest that the moving features are compressions established by slow and fast mode magneto-acoustic MHD waves. We illustrate the situation with a simple model in which the slowest moving feature is a slow-mode wave, and the fastest feature is a fast-mode wave. In the model the beam has Lorentz factor about 3.5 in the frame of the host galaxy, and the fast mode wave has Lorentz factor about 1.6 in the frame of the beam. This gives a maximum apparent speed for the moving features 10c. In this model the Lorentz factor of the pattern in the galaxy frame is approximately 3 times larger than that of the beam itself.
  • We have obtained milliarcsecond-scale spectral index distributions for a sample of 190 extragalactic radio jets through the MOJAVE project. The sources were observed in 2006 at 8.1, 8.4, 12.1, and 15.4 GHz, and we study the four-frequency spectrum in individual jet features. We performed detailed simulations to study the effects of image alignment and (u,v)-plane coverage on the spectral index maps to verify our results. The core spectral indices are on average flat, while the jet spectrum is in general steep with a mean index of -1.04. A simple power-law fit is often inadequate for the core regions, as expected if the cores are partially self-absorbed. The overall jet spectrum steepens at a rate of about -0.001 to -0.004 per deprojected parsec when moving further out from the core with flat spectrum radio quasars having significantly steeper spectra than the BL Lac objects. The spectrum in both types of objects flattens at the locations of the jet components indicating particle acceleration or density enhancements along the jet. The mean spectral index at the component locations of -0.81 corresponds to a power-law index of ~2.6 for the electron energy distribution. We find a significant trend that jet components with linear polarization parallel to the jet (magnetic field perpendicular to the jet) have flatter spectra, as expected for transverse shocks. Compared to quasars, BL Lacs have more jet components with perpendicular magnetic field alignment, which may explain their generally flatter spectra. The overall steepening of the spectra with distance can be explained with radiative losses if the jets are collimating or with the evolution of the high-energy cutoff in the electron spectrum if the jets are conical. This interpretation is supported by a significant correlation with the age of the component and the spectral index, with older components having steeper spectra. (abridged)
  • Previously unremarkable, the extragalactic radio source GB 1310+487 showed a gamma-ray flare on 2009 November 18, reaching a daily flux of ~10^-6 photons/cm^2/s at energies E>100 MeV and becoming one of the brightest GeV sources for about two weeks. Its optical spectrum is not typical for a blazar, instead, it resembles those of narrow emission-line galaxies. We investigate changes of the object's radio-to-GeV spectral energy distribution (SED) during and after the prominent GeV flare with the aim to determine the nature of the object and constrain the origin of the variable high-energy emission. The data collected by the Fermi and AGILE satellites at gamma-ray energies, Swift at X-ray and ultraviolet, Kanata, NOT, and Keck telescopes at optical, OAGH and WISE at infrared, and IRAM 30m, OVRO 40m, Effelsberg 100m, RATAN-600, and VLBA at radio, are analysed together to trace the SED evolution on timescales of months. The gamma-ray/radio-loud narrow-line active galactic nucleus (AGN) is located at redshift z=0.638. It is shining through an unrelated foreground galaxy at z=0.500. The AGN light is likely amplified by a factor of a few because of gravitational lensing. The AGN SED shows a two-humped structure typical of blazars and gamma-ray-loud NLSy1 galaxies, with the high-energy (inverse-Compton) emission dominating by more than an order of magnitude over the low-energy (synchrotron) emission during gamma-ray flares. The difference between the two SED humps is smaller during the low-activity state. Fermi observations reveal a strong correlation between the gamma-ray flux and spectral index, with the hardest spectrum observed during the brightest gamma-ray state. If the gamma-ray flux is a mixture of synchrotron self-Compton (SSC) and external Compton (EC) emission, the observed GeV spectral variability may result from varying relative contributions of these two emission components.
  • We report that active galactic nucleus (AGN) jets are causally connected on parsec scales, based on 15 GHz Very Long Baseline Array (VLBA) data from a sample of 133 AGN jets. This result is achieved through a new method for measuring the product of the jet Lorentz factor and the intrinsic opening angle Gamma*theta_j from measured apparent opening angles in flux density limited samples of AGN jets. The Gamma*theta_j parameter is important for jet physics because it is related to the jet-frame sidewise expansion speed and causal connection between the jet edges and its symmetry axis. Most importantly, the standard model of jet production requires that the jet be causally connected with its symmetry axis, implying that Gamma*theta_j < 1. When we apply our method to the MOJAVE flux density limited sample of radio loud objects, we find Gamma*theta_j = 0.2, implying that AGN jets are causally connected. We also find evidence that AGN jets viewed very close to the line of sight effectively have smaller intrinsic opening angles compared with jets viewed more off-axis, which is consistent with Doppler beaming and a fast inner spine/slow outer sheath velocity field. Notably, gamma-ray burst (GRB) jets have a typical Gamma*theta_j that is two orders of magnitude higher, suggesting that different physical mechanisms are at work in GRB jets compared to AGN jets. A useful application of our result is that a jet's beaming parameters can be derived. Assuming Gamma*theta_j is approximately constant in the AGN jet population, an individual jet's Doppler factor and Lorentz factor (and therefore also its viewing angle) can be determined using two observable quantities: apparent jet opening angle and the apparent speed of jet components.
  • We report the detection of a new TeV gamma-ray source, VER J0521+211, based on observations made with the VERITAS imaging atmospheric Cherenkov telescope array. These observations were motivated by the discovery of a cluster of >30GeV photons in the first year of Fermi-LAT observations. VER J0521+211 is relatively bright at TeV energies, with a mean photon flux of 1.93 +/- 0.13_stat +/- 0.78_sys 10^-11 cm-2 s-1 above 0.2 TeV during the period of the VERITAS observations. The source is strongly variable on a daily timescale across all wavebands, from optical to TeV, with a peak flux corresponding to ~0.3 times the steady Crab Nebula flux at TeV energies. Follow-up observations in the optical and X-ray bands classify the newly-discovered TeV source as a BL Lac-type blazar with uncertain redshift, although recent measurements suggest z=0.108. VER J0521+211 exhibits all the defining properties of blazars in radio, optical, X-ray, and gamma-ray wavelengths.
  • We report observations of Faraday rotation measures (RMs) for a sample of 191 extragalactic radio jets observed within the MOJAVE program. Multifrequency VLBA observations were carried out over twelve epochs in 2006 at four frequencies between 8 and 15 GHz. We detect parsec-scale Faraday RMs in 149 sources and find the quasars to have larger RMs on average than BL Lac objects. The median core RMs are significantly higher than in the jet components. This is especially true for quasars where we detect a significant negative correlation between the magnitude of the RM and the de-projected distance from the core. We perform detailed simulations of the observational errors of total intensity, polarization and Faraday rotation, and concentrate on the errors of transverse Faraday RM gradients in unresolved jets. Our simulations show that the finite image restoring beam size has a significant effect on the observed RM gradients, and spurious gradients can occur due to noise in the data if the jet is less than two beams wide in polarization. We detect significant transverse RM gradients in four sources (0923+392, 1226+023, 2230+114 and 2251+158). In 1226+023 the RM is for the first time seen to change sign from positive to negative over the transverse cuts, which supports the presence of a helical magnetic field in the jet. In this source we also detect variations in the jet RM over a time scale of three months, which are difficult to explain with external Faraday screens and suggest internal Faraday rotation. By comparing fractional polarization changes in jet components between the four frequency bands to depolarization models we find that an external purely random Faraday screen viewed through only a few lines of sight can explain most of our polarization observations but in some sources, such as 1226+023 and 2251+158, internal Faraday rotation is needed.
  • We report on the detection of a very rapid TeV gamma-ray flare from BL Lacertae on 2011 June 28 with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The flaring activity was observed during a 34.6-minute exposure, when the integral flux above 200 GeV reached $(3.4\pm0.6) \times 10^{-6} \;\text{photons}\;\text{m}^{-2}\text{s}^{-1}$, roughly 125% of the Crab Nebula flux measured by VERITAS. The light curve indicates that the observations missed the rising phase of the flare but covered a significant portion of the decaying phase. The exponential decay time was determined to be $13\pm4$ minutes, making it one of the most rapid gamma-ray flares seen from a TeV blazar. The gamma-ray spectrum of BL Lacertae during the flare was soft, with a photon index of $3.6\pm 0.4$, which is in agreement with the measurement made previously by MAGIC in a lower flaring state. Contemporaneous radio observations of the source with the Very Long Baseline Array (VLBA) revealed the emergence of a new, superluminal component from the core around the time of the TeV gamma-ray flare, accompanied by changes in the optical polarization angle. Changes in flux also appear to have occurred at optical, UV, and GeV gamma-ray wavelengths at the time of the flare, although they are difficult to quantify precisely due to sparse coverage. A strong flare was seen at radio wavelengths roughly four months later, which might be related to the gamma-ray flaring activities. We discuss the implications of these multiwavelength results.
  • We present more than three years of observations at different frequencies, from radio to high-energy gamma-rays, of the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z=0.585). This source is the first NLS1 detected at energies above 100 MeV and therefore can be considered the prototype of this emerging new class of gamma-ray emitting active galactic nuclei (AGN). The observations performed from 2008 August 1 to 2011 December 31 confirmed that PMN J0948+0022 generates a powerful relativistic jet, able to develop an isotropic luminosity at gamma-rays of the order of 10^48 erg s^-1, at the level of powerful quasars. The evolution of the radiation emission of this source in 2009 and 2010 followed the canonical expectations of relativistic jets, with correlated multiwavelength variability (gamma-rays followed by radio emission after a few months), but it was difficult to retrieve a similar pattern in the light curves of 2011. The comparison of gamma-ray spectra before and including 2011 data suggested that there was a softening of the high-energy spectral slope. We selected five specific epochs to be studied by modelling the broad-band spectrum, characterised by an outburst at gamma-rays or very low/high flux at other wavelengths. The observed variability can largely be explained either by changes in the injected power, the bulk Lorentz factor of the jet or the electron spectrum. The characteristic time scale of doubling/halving flux ranges from a few days to a few months, depending on the frequency and the sampling rate. The shortest doubling time scale at gamma-rays is 2.3+-0.5 days. These small values underline the need of highly-sampled multiwavelength campaigns to better understand the physics of these sources.
  • (Abridged) We analyze blazar jet apparent speeds and accelerations from the RDV series of astrometric and geodetic VLBI experiments. From these experiments, we have produced and analyzed 2753 global VLBI images of 68 sources at 8 GHz with a median beam size of 0.9 milliarcseconds (mas), and a median of 43 epochs per source. From this sample, we analyze the motions of 225 jet components in 66 sources. The distribution of the fastest measured apparent speed in each source has a median of 8.3c and a maximum of 44c. Sources in the 2FGL Fermi LAT catalog display higher apparent speeds than those that have not been detected. On average, components farther from the core in a given source have significantly higher apparent speeds than components closer to the core. We measure accelerations of components in orthogonal directions parallel and perpendicular to their average velocity vector. Parallel accelerations have significantly larger magnitudes than perpendicular accelerations, implying observed accelerations are predominantly due to changes in the Lorentz factor (bulk or pattern) rather than projection effects from jet bending. Positive parallel accelerations are significantly more common than negative ones, so the Lorentz factor (bulk or pattern)tends to increase on the scales observed here. Observed parallel accelerations correspond to modest source frame increases in the bulk or pattern Lorentz factor.
  • We have investigated a frequency-dependent shift in the absolute position of the optically thick apparent origin of parsec-scale jets ("core shift" effect) to probe physical conditions in ultra-compact relativistic outflows in AGN. We used multi-frequency Very Long Baseline Array (VLBA) observations of 191 sources carried out in 12 epochs in 2006 within the MOJAVE program. The observations were performed at 8.1, 8.4, 12.1, and 15.4 GHz. We implemented a method of determining the core shift vector based on (i) image registration by two-dimensional normalized cross-correlation and (ii) model-fitting the source brightness distribution to take into account a non-zero core component offset from the phase center. The 15.4-8.1, 15.4-8.4, and 15.4-12.1 GHz core shift vectors are derived for 163 sources, and have median values of 0.128, 0.125, and 0.088 mas, respectively, compared to the typical measured errors of 0.050, 0.051, 0.035 mas. The effect occurs predominantly along the jet direction, with departures smaller than 45 deg from the median jet position angle in over 80% of the cases. Despite the moderate ratio of the observed frequencies (<2), core shifts significantly different from zero (>2sigma) are detected for about 55% of the sources. These shifts are even better aligned with the jet direction, deviating from the latter by less than 30 deg in over 90% of the cases. There is an indication that the core shift decreases with increasing redshift. Magnetic fields in the jet at a distance of 1 parsec from the central black hole, calculated from the obtained core shifts, are found to be systematically stronger in quasars (median B1~0.9 G) than those in BL Lacs (median B1~0.4 G). We also constrained the absolute distance of the core from the apex of the jet at 15 GHz as well as the magnetic field strength in the 15 GHz core region.
  • We investigate statistical and individual astrophysical properties of active galactic nuclei (AGNs), such as parsec-scale flux density, core dominance, angular and linear sizes, maximum observed brightness temperatures of VLBI core components, spectral index distributions for core and jet components, and evolution of brightness temperature along the jets. Furthermore, we statistically compare core flux densities and brightness temperature as well as jet spectral indices of gamma-ray bright and weak sources. We used 19 very long baseline interferometry (VLBI) observing sessions carried out simultaneously at 2.3 and 8.6 GHz with the participation of 10 Very Long Baseline Array (VLBA) stations and up to 10 additional geodetic telescopes. The observations span the period 1998-2003. We present here single-epoch results from high-resolution radio observations of 370 AGNs. Our VLBI images at 2.3 and 8.6 GHz as well as Gaussian models are presented and analyzed. At least one-fourth of the cores are completely unresolved on the longest baselines of the global VLBI observations. The VLBI core components are partially opaque with the median value of spectral index of alpha_core=0.3, while the jet features are usually optically thin alpha_jet=-0.7. The spectral index typically decreases along the jet ridge line owing to the spectral aging, with a median value of -0.05 mas^-1. Brightness temperatures are found to be affected by Doppler boosting and reaching up to \sim10^13 K with a median of \sim2.5x10^11 K at both frequencies. The brightness temperature gradients along the jets typically follow a power law T_b\simr^-2.2 at both frequencies. 147 sources (40%) positionally associated with gamma-ray detections from the Fermi LAT Second Source Catalog have higher core flux densities and brightness temperatures, and are characterized by the less steep radio spectrum of the optically thin jet emission.
  • We have conducted a survey of Faraday rotation in a sample of 191 compact radio-loud AGNs as part of the MOJAVE (Monitoring of Jets in Active galactic nuclei with VLBA Experiments) project. The observations were carried out with the VLBA at 8.1, 8.4, 12.1 and 15.3 GHz over 12 epochs in 2006. We detect sufficiently strong linear polarization in 159 out of 211 observations to calculate the rotation measure values, resulting in a large enough sample for statistical analysis of the Faraday rotation in blazars. These Faraday rotation measures can be used to study the intrinsic magnetic field order and orientation in parsec-scale blazar jets. Our sample includes 119 sources listed in the 1FGL or 2FGL catalogs and we detect rotation measure values in 111 out of 131 maps. Of the 72 sources that are not in the gamma-ray catalogs we detect RM in 48 out of 80 maps. The median RM values of the LAT-detected sources do not differ significantly from the non-LAT-detected sources. Nine of the sources in our sample have resolved enough jets to study the transverse Faraday rotation structure, and we detect significant transverse rotation measure gradients in four sources. In two of these (3C~273 and 3C~454.3) there is additional evidence to support helical magnetic field in the parsec-scale jets. The two others (0923+392 and 2230+114) require further observations to identify the nature of the gradient. It is interesting that three of the four sources with significant rotation measure gradients are sources that have shown large gamma-ray flares.
  • We report about the multiwavelength campaign on the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z = 0.5846) performed in 2010 July-September and triggered by high activity as measured by Fermi/LAT. The peak luminosity in the 0.1-100 GeV energy band exceeded, for the first time in this type of source, the value of 10^48 erg/s, a level comparable to the most powerful blazars. The comparison of the spectral energy distribution of the NLS1 PMN J0948+0022 with that of a typical blazar - like 3C 273 - shows that the power emitted at gamma rays is extreme.
  • The Compact Symmetric Objects (CSOs) are small (<1 kiloparsec) and powerful extragalactic radio sources showing emission on both sides of an active galactic nucleus and no signs of strong relativistic beaming. They may be young radio sources, progenitors of large FRII radio galaxies. We aim to study the statistical properties of CSOs by constructing and investigating a new large sample of CSO candidates on the basis of dual-frequency, parsec-scale morphology. For the candidate selection we utilized VLBI data for 4170 extragalactic objects obtained simultaneously at 2.3 and 8.6 GHz (S and X band) within the VLBA Calibrator Survey 1-6 and the Research and Development - VLBA projects. Properties of their broad-band radio spectra were characterized by using RATAN-600 observations. Numerical modeling was applied in an attempt to explain the observed effects. A sample of 64 candidate CSOs is identified. The median two-point S-X band spectral index of parsec-scale hotspots is found to be -0.52; with the median brightness temperature ~10^9 K at X band. Statistical analysis reveals a systematic difference between positions of brightest CSO components (associated with hotspots) measured in the S and X bands. The distance between these components is found to be on average 0.32+/-0.06 mas greater at 8.6 GHz than at 2.3 GHz. This difference in distances cannot be explained by different resolutions at the S and X bands. It is a manifestation of spectral index gradients across CSO components, which may potentially provide important physical information about them. Despite our detailed numerical modeling of a CSO hotspot, the model was not able to reproduce the magnitude of the observed positional difference. A more detailed modeling may shed light on the origin of the effect.
  • M. L. Lister, M. Aller, H. Aller, T. Hovatta, K. I. Kellermann, Y. Y. Kovalev, E. T. Meyer, A. B. Pushkarev, E. Ros, M. Ackermann, E. Antolini, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, K. Bechtol, R. Bellazzini, B. Berenji, R. D. Blandford, E. D. Bloom, M. Boeck, E. Bonamente, A. W. Borgland, J. Bregeon, M. Brigida, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, P. A. Caraveo, J. M. Casandjian, E. Cavazzuti, C. Cecchi, C. S. Chang, E. Charles, A. Chekhtman, C. C. Cheung, J. Chiang, S. Ciprini, R. Claus, J. Cohen-Tanugi, J. Conrad, S. Cutini, F. de Palma, C. D. Dermer, E. do Couto e Silva, P. S. Drell, A. Drlica-Wagner, C. Favuzzi, S. J. Fegan, E. C. Ferrara, J. Finke, W. B. Focke, P. Fortin, Y. Fukazawa, P. Fusco, F. Gargano, D. Gasparrini, N. Gehrels, S. Germani, N. Giglietto, F. Giordano, M. Giroletti, T. Glanzman, G. Godfrey, I. A. Grenier, S. Guiriec, D. Hadasch, M. Hayashida, E. Hays, D. Horan, R. E. Hughes, G. Jóhannesson, A. S. Johnson, M. Kadler, H. Katagiri, J. Kataoka, J. Knödlseder, M. Kuss, J. Lande, F. Longo, F. Loparco, B. Lott, M. N. Lovellette, P. Lubrano, G. M. Madejski, M. N. Mazziotta, W. McConville, J. E. McEnery, J. Mehault, P. F. Michelson, T. Mizuno, C. Monte, M. E. Monzani, A. Morselli, I. V. Moskalenko, S. Murgia, M. Naumann-Godo, S. Nishino, P. L. Nolan, J. P. Norris, E. Nuss, M. Ohno, T. Ohsugi, A. Okumura, N. Omodei, E. Orlando, M. Ozaki, D. Paneque, D. Parent, M. Pesce-Rollins, M. Pierbattista, F. Piron, G. Pivato, S. Rainò, A. Readhead, A. Reimer, O. Reimer, J. L. Richards, S. Ritz, H. F.-W. Sadrozinski, C. Sgrò, M. S. Shaw, E. J. Siskind, G. Spandre, P. Spinelli, H. Takahashi, T. Tanaka, J. G. Thayer, J. B. Thayer, D. J. Thompson, G. Tosti, A. Tramacere, E. Troja, T. L. Usher, J. Vandenbroucke, V. Vasileiou, G. Vianello, V. Vitale, A. P. Waite, P. Wang, B. L. Winer, K. S. Wood, S. Zimmer
    Aug. 31, 2011 astro-ph.CO, astro-ph.HE
    We investigate the Fermi LAT gamma-ray and 15 GHz VLBA radio properties of a joint gamma-ray- and radio-selected sample of AGNs obtained during the first 11 months of the Fermi mission (2008 Aug 4 - 2009 Jul 5). Our sample contains the brightest 173 AGNs in these bands above declination -30 deg. during this period, and thus probes the full range of gamma-ray loudness (gamma-ray to radio band luminosity ratio) in the bright blazar population. The latter quantity spans at least four orders of magnitude, reflecting a wide range of spectral energy distribution (SED) parameters in the bright blazar population. The BL Lac objects, however, display a linear correlation of increasing gamma-ray loudness with synchrotron SED peak frequency, suggesting a universal SED shape for objects of this class. The synchrotron self-Compton model is favored for the gamma-ray emission in these BL Lacs over external seed photon models, since the latter predict a dependence of Compton dominance on Doppler factor that would destroy any observed synchrotron SED peak - gamma-ray loudness correlation. The high-synchrotron peaked (HSP) BL Lac objects are distinguished by lower than average radio core brightness temperatures, and none display large radio modulation indices or high linear core polarization levels. No equivalent trends are seen for the flat-spectrum radio quasars (FSRQ) in our sample. Given the association of such properties with relativistic beaming, we suggest that the HSP BL Lacs have generally lower Doppler factors than the lower-synchrotron peaked BL Lacs or FSRQs in our sample.
  • The effect of a frequency dependent shift of the VLBI core position (known as the "core shift") was predicted more than three decades ago and has since been observed in a few sources, but often within a narrow frequency range. This effect has important astrophysical and astrometric applications. To achieve a broader understanding of the core shift effect and the physics behind it, we conducted a dedicated survey with NRAO's Very Long Baseline Array (VLBA). We used the VLBA to image 20 pre-selected sources simultaneously at nine frequencies in the 1.4-15.4 GHz range. The core position at each frequency was measured by referencing it to a bright, optically thin feature in the jet. A significant core shift has been successfully measured in each of the twenty sources observed. The median value of the core shift is found to be 1.21 mas if measured between 1.4 and 15.4 GHz, and 0.24 mas between 5.0 and 15.4 GHz. The core position, r, as a function of frequency, n, is found to be consistent with an r n^-1 law. This behavior is predicted by the Blandford & Koenigl model of a purely synchrotron self-absorbed conical jet in equipartition. No systematic deviation from unity of the power law index in the r(n) relation has been convincingly detected. We conclude that neither free-free absorption nor gradients in pressure and/or density in the jet itself and in the ambient medium surrounding the jet play a significant role in the sources observed within the 1.4-15.4 GHz frequency range. These results support the interpretation of the parsec-scale core as a continuous Blandford-Koenigl type jet with smooth gradients of physical properties along it.
  • We have compared the parsec-scale jet linear polarization properties of the Fermi LAT-detected and non-detected sources in the complete flux-density-limited (MOJAVE-1) sample of highly beamed AGN. Of the 123 MOJAVE sources, 30 were detected by the LAT during its first three months of operation. We find that during the era since the launch of Fermi, the unresolved core components of the LAT-detected jets have significantly higher median fractional polarization at 15 GHz. This complements our previous findings that these LAT sources have higher apparent jet speeds, brightness temperatures and Doppler factors, and are preferentially found in higher activity states.
  • We report the discovery of high-energy (E>100 MeV) gamma-ray emission from NGC 1275, a giant elliptical galaxy lying at the center of the Perseus cluster of galaxies, based on observations made with the Large Area Telescope (LAT) of the Fermi Gamma ray Space Telescope. The positional center of the gamma-ray source is only ~3' away from the NGC 1275 nucleus, well within the 95% LAT error circle of ~5'.The spatial distribution of gamma-ray photons is consistent with a point source. The average flux and power-law photon index measured with the LAT from 2008 August 4 to 2008 December 5 are F_gamma = (2.10+-0.23)x 10^{-7} ph (>100 MeV) cm^{-2} s^{-1} and Gamma = 2.17+-0.05, respectively. The measurements are statistically consistent with constant flux during the four-month LAT observing period.Previous EGRET observations gave an upper limit of F_gamma < 3.72x 10 ^{-8} ph (>100 MeV) cm^{-2} s^{-1} to the gamma-ray flux from NGC 1275. This indicates that the source is variable on timescales of years to decades, and therefore restricts the fraction of emission that can be produced in extended regions of the galaxy cluster. Contemporaneous and historical radio observations are also reported. The broadband spectrum of NGC 1275 is modeled with a simple one-zone synchrotron/synchrotron self-Compton model and a model with a decelerating jet flow.