Context.
The magnetic field strength in interstellar clouds can be estimated indirectly from measurements of dust polarization by assuming that turbulent kinetic energy is comparable to the ...fluctuating magnetic energy, and using the spread of polarization angles to estimate the latter. The method developed by Davis (1951, Phys. Rev., 81, 890) and by Chandrasekhar and Fermi (1953, ApJ, 118, 1137) (DCF) assumes that incompressible magnetohydrodynamic (MHD) fluctuations induce the observed dispersion of polarization angles, deriving
B
∝ 1∕
δθ
(or, equivalently,
δθ
∝
M
A
, in terms of the Alfvénic Mach number). However, observations show that the interstellar medium is highly compressible. Recently, two of us (ST) relaxed the incompressibility assumption and derived instead B ∝ 1/√δθ (equivalently,
δθ
∝
M
A
2
).
Aims.
We explored what the correct scaling is in compressible and magnetized turbulence through theoretical arguments, and tested the assumptions and the accuracy of the two methods with numerical simulations.
Methods.
We used 26 magnetized, isothermal, ideal-MHD numerical simulations without self-gravity and with different types of forcing. The range of
M
A
and sonic Mach numbers
M
s
explored are 0.1 ≤
M
A
≤ 2.0 and 0.5 ≤
M
s
≤ 20. We created synthetic polarization maps and tested the assumptions and accuracy of the two methods.
Results.
The synthetic data have a remarkable consistency with the δθ ∝ M
A
2
scaling, which is inferred by ST, while the DCF scaling failed to follow the data. Similarly, the assumption of ST that the turbulent kinetic energy is comparable to the root-mean-square (rms) of the coupling term of the magnetic energy between the mean and fluctuating magnetic field is valid within a factor of two for all
M
A
(with the exception of solenoidally driven simulations at high
M
A
, where the assumption fails by a factor of 10). In contrast, the assumption of DCF that the turbulent kinetic energy is comparable to the rms of the second-order fluctuating magnetic field term fails by factors of several to hundreds for sub-Alfvénic simulations. The ST method shows an accuracy better than 50% over the entire range of
M
A
explored; DCF performs adequately only in the range of
M
A
for which it has been optimized through the use of a “fudge factor”. For low
M
A
, it is inaccurate by factors of tens, since it omits the magnetic energy coupling term, which is of first order and corresponds to compressible modes. We found no dependence of the accuracy of the two methods on
M
s
.
Conclusions.
The assumptions of the ST method reflect better the physical reality in clouds with compressible and magnetized turbulence, and for this reason the method provides a much better estimate of the magnetic field strength over the DCF method. Even in super-Alfvénic cases where DCF might outperform ST, the ST method still provides an adequate estimate of the magnetic field strength, while the reverse is not true.
If a single line of sight (LOS) intercepts multiple dust clouds with different spectral energy distributions and magnetic field orientations, then the frequency scaling of each of the Stokes
Q
and
U
...parameters of the thermal dust emission may be different, a phenomenon we refer to as LOS frequency decorrelation. We present first evidence for LOS frequency decorrelation in
Planck
data using independent measurements of neutral-hydrogen (H
I
) emission to probe the 3D structure of the magnetized interstellar medium (ISM). We use H
I
-based measurements of the number of clouds per LOS and the magnetic field orientation in each cloud to select two sets of sightlines: (i) a target sample of pixels that are likely to exhibit LOS frequency decorrelation and (ii) a control sample of pixels that lack complex LOS structure. We test the null hypothesis that LOS frequency decorrelation is not detectable in
Planck
353 and 217 GHz polarization data at high Galactic latitudes. We reject the null hypothesis at high significance based on data that show that the combined effect of polarization angle variation with frequency and depolarization are detected in the target sample. This detection is robust against the choice of cosmic microwave background (CMB) map and map-making pipeline. The observed change in polarization angle due to LOS frequency decorrelation is detectable above the
Planck
noise level. The probability that the detected effect is due to noise alone ranges from 5 × 10
−2
to 4 × 10
−7
, depending on the CMB subtraction algorithm and treatment of residual systematic errors; correcting for residual systematic errors consistently increases the significance of the effect. Within the target sample, the LOS decorrelation effect is stronger for sightlines with more misaligned magnetic fields, as expected. With our sample, we estimate that an intrinsic variation of ~15% in the ratio of 353 to 217 GHz polarized emission between clouds is sufficient to reproduce the measured effect. Our finding underlines the importance of ongoing studies to map the three-dimensional structure of the magnetized and dusty ISM that could ultimately help component separation methods to account for frequency decorrelation effects in CMB polarization studies.
In a companion paper we presented a statistical model for the blazar population, consisting of distributions for the unbeamed radio luminosity function and the Lorentz factor distribution of each of ...the BL Lac and Flat Spectrum Radio Quasar (FSRQ) classes. Our model has been optimized so that it reproduces the MOJAVE (Monitoring Of Jets in Active galactic nuclei with VLBA Experiments) distributions of apparent speeds and redshifts when the appropriate flux limit is applied and a uniform distribution of jet viewing angles is assumed for the population. Here we use this model to predict the Doppler factor distribution for various flux-limited samples for which Doppler factors have been estimated in a variety of ways (equipartition, variability + equipartition, inverse Compton dominance) on a blazar-by-blazar basis. By comparing the simulated and data-estimated Doppler factor distributions in each case, we evaluate the different methods of estimating blazar Doppler factors. We find that the variability Doppler factors assuming equipartition are the ones in the best agreement with our statistical model, whereas the inverse Compton Doppler factor method is only suitable for FSRQs. In the case of variability Doppler factors, we find that while random errors are relatively low (∼30 per cent), uncertainties are dominated by systematic effects. In the case of inverse Compton Doppler factors, random errors appear to dominate, but are significantly larger (∼60 per cent).
We present a practical implementation of a Monte Carlo method to estimate the significance of cross-correlations in unevenly sampled time series of data, whose statistical properties are modelled ...with a simple power-law power spectral density. This implementation builds on published methods; we introduce a number of improvements in the normalization of the cross-correlation function estimate and a bootstrap method for estimating the significance of the cross-correlations. A closely related matter is the estimation of a model for the light curves, which is critical for the significance estimates. We present a graphical and quantitative demonstration that uses simulations to show how common it is to get high cross-correlations for unrelated light curves with steep power spectral densities. This demonstration highlights the dangers of interpreting them as signs of a physical connection. We show that by using interpolation and the Hanning sampling window function we are able to reduce the effects of red-noise leakage and to recover steep simple power-law power spectral densities. We also introduce the use of a Neyman construction for the estimation of the errors in the power-law index of the power spectral density. This method provides a consistent way to estimate the significance of cross-correlations in unevenly sampled time series of data.
RoboPol: AGN polarimetric monitoring data Blinov, D; Kiehlmann, S; Pavlidou, V ...
Monthly notices of the Royal Astronomical Society,
03/2021, Letnik:
501, Številka:
3
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
We present uniformly reprocessed and re-calibrated data from the RoboPol programme of optopolarimetric monitoring of active galactic nuclei (AGNs), covering observations between 2013, when ...the instrument was commissioned, and 2017. In total, the data set presented in this paper includes 5068 observations of 222 AGN with Dec. > −25○. We describe the current version of the RoboPol pipeline that was used to process and calibrate the entire data set, and we make the data publicly available for use by the astronomical community. Average quantities summarizing optopolarimetric behaviour (average degree of polarization, polarization variability index) are also provided for each source we have observed and for the time interval we have followed it.
Magnetohydrodynamic (MHD) turbulence is a cross-field process relevant to many systems. A prerequisite for understanding these systems is to constrain the role of MHD turbulence, and in particular, ...the energy exchange between kinetic and magnetic forms. The energetics of strongly magnetized and compressible turbulence has so far resisted attempts to understand them. Numerical simulations reveal that kinetic energy can be orders of magnitude higher than fluctuating magnetic energy. We solved this lack-of-balance puzzle by calculating the energetics of compressible and sub-Alfvénic turbulence based on the dynamics of coherent cylindrical fluid parcels. Using the MHD Lagrangian, we proved analytically that the bulk of the magnetic energy transferred to kinetic energy is the energy that is stored in the coupling between the ordered and fluctuating magnetic field. The analytical relations are in strikingly good agreement with numerical data, up to second-order terms.
Abstract
We use results of our 3 yr polarimetric monitoring programme to investigate the previously suggested connection between rotations of the polarization plane in the optical emission of blazars ...and their gamma-ray flares in the GeV band. The homogeneous set of 40 rotation events in 24 sources detected by RoboPol is analysed together with the gamma-ray data provided by Fermi-LAT. We confirm that polarization plane rotations are indeed related to the closest gamma-ray flares in blazars and the time lags between these events are consistent with zero. Amplitudes of the rotations are anticorrelated with amplitudes of the gamma-ray flares. This is presumably caused by higher relativistic boosting (higher Doppler factors) in blazars that exhibit smaller amplitude polarization plane rotations. Moreover, the time-scales of rotations and flares are marginally correlated.
Context. We investigate the optical polarization properties of high-energy BL Lac objects using data from the RoboPol blazar monitoring program and the Nordic Optical Telescope. Aims. We wish to ...understand if there are differences between the BL Lac objects that have been detected with the current-generation TeV instruments and those objects that have not yet been detected. Methods. We used a maximum-likelihood method to investigate the optical polarization fraction and its variability in these sources. In order to study the polarization position angle variability, we calculated the time derivative of the electric vector position angle (EVPA) change. We also studied the spread in the Stokes Q/I−U/I plane and rotations in the polarization plane. Results. The mean polarization fraction of the TeV-detected BL Lacs is 5%, while the non-TeV sources show a higher mean polarization fraction of 7%. This difference in polarization fraction disappears when the dilution by the unpolarized light of the host galaxy is accounted for. The TeV sources show somewhat lower fractional polarization variability amplitudes than the non-TeV sources. Also the fraction of sources with a smaller spread in the Q/I−U/I plane and a clumped distribution of points away from the origin, possibly indicating a preferred polarization angle, is larger in the TeV than in the non-TeV sources. These differences between TeV and non-TeV samples seem to arise from differences between intermediate and high spectral peaking sources instead of the TeV detection. When the EVPA variations are studied, the rate of EVPA change is similar in both samples. We detect significant EVPA rotations in both TeV and non-TeV sources, showing that rotations can occur in high spectral peaking BL Lac objects when the monitoring cadence is dense enough. Our simulations show that we cannot exclude a random walk origin for these rotations. Conclusions. These results indicate that there are no intrinsic differences in the polarization properties of the TeV-detected and non-TeV-detected high-energy BL Lac objects. This suggests that the polarization properties are not directly related to the TeV-detection, but instead the TeV loudness is connected to the general flaring activity, redshift, and the synchrotron peak location.
In order to determine the location of the gamma-ray emission site in blazars, we investigate the time-domain relationship between their radio and gamma-ray emission. Light curves for the brightest ...detected blazars from the first 3 yr of the mission of the Fermi Gamma-ray Space Telescope are cross-correlated with 4 yr of 15 GHz observations from the Owens Valley Radio Observatory 40 m monitoring programme. The large sample and long light-curve duration enable us to carry out a statistically robust analysis of the significance of the cross-correlations, which is investigated using Monte Carlo simulations including the uneven sampling and noise properties of the light curves. Modelling the light curves as red noise processes with power-law power spectral densities, we find that only one of 41 sources with high-quality data in both bands shows correlations with significance larger than 3σ (AO 0235+164), with only two more larger than even 2.25σ (PKS 1502+106 and B2 2308+34). Additionally, we find correlated variability in Mrk 421 when including a strong flare that occurred in 2012 July–September. These results demonstrate very clearly the difficulty of measuring statistically robust multiwavelength correlations and the care needed when comparing light curves even when many years of data are used. This should be a caution. In all four sources, the radio variations lag the gamma-ray variations, suggesting that the gamma-ray emission originates upstream of the radio emission. Continuous simultaneous monitoring over a longer time period is required to obtain high significance levels in cross-correlations between gamma-ray and radio variability in most blazars.