The first IRAM/PdBI polarimetric millimeter survey of active galactic nuclei
II. Activity and properties of individual sources^{⋆,}^{⋆⋆}
^{1} Seoul National University, Department of Physics and Astronomy, 599 Gwanakro, Gwanakgu, 151742 Seoul, South Korea
email: trippe@astro.snu.ac.kr
^{2} Institut de Radioastronomie Millimétrique (IRAM), 300 rue de la Piscine, 38406 Saint Martin d’Hères, France
Received: 2 December 2011
Accepted: 15 February 2012
We present an analysis of the linear polarization of six active galactic nuclei – 0415+379 (3C 111), 0507+179, 0528+134 (OG+134), 0954+658, 1418+546 (OQ+530), and 1637+574 (OS+562). Our targets were monitored from 2007 to 2011 in the observatoryframe frequency range 80–253 GHz, corresponding to a restframe frequency range 88–705 GHz. We find average degrees of polarization m_{L} ≈ 2 − 7%; this indicates that the polarization signals are effectively averaged out by the emitter geometries. From a comparison of the fluctuation rates in flux and degree of polarization we conclude that the spatial scales relevant for polarized emission are of the same order of, but probably not smaller than, the spatial scales relevant for the emission of the total flux. We see indication for fairly strong shocks and/or complex, variable emission region geometries in our sources, with compression factors ≲ 0.9 and/or changes in viewing angles by ≳ 10°. An analysis of correlations between source fluxes and polarization parameter points out special cases: the presence of (at least) two distinct emission regions with different levels of polarization (for 0415+379) as well as emission from a single, predominant component (for 0507+179 and 1418+546). Regarding the evolution of flux and polarization, we find good agreement between observations and the signal predicted by “oblique shock in jet” scenarios in one source (1418+546). We attempt to derive rotation measures for all sources, leading to actual measurements for two AGN and upper limits for three sources. We derive values of RM = (−39 ± 1_{stat} ± 13_{sys}) × 10^{3} rad m^{2} and RM = (42 ± 1_{stat} ± 11_{sys}) × 10^{4} rad m^{2} for 1418+546 and 1637+574, respectively; these are the highest values reported to date for AGN. These values indicate magnetic field strengths of the order ~10^{4} G. For 0415+379, 0507+179, and 0954+658 we derive upper limits RM < 1.7 × 10^{4} rad m^{2}. From the relation RM ∝ ν^{a} we find a = 1.9 ± 0.3 for 1418+546, in good agreement with a = 2 as expected for a spherical or conical outflow.
Key words: galaxies: active / quasars: general / radiation mechanisms: nonthermal / polarization / techniques: polarimetric
This study is based on observations carried out with the IRAM Plateau de Bure Interferometer. IRAM is supported by INSU/CNRS (France), MPG (Germany), and IGN (Spain).
Table 4 is available in electronic form at http://www.aanda.org
© ESO, 2012
1. Introduction
Active galactic nuclei (AGN) have been studied extensively in the wavelength range from cmradio to γ radiation in the last decades (see, e.g., Kembhavi & Narlikar 1999; or Krolik 1999, and references therein for a review). There is overwhelming observational evidence that the main source of their emission is accretion onto supermassive black holes (SMBH) with masses M_{•} ≃ 10^{6...9} M_{⊙} (e.g., Ferrarese & Ford 2005, and references therein). Important information on the physics of active galacic nuclei is encoded in their linear polarization m_{L}. Degrees and angles of polarization provide details on synchrotron emission, the geometry of emission regions, strength and orientation of magnetic fields, and (via Faraday rotation and/or depolarization) on particle densities and matter distributions of the surrounding or outflowing matter (see, e.g., Saikia & Salter 1988, and references therein).
Physical properties and observed frequency ranges for our six target AGN.
In Trippe et al. (2010, hereafter Paper I), we presented the results of our polarimetric AGN survey with the IRAM Plateau de Bure Interferometer (PdBI; e.g., Winters & Neri 2011)^{1}. This survey observed 86 sources in the (observatory frame) frequency range 80–267 GHz. Sources were not resolved spatially, all properties observed were sourceintegrated ones. Paper I provided a statistical analysis of the entire survey sample. The main results of this analysis were

For 73 out of the 86 AGN observed linear polarization was detected. The median degree of polarization was ≈ 4%, the highest value ever measured was ≈ 19%. These numbers are far below the value ≳ 60% expected from “ideal” synchrotron emission from optically thin sources, pointing toward a substantial averaging of polarization due to limited spatial resolution and nonuniform Bfield geometries, and/or a substantial contribution by optically thick emission from galactic cores.

On average, BL Lacertae objects show a higher degree of polarization (m_{L} ≃ 7%) than Seyfert galaxies or QSOs (m_{L} ≃ 3% and m_{L} ≃ 5%, respectively). This “polarization level sequence” BL Lac → QSOs → Seyferts can be understood approximately in the context of the viewing angle unification scheme of AGN: a more direct view into a more active nucleus permits observation of more concentrated emission. This reduces the impact of averaging of the polarization through source geometry.

We did not find evidence for correlations of polarization with restframe frequency or redshift z. The absence of a m_{L} − z correlation indicates that there has been no noticeable change of polarization properties since z ≃ 2.4.

There is a trend toward higher degrees of polarization in sources with flat spectral mm/radio indices that indicate partially optically thick, coredominated emission. This is somewhat unexpected because in general, optically thin AGN outflows are more polarized than optically thick cores. A possible explanation is – again – provided by source geometry: because the jet emission is distributed over much larger spatial areas than the nuclear emission, averaging of polarization through nonuniform magnetic field geometries could be more efficient in the outflows than in the cores, resulting in observation of the highest polarization levels in the most coredominated sources.
In this article, Paper II, we focus on the linear polarization properties of AGN that have been observed multiple times and resolved in time and frequency. We selected sources from the sample analyzed in Paper I (1) for which significant polarization was detected at least 20 times, or (2) that had been observed at least 25 times and polarization had been detected at least 15 times. When selecting targets, these limits proved to be the minimum numbers necessary for reliable statistical statements. Observations were obtained between February 2007 and January 2011 at frequencies 82–253 GHz. Eventually, this left us with a sample of six targets: 0415+379, 0507+179, 0528+134, 0954+658, 1418+546, and 1637+574. We present the basic physical properties of our target sources in Table 1.
2. Observations and data processing
2.1. Polarization
In January 2007 (January 2008 for the 2 mmband), the antennas of the PdBI were equipped with dual linear polarization Cassegrain focus receivers. Since then, it is possible to observe both orthogonal polarizations – “horizontal” (H) and “vertical” (V) with respect to the antenna frame – simultaneously. Observations can be carried out (nonsimultaneously) in three atmospheric windows located around wavelengths of 1.3 mm, 2 mm, and 3 mm. Each of these bands covers a continuous range of frequencies; their ranges are 201–267 GHz for the 1.3mm band, 129–174 GHz for the 2mm band, and 80–116 GHz for the 3mm band. Within a given band, any frequency is available for observations. Receivers covering a fourth wavelength band centered at 0.8 mm began operations at the end of 2010 but did not provide polarization data as of January 2011.
At the time of the observations presented here, the PdBI was not yet equipped for observations of all Stokes parameters. We collected linear polarization data on point sources via Earth rotation polarimetry, i.e. by monitoring the fluxes in the H and V channels as functions of parallactic angle ψ. For assessing the polarization of a source, we calculated the parameter (1)from the fluxes H(ψ) and V(ψ). The second equality means that q(ψ) provides full information on linear polarization (see, e.g., Sault et al. 1996; Thompson et al. 2001; but also Beltrán et al. 2004) if a sufficient range of ψ is observed.
Fig. 1 Linear polarization state as function of time t and observation frequency ν. We give the result of each polarization measurement via a black bar in the t − ν plane. The length of each bar corresponds to the degree of linear polarization according to the scales in the upper left corners of the corresponding panels (in units of %). The orientation of the bar indicates the polarization angle; angles are counted counterclockwise from the 12o’clock position. Points in the t − ν plane where a measurement was carried out but no polarization was detected are indicated by gray bars. The length of each gray bar corresponds to a 3σ upper limit on the polarization level. As upper limits do not have meaningful polarization angles, all gray bars are fixed at the 12o’clock position. 

Open with DEXTER 
Fig. 2 3mm flux densities of our target AGN as function of time. Please note the different axis scales. The errors given are dominated by a ≈ 10% systematic uncertainty of the flux scale calibration; therefore error bars tend to be smaller for lower fluxes. 

