We present the third Fermi Large Area Telescope (LAT) source catalog (3FGL) of sources in the 100 MeV-300 GeV range. Based on the first 4 yr of science data from the Fermi Gamma-ray Space Telescope mission, it is the deepest yet in this energy range. Relative to the Second Fermi LAT catalog, the 3FGL catalog incorporates twice as much data, as well as a number of analysis improvements, including improved calibrations at the event reconstruction level, an updated model for Galactic diffuse γ-ray emission, a refined procedure for source detection, and improved methods for associating LAT sources with potential counterparts at other wavelengths. The 3FGL catalog includes 3033 sources above 4σ significance, with source location regions, spectral properties, and monthly light curves for each. Of these, 78 are flagged as potentially being due to imperfections in the model for Galactic diffuse emission. Twenty-five sources are modeled explicitly as spatially extended, and overall 238 sources are considered as identified based on angular extent or correlated variability (periodic or otherwise) observed at other wavelengths. For 1010 sources we have not found plausible counterparts at other wavelengths. More than 1100 of the identified or associated sources are active galaxies of the blazar class; several other classes of non-blazar active galaxies are also represented in the 3FGL. Pulsars represent the largest Galactic source class. From source counts of Galactic sources we estimate that the contribution of unresolved sources to the Galactic diffuse emission is ∼3% at 1 GeV.

We present the second catalog of high-energy γ -ray sources detected by the Large Area Telescope (LAT), the primary science instrument on the Fermi Gamma-ray Space Telescope (Fermi), derived from data taken during the first 24 months of the science phase of the mission, which began on 2008 August 4. Source detection is based on the average flux over the 24 month period. The second Fermi-LAT catalog (2FGL) includes source location regions, defined in terms of elliptical fits to the 95% confidence regions and spectral fits in terms of power-law, exponentially cutoff power-law, or log-normal forms. Also included are flux measurements in five energy bands and light curves on monthly intervals for each source. Twelve sources in the catalog are modeled as spatially extended. We provide a detailed comparison of the results from this catalog with those from the first Fermi-LAT catalog (1FGL). Although the diffuse Galactic and isotropic models used in the 2FGL analysis are improved compared to the 1FGL catalog, we attach caution flags to 162 of the sources to indicate possible confusion with residual imperfections in the diffuse model. The 2FGL catalog contains 1873 sources detected and characterized in the 100 MeV to 100 GeV range of which we consider 127 as being firmly identified and 1171 as being reliably associated with counterparts of known or likely γ -ray-producing source classes.

We present a catalog of high-energy gamma-ray sources detected by the Large Area Telescope (LAT), the primary science instrument on the Fermi Gamma-ray Space Telescope (Fermi), during the first 11 months of the science phase of the mission, which began on 2008 August 4. The First Fermi-LAT catalog (1FGL) contains 1451 sources detected and characterized in the 100 MeV to 100 GeV range. Source detection was based on the average flux over the 11 month period, and the threshold likelihood Test Statistic is 25, corresponding to a significance of just over 4σ. The 1FGL catalog includes source location regions, defined in terms of elliptical fits to the 95% confidence regions and power-law spectral fits as well as flux measurements in five energy bands for each source. In addition, monthly light curves are provided. Using a protocol defined before launch we have tested for several populations of gamma-ray sources among the sources in the catalog. For individual LAT-detected sources we provide firm identifications or plausible associations with sources in other astronomical catalogs. Identifications are based on correlated variability with counterparts at other wavelengths, or on spin or orbital periodicity. For the catalogs and association criteria that we have selected, 630 of the sources are unassociated. Care was taken to characterize the sensitivity of the results to the model of interstellar diffuse gamma-ray emission used to model the bright foreground, with the result that 161 sources at low Galactic latitudes and toward bright local interstellar clouds are flagged as having properties that are strongly dependent on the model or as potentially being due to incorrectly modeled structure in the Galactic diffuse emission.

The third catalog of active galactic nuclei (AGNs) detected by the Fermi-LAT (3LAC) is presented. It is based on the third Fermi-LAT catalog (3FGL) of sources detected between 100 MeV and 300 GeV with a Test Statistic greater than 25, between 2008 August 4 and 2012 July 31. The 3LAC includes 1591 AGNs located at high Galactic latitudes (>  | | b 10), a 71% increase over the second catalog based on 2 years of data. There are 28 duplicate associations, thus 1563 of the 2192 high-latitude gamma-ray sources of the 3FGL catalog are AGNs. Most of them (98%) are blazars. About half of the newly detected blazars are of unknown type, i.e., they lack spectroscopic information of sufficient quality to determine the strength of their emission lines. Based on their gamma-ray spectral properties, these sources are evenly split between flat-spectrum radio quasars (FSRQs) and BL Lacs. The most abundant detected BL Lacs are of the high-synchrotron-peaked (HSP) type. About 50% of the BL Lacs have no measured redshifts. A few new rare outliers (HSP-FSRQs and high-luminosity HSP BL Lacs) are reported. The general properties of the 3LAC sample confirm previous findings from earlier catalogs. The fraction of 3LAC blazars in the total population of blazars listed in BZCAT remains non-negligible even at the faint ends of the BZCAT-blazar radio, optical, and X-ray flux distributions, which hints that even the faintest known blazars could eventually shine in gamma-rays at LAT-detection levels. The energy-flux distributions of the different blazar populations are in good agreement with extrapolation from earlier catalogs.

