Search

We present an energy-dependent analysis for the type-C quasi-periodic oscillations (QPOs) observed in the black hole X-ray binary Swift J1727.8-1613 using Insight-HXMT observations. We find that the QPO fractional rms at energies above 40 keV is significantly higher than that below 20 keV. This is the first report of a high energy (HE)-rms excess in the rms spectrum of a black hole X-ray binary. In the high energy band, an extra hard component is observed in additional to the standard thermal Comptonization component at similar energy band. The value of the QPO HE-rms excess is not only correlated with the disk parameters and the photon index of the standard Comptonization component, but also exhibits a moderate positive correlation with the flux of the additional hard spectral component. No features in the QPO phase-lag spectra are seen corresponding to the additional hard component. We propose that the additional hard component in the spectrum may originate from jet emission and the associated QPO HE-rms excess can be explained by the precession of the jet base.

We present a comprehensive timing analysis on the Be/X-ray binary pulsar RX J0440.9+4431 using observations from \textit{NICER} and \textit{Insight}-HXMT during the 2022--2023 outburst. The power density spectrum (PDS) of RX J0440.9+4431 exhibits typical aperiodic variability in X-ray flux across a wide frequency range. During a super-critical accretion state, we detect quasi-periodic oscillations (QPOs) at 0.2--0.5\,Hz in the light curves of five pulses for RX J0440.9+4431. The observed QPOs manifest during flares, while the flares appear at the peaks of the pulse profiles on a timescale of seconds and are primarily caused by an increase in hard photons. These flares can be explained by increased material ingestion in the accretion column at a fixed phase, primarily generating hard photons. Alternatively, an increase in accretion rate, independent of phase, may result in highly beamed hard photons within the accretion column, causing the flares. We argue the origin of QPOs to instabilities within the accretion flow. Additionally, we find that the break frequencies in the noise power spectra align well with $\propto L_{\mathrm{x}}^{3 / 7}$ across three orders of magnitude in the luminosity, which points to a relatively strong magnetic field in RX J0440.9+4431, estimated to be \textasciitilde$10^{13}$\,G., Comment: 8 pages, 7 figures. Accepted in MNRAS

We report the phase-resolved spectral results of the first Galactic Pulsating Ultra-Luminous X-ray source (PULX) Swift J0243.6+6124, modeling at its 2017-2018 outburst peak using data collected by the Hard X-ray Modulation Telescope (Insight-HXMT). The broad energy coverage of Insight-HXMT allows us to obtain more accurate spectral continuum to reduce the coupling of broad iron line profiles with other components. We use three different continuum spectrum models but obtain similar iron line results. For the first time, we detected the pulse characteristics of the broad iron line in a PULX. The variation in width and intensity of this iron line with $\sigma \sim 1.2-1.5$\,keV has a phase offset of about 0.25 from the pulse phase. We suggest that the uneven irradiation of the thick inner disk by the accretion column produces the modulated variation of the broad iron line. In addition, the non-pulsed narrow line is suggested to come from the outer disk region.

RX J0440.9+4431, a Be/X-ray binary, had its brightest outburst in 2022 since its discovery, with a peak X-ray flux of 2.25 Crab (as recorded by Swift/BAT, 15-50 keV). We analyze the timing properties of this giant outburst using data from Insight-HXMT and NICER, focusing on the evolution of the pulse profile and pulse fraction. We observe that when the luminosity reached around ~ 3*10^{37} er s^{-1}, a transition from double-peaked to single-peaked pulse profiles occurred across the energy range, with the peak of the low-energy profile aligning gradually with the peak of the high-energy profile. This change indicates a transition from subcritical to supercritical accretion. Additionally, we found a concave in the pulse fraction as a function of energy around 20-30 keV throughout the entire outburst period. Compared to the low luminosity, the concave becomes weaker in high luminosities, and overall, the pulse fraction is higher. We propose that this concave could be caused by the scattering of high-energy photons by the atmosphere of a neutron star, leading to a dilution of the pulse fraction. As the accretion reaches the supercritical state, the accretion column height increases, resulting in a larger direct component of strongly beamed X-ray flux, and an elevated pulse fraction., Comment: 15 pages, 8 figures, accepted for publication in MNRAS

