Spherical collapse predicts that a single value of the turnaround density (average matter density within the scale on which a structure detaches from the Hubble flow) characterizes all cosmic structures at the same redshift. It has been recently shown by Korkidis et al. that this feature persists in complex non-spherical galaxy clusters identified in N-body simulations. Here we show that the low-redshift evolution of the turnaround density constrains the cosmological parameters, and that it can be used to derive a local constraint on $\Omega_\Lambda$ alone, independent of $\Omega_m$. The turnaround density thus provides a promising new way to exploit upcoming large cosmological datasets.

We use N-body simulations to examine whether a characteristic turnaround radius, as predicted from the spherical collapse model in a $\rm {\Lambda CDM}$ Universe, can be meaningfully identified for galaxy clusters, in the presence of full three-dimensional effects. We use The Dark Sky Simulations and Illustris-TNG dark-matter--only cosmological runs to calculate radial velocity profiles around collapsed structures, extending out to many times the virial radius $R_{200}$. There, the turnaround radius can be unambiguously identified as the largest non-expanding scale around a center of gravity. We find that: (a) Indeed, a single turnaround scale can meaningfully describe strongly non-spherical structures. (b) For halos of masses $M_{200}>10^{13}M_\odot$, the turnaround radius $R_{ta}$ scales with the enclosed mass $M_{ta}$ as $M_{ta}^{1/3}$, as predicted by the spherical collapse model. (c) The deviation of $R_{ta}$ in simulated halos from the spherical collapse model prediction is insensitive to halo asphericity. Rather, it is sensitive to the tidal forces due to massive neighbors when such are present. (d) Halos exhibit a characteristic average density within the turnaround scale. This characteristic density is dependent on cosmology and redshift. For the present cosmic epoch and for concordance cosmological parameters ($\Omega_m \sim 0.7$; $\Omega_\Lambda \sim 0.3$) turnaround structures exhibit an average matter density contrast with the background Universe of $\delta \sim 11$. Thus $R_{ta}$ is equivalent to $R_{11}$ -- in a way analogous to defining the "virial" radius as $R_{200}$ -- with the advantage that $R_{11}$ is shown in this work to correspond to a kinematically relevant scale in N-body simulations.

After three years of polarimetric monitoring of blazars, the RoboPol project has uncovered several key characteristics of polarimetric rotations in the optical for these most variable sources. The most important of these is that polarization properties of the synchrotron emission in the optical appear to be directly linked with gamma-ray activity. In this paper, we discuss the evidence for this connection, as well as the broader features of polarimetric behavior in blazars that are key in making progress with theoretical modeling of blazar emission.

The atomic interstellar medium (ISM) is observed to be full of linear structures that are referred to as "fibers". Fibers exhibit similar properties to linear structures found in molecular clouds known as striations. Suggestive of a similar formation mechanism, both striations and fibers appear to be ordered, quasi-periodic, and well-aligned with the magnetic field. The prevailing formation mechanism for striations involves the excitation of fast magnetosonic waves. Based on this theoretical model, and through a combination of velocity centroids and column density maps, Tritsis et al. (2018) developed a method for estimating the plane-of-sky (POS) magnetic field from molecular cloud striations. We apply this method in two H\textsc{I} clouds with fibers along the same line-of-sight (LOS) towards the ultra-high-energy cosmic-ray (UHECR) hotspot, at the boundaries of Ursa Major. For the cloud located closer to Earth, where Zeeman observations from the literature were also available, we find general agreement in the distributions of the LOS and POS components of the magnetic field. We find relatively large values for the total magnetic field (ranging from $\sim$$\rm{10}$ to $\sim$$\rm{20} ~\rm{\mu G}$) and an average projection angle with respect to the LOS of $\sim$ 50$^\circ$. For the cloud located further away, we also find a large value for the POS component of the magnetic field of $15^{+8}_{-3}~\rm{\mu G}$. We discuss the potential of our new magnetic-field tomography method for large-scale application. We consider the implications of our findings for the accuracy of current reconstructions of the Galactic magnetic field and on the propagation of UHECR through the ISM.

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Energetic particles in our own galaxy, galactic cosmic rays (GCRs), engage in a complex interplay with the interstellar medium and magnetic fields in the galaxy, giving rise to many of its key characteristics. This White Paper is the first of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. The focus of this white paper is galactic cosmic rays, supernova remnants, protostellar jets and superbubbles, and colliding wind binaries.

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Though they play a key role in cultivating the cosmological environment and/or enabling our studies of it, there is still much we do not know about AGNs and GRBs, particularly the avenue in which and through which they supply radiation and energetic particles, namely their jets. This White Paper is the second of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. The focus of this white paper is active galactic nuclei and gamma-ray bursts.

The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is measuring polarization. The capability to measure polarization is being realized currently at X-ray energies (e.g. with IXPE), and sensitive gamma-ray telescopes capable of measuring polarization, such as AMEGO, AdEPT, e-ASTROGAM, etc., are being developed. These future gamma-ray telescopes will probe the radiation mechanisms and magnetic fields of relativistic jets from active galactic nuclei at spatial scales much smaller than the angular resolution achieved with continuum observations of the instrument. In this white paper, we discuss the scientific potentials of high-energy polarimetry, especially gamma-ray polarimetry, including the theoretical implications, and observational technology advances being made. In particular, we will explore the primary scientific opportunities and wealth of information expected from synergy of multi-wavelength polarimetry that will be brought to multi-messenger astronomy.

Locating the sources of ultra-high energy cosmic rays (UHECRs) still remains a difficult puzzle for modern astrophysics. A major hurdle in the search for the sources is the fact that UHECRs are deflected by the Galactic magnetic field (GMF). Current knowledge of the GMF is limited, as most experimental measurements track line-of-sight--integrated quantities that are used to obtain best-fit parameters for global models including a large random component. The advent, however, of the Gaia era, with measurements of $\sim 10^9$ stellar parallaxes, in combination with upcoming large polarimetric surveys, make, for the first time, a 3D measurement of the GMF possible in principle. Such measurements can then be used to attempt a reconstruction of the trajectories of individual UHECRs through the Galaxy, in order to correct for their deflection. Motivated by these developments, in the present work, we study the limits of such a correction, by examining how its effectiveness depends on the uncertainty of any such future magnetic field measurements. To that end, we simulate attempts to reconstruct the trajectory of the cosmic ray by using hypothetical measurements of the GMF, based on values received from two recently updated GMF models. To simulate the uncertainty of a 3D measurement, random errors to these values are introduced separately for the plane-of-the-sky (POS) magnitude, the line-of-sight (LOS) magnitude and the POS direction. Our results highlight the conditions under which an effective correction is achievable. We find that the effectiveness is dependent on the particle rigidity and arrival direction, and can vary significantly depending on the GMF model used.

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.

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.

Despite the influx of unprecedented-quality data from the Fermi Gamma-Ray Space Telescope that have been collected over nine years of operation, the contribution of normal star-forming galaxies to the extragalactic gamma-ray background (EGRB) remains poorly constrained. Different estimates are discrepant both their underlying physical assumptions and their results. With several detections and many upper limits for the gamma-ray fluxes of nearby starforming galaxies now available, estimates that rely on empirical scalings between gamma-ray and longer-wavelength luminosities have become possible and increasingly popular. In this paper we examine factors that can bias such estimates, including: a) possible sources of nontrivial redshift dependence; b) dependence on the choice of star-formation tracer; c) uncertainties in the slope and normalisation of empirical scalings. We find that such biases can be significant, pointing towards the need for more sophisticated models for the star-forming galaxy contribution to the gamma-ray background, implementing more, and more confident, physics in their buildup. Finally, we show that there are large regions of acceptable parameter space in observational inputs that significantly overproduce the gamma-ray background, implying that the observed level of the background can yield significant constraints on models of the average cosmic gamma-ray emissivity associated with star formation.

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.

In this paper we develop a new semianalytical approach to quantifying the density profile of outer dark matter halos, motivated by the remarkable universality those profiles, for a wide range of dynamical parameters of the inner halos. We show that our minimalist model is robust under significant variations of its input parameters around the currently known values and we identify the turnaround radius as the most important scale of the problem. Based on that observation, we argue that the turnaround radius accurately represents the transition to the external density profile, and we provide a universal test based on geometrical characteristics of the outer profile as a proposal for measuring the turnaround radius of structures.

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).

The latest observations of extensive air showers (EAS) induced by ultra-high-energy cosmic rays (UHECR) appear to indicate, prima facie, a transition to heavy primaries at the highest energies. However, this interpretation, which is based on extrapolations of the Standard Model (SM) to ultra-LHC energies, is strained from both astrophysical and particle phenomenology perspectives. We consider the alternative that after some energy threshold, the first collision of the primary in the atmosphere results in a state, the decay of which leads to a considerably increased shower particle multiplicity, so that light-primary EAS appear heavy-like. We show that a minimal implementation of such a model yields predictions for the average EAS depth and shower-to-shower fluctuations that are consistent with each other, and an excellent fit to Auger data. If such an effect indeed takes place, we predict that: (a) the center-of-momentum (CM) energy threshold for the effect is of order 50 TeV; (b) the probability with which the effect occurs is high, and it will be detected easily by next-generation accelerators; (c) the increase in multiplicity compared to the SM prediction grows with CM energy roughly as $\sim E_{\rm CM}$; (d) the cosmic-ray composition at the highest energies is light. Remarkably, if the latter is confirmed electromagnetically this would necessitate the existence of new physics by these energies.

