The radiation-loud and radiation-quiet anti-glitches in magnetar 1E 2259+586 are discussed together. It is shown that the 2019 radiation-quiet anti-glitch in 1E 2259+586 is not constraining, because its amplitude is 2.5 smaller than the previous radiation-loud anti-glitch. We note that magnetars with high X-ray luminosity and low rotational energy loss rate are more likely to show anti-glitches. Anti-glitchers and radio-emitters of magnetars occupy opposite corner on the luminosity vs. rotational energy loss rate diagram. By introducing the X-ray luminosity parameter and the activity parameter, we try to characterize this trend quantitatively. A general picture of magnetar anti-glitch is given.

The multi-messenger observation of compact binary coalescence promises great scientific treasure. However, a synthetic observation from both gravitational wave and electromagnetic channels remains challenging. Relying on the day-to-week long macronova emission, GW170817 remains the only event with successful electromagnetic followup. In this manuscript, we explore the possibility of using the early stage X-ray afterglow to search for the electromagnetic counterpart of gravitational wave events. Two algorithms, the sequential observation and the local optimization are considered and applied to three simulated events. We consider the proposed Einstein probe as a candidate X-ray telescope. Benefiting from the large field of view and high sensitivity, we find that the sequential observation algorithm not only is easy to implement, but also promises a good chance of actual detection.

The outburst and spin-down glitch of magnetars are modeled from the magnetospheric point of view. We try to answer the following four questions: (1) Which pulsar on the period and peirod-derivative diagram are more likely to show magnetar outburst? (2) Which outburst will make the glitch that triggered the outburst to become a spin-down glitch? (3) Can we model the outburst and spin-down glitch in PSR J1119$-$6127 simultaneously? (4) Why the torque variation is delayed compared with the peak of the X-ray luminosity in 1E 1048.1$-$5937 and PSR J1119$-$6127? It is found that both the global and local twisted magnetic field will affect the radiation and timing behaviors of magnetars. Especially, the delay of torque variations may due to the combined effect of increasing twist in the j-bundle and untwisting of the global magnetosphere. A toy model is build for magnetar outburst and torque variations. It can catch the general trend of magnetar outburst: decaying flux, shrinking hot spot, and torque variations.

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

Recently, a 16-day periodicity in one fast radio burst is reported. We propose that this 16-day periodicity may be due to forced precession of the neutron star by a fallback disk. When the rotation axis is misaligned with the normal direction of the disk plane, the neutron star will precess. The eccentricity of the neutron star may be due to rotation or strong magnetic field etc. We found that the 16-day period may be understood using typical masses of the fallback disk. Polarization observations and information about the neutron star rotation period may help to discriminate different models. The possible precession observations in pulsars, magnetars, and fast radio bursts may be understood together using forced precession by a fallback disk.

Assuming the wind-fed accretion magnetars in long period X-ray pulsars, we calculated the rotational evolution of the neutron stars. Our calculations considered the effects of the magnetic field decay in magnetars. The results show that wind-fed accretion magnetars can evolve to the long period X-ray pulsars with a spin period much longer than 1000 s. The spin-down trend observed in 4U 2206+54 like sources is expected when the young X-ray binary systems are on the way to their equilibrium period. Detailed calculations showed that its spin-down may be affected by accretion with outflows or accretion while spin-down. Due to the magnetic field decay in magnetars, wind-fed accretion magnetars will have a decreasing equilibrium period for a constant mass accretion rate. For 2S 0114+65, the spin-up rate due to magnetic field decay is one order of magnitude smaller than the observations. The spin-up rate of 2S 0114+65 may be attributed to the formation of a transient disk during wind accretion. The slowest X-ray pulsar AX J1910.7+0917 would be a link source between 4U 2206+54 and 2S 0114+65.

The magnetic field of magnetars may be twisted compared with that of normal pulsars. Previous works mainly discussed magnetic energy release in the closed field line regions of magnetars. For a twisted magnetic field, the field lines will inflate in the radial direction. Similar to normal pulsars, the idea of light cylinder radius is introduced. More field lines will cross the light cylinder and become open for a twisted magnetic field. Therefore, magnetars may have a large polar cap, which may correspond to the hot spot during outburst. Particle flow in the open field line regions will result in the untwisting of the magnetic field. Magnetic energy release in the open field line regions can be calculated. The model calculations can catch the general trend of magnetar outburst: decreasing X-ray luminosity, shrinking hot spot etc. For magnetic energy release in the open field line regions, the geometry will be the same for different outburst in one magnetar.

Using the archive data from the Rossi X-ray Timing Explorer ({\sl RXTE}), we have studied the evolution of the X-ray profile of the Crab pulsar in a time span of 11 years. The X-ray profiles, as characterized by a few parameters, changed slightly but significantly in these years: the separation of the two peaks increased with a rate $0.88\pm0.20\,\textordmasculine$\,per century, the flux ratio of the second pulse to the first pulse decreased with $(3.64\pm0.86)\times10^{-2}$\,per century, and the pulse widths of the two pulses descended with $1.44\pm0.15\,\textordmasculine$\, and $1.09\pm0.73\,\textordmasculine$\,per century, respectively. The evolutionary trends of the above parameters are similar to the radio results, but the values are different. We briefly discussed the constraints of these X-ray properties on the geometry of the emission region of this pulsar.