MAXI J181-195 is a newly discovered accreting millisecond X-ray pulsar that went outburst in June 2022. Through timing analysis with NICER and NuSTAR observations, we find a transient modulation at ~2.5 Hz during the decay period of MAXI J1816-195. The modulation is strongly correlated with a spectral hardening, and its fractional rms amplitude increases with energy. These results suggest that the modulation is likely to be produced in an unstable corona. In addition, the presence of the modulation during thermonuclear bursts indicates that it may originate from a disk-corona where the optical depth is likely the main factor affecting the modulation, rather than temperature. Moreover, we find significant reflection features in the spectra observed simultaneously by NICER and NuSTAR, including a relativistically broadened Fe-K line around 6-7 keV, and a Compton hump in the 10-30 keV energy band. The radius of the inner disc is constrained to be Rin = (1.04-1.23) RISCO based on reflection modeling of the broadband spectra. Assuming that the inner disc is truncated at the magnetosphere radius, we estimate that the magnetic field strength is < 4.67 * 10e8 G., Comment: 12 pages, 13 figures

We investigate the dynamical properties of low frequency quasi-periodic oscillations (QPOs) observed from the black hole X-ray binary MAXI J1820+070 during the early part of its 2018 outburst, when the system was in a bright hard state. To this aim, we use a series of observations from the Hard X-ray Modulation Telescope Insight-HXMT, and apply a wavelet decomposition (weighted wavelet Z-transforms) to the X-ray light-curve. We find that the QPO phenomenon is intermittent within each individual observation, with some sub-intervals where the oscillation is strongly detected (high root-mean-square amplitude) and others where it is weak or absent. The average life time of individual QPO segments is ~ 5 oscillation cycles, with a 3 sigma tail up to ~ 20 cycles. There is no substantial difference between the energy spectra during intervals with strong and weak/absent QPOs. We discuss two possible reasons for the intermittent QPO strength, within the precessing jet model previously proposed for MAXI J1820+070. In the rigid precession model, intermittent QPOs are predicted to occur with a coherence Q ~ a few when the disk alignment time-scale is only a few times the precession time-scale. Alternatively, we suggest that changes in oscillation amplitude can be caused by changes in the jet speed. We discuss a possible reason for the intermittent QPO strength, within the precessing jet model previously proposed for MAXI J1820+070: we suggest that changes in oscillation amplitude are caused by changes in the jet speed. We argue that a misaligned, precessing jet scenario is also consistent with other recent observational findings that suggest an oscillation of the Compton reflection component in phase with the QPOs., Comment: 8 pages, 4 figures

We report the X-ray timing results of the black hole candidate MAXI J1820+070 during its 2018 outburst using the Hard X-ray Modulation Telescope (Insight-HXMT) and Neutron Star Interior Composition Explorer Mission (NICER) observations. Low frequency quasi-periodic oscillations (LFQPOs) are detected in the low/hard state and the hard intermediate state, which lasted for about 90 days. Thanks to the large effective area of Insight-HXMT at high energies and NICER at low energies, we are able to present the energy dependence of the LFQPO characteristics and phase lags from 0.2 keV to 200 keV, which has never been explored by previous missions. We find that the centroid frequency of the LFQPOs do not change significantly with energy, while the full width at half maximum (FWHM) and fractional rms show a complex evolution with energy. The LFQPO phase lags at high energies and low energies show consistent energy-dependence relations taking the ~2 keV as reference. Our results suggest that the LFQPOs from high energy come from the LT precession of the relativistic jet, while the low-energy radiation is mainly from the perpendicular innermost regions of the accretion disk.