The Fermi-LAT revealed that the census of the gamma-ray sky is dominated by blazars. Looking for a possible connection between radio and gamma-ray emission is a central issue for understanding the blazar physics, and various works were dedicated to this topic. However, while a strong and significant correlation was found between radio and gamma-ray emission in the 0.1-100 GeV energy range, the connection between radio and very high energy (VHE, E>0.1 TeV) emission is still elusive. The main reason is the lack of a homogeneous VHE sky coverage, due to the operational mode of the imaging atmospheric Cherenkov telescopes. With the present work we aim to quantify and assess the significance of the possible connection between high-resolution radio emission, on milliarcsecond scale, and GeV-TeV gamma-ray emission in blazars. For achieving our goal we extract two large and unbiased blazar samples from the 1FHL and 2FHL Fermi catalogs, above 10 GeV and 50 GeV, respectively. To investigate how the correlation evolves as the gamma-ray energy increases, we perform the same analysis by using the 0.1-300 GeV 3FGL gamma-ray energy fluxes. When we consider the 0.1-300 GeV gamma-ray energy range, we find a strong and significant correlation for all of the blazar sub-classes. Conversely, when we consider the gamma-ray emission above 10 GeV the correlation with the radio emission vanishes, with the exception of the blazar sub-class of high synchrotron peaked objects.

The RoboPol instrument and the relevant program was developed in order to conduct a systematic study of the optical polarisation variability of blazars. Driven by the discovery that long smooth rotations of the optical polarisation plane can be associated with the activity in other bands and especially in gamma rays, the program was meant to investigate the physical mechanisms causing them and quantify the optical polarisation behaviour in blazars. Over the first three nominal observing seasons (2013, 2014 and 2015) RoboPol detected 40 rotations in 24 blazars by observing a gamma-ray-loud and gamma-ray-quite unbiassed sample of blazars, providing a reliable set of events for exploring the phenomenon. The obtain datasets provided the ground for a systematic quantification of the variability of the optical polarisation in such systems. In the following after a brief review of the discoveries that relate to the gamma-ray loudness of the sources we move on to discuss a simple jet model that explains the observed dichotomy in terms of polarisation between gamma-ray-loud and quite sources and the dependence of polarisation and the stability of the polarisation angle on the synchrotron peak frequency.

Black holes, anywhere in the stellar-mass to supermassive range, are often associated with relativistic jets. Models suggest that jet production may be a universal process common in all black hole systems regardless of their mass. Although in many cases observations support such hypotheses for microquasars and Seyfert galaxies, little is known on whether boosted blazar jets also comply with such universal scaling laws. We use uniquely rich multiwavelength radio light curves from the F-GAMMA program and the most accurate Doppler factors available to date to probe blazar jets in their emission rest frame with unprecedented accuracy. We identify for the first time a strong correlation between the blazar intrinsic broad-band radio luminosity and black hole mass, which extends over $\sim$ 9 orders of magnitude down to microquasars scales. Our results reveal the presence of a universal scaling law that bridges the observing and emission rest frames in beamed sources and allows us to effectively constrain jet models. They consequently provide an independent method for estimating the Doppler factor, and for predicting expected radio luminosities of boosted jets operating in systems of intermediate or tens-of-solar mass black holes, immediately applicable to cases as those recently observed by LIGO.

The Fermi Large Area Telescope (LAT) revealed that blazars, representing the most extreme radio-loud active galactic nuclei (AGN) population, dominate the census of the gamma-ray sky, and a significant correlation was found between radio and gamma-ray emission in the 0.1-100 GeV energy range. However, the possible connection between radio and very high energy (VHE, E>0.1 TeV) emission still remains elusive, owing to the lack of a homogeneous coverage of the VHE sky. The main goal of this work is to quantify and assess the significance of a possible connection between the radio emission on parsec scale measured by the very long baseline interferometry (VLBI) and GeV-TeV gamma-ray emission in blazars, which is a central issue for understanding the blazar physics and the emission processes. We investigate the radio VLBI and high energy gamma-ray emission by using two large and unbiased AGN samples extracted from the first and second Fermi-LAT catalogs of hard gamma-ray sources detected above 10 GeV (1FHL) and 50 GeV (2FHL). For comparison, we perform the same correlation analysis by using the 0.1-300 GeV gamma-ray energy flux provided by the third Fermi-LAT source catalog. We find that the correlation strength and significance depend on the gamma-ray energy range with a different behavior among the blazar sub-classes. Overall, the radio and gamma-ray emission above 10 GeV turns out to be uncorrelated for the full samples and for all of the blazar sub-classes with the exception of high synchrotron peaked (HSP) objects, which show a strong and significant correlation. On the contrary, when 0.1-300 GeV gamma-ray energies are considered, a strong and significant correlation is found for the full blazar sample as well as for all of the blazar sub-classes. We interpret and explain this correlation behavior within the framework of the blazar spectral energy distribution properties.

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.

Blazar population models have shown that the inverse-Compton and variability Doppler factor estimates yield consistent results at the population level for flat spectrum radio quasars (FSRQs). The two methods, however, are inconsistent when compared on a source-by-source basis. In this work, we attempt to understand the source of the discrepancy by tracing the potential sources of systematic and statistical error for the inverse-Compton Doppler factors. By eliminating these sources of error, we provide stronger constrains on the value of the Doppler factor in blazar jets. We re-estimate the inverse-Compton Doppler factor for 11 sources that meet certain criteria for their synchrotron peak frequency and the availability of Doppler factor estimates in the literature. We compare these estimates with the average of two different estimates of the variability Doppler factor obtained using various datasets and methodologies to identify any discrepancies and, in each case, trace their sources in the methodology or assumptions adopted. We identify three significant sources of error for the inverse-Compton Doppler factors: a) contamination of the X-ray flux by non-synchrotron self-Compton emission; b) radio observations at frequencies other than the synchrotron turnover frequency; c) non-simultaneity between radio and X-ray observations. We discuss key aspects in the correct application of the inverse-Compton method in light of these potential errors. We are able to constrain the Doppler factor of 3C273, 3C345, 3C454.3, PKS1510-089, and PKS1633+382 effectively, since all available estimates from both methods converge to the same values for these five sources.

Blazars are known to show periods of quiescence followed by outbursts visible throughout the electromagnetic spectrum. We present a novel maximum likelihood approach to capture this bimodal behavior by examining blazar radio variability in the flux-density domain. We separate quiescent and flaring components of a source's light curve by modeling its flux-density distribution as a series of "off" and "on" states. Our modeling allows us to extract information regarding the flaring ratio, duty cycle, and the modulation index in the "off"-state, in the "on"-state, as well as throughout the monitoring period of each blazar. We apply our method to a flux-density-limited subsample from the Owens Valley Radio observatory's 15 GHz blazar monitoring program, and explore differences in the variability characteristics between BL Lacs and FSRQs as well as between $\gamma$-ray detected and non-detected sources. We find that: (1) BL Lacs are more variable and have relatively larger outbursts than the FSRQs, (2) unclassified blazar candidates in our sample show similar variability characteristics as the FSRQs, and (3) $\gamma$-ray detected differ from the $\gamma$-ray non-detected sources in all their variability properties, suggesting a link between the production of $\gamma$-rays and the mechanism responsible for the radio variability. Finally, we fit distributions for blazar flaring ratios, duty cycles, and on- and off- modulation indices that can be used in population studies of variability-dependent blazar properties.

The radio light-curve of J1415+1320 (PKS 1413+135) shows time-symmetric and recurring U-shaped features across the cm-wave and mm-wave bands. The symmetry of these features points to lensing by an intervening object as the cause. U-shaped events in radio light curves in the cm-wave band have previously been attributed to Extreme scattering events (ESE). ESEs are thought to be the result of lensing by compact plasma structures in the Galactic interstellar medium, but the precise nature of these plasma structures remains unknown. Since the strength of a plasma lens evolves with wavelength $\lambda$ as $\lambda^2$, the presence of correlated variations at over a wide wavelength range casts doubt on the canonical ESE interpretation for J1415+1320. In this paper, we critically examine the evidence for plasma lensing in J1415+1320. We compute limits on the lensing strength, and the associated free-free opacity of the putative plasma lenses. We compare the observed and model ESE light curves, and also derive a lower limit on the lens distance based on the effects of parallax due to the Earth's orbit around the Sun. We conclude that plasma lensing is not a viable interpretation for J1415+1320's light curves and that symmetric U-shaped features in the radio light curves of extragalactic sources do not present {\em prima facie} evidence for ESEs. The methodology presented here is generic enough to be applicable to any plasma lensing candidate.

We report the discovery of a rare new form of long-term radio variability in the light-curves of active galaxies (AG) --- Symmetric Achromatic Variability (SAV) --- a pair of opposed and strongly skewed peaks in the radio flux density observed over a broad frequency range. We propose that SAV arises through gravitational milli-lensing when relativistically moving features in AG jets move through gravitational lensing caustics created by $10^3-10^6 \;{\rm M}_{\odot}$ subhalo condensates or black holes located within intervening galaxies. The lower end of this mass range has been inaccessible with previous gravitational lensing techniques. This new interpretation of some AG variability can easily be tested and if it passes these tests, will enable a new and powerful probe of cosmological matter distribution on these intermediate mass scales, as well as provide, for the first time, micro-arcsecond resolution of the nuclei of AG --- a factor of 30--100 greater resolution than is possible with ground-based millimeter VLBI.

In any flux-density limited sample of blazars, the distribution of the timescale modulation factor $\Delta t'/\Delta t$, which quantifies the change in observed timescales compared to the rest-frame ones due to redshift and relativistic compression follows an exponential distribution with a mean depending on the flux limit of the sample. In this work we produce the mathematical formalism that allows us to use this information in order to uncover the underlining rest-frame probability density function of measurable timescales of blazar jets. We extensively test our proposed methodology using a simulated FSRQ population with a 1.5 Jy flux-density limit in the simple case (where all blazars share the same intrinsic timescale), in order to identify limits of applicability and potential biases due to observational systematics and sample selection. We find that for monitoring with time intervals between observations longer than $\sim$30\% of the intrinsic timescale under investigation the method loses its ability to produce robust results. For time intervals of $\sim$3\% of the intrinsic timescale the error of the method is as low as 1\% in recovering the intrinsic rest-frame timescale. We applied our method to rotations of the optical polarization angle of blazars observed by RoboPol. We found that the intrinsic timescales of the longest-duration rotation event in each blazar follows a narrow distribution, well-described by a normal distribution with mean 87 days and standard deviation 5 days. We discuss possible interpretations of this result.

