The hygroscopic growth of atmospheric submicrometer particle size distributions (diameter D-p ranging from 22 to 900 nm) was studied at a rural/suburban site in the North China Plain within the framework of the international Campaigns of Air Quality Research in Beijing and Surrounding Region 2006 (CAREBeijing-2006) research project. The goal was to characterize the regional aerosol in the polluted northeastern plain in China. Size descriptive hygroscopic growth factors (DHGFs) were determined as a function of relative humidity (RH) by relating the particle number size distribution at a dry condition ( 100 nm), the DHGF are substantially higher than in the Aitken particle mode (D-p < 100 nm) as a result of different chemical composition. The size-dependent behavior of the DHGF highlights the relevance of particulate sulfate production over the North China Plain, accomplished by secondary formation from the gas phase and, potentially, liquid phase processes in convective clouds. Furthermore, all results concerning the DHGF show a significant dependency on meteorological air masses. The hygroscopic growth of accumulation mode particles correlates significantly with the PM1-mass fraction of sulfate ions determined by chemical analysis. Finally, this investigation provides a parameterization of the hygroscopic growth of 250-nm particles, which might be useful when predicting visibility and radiative forcing and performing atmospheric aerosol model validations.

The type of Polar stratospheric clouds (PSCs) as well as their temporal and spatial extent are important for the occurrence of heterogeneous reactions in the polar stratosphere. The formation of PSCs depends strongly on temperature. However, the mechanisms of the formation of solid PSCs are still poorly understood. Recent satellite studies of Antarctic PSCs have shown that their formation can be associated with deep-tropospheric clouds which have the ability to cool the lower stratosphere radiatively and/or adiabatically. In the present study, lidar measurements aboard the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellite were used to investigate whether the formation of Arctic PSCs can be associated with deep-tropospheric clouds as well. Deep-tropospheric cloud systems have a vertical extent of more than 6.5 km with a cloud top height above 7 km altitude. PSCs observed by CALIPSO during the Arctic winter 2007/2008 were classified according to their type (STS, NAT, or ice) and to the kind of underlying tropospheric clouds. Our analysis reveals that 172 out of 211 observed PSCs occurred in connection with tropospheric clouds. 72% of these 172 observed PSCs occurred above deep-tropospheric clouds. We also find that the type of PSC seems to be connected to the characteristics of the underlying tropospheric cloud system. During the Arctic winter 2007/2008 PSCs consisting of ice were mainly observed in connection with deep-tropospheric cloud systems while no ice PSC was detected above cirrus. Furthermore, we find no correlation between the occurrence of PSCs and the top temperature of tropospheric clouds. Thus, our findings suggest that Arctic PSC formation is connected to adiabatice cooling, i.e. dynamic effects rather than radiative cooling.

The Department of Meteorology at Stockholm University operates the Esrange Rayleigh/Raman lidar at Esrange(68° N, 21° E) near the Swedish city of Kiruna. This paper describes the design and first measurements of the newpure rotational-Raman channel of the Esrange lidar. The Esrange lidar uses a pulsed Nd:YAG solid-state laser operating at 532 nm as light source with a repetition rate of 20 Hz and a pulse energy of 350 mJ. The minimum vertical resolution is 150m and the integration time for one profile is 5000 shots. The newly implemented channel allows for measurements of atmospheric temperature at altitudes below 35 km and is currently optimized for temperature measurements between 180 and 200 K. This corresponds to conditions in the lower Arctic stratosphere during winter. In addition to the temperature measurements, the aerosol extinction coefficientand the aerosol backscatter coefficient at 532 nm can be measured in dependently. Our filter-based design minimizes the systematic error in the obtained temperature profile to less than 0.51 K. By combining rotational-Raman measurements (5–35 km height) and the integration technique (30–80 kmheight), the Esrange lidar is now capable of measuring atmospheric temperature profiles from the upper troposphere up to the mesosphere. With the improved setup, the system can be used to validate current lidar-based polar stratospheric cloud classification schemes. The new capability of the instrument measuring temperature and aerosol extinction furthermore enables studies of the thermal structure and variability of the upper troposphere/lower stratosphere. Although several lidars are operated at polar latitudes, there are few instruments that are capable of measuring temperature profiles in the troposphere, stratosphere, and mesosphere, as well as aerosols extinction in the troposphere and lower stratospherewith daylight capability.

