The microquasar V404 Cygni (also known as GS 2023+338) was previously reported to have possible GeV gamma-ray emission in two days during its 2015 outburst. In order to provide more detailed information at the high energy range for this black hole binary system, we conduct detailed analysis to the data obtained with the Large Area Telescope (LAT) onboard the Fermi Gamma Ray Space Telescope (Fermi). Both LAT database and source catalog used are the latest. From the analysis, we can not confirm the previous detection. Instead, we find one possible detection (~4 sigma) of the source at the end of the outburst in the time period of 2015 Aug. 17-19, and one convincing detection (~7 sigma) in 2016 Aug. 23-25. The latter shows that the gamma-ray emission of the source is soft with photon index $\Gamma\sim 2.9$, mostly detected below ~1.3 GeV with Fermi LAT. As gamma-ray emission from microquasars is likely associated with their jet activity, we discuss the results by comparing with those well studied cases, namely Cyg X-3 and Cyg X-1. The detection establishes V404 Cygni as one of four microquasars with detectable gamma-ray emission, and adds interesting features to the small group, or in a more general context to X-ray binaries with jets.

We detail the sensitivity of the liquid xenon (LXe) DARWIN observatory to solar neutrinos via elastic electron scattering. We find that DARWIN will have the potential to measure the fluxes of five solar neutrino components: $pp$, $^7$Be, $^{13}$N, $^{15}$O and $pep$. The precision of the $^{13}$N, $^{15}$O and $pep$ components is hindered by the double-beta decay of $^{136}$Xe and, thus, would benefit from a depleted target. A high-statistics observation of $pp$ neutrinos would allow us to infer the values of the weak mixing angle, $\sin^2\theta_w$, and the electron-type neutrino survival probability, $P_e$, in the electron recoil energy region from a few keV up to 200 keV for the first time, with relative precision of 5% and 4%, respectively, at an exposure of 300 ty. An observation of $pp$ and $^7$Be neutrinos would constrain the neutrino-inferred solar luminosity down to 0.2%. A combination of all flux measurements would distinguish between the high (GS98) and low metallicity (AGS09) solar models with 2.1-2.5$\sigma$ significance, independent of external measurements from other experiments or a measurement of $^8$B neutrinos through coherent elastic neutrino-nucleus scattering in DARWIN. Finally, we demonstrate that with a depleted target DARWIN may be sensitive to the neutrino capture process of $^{131}$Xe.

This paper presents a class of novel high-order accurate discontinuous Galerkin (DG) schemes for the compressible Euler equations under gravitational fields. A notable feature of these schemes is that they are well-balanced for a general hydrostatic equilibrium state, and at the same time, provably preserve the positivity of density and pressure. In order to achieve the well-balanced and positivity-preserving properties simultaneously, a novel DG spatial discretization is carefully designed with suitable source term reformulation and a properly modified Harten-Lax-van Leer contact (HLLC) flux. Based on some technical decompositions as well as several key properties of the admissible states and HLLC flux, rigorous positivity-preserving analyses are carried out. It is proven that the resulting well-balanced DG schemes, coupled with strong stability preserving time discretizations, satisfy a weak positivity property, which implies that one can apply a simple existing limiter to effectively enforce the positivity-preserving property, without losing high-order accuracy and conservation. The proposed methods and analyses are applicable to the Euler system with general equation of state. Extensive one- and two-dimensional numerical tests demonstrate the desired properties of these schemes, including the exact preservation of the equilibrium state, the ability to capture small perturbation of such state, the robustness for solving problems involving low density and/or low pressure, and good resolution for smooth and discontinuous solutions.

4FGL J0935.3+0901 is a gamma-ray source detected by the Large Area Telescope (LAT) onboard the Fermi Gamma-Ray Space Telescope. We have conducted detailed analysis of the LAT data for this source and multi-wavelength studies of the source field. Its gamma-ray emission can be described with a power law ($\Gamma=2.0\pm0.2$) with an exponential cutoff ($E_c=2.9\pm1.6$ GeV), while the flux shows significant long-term variations. From analysis of archival Neil Gehrels Swift X-ray Telescope data, we find only one X-ray source in the LAT's 2$\sigma$ error region. Within a 3.7" radius error circle of the X-ray source, there is only one optical object down to $r'\sim$23 mag. Time-resolved photometry of the optical object indicates a likely 2.5~hr periodic modulation, while its spectrum shows double-peaked hydrogen and helium emission lines (similar to those seen in accretion discs in low-mass X-ray binaries). Combining these results, we conclude that we have discovered a compact X-ray emitting binary in likely association with 4FGL J0935.3+0901, i.e., a millisecond pulsar (MSP) binary. We discuss the implication of the optical spectral features: this binary could be a transitional MSP system at a sub-luminous disc state, although the other possibility, the binary in a rotation-powered state showing the optical emission lines due to intrabinary interaction processes, can not be excluded. Further observational studies will help determine detailed properties of this candidate MSP binary and thus clarify its current state.

