• ### Search for Spatial Extension in High-Latitude Sources Detected by the Fermi Large Area Telescope(1804.08035)

April 21, 2018 astro-ph.HE
We present a search for spatial extension in high-latitude ($|b|>5^\circ$) sources in the Fermi Large Area Telescope Third Source Catalog (3FGL) and Third Catalog of Hard Sources (3FHL) using the latest Pass 8 software and data set. The result is the Fermi High-Latitude Extended Sources Catalog (FHES) that provides source extensions (or upper limits thereof) and likelihood profiles for the tested source morphologies, i.e., a two-dimensional Gaussian, a disk, and a two-dimensional disk superimposed with a point source (halo). We find 24 extended sources of which 19 sources were not previously characterized as extended. These sources include ones with potential associations with supernova remnants and star forming regions. We also found extended $\gamma$-ray emission in the vicinity of the Cen A radio lobes, and, for the first time at GeV energies, spatially coincident with the radio emission of the SNR CTA 1, and from the Crab Nebula. We also search for halos around active galactic nuclei, which is predicted in the scenario where the $\gamma$ rays that produce the $e^+e^-$ pairs in interactions with background radiation fields induce an electromagnetic cascade. We do not find evidence for extension in individual sources or in stacked source samples. This enables us to place limits on the flux of the extended source components, which are used to derive constraints on the intergalactic magnetic field (IGMF) using spectral and spatial templates of the cascades. This allows us to constrain the IGMF to be stronger than $3\times10^{-16}\,$G for a coherence length $\lambda \gtrsim 10\,$kpc even when conservative assumptions on the source duty cycle are made. This improves previous limits by several orders of magnitude.
• ### Fermi-LAT observations of the LIGO/Virgo event GW170817(1710.05450)

Oct. 16, 2017 astro-ph.HE
We present the Fermi Large Area Telescope (LAT) observations of the binary neutron star merger event GW170817 and the associated short gamma-ray burst (SGRB) GRB\,170817A detected by the Fermi Gamma-ray Burst Monitor. The LAT was entering the South Atlantic Anomaly at the time of the LIGO/Virgo trigger ($t_{\rm GW}$) and therefore cannot place constraints on the existence of high-energy (E $>$ 100 MeV) emission associated with the moment of binary coalescence. We focus instead on constraining high-energy emission on longer timescales. No candidate electromagnetic counterpart was detected by the LAT on timescales of minutes, hours, or days after the LIGO/Virgo detection. The resulting flux upper bound (at 95\% C.L.\/) from the LAT is $4.5\times$10$^{-10}$ erg cm$^{-2}$ s$^{-1}$ in the 0.1--1 GeV range covering a period from T0 + 1153 s to T0 + 2027 s. At the distance of GRB\,170817A, this flux upper bound corresponds to a luminosity upper bound of 9.7$\times10^{43}$ erg s$^{-1}$, which is 5 orders of magnitude less luminous than the only other LAT SGRB with known redshift, GRB\,090510. We also discuss the prospects for LAT detection of electromagnetic counterparts to future gravitational wave events from Advanced LIGO/Virgo in the context of GW170817/GRB\,170817A.
• ### VERITAS and Fermi-LAT observations of TeV gamma-ray sources from the second HAWC catalog(1708.05744)

Aug. 18, 2017 astro-ph.HE
The HAWC observatory recently published their second source catalog with 39 very high energy gamma-ray sources based on 507 days of exposure time. Among these, there were 16 sources that are more than one degree away from any known TeV source. We studied 13 of these 16 sources with VERITAS and \textit{Fermi}-LAT data. VERITAS, an array of four imaging atmospheric Cherenkov telescopes observing gamma rays with energies higher than 85 GeV, can provide a more detailed image of the source with much shorter exposure time and with better angular resolution. With \textit{Fermi}-LAT data, we searched for the counterparts at lower energies (E>10 GeV). The exposure of VERITAS varies among the 13 sources, and we used eight years of \textit{Fermi}-LAT data. VERITAS found weak gamma-ray emission in the region of PWN~DA~495 coinciding with 2HWC J1953+294 in this follow-up study. We will discuss the results especially focusing on the PWN~DA~495 region and the SNR~G~54.1+0.3 region where \textit{Fermi}-LAT detected a GeV counterpart of SNR~G~54.1+0.3, a known TeV source detected by both VERITAS and HAWC.
• ### Characterizing the population of pulsars in the Galactic bulge with the $Fermi$ Large Area Telescope(1705.00009)

