We have used the Nanshan 25-m Radio Telescope at Xinjiang Astronomical Observatory to obtain timing observations of 87 pulsars from 2002 July to 2014 March. Using the "Cholesky" timing analysis method we have determined positions and proper motions for 48 pulsars, 24 of which are improved positions compared to previously published values. We also present the first published proper motions for nine pulsars and improved proper motions for 21 pulsars using pulsar timing and position comparison method. The pulsar rotation parameters are derived and are more accurate than previously published values for 36 pulsars. Glitches are detected in three pulsars: PSRs J1722$-$3632, J1852$-$0635 and J1957+2831. For the first two, the glitches are large, with $\Delta\nu_g/\nu > 10^{-6}$, and they are the first detected glitches in these pulsars. PSR J1722$-$3632 is the second oldest pulsar with large glitch. For the middle-age pulsars ($\tau_c > 10^5$~yr), the calculated braking indices, $|n|$, are strongly correlated with $\tau_c$ and the numbers of positive and negative values of $n$ are almost equal. For young pulsars ($\tau_c < 10^5$~yr), there is no correlation between $|n|$ and $\tau_c$ and most have $n>0$.

We report the discovery of a binary millisecond pulsar (namely PSR J1641+3627F or M13F) in the globular cluster M13 (NGC 6205) and timing solutions of M13A to F using observations made with the Five-hundred-metre Aperture Spherical radio Telescope (FAST). PSR J1641+3627F has a spin period of 3.00 ms and an orbital period of 1.4 days. The most likely companion mass is 0.16 M$_{\odot}$. M13A to E all have short spin periods and small period derivatives. We also confirm that the binary millisecond pulsar PSR J1641$+$3627E (also M13E) is a black widow with a companion mass around 0.02 M$_{\odot}$. We find that all the binary systems have low eccentricities compared to those typical for globular cluster pulsars and that they decrease with distance from the cluster core. This is consistent with what is expected as this cluster has a very low encounter rate per binary.

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.

Rotating radio transients (RRATs) are peculiar astronomical objects whose emission mechanism remains under investigation. In this paper, we present observations of three RRATs, J1538+2345, J1854+0306 and J1913+1330, observed with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Specifically, we analyze the mean pulse profiles and temporal flux density evolutions of the RRATs. Owing to the high sensitivity of FAST, the derived burst rates of the three RRATs are higher than those in previous reports. RRAT J1854+0306 exhibited a time-dynamic mean pulse profile, whereas RRAT J1913+1330 showed distinct radiation and nulling segments on its pulse intensity trains. The mean pulse profile variation with frequency is also studied for RRAT J1538+2345 and RRAT J1913+1330, and the profiles at different frequencies could be well fitted with a cone-core model and a conal-beam model, respectively.

PSR B0919+06 is known for its abnormal emission phenomenon, where the pulse emission window occasionally shifts progressively in longitude and returns afterwards. The physical mechanism behind this phenomenon is still under investigation. In this paper, we present our ultra-wideband observation of this pulsar using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), with simultaneous measurements in the frequency ranges 280-780 MHz and 1250-1550 MHz. We have identified three abnormal events, each of which becomes less apparent as the frequency decreases. At 1400 MHz, the averaged profile slightly shifted after the first and third abnormal events, implying a relationship between abnormal event and profile variation. We also found a linear trend in the left-edge position of the averaged profiles between the first and third event as well as after the third event, suggesting the existence of a slow-drifting mode between the two major events. The second event has a comparatively small shift in phase and is thus categorized as a `small flare state'. During the third event, a sequence of approximately nine pulses was seen to significantly weaken in all frequency bands, likely associated with the pseudo-nulling observed at 150 MHz. A three-component de-composition analysis of the normal averaged profiles shows that the trailing component is dominant at our observing frequencies, while the centre component has a comparatively steeper spectrum. We found the overall flux density in an abnormal event to slightly differ from that in an ordinary state, and the difference shows a frequency dependence. A comparison of the normal, abnormal and dimmed averaged profile indicates that the leading component is likely to be stable in all states.

We presented timing results of PSR B2035+36 using ~9-yr observations with the Nanshan 25-m radio telescope. PSR B2035+36 was reported to exhibit significant changes in pulse profile correlated with spin-down state variations. We found that the pulsar underwent a glitch with a jump in the frequency of $\Delta{\nu}\sim12.4(5)\, \rm nHz$ around MJD 52950. Unusually, the spin-down rate increased persistently over $800$ days after the glitch, and the average spin-down rate of post-glitch was about $9.6\%$ larger than pre-glitch. Accompanied with the glitch activity, pulse profile became narrower. The pulsar began to switch between two emission modes after the glitch, with pulse width ($W_{\ 50 \rm mean}$) of $8.5(7)^{\circ}$ and $3.7(3)^{\circ}$, respectively. Besides that, the relatively narrow pulse profile gradually became dominant. All of the observations indicate that there should be connection between magnetospheric behavior and glitch activity. We discuss one possibility of magnetosphere fluctuation triggered by glitch event.

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.

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.

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.

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.