Search

In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes of astrophysical $\gamma$-ray background while large amount of dark matter. By analyzing more than 700 days observational data at LHAASO, no significant dark matter signal from 1 TeV to 1 EeV is detected. Accordingly we derive the most stringent constraints on the ultra-heavy dark matter annihilation cross-section up to EeV. The constraints on the lifetime of dark matter in decay mode are also derived., Comment: 17 pages, 12 figures, accepted by PRL

The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To ensure the reliability of the LHAASO-KM2A data, a three-level quality control system has been established. It is used to monitor the status of detector units, stability of reconstructed parameters and the performance of the array based on observations of the Crab Nebula and Moon shadow. This paper will introduce the control system and its application on the LHAASO-KM2A data collected from August 2021 to July 2023. During this period, the pointing and angular resolution of the array were stable. From the observations of the Moon shadow and Crab Nebula, the results achieved using the two methods are consistent with each other. According to the observation of the Crab Nebula at energies from 25 TeV to 100 TeV, the time averaged pointing errors are estimated to be $-0.003^{\circ} \pm 0.005^{\circ}$ and $0.001^{\circ} \pm 0.006^{\circ}$ in the R.A. and Dec directions, respectively., Comment: 15 pages, 9 figures

The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) is compatible with NGC 4278 within $\sim0.03$ degree. Variation analysis shows an indication of the variability at a few months level in the TeV band, which is consistent with low frequency observations. Based on these observations, we report the detection of TeV $\gamma$-ray emissions from this low-luminosity AGN NGC 4278. The observations by LHAASO-WCDA during active period has a significance level of 8.8\,$\sigma$ with best-fit photon spectral index $\varGamma=2.56\pm0.14$ and a flux $f_{1-10\,\rm{TeV}}=(7.0\pm1.1_{\rm{sta}}\pm0.35_{\rm{syst}})\times10^{-13}\,\rm{photons\,cm^{-2}\,s^{-1}}$, or approximately $5\%$ of the Crab Nebula. The discovery of VHE from NGC 4278 indicates that the compact, weak radio jet can efficiently accelerate particles and emit TeV photons., Comment: 11 pages, 5 figures

KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with large altitude difference (30 m) and huge coverage (1.3 km^2). In this paper, the design of the KM2A simulation code G4KM2A based on Geant4 is introduced. The process of G4KM2A is optimized mainly in memory consumption to avoid memory overffow. Some simpliffcations are used to signiffcantly speed up the execution of G4KM2A. The running time is reduced by at least 30 times compared to full detector simulation. The particle distributions and the core/angle resolution comparison between simulation and experimental data of the full KM2A array are also presented, which show good agreement.

Double neutron star (DNS) systems offer excellent opportunities to test gravity theories. We report the timing results of PSR J1901+0658, the first pulsar discovered in the FAST Galactic Plane Pulsar Snapshot (GPPS) Survey. Based on timing observations by FAST over 5 yr, we obtain the phase-coherent timing solutions and derive the precise measurements of its position, spin parameters, orbital parameters, and dispersion measure. It has a period of 75.7 ms, a period derivative of 2.169(6)$\times 10^{-19}$ s s$^{-1}$, and a characteristic age of 5.5 Gyr. This pulsar is in an orbit with a period of 14.45 d and an eccentricity of 0.366. One post-Keplerian parameter, periastron advance, has been well-measured as being 0.00531(9) deg yr$^{-1}$, from which the total mass of this system is derived to be 2.79(7) M$_{\odot}$. The pulsar has the mass upper limit of 1.68 M$_{\odot}$, so the lower limit for the companion mass is 1.11 M$_{\odot}$. Because PSR J1901+0658 is a partially recycled pulsar in an eccentric binary orbit with such a large companion mass, it should be in a DNS system according to the evolution history of the binary system., Comment: 6 pages, 6 figures

We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at $3.67 \pm 0.05 \pm 0.15$ PeV. Below the knee, the spectral index is found to be -$2.7413 \pm 0.0004 \pm 0.0050$, while above the knee, it is -$3.128 \pm 0.005 \pm 0.027$, with the sharpness of the transition measured with a statistical error of 2%. The mean logarithmic mass of cosmic rays is almost heavier than helium in the whole measured energy range. It decreases from 1.7 at 0.3 PeV to 1.3 at 3 PeV, representing a 24% decline following a power law with an index of -$0.1200 \pm 0.0003 \pm 0.0341$. This is equivalent to an increase in abundance of light components. Above the knee, the mean logarithmic mass exhibits a power law trend towards heavier components, which is reversal to the behavior observed in the all-particle energy spectrum. Additionally, the knee position and the change in power-law index are approximately the same. These findings suggest that the knee observed in the all-particle spectrum corresponds to the knee of the light component, rather than the medium-heavy components., Comment: 8 pages, 3 figures