April 28, 2017 astro-ph.HE
An excess of $\gamma$-ray emission from the Galactic Center (GC) region with respect to predictions based on a variety of interstellar emission models and $\gamma$-ray source catalogs has been found by many groups using data from the $Fermi$ Large Area Telescope (LAT). Several interpretations of this excess have been invoked. In this paper we test the interpretation that the excess is caused by an unresolved population of $\gamma$-ray pulsars located in the Galactic bulge. We use cataloged LAT sources to derive criteria that efficiently select pulsars with very small contamination from blazars. We search for point sources in the inner $40^\circ\times40^\circ$ region of the Galaxy, derive a list of approximately 400 sources, and apply pulsar selection criteria to extract pulsar candidates among our source list. We also derive the efficiency of these selection criteria for $\gamma$-ray pulsars as a function of source energy flux and location. We demonstrate that given the observed spatial and flux distribution of pulsar candidates, a model that includes a population with about 2.7 $\gamma$-ray pulsars in the Galactic disk (in our $40^\circ\times40^\circ$ analysis region) for each pulsar in the Galactic bulge is preferred at the level of 7~standard deviations with respect to a disk-only model. The properties of these disk and bulge pulsar populations are consistent with the population of known $\gamma$-ray pulsars as well as with the spatial profile and energy spectrum of the GC excess. Finally, we show that the dark matter interpretation of the GC excess is strongly disfavored since a distribution of dark matter is not able to mimic the observed properties of the population of sources detected in our analysis.
• ### Search for Cosmic-Ray Electron and Positron Anisotropies with Seven Years of Fermi Large Area Telescope Data(1703.01073)

March 3, 2017 astro-ph.HE
The Large Area Telescope on board the Fermi Gamma-ray Space Telescope has collected the largest ever sample of high-energy cosmic-ray electron and positron events since the beginning of its operation. Potential anisotropies in the arrival directions of cosmic-ray electrons or positrons could be a signature of the presence of nearby sources. We use almost seven years of data with energies above 42 GeV processed with the Pass~8 reconstruction. The present data sample can probe dipole anisotropies down to a level of $10^{-3}$. We take into account systematic effects that could mimic true anisotropies at this level. We present a detailed study of the event selection optimization of the cosmic-ray electrons and positrons to be used for anisotropy searches. Since no significant anisotropies have been detected on any angular scale, we present upper limits on the dipole anisotropy. The present constraints are among the strongest to date probing the presence of nearby young and middle-aged sources.
• ### Observations of M31 and M33 with the Fermi Large Area Telescope: a galactic center excess in Andromeda?(1702.08602)

Feb. 28, 2017 astro-ph.HE
The Fermi LAT has opened the way for comparative studies of cosmic rays (CRs) and high-energy objects in the Milky Way (MW) and in other, external, star-forming galaxies. Using 2 yr of observations with the Fermi LAT, local Group galaxy M31 was detected as a marginally extended gamma-ray source, while only an upper limit (UL) has been derived for the other nearby galaxy M33. We revisited the gamma-ray emission in the direction of M31 and M33 using more than 7 yr of LAT Pass 8 data in the energy range 0.1-100 GeV, presenting detailed morphological and spectral analyses. M33 remains undetected and we computed an UL for it. This revised UL remains consistent with the observed correlation between gamma-ray luminosity and star-formation rate tracers and implies an average CR density in M33 that is at most half of that of the MW. M31 is detected with a significance of nearly 10 sigma and to be extended with 4 sigma. Its spectrum is consistent with a power law. The spatial distribution of the emission is consistent with a uniform disk with a radius of 0.4 deg and no offset from the center of M31, but nonuniform intensity distributions cannot be excluded. The flux from M31 appears confined to the inner regions of the galaxy and does not fill the disk or extend far from it. The gamma-ray signal is not correlated with regions rich in gas or star-formation activity suggesting that the emission is not interstellar in origin, unless the energetic particles radiating in gamma rays do not originate in recent star formation. Alternative and nonexclusive interpretations are that the emission results from a population of millisecond pulsars dispersed in the bulge and disk of M31 by disrupted globular clusters or from the decay or annihilation of dark matter particles, similar to what has been proposed to account for the so-called Galactic Center excess found in Fermi-LAT observations of the MW.
• The supernova remnant (SNR) W49B originated from a core-collapse supernova that occurred between one and four thousand years ago, and subsequently evolved into a mixed-morphology remnant, which is interacting with molecular clouds (MC). $\gamma$-ray observations of SNR/MC associations are a powerful tool to constrain the origin of Galactic cosmic-rays, as they can probe the acceleration of hadrons through their interaction with the surrounding medium and subsequent emission of non-thermal photons. The detection of a $\gamma$-ray source coincident with W49B at very high energies (VHE; E > 100 GeV) with the H.E.S.S. Cherenkov telescopes is reported together with a study of the source with 5 years of Fermi-LAT high energy $\gamma$-ray (0.06 - 300 GeV) data. The smoothly-connected combined source spectrum, measured from 60 MeV to multi-TeV energies, shows two significant spectral breaks at $304\pm20$ MeV and $8.4_{-2.5}^{+2.2}$ GeV, the latter being constrained by the joint fit from the two instruments. The detected spectral features are similar to those observed in several other SNR/MC associations and are found to be indicative of $\gamma$-ray emission produced through neutral-pion decay.
• ### Fermi-LAT Observations of the LIGO event GW150914(1602.04488)

