We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities.

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.

We report here on initial results from the Thousand Pulsar Array (TPA) programme, part of the Large Survey Project "MeerTime" on the MeerKAT telescope. The interferometer is used in tied-array mode in the band from 856 to 1712~MHz, and the wide band coupled with the large collecting area and low receiver temperature make it an excellent telescope for the study of radio pulsars. The TPA is a 5 year project which aims to observe (a) more than 1000 pulsars to obtain high-fidelity pulse profiles, (b) some 500 of these pulsars over multiple epochs, (c) long sequences of single-pulse trains from several hundred pulsars. The scientific outcomes from the programme will include determination of pulsar geometries, the location of the radio emission within the pulsar magnetosphere, the connection between the magnetosphere and the crust and core of the star, tighter constraints on the nature of the radio emission itself as well as interstellar medium studies. First results presented here include updated dispersion measures, 26 pulsars with Faraday rotation measures derived for the first time and a description of interesting emission phenomena observed thus far.

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.

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.

Motivated by the recent discoveries that six Ultraluminous X-ray Sources (ULXs) are powered by highly super-Eddington X-ray pulsars, we searched for additional pulsating ULX (PULX) candidates by identifying sources that exhibit long-term flux variability of at least an order of magnitude (a common feature seen in the 6 known PULXs, which may potentially be related to transitions to the propeller regime). Expanding on previous studies, we used the available fluxes from XMM-Newton, Swift and Chandra, along with carefully computed upper limits in cases of a non-detection, to construct long-term lightcurves for a sample of 296 ULXs selected from the XMM-Newton archive. Among these 296, we find 25 sources showing flux variability larger than a factor of 10, of which 17 show some evidence for (or are at least consistent with) exhibiting bi-modal flux distributions, as would be expected for sources undergoing propeller transitions. These sources are excellent candidates for continued monitoring programs to further test for this behaviour. There are 3 sources in our final sample with fluxes similar to NGC 5907 ULX1, currently the faintest known PULX, which would also be good targets for deeper observations with current facilities to search for pulsations. For the rest of the PULX candidates identified here, the next generation of X-ray telescopes (such as Athena) may be required to determine their nature owing to their lower peak fluxes.

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 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.

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.

SPIDER is a balloon-borne instrument designed to map the polarization of the millimeter-wave sky at large angular scales. SPIDER targets the B-mode signature of primordial gravitational waves in the cosmic microwave background (CMB), with a focus on mapping a large sky area with high fidelity at multiple frequencies. SPIDER's first longduration balloon (LDB) flight in January 2015 deployed a total of 2400 antenna-coupled Transition Edge Sensors (TESs) at 90 GHz and 150 GHz. In this work we review the design and in-flight performance of the SPIDER instrument, with a particular focus on the measured performance of the detectors and instrument in a space-like loading and radiation environment. SPIDER's second flight in December 2018 will incorporate payload upgrades and new receivers to map the sky at 285 GHz, providing valuable information for cleaning polarized dust emission from CMB maps.

We describe the construction and characterization of the 280 GHz bolometric focal plane units (FPUs) to be deployed on the second flight of the balloon-borne SPIDER instrument. These FPUs are vital to SPIDER's primary science goal of detecting or placing an upper limit on the amplitude of the primordial gravitational wave signature in the cosmic microwave background (CMB) by constraining the B-mode contamination in the CMB from Galactic dust emission. Each 280 GHz focal plane contains a 16 x 16 grid of corrugated silicon feedhorns coupled to an array of aluminum-manganese transition-edge sensor (TES) bolometers fabricated on 150 mm diameter substrates. In total, the three 280 GHz FPUs contain 1,530 polarization sensitive bolometers (765 spatial pixels) optimized for the low loading environment in flight and read out by time-division SQUID multiplexing. In this paper we describe the mechanical, thermal, and magnetic shielding architecture of the focal planes and present cryogenic measurements which characterize yield and the uniformity of several bolometer parameters. The assembled FPUs have high yields, with one array as high as 95% including defects from wiring and readout. We demonstrate high uniformity in device parameters, finding the median saturation power for each TES array to be ~3 pW at 300 mK with a less than 6% variation across each array at one standard deviation. These focal planes will be deployed alongside the 95 and 150 GHz telescopes in the SPIDER-2 instrument, slated to fly from McMurdo Station in Antarctica in December 2018.