In the era of large language models, Mixture-of-Experts (MoE) is a promising architecture for managing computational costs when scaling up model parameters. However, conventional MoE architectures like GShard, which activate the top-$K$ out of $N$ experts, face challenges in ensuring expert specialization, i.e. each expert acquires non-overlapping and focused knowledge. In response, we propose the DeepSeekMoE architecture towards ultimate expert specialization. It involves two principal strategies: (1) finely segmenting the experts into $mN$ ones and activating $mK$ from them, allowing for a more flexible combination of activated experts; (2) isolating $K_s$ experts as shared ones, aiming at capturing common knowledge and mitigating redundancy in routed experts. Starting from a modest scale with 2B parameters, we demonstrate that DeepSeekMoE 2B achieves comparable performance with GShard 2.9B, which has 1.5 times the expert parameters and computation. In addition, DeepSeekMoE 2B nearly approaches the performance of its dense counterpart with the same number of total parameters, which set the upper bound of MoE models. Subsequently, we scale up DeepSeekMoE to 16B parameters and show that it achieves comparable performance with LLaMA2 7B, with only about 40% of computations. Further, our preliminary efforts to scale up DeepSeekMoE to 145B parameters consistently validate its substantial advantages over the GShard architecture, and show its performance comparable with DeepSeek 67B, using only 28.5% (maybe even 18.2%) of computations.

The rapid development of open-source large language models (LLMs) has been truly remarkable. However, the scaling law described in previous literature presents varying conclusions, which casts a dark cloud over scaling LLMs. We delve into the study of scaling laws and present our distinctive findings that facilitate scaling of large scale models in two commonly used open-source configurations, 7B and 67B. Guided by the scaling laws, we introduce DeepSeek LLM, a project dedicated to advancing open-source language models with a long-term perspective. To support the pre-training phase, we have developed a dataset that currently consists of 2 trillion tokens and is continuously expanding. We further conduct supervised fine-tuning (SFT) and Direct Preference Optimization (DPO) on DeepSeek LLM Base models, resulting in the creation of DeepSeek Chat models. Our evaluation results demonstrate that DeepSeek LLM 67B surpasses LLaMA-2 70B on various benchmarks, particularly in the domains of code, mathematics, and reasoning. Furthermore, open-ended evaluations reveal that DeepSeek LLM 67B Chat exhibits superior performance compared to GPT-3.5.

In this paper, we present an innovative process-oriented math process reward model called \textbf{Math-Shepherd}, which assigns a reward score to each step of math problem solutions. The training of Math-Shepherd is achieved using automatically constructed process-wise supervision data, breaking the bottleneck of heavy reliance on manual annotation in existing work. We explore the effectiveness of Math-Shepherd in two scenarios: 1) \textit{Verification}: Math-Shepherd is utilized for reranking multiple outputs generated by Large Language Models (LLMs); 2) \textit{Reinforcement Learning}: Math-Shepherd is employed to reinforce LLMs with step-by-step Proximal Policy Optimization (PPO). With Math-Shepherd, a series of open-source LLMs demonstrates exceptional performance. For instance, the step-by-step PPO with Math-Shepherd significantly improves the accuracy of Mistral-7B (77.9\%$\to$84.1\% on GSM8K and 28.6\%$\to$33.0\% on MATH). The accuracy can be further enhanced to 89.1\% and 43.5\% on GSM8K and MATH with the verification of Math-Shepherd, respectively. We believe that automatic process supervision holds significant potential for the future evolution of LLMs., Comment: Add Step-by-Step reinforcement learning results

For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE; $E_\gamma \geq 100$~TeV) $\gamma$-rays. In this context, the historical SNR Cassiopeia A (Cas A) is considered one of the most promising target for UHE observations. This paper presents the observation of Cas A and its vicinity by the LHAASO KM2A detector. The exceptional sensitivity of LHAASO KM2A in the UHE band, combined with the young age of Cas A, enabled us to derive stringent model-independent limits on the energy budget of UHE protons and nuclei accelerated by Cas A at any epoch after the explosion. The results challenge the prevailing paradigm that Cas A-type SNRs are major suppliers of PeV CRs in the Milky Way., Comment: 11 pages, 3 figures, Accepted by the APJL

The highest energy gamma-rays from gamma-ray bursts (GRBs) have important implications for their radiation mechanism. Here we report for the first time the detection of gamma-rays up to 13 TeV from the brightest GRB 221009A by the Large High Altitude Air-shower Observatory (LHAASO). The LHAASO-KM2A detector registered more than 140 gamma-rays with energies above 3 TeV during 230$-$900s after the trigger. The intrinsic energy spectrum of gamma-rays can be described by a power-law after correcting for extragalactic background light (EBL) absorption. Such a hard spectrum challenges the synchrotron self-Compton (SSC) scenario of relativistic electrons for the afterglow emission above several TeV. Observations of gamma-rays up to 13 TeV from a source with a measured redshift of z=0.151 hints more transparency in intergalactic space than previously expected. Alternatively, one may invoke new physics such as Lorentz Invariance Violation (LIV) or an axion origin of very high energy (VHE) signals., Comment: 49pages, 11figures

The axion is well motivated in physics. It solves the strong charge conjugation-parity reversal problem CP in fundamental physics and the dark matter problem in astronomy. Its interaction with the electromagnetic field has been expected but never detected experimentally. Such particles may convert to radio waves in the environment with a strong magnetic field. Inspired by the idea, various research groups have been working on theoretical modeling and radio data analysis to search for the signature of radio signals generated by the axion conversion in the magnetosphere of compact stars, where the surface magnetic field as strong as $10^{13}$-$10^{14}$ G is expected. In this work, we calculate the observational properties of the axion-induced radio signals (AIRSs) in the neutron star magnetosphere, where both the total intensity and polarization properties of radio emission are derived. Based on the ray tracing method, assuming 100% linear polarization of radio waves generated in each conversion, we compute the polarization emission profile concerning different viewing angles. We note that plasma and general relativistic effects are important for the polarization properties of AIRSs. Our work suggests that AIRSs can be identified by the narrow bandwidth and distinct polarization features., Comment: 15 pages, 7 figures. Published in Physical Review D