Feb. 23, 2016 astro-ph.HE
The Fermi Large Area Telescope (LAT) has an instantaneous field of view covering $\sim 1/5$ of the sky and completes a survey of the full sky every ~3 hours. It provides a continuous, all-sky survey of high-energy gamma-rays, enabling searches for transient phenomena over timescales from milliseconds to years. Among these phenomena could be electromagnetic counterparts to gravitational wave sources. In this paper, we present a detailed study of the LAT observations relevant to Laser Interferometer Gravitational-wave Observatory (LIGO) event GW150904 (Abbott et al. 2016), which is the first direct detection of gravitational waves and has been interpreted as due to coalescence of two stellar-mass black holes. The localization region for GW150904 was outside the LAT field of view at the time of the gravitational wave signal. However, as part of routine survey observations, the LAT observed the entire LIGO localization region within ~70 minutes of the trigger, and thus enabled a comprehensive search for a gamma-ray counterpart to GW150904. The study of the LAT data presented here did not find any potential counterparts to GW150904, but it did provide limits on the presence of a transient counterpart above 100 MeV on timescales of hours to days over the entire GW150904 localization region.
• ### Searching for Dark Matter Annihilation from Milky Way Dwarf Spheroidal Galaxies with Six Years of Fermi-LAT Data(1503.02641)

Nov. 3, 2015 hep-ex, astro-ph.HE
The dwarf spheroidal satellite galaxies (dSphs) of the Milky Way are some of the most dark matter (DM) dominated objects known. We report on gamma-ray observations of Milky Way dSphs based on 6 years of Fermi Large Area Telescope data processed with the new Pass 8 event-level analysis. None of the dSphs are significantly detected in gamma rays, and we present upper limits on the DM annihilation cross section from a combined analysis of 15 dSphs. These constraints are among the strongest and most robust to date and lie below the canonical thermal relic cross section for DM of mass $\lesssim$ 100 GeV annihilating via quark and $\tau$-lepton channels.
• ### The Spectrum and Morphology of the Fermi Bubbles(1407.7905)

