Context: The Crab pulsar is a bright $\gamma$-ray source detected at photon energies up to $\sim$1~TeV. Its phase-averaged and phase-resolved $\gamma$-ray spectra below 10 GeV exhibit exponential cutoffs while those above 10~GeV apparently follow simple power-laws. Aims: We re-visit the $\gamma$-ray properties of the Crab pulsar with 10-year \emph{Fermi} Large Area Telescope (LAT) data in the range of 60~MeV--500~GeV. With the phase-resolved spectra, we investigate the origins and mechanisms responsible for the emissions. Methods: The phaseograms are reconstructed for different energy bands and further analysed using a wavelet decomposition. The phase-resolved energy spectra are combined with the observations of ground-based instruments like MAGIC and VERITAS to achieve a larger energy converage. We fit power-law models to the overlapping energy spectra from 10~GeV to $\sim$1~TeV. We include in the fit a relative cross-calibration of energy scales between air-shower based gamma-ray telescopes with the orbital pair-production telescope of the Fermi mission. Results: We confirm the energy-dependence of the $\gamma$-ray pulse shape, and equivalently, the phase-dependence of the spectral shape for the Crab pulsar. A relatively sharp cutoff at a relatively high energy of $\sim$8~GeV is observed for the bridge-phase emission. The $E>$10~GeV spectrum observed for the second pulse peak is harder than those for other phases. Conclusions: In view of the diversity of phase-resolved spectral shapes of the Crab pulsar, we tentatively propose a multi-origin scenario where the polar-cap, outer-gap and relativistic-wind regions are involved.

Context: The Crab pulsar and its nebula are the origin of relativistic electrons which can be observed through their synchrotron and inverse Compton emission. The transition between synchrotron-dominated and inverse-Compton-dominated emissions takes place at $\approx 10^9$ eV. Aims: The short-term (weeks to months) flux variability of the synchrotron emission from the most energetic electrons is investigated with data from ten years of observations with the Fermi Large Area Telescope (LAT) in the energy range from 60 MeV to 600 MeV. Methods: The off-pulse light-curve has been reconstructed from phase-resolved data. The corresponding histogram of flux measurements is used to identify distributions of flux-states and the statistical significance of a lower-flux component is estimated with dedicated simulations of mock light-curves. The energy spectra for different flux states are reconstructed. Results: We confirm the presence of flaring-states which follow a log-normal flux distribution. Additionally, we discover a low-flux state where the flux drops to as low as 18.4% of the intermediate-state average flux and stays there for several weeks. The transition time is observed to be as short as 2 days. The energy spectrum during the low-flux state resembles the extrapolation of the inverse-Compton spectrum measured at energies beyond several GeV energy, implying that the high-energy part of the synchrotron emission is dramatically depressed. Conclusions: The low-flux state found here and the transition time of at most 10 days indicate that the bulk ($>75$%) of the synchrotron emission above $10^8$ eV originates in a compact volume with apparent angular size of $\theta\approx0.4" t_\mathrm{var}/(5 \mathrm{d})$. We tentatively infer that the so-called inner knot feature is the origin of the bulk of the $\gamma$-ray emission.

The Crab Nebula is a bright emitter of non-thermal radiation across the entire accessible range of wavelengths. The spatial and spectral structures of the synchrotron nebula are well-resolved from radio to hard X-ray emission. The un-pulsed emission at GeV to TeV energies is mostly produced via inverse-Compton scattering of energetic electrons with the synchrotron-emitted photons. The spatial structure observed at these energies provides insights into the distribution of electrons and indirectly constrains the so-far unknown structure of the magnetic field in the nebula. Analyzing the LAT data accumulated over $\sim$9.1 years with a properly refined model for the Crab pulsar's spectrum, we determined the 68\% containment radius ($R_{68}$) of the Crab Nebula to be $({0.0330\pm0.0025_{stat}}{^{+0.0012}_{-0.0075}}_{sys})^\circ$ (${1.98'\pm0.15'_{stat}}{^{+0.07'}_{-0.45'}}_{sys}$) in the 5--500 GeV band. The estimated systematic uncertainty is based on two factors: (1) different analysis methods, morphological models and event types, and (2) the point-spread-function evaluated with observations of Mkn 421. When comparing the \emph{Fermi} LAT and H.E.S.S. results on the spatial extension, we find evidence for an energy-dependent shrinking of the Crab Nebula's $\gamma$-ray extension ($R_{68}\propto E_\mathrm{IC}^{-\alpha}$ where $\alpha=0.155\pm0.035_{stat}{-0.037}_{sys}$).

PSR J1119-6127 is a radio pulsar that behaved magnetar like bursts, and we performed a comprehensive investigation of this pulsar using the archival high-energy observations obtained after its outburst in 2016 July. After the 2016 outburst, specific regions on the neutron star surface were heated up to > 0.3 and > 1 keV from about 0.2 keV. A hard non-thermal spectral component with a photon index < 0.5 related to the magnetospheric emission can be resolved from the NuSTAR spectra above 10 keV. We find that the thermal emitting regions did not cool down and gradually shrank by about (20-35)% four months after the outburst. Hard X-ray pulsations were detected with NuSTAR immediately after the outburst, at 5 sigma confidence level and with a background-subtracted pulsed fraction of (40+-10)%. However, the signal became undetectable after a few days. Using Fermi data, we found that the gamma-ray emission in 0.5-300 GeV was suppressed along with the disappearance of the radio pulsations. This is likely caused by a reconfiguration of the magnetic field. We also discovered that the timing noise evolved dramatically, and the spin-down rate significantly increased after the 2016 glitch. We proposed that post-outburst temporal and spectral behaviors from radio to gamma-ray bands were caused by changes of the magnetosphere structure, pair plasma injection and the shrinking emission sites on the neutron star.

Supernova remnant (SNR)\,W28 is well known for its classic hadronic scenario, in which the TeV cosmic rays (CRs) released at the early stage of this mid-aged SNR are illuminating nearby molecular clouds (MCs). Overwhelming evidences have shown that the northeast of the SNR (W28-North) has already encountered with the MC clumps. Through this broken shell -- W28-North, we believe the CRs with energy down to $<$1\,GeV to be able to be injected into nearby MCs. To further testify this hadronic scenario, we first analyse the 9 years Fermi-LAT data in/around W28 with energy down to 0.3\,GeV. Our Fermi-LAT analysis display a 10-200 GeV skymap which spatially matches well with the known TeV sources -- HESS\,J1801-233 (W28-North), HESS\,J1800-240\,A,\,B\,\&\,C (240\,A\,B\,\&\,C). At low energy band, we has discovered a 0.5-1\,GeV blob located to the south of 240\,B\,\&\,C, and a low flux of 0.3-1\,GeV at 240\,A. A hadronic model is build to explain our analysis results and previous multi-wavelength observations of W28. Our model consists of three CR sources: The run-away CRs escaped from a strong shock; The leaked GeV CRs from the broken shell -- W28-North; And the local CR sea. Through modelling the SNR evolution, CR acceleration\,\&\,releasing, we have explained the GeV-TeV emission in/around SNR\,W28 (except for 240\,A) in one model. Both the damping of the magnetic waves by the neutrals and the decreased acceleration efficiency are taken into account in our model due to the mid-age of SNR W28.

While the Kes 73/1E 1841-045 system had been confirmed as an extended GeV source, whether its morphology depends on the photon energy or not deserves our further investigation. Adopting data collected by Fermi Large Area Telescope (LAT) again, we look into the extensions of this source in three energy bands individually: 0.3-1 GeV, 1-3 GeV and 3-200 GeV. We find that the 0.3-1 GeV morphology is point-like and is quite different from those in the other two bands, although we cannot robustly reject a unified morphology for the whole LAT band.

The supernova remnant Kes 73 and/or the magnetar 1E 1841-045 at its center can deposit a large amount of energy to the surroundings and is potentially responsible for particle acceleration. Using the data taken with the Fermi Large Area Telescope (LAT), we confirmed the presence of an extended source whose centroid position is highly consistent with this magnetar/supernova-remnant pair. Its emission is intense from 100 MeV to $>$100 GeV. Its LAT spectrum can be decoupled into two components which are respectively governed by two different mechanisms. According to the young age of this system, the magnetar is seemingly a necessary and sufficient source for the downward-curved spectrum below 10 GeV, as the observed $<$10 GeV flux is too high for the supernova remnant to account for. On the other hand, the supernova remnant is reasonably responsible for the hard spectrum above 10 GeV. Further studies of this region in the TeV regime is required, so that we can perform physically meaningful comparisons of the $>$10 GeV spectrum and the TeV spectrum.

We present observations at 7 mm that fully resolve the two circumstellar disks, and a reanalyses of archival observations at 3.5 cm that resolve along their major axes the two ionized jets, of the class I binary protostellar system L1551 NE. We show that the two circumstellar disks are better fit by a shallow inner and steep outer power-law than a truncated power-law. The two disks have very different transition radii between their inner and outer regions of $\sim$18.6 AU and $\sim$8.9 AU respectively. Assuming that they are intrinsically circular and geometrically thin, we find that the two circumstellar disks are parallel with each other and orthogonal in projection to their respective ionized jets. Furthermore, the two disks are closely aligned if not parallel with their circumbinary disk. Over an interval of $\sim$10 yr, source B (possessing the circumsecondary disk) has moved northwards with respect to and likely away from source A, indicating an orbital motion in the same direction as the rotational motion of their circumbinary disk. All the aforementioned elements therefore share the same axis for their angular momentum, indicating that L1551 NE is a product of rotationally-driven fragmentation of its parental core. Assuming a circular orbit, the relative disk sizes are compatible with theoretical predictions for tidal truncation by a binary system having a mass ratio of $\sim$0.2, in agreement with the reported relative separations of the two protostars from the center of their circumbinary disk. The transition radii of both disks, however, are a factor of $\gtrsim$1.5 smaller than their predicted tidally-truncated radii.

Either bulk rotation or local turbulence is widely invoked to drive fragmentation in collapsing cores so as to produce multiple star systems. Even when the two mechanisms predict different manners in which the stellar spins and orbits are aligned, subsequent internal or external interactions can drive multiple systems towards or away from alignment thus masking their formation process. Here, we demonstrate that the geometrical and dynamical relationship between the binary system and its surrounding bulk envelope provide the crucial distinction between fragmentation models. We find that the circumstellar disks of the binary protostellar system L1551 IRS 5 are closely parallel not just with each other but also with their surrounding flattened envelope. Measurements of the relative proper motion of the binary components spanning nearly 30 yr indicate an orbital motion in the same sense as the envelope rotation. Eliminating orbital solutions whereby the circumstellar disks would be tidally truncated to sizes smaller than are observed, the remaining solutions favor a circular or low-eccentricity orbit tilted by up to $\sim$25$^\circ$ from the circumstellar disks. Turbulence-driven fragmentation can generate local angular momentum to produce a coplanar binary system, but which bears no particular relationship with its surrounding envelope. Instead, the observed properties conform with predictions for rotationally-driven fragmentation. If the fragments were produced at different heights or on opposite sides of the midplane in the flattened central region of a rotating core, the resulting protostars would then exhibit circumstellar disks parallel with the surrounding envelope but tilted from the orbital plane as is observed.

The region around SGR 1806-20 and its host stellar cluster Cl* 1806-20 is a potentially important site of particle acceleration. The soft $\gamma-$ray repeater and Cl* 1806-20, which also contains several very massive stars including a luminous blue variable hypergiant LBV 1806-20, are capable of depositing a large amount of energy to the surroundings. Using the data taken with the \emph{Fermi} Large Area Telescope (LAT), we identified an extended LAT source to the south-west of Cl* 1806-20. The centroid of the 1-50~GeV emission is consistent with that of HESS J1808-204 (until now unidentified). The LAT spectrum is best-fit by a broken power-law with the break energy $E_\mathrm{b}=297\pm15$ MeV. The index above $E_\mathrm{b}$ is $2.60\pm0.04$, and is consistent with the flux and spectral index above 100 GeV for HESS J1808-204, suggesting an association between the two sources. Meanwhile, the interacting supernova remnant SNR G9.7-0.0 is also a potential contributor to the LAT flux. A tentative flux enhancement at the MeV band during a 45-day interval (2011 Jan 21 - 2011 Mar 7) is also reported. We discuss possible origins of the extended LAT source in the context of both leptonic and hadronic scenarios.

We report the results from a detailed $\gamma-$ray investigation in the field of two "dark accelerators", HESS J1745-303 and HESS J1741-302, with $6.9$ years of data obtained by the Fermi Large Area Telescope. For HESS J1745-303, we found that its MeV-GeV emission is mainly originated from the "Region A" of the TeV feature. Its $\gamma-$ray spectrum can be modeled with a single power-law with a photon index of $\Gamma\sim2.5$ from few hundreds MeV to TeV. Moreover, an elongated feature, which extends from "Region A" toward northwest for $\sim1.3^{\circ}$, is discovered for the first time. The orientation of this feature is similar to that of a large scale atomic/molecular gas distribution. For HESS J1741-302, our analysis does not yield any MeV-GeV counterpart for this unidentified TeV source. On the other hand, we have detected a new point source, Fermi J1740.1-3013, serendipitously. Its spectrum is apparently curved which resembles that of a $\gamma-$ray pulsar. This makes it possibly associated with PSR B1737-20 or PSR J1739-3023.

We explore very fine scales of scalar dissipation in turbulent mixing, below Kolmogorov and around Batchelor scales, by performing direct numerical simulations at much finer grid resolution than is usually adopted in the past. We consider the resolution in terms of a local, fluctuating Batchelor scale and study the effects on the tails of the probability density function and multifractal properties of the scalar dissipation. The origin and importance of these very fine-scale fluctuations are discussed. One conclusion is that they are unlikely to be related to the most intense dissipation events.