Polar stratospheric clouds (PSCs) play a key role in heterogeneous chemistry and ozone depletion in the lower stratosphere. The type of PSC as well as their temporal and spatial extent are important for the occurrence of heterogeneous reactions and, thus, ozone depletion. In this study a combination of ground-based and spaceborne lidar measurements were used together with microphysical box model simulations along back trajectories to investigate the formation and alteration of Arctic PSCs. The measurements were made by the Rayleigh/Mie/Raman lidar system at Esrange and by the Cloud-Aerosol Lidar with Orthogonal Polarization aboard the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellite. Between 20 and 23 January 2008 PSCs composed of liquid particles were observed by CALIPSO between Greenland and the western side of the Scandinavian Mountains. Between 21 and 23 January 2008 the Esrange lidar observed a PSC composed of distinct layers of liquid and solid particles on the eastern side of the mountain range. Microphysical box model simulations along air parcel back trajectories indicate that liquid particles had formed at least 40 h before the observation at Esrange. Furthermore, the model indicates a high HNO(3) uptake into the liquid layer between 10 and 20 h before the observation. The PSC was formed when the air mass was over Greenland. On two occasions during these 20 h, CALIPSO observed PSCs when its measurement tracks crossed the air parcel back trajectory ending at the location of the Esrange lidar. Backscatter ratios calculated from the output of the box model simulation indicate good agreement with the values observed with the Esrange lidar and by CALIPSO. The box model simulations along the back trajectories from Esrange to the CALIPSO ground track and beyond provide us with the unique opportunity to relate ground-based and spaceborne lidar measurements that were not performed at the same spatial location and time. Furthermore, possible differences in the observations from ground and space can be traced to temporal and/or geographically induced changes in particle microphysics within the measured PSCs.

Lidar measurements of polar stratospheric clouds (PSCs) are commonly analyzed in classification schemes that apply the backscatter ratio and the particle depolarization ratio. This similarity of input data suggests comparable results of different classification schemes—despite measurements being performed with a variety of mostly custom-made instruments. Based on a time series of 16 years of lidar measurements at Esrange (68°N, 21°E), Sweden, we show that PSC classification differs substantially depending on the applied scheme. The discrepancies result from varying threshold values of lidar-derived parameters used to define certain PSC types. The resulting inconsistencies could impact the understanding of long-term PSC observations documented in the literature. We identify two out of seven considered classification schemes that are most likely to give reliable results and should be used in future lidar-based studies. Using polarized backscatter ratios gives the advantage of increased contrast for observations of weakly backscattering and weakly depolarizing particles. Improved confidence in PSC classification can be achieved by a more comprehensive consideration of the effect of measurement uncertainties. The particle depolarization ratio is the key to a reliable identification of different PSC types. Hence, detailed information on the calibration of the polarization-sensitive measurement channels should be provided to assess the findings of a study. Presently, most PSC measurements with lidar are performed at 532 nm only. The information from additional polarization-sensitive measurements in the near infrared could lead to an improved PSC classification. Coincident lidar-based temperature measurements at PSC level might provide useful information for an assessment of PSC classification.