Ultraluminous x-ray sources (ULXs) in nearby galaxies shine brighter than any X-ray source in our Galaxy. ULXs are usually modeled as stellar-mass black holes (BHs) accreting at very high rates or intermediate-mass BHs. We present observations showing that NGC 5907 ULX is instead an x-ray accreting neutron star (NS) with a spin period evolving from 1.43 s in 2003 to 1.13 s in 2014. It has an isotropic peak luminosity of ∼1000 times the Eddington limit for a NS at 17.1 Mpc. Standard accretion models fail to explain its luminosity, even assuming beamed emission, but a strong multipolar magnetic field can describe its properties. These findings suggest that other extreme ULXs (x-ray luminosity ≥1041 erg s −1 ) might harbor NSs.Ultraluminous x-ray sources (ULXs) are observed in off-nucleus regions of nearby galaxiesand have x-ray luminosities in excess of a few 10 39 erg s −1 , which is the Eddington luminosity (L Edd ) for a black hole (BH) of 10 M (1). The L Edd sets an upper limit on the accretion luminosity (L acc ) of a compact object steadily accreting, since for L acc > L Edd accretion will be halted by radiation forces. For spherical accretion of fully ionized hydrogen, the limit can be written as, where σ T is the Thomson scattering cross section, m p is the proton mass, and M/M is the compact object mass in solar masses; for a 1.4 M neutron star (NS), the maximum accreting luminosity is ∼2×10 38 erg s −1 .The high luminosity of ULXs has thus been explained as accretion at or above the Eddington luminosity onto BHs of stellar origin (<80-100 M ), or onto intermediate-mass (10BHs (2, 3). However, if the emission of ULXs were beamed over a fraction b < 1 of the sky, their true luminosity, and thus also the compact object mass required not to exceed L Edd , would be reduced by the same factor. This possibility, together with the recent identification of two accreting NSs associated with the ∼10 40 erg s −1 M82 X-2 (4) and NGC 7793 P13 (5, 6) x-ray sources, have brought support to the view that most low-luminosity ULXs likely host a NS (7) 2 or a stellar-mass BH (8). For the most extreme ULXs with x-ray luminosity exceeding a few ×10 40 erg s −1 , BHs with masses in excess of 100 M are still commonly considered (9, 10).Despite several searches for coherent x-ray pulsations,no other ultraluminous x-ray source has been found to host a NS so far (11).Within the framework of "Exploring the X-ray Transient and variable Sky", EXTraS (12) Fig. 1 and Table 1). In all cases, a strong first period derivative term is present (see Table 1). The pulse shape is nearly sinusoidal, while the pulsed fraction (the semi-amplitude of the sinusoid divided by the average count rate)is energy dependent and increases from about 12% at low energies (<2.5 keV) to ∼20% in the hard band (>7 keV; Fig. 1).To derive constraints on the orbital period (P orb ), we applied a likelihood analysis to the two 2014 NuSTAR observations (see supplementary online text), which have the longest baseline. 3By assuming a circular orbit (as in the case of M...

We apply a number of statistical and machine learning techniques to classify and rank gamma-ray sources from the Third Fermi Large Area Telescope Source Catalog (3FGL), according to their likelihood of falling into the two major classes of gamma-ray emitters: pulsars (PSR) or active galactic nuclei (AGNs). Using 1904 3FGL sources that have been identified/associated with AGNs (1738) and PSR (166), we train (using 70% of our sample) and test (using 30%) our algorithms and find that the best overall accuracy (>96%) is obtained with the Random Forest (RF) technique, while using a logistic regression (LR) algorithm results in only marginally lower accuracy. We apply the same techniques on a subsample of 142 known gamma-ray pulsars to classify them into two major subcategories: young (YNG) and millisecond pulsars (MSP). Once more, the RF algorithm has the best overall accuracy (∼90%), while a boosted LR analysis comes a close second. We apply our two best models (RF and LR) to the entire 3FGL catalog, providing predictions on the likely nature of unassociated sources, including the likely type of pulsar (YNG or MSP). We also use our predictions to shed light on the possible nature of some gamma-ray sources with known associations (e.g., binaries, supernova remnants/pulsar wind nebulae). Finally, we provide a list of plausible X-ray counterparts for some pulsar candidates, obtained using Swift, Chandra, and XMM. The results of our study will be of interest both for in-depth follow-up searches (e.g., pulsar) at various wavelengths and for broader population studies.

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground-and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline, and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

We report on the long term X-ray monitoring of the outburst decay of the low magnetic field magnetar SGR 0418+5729 , using all the available X-ray data obtained with RXTE, Swift, Chandra, and XMM-Newton observations, from the discovery of the source in June 2009, up to August 2012. The timing analysis allowed us to obtain the first measurement of the period derivative of SGR 0418+5729: P = 4(1) × 10 −15 s s −1 , significant at ∼ 3.5σ confidence level. This leads to a surface dipolar magnetic field of B dip 6 × 10 12 Gauss. This measurement confirms SGR 0418+5729 as the lowest magnetic field magnetar. Following the flux and spectral evolution from the beginning of the outburst up to ∼1200 days, we observe a gradual cooling of the tiny hot spot responsible for the X-ray emission, from a temperature of ∼ 0.9 to 0.3 keV. Simultaneously, the X-ray flux decreased by about 3 orders of magnitude: from about 1.4 × 10 −11 to 1.2 × 10 −14 erg s −1 cm −2 . Deep radio, millimeter, optical and gamma-ray observations did not detect the source counterpart, implying stringent limits on its multiband emission, as well as constraints on the presence of a fossil disk. By modeling the magneto-thermal secular evolution of SGR 0418+5729 , we infer a realistic age of ∼550 kyr, and a dipolar magnetic field at birth of ∼ 10 14 G. The outburst characteristics suggest the presence of a thin twisted bundle with a small heated spot at its base. The bundle untwisted in the first few months following the outburst, while the hot spot decreases in temperature and size. We estimate the outburst rate of low magnetic field magnetars to be about one per year per galaxy, and we briefly discuss the consequences of such result in several other astrophysical contexts.