July 29, 2014 astro-ph.HE
The Fermi bubbles are two large structures in the gamma-ray sky extending to $55^\circ$ above and below the Galactic center. We analyze 50 months of Fermi Large Area Telescope data between 100 MeV and 500 GeV above $10^\circ$ in Galactic latitude to derive the spectrum and morphology of the Fermi bubbles. We thoroughly explore the systematic uncertainties that arise when modeling the Galactic diffuse emission through two separate approaches. The gamma-ray spectrum is well described by either a log parabola or a power law with an exponential cutoff. We exclude a simple power law with more than 7$\sigma$ significance. The power law with an exponential cutoff has an index of $1.9 \pm 0.2$ and a cutoff energy of $110\pm 50$ GeV. We find that the gamma-ray luminosity of the bubbles is $4.4^{+2.4}_{-0.9} \times 10^{37}$ erg s$^{-1}$. We confirm a significant enhancement of gamma-ray emission in the south-eastern part of the bubbles, but we do not find significant evidence for a jet. No significant variation of the spectrum across the bubbles is detected. The width of the boundary of the bubbles is estimated to be $3.4^{+3.7}_{-2.6}$ deg. Both inverse Compton (IC) models and hadronic models including IC emission from secondary leptons fit the gamma-ray data well. In the IC scenario, the synchrotron emission from the same population of electrons can also explain the WMAP and Planck microwave haze with a magnetic field between 5 and 20 $\mu$G.
• Previous observations with HESS have revealed the existence of an extended very-high-energy (VHE; E>100 GeV) gamma-ray source, HESS J1834-087, coincident with the SNR W41. The origin of the gamma-ray emission has been further investigated with HESS and the Fermi-LAT. The gamma-ray data provided by 61h (HESS) and 4 yrs (Fermi LAT) of observations cover over 5 decades in energy (1.8GeV - 30TeV). The morphology and spectrum of the TeV and GeV sources have been studied and multi-wavelength data have been used to investigate the origin of the observed emission. The TeV source can be modeled with a sum of two components: one point-like and one significantly extended (sig_TeV = 0.17{\deg}), both centered on SNR W41 and exhibiting spectra described by a power law of index 2.6. The GeV source detected with Fermi is extended (sig_GeV =0.15{\deg}) and morphologically matches the VHE emission. Its spectrum can be described by a power-law with index 2.15 and joins smoothly the one of the whole TeV source. A break appears in the spectra around 100 GeV. Two main scenarios are proposed to explain the emission: a pulsar wind nebula (PWN) or the interaction of SNR W41 with a molecular cloud. X-ray observations suggest the presence of a point-like source (pulsar candidate) near the center of the SNR and non-thermal X-ray diffuse emission which could arise from a potential PWN. The PWN scenario is supported by the match of of the TeV and GeV positions with the putative pulsar. However, the overall spectrum is reproduced by a 1-zone leptonic model only if an excess of low-energy electrons is injected by a high spin-down power pulsar. This low-energy component is not needed if the point-like TeV source is unrelated to the extended GeV and TeV sources. The interacting SNR scenario is supported by the spatial coincidence between the gamma-ray sources, the detection of OH maser lines and the hadronic modeling.
• ### Fermi-LAT Observations of the Gamma-ray Burst GRB 130427A(1311.5623)

Dec. 9, 2013 astro-ph.HE
The observations of the exceptionally bright gamma-ray burst (GRB) 130427A by the Large Area Telescope aboard the Fermi Gamma-ray Space Telescope provide constraints on the nature of such unique astrophysical sources. GRB 130427A had the largest fluence, highest-energy photon (95 GeV), longest $\gamma$-ray duration (20 hours), and one of the largest isotropic energy releases ever observed from a GRB. Temporal and spectral analyses of GRB 130427A challenge the widely accepted model that the non-thermal high-energy emission in the afterglow phase of GRBs is synchrotron emission radiated by electrons accelerated at an external shock.
• ### Search for Gamma-ray Spectral Lines with the Fermi Large Area Telescope and Dark Matter Implications(1305.5597)

Sept. 23, 2013 astro-ph.HE
Weakly Interacting Massive Particles (WIMPs) are a theoretical class of particles that are excellent dark matter candidates. WIMP annihilation or decay may produce essentially monochromatic gamma rays detectable by the Fermi Large Area Telescope (LAT) against the astrophysical gamma-ray emission of the Galaxy. We have searched for spectral lines in the energy range 5--300 GeV using 3.7 years of data, reprocessed with updated instrument calibrations and an improved energy dispersion model compared to the previous Fermi-LAT Collaboration line searches. We searched in five regions selected to optimize sensitivity to different theoretically-motivated dark matter density distributions. We did not find any globally significant lines in our a priori search regions and present 95% confidence limits for annihilation cross sections of self-conjugate WIMPs and decay lifetimes. Our most significant fit occurred at 133 GeV in our smallest search region and had a local significance of 3.3 standard deviations, which translates to a global significance of 1.5 standard deviations. We discuss potential systematic effects in this search, and examine the feature at 133 GeV in detail. We find that both the use of reprocessed data and of additional information in the energy dispersion model contribute to the reduction in significance of the line-like feature near 130 GeV relative to significances reported in other works. We also find that the feature is narrower than the LAT energy resolution at the level of 2 to 3 standard deviations, which somewhat disfavors the interpretation of the 133 GeV feature as a real WIMP signal.
• ### Determination of the Point-Spread Function for the Fermi Large Area Telescope from On-orbit Data and Limits on Pair Halos of Active Galactic Nuclei(1309.5416)

