Search

We present detailed timing and spectral analyses of the transient X-ray pulsar RX J0209.6$-$7427 in the Small Magellanic Cloud during its 2019 giant outburst. With a better known distance than most galactic X-ray pulsars, its peak luminosity is determined to be $(1.11\pm0.06)\times 10^{39}\, \rm erg\ s^{-1}$; it is thus a {\it bonda fide} pulsating ultraluminous X-ray source (PULX). Owing to the broad energy band of \textit{Insight}-HXMT, its pulsed X-ray emission was detected from 1 keV up to the 130$-$180 keV band, which is the highest energy emission detected from any PULXs outside the Milky Way. This allows us to conclude that its main pulsed X-ray emission is from the "fan beam" of the accretion column, and its luminosity is thus intrinsic. We also estimate its magnetic field of (4.8$-$8.6)$\times10^{12}$ G or (1.7$-$2.2)$\times10^{13}$ G, from its spin evolution or transition in the accretion column structure during the outburst; we suggest that the two values of the magnetic field strength correspond to the dipole and multipole magnetic fields of the neutron star, similar to the recent discovery in the Galactic PULX Swift J0243.6+6124. Therefore, the nature of the neutron star and its ULX emission can be understood within the current theoretical frame of accreting neutron stars. This may have implications for understanding the nature of those farther away extragalactic PULXs., Comment: 31 pages, 14 figures, 5 tables. Accepted by ApJ

The fast transitions between different types of quasi-periodic oscillations (QPOs) are generally observed in black hole transient sources (BHTs). We present a detailed study on the timing and spectral properties of the transitions of type-B QPOs in MAXI~J1348--630, observed by \emph{Insight}-HXMT. The fractional rms variability--energy relationship and energy spectra reveal that type-B QPOs probably originate from jet precession. Compared to weak power-law dominated power spectrum, when type-B QPO is present, the corresponding energy spectrum shows an increase in Comptonization component and the need for {\tt\string xillverCp} component, and a slight increase of height of the corona when using {\tt\string relxilllp} model. Therefore, we suggest that a coupled inner disk-jet region is responsible for the observed type-B QPOs transitions. The time scale for the appearance/disappearance of type-B QPOs is either long or short (seconds), which may indicate an instability of disk-jet structure. For these phenomena, we give the hypothesis that the Bardeen-Petterson effect causes disk-jet structure to align with BH spin axis, or that the disappearance of small-scale jets bound by the magnetic flux tubes lead to the disappearance of type-B QPOs. We observed three events regarding the B/C transitions, one of which occurred in a short time from $\sim 9.2$ Hz (C) to $\sim 4.8$ Hz (B). The energy spectral analysis for the other two transitions shows that when type-C QPO is present, the Comptonization flux is higher, the spectrum is harder and the inner radius of disk changes insignificantly. We suggest that type-C QPOs probably originate from relatively stronger jets or corona., Comment: 26 pages, 9 figures

We studied the transient Galactic black hole candidate MAXI J0637-430 with data from Insight-HXMT, Swift and XMM-Newton. The broad-band X-ray observations from Insight-HXMT help us constrain the power-law component. MAXI J0637-430 is located at unusually high Galactic latitude; if it belongs to the Galactic thick disk, we suggest a most likely distance <7 kpc. Compared with other black hole transients, MAXI J0637-430 is also unusual for other reasons: a fast transition to the thermal dominant state at the start of the outburst; a low peak temperature and luminosity (we estimate them at ~ 0.7 keV and <0.1 times Eddington, respectively); a short decline timescale; a low soft-to-hard transition luminosity (<0.01 times Eddington). We argue that such properties are consistent with a small binary separation, short binary period (P ~ 2 hr), and low-mass donor star (M2 ~ 0.2 M_sun). Moreover, spectral modelling shows that a single disk-blackbody component is not a good fit to the thermal emission. Soft spectral residuals, and deviations from the standard L_disk ~ T^4 in relation, suggest the need for a second thermal component. We propose and discuss various scenarios for such component, in addition to those presented in previous studies of this source. For example, a gap in the accretion disk between a hotter inner ring near the innermost stable orbit, and a cooler outer disk. Another possibility is that the second thermal component is the thermal plasma emission from an ionized outflow.

We present a detailed study of the high mass X-ray binary Vela X-1, using observations performed by Insight-HXMT in 2019 and 2020, concentrating on timing analysis and spectral studies including pulse phase-resolved spectroscopy. The cyclotron line energy is found to be at ~21-27 keV and 43-50 keV for the fundamental and first harmonic, respectively. We present the evolution of spectral parameters and find that two line centroid energy ratio E2/E1 evolved from ~2 before MJD 58900 to ~1.7 after that. The harmonic cyclotron line energy has no relation to the luminosity but the fundamental line energy shows a positive correlation with X-ray luminosity, suggesting that Vela X-1 is located in the sub-critical accreting regime. In addition, the pulse phase-resolved spectroscopy in Vela X-1 is performed. Both the CRSF and continuum parameters show strong variability over the pulse phase with the ratio of two line energies about 2 near the peak phases, and down to ~1.6 around off-peak phases. Long-term significant variations of the absorption column density and its evolution over the pulse phase may imply the existence of the clumpy wind structure near the neutron star., Comment: 10 pages, MNRAS in press