The RoboPol program has been monitoring the $R$-band linear polarisation parameters of an unbiased sample of 60 gamma-ray-loud blazars and a "control" sample of 15 gamma-ray-quite ones. The prime drive for the program has been the systematic study of the temporal behaviour of the optical polarisation and particularly the potential association of smooth and long rotations of the polarisation angle with flaring activity at high energies. Here we present the program and discuss a list of selected topics from our studies of the first three observing seasons (2013--2015) both in the angle and in the amplitude domain.

Detection of blazar pc scale counter-jets is difficult, but it can provide invaluable insight into the relativistic effects, radiative processes and the complex mechanisms of jet production, collimation and accelation in blazars. We build on recent populations models (optimized using the MOJAVE apparent velocity and redshift distributions) in order to derive the distribution of jet-to-counter-jet ratios and the flux densities of the counter-jet at different frequencies, in an effort to set minimum sensitivity limits required for existing and future telescope arrays in order to detect these elusive counter-jets. We find that: for the BL Lacs $5\%$ of their counter-jets have a flux-density higher than 100mJy, $15\%$ are higher than 10 mJy, and $32\%$ have higher flux-density than 1 mJy, whereas for the FSRQs $8\%$ have a flux-density higher than 10mJy, $17\%$ are higher than 1 mJy, and $32\%$ are higher than 0.1 mJy (at 15 GHz). Future telescopes like the SKA and newly operating like e-MERLIN and JVLA may detect up to $99\%$ of the BL Lac and $77\%$ of the FSRQ counter-jets. Sources with both low apparent velocity and a low Doppler factor make prime candidates for counter-jet detection. Combining our findings with literature values we have identified five such counter-jet detection candidates. Finally, we discuss possible effects beyond relativistic deboosting that may complicate the detection of counter-jets and that need to be accounted for in the interpretation of detections.

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 investigate the optical polarization properties of high-energy BL Lac objects using data from the RoboPol blazar monitoring program and the Nordic Optical Telescope. We wish to understand if there are differences in the BL Lac objects that are detected with the current-generation TeV instruments compared to those that have not yet been detected. 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 seems 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. 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. (Abridged)

To fully exploit the scientific potential of the Fermi mission, we initiated the F-GAMMA program. Between 2007 and 2015 it was the prime provider of complementary multi-frequency monitoring in the radio regime. We quantify the radio variability of gamma-ray blazars. We investigate its dependence on source class and examine whether the radio variability is related to the gamma-ray loudness. Finally, we assess the validity of a putative correlation between the two bands. The F-GAMMA monitored monthly a sample of about 60 sources at up to twelve radio frequencies between 2.64 and 228.39 GHz. We perform a time series analysis on the first 2.5-year dataset to obtain variability parameters. A maximum likelihood analysis is used to assess the significance of a correlation between radio and gamma-ray fluxes. We present light curves and spectra (coherent within ten days) obtained with the Effelsberg 100-m and IRAM 30-m telescopes. All sources are variable across all frequency bands with amplitudes increasing with frequency up to rest frame frequencies of around 60 - 80 GHz as expected by shock-in-jet models. Compared to FSRQs, BL Lacs show systematically lower variability amplitudes, brightness temperatures and Doppler factors at lower frequencies, while the difference vanishes towards higher ones. The time scales appear similar for the two classes. The distribution of spectral indices appears flatter or more inverted at higher frequencies for BL Lacs. Evolving synchrotron self-absorbed components can naturally account for the observed spectral variability. We find that the Fermi-detected sources show larger variability amplitudes as well as brightness temperatures and Doppler factors, than non-detected ones. Flux densities at 86.2 and 142.3 GHz correlate with 1 GeV fluxes at a significance level better than 3sigma, implying that gamma rays are produced very close to the mm-band emission region.

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.

We compare the variability properties of very high energy gamma-ray emitting BL Lac objects in the optical and radio bands. We use the variability information to distinguish multiple emission components in the jet, to be used as a guidance for spectral energy distribution modelling. Our sample includes 32 objects in the Northern sky that have data for at least 2 years in both bands. We use optical R-band data from the Tuorla blazar monitoring program and 15 GHz radio data from the Owens Valley Radio Observatory blazer monitoring program. We estimate the variability amplitudes using the intrinsic modulation index, and study the time-domain connection by cross-correlating the optical and radio light curves assuming power law power spectral density. Our sample objects are in general more variable in the optical than radio. We find correlated flares in about half of the objects, and correlated long-term trends in more than 40% of the objects. In these objects we estimate that at least 10%-50% of the optical emission originates in the same emission region as the radio, while the other half is due to faster variations not seen in the radio. This implies that simple single-zone spectral energy distribution models are not adequate for many of these objects.

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.

Low-frequency radio arrays are opening a new window for the study of the sky, both to study new phenomena and to better characterize known source classes. Being flat-spectrum sources, blazars are so far poorly studied at low radio frequencies. We characterize the spectral properties of the blazar population at low radio frequency compare the radio and high-energy properties of the gamma-ray blazar population, and search for radio counterparts of unidentified gamma-ray sources. We cross-correlated the 6,100 deg^2 Murchison Widefield Array Commissioning Survey catalogue with the Roma blazar catalogue, the third catalogue of active galactic nuclei detected by Fermi-LAT, and the unidentified members of the entire third catalogue of gamma-ray sources detected by \fermilat. When available, we also added high-frequency radio data from the Australia Telescope 20 GHz catalogue. We find low-frequency counterparts for 186 out of 517 (36%) blazars, 79 out of 174 (45%) gamma-ray blazars, and 8 out of 73 (11%) gamma-ray blazar candidates. The mean low-frequency (120--180 MHz) blazar spectral index is $\langle \alpha_\mathrm{low} \rangle=0.57\pm0.02$: blazar spectra are flatter than the rest of the population of low-frequency sources, but are steeper than at $\sim$GHz frequencies. Low-frequency radio flux density and gamma-ray energy flux display a mildly significant and broadly scattered correlation. Ten unidentified gamma-ray sources have a (probably fortuitous) positional match with low radio frequency sources. Low-frequency radio astronomy provides important information about sources with a flat radio spectrum and high energy. However, the relatively low sensitivity of the present surveys still misses a significant fraction of these objects. Upcoming deeper surveys, such as the GaLactic and Extragalactic All-Sky MWA (GLEAM) survey, will provide further insight into this population.

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 demonstrate that, in the context of the $\Lambda$CDM model, it is in principle possible to measure the value of the cosmological constant by tracing, across cosmic time, the evolution of the turnaround radius of cosmic structures. The novelty of the presented method is that it is local, in the sense that it uses the effect of the cosmological constant on the relatively short scales of cosmic structures and not on the dynamics of the Universe at its largest scales. In this way, it can provide an important consistency check for the standard cosmological model and can give signs of new physics, beyond $\Lambda$CDM.

In $\Lambda$CDM cosmology, structure formation is halted shortly after dark energy dominates the mass/energy budget of the Universe. A manifestation of this effect is that in such a cosmology the turnaround radius has an upper bound. Recently, a new, local, test for the existence of dark energy in the form of a cosmological constant was proposed based on this turnaround bound. Before designing an experiment that, through high-precision determination of masses and turnaround radii, will challenge $\Lambda$CDM cosmology, we have to answer two important questions: first, when turnaround-scale structures are predicted to be close enough to their maximum size, so that a possible violation of the bound may be observable. Second, which is the best mass scale to target for possible violations of the bound. Using the Press-Schechter formalism, we find that turnaround structures have in practice already stopped forming, and consequently, the turnaround radius of structures must be very close to the maximum value today. We also find that the mass scale of $\sim 10^{13} M_\odot$ characterizes turnaround structures that start to form in a statistically important number density today. This mass scale also separates turnaround structures with different cosmological evolution: smaller structures are no longer readjusting their mass distribution inside the turnaround scale, they asymptotically approach their ultimate abundance from higher values, and they are common enough to have, at some epoch, experienced major mergers with structures of comparable mass; larger structures exhibit the opposite behavior. We call this mass scale the transitional mass scale and we argue that it is the optimal for the purpose outlined above. As a corollary result, we explain the different accretion behavior of small and larger structures observed in already conducted numerical simulations.

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 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 ($\sim 30\%$), uncertainties are dominated by systematic effects. In the case of inverse Compton Doppler factors, random errors appear to dominate, but are significantly larger ($\sim 60\%$).

Observations of blazar jets are shrouded in relativistic effects, thus hindering our understanding of their intrinsic properties and dominant physical processes responsible for their generation, evolution,radiation and particle emission. In this work we focus on extracting information about timescales in the jet rest frame using a population-modeling approach. We employ Monte Carlo simulations to derive a simple population model for the intrinsic unbeamed luminosities and the Lorentz $\Gamma$ of blazar jets that adequately describe the observed redshift and apparent superluminal speed distributions for flux-limited blazar samples. We derive separate models for BL Lacs and Flat Spectrum Radio Quasars. We then use these models to compute the predicted distribution of Doppler factors in each blazar class, and address the following questions: (a) What is the relativistically induced spread in observed timescales (e.g., event duration, time lags) in a flux-limited sample of relativistic jets? (b) Could differences between BL Lacs and FSRQs observed in the time domain be attributed to differences in beaming between the two populations? (c) Is there a statistically preferred amount of beaming in a flux-limited sample? How large are statistical deviations from that preferred value? We use our findings to propose promising approaches in phenomenological studies of timescales in blazar jets.

Recent Planck results have shown that radiation from the cosmic microwave background passes through foregrounds in which aligned dust grains produce polarized dust emission, even in regions of the sky with the lowest level of dust emission. One of the most commonly used ways to remove the dust foreground is to extrapolate the polarized dust emission signal from frequencies where it dominates (e.g., ~ 350 GHz) to frequencies commonly targeted by cosmic microwave background experiments (e.g., ~150 GHz). In this paper, we describe an interstellar medium effect that can lead to decorrelation of the dust emission polarization pattern between different frequencies due to multiple contributions along the line of sight. Using a simple 2-cloud model we show that there are two conditions under which this decorrelation can be large: (a) the ratio of polarized intensities between the two clouds changes between the two frequencies; (b) the magnetic fields between the two clouds contributing along a line of sight are significantly misaligned. In such cases, the 350 GHz polarized sky map is not predictive of that at 150 GHz. We propose a possible correction for this effect, using information from optopolarimetric surveys of dichroicly absorbed starlight.

