Millisecond pulsars (MSPs) are a growing class of gamma-ray emitters. Pulsed gamma-ray signals have been detected from more than 40 MSPs with the Fermi Large Area Telescope (LAT). The wider radio beams and more compact magnetospheres of MSPs enable studies of emission geometries over a broader range of phase space than non-recycled radio-loud gamma-ray pulsars. We have modeled the gamma-ray light curves of 40 LAT-detected MSPs using geometric emission models assuming a vacuum retarded-dipole magnetic field. We modeled the radio profiles using a single-altitude hollow-cone beam, with a core component when indicated by polarimetry; however, for MSPs with gamma-ray and radio light curve peaks occurring at nearly the same rotational phase we assume that the radio emission is co-located with the gamma rays and caustic in nature. The best-fit parameters and confidence intervals are determined using a maximum likelihood technique. We divide the light curves into three model classes, with gamma-ray peaks trailing (Class I), aligned (Class II) or leading (Class III) the radio peaks. Outer gap and slot gap (two-pole caustic) models best fit roughly equal numbers of Class I and II, while Class III are exclusively fit with pair-starved polar cap models. Distinguishing between the model classes based on typical derived parameters is difficult. We explore the evolution of magnetic inclination angle with period and spin-down power, finding possible correlations. While the presence of significant off-peak emission can often be used as a discriminator between outer gap and slot gap models, a hybrid model may be needed.

We report on the discovery of gamma-ray pulsations from five millisecond pulsars (MSPs) using the Fermi Large Area Telescope (LAT) and timing ephemerides provided by various radio observatories. We also present confirmation of the gamma-ray pulsations from a sixth source, PSR J2051-0827. Five of these six MSPs are in binary systems: PSRs J1713+0747, J1741+1351, J1600-3053 and the two black widow binary pulsars PSRs J0610-2100 and 2051-0827. The only isolated MSP is the nearby PSR J1024-0719, which is also known to emit X-rays. We present X-ray observations in the direction of PSRs J1600-3053 and J2051-0827. While the latter is firmly detected, we an only give upper limits for the X-ray flux of the former. There are no dedicated X-ray observations available for the other 3 objects. The MSPs mentioned above, together with most of the MSPs detected by Fermi, are used to put together a sample of 30 gamma-ray MSPs which is used to study the morphology and phase connection of radio and gamma-ray pulse profiles. We show that MSPs with pulsed gamma-ray emission which is phase aligned with the radio emission present the steepest radio spectra and the largest magnetic fields at the light cylinder among all MSPs. As well, we also observe a trend towards very low, or undetectable, radio linear polarisation levels. These properties could be attributed to caustic radio emission produced at a range of different altitudes in the magnetosphere. We note that most of these characteristics are also observed in the Crab pulsar, the only other radio pulsar known to exhibit phase-aligned radio and gamma-ray emission.

We report the discovery of nine previously unknown gamma-ray pulsars in a blind search of data from the Fermi Large Area Telescope (LAT). The pulsars were found with a novel hierarchical search method originally developed for detecting continuous gravitational waves from rapidly rotating neutron stars. Designed to find isolated pulsars spinning at up to kHz frequencies, the new method is computationally efficient, and incorporates several advances, including a metric-based gridding of the search parameter space (frequency, frequency derivative and sky location) and the use of photon probability weights. The nine pulsars have spin frequencies between 3 and 12 Hz, and characteristic ages ranging from 17 kyr to 3 Myr. Two of them, PSRs J1803-2149 and J2111+4606, are young and energetic Galactic-plane pulsars (spin-down power above 6e35 erg/s and ages below 100 kyr). The seven remaining pulsars, PSRs J0106+4855, J0622+3749, J1620-4927, J1746-3239, J2028+3332, J2030+4415, J2139+4716, are older and less energetic; two of them are located at higher Galactic latitudes (|b| > 10 deg). PSR J0106+4855 has the largest characteristic age (3 Myr) and the smallest surface magnetic field (2e11 G) of all LAT blind-search pulsars. PSR J2139+4716 has the lowest spin-down power (3e33 erg/s) among all non-recycled gamma-ray pulsars ever found. Despite extensive multi-frequency observations, only PSR J0106+4855 has detectable pulsations in the radio band. The other eight pulsars belong to the increasing population of radio-quiet gamma-ray pulsars.