Open with DEXTER 
Owing to the nature of polarized light and because the PdBI antenna receivers are located in the Cassegrain foci, observing a polarized target results in q(ψ) being a cosinusoidal signal with a period of 180°. The functional form of q(ψ) is thus (2)Here m_{L} is the fraction of linear polarization (ranging from 0 to 1; in the following, we will express m_{L} in units of %) and χ is the polarization angle (ranging from 0° to 180°). A more detailed discussion of the methodology is provided in Paper I. Our observational results are summarized in Table 4.
2.2. Flux densities
Source fluxes are initially recorded as antenna temperatures. Antenna temperatures are converted into physical flux densities (quoted in units of Jansky throughout this paper) by using empirical antenna efficiencies as conversion factors. Those factors are functions of frequencies and are located in the range from ≈ 22 Jy K^{1} (for the 3mm band) to ≈37 Jy K^{1} (for the 1.3mm band). To estimate the systematic uncertainties of our dataset, we analyzed observations of the radio continuum star MWC 349A (e.g. Tafoya et al. 2004) that were obtained and calibrated like the AGN measurements. Owing to its wellknown flux properties MWC 349 serves as a “standard candle” for the PdBI. From the systematic scatter of the MWC 349 lightcurves we conclude that the systematic relative uncertainties of our AGN observations are ≈ 10% for the 2mm and 3mm bands and ≈ 16% for the 1.3mm band. For more details regarding AGN flux monitoring with the PdBI, please see Trippe et al. (2011).
We note that the systematic accuracy of our flux data is limited by the automatic monitoring and calibration. When calibrating observations individually and interactively while using MWC 349 as reference, systematic errors can be reduced to ≲ 5%; Krips et al. (in prep.) are going to discuss the recent developments regarding the PdBI flux calibration and the resulting improved performance.
3. Analysis
3.1. Variability of polarization and flux
Our analysis focuses on the properties of AGN for which polarization has been observed multiple ( ≥ 15) times. We take into account that the polarization is, a priori, a function of two independent parameters: time and frequency. Accordingly, we may present the polarization of our six target AGN in twodimensional diagrams with the two dimensions being observing time and observing frequency, respectively. The result is shown in Fig. 1. The diagrams provide degrees of polarization and polarization angles for each measurement that detected significant polarization and upper limits on polarization levels otherwise.
Estimates on the variability of fluxes and polarization levels.
In Fig. 2, we present the 3mm lightcurves of our target sources. To prevent systematic errors caused by spectral slope effects, we limit our variability analysis to flux measurements obtained in the 3mm band.
To assess variabilities quantitatively, we define the redshiftcorrected fluctuation rate of a timedependent parameter x like (3)where z is the redshift, x_{i} is the ith (i = 1,2,3,...) value of the dataset, t_{i} the time when x_{i} was obtained, and N the total number of measurements; ... denotes the absolute value of the enclosed function. By construction, ξ_{x} is a measure of the rapidity of variations of the parameter x, taking explicitly into account the time axis of the dataset. Thus our parametrization provides additional information compared to commonly used, simpler variability estimators such as the σ_{V} parameter proposed by Hook et al. (1994). We note that our ξ is similar to a parameter introduced by Villforth et al. (2009), which they named “violence”.
Fig. 3 Degree of linear polarization m_{L} as function of rescaled flux density . Please note the different axis scales. Errorbars denote 1σ errors. In each diagram we give the Pearson correlation coefficient r. We also quote (in units of %) for each value of r the significance s = 1 − p with p being the falsealarm probability for the null hypothesis “the true value of r is zero”. 

Open with DEXTER 
Fig. 4 Rescaled linearly polarized flux P^{ ∗ } vs. rescaled flux density . Please note the different axis scales. Errorbars denote 1σ errors. We quote correlation coefficients r and significance levels s like in Fig. 3. 

Open with DEXTER 
Fig. 5 Polarization angle χ vs. rescaled flux density . Please note the different axis scales. Polarization angles are restricted to the interval [0°,180°] ; for 0415+379 and 0954+658 these intervals have been shifted to resolve 0°/180° ambiguities. Errorbars denote 1σ errors. We quote correlation coefficients r and significance levels s like in Fig. 3. 

Open with DEXTER 
Fig. 6 Degree of linear polarization m_{L} as function of observing frequency ν. Squares with errorbars denote polarization measurements; diamonds indicate 3σ upper limits. Please note the different m_{L} axis scales. Errorbars denote 1σ errors. We quote correlation coefficients r and significance levels s like in Fig. 3. 

Open with DEXTER 
Fig. 7 Polarization angle χ as function of the square of the restframe wavelength λ_{0}. Polarization angles are restricted to the interval [0°,180°] ; for 0954+658 and 1637+574 these intervals have been shifted to resolve 0°/180° ambiguities. The selection of data is limited to the epoch ranges quoted on top of each plot to maximize the significance of a rotation measure signal. We quote correlation coefficients r and significance levels s like in Fig. 3. For 0415+379, 0507+179, and 0954+658, we derive 3σ upper limits on RM. For 1418+546 and 1637+574, we quote the rotation measures together with their statistical and systematic errors. Rotation measures are derived via linear fits (gray lines) to the data. 