We report five new fast radio bursts (FRBs) discovered from the Galactic Plane Pulsar Snapshot (GPPS) survey by the Five-hundred-meter Aperture Spherical radio Telescope (FAST): FRB\,20210126, FRB\,20210208, FRB\,20210705, FRB\,20211005 and FRB\,20220306. To date, no repeating bursts from these FRB sources have been detected in the follow-up monitoring observations, leading to their classification as potential one-off events. We obtain the basic parameters for these bursts, including position, dispersion measure (DM), pulse width, spectral index, scattering time-scale, etc. The fluences and flux densities are generally lower in comparison to the values observed in one-off bursts discovered by other telescopes. Among the observed bursts, polarization data for 4 bursts were recorded during observations. Consequently, we obtain polarization profiles and Faraday rotation measures (RMs) for these bursts., Comment: 7 pages, 5 figures, 2 tables. Published in MNRAS, Volume 526, Issue 2, December 2023

Pulsars radiate radio signals when they rotate. However, some old pulsars often stop radiating for some periods. The underlying mechanism remains unknown, while the magnetosphere during nulling phases is hard to probe due to the absence of emission measurement. Here we report the detection and accurate polarization measurements of sporadic weak narrow dwarf pulses detected in the ordinary nulling state of pulsar B2111+46 via the Five-Hundred-Meter Aperture Spherical radio Telescope (FAST). Further analysis shows that their polarization angles follow the average polarization angle curve of normal pulses, suggesting no change of magnetic field structure in the emission region in the two emission states. Whereas radio emission of normal individual pulses is radiated by a thunderstorm of particles produced by copious discharges in regularly formed gaps, dwarf pulses are produced by one or a few raindrops of particles generated by pair production in a fragile gap of this near-death pulsar., Comment: 37 pages, 15 figures. Authors version. Great FAST pulsar data. Share full-text access: https://rdcu.be/djLSL

We present the first catalog of very-high energy and ultra-high energy gamma-ray sources detected by the Large High Altitude Air Shower Observatory (LHAASO). The catalog was compiled using 508 days of data collected by the Water Cherenkov Detector Array (WCDA) from March 2021 to September 2022 and 933 days of data recorded by the Kilometer Squared Array (KM2A) from January 2020 to September 2022. This catalog represents the main result from the most sensitive large coverage gamma-ray survey of the sky above 1 TeV, covering declination from $-$20$^{\circ}$ to 80$^{\circ}$. In total, the catalog contains 90 sources with an extended size smaller than $2^\circ$ and a significance of detection at $> 5\sigma$. Based on our source association criteria, 32 new TeV sources are proposed in this study. Among the 90 sources, 43 sources are detected with ultra-high energy ($E > 100$ TeV) emission at $> 4\sigma$ significance level. We provide the position, extension, and spectral characteristics of all the sources in this catalog., Comment: 40 pages, 13 figures, 4 tables

Timing observations are crucial for determining the basic parameters of newly discovered pulsars. Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) with the L-band 19-beam receiver covering the frequency range of 1.0--1.5 GHz, the FAST Galactic Plane Pulsar Snapshot (GPPS) Survey has discovered more than 600 faint pulsars with flux densities of only a few or a few tens of $\mu$Jy at 1.25 GHz. To obtain accurate position, spin parameters and dispersion measure of a pulsar, and to calculate derived parameters such as the characteristic age and surface magnetic field, we collect available FAST pulsar data obtained either through targeted follow-up observations or through coincidental survey observations with one of the 19 beams of the receiver. From these data we obtain time of arrival (TOA) measurements for 30 newly discovered pulsars as well as for 13 known pulsars. We demonstrate that the TOA measurements acquired by the FAST from any beams of the receiver in any observation mode (e.g. the tracking mode or the snapshot mode) can be combined to get timing solutions. We update the ephemerides of 13 previously known pulsars and obtain the first phase-coherent timing results for 30 isolated pulsars discovered in the FAST GPPS Survey. Notably, PSR J1904+0853 is an isolated millisecond pulsar, PSR J1906+0757 is a disrupted recycled pulsar, and PSR J1856+0211 has a long period of 9.89 s that can constrain pulsar death lines. Based on these timing solutions, all available FAST data have been added together to obtain the best pulse profiles for these pulsars.

The diffuse Galactic $\gamma$-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this work we report the measurements of diffuse $\gamma$-rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer array of the Large High Altitude Air Shower Observatory (LHAASO). Diffuse emissions from the inner ($15^{\circ}10$~TeV). The energy spectrum in the inner Galaxy regions can be described by a power-law function with an index of $-2.99\pm0.04$, which is different from the curved spectrum as expected from hadronic interactions between locally measured cosmic rays and the line-of-sight integrated gas content. Furthermore, the measured flux is higher by a factor of $\sim3$ than the prediction. A similar spectrum with an index of $-2.99\pm0.07$ is found in the outer Galaxy region, and the absolute flux for $10\lesssim E\lesssim60$ TeV is again higher than the prediction for hadronic cosmic ray interactions. The latitude distributions of the diffuse emission are consistent with the gas distribution, while the longitude distributions show clear deviation from the gas distribution. The LHAASO measurements imply that either additional emission sources exist or cosmic ray intensities have spatial variations., Comment: 12 pages, 8 figures, 5 tables; accepted for publication in Physical Review Letters; source mask file provided as ancillary file