Using the Parkes radio telescope we have carried out deep observations of eleven unassociated gamma-ray sources. Periodicity searches of these data have discovered two millisecond pulsars, PSR J1103-5403 (1FGL J1103.9-5355) and PSR J2241-5236 (1FGL J2241.9-5236), and a long period pulsar, PSR J1604-44 (1FGL J1604.7-4443). In addition we searched for but did not detect any radio pulsations from six gammaray pulsars discovered by the Fermi satellite to a level of - 0.04 mJy (for pulsars with a 10% duty cycle). Timing of the millisecond pulsar PSR J1103-5403 has shown that its position is 9' from the centroid of the gamma-ray source. Since these observations were carried out, independent evidence has shown that 1FGL J1103.9-5355 is associated with the flat spectrum radio source PKS 1101-536. It appears certain that the pulsar is not associated with the gamma-ray source, despite the seemingly low probability of a chance detection of a radio millisecond pulsar. We consider that PSR J1604-44 is a chance discovery of a weak, long period pulsar and is unlikely to be associated with 1FGL J1604.7-4443. PSR J2241-5236 has a spin period of 2.2 ms and orbits a very low mass companion with a 3.5 hour orbital period. The relatively high flux density and low dispersion measure of PSR J2241-5236 makes it an excellent candidate for high precision timing experiments. The gamma-rays of 1FGL J2241.9-5236 have a spectrum that is well modelled by a power law with exponential cutoff, and phasebinning with the radio ephemeris results in a multi-peaked gamma-ray pulse profile. Observations with Chandra have identified a coincident X-ray source within 0.1" of the position of the pulsar obtained by radio timing

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.

We searched for radio pulsars in 25 of the non-variable, unassociated sources in the Fermi LAT Bright Source List with the Green Bank Telescope at 820 MHz. We report the discovery of three radio and gamma-ray millisecond pulsars (MSPs) from a high Galactic latitude subset of these sources. All of the pulsars are in binary systems, which would have made them virtually impossible to detect in blind gamma-ray pulsation searches. They seem to be relatively normal, nearby (<=2 kpc) millisecond pulsars. These observations, in combination with the Fermi detection of gamma-rays from other known radio MSPs, imply that most, if not all, radio MSPs are efficient gamma-ray producers. The gamma-ray spectra of the pulsars are power-law in nature with exponential cutoffs at a few GeV, as has been found with most other pulsars. The MSPs have all been detected as X-ray point sources. Their soft X-ray luminosities of ~10^{30-31} erg/s are typical of the rare radio MSPs seen in X-rays.

The intermediate-frequency peaked BL Lacertae (IBL) object 3C 66A is detected during 2007 - 2008 in VHE (very high energy: E > 100 GeV) gamma-rays with the VERITAS stereoscopic array of imaging atmospheric Cherenkov telescopes. An excess of 1791 events is detected, corresponding to a significance of 21.2 standard deviations (sigma), in these observations (32.8 hours live time). The observed integral flux above 200 GeV is 6% of the Crab Nebula's flux and shows evidence for variability on the time-scale of days. The measured energy spectrum is characterized by a soft power law with photon index Gamma = 4.1 +- 0.4_stat +- 0.6_sys. The radio galaxy 3C 66B is excluded as a possible source of the VHE emission.

(Abridged) We have conducted a detailed investigation of the broad-band spectral properties of the \gamma-ray selected blazars of the Fermi LAT Bright AGN Sample (LBAS). By combining our accurately estimated Fermi gamma-ray spectra with Swift, radio, infra-red, optical and other hard X-ray/gamma-ray data, collected within three months of the LBAS data taking period, we were able to assemble high-quality and quasi-simultaneous Spectral Energy Distributions (SED) for 48 LBAS blazars.The SED of these gamma-ray sources is similar to that of blazars discovered at other wavelengths, clearly showing, in the usual Log $\nu $ - Log $\nu$ F$_\nu$ representation, the typical broad-band spectral signatures normally attributed to a combination of low-energy synchrotron radiation followed by inverse Compton emission of one or more components. We have used these SEDs to characterize the peak intensity of both the low and the high-energy components. The results have been used to derive empirical relationships that estimate the position of the two peaks from the broad-band colors (i.e. the radio to optical and optical to X-ray spectral slopes) and from the gamma-ray spectral index. Our data show that the synchrotron peak frequency $\nu_p^S$ is positioned between 10$^{12.5}$ and 10$^{14.5}$ Hz in broad-lined FSRQs and between $10^{13}$ and $10^{17}$ Hz in featureless BL Lacertae objects.We find that the gamma-ray spectral slope is strongly correlated with the synchrotron peak energy and with the X-ray spectral index, as expected at first order in synchrotron - inverse Compton scenarios. However, simple homogeneous, one-zone, Synchrotron Self Compton (SSC) models cannot explain most of our SEDs, especially in the case of FSRQs and low energy peaked (LBL) BL Lacs. (...)

