TeV Halos, extended regions of TeV gamma-ray emission around middle-aged pulsars, have recently been established as a new source class in gamma-ray astronomy. These halos have been attributed to relativistic electrons and positrons that have left the acceleration region close to the pulsar and are diffusing in the surrounding medium. Measuring the morphology of TeV Halos enables for the first time a direct measurement of the electron diffusion on scales of tens of parsecs. There are hints that the presence of relativistic particles affects the diffusion rate in the pulsars' surroundings. Understanding electron diffusion is necessary to constrain the origins of the apparent `excess' of cosmic-ray positrons at tens of GeV. TeV Halos can also be used to find mis-aligned pulsars, as well as study certain properties of the Galaxy's pulsar population. Future VHE gamma-ray instruments will detect more of those TeV Halos and determine how much pulsars contribute to the observed cosmic-ray electron and positron fluxes, and how they affect diffusion in their environments.

The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1~TeV-30~TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected fourteen new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected GeV gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region.

The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5 {\sigma} using the MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357 (z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies with the Large Area Telescope on board Fermi. Multi-wavelength observations suggest a subdivision of the high state into two distinct flux states. In the band covered by MAGIC, the variability time scale is estimated to be 6.4 +/- 1.9 days. Modeling the broadband spectral energy distribution with an external Compton model, the location of the emitting region is understood as originating in the jet outside the broad line region (BLR) during the period of high activity, while being partially within the BLR during the period of low (typical) activity. The observed VHE spectrum during the highest activity is used to probe the extragalactic background light at an unprecedented distance scale for ground-based gamma-ray astronomy.

The High Altitude Water Cherenkov (HAWC) observatory is a wide field-of-view detector sensitive to gamma rays of 100 GeV to a few hundred TeV. Located in central Mexico at 19 degrees North latitude and 4100 m above sea level, HAWC will observe gamma rays and cosmic rays with an array of water Cherenkov detectors. The full HAWC array is scheduled to be operational in Spring 2015. In this paper, we study the HAWC sensitivity to the gamma-ray signatures of high-mass (multi- TeV) dark matter annihilation. The HAWC observatory will be sensitive to diverse searches for dark matter annihilation, including annihilation from extended dark matter sources, the diffuse gamma-ray emission from dark matter annihilation, and gamma-ray emission from non-luminous dark matter subhalos. Here we consider the HAWC sensitivity to a subset of these sources, including dwarf galaxies, the M31 galaxy, the Virgo cluster, and the Galactic center. We simulate the HAWC response to gamma rays from these sources in several well-motivated dark matter annihilation channels. If no gamma-ray excess is observed, we show the limits HAWC can place on the dark matter cross-section from these sources. In particular, in the case of dark matter annihilation into gauge bosons, HAWC will be able to detect a narrow range of dark matter masses to cross-sections below thermal. HAWC should also be sensitive to non-thermal cross-sections for masses up to nearly 1000 TeV. The constraints placed by HAWC on the dark matter cross-section from known sources should be competitive with current limits in the mass range where HAWC has similar sensitivity. HAWC can additionally explore higher dark matter masses than are currently constrained.

We report on radio-to-gamma-ray observations during 2011 May-September of PMN J0948+0022, the first narrow-line Seyfert 1 (NLSy1) galaxy detected in gamma-rays by Fermi-LAT. Strong variability was observed in gamma-rays, with two flaring periods peaking on 2011 June 20 and July 28. The variability observed in optical and near-infrared seems to have no counterpart in gamma-rays. This different behaviour could be related to a bending and inhomogeneous jet or a turbulent extreme multi-cell scenario. The radio spectra showed a variability pattern typical of relativistic jets. The XMM spectrum shows that the emission from the jet dominates above 2 keV, while a soft X-ray excess is evident in the low-energy part of the X-ray spectrum. Models where the soft emission is partly produced by blurred reflection or Comptonisation of the thermal disc emission provide good fits to the data. The X-ray spectral slope is similar to that found in radio-quiet NLSy1, suggesting that a standard accretion disc is present, as expected from the high accretion rate. Except for the soft X-ray excess, unusual in jet-dominated AGNs, PMN J0948+0022 shows all characteristics of the blazar class.

We present more than three years of observations at different frequencies, from radio to high-energy gamma-rays, of the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z=0.585). This source is the first NLS1 detected at energies above 100 MeV and therefore can be considered the prototype of this emerging new class of gamma-ray emitting active galactic nuclei (AGN). The observations performed from 2008 August 1 to 2011 December 31 confirmed that PMN J0948+0022 generates a powerful relativistic jet, able to develop an isotropic luminosity at gamma-rays of the order of 10^48 erg s^-1, at the level of powerful quasars. The evolution of the radiation emission of this source in 2009 and 2010 followed the canonical expectations of relativistic jets, with correlated multiwavelength variability (gamma-rays followed by radio emission after a few months), but it was difficult to retrieve a similar pattern in the light curves of 2011. The comparison of gamma-ray spectra before and including 2011 data suggested that there was a softening of the high-energy spectral slope. We selected five specific epochs to be studied by modelling the broad-band spectrum, characterised by an outburst at gamma-rays or very low/high flux at other wavelengths. The observed variability can largely be explained either by changes in the injected power, the bulk Lorentz factor of the jet or the electron spectrum. The characteristic time scale of doubling/halving flux ranges from a few days to a few months, depending on the frequency and the sampling rate. The shortest doubling time scale at gamma-rays is 2.3+-0.5 days. These small values underline the need of highly-sampled multiwavelength campaigns to better understand the physics of these sources.