Sept. 21, 2013 astro-ph.HE
The Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope is a pair-conversion telescope designed to detect photons with energies from 20 MeV to > 300 GeV. The pre-launch response functions of the LAT were determined through extensive Monte Carlo simulations and beam tests. The point-spread function (PSF) characterizing the angular distribution of reconstructed photons as a function of energy and geometry in the detector is determined here from two years of on-orbit data by examining the distributions of \gamma rays from pulsars and active galactic nuclei (AGN). Above 3 GeV, the PSF is found to be broader than the pre-launch PSF. We checked for dependence of the PSF on the class of \gamma-ray source and observation epoch and found none. We also investigated several possible spatial models for pair-halo emission around BL Lac AGN. We found no evidence for a component with spatial extension larger than the PSF and set upper limits on the amplitude of halo emission in stacked images of low and high redshift BL Lac AGN and the TeV blazars 1ES0229+200 and 1ES0347-121.
• ### The Fermi All-sky Variability Analysis: A list of flaring gamma-ray sources and the search for transients in our Galaxy(1304.6082)

April 22, 2013 astro-ph.HE
In this paper we present the Fermi All-sky Variability Analysis (FAVA), a tool to systematically study the variability of the gamma-ray sky measured by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. For each direction on the sky, FAVA compares the number of gamma rays observed in a given time window to the number of gamma rays expected for the average emission detected from that direction. This method is used in weekly time intervals to derive a list of 215 flaring gamma-ray sources. We proceed to discuss the 27 sources found at Galactic latitudes smaller than 10 degrees and show that, despite their low latitudes, most of them are likely of extragalactic origin.
• ### Impulsive and Long Duration High-Energy Gamma-ray Emission From the Very Bright 2012 March 7 Solar Flares(1304.5559)

April 19, 2013 astro-ph.HE
The Fermi Large Area Telescope (LAT) observed two bright X-class solar flares on 2012 March 7, and detected gamma-rays up to 4 GeV. We detected gamma-rays both during the impulsive and temporally-extended emission phases, with emission above 100 MeV lasting for approximately 20 hours. Accurate localization of the gamma-ray production site(s) coincide with the solar active region from which X-ray emissions associated with these flares originated. Our analysis of the >100 MeV gamma-ray emission shows a relatively rapid monotonic decrease in flux during the first hour of the impulsive phase, and a much slower, almost monotonic decrease in flux for the next 20 hours. The spectra can be adequately described by a power law with a high energy exponential cutoff, or as resulting from the decay of neutral pions produced by accelerated protons and ions with an isotropic power-law energy distribution. The required proton spectrum has a number index ~3, with minor variations during the impulsive phase, while during the temporally extended phase the spectrum softens monotonically, starting with index ~4. The >30 MeV proton flux and spectra observed near the Earth by the GOES satellites also show a monotonic flux decrease and spectral softening during the extended phase, but with a harder spectrum, with index ~3. Based on the Fermi-LAT and GOES observations of the flux and spectral evolution of these bright flares, we explore the relative merits of prompt and continuous acceleration scenarios, hadronic and leptonic emission processes, and acceleration at the solar corona by the fast Coronal Mass Ejections (CME) as explanations for the observations. We conclude that the most likely scenario is continuous acceleration of protons in the solar corona which penetrate the lower solar atmosphere and produce pions that decay into gamma-rays.
• ### High-Energy Gamma-Ray Emission From Solar Flares: Summary of Fermi LAT Detections and Analysis of Two M-Class Flares(1304.3749)

April 12, 2013 astro-ph.HE
We present the detections of 19 solar flares detected in high-energy gamma rays (above 100 MeV) with the Fermi Large Area Telescope (LAT) during its first four years of operation. Interestingly, all flares are associated with fairly fast Coronal Mass Ejections (CMEs) and are not all powerful X-ray flares. We then describe the detailed temporal, spatial and spectral characteristics of the first two long-lasting events: the 2011 March 7 flare, a moderate (M3.7) impulsive flare followed by slowly varying gamma-ray emission over 13 hours, and the 2011 June 7 M2.5 flare, which was followed by gamma-ray emission lasting for 2 hours. We compare the Fermi-LAT data with X-ray and proton data measurements from GOES and RHESSI. We argue that a hadronic origin of the gamma rays is more likely than a leptonic origin and find that the energy spectrum of the proton distribution softens after the 2011 March 7 flare, favoring a scenario with continuous acceleration at the flare site. This work suggests that proton acceleration in solar flares is more common than previously thought, occurring for even modest X-ray flares, and for longer durations.
• ### The First Fermi LAT Gamma-Ray Burst Catalog(1303.2908)