We search through $\gamma$-ray data obtained with the Large Area Telescope (LAT) onboard the Fermi Gamma-Ray Space Telescope and find 24 blazars (or candidates) that have a single clear flare event in their 9.5 year long-term light curves. We define these events as loner flares since each flare stands out significantly above the relatively stable, low-flux light curve. We analyze the LAT data in detail for these 24 sources. The flares in ten of them are primarily due to a single sharp peak, for which we study by fitting with two different analytic functions. The time durations thus determined for the sharp peaks are in a range of 4-25 days. The $\gamma$-ray spectra of the 24 blazar sources can be described with a power-law or a log-parabola function. We obtain their spectral properties in the flaring and quiescent states, and find that in the flares 16 of the sources have harder emission and three have softer emission while the other five keep the same emission. We discuss a possible correlation between the differences in photon index in the quiescent and flaring states and photon indices in quiescence. In addition, the sharp peak flares seem to have a tendency of having long time durations and hard emission, possibly related to their physical origin in a blazar jet. Studies of more similar flares will help establish these possible features.

In our search for gamma-ray emission from Be X-ray binaries from analysis of the data obtained with the Large Area Telescope (LAT) on board the Fermi Gamma-Ray Space Telescope, we find likely detection of GRO J1008-57. The binary has an orbital period of 249.48 days, and it is only significantly detected in its orbital phase 0.8--0.9 (>4 sigma). Further light curve analysis indicates that the detection is probably largely due to an emitting event in one orbital cycle around year 2012--2013, following a giant X-ray outburst of the source. This property of having occasional gamma-ray emitting events is similar to that seen in another high-mass X-ray binary 4U 1036-56. However, models considering possible gamma-ray emission from an accreting neutron star have difficulty in explaining the observed 10^{34} erg s^{-1} luminosity of the source, unless the distance was largely over-estimated. Further observational studies are required, in order to more clearly establish the high-energy emission properties of GRO J1008-57 or similar high-mass X-ray binaries and find clues for understanding how gamma-ray emission is produced from them.

The data from the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope have recently been updated. We thus re-analyze the LAT data for the supernova remnant (SNR) SN 1006. Two parts of gamma-ray emission from the region is clearly resolved, which correspond to the north-east (NE) and south-west (SW) limbs of the SNR. The former was detected in the previous LAT data (Xing et al. 2016), but the latter is newly detected in this work. The detection of the two limbs are at a 4 sigma significance level, and the spectral results for the NE limb is consistent with those obtained in previous detection analyses. We construct the broadband spectral energy distribution (SED) for the SW limb. Different scenarios are considered for the SED in gamma-ray energies. We conclude very similar to that of the NE limb, the high-energy and very high-energy emission from the SW limb is likely dominated by the leptonic process, in which high-energy electrons accelerated from the shell region of the SNR inverse-Compton scatter background photons to gamma-rays.

We report our analysis of the data obtained with the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope (Fermi) for the SS 433/W50 region. The total data length is ten years. Different from the previous results reported for this region (for which an old version of database was used), we show that excess emission is detected around the w1 region in the western lobe of the jets from SS 433. The region is bright at X-rays due to the interaction between the jet and the ambient medium. This detection also matches well the recent results of the very-high-energy detection of SS 433/W50 with the High Altitude Water Cherenkov (HAWC). However, the eastern regions that are slightly brighter in HAWC's observation are not detected in the Fermi data. Constructing the broad-band spectral energy distribution (SED) for the western region, we compare with the HAWC results for the eastern regions and discuss the possible origin of the emission. A leptonic scenario can provide a fit to the Fermi GeV spectrum and HAWC TeV detection, where the former and latter are due to the synchrotron radiation and inverse-Compton process respectively. However, the model can not explain the X-ray and radio emission from the region simultaneously, which thus requires further observational and theoretical studies of the region in order to clarify the reasons.