Open with DEXTER 
In Table 2, we present fluctuation rates for the degree of linear polarization m_{L}, ξ_{m}, as well as for the 3mm flux densities S_{ν}, ξ_{S}. In addition, we show normalized fluctuation rates ξ_{m}/ ⟨ m_{L} ⟩ and ξ_{S}/ ⟨ S_{ν} ⟩ , with ⟨ ... ⟩ denoting the time average of the corresponding parameter. Finally, we also present the values of the fluctuation rate ratio (4)For the two fluctuation rates, ξ_{m} and ξ_{S}, we calculated statistical errors numerically. Because all data points are affected by measurement errors, the calculated fluctuation rates deviate from zero even for signals with zero “true” variability. Our error analysis is based on a Monte Carlo scheme: we repeated each calculation 1000 times with different realizations of a given data set. The measurement error is defined as the central 68% (1σ) interval of the resulting distribution of ξ_{x} values. A realization of a data set is drawn by adding random, Gaussian errors – which are in accord with the individual statistical uncertainties – to all elements of a data set. For the ξ_{m}, errors are small compared to the parameter values and are quoted explicitly in Table 2. For the ξ_{S}, uncertainty ranges are asymmetric and reach the magnitudes of the parameter values; we therefore do not quote errors individually but restrict ourselves to the statement that the ξ_{S} values are accurate to factors ≈ 2.
3.2. Correlations between flux and polarization
The interplay between source polarizations and fluxes promises new insights into the emission mechanisms. We therefore study the relation between fluxes S_{ν} on the one hand and polarization levels m_{L}, polarized fluxes P, and polarization angles χ on the other hand. Our observational procedures measure m_{L} according to Eq. (2) but not P. Thus we calculate linearly polarized fluxes like P = m_{L}S_{ν} by combining those – initially separate – polarization and flux measurements obtained simultaneously at the same frequency. Because flux densities are a function of observation frequency ν like S_{ν} ∝ ν^{ − α}, we apply spectral rescaling: (5)This procedure means that we scale all observed flux densities to an observing frequency ν = 90 GHz. We used the spectral index information from Paper I (see also Table 1). Accordingly, we computed rescaled polarized fluxes . For the special case of 0415+379, we included only 3mm flux densities into our calculations. For this AGN our flux data are affected by a combination of rapid flux variability and irregular sampling that introduce additional systematic errors when combining several wavelength bands; see also Trippe et al. (2011). For the remaining five sources we included all data available regardless of frequency. We note that rescaling is based on spectral slopes derived from an ensemble of flux measurements. Accordingly, rescaled parameters can only be used when analyzing the same ensemble of flux data collectively, though not for the analysis of individual or differences between individual measurements as in Sect. 3.1. Our subsequent correlation analysis is based on the rescaled parameters and P^{ ∗ }.
In general, we quantify the correlation between two parameters x, y by means of the Pearson correlation coefficient (6)(e.g., Snedecor & Cochran 1980). By construction, r values of ( − 1) +1 correspond to perfect (anti)correlation, a value of 0 means absence of any correlation. For each AGN, we computed correlation coefficients with x being and y being either m_{L} or P^{ ∗ } or χ, thus providing^{2}, , and . We present the distributions and correlations of our parameters in Figs. 3–5.
Inspecting the distributions in Figs. 3–5 shows that the presence or absence of correlations is far from obvious due to the substantial scatter within our data. To permit quantitative conclusions, we calculated significances by means of a permutation test. For each correlation coefficient r(x,y) we computed a sample of modified coefficients r(x,y′) with y′ being a copy of the initial parameter vector y with its elements rearranged in random order. Because the pairs (x,y′) are drawn randomly, the expected value for r(x,y′) is zero. From a total of 10 000 random realizations we determined the fraction p for which The fraction p is thus the false alarm probability, i.e. the probability that a value r(x,y) ≠ 0 is produced by a statistical fluctuation. Accordingly, our randomization test probes the null hypothesis “the parameters x and y are intrinsically uncorrelated”. In the following discussion as well as in Figs. 3–5, we quote significance levels s = 1 − p in units of %.
3.3. Relations between polarization and frequency
Because our observations are resolved in frequency, our data permit probing correlations between polarization properties and frequency. In Fig. 6, we present the degree of linear polarization m_{L} as function of observing frequency ν. We quote Pearson correlation coefficients and their significances computed analogously to the procedure outlined in Sect. 3.2. As already noted in Paper I (see especially the top panel of Fig. 6 therein), a positive correlation – if any – may be expected already because of observational bias: measurements obtained at higher frequencies tend to suffer from reduced signaltonoise (S/N) ratios. This effect reduces the probability of detecting low values of m_{L} at high frequencies, resulting in an apparent increase of the average m_{L} with frequency.
A more promising approach is the analysis of the polarization angle χ. Other than for m_{L}, we are not aware of any mechanism that could introduce a bias affecting the observed values of χ. This is a convenient situation because it permits an analysis of Faraday rotation. For linearly polarized light propagating through a plasma permeated by a magnetic field, the polarization angle is related to the restframe wavelength^{3}λ_{0} like (7)where RM is the rotation measure (8)with B_{} being the strength of the magnetic field parallel to the line of sight (l.o.s.), n_{e} being the electron number density, and l being the coordinate directed along the l.o.s. (e.g., Wilson et al. 2010). Thus, deriving RM permits constraining the l.o.s.integrated magnetic field strength and particle density of an active galaxy.
In Fig. 7 we present for each source observed polarization angles χ as function of squared restframe wavelength λ_{0}. Our data are affected by strong temporal variability. Where possible, we therefore limit our analysis to time windows for which the polarization appears coherent (compare Fig. 1, Table 4). We selected those time windows such that a minimum number of data points is excluded and the significance of a potential rotation measure signal is enhanced. We note that this “windowing” does not bias our results because it applies a selection in time but not in polarization angle or frequency. We quote Pearson correlation coefficients and their significance levels computed analogously to the procedure outlined in Sect. 3.2. For all sources except 0528+134 the observed polarization angles (after windowing) are distributed over ranges notably smaller than the permitted [0,π] intervals. For these AGN we attempted to derive rotation measures via linear fits to the data in accord with Eq. (7). Subsequently, we assumed that the observed distributions of χ vs. are superpositions of two principal effects: (1) a linear trend due to Faraday rotation, and (2) scatter due to temporal variability of the AGN polarization that is independent of wavelength. In this case the RM is given by the slope of the straight line fitting the data best. We claim a successful rotation measure measurement if, and only if, the data show a linear correlation r between χ vs. that is significantly different from zero.
We note that the statistical errors of the RM appear to be very small when regarding the scatter of the data. This is a consequence of our “twoeffect assumption” outlined above. Our RM measurements used a χ^{2} minimization algorithm^{4} that fits linear models to the data. For models with few parameters, the 1σ confidence intervals of the parameters are given by variations of χ^{2} around its minimum of order unity. In our case however, the scatter of the data is much larger than the statistical error of any individual data point; therefore even the bestfitting solutions have χ^{2} in the order of several thousands. From this it follows that the statistical errors we derive are meaningful if, and only if, our initial “twoeffect assumption” is applicable. Obviously, the accuracy of our results is limited by systematic errors due to intrinsic temporal variability; statistical errors due to random measurement errors are hardly relevant. We therefore eventually quote for each RM value two errors: the statistical fit error and a systematic error. The systematic error essentially quantifies the uncertainty due to the temporal variability of χ. We derived systematic errors from the confidence levels s of the Pearson correlation coefficients, expressed in units of Gaussian σ. For a significance level s that corresponds to n × σ in Gaussian terms (with n being a positive real number), we quote RM/n as an estimate of the systematic error. Our results are summarized in Table 3.
Rotation measures RM for our six target sources.
4. Results
4.1. 0415+379 (3C 111)
For this Seyfert 1 galaxy located at z ≈ 0.05, 17 out of 34 polarimetric observations, obtained between March 2007 and January 2010, detected linearly polarized emission. The average degree of polarization is ⟨ m_{L} ⟩ ≈ 2.3%. Initially, low polarization levels below 2% and polarization angles around 0° (modulo 0°/180° ambiguities) are predominant. In August 2007, the polarization changed to m_{L} ≈ 4 − 6% and χ ≈ 100°. After April 2008, polarization could only be detected at one occasion (April 2009, m_{L} ≈ 1%,χ ≈ 120°). The degree of polarization fluctuates at (normalized) rates of ξ_{m} ≈ 0.09% d^{1} and ξ_{m}/ ⟨ m_{L} ⟩ ≈ 0.04 d^{1}, respectively.
The 3mm flux density was highly variable, with the average value being ⟨ S_{ν} ⟩ ≈ 7.5 Jy. In mid2007, flux measurements commenced at values S_{ν} ≈ 10 Jy, rising to S_{ν} ≈ 15 Jy shortly thereafter. Since then, the flux declined to S_{ν} ≈ 2 Jy, with a shortlived intermediate high of S_{ν} ≈ 7 Jy observed at the end of 2008. The flux density fluctuates with ξ_{S} ≈ 0.7 Jy d^{1}, the normalized fluctuation rate being ξ_{S}/ ⟨ S_{ν} ⟩ ≈ 0.1 d^{1}. From a comparison of the fluctuation rates in flux and polarization, one obtains a fluctuation rate ratio of ρ ≈ 0.4.
Degree of polarization and rescaled flux show a strong (r = −0.9) anticorrelation. With s = 99.91% (corresponding to ≈ 3.3σ in terms of a Gaussian distribution), this anticorrelation is significant. An anticorrelation (r = −0.51) is apparent also when comparing rescaled polarized flux P^{∗} and ; however, as s ≈ 94% (corresponding to ≈ 1.9σ), this signal is statistically insignificant. Additionally, polarization angle χ and show a strong (r = 0.91) correlation, which is significant at a level of s = 99.91% (≈3.3σ).
A search for a correlation between m_{L} and observing frequency ν provides a null result (r = 0.37, s = 92%). When comparing polarization angle χ versus squared restframe wavelength , we formally derive a rotation measure RM = (−4.0 ± 2.3) × 10^{3} rad m^{2} (statistical error only). As s = 71%, this signal is insignificant. Accordingly, we quote a 3σ upper limit on RM given by the absolute value of the RM formally derived plus three times its statistical error. Eventually, we find RM < 1.1 × 10^{4} rad m^{2}.
4.2. 0507+179
This quasar, located at z ≈ 0.42, shows significant polarization in 21 out of 22 measurements obtained between September 2007 and January 2011. The average polarization level is ⟨ m_{L} ⟩ ≈ 7.4%. In 2007, the degree of polarization is m_{L} ≈ 6%, declines to m_{L} ≈ 3% in the first half of 2008, rises again to levels up to m_{L} ≈ 11%, and decreases again to m_{L} ≈ 8%. In 2007 and 2008, the polarization angle varies in the range χ ≈ 50...110°, remaining around χ ≈ 100° thereafter. The degree of polarization fluctuates at rates of ξ_{m} ≈ 0.4 % d^{1} and ξ_{m}/ ⟨ m_{L} ⟩ ≈ 0.05 d^{1}, respectively.
The 3mm flux evolved smoothly, with values S_{ν} ≈ 0.7 Jy in 2007, peaking at S_{ν} ≈ 1.5 Jy at the beginning of 2009, and decreasing again to S_{ν} ≈ 0.5 Jy. The average flux density is found to be ⟨ S_{ν} ⟩ ≈ 1 Jy. It fluctuates with ξ_{S} ≈ 0.1 Jy d^{1}, corresponding to a normalized fluctuation rate of ξ_{S}/ ⟨ S_{ν} ⟩ ≈ 0.1 d^{1}. The polarization – flux fluctuation rate ratio is ρ ≈ 0.7.
Polarization level and rescaled flux density appear to be uncorrelated (r = 0.22, s ≈ 79%). When comparing rescaled polarized flux and rescaled flux density, we find a positive correlation (r = 0.73), which is statistically significant at a level of s = 99.96%, corresponding to ≈ 3.5σ. Polarization angle and rescaled flux are uncorrelated (r = −0.08, s ≈ 61%).
From a comparison of m_{L} and observing frequency ν we do not find any significant correlation (r = 0.3, s ≈ 91%). We find a slightly positive correlation (r = 0.3) between polarization angle χ and the square of the restframe wavelength . We can formally derive a rotation measure of RM = (9.3 ± 1.6) × 10^{3} rad m^{2} (statistical error only). However, as the significance of the correlation is only s = 60% we have to reject this result as insignificant. Accordingly, we quote a 3σ upper limit of RM < 1.4 × 10^{4} rad m^{2}.
4.3. 0528+134 (OG+134)
For this highpolarization quasar located at z ≈ 2.1, 18 out of 32 polarimetric observations, obtained between February 2007 and December 2008, detected linearly polarized emission. The average degree of polarization is ⟨ m_{L} ⟩ ≈ 3.2%. During the observation timeline, polarization levels varied in the range m_{L} ≈ 2...6%, whereas the polarization angles varied in the range χ ≈ 0...140°. The degree of polarization fluctuates at rates of ξ_{m} ≈ 0.3% d^{1} and ξ_{m}/ ⟨ m_{L} ⟩ ≈ 0.1 d^{1}, respectively.
The 3mm flux density was at its peak of S_{ν} ≈ 5 Jy when observations commenced in 2007 and decreased steadily to levels of S_{ν} ≈ 1 Jy in 2010. The average of the flux density is ⟨ S_{ν} ⟩ ≈ 3.5 Jy. The flux density fluctuates with ξ_{S} ≈ 1.1 Jy d^{1}, the normalized fluctuation rate being ξ_{S}/ ⟨ S_{ν} ⟩ ≈ 0.3 d^{1}. By comparing the fluctuation rates in flux and polarization, we obtain a fluctuation rate ratio of ρ ≈ 0.3.
The degree of polarization and the rescaled flux do not show a correlation (r = −0.21, s ≈ 72%). When analyzing rescaled polarized flux P^{∗} versus , we see a positive correlation (r = 0.73). However, as s = 99.45%, being equivalent to ≈ 2.8σ, this signal is only marginally significant. For the correlation between polarization angle and rescaled flux we find r = −0.37 and s ≈ 87%.
Polarization level and observing frequency ν are uncorrelated (r = −0.01). A comparison of polarization angle χ and squared restframe wavelength finds that the values observed for χ occupy the entire range [0,π] ; evidently, strong temporal variability renders our rotation measure analysis inconclusive even when using “windowing”.
4.4. 0954+658
This BL Lacertae object, located at z ≈ 0.37, shows significant polarization in 35 out of 37 measurements obtained between March 2007 and January 2011. The average polarization level is ⟨ m_{L} ⟩ ≈ 7.4%. 0954+658 shows by far the highest degrees of polarization among our sources, ranging up to m_{L} ≈ 18% in 2007. In 2009, the polarization drops to m_{L} ≈ 5% and increases again up to m_{L} ≈ 10% in 2010/2011. The polarization angles remain around χ ≈ 0° for most of the observing timeline; the only notable exception occurs in 2009 when values of χ ≈ 10...30° occur. The degree of polarization fluctuates at a rate of ξ_{m} ≈ 0.5% d^{1} and a normalized rate ξ_{m}/ ⟨ m_{L} ⟩ ≈ 0.07 d^{1}, respectively.
The average of the 3mm flux density is ⟨ S_{ν} ⟩ ≈ 1.2 Jy. The flux varies remains within the range S_{ν} ≈ 0.5...1.5 Jy throughout the observing timeline, the only exception being a short peak with S_{ν} ≈ 2.5 Jy occurring mid2008. The flux fluctuation rate is ξ_{S} ≈ 0.1 Jy d^{1}, corresponding to a normalized fluctuation rate of ξ_{S}/ ⟨ S_{ν} ⟩ ≈ 0.1 d^{1}. The fluctuation rate ratio of polarization and flux is ρ ≈ 0.8.
Polarization level m_{L} and rescaled flux are uncorrelated (r = −0.31, s ≈ 96%). Likewise, rescaled polarized flux P^{∗} and are uncorrelated (r = 0.2, s ≈ 85%). Polarization angle χ and rescaled flux are uncorrelated (r = 0.36, s ≈ 97%).
A correlation analysis of m_{L} versus observing frequency ν unveils a positive correlation with r = 0.52. As s ≈ 99.9% ( ≈ 3.2σ), the correlation is statistically significant. We find a slightly negative correlation (r = −0.21) between polarization angle χ and the square of the restframe wavelength . Regarding the polarization angle, we can formally derive a rotation measure of RM = (−14 ± 1) × 10^{3} rad m^{2} (statistical error only). However, as the significance of the correlation is s ≈ 99.2% ( ≈ 2.6σ), the correlation falls short of being statistically significant. Accordingly, we quote a 3σ upper limit of RM < 1.7 × 10^{4} rad m^{2}.
4.5. 1418+546 (OQ+530)
For this BL Lacertae object located at z ≈ 0.15, 52 out of 76 polarimetric observations, obtained between January 2008 and January 2011, detected linearly polarized emission. The average degree of polarization is ⟨ m_{L} ⟩ ≈ 5.0%. Polarization levels vary in the range m_{L} ≈ 3...10% with a tendency for the highest values to occur at the beginning and the lowest values to be found at the end of the observed period. The degree of polarization fluctuates at rates of ξ_{m} ≈ 0.3% d^{1} and ξ_{m}/ ⟨ m_{L} ⟩ ≈ 0.06 d^{1}, respectively. The angles of polarization are located in the range χ ≈ 110...150°. From inspecting the corresponding time–frequency diagram (Fig. 1) we note a “swing pattern” for the polarization angles, with the highest values of χ to be found around the beginning of 2010.
The 3mm lightcurve shows a doublepeak morphology, with values ranging from S_{ν} ≈ 0.5 Jy at the beginnings of 2009 and 2011 to S_{ν} ≈ 1.3 Jy at the beginnings of 2008 and 2010. The average value of the flux is ⟨ S_{ν} ⟩ ≈ 0.8 Jy. The flux density fluctuates with ξ_{S} ≈ 0.04 Jy d^{1}, the normalized fluctuation rate being ξ_{S}/ ⟨ S_{ν} ⟩ ≈ 0.05 d^{1}. From comparison of the fluctuation rates in flux and polarization, we obtain a fluctuation rate ratio of ρ ≈ 1.1.
When comparing m_{L} and , we find a positive correlation (r = 0.57), which is significant as s ≈ 99.98% (i.e. ≈ 3.7σ). An analysis of rescaled polarized flux P^{∗} versus indicates a strong (r = 0.75) and significant (s ≈ 99.99%, corresponding to ≈ 3.9σ) positive correlation. Polarization angle and flux are uncorrelated (r = 0.31, s ≈ 97%).
Polarization level and observing frequency ν show a positive correlation (r = 0.48), which appears to be significant (s ≈ 99.6%, i.e. ≈ 2.9σ). We find a significant – with s ≈ 99.6%, meaning ≈2.9σ – negative correlation (r = −0.35) between polarization angle χ and the square of the restframe wavelength . We are able to derive a rotation measure of RM = (−39 ± 1_{stat} ± 13_{sys}) × 10^{3} rad m^{2}.
4.6. 1637+574 (OS+562)
This Seyfert 1 galaxy, located at z ≈ 0.75, shows significant polarization in 17 out of 26 measurements obtained between February 2007 and August 2009. The average polarization level is ⟨ m_{L} ⟩ ≈ 3.0%. Degrees of polarization are found in the range m_{L} ≈ 1...5%, with the lowest values occurring at the beginning of the observed period. The angles of polarization swing from χ ≈ 150° to χ ≈ 90°. The degree of polarization fluctuates at a rate of ξ_{m} ≈ 0.2% d^{1} and a normalized rate ξ_{m}/ ⟨ m_{L} ⟩ ≈ 0.05 d^{1}, respectively.
The 3mm flux density is initially observed at values S_{ν} ≈ 1.2 Jy, increases up to S_{ν} ≈ 2 Jy in 2008 and decreases again to S_{ν} ≈ 1 Jy in 2009. The average flux is ⟨ S_{ν} ⟩ ≈ 1.2 Jy. The flux fluctuation rate is ξ_{S} ≈ 0.17 Jy d^{1}, corresponding to a normalized fluctuation rate of ξ_{S}/ ⟨ S_{ν} ⟩ ≈ 0.12 d^{1}. The fluctuation rate ratio of polarization and flux is ρ ≈ 0.4.
Polarization level m_{L} and rescaled flux are uncorrelated (r = 0.1). When comparing rescaled polarized flux P^{∗} and , we find a positive correlation with r = 0.55; however, as s ≈ 97%, the correlation is insignificant. Likewise, polarization angle and flux are uncorrelated (r = −0.22, s ≈ 76%).
An analysis of m_{L} versus observing frequency finds a positive correlation (r = 0.46), which is insignificant as s ≈ 97%. We find a strong, positive correlation (r = 0.78) between polarization angle χ and the square of the restframe wavelength . As s ≈ 99.98%, corresponding to ≈ 3.7σ, the correlation is significant. The rotation measure we derive is RM = (42 ± 1_{stat} ± 11_{sys}) × 10^{4} rad m^{2}.
5. Discussion
5.1. Polarization levels
We present the polarization properties of six mm/radioluminous – with S_{ν} ≳ 0.8 Jy in the 3mm band – AGN. Our sample is a mix of Seyfert galaxies, BL Lacertae objects, and quasars. An inspection of their mm/radio spectral indices (as listed in Table 1) indicates that the radiation of two sources – 0415+379 and 0528+134 – is optically thin and outflowdominated, whereas the remaining four AGN show optically thick, coredominated emission when applying the conventional division between the two cases placed at α = 0.5 (e.g., Krolik 1999).
The average degrees of polarization m_{L} we observe are located between ≈ 2% and ≈ 7%. This is substantially lower than the values one can derive from the theory of synchrotron emission from homogeneous and isotropic ensembles of relativistic electrons moving in uniform magnetic fields (e.g., Ginzburg & Syrovatskii 1965; Pacholczyk 1970; see also the corresponding discussion in Paper I). Given the range of α observed, we may expect theoretical polarization levels from m_{L} ≈ 13% (for α = 0.3, optically thick emission) to m_{L} ≈ 74% (for α = 0.9, optically thin emission). The most probable explanation is provided by geometric averaging: our observations do not resolve our targets spatially, whereas extended outflows (jets or halos) have been observed in all of them. More specifically, Lee et al. (2008) have been able to map four sources – 0415+379, 0528+134, 0954+658, and 1637+574 – using global VLBI observations at 86 GHz. In addition, Sokolovsky et al. (2011) observed an optically thick core as well as an optically thin jet in VLBA radio maps obtained in the range 1.4–15 GHz. These maps show jets for all sources observed, indicating that probably all our targets show optically thin outflows if observed with sufficient angular resolution. Accordingly we have to assume that any of our observations takes the average of multiple emission regions with different optical depths and different magnetic field orientations. Assuming we observe N emission zones with randomly distributed magnetic field orientations in any given source, we may expect sourceintegrated degrees of polarization . From our observations (see Tables 1, 2) we can then estimate values between N ≈ 3 (for 0507+179 and 0954+658, assuming optically thick emission) and N ≈ 1000 (for 0415+379, assuming optically thin emission). The variation in N mirrors the variation in relevant emission region sizes: AGN radio cores are placed at spatial scales of parsecs whereas jets occur on scales of many kiloparsecs ( ≈ 80 kpc for the radio jet of 0415+379; Liu & Xie 1992). However, this simple “static” picture is complicated by any dynamical evolution of the emission regions; we discuss the potential signatures of shocks in the following subsections.
An independent check of our results is provided by Agudo et al. (2010). This study provides singleepoch measurements of linear polarization obtained between 2005 and 2006 at 86 GHz for five of our targets. The authors found degrees of polarization of < 1.5% (i.e., an upper limit only), 2.6 ± 0.5%, 5.2 ± 0.5%, 2.8 ± 0.6%, and < 1.5% for 0415+379, 0528+134, 0954+658, 1418+546, and 1637+574, respectively. In addition, Agudo et al. (2010) derived polarization angles of 91 ± 6°, 8 ± 3°, and 108 ± 6° for 0528+134, 0954+658, and 1418+546, respectively. These values may be compared with our results given in Tables 2 and 4 as well as Fig. 1; taking into account the temporal variability, both datasets agree well. This indicates that neither Agudo et al. (2010) nor our study seem to have observed a source in a temporary atypical or extreme state. This reduces the probability that any of the two studies is affected by temporal selection effects.
5.2. Variability
Variability of AGN polarization is wellknown in the cm/radio regime (e.g., Aller et al. 1985) and has been probed in the mm/radio regime for a few selected sources (e.g., Jorstad et al. 2005) as well as in the optical (e.g., Villforth et al. 2009). Our variability analysis (see Eqs. (3), (4) and Table 2) treats variability in source flux S_{ν} and degree of polarization m_{L} jointly (see also Figs. 1, 2). Both parameters fluctuate at relative rates of few per cent per day. For the flux densities, those fluctuation rates have been observed in cm/radio as well as mm/radio bands and have been discussed in the context of potential intraday variability (e.g., Fuhrmann et al. 2008). Rapid flux variability is the basis of the famous historical statement that the size of the (relevant parts of the) individual emission regions cannot exceed a few light days (Smith & Hoffleit 1963). Our analysis makes it possible to extend this statement to the observed polarization: because the fluctuation rates of flux density and degree of polarization are compatible, the spatial scales of the corresponding emission regions should be compatible as well. We may quantify this by means of the fluctuation rate ratios ρ as given by Eq. (4). Our analysis finds values ρ ≈ 0.3 − 1.1, permitting two principal conclusions:

1.
Within their uncertainties of factors ≈ 2, the ρ values are on the order of unity for all our six target sources. This suggests that ratios ρ ≈ 1 are indeed characteristic for AGN emission (at least in mm/radio) – pointing toward a strong structural correlation between total and polarized flux.

2.
Even though all ratios ρ are consistent with being on the order of unity, the individual values tend toward numbers somewhat smaller than one^{5}. This indicates that the size scales relevant for the polarized emission – like turbulence or eddy scales of the magnetic fields involved – are compatible with, but not smaller than, the relevant size scales for the overall emission.
Fig. 8 Constraints on the values and the variability of the shock parameter δ. The shock parameter δ relates to the factor μ = δ/(2 − δ) by which linear polarization is reduced by shock compression. Our observations provide Δμ = 2(δ_{1} − δ_{2})/(2 − δ_{1})(2 − δ_{2}) ≈ σ_{m}/ ⟨ m_{L} ⟩ for two arbitrary points in time “1”, “2”. The grayshaded area marks the δ_{1} − δ_{2} plane range consistent with the observed σ_{m}/ ⟨ m_{L} ⟩ ≈ 0.33 − 0.65 for our sample of sources. 

Open with DEXTER 
Rapid – referring here to timescales shorter than a few months – variability of polarization is supposedly caused by shocks propagating through magnetized plasma (e.g., Hughes et al. 1985). In this case, the degree of linear polarization of synchrotron radiation from a partially compressed plasma is reduced by a factor μ compared to the case without compression (Hughes et al. 1985; Cawthorne & Wardle 1988), where (9)with the shock parameter (10)Here 0 ≤ k ≤ 1 is the compression factor, meaning the factor by which the length of the affected region is reduced by compression. The angle ϵ is the angle between the line of sight and the plane of compression in the frame of the emitter. From Eq. (10) follows 0 ≤ δ ≤ 1. Additional information from, e.g., timeresolved radio interferometric maps, is not available, making it impossible to disentangle k and ϵ in δ; accordingly, δ is our only observable. As discussed in Sect. 5.1, we probably observe polarization from multiple emission regions; in this case, any attempt to constrain the parameters of Eqs. (9) and (10) actually provides sourceaveraged, “typical” values. Assuming we observe different values for μ at two points of time “1” and “2”, we may write (11)Our observations do not provide the absolute, intrinsic degrees of polarization – meaning the polarization levels corresponding to the case when shocks are absent – of individual emission regions. Instead, we observe sourceaveraged (except, possibly, 1418+546; see Sect. 5.3) polarization levels that fluctuate with time. From this we attempt to constrain at least average, “typical” shock parameters δ. For a given source, we use the time average of the degree of polarization ⟨ m_{L} ⟩ as an approximate reference. We then use the normalized standard deviation of the degree of polarization σ_{m}/ ⟨ m_{L} ⟩ as an estimate for Δμ, meaning (12)Evidently, Δμ depends on the individual values of δ_{1} and δ_{2} as well as on their difference δ_{1} − δ_{2} (see Eq. (11)). Accordingly, it is not possible to derive single values for δ_{1,2} but ranges in a plane spanned by δ_{1} and δ_{2}. In Fig. 8 we present such a plane with the ranges marked that are permitted for σ_{m}/ ⟨ m_{L} ⟩ ≈ 0.33 − 0.65 as observed for our sources (compare Table 2). Our result points toward the occurrence of fairly violent shocks and/or multiple shocks with varying geometries within our sources. According to Eq. (10) and regarding, e.g., the pair δ_{1} ≈ 1 and δ_{2} ≈ 0.8 in Fig. 8, our observations indicate the occurrence of k ≲ 0.9, or of variations in ϵ by ≳ 10°, or combinations of both.
5.3. Correlations between flux and polarization
Immediate correlations between rescaled (Eq. (5)) polarization parameters and source flux (Figs. 3–5) can occur in several varieties. Because we might observe multiple emission components within the same source, we can write the observed linear polarization vector P as (13)with twodimensional vectors P = (Q,U) with Stokes parameters Q and U. To separate out the variability in source flux S_{ν}, we also write a normalized polarization (14)with the Stokes parameter I being identical to S_{ν}; the index i indicates the contribution from the ith emission component. Here p = (Q/I,U/I) ≡ (q,u), which relates to our observed parameters like and χ = 0.5tan^{1}(u/q). In general, one may assume that all parameters involved can vary independently, leaving the polarization parameters m_{L}, χ and rescaled fluxes uncorrelated; indeed, this is what our data suggest for most of our targets. Nevertheless, we are able to probe at least two special cases:
 1.