Strangeon star model has passed various observational tests, such as the massive pulsars and the tidal deformability during binary mergers. Pulsar glitch, as a useful probe for studying the interior structure of pulsars, has also been studied in strangeon star model in our previous papers, including the recovery coefficient, the waiting time of glitches and glitch activity. In this paper, the recovery process of a glitch is described in the strangeon star model, based on the starquake picture established in Paper I. After the starquake, the inner motion of the stellar matter would reduce the tangential pressure in the cracked places at the equatorial plane. The recovery (increase) of the tangential pressure would be achieved by a viscous flow towards the cracked places at equatorial plane, which leads to the exponential recovery of the spin frequency. A uniform viscous flow can reproduce the single exponential decay observed in some glitches, and the viscous time-scale $\tau$ and the depth $h$ of the cracking place below the surface can be fitted by the recovery data. It is found that $h$ increases with glitch size $\Delta\nu/\nu$, which is expected in the glitch scenario of strangeon stars. The magnitude of the recovery predicted in this recovery model is also consistent with that derived from observations. The single exponential decay reproduced by a uniform viscous flow can be generalized to two or more exponentials by the multi-component of viscous flows., Comment: 7 pages, 4 figures, 1 table. Accepted for publication in MNRAS

More than 100 rotating radio transients (RRATs) have been discovered since 2006. However, it is unclear whether RRATs radiate in the nulling states. PSR J0628+0909 has been classified as an RRAT. In this paper, we study the single pulses and integrated pulse profile of PSR J0628+0909 to check whether we can detect pulsed radio emission in the nulling states. We also aim to study the polarization of the RRAT and its relationship to the general pulsar population. We used the Five-hundred-meter Aperture Spherical radio Telescope (FAST) to observe PSR J0628+0909 in the frequency range from 1.0 to 1.5 GHz. We searched for strong single pulses and looked for pulsed emission in the RRAT nulling states. Polarisation profiles, the single-pulse energy distribution, and waiting-time statistics were measured. The Faraday rotation measure and dispersion measure values are updated with the current observation. The single-pulse polarisation behaviours show great diversity, similar to the case of pulsars. Based on the integrated pulse profile and single-pulse energy statistics, we argue that continuous pulsar-like emission exists in addition to the transient-like burst emission for PSR J0628+0909. We find that the pulse waiting time is not correlated with the pulse energy and conclude that the strong transient emission of RRAT is not generated by the energy store-release mechanism., Comment: 10 pages, 9 figures, accepted by MNRAS

The Large High Altitude Air Shower Observatory (LHAASO) has three sub-arrays, KM2A, WCDA and WFCTA. The flux variations of cosmic ray air showers were studied by analyzing the KM2A data during the thunderstorm on 10 June 2021. The number of shower events that meet the trigger conditions increases significantly in atmospheric electric fields, with maximum fractional increase of 20%. The variations of trigger rates (increases or decreases) are found to be strongly dependent on the primary zenith angle. The flux of secondary particles increases significantly, following a similar trend with that of the shower events. To better understand the observed behavior, Monte Carlo simulations are performed with CORSIKA and G4KM2A (a code based on GEANT4). We find that the experimental data (in saturated negative fields) are in good agreement with simulations, assuming the presence of a uniform upward electric field of 700 V/cm with a thickness of 1500 m in the atmosphere above the observation level. Due to the acceleration/deceleration and deflection by the atmospheric electric field, the number of secondary particles with energy above the detector threshold is modified, resulting in the changes in shower detection rate., Comment: 18 pages, 11 figures

Fast radio bursts (FRBs) are highly dispersed millisecond-duration radio bursts. Recent observations of a Galactic FRB suggest that at least some FRBs originate from magnetars, but the origin of cosmological FRBs is still not settled. Here we report the detection of 1863 bursts in 82 hr over 54 days from the repeating source FRB~20201124A. These observations show irregular short-time variation of the Faraday rotation measure (RM), which probes the density-weighted line-of-sight magnetic field strength, of individual bursts during the first 36 days, followed by a constant RM. We detected circular polarisation in more than half of the burst sample, including one burst reaching a high fractional circular polarisation of 75%. Oscillations in fractional linear and circular polarisations as well as polarisation angle as a function of wavelength were detected. All of these features provide evidence for a complicated, dynamically evolving, magnetised immediate environment within about an astronomical unit (au; Earth-Sun distance) of the source. Our optical observations of its Milky-Way-sized, metal-rich host galaxy reveal a barred spiral, with the FRB source residing in a low stellar density, interarm region at an intermediate galactocentric distance. This environment is inconsistent with a young magnetar engine formed during an extreme explosion of a massive star that resulted in a long gamma-ray burst or superluminous supernova., Comment: 28 pages, 10 figures, updated to match the published version