In this paper, wavelet analysis is used to study spectral-timing properties of MAXI J1535-571 observed by Insight-HXMT. The low-frequency quasi-periodic oscillations (QPOs) are detected in nine observations. Based on wavelet analysis, the time intervals with QPO and non-QPO are isolated separately, and the corresponding spectra with QPO and non-QPO are analyzed. We find that the spectra with QPO (hereafter QPO spectra) are softer than those without QPO (hereafter non-QPO spectra) in the hard intermediate state (HIMS). While in the soft intermediate state (SIMS), the QPO spectra are slightly harder. The disk temperature of QPO regime is slightly lower during HIMS, but becomes higher during SIMS. The cutoff energies of QPO spectra and non-QPO spectra do not show significant differences. The flux ratio of the disk to total flux is higher for the time intervals with non-QPO than that of QPO regime. We propose that these differences in the spectral properties between QPO and non-QPO regimes could be explained in the scenario of Lense-Thirring precession, and the reversal of the QPO/non-QPO behavior between HIMS and SIMS may be associated with appearance/disappearance of a type-B QPO which might origin from the precession of the jet., Comment: 13 pages, 9 figures in main part, and 3 figures in the appendix; MNRAS in press, the first submission in 2021-06

Multiple isotopes of samarium element are the isotopes produced by the s process, and 154Sm is produced by the r process. In addition, 144Sm is p nuclei in nuclear astrophysics. The measurement of these can help us to better understand the results of relevant photonuclear reaction experiments. On the other hand, 149Sm is a 235U fission product with a 1% yield, its cross sections are important to reactor neutronics. In this work, the neutron capture yield of the natural samarium target was measured at the back-streaming white neutron beamline (Back-n) of the China Spallation Neutron Source (CSNS), and the resonance parameters were analyzed by SAMMY code. The resonance peaks and the neutron separation energies contributed by the different isotopes are considered individually. The results of the capture yield found signs of the possibility of two resonance peaks at 8 eV, which awaits further experimental examination. Cross-section was calculated according to resonance parameters and was compared with other experimental results and evaluation databases of ENDF/B-VIII.0 and CENDL-3.2. A clear difference between ENDF/B VIII.0 and CENDL-3.2 database appears at 23.2 eV, the experimental result at this energy is smaller than data of ENDF/B VIII.0 database but CENDL-3.2 database. Most of the controversial experimental results invariably come from the samarium 149 isotope., Comment: 10 pages,7 figures

Here we present a detailed study of the broadband noise in the power density spectra of the black hole X-ray binary MAXI J1820+070 during the hard state of its 2018 outburst, using the Hard X-ray Modulation Telescope (Insight-HXMT) observations. The broadband noise shows two main humps, which might separately correspond to variability from a variable disk and two Comptonization regions. We fitted the two humps with multiple Lorentzian functions and studied the energy-dependent properties of each component up to 100--150 keV and their evolution with spectral changes. The lowest frequency component is considered as the sub-harmonic of QPO component and shows different energy dependence compared with other broadband noise components. We found that although the fractional rms of all the broadband noise components mainly decrease with energy, their rms spectra are different in shape. Above $\sim$ 20--30 keV, the characteristic frequencies of these components increase sharply with energy, meaning that the high-energy component is more variable on short timescales. Our results suggest that the hot inner flow in MAXI J1820+070 is likely to be inhomogeneous. We propose a geometry with a truncated accretion disk, two Comptonization regions.

We present the catalog of Interplanetary Network (IPN) localizations for 199 short-duration gamma-ray bursts (sGRBs) detected by the Konus-Wind (KW) experiment between 2011 January 1 and 2021 August 31, which extends the initial sample of IPN localized KW sGRBs (arXiv:1301.3740) to 495 events. We present the most comprehensive IPN localization data on these events, including probability sky maps in HEALPix format., Comment: Published in ApJS

Using the observations of the high-energy (HE) detector of the Hard X-ray Modulation Telescope (Insight-HXMT) for Scorpius X-1 in 2018, we search for hard X-ray tails in the hard X-ray spectra in ~30-200 keV. The hard X-ray tails are found throughout the Z-track on the hardness-intensity diagram and they harden and fade away from the horizontal branch (HB), through the normal branch (NB), to the flaring branch (FB). Comparing the hard X-ray spectra between Cyg X-1 and Sco X-1, it is concluded that the hard X-ray spectrum of Cyg X-1 shows high-energy cutoff, implying a hot corona in it, but the high-energy cutoff does not reveal in the hard X-ray spectrum of Sco X-1. Jointly fitting the HE spectrum with the medium-energy and low-energy spectra of Sco X-1 in ~2-200 keV, it is suggested that the upscattering Comptonization of the neutron star (NS) emission photons by the energetic free-falling electrons onto the NS or by the hybrid electrons in the boundary layer between the NS and the accretion disk could be responsible for the hard X-ray tails of Sco X-1 on the HB and NB, but neither of the two mechanisms can be responsible for the hard X-ray tail on the FB. Some possible origins for the peculiar hard X-ray tail of FB are argued., Comment: 26 pages, 7 figures