We apply theoretical spin-down models of magnetospheric evolution and magnetic field decay to simulate the possible evolution of PSR J0250+5854, which is the slowest-spinning radio pulsar detected to date. Considering the alignment of inclination angle in a 3-D magnetosphere, it is possible that PSR J0250+5854 has a high magnetic field comparable with magnetars or/and high magnetic field pulsars, if a small inclination angle is considered. Our calculations show that similar long-period pulsars tend to have a relatively low period derivative in this case. In another case of magnetic field decay, calculations also show a possible connection between PSR J0250+5854 and high dipole-magnetic field magnetars. The evolutionary path indicates a relatively high spin-down rate for similar long-period pulsars.

The magnetic dipole field of ultraluminous X-ray pulsars may not be very high. However, it is too early to say that they are not magnetars. The existence of low magnetic field magnetars should be taken into consideration.

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.

Possible manifestations of accreting magnetars are discussed. It is shown that the four ultra-luminous X-ray pulsars can be understood in the accreting low magnetic field magnetar scenario. The NGC300 ULX1 pulsar may have a higher dipole magnetic field than other sources. General constraint on their mass accretion rate confirmed their super-Eddington nature. Lower limits on their beaming factor are obtained. They do not seem to have strong beaming. The duty cycle of the ULX burst state can also be constrained by their timing observations. ULX pulsars may be in accretion equilibrium in the long run. During the outburst, they will spin up, and run from the previous equilibrium state to the new equilibrium state. It is proposed that the slowest puslation X-ray pulsar AX J1910.7+0917 may be an accreting magnetar with a low mass accretion rate. ULX pulsars, slow pulsation X-ray pulsars may all be accreting magnetars with different accretion rates. Seven possible signatures of an accreting magnetar are summarized.

The gamma-ray burst GRB 170817A associated with GW170817 is subluminous and subenergetic compared with other typical short GRBs. It may be due to a relativistic jet viewed off-axis, or a structured jet, or cocoon emission. Giant flares from magnetars may possibly be ruled out. However, the luminosity and energetics of GRB 170817A is coincident with that of magnetar giant flares. After the coalescence of the binary neutron star, a hypermassive neutron star may be formed. The hypermassive neutron star may have magnetar-strength magnetic field. During the collapse of the hypermassive neutron star, the magnetic field energy will also be released. This giant-flare-like event may explain the the luminosity and energetics of GRB 170817A. Bursts with similar luminosity and energetics are expected in future neutron star-neutron star or neutron star-black hole mergers.

Calvera is a $59\,\mathrm{ms}$ isolated pulsar, being unique due to its non-detection in radio, optical and gamma-rays but the purely thermal emission in soft X-rays. It is suggested that Calvera could be an ordinary middle-aged pulsar with significant magnetospheric activity at a large distance. Alternatively, it is proposed in this paper that Calvera is a low-mass strangeon star with inactive magnetosphere (dead). In this scenario, we jointly fit the spectra obtained by the {\it XMM-Newton} Observatory and the {\it Chandra} X-ray Observatory with the strangeon star atmosphere model. The spectral model is successful in explaining the radiation properties of Calvera and X-ray Dim Isolated Neutron Stars, both showing similar observation features. Within the dead pulsar picture, Calvera might be of high temperature at $0.67\,\mathrm{keV}$, possessing a small stellar radius $R\la4\,\mathrm{km}$ and a presumably small magnetic field $B\la10^{11}\,\mathrm{G}$ and is probably braked by the fall-back disk accretion. Future advanced facilities may provide unique opportunities to know the real nature of Calvera.

We study the X-ray phases of the Crab pulsar utilizing the 11-year observations from the Rossi X-ray Timing Explorer, 6-year radio observations from the Nanshan Telescope, and the ephemeris from Jodrell Bank Observatory. It is found that the X-ray phases in different energy bands and the radio phases from Nanshan Telescope show similar behaviors, including long-time evolution and short-time variations. Such strong correlations between the X-ray and radio phases imply that the radio and X-ray timing noises are both generated from the pulsar spin that cannot be well described by the the monthly ephemeris from the Jodrell Bank observatory. When using the Nanshan phases as references to study the X-ray timing noise, it has a significantly smaller variation amplitude and shows no long-time evolution, with a change rate of $(-1.1\pm1.1)\times10^{-7}$ periods per day. These results show that the distance of the X-ray and radio emission regions on the Crab pulsar has no detectable secular change, and it is unlikely that timing-noises resulted from any unique physical processes in the radio or X-ray emitting regions. The similar behaviors of the X-ray and radio timing noises also imply that the variation of the interstellar medium is not the origin of the Crab pulsar's timing noises, which is consistent with the results obtained from the multi-frequency radio observations of PSR B1540-06.

Pulsar death line can be defined in $P-\dot{P}$ diagram. Traditionally, radio-loud pulsars are supposed to locate above the death line where is the radio-loud region. With the development of observational equipment, the observational properties of the neutron star are remarkably diverse. In $P-\dot{P}$ diagram, some of the special sources are radio-quiet but lie above the death line. From the definition of the pulsar death line, different equation of states for neutron star or strange star results in different death lines. We discuss the influence of the equation of state on the pulsar death line and the possible link between different neutron star groups. The results show that central compact objects would be the small mass of self-bound strange stars, and rotating radio transients might be old pulsars on the verge of death. We suggest that PSR J2144-3933 is likely to be a large mass pulsar, which would be larger than $2.0\rm{M_{\odot}}$. Multiple observational facts would help us to reveal the nature of pulsars.