We revisit the relation between magnetic-field strength ($B$) and gas density ($\rho$) for contracting interstellar clouds and fragments (or, cores), which is central in observationally determining the dynamical importance of magnetic fields in cloud evolution and star formation. Recently, it has been claimed that a relation $B \propto \rho^{2/3} $ is statistically preferred over $B \propto \rho^{1/2}$ in molecular clouds, when magnetic field detections and nondetections from Zeeman observations are combined. This finding has unique observational implications on cloud and core geometry: The relation $B \propto \rho^{2/3} $ can only be realized under spherical contraction. However, no indication of spherical geometry can be found for the objects used in the original statistical analysis of the $B-\rho$ relation. We trace the origin of the inconsistency to simplifying assumptions in the statistical model used to arrive at the $B\propto \rho^{2/3}$ conclusion and to an underestimate of observational uncertainties in the determination of cloud and core densities. We show that, when these restrictive assumptions are relaxed, $B \propto \rho^{1/2}$ is the preferred relation for the (self-gravitating) molecular-cloud data, as theoretically predicted four decades ago.

Emission from blazar jets in the ultraviolet, optical, and infrared is polarized. If these low-energy photons were inverse-Compton scattered, the upscattered high-energy photons retain a fraction of the polarization. Current and future X-ray and gamma-ray polarimeters such as INTEGRAL-SPI, PoGOLITE, X-Calibur, Gamma-Ray Burst Polarimeter, GEMS-like missions, ASTRO-H, and POLARIX have the potential to discover polarized X-rays and gamma-rays from blazar jets for the first time. Detection of such polarization will open a qualitatively new window into high-energy blazar emission; actual measurements of polarization degree and angle will quantitatively test theories of jet emission mechanisms. We examine the detection prospects of blazars by these polarimetry missions using examples of 3C 279, PKS 1510-089, and 3C 454.3, bright sources with relatively high degrees of low-energy polarization. We conclude that while balloon polarimeters will be challenged to detect blazars within reasonable observational times (with X-Calibur offering the most promising prospects), space-based missions should detect the brightest blazars for polarization fractions down to a few percent. Typical flaring activity of blazars could boost the overall number of polarimetric detections by nearly a factor of five to six purely accounting for flux increase of the brightest of the comprehensive, all-sky, Fermi-LAT blazar distribution. The instantaneous increase in the number of detections is approximately a factor of two, assuming a duty cycle of 20% for every source. The detectability of particular blazars may be reduced if variations in the flux and polarization fraction are anticorrelated. Simultaneous use of variability and polarization trends could guide the selection of blazars for high-energy polarimetric observations.

Since mid-2007 we have carried out a dedicated long-term monitoring programme at 15 GHz using the Owens Valley Radio Observatory 40 meter telescope. One of the main goals of this programme is to study the relation between the radio and gamma-ray emission in blazars and to use it as a tool to locate the site of high energy emission. Using this large sample of objects we are able to characterize the radio variability, and study the significance of correlations between the radio and gamma-ray bands. We find that the radio variability of many sources can be described using a simple power law power spectral density, and that when taking into account the red-noise characteristics of the light curves, cases with significant correlation are rare. We note that while significant correlations are found in few individual objects, radio variations are most often delayed with respect to the gamma-ray variations. This suggests that the gamma-ray emission originates upstream of the radio emission. Because strong flares in most known gamma-ray-loud blazars are infrequent, longer light curves are required to settle the issue of the strength of radio-gamma cross-correlations and establish confidently possible delays between the two. For this reason continuous multiwavelength monitoring over a longer time period is essential for statistical tests of jet emission models.

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.

Our intuitive understanding of cosmic structure formation works best in scales small enough so that bound, relaxed gravitating systems are no longer adjusting their radius; and large enough so that space and matter follow the average expansion of the Universe. Yet one of the most robust predictions of $\Lambda$CDM cosmology concerns the scale that separates these limits: the turnaround radius, which is the non-expanding shell furthest away from the center of a bound structure. The maximum possible value of the turnaround radius within the framework of the $\Lambda$CDM model is, for a given mass $M$, equal to $(3GM/\Lambda c^2)^{1/3}$, with $G$ Newton's constant and $c$ the speed of light, independently of cosmic epoch, exact nature of dark matter, or baryonic effects. We discuss the possible use of this prediction as an observational test for $\Lambda$CDM cosmology.

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 modeled 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.

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 years of the mission of the Fermi Gamma-ray Space Telescope are cross-correlated with 4 years of 15GHz observations from the OVRO 40-m monitoring program. 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. Modeling 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-sigma (AO 0235+164), with only two more larger than even 2.25-sigma (PKS 1502+106 and B2 2308+34). Additionally, we find correlated variability in Mrk 421 when including a strong flare that occurred in July-September 2012. 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.

We present a study of the very high energy (VHE; E > 100 GeV) gamma-ray emission of the blazar PKS 1424+240 observed with the MAGIC telescopes. The primary aim of this paper is the multiwavelength spectral characterization and modeling of this blazar, which is made particularly interesting by the recent discovery of a lower limit of its redshift of z > 0.6 and makes it a promising candidate to be the most distant VHE source. The source has been observed with the MAGIC telescopes in VHE gamma rays for a total observation time of ~33.6 h from 2009 to 2011. The source was marginally detected in VHE gamma rays during 2009 and 2010, and later, the detection was confirmed during an optical outburst in 2011. The combined significance of the stacked sample is ~7.2 sigma. The differential spectra measured during the different campaigns can be described by steep power laws with the indices ranging from 3.5 +/- 1.2 to 5.0 +/- 1.7. The MAGIC spectra corrected for the absorption due to the extragalactic background light connect smoothly, within systematic errors, with the mean spectrum in 2009-2011 observed at lower energies by the Fermi-LAT. The absorption-corrected MAGIC spectrum is flat with no apparent turn down up to 400 GeV. The multiwavelength light curve shows increasing flux in radio and optical bands that could point to a common origin from the same region of the jet. The large separation between the two peaks of the constructed non-simultaneous spectral energy distribution also requires an extremely high Doppler factor if an one zone synchrotron self-Compton model is applied. We find that a two-component synchrotron self-Compton model describes the spectral energy distribution of the source well, if the source is located at z~0.6.

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.

For a general dark-energy equation of state, we estimate the maximum possible radius of massive structures that are not destabilized by the acceleration of the cosmological expansion. A comparison with known stable structures constrains the equation of state. The robustness of the constraint can be enhanced through the accumulation of additional astrophysical data and a better understanding of the dynamics of bound cosmic structures.

Radio emission in blazars -- the aligned subset of Active Galactic Nuclei (AGN) -- is produced by synchrotron electrons moving relativistically in their jet's magnetic field. Under the assumption of some degree of uniformity of the field, the emission can be highly polarized -- linearly and circularly. In the radio regime, the observed variability is in most of the cases attributed to flaring events undergoing opacity evolution, i.e. transitions from optically thick to thin emission (or vice versa). These transistions have a specific signature in the polarization parameter space (angle and magnitude) which can be traced with high cadence polarization monitoring and provide us with a unique probe of the microphysics of the emitting region. Here we present the full Stokes analysis of radio emission from blazars observed in the framework of the F-GAMMA program and discuss the case study of PKS\,1510$-$089 which has shown a prominent polarization event around MJD 55900.

The exact location of the gamma-ray emitting region in blazars is still controversial. In order to attack this problem we present first results of a cross-correlation analysis between radio (11 cm to 0.8 mm wavelength, F-GAMMA program) and gamma-ray (0.1-300 GeV) ~ 3.5 year light curves of 54 Fermi-bright blazars. We perform a source stacking analysis and estimate significances and chance correlations using mixed source correlations. Our results reveal: (i) the first highly significant multi-band radio and gamma-ray correlations (radio lagging gamma rays) when averaging over the whole sample, (ii) average time delays (source frame: 76+/-23 to 7+/-9 days), systematically decreasing from cm to mm/sub-mm bands with a frequency dependence tau_r,gamma (nu) ~ nu^-1, in good agreement with jet opacity dominated by synchrotron self-absorption, (iii) a bulk gamma-ray production region typically located within/upstream of the 3 mm core region (tau_3mm,gamma=12+/-8 days), (iv) mean distances between the region of gamma-ray peak emission and the radio "tau=1 photosphere" decreasing from 9.8+/-3.0 pc (11 cm) to 0.9+/-1.1 pc (2 mm) and 1.4+/-0.8 pc (0.8 mm), (v) 3 mm/gamma-ray correlations in 9 individual sources at a significance level where one is expected by chance (probability: 4 times 10^-6), (vi) opacity and "time lag core shift" estimates for quasar 3C 454.3 providing a lower limit for the distance of the bulk gamma-ray production region from the supermassive black hole (SMBH) of ~ 0.8-1.6 pc, i.e. at the outer edge of the Broad Line Region (BLR) or beyond. A 3 mm tau=1 surface at ~ 2-3 pc from the jet-base (i.e. well outside the "canonical BLR") finally suggests that BLR material extends to several pc distances from the SMBH.

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 use optical data from the Palomar Transient Factory (PTF) and the Catalina Real-Time Transient Survey (CRTS) to study the variability of gamma-ray detected and non-detected objects in a large population of active galactic nuclei (AGN) selected from the Candidate Gamma-Ray Blazar Survey and Fermi Gamma-Ray Space Telescope catalogs. Our samples include 714 sources with PTF data and 1244 sources with CRTS data. We calculate the intrinsic modulation index to quantify the optical variability amplitude in these samples. We find the gamma-ray detected objects to be more variable than the non-detected ones. The flat spectrum radio quasars (FSRQs) are more variable than the BL Lac objects in our sample, but the significance of the difference depends on the sample used. When dividing the objects based on their synchrotron peak frequency, we find the low synchrotron peaked (LSP) objects to be significantly more variable than the high synchrotron peaked (HSP) ones, explaining the difference between the FSRQs and BL Lacs. This could be due to the LSPs being observed near their electron energy peak, while in the HSPs the emission is caused by lower energy electrons, which cool more slowly. We also find a significant correlation between the optical and gamma-ray fluxes that is stronger in the HSP BL Lacs than in the FSRQs. The FSRQs in our sample are also more Compton dominated than the HSP BL Lacs. These findings are consistent with models where the gamma-ray emission of HSP objects is produced by the synchrotron self-Compton mechanism, while the LSP objects need an additional external Compton component that increases the scatter in the flux-flux correlation.

