The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope routinely detects the highly dust-absorbed, reddened, and MeV-peaked flat spectrum radio quasar PKS 1830-211 (z=2.507). Its apparent isotropic gamma-ray luminosity (E>100 MeV) averaged over $\sim$ 3 years of observations and peaking on 2010 October 14/15 at 2.9 X 10^{50} erg s^{-1}, makes it among the brightest high-redshift Fermi blazars. No published model with a single lens can account for all of the observed characteristics of this complex system. Based on radio observations, one expects time delayed variability to follow about 25 days after a primary flare, with flux about a factor 1.5 less. Two large gamma-ray flares of PKS 1830-211 have been detected by the LAT in the considered period and no substantial evidence for such a delayed activity was found. This allows us to place a lower limit of about 6 on the gamma rays flux ratio between the two lensed images. Swift XRT observations from a dedicated Target of Opportunity program indicate a hard spectrum and with no significant correlation of X-ray flux with the gamma-ray variability. The spectral energy distribution can be modeled with inverse Compton scattering of thermal photons from the dusty torus. The implications of the LAT data in terms of variability, the lack of evident delayed flare events, and different radio and gamma-ray flux ratios are discussed. Microlensing effects, absorption, size and location of the emitting regions, the complex mass distribution of the system, an energy-dependent inner structure of the source, and flux suppression by the lens galaxy for one image path may be considered as hypotheses for understanding our results.

Current theories predict relativistic hadronic particle populations in clusters of galaxies in addition to the already observed relativistic leptons. In these scenarios hadronic interactions give rise to neutral pions which decay into $\gamma$ rays, that are potentially observable with the Large Area Telescope (LAT) on board the Fermi space telescope. We present a joint likelihood analysis searching for spatially extended $\gamma$-ray emission at the locations of 50 galaxy clusters in 4 years of Fermi-LAT data under the assumption of the universal cosmic-ray model proposed by Pinzke & Pfrommer (2010). We find an excess at a significance of $2.7\sigma$ which upon closer inspection is however correlated to individual excess emission towards three galaxy clusters: Abell 400, Abell 1367 and Abell 3112. We discuss these cases in detail and conservatively attribute the emission to unmodeled background (for example, radio galaxies within the clusters). Through the combined analysis of 50 clusters we exclude hadronic injection efficiencies in simple hadronic models above 21% and establish limits on the cosmic-ray to thermal pressure ratio within the virial radius, $R_{200}$, to be below 1.2-1.4% depending on the morphological classification. In addition we derive new limits on the $\gamma$-ray flux from individual clusters in our sample.

Recent accurate measurements of cosmic-ray (CR) species by ATIC-2, CREAM, and PAMELA reveal an unexpected hardening in the proton and He spectra above a few hundred GeV, a gradual softening of the spectra just below a few hundred GeV, and a harder spectrum of He compared to that of protons. These newly-discovered features may offer a clue to the origin of high-energy CRs. We use the ${\it Fermi}$ Large Area Telescope observations of the $\gamma$-ray emission from the Earth's limb for an indirect measurement of the local spectrum of CR protons in the energy range $\sim 90~$GeV-$6~$TeV (derived from a photon energy range $15~$GeV-$1~$TeV). Our analysis shows that single power law and broken power law spectra fit the data equally well and yield a proton spectrum with index $2.68 \pm 0.04$ and $2.61 \pm 0.08$ above $\sim 200~$GeV, respectively.

The dwarf spheroidal satellite galaxies of the Milky Way are some of the most dark-matter-dominated objects known. Due to their proximity, high dark matter content, and lack of astrophysical backgrounds, dwarf spheroidal galaxies are widely considered to be among the most promising targets for the indirect detection of dark matter via gamma rays. Here we report on gamma-ray observations of 25 Milky Way dwarf spheroidal satellite galaxies based on 4 years of Fermi Large Area Telescope (LAT) data. None of the dwarf galaxies are significantly detected in gamma rays, and we present gamma-ray flux upper limits between 500 MeV and 500 GeV. We determine the dark matter content of 18 dwarf spheroidal galaxies from stellar kinematic data and combine LAT observations of 15 dwarf galaxies to constrain the dark matter annihilation cross section. We set some of the tightest constraints to date on the the annihilation of dark matter particles with masses between 2 GeV and 10 TeV into prototypical Standard Model channels. We find these results to be robust against systematic uncertainties in the LAT instrument performance, diffuse gamma-ray background modeling, and assumed dark matter density profile.

Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial pulse up to 2.5 s is possibly the brightest well-isolated pulse observed to date. A fine time resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly correlated power-law behavior is observed between the luminosity and the spectral peak energy that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the existing models to account for all of the observed spectral and temporal behaviors simultaneously.

Based on the experience gained during the four and a half years of the mission, the Fermi -LAT collaboration has undertaken a comprehensive revision of the event-level analysis going under the name of Pass 8. Although it is not yet finalized, we can test the improvements in the new event reconstruction with the special case of the prompt phase of bright Gamma-Ray Bursts (GRBs), where the signal to noise ratio is large enough that loose selection cuts are sufficient to identify gamma- rays associated with the source. Using the new event reconstruction, we have re-analyzed ten GRBs previously detected by the LAT for which an x-ray/optical follow-up was possible and found four new gamma rays with energies greater than 10 GeV in addition to the seven previously known. Among these four is a 27.4 GeV gamma-ray from GRB 080916C, which has a redshift of 4.35, thus making it the gamma ray with the highest intrinsic energy (147 GeV) detected from a GRB. We present here the salient aspects of the new event reconstruction and discuss the scientific implications of these new high-energy gamma rays, such as constraining extragalactic background light models, Lorentz invariance violation (LIV) tests, the prompt emission mechanism and the bulk Lorentz factor of the emitting region.