We receive radiation from a single predominant emission component with an arbitrary but approximately fixed degree and angle of polarization. Any change in source flux should be accompanied by a proportional change in polarized flux. In this case we would observe a strong correlation between flux density and polarized flux. The flux would not correlate with degree or angle of polarization.
 2.
We receive radiation from two (or more) emission components with each component having an arbitrary but approximately constant q_{i} and u_{i}. The observed polarization signal is a superposition of the individual components. Any increase (or decrease) in the luminosity of one component increases (decreases) its contribution to the average signal observed. In this scenario, we would observe a tight correlation between flux and degree of polarization as well as between flux and polarization angle^{6}.
For 0507+179, we observe a signature in agreement with scenario 1: while the source flux increases by a factor ≈3, the polarized flux increases by a factor ≈ 3, too. There is no significant correlation between flux on the one hand and degree or angle of polarization on the other hand. Radio interferometric VLBA maps obtained by Fey & Charlot (1997) at 8.6 GHz show a dominating core and a fainter jet component about 3 mas west of the core. However, as this map was obtained in 1995, any comparison between the data of Fey & Charlot (1997) and our results has to be taken with care. We note that our “onecomponent scenario” does not contradict the discussion in Sect. 5.1. The scenarios are reconciled when assuming largescale regions with coordinated luminosity variations – like shock fronts – that have substructure defined by, e.g., characteristic ordering scales of magnetic fields.
In 1418+546, we observe a more complex signature. While the source flux increases by a factor ≈ 3, the polarized flux inceases by a factor ≈ 8, resulting in a significant correlation between the two parameters in agreement with scenario 1. Additionally, there is also a significant, though more scattered, positive correlation between flux and degree of polarization, with m_{L} rising by a factor ≈ 4. There is no significant correlation between flux and polarization angle. In essence, we seem to observe a single, dominating emission component that tends toward higher polarization levels at higher flux densities. We note the agreement of the behavior of 1418+546 with the signatures expected from recent model calculations based on “shock in jet” scenarios (Hughes et al. 2011). Regarding the lightcurve, our observations indeed appear to cover two outbursts occurring within about two years. Our observations find (a) a correlated variation of flux and polarization, (b) changes of the polarization angle by tens of degrees, and (c) no correlation between flux and polarization angle. We see the best agreement of our data with the scenario of an oblique shock front observed at an angle > 60° to the jet axis (see Fig. 7f of Hughes et al. 2011). We also note that the calculations of Hughes et al. (2011) assume a spectral index α = 0.25; this may be compared to the observational result α = 0.2 ± 0.1 (see Table 1 and Paper I). Overall, our observations appear to provide new support for the “oblique shock” scenario.
5.4. Rotation measures
We analyzed the rotation measures of our target AGN according to the procedure outlined in Sect. 3.3. We are able to derive statistically significant rotation measures for 1418+546 and 1637+574 (see Fig. 7 and Sects. 4.5, 4.6). We note that for these two sources we included all data available into our calculation because “time windowing” did not increase the significance of the result or reduce the scatter caused by temporal variability. Our calculations assume that any Faraday rotation detected is caused by matter in or within the vicinity of the emitters. A priori, there are two additional sources of Faraday rotation that might distort our analysis: the terrestrial ionosphere and Galactic interstellar matter. For the frequency range of our observations, ν > 80 GHz, the ionosphere is able to rotate the polarization angle of electromagnetic radiation by less than 0.01° (Thompson et al. 2001), meaning its influence may be neglected. Galactic interstellar matter introduces rotation measures RM ≲ 200 at radio frequencies (Taylor et al. 2009). This number is about one order of magnitude smaller than the statistical errors of our rotation measure values, which renders the influence of interstellar matter irrelevant.
We find RM = (−39 ± 1_{stat} ± 13_{sys}) × 10^{3} rad m^{2} and RM = (42 ± 1_{stat} ± 11_{sys}) × 10^{4} rad m^{2} for 1418+546 and 1637+574, respectively. For 0415+379, 0507+179, and 0954+658 we derive 3σ upper limits on RM of 1.1 × 10^{4} rad m^{2}, 1.4 × 10^{4} rad m^{2}, and 1.7 × 10^{4} rad m^{2}, respectively. Evidently, our observations are subject to strong selection effects. As already noted by Agudo et al. (2010), Eq. (7) implies that, at mmwavelengths, changes in polarization angle as small as ≈5° require RM ≈ 7 × 10^{3} rad m^{2}. In addition, our polarization angle data show substantial scatter caused by rapid temporal variability. Therefore rotation measures smaller than a few 10^{4} rad m^{2} probably remain unnoticed simply because they do not imprint a statistically significant signal into a diagram. This is illustrated by the case of 1418+546, where an RM ≈ − 3.9 × 10^{4} rad m^{2} signal still comes with a statistical falsealarm probability of ≈ 0.4%, barely corresponding to a 3σ significance threshold. Rotation measure values derived from cm/radio observations are usually located in ranges from a few 10 rad m^{2} to a few 1000 rad m^{2} at most (e.g., O’Sullivan & Gabuzda 2009; Taylor et al. 2009). At mm/radio frequencies, AGN polarization and Faraday rotations are poorly probed, with the most comprehensive study to date probably being the work by Jorstad et al. (2007). Those studies observe values up to RM ≈ 3.4 × 10^{4} rad m^{2} (for CTA 102; Jorstad et al. 2007). Accordingly, we report in this study the highest rotation measures observed to date in AGN.
From Eq. (8) we may estimate the product of magnetic field strength and electron density averaged along the line of sight ⟨ B_{}n_{e} ⟩ . For emission region sizes on the order of one parsec (compare, e.g., O’Sullivan & Gabuzda 2009, for the case of 1418+546) and RM of order 10^{5} rad m^{2}, we find ⟨ B_{}n_{e} ⟩ ~ 0.1 G cm^{3}. For realistic electron densities n_{e} ~ 10^{3} cm^{3} (e.g., Koski 1978; Bennert et al. 2009) this corresponds to average magnetic field strengths B_{} ~ 10^{4} G.
Observing increasing rotation measures when going from cm/radio to mm/radio frequencies may be expected from Eq. (8) for optically thick sources. Higher frequencies probe regions closer to the central engine with higher electron densities and stronger magnetic fields. For realistic emission region geometries, one may assume n_{e} ∝ l^{ − a} and B_{} ∝ l^{1}, with n_{e} being the electron density, B_{} being the magnetic field parallel to the line of sight, l being the distance from the core where the observed emission predominantly originates, and a being a positive real number (Jorstad et al. 2007). For outflows with spherical or conical geometries, a = 2. From Eq. (8) one derives immediately RM ∝ l^{ − a}. For optically thick sources, distance l and frequency ν are related like l ∝ ν^{1}, leading to a “core shift” effect. For the specific case of 0507+179, Sokolovsky et al. (2011) explicitly observed a core shift for the frequency range 1.4–15 GHz. Combining the various relations leads to (15)when inserting values for a in agreement with observations, i.e. a ≈ 1−3 (O’Sullivan & Gabuzda 2009), one may indeed expect RM of order 10^{5} rad m^{2}. This is additionally supported by observations of Sagittarius A*, the radio source commonly equated with the supermassive black hole at the center of the Milky Way (e.g., Genzel et al. 2010, and references therein). Here rotation measures up to RM ≈ 5 × 10^{5} rad m^{2} at frequencies around 340 GHz have been reported (Marrone et al. 2006; Macquart et al. 2006).
For 1418+546 it is possible to compare our result with RM observations obtained at cm/radio frequencies and from these to derive a. O’Sullivan & Gabuzda (2009) report RM = −501 ± 48 rad m^{2} at observed frequencies around 15 GHz. This can be related to our result, RM = (−39 ± 1_{stat} ± 13_{sys}) × 10^{3} rad m^{2} at ν ≈ 150 GHz. From Eq. (15) we obtain a = 1.9 with a formal statistical error of only ± 0.04. From the systematic error in RM and the uncertainties in defining the two reference frequencies used in our calculation we expect an additional systematic error of approximately 0.3. Eventually, we obtain a = 1.9 ± 0.3 with the error being predominantly of systematic nature. This may be compared to the theoretical value a = 2 for spherical or conical outflows. Our rotation measure analysis therefore indicates that for the specific case of 1418+546 the assumption of a simple spherical or conical outflow geometry is sufficient.
6. Summary and conclusions
We have analyzed the linear polarization of six AGN – 0415+379 (3C 111), 0507+179, 0528+134 (OG+134), 0954+658, 1418+546 (OQ+530), and 1637+574 (OS+562) – in the observatory frame frequency range 80–253 GHz, corresponding to a range of rest frame frequencies 88–705 GHz. Our data were obtained in the years 2007–2011 and provide information resolved in time and frequency for all targets. Our sample includes two optically thin, outflowdominated and four optically thick, coredominated emitters. This study arrives at the following principal conclusions:

1.
In agreement with the results of Paper I, the observed degrees of polarization m_{L} are lower than the levels expected from the theory of synchrotron emission by up to one order of magnitude. This indicates that the polarization signal is effectively averaged out. Assuming radiation from N emission regions with random orientations and geometries, we expect . Our results point to values from N ≈ 3 (for two optically thick emitters) to N ≈ 1000 (0415+379, an optically thin emitter).

2.
Degrees of polarization as well as fluxes vary with relative rates of a few per cent per day. We find fluctuation rate ratios ρ (Eq. (4)) on the order of unity though typically slightly lower than one. This indicates that the size scales relevant for the polarized emission are compatible with, but probably not smaller than, the size scales relevant for the overall emission.

3.
The degrees of polarization in our sample fluctuate with normalized standard deviations σ_{m}/ ⟨ m_{L} ⟩ ≈ 0.33−0.65. We identify this variability approximately with sourceaveraged fluctuations expected from shock compression of magnetized plasmas. Our calculations point toward fairly strong shocks and/or complex, variable shock geometries, with compression factors k ≲ 0.9 and/or changes in viewing angles ≳ 10°.

4.
A search for correlations between source fluxes and polarization parameters identifies two special cases. For 0415+379, we observe a strong anticorrelation between flux and degree of polarization. This indicates the presence of (at least) two distinct emission regions with different levels of polarization. For 0507+179 and 1418+546, we find a positive correlation between flux and polarized flux. Here the data point toward emission from a single, predominant component.