Although Galactic cosmic rays (protons and nuclei) are widely believed to be dominantly accelerated by the winds and supernovae of massive stars, definitive evidence of this origin remains elusive nearly a century after their discovery [1]. The active regions of starburst galaxies have exceptionally high rates of star formation, and their large size, more than 50 times the diameter of similar Galactic regions, uniquely enables reliable calorimetric measurements of their potentially high cosmic-ray density [2]. The cosmic rays produced in the formation, life, and death of their massive stars are expected to eventually produce diffuse gamma-ray emission via their interactions with interstellar gas and radiation. M 82, the prototype small starburst galaxy, is predicted to be the brightest starburst galaxy in gamma rays [3, 4]. Here we report the detection of >700 GeV gamma rays from M 82. From these data we determine a cosmic-ray density of 250 eV cm-3 in the starburst core of M 82, or about 500 times the average Galactic density. This result strongly supports that cosmic-ray acceleration is tied to star formation activity, and that supernovae and massive-star winds are the dominant accelerators.

The recent detection by the Fermi gamma-ray space telescope of high-energy gamma-rays from the radio galaxy NGC 1275 makes the observation of the very high energy (VHE: E > 100 GeV) part of its broadband spectrum particularly interesting, especially for the understanding of active galactic nuclei (AGN) with misaligned multi-structured jets. The radio galaxy NGC 1275 was recently observed by VERITAS at energies above 100 GeV for about 8 hours. No VHE gamma-ray emission was detected by VERITAS from NGC 1275. A 99% confidence level upper limit of 2.1% of the Crab Nebula flux level is obtained at the decorrelation energy of approximately 340 GeV, corresponding to 19% of the power-law extrapolation of the Fermi Large Area Telescope (LAT) result.

We present evidence that the very-high-energy (VHE, E > 100 GeV) gamma-ray emission coincident with the supernova remnant IC 443 is extended. IC 443 contains one of the best-studied sites of supernova remnant/molecular cloud interaction and the pulsar wind nebula CXOU J061705.3+222127, both of which are important targets for VHE observations. VERITAS observed IC 443 for 37.9 hours during 2007 and detected emission above 300 GeV with an excess of 247 events, resulting in a significance of 8.3 standard deviations (sigma) before trials and 7.5 sigma after trials in a point-source search. The emission is centered at 06 16 51 +22 30 11 (J2000) +- 0.03_stat +- 0.08_sys degrees, with an intrinsic extension of 0.16 +- 0.03_stat +- 0.04_sys degrees. The VHE spectrum is well fit by a power law (dN/dE = N_0 * (E/TeV)^-Gamma) with a photon index of 2.99 +- 0.38_stat +- 0.3_sys and an integral flux above 300 GeV of (4.63 +- 0.90_stat +- 0.93_sys) * 10^-12 cm^-2 s^-1. These results are discussed in the context of existing models for gamma-ray production in IC 443.

The VERITAS collaboration reports the detection of very-high-energy (VHE) gamma-ray emission from the high-frequency-peaked BL Lac object 1ES 1218+304 located at a redshift of z=0.182. A gamma-ray signal was detected with a statistical significance of 10.4 standard deviations (10.4 sigma) for the observations taken during the first three months of 2007, confirming the discovery of this object made by the MAGIC collaboration. The photon spectrum between ~160 GeV and ~1.8 TeV is well described by a power law with an index of Gamma = 3.08 +/- 0.34_stat +/- 0.2_sys. The integral flux is Phi(E > 200 GeV) = (12.2 +/- 2.6) X 10^-12 cm^-2 s^-1, which corresponds to ~6% of that of the Crab Nebula. The light curve does not show any evidence for VHE flux variability. Using lower limits on the density of the extragalactic background light in the near to mid-infrared we are able to limit the range of intrinsic energy spectra for 1ES 1218+304. We show that the intrinsic photon spectrum has an index that is harder than Gamma = 2.32 +/- 0.37_stat. When including constraints from the spectra of 1ES 1101-232 and 1ES 0229+200, the spectrum of 1ES 1218+304 is likely to be harder than Gamma = 1.86 +/- 0.37_stat.