We present a spectral study of the black hole candidate MAXI J1348-630 during its 2019 outburst, based on monitoring observations with Insight-HXMT and Swift. Throughout the outburst, the spectra are well fitted with power-law plus disk-blackbody components. In the soft-intermediate and soft states, we observed the canonical relation L ~ T_in^4 between disk luminosity L and peak colour temperature T_in, with a constant inner radius R_in (traditionally identified with the innermost stable circular orbit). At other stages of the outburst cycle, the behaviour is more unusual, inconsistent with the canonical outburst evolution of black hole transients. In particular, during the hard rise, the apparent inner radius is smaller than in the soft state (and increasing), and the peak colour temperature is higher (and decreasing). This anomalous behaviour is found even when we model the spectra with self-consistent Comptonization models, which take into account the up-scattering of photons from the disk component into the power-law component. To explain both those anomalous trends at the same time, we suggest that the hardening factor for the inner disk emission was larger than the canonical value of ~1.7 at the beginning of the outburst. A more physical trend of radii and temperature evolution requires a hardening factor evolving from ~3.5 at the beginning of the hard state to ~1.7 in the hard intermediate state. This could be evidence that the inner disk was in the process of condensing from the hot, optically thin medium and had not yet reached a sufficiently high optical depth for its emission spectrum to be described by the standard optically-thick disk solution., Comment: 20 pages, 10 figures, 4 tables, accepted by ApJ

One month after launching Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM), a bright thermonuclear X-ray burst from 4U~0614+09, was observed on January 24, 2021. We report the time-resolved spectroscopy of the burst and a burst oscillation detection at 413 Hz with a fractional amplitude 3.4\% (rms). This coincides with the burst oscillation previously discovered with \textit{Swift}/BAT \citep{Strohmayer2008}, and therefore confirms the spin frequency of this source. This burst is the brightest one in the normal bursts (except the superburst) ever detected from 4U~0614+09, which leads to an upper limit of distance estimation as 3.1 kpc. The folded light curve during the burst oscillation shows a multi-peak structure, which is the first case observed during a single burst oscillation in nonpulsating sources. The multi-peak profile could be due to additional harmonics of the burst oscillation, which is corresponding to several brighter/fainter spots at the stellar surface.

As the main detector of the GECAM satellite, the calibration of the energy response and detection efficiency of the GRD detector is the main content of the ground-based calibration. The calibration goal requires the calibrated energy points to sample the full energy range (8 keV-2 MeV) as much as possible. The low energy band (8-160 keV) is calibrated with the X-ray beam, while the high energy band (>160 keV) with radioactive sources. This article mainly focuses on the calibration of the energy response and detection efficiency in the 8-160 keV with a refined measurement around the absorption edges of the lanthanum bromide crystal. The GRD performances for different crystal types, data acquisition modes, working modes, and incident positions are also analyzed in detail. We show that the calibration campaign is comprehensive, and the calibration results are generally consistent with simulations as expected., Comment: 16 pages,46 figures

The GECAM mission consists of two identical microsatellites (GECAM-A and GECAM-B). Each satellite is equipped with 25 gamma-ray detectors (GRD) and 8 charged particle detectors (CPD). The main scientific objective of the GECAM mission is to detect gamma-ray bursts (GRBs) associated with the gravitational wave events produced by the merging of binary compact stars. After the launch on Dec. 10, 2020 , we carried out a series of on orbit tests. This paper introduces the test results of the GECAM-B satellite. According to the in-flight performance, the energy band for gamma-ray detection of GECAM-B is from about 7 keV to 3.5 MeV. GECAM-B can achieve prompt localization of GRBs. For the first time, GECAM-B realized a quasi-real-time transmission of trigger information using Beidou-3 RDSS. Keywords GECAM, gamma-ray burst, gravitational wave, GRD, CPD, Comment: 16 pages, 8 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

In the standard model of magnetic reconnection, both ions and electrons couple to the newly reconnected magnetic field lines and are ejected away from the reconnection diffusion region in the form of bidirectional burst ion and electron jets. Recent observations propose a new model: electron only magnetic reconnection without ion coupling in electron scale current sheet. Based on the data from Magnetospheric Multiscale (MMS) Mission, we observe a long extension inner electron diffusion region (EDR) at least 40 di away from the X line at the terrestrial Magnetopause, implying that the extension of EDR is much longer than the prediction of the theory and simulations. This inner EDR is embedded in an ion scale current sheet (the width of 4 di, di is ion inertial length). However, such ongoing magnetic reconnection was not accompanied with burst ion outflow, implying the presence of electron only reconnection in ion scale current sheet. Our observations present new challenge for understanding the model of standard magnetic reconnection and electron only reconnection model in electron scale current sheet., Comment: 13 pages, 4 figures, accepted by ApJ