We report the evolution of the X-ray pulsations of EXO 2030+375 during its 2021 outburst using the observations from \textit{Insight}-HXMT. Based on the accretion torque model, we study the correlation between the spin frequency derivatives and the luminosity. Pulsations can be detected in the energy band of 1--160 keV. The pulse profile evolves significantly with luminosity during the outburst, leading to that the whole outburst can be divided into several parts with different characteristics. The evolution of the pulse profile reveals the transition between the super-critical (fan-beam dominated) and the sub-critical accretion (pencil-beam dominated) mode. From the accretion torque model and the critical luminosity model, based on a distance of 7.1 kpc, the inferred magnetic fields are $(0.41-0.74) \times 10^{12}$ G and $(3.48-3.96) \times 10^{12}$ G, respectively, or based on a distance of 3.6 kpc, the estimated magnetic fields are $(2.4-4.3) \times 10^{13}$ G and $(0.98-1.11)\times 10^{12}$ G, respectively. Two different sets of magnetic fields both support the presence of multipole magnetic fields of the NS.

A double-peak X-ray burst from the Galactic magnetar SGR J1935+2154 was discovered as associated with the two radio pulses of FRB 200428 separated by 28.97+-0.02 ms. Precise measurements of the timing and spectral properties of the X-ray bursts are helpful for understanding the physical origin of fast radio bursts (FRBs). In this paper, we have reconstructed some information about the hard X-ray events, which were lost because the High Energy X-ray Telescope (HE) onboard the Insight-HXMT mission was saturated by this extremely bright burst, and used the information to improve the temporal and spectral analyses of the X-ray burst. The arrival times of the two X-ray peaks by fitting the new Insight-HXMT/HE lightcurve with multi-Gaussian profiles are 2.77+-0.45 ms and 34.30+-0.56 ms after the first peak of FRB 200428, respectively, while these two parameters are 2.57+-0.52 ms and 32.5+-1.4 ms if the fitting profile is a fast rise and exponential decay function. The spectrum of the two X-ray peaks could be described by a cutoff power-law with cutoff energy ~60 keV and photon index ~1.4, the latter is softer than that of the underlying bright and broader X-ray burst when the two X-ray peaks appeared., Comment: 11 pages, 7 figures

We present a systematic analysis of type C quasi-periodic oscillation (QPO) observations of H 1743--322 throughout the Rossi X-ray Timing Explorer (RXTE) era. We find that, while different outbursts have significant flux differences, they show consistent positive correlations between the QPO fractional root-mean-square (rms) amplitude and non-thermal fraction of the emission, which indicate an independence of the intrinsic QPO rms on individual outburst brightness in H 1743--322. However, the dependence of the QPO rms on frequency is different between the outburst rise and decay phases, where QPO fractional rms of the decay phase is significantly lower than that of the rise phase at low frequencies. The spectral analysis also reveals different ranges of coronal temperature between the two outburst stages. A semi-quantitative analysis shows that the Lense-Thirring precession model could be responsible for the QPO rms differences, requiring a variable coronal geometric shape. However, the variable-Comptonization model could also account for the findings. The fact that the rms differences and the hysteresis traces in the hardness-intensity diagram (HID) accompany each other indicates a connection between the two phenomena. By correlating the findings with QPO phase lags and the quasi-simultaneous radio flux previously published, we propose there could be corona-jet transitions in H 1743--322 similar to those that have been recently reported in GRS 1915+105., Comment: 21 pages, 12 figures

Using wavelet analysis and power density spectrum, we investigate two transient quasi-periodic oscillations (QPOs) observed in MAXI J1535$-$571 observed with Insight-HXMT. The transient QPOs have a centroid frequency of $\sim 10$ Hz with a FWHM $\sim 0.6$ Hz and an rms amplitude $\sim 14\%$. Energy spectra of QPO and non-QPO regimes are also separated and analyzed, and the spectra become softer with higher $E_{cut}$ in the non-QPO regime compared to the QPO regime. Our results suggest that the transient QPOs detected in MJD 58016 and 58017 are still the type-C QPO, and the source remains in its HIMS. The duration of all type-C QPO signals based on wavelet is positively correlated with the mean count rate above $\sim 10$ keV, implying appearance of QPOs in different time scales should be coupled with the corona. The transient QPO properties could be related to the jet or flares, perhaps the partial ejection of the corona is responsible for the disappearance of the type-C QPO., Comment: 10 pages, MNRAS in press