We report the probable identification of the X-ray counterpart to the gamma-ray pulsar PSR J2021+4026 using imaging with the Chandra X-ray Observatory ACIS and timing analysis with the Fermi satellite. Given the statistical and systematic errors, the positions determined by both satellites are coincident. The X-ray source position is R.A. 20h21m30.733s, Decl. +40 deg 26 min 46.04sec (J2000) with an estimated uncertainty of 1.3 arsec combined statistical and systematic error. Moreover, both the X-ray to gamma-ray and the X-ray to optical flux ratios are sensible assuming a neutron star origin for the X-ray flux. The X-ray source has no cataloged infrared-to-visible counterpart and, through new observations, we set upper limits to its optical emission of i' >23.0 mag and r' > 25.2mag. The source exhibits an X-ray spectrum with most likely both a powerlaw and a thermal component. We also report on the X-ray and visible light properties of the 43 other sources detected in our Chandra observation.

Single-zone synchrotron self-Compton and external Compton models are widely used to explain broad-band Spectral Energy Distributions (SEDs) of blazars from infrared to gamma-rays. These models bear obvious similarities to the homogeneous synchrotron cloud model which is often applied to explain radio emission from individual components of parsec-scale radio jets. The parsec-scale core, typically the brightest and most compact feature of blazar radio jet, could be the source of high-energy emission. We report on ongoing work to test this hypothesis by deriving the physical properties of parsec-scale radio emitting regions of twenty bright Fermi blazars using dedicated 5-43 GHz VLBA observations and comparing these parameters to results of SED modeling.

We established a sample of millimeter-wave and gamma-ray bright active galactic nuclei matching the WMAP catalog with the EGRET catalog, highest energy photons and the Fermi bright source list. We have monitored over 80 of these objects in the near infrared, obtaining over 2000 JHKs data points directly comparable with Fermi data. We present examples of correlated near infrared and gamma-ray activity of known blazars and recently identified sources.

Following the recent discovery of gamma rays from the radio-loud narrow-line Seyfert 1 galaxy PMN J0948+0022 (z=0.5846), we started a multiwavelength campaign from radio to gamma rays, which was carried out between the end of March and the beginning of July 2009. The source displayed activity at all the observed wavelengths: a general decreasing trend from optical to gamma-ray frequencies was followed by an increase of radio emission after less than two months from the peak of the gamma-ray emission. The largest flux change, about a factor of about 4, occurred in the X-ray band. The smallest was at ultraviolet and near-infrared frequencies, where the rate of the detected photons dropped by a factor 1.6-1.9. At optical wavelengths, where the sampling rate was the highest, it was possible to observe day-scale variability, with flux variations up to a factor of about 3. The behavior of PMN J0948+0022 observed in this campaign and the calculated power carried out by its jet in the form of protons, electrons, radiation and magnetic field are quite similar to that of blazars, specifically of flat-spectrum radio quasars. These results confirm the idea that radio-loud narrow-line Seyfert 1 galaxies host relativistic jets with power similar to that of average blazars.