The High Energy X-ray telescope (HE) on-board the Hard X-ray Modulation Telescope (\textit{Insight}-HXMT) can serve as a wide Field of View (FOV) gamma-ray monitor with high time resolution ($\mu$s) and large effective area (up to thousands cm$^2$). We developed a pipeline to search for Gamma-Ray Bursts (GRBs), using the traditional signal-to-noise ratio (SNR) method for blind search and the coherent search method for targeted search. By taking into account the location and spectrum of the burst and the detector response, the targeted coherent search is more powerful to unveil weak and sub-threshold bursts, especially those in temporal coincidence with Gravitational Wave (GW) events. Based on the original method in literature, we further improved the coherent search to filter out false triggers caused by spikes in light curves, which are commonly seen in gamma-ray instruments (e.g. \textit{Fermi}/GBM, \textit{POLAR}). We show that our improved targeted coherent search method could eliminate almost all false triggers caused by spikes. Based on the first two years of \textit{Insight}-HXMT/HE data, our targeted search recovered 40 GRBs, which were detected by either \textit{Swift}/BAT or \textit{Fermi}/GBM but too weak to be found in our blind search. With this coherent search pipeline, the GRB detection sensitivity of \textit{Insight}-HXMT/HE is increased to about 1.5E-08 erg/cm$^2$ (200 keV$-$3 MeV). We also used this targeted coherent method to search \textit{Insight}-HXMT/HE data for electromagnetic (EM) counterparts of LIGO-Virgo GW events (including O2 and O3a runs). However, we did not find any significant burst associated with GW events., Comment: 12 pages, 14 figures, 5 tables; accepted for publication in MNRAS

Cyclotron line scattering features are detected in a few tens of X-ray pulsars (XRPs) and used as direct indicators of a strong magnetic field at the surface of accreting neutron stars (NSs). In a few cases, cyclotron lines are known to be variable with accretion luminosity of XRPs. It is accepted that the observed variations of cyclotron line scattering features are related to variations of geometry and dynamics of accretion flow above the magnetic poles of a NS. A positive correlation between the line centroid energy and luminosity is typical for sub-critical XRPs, where the accretion results in hot spots at the magnetic poles. The negative correlation was proposed to be a specific feature of bright super-critical XRPs, where radiation pressure supports accretion columns above the stellar surface. Cyclotron line in spectra of Be-transient X-ray pulsar GRO J1008-57 is detected at energies from $\sim 75 -90$ keV, the highest observed energy of cyclotron line feature in XRPs. We report the peculiar relation of cyclotron line centroid energies with luminosity in GRO J1008-57 during the Type II outburst in August 2017 observed by Insight-HXMT. The cyclotron line energy was detected to be negatively correlated with the luminosity at $3.2\times 10^{37}\,\ergs

We report the orbital X-ray variability of high mass X-ray binary (HMXB) GX301--2. GX301--2 undergone a spin up process in 2018--2020 with the period evolving from $\sim$ 685 s to 670 s. The energy resolved pulse-profiles of the pulsar in 1--60 keV varied from single peaked and sinusoidal shapes to multi-peaked across different orbital phases. Pulse fractions evolving over orbit had negative correlations with the X-ray flux. The broad-band X-ray energy spectrum of the pulsar can be described with a partial covering negative positive cutoff power-law continuum model. Near the periastron passage of the pulsar we found a strong variation in the additional column density ($NH_{2}$), which correlated with variation of the flux. Curves of growth for both Fe K$\alpha$ and Fe K$\beta$ lines were plotted to investigate the distribution of matter around neutron star. We have also found the evidence of two cyclotron absorption lines in the phase-averaged spectra in GX301--2, with one line of 30--42 keV and the other line varying in 48--56 keV. Both two line's centroid energies show the similar relationship with X-ray luminosity: positive correlation in lower luminosity range, and negative relation above a critical luminosity of $10^{37}$ erg s$^{-1}$. We estimated the surface magnetic field of the neutron star in GX301--2 at ~$(0.5-2)\times 10^{13}$ G. Two cyclotron line energies have a nearly fixed ratio of ~1.63 while having a low strength ratio (~0.05), suggesting that these two features may actually be one line., Comment: 14 pages, 14 figures and 1 table in the main content, and 1 table and 1 figure in the appendix, accept for publication in MNRAS