We investigated the quasi-periodic oscillation (QPO) features in the accretion-powered X-ray pulsar Cen X-3 observed by Insight-HXMT. For two observations in 2020 when Cen X-3 was in an extremely soft state, the power density spectrum revealed the presence of obvious QPO features at $\sim$40 mHz with an averaged fractional rms amplitude of $\sim9\%$. We study the mHz QPO frequency and rms amplitude over orbital phases, and find that the QPO frequency is $\sim$33-39 mHz at the orbital phase of 0.1-0.4, increasing to $\sim$37-43 mHz in the orbital phase of 0.4-0.8, but has no strong dependence on X-ray intensity. We also carried out an energy-dependent QPO analysis, the rms amplitude of the mHz QPOs have a decreasing trend as the energy increases from 2 to 20 keV. In addition, the QPO time-lag analysis shows that the time delay is $\sim 20$ ms (a hard lag) in the range of $\sim$5-10 keV, and becomes negative (time lag of $-(20-70)$ ms) above $\sim 10$ keV. The different QPO theoretical models are summarized and discussed. In the end, we suggest that these energy-dependent timing features as well as the origin of mHz QPOs in Cen X-3 may be ascribed to an instability when the accretion disk is truncated near the corotation radius., Comment: 9 pages, MNRAS in press

Cygnus X-1, as the first discovered black hole binary, is a key source for understanding the mechanisms of state transitions, and the scenarios of accretion in extreme gravity fields. We present a spectral-timing analysis of observations taken with the Insight-HXMT mission, focusing on the spectral-state dependent timing properties in the broad energy range of 1--150 keV, thus extending previous RXTE-based studies to both lower and higher energies. Our main results are the following: a) We successfully use a simple empirical model to fit all spectra, confirming that the reflection component is stronger in the soft state than in the hard state; b) The evolution of the total fractional root mean square (rms) depends on the selected energy band and the spectral shape, which is a direct result of the evolution of the power spectral densities (PSDs); c) In the hard/intermediate state, we see clear short-term variability features and a positive correlation between central frequencies of the variability components and the soft photon index $\Gamma_1$, also at energies above 15 keV. The power spectrum is dominated by red noise in the soft state instead. These behaviors can be traced to at least 90 keV; d) The coherence and the phase-lag spectra show different behaviors dependent on different spectral shapes., Comment: 18 pages, 16 figures, 1 table. accepted for publication in Astronomy & Astrophysics

We report results on the timing analysis of the 2020 giant outburst of 1A 0535+262, using broadband data from Insight-HXMT. The analysis of the pulse profile evolution from the sub-critical luminosity to super-critical luminosity regime is presented for the first time. We found that the observed pulse profile exhibits a complex dependence on both energy and luminosity.A dip structure at the energy of the cyclotron resonant scattering features (CRSFs) is found for the first time in the pulse fraction-energy relation of 1A 0535+262, when the outburst evolves in a luminosity range from 4.8 $\times 10^{37}$ erg s$^{-1}$ to 1.0 $\times 10^{38}$ erg s$^{-1}$. The observed structure is luminosity dependent and appears around the source critical luminosity ($\sim$ 6.7 $\times 10^{37}$ erg s$^{-1}$).

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

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.

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

We execute a detailed spectral-timing study of Cygnus X-1 in the low/hard, intermediate and high/soft states with the Hard X-ray Modulation Telescope observations. The broad band energy spectra fit well with the "truncated disk model" with the inner boundary of the accretion disk stays within $\sim$10 \textit R$_{\rm g}$ and moves inward as the source softens. Through studying of PDS, rms and Fourier-frequency component resolved spectroscopy, we find that the X-ray variations are generated in two different regions for each state. We discover that the major role that contributes to the X-ray variation is the hot corona rather than the accretion disk. We suggest a scenario with different corona geometry for each state based on the truncated disk geometry in which the corona wraps up the disk to form a sandwich geometry in the low/hard state, and then gradually moves away from the disk in direction that is perpendicular to the disk until forming a jet-like geometry in the high/soft state., Comment: 13 pages, 12 figures, and 3 tables, accepted for publication in ApJ

