It has not been shown so far whether the diffuse Galactic polarized emission at frequencies relevant for cosmic microwave background (CMB) studies originates from nearby or more distant regions of our Galaxy. This questions previous attempts that have been made to constrain magnetic field models at local and large scales. The scope of this work is to investigate and quantify the contribution of the dusty and magnetized local interstellar medium to the observed emission that is polarized by thermal dust. We used stars as distance candles and probed the line-of-sight submillimeter polarization properties by comparing the emission that is polarized by thermal dust at submillimeter wavelengths and the optical polarization caused by starlight. We provide statistically robust evidence that at high Galactic latitudes ($|b| \geq 60^\circ$), the $353$ GHz polarized sky as observed by \textit{Planck} is dominated by a close-by magnetized structure that extends between $200$ and $300$ pc and coincides with the shell of the Local Bubble. Our result will assist modeling the magnetic field of the Local Bubble and characterizing the CMB Galactic foregrounds.

Galactic dust emission is polarized at unexpectedly high levels, as revealed by Planck. The origin of the observed $\simeq 20\%$ polarization fractions can be identified by characterizing the properties of optical starlight polarization in a region with maximally polarized dust emission. We measure the R-band linear polarization of 22 stars in a region with a submillimeter polarization fraction of $\simeq 20%$. A subset of 6 stars is also measured in the B, V and I bands to investigate the wavelength dependence of polarization. We find that starlight is polarized at correspondingly high levels. Through multiband polarimetry we find that the high polarization fractions are unlikely to arise from unusual dust properties, such as enhanced grain alignment. Instead, a favorable magnetic field geometry is the most likely explanation, and is supported by observational probes of the magnetic field morphology. The observed starlight polarization exceeds the classical upper limit of $\left[p_V/E\left(B-V\right)\right]_{\rm max} = 9$%mag$^{-1}$ and is at least as high as 13%mag$^{-1}$ that was inferred from a joint analysis of Planck data, starlight polarization and reddening measurements. Thus, we confirm that the intrinsic polarizing ability of dust grains at optical wavelengths has long been underestimated.

Open clusters (OCs) are popular tracers of the structure and evolutionary history of the Galactic disk. The OC population is often considered to be complete within 1.8 kpc of the Sun. The recent Gaia Data Release 2 (DR2) allows the latter claim to be challenged. We perform a systematic search for new OCs in the direction of Perseus using precise and accurate astrometry from Gaia DR2. We implement a coarse-to-fine search method. First, we exploit spatial proximity using a fast density-aware partitioning of the sky via a k-d tree in the spatial domain of Galactic coordinates, (l, b). Secondly, we employ a Gaussian mixture model in the proper motion space to quickly tag fields around OC candidates. Thirdly, we apply an unsupervised membership assignment method, UPMASK, to scrutinise the candidates. We visually inspect colour-magnitude diagrams to validate the detected objects. Finally, we perform a diagnostic to quantify the significance of each identified overdensity in proper motion and in parallax space We report the discovery of 41 new stellar clusters. This represents an increment of at least 20% of the previously known OC population in this volume of the Milky Way. We also report on the clear identification of NGC 886, an object previously considered an asterism. This letter challenges the previous claim of a near-complete sample of open clusters up to 1.8 kpc. Our results reveal that this claim requires revision, and a complete census of nearby open clusters is yet to be found.

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.

The availability of large datasets with stellar distance and polarization information will enable a tomographic reconstruction of the (plane-of-the-sky-projected) interstellar magnetic field in the near future. We demonstrate the feasibility of such a decomposition within a small region of the diffuse ISM. We combine measurements of starlight (R-band) linear polarization obtained using the RoboPol polarimeter with stellar distances from the second Gaia data release. The stellar sample is brighter than 17 mag in the R band and reaches out to several kpc from the Sun. HI emission spectra reveal the existence of two distinct clouds along the line of sight. We decompose the line-of-sight-integrated stellar polarizations to obtain the mean polarization properties of the two clouds. The two clouds exhibit significant differences in terms of column density and polarization properties. Their mean plane-of-the-sky magnetic field orientation differs by 60 degrees. We show how our tomographic decomposition can be used to constrain our estimates of the polarizing efficiency of the clouds as well as the frequency dependence of the polarization angle of polarized dust emission. We also demonstrate a new method to constrain cloud distances based on this decomposition. Our results represent a preview of the wealth of information that can be obtained from a tomographic map of the ISM magnetic field.