The X9.3 flare SOL20170906T11:55 was observed by the CsI detector aboard the first Chinese X-ray observatory Hard X-ray Modulation telescope (Insight-HXMT). By using wavelets method, we report about 22 s quasiperiodic pulsations(QPPs) during the impulsive phase. And the spectra from 100 keV to 800 keV showed the evolution with the gamma-ray flux, of a power-law photon index from $\sim 1.8$ before the peak, $\sim 2.0$ around the flare peak, to $\sim 1.8$ again. The gyrosynchrotron microwave spectral analysis reveals a $36.6 \pm 0.6 \arcsec$ radius gyrosynchrotron source with mean transverse magnetic field around 608.2 Gauss, and the penetrated $\ge$ 10 keV non-thermal electron density is about $10^{6.7} \mathrm{cm}^{-3}$ at peak time. The magnetic field strength followed the evolution of high-frequency radio flux. Further gyrosynchrotron source modeling analysis implies that there exists a quite steady gyrosynchrotron source, the non-thermal electron density and transverse magnetic field evolution are similar to higher-frequency light curves. The temporally spectral analysis reveals that those non-thermal electrons are accelerated by repeated magnetic reconnection, likely from a lower corona source., Comment: 14 pages, 9 figures, accept for publication by ApJ

We present a detailed spectral-timing analysis of the Kilohertz quasi-periodic oscillations (kHz QPOs) in Sco X-1 using the data of Rossi X-ray Timing Explorer ($RXTE$) and the Hard X-ray Modulation Telescope ($Insight$-HXMT). The energy band with detectable kHz QPOs is studied for the first time: on the horizontal branch, it is $\sim$ 6.89--24.01 keV and $\sim$ 8.68--21.78 keV for the upper and lower kHz QPOs detected by $RXTE$, and $\sim$ 9--27.5 keV for the upper kHz QPOs by $Insight$-HXMT; on the lower normal branch, the energy band is narrower. The fractional root mean square (rms) of the kHz QPOs increases with energy at lower energy, reaches a plateau at about 16 keV and 20 keV for the lower and upper peaks, and then levels off though with a large uncertainty. The simulation of the deadtime effect of $RXTE$/PCA shows that the deadtime does not affect much the search of the kHz QPOs but makes the rms amplitude underestimated. No significant QPO is detected below $\sim$ 6 keV as shown by the $RXTE$ data, implying that the kHz QPOs do not originate from the black body emission of the accretion disk and neutron star surface. In addition, with the combined analysis of the energy spectra and the absolute rms spectra of kHz QPOs, we suggest that the kHz QPOs in Sco X-1 originate from the Comptonization of the inner part of the transition layer, where the rotation sets the frequency and the inward bulk motion makes the spectrum harder.

A novel framework of intelligent reflecting surface (IRS)-aided multiple-input single-output (MISO) non-orthogonal multiple access (NOMA) network is proposed, where a base station (BS) serves multiple clusters with unfixed number of users in each cluster. The goal is to maximize the sum rate of all users by jointly optimizing the passive beamforming vector at the IRS, decoding order, power allocation coefficient vector and number of clusters, subject to the rate requirements of users. In order to tackle the formulated problem, a three-step approach is proposed. More particularly, a long short-term memory (LSTM) based algorithm is first adopted for predicting the mobility of users. Secondly, a K-means based Gaussian mixture model (K-GMM) algorithm is proposed for user clustering. Thirdly, a deep Q-network (DQN) based algorithm is invoked for jointly determining the phase shift matrix and power allocation policy. Simulation results are provided for demonstrating that the proposed algorithm outperforms the benchmarks, while the throughput gain of 35% can be achieved by invoking NOMA technique instead of orthogonal multiple access (OMA).

