Theoretical simulations and observations at different angular resolutions have shown that magnetic fields have a central role in massive star formation. Like in low-mass star formation, the magnetic field in massive young stellar objects can either be oriented along the outflow axis or randomly. Measuring the magnetic field at milliarcsecond resolution (10-100 au) around a substantial number of massive young stellar objects permits determining with a high statistical significance whether the direction of the magnetic field is correlated with the orientation of the outflow axis or not. In late 2012, we started a large VLBI campaign with the European VLBI Network to measure the linearly and circularly polarized emission of 6.7 GHz methanol masers around a sample of massive star-forming regions. This paper focuses on the first seven observed sources, G24.78+0.08, G25.65+1.05, G29.86-0.04, G35.03+0.35, G37.43+1.51, G174.20-0.08, and G213.70-12.6. For all these sources, molecular outflows have been detected in the past. We detected a total of 176 methanol masing cloudlets toward the seven massive star-forming regions, 19% of which show linearly polarized emission. The methanol masers around the massive young stellar object MM1 in G174.20-0.08 show neither linearly nor circularly polarized emission. The linear polarization vectors are well ordered in all the other massive young stellar objects. We measured significant Zeeman splitting toward both A1 and A2 in G24.78+0.08, and toward G29.86-0.04 and G213.70-12.6. By considering all the 19 massive young stellar objects reported in the literature for which both the orientation of the magnetic field at milliarcsecond resolution and the orientation of outflow axes are known, we find evidence that the magnetic field (on scales 10-100 au) is preferentially oriented along the outflow axes.

In this work we present a polarization analysis at radio frequencies of Markarian 421, one of the closest (z=0.03) TeV blazars. The observations were obtained, both in total and in polarized intensity, with the Very Long Baseline Array (VLBA) at 15, 24, and 43 GHz throughout 2011, with one observation per month (for a total of twelve epochs). We investigate the magnetic field topology and the polarization structure on parsec scale and their evolution with time. We detect polarized emission both in the core and in the jet region, and it varies with frequency, location and time. In the core region we measure a mean fractional polarization of about 1-2%, with a peak of about 4% in March at 43 GHz; the polarization angle is almost stable at 43 GHz, but it shows significant variability in the range 114-173 deg at 15 GHz. In the jet region the polarization properties show a more stable behavior; the fractional polarization is about 16% and the polarization angle is nearly perpendicular to the jet axis. The higher EVPA variability observed at 15 GHz is due both to a variable Faraday rotation effect and to opacity. The residual variability observed in the intrinsic polarization angle, together with the low degree of polarization in the core region, could be explained with the presence of a blend of variable cross-polarized subcomponents within the beam.

We analyzed data in polarized intensity obtained with the Very Long Baseline Array (VLBA) at twelve epochs (one observation per month from January to December 2011) at 15, 24, and 43 GHz. For the absolute orientation of the electric vector position angles (EVPA) we used the D-terms method. We also used gamma-ray data from the Fermi Large Area Telescope on weekly time bins throughout 2011. The source shows polarized emission, and its properties vary with time, frequency, and location along the jet. The core mean polarization fraction is generally between 1% and 2%, with a 4% peak at 43 GHz in March; the polarization angle is variable, mainly at 15 GHz, where it changes frequently, and less so at 43 GHz, where it oscillates in the range 114 - 173 deg. The jet polarization properties are more stable, with a fractional polarization of around 16% and a polarization angle nearly perpendicular to the jet axis. The average flux and photon index at gamma-ray energies are (17.7+-0.5)x10^8 ph cm^-2 s^-1 and the photon index is 1.77+-0.02. The gamma-ray light curve shows variability, with a main peak that appears to be associated with the peak in the core polarized emission at 43 GHz, as well as with the total intensity light curve. A discrete correlation function analysis yields a correlation coefficient of 0.54 at zero delay, with a significance level above 99.7%. We accurately determine the polarization properties of Mrk 421, both in the core and in the jet region. The radio and gamma-ray light curves are correlated. The observed EVPA variability at 15 GHz is partly due to opacity and partly to a variable Faraday rotation effect. To explain the residual variability of the intrinsic polarization angle and the low degree of polarization in the core region, we invoke a blend of variable cross-polarized subcomponents with different polarization properties within the beam.

We report the results obtained for the AGN Markarian 421 by model-fitting the data in the visibility plane, studing the proper motion of jet components, the light curve, and the spectral index of the jet features. We compare the radio data with optical light curves obtained at the Steward Observatory, considering also the optical polarization information. Mrk 421 has a bright nucleus and a one-sided jet extending towards the north-west for a few parsecs. The model-fits show that brightness distribution is well described using 6-7 circular Gaussian components, four of which are reliably identified at all epochs; all components are effectively stationary except for component D, at ~0.4 mas from the core, whose motion is however subluminal. Analysis of the light curve shows two different states, with the source being brighter and more variable in the first half of 2011 than in the second half. The highest flux density is reached in February. A comparison with the optical data reveals an increase of the V magnitude and of the fractional polarization simultaneous with the enhancement of the radio activity.