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.

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 recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

In this paper we present the science potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies of strongly magnetized objects. We will focus on the physics and astrophysics of strongly magnetized objects, namely magnetars, accreting X-ray pulsars, and rotation powered pulsars. We also discuss the science potential of eXTP for QED studies. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Sciences, the eXTP mission is expected to be launched in the mid 2020s.

Aql X--1 is one of the most prolific low mass X-ray binary transients (LMXBTs) showing outbursts almost annually. We present the results of our spectral analyses of RXTE/PCA observations of the 2000 and the 2011 outbursts. We investigate the spectral changes related to the changing disk-magnetosphere interaction modes of Aql X--1. The X-ray light curves of the outbursts of LMXBTs typically show phases of fast rise and exponential decay. The decay phase shows a "knee" where the flux goes from the slow decay to the rapid decay stage. We assume that the rapid decay corresponds to a weak propeller stage at which a fraction of the inflowing matter in the disk accretes onto the star. We introduce a novel method for inferring, from the light curve, the fraction of the inflowing matter in the disk that accretes onto the NS depending on the fastness parameter. We determine the fastness parameter range within which the transition from the accretion to the partial propeller stage is realized. This fastness parameter range is a measure of the scale-height of the disk in units of the inner disk radius. We applied the method to a sample of outbursts of Aql X--1 with different maximum flux and duration times. We show that different outbursts with different maximum luminosity and duration follow a similar path in the parameter space of accreted/inflowing mass flux fraction versus fastness parameter.

We present the X-ray light curves of the last two outbursts --2014 & 2016-- of the well known accreting millisecond X-ray pulsar (AMXP) Aquila X-1 using the monitor of all sky X-ray image (MAXI) observations in the $2-20$ keV band. After calibrating the (MAXI) count rates to the all-sky monitor (ASM) level, we report that the 2016 outburst is the most energetic event of Aql X-1, ever observed from this source. We show that 2016 outburst is a member of the long-high class according to the classification presented by G\"ung\"or et al. with $\sim68$ cnt/s maximum flux and $\sim60$ days duration time and the previous outburst, 2014, belongs to the short-low class with $\sim 25$ cnt/s maximum flux and $\sim 30$ days duration time. In order to understand differences between outbursts, we investigate the possible dependence of the peak intensity to the quiescent duration leading to the outburst and find that the outbursts following longer quiescent episodes tend to reach higher peak energetic.

The recent detection of pulsations from the ultra luminous X-ray source (ULX) NuSTAR J095551+6940.8 in M82 by Bachetti et al. indicates that the object is an accreting neutron star in a high mass X-ray binary (HMXB) system. The super-Eddington luminosity of the object implies that the magnetic field is sufficiently strong to suppress the scattering cross-section unless its beam is viewed at a favourable angle. We show that the torque equilibrium condition for the pulsar indicates the dipole magnetic field of the neutron star is $6.7 \times 10^{13}$ G, two orders of magnitude higher than that estimated by Bachetti et al., and further point to the possibility that even stronger magnetic fields could well be in the higher multipoles. This supports the recent view that magnetars descent from HMXBs if the magnetic field decays an order of magnitude during the process of transition.

Neutron stars (NSs) are thought to be excellent laboratories for determining the equation of state (EoS) of cold dense matter. Their strong gravity suggests that they can also be used to constrain gravity models. The mass and radius (M-R) of a NS both depend on the choice of EoS and gravity, meaning that NSs cannot be simultaneously good laboratories for both of these questions. A measurement of M-R would constrain the less well known physics input. The assumption that M-R measurements can be used to constrain EoS-presumes general relativity (GR) is the ultimate model of gravity in the classical regime. We calculate the radial profile of compactness and curvature (square root of the full contraction of the Weyl tensor) within a NS and determine the domain not probed by the Solar System tests of GR. We find that, except for a tiny sphere of radius less than a millimeter at the center, the curvature is several orders of magnitude above the values present in Solar System tests. The compactness is beyond the solar surface value for r>10 m, and increases by 5 orders of magnitude towards the surface. With the density being only an order of magnitude higher than that probed by nuclear scattering experiments, our results suggest that the employment of GR as the theory of gravity describing the hydrostatic equilibrium of NSs is a rather remarkable extrapolation from the regime of tested validity, as opposed to that of EoS models. Our larger ignorance of gravity within NSs suggests that a measurement of M-R constrains gravity rather than EoS, and given that EoS has yet to be determined by nucleon scattering experiments, M-R measurements cannot tightly constrain the gravity models either. Near the surface the curvature and compactness attain their largest values, while EoS in this region is fairly well known. This renders the crust as the best site to look for deviations from GR.

We present a broad classification of all outbursts detected with the All-Sky Monitor (ASM) on the Rossi X-Ray Timing Explorer (RXTE) and the Monitor of All Sky X-Ray Image (MAXI) of Aql X-1. We identify three types of outbursts; long-high, medium-low, and short-low, based on the duration and maximum flux. We analyse the trends in the "phase-space" of flux-derivative versus flux to demonstrate the differences in the three identified outburst types. We present a spectral analysis of the observations of Aql X-1 performed by the Proportional Counter Array (PCA) onboard RXTE during the 2000 and 2011 outbursts of the long-high class and the 2010 outburst of the medium-low class. We model the source spectrum with a hybrid thermal/non-thermal hot plasma emission model (EQPAIR in XSPEC, Coppi 2000) together with a Gaussian component to model the Fe K_alpha emission line. We construct time histories of the source flux, the optical depth of the corona (tau), the seed photon temperature (kT_bb) and the hard state compactness (l_h) for these three outbursts. We show that the physical parameters of either classes reach the same values throughout the outbursts, the only difference being the maximum flux. We discuss our results in the terms of modes of interaction of the star with the disc and size of the disc kept hot by irradiation. We conclude that irradiation is the dominant physical process leading to the different classes of outbursts.

We investigate the effect of a strong magnetic field on the structure of neutron stars in a model with perturbative $f(R)$ gravity. The effect of an interior strong magnetic field of about $10^{17 \sim 18}$ G on the equation of state is derived in the context of a quantum hadrodynamics (QHD) model. We solve the modified spherically symmetric hydrostatic equilibrium equations derived for a gravity model with $f(R)=R+\alpha R^2$. Effects of both the finite magnetic field and the modified gravity are detailed for various values of the magnetic field and the perturbation parameter $\alpha$ along with a discussion of their physical implications. We show that there exists a parameter space of the modified gravity and the magnetic field strength, in which even a soft equation of state can accommodate a large ($> 2$ M$_\odot$) maximum neutron star mass through the modified mass-radius relation.