The Crab pulsar and the surrounding nebula powered by the pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind is a bright source of gamma-rays carrying crucial information about this complex conglomerate. We report the detection of $\gamma$-rays with a spectrum showing gradual steepening over three energy decades, from $5\times 10^{-4}$ to $1.1$ petaelectronvolt (PeV). The ultra-high-energy photons exhibit the presence of a PeV electron accelerator (a pevatron) with an acceleration rate exceeding 15% of the absolute theoretical limit. Assuming that unpulsed $\gamma$-rays are produced at the termination of the pulsar's wind, we constrain the pevatron's size, between $0.025$ and $0.1$ pc, and the magnetic field $\approx 110 \mu$G. The production rate of PeV electrons, $2.5 \times 10^{36}$ erg $\rm s^{-1}$, constitutes 0.5% of the pulsar's spin-down luminosity, although we do not exclude a non-negligible contribution of PeV protons to the production of the highest energy $\gamma$-rays., Comment: 43 pages, 13 figures, 2 tables; Published in Science

The Wide Field-of-View Cherenkov Telescope Array (WFCTA) and the Water Cherenkov Detector Arrays (WCDA) of LHAASO are designed to work in combination for measuring the energy spectra of various cosmic ray species over a very wide energy range from a few TeV to 10 PeV. The energy calibration of WCDA can be achieved with a proven technique of measuring the westward shift of the Moon shadow of galactic cosmic rays due to the geomagnetic field. This deflection angle $\Delta$ is inversely proportional to the energy of the cosmic rays. The precise measurements of the shifts by WCDA allows us to calibrate its energy scale for energies as high as 35 TeV. The energy scale measured by WCDA can be used to cross calibrate the energy reconstructed by WFCTA, which spans the whole energy range up to 10 PeV. In this work, we will demonstrate the feasibility of the method using the data collected from April 2019 to January 2020 by the WFCTA array and WCDA-1 detector, the first of the three water Cherenkov ponds, already commissioned at LHAASO site.

The focal plane camera is the core component of the Wide Field-of-view Cherenkov/fluorescence Telescope Array (WFCTA) of the Large High-Altitude Air Shower Observatory (LHAASO). Because of the capability of working under moonlight without aging, silicon photomultipliers (SiPM) have been proven to be not only an alternative but also an improvement to conventional photomultiplier tubes (PMT) in this application. Eighteen SiPM-based cameras with square light funnels have been built for WFCTA. The telescopes have collected more than 100 million cosmic ray events and preliminary results indicate that these cameras are capable of working under moonlight. The characteristics of the light funnels and SiPMs pose challenges (e.g. dynamic range, dark count rate, assembly techniques). In this paper, we present the design features, manufacturing techniques and performances of these cameras. Finally, the test facilities, the test methods and results of SiPMs in the cameras are reported here., Comment: 45 pages, 21 figures, article

Fast radio bursts (FRBs) are millisecond-duration radio transients of unknown origin. Two possible mechanisms that could generate extremely coherent emission from FRBs invoke neutron star magnetospheres or relativistic shocks far from the central energy source. Detailed polarization observations may help us to understand the emission mechanism. However, the available FRB polarization data have been perplexing, because they show a host of polarimetric properties, including either a constant polarization angle during each burst for some repeaters or variable polarization angles in some other apparently one-off events. Here we report observations of 15 bursts from FRB 180301 and find various polarization angle swings in seven of them. The diversity of the polarization angle features of these bursts is consistent with a magnetospheric origin of the radio emission, and disfavours the radiation models invoking relativistic shocks., Comment: Published online in Nature on 29 Oct, 2020

As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the first observation on the Crab Nebula, a standard candle in very high energy gamma-ray astronomy. We detect gamma-ray signals from the Crab Nebula in both energy ranges of 10$-$100 TeV and $>$100 TeV with high significance, by analyzing the KM2A data of 136 live days between December 2019 and May 2020. With the observations, we test the detector performance including angular resolution, pointing accuracy and cosmic ray background rejection power. The energy spectrum of the Crab Nebula in the energy range 10-250 TeV fits well with a single power-law function dN/dE =(1.13$\pm$0.05$_{stat}$$\pm$0.08$_{sys}$)$\times$10$^{-14}$$\cdot$(E/20TeV)$^{-3.09\pm0.06_{stat}\pm0.02_{sys}}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$. It is consistent with previous measurements by other experiments. This opens a new window of gamma-ray astronomy above 0.1 PeV through which ultrahigh-energy gamma-ray new phenomena, such as cosmic PeVatrons, might be discovered., Comment: 13 pages, 15 figures,submitted to CPC

The state of supranuclear matter in compact stars remains puzzling, and it is argued that pulsars could be strangeon stars. The consequences of merging double strangeon stars are worth exploring, especially in the new era of multi-messenger astronomy. To develop the "strangeon kilonova" scenario proposed in Paper I, we make a qualitative description about the evolution of ejecta and light curves for merging double strangeon stars. In the hot environment of the merger, the strangeon nuggets ejected by tidal disruption and hydrodynamical squeezing would suffer from evaporation, in which process particles, such as strangeons, neutrons and protons, are emitted. Taking into account both the evaporation of strangeon nuggets and the decay of strangeons, most of the strangeon nuggets would turn into neutrons and protons, within dozens of milliseconds after being ejected. The evaporation rates of different particles depend on temperature, and we find that the ejecta could end up with two components, with high and low opacity respectively. The high opacity component would be in the directions around the equatorial plane, and the low opacity component would be in a broad range of angular directions. The bolometric light curves show that even if the total ejected mass would be as low as $\sim 10^{-4} M_\odot$, the spin-down power of the long-lived remnant would account for the whole emission of kilonova AT2017gfo associated with GW 170817. The detailed picture of merging double strangeon stars is expected to be tested by future numerical simulations.