Since September 2005, the Whipple 10m Gamma-ray Telescope has been operated primarily as a blazar monitor. The five Northern Hemisphere blazars that have already been detected at the Whipple Observatory, Markarian 421, H1426+428, Markarian 501, 1ES 1959+650 and 1ES 2344+514, are monitored routinely each night that they are visible. We report on the Markarian 421 observations taken from November 2005 to June 2006 in the gamma-ray, X-ray, optical and radio bands. During this time, Markarian 421 was found to be variable at all wavelengths probed. Both the variability and the correlations among different energy regimes are studied in detail here. A tentative correlation, with large spread, was measured between the X-ray and gamma-ray bands, while no clear correlation was evident among the other energy bands. In addition to this, the well-sampled spectral energy distribution of Markarian 421 (1101+384) is presented for three different activity levels. The observations of the other blazar targets will be reported separately.

The high-frequency-peaked BL-Lacertae object \objectname{1ES 0806+524}, at redshift z=0.138, was observed in the very-high-energy (VHE) gamma-ray regime by VERITAS between November 2006 and April 2008. These data encompass the two-, and three-telescope commissioning phases, as well as observations with the full four-telescope array. \objectname{1ES 0806+524} is detected with a statistical significance of 6.3 standard deviations from 245 excess events. Little or no measurable variability on monthly time scales is found. The photon spectrum for the period November 2007 to April 2008 can be characterized by a power law with photon index $3.6 \pm 1.0_{\mathrm{stat}} \pm 0.3_{\mathrm{sys}}$ between $\sim$300 GeV and $\sim$700 GeV. The integral flux above 300 GeV is $(2.2\pm0.5_{\mathrm{stat}}\pm0.4_{\mathrm{sys}})\times10^{-12}\:\mathrm{cm}^{2}\:\mathrm{s}^{-1}$ which corresponds to 1.8% of the Crab Nebula flux. Non contemporaneous multiwavelength observations are combined with the VHE data to produce a broadband spectral energy distribution that can be reasonably described using a synchrotron-self Compton model.

VERITAS, an Imaging Atmospheric Cherenkov Telescope (IACT) system for gammma-ray astronomy in the GeV-TeV range, has recently completed its first season of observations with a full array of four telescopes. A number of astrophysical gamma-ray sources have been detected, both galactic and extragalactic, including sources previously unknown at TeV energies. We describe the status of the array and some highlight results, and assess the technical performance, sensitivity and shower reconstruction capabilities.

We report the detection of very high-energy gamma-ray emission from the intermediate-frequency-peaked BL Lacertae object W Comae (z=0.102) by VERITAS. The source was observed between January and April 2008. A strong outburst of gamma-ray emission was measured in the middle of March, lasting for only four days. The energy spectrum measured during the two highest flare nights is fit by a power-law and is found to be very steep, with a differential photon spectral index of Gamma = 3.81 +- 0.35_stat +- 0.34_syst. The integral photon flux above 200GeV during those two nights corresponds to roughly 9% of the flux from the Crab Nebula. Quasi-simultaneous Swift observations at X-ray energies were triggered by the VERITAS observations. The spectral energy distribution of the flare data can be described by synchrotron-self-Compton (SSC) or external-Compton (EC) leptonic jet models, with the latter offering a more natural set of parameters to fit the data.

LS I +61 303 is one of only a few high-mass X-ray binaries currently detected at high significance in very high energy gamma-rays. The system was observed over several orbital cycles (between September 2006 and February 2007) with the VERITAS array of imaging air-Cherenkov telescopes. A signal of gamma-rays with energies above 300 GeV is found with a statistical significance of 8.4 standard deviations. The detected flux is measured to be strongly variable; the maximum flux is found during most orbital cycles at apastron. The energy spectrum for the period of maximum emission can be characterized by a power law with a photon index of Gamma=2.40+-0.16_stat+-0.2_sys and a flux above 300 GeV corresponding to 15-20% of the flux from the Crab Nebula.

The multiwavelength observation of the nearby radio galaxy M87 provides a unique opportunity to study in detail processes occurring in Active Galactic Nuclei from radio waves to TeV gamma-rays. Here we report the detection of gamma-ray emission above 250 GeV from M87 in spring 2007 with the VERITAS atmospheric Cherenkov telescope array and discuss its correlation with the X-ray emission. The gamma-ray emission is measured to be point-like with an intrinsic source radius less than 4.5 arcmin. The differential energy spectrum is fitted well by a power-law function: dPhi/dE=(7.4+-1.3_{stat}+-1.5_{sys})(E/TeV)^{-2.31+-0.17_{stat}+-0.2_{sys}} 10^{-9}m^{-2}s^{-1}TeV^{-1}. We show strong evidence for a year-scale correlation between the gamma-ray flux reported by TeV experiments and the X-ray emission measured by the ASM/RXTE observatory, and discuss the possible short-time-scale variability. These results imply that the gamma-ray emission from M87 is more likely associated with the core of the galaxy than with other bright X-ray features in the jet.

