Chang'E 4 is the first mission to the far side of the Moon and consists of a lander, a rover, and a relay spacecraft. Lander and rover were launched at 18:23 UTC on December 7, 2018 and landed in the von K\'arm\'an crater at 02:26 UTC on January 3, 2019. Here we describe the Lunar Lander Neutron \& Dosimetry experiment (LND) which is part of the Chang'E 4 Lander scientific payload. Its chief scientific goal is to obtain first active dosimetric measurements on the surface of the Moon. LND also provides observations of fast neutrons which are a result of the interaction of high-energy particle radiation with the lunar regolith and of their thermalized counterpart, thermal neutrons, which are a sensitive indicator of subsurface water content.

This paper studies the secondary's rotation in a synchronous binary asteroid system in which the secondary enters the 1:1 spin-orbit resonance. The model used is the planar full two-body problem composed of a spherical primary plus a tri-axial ellipsoid secondary. Compared with classical spin-orbit work, there are two differences: (1) Influence on the mutual orbit from the secondary's rotation is considered; (2) Instead of the Hamiltonian approach, the approach of periodic orbits is adopted. Our studies find: (1) Genealogy of the two families of periodic orbits is same as that of the families around triangular libration points in the restricted three-body problem. That is, the long-period family terminates onto a short-period orbit travelling N times; (2) In the limiting case where the secondary's mass is negligible, our results can be reduced to the classical spin-orbit theory, by equating the long-period orbit with the free libration, and by equating the short-period orbit with the forced libration caused by orbit eccentricity. However, the two models show obvious differences when the secondary's mass is non-negligible. (3) By studying the stability of periodic orbits, for a specific binary asteroid system, we are able to obtain the maximum libration amplitude of the secondary (which is usually less than 90 degree), and the maximum mutual orbit eccentricity which does not break the secondary's synchronous state. We also find the anti-correlation between the secondary's libration amplitude and the orbit eccentricity. The (65803) Didymos system is taken as an example to show the results.

The Andromeda (M31) and Triangulum (M33) galaxies are the closest Local Group galaxies to the Milky Way, being only 785 and 870 kpc away. These two galaxies provide an independent view of high-energy processes that are often obscured in our own Galaxy, including possible signals of dark matter (DM) particle interactions. The Fermi Large Area Telescope (Fermi-LAT) preliminary eight year list of sources includes both M31, which is detected as extended with a size of about 0.4$^\circ$, and M33, which is detected as a point-like source. The spatial morphology of M31 $\gamma$-ray emission could trace a population of unresolved sources and energetic particles originating in sources not related to massive star formation. Alternatively, the $\gamma$-ray emission could also be an indication of annihilation or decay of DM particles. We investigate these two possibilities using almost 10 years of data from the Fermi LAT. An interpretation that involves only a DM $\gamma$-ray emission is in tension with the current limits from other searches, such as those targeting Milky Way dwarf spheroidal galaxies. When we include a template of astrophysical emission, tuned on $\gamma$-ray data or from observations of these galaxies in other wavelengths, we do not find any significant evidence for a DM contribution and we set limits for the annihilation cross section that probe the thermal cross section for DM masses up to a few tens of GeV in the $b\bar{b}$ and $\tau^+\tau^-$ channels. For models where the DM substructures have masses above $10^{-6}$ solar masses our limits probe the DM interpretation of the Fermi LAT Galactic center excess. We provide also the lower limit for the DM decay time assuming the same spatial models of the DM distribution in M31 and M33.

We naturally obtain a power-law total mass-density profile for the inner regions (within several effective radius) of early-type galaxies(ETGs) from space-time structures, which are described by the static spherically-symmetric solutions of Einstein equations with the perfect fluid matters. We point out that, in the inner region of an early type galaxy, the average value of the equation of state $\omega$ determines the power-law index $\gamma$ of the total mass-density profile given by \cite{2006ApJ...649..599K}. We then use 119 galaxy-galaxy strong gravitational lensing systems from SLACS, BELLS, BELLS GALLERY and SL2S, to constrain the equation of state $\omega$ for the inner regions of ETGs. We find that the average value of $\omega$ for totally 119 ETGs is $\omega=-0.315\pm0.073$ (68\% CL), which is very close to $-1/3$. If $\omega=-1/3$ is the best value to describe the equation of state of the inner regions of the dark matter in galaxies, the deviation from $-1/3$ could represent the mixture of baryonic matter and dark matter in galaxies.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

We present a multi-wavelength study of the unassociated Fermi-LAT source, 3FGL J0954.8-3948, which is likely the gamma-ray counterpart of a 9.3-hour binary in the field. With more than 9 years of Pass 8 LAT data, we updated the gamma-ray spectral properties and the LAT localization of the gamma-ray source. While the binary lies outside the cataloged 95% error ellipse, the optimized LAT ellipse is 0.1 degrees closer and encloses the binary. The system is likely spectrally hard in X-rays (photon index ~ 1.4) with orbital modulations detected in optical, UV, and possibly X-rays. A steep spectrum radio counterpart (spectral index ~ -1.6) is also found in the TIFR GMRT Sky Survey (TGSS), implying that it is a pulsar system. We obtained a series of SOAR and Gemini spectroscopic observations in 2017/2018, which show a low-mass secondary orbiting in a close circular orbit with K2 = 272 km/s under strong irradiation by the primary compact object. All the observations as well as the modelling of the X/gamma-ray high-energy emission suggest that 3FGL J0954.8-3948 is a redback millisecond pulsar in a rotation-powered state.

A new measurement of a spatially extended gamma-ray signal from the center of the Andromeda galaxy (M31) has been recently published by the Fermi-LAT collaboration, reporting that the emission broadly resembles the so-called Galactic center excess (GCE) of the Milky Way (MW). At the same time, evidence is accumulating on a millisecond pulsar (MSPs) origin for the GCE. These elements prompt us to compare the mentioned observations with what is, perhaps, the simplest model for an MSP population, solely obtained by rescaling of the MSP luminosity function determined in the local MW disk via the respective stellar mass of the systems. It is remarkable that without free fitting parameters, this model can account for both the energetics and the morphology of the GCE within uncertainties. For M31, the estimated luminosity due to primordial MSPs is expected to contribute only about a quarter of the detected emission, although a dominant contribution cannot be excluded given the large uncertainties. If correct, the model predicts that the M31 disk emission due to MSP is not far below the present upper bound. We also discuss a few refinements of this simple model. In particular, we use the correlation between globular cluster gamma-ray luminosity and stellar encounter rate to gauge the dynamical MSP formation in the bulge. This component is expected to contribute to the GCE only at a level $\lesssim 5\%$, but it may be of some importance in explaining the signal's morphology in the inner region of the Galaxy. We also comment on some effects which may lead to violations of the simple scaling used, on alternative models, and on future perspectives for improved diagnostics.

For the first time, we obtain the analytical form of black hole space-time metric in dark matter halo for the stationary situation. Using the relation between the rotation velocity (in the equatorial plane) and the spherical symmetric space-time metric coefficient, we obtain the space-time metric for pure dark matter. By considering the dark matter halo in spherical symmetric space-time as part of the energy-momentum tensors in the Einstein field equation, we then obtain the spherical symmetric black hole solutions in dark matter halo. Utilizing Newman-Jains method, we further generalize spherical symmetric black holes to rotational black holes. As examples, we obtain the space-time metric of black holes surrounded by Cold Dark Matter and Scalar Field Dark Matter halos, respectively. Our main results regarding the interaction between black hole and dark matter halo are as follows: (i) For both dark matter models, the density profile always produces "cusp" phenomenon in small scale in the relativity situation; (ii) Dark matter halo makes the black hole horizon to increase but the ergosphere to decrease, while the magnitude is small; (iii) Dark matter does not change the singularity of black holes. These results are useful to study the interaction of black hole and dark matter halo in stationary situation. Particularly, the "cusp" produced in the $0\sim 1$ kpc scale would be observable in the Milky Way. Perspectives on future work regarding the applications of our results in astrophysics are also briefly discussed.

The spatial extension of a gamma-ray source is an essential ingredient to determine its spectral properties as well as its potential multi-wavelength counterpart. The capability to spatially resolve gamma-ray sources is greatly improved by the newly delivered Fermi-Large Area Telescope (LAT) Pass 8 event-level analysis which provides a greater acceptance and an improved point spread function, two crucial factors for the detection of extended sources. Here, we present a complete search for extended sources located within 7 degrees from the Galactic plane, using 6 years of LAT data above 10 GeV. We find 46 extended sources and provide their morphological and spectral characteristics. This constitutes the first catalog of hard LAT extended sources, named the Fermi Galactic Extended Source Catalog, which allows a thorough study of the properties of the Galactic plane in the sub-TeV domain.

HESS J0632+057 is the only gamma-ray binary that has been detected at TeV energies, but not at GeV energies yet. Based on nearly nine years of Fermi Large Area Telescope (LAT) Pass 8 data, we report here on a deep search for the gamma-ray emission from HESS J0632+057 in the 0.1-300 GeV energy range. We find a previously unknown gamma-ray source, Fermi J0632.6+0548, spatially coincident with HESS J0632+057. The measured flux of Fermi J0632.6+0548 is consistent with the previous flux upper limit on HESS J0632+057 and shows variability that can be related to the HESS J0632+057 orbital phase. We propose that Fermi J0632.6+0548 is the GeV counterpart of HESS J0632+057. Considering the Very High Energy (VHE) spectrum of HESS J0632+057, a possible spectral turnover above 10 GeV may exist in Fermi J0632.6+0548, as appears to be common in other established gamma-ray binaries.

We report our recent Swift, NuSTAR, and XMM-Newton X-ray and Lijiang optical observations on PSR J2032+4127/MT91 213, the gamma-ray binary candidate with a period of 45-50 years. The coming periastron of the system was predicted to be in November 2017, around which high-energy flares from keV to TeV are expected. Recent studies with Chandra and Swift X-ray observations taken in 2015/16 showed that its X-ray emission has been brighter by a factors of ~10 than that before 2013, probably revealing some on-going activities between the pulsar wind and the stellar wind. Our new Swift/XRT lightcurve shows no strong evidence of a single vigorous brightening trend, but rather several strong X-ray flares on weekly to monthly timescales with a slowly brightening baseline, namely the low state. The NuSTAR and XMM-Newton observations taken during the flaring and the low states, respectively, show a denser environment and a softer power-law index during the flaring state, implying that the pulsar wind interacted with stronger stellar winds of the companion to produce the flares. These precursors would be crucial in studying the predicted giant outburst from this extreme gamma-ray binary during the periastron passage in late 2017.

We present a multi-wavelength study of the unidentified Fermi object, 3FGL J0212.1+5320. Within the 95% error ellipse, Chandra detects a bright X-ray source (i.e., F(0.5-7keV) = 1.4e-12 erg/cm^2/s), which has a low-mass optical counterpart (M < 0.4 Msun and T ~ 6000 K). A clear ellipsoidal modulation is shown in optical/infrared at 20.87 hours. The gamma-ray properties of 3FGL J0212.1+5320 are all consistent with that of a millisecond pulsar, suggesting that it is a gamma-ray redback millisecond pulsar binary with a low-mass companion filling > 64% of the Roche-lobe. If confirmed, it will be a redback binary with one of the longest orbital periods known. Spectroscopic data taken in 2015 from the Lijiang observatory show no evidence of strong emission lines, revealing that the accretion is currently inactive (the rotation-powered pulsar state). This is consistent with the low X-ray luminosities (Lx ~ 10^32 erg/s) and the possible X-ray modulation seen by Chandra and Swift. Considering that the X-ray luminosity and the high X-ray-to-gamma-ray flux ratio (8%) are both comparable to that of the two known gamma-ray transitional millisecond pulsars, we suspect that 3FGL J0212.1+5320 could be a potential target to search for future transition to the accretion active state.

AE Aquarii (AE Aqr) is a cataclysmic binary hosting one of the fastest rotating (P$_{\rm spin}$ = 33.08 s) white dwarfs known. Based on seven years of Fermi Large Area Telescope (LAT) Pass 8 data, we report on a deep search for gamma-ray emission from AE Aqr. Using X-ray observations from ASCA, XMM-Newton, Chandra, Swift, Suzaku, and NuSTAR, spanning 20 years, we substantially extend and improve the spin ephemeris of AE Aqr. Using this ephemeris, we searched for gamma-ray pulsations at the spin period of the white dwarf. No gamma-ray pulsations were detected above 3 $\sigma$ significance. Neither phase-averaged gamma-ray emission nor gamma-ray variability of AE Aquarii is detected by Fermi-LAT. We impose the most restrictive upper limit to the gamma-ray flux from AE Aqr to date: $1.3\times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$ in the 100 MeV-300 GeV energy range, providing constraints on models.

We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first 7 years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interactions of cosmic rays with the lunar surface. The results of the present analysis can be explained in the framework of this model, where the production of gamma rays is due to the interactions of cosmic-ray proton and helium nuclei with the surface of the Moon. Finally, we have used our simulation to derive the cosmic-ray proton and helium spectra near Earth from the Moon gamma-ray data.

Most of the celestial gamma rays detected by the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point and extended source studies rely on the modeling of this diffuse emission for accurate characterization. We describe here the development of the Galactic Interstellar Emission Model (GIEM) that is the standard adopted by the LAT Collaboration and is publicly available. The model is based on a linear combination of maps for interstellar gas column density in Galactocentric annuli and for the inverse Compton emission produced in the Galaxy. We also include in the GIEM large-scale structures like Loop I and the Fermi bubbles. The measured gas emissivity spectra confirm that the cosmic-ray proton density decreases with Galactocentric distance beyond 5 kpc from the Galactic Center. The measurements also suggest a softening of the proton spectrum with Galactocentric distance. We observe that the Fermi bubbles have boundaries with a shape similar to a catenary at latitudes below 20 degrees and we observe an enhanced emission toward their base extending in the North and South Galactic direction and located within 4 degrees of the Galactic Center.

To uniformly determine the properties of supernova remnants (SNRs) at high energies, we have developed the first systematic survey at energies from 1 to 100 GeV using data from the Fermi Large Area Telescope. Based on the spatial overlap of sources detected at GeV energies with SNRs known from radio surveys, we classify 30 sources as likely GeV SNRs. We also report 14 marginal associations and 245 flux upper limits. A mock catalog in which the positions of known remnants are scrambled in Galactic longitude, allows us to determine an upper limit of 22% on the number of GeV candidates falsely identified as SNRs. We have also developed a method to estimate spectral and spatial systematic errors arising from the diffuse interstellar emission model, a key component of all Galactic Fermi LAT analyses. By studying remnants uniformly in aggregate, we measure the GeV properties common to these objects and provide a crucial context for the detailed modeling of individual SNRs. Combining our GeV results with multiwavelength (MW) data, including radio, X-ray, and TeV, demonstrates the need for improvements to previously sufficient, simple models describing the GeV and radio emission from these objects. We model the GeV and MW emission from SNRs in aggregate to constrain their maximal contribution to observed Galactic cosmic rays.

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

We report the discovery of PSR J1906+0722, a gamma-ray pulsar detected as part of a blind survey of unidentified Fermi Large Area Telescope (LAT) sources being carried out on the volunteer distributed computing system, Einstein@Home. This newly discovered pulsar previously appeared as the most significant remaining unidentified gamma-ray source without a known association in the second Fermi-LAT source catalog (2FGL) and was among the top ten most significant unassociated sources in the recent third catalog (3FGL). PSR J1906+0722 is a young, energetic, isolated pulsar, with a spin frequency of $8.9$ Hz, a characteristic age of $49$ kyr, and spin-down power $1.0 \times 10^{36}$ erg s$^{-1}$. In 2009 August it suffered one of the largest glitches detected from a gamma-ray pulsar ($\Delta f / f \approx 4.5\times10^{-6}$). Remaining undetected in dedicated radio follow-up observations, the pulsar is likely radio-quiet. An off-pulse analysis of the gamma-ray flux from the location of PSR J1906+0722 revealed the presence of an additional nearby source, which may be emission from the interaction between a neighboring supernova remnant and a molecular cloud. We discuss possible effects which may have hindered the detection of PSR J1906+0722 in previous searches and describe the methods by which these effects were mitigated in this survey. We also demonstrate the use of advanced timing methods for estimating the positional, spin and glitch parameters of difficult-to-time pulsars such as this.

In this work, we investigate the rotational dynamics of the ginger-shaped near-Earth asteroid 4179 Toutatis, which was closely observed by Chang'e-2 at a distance of $770\pm120~$ meters from the asteroid's surface during the outbound flyby \citep{Huang2013} on 13 December 2012. A sequence of high-resolution images was acquired during the flyby mission. In combination with ground-based radar observations collected over the last two decades, we analyze these flyby images and determine the orientation of the asteroid at the flyby epoch. The 3-1-3 Euler angles of the conversion matrix from the J2000 ecliptic coordinate system to the body-fixed frame are evaluated to be $-20.1^\circ\pm1^\circ$, $27.6^\circ\pm1^\circ$ and $42.2^\circ\pm1^\circ$, respectively. The least-squares method is utilized to determine the rotational parameters and spin state of Toutatis. The characteristics of the spin-state parameters and angular momentum variations are extensively studied using numerical simulations, which confirm those reported by \citet{Takahashi2013}. The large amplitude of Toutatis' precession is assumed to be responsible for its tumbling attitude as observed from Earth. Toutatis' angular momentum orientation is determined to be described by $\lambda_{H}=180.2^{+0.2^\circ}_{-0.3^\circ}$ and $\beta_{H}=-54.75^{+0.15^\circ}_{-0.10^\circ}$, implying that it has remained nearly unchanged for two decades. Furthermore, using Fourier analysis to explore the change in the orientation of Toutatis' axes, we reveal that the two rotational periods are 5.38 and 7.40 days, respectively, consistent with the results of the former investigation. Hence, our investigation provides a clear understanding of the state of the rotational dynamics of Toutatis.

Context: GeV gamma-ray pulsations from over 140 pulsars have been characterized using the Fermi Large Area Telescope, enabling improved understanding of the emission regions within the neutron star magnetospheres, and the contributions of pulsars to high energy electrons and diffuse gamma rays in the Milky Way. The first gamma-ray pulsars to be detected were the most intense and/or those with narrow pulses. Aims: As the Fermi mission progresses, progressively fainter objects can be studied. In addition to more distant pulsars (thus probing a larger volume of the Galaxy), or ones in high background regions (thus improving the sampling uniformity across the Galactic plane), we detect pulsars with broader pulses or lower luminosity. Adding pulsars to our catalog with inclination angles that are rare in the observed sample, and/or with lower spindown power, will reduce the bias in the currently known gamma-ray pulsar population. Methods: We use rotation ephemerides derived from radio observations to phase-fold gamma rays recorded by the Fermi Large Area Telescope, to then determine the pulse profile properties. Spectral analysis provides the luminosities and, when the signal-to-noise ratio allows, the cutoff energies. We constrain the pulsar distances by different means in order to minimize the luminosity uncertainties. Results: We present six new gamma-ray pulsars with an eclectic mix of properties. Three are young, and three are recycled. They include the farthest, the lowest power, two of the highest duty-cycle pulsars seen, and only the fourth young gamma-ray pulsar with a radio interpulse. We discuss the biases existing in the current gamma-ray pulsar catalog, and steps to be taken to mitigate the bias.

Millisecond pulsars (MSPs) are a growing class of gamma-ray emitters. Pulsed gamma-ray signals have been detected from more than 40 MSPs with the Fermi Large Area Telescope (LAT). The wider radio beams and more compact magnetospheres of MSPs enable studies of emission geometries over a broader range of phase space than non-recycled radio-loud gamma-ray pulsars. We have modeled the gamma-ray light curves of 40 LAT-detected MSPs using geometric emission models assuming a vacuum retarded-dipole magnetic field. We modeled the radio profiles using a single-altitude hollow-cone beam, with a core component when indicated by polarimetry; however, for MSPs with gamma-ray and radio light curve peaks occurring at nearly the same rotational phase we assume that the radio emission is co-located with the gamma rays and caustic in nature. The best-fit parameters and confidence intervals are determined using a maximum likelihood technique. We divide the light curves into three model classes, with gamma-ray peaks trailing (Class I), aligned (Class II) or leading (Class III) the radio peaks. Outer gap and slot gap (two-pole caustic) models best fit roughly equal numbers of Class I and II, while Class III are exclusively fit with pair-starved polar cap models. Distinguishing between the model classes based on typical derived parameters is difficult. We explore the evolution of magnetic inclination angle with period and spin-down power, finding possible correlations. While the presence of significant off-peak emission can often be used as a discriminator between outer gap and slot gap models, a hybrid model may be needed.

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

Over thirteen times more gamma-ray pulsars have now been studied with the Large Area Telescope on NASA's Fermi satellite than the ten seen with the Compton Gamma-Ray Observatory in the nineteen-nineties. The large sample is diverse, allowing better understanding both of the pulsars themselves and of their roles in various cosmic processes. Here we explore the prospects for even more gamma-ray pulsars as Fermi enters the 2nd half of its nominal ten-year mission. New pulsars will naturally tend to be fainter than the first ones discovered. Some of them will have unusual characteristics compared to the current population, which may help discriminate between models. We illustrate a vision of the future with a sample of six pulsars discovered after the 2nd Fermi Pulsar Catalog was written.

Pulsar wind nebulae (PWNe) have been established as the most populous class of TeV gamma-ray emitters. Since launch, the Fermi Large Area Telescope (LAT)identified five high-energy (100MeV <E< 100 GeV) gamma-ray sources as PWNe, and detected a large number of PWNe candidates, all powered by young and energetic pulsars. The wealth of multi-wavelength data available and the new results provided by Fermi-LAT give us an opportunity to find new PWNe and to explore the radiative processes taking place in known ones. The TeV gamma-ray unidentifiedsources (UNIDs) are the best candidates for finding new PWNe. Using 45 months of Fermi-LAT data for energies above 10 GeV, an analysis was performed near the position of 58TeV PWNe and UNIDs within 5deg of the Galactic Plane to establish new constraints on PWNe properties and find new clues on the nature of UNIDs. Of the 58 sources, 30 were detected, and this work provides their gamma-rayfluxes for energies above 10 GeV. The spectral energy distributions (SED) andupper limits, in the multi-wavelength context, also provide new information on the source nature and can help distinguish between emission scenarios, i.e. between classification as a pulsar candidate or as a PWN candidate. Six new GeV PWNe candidates are described in detail and compared with existing models. A population study of GeV PWNe candidates as a function of the pulsar/PWN system characteristics is presented.

We report on the discovery of gamma-ray pulsations from five millisecond pulsars (MSPs) using the Fermi Large Area Telescope (LAT) and timing ephemerides provided by various radio observatories. We also present confirmation of the gamma-ray pulsations from a sixth source, PSR J2051-0827. Five of these six MSPs are in binary systems: PSRs J1713+0747, J1741+1351, J1600-3053 and the two black widow binary pulsars PSRs J0610-2100 and 2051-0827. The only isolated MSP is the nearby PSR J1024-0719, which is also known to emit X-rays. We present X-ray observations in the direction of PSRs J1600-3053 and J2051-0827. While the latter is firmly detected, we an only give upper limits for the X-ray flux of the former. There are no dedicated X-ray observations available for the other 3 objects. The MSPs mentioned above, together with most of the MSPs detected by Fermi, are used to put together a sample of 30 gamma-ray MSPs which is used to study the morphology and phase connection of radio and gamma-ray pulse profiles. We show that MSPs with pulsed gamma-ray emission which is phase aligned with the radio emission present the steepest radio spectra and the largest magnetic fields at the light cylinder among all MSPs. As well, we also observe a trend towards very low, or undetectable, radio linear polarisation levels. These properties could be attributed to caustic radio emission produced at a range of different altitudes in the magnetosphere. We note that most of these characteristics are also observed in the Crab pulsar, the only other radio pulsar known to exhibit phase-aligned radio and gamma-ray emission.

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

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

The contribution of unresolved sources to the diffuse gamma-ray background could induce anisotropies in this emission on small angular scales. We analyze the angular power spectrum of the diffuse emission measured by the Fermi LAT at Galactic latitudes |b| > 30 deg in four energy bins spanning 1 to 50 GeV. At multipoles \ell \ge 155, corresponding to angular scales \lesssim 2 deg, angular power above the photon noise level is detected at >99.99% CL in the 1-2 GeV, 2-5 GeV, and 5-10 GeV energy bins, and at >99% CL at 10-50 GeV. Within each energy bin the measured angular power takes approximately the same value at all multipoles \ell \ge 155, suggesting that it originates from the contribution of one or more unclustered source populations. The amplitude of the angular power normalized to the mean intensity in each energy bin is consistent with a constant value at all energies, C_P/<I>^2 = 9.05 +/- 0.84 x 10^{-6} sr, while the energy dependence of C_P is consistent with the anisotropy arising from one or more source populations with power-law photon spectra with spectral index \Gamma_s = 2.40 +/- 0.07. We discuss the implications of the measured angular power for gamma-ray source populations that may provide a contribution to the diffuse gamma-ray background.

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

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