Fast radio bursts (FRBs) are mysterious millisecond-duration radio transients of unknown origin observed at extragalactic distances. It has been long speculated that magnetars are the engine powering repeating bursts from FRB sources, but no convincing evidence has been collected so far\cite{sun19}. Recently, the Galactic magnetar SGR J1935+2154 entered an active phase by emitting intense soft Gamma-ray bursts. One FRB-like event with two peaks (FRB 200428) and a luminosity slightly lower than the faintest extragalactic FRBs was detected from the source, in association with a soft Gamma-ray / hard X-ray flare. Here we report an eight-hour targeted radio observational campaign comprising four sessions and assisted by multi-wavelength (optical and hard X-rays) data. During the third session, 29 soft Gamma-ray repeater (SGR) bursts were detected in Gamma-ray energies. Throughout the observing period, we detected no single dispersed pulsed emission coincident with the arrivals of SGR bursts, but unfortunately we were not observing when the FRB was detected. The non-detection places a fluence upper limit that is eight orders of magnitude lower than the fluence of FRB 200428. Our results suggest that FRB -- SGR burst associations are rare. FRBs may be highly relativistic and geometrically beamed, or FRB-like events associated with SGR bursts may have narrow spectra and characteristic frequencies outside the observed band. It is also possible that the physical conditions required to achieve coherent radiation in SGR bursts are difficult to satisfy, and that only under extreme conditions could an FRB be associated with an SGR burst., Comment: 27 pages, 7 figures, Published in Nature 2020, Vol 587