Using isochronous mass spectrometry at the experimental storage ring CSRe in Lanzhou, the masses of $^{82}$Zr and $^{84}$Nb were measured for the first time with an uncertainty of $\sim 10$ keV, and the masses of $^{79}$Y, $^{81}$Zr, and $^{83}$Nb were re-determined with a higher precision. %The latter differ significantly from their literature values. The latter are significantly less bound than their literature values. Our new and accurate masses remove the irregularities of the mass surface in this region of the nuclear chart. Our results do not support the predicted island of pronounced low $\alpha$ separation energies for neutron-deficient Mo and Tc isotopes, making the formation of Zr-Nb cycle in the $rp$-process unlikely. The new proton separation energy of $^{83}$Nb was determined to be 490(400)~keV smaller than that in the Atomic Mass Evaluation 2012. This partly removes the overproduction of the $p$-nucleus $^{84}$Sr relative to the neutron-deficient molybdenum isotopes in the previous $\nu p$-process simulations.

The mean white dwarf (WD) mass in the Galactic bulge cataclysmic variables (CVs) was measured by applying the the shock temperature-WD mass correlation of magnetic cataclysmic variables (mCVs) to the Galactic Bulge X-ray Emission (GBXE) spectra. However, the resulting mean WD mass is lower than that of the local CVs. This discrepancy could be explained by the dominating sources in the GBXE are non-mCVs instead of mCVs. In this work, we conduct an thorough investigation on the X-ray spectra of local DNe from \suzaku\ archives, and derive semi-empirical correlations between the shock temperature $T_{\rm max}$, the flux ratio of \fec\ to \feb\ lines, and WD mass for quiescent, non-magnetic CVs. By applying these correlations to the GBXE, we derive the average WD mass of CVs in the Galactic bulge to be $0.81\pm 0.07M_\odot$. This value is consistent with previous optical measurements of WD mass in local CVs.

We report the results from our analysis of 8 years of the data for the gamma-ray binary LS I +61 303, obtained with the Large Area Telescope onboard the Fermi Gamma-Ray Space Telescope. We find a significant dip around the binary's periastron in the superorbital light curves, and by fitting the light curves with a sinusoidal function, clear phase shifts are obtained. The superorbital modulation seen in the binary has been long known and different scenarios have been proposed. Based on our results, we suggest that the circumstellar disk around the Be companion of this binary has an elliptical shape and the major axis of the disk rotates at the superorbital period of 1667 days. As a result, the density of the ambient material around the compact star of the binary changes along the binary orbit over the superorbital period, causing the phase shifts in the modulation, and around periastron, the compact star probably enters the disk, causing the appearance of the dip. Numerical simulations may be conducted in order to study the detailed physical processes and verify our suggestion.

Following our work presented in Dai et al. (2016), we report our detailed data analysis for another 38 Fermi gamma-ray un-associated sources. These sources are selected from the Fermi Large Area Telescope (LAT) third source catalog on the basis of the properties of known gamma-ray millisecond pulsars (MSPs) and for the purpose of finding likely candidate MSPs. From our analysis of the LAT data, we identify that among the 38 sources, 28 of them are single point-like sources with clean background and their spectra show significant curvature. We also conduct analysis of archival X-ray data available for 24 of the 28 sources. In the fields of 10 sources, there are at least one X-ray object, and in those of the other 14 sources, no X-ray object is detected but probably due to the X-ray observations being short. We discuss the possible MSP nature for these sources. Six of them(J0514.6-4406, J1035.7-6720, J1624.2-4041, J1744.1-7619, J1946.4-5403, and J2039.6-5618) are most likely associated with pulsars because of multi-wavelength identifications including direct radio or gamma-ray detection of pulsations. To firmly establish the associations or verify the MSP nature for other sources, deep X-ray and/or optical observations are needed.

We report our analysis of the Fermi Large Area Telescope data for 39 millisecond pulsars (MSPs) listed in the second $\gamma$-ray pulsar catalog. Spectra of the pulsars are obtained. We fit the spectra with a function of a power law with exponential cutoff, and find the best-fit parameters of photon index $\Gamma = 1.54^{+0.10}_{-0.11}$ and cutoff energy $E_{c} = 3.70^{+0.95}_{-0.70}$ GeV. This spectral shape, which includes the intrinsic differences in the spectra of the MSPs, can be used for finding candidate MSPs and unidentified types of sources detected by Fermi at high Galactic latitudes. In one of the MSPs PSR J0614-3329, we find significant pulsed emission upto 200 GeV. The result has thus added this MSP to the group of the Crab and Vela pulsars that have been detected with >50 GeV pulsed emission. Comparing the $\gamma$-ray spectrum of PSR J0614-3329 with those of the Crab and Vela pulsars, we discuss possible emission mechanisms for the very high-energy component.

We report the results from our detailed analysis of the Fermi Large Area Telescope (LAT) data for the pulsar binary PSR B1259-63/LS 2883. During the GeV flares that occurred when the pulsar was in the periastron passages, we have detected a 5--300 GeV component at $\simeq 5\sigma$ in emission from the binary. The detection verifies the presence of the component that has been marginally found in the previous studies of the binary. Furthermore, we have discovered that this component was present even in the quiescent state of the binary, specifically the mean anomaly phase 0.7--0.9. The component can be described by a power law with photon index $\Gamma\sim 1.4$, and the flux in the flares is approximately one order of magnitude higher than that in quiescence. We discuss the origin of this component. It likely arises from the inverse-Compton process: high-energy particles from the shock, due to the interaction between the winds from the pulsar and massive companion, scatter the seed photons from the companion to GeV/TeV energies. Based on this scenario, model fits to the broad-band X-ray--to--TeV spectra of the binary in the flaring and quiescent states are provided.

We report the likely detection of gamma-ray emission from the northeast shell region of the historical supernova remnant (SNR) SN 1006. Having analyzed 7 years of Fermi LAT Pass 8 data for the region of SN 1006, we found a GeV gamma-ray source detected with 4 sigma significance. Both the position and spectrum of the source match those of HESS J1504-418 respectively, which is TeV emission from SN 1006. Considering the source as the GeV gamma-ray counterpart to SN~1006, the broadband spectral energy distribution is found to be approximately consistent with the leptonic scenario that has been proposed for the TeV emission from the SNR. Our result has likely confirmed the previous study of the SNRs with TeV shell-like morphology: SN 1006 is one of them sharing very similar peak luminosity and spectral shape.

We report our detailed data analysis for 39 $\gamma$-ray sources selected from the 992 unassociated sources in the \textit{Fermi} Large Area Telescope (LAT) third source catalog. The selection criteria, which were set for finding candidate millisecond pulsars (MSPs), are non-variables with curved spectra and $>$5$^{\circ}$ Galactic latitudes. From our analysis, 24 sources were found to be point-like sources not contaminated by background or nearby unknown sources. Three of them, J1544.6$-$1125, J1625.1$-$0021, and J1653.6$-$0158, have been previously studied, indicating that they are likely MSPs. The spectra of J0318.1+0252 and J2053.9+2922 do not have properties similar to that of known $\gamma$-ray MSPs, and we thus suggest that they are not MSPs. Analysis of archival X-ray data for most of the 24 sources were also conducted. Four sources were found with X-ray objects in their error circles, and 16 with no detection. The ratios between the $\gamma$-ray fluxes and X-ray fluxes or flux upper limits are generally lower than those of the known $\gamma$-ray MSPs, suggesting that if the $\gamma$-ray sources are MSPs, none of the X-ray objects are the counterparts. Deep X-ray or radio observations of these sources are needed in order to identify their MSP nature.

Hadron emission from supernova remnant (SNR)-molecular cloud (MC) association systems has been widely regarded as a probe of the shock-accelerated cosmic-ray protons. Here, we report on the detection of a {\gamma}-ray emission source, with a significance of 24{\sigma} in 0.2-300 GeV, projected to lie to the northwest of the thermal composite SNR Kesteven 41, using 5.6 years of Fermi-Large Area Telescope (LAT) observation data. No significant long-term variability in the energy range 0.2--300 GeV is detected around this source. The 3{\sigma} error circle, 0.09 degree; in radius, covers the 1720MHz OH maser and is essentially consistent with the location of the V_{LSR} ~-50 km/s MC with which the SNR interacts. The source emission has an exponential cutoff power-law spectrum with a photon index of 1.9+/-0.1 and a cutoff energy of 4.0+/-0.9 GeV, and the corresponding 0.2-300 GeV luminosity is ~1.3*10^36 erg/s at a distance of 12 kpc. There is no radio pulsar in the 3{\sigma} circle responsible for the high {\gamma}-ray luminosity. While the inverse Compton scattering scenario would lead to a difficulty in the electron energy budget, the source emission can naturally be explained by the hadronic interaction between the relativistic protons accelerated by the shock of SNR Kesteven 41 and the adjacent northwestern MC. In this paper, we present a list of Galactic thermal composite SNRs detected at GeV {\gamma}-ray energies by Fermi-LAT is presented.

In de Rham-Gabadadze-Tolley (dRGT) massive gravity and bigravity, a non-minimal matter coupling involving both metrics generically re-introduces the Boulware--Deser (BD) ghost. A non-minimal matter coupling via a simple, yet specific composite metric has been proposed, which eliminates the BD ghost below the strong coupling scale. Working explicitly in the metric formulation and for arbitrary spacetime dimensions, we show that this composite metric is the unique consistent non-minimal matter coupling below the strong coupling scale, which emerges out of two diagnostics, namely, absence of Ostrogradski ghosts in the decoupling limit and absence of the BD ghost from matter quantum loop corrections.

We report our discovery of orbitally modulated $\gamma$-ray emission from the black widow system PSR J1311-3430. We analyze the \textit{Fermi} Large Area Telescope data during the offpulse phase interval of the pulsar, and find the orbital modulation signal at a $\sim$3$\sigma$ confidence level. Further spectral analysis shows no significant differences for the spectra obtained during the bright and faint orbital phase ranges. A simple sinusoid-like function can describe the modulation. Given these properties, we suggest that the intrabinary $\gamma$-ray emission arises from the region close to the companion and the modulation is caused by the occultation of the emitting region by the companion, similar to that is seen in the transitional millisecond pulsar binary (MSP) PSR J1023+0038. Considering the X-ray detection of intrabinary shock emission from eclipsing MSP binaries recently reported, this discovery further suggests the general existence of intrabinary $\gamma$-ray emission from them.

We report on the likely detection of gamma-ray emission from the supernova remnant (SNR) Kesteven 27 (Kes 27). We analyze 5.7 yr Fermi Large Area Telescope data of the SNR region and find an unresolved source at a position consistent with the radio brightness peak and the X-ray knot of Kes 27, which is located in the eastern region of the SNR and caused by the interaction with a nearby HI cloud. The source's emission is best fit with a power-law spectrum with a photon index of 2.5 +/- 0.1 and a >0.2 GeV luminosity of 5.8*10^34 erg s^-1 assuming a distance of 4.3 kpc, as derived from radio observations of the nearby HI cloud. Comparing the properties of the source with that of other SNRs that are known to be interacting with nearby high-density clouds, we discuss the origin of the source's emission. The spectral energy distribution of the source can be described by a hadronic model that considers the interaction of energetic protons, escaping from the shock front of Kes 27, with a high-density cloud.

We report the Fermi Large Area Telescope (LAT) detection of a $\gamma$-ray source at the position of SAX J1808.4$-$3658. This transient low-mass X-ray binary contains an accreting millisecond puslar, which is only seen during its month-long outbursts and likely switches to be rotation powered during its quiescent state. Emission from the $\gamma$-ray source can be described by a power law with an exponential cutoff, the characteristic form for pulsar emission. Folding the source's 2.0--300 GeV photons at the binary orbital period, a weak modulation is seen (with an H-test value of $\sim$17). In addition, three sets of archival XMM-Newton data for the source field are analyzed, and we find only one X-ray source with 3--4$\sigma$ flux variations in the 2$\sigma$ error circle of the $\gamma$-ray source. However based on the X-ray properties, this X-ray source is not likely a background AGN, the major class of Fermi sources detected by LAT. These results support the possible association between the $\gamma$-ray source and SAX J1808.4$-$3658 and thus the scenario that the millisecond pulsar is rotation powered in the quiescent state. Considering a source distance of 3.5 kpc for SAX J1808.4$-$3658, the 0.1--300 GeV luminosity is 5.7$\times 10^{33}$ erg s$^{-1}$, implying a $\gamma$-ray conversion efficiency of 63\% for the pulsar in this binary.

Because of the disappearance of its accretion disk since the time period around 2012 November--December, XSS J12270-4859 has recently been identified as, in addition to PSR J1023+0038, another transitional millisecond pulsar binary. We have carried out detailed analysis of the Fermi Large Area Telescope data for the binary. While both spectra before and after the disk-disappearance transition are well described by an exponentially cut-off power law, typical for pulsars' emission in the Fermi's 0.2-300 GeV, a factor of 2 flux decrease related to the transition is detected. A weak orbital modulation is seen, but only detectable in the after-transition data, same to that found at X-rays. In the long-term light curve of the source before the transition, a factor of 3 flux variations are seen. Comparing to the properties of J1023+0038, we disucss the implications from these results. We suggest that since the modulation is aligned with that at X-rays in orbital phase, it possibly arises due to the occultation of the gamma-ray emitting region by the companion. The origin of the variations in the long-term light curve is not clear, because the source field also contains unidentified radio or X-ray sources and their contamination can not be excluded. Multi-wavelength observations of the source field will help identify the origin of the variations by detecting any related flux changes from the in-field sources.

We extend the positivity-preserving method of Zhang & Shu (2010, JCP, 229, 3091-3120) to simulate the advection of neutral particles in phase space using curvilinear coordinates. The ability to utilize these coordinates is important for non-equilibrium transport problems in general relativity and also in science and engineering applications with specific geometries. The method achieves high-order accuracy using Discontinuous Galerkin (DG) discretization of phase space and strong stability-preserving, Runge-Kutta (SSP-RK) time integration. Special care in taken to ensure that the method preserves strict bounds for the phase space distribution function $f$; i.e., $f\in[0,1]$. The combination of suitable CFL conditions and the use of the high-order limiter proposed in Zhang & Shu (2010) is sufficient to ensure positivity of the distribution function. However, to ensure that the distribution function satisfies the upper bound, the discretization must, in addition, preserve the divergence-free property of the phase space flow. Proofs that highlight the necessary conditions are presented for general curvilinear coordinates, and the details of these conditions are worked out for some commonly used coordinate systems (i.e., spherical polar spatial coordinates in spherical symmetry and cylindrical spatial coordinates in axial symmetry, both with spherical momentum coordinates). Results from numerical experiments --- including one example in spherical symmetry adopting the Schwarzschild metric --- demonstrate that the method achieves high-order accuracy and that the distribution function satisfies the maximum principle.

The Fermi source 2FGL~J0523.3$-$2530 has recently been identified as a candidate millisecond pulsar binary with an orbital period of 16.5 hrs. We have carried out detailed studies of the source's emission properties by analyzing data taken with the Fermi Large Area Telescope in the 0.2--300 GeV energy range. Long-term, yearly variability from the source has been found, with a factor of 4 flux variations in 1--300 GeV. From spectral analysis, we find an extra spectral component at 2--3 GeV that causes the source brightening. While no orbital modulations have been found from the Fermi data over the whole period of 2008--2014, orbital modulation in the source's $>$2 GeV emission is detected during the last 1.5 yrs of the Fermi observation. Our results support the millisecond pulsar binary nature of 2FGL~J0523.3$-$2530. Multi-wavelength observations of the source are warranted in order to find any correlated flux variations and thus help determine the origin of the long-term variability, which currently is not understood.

We report the results from our analysis of the \textit{Fermi} Large Area Telescope data of the \textit{Fermi} unassociated source 2FGL J1906.5$+$0720, which is a high-ranked candidate pulsar. In order to better study our target, we first update the ephemeris for PSR J1907$+$0602, which is used to help remove any possible contamination due to strong emission from this nearby pulsar. From our analysis, 2FGL J1906.5$+$0720 is confirmed to have a significant low energy cutoff at $\sim$ 1~GeV in its emission (14$\sigma$--18$\sigma$ significance), consistent with those seen in young pulsars. We search for pulsations but no spin frequency signals are found in a frequency range of 0.5--32 Hz. No single model can fully describe the source's overall \fermi\ $\gamma$-ray spectrum, and the reason for this is the excess emission detected at energies of $\geq$4 GeV. The high-energy component possibly indicates emission from a pulsar wind nebula, when considering 2FGL J1906.5$+$0720 as a young pulsar. We conclude that 2FGL J1906.5$+$0720 is likely a pulsar based on the emission properties we have obtained, and observations at other energies are needed in order to confirm its pulsar nature.

Recent observations strongly suggest that the millisecond pulsar binary PSR J1023+0038 has developed an accretion disk since 2013 June. We present a multi-wavelength analysis of PSR J1023+0038, which reveals that 1) its gamma-rays suddenly brightened within a few days in June/July 2013 and has remained at a high gamma-ray state for several months; 2) both UV and X-ray fluxes have increased by roughly an order of magnitude, and 3) the spectral energy distribution has changed significantly after the gamma-ray sudden flux change. Time variabilities associated with UV and X-rays are on the order of 100-500 seconds and 50-100 seconds, respectively. Our model suggests that a newly formed accretion disk due to the sudden increase of the stellar wind could explain the changes of all these observed features. The increase of UV is emitted from the disk, and a new component in gamma-rays is produced by inverse Compton scattering between the new UV component and pulsar wind. The increase of X-rays results from the enhancement of injection pulsar wind energy into the intra-binary shock due to the increase of the stellar wind. We also predict that the radio pulses may be blocked by the evaporated winds from the disk and the pulsar is still powered by rotation.

We report on the results from our $\gamma$-ray analysis of the supernova remnant (SNR) RCW 103 region. The data were taken with the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. An extended source is found at a position consistent with that of RCW 103, and its emission was only detected above 1 GeV (10$\sigma$ significance), having a power-law spectrum with a photon index of 2.0$\pm$0.1. We obtain its 1--300 GeV spectrum, and the total flux gives a luminosity of 8.3$\times 10^{33}$ erg s$^{-1}$ at a source distance of 3.3 kpc. Given the positional coincidence and property similarities of this source with other SNRs, we identify it as the likely Fermi $\gamma$-ray counterpart to RCW 103. Including radio measurements of RCW 103, the spectral energy distribution (SED) is modeled by considering emission mechanisms based on both hadronic and leptonic scenarios. We find that models in the two scenarios can reproduce the observed SED, while in the hadronic scenario the existence of SNR--molecular-cloud interaction is suggested as a high density of the target protons is required.

We report our search for \gamma-ray emission in the energy range from 100 MeV to 300 GeV from four Accreting Millisecond Pulsars (AMPs), SAX J1808.4-3658, IGR J00291+5934, XTE J1814-338, and XTE J0929-314. The data are from four-year observations carried out by Large Area Telescope (LAT) onboard the Fermi \gamma-ray Space Telescope. The AMPs were not detected, and their \gamma-ray luminosity upper limits we obtain are 5.1*10^33 ergs/s for SAX J1808.4-3658, 2.1*10^33 ergs/s for IGR J00291+5934, 1.2*10^34 ergs/s for XTE J1814-338, and 2.2*10^33 ergs/s for XTE J0929-314. We compare our results with \gamma-ray irradiation luminosities required for producing optical modulations seen from the companions in the AMPs, which has been suggested by Takata et al. (2012), and our upper limits have excluded \gamma-ray emission as the heating source in these systems except XTE J0929-314, the upper limit of which is not deep enough. Our results also do not support the model proposed by Takata et al. (2012) that relatively strong \gamma-ray emission could arise from the outer gap of a high-mass neutron star controlled by the photon-photon pair-creation for the AMPs. Two AMPs, SAX J1808.4-3658 and IGR J00291+5934, have the measurements of their spin-down rates, and we derive the upper limits of their \gamma-ray conversion efficiencies, which are 57% and 3%, respectively. We discuss the implications to the AMP systems by comparing the efficiency upper limit values with that of 20 \gamma-ray millisecond pulsars (MSP) detected by Fermi and the newly discovered transitional MSP binary J1023+0038.

We analyzed the observations of SGR 1806-20 performed with the \textit{Rossi X-ray Timing Explorer} (RXTE) during its 2004 giant flare. We studied the phase evolution of the sub-pulses identified in the X-ray waveform and found that the sub-pulses varied in phase with time and then gradually settled, which might indicate drifts of the emission regions in relative to the neutron star surface, or changes in the local emission geometry before the magnetic field became stable. The characteristic e-folding timescale of the phase drifts measured starting about 15 s following the initial flux spike are in the range between 37 s and 84 s. This leads to the first measurements of the characteristic timescale for the magnetic field of the neutron star to settle after a field reconfiguration during the giant flare.