The wind braking model and its applications of magnetars are discussed. The decreasing torque of magnetars during outbursts, anti-glitch, and anti-correlations between radiation and timing are understandable in the wind braking model. Recent timing observations of magnetars are also consistent with the previous modeling. A magnetism-powered wind nebula and a braking index smaller than three are the two predictions. Besides isolated magnetars, there may also be accreting magnetars in binary systems and magnetars accreting from fallback disks. Observationally, ultra-luminous X-ray pulsars may be accreting magnetars, while super-slow magnetars may be magnetars with fallback disks in the past. Many works are needed for both isolated magnetars and accreting magnetars.

The coupled evolution of pulsar rotation and inclination angle in the wind braking model is calculated. The oblique pulsar tends to align. The pulsar alignment will affect its spin-down behavior. As a pulsar evolves from the magneto-dipole radiation dominated case to the particle wind dominated case, the braking index will first increase and then decrease. In the early time, the braking index may be larger than 3. And during the following long time, the braking index will be always smaller than 3. The minimum braking index is about one. This can explain the existence of high braking index larger than 3, and low braking index of pulsars simultaneously. The pulsar braking index is expected to evolve from larger than three to about one. A general trend is that the pulsar braking index will evolve from the Crab-like case to the Vela-like case.

X-ray dim isolated neutron stars (XDINSs) are characterized by Planckian spectra in X-ray bands, but show optical/ultraviolet(UV) excesses which are the measured photometry exceeding that is extrapolated from X-ray spectra. To solve this problem, a radiative model of bremsstrahlung emission from a plasma atmosphere is established in the regime of strangeon star. A strangeon star atmosphere could simply be regarded as the upper layer of a normal neutron star. This plasma atmosphere, formed and maintained by the ISM-accreted matter due to the so-called strangeness barrier, is supposed to be of two-temperature. All the seven XDINS spectra could be well fitted by the radiative model, from optical/UV to X-ray bands. The fitted radiation radii of XDINSs are from 7 to 13km, while the modelled electron temperatures are between 50 and 250eV, except RX J0806.4-4123 with a radiation radius ~3.5km, indicating that this source could be a low-mass strangeon star candidate. This strangeon star model could further be tested by soft X-ray polarimetry, such as the Lightweight Asymmetry and Magnetism Probe which is expected to work on Chinese space station around 2020.

Magnetars may have strong surface dipole field. Observationally, two magnetars may have passive fallback disks. In the presence of a fallback disk, the rotational evolution of magnetars may be changed. In the self-similar fallback disk model, it is found that: (1) When the disk mass is significantly smaller than $10^{-6} \,\rm M_{\odot}$, the magnetar is unaffected by the fallback disk and it will be a normal magnetar. (2) When the disk mass is large, but the magnetar's surface dipole field is about or below $10^{14} \,\rm G$, the magnetar will also be a normal magnetar. A magnetar plus a passive fallback disk system is expected. This may correspond to the observations of magnetars 4U 0142$+$61, and 1E 2259$+$586. (3) When the disk mass is large, and the magnetar's surface dipole field is as high as $4\times 10^{15} \,\rm G$, the magnetar will evolve from the ejector phase to the propeller phase, and then enter rotational equilibrium. The magnetar will be slowed down quickly in the propeller phase. The final rotational period can be as high $2\times 10^4 \,\rm s$. This may correspond to the super-slow magnetar in the supernova remnant RCW 103. Therefore, the three kinds of magnetars can be understood in a unified way.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

We propose that the magnetosphere reconfiguration induced by magnetic twists in the closed field line region can account for the mode-switching of intermittent pulsars. We carefully investigate the properties of axisymmetric force-free pulsar magnetospheres with magnetic twists in closed field line region around the polar caps. The magnetosphere with twisted closed lines leads to enhanced spin-down rates. The enhancement in spin-down rate depends on the size of region with twisted closed lines. Typically, it is increased by a factor of $\sim2$, which is consistent with the intermittent pulsars' spin down behavior during the `off' and `on' states. We find there is a threshold of maximal twist angle $\Delta\phi_{\rm thres}\sim1$. The magnetosphere is stable only if the closed line twist angle is less than $\Delta\phi_{\rm thres}$. Beyond this value, the magnetosphere becomes unstable and gets untwisted. The spin-down rate would reduce to its off-state value. The quasi-periodicity in spin-down rate change can be explained by long-term activities in star's crust and the untwisting induced by MHD instability. The estimated duration time of on-state is about one week, consistent with observations. Due to the MHD instability, there exists an upper limit for the spin down ratio ($f\sim3$) between the on-state and the off-state, if the Y-point remains at the light cylinder.

Recently, the magnetar Swift J1834.9$-$0846 is reported to have a possible wind nebula. It is shown that both the magnetar and its wind nebula are understandable in the wind braking scenario. The magnetar's rotational energy loss rate is not enough to power the particle luminosity. The required particle luminosity should be about $10^{36} \,\rm erg \,s^{-1}$ to $10^{38} \,\rm erg \,s^{-1}$. It is obtained in three different approaches: considering wind braking of Swift J1834.9$-$0846; the spectral and spatial observations of the wind nebula; and an empirical upper bound on wind nebula X-ray luminosity. The nebula magnetic field is about $10^{-4} \,\rm G$. The possible wind nebula of Swift J1834.9$-$0846 should be a magnetar wind nebula. It is powered by the magnetic energy release of the magnetar.

Ultra-luminous X-ray sources are usually believed to be black holes with mass about $10^{2-3}M_{\odot}$. However, the recent discovery of NuSTAR J095551+6940.8 in M82 by Bachetti et al. shows that it holds the spin period $P=1.37\rm\,s$ and period {\bf derivative} $\dot{P}\approx-2\times10^{-10}\rm\,s\,s^{-1}$, which provides a strong evidence that some ultra-luminous X-ray sources could be neutron stars. We obtain that the source may be an evolved magetar, according to our simulation by employing the model of accretion induced the polar magnetic field decay and standard spin-up torque of an accreting neutron star. The results show that NuSTAR J095551+6940.8 is still in the spin-up process, and the polar magnetic field decays to about $4.5\times10^{12}\rm\,G$ after accreting $\sim 10^{-2.5}$\ms, while the strong magnetic field exists in the out-polar region, which could be responsible for the observed low field magnetar. The ultra luminosity of the source can be explained by the beaming effort and two kinds of accretion--radial random accretion and disk accretion. Since the birth rate of magnetars is about ten percent of the normal neutron stars, we guess that several ultra-luminous X-ray sources should share the similar properties to that of NuSTAR J095551+6940.8.

Eight pulsars have low braking indices, which challenge the magnetic dipole braking of pulsars. 222 pulsars and 15 magnetars have abnormal distribution of frequency second derivatives, which also make contradiction with classical understanding. How neutron star magnetospheric activities affect these two phenomena are investigated by using the wind braking model of pulsars. It is based on the observational evidence that pulsar timing is correlated with emission and both aspects reflect the magnetospheric activities. Fluctuations are unavoidable for a physical neutron star magnetosphere. Young pulsars have meaningful braking indices, while old pulsars' and magnetars' fluctuation item dominates their frequency second derivatives. It can explain both the braking index and frequency second derivative of pulsars uniformly. The braking indices of eight pulsars are the combined effect of magnetic dipole radiation and particle wind. During the life time of a pulsar, its braking index will evolve from three to one. Pulsars with low braking index may put strong constraint on the particle acceleration process in the neutron star magnetosphere. The effect of pulsar death should be considered during the long term rotational evolution of pulsars. An equation like the Langevin equation for Brownian motion was derived for pulsar spin-down. The fluctuation in the neutron star magnetosphere can be either periodic or random, which result in anomalous frequency second derivative and they have similar results. The magnetospheric activities of magnetars are always stronger than those of normal pulsars.

Pulsars are good clocks in the universe. One fundamental question is that why they are good clocks? This is related to the braking mechanism of pulsars. Nowadays pulsar timing is done with unprecedented accuracy. More pulsars have braking indices measured. The period derivative of intermittent pulsars and magnetars can vary by a factor of several. However, during pulsar studies, the magnetic dipole braking in vacuum is still often assumed. It is shown that the fundamental assumption of magnetic dipole braking (vacuum condition) does not exist and it is not consistent with the observations. The physical torque must consider the presence of the pulsar magnetosphere. Among various efforts, the wind braking model can explain many observations of pulsars and magnetars in a unified way. It is also consistent with the up-to-date observations. It is time for a paradigm shift in pulsar studies: from magnetic dipole braking to wind braking. As one alternative to the magnetospheric model, the fallback disk model is also discussed.

PSR B0540-69 has braking index measurement in its persistent state: n=2.129 \pm 0.012. Recently, it is reported to have spin-down state changes: a suddenly 36 % increase in the spin-down rate. Combining the persistent state braking index measurement and different spin-down states, PSR B0540-69 is more powerful than intermittent pulsars in constraining pulsar spin-down models. The pulsar wind model is applied to explain the variable timing behavior of PSR B0540-69. The persistent state braking index of PSR B0540-69 is the combined effect of magnetic dipole radiation and particle wind. The particle density reflects the magnetospheric activity in real-time and may be responsible for the changing spin-down behavior. Corresponding to the 36 % increase in the spin-down rate of PSR B0540-69, the relative increase in the particle density is 88 % in the vacuum gap model. And the model calculated braking index in the new state is n=1.79. Future braking index observation of PSR B0540-69 in the new spin-down state will be very powerful in distinguishing between different pulsar spin-down models and different particle acceleration models in the wind braking scenario. The variable timing behavior of PSR J1846-0258 is also understandable in the pulsar wind model.

Magnetars are a special kind of neutron stars. There may also be accreting magnetars. From the studies of isolated magnetars, it is known that a neutron star with a strong dipole field only is not a magnetar. Super-slow X-ray pulsars may just be accreting high magnetic field neutron stars. The ultra-luminous X-ray pulsar NuSTAR J095551+6940.8 may be an accreting magnetar. It may be an accreting low magnetic field magnetar with multipole field of 10^14 G and dipole field of 10^12 G. This point of view is consistent with the study of isolated magnetars. An ultra-luminous X-ray pulsar phase in the binary evolution may result in massive millisecond pulsars.

The detection of 1.37$\, $s pulsations from NuSTAR J095551+6940.8 implies the existence of an accreting pulsar, which challenges the conventional understanding of ultraluminous X-ray source. This kind of sources are proposed to be massive pulsars in this paper. Considering the general relativistic effect, stronger gravity of massive pulsars could lead to a larger maximum luminosity, scaled as the Eddington luminosity. The pseudo-Newtonian potential is employed to simulate the gravitational field in general relativity, and the Eddington luminosity is calculated for self-bound stars (quark star and quark-cluster star) and for the Tolman IV solution. It is found that, for a massive pulsar with radius close to the Schwarszchild radius, the Eddington luminosity could be as high as $2\times10^{39} \, {\rm erg\, s}^{-1}$. It is able to account for the X-ray luminosity of NuSTAR J095551+6940.8 with reasonable beaming factor. It is also suggested that massive pulsar-like compact stars could form via this super-Eddington phase of ultraluminous X-ray source.

The pulsar wind model is updated by considering the effect of particle density and pulsar death. It can describe both the short term and long term rotational evolution of pulsars consistently. It is applied to model the rotational evolution of the Crab pulsar. The pulsar is spun down by a combination of magnetic dipole radiation and particle wind. The parameters of the Crab pulsar, including magnetic field, inclination angle, and particle density are calculated. The primary particle density in acceleration region is about 10^3 times the Goldreich-Julian charge density. The lower braking index between glitches is due to a larger outflowing particle density. This may be glitch induced magnetospheric activities in normal pulsars. Evolution of braking index and the Crab pulsar in P-Pdot diagram are calculated. The Crab pulsar will evolve from magnetic dipole radiation dominated case towards particle wind dominated case. Considering the effect of pulsar "death", the Crab pulsar (and other normal pulsars) will not evolve to the cluster of magnetars but downwards to the death valley. Different acceleration models are also considered. Applications to other sources are also discussed, including pulsars with braking index measured, and the magnetar population.

The basic concepts of wind braking of magnetars, and its applications to various timing events of magnetars are outlined. The case of the Galactic centre magnetar SGR J1745-2900 is especially discussed. Two predictions for SGR J1745-2900 are made (the period derivative will first increase then decrease; and possible outbursts in the near future).

One ultraluminous X-ray source in M82 is identified as an accreting neutron star recently (named as NuSTAR J095551+6940.8). It has a super-Eddington luminosity and is spinning up. For an aged magnetar, it is more likely to be a low magnetic field magnetar. An accreting low magnetic field magnetar may explain both the super-Eddington luminosity and the rotational behavior of this source. Considering the effect of beaming, the spin-up rate is understandable using the traditional form of accretion torque. The transient nature, spectral properties of M82 X-2 are discussed. The theoretical period range of accreting magnetar is provided. Three observational appearances of accreting magnetars are summarized.

SGR J1745$-$2900 is a magnetar near the Galactic center. X-ray observations of this source found a decreasing X-ray luminosity accompanied by an enhanced spindown rate. This negative correlation between X-ray luminosity and spindown rate is hard to understand. The wind braking model of magnetars is employed to explain this puzzling spindown behavior. During the release of magnetic energy of magnetars, a system of particles may be generated. Some of these particles remain trapped in the magnetosphere and may contribute to the X-ray luminosity. The rest of the particles can flow out and take away the rotational energy of the central neutron star. A smaller polar cap angle will cause the decrease of X-ray luminosity and enhanced spindown rate of SGR J1745$-$2900. This magnetar is expected to have a maximum spindown rate shortly.

Intermittent pulsars are part-time radio pulsars. They have higher slow down rate in the on state (radio-loud) than in the off state (radio-quiet). This gives the evidence that particle wind may play an important role in pulsar spindown. The effect of particle acceleration is included in modeling the rotational energy loss rate of the neutron star. Applying the pulsar wind model to the three intermittent pulsars (PSR B1931+24, PSR J1841-0500, and PSR J1832+0029), their magnetic field and inclination angle are calculated simultaneously. The theoretical braking indices of intermittent pulsars are also given. In the pulsar wind model, the density of the particle wind can always be the Goldreich-Julian density. This may ensure that different on states of intermittent pulsars are stable. The duty cycle of particle wind can be determined from timing observations. It is consistent with the duty cycle of the on state. Inclination angle and braking index observations of intermittent pulsars may help to test different models of particle acceleration. At present, the inverse Compton scattering induced space charge limited flow with field saturation model can be ruled out.

Glitch (sudden spinup) is a common phenomenon in pulsar observations. However, the physical mechanism of glitch is still a matter of debate because it depends on the puzzle of pulsar's inner structure, i.e., the equation of state of dense matter. Some pulsars (e.g., Vela-like) show large glitches ({\Delta}{\nu}/{\nu}~10^-6) but release negligible energy, whereas the large glitches of AXPs/SGRs (anomalous X-ray pulsars/soft gamma repeaters) are usually (but not always) accompanied with detectable energy releases manifesting as X-ray bursts or outbursts. We try to understand this aspect of glitches in a starquake model of solid quark stars. There are actually two kinds of glitches in this scenario: bulk-invariable (Type I) and bulk-variable (Type II) ones. The total stellar volume changes (and then energy releases) significantly for the latter but not for the former. Therefore, glitches accompanied with X-ray bursts (e.g., that of AXP/SGRs) could originate from Type II starquakes induced probably by accretion, while the others without evident energy release (e.g., that of Vela pulsar) would be the result of Type I starquakes due to, simply, a change of stellar ellipticity.

Magnetars and many of the magnetar-related objects are summarized together and discussed. It is shown that there is an abuse of language in the use of "magnetar". Anomalous X-ray pulsars and soft gamma-ray repeaters are well-known magnetar candidates. The current so called anti-magnetar (for central compact objects), accreting magnetar (for superslow X-ray pulsars in high mass X-ray binaries), and millisecond magnetar (for the central engine of some gamma-ray bursts), they may not be real magnetars in present understandings. Their observational behaviors are not caused by the magnetic energy. Many of them are just neutron stars with strong surface dipole field. A neutron star plus strong dipole field is not a magnetar. The characteristic parameters of the neutron stars for the central engine of some gamma-ray bursts are atypical from the neutron stars in the Galaxy. Possible signature of magnetic activities in accreting systems are discussed, including repeated bursts and a hard X-ray tail. China's future hard X-ray modulation telescope may contribute to finding some magnetic activities in accreting neutron star systems.

Hard X-rays above 10 keV are detected from several anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs), and different models have been proposed to explain the physical origin within the frame of either magnetar model or fallback disk system. Using data from $\it Suzaku$ and INTEGRAL, we study the soft and hard X-ray spectra of four AXPs/SGRs, 1RXS J170849-400910, 1E 1547.0-5408, SGR 1806-20 and SGR 0501+4516. It is found that the spectra could be well reproduced by bulk-motion Comptonization (BMC) process as was first suggested by Tr$\"u$mper et al., showing that the accretion scenario could be compatible with X-ray emission of AXPs/SGRs. HXMT simulations for BMC model show that the spectra would have discrepancies from power-law, especially the cutoff at $\sim$ 200 keV. Thus future observations are promising to distinguish different models for the hard tail and may help us understand the nature of AXPs/SGRs.

The very small braking index of PSR J1734-3333, $n=0.9\pm0.2$, challenges the current theories of braking mechanisms in pulsars. We present a possible interpretation that this pulsar is surrounded by a fall-back disk and braked by it. A modified braking torque is proposed based on the competition between the magnetic energy density of a pulsar and the kinetic energy density of a fall-back disk. With this torque, a self-similar disk can fit all the observed parameters of PSR J1734-3333 with natural initial parameters. In this regime, the star will evolve to the region having anomalous X-ray pulsars and soft gamma repeaters in the $P-\dot{P}$ diagram in about 20000 years and stay there for a very long time. The mass of the disk around PSR J1734-3333 in our model is about $10M_{\oplus}$, similar to the observed mass of the disk around AXP 4U 0142+61.

The anti-glitch of magnetar 1E 2259+586 is analyzed theoretically. An enhanced particle wind during the observational interval will taken away additional rotational energy of the neutron star. This will result in a net spin-down of the magnetar, i.e., an anti-glitch. In the wind braking scenario of anti-glitch, there are several predictions: (1) A radiative event will always accompany the anti-glitch, (2) Decrease/variation of braking index after anti-glitch, (3) Anti-glitch is just a period of enhanced spin-down. If there are enough timing observations, a period of enhanced spin-down is expected instead of an anti-glitch. Applications to previous timing events of SGR 1900+14, and PSR J1846-0258 are also included. It is shown that current timing events of 1E 2259+586, SGR 1900+14, and PSR J1846-0258 can be understood safely in the wind braking model. The enhanced spin-down and absence of an anti-glitch before the giant flare of SGR 1806-20 is consistent with the wind braking scenario.

Magnetars are proposed to be peculiar neutron stars powered by their super strong magnetic field. Observationally, anomalous X-ray pulsars and soft gamma-ray repeaters are believed to be magnetar candidates. While more and more multiwave observations of magnetars are available, unfortunately, we see accumulating failed predictions of the traditional magnetar model. These challenges urge rethinking of magnetar. Wind braking of magnetars is one of the alternative modelings. The release of magnetic energy may generate a particle outflow (i.e., particle wind), that results in both an anomalous X-ray luminosity and significantly high spindown rate. In this wind braking scenario, only strong multipole field is necessary for a magnetar (a strong dipole field is no longer needed). Wind braking of magnetars may help us to understand their multiwave radiation properties, including (1) Non-detection of magnetars in Fermi-LAT observations, (2) The timing behaviors of low magnetic field magnetars, (3) The nature of anti-glitches, (4) The criterion for magnetar's radio emission, etc. In the wind braking model of magentars, timing events of magnetars should always be accompanied by radiative events. It is worth noting that the wind engine should be the central point in the research since other efforts with any reasonable energy mechanism may also reproduce the results.

The anomalous X-ray pulsar 4U 0142+61 was studied by the INTEGRAL observations. The hard X-ray spectrum of 18 -- 500 keV for 4U 0142+61 was derived using near 9 years of INTEGRAL/IBIS data. We obtained the average hard X-ray spectrum of 4U 0142+61 with all available data. The spectrum of 4U 0142+61 can be fitted with a power-law with an exponential high energy cutoff. This average spectrum is well fitted with a power-law of $\Gamma\sim 0.51\pm 0.11$ plus a cutoff energy at $128.6\pm 17.2$ keV. The hard X-ray flux of the source from 20 -- 150 keV showed no significant variations (within 20$\%$) from 2003 -- 2011. The spectral profiles have some variability in nine years: photon index varied from 0.3 -- 1.5, and cutoff energies of 110 -- 250 keV. The detection of the high energy cutoff around 130 keV shows some constraints on the radiation mechanisms of magnetars and possibly probes the differences between magnetar and accretion models for these special class of neutron stars. Future HXMT observations could provide stronger constraints on the hard X-ray spectral properties of this source and other magnetar candidates.

The different timing results of the magnetar Swift J1822.3-1606 is analyzed and understood theoretically. It is pointed that different timing solutions are caused not only by timing noise, but also that the period derivative is decreasing after outburst. Both the decreasing period derivative and the large timing noise may be originated from wind braking of the magnetar. Future timing of Swift J1822.3-1606 will help us make clear whether its period derivative is decreasing with time or not.

(adapted)Considering recent observations challenging the traditional magnetar model, we explore the wind braking of magnetars. There is evidence for strong multipole magnetic fields in active magnetars, but the dipole field inferred from spin down measurements may be strongly biased by a particle wind. Recent challenging observations of magnetars may be explained naturally in the wind braking scenario: (1) The supernova energies of magnetars are of normal value; (2) The non-detection in Fermi observations of magnetars; (3) The problem posed by the low-magnetic field soft gamma-ray repeaters; (4) The relation between magnetars and high magnetic field pulsars; (5) A decreasing period derivative during magnetar outbursts. Transient magnetars may still be magnetic dipole braking. This may explain why low luminosity magnetars are more likely to have radio emissions. In the wind braking scenario, magnetars are neutron stars with strong multipole field. For some sources, a strong dipole field may be no longer needed. A magnetism-powered pulsar wind nebula and a braking index smaller than three are the two predictions of the wind braking model.

The magnetar Swift J1834.9-0846 is observed using Nanshan 25 meter radio telescope. No pulsed radio emission is detected. The upper limit on pulsed radio emission from this source is 0.5 mJy. According to the fundamental plane for radio magnetars, this source should have radio emission. Therefore, our results put constraints on the existence of a fundamental plane of magnetar radio emission. We argue that a magnetar's ability to emit radio emission may have little to do with the spin down luminosity and it is related to the magnetar X-ray luminosity. The only necessary condition is a relatively low X-ray luminosity.

The concept of a "magnetar" was proposed mainly because of two factors. First, the X-ray luminosity of Anomalous X-ray Pulsars (AXPs) and Soft Gamma-Ray Repeaters (SGRs) is larger than the rotational energy loss rate, and second, the magnetic field strength calculated from "normal method" is super strong. It is proposed that the radiation energy of magnetar comes from its magnetic fields. Here it is argued that the magnetic field strength calculated through the normal method is incorrect when X-ray luminosity is larger than rotational energy loss rate, because the wind braking is not taken into account. Besides, the "anti-magnetar" and some other X-ray and radio observations are difficult to understand with a magnetar model. Instead of the magnetar, we propose a "quarctar", which is a crusted quark star in an accretion disk, to explain the observations. In this model, the persistent X-ray emission, burst luminosity, spectrum of AXPs and SGRs can be understood naturally. The radio-emitting AXPs, which are challenging the magnetar, can also be explained by the quarctar model.

The key point of studying AXPs/SGRs (anomalous X-ray pulsars/soft gamma-ray repeaters) is relevant to the energy budget. Historically, rotation was thought to be the only free energy of pulsar until the discovery of accretion power in X-ray binaries. AXPs/SGRs could be magnetars if they are magnetism-powered, but would alternatively be quark-star/fallback-disk systems if more and more observations would hardly be understood in the magnetar scenario.

We have analyzed the physical implications of Fermi observations of magnetars. Observationally, no significant detection is reported in Fermi observations of all magnetars. Then there are conflicts between outer gap model in the case of magnetars and Fermi observations. One possible explanation is that magnetars are wind braking instead of magnetic dipole braking. In the wind braking scenario, magnetars are neutron stars with strong multipole field. A strong dipole field is no longer required. A magnetism-powered pulsar wind nebula and a braking index smaller than three are the two predictions of wind braking of magnetars. Future deeper Fermi observations will help us make clear whether they are wind braking or magnetic dipole braking. It will also help us to distinguish between the magnetar model and the accretion model for AXPs and SGRs.

The spin down behaviors of SGR 0418+5729 are investigated. The pulsar spin down model of Contopoulos & Spitkovsky (2006) is applied to SGR 0418+5729. It is shown that SGR 0418+5729 lies below the pulsar death line and its rotation-powered magnetospheric activities may therefore have stopped. The compact star is now spun down by the magnetic dipole moment perpendicular to its rotation axis. Our calculations show that under these assumption there is the possibility of SGR 0418+5729 having a strong dipole magnetic field, if there is a small magnetic inclination angle. Its dipole magnetic field may be much higher than the characteristic magnetic field. Therefore, SGR 0418+5729 may be a normal magnetar instead of a low magnetic field magnetar.

From 2000 to 2010, monitoring of radio emission from the Crab pulsar at Xinjiang Observatory detected a total of nine glitches. The occurrence of glitches appears to be a random process as described by previous researches. A persistent change in pulse frequency and pulse frequency derivative after each glitch was found. There is no obvious correlation between glitch sizes and the time since last glitch. For these glitches $\Delta\nu_{p}$ and $\Delta\dot{\nu}_{p}$ span two orders of magnitude. The pulsar suffered the largest frequency jump ever seen on MJD 53067.1. The size of the glitch is $\sim$ 6.8 $\times 10^{-6}$ Hz, $\sim$ 3.5 times that of the glitch occured in 1989 glitch, with a very large permanent changes in frequency and pulse frequency derivative and followed by a decay with time constant $\sim$ 21 days. The braking index presents significant changes. We attribute this variation to a varying particle wind strength which may be caused by glitch activities. We discuss the properties of detected glitches in Crab pulsar and compare them with glitches in the Vela pulsar.

Anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) are enigmatic pulsar-like objects. The energy budget is the fundamental problem in their studies. In the magnetar model, they are supposed to be powered by the extremely strong magnetic fields (>~ 10^14 G) of neutron stars. Observations for and against the magnetar model are both summarized. Considering the difficulties encountered by the magnetar model to comfortably understand more and more observations, one may doubt that AXPs and SGRs are really magnetars. If they are not magnetar candidates (including magnetar-based models), then they must be "quark star/fallback disk" systems.

The X-ray dim isolated neutron stars (XDINSs) are peculiar pulsar-like objects, characterized by their very well Planck-like spectrum. In studying their spectral energy distributions, the optical/UV excess is a long standing problem. Recently, Kaplan et al. (2011) have measured the optical/UV excess for all seven sources, which is understandable in the resonant cyclotron scattering (RCS) model previously addressed. The RCS model calculations show that the RCS process can account for the observed optical/UV excess for most sources . The flat spectrum of RX J2143.0+0654 may due to contribution from bremsstrahlung emission of the electron system in addition to the RCS process.

Anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) are magnetar candidates, i.e., neutron stars powered by strong magnetic field. If they are indeed magnetars, they will emit high-energy gamma-rays which are detectable by Fermi-LAT according to the outer gap model. However, no significant detection is reported in recent Fermi-LAT observations of all known AXPs and SGRs. Considering the discrepancy between theory and observations, we calculate the theoretical spectra for all AXPs and SGRs with sufficient observational parameters. Our results show that most AXPs and SGRs are high-energy gamma-ray emitters if they are really magnetars. The four AXPs 1E 1547.0-5408, XTE J1810-197, 1E 1048.1-5937, and 4U 0142+61 should have been detected by Fermi-LAT. Then there is conflict between out gap model in the case of magnetars and Fermi observations. Possible explanations in the magnetar model are discussed. On the other hand, if AXPs and SGRs are fallback disk systems, i.e., accretion-powered for the persistent emissions, most of them are not high-energy gamma-ray emitters. Future deep Fermi-LAT observations of AXPs and SGRs will help us make clear whether they are magnetars or fallback disk systems.

Significant research in compact stars is currently focused on two kinds of enigmatic sources: anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs). Although AXPs and SGRs are popularly thought to be magnetars, other models (e.g. the accretion model) to understand the observations can still not be ruled out. It is worth noting that a non-detection in a Fermi/LAT observation of AXP 4U 0142+61 has been reported recently by Sasmaz Mus & Gogus. We propose here that Fermi/LAT observations may distinguish between the magnetar model and the accretion model for AXPs and SGRs. We explain how this null observation of AXP 4U 0142+61 favors the accretion model. Future Fermi/LAT observations of AXP 1E 1547.0-5408 and AXP 1E 1048.1-5937 are highly recommended.

Resonant cyclotron scattering (RCS) in pulsar magnetospheres is considered. The photon diffusion equation (Kompaneets equation) for RCS is derived. The photon system is modeled three dimensionally. Numerical calculations show that there exist not only up scattering but also down scattering of RCS, depending on the parameter space. RCS's possible applications to the spectra energy distributions of magnetar candidates and radio quiet isolated neutron stars (INSs) are point out. The optical/UV excess of INSs may caused by the down scattering of RCS. The calculations for RX J1856.5-3754 and RX J0720.4-3125 are presented and compared with their observational data. In our model, the INSs are proposed to be normal neutron stars, although the quark star hypothesis is still possible. The low pulsation amplitude of INSs is a natural consequence in the RCS model.

A phase transition from paramagnetism to ferromagnetism in neutron star interior is explored. Since there is $^3$P$_2$ neutron superfluid in neutron star interior, it can be treated as a system of magnetic dipoles. Under the presence of background magnetic field, the magnetic dipoles tend to align in the same direction. Below a critical temperature, there is a phase transition from paramagnetism to ferromagnetism. And this gives a convenient explanation of the strong magnetic field of magnetars. In our point of view, there is an upper limit for the magnetic field strength of magnetars. The maximum field strength of magnetars is about $(3.0-4.0)\times 10^{15}$ G. This can be tested directly by further investigations. Magnetars are instable due to the ultra high Fermi energy of electrons. The Landau column becomes a very long cylinder along the magnetic field, but it is very narrow and the Fermi energy of electron gas is given as $E_F(e)\approx40(B/B_{cr})^{1/4}$ when $B\gg B_{cr}$. $E_F(e)\approx90MeV$ When $B\sim10^{15}$ G. Hence, the electron capture process $e^-+p\to n+\nu_e$ will be happen rapidly. Thus the $^3$P$_2$ Cooper pairs will be destroyed quickly by the outgoing neutrons with high energy. It will cause the isotropic superfluid disappear and then the magnetic field induced by the $^3$P$_2$ Cooper pairs will be also disappear. These energy will immediately be transmitted into thermal energy and then transformed into the radiation energy with X-ray - soft $\gamma$-ray. We may get a conclusion that the activity of magnetars originates from instability caused by the high Fermi energy of electrons in extra strong magnetic field.

There is ${}^3P_2$ neutron superfluid region in NS (neutron star) interior. For a rotating NS, the ${}^3P_2$ superfluid region is like a system of rotating magnetic dipoles. It will give out electromagnetic radiation, which may provides a new heating mechanism of NSs. This heating mechanism plus some cooling agent may give sound explanation to NS glitches.

According to recent observational and theoretical progresses, the DUrca process (direct Urca process) may be excluded from the category of neutron star cooling mechanisms. This result combined with the latest nuclear symmetry energy experiments, will provide us an independent way of testing the EOS (equation of state) for supranormal density. For example, soft EOSs such as FPS will probably be excluded.

In this paper we present a new result, namely that the primal magnetic field of the collapsed core during a supernova explosion will, as a result of the conservation of magnetic flux, receive a massive boost to more than 90 times its original value by the Pauli paramagnetization of the highly degenerate relativistic electron gas just after the formation of the neutron star. Thus, the observed super-strong magnetic field of neutron stars may originate from the induced Pauli paramagnetization of the highly degenerate relativistic electron gas in the interior of the neutron star. We therefore have an apparently natural explanation for the surface magnetic field of a neutron star.