We have used sensitive LOw Frequency ARray (LOFAR) observations of PSR B0809+74 at 15--62 MHz to study the anomalously intensive pulses, first reported by Ulyanov et al. (2006) at 18--30MHz. Similarly to Ulyanov et al., we found that the spectra of strong pulses consist of distinct bright patches. Moreover, these spectral patches were spotted to drift upwards in frequency over the course of several pulse sequences. We established that this drift is not pulsar-intrinsic, but is caused by the broadband ~20 second-long enhancements of recorded signal, which influenced the dispersed tracks of several pulses at once. We speculate on the cause of such enhancements (i.e. propagation or telescope-related) and the ramifications they bring to the single-pulse studies at the very low radio frequencies. Depending on the origin, the phenomenon may also affect the analysis of highly dispersed single pulses at higher radio frequencies, e.g. Fast Radio Bursts.

We present a new upper limit on CMB circular polarization from the 2015 flight of SPIDER, a balloon-borne telescope designed to search for $B$-mode linear polarization from cosmic inflation. Although the level of circular polarization in the CMB is predicted to be very small, experimental limits provide a valuable test of the underlying models. By exploiting the non-zero circular-to-linear polarization coupling of the HWP polarization modulators, data from SPIDER's 2015 Antarctic flight provide a constraint on Stokes $V$ at 95 and 150 GHz from $33<\ell<307$. No other limits exist over this full range of angular scales, and SPIDER improves upon the previous limit by several orders of magnitude, providing 95% C.L. constraints on $\ell (\ell+1)C_{\ell}^{VV}/(2\pi)$ ranging from 141 $\mu K ^2$ to 255 $\mu K ^2$ at 150 GHz for a thermal CMB spectrum. As linear CMB polarization experiments become increasingly sensitive, the techniques described in this paper can be applied to obtain even stronger constraints on circular polarization.

We analyze sunspots rotation and magnetic transients in NOAA AR 11429 during two X-class (X5.4 and X1.3) flares using the data from the Helioseismic and Magnetic Imager on board the \emph{Solar Dynamics Observatory}. A large leading sunspot with positive magnetic polarity rotated counterclockwise. As expected, the rotation was significantly affected by the two flares. The magnetic transients induced by the flares were clearly evident in the sunspots with negative polarity. They were moving across the sunspots with speed of order $3-7\ \rm km \ s^{-1}$. Furthermore, the trend of magnetic flux evolution of these sunspots exhibited changes associated with the flares. These results may shed light on the understanding of the evolution of sunspots.

We describe 280 GHz bolometric detector arrays that instrument the balloon-borne polarimeter SPIDER. A primary science goal of SPIDER is to measure the large-scale B-mode polarization of the cosmic microwave background in search of the cosmic-inflation, gravitational-wave signature. 280 GHz channels aid this science goal by constraining the level of B-mode contamination from galactic dust emission. We present the focal plane unit design, which consists of a 16$\times$16 array of conical, corrugated feedhorns coupled to a monolithic detector array fabricated on a 150 mm diameter silicon wafer. Detector arrays are capable of polarimetric sensing via waveguide probe-coupling to a multiplexed array of transition-edge-sensor (TES) bolometers. The SPIDER receiver has three focal plane units at 280 GHz, which in total contains 765 spatial pixels and 1,530 polarization sensitive bolometers. By fabrication and measurement of single feedhorns, we demonstrate 14.7$^{\circ}$ FHWM Gaussian-shaped beams with $<$1% ellipticity in a 30% fractional bandwidth centered at 280 GHz. We present electromagnetic simulations of the detection circuit, which show 94% band-averaged, single-polarization coupling efficiency, 3% reflection and 3% radiative loss. Lastly, we demonstrate a low thermal conductance bolometer, which is well-described by a simple TES model and exhibits an electrical noise equivalent power (NEP) = 2.6 $\times$ 10$^{-17}$ W/$\sqrt{\mathrm{Hz}}$, consistent with the phonon noise prediction.