Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 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 $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.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.

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.

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.

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.

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.

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.

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.

LS 5039 is a well-known $\gamma$-ray binary system which consists of an unknown compact object and a massive companion O star. It shows rather stable emissions at high energies over years and hence serves as an ideal laboratory to investigate the emission mechanism for such peculiar systems which emit prominent $\gamma$-rays. To this end, we take the orbital phase resolved energy spectrum as observed by \fermi over 10 years. We divide the orbit into four orbital phases, each with an orbital phase range of 0.25, centered at 0.00, 0.25, 0.50 and 0.75 respectively, where the phase 0.0 is the periastron and phase 0.5 is the apastron. The phases around 0.25 and 0.75 are symmetric and hence are supposed to have identical local acceleration environment. The spectral analysis shows that, the \fermi spectra are largely different from these two symmetric orbital phases: the emission from orbital phase 0.25 turns out to be significantly stronger than that from 0.75. This result does not fit a scenario that $\gamma$-rays are Doppler boosted emission from bow shock tails if LS 5039 has a shock configuration similar to PSR B1259-63, and indicates that the inverse Compton scatterings between the shock accelerated plasma and the stellar particle environment is the underline procedure. We also find that the previous report for a disappearance of the orbital modulation at 3--20 GeV is due to the similar spectral turn-over energies of the different orbital phases. The spectral properties of periastron and apastron regions are addressed in the context of the measurements in phase regions around 0.25 and 0.75.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

The hard X-ray transient source Swift J1845.7-0037 was discovered in 2012 by Swift/BAT. However, at that time no dedicated observations of the source were performed. On Oct 2019 the source became active again, and X-ray pulsations with a period of ~199s were detected with Swift/XRT. This triggered follow-up observations with NuSTAR. Here we report on the timing and spectral analysis of the source properties using NuSTAR and Swift/XRT. The main goal was to confirm pulsations and search for possible cyclotron lines in the broadband spectrum of the source to probe its magnetic field. Despite highly significant pulsations with period of 207.379(2) were detected, no evidence for a cyclotron line was found in the spectrum of the source. We therefore discuss the strength of the magnetic field based on the source flux and the detection of the transition to the "cold-disc" accretion regime during the 2012 outburst. Our conclusion is that, most likely, the source is a highly magnetized neutron star with B 1e13G at a large distance of d~10 kpc. The latter one consistent with the non-detection of a cyclotron line in the NuSTAR energy band.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

PSR B1259-63 is a $\gamma$-ray emitting high mass X-ray binary system, in which the compact object is a millisecond pulsar. The system has an orbital period of 1236.7 d and shows peculiar $\gamma$-ray flares when the neutron star moves out of the stellar disk of the companion star. The $\gamma$-ray flare events were firstly discovered by using Fermi-LAT around the 2010 periastron passage, which was repeated for the 2014 and 2017 periastron passages. We analyze the Fermi-LAT data for all the three periastron passages and found that in each flare the energy spectrum can be represented well by a simple power law. The $\gamma$-ray light curves show that in 2010 and 2014 after each periastron there are two main flares, but in 2017 there are four flares including one precursor about 10 d after the periastron passage. The first main flares in 2010 and 2014 are located at around 35 d after the periastron passage, and the main flare in 2014 is delayed by roughly 1.7 d with respect to that in 2010. In the 2017 flare, the source shows a precursor about 10 d after the periastron passage, but the following two flares become weaker and lag behind those in 2010 by roughly 5 d. The strongest flares in 2017 occurred 58 d and 70 d after the periastron passage. These results challenge the previous models.

We discuss a counterintuitive phenomenon of classical general relativity, in which a significant fraction of the radiation emitted by a collapsing object and detected by a distant observer may be blueshifted rather than redshifted. The key-point is that when the radiation propagates inside the collapsing body it is blueshifted, and this time interval may be sufficiently long for the effect to be larger than the later redshift due to the propagation in the vacuum exterior, from the surface of the body to the distant observer. Unfortunately, the phenomenon can unlikely have direct observational implications, but it is interesting by itself as a pure relativistic effect.

In a previous paper, one of us has described a code to compute the thermal spectrum of geometrically thin and optically thick accretion disks around generic stationary and axisymmetric black holes, which are not necessarily of the Kerr type. As the structure of the accretion disk and the propagation of electromagnetic radiation from the disk to the distant observer depend on the background metric, the analysis of the thermal spectrum of thin disks can be used to test the actual nature of black hole candidates. In this paper, we consider the 10 stellar-mass black hole candidates for which the spin parameter has been already estimated from the analysis of the disk's thermal spectrum and under the assumption of the Kerr background, and we translate the measurements reported in the literature into constraints on the spin parameter--deformation parameter plane. The analysis of the disk's thermal spectrum can be used to estimate only one parameter of the geometry close to the compact object, and therefore it is not possible to get independent measurements of both the spin and the deformation parameters. The constraints obtained here will be used in combination with other measurements in future work, with the final goal to break the degeneracy between the spin and possible deviations from the Kerr solution and thus test the Kerr black hole hypothesis.

In general relativity, gravitational collapse of matter fields ends with the formation of a spacetime singularity, where the matter density becomes infinite and standard physics breaks down. According to the weak cosmic censorship conjecture, singularities produced in the gravitational collapse cannot be seen by distant observers and must be hidden within black holes. The validity of this conjecture is still controversial and at present we cannot exclude that naked singularities can be created in our Universe from regular initial data. In this paper, we study the radiation emitted by a collapsing cloud of dust and check whether it is possible to distinguish the birth of a black hole from the one of a naked singularity. In our simple dust model, we find that the properties of the radiation emitted in the two scenarios is qualitatively similar. That suggests that observational tests of the cosmic censorship conjecture may be very difficult, even in principle.

General relativity makes clear predictions about the spacetime geometry around black holes. In the near future, new facilities will have the capability to explore the metric around SgrA$^*$, the supermassive black hole candidate at the Center of our Galaxy, and open a new window to test the Kerr black hole hypothesis. In this paper, we compute light curves and images associated with compact emission regions (hot spots) orbiting around Kerr and non-Kerr black holes. We study how the analysis of the properties of the radiation emitted by a hot spot can be used to test the Kerr nature of SgrA$^*$. We find that the sole observation of the hot spot light curve can at most constrain a combination of the black hole spin and of possible deviations from the Kerr solution. This happens because the same orbital frequency around a Kerr black hole can be found for a non-Kerr object with a different spin parameter. Second order corrections in the light curve due to the background geometry are typically too small to be identified. While the observation of the hot spot centroid track can potentially bound possible deviations from the Kerr solution, that is out of reach for the near future VLTI instrument GRAVITY. The Kerr black hole hypothesis could really be tested in the case of the discovery of a radio pulsar in a compact orbit around SgrA$^*$. Radio observations of such a pulsar would provide precise estimates of the mass and the spin of SgrA$^*$, and the combination of these measurements (probing the weak field) with the hot spot light curve information (probing the strong field) may constrain/find possible deviations from the Kerr solution with quite good precision.