PASIPHAE (the Polar-Areas Stellar Imaging in Polarization High-Accuracy Experiment) is an optopolarimetric survey aiming to measure the linear polarization from millions of stars, and use these to create a three-dimensional tomographic map of the magnetic field threading dust clouds within the Milky Way. This map will provide invaluable information for future CMB B-mode experiments searching for inflationary gravitational waves, providing unique information regarding line-of-sight integration effects. Optical polarization observations of a large number of stars at known distances, tracing the same dust that emits polarized microwaves, can map the magnetic field between them. The Gaia mission is measuring distances to a billion stars, providing an opportunity to produce a tomographic map of Galactic magnetic field directions, using optical polarization of starlight. Such a map will not only boost CMB polarization foreground removal, but it will also have a profound impact in a wide range of astrophysical research, including interstellar medium physics, high-energy astrophysics, and evolution of the Galaxy. Taking advantage of the novel technology implemented in our high-accuracy Wide-Area Linear Optical Polarimeters (WALOPs) currently under construction at IUCAA, India, we will engage in a large-scale optopolarimetric program that can meet this challenge: a survey of both northern and southern Galactic polar regions targeted by CMB experiments, covering over 10,000 square degrees, which will measure linear optical polarization of over 360 stars per square degree (over 3.5 million stars, a 1000-fold increase over the state of the art). The survey will be conducted concurrently from the South African Astronomical Observatory in Sutherland, South Africa in the southern hemisphere, and the Skinakas Observatory in Crete, Greece, in the north.

We conduct a small-scale pathfinding survey designed to identify the average polarization properties of the diffuse ISM locally at the lowest dust content regions. We perform deep optopolarimetric surveys within three $\sim 15' \times 15'$ regions located at $b > 48^\circ$, using the RoboPol instrument. The observed samples of stars are photometrically complete to $\sim$16 mag in the R-band. The selected regions exhibit low dust emission at 353 GHz and low total reddening compared to the majority of high-latitude sightlines. We measure the level of systematic uncertainty for all observing epochs and find it to be 0.1\% in fractional linear polarization, $p$. The majority of individual stellar measurements are non-detections. However, our survey strategy enables us to locate the mean fractional linear polarization $p_{mean}$ in each of the three regions. The region with lowest dust content yields $p_{mean}=(0.054 \pm 0.038) \%$, not significantly different from zero. We find significant detections for the remaining two regions of: $p_{mean}=(0.113 \pm 0.036) \%$ and $p_{mean}=(0.208 \pm 0.044) \%$. Using a Bayesian approach we provide upper limits on the intrinsic spread of the small-scale distributions of $q$ and $u$. At the detected $p_{mean}$ levels, the determination of the systematic uncertainty is critical for the reliability of the measurements. We verify the significance of our detections with statistical tests, accounting for all sources of uncertainty. Using publicly available HI emission data, we identify the velocity components that most likely account for the observed $p_{mean}$ and find their morphologies to be misaligned with the orientation of the mean plane-of-sky magnetic field at a spatial resolution of 10$\arcmin$. We find indications that the standard upper envelope of $p$ with reddening underestimates the maximum $p$ at very low E(B-V) ($\leq 0.01$ mag).

Filaments in Herschel molecular cloud images are found to exhibit a "characteristic width". This finding is in tension with spatial power spectra of the data, which show no indication of this characteristic scale. We demonstrate that this discrepancy is a result of the methodology adopted for measuring filament widths. First, we perform the previously used analysis technique on artificial scale-free data, and obtain a peaked width distribution of filament-like structures. Next, we repeat the analysis on three Herschel maps and reproduce the narrow distribution of widths found in previous studies $-$ when considering the average width of each filament. However, the distribution of widths measured at all points along a filament spine is broader than the distribution of mean filament widths, indicating that the narrow spread (interpreted as a "characteristic" width) results from averaging. Furthermore, the width is found to vary significantly from one end of a filament to the other. Therefore, the previously identified peak at 0.1 pc cannot be understood as representing the typical width of filaments. We find an alternative explanation by modelling the observed width distribution as a truncated power-law distribution, sampled with uncertainties. The position of the peak is connected to the lower truncation scale and is likely set by the choice of parameters used in measuring filament widths. We conclude that a "characteristic" width of filaments is not supported by the available data.