5.
For 1418+546 we find (a) a correlated variation of flux and polarization, (b) changes of the polarization angle by tens of degrees, and (c) no correlation between flux and polarization angle. Comparison with the “oblique shock” models by Hughes et al. (2011) indicates good agreement with the presence of an oblique shock front observed at an angle > 60° to the jet axis.

6.
We derive rotation measures of RM = (−39 ± 1_{stat} ± 13_{sys}) × 10^{3} rad m^{2} and RM = (42 ± 1_{stat} ± 11_{sys}) × 10^{4} rad m^{2} for 1418+546 and 1637+574, respectively. These are the highest values reported to date for AGN except for Sagittarius A*. Our results point to magnetic field strengths on the order of ~10^{4} G. For 0415+379, 0507+179, and 0954+658 we are able to derive upper limits RM < 1.7 × 10^{4} rad m^{2}.

7.
Comparing our rotation measure values with those obtained by other studies at lower frequencies, we are able to constrain the emission region geometry of 1418+546. Exploiting the relation RM ∝ ν^{a} we find a = 1.9 ± 0.3, in good agreement with a = 2 as expected for a spherical or conical outflow.

8.
We do not see a relation between the observed polarization properties of our sources and their physical parameters such as black hole masses or kpcmorphology, as summarized in Table 1.
This paper is the second of a set of two articles. Earth rotation polarimetry of AGN is a byproduct of regular PdBI operations; accordingly, the database is going to grow steadily in the future. We therefore expect that a new study – analogous to the one outlined in this article as well as in Paper I – carried out in a few years from now promises new insights already because of substantially improved statistics. In addition, the PdBI crew currently evaluates the implementation of an observing mode that permits observations of all four Stokes parameters. Evidently, such an observing mode would permit additional studies if implemented eventually.
Our work demonstrates the power of mm/radio polarization studies of AGN based on data resolved in time and frequency. Given the strong temporal variability of AGN polarization, it becomes clear that more accurate measurements of differential parameters such as rotation measures require simultaneous multifrequency observations. New facilities such as the Korean VLBI Network (KVN; e.g., Kim et al. 2010), which permits singlepolarization observations simultaneously at four frequencies – namely 22, 43, 86, and 129 GHz – or dualpolarization
observations simultaneously at two of the aforementioned frequencies, will be important tools for new discoveries.
Acknowledgments
We are grateful to the entire IRAM and PdBI staff, whose hard work over many years made this study possible. Our work has made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. We have also made use of data from the MOJAVE database that is maintained by the MOJAVE team (Lister et al. 2009). For data processing and analysis we applied the software package GILDAS, developed and maintained by the GILDAS team (http://www.iram.fr/IRAMFR/GILDAS/), as well as the software package DPUSER (http://www.mpe.mpg.de/~ott/dpuser/index.html), developed and maintained by Thomas Ott at MPE Garching. Last but not least, we are grateful to the referee, E. Ros, whose review improved the quality of our paper.
References
 Agudo, I., Thum, C., Wiesemeyer, H., & Krichbaum, T. P. 2010, ApJSS, 189, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Aller, H. D., Aller, M. F., Latimer, G. E., & Hodge, P. E. 1985, ApJSS, 59, 513 [NASA ADS] [CrossRef] [Google Scholar]
 Beltrán, M. T., Gueth, F., Guilloteau, S., & Dutrey, A. 2004, A&A, 416, 631 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bennert, N., Barvainis, R., Henkel, C., & Antonucci, R. 2009, A&A, 459, 55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cawthorne, T. V., & Wardle, J. F. C. 1988, ApJ, 332, 696 [NASA ADS] [CrossRef] [Google Scholar]
 D’Elia, V., Padovani, P., & Landt, H. 2003, MNRAS, 339, 1081 [NASA ADS] [CrossRef] [Google Scholar]
 Falomo, R., Carangelo, N., & Treves, A. 2003, MNRAS, 343, 505 [NASA ADS] [CrossRef] [Google Scholar]
 Fan, Z.H., & Cao, X. 2004, ApJ, 602, 103 [NASA ADS] [CrossRef] [Google Scholar]
 Ferrarese, L., & Ford, H. 2005, SSR, 116, 523 [NASA ADS] [CrossRef] [Google Scholar]
 Fey, A. L., & Charlot, P. 1997, ApJSS, 111, 95 [NASA ADS] [CrossRef] [Google Scholar]
 Fuhrmann, L., Krichbaum, T. P., Witzel, A., et al. 2008, A&A, 490, 1019 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Genzel, R., Eisenhauer, F., & Gillessen, S. 2010, Rev. Mod. Phys., 82, 3121 [NASA ADS] [CrossRef] [Google Scholar]
 Ghisellini, G., Tavecchio, F., & Ghirlanda, G. 2009, MNRAS, 399, 2041 [NASA ADS] [CrossRef] [Google Scholar]
 Ginzburg, V. L., & Syrovatskii, S. I. 1965, ARA&A, 3, 297 [NASA ADS] [CrossRef] [Google Scholar]
 Grossberger, C., Kadler, M., Wilms, J., et al. 2011, Proc. 8th INTEGRAL/BART Workshop 2011 [arXiv:1110.1197] [Google Scholar]
 Hook, I. M., McMahon, R. G., Boyle, B. J., & Irwin, M. J. 1994, MNRAS, 268, 305 [NASA ADS] [CrossRef] [Google Scholar]
 Hughes, P. A., Aller, H. D., & Aller, M. F. 1985, ApJ, 298, 301 [NASA ADS] [CrossRef] [Google Scholar]
 Hughes, P. A., Aller, M. F., & Aller, H. D. 2011, ApJ, 735, 81 [NASA ADS] [CrossRef] [Google Scholar]
 Jorstad, S. G., Marscher, A. P., Lister, M. L., et al. 2005, AJ, 130, 1418 [NASA ADS] [CrossRef] [Google Scholar]
 Jorstad, S. G., Marscher, A. P., Stevens, J. A., et al. 2007, AJ, 134, 799 [NASA ADS] [CrossRef] [Google Scholar]
 Kembhavi, A. K., & Narlikar, J. V. 1999, Quasars and Active Galactic Nuclei (Cambridge University Press) [Google Scholar]
 Kim, J. H., Cho, S.H., Oh, C. S., & Byun, D.Y., et al. 2010, ApJSS, 188, 209 [NASA ADS] [CrossRef] [Google Scholar]
 Koski, A. T. 1978, ApJ, 223, 56 [NASA ADS] [CrossRef] [Google Scholar]
 Krolik, J. H. 1999, Active Galactic Nuclei (Princeton University Press) [Google Scholar]
 Lee, S.S., Lobanov, A. P., Krichbaum, T. P., et al. 2008, AJ, 136, 159 [NASA ADS] [CrossRef] [Google Scholar]
 Lister, M. L., Aller, H. D., Aller, M. F., et al. 2009, AJ, 137, 3718 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, F. K., & Xie, G. Z. 1992, A&ASS, 95, 249 [NASA ADS] [Google Scholar]
 Liu, Y., Jiang, D. R., & Gu, M. F. 2006, ApJ, 637, 669 [NASA ADS] [CrossRef] [Google Scholar]
 Macquart, J.P., Bower, G. C., Wright, M. C. H., Backer, D. C., & Falcke, H. 2006, ApJ, 646, L111 [NASA ADS] [CrossRef] [Google Scholar]
 Marchesini, D., Celotti, A., & Ferrarese, L. 2004, MNRAS, 351, 733 [NASA ADS] [CrossRef] [Google Scholar]
 Marrone, D. P., Moran, J. M., Zhao, J.H., & Rao, R. 2006, ApJ, 640, 308 [NASA ADS] [CrossRef] [Google Scholar]
 O’Sullivan, S. P., & Gabuzda, D. C. 2009, MNRAS, 393, 429 [NASA ADS] [CrossRef] [Google Scholar]
 Pacholczyk, A. G. 1970, Radio Astrophysics (Freeman) [Google Scholar]
 Saikia, D. J., & Salter, C. J. 1988, ARA&A, 26, 93 [NASA ADS] [CrossRef] [Google Scholar]
 Sault, R. J., Hamaker, J. P., & Bregman, J. D. 1996, A&AS, 117, 149 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Smith, H. J., & Hoffleit, D. 1963, Nature, 198, 650 [NASA ADS] [CrossRef] [Google Scholar]
 Snedecor, G. W., & Cochran, W. G. 1980, Statistical Methods (Iowa State University Press) [Google Scholar]
 Sokolovsky, K. V., Kovalev, Y. Y., Pushkarev, A. B., & Lobanov, A. P. 2011, A&A, 532, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tafoya, D., Gómez, Y., & Rodríguez, L. F. 2004, ApJ, 610, 827 [NASA ADS] [CrossRef] [Google Scholar]
 Taylor, G. B., Gugliucci, N. E., Fabian, A. C., et al. 2006, MNRAS, 368, 1500 [NASA ADS] [CrossRef] [Google Scholar]
 Taylor, A. R., Stil, J. M., & Sunstrom, C. 2009, ApJ, 702, 1230 [NASA ADS] [CrossRef] [Google Scholar]
 Thompson, A. R., Moran, J. M., & Swenson, G. W. 2001, Interferometry and Synthesis in Radio Astronomy, 2nd edn. (WileyInterscience) [Google Scholar]
 Trippe, S., Neri, R., Krips, M., et al. 2010, A&A, 515, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Trippe, S., Krips, M., Piétu, V., et al. 2011, A&A, 533, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Villforth, C., Nilsson, K., Østensen, R., et al. 2009, MNRAS, 397, 1893 [NASA ADS] [CrossRef] [Google Scholar]
 Wilson, T. L., Rohlfs, K., & Hüttemeister, S. 2010, Tools of Radio Astronomy, 5th edn. (Springer) [Google Scholar]
 Winters, J. M., & Neri, R. 2011, An Introduction to the IRAM Plateau de Bure Interferometer, IRAM manual, version 4.200 [Google Scholar]
Online material
Observations journal of our AGN linear polarization study.
All Tables
All Figures
Fig. 1 Linear polarization state as function of time t and observation frequency ν. We give the result of each polarization measurement via a black bar in the t − ν plane. The length of each bar corresponds to the degree of linear polarization according to the scales in the upper left corners of the corresponding panels (in units of %). The orientation of the bar indicates the polarization angle; angles are counted counterclockwise from the 12o’clock position. Points in the t − ν plane where a measurement was carried out but no polarization was detected are indicated by gray bars. The length of each gray bar corresponds to a 3σ upper limit on the polarization level. As upper limits do not have meaningful polarization angles, all gray bars are fixed at the 12o’clock position. 

Open with DEXTER  
In the text 
Fig. 2 3mm flux densities of our target AGN as function of time. Please note the different axis scales. The errors given are dominated by a ≈ 10% systematic uncertainty of the flux scale calibration; therefore error bars tend to be smaller for lower fluxes. 

Open with DEXTER  
In the text 
Fig. 3 Degree of linear polarization m_{L} as function of rescaled flux density . Please note the different axis scales. Errorbars denote 1σ errors. In each diagram we give the Pearson correlation coefficient r. We also quote (in units of %) for each value of r the significance s = 1 − p with p being the falsealarm probability for the null hypothesis “the true value of r is zero”. 

Open with DEXTER  
In the text 
Fig. 4 Rescaled linearly polarized flux P^{ ∗ } vs. rescaled flux density . Please note the different axis scales. Errorbars denote 1σ errors. We quote correlation coefficients r and significance levels s like in Fig. 3. 

Open with DEXTER  
In the text 
Fig. 5 Polarization angle χ vs. rescaled flux density . Please note the different axis scales. Polarization angles are restricted to the interval [0°,180°] ; for 0415+379 and 0954+658 these intervals have been shifted to resolve 0°/180° ambiguities. Errorbars denote 1σ errors. We quote correlation coefficients r and significance levels s like in Fig. 3. 

Open with DEXTER  
In the text 
Fig. 6 Degree of linear polarization m_{L} as function of observing frequency ν. Squares with errorbars denote polarization measurements; diamonds indicate 3σ upper limits. Please note the different m_{L} axis scales. Errorbars denote 1σ errors. We quote correlation coefficients r and significance levels s like in Fig. 3. 

Open with DEXTER  
In the text 
Fig. 7 Polarization angle χ as function of the square of the restframe wavelength λ_{0}. Polarization angles are restricted to the interval [0°,180°] ; for 0954+658 and 1637+574 these intervals have been shifted to resolve 0°/180° ambiguities. The selection of data is limited to the epoch ranges quoted on top of each plot to maximize the significance of a rotation measure signal. We quote correlation coefficients r and significance levels s like in Fig. 3. For 0415+379, 0507+179, and 0954+658, we derive 3σ upper limits on RM. For 1418+546 and 1637+574, we quote the rotation measures together with their statistical and systematic errors. Rotation measures are derived via linear fits (gray lines) to the data. 

Open with DEXTER  
In the text 
Fig. 8 Constraints on the values and the variability of the shock parameter δ. The shock parameter δ relates to the factor μ = δ/(2 − δ) by which linear polarization is reduced by shock compression. Our observations provide Δμ = 2(δ_{1} − δ_{2})/(2 − δ_{1})(2 − δ_{2}) ≈ σ_{m}/ ⟨ m_{L} ⟩ for two arbitrary points in time “1”, “2”. The grayshaded area marks the δ_{1} − δ_{2} plane range consistent with the observed σ_{m}/ ⟨ m_{L} ⟩ ≈ 0.33 − 0.65 for our sample of sources. 

Open with DEXTER  
In the text 