We present the design and performance of RoboPol, a four-channel optical polarimeter operating at the Skinakas Observatory in Crete, Greece. RoboPol is capable of measuring both relative linear Stokes parameters $q$ and $u$ (and the total intensity $I$) in one sky exposure. Though primarily used to measure the polarization of point sources in the R-band, the instrument features additional filters (B, V and I), enabling multi-wavelength imaging polarimetry over a large field of view (13.6' $\times$ 13.6'). We demonstrate the accuracy and stability of the instrument throughout its five years of operation. Best performance is achieved within the central region of the field of view and in the R band. For such measurements the systematic uncertainty is below 0.1% in fractional linear polarization, $p$ (0.05% maximum likelihood). Throughout all observing seasons the instrumental polarization varies within 0.1% in $p$ and within 1$^\circ$ in polarization angle.

We use results of our 3 year polarimetric monitoring program 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 {\em RoboPol} is analysed together with the gamma-ray data provided by {\em 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.

We present average R-band optopolarimetric data, as well as variability parameters, from the first and second RoboPol observing season. We investigate whether gamma- ray--loud and gamma-ray--quiet blazars exhibit systematic differences in their optical polarization properties. We find that gamma-ray--loud blazars have a systematically higher polarization fraction (0.092) than gamma-ray--quiet blazars (0.031), with the hypothesis of the two samples being drawn from the same distribution of polarization fractions being rejected at the 3{\sigma} level. We have not found any evidence that this discrepancy is related to differences in the redshift distribution, rest-frame R-band lu- minosity density, or the source classification. The median polarization fraction versus synchrotron-peak-frequency plot shows an envelope implying that high synchrotron- peaked sources have a smaller range of median polarization fractions concentrated around lower values. Our gamma-ray--quiet sources show similar median polarization fractions although they are all low synchrotron-peaked. We also find that the random- ness of the polarization angle depends on the synchrotron peak frequency. For high synchrotron-peaked sources it tends to concentrate around preferred directions while for low synchrotron-peaked sources it is more variable and less likely to have a pre- ferred direction. We propose a scenario which mediates efficient particle acceleration in shocks and increases the helical B-field component immediately downstream of the shock.

We present a new set of optical polarization plane rotations in blazars, observed during the third year of operation of RoboPol. The entire set of rotation events discovered during three years of observations is analysed with the aim of determining whether these events are inherent in all blazars. It is found that the frequency of the polarization plane rotations varies widely among blazars. This variation cannot be explained either by a difference in the relativistic boosting or by selection effects caused by a difference in the average fractional polarization. We conclude that the rotations are characteristic of a subset of blazars and that they occur as a consequence of their intrinsic properties.

The stages before the formation of stars in molecular clouds are poorly understood. Insights can be gained by studying the properties of quiescent clouds, such as their magnetic field structure. The plane-of-the-sky orientation of the field can be traced by polarized starlight. We present the first extended, wide-field ($\sim$10 $\rm deg^2$) map of the Polaris Flare cloud in dust-absorption induced optical polarization of background stars, using the RoboPol polarimeter at the Skinakas Observatory. This is the first application of the wide-field imaging capabilities of RoboPol. The data were taken in the R-band and analysed with the automated reduction pipeline of the instrument. We present in detail optimizations in the reduction pipeline specific to wide-field observations. Our analysis resulted in reliable measurements of 641 stars with median fractional linear polarization 1.3%. The projected magnetic field shows a large scale ordered pattern. At high longitudes it appears to align with faint striations seen in the Herschel-SPIRE map of dust emission (250 $\mu m$), while in the central 4-5 deg$^2$ it shows an eddy-like feature. The overall polarization pattern we obtain is in good agreement with large scale measurements by Planck of the dust emission polarization in the same area of the sky.

We present first results on polarization swings in optical emission of blazars obtained by RoboPol, a monitoring program of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events. A possible connection of polarization swing events with periods of high activity in gamma rays is investigated using the dataset obtained during the first season of operation. It was found that the brightest gamma-ray flares tend to be located closer in time to rotation events, which may be an indication of two separate mechanisms responsible for the rotations. Blazars with detected rotations have significantly larger amplitude and faster variations of polarization angle in optical than blazars without rotations. Our simulations show that the full set of observed rotations is not a likely outcome (probability $\le 1.5 \times 10^{-2}$) of a random walk of the polarization vector simulated by a multicell model. Furthermore, it is highly unlikely ($\sim 5 \times 10^{-5}$) that none of our rotations is physically connected with an increase in gamma-ray activity.

We present measurements of rotations of the optical polarization of blazars during the second year of operation of RoboPol, a monitoring programme of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events, and we analyse the large set of rotation events discovered in two years of observation. We investigate patterns of variability in the polarization parameters and total flux density during the rotation events and compare them to the behaviour in a non-rotating state. We have searched for possible correlations between average parameters of the polarization-plane rotations and average parameters of polarization, with the following results: (1) there is no statistical association of the rotations with contemporaneous optical flares; (2) the average fractional polarization during the rotations tends to be lower than that in a non-rotating state; (3) the average fractional polarization during rotations is correlated with the rotation rate of the polarization plane in the jet rest frame; (4) it is likely that distributions of amplitudes and durations of the rotations have physical upper bounds, so arbitrarily long rotations are not realised in nature.

We investigate the effectiveness of blind surveys for radio sources and galaxy cluster thermal Sunyaev-Zel'dovich effects (TSZEs) using the four-pair, beam-switched OCRA-f radiometer on the 32-m radio telescope in Poland. The predictions are based on mock maps that include the cosmic microwave background, TSZEs from hydrodynamical simulations, and unresolved radio sources. We estimate the effects of source clustering towards galaxy clusters from NVSS source counts around Planck-selected cluster candidates, and include appropriate correlations in our mock maps. The study allows us to quantify the effects of halo line-of-sight alignments, source confusion, and telescope angular resolution on the detections of TSZEs. We perform a similar analysis for the planned 100-m Hevelius radio telescope (RTH) equipped with a 49-beam radio camera. We find that RT32/OCRA-f will be suitable for small-field blind radio source surveys, and will detect $33^{+17}_{-11}$ new radio sources brighter than 0.87 mJy at 30~GHz in a 1 deg$^2$ field at $>5\sigma$ CL during a one-year, non-continuous, observing campaign, taking account of Polish weather conditions. It is unlikely that any galaxy cluster will be detected at $3\sigma$ CL in such a survey. A $60$-deg$^2$ survey, with field coverage of $2^2$ beams per pixel, at 15 GHz with the RTH, would find <1.5 galaxy clusters per year brighter than 60 $\mu$Jy (at $3\sigma$ CL), and would detect about $3.4 \times 10^4$ point sources brighter than 1 mJy at $5\sigma$ CL, with confusion causing flux density errors $\lesssim 2\%$ (20%) in 68% (95%) of the detected sources. A primary goal of the planned RTH will be a wide-area ($\pi$~sr) radio source survey at 15 GHz. This survey will detect nearly $3 \times 10^5$ radio sources at $5\sigma$ CL down to 1.3 mJy, and tens of galaxy clusters, in one year of operation with typical weather conditions. ABRIDGED

We report the polarized optical light curve of a gamma-ray burst afterglow obtained using the RoboPol instrument. Observations began 655 seconds after the initial burst of gamma-rays from GRB131030A, and continued uninterrupted for 2 hours. The afterglow displayed a low, constant fractional linear polarization of $p = (2.1 \pm 1.6)\,\%$ throughout, which is similar to the interstellar polarization measured on nearby stars. The optical brightness decay is consistent with a forward-shock propagating in a medium of constant density, and the low polarization fraction indicates a disordered magnetic field in the shock front. This supports the idea that the magnetic field is amplified by plasma instabilities on the shock front. These plasma instabilities produce strong magnetic fields with random directions on scales much smaller than the total observable region of the shock, and the resulting randomly-oriented polarization vectors sum to produce a low net polarization over the total observable region of the shock.

We present first results from RoboPol, a novel-design optical polarimeter operating at the Skinakas Observatory in Crete. The data, taken during the May - June 2013 commissioning of the instrument, constitute a single-epoch linear polarization survey of a sample of gamma-ray - loud blazars, defined according to unbiased and objective selection criteria, easily reproducible in simulations, as well as a comparison sample of, otherwise similar, gamma-ray - quiet blazars. As such, the results of this survey are appropriate for both phenomenological population studies and for tests of theoretical population models. We have measured polarization fractions as low as $0.015$ down to $R$ magnitude of 17 and as low as $0.035$ down to 18 magnitude. The hypothesis that the polarization fractions of gamma-ray - loud and gamma-ray - quiet blazars are drawn from the same distribution is rejected at the $10^{-3}$ level. We therefore conclude that gamma-ray - loud and gamma-ray - quiet sources have different optical polarization properties. This is the first time this statistical difference is demonstrated in optical wavelengths. The polarization fraction distributions of both samples are well-described by exponential distributions with averages of $\langle p \rangle =6.4 ^{+0.9}_{-0.8}\times 10^{-2}$ for gamma-ray--loud blazars, and $\langle p \rangle =3.2 ^{+2.0}_{-1.1}\times 10^{-2}$ for gamma-ray--quiet blazars. The most probable value for the difference of the means is $3.4^{+1.5}_{-2.0}\times 10^{-2}$. The distribution of polarization angles is statistically consistent with being uniform.

We describe the data reduction pipeline and control system for the RoboPol project. The RoboPol project is monitoring the optical $R$-band magnitude and linear polarization of a large sample of active galactic nuclei that is dominated by blazars. The pipeline calibrates and reduces each exposure frame, producing a measurement of the magnitude and linear polarization of every source in the $13'\times 13'$ field of view. The control system combines a dynamic scheduler, real-time data reduction, and telescope automation to allow high-efficiency unassisted observations.

We present 30 GHz Sunyaev Zel'dovich observations of a statistically complete sample of galaxy clusters with OCRA-p. The clusters are the 18 most X-ray luminous clusters at z > 0.2 in the ROSAT Brightest Cluster Sample. We correct for contaminant radio sources via supplementary observations with the Green Bank Telescope, also at 30 GHz, and remove a cluster that is contaminated by an unresolved X-ray source. All 17 remaining clusters have central SZ effects with Comptonisation parameter y_0 exceeding 1.9x10^-4, and 13 are detected at significance > 3 sigma. We use our data to examine scalings between y_0 and X-ray temperature, X-ray luminosity, and the X-ray mass proxy Y_X, and find good agreement with predictions from self-similar models of cluster formation, with an intrinsic scatter in y_0 of about 25%. We also comment on the success of the observations in the face of the contaminant source population, and the implications for upcoming cm-wave surveys.

Knowledge of the population of radio sources in the range ~2-200 GHz is important for understanding their effects on measurements of the Cosmic Microwave Background power spectrum. We report measurements of the 30 GHz flux densities of 605 radio sources from the Combined Radio All-sky Targeted Eight-GHz Survey (CRATES), which have been made with the One Centimetre Receiver Array prototype (OCRA-p) on the Torun 32-m telescope. The flux densities of sources that were also observed by WMAP and previous OCRA surveys are in broad agreement with those reported here, however a number of sources display intrinsic variability. We find a good correlation between the 30 GHz and Fermi gamma-ray flux densities for common sources. We examine the radio spectra of all observed sources and report a number of Gigahertz-peaked and inverted spectrum sources. These measurements will be useful for comparison to those from the Low Frequency Instrument of the Planck satellite, which will make some of its most sensitive observations in the region covered here.

Small angular scale (high l) studies of cosmic microwave background anisotropies require accurate knowledge of the statistical properties of extragalactic sources at cm-mm wavelengths. We have used a 30 GHz dual-beam receiver (OCRA-p) on the Torun 32-m telescope to measure the flux densities of 121 sources in VSA fields selected at 15 GHz with the Ryle Telescope. We have detected 57 sources above a limiting flux density of 5mJy, of which 31 sources have a flux density greater than 10mJy, which is our effective completeness limit. From these measurements we derive a surface density of sources above 10mJy at 30 GHz of 2.0+/-0.4 per square degree. This is consistent with the surface density obtained by Mason et al. (2009) who observed a large sample of sources selected at a much lower frequency (1.4 GHz). We have also investigated the dependence of the spectral index distribution on flux density by comparing our results with those for sources above 1 Jy selected from the WMAP 22 GHz catalogue. We conclude that the proportion of steep spectrum sources increases with decreasing flux density, qualitatively consistent with the predictions of deZotti et al. (2005). We find no evidence for an unexpected population of sources whose spectra rise towards high frequencies, which would affect our ability to interpret current high resolution CMB observations at 30 GHz and above.

We report the results of a survey of 442 planetary nebulae at 30 GHz. The purpose of the survey is to develop a list of planetary nebulae as calibration sources which could be used for high frequency calibration in future. For 41 PNe with sufficient data, we test the emission mechanisms in order to evaluate whether or not spinning dust plays an important role in their spectra at 30 GHz. The 30-GHz data were obtained with a twin-beam differencing radiometer, OCRA-p, which is in operation on the Torun 32-m telescope. Sources were scanned both in right ascension and declination. We estimated flux densities at 30 GHz using a free-free emission model and compared it with our data. The primary result is a catalogue containing the flux densities of 93 planetary nebulae at 30 GHz. Sources with sufficient data were compared with a spectral model of free-free emission. The model shows that free-free emission can generally explain the observed flux densities at 30 GHz thus no other emission mechanism is needed to account for the high frequency spectra.

Almost 50 radio SNR objects were selected from Green's catalogue for continuum and polarization observations. We present preliminary results of this observational campaign at 4.7 GHz carried out with Torun 32-metre radio telescope.

To measure the 30-GHz flux densities of the 293 sources in the Caltech-Jodrell Bank flat-spectrum (CJF) sample. The measurements are part of an ongoing programme to measure the spectral energy distributions of flat spectrum radio sources and to correlate them with the milliarcsecond structures from VLBI and other measured astrophysical properties. The 30-GHz data were obtained with a twin-beam differencing radiometer system mounted on the Torun 32-m telescope. The system has an angular resolution of 1.2 arcmin. Together with radio spectral data obtained from the literature, the 30-GHz data have enabled us to identify 42 of the CJF sources as Giga-hertz Peaked Spectrum (GPS) sources. Seventeen percent of the sources have rising spectra (alpha > 0) between 5 and 30 GHz.

We present 30 GHz Sunyaev Zel'dovich (SZ) observations of a sample of four galaxy clusters with a prototype of the One Centimetre Receiver Array (OCRA-p) which is mounted on the Torun 32-m telescope. The clusters (Cl0016+16, MS0451.6-0305, MS1054.4-0321 and Abell 2218) are popular SZ targets and serve as commissioning observations. All four are detected with clear significance (4-6 sigma) and values for the central temperature decrements are in good agreement with measurements reported in the literature. We believe that systematic effects are successfully suppressed by our observing strategy. The relatively short integration times required to obtain these results demonstrate the power of OCRA-p and its successors for future SZ studies.

We present multi-frequency observational results for the Compact Steep-Spectrum (CSS) quasar 3C309.1. The observations were carried out with the VSOP at 1.6 and 4.8GHz and the VLBA at 15.4GHz. The source has a distorted, one-sided radio jet. Relativistics effects and physical properties are discussed. Comparing the predicted and observed X-ray fluxes in the framework of the synchrotron self-Compton model we derive the beaming factor for 3C309.1. The complex structure of the jet may be explained by a relativistic helical flow in a heterogeneous, clumpy ISM.

The so-called HYbrid MOrphology Radio Sources (HYMORS) are a class of objects that appear to have a mixed Fanaroff-Riley (FR) morphology in a single object; i.e. a HYMORS has an FR I-type lobe on one side of its nucleus and an FR II-type lobe on the other side. Because of this unique feature and given that the origin of the FR morphological dichotomy is still unclear, HYMORS may possibly play a crucial role in our understanding of the FR-dichotomy. As the number of known HYMORS is quite small, we aimed to increase that number by inspecting a few areas of the sky covered by the VLA FIRST survey and by selecting 21 HYMORS candidates based on the morphology shown in the FIRST images. They were observed with the VLA in B-conf. at 4.9 GHz. Three objects from the initial sample turned out to be actual HYMORS and two others very likely to fulfill the criteria. These five were subsequently re-observed with the VLA in A-conf. at 1.4 GHz. Our results provide strong support to the findings of Gopal-Krishna & Wiita (2000), namely that there are two different kinds of jets in HYMORS; consequently, the existence of FR-dichotomy as a whole is difficult to reconcile with the class of explanations that posit fundamental differences in the central engine.

We report early results of an observational campaign targeted on a sample of compact steep spectrum sources selected from the FIRST survey which are significantly weaker than those investigated before. The selection criteria and procedure are given in detail. We present here an assortment of MERLIN and VLBI observations and make some general comments based on the morphologies of the sources presented.

Compact Steep Spectrum (CSS) sources are powerful extragalactic radio sources with angular dimensions of the order of a few arcseconds or less. Such a compactness is apparently linked to the youth of these objects. The majority of CSSs investigated so far have been known since the early 1980s. This paper is the first in a series where we report the results of an observational campaign targeted on a completely new sample of CSSs which are significantly weaker than those investigated before. The ultimate goal of that campaign is to find out how ``weak'' CSSs compare to ``strong'', classical ones, especially with regard to the morphologies. Here we present an analysis of morphological and physical properties of five relatively large sources based on MERLIN observations at 1.6 and 5 GHz.

We obtained EVN observations of the X-ray transient source GRS1915+105 in its radio-loud state. In 2 images obtained before a major flare, the source is resolved, showing two opposite asymmetric jets. This indicates that steady jets exist also before a major flare, with an estimated velocity of 0.2-0.6c.