Context. A recent observational census of Kuiper belt objects (KBOs) has unveiled anomalous orbital structures. This has led to the hypothesis that an additional $\sim5-10~m_{\oplus}$ planet exists. This planet, known as Planet 9, occupies an eccentric and inclined orbit at hundreds of astronomical units. However, the KBOs under consideration have the largest known semimajor axes at $a>250$ AU; thus they are very difficult to detect.// Aims. In the context of the proposed Planet 9, we aim to measure the mean plane of the Kuiper belt at $a>50$ AU. In a comparison of the expected and observed mean planes, some constraints would be put on the mass and orbit of this undiscovered planet. Methods. We adopted and developed the theoretical approach of Volk & Malhotra (2017) to the relative angle $\delta$ between the expected mean plane of the Kuiper belt and the invariable plane determined by the eight known planets. Numerical simulations were constructed to validate our theoretical approach. Then similar to Volk & Malhotra (2017), we derived the angle $\delta$ for the real observed KBOs with $100<a<200$ AU, and the measurement uncertainties were also estimated. Finally, for comparison, maps of the theoretically expected $\delta$ were created for different combinations of possible Planet 9 parameters. Results. The expected mean plane of the Kuiper belt nearly coincides with the said invariable plane interior to $a=90$ AU. But these two planes deviate noticeably from each other at $a>100$ AU owing to the presence of Planet 9 because the relative angle $\delta$ could be as large as $\sim10^{\circ}$. Using the $1\sigma$ upper limit of $\delta<5^{\circ}$ deduced from real KBO samples as a constraint, we present the most probable parameters of Planet 9.

The conversion between photons and axionlike particles (ALPs) in the Milky Way magnetic field could result in the detectable oscillation phenomena in $\gamma$-ray spectra of Galactic sources. In this work, the GeV (Fermi-LAT) and TeV (MAGIC/VERITAS/H.E.S.S.) data of three bright supernova remnants (SNRs, ie. IC443, W51C and W49B) have been adopted together to search such the oscillation effect. Different from our previous analysis of the sole Fermi-LAT data of IC443, we do not find any reliable signal for the photon-ALP oscillation in the joint broadband spectrum of each SNR. The reason for the inconsistence is that in this work we use the latest revision (P8R3) of Fermi-LAT data, updated diffuse emission templates and the new version of the source catalog (4FGL), which lead to some modification of the GeV spectrum of IC443. Then we set constraints on ALP parameters based on the combined analysis of all the three sources. Though these constraints are somewhat weaker than limits from the CAST experiment and globular clusters, they are supportive of and complementary to these other results.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

We aim to determine the relative angle between the total angular momentum of the minor planets and that of the Sun-planets system, and to improve the orientation of the invariable plane of the solar system. By utilizing physical parameters available in public domain archives, we assigned reasonable masses to 718041 minor planets throughout the solar system, including near-Earth objects, main belt asteroids, Jupiter trojans, trans-Neptunian objects, scattered-disk objects, and centaurs. Then we combined the orbital data to calibrate the angular momenta of these small bodies, and evaluated the specific contribution of the massive dwarf planets. The effects of uncertainties on the mass determination and the observational incompleteness were also estimated. We determine the total angular momentum of the known minor planets to be $1.7817\times10^{46}$ g $\cdot$ cm$^2$ $\cdot$ s$^{-1}$. The relative angle $\alpha$ between this vector and the total angular momentum of the Sun-planets system is calculated to be 14.74 deg.. By excluding the dwarf planets Eris, Pluto, and Haumea, which have peculiar angular momentum directions, $\alpha$ drops sharply to 1.76 deg.; a similar result applies to each individual minor planet group. This suggests that, without these three most massive bodies, the plane perpendicular to the total angular momentum of the minor planets would be close to the invariable plane of the solar system. On the other hand, the inclusion of Eris, Haumea, and Makemake can produce a difference of 1254 mas in the inclination of the invariable plane, which is much larger than the difference of 9 mas induced by Ceres, Vesta, and Pallas as found previously. By taking into account the angular momentum contributions from all minor planets, including unseen ones, the orientation improvement of the invariable plane is larger than 1000 mas in inclination with a $1\sigma$ error of 50-140 mas.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

In order to search for the line-like signals in the Fermi-LAT data, we have analyzed totally 49152 regions of interest (ROIs) that cover the whole sky. No ROI displays a line signal with test statistic (TS) value above 25, while for 50 ROIs weak line-like excesses with ${\rm TS}>16$ are presented. The intrinsic significances of these potential signals are further reduced by the large trial factor introduced in such kind of analysis. For the largest TS value of 24.3 derived in our analysis, the corresponding global significance is only $0.54\sigma$. We thus do not find any significant line-like signal and set up constraints on the cross section of dark matter annihilating to gamma-ray lines, $\left<\sigma{v}\right>_{\gamma\gamma}$.

In this work, we search for $\gamma$-ray emission from the directions of some nearby Milky Way dwarf spheroidal galaxies (dSphs) and candidates with the publicly-available Pass 8 data of Fermi-LAT. Our sample includes 12 sources with the distances $<50$ kpc. Very weak $\gamma$-ray excesses ($\sim 2\sigma$) are found in some dSphs/candidates, consistent with that reported in the previous literature. Intriguingly, the peak test statistic (TS) value of the weak emission from the direction of Reticulum II rises continually. If interpreted as dark matter (DM) annihilation, the peak TS value is 13.5 for the annihilation channel of $\chi\chi \rightarrow \tau^{+}\tau^{-}$ and the DM mass of $m_\chi \sim 16$ GeV. The combination of all these nearby sources yields a more significant (with local significance $> 4\sigma$) $\gamma$-ray signal.

Axion-like particles (ALPs) can oscillate to photons and vice versa in electromagnetic fields. The photon-ALP oscillation provides an attractive solution to the apparent transparency of the Universe to TeV photons. The allowed parameter regions for the ALP mass $m_{\rm a}\leq 10^{-7}$ eV have been tightly constrained by the Fermi-LAT and H.E.S.S observations of some extragalactic sources. In this work we show for the first time that the H.E.S.S observations of some TeV sources in the Galactic plane exclude the highest ALP mass region (i.e., $m_{\rm a}\sim {\rm a ~few\times 10^{-7}}$ eV) that accounts for the TeV transparency of the Universe. The upcoming CTA observations of the Galactic TeV sources are shown to be able to improve the constraints significantly.

Axionlike-particles (ALPs) are one promising type of dark matter candidate particle that may generate detectable effects on $\gamma$-ray spectra other than the canonical weakly interacting massive particles. In this work we search for such oscillation effects in the spectra of supernova remnants caused by the photon-ALP conversion, using the Fermi Large Area Telescope data. Three bright supernova remnants, IC443, W44, and W51C, are analyzed. The inclusion of photon-ALP oscillations yields an improved fit to the $\gamma$-ray spectrum of IC443, which gives a statistical significance of $4.2\sigma$ in favor of such spectral oscillation. However, the best-fit parameters of ALPs ($m_{a}=6.6\,{\rm neV}$, $g_{a\gamma}=13.4 \times 10^{-11}\,{\rm GeV}^{-1}$) are in tension with the upper bound ($g_{a\gamma}< 6.6 \times 10^{-11}\,{\rm GeV}^{-1}$) set by the CAST experiment. It is difficult to explain the results using the systematic uncertainties of the flux measurements. We speculate that the "irregularity" displayed in the spectrum of IC443 may be due to the superposition of the emission from different parts of the remnant.

Recently, gamma-ray emission in the direction of Coma, with a TS value of $\sim 40$, has been reported. In this work we will discuss the possibility of such a residual emission coming from dark matter annihilation. Our results show that the gamma-ray emission within the Coma region is also spatially correlated to the mass distribution derived from weak gravitational lensing measurements very well. However the dark matter models are not supported by the spectral analysis results and constraints by observations of other targets. Thus we derive the upper limits of the dark matter annihilation cross section according to the observation of the Coma region.

High-redshift blazars are one of the most powerful monsters in the universe and $\gamma$-ray variability carries crucial information of their relativistic jets. In this work we present results of the first systematical temporal analysis of {\it Fermi}-LAT data of all known seven $\gamma$-ray blazars beyond redshift 3. Significant long-term $\gamma$-ray variability is found from five sources in monthly $\gamma$-ray light curves, in which three of them are reported for the first time. Furthermore, intraday $\gamma$-ray variations are detected from NVSS J053954$-$283956 and NVSS J080518$+$614423. Doubling variability timescale of the former source is limited as short as $\lesssim$ 1 hour (at the source frame). Together with variability amplitude over one order of magnitude, NVSS J053954$-$283956 is the most distant $\gamma$-ray flaring blazar so far. Meanwhile, intraday optical variability of NVSS J163547$+$362930 is found based on archival PTF/iPTF light curve. Benefited from multiwavelength activity of these sources, constraints of their Doppler factors as well as locations of $\gamma$-ray radiation region and indications for the SDSS high redshift jetted active galactic nuclei deficit are discussed.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

Milky Way-like galaxies are predicted to host a very large number of dark matter subhalos. Some massive and nearby subhalos could generate detectable gamma-rays, appearing as unidentified, spatially-extended and stable gamma-ray sources. We search for such sources in the third Fermi Large Area Telescope source List (3FGL) and report the identification of a new candidate, 3FGL J1924.8-1034. With the Fermi-LAT Pass 8 data, we find that 3FGL J1924.8-1034 is spatially-extended at a high confidence level of $5.4\sigma$, with a best-fit extension radius of $\sim0.15^{\circ}$. No significant variability has been found and its gamma-ray spectrum is well fitted by the dark matter annihilation into $b\bar{b}$ with a mass of $\sim 43$ GeV. All these facts make 3FGL J1924.8-1034 a possible dark matter subhalo candidate. However, due to the limited angular resolution, the possibility of that the spatial extension of 3FGL J1924.8-1034 is caused by the contamination from the other un-resolved point source can not be ruled out.

In this work, we simulate a set of realizations of local volume dark matter subhalo population based on the distributions and relations derived from Via Lactea II N-body simulation. We calculate the J-factors of these subhalos, and find that the low mass subhalos contribute a lot to the total J-factors. Combining with 91 months of the Fermi LAT observation, we constrain on the cross section of dark matter annihilating directly to two gamma rays. This is the first work combining numerical simulation results and Fermi LAT observations to constrain dark matter cross section to gamma-ray line with subhalo population. Though the constraints derived from subhalo population are weaker than those from Fermi LAT observation of the Galactic center, they are supports of and complementary to these other results.

With 91 months of the publicly available Fermi-LAT Pass 8 data, we analyze the gamma-ray emission from the Milky Way satellites to search for potential line signals due to the annihilation of dark matter particles into double photons. The searched targets include a sample of dwarf spheroidal galaxies, the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC). No significant line emission has been found neither in the stacked dwarf galaxy sample nor in the direction of LMC/SMC. The corresponding upper limits on the cross section of DM annihilation into two photons are derived. Compared with results of previous gamma-ray line searches with the Pass 7 data, the current constraints on the line emission from dwarf spheroidal galaxies has been significantly improved in a wide energy range. With the rapid increase of the sample of dwarf spheroidal galaxies (candidates), we expect that the sensitivity of gamma ray line searches will be significantly improved in the near future.

This paper is to explore a model of the ABS Algorithms dealing with the solution of a class of systems of linear stochastic equations $A\xi=\eta$ when $\eta$ is a $m$-dimensional normal distribution. It is shown that the stepsize $\alpha_i$ is distributed as $N(u_i,\sigma_i)$ (being $u_i$ the expected value of $\alpha_i$ and $\sigma_i$ its variance) and the approximation to the solutions $\xi_{i}$ is distributed as $N_n(U_i,\Sigma_i)$ (being $U_i$ the expected value of $\xi_i$ and $\Sigma_i$ its variance), for this algorithm model.

We consider the energy release associated with first-order transition by Gibbs construction and present such energy release as an accumulation of a series of tiny binding energy differences between over-compressed states and stable ones. Universal formulae for the energy release from one homogeneous phase to the other is given. We find the energy release per converted particle varies with number density. As an example, the deconfinement phase transition at supranuclear densities is discussed in detail. The mean energy release per converted baryon is of order 0.1MeV in RMF theory and MIT bag descriptions for hadronic matter and strange quark matter for a wider parameter region.

We investigated the thermal evolution of rotating strange stars with the deconfinement heating due to magnetic braking. We consider the stars consisting of either normal quark matter or color-flavor-locked phase. Combining deconfinement heating with magnetic field decay, we find that the thermal evolution curves are identical to pulsar data.

We investigated the non-leptonic weak interaction in magnetic field. We discussed an improvement of previous method to analytical work out the rate for weak field case.Our result easily goes over to field-free limit.Then we calculated the reaction rate in strong magnetic field where the charged particles are confined to the lowest Landau level. A strong magnetic field strongly suppressed the rate,which will be foreseen to affect viscous dynamics in SQM .We also derived a few approximation formulae under given conditions that can be conveniently applied.

We study the non-linear evolution of a dust ellipsoid,embedded in a Friedmann flat background universe, in order to determine the evolution of the density of the ellipsoid as the perturbation to it related detaches from general expansion and begins to collapse. We show that while the growth rate of the density contrast of a mass element is enhanced by the shear in agreement with Hoffman (1986a), the angular momentum acquired by the ellipsoid has the right magnitude to counterbalance the effect of the shear. The result confirms the previrialization conjecture (Peebles & Groth 1976; Davis & Peebles 1977; Peebles 1990) by showing that initial asphericities and tidal interactions begin to slow down the collapse after the system has broken away from the general expansion.

This paper is to explore a model of the ABS Algorithms for dealing with a class of systems of linear stochastic equations A xi=eta satisfying eta sim N_m(v, I_{m}). It is shown that the iteration step alpha_{i} is N(V,\pi) and approximation solutions is xi_{i} \sim N_n(U,\Sigma) for this algorithm model. And some properties of (V,\pi)$ and $(U,\Sigma) are given.

We present a review and bibliography of the main results obtained during a research on ABS (Abaffy, Broyden, Spedicato) methods.

We present the main results obtained during a research on ABS methods in the framework of the project Analisi Numerica e Matematica Computazionale.