Observations of the Crab Nebula have proven to be the best tool to calibrate and to characterize the performance of a Cherenkov telescope. Scientifically, it is interesting to measure the energy spectrum of the Crab Nebula close to the inverse-Compton peak where a deviation is expected from the power law seen at energies above 300 GeV. Additionally, it is important to search for pulsed emission from the Crab Pulsar at energies beyond the 10 GeV upper limit of the EGRET pulsar detection. Since current models predict a cut-off in pulsed emission between 10 and 100 GeV, measurements at energies close to this range may help to discriminate between them. We observed the Crab extensively in the 2006-2007 season during the VERITAS 2- and 3-telescope commissioning phases. Using this data set we reconstructed a preliminary energy spectrum of the signal from the Crab Nebula. We also measured the optical pulsed signal to validate our GPS time-stamping and barycentering techniques and obtained an upper limit for the pulsed emission at gamma-ray energies.

Gamma-ray burst (GRB) observations at very high energies (VHE, E > 100 GeV) can impose tight constraints on some GRB emission models. Many GRB afterglow models predict a VHE component similar to that seen in blazars and plerions, in which the GRB spectral energy distribution has a double-peaked shape extending into the VHE regime. VHE emission coincident with delayed X-ray flare emission has also been predicted. GRB follow-up observations have had high priority in the observing program at the Whipple 10m Gamma-ray Telescope and GRBs will continue to be high priority targets as the next generation observatory, VERITAS, comes on-line. Upper limits on the VHE emission, at late times (>~4 hours), from seven GRBs observed with the Whipple Telescope are reported here.

We report on observations of the sky region around the unidentified TeV gamma-ray source TeV J2032+4130 carried out with the Whipple Observatory 10 m atmospheric Cherenkov telescope for a total of 65.5 hrs between 2003 and 2005. The standard two-dimensional analysis developed by the Whipple collaboration for a stand-alone telescope reveals an excess in the field of view at a pre-trials significance level of 6.1 standard deviations. The measured position of this excess is alpha(2000) =20 h 32 m 27 s, delta(2000) = 41 deg 39 min 17 s. The estimated integral flux for this gamma-ray source is about 8% of the Crab-Nebula flux. The data are consistent with a point-like source. Here we present a detailed description of the standard two-dimensional analysis technique used for the analysis of data taken with the Whipple Observatory 10 m telescope and the results for the TeV J2032+4130 campaign. We include a short discussion of the physical mechanisms that may be responsible for the observed gamma-ray emission, based on possible association with known astrophysical objects, in particular Cygnus OB2.

The first atmospheric Cherenkov telescope of VERITAS (the Very Energetic Radiation Imaging Telescope Array System) has been in operation since February 2005. We present here a technical description of the instrument and a summary of its performance. The calibration methods are described, along with the results of Monte Carlo simulations of the telescope and comparisons between real and simulated data. The analysis of TeV $\gamma$-ray observations of the Crab Nebula, including the reconstructed energy spectrum, is shown to give results consistent with earlier measurements. The telescope is operating as expected and has met or exceeded all design specifications.

In April-May 2003, the blazar 1ES 1959+650 showed an increased level of X-ray activity. This prompted a multiwavelength observation campaign with the Whipple 10 m gamma-ray telescope, the Rossi X-ray Timing Explorer, the Bordeaux Optical Observatory, and the University of Michigan Radio Astrophysical Observatory. We present the multiwavelength data taken from May 2, 2003 to June 7, 2003 and compare the source characteristics with those measured during observations taken during the years 2000 and 2002. The X-ray observations gave a data set with high signal-to-noise light curves and energy spectra; however, the gamma-ray observations did not reveal a major TeV gamma-ray flare. Furthermore, we find that the radio and optical fluxes do not show statistically significant deviations from those measured during the 2002 flaring periods. While the X-ray flux and X-ray photon index appear correlated during subsequent observations, the apparent correlation evolved significantly between the years 2000, 2002, and 2003. We discuss the implications of this finding for the mechanism that causes the flaring activity.

A calibrated laser pulse propagating through the atmosphere produces a flash of Rayleigh scattered light with an intensity that can be calculated very accurately when atmospheric conditions are good. This is used in a technique developed for the absolute calibration of ultra high energy cosmic ray fluorescence telescopes, and it can also be applied to imaging atmospheric Cherenkov telescopes (IACTs). In this paper we present the absolute calibration system being constructed and tested for the VERITAS project.