Millisecond pulsars, old neutron stars spun-up by accreting matter from a companion star, can reach high rotation rates of hundreds of revolutions per second. Until now, all such "recycled" rotation-powered pulsars have been detected by their spin-modulated radio emission. In a computing-intensive blind search of gamma-ray data from the Fermi Large Area Telescope (with partial constraints from optical data), we detected a 2.5-millisecond pulsar, PSR J1311-3430. This unambiguously explains a formerly unidentified gamma-ray source that had been a decade-long enigma, confirming previous conjectures. The pulsar is in a circular orbit with an orbital period of only 93 minutes, the shortest of any spin-powered pulsar binary ever found.

The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009.

Forty six gamma-ray pulsars were reported in the First Fermi Large Area Telescope (LAT) Catalog of Gamma-ray Pulsars. Over forty more have been seen since then. A simple but effective figure-of-merit for gamma-detectability is sqrt(Edot)/d^2, where Edot is the pulsar spindown power and d the distance. We are tracking down the best gamma-ray candidates not yet seen. We present the timing and spectral analysis results of some new high spindown power, nearby gamma-ray pulsars. We also update some population distribution plots in preparation for the 2nd Fermi LAT gamma-ray Pulsar Catalog.

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.

We report the discovery of two millisecond pulsars in a search for radio pulsations at the positions of \emph{Fermi Large Area Telescope} sources with no previously known counterparts, using the Nan\c{c}ay radio telescope. The two millisecond pulsars, PSRs J2017+0603 and J2302+4442, have rotational periods of 2.896 and 5.192 ms and are both in binary systems with low-eccentricity orbits and orbital periods of 2.2 and 125.9 days respectively, suggesting long recycling processes. Gamma-ray pulsations were subsequently detected for both objects, indicating that they power the associated \emph{Fermi} sources in which they were found. The gamma-ray light curves and spectral properties are similar to those of previously-detected gamma-ray millisecond pulsars. Detailed modeling of the observed radio and gamma-ray light curves shows that the gamma-ray emission seems to originate at high altitudes in their magnetospheres. Additionally, X-ray observations revealed the presence of an X-ray source at the position of PSR J2302+4442, consistent with thermal emission from a neutron star. These discoveries along with the numerous detections of radio-loud millisecond pulsars in gamma rays suggest that many \emph{Fermi} sources with no known counterpart could be unknown millisecond pulsars.

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.