GRO J1008-57, as a Be/X-ray transient pulsar, is considered to have the highest magnetic field in known neutron star X-ray binary systems. Observational data of the X-ray outbursts in GRO J1008-57 from 2017 to 2020 were collected by the Insight-HXMT satellite. In this work, the spin period of the neutron star in GRO J1008-57 was determined to be about 93.28 seconds in August 2017, 93.22 seconds in February 2018, 93.25 seconds in June 2019 and 93.14 seconds in June 2020. GRO J1008-57 evolved in the spin-up process with a mean rate of $-(2.10\pm 0.05)\times$10$^{-4}$ s/d from 2009 -- 2018, and turned into a spin down process with a rate of $(6.7\pm 0.6)\times$10$^{-5}$ s/d from Feb 2018 to June 2019. During the type II outburst of 2020, GRO J1008-57 had the spin-up torque again. During the torque reversals, the pulse profiles and continuum X-ray spectra did not change significantly, and the cyclotron resonant scattering feature around 80 keV was only detected during the outbursts in 2017 and 2020. Based on the observed mean spin-up rate, we estimated the inner accretion disk radius in GRO J1008-57 (about 1 - 2 times of the Alfv\'{e}n radius) by comparing different accretion torque models of magnetic neutron stars. During the spin-down process, the magnetic torque should dominate over the matter accreting inflow torque, and we constrained the surface dipole magnetic field $B\geq 6\times 10^{12}$ G for the neutron star in GRO J1008-57, which is consistent with the magnetic field strength obtained by cyclotron line centroid energy., Comment: 10 pages, 5 figures, and 3 tables, the Journal of High Energy Astrophysics in press

In this paper, we presented a detailed timing analysis of a prominent outburst of 4U 0115+63 detected by \textit{Insight}-HXMT in 2017 August. The spin period of the neutron star was determined to be $3.61398\pm 0.00002$ s at MJD 57978. We measured the period variability and extract the orbital elements of the binary system. The angle of periastron evolved with a rate of $0.048\pm0.003$ $yr^{-1}$. The light curves are folded to sketch the pulse profiles in different energy ranges. A multi-peak structure in 1-10 keV is clearly illustrated. We introduced wavelet analysis into our data analysis procedures to study QPO signals and perform a detailed wavelet analysis in many different energy ranges. Through the wavelet spectra, we report the discovery of a QPO at the frequency $\sim 10$ mHz. In addition, the X-ray light curves showed multiple QPOs in the period of $\sim 16-32 $ s and $\sim 67- 200 $ s. We found that the $\sim100$ s QPO was significant in most of the observations and energies. There exist positive relations between X-ray luminosity and their Q-factors and S-factors, while the QPO periods have no correlation with X-ray luminosity. In wavelet phase maps, we found that the pulse phase of $\sim 67- 200 $ s QPO drifting frequently while the $\sim 16-32 $ s QPO scarcely drifting. The dissipation of oscillations from high energy to low energy was also observed. These features of QPOs in 4U 0115+63 provide new challenge to our understanding of their physical origins., Comment: 14 pages, 9 figures, and 6 tables. This work has been submitted to MNRAS after the referee's report

Based on Insight-HXMT data, we report on the pulse fraction evolution during the 2017-2018 outburst of the newly discovered first Galactic ultraluminous X-ray source (ULX) Swift J0243.6+6124. The pulse fractions of 19 observation pairs selected in the rising and fading phases with similar luminosity are investigated. The results show a general trend of the pulse fraction increasing with luminosity and energy at super-critical luminosity. However, the relative strength of the pulsation between each pair evolves strongly with luminosity. The pulse fraction in the rising phase is larger at luminosity below $7.71\times10^{38}$~erg~s$^{-1}$, but smaller at above. A transition luminosity is found to be energy independent. Such a phenomena is firstly confirmed by Insight-HXMT observations and we speculate it may have relation with the radiation pressure dominated accretion disk.

We report on the Insight-HXMT observations of the new black hole X-ray binary MAXI J1820+070 during its 2018 outburst. Detailed spectral analysis via the continuum fitting method shows an evolution of the inferred spin during its high soft sate. Moreover, the hardness ratio, the non-thermal luminosity and the reflection fraction also undergo an evolution, exactly coincident to the period when the inferred spin transition takes place. The unphysical evolution of the spin is attributed to the evolution of the inner disc, which is caused by the collapse of a hot corona due to condensation mechanism or may be related to the deceleration of a jet-like corona. The studies of the inner disc radius and the relation between the disc luminosity and the inner disc radius suggest that, only at a particular epoch, did the inner edge of the disc reach the innermost stable circular orbit and the spin measurement is reliable. We then constrain the spin of MAXI J1820+070 to be a*=0.2^{+0.2}_{-0.3}. Such a slowly spinning black hole possessing a strong jet suggests that its jet activity is driven mainly by the accretion disc rather than by the black hole spin., Comment: 14 pages, 13 figures, 5 tables, accepted for publication in MNRAS

We investigate the absorption and emission features in observations of GX 301-2 detected with Insight-HXMT/LE in 2017-2019. At different orbital phases, we found prominent Fe Kalpha, Kbeta and Ni Kalpha lines, as well as Compton shoulders and Fe K-shell absorption edges. These features are due to the X-ray reprocessing caused by the interaction between the radiation from the source and surrounding accretion material. According to the ratio of iron lines Kalpha and Kbeta, we infer the accretion material is in a low ionisation state. We find an orbital-dependent local absorption column density, which has a large value and strong variability around the periastron. We explain its variability as a result of inhomogeneities of the accretion environment and/or instabilities of accretion processes. In addition, the variable local column density is correlated with the equivalent width of the iron Kalpha lines throughout the orbit, which suggests that the accretion material near the neutron star is spherically distributed., Comment: 10 pages, 5 figures, 2 tables, accepted for publication in MNRAS

We present the analysis of the brightest flare that was recorded in the \emph{Insight}-HMXT data set, in a broad energy range (2$-$200 keV) from the microquasar GRS~1915+105 during an unusual low-luminosity state. This flare was detected by \emph{Insight}-HXMT among a series of flares during 2 June 2019 UTC 16:37:06 to 20:11:36, with a 2-200 keV luminosity of 3.4$-$7.27$\times10^{38}$ erg s$^{-1}$. Basing on the broad-band spectral analysis, we find that the flare spectrum shows different behaviors during bright and faint epochs. The spectrum of the flare can be fitted with a model dominated by a power-law component. Additional components show up in the bright epoch with a hard tail and in the faint epoch with an absorption line $\sim$ 6.78 keV. The reflection component of the latter is consistent with an inner disk radius $\sim$ 5 times larger than that of the former. These results on the giant flare during the "unusual" low-luminosity state of GRS~1915+105 may suggest that the source experiences a possible fast transition from a jet-dominated state to a wind-dominated state. We speculate that the evolving accretion disk and the large-scale magnetic field may play important roles in this peculiar huge flare.

After 34 years, the black-hole candidate EXO 1846-031 went into outburst again in 2019. We investigate its spectral properties in the hard intermediate and the soft states with NuSTAR and Insight-HXMT. A reflection component has been detected in the two spectral states but possibly originating from different illumination spectra: in the intermediate state, the illuminating source is attributed to a hard coronal component, which has been commonly observed in other X-ray binaries, whereas in the soft state the reflection is probably produced by the disk self-irradiation. Both cases support EXO 1846-031 as a low inclination system of ~40 degrees. An absorption line is clearly detected at ~7.2 keV in the hard intermediate state, corresponding to a highly ionized disk wind (log {\xi} > 6.1) with a velocity up to 0.06c. Meanwhile, quasi-simultaneous radio emissions have been detected before and after the X-rays, implying the co-existence of disk winds and jets in this system. Additionally, the observed wind in this source is potentially driven by magnetic forces. The absorption line disappeared in the soft state and a narrow emission line appeared at ~6.7 keV on top of the reflection component, which may be evidence for disk winds, but data with the higher spectral resolution are required to examine this., Comment: Accepted for publication in ApJ

We systematically analyze three GRB samples named as radio-loud, radio-quiet and radio-none afterglows, respectively. It is shown that dichotomy of the radio-loud afterglows is not necessary. Interestingly, we find that the intrinsic durations ($T_{int}$), isotropic energies of prompt gamma-rays ($E_{\gamma, iso}$) and redshifts ($z$) of their host galaxies are log-normally distributed for both the radio-loud and radio-quiet samples except those GRBs without any radio detections. Based on the distinct distributions of $T_{int}$, $E_{\gamma, iso}$, the circum-burst medium density ($n$) and the isotropic equivalent energy of radio afterglows ($L_{\nu,p}$), we confirm that the GRB radio afterglows are really better to be divided into the dim and the bright types. However, it is noticeable that the distributions of flux densities ($F_{host}$) from host galaxies of both classes of radio afterglows are intrinsically quite similar. Meanwhile, we point out that the radio-none sample is also obviously different from the above two samples with radio afterglows observed, according to the cumulative frequency distributions of the $T_{int}$ and the $E_{\gamma, iso}$, together with correlations between $T_{int}$ and $z$. In addition, a positive correlation between $E_{\gamma, iso}$ and $L_{\nu,p}$ is found in the radio-loud samples especially for the supernova-associated GRBs. Besides, we also find this positive correlation in the radio-quiet sample. A negative correlation between $T_{int}$ and $z$ is confirmed to hold for the radio-quiet sample too. The dividing line between short and long GRBs in the rest frame is at $T_{int}\simeq$1 s. Consequently, we propose that the radio-loud, the radio-quiet and the radio-none GRBs could be originated from different progenitors., Comment: 20 pages, 30 figures, submitted to MNRAS

A link between magnetars and fast radio burst (FRB) sources has finally been established. In this context, one of the open issues is whether/which sources of extra galactic FRBs exhibit X/gamma-ray outbursts and whether it is correlated with radio activity. We aim to constrain possible X/gamma-ray burst activity from one of the nearest extragalactic FRB sources currently known over a broad energy range, by looking for bursts over a range of timescales and energies that are compatible with being powerful flares from extragalactic magnetars. We followed up the as-yet nearest extragalactic FRB source at a mere 149 Mpc distance, the periodic repeater FRB180916.J0158+65, during the active phase on February 4-7, 2020, with the Insight-Hard X-ray Modulation Telescope (HXMT). Taking advantage of the combination of broad band, large effective area, and several independent detectors available, we searched for bursts over a set of timescales from 1 ms to 1.024 s with a sensitive algorithm, that had previously been characterised and optimised. Moreover, through simulations we studied the sensitivity of our technique in the released energy-duration phase space for a set of synthetic flares and assuming different energy spectra. We constrain the possible occurrence of flares in the 1-100 keV energy band to E<10^46 erg for durations <0.1 s over several tens of ks exposure. We can rule out the occurrence of giant flares similar to the ones that were observed in the few cases of Galactic magnetars. The absence of reported radio activity during our observations does not allow us to make any statements on the possible simultaneous high-energy emission., Comment: 10 pages, 3 figures, accepted by A&A

With the twelve-year long observations by {\sl Fermi}-LAT, we discover two pairs of spin-down state transitions of PSR J1124--5916, making it the second young pulsar detected to have such behaviors. PSR J1124--5916 shows mainly two states according to its spin-down rate evolution, the normal spin-down state and the low spin-down state. In about 80\% of the observation time, the pulsar is in the normal spin-down state, in which the spin-down rate decreases linearly and gives a braking index of $1.98\pm0.04$. The two transitions to the low spin-down state are in MJD 55183--55803 and MJD 56114--56398 respectively, with fractional amplitudes both $\sim0.4\%$. No significant difference between the $\gamma$-ray profiles of the two spin-down states is detected, which is similar to PSR B0540-69, the other young pulsar with state transition detected., Comment: 13 pages, 2 figures

We report on the observation of the accreting pulsar GRO J1008-57 performed by Insight-HXMT at the peak of the source's 2017 outburst. Pulsations are detected with a spin period of 93.283(1) s. The pulse profile shows double peaks at soft X-rays, and only one peak above 20 keV. The spectrum is well described by the phenomenological models of X-ray pulsars. A cyclotron resonant scattering feature is detected with very high statistical significance at a centroid energy of $E_{\rm cyc}=90.32_{-0.28}^{+0.32}$ keV, for the reference continuum and line models, HIGHECUT and GABS respectively. Detection is very robust with respect to different continuum models. The line energy is significantly higher than what suggested from previous observations, which provided very marginal evidence for the line. This establishes a new record for the centroid energy of a fundamental cyclotron resonant scattering feature observed in accreting pulsars. We also discuss the accretion regime of the source during the Insight-HXMT observation., Comment: 8 pages, 3 figures, accepted for publication in ApJL

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Some radio pulses have been observed from some magnetars, however, no FRB-like events had been detected in association any magnetar burst, including one giant flare. Recently, a pair of FRB-like bursts (FRB 200428 hereafter) separated by milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR J1935+2154. The burst showed two hard peaks with a separation of 34 ms, broadly consistent with that of the two bursts in FRB 200428. The delay time between the double radio and X-ray peaks is about 8.57 s, fully consistent with the dispersion delay of FRB 200428. We thus identify the non-thermal X-ray burst is associated with FRB 200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB 200428 share the same physical origin in an explosive event from SGR J1935+2154., Comment: 24 pages, 9 figures, 6 tables; initial submission to a journal on May 9th, 2020. Significant changes include updated localization and detailed spectral evolution of the X-ray burst, and better determination of the two narrow X-ray peaks corresponding to the two radio pulses. Conclusions are strengthened. Nature Astronomy online on Feb. 18, 2021

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc., Comment: 8 pages, 5 figures, accepted by JHEAP

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet., Comment: Accepted for publication in ApJ

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches., Comment: 25 pages, 8 figures

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB., Comment: 15 pages, 3 figures, 6 tables, accepted by A&A

In Dec. 2011 PSR B0540-69 experienced a spin-down rate transition (SRT), after which the spin-down power of the pulsar increased by ~36%. About 1000 days after the SRT, the X-ray luminosity of the associated pulsar wind nebula (PWN) was found to brighten by 32+/-8%. After the SRT, the braking index n of PSR B0540-69 changes from n=2.12 to n=0.03 and then keeps this value for about five years before rising to n=0.9 in the following years. We find that most of the current models have difficulties in explaining the measured braking index. One exceptive model of the braking index evolution is the increasing dipole magnetic field of PSR B0540-69. We suggest that the field increase may result from some instabilities within the pulsar core that enhance the poloidal component at the price of toroidal component of the magnetic field. The increasing dipole magnetic field will result in the X-ray brightening of the PWN. We fit the PWN X-ray light curve by two models: one assumes a constant magnetic field within the PWN during the brightening and the other assumes an enhanced magnetic field proportional to the energy density of the PWN. It appears that the two models fit the data well, though the later model seems to fit the data a bit better. This provides marginal observational evidence that magnetic field in the PWN is generated by the termination shock. Future high-quality and high-cadence data are required to draw a solid conclusion., Comment: 7 pages, 2 figures, the published version

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate., Comment: 11 pages, 8 figures, accepted for publication in MNRAS

In the current unmanned aircraft systems (UASs) for sensing services, unmanned aerial vehicles (UAVs) transmit their sensory data to terrestrial mobile devices over the unlicensed spectrum. However, the interference from surrounding terminals is uncontrollable due to the opportunistic channel access. In this paper, we consider a cellular Internet of UAVs to guarantee the Quality-of-Service (QoS), where the sensory data can be transmitted to the mobile devices either by UAV-to-Device (U2D) communications over cellular networks, or directly through the base station (BS). Since UAVs' sensing and transmission may influence their trajectories, we study the trajectory design problem for UAVs in consideration of their sensing and transmission. This is a Markov decision problem (MDP) with a large state-action space, and thus, we utilize multi-agent deep reinforcement learning (DRL) to approximate the state-action space, and then propose a multi-UAV trajectory design algorithm to solve this problem. Simulation results show that our proposed algorithm can achieve a higher total utility than policy gradient algorithm and single-agent algorithm., Comment: 33 pages, 12 figures

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.