Context. Fomalhaut is a young (2 +/- 1 x 10(8) years), nearby (7.7 pc), 2 M-circle dot star that is suspected to harbor an infant planetary system, interspersed with one or more belts of dusty debris. Aims. We present far-infrared images obtained with the Herschel Space Observatory with an angular resolution between 5.7 '' and 36.7 '' at wavelengths between 70 mu m and 500 mu m. The images show the main debris belt in great detail. Even at high spatial resolution, the belt appears smooth. The region in between the belt and the central star is not devoid of material; thermal emission is observed here as well. Also at the location of the star, excess emission is detected. We aim to construct a consistent image of the Fomalhaut system. Methods. We use a dynamical model together with radiative-transfer tools to derive the parameters of the debris disk. We include detailed models of the interaction of the dust grains with radiation, for both the radiation pressure and the temperature determination. Comparing these models to the spatially resolved temperature information contained in the images allows us to place strong constraints on the presence of grains that will be blown out of the system by radiation pressure. We use this to derive the dynamical parameters of the system. Results. The appearance of the belt points toward a remarkably active system in which dust grains are produced at a very high rate by a collisional cascade in a narrow region filled with dynamically excited planetesimals. Dust particles with sizes below the blow-out size are abundantly present. The equivalent of 2000 one-km-sized comets are destroyed every day, out of a cometary reservoir amounting to 110 Earth masses. From comparison of their scattering and thermal properties, we find evidence that the dust grains are fluffy aggregates, which indicates a cometary origin. The excess emission at the location of the star may be produced by hot dust with a range of temperatures, but may also be due to gaseous free-free emission from a stellar wind.

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

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.

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 four years of Fermi-LAT data under the assumption of the universal cosmic-ray (CR) model proposed by Pinzke & Pfrommer. We find an excess at a significance of 2.7 sigma, which upon closer inspection, however, is correlated to individual excess emission toward three galaxy clusters: A400, A1367, and A3112. We discuss these cases in detail and conservatively attribute the emission to unmodeled background systems (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 CR to thermal pressure ratio within the virial radius, R-200, to be below 1.25%-1.4% depending on the morphological classification. In addition, we derive new limits on the gamma-ray flux from individual clusters in our sample.

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 the use both of reprocessed data and of additional

We present the detections of 18 solar flares detected in high-energy gamma-rays (above 100 MeV) with the Fermi Large Area Telescope (LAT) during its first 4 yr of operation. This work suggests that particle acceleration up to very high energies in solar flares is more common than previously thought, occurring even in modest flares, and for longer durations. Interestingly, all these flares are associated with fairly fast coronal mass ejections (CMEs). 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 hr, and the 2011 June 7 M2.5 flare, which was followed by gamma-ray emission lasting for 2 hr. We compare the Fermi LAT data with X-ray and proton data measurements from GOES and RHESSI. We argue that the gamma-rays are more likely produced through pion decay than electron bremsstrahlung, and we find that the energy spectrum of the proton distribution softens during the extended emission of the 2011 March 7 flare. This would disfavor a trapping scenario for particles accelerated during the impulsive phase of the flare and point to a continuous acceleration process at play for the duration of the flares. CME shocks are known for accelerating the solar energetic particles (SEPs) observed in situ on similar timescales, but it might be challenging to explain the production of gamma-rays at the surface of the Sun while the CME is halfway to the Earth. A stochastic turbulence acceleration process occurring in the solar corona is another likely scenario. Detailed comparison of characteristics of SEPs and gamma-ray-emitting particles for several flares will be helpful to distinguish between these two possibilities.

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 Fermi Large Area Telescope observations of the gamma-ray emission from Earth's limb for an indirect measurement of the local spectrum of CR protons in the energy range similar to 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 +/- 0.04 and 2.61 +/- 0.08 above similar to 200 GeV, respectively.

The Fermi Large Area Telescope (LAT) has an instantaneous field of view (FoV) covering similar to 1/5 of the sky and it completes a survey of the entire sky in high-energy gamma-rays every 3 hr. It enables searches for transient phenomena over timescales from milliseconds to years. Among these phenomena could be electromagnetic counterparts to gravitational wave (GW) sources. In this paper, we present a detailed study of the LAT observations relevant to Laser Interferometer Gravitational-wave Observatory (LIGO) event GW150914, which is the first direct detection of gravitational waves and has been interpreted as being due to the coalescence of two stellar-mass black holes. The localization region for GW150914 was outside the LAT FoV at the time of the GW signal. However, as part of routine survey observations, the LAT observed the entire LIGO localization region within similar to 70 minutes of the trigger and thus enabled a comprehensive search for a.-ray counterpart to GW150914. The study of the LAT data presented here did not find any potential counterparts to GW150914, but it did provide limits on the presence of a transient counterpart above 100 MeV on timescales of hours to days over the entire GW150914 localization region.

Dark matter in the Milky Way may annihilate directly into. rays, producing a monoenergetic spectral line. Therefore, detecting such a signature would be strong evidence for dark matter annihilation or decay. We search for spectral lines in the Fermi Large Area Telescope observations of the Milky Way halo in the energy range 200 MeV-500 GeV using analysis methods from our most recent line searches. The main improvements relative to previous works are our use of 5.8 years of data reprocessed with the Pass 8 event-level analysis and the additional data resulting from the modified observing strategy designed to increase exposure of the Galactic center region. We search in five sky regions selected to optimize sensitivity to different theoretically motivated dark matter scenarios and find no significant detections. In addition to presenting the results from our search for lines, we also investigate the previously reported tentative detection of a line at 133 GeV using the new Pass 8 data.

We report for the first time a gamma-ray and multiwavelength nearly periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope we have discovered an apparent quasi-periodicity in the gamma-ray flux (E > 100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/- 0.08 year period gamma-ray cycle is strengthened by correlated oscillations observed in radio and optical fluxes, through data collected in the Owens Valley Radio Observatory, Tuorla, Katzman Automatic Imaging Telescope, and Catalina Sky Survey monitoring programs and Swift-UVOT. The optical cycle appearing in similar to 10 years of data has a similar period, while the 15 GHz oscillation is less regular than seen in the other bands. Further long-term multiwavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity.

Satellite galaxies of the Milky Way are among the most promising targets for dark matter searches in gamma rays. We present a search for dark matter consisting of weakly interacting massive particles, applying a joint likelihood analysis to 10 satellite galaxies with 24 months of data of the Fermi Large Area Telescope. No dark matter signal is detected. Including the uncertainty in the dark matter distribution, robust upper limits are placed on dark matter annihilation cross sections. The 95% confidence level upper limits range from about 10(-26) cm(3) s(-1) at 5 GeV to about 5 x 10(-23) cm(3) s(-1) at 1 TeV, depending on the dark matter annihilation final state. For the first time, using gamma rays, we are able to rule out models with the most generic cross section (similar to 3 x 10(-26) cm(3) s(-1) for a purely s-wave cross section), without assuming additional boost factors.

We present results from gamma-ray observations of the Coma cluster incorporating six years of Fermi-LAT data and the newly released "Pass 8" event-level analysis. Our analysis of the region reveals low-significance residual structures within the virial radius of the cluster that are too faint for a detailed investigation with the current data. Using a likelihood approach that is free of assumptions on the spectral shape we derive upper limits on the gamma-ray flux that is expected from energetic particle interactions in the cluster. We also consider a benchmark spatial and spectral template motivated by models in which the observed radio halo is mostly emission by secondary electrons. In this case, the median expected and observed upper limits for the flux above 100MeV are 1.7 x 10(-9) ph cm(-2) s(-1) and 5.2 x 10(-9) ph cm(-2) s(-1) respectively (the latter corresponds to residual emission at the level of 1.8 sigma). These bounds are comparable to or higher than predicted levels of hadronic gamma-ray emission in cosmic-ray (CR) models with or without reacceleration of secondary electrons, although direct comparisons are sensitive to assumptions regarding the origin and propagation mode of CRs and magnetic field properties. The minimal expected.-ray flux from radio and star-forming galaxies within the Coma cluster is roughly an order of magnitude below the median sensitivity of our analysis.

The Fermi Large Area Telescope (LAT) Collaboration has recently released a catalog of 360 sources detected above 50 GeV (2FHL). This catalog was obtained using 80 months of data re-processed with Pass 8, the newest event-level analysis, which significantly improves the acceptance and angular resolution of the instrument. Most of the 2FHL sources at high Galactic latitude are blazars. Using detailed Monte Carlo simulations, we measure, for the first time, the source count distribution, dN= dS, of extragalactic.-ray sources at E > 50 GeV and find that it is compatible with a Euclidean distribution down to the lowest measured source flux in the 2FHL (8 x 10(-12) ph cm(-2) s(-1)). We employ a one-point photon fluctuation analysis to constrain the behavior of dN= dS below the source detection threshold. Overall, the source count distribution is constrained over three decades in flux and found compatible with a broken power law with a break flux, Sb, in the range [8 x 10-12; 1.5 x 10-11] ph cm(-2) s(-1) and power-law indices below and above the break of a 2. [1.60; 1.75] and a 1 +/- 2.49 +/- 0.12, respectively. Integration of dN= dS shows that point sources account for at least 86_16 -14 % of the total extragalactic gamma-ray background. The simple form of the derived source count distribution is consistent with a single population (i. e., blazars) dominating the source counts to the minimum flux explored by this analysis. We estimate the density of sources detectable in blind surveys that will be performed in the coming years by the Cherenkov Telescope Array.

The gamma-ray sky can be decomposed into individually detected sources, diffuse emission attributed to the interactions of Galactic cosmic rays with gas and radiation fields, and a residual all-sky emission component commonly called the isotropic diffuse gamma-ray background (IGRB). The IGRB comprises all extragalactic emissions too faint or too diffuse to be resolved in a given survey, as well as any residual Galactic foregrounds that are approximately isotropic. The first IGRB measurement with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) used 10 months of sky-survey data and considered an energy range between 200 MeV and 100 GeV. Improvements in event selection and characterization of cosmic-ray backgrounds, better understanding of the diffuse Galactic emission (DGE), and a longer data accumulation of 50 months allow for a refinement and extension of the IGRB measurement with the LAT, now covering the energy range from 100 MeV to 820 GeV. The IGRB spectrum shows a significant high-energy cutoff feature and can be well described over nearly four decades in energy by a power law with exponential cutoff having a spectral index of 2.32 +/- 0.02 and a break energy of (279 +/- 52) GeV using our baseline DGE model. The total intensity attributed to the IGRB is (7.2 +/- 0.6) x 10(-6) cm(-2) s(-1) sr(-1) above 100 MeV, with an additional +15%/-30% systematic uncertainty due to the Galactic diffuse foregrounds.

The region around the Galactic Center (GC) is now well established to be brighter at energies of a few GeV than what is expected from conventional models of diffuse gamma-ray emission and catalogs of known gamma-ray sources. We study the GeV excess using 6.5 yr of data from the Fermi Large Area Telescope. We characterize the uncertainty of the GC excess spectrum and morphology due to uncertainties in cosmic-ray source distributions and propagation, uncertainties in the distribution of interstellar gas in the Milky Way, and uncertainties due to a potential contribution from the Fermi bubbles. We also evaluate uncertainties in the excess properties due to resolved point sources of gamma rays. The GC is of particular interest, as it would be expected to have the brightest signal from annihilation of weakly interacting massive dark matter (DM) particles. However, control regions along the Galactic plane, where a DM signal is not expected, show excesses of similar amplitude relative to the local background. Based on the magnitude of the systematic uncertainties, we conservatively report upper limits for the annihilation cross-section as a function of particle mass and annihilation channel.

Galaxy clusters are one of the prime sites to search for dark matter (DM) annihilation signals. Depending on the substructure of the DM halo of a galaxy cluster and the cross sections for DM annihilation channels, these signals might be detectable by the latest generation of gamma-ray telescopes. Here we use three years of Fermi-Large Area Telescope data, which are the most suitable for searching for very extended emission in the vicinity of the nearby Virgo galaxy cluster. Our analysis reveals statistically significant extended emission which can be well characterized by a uniformly emitting disk profile with a radius of 3 degrees that moreover is offset from the cluster center. We demonstrate that the significance of this extended emission strongly depends on the adopted interstellar emission model (IEM) and is most likely an artifact of our incomplete description of the IEM in this region. We also search for and find new point source candidates in the region. We then derive conservative upper limits on the velocity-averaged DM pair annihilation cross section from Virgo. We take into account the potential gamma-ray flux enhancement due to DM sub-halos and its complex morphology as a merging cluster. For DM annihilating into b (b) over bar, assuming a conservative sub-halo model setup, we find limits that are between 1 and 1.5 orders of magnitude above the expectation from the thermal cross section for m(DM) <= 100 GeV. In a more optimistic scenario, we exclude <sigma v > similar to 3 x 10(-26)cm(3)s(-1) for m(DM)less than or similar to 40 GeV for the same channel. Finally, we derive upper limits on the gamma-ray-flux produced by hadronic cosmic-ray interactions in the inter cluster medium. We find that the volume-averaged cosmic-ray-to-thermal pressure ratio is less than similar to 6%.

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 Large Area Telescope at Galactic latitudes vertical bar b vertical bar > 30 degrees in four energy bins spanning 1-50 GeV. At multipoles l >= 155, corresponding to angular scales less than or similar to 2 degrees, angular power above the photon noise level is detected at >99.99% confidence level in the 1-2 GeV, 2-5 GeV, and 5-10 GeV energy bins, and at >99% confidence level at 10-50 GeV. Within each energy bin the measured angular power takes approximately the same value at all multipoles l >= 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.

A classical nova results from runaway thermonuclear explosions on the surface of a white dwarf that accretes matter from a low-mass main-sequence stellar companion. In 2012 and 2013, three novae were detected in gamma rays and stood in contrast to the first gamma-ray-detected nova V407 Cygni 2010, which belongs to a rare class of symbiotic binary systems. Despite likely differences in the compositions and masses of their white dwarf progenitors, the three classical novae are similarly characterized as soft-spectrum transient gamma-ray sources detected over 2- to 3-week durations. The gamma-ray detections point to unexpected high-energy particle acceleration processes linked to the mass ejection from thermonuclear explosions in an unanticipated class of Galactic gamma-ray sources.

The Fermi Large Area Telescope (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 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 upper limit of 2.0 x 10(-12) erg cm(-2) s(-1) on the 0.1-100 GeV energy flux (95% confidence level). This revised upper limit 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. Its spectrum is consistent with a power law with photon index Gamma = 2.4 +/- 0.1(stat) (vertical bar) (syst) and a 0.1-100 GeV energy flux of (5.6 +/- 0.6(stat vertical bar syst)) x 10(-12) erg cm(-1) s(-1). M31 is detected to be extended with a 4 sigma significance. The spatial distribution of the emission is consistent with a uniform-brightness disk with a radius of 0 degrees.4 and no offset from the center of the galaxy, 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 of the galaxy or extend far from it. The gamma-ray signal is not correlated with regions rich in gas or star formation activity, which suggests 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.

We search for evidence of dark matter (DM) annihilation in the isotropic gamma-ray background (IGRB) measured with 50 months of Fermi Large Area Telescope (LAT) observations. An improved theoretical description of the cosmological DM annihilation signal, based on two complementary techniques and assuming generic weakly interacting massive particle (WIMP) properties, renders more precise predictions compared to previous work. More specifically, we estimate the cosmologically-induced gamma-ray intensity to have an uncertainty of a factor similar to 20 in canonical setups. We consistently include both the Galactic and extragalactic signals under the same theoretical framework, and study the impact of the former on the IGRB spectrum derivation. We find no evidence for a DM signal and we set limits on the DM-induced isotropic gamma-ray signal. Our limits are competitive for DM particle masses up to tens of TeV and, indeed, are the strongest limits derived from Fermi LAT data at TeV energies. This is possible thanks to the new Fermi LAT IGRB measurement, which now extends up to an energy of 820 GeV. We quantify uncertainties in detail and show the potential this type of search offers for testing the WIMP paradigm with a complementary and truly cosmological probe of DM particle signals.

Dark matter particle annihilation or decay can produce monochromatic gamma-ray lines and contribute to the diffuse gamma-ray background. Flux upper limits are presented for gamma-ray spectral lines from 7 to 200 GeV and for the diffuse gamma-ray background from 4.8 GeV to 264 GeV obtained from two years of Fermi Large Area Telescope data integrated over most of the sky. We give cross-section upper limits and decay lifetime lower limits for dark matter models that produce gamma-ray lines or contribute to the diffuse spectrum, including models proposed as explanations of the PAMELA and Fermi cosmic-ray data.

We present a detailed statistical analysis of the correlation between radio and gamma-ray emission of the active galactic nuclei (AGNs) detected by Fermi during its first year of operation, with the largest data sets ever used for this purpose. We use both archival interferometric 8.4 GHz data (from the Very Large Array 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 centimeter radio and the broadband (E > 100 MeV) gamma-ray energy flux is very high for the whole AGN sample, with a probability of < 10(-7) for the correlation appearing by chance. Using the OVRO data, we find that concurrent data improve the significance of the correlation from 1.6 x 10(-6) to 9.0 x 10(-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 < 10(-7)) for both flat spectrum radio quasars and BL Lac objects 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.

The second catalog of active galactic nuclei (AGNs) detected by the Fermi Large Area Telescope (LAT) in two years of scientific operation is presented. The second LAT AGN catalog (2LAC) includes 1017 γ-ray sources located at high Galactic latitudes (|b| > 10°) that are detected with a test statistic (TS) greater than 25 and associated statistically with AGNs. However, some of these are affected by analysis issues and some are associated with multiple AGNs. Consequently, we define a Clean Sample which includes 886 AGNs, comprising 395 BL Lacertae objects (BL Lac objects), 310 flat-spectrum radio quasars (FSRQs), 157 candidate blazars of unknown type (i.e., with broadband blazar characteristics but with no optical spectral measurement yet), 8 misaligned AGNs, 4 narrow-line Seyfert 1 (NLS1s), 10 AGNs of other types, and 2 starburst galaxies. Where possible, the blazars have been further classified based on their spectral energy distributions (SEDs) as archival radio, optical, and X-ray data permit. While almost all FSRQs have a synchrotron-peak frequency <1014 Hz, about half of the BL Lac objects have a synchrotron-peak frequency >1015 Hz. The 2LAC represents a significant improvement relative to the first LAT AGN catalog (1LAC), with 52% more associated sources. The full characterization of the newly detected sources will require more broadband data. Various properties, such as γ-ray fluxes and photon power-law spectral indices, redshifts, γ-ray luminosities, variability, and archival radio luminosities and their correlations are presented and discussed for the different blazar classes. The general trends observed in 1LAC are confirmed.

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 similar to 100 MeV and similar to 100 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 similar to 10% 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 (W-CO) at a 1 degrees x 1 degrees pixel level. The correlation is found to be linear over a W-CO range of similar to 10-fold when divided in three regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The W-CO-to-mass conversion factor, X-CO, is found to be similar to 2.3 x 10(20) cm(-2) (K km s(-1))(-1) for the high-longitude part of Orion A (l > 212 degrees), similar to 1.7 times higher than similar to 1.3 x 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 nonlinear relation between H-2 and CO densities in the diffuse molecular gas. W-CO decreases faster than the H-2 column density in the region making the gas darker to W-CO.

The Fermi Large Area Telescope (LAT) First Source Catalog (1FGL) provided spatial, spectral, and temporal properties for a large number of gamma-ray sources using a uniform analysis method. After correlating with the most-complete catalogs of source types known to emit gamma rays, 630 of these sources are unassociated (i.e., have no obvious counterparts at other wavelengths). Here, we employ two statistical analyses of the primary gamma-ray characteristics for these unassociated sources in an effort to correlate their gamma-ray properties with the active galactic nucleus (AGN) and pulsar populations in 1FGL. Based on the correlation results, we classify 221 AGN-like and 134 pulsar-like sources in the 1FGL unassociated sources. The results of these source classifications appear to match the expected source distributions, especially at high Galactic latitudes. While useful for planning future multiwavelength follow-up observations, these analyses use limited inputs, and their predictions should not be considered equivalent to probable source classes for these sources. We discuss multiwavelength results and catalog cross-correlations to date, and provide new source associations for 229 Fermi-LAT sources that had no association listed in the 1FGL catalog. By validating the source classifications against these new associations, we find that the new association matches the predicted source class in similar to 80% of the sources.

The Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope is a pair-conversion telescope designed to detect photons with energies from approximate to 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 (AGNs). 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 AGNs. 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 AGNs and the TeV blazars 1ES0229 + 200 and 1ES0347-121.

The Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope is a pair-conversion telescope designed to survey the gamma-ray sky from 20 MeV to several hundreds of GeV. In this energy band there are no astronomical sources with sufficiently well known and sharp spectral features to allow an absolute calibration of the LAT energy scale. However, the geomagnetic cutoff in the cosmic ray electron-plus-positron (CRE) spectrum in low Earth orbit does provide such a spectral feature. The energy and spectral shape of this cutoff can be calculated with the aid of a numerical code tracing charged particles in the Earth's magnetic field. By comparing the cutoff value with that measured by the LAT in different geomagnetic positions, we have obtained several calibration points between similar to 6 and similar to 13 GeV with an estimated uncertainty of similar to 2%. An energy calibration with such high accuracy reduces the systematic uncertainty in LAT measurements of, for example, the spectral cutoff in the emission from gamma ray pulsars.

We present a catalog of gamma-ray sources at energies above 10 GeV based on data from the Large Area Telescope (LAT) accumulated during the first 3 yr of the Fermi Gamma-ray Space Telescope mission. The first Fermi catalog of > 10 GeV sources (1FHL) has 514 sources. For each source we present location, spectrum, a measure of variability, and associations with cataloged sources at other wavelengths. We found that 449 (87%) could be associated with known sources, of which 393 (76% of the 1FHL sources) are active galactic nuclei. Of the 27 sources associated with known pulsars, we find 20 (12) to have significant pulsations in the range > 10 GeV (> 25 GeV). In this work we also report that, at energies above 10 GeV, unresolved sources account for 27% +/- 8% of the isotropic. gamma-ray background, while the unresolved Galactic population contributes only at the few percent level to the Galactic diffuse background. We also highlight the subset of the 1FHL sources that are best candidates for detection at energies above 50-100 GeV with current and future ground-based gamma-ray observatories.

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 Earth's shadow, which is offset in opposite directions for opposite charges due to 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 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 GeVare consistent with a fraction that is continuing to rise with energy.

The Geostationary Operational Environmental Satellite (GOES) M2-class solar flare, SOL2010-06-12T00: 57, was modest in many respects yet exhibited remarkable acceleration of energetic particles. The flare produced an similar to 50 s impulsive burst of hard X-and gamma-ray emission up to at least 400 MeV observed by the Fermi Gamma-ray Burst Monitor and Large Area Telescope experiments. The remarkably similar hard X-ray and high-energy gamma-ray time profiles suggest that most of the particles were accelerated to energies greater than or similar to 300 MeV with a delay of similar to 10 s from mildly relativistic electrons, but some reached these energies in as little as similar to 3 s. The gamma-ray line fluence from this flare was about 10 times higher than that typically observed from this modest GOES class of X-ray flare. There is no evidence for time-extended >100 MeV emission as has been found for other flares with high-energy gamma-rays.

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 report an analysis of the interstellar gamma-ray emission from the Chamaeleon, R Coronae Australis (R CrA), and Cepheus and Polaris flare regions with the Fermi Large Area Telescope. They are among the nearest molecular cloud complexes, within similar to 300 pc from the solar system. The gamma-ray emission produced by interactions of cosmic rays (CRs) and interstellar gas in those molecular clouds is useful to study the CR densities and distributions of molecular gas close to the solar system. The obtained gamma-ray emissivities above 250 MeV are (5.9 +/- 0.1(stat-1.0sys)(+0.9)) x 10(-27) photons s(-1) sr(-1) H-atom(-1), (10.2 +/- 0.4(stat-1.7sys)(+1.2)) x 10(-27) photons s(-1) sr(-1) H-atom(-1), and (9.1 +/- 0.3(stat-0.6sys)(+1.5)) x 10(-27) photons s(-1) sr(-1) H-atom(-1) for the Chamaeleon, R CrA, and Cepheus and Polaris flare regions, respectively. Whereas the energy dependences of the emissivities agree well with that predicted from direct CR observations at the Earth, the measured emissivities from 250 MeV to 10 GeV indicate a variation of the CR density by similar to 20% in the neighborhood of the solar system, even if we consider systematic uncertainties. The molecular mass calibrating ratio, X-CO = N(H-2)/W-CO, is found to be (0.96 +/- 0.06(stat-0.12sys)(+0.15)) x 10(20) H-2-molecule cm(-2) (K km s(-1))(-1), (0.99 +/- 0.08(stat-0.10sys)(+0.18)) x 10(20) H-2-molecule cm(-2) (K km s(-1))(-1), and (0.63 +/- 0.02(stat-0.07sys)(+0.09)) x 10(20) H-2-molecule cm(-2) (K km s(-1))(-1) for the Chamaeleon, R CrA, and Cepheus and Polaris flare regions, respectively, suggesting a variation of X-CO in the vicinity of the solar system. From the obtained values of X-CO, the masses of molecular gas traced by W-CO in the Chamaeleon, R CrA, and Cepheus and Polaris flare regions are estimated to be similar to 5 x 10(3)M(circle dot), similar to 10(3)M(circle dot), and similar to 3.3 x 10(4)M(circle dot), respectively. A comparable amount of gas not traced well by standard Hi and CO surveys is found in the regions investigated.