The high-frequency-peaked BL Lacertae object RGB J0710+591 was observed in the very high-energy (VHE; E > 100 GeV) wave band by the VERITAS array of atmospheric Cherenkov telescopes. The observations, taken between 2008 December and 2009 March and totaling 22.1 hr, yield the discovery of VHE gamma rays from the source. RGB J0710+591 is detected at a statistical significance of 5.5 standard deviations (5.5{\sigma}) above the background, corresponding to an integral flux of (3.9 +/- 0.8) x 10-12 cm-2 s-1 (3% of the Crab Nebula's flux) above 300 GeV. The observed spectrum can be fit by a power law from 0.31 to 4.6 TeV with a photon spectral index of 2.69 +/- 0.26stat +/- 0.20sys. These data are complemented by contemporaneous multiwavelength data from the Fermi Large Area Telescope, the Swift X-ray Telescope, the Swift Ultra-Violet and Optical Telescope, and the Michigan-Dartmouth-MIT observatory. Modeling the broadband spectral energy distribution (SED) with an equilibrium synchrotron self-Compton model yields a good statistical fit to the data. The addition of an external-Compton component to the model does not improve the fit nor brings the system closer to equipartition. The combined Fermi and VERITAS data constrain the properties of the high-energy emission component of the source over 4 orders of magnitude and give measurements of the rising and falling sections of the SED.

The first published Fermi large area telescope (Fermi-LAT) measurement of the isotropic diffuse gamma-ray emission is in good agreement with a single power law, and is not showing any signature of a dominant contribution from dark matter sources in the energy range from 20 to 100 GeV. We use the absolute size and spectral shape of this measured flux to derive cross section limits on three types of generic dark matter candidates: annihilating into quarks, charged leptons and monochromatic photons. Predicted gamma-ray fluxes from annihilating dark matter are strongly affected by the underlying distribution of dark matter, and by using different available results of matter structure formation we assess these uncertainties. We also quantify how the dark matter constraints depend on the assumed conventional backgrounds and on the Universe's transparency to high-energy gamma-rays. In reasonable background and dark matter structure scenarios (but not in all scenarios we consider) it is possible to exclude models proposed to explain the excess of electrons and positrons measured by the Fermi-LAT and PAMELA experiments. Derived limits also start to probe cross sections expected from thermally produced relics (e.g. in minimal supersymmetry models) annihilating predominantly into quarks. For the monochromatic gamma-ray signature, the current measurement constrains only dark matter scenarios with very strong signals.

Dark matter (DM) particle annihilation or decay can produce monochromatic $\gamma$-rays readily distinguishable from astrophysical sources. $\gamma$-ray line limits from 30 GeV to 200 GeV obtained from 11 months of Fermi Large Area Space Telescope data from 20-300 GeV are presented using a selection based on requirements for a $\gamma$-ray line analysis, and integrated over most of the sky. We obtain $\gamma$-ray line flux upper limits in the range $0.6-4.5\times 10^{-9}\mathrm{cm}^{-2}\mathrm{s}^{-1}$, and give corresponding DM annihilation cross-section and decay lifetime limits. Theoretical implications are briefly discussed.

We report on the observations of 14 dwarf spheroidal galaxies with the Fermi Gamma-Ray Space Telescope taken during the first 11 months of survey mode operations. The Fermi telescope provides a new opportunity to test particle dark matter models through the expected gamma-ray emission produced by pair annihilation of weakly interacting massive particles (WIMPs). Local Group dwarf spheroidal galaxies, the largest galactic substructures predicted by the cold dark matter scenario, are attractive targets for such indirect searches for dark matter because they are nearby and among the most extreme dark matter dominated environments. No significant gamma-ray emission was detected above 100 MeV from the candidate dwarf galaxies. We determine upper limits to the gamma-ray flux assuming both power-law spectra and representative spectra from WIMP annihilation. The resulting integral flux above 100 MeV is constrained to be at a level below around 10^-9 photons cm^-2 s^-1. Using recent stellar kinematic data, the gamma-ray flux limits are combined with improved determinations of the dark matter density profile in 8 of the 14 candidate dwarfs to place limits on the pair annihilation cross-section of WIMPs in several widely studied extensions of the standard model. With the present data, we are able to rule out large parts of the parameter space where the thermal relic density is below the observed cosmological dark matter density and WIMPs (neutralinos here) are dominantly produced non-thermally, e.g. in models where supersymmetry breaking occurs via anomaly mediation. The gamma-ray limits presented here also constrain some WIMP models proposed to explain the Fermi and PAMELA e^+e^- data, including low-mass wino-like neutralinos and models with TeV masses pair-annihilating into muon-antimuon pairs. (Abridged)

Geminga is the second brightest persistent source in the GeV gamma-ray sky. Discovered in 1975 by SAS-2 mission, it was identified as a pulsar only in the 90s, when ROSAT detected the 237 ms X-ray periodicity, that was later also found by EGRET in gamma rays. Even though Geminga has been one of the most intensively studied isolated neutron star during the last 30 years, its interest remains intact especially at gamma-ray energies, where instruments like the Large Area Telescope (LAT) aboard the Fermi mission will provide an unprecedented view of this pulsars. We will report on the preliminary results obtained on the analysis of the first year of observations. We have been able to do precise timing of Geminga using solely gamma rays, producing a timing solution and allowing a deep study of the evolution of the light curve with energy. We have also measured and studied the high-energy cutoff in the phase-averaged spectrum and produced a detailed study of the spectral evolution with phase.

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

We describe a comprehensive pulsar monitoring campaign for the Large Area Telescope (LAT) on the {\em Fermi Gamma-ray Space Telescope} (formerly GLAST). The detection and study of pulsars in gamma rays give insights into the populations of neutron stars and supernova rates in the Galaxy, into particle acceleration mechanisms in neutron star magnetospheres, and into the "engines" driving pulsar wind nebulae. LAT's unprecedented sensitivity between 20 MeV and 300 GeV together with its 2.4 sr field-of-view makes detection of many gamma-ray pulsars likely, justifying the monitoring of over two hundred pulsars with large spin-down powers. To search for gamma-ray pulsations from most of these pulsars requires a set of phase-connected timing solutions spanning a year or more to properly align the sparse photon arrival times. We describe the choice of pulsars and the instruments involved in the campaign. Attention is paid to verifications of the LAT pulsar software, using for example giant radio pulses from the Crab and from PSR B1937+21 recorded at Nan\c{c}ay, and using X-ray data on PSR J0218+4232 from XMM-Newton. We demonstrate accuracy of the pulsar phase calculations at the microsecond level. Data Table 1 is only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ .

We have measured the gamma-ray fluxes of the blazars Mrk 421 and Mrk 501 in the energy range between 50 and 350 GeV (1.2 to 8.3 x 10^25 Hz). The detector, called CELESTE, used first 40, then 53 heliostats of the former solar facility "Themis" in the French Pyrenees to collect Cherenkov light generated in atmospheric particle cascades. The signal from Mrk 421 is often strong. We compare its flux with previously published multi-wavelength studies and infer that we are straddling the high energy peak of the spectral energy distribution. The signal from Mrk 501 in 2000 was weak (3.4 sigma). We obtain an upper limit on the flux from 1ES 1426+428 of less than half that of the Crab flux near 100 GeV. The data analysis and understanding of systematic biases have improved compared to previous work, increasing the detector's sensitivity.

If dark matter is made of neutralinos, annihilation of such Majorana particles should produce high energy cosmic rays, especially in galaxy halo high density regions like galaxy centres. M31 (Andromeda) is our nearest neighbour spiral galaxy, and both its high mass and its low distance make it a source of interest for the indirect search for dark matter through gamma-ray detection. The ground based atmospheric Cherenkov telescope CELESTE observed M31 from 2001 to 2003, in the mostly unexplored energy range 50-500 GeV. These observations provide an upper limit on the flux above 50 GeV around $10^{-10}\rm{cm}^{-2}\rm{s}^{-1}$ in the frame of supersymmetric dark matter, and more generally on any gamma emission from M31.

We have converted the former solar electrical plant THEMIS (French Pyrenees) into an atmospheric Cherenkov detector called CELESTE, which records gamma rays above 30 GeV (7E24 Hz). Here we present the first sub-100 GeV detection by a ground based telescope of a gamma ray source, the Crab nebula, in the energy region between satellite measurements and imaging atmospheric Cherenkov telescopes. At our analysis threshold energy of 60 +/- 20 GeV we measure a gamma ray rate of 6.1 +/- 0.8 per minute. Allowing for 30% systematic uncertainties and a 30% error on the energy scale yields an integral gamma ray flux of I(E>60 GeV) = 6.2^{+5.3}_{-2.3} E-6 photons m^-2 s^-1. The analysis methods used to obtain the gamma ray signal from the raw data are detailed. In addition, we determine the upper limit for pulsed emission to be <12% of the Crab flux at the 99% confidence level, in the same energy range. Our result indicates that if the power law observed by EGRET is attenuated by a cutoff of form e^{-E/E_0} then E_0 < 26 GeV. This is the lowest energy probed by a Cherenkov detector and leaves only a narrow range unexplored beyond the energy range studied by EGRET.

The BL Lac object Mrk 501 went into a very high state of activity during 1997, both in VHE gamma-rays and X-rays. We present here results from observations at energies above 250 GeV carried out between March and October 1997 with the CAT Cerenkov imaging Telescope. The average differential spectrum between 30 GeV and 13 TeV shows significant curvature and is well represented by phi_0 * E_TeV^{-(alpha + beta*log10(E_TeV))}, with: phi_0 = 5.19 +/- 0.13 {stat} +/- 0.12 {sys-MC} +1.66/-1.04 {sys-atm} * 10^-11 /cm^2/s/TeV alpha = 2.24 +/- 0.04 {stat} +/- 0.05 {sys} beta = 0.50 +/- 0.07 {stat} (negligible systematics). The TeV spectral energy distribution of Mrk 501 clearly peaks in the range 500 GeV-1 TeV. Investigation of spectral variations shows a significant hardness-intensity correlation with no measurable effect on the curvature. This can be described as an increase of the peak TeV emission energy with intensity. Simultaneous and quasi-simultaneous CAT VHE gamma-ray and BeppoSAX hard X-ray detections for the highest recorded flare on 16th April and for lower-activity states of the same period show correlated variability with a higher luminosity in X-rays than in gamma-rays. The observed spectral energy distribution and the correlated variability between X-rays and gamma-rays, both in amplitude and in hardening of spectra, favour a two-component emission scheme where the low and high energy components are attributed to synchrotron and inverse Compton (IC) radiation, respectively.

A new method of shower-image analysis is presented which appears very powerful as applied to those Cherenkov Imaging Telescopes with very high definition imaging capability. It provides hadron rejection on the basis of a single cut on the image shape, and simultaneously determines the energy of the electromagnetic shower and the position of the shower axis with respect to the detector. The source location is also reconstructed for each individual gamma-ray shower, even with one single telescope, so for a point source the hadron rejection can be further improved. As an example, this new method is applied to data from the CAT (Cherenkov Array at Themis) imaging telescope, which has been operational since Autumn, 1996.

The CAT (Cherenkov Array at Themis) imaging telescope, equipped with a very-high-definition camera (546 fast phototubes with 0.12 degrees spacing surrounded by 54 larger tubes in two guard rings) started operation in Autumn 1996 on the site of the former solar plant Themis (France). Using the atmospheric Cherenkov technique, it detects and identifies very high energy gamma-rays in the range 250 GeV to a few tens of TeV. The instrument, which has detected three sources (Crab nebula, Mrk 421 and Mrk 501), is described in detail.

The CELESTE experiment will be an Atmospheric Cherenkov detector designed to bridge the gap in energy sensitivity between current satellite and ground-based gamma-ray telescopes, 20 to 300 GeV. We present test results made at the former solar power plant, Themis, in the French Pyrenees. The tests confirm the viability of using a central tower heliostat array for Cherenkov wavefront sampling.

The CAT imaging telesope, which uses the Atmospheric Cherenkov Technique to observe TeV gamma-ray sources, has been operating since September 1996. Located in southern France, it features a fine grained camera consisting of 546 PMT with a pixel size of 0.12 degrees. Strong gamma-ray signals from Mkn 501 above 200 GeV were detected during observations made in April 1997. The high signal/noise ratio for these data enables direct characterization of the gamma-ray images as recorded in the camera. These results are used to validate the Monte-Carlo simulations of the response of the telescope to gamma-ray showers and to improve its performances. Observational results on Mkn 180 and the Crab nebula spectrum are presented.

The CAT Imaging Telescope has observed the BL Lac object Markarian 501 between March and August 1997. We report here on the variability over this time including several large flares. We present also preliminary spectra for all these data, for the low emission state, and for the largest flare.

The VHE gamma-ray imaging telescope CAT started taking data in October 1996. Located at the Themis solar site in southern France, it features a 17.7 m^2 Davies-Cotton mirror equipped with 600 PMT camera at the focal plane. The mechanics and optics, the PMTs and the electronics are presented. The performance, based on the first 7 months of operation, is discribed.

The CAT imaging telescope, recently built on the site of the former solar plant Themis (French Pyrenees), observed gamma-rays from the Crab nebula from October 1996 to March 1997. This steady source, often considered as the standard candle of very-high-energy gamma-ray astronomy, is used as a test-beam to probe the performances of the new telescope, particularly its energy threshold (220 GeV at 20 degrees zenith angle) and the stability of its response. Due to the fine-grain camera, an accurate analysis of the longitudinal profiles of shower images is performed, yielding the source position in two dimensions for each individual shower.

The CAT imaging telescope on the site on the former solar plant Themis has been observing gamma-rays from Mrk501 above 220 GeV in March and April 1997. This source is shown to be highly variable and the light curve is presented. The detected gamma-ray rate for the most intense flare is in excess of 10 per minute.