March 12, 2013 astro-ph.HE
In three years of observations since the beginning of nominal science operations in August 2008, the Large Area Telescope (LAT) on board the Fermi Gamma Ray Space Telescope has observed high-energy (>20 MeV) \gamma-ray emission from 35 gamma-ray bursts (GRBs). Among these, 28 GRBs have been detected above 100 MeV and 7 GRBs above ~ 20 MeV. The first Fermi-LAT catalog of GRBs is a compilation of these detections and provides a systematic study of high-energy emission from GRBs for the first time. To generate the catalog, we examined 733 GRBs detected by the Gamma-Ray Burst Monitor (GBM) on Fermi and processed each of them using the same analysis sequence. Details of the methodology followed by the LAT collaboration for GRB analysis are provided. We summarize the temporal and spectral properties of the LAT-detected GRBs. We also discuss characteristics of LAT-detected emission such as its delayed onset and longer duration compared to emission detected by the GBM, its power-law temporal decay at late times, and the fact that it is dominated by a power-law spectral component that appears in addition to the usual Band model.
• ### The Fermi Large Area Telescope On Orbit: Event Classification, Instrument Response Functions, and Calibration(1206.1896)

The Fermi Large Area Telescope (Fermi-LAT, hereafter LAT), the primary instrument on the Fermi Gamma-ray Space Telescope (Fermi) mission, is an imaging, wide field-of-view, high-energy \gamma-ray telescope, covering the energy range from 20 MeV to more than 300 GeV. During the first years of the mission the LAT team has gained considerable insight into the in-flight performance of the instrument. Accordingly, we have updated the analysis used to reduce LAT data for public release as well as the Instrument Response Functions (IRFs), the description of the instrument performance provided for data analysis. In this paper we describe the effects that motivated these updates. Furthermore, we discuss how we originally derived IRFs from Monte Carlo simulations and later corrected those IRFs for discrepancies observed between flight and simulated data. We also give details of the validations performed using flight data and quantify the residual uncertainties in the IRFs. Finally, we describe techniques the LAT team has developed to propagate those uncertainties into estimates of the systematic errors on common measurements such as fluxes and spectra of astrophysical sources.
• ### Gamma-ray flares from the Crab Nebula(1011.3855)

Jan. 10, 2011 astro-ph.HE
A young and energetic pulsar powers the well-known Crab Nebula. Here we describe two separate gamma-ray (photon energy >100 MeV) flares from this source detected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. The first flare occurred in February 2009 and lasted approximately 16 days. The second flare was detected in September 2010 and lasted approximately 4 days. During these outbursts the gamma-ray flux from the nebula increased by factors of four and six, respectively. The brevity of the flares implies that the gamma rays were emitted via synchrotron radiation from PeV (10^15 eV) electrons in a region smaller than 1.4 10^-2 pc. These are the highest energy particles that can be associated with a discrete astronomical source, and they pose challenges to particle acceleration theory.
• ### Fermi Large Area Telescope Observations of Misaligned AGN(1007.1624)

July 9, 2010 astro-ph.HE
We report the detection of high energy gamma-ray emission from the young and energetic pulsar PSR B1509$-$58 and its pulsar wind nebula (PWN) in the composite supernova remnant SNR G320.4-1.2 (aka MSH 15-52). Using 1 year of survey data with the Fermi-Large Area Telescope (LAT), we detected pulsations from PSR B1509-58 up to 1 GeV and extended gamma-ray emission above 1 GeV spatially coincident with the PWN. The pulsar light curve presents two peaks offset from the radio peak by phases 0.96 $\pm$ 0.01 and 0.33 $\pm$ 0.02. New constraining upper limits on the pulsar emission are derived below 1 GeV and confirm a severe spectral break at a few tens of MeV. The nebular spectrum in the 1 - 100 GeV energy range is well described by a power-law with a spectral index of (1.57 $\pm$ 0.17 $\pm$ 0.13) and a flux above 1 GeV of (2.91 $\pm$ 0.79 $\pm$ 1.35) 10^{-9} cm^{-2} s^{-1}. The first errors represent the statistical errors on the fit parameters, while the second ones are the systematic uncertainties. The LAT spectrum of the nebula connects nicely with Cherenkov observations, and indicates a spectral break between GeV and TeV energies.