Electron beam polarization in the bubble regime of the interaction between a high-intensity laser and a longitudinally pre-polarized plasma is investigated by means of the Thomas-Bargmann-Michel-Telegdi equation. Using a test-particle model, the dependence of the accelerated electron polarization on the bubble geometry is analyzed in detail. Tracking the polarization dynamics of individual electrons reveals that although the spin direction changes during both the self-injection process and acceleration phase, the former has the biggest impact. For nearly spherical bubbles, the polarization of electron beam persists after capture and acceleration in the bubble. By contrast, for aspherical bubble shapes, the electron beam becomes rapidly depolarized, and the net polarization direction can even reverse in the case of a oblate spheroidal bubble. These findings are confirmed via particle-in-cell simulations., Comment: 8pages, 4 figures

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

We investigate systematically four outbursts of black hole system GX 339-4 observed by the Rossi X-ray Timing Explorer (RXTE) in both spectral and timing domains and find that these outbursts have some common properties although they experience different 'q' tracks in the hardness-intensity diagram (HID). While the spectral indices are around 1.5 in low hard state (LHS), 2.4 in soft intermediate state (SIMS) and high soft state (HSS), the spectral parameters of thermal, non-thermal and reflection components vary significantly in transitions from LHS to HIMS. Also the quasi periodic oscillation (QPO) shows a peculiar behavior during the state transition between LHS and HIMS: the RMS drop of type C fundamental QPO is accompanied with showing-up of the second harmonic. Interestingly, the QPO RMS is found to have a similar linear relationship with the non-thermal fraction of emission in different outbursts. These findings provide more clues to our understanding the outburst of the black hole X-ray binary system., Comment: 13 papes, 12 figures

We report on a detailed spectral analysis of the transient X-ray pulsar 1A~0535+262, which underwent the brightest giant outburst ever recorded for this source from November to December 2020 with a peak luminosity of $1.2$ $\times10^{38}\ \rm erg\ s^{-1}$. Thanks to the unprecedented energy coverage and high cadence observations provided by Insight-HXMT, we were able to find for the first time evidence for a transition of the accretion regime. At high luminosity, above the critical luminosity $6.7\times10^{37}$ erg s$^{-1}$, the cyclotron absorption line energy anti-correlates with luminosity. Below the critical luminosity, a positive correlation is observed. The 1A~0535+262 becomes, therefore, the second source after V~0332+53, which clearly shows an anti-correlation above and transition between correlation and anti-correlation around the critical luminosity. The evolution of both the observed CRSF line energy and broadband X-ray continuum spectrum throughout the outburst exhibits significant differences during the rising and fading phases: that is, for a similar luminosity the spectral parameters take different values which results in hysteresis patterns for several spectral parameters including the cyclotron line energy. We argue that, similarly to V~0332+53, these changes might be related to different geometry of the emission region in rising and declining parts of the outburst, probably due to changes in the accretion disk structure and its interaction with the magnetosphere of the neutron star.

We report the energy-resolved broadband timing analysis of the black hole X-ray transient MAXI J1631-479 during its 2019 outburst from February 11 to April 9, using data from the Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which caught the source from its hard intermediate state to the soft state. Thanks to the large effective area of Insight-HXMT at high energies, we are able to present the energy dependence of fast variability up to ~100 keV. Type-C quasi-periodic oscillations (QPOs) with frequency varying between 4.9 Hz and 6.5 Hz are observed in the 1-100 keV energy band. While the QPO fractional rms increases with photon energy from 1 keV to ~10 keV and remains more or less constant from ~10 keV to ~100 keV, the rms of the flat-top noise first increases from 1 keV to ~8 keV then drops to less than 0.1% above ~30 keV. We suggest that the disappearance of the broadband variability above 30 keV could be caused by the non-thermal acceleration in the Comptonizing plasma. At the same time, the QPOs could be produced by the precession of either a small-scale jet or a hot inner flow model., Comment: To be published in The Astrophysical Journal

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

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.

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

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

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

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