We present results from four years of twice-weekly 15 GHz radio monitoring of about 1500 blazars with the Owens Valley Radio Observatory 40 m telescope. Using the intrinsic modulation index to measure variability amplitude, we find that, with $>\!6\sigma$ significance, the radio variability of radio-selected gamma-ray-loud blazars is stronger than that of gamma-ray-quiet blazars. Our extended data set also includes at least 21 months of data for all AGN with `clean' associations in the Fermi Large Area Telescope First AGN catalogue, 1LAC. With these additional data we examine the radio variability properties of a gamma-ray-selected blazar sample. Within this sample, we find no evidence for a connection between radio variability amplitude and optical classification. In contrast, for our radio-selected sample we find that the BL Lac object subpopulation is more variable than the flat spectrum radio quasar (FSRQ) subpopulation. Radio variability is found to correlate with the synchrotron peak frequency, with low- and intermediate-synchrotron-peaked blazars varying less than high-synchrotron-peaked ones. We find evidence for a significant negative correlation between redshift and radio variability among bright FSRQs.

The intergalactic magnetic field (IGMF) may leave an imprint on the angular anisotropy of the extragalactic gamma-ray background through its effect on electromagnetic cascades triggered by interactions between very high energy photons and the extragalactic background light. A strong IGMF will deflect secondary particles produced in these cascades and will thus tend to isotropize lower energy cascade photons, thereby inducing a modulation in the anisotropy energy spectrum of the gamma-ray background. Here we present a simple, proof-of-concept calculation of the magnitude of this effect and demonstrate that current Fermi data already seem to prefer non-negligible IGMF values. The anisotropy energy spectrum of the Fermi gamma-ray background could thus be used as a probe of the IGMF strength.

The total anisotropy of a diffuse background composed of two or more sources, such as the Fermi-LAT--measured gamma-ray background, is set by the anisotropy of each source population and the contribution of each population to the total intensity. The total anisotropy as a function of energy (the anisotropy energy spectrum) will modulate as the relative contributions of the sources change, implying that the anisotropy energy spectrum also encodes the intensity spectrum of each source class. We develop techniques, applicable to any such diffuse background, for unraveling the intensity spectrum of each component source population given a measurement of the total intensity spectrum and the total anisotropy energy spectrum, without introducing \emph{a priori} assumptions about the spectra of the source classes. We demonstrate the potential of these methods by applying them to example scenarios for the composition of the Fermi-LAT gamma-ray background consistent with current data and feasible within 10 years of observation.

Blazars are powerful, variable emitters from radio to gamma-ray wavelengths. Even though the general picture of synchrotron emission at low energies and inverse Compton at the high energy end is well established, many important aspects of these remarkable objects are still not well understood. For example, even the location of the gamma-ray emission region is still not clearly established, with some theories locating it close to the black hole/accretion disk while others place it at parsec scales in the radio jet. Since mid-2007 we have carried out a large scale monitoring program at 15 GHz using the OVRO 40 m telescope. We are currently observing about 1700 blazars twice per week. The sample includes all the Fermi-LAT detected blazars north of declination -20 degrees. Here, we study the existence of correlated variability between these two bands for 86 sources bright enough to be detected weekly by LAT. The existence of correlated variability can be interpreted as an indication of a related spatial locations for the radio and gamma-ray emission, making the evaluation of its statistical significance a key goal of our program. A study of the statistical significance of these cross-correlations is presented along with a discussion of the Monte Carlo simulations used to evaluate them. More information about the conditions on the radio emission zone can be obtained through polarization monitoring which tells us about the configuration of the magnetic fields in this region. To study radio polarization variability we are building KuPol, a radio polarization receiver for the 12 to 18 GHz band that will replace the current total power receiver at the OVRO 40 meter telescope.

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and gamma-ray bands. Our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi-LAT, whereas ~40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-mm spectral slope of blazars is quite flat up to ~70GHz, above which it steepens to <\alpha>~-0.65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (\nupS) in the SED of FSRQs is the same in all the blazar samples with <\nupS>=10^13.1 Hz, while the mean inverse-Compton peak frequency, <\nupIC>, ranges from 10^21 to 10^22 Hz. The distributions of \nupS and of \nupIC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs and strongly depend on the selection method. The Compton dominance of blazars ranges from ~0.2 to ~100, with only FSRQs reaching values >3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at ~7 or more, and radio-selected blazars at values ~1, thus implying that the assumption that the blazar power is dominated by high-energy emission is a selection effect. Simple SSC models cannot explain the SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and \nupS predicted by the blazar sequence.

A large sample of known and likely gamma-ray blazars has been monitored twice per week since late 2007 at 15 GHz with the Owens Valley Radio Observatory (OVRO) 40-meter Telescope. The sample contains about 1700 sources, including the initial sample of 1158 sources above declination -20 degrees from the Candidate Gamma-Ray Blazar Survey (CGRaBS) plus all the blazars associated with Fermi-LAT detections as released in the Fermi AGN catalogs. Using statistical likelihood analyses, we compare the variability amplitude for various sub-populations within our sample. These include comparisons of gamma-ray-loud versus quiet objects, BL Lac objects versus flat-spectrum radio quasars, and a study of the variability amplitude trend with redshift. To learn about the location of the gamma-ray emission region we study the significance of peaks in the radio/gamma-ray cross-correlation using Monte Carlo simulations. First results for 52 sources with data from both the high-confidence Fermi Large Area Telescope Bright AGN Sample and the first 2 years of our monitoring program are presented. We find that assuming a power spectral density with power law slope of -2 at 15 GHz and -1.5 at gamma-ray energies, 7 of our objects show cross-correlations at the 3sigma level. We are now studying the physical significance of these correlations by further exploring the range of power law slopes that are consistent with the data. An extension of this to a larger sample and longer light curves is underway and preliminary results are presented. We also describe KuPol, the new digital Ku-band receiver being constructed for the 40-meter telescope. This new receiver will provide total intensity and linear polarization measurements over the 12-18 GHz band, with 16 MHz spectral resolution. The polarization data will provide important clues about the magnetic field configuration in the radio emission region.

Whether a correlation exists between the radio and gamma-ray flux densities of blazars is a long-standing question, and one that is difficult to answer confidently because of various observational biases which may either dilute or apparently enhance any intrinsic correlation between radio and gamma-ray luminosities. We introduce a novel method of data randomization to evaluate quantitatively the effect of these biases and to assess the intrinsic significance of an apparent correlation between radio and gamma-ray flux densities of blazars. The novelty of the method lies in a combination of data randomization in luminosity space (to ensure that the randomized data are intrinsically, and not just apparently, uncorrelated) and significance assessment in flux space (to explicitly avoid Malmquist bias and automatically account for the limited dynamical range in both frequencies). The method is applicable even to small samples that are not selected with strict statistical criteria. For larger samples we describe a variation of the method in which the sample is split in redshift bins, and the randomization is applied in each bin individually; this variation is designed to yield the equivalent to luminosity-function sampling of the underlying population in the limit of very large, statistically complete samples. We show that for a smaller number of redshift bins, the method yields a worse significance, and in this way it is conservative in that it does not assign a stronger, artificially enhanced significance. We demonstrate how our test performs as a function of number of sources, strength of correlation, and number of redshift bins used, and we show that while our test is robust against common-distance biases and associated false positives for uncorrelated data, it retains the power of other methods in rejecting the null hypothesis of no correlation for correlated data.

The F-GAMMA program is a coordinated effort to investigate the physics of Active Galactic Nuclei (AGNs) via multi-frequency monitoring of Fermi blazars. In the current study we show and discuss the evolution of broad-band radio spectra, which are measured at ten frequencies between 2.64 and 142 GHz using the Effelsberg 100-m and the IRAM 30-m telescopes. It is shown that any of the 78 sources studied can be classified in terms of their variability characteristics in merely 5 types of variability. It is argued that these can be attributed to only two classes of variability mechanisms. The first four types are dominated by spectral evolution and can be described by a simple two-component system composed of: (a) a steep quiescent spectral component from a large scale jet and (b) a time evolving flare component following the "Shock-in-Jet" evolutionary path. The fifth type is characterised by an achromatic change of the broad band spectrum, which could be attributed to a different mechanism, likely involving differential Doppler boosting caused by geometrical effects. Here we present the classification, the assumed physical scenario and the results of calculations that have been performed for the spectral evolution of flares.

The contribution of unresolved sources to the diffuse gamma-ray background could induce anisotropies in this emission on small angular scales. We analyze the angular power spectrum of the diffuse emission measured by the Fermi LAT at Galactic latitudes |b| > 30 deg in four energy bins spanning 1 to 50 GeV. At multipoles \ell \ge 155, corresponding to angular scales \lesssim 2 deg, angular power above the photon noise level is detected at >99.99% CL in the 1-2 GeV, 2-5 GeV, and 5-10 GeV energy bins, and at >99% CL at 10-50 GeV. Within each energy bin the measured angular power takes approximately the same value at all multipoles \ell \ge 155, suggesting that it originates from the contribution of one or more unclustered source populations. The amplitude of the angular power normalized to the mean intensity in each energy bin is consistent with a constant value at all energies, C_P/<I>^2 = 9.05 +/- 0.84 x 10^{-6} sr, while the energy dependence of C_P is consistent with the anisotropy arising from one or more source populations with power-law photon spectra with spectral index \Gamma_s = 2.40 +/- 0.07. We discuss the implications of the measured angular power for gamma-ray source populations that may provide a contribution to the diffuse gamma-ray background.

Spectral energy distributions (SEDs) and radio continuum spectra are presented for a northern sample of 104 extragalactic radio sources, based on the Planck Early Release Compact Source Catalogue (ERCSC) and simultaneous multifrequency data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous observations ranging from radio to gamma-rays. This is the first extensive frequency coverage in the radio and millimetre domains for an essentially complete sample of extragalactic radio sources, and it shows how the individual shocks, each in their own phase of development, shape the radio spectra as they move in the relativistic jet. The SEDs presented in this paper were fitted with second and third degree polynomials to estimate the frequencies of the synchrotron and inverse Compton (IC) peaks, and the spectral indices of low and high frequency radio data, including the Planck ERCSC data, were calculated. SED modelling methods are discussed, with an emphasis on proper, physical modelling of the synchrotron bump using multiplecomponents. Planck ERCSC data also suggest that the original accelerated electron energy spectrum could be much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shock. Furthermore in many cases the Planck data indicate that gamma-ray emission must originate in the same shocks that produce the radio emission.

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above $3\times 10^{18}$ eV, for all zenith angles between 0$^\circ$ and 60$^\circ$, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance.

Since late 2007, we have been monitoring a large sample of known and likely gamma-ray-loud blazars at 15 GHz twice per week with the Owens Valley Radio Observatory (OVRO) 40 m Telescope. Our initial sample included the 1158 sources above declination -20 degrees from the Candidate Gamma-Ray Blazar Survey (CGRaBS), and we have since added nearly 400 more sources, including all blazars associated with Fermi Large Area Telescope (LAT) detections in the First AGN Catalog (1LAC). Here, we describe the new sample and present results for 2008 through early 2011. Using statistical likelihood analyses, we compare the variability amplitude for various sub-populations within our sample. These include comparisons of gamma-ray-loud versus gamma-ray-quiet objects, BL Lac objects versus flat-spectrum radio quasars, and a study of the variability amplitude trend with redshift. We also describe KuPol, the new digital Ku-band receiver being constructed for the 40 m telescope. This new receiver will provide total intensity and linear polarization measurements over the 12-18 GHz band, with 16 MHz spectral resolution.

We report about the multiwavelength campaign on the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z = 0.5846) performed in 2010 July-September and triggered by high activity as measured by Fermi/LAT. The peak luminosity in the 0.1-100 GeV energy band exceeded, for the first time in this type of source, the value of 10^48 erg/s, a level comparable to the most powerful blazars. The comparison of the spectral energy distribution of the NLS1 PMN J0948+0022 with that of a typical blazar - like 3C 273 - shows that the power emitted at gamma rays is extreme.

We present a detailed statistical analysis of the correlation between radio and gamma-ray emission of the Active Galactic Nuclei (AGN) detected by Fermi during its first year of operation, with the largest datasets ever used for this purpose. We use both archival interferometric 8.4 GHz data (from the VLA and ATCA, for the full sample of 599 sources) and concurrent single-dish 15 GHz measurements from the Owens Valley Radio Observatory (OVRO, for a sub sample of 199 objects). Our unprecedentedly large sample permits us to assess with high accuracy the statistical significance of the correlation, using a surrogate-data method designed to simultaneously account for common-distance bias and the effect of a limited dynamical range in the observed quantities. We find that the statistical significance of a positive correlation between the cm radio and the broad band (E>100 MeV) gamma-ray energy flux is very high for the whole AGN sample, with a probability <1e-7 for the correlation appearing by chance. Using the OVRO data, we find that concurrent data improve the significance of the correlation from 1.6e-6 to 9.0e-8. Our large sample size allows us to study the dependence of correlation strength and significance on specific source types and gamma-ray energy band. We find that the correlation is very significant (chance probability <1e-7) for both FSRQs and BL Lacs separately; a dependence of the correlation strength on the considered gamma-ray energy band is also present, but additional data will be necessary to constrain its significance.

The spectral shapes of the contributions of different classes of unresolved gamma-ray emitters can provide insight into their relative contributions to the extragalactic gamma-ray background (EGB) and the natures of their spectra at GeV energies. We calculate the spectral shapes of the contributions to the EGB arising from BL Lacertae type objects (BL Lacs) and flat-spectrum radio quasars (FSRQs) assuming blazar spectra can be described as broken power laws. We fit the resulting total blazar spectral shape to the Fermi Large Area Telescope measurements of the EGB, finding that the best-fit shape reproduces well the shape of the Fermi EGB for various break scenarios. We conclude that a scenario in which the contribution of blazars is dominant cannot be excluded on spectral grounds alone, even if spectral breaks are shown to be common among Fermi blazars. We also find that while the observation of a featureless (within uncertainties) power-law EGB spectrum by Fermi does not necessarily imply a single class of contributing unresolved sources with featureless individual spectra, such an observation and the collective spectra of the separate contributing populations determine the ratios of their contributions. As such, a comparison with studies including blazar gamma-ray luminosity functions could have profound implications for the blazar contribution to the EGB, blazar evolution, and blazar gamma-ray spectra and emission.

Pulsars emerge in the Fermi era as a sizable population of gamma-ray sources. Millisecond pulsars (MSPs) constitute an older subpopulation whose sky distribution extends to high Galactic latitudes, and it has been suggested that unresolved members of this class may contribute a significant fraction of the measured large-scale isotropic gamma-ray background (IGRB). We investigate the possible energy-dependent contribution of unresolved MSPs to the anisotropy of the Fermi-measured IGRB. For observationally-motivated MSP population models, we show that the preliminary Fermi anisotropy measurement places an interesting constraint on the abundance of MSPs in the Galaxy and the typical MSP flux, about an order of magnitude stronger than constraints on this population derived from the intensity of the IGRB alone. We also examine the possibility of a MSP component in the IGRB mimicking a dark matter signal in anisotropy-based searches, and conclude that the energy dependence of an anisotropy signature would distinguish MSPs from all but very light dark matter candidates.

We report on a multiwavelength campaign on the radio-loud Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z=0.5846) performed in 2010 July-September and triggered by a high-energy gamma-ray outburst observed by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope. The peak flux in the 0.1-100 GeV energy band exceeded, for the first time in this type of source, the value of 10^-6 ph cm^-2 s^-1, corresponding to an observed luminosity of 10^48 erg s^-1. Although the source was too close to the Sun position to organize a densely sampled follow-up, it was possible to gather some multiwavelength data that confirmed the state of high activity across the sampled electromagnetic spectrum. The comparison of the spectral energy distribution of the NLS1 PMN J0948+0022 with that of a typical blazar - like 3C 273 - shows that the power emitted at gamma rays is extreme.

The origin of the extragalactic gamma-ray background is a pressing cosmological mystery. The Fermi Gamma-Ray Space Telescope has recently measured the intensity and spectrum of this background; both are substantially different from previous measurements. We present a novel calculation of the gamma-ray background from normal star-forming galaxies. Contrary to long-standing expectations, we find that numerous but individually faint normal galaxies may comprise the bulk of the Fermi signal, rather than rare but intrinsically bright active galaxies. This result has wide-ranging implications, including: the possibility to probe the cosmic star formation history with gamma rays; the ability to infer the cosmological evolution of cosmic rays and galactic magnetic fields; and an increased likelihood of identifying subdominant components from rare sources (e.g., dark matter clumps) through their large anisotropy.

The Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope provides an unprecedented opportunity to study gamma-ray blazars. To capitalize on this opportunity, beginning in late 2007, about a year before the start of LAT science operations, we began a large-scale, fast-cadence 15 GHz radio monitoring program with the 40-m telescope at the Owens Valley Radio Observatory (OVRO). This program began with the 1158 northern (declination>-20 deg) sources from the Candidate Gamma-ray Blazar Survey (CGRaBS) and now encompasses over 1500 sources, each observed twice per week with a ~4 mJy (minimum) and 3% (typical) uncertainty. Here, we describe this monitoring program and our methods, and present radio light curves from the first two years (2008 and 2009). As a first application, we combine these data with a novel measure of light curve variability amplitude, the intrinsic modulation index, through a likelihood analysis to examine the variability properties of subpopulations of our sample. We demonstrate that, with high significance (7-sigma), gamma-ray-loud blazars detected by the LAT during its first 11 months of operation vary with about a factor of two greater amplitude than do the gamma-ray-quiet blazars in our sample. We also find a significant (3-sigma) difference between variability amplitude in BL Lacertae objects and flat-spectrum radio quasars (FSRQs), with the former exhibiting larger variability amplitudes. Finally, low-redshift (z<1) FSRQs are found to vary more strongly than high-redshift FSRQs, with 3-sigma significance. These findings represent an important step toward understanding why some blazars emit gamma-rays while others, with apparently similar properties, remain silent.

The energy-dependence of the anisotropy (the anisotropy energy spectrum) of the large-scale diffuse gamma-ray background can reveal the presence of multiple source populations. Annihilating dark matter in the substructure of the Milky Way halo could give rise to a modulation in the anisotropy energy spectrum of the diffuse gamma-ray emission measured by Fermi, enabling the detection of a dark matter signal. We determine the detectability of a dark-matter-induced modulation for scenarios in which unresolved blazars are the primary contributor to the measured emission above ~1 GeV and find that in some scenarios pair-annihilation cross sections of order the value expected for thermal relic dark matter can produce a detectable feature. We anticipate that the sensitivity of this technique to specific dark matter models could be improved by tailored likelihood analysis methods.

Data collected by the Pierre Auger Observatory through 31 August 2007 showed evidence for anisotropy in the arrival directions of cosmic rays above the Greisen-Zatsepin-Kuz'min energy threshold, \nobreak{$6\times 10^{19}$eV}. The anisotropy was measured by the fraction of arrival directions that are less than $3.1^\circ$ from the position of an active galactic nucleus within 75 Mpc (using the V\'eron-Cetty and V\'eron $12^{\rm th}$ catalog). An updated measurement of this fraction is reported here using the arrival directions of cosmic rays recorded above the same energy threshold through 31 December 2009. The number of arrival directions has increased from 27 to 69, allowing a more precise measurement. The correlating fraction is $(38^{+7}_{-6})%$, compared with $21%$ expected for isotropic cosmic rays. This is down from the early estimate of $(69^{+11}_{-13})%$. The enlarged set of arrival directions is examined also in relation to other populations of nearby extragalactic objects: galaxies in the 2 Microns All Sky Survey and active galactic nuclei detected in hard X-rays by the Swift Burst Alert Telescope. A celestial region around the position of the radiogalaxy Cen A has the largest excess of arrival directions relative to isotropic expectations. The 2-point autocorrelation function is shown for the enlarged set of arrival directions and compared to the isotropic expectation.

The F -GAMMA program is a coordinated effort of several observing facilities to understand the AGN/blazar phenomenon via a multi-frequency monitoring approach, especially in the era of Fermi-GST. Some 60 prominent sources are monitored monthly with the Effelsberg 100-m telescope, the IRAM 30-m telescope and more frequently but in a less uniform fashion at the APEX 12-m telescope, covering from 2.64 to 345 GHz. The program has been running since January 2007 and here some of its findings are summarized. (a) There are two major variability patterns that the spectra of sources follow, one spectral-evolution-dominated and one achromatic. (b) The FSRQs show higher brightness temperatures indicative of larger Doppler factors at play and (c) a statistically significant radio-gamma-ray correlation has been found with a method recently suggested by Pavlidou et al. (in prep.).

Radio and gamma-ray emission from active galactic nuclei (AGN) are thought to share a common origin, related to the ejection phenomena in the vicinity of supermassive black holes. Thanks to its sensitivity, surveying capability, and broad energy range, the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope has permitted us to discover and characterize a huge number of extragalactic gamma-ray sources. Similarly to what was found by EGRET, these sources are typically associated with blazars, characterized by significant radio emission and flat spectrum. The radio luminosity distribution is extended over 7 orders of magnitudes, with flat spectrum radio quasars clustered at higher powers and BL Lacs more scattered; the average spectral index is consistent with alpha=0, although a few remarkable sources have alpha>0.5. A comparison of the radio flux density and the gamma-ray photon flux is presented, although claims on its significance require a detailed discussion and Monte Carlo simulations which will be presented in a future paper.

We present results from the first two years of our fast-cadence 15 GHz gamma-ray blazar monitoring program, part of the F-GAMMA radio monitoring project. Our sample includes the 1158 blazars north of -20 degrees declination from the Candidate Gamma-Ray Blazar Survey (CGRaBS), which encompasses a significant fraction of the extragalactic sources detected by the Fermi Gamma-ray Space Telescope. We introduce a novel likelihood analysis for computing a time series variability amplitude statistic that separates intrinsic variability from measurement noise and produces a quantitative error estimate. We use this method to characterize our radio light curves. We also present results indicating a statistically significant correlation between simultaneous average 15 GHz radio flux density and gamma-ray photon flux.

The F-GAMMA-project is the coordinated effort of several observatories to understand the AGN phenomenon and specifically blazars via multi-frequency monitoring in collaboration with the {\sl Fermi}-GST satellite since January 2007. The core observatories are: the Effelsberg 100-m, the IRAM 30-m and the OVRO 40-m telescope covering the range between 2.6 and 270 GHz. Effelsberg and IRAM stations do a monthly monitoring of the cm to sub-mm radio spectra of 60 selected blazars whereas the OVRO telescope is observing roughly 1200 objects at 15 GHz with a dense sampling of 2 points per week. The calibration uncertainty even at high frequencies, is of a few percent. 47% of the Effelsberg/Pico Veleta sample is included in the LBAS list. An update of the monitored sample is currently underway.

High-energy photons from cosmological emitters suffer attenuation due to pair production interactions with the extragalactic background light (EBL). The collective emission of any high-energy emitting cosmological population will exhibit an absorption feature at the highest energies. We calculate this absorption feature in the collective emission of blazars for various models of the blazar gamma-ray luminosity function (GLF) and the EBL. We find that models of the blazar GLF that predict higher relative contributions of high-redshift blazars to the blazar collective spectrum result in emission that is more susceptible to attenuation by the EBL, and hence result in more prominent absorption features, allowing for better differentiation amongst EBL models. We thus demonstrate that observations of such an absorption feature will contain information regarding both the blazar GLF and the EBL, and we discuss tests for EBL models and the blazar GLF that will become possible with upcoming Fermi observations.

Star-forming galaxies are predicted to contribute considerably to the cosmic gamma-ray background (CGB) as they are confirmed gamma-ray emitters and are the most numerous population of gamma-ray sources, although individually faint. Even though the Fermi Gamma-ray Space Telescope will be able to resolve few star-forming galaxies individually, their fractional contribution to the CGB should become far more significant than it was for past measurements of the CGB as many of the brighter, formerly unresolved sources will be resolved out. Thus, the clustering feature of galaxies imprinted on the CGB might be detectable by Fermi. In anticipation of such measurements, we calculate the predicted angular auto-power and cross-power spectra of the CGB from normal galaxies. We find that the amplitude of the auto-power spectrum is smaller than that for other sources such as blazars and dark-matter annihilation; the shape is also characteristic. We also show that the cross-power spectrum with galaxy surveys features larger amplitude. Fermi should be able to detect the correlation signature in both the auto-power and cross-power spectra at angular scales of ~1-10 deg after 5-yr of operation. Such a detection would be valuable in confirming the level of the star-forming galaxy contribution to the CGB, and more importantly, in serving as a tool in the effort to discriminate between possible origins of the CGB.

The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic rays. It combines a surface array to measure secondary particles at ground level together with a fluorescence detector to measure the development of air showers in the atmosphere above the array. The fluorescence detector comprises 24 large telescopes specialized for measuring the nitrogen fluorescence caused by charged particles of cosmic ray air showers. In this paper we describe the components of the fluorescence detector including its optical system, the design of the camera, the electronics, and the systems for relative and absolute calibration. We also discuss the operation and the monitoring of the detector. Finally, we evaluate the detector performance and precision of shower reconstructions.

Studies of the cosmic ray energy spectrum at the highest energies with the Pierre Auger Observatory.

Studies of the correlations of ultra-high energy cosmic ray directions with extra-Galactic objects, of general anisotropy, of photons and neutrinos, and of other astrophysical effects, with the Pierre Auger Observatory. Contributions to the 31st ICRC, Lodz, Poland, July 2009.

Technical reports on operations and features of the Pierre Auger Observatory, including ongoing and planned enhancements and the status of the future northern hemisphere portion of the Observatory. Contributions to the 31st International Cosmic Ray Conference, Lodz, Poland, July 2009.

Dark matter annihilation in Galactic substructure will produce diffuse gamma-ray emission of remarkably constant intensity across the sky, making it difficult to disentangle this Galactic dark matter signal from the extragalactic gamma-ray background. We show that if Galactic dark matter contributes a modest fraction of the measured emission in an energy range accessible to the Fermi Gamma-ray Space Telescope, the energy dependence of the angular power spectrum of the total measured emission could be used to confidently identify gamma rays from Galactic dark matter substructure.

Reports on the atmospheric monitoring, calibration, and other operating systems of the Pierre Auger Observatory. Contributions to the 31st International Cosmic Ray Conference, Lodz, Poland, July 2009.

Studies of the composition of the highest energy cosmic rays with the Pierre Auger Observatory, including examination of hadronic physics effects on the structure of extensive air showers. Submissions to the 31st ICRC, Lodz, Poland (July 2009).

We introduce a new method for testing departure from isotropy of points on a sphere based on an enhanced form of the two-point correlation function that we named 2pt+. This method uses information from the two extra variables that define the vector between two points on a sphere. We show that this is a powerful method to test departure from isotropy of a distribution of points on a sphere especially when the number of events is small. We apply the method to a few examples in astronomy and discuss the relevance for limited datasets, such as the case of ultra-high energy cosmic rays.

From direct observations of the longitudinal development of ultra-high energy air showers performed with the Pierre Auger Observatory, upper limits of 3.8%, 2.4%, 3.5% and 11.7% (at 95% c.l.) are obtained on the fraction of cosmic-ray photons above 2, 3, 5 and 10 EeV (1 EeV = 10^18 eV) respectively. These are the first experimental limits on ultra-high energy photons at energies below 10 EeV. The results complement previous constraints on top-down models from array data and they reduce systematic uncertainties in the interpretation of shower data in terms of primary flux, nuclear composition and proton-air cross-section.

The identification of point sources poses a great challenge for the high energy community. We present a new approach to evaluate the likelihood of a set of sources being a Galactic population based on the simple assumption that galaxies similar to the Milky Way host comparable populations of gamma-ray emitters. We propose a luminosity constraint on Galactic source populations which complements existing approaches by constraining the abundance and spatial distribution of any objects of Galactic origin, rather than focusing on the properties of a specific candidate emitter. We use M31 as a proxy for the Milky Way, and demonstrate this technique by applying it to the unidentified EGRET sources. We find that it is highly improbable that the majority of the unidentified EGRET sources are members of a Galactic halo population (e.g., dark matter subhalos), but that current observations do not provide any constraints on all of these sources being Galactic objects if they reside entirely in the disk and bulge. Applying this method to upcoming observations by the Fermi Gamma-ray Space Telescope has the potential to exclude association of an even larger number of unidentified sources with any Galactic source class.

The large majority of EGRET point sources remain without an identified low-energy counterpart, and a large fraction of these sources are most likely extragalactic. Whatever the nature of the extragalactic EGRET unidentified sources, faint unresolved objects of the same class must have a contribution to the diffuse extragalactic gamma-ray background (EGRB). Understanding this component of the EGRB, along with other guaranteed contributions from known sources, is essential if we are to use this emission to constrain exotic high-energy physics. Here, we follow an empirical approach to estimate whether a potential contribution of unidentified sources to the EGRB is likely to be important, and we find that it is. Additionally, we show how upcoming GLAST observations of EGRET unidentified sources, as well as of their fainter counterparts, can be combined with GLAST observations of the Galactic and extragalactic diffuse backgrounds to shed light on the nature of the EGRET unidentified sources even without any positional association of such sources with low-energy counterparts.

The spectral shape of the unresolved emission from different classes of gamma-ray emitters can be used to disentangle the contributions from these populations to the extragalactic gamma-ray background (EGRB). We present a calculation of the unabsorbed spectral shape of the unresolved blazar contribution to the EGRB starting from the spectral index distribution (SID) of resolved EGRET blazars derived through a maximum-likelihood analysis accounting for measurement errors. In addition, we explicitly calculate the uncertainty in this theoretically predicted spectral shape, which enters through the spectral index distribution parameters. We find that: (a) the unresolved blazar emission spectrum is only mildly convex, and thus, even if blazars are shown by GLAST to be a dominant contribution to the ERGB at lower energies, they may be insufficient to explain the EGRB at higher energies; (b) the theoretically predicted unresolved spectral shape involves significant uncertainties due to the limited constraints provided by EGRET data on the SID parameters, which are comparable to the statistical uncertainties of the observed EGRET EGRB at high energies; (c) the increased number statistics which will be provided by GLAST will be sufficient to reduce this uncertainty by at least a factor of three.

The intrinsic distribution of spectral indices in GeV energies of gamma-ray--loud blazars is a critical input in determining the spectral shape of the unresolved blazar contribution to the diffuse extragalactic gamma-ray background, as well as an important test of blazar emission theories. We present a maximum-likelihood method of determining the intrinsic spectral index distribution (ISID) of a population of gamma-ray emitters which accounts for error in measurement of individual spectral indices, and we apply it to EGRET blazars. We find that the most likely Gaussian ISID for EGRET blazars has a mean of 2.27 and a standard deviation of 0.20. We additionally find some indication that FSRQs and BL Lacs may have different ISIDs (with BL Lacs being harder). We also test for spectral index hardening associated with blazar variability for which we find no evidence. Finally, we produce simulated GLAST spectral index datasets and perform the same analyses. With improved statistics due to the much larger number of resolvable blazars, GLAST data will help us determine the ISIDs with much improved accuracy. Should any difference exist between the ISIDs of BL Lacs and FSRQs or between the ISIDs of blazars in the quiescent and flaring states, GLAST data will be adequate to separate these ISIDs at a significance better than 3 sigma.

We present an analytical description of the energetics of the population of cosmic accretion shocks, for a concordance cosmology. We calculate how the shock-processed accretion power and mass current are distributed among different shock Mach numbers, and how they evolve with cosmic time. We calculate the cumulative energy input of cosmic accretion shocks of any Mach number to the intergalactic medium as a function of redshift, and we compare it with the energy output of supernova explosions as well as with the energy input required to reionize the universe. In addition, we investigate and quantify the effect of environmental factors, such as local clustering properties and filament preheating on the statistical properties of these shocks. We find that the energy processed by accretion shocks is higher than the supernova energy output for z<3 and that it becomes more than an order of magnitude higher in the local universe. The energy processed by accretion shocks alone becomes comparable to the energy required to reionize the universe by z~3.5. Finally, we establish both qualitative and quantitatively that both local clustering as well as filament compression and preheating are important factors in determining the statistical properties of the cosmic accretion shock population.

We present a new analytic tool for the study of cosmic structure formation, a double distribution of the number density of dark matter halos with respect to both halo mass and local over- (or under-) density. The double distribution provides a statistical treatment of the properties of matter surrounding collapsed objects, and can be used to provide analytical insight into environmental effects on hierarchical structure formation. We apply this new tool to the case of cosmic accretion shocks. We investigate and quantify the effect of environmental factors on the statistical properties of these shocks. For this purpose, we explore two different models. The first control model uses a Press-Schechter mass function to describe the population of collapsed structures, and assumes that all objects accrete gas of the same density and temperature. The second model treats the accreted material as a multi-temperature, multi-density medium with densities and temperatures derived from the double distribution. We find that the shock environment significantly alters the physical impact of cosmic accretion shocks on the intergalactic medium, as well as the cosmic history of their properties.

The large majority of EGRET point sources remain to this day without an identified low-energy counterpart. Whatever the nature of the EGRET unidentified sources, faint unresolved objects of the same class must have a contribution to the diffuse gamma-ray background: if most unidentified objects are extragalactic, faint unresolved sources of the same class contribute to the background, as a distinct extragalactic population; on the other hand, if most unidentified sources are Galactic, their counterparts in external galaxies will contribute to the unresolved emission from these systems. Understanding this component of the gamma-ray background, along with other guaranteed contributions from known sources, is essential in any attempt to use gamma-ray observations to constrain exotic high-energy physics. Here, we follow an empirical approach to estimate whether a potential contribution of unidentified sources to the extragalactic gamma-ray background is likely to be important, and we find that it is. Additionally, we comment on how the anticipated GLAST measurement of the diffuse gamma-ray background will change, depending on the nature of the majority of these sources.

The third EGRET catalog contains a large number of unidentified sources. This subset of objects is expected to include known gamma-ray emitters of Galactic origin such as pulsars and supernova remnants, in addition to an extragalactic population of blazars. However, current data allows the intriguing possibility that some of these objects may represent a new class of yet undiscovered gamma-ray sources. Many theoretically motivated candidate emitters (e.g. clumps of annihilating dark matter particles) have been suggested to account for these detections. We take a new approach to determine to what extent this population is Galactic and to investigate the nature of the possible Galactic component. By assuming that galaxies similar to the Milky Way should host comparable populations of objects, we constrain the allowed Galactic abundance and distribution of various classes of gamma-ray sources using the EGRET data set. We find it is highly improbable that a large number of the unidentified sources are members of a Galactic halo population, but that a distribution of the sources entirely in the disk and bulge is plausible. Finally, we discuss the additional constraints and new insights that GLAST will provide.

We present a double distribution function of dark matter halos, with respect to both object mass and local over- (or under-) density. This analytical tool provides a statistical treatment of the properties of matter surrounding collapsed objects, and can be used to study environmental effects on hierarchical structure formation. The size of the "local environment" of a collapsed object is defined to depend on the mass of the object. The Press-Schechter mass function is recovered by integration of our double distribution over the density contrast. We also present a detailed treatment of the evolution of overdensities and underdensities in Einstein-deSitter and flat LCDM universes, according to the spherical evolution model. We explicitly distinguish between true and linearly extrapolated overdensities and provide conversion relations between the two quantities.

The diffuse extragalactic gamma-ray background (EGRB) above 100 MeV encodes unique information about high-energy processes in the universe. Numerous sources for the EGRB have been proposed, but the two systems which are certain to make some contribution are active galaxies (blazars) as well as normal galaxies. In this paper, we evaluate the contribution to the background from both sources. The active galaxy contribution arises from unresolved blazars. We compute this contribution using the Stecker-Salamon model. For normal galaxies, the emission is due to cosmic-ray interactions with diffuse gas. Our key assumption here is that the cosmic-ray flux in a galaxy is proportional to the supernova rate and thus the massive star formation rate, quantified observationally by the cosmic star formation rate (CSFR). In addition, the existence of stars today requires a considerably higher ISM mass in the past. Using the CSFR to compute both these effects, we find that normal galaxies are responsible for a significant portion (\sim 1/3) of the EGRB near 1 GeV, but make a smaller contribution at other energies. Finally, we present a "minimal" two-component model which includes contributions from both normal galaxies and blazars. We show that the spectrum of the diffuse radiation is a key constraint on this model: while neither the blazar spectra, nor the galactic spectra, are separately optimal fits to the observed spectrum, the combined emission provides an excellent fit. We close by noting key observational tests of this two-component model, which can be probed by future gamma-ray observatories such as GLAST.

Self-interacting dark matter (SIDM) was introduced by Spergel & Steinhardt to address possible discrepancies between collisionless dark matter simulations and observations on scales of less than 1 Mpc. We examine the case in which dark matter particles not only have strong self-interactions but also have strong interactions with baryons. The presence of such interactions will have direct implications for nuclear and particle astrophysics. Among these are a change in the predicted abundances from big bang nucleosynthesis (BBN) and the flux of gamma-rays produced by the decay of neutral pions which originate in collisions between dark matter and Galactic cosmic rays (CR). From these effects we constrain the strength of the baryon--dark matter interactions through the ratio of baryon - dark matter interaction cross section to dark matter mass, $s$. We find that BBN places a weak upper limit to this ratio $< 10^8 cm^2/g$. CR-SIDM interactions, however, limit the possible DM-baryon cross section to $< 5 \times 10^{-3} cm^2/g$; this rules out an energy-independent interaction, but not one which falls with center-of-mass velocity as $s \propto 1/v$ or steeper.

Diffuse gamma-ray radiation in galaxies is produced by cosmic ray interactions with the interstellar medium. With the completion of EGRET observations, the only extragalactic object from which there has been a positive detection of diffuse gamma-ray emission is the Large Magellanic Cloud. We systematically estimate the expected diffuse gamma-ray flux from Local Group galaxies, and determine their detectability by new generation gamma-ray observatories such as GLAST. For each galaxy, the expected gamma-ray flux depends only on its total gas content and its cosmic ray flux. We present a method for calculating cosmic ray flux in these galaxies in terms of the observed rate of supernova explosions, where cosmic ray acceleration is believed to take place. The difficulty in deriving accurate supernova rates from observational data is a dominant uncertainty in our calculations. We estimate the gamma-ray flux for Local Group galaxies and find that our predictions are consistent with the observations for the LMC and with the observational upper limits for the Small Magellanic Cloud and M31. Both the Andromeda galaxy, with a flux of $\sim 1.0 \times 10^{-8}$ photons sec$^{-1}$ cm$^{-2}$ above 100 MeV, and the SMC, with a flux of $\sim 1.7 \times 10^{-8}$ photons sec$^{-1}$ cm$^{-2}$ above 100 MeV, are expected to be observable by GLAST. M33 is at the limit of detectability with a flux of $\sim 0.11 \times 10^{-8}$ sec$^{-1}$ cm$^{-2}$. Other Local Group galaxies are at least two orders of magnitude below GLAST sensitivity.