Cosmic rays are particles (mostly protons) accelerated to relativistic speeds. Despite wide agreement that supernova remnants (SNRs) are the sources of galactic cosmic rays, unequivocal evidence for the acceleration of protons in these objects is still lacking. When accelerated protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma rays. This offers a compelling way to detect the acceleration sites of protons. The identification of pion-decay gamma rays has been difficult because high-energy electrons also produce gamma rays via bremsstrahlung and inverse Compton scattering. We detected the characteristic pion-decay feature in the gamma-ray spectra of two SNRs, IC 443 and W44, with the Fermi Large Area Telescope. This detection provides direct evidence that cosmic-ray protons are accelerated in SNRs.

We have performed an analysis of the diffuse gamma-ray emission with the Fermi Large Area Telescope in the Milky Way Halo region searching for a signal from dark matter annihilation or decay. In the absence of a robust dark matter signal, constraints are presented. We consider both gamma rays produced directly in the dark matter annihilation/decay and produced by inverse Compton scattering of the e+e- produced in the annihilation/decay. Conservative limits are derived requiring that the dark matter signal does not exceed the observed diffuse gamma-ray emission. A second set of more stringent limits is derived based on modeling the foreground astrophysical diffuse emission using the GALPROP code. Uncertainties in the height of the diffusive cosmic-ray halo, the distribution of the cosmic-ray sources in the Galaxy, the index of the injection cosmic-ray electron spectrum and the column density of the interstellar gas are taken into account using a profile likelihood formalism, while the parameters governing the cosmic-ray propagation have been derived from fits to local cosmic-ray data. The resulting limits impact the range of particle masses over which dark matter thermal production in the early Universe is possible, and challenge the interpretation of the PAMELA/Fermi-LAT cosmic ray anomalies as annihilation of dark matter.

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 blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

We report on the gamma-ray observations of giant molecular clouds Orion A and B with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. The gamma-ray emission in the energy band between \sim100 MeV and \sim100 GeV is predicted to trace the gas mass distribution in the clouds through nuclear interactions between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production cross-section for the nuclear interaction is known to \sim10% precision which makes the LAT a powerful tool to measure the gas mass column density distribution of molecular clouds for a known CR intensity. We present here such distributions for Orion A and B, and correlate them with those of the velocity integrated CO intensity (WCO) at a 1{\deg} \times1{\deg} pixel level. The correlation is found to be linear over a WCO range of ~10 fold when divided in 3 regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The Wco-to-mass conversion factor, Xco, is found to be \sim2.3\times10^20 cm-2(K km s-1)-1 for the high-longitude part of Orion A (l > 212{\deg}), \sim1.7 times higher than \sim1.3 \times 10^20 found for the rest of Orion A and B. We interpret the apparent high Xco in the high-longitude region of Orion A in the light of recent works proposing a non-linear relation between H2 and CO densities in the diffuse molecular gas. Wco decreases faster than the H2 column density in the region making the gas "darker" to Wco.

We measured separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope. Because the instrument does not have an onboard magnet, we distinguish the two species by exploiting the Earth's shadow, which is offset in opposite directions for opposite charges due to the Earth's magnetic field. We estimate and subtract the cosmic-ray proton background using two different methods that produce consistent results. We report the electron-only spectrum, the positron-only spectrum, and the positron fraction between 20 GeV and 200 GeV. We confirm that the fraction rises with energy in the 20-100 GeV range. The three new spectral points between 100 and 200 GeV are consistent with a fraction that is continuing to rise with energy.

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.

Numerical simulations based on the Lambda-CDM model of cosmology predict a large number of as yet unobserved Galactic dark matter satellites. We report the results of a Large Area Telescope (LAT) search for these satellites via the gamma-ray emission expected from the annihilation of weakly interacting massive particle (WIMP) dark matter. Some dark matter satellites are expected to have hard gamma-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We sought to identify LAT sources with these characteristics, focusing on gamma-ray spectra consistent with WIMP annihilation through the $b \bar b$ channel. We found no viable dark matter satellite candidates using one year of data, and we present a framework for interpreting this result in the context of numerical simulations to constrain the velocity-averaged annihilation cross section for a conventional 100 GeV WIMP annihilating through the $b \bar b$ channel.

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.

We investigate the Fermi LAT gamma-ray and 15 GHz VLBA radio properties of a joint gamma-ray- and radio-selected sample of AGNs obtained during the first 11 months of the Fermi mission (2008 Aug 4 - 2009 Jul 5). Our sample contains the brightest 173 AGNs in these bands above declination -30 deg. during this period, and thus probes the full range of gamma-ray loudness (gamma-ray to radio band luminosity ratio) in the bright blazar population. The latter quantity spans at least four orders of magnitude, reflecting a wide range of spectral energy distribution (SED) parameters in the bright blazar population. The BL Lac objects, however, display a linear correlation of increasing gamma-ray loudness with synchrotron SED peak frequency, suggesting a universal SED shape for objects of this class. The synchrotron self-Compton model is favored for the gamma-ray emission in these BL Lacs over external seed photon models, since the latter predict a dependence of Compton dominance on Doppler factor that would destroy any observed synchrotron SED peak - gamma-ray loudness correlation. The high-synchrotron peaked (HSP) BL Lac objects are distinguished by lower than average radio core brightness temperatures, and none display large radio modulation indices or high linear core polarization levels. No equivalent trends are seen for the flat-spectrum radio quasars (FSRQ) in our sample. Given the association of such properties with relativistic beaming, we suggest that the HSP BL Lacs have generally lower Doppler factors than the lower-synchrotron peaked BL Lacs or FSRQs in our sample.

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.