This is a report on the findings of the extragalactic science working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper was commissioned by the American Physical Society, and the full white paper can be found on astro-ph (arXiv:0810.0444). This detailed section discusses extragalactic science topics including active galactic nuclei, cosmic ray acceleration in galaxies, galaxy clusters and large scale structure formation shocks, and the study of the extragalactic infrared and optical background radiation. The scientific potential of ground based gamma-ray observations of Gamma-Ray Bursts and dark matter annihilation radiation is covered in other sections of the white paper.

We report observations with the Chandra X-ray Observatory of a field in the gamma$-Cygni supernova remnant (SNR78.2+2.1) centered on the cataloged location of the unidentified, bright gamma-ray source 3EG J2020+4017. In this search for an X-ray counterpart to the gamma-ray source, we detected 30 X-ray sources. Of these, we found 17 strong-candidate counterparts in optical (visible through near-infrared) cataloged and an additional 3 through our optical observations. Based upon colors and (for several objects) optical spectra, nearly all the optically identified objects appear to be reddened main-sequence stars. None of the X-ray sources with an optical counterpart is a plausible X-ray counterpart to 3EG J2020+4017 --if that gamma-ray source is a spin-powered pulsar. Many of the 10 X-ray sources lacking optical counterparts are likely (extragalactic) active galactic nuclei, based upon the sky density of such sources. Although one of the 10 optically unidentified X-ray sources could be the gamma-ray source, there is no auxiliary evidence supporting such an identification.

We present the 90 micron luminosity function of the Final Analysis of the European Large Area ISO Survey (ELAIS), extending the sample size of our previous analysis (paper IV) by about a factor of 4. Our sample extends to z=1.1, around 50 times the comoving volume of paper IV, and 10^{7.7} < h^{-2}L/Lsun < 10^{12.5}. From our optical spectroscopy campaigns of the northern ELAIS 90 mircon survey (7.4 deg^2 in total, to S(90um)>70mJy), we obtained redshifts for 61% of the sample (151 redshifts) to B<21 identified at 7 microns, 15 microns, 20cm or with bright (B<18.5) optical identifications. The selection function is well-defined, permitting the construction of the 90 micron luminosity function of the Final Analysis catalogue in the ELAIS northern fields, which is in excellent agreement with our Preliminary Analysis luminosity function in the ELAIS S1 field from paper IV. The luminosity function is also in good agreement with the IRAS-based prediction of Serjeant & Harrison (2004), which if correct requires luminosity evolution of (1+z)^{3.4 +/- 1.0} for consistency with the source counts. This implies an evolution in comoving volume averaged star formation rate at z<~1 consistent with that derived from rest-frame optical and ultraviolet surveys.

Spectroscopy of the Crab nebula along different slit directions reveals the 3 dimensional structure of the optical nebula. On the basis of the linear radial expansion result first discovered by Trimble (1968), we make a 3D model of the optical emission. Results from a limited number of slit directions suggest that optical lines originate from a complicated array of wisps that are located in a rather thin shell, pierced by a jet. The jet is certainly not prominent in optical emission lines, but the direction of the piercing is consistent with the direction of the X-ray and radio jet. The shell's effective radius is ~ 79 seconds of arc, its thickness about a third of the radius and it is moving out with an average velocity 1160 km/s.

We present a stroboscopic system developed for optical observations of pulsars and its application in the CLYPOS survey. The stroboscopic device is connected to a GPS clock and provides absolute timing to the stroboscopic shutter relative to the pulsar's radio ephemerides. By changing the phase we can examine the pulsar's light curve. The precisely timed stroboscope in front of the CCD camera can perform highly accurate time resolved pulsar photometry and offers the advantages of CCD cameras, which are high quantum efficiency as well as relatively large field of view, which is important for flux calibrations. CLYPOS (Cananea Ljubljana Young Pulsar Optical Survey) is an extensive search for optical counterparts of about 30 northern hemisphere radio pulsars. It is a collaboration between the INAOE, Mexico and the Faculty of Mathematics and Physics of the University of Ljubljana. Stroboscopic observations were done between December 1998 and November 2000 at the 2.12 m telescope of the Observatory Guillermo Haro in Cananea, Sonora. The first results of the survey are presented. Analyzed data indicate that there is no optical counterpart brighter than ~22.

Optical seeing measurements carried out at Sierra Negra, the site of the Large Millimeter Telescope, are reported. The site, one of the highest peaks of Central Mexico, offers good coverage of Northern and Southern hemispheres and we have undertaken several campaigns to investigate the astronomical potential of the site in the optical. Here we report on our campaign aimed at establishing the seeing quality of the site. We present data of the first three campaigns of optical seeing monitoring covering from February 2000 to May 2002, carried out with a Differential Image Motion Monitor. The results clearly indicate a sub-arcsec seeing, better statistics during the dry season and no dependence with the time of night. We find no dependence of our results with the integration time used.

Optical photometry for the pulsar PSR0531+21 has been extended with new observations that strengthen evidence for a previously observed 60 seconds periodicity. This period is found to be increasing with time at approximately the same rate as the rotational period of the pulsar. The observed period and its time dependence fit a simple free precession model.

Photometric data of the Crab pulsar, obtained in stroboscopic mode over a period of more than eight years, are presented here. The applied Fourier analysis reveals a faint 60 second modulation of the pulsar's signal, which we interpret as a free precession of the pulsar.

The EGRET telescope has repeatedly observed 3EG J1835+5918 as a bright and steady source of high-energy gamma-ray radiation which has not yet been indentified. EGRET data from CGRO observation cycle 1 to 7 have been reanalysed above 100 MeV and above 1 GeV. The gamma-ray source location represents the latest and probably final positional assessment based on EGRET data. We especially adress the question of flux variability and spectral variability. The results of a X-ray/optical identification campaign towards 3EG J1835+5918 are given. The one object which might be associated with the gamma-ray source 3EG J1835+5918 has the characteristics of an isolated neutron star and possible of a radio-quiet pulsar.

The high redshift (z=0.997) blazar B2 1308+326 was observed contemporaneously at x-ray, optical and radio wavelengths in June 1996. The x-ray observations were performed with ASCA. The ASCA results were found to be consistent with reanalysed data from two earlier ROSAT observations. The combined ASCA and ROSAT data reveal an x-ray spectrum that is best fit by a broken power law with absorber model. The break in the x-ray spectrum is interpreted, to be the emerging importance of inverse Compton (IC) emission which dominates the ASCA spectrum. The faint optical state reported for these observations (m_V=18.3+/-0.25) is incompatible with the high synchrotron flux previously detected by ROSAT. The IC emission detected by both ROSAT and ASCA was not significantly affected by the large change in the synchrotron component. MgII emission was detected with an equivalent width (EW) significantly different from previously reported values. Absorption at a level of in excess of the Galactic value was detected, indicating the possible presence of a foreground absorber. A gravitational microlensing scenario cannot therefore be ruled out for this blazar. B2 1308+326 could be a typical radio-selected BL Lac in terms of peak synchrotron frequency and optical and radio variability but its high bolometric luminosity, variable line emission and high Doppler boost factor make it appear more like a quasar than a BL Lac. It is suggested that B2 1308+326 be considered as the prototype of this class of composite source.

A dedicated stroboscopic device was used to obtain optical spectra of the Crab main-pulse and inter-pulse as well as the spectrum of the underlying nebula when the pulsar is turned off. As the nebular emission is very inhomogeneous, our ability to effectively subtract the nebular background signal is crucial. No spectral lines intrinsic to the pulsar are detected. The main-pulse and the inter-pulse behave as power laws, both with the same de-reddened index Alpha = +0.2 +- 0.1. This value was obtained by subtracting the nebular spectrum at the exact position of the pulsar. The underlying nebula is redder, Alpha = -0.4 +- 0.1. Its emission lines are split into approaching (sim. -1200 km/s) and receding (sim. +600 km/s) components. The strength of emission line components and the flux in nebular continuum vary on arcsec scale. The nebular line and continuum intensities along the N-S slit are given.

We present the catalogue, mask, redshift data and selection function for the PSCz survey of 15411 IRAS galaxies across 84% of the sky. Most of the IRAS data is taken from the Point Source Catalog, but this has been supplemented and corrected in various ways to improve the completeness and uniformity. We quantify the known imperfections in the catalogue, and we assess the overall uniformity, completeness and data quality. We find that overall the catalogue is complete and uniform to within a few percent at high latitudes and 10% at low latitudes. Ancillary information, access details, guidelines and caveats for using the catalogue are given.

We report on the first firm detection of pulsed gamma-ray emission from PSR B1509-58 in the 0.75-30 MeV energy range in CGRO COMPTEL data collected over more than 6 years. The modulation significance in the 0.75-30 MeV pulse-phase distribution is 5.4 sigma and the lightcurve is similar to the lightcurves found earlier between 0.7 and 700 keV: a single broad asymmetric pulse reaching its maximum 0.38 +/- 0.03 in phase after the radio peak, compared to the offset of 0.30 found in the CGRO BATSE soft gamma-ray data, and 0.27 +/- 0.01 for RXTE (2-16 keV), compatible with ASCA (0.7-2.2 keV). Spectral analysis based on the excess counts in the broad pulse of the lightcurve shows that extrapolation of the OSSE power-law spectral fit with index -1.68 describes our data well up to 10 MeV. Above 10 MeV the spectrum breaks abruptly. The precise location of the break/bend between 10 and 30 MeV depends on the interpretation of the structure in the lightcurve measured by COMPTEL and EGRET above 10 MeV.

We report the COMPTEL detection of pulsed $\gamma$-emission from PSR B1509-58 up to at least 10 MeV using data collected over more than 6 years. The 0.75-10 MeV lightcurve is broad and reaches its maximum near radio-phase 0.38, slightly beyond the maximum found at hard X-rays/ soft $\gamma$-rays. In the 10-30 MeV energy range a strong source is present in the skymap positionally consistent with the pulsar, but we do not detect significant pulsed emission. However, the lightcurve is consistent with the pulse shape changing from a single broad pulse into a double-peak morphology. Our results significantly constrain pulsar modelling.

We present a search for high energy gamma-ray emission from 9 nearby starburst galaxies and M31 with the EGRET instrument aboard CGRO. Though the diffuse gamma-ray emission from starburst galaxies was suspected to be detectable, we find no emission from NGC 253, M82 nor from the average of all 9 galaxies. The 2 sigma upper limit for the EGRET flux above 100 MeV for the averaged survey observations is 1.8 x 10-8 ph cm-2 s-1. From a model of the expected radio and gamma-ray emission, we find that the magnetic field in the nuclei of these galaxies is > 25 micro Gauss, and the ratio of proton and electron densities is < 400. The EGRET limits indicate that the rate of massive star formation in the survey galaxies is only about an order of magnitude higher than in the Milky Way. The upper limit to the gamma-ray flux above 100 MeV for M31 is 1.6 x 10-8 ph cm-2 s-1. At the distance of M31, the Milky Way flux would be over twice this value, indicating higher gamma-ray emissivities in our Galaxy. Therefore, since the supernova rate of the Milky Way is higher than in M31, our null detection of M31 supports the theory of the supernova origin of cosmic rays in galaxies.

We present a detailed analysis of the high energy gamma-ray source 2EG J0008+7307. The source has a steady flux and a hard spectrum, softening above 2 GeV. The properties of the gamma-ray source are suggestive of emission from a young pulsar in the spatially coincident CTA 1 supernova remnant, which has recently been found to have a non-thermal X-ray plerion. Our 95% uncertainty contour around the >1 GeV source position includes the point-like X-ray source at the centre of the plerion. We propose that this object is a young pulsar and is the most likely counterpart of 2EG J0008+7307.