We present timing solutions and analyses of 11 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These pulsars were discovered using an ultra-wide bandwidth receiver in drift-scan observations made during the commissioning phase of FAST, and were then confirmed and timed using the 64-m Parkes Radio Telescope. Each pulsar has been observed over a span of at least one year. Highlighted discoveries include PSR J0344-0901, which displays mode-changing behaviour and may belong to the class of so-called `swooshing' pulsars (alongside PSRs B0919+06 and B1859+07); PSR J0803-0942, whose emission is almost completely linearly polarised; and PSRs J1900-0134 and J1945+1211, whose well defined polarisation angle curves place stringent constraints on their emission geometry. We further discuss the detectability of these pulsars by earlier surveys, and highlight lessons learned from our work in carrying out confirmation and monitoring observations of pulsars discovered by a highly sensitive telescope, many of which may be applicable to next-generation pulsar surveys. This paper marks one of the first major releases of FAST-discovered pulsars, and paves the way for future discoveries anticipated from the Commensal Radio Astronomy FAST Survey (CRAFTS)., Comment: 18 pages, 8 figures, 5 tables. Accepted for publication by MNRAS. Correcting error in the caption of Figure 8

We present our piggyback search for fast radio bursts using the Nanshan 26m Radio Telescope and the Kunming 40m Radio Telescope. The observations are performed in the L-band from 1380 MHz to 1700 MHz at Nanshan and S-band from 2170 MHz to 2310 MHz at Kunming. We built the \textsc{Roach2}-based FFT spectrometer and developed the real-time transient search software. We introduce a new radio interference mitigation technique named \emph{zero-DM matched filter} and give the formula of the signal-to-noise ratio loss in the transient search. Though we have no positive detection of bursts in about 1600 and 2400 hours data at Nanshan and Kunming respectively, an intriguing peryton was detected at Nanshan, from which hundreds of bursts were recorded. Perytons are terrestrial radio signals that mimic celestial fast radio bursts. They were first reported at Parkes and identified as microwave oven interferences later. The bursts detected at Nanshan show similar frequency swept emission and have double-peaked profiles. They appeared in different sky regions in about tens of minutes observations and the dispersion measure index is not exactly 2, which indicates the terrestrial origin. The peryton differs drastically from the known perytons detected at Parkes, because it appeared in a precise period of $p=1.71287\pm 0.00004$ s. Its origin remains unknown., Comment: 17 pages, 13 figures

We apply theoretical spin-down models of magnetospheric evolution and magnetic field decay to simulate the possible evolution of PSR J0250+5854, which is the slowest-spinning radio pulsar detected to date. Considering the alignment of inclination angle in a 3-D magnetosphere, it is possible that PSR J0250+5854 has a high magnetic field comparable with magnetars or/and high magnetic field pulsars, if a small inclination angle is considered. Our calculations show that similar long-period pulsars tend to have a relatively low period derivative in this case. In another case of magnetic field decay, calculations also show a possible connection between PSR J0250+5854 and high dipole-magnetic field magnetars. The evolutionary path indicates a relatively high spin-down rate for similar long-period pulsars., Comment: 8 pages, 4 figures

The state of supranuclear matter in compact star remains puzzling, and it is argued that pulsars could be strangeon stars. What if binary strangeon stars merge? This kind of merger could result in the formation of a hyper-massive strangeon star, accompanied by bursts of gravitational waves and electromagnetic radiation (and even strangeon kilonova explained in the paper). The tidal polarizability of binary strangeon stars is different from that of binary neutron stars, because a strangeon star is self-bound on surface by fundamental strong force while a neutron star by the gravity, and their equations of state are different. Our calculation shows that the tidal polarizability of merging binary strangeon stars is favored by GW170817. Three kinds of kilonovae (i.e., of neutron, quark and strangeon) are discussed, and the light curve of the kilonova AT 2017gfo following GW170817 could be explained by considering the decaying strangeon nuggets and remnant star spin-down. Additionally, the energy ejected to the fireball around the nascent remnant strangeon star, being manifested as a Gamma-ray burst (GRB), is calculated. It is found that, after a promote burst, an X-ray plateau could follow in a timescale of $10^{2-3}$ s. Certainly, the results could be tested also by further observational synergies between gravitational wave detectors (e.g., aLIGO) and X-ray telescopes (e.g., Chinese HXMT and eXTP), and especially if the detected gravitational wave form is checked by peculiar equation of state provided by the numerical relativistical simulation., Comment: Research in Astronomy and Astrophysics, 18, 24 (2018)

Pulsar-like compact stars provide us a unique laboratory to explore properties of dense matter at supra-nuclear densities. One of the models for pulsar-like stars is that they are totally composed of "strangeons", and in this paper we studied the pulsar glitches in a strangeon star model. Strangeon stars would be solidified during cooling, and the solid stars would be natural to have glitches as the result of starquakes. Based on the starquake model established before, we proposed that when the starquake occurs, the inner motion of the star which changes the moment of inertia and has impact on the glitch sizes, is divided into plastic flow and elastic motion. The plastic flow which is induced in the fractured part of the outer layer, would move tangentially to redistribute the matter of the star and would be hard to recover. The elastic motion, on the other hand, changes its shape and would recover significantly. Under this scenario, we could understand the behaviors of glitches without significant energy releasing, including the Crab and the Vela pulsars, in an uniform model. We derive the recovery coefficient as a function of glitch size, as well as the time interval between two successive glitches as the function of the released stress. Our results show consistency with observational data under reasonable ranges of parameters. The implications on the oblateness of the Crab and the Vela pulsars are discussed., Comment: MNRAS, accepted

We report identification of the optical counterpart to the companion of the millisecond pulsar J2317+1439. At the timing position of the pulsar, we find an object with $g=22.96\pm0.05$, $r=22.86\pm0.04$ and $i=22.82\pm0.05$. The magnitudes and colors of the object are consistent with it being a white dwarf. By comparing with white dwarf cooling models, we estimate that it has a mass of $0.39^{+0.13}_{-0.10}$ M$_{\odot}$, an effective temperature of $8077^{+550}_{-470}$ K and a cooling age of $10.9\pm0.3$ Gyr. Combining our results with published constraints on the orbital parameters obtained through pulsar timing, we estimate the pulsar mass to be $3.4^{+1.4}_{-1.1}$ M$_{\odot}$. Although the constraint on the pulsar mass is still weak, there is a significant possibility that the pulsar could be more massive than two solar mass., Comment: 7 pages, 6 figures, accepted for publication in ApJ

The emission geometries, e.g. the emission region height, the beam shape, and radius-to-frequency mapping, are important predictions of pulsar radiation model. The multi-band radio observations carry such valuable information. In this paper, we study two bright pulsars, (PSRs B0329+54 and B1642-03) and observe them in high frequency (2.5 GHz, 5 GHz, and 8 GHz). The newly acquired data together with historical archive provide an atlas of multi-frequency profiles spanning from 100 MHz to 10 GHz. We study the frequency evolution of pulse profiles and the radiation regions with the these data. We firstly fit the pulse profiles with Gaussian functions to determine the phase of each component, and then calculate the radiation altitudes of different emission components and the radiation regions. We find that the inverse Compton scattering (ICS) model can reproduce the radiation geometry of these two pulsars. But for PSR B0329+54 the radiation can be generated in either annular gap (AG) or core gap (CG), while the radiation of PSR B1642-03 can only be generated in the CG. This difference is caused by the inclination angle and the impact angle of these two pulsars. The relation of beaming angle (the angle between the radiation direction and the magnetic axis) and the radiation altitudes versus frequency is also presented by modelling the beam-frequency evolution in the ICS model. The multi-band pulse profiles of these two pulsars can be described well by the ICS model combined with the CG and AG.

Magnetars may have strong surface dipole field. Observationally, two magnetars may have passive fallback disks. In the presence of a fallback disk, the rotational evolution of magnetars may be changed. In the self-similar fallback disk model, it is found that: (1) When the disk mass is significantly smaller than $10^{-6} \,\rm M_{\odot}$, the magnetar is unaffected by the fallback disk and it will be a normal magnetar. (2) When the disk mass is large, but the magnetar's surface dipole field is about or below $10^{14} \,\rm G$, the magnetar will also be a normal magnetar. A magnetar plus a passive fallback disk system is expected. This may correspond to the observations of magnetars 4U 0142$+$61, and 1E 2259$+$586. (3) When the disk mass is large, and the magnetar's surface dipole field is as high as $4\times 10^{15} \,\rm G$, the magnetar will evolve from the ejector phase to the propeller phase, and then enter rotational equilibrium. The magnetar will be slowed down quickly in the propeller phase. The final rotational period can be as high $2\times 10^4 \,\rm s$. This may correspond to the super-slow magnetar in the supernova remnant RCW 103. Therefore, the three kinds of magnetars can be understood in a unified way., Comment: 14 pages, 4 figures, accepted for ApJ

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025., Comment: 16 pages, 16 figures. Oral talk presented at SPIE Astronomical Telescopes and Instrumentation, June 26 to July 1, 2016, Edingurgh, UK

Detailed structures of the Vela pulsar (PSR B0833-45, with a period of $89.33$ milliseconds) are predicted by adopting a recently-constructed unified treatment of all parts of neutron stars: the outer crust, the inner crust and the core based on modern microscopic Brueckner-Hartree-Fock calculations. To take the pulsar mass ranging from $1.0M_{\odot}$ to $2.0M_{\odot}$, we calculate the central density, the core/crust radii, the core/crust mass, the core/crustal thickness, the moment of inertia, and the crustal moment of inertia. Among them, the crustal moment of inertia could be effectively constrained from the accumulated glitch observations, which has been a great debate recently, known as "glitch crisis". Namely, superfluid neutrons contained in the inner crust, which are regarded as the origin of the glitch in the standard two-component model, could be largely entrained in the nuclei lattices, then there may not be enough superfluid neutrons ($\sim 4/5$ less than the previous value) to trigger the large glitches ($\Delta \nu/\nu_0 \sim 10^{-6}$) in the Vela pulsar. We then provide the first analysis of the crisis based on the microscopic basis, by confronting the glitch observations with the theoretical calculations for the crustal moment of inertia. We find that despite some recent opposition to the crisis argument, the glitch crisis is still present, which means that besides the crust superfluid neutrons, core neutrons might be necessary for explaining the large glitches of the Vela pulsar., Comment: 28 pages, 6 figures, 4 tables

Pulsar glitches, i.e. the sudden spin-ups of pulsars, have been detected for most pulsars that we known. The mechanism giving rise to this kind of phenomenon is uncertain, although a large data set has been built. In the framework of star-quake model, based on~\cite{Baym1971}, the glitch-sizes (the relative increases of spin-frequencies during glitches) $\Delta \Omega/\Omega$ depend on the released energies during glitches, with less released energies corresponding to smaller glitch sizes. On the other hand, as one of dark matter candidates, our Galaxy might be filled with the so called strange nuggets (SNs) which are the relics from the early Universe. In this case the collisions between pulsars and SNs is inevitable, and these collisions would lead to glitches when enough elastic energy has been accumulated during the spin-down process. The SNs-triggered glitches could release less energy, because the accumulated elastic energy would be less than that in the scenario of glitches without SNs. Therefore, if a pulsar is hit frequently by SNs, it would tend to have more small size glitches, whose values of $\Delta \Omega/\Omega$ are smaller than that in standard star-quake model (with larger amounts of released energy). Based on the assumption that in our Galaxy the distribution of SNs is similar to that of dark matter, as well as on the glitch data in ATNF Pulsar Catalogue and Jodrell Bank glitch table, we find that in our Galaxy the incidences of small size glitches exhibit tendencies consistent with the collision rates of pulsars and strange nuggets. Further test of this scenario is expected by detecting more small glitches (e.g., by the Square Kilometre Array)., Comment: 11 pages, 1 figure, 1 table. Accepted by RAA

We present high signal-to-noise ratio, multi-frequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array (PPTA) project. The pulsars are observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20 cm multibeam receiver. Observations spanning approximately six years have been carefully calibrated and summed to produce high S/N profiles. This allows us to study the individual profile components and in particular how they evolve with frequency. We also identify previously undetected profile features. For many pulsars we show that pulsed emission extends across almost the entire pulse profile. The pulse component widths and component separations follow a complex evolution with frequency; in some cases these parameters increase and in other cases they decrease with increasing frequency. The evolution with frequency of the polarization properties of the profile is also non-trivial. We provide evidence that the pre- and post-cursors generally have higher fractional linear polarization than the main pulse. We have obtained the spectral index and rotation measure for each pulsar by fitting across all three observing bands. For the majority of pulsars, the spectra follow a single power-law and the position angles follow a $\lambda^2$ relation, as expected. However, clear deviations are seen for some pulsars. We also present phase-resolved measurements of the spectral index, fractional linear polarization and rotation measure. All these properties are shown to vary systematically over the pulse profile., Comment: 41 pages, 28 figures, accepted for publication in MNRAS

Pulsar-like compact stars usually have strong magnetic fields, with the strength from $\sim 10^8$ to $\sim 10^{12}$ Gauss on surface. How such strong magnetic fields can be generated and maintained is still an unsolved problem, which is, in principle, related to the interior structure of compact stars, i.e., the equation of state of cold matter at supra-nuclear density. In this paper we are trying to solve the problem in the regime of solid quark-cluster stars. Inside quark-cluster stars, the extremely low ratio of number density of electrons to that of baryons $n_e/n_b$ and the screening effect from quark-clusters could reduce the long-range Coulomb interaction between electrons to short-range interaction. In this case, the Stoner's model could apply, and we find that the condition for ferromagnetism is consistent with that for validity of Stoner's model. Under the screened Coulomb repulsion, the electrons inside the stars could spontaneously magnetized and become ferromagnetic, and hence would contribute non-zero net magnetic momentum to the whole star. We conclude that, for most cases in solid quark-cluster stars, the amount of net magnetic momentum, which is proportional to the amount of unbalanced spins $\xi=(n_+-n_-)/n_e$ and depends on the number density of electrons $n_e=n_++n_-$, could be significant with non-zero $\xi$. The net magnetic moments of electron system in solid quark-cluster stars could be large enough to induce the observed magnetic fields for pulsars with $B\sim 10^{11}$ to $\sim 10^{13}$ Gauss., Comment: 7 pages, 1 figure. Accepted by Chinese Physics C