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.

The third Fermi-LAT catalog (3FGL) presented the data of the first four years of observations from the Fermi Gamma-ray Space Telescope mission. There are 3034 sources, 1010 of which still remain unidentified. Identifying and classifying gamma-ray emitters is of high significance with regard to studying high-energy astrophysics. We demonstrate that optical polarimetry can be an advantageous and practical tool in the hunt for counterparts of the unidentified gamma-ray sources (UGSs). Using data from the RoboPol project, we validated that a significant fraction of active galactic nuclei (AGN) associated with 3FGL sources can be identified due to their high optical polarization exceeding that of the field stars. We performed an optical polarimetric survey within $3\sigma$ uncertainties of four unidentified 3FGL sources. We discovered a previously unknown extragalactic object within the positional uncertainty of 3FGL J0221.2+2518. We obtained its spectrum and measured a redshift of $z=0.0609\pm0.0004$. Using these measurements and archival data we demonstrate that this source is a candidate counterpart for 3FGL J0221.2+2518 and most probably is a composite object: a star-forming galaxy accompanied by AGN. We conclude that polarimetry can be a powerful asset in the search for AGN candidate counterparts for unidentified Fermi sources. Future extensive polarimetric surveys at high galactic latitudes (e.g., PASIPHAE) will allow the association of a significant fraction of currently unidentified gamma-ray sources.

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 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 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 the first results of an [O III] 5007 A interference filter survey for Planetary Nebulae (PNe) in the Galactic bulge. Covering (at first) the 66 per cent of the survey area, we detected a total of 90 objects, including 25 new PNe, 57 known PNe and 8 known PNe candidates. Deep Halpha+[N II] CCD images have been obtained as well as low resolution spectra for the newly discovered PNe. Their spectral signature suggests that the detected emission originates from a photoionized nebula. In addition, absolute line fluxes have been measured and the electron densities are given. Accurate optical positions and optical diameters were also determined.

Deep Halpha+[N II] CCD images have been obtained in the area of the Pegasus Constellation. The resulting mosaic covers an extent of ~7.5 degr x ~8.5 degr and filamentary and diffuse emission is discovered. Several long filaments (up to ~1 degr) are found within the field, while diffuse emission is present mainly in the central and northern areas. The filaments show variations in the intensity along their extent suggesting inhomogeneous interstellar clouds. Faint soft X-ray emission is also detected in the ROSAT All-Sky Survey. It is mainly concentrated in the central areas of our field and overlaps the optical emission. The low ionization images of [S II] of selected areas mainly show faint diffuse emission, while in the medium ionization images of [O III] diffuse and faint filamentary structures are present. Spectrophotometric observations have been performed on the brightest filaments and indicate emission from photoionized or shock-heated gas. The sulfur line ratios indicate electron densities below ~600 cm^{-3}, while the absolute Halpha emission lies in the range of 1.1 - 8.8 x 10^{-17} erg s^{-1} cm^{-2} arcsec^{-2}. The detected optical line emission could be part of a single or multiple supernova explosions.

We performed deep optical observations of the area of the new supernova remnant G 69.4+1.2 in the emission lines of [O III], Halpha+[N II] and [S II]. The low ionization images reveal diffuse and filamentary emission in the central and south, south-west areas of our field. Estimates of the [S II]/Halpha ratio suggest that the detected emission in these areas originates from shock heated gas, while the strong extended source in the north must be an HII region. The medium ionization image of [O III] shows a single filament close to the field center. Emission from [O III] is not detected elsewhere in the field but only in the north from LBN 069.96+01.35. Deep long-slit spectra taken at the position of the [O III] filament suggest shock velocities ~120 km/s, while in other areas velocities around 50 km/s are expected. The sulfur lines ratio indicates electron densities less than 120 cm^{-3}. The absolute Halpha flux is ~5 x 10^{-17} erg s^{-1} cm^{-2} arcsec^{-2}. The optical emission is very well correlated with the radio emission, especially in the south west. The soft X-ray emission detected in the ROSAT All-Sky survey shows a satisfactory degree of correlation with the optical data in the south-west suggesting their association.

Deep optical CCD images of the supernova remnant G 17.4-2.3 were obtained and faint emission has been discovered. The images, taken in the emission lines of Halpha+[N II], [S II] and [O III], reveal filamentary structures in the east, south-east area, while diffuse emission in the south and central regions of the remnant is also present. The radio emission in the same area is found to be well correlated with the brightest optical filament. The flux calibrated images suggest that the optical filamentary emission originates from shock-heated gas ([S II]/Halpha > 0.4), while the diffuse emission seems to originate from an HII region ([S II]/Halpha < 0.3). Furthermore, deep long-slit spectra were taken at the bright [O III] filament and clearly show that the emission originates from shock heated gas. The [O III] flux suggests shock velocities into the interstellar "clouds" greater than 100 km/s, while the [S II] 6716/6731 ratio indicates electron densities ~240 cm^{-3}. Finally, the Halpha emission has been measured to be between 7 to 20 x 10^{-17} erg s^{-1} cm^{-2} arcsec^{-2}.

We present the first calibrated CCD images of the faint supernova remnant G 114.3+0.3 in the emission lines of [OII], [OIII], Halpha+[NII] and [SII]. The deep low ionization CCD images reveal diffuse emission in the south and central areas of the remnant. These are correlated with areas of intense radio emission, while estimates of the [SII]/Halpha ratio suggest that the detected emission originates from shock heated gas. In the medium ionization image of [OIII] we discovered a thin filament in the south matching very well the outer radio contours. This filament is not continuous over its total extent but shows variations in the intensity, mainly in the south-west suggesting inhomogeneous interstellar clouds. Deep long-slit spectra were also taken along the [OIII] filament clearly identifying the observed emission as emission from shock heated gas. The Halpha emission is a few times 10^(-17) erg s^(-1) cm^(-2) arcsec^(-2), while the variations seen in the [OIII] flux suggest shock velocities into the interstellar clouds around or below 100 km/s. The sulfur line ratio approaches the low density limit implying electron densities less than ~500 cm^(-3).

We present the first CCD mosaic ever produced for the supernova remnant G 65.3+5.7 in [OII] 3727 at a moderate angular resolution. The remnant was observed with the 0.3 m wide-field telescope at Skinakas Observatory, Crete, Greece. Low resolution spectroscopy was performed at selected areas around this extended remnant. The spectral data are under analysis, while a first processing of the imaging observations produced impressive results.

We present results from a first analysis of optical data on G 69.4+2.4 obtained with the 0.3 m and 1.3 m telescopes at Skinakas Observatory. Filamentary emission is detected in the [OIII] 5007 A and Halpha+[NII] images oriented in the south-west, north-east direction. The Halpha+[NII] images reveal a complex network of filaments, while the [OIII] images show a single filament. Low resolution spectra were also obtained which identify the observed emission as emission from shock heated gas.

Wide field CCD observations of three, less known, supernova remnants. Filamentary and diffuse emission is discovered from the supernova remnant G67.7+1.8. Emission can also be seen in the ROSAT All Sky Survey data which indicate an extended hard X-ray source. The optical radiation from the G31.5-0.6 area shows weak sulfur emission and is probably not related to the remnant. Weak sulfur emission is also observed from W51C or G49.2-0.7. However, a small patch of emission in the south-east seems to emit strong sulfur emission. Results from long slit spectra are also presented.

Deep CCD exposures of the peculiar supernova remnant CTB 80 in the light of major optical lines have been obtained. These images reveal significant shock heated emission in the area of the remnant. The sulfur line image shows emission in the north along the outer boundary of the IRAS and HI shells. The comparison between the [OIII] and [OII] line images further suggest the presence of significant inhomogeneities in the interstellar medium. The flux calibrated images do not indicate the presence of incomplete recombination zones, and we estimate that the densities of the preshock clouds should not exceed a few atoms per cm^3. The area covered by the optical radiation along with the radio emission at 1410 MHz suggest that CTB 80 occupies a larger angular extent than was previously known.

A lightcurve of the eclipsing binary CM Draconis has been analyzed for the presence of transits of planets of size >= 2.5 Earth-radii (Re), with periods of 60 days or less, and in co-planar orbits around the binary system. About 400 million model lightcurves, representing transits from planets with periods ranging from 7 to 60 days, have been matched/correlated against these data. This process we call the "transit detection algorithm" or TDA. The resulting `transit-statistics' for each planet candidate allow the quantification of detection probabilities, and of false alarm rates. Our current lightcurve of CM Dra has a coverage of 1014 hours with 26,043 individual points, at a photometric precision between 0.2% and 0.7%. Planets significantly larger then 3Re would constitute a `supra-noise' detection, and for periods of 60 days or less, they would have been detected with a probability of 90%. `Subnoise' detections of smaller planets are more constrained. For example, 2.5 Re planets with 10-day periods or less would have been detected with an 80% probability. The necessity for predicted observations is illustrated with the nine top planet candidates that emerged from our TDA analysis. They are the planet candidates with the highest transit-statistics from the 1994-1998 observing seasons and, for them, transits for the 1999 observing season were predicted. Of the seven candidates that were then observationally tested in 1999, all were ruled out except one, which needs further observational confirmation. We conclude that the photometric transit method is a viable way to search for relatively small, inner extrasolar planets with moderate-sized telescopes using CCD photometry with a matching-filter analysis.

The opacity of spiral galaxies is examined by modelling the dust and stellar content of individual galaxies. The model is applied to five late-type spiral galaxies (NGC 4013, IC 2531, UGC 1082, NGC 5529 and NGC 5907). Having analyzed a total of seven galaxies thus far, the five galaxies mentioned above plus UGC 2048 and NGC 891 presented in (Xilouris et al. 1997, 1998), we are able to draw some general conclusions, the most significant of which are: 1) The face-on central optical depth is less than one in all optical bands indicating that typical spiral galaxies like the ones that we have modelled would be completely transparent if they were to be seen face-on. 2) The dust scaleheight is about half that of the stars, which means that the dust is more concentrated near the plane of the disk. 3) The dust scalelength is about 1.4 times larger than that of the stars and the dust is more radially extended than the stars. 4) The dust mass is found to be about an order of a magnitude more than previously measured using the IRAS fluxes, indicating the existence of a cold dust component. The gas-to-dust mass ratio calculated is close to the value derived for our Galaxy. 5) The derived extinction law matches quite well the Galactic extinction law, indicating a universal dust behaviour.