Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

We present here the analysis of giant micropulses from the Vela pulsar. A total of 4187 giant micropulses with peak flux density $>$2.5 Jy were detected during almost 4 hours of observations carried out with the Yunnan 40-m radio telescope at 6800 MHz. Nine of the giant micropulses arrived approximately 3 to 4 ms earlier than the peak of average pulse profile, longer than that at lower frequencies. The remaining giant micropulses were clustered into three distributions which correspond to three main emission regions, including four occurring on the trailing edge of averaged profile.

There are many stars that are rotating spheroids in the Universe, and studying them is of very important significance. Since the times of Newton, many astronomers and physicists have researched gravitational properties of stars by considering the moment equations derived from Eulerian hydrodynamic equations. In this paper we study the scattering of spinors of the Dirac equation, and in particular investigate the scattering issue in the limit case of rotating Maclaurin spheroids. Firstly we give the metric of a rotating ellipsoid star, then write the Dirac equation under this metric, and finally derive the scattering solution to the Dirac equation and establish a relation between differential scattering cross-section, $\sigma$, and stellar matter density, $\mu$. It is found that the sensitivity of $\sigma$ to the change in $\mu$ is proportional to the density $\mu$. Because of weak gravitational field and constant mass density, our results are reasonable. The results can be applied to white dwarfs, main sequence stars, red giants, supergiant stars and so on, as long as their gravitational fields are so weak that they can be treated in the Newtonan approximations, and the fluid is assumed to be incompressible. Notice that we take the star's matter density to be its average density and the star is not taken to be compact. Obviously our results cannot be used to study neutron stars and black holes. In particular, our results are suitable for white dwarfs, which have average densities of about $10^{5}-10^{6}$\,g~cm$^{-3}$, corresponding to a range of mass of about $0.21-0.61 M_{\bigodot}$ and a range of radius of about $6000-10000$\,km.

We present long-term observations of two intermittent pulsars, PSRs~J1832+0029 and J1841$-$0500 using the Parkes 64\,m radio telescope. The radio emission for these pulsars switches "off" for year-long durations. Our new observations have enabled us to improve the determination of the on-off timescales and the spin down rates during those emission states. In general our results agree with previous studies of these pulsars, but we now have significantly longer data spans. We have identified two unexpected signatures in the data. Weak emission was detected in a single observation of PSR~J1832$+$0029 during an "off" emission state. For PSR~J1841$-$0500, we identified a quasi-periodic fluctuation in the intensities of the detectable single pulses, with a modulation period between 21 and 36 pulse periods.

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

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.

By analysing the data acquired from the Parkes 64-m radio telescope at 1369 MHz, we report on the phase-stationary non-drift amplitude modulation observed in PSR J1048-5832. The high-sensitivity observations revealed that the central and trailing components of the pulse profile of this pulsar switch between a strong mode and a weak mode periodically. However, the leading component remains unchanged. Polarization properties of the strong and weak modes are investigated. Considering the similarity to mode changing, we argue that the periodic amplitude modulation in PSR J1048$-$5832 is periodic mode changing. The fluctuation spectral analysis showed that the modulation period is very short (~2.1 s or 17 P1), where P1 is the rotation period of the pulsar. We find that this periodic amplitude modulation is hard to explain by existing models that account for the periodic phenomena in pulsars like subpulse drifting.

According to the recycling model, neutron stars in low-mass X-ray binaries were spun up to millisecond pulsars (MSPs), which indicates that all MSPs in the Galactic plane ought to be harbored in binaries. However, about $20\%$ Galactic field MSPs are found to be solitary. To interpret this problem, we assume that the accreting neutron star in binaries may collapse and become a strange star when it reaches some critical mass limit. Mass loss and a weak kick induced by asymmetric collapse during the phase transition (PT) from neutron star to strange star can result in isolated MSPs. In this work, we use a population-synthesis code to examine the PT model. The simulated results show that a kick velocity of $\sim60~{\rm km~s}^{-1}$ can produce $\sim6\times10^3$ isolated MSPs and birth rate of $\sim6.6\times10^{-7} {\rm ~yr}^{-1}$ in the Galaxy, which is approximately in agreement with predictions from observations. For the purpose of comparisons with future observation, we also give the mass distributions of radio and X-ray binary MSPs, along with the delay time distribution.

We describe an ultra-wide-bandwidth, low-frequency receiver ("UWL") recently installed on the Parkes radio telescope. The receiver system provides continuous frequency coverage from 704 to 4032 MHz. For much of the band (~60%) the system temperature is approximately 22K and the receiver system remains in a linear regime even in the presence of strong mobile phone transmissions. We discuss the scientific and technical aspects of the new receiver including its astronomical objectives, as well as the feed, receiver, digitiser and signal-processor design. We describe the pipeline routines that form the archive-ready data products and how those data files can be accessed from the archives. The system performance is quantified including the system noise and linearity, beam shape, antenna efficiency, polarisation calibration and timing stability.

We build a conceptual coupled model of the climate and tidal evolution of the Earth-Moon system to find the influence of the former on the latter. An energy balance model is applied to calculate steady-state temperature field from the mean annual insolation as a function of varying astronomical parameters. A harmonic oscillator model is applied to integrate the lunar orbit and Earth's rotation with the tidal torque dependent on the dominant natural frequency of ocean. An ocean geometry acts as a bridge between temperature and oceanic frequency. On assumptions of a fixed hemispherical continent and an equatorial circular lunar orbit, considering only the 41 kyr periodicity of Earth's obliquity $\varepsilon$ and the $M_2$ tide, simulations are performed near tidal resonance for $10^6$ yr. It is verified that the climate can influence the tidal evolution via ocean. Compared with the tidal evolution with constant $\varepsilon$, that with varying $\varepsilon$ is slowed down; the Earth-Moon distance oscillates in phase with $\varepsilon$ before the resonance maximum but exactly out of phase after that; the displacement of the oscillation is in positive correlation with the difference between oceanic frequency and tidal frequency.

PSR J1825$-$0935 (PSR B1822$-$09) switches between radio-quiet (Q-mode) and radio-bright (B-mode) modes. The Q-mode is known to have a periodic fluctuation that modulates both the interpulse and the main pulse with the same period. Earlier investigators argued that the periodic Q-mode modulation is associated with drifting subpulses. We report on single-pulse observations of PSR J1825$-$0935 that were made using the Parkes 64-m radio telescope with a central frequency of 1369 MHz. The high-sensitivity observations revealed that the periodic Q-mode modulation is in fact a periodic longitude-stationary intensity modulation occurring in the interpulse and the main pulse. The fluctuation spectral analysis showed that the modulation period is about $43 P_1$, where $P_1$ is the rotation period of the pulsar. Furthermore, we confirm that the modulation patterns in the interpulse and the main pulse are phase-locked. Specifically, the intensities of the interpulse and the immediately following main pulse are more highly correlated than for the main pulse and interpulse at any other lag. Polarization properties of the strong and weak Q-mode states are different, even for the trailing part of the main pulse which does not show the periodic intensity modulation.

The Beijing-Arizona Sky Survey (BASS) is an imaging survey and uses the 2.3 m Bok telescope at Kitt Peak. In order to tap the astrometric potential of the Bok telescope and improve the astrometry of BASS, a distortion solution for the Bok telescope is made. In the past, we used a single lookup table to correct all the positional errors. However, this method can not be applied to the reduction of the observation for the Bok telescope where four CCD chips are equipped with. Then quite different from our previous method, two third-order polynomials were used to fit the lookup table. By using the polynomial Geometric Distortion (called GD hereafter) correction the astrometry of BASS is improved greatly. Moreover, an additional lookup table correction is found to be more effective to obtain a final GD. The results show that the positional measurement precision of the appropriate bright stars is estimated at about 20 mas and even better in each direction. Besides, the relative positions of the chips are measured. The change of the inter-chip gaps in horizontal or vertical is no more than 5 pixel between 2016 and 2017 and the change of the roll angle is no more than 0.1 degree.

With the constraint from gravitational wave emission of a binary merger system (GW170817) and two-solar-mass pulsar observations, we investigate the r-mode instability windows of strange stars with unpaired and color-flavor-locked phase strange quark matter. Shear viscosities due to surface rubbing and electron-electron scattering are taken into account in this work. The results show that the effects of the equation of state of unpaired strange quark matter are only dominant at low temperature, but do not have significant effects on strange stars in the color-flavor-locked phase. A color-flavor-locked phase strange star, which is surrounded by an insulating nuclear crust, seems to be consistent with observational data of young pulsars. We find that an additional enhanced dissipation mechanisms might exist in SAX J1808.4-3658. Fast spinning young pulsar PSR J0537-6910 is a primary source for detecting gravitational waves from a rotating strange star, and young pulsars might be strange stars with color-flavor-locked phase strange quark matter.

The ANtarctic Impulsive Transient Antenna (ANITA) NASA long-duration balloon payload completed its fourth flight in December 2016, after 28 days of flight time. ANITA is sensitive to impulsive broadband radio emission from interactions of ultra-high-energy neutrinos in polar ice (Askaryan emission). We present the results of two separate blind analyses searching for signals from Askaryan emission in the data from the fourth flight of ANITA. The more sensitive analysis, with a better expected limit, has a background estimate of $0.64^{+0.69}_{-0.45}$ and an analysis efficiency of $82\pm2\%$. The second analysis has a background estimate of $0.34^{+0.66}_{-0.16}$ and an analysis efficiency of $71\pm6\%$. Each analysis found one event in the signal region, consistent with the background estimate for each analysis. The resulting limit further tightens the constraints on the diffuse flux of ultra-high-energy neutrinos at energies above $10^{19.5}$ eV.

A new method based on parameters iteration technique has been developed to determine the optimal subreflector position for shaped Cassegrain antennas to improve the electromagnetic (EM) performance distorted by gravity. Both the features of shaped surface and the relationship between optical path difference (OPD) and far field beam pattern are employed. By describing the shaped dual-reflector surface as a standard discrete parabola set, we can utilize the optical features of the standard Cassegrain system in the classical OPD relationship. Then, the actual far field beam pattern is expressed as the synthesis of ideal beam and error beam by decomposing subreflector adjustment parameters using mechanical-electromagnetic-field-coupling-model (MEFCM). Furthermore, a numerical method for determining optimal subreflector position is presented. The proposed method is based on the iteration technique of subreflector adjustment parameters, and the optimal far field pattern is used as the iteration goal. The numerical solution of optimal adjustment parameters can be obtained rapidly. Results of a 25 m Shaped Cassegrain antenna demonstrate that the adjustment of the subreflector to the optimal position determined by the proposed method can improve the EM performance effectively.

One large glitch was detected in PSR B1737$-$30 using data spanning from MJD 57999 to 58406 obtained with the newly built Shanghai Tian Ma Radio Telescope (TMRT). The glitch took place at the time around MJD 58232.4 when the pulsar underwent an increase in the rotation frequency of $\Delta \nu$ about 1.38$\times 10^{-6}$ Hz, corresponding to a fractional step change of $\Delta \nu / \nu$ $\thicksim$ 8.39$\times 10^{-7}$. Post$\textrm{-}$glitch $\nu$ gradually decreased to the pre$\textrm{-}$glitch value. The frequency derivative was observed to undergo a step change of about $-$9$\times 10^{-16}$ s$^{-2}$. Since July 1987, there are 36 glitches already reported in PSR B1737$-$30 including this one. According to our analysis, the glitch size distribution is well described by the power law with index of 1.13. The distribution of the interval between two adjacent glitches (waiting time $\Delta T$) follows a Poissonian probability density function. For PSR B1737$-$30, the interval is prone to be long after a large glitch. But no correlation is found between glitch size and the interval since previous glitch.

Timing analysis of PSR J1705$-$1906 using data from Nanshan 25-m and Parkes 64-m radio telescopes, which span over fourteen years, shows that the pulsar exhibits significant proper motion, and rotation instability. We updated the astrometry parameters and the spin parameters of the pulsar. In order to minimize the effect of timing irregularities on measuring its position, we employ the Cholesky method to analyse the timing noise. We obtain the proper motion of $-$77(3) \,mas\,yr$^{-1}$ in right ascension and $-$38(29) \,mas\,yr$^{-1}$ in declination. The power spectrum of timing noise is analyzed for the first time, which gives the spectral exponent $\alpha=-5.2$ for the power-law model indicating that the fluctuations in spin frequency and spin-down rate dominate the red noise. We detect two small glitches from this pulsar with fractional jump in spin frequency of $\Delta \nu/\nu\sim2.9\times10^{-10}$ around MJD~55199 and $\Delta \nu/\nu\sim2.7\times10^{-10}$ around MJD~55953. Investigations of pulse profile at different time segments suggest no significant changes in the pulse profiles around the two glitches.

Based on the model of the accretion-induced magnetic field decay of a neutron star (NS), millisecond pulsars (MSPs) will obtain their minimum magnetic field when the NS magnetosphere radius shrinks to the stellar surface during the binary accretion phase. We find that this minimummagnetic field is related to the accretion rate Mdot as Bmin ~2.0*10^7 G( Mdot/Mdot_min)^1/2, where Mdot_min = 4.6*10^15 g/s is the average minimum accretion rate required for MSP formation and is constrained by the long-term accretion time, which corresponds to the companion lifetime, being less than the Hubble time. The value of Bmin is consistent with that of observed radio MSPs and accreting MSPs in low-mass X-ray binaries, which can be found the illustrated case of the minimum and present field strength of SAX J1808.4-3658. The prediction of the minimum magnetic field of MSPs would be the lowest field strength of NSs in the Universe, which could constrain the evolution mechanism of the magnetic field of accreting NSs.

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.

We report on single-pulse observations of the Galactic Center magnetar PSR J1745$-$2900 that were made using the Parkes 64-m radio telescope with a central frequency of 3.1 GHz at five observing epochs between 2013 July and August. The shape of the integrated pulse profiles was relatively stable across the five observations, indicating that the pulsar was in a stable state between MJDs 56475 and 56514. This extends the known stable state of this pulsar to 6.8 months. Short term pulse shape variations were also detected. It is shown that this pulsar switches between two emission modes frequently and that the typical duration of each mode is about ten minutes. No giant pulses or subpulse drifting were observed. Apparent nulls in the pulse emission were detected on MJD 56500. Although there are many differences between the radio emission of magnetars and normal radio pulsars, they also share some properties. The detection of mode changing and pulse nulling in PSR J1745$-$2900 suggests that the basic radio emission process for magnetars and normal pulsars is the same.

Timing analysis provides information about the dynamics of matter accreting on to neutron stars and black holes, and hence is crucial for studying the physics of the accretion flow around these objects. It is difficult, however, to associate the different variability components with each of the spectral components of the accretion flow. We apply several new methods to two Rossi X-ray Timing Explorer observations of the black hole binary GRS 1915+105 during its heartbeat state to explore the origin of the X-ray variability and the interactions of the accretion-flow components. We offer a promising window into the disc--corona interaction through analysing the formation regions of the disc aperiodic variabilities with different time-scales via comparing the corresponding transition energies of the amplitude-ratio spectra. In a previous paper, we analysed the Fourier power density as a function of energy and frequency to study the origin of the aperiodic variability, and combined that analysis with the phase lag as a function of frequency to derive a picture of the disc--corona interaction in this source. We here, for the first time, investigate the phase lag as a function of energy and frequency, and display some interesting details of the disc--corona interaction. Besides, the results from the shape of amplitude-ratio spectrum and from several other aspects suggest that the quasi-periodic oscillation originates from the corona.

The low braking-index pulsar PSR J1734$-$3333 could be born with superhigh internal magnetic fields $B_{\rm in}\sim10^{15}-10^{16}$ G, and undergo a supercritical accretion soon after its formation in a supernova explosion. The buried multipole magnetic fields will merger into a dipole magnetic field. Since the magnetic flow transfers from the core to the crust of the pulsar, its surface dipole field grows quickly at a power-law form assumed until it saturates at the level of internal dipole field. The increase in surface dipole magnetic field results in the observed low braking index of $n=0.9(2)$. Keeping an average field growth index $\varepsilon=1.34(6)$, this pulsar will become a magnetar with surface dipole magnetic field at the equator $B_{\rm d}\sim 2.6(1)\times 10^{14}$\,G and $\sim 5.3(2)\times 10^{14}$\,G after next 50\,kyrs and 100\,kys, respectively.

In this work, we reinvestigate the electron fraction $Y_{e}$ and electron Fermi energy $E_{F}(e)$ of neutron stars, based on our previous work of Li et al.(2016), in which we firstly deduced a special solution to $E_{F}(e)$, and then obtained several useful analytical formulae for $Y_{\rm e}$ and matter density $\rho$ within classical models and the relativistic mean field(RMF) theory using numerically fitting. The advantages of this work include the following aspects:(1) The linear functions are substituted for the nonlinear exponential functions used in the previous work. This method may be more simple, and closer to realistic equation of state\,(EoS) of a neutron star(NS), because there are linear or quasi-linear relationships between number fractions of leptons and matter density, which can be seen by solving NS EoS; (2)we introduce a dimensionless variable $\varrho$\,($\varrho=\rho/\rho_0$, $\rho_{0}$ is the standard saturated nuclear density), which greatly reduces the scope of the fitting coefficients;(3)we present numerical errors including absolute and relative deviations between the data and fit. By numerically simulating, we have obtained several analytical formulae for $Y_{e}$ and $\rho$ for both APR98 and RMF models. Combining these analytical formulae with the special solution, we can calculate the value of $E_{\rm F}(e)$ for any given matter density. Since $Y_e$ and $E_{ F}(e)$ are important in assessing cooling rate of a NS and the possibility of kaon/pion condensation in the NS interior, this study could be useful in the future study on the thermal evolution of a NS.

In this work, we interpreted the high braking index of PSR J1640$-$4631 with a combination of the magneto-dipole radiation and dipole magnetic field decay models. By introducing a mean rotation energy conversion coefficient $\overline{\zeta}$, the ratio of the total high-energy photon energy to the total rotation energy loss in the whole life of the pulsar, and combining the pulsar's high-energy and timing observations with reliable nuclear equation of state, we estimate the pulsar's initial spin period, $P_{0}\sim (17-44)$ ms, corresponding to the moment of inertia $I\sim (0.8-2.1)\times 10^{45}$ g cm$^{2}$. Assuming that PSR J1640$-$4631 has experienced a long-term exponential decay of the dipole magnetic field, we calculate the true age $t_{\rm age}$, the effective magnetic field decay timescale $\tau_{D}$, and the initial surface dipole magnetic field at the pole $B_{p}(0)$ of the pulsar to be $(2900-3100)$ yrs, $1.07(2)\times10^{5}$ yrs, and $(1.84-4.20)\times10^{13}$ G, respectively. The measured braking index of $n=3.15(3)$ for PSR J1640$-$4631 is attributed to its long-term dipole magnetic field decay and a low magnetic field decay rate, $dB_{\rm p}/dt\sim -(1.66-3.85)\times10^{8}$ G yr$^{-1}$. Our model can be applied to both the high braking index ($n>3$) and low braking index ($n<3$) pulsars, tested by the future polarization, timing, and high-energy observations of PSR J1640$-$4631.

We report detailed investigation of nulling and drifting behavior of two pulsars PSRs J1741-0840 and J1840-0840 observed from the Giant Meterwave Radio Telescope at 625 MHz. PSR J1741-840 was found to show nulling fraction (NF) of around 30(5)% while PSR J1840-0840 was shown to have NF of around 50(6)%. We measured drifting behavior from different profile components in PSR J1840-0840 for the first time with the leading component showing drifting with 13.5(0.7) periods while the weak trailing component showed drifting of around 18(1) periods. Large nulling do hamper accuracy of these quantities derived using standard Fourier techniques. A more accurate comparison was drawn from driftband slopes, measured after sub-pulse modeling. These measurements revealed interesting sporadic and irregular drifting behavior in both pulsars. We conclude that the previously reported different drifting periodicities in the trailing component of PSR J1741-0840 is likely due to the spread in these driftband slopes. We also find that both components of PSR J1840-0840 show similar driftband slopes within the uncertainties. Unique nulling-drifting interaction is identified in PSR J1840-0840 where, in most occasions, the pulsar tends to start nulling after what appears to be an end of a driftband. Similarly, when the pulsar switches back to an emission phase, in most occasions it starts at the beginning of a new driftband in both components. Such behaviors have not been detected in any other pulsars to our knowledge. We also found that PSR J1741-0840 seems to have no memory of its previous burst phase while PSR J1840-0840 clearly exhibits memory of its previous state even after longer nulls for both components. We discuss possible explanations for these intriguing nulling-drifting interactions seen in both pulsars based on various pulsar nulling models.

The statistical uncertainties of 13 model parameters in the Weizs\"acker-Skyrme (WS*) mass model are investigated for the first time with an efficient approach, and the propagated errors in the predicted masses are estimated. The discrepancies between the predicted masses and the experimental data, including the new data in AME2016, are almost all smaller than the model errors. For neutron-rich heavy nuclei, the model errors increase considerably, and go up to a few MeV when the nucleus approaches the neutron drip line. The most sensitive model parameter which causes the largest statistical error is analyzed for all bound nuclei. We find that the two coefficients of symmetry energy term significantly influence the mass predictions of extremely neutron-rich nuclei, and the deformation energy coefficients play a key role for well-deformed nuclei around the $\beta$-stability line.

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.

We derive the Sun's offset from the local mean Galactic plane($z_\odot$) using the observed $z$ distribution of young pulsars. Pulsar distances are obtained from measurements of annual parallax, HI absorption spectra or associations where available and otherwise from the observed pulsar dispersion and a model for the distribution of free electrons in the Galaxy. We fit the cumulative distribution function for a ${\rm sech}^2(z)$ distribution function, representing an isothermal self-gravitating disk, with uncertainties being estimated using the bootstrap method. We take pulsars having characteristic age $\tau_c<10^{6.5}$~yr and located within 4.5~kpc of the Sun, omitting those within the local spiral arm and those significantly affected by the Galactic warp, and solve for $z_\odot$ and the scale height, $H$, for different cutoffs in $\tau_c$. We compute these quantities using just the independently determined distances, and these together with DM-based distances separately using the YMW16 and NE2001 Galactic electron density models. We find that an age cutoff at $10^{5.75}$~yr with YMW16 DM-distances gives the best results with a minimum uncertainty in $z_\odot$ and an asymptotically stable value for $H$ showing that, at this age and below, the observed pulsar $z$-distribution is dominated by the dispersion in their birth locations. From this sample of 115 pulsars, we obtain $z_\odot=13.4\pm$4.4~pc and $H=56.9\pm$6.5~pc, similar to estimated scale heights for OB stars and open clusters. Consistent results are obtained using the independent-only distances and using the NE2001 model for the DM-based distances.

Pulsar positions can be measured with high precision using both pulsar timing methods and very-long-baseline interferometry (VLBI). Pulsar timing positions are referenced to a solar-system ephemeris, whereas VLBI positions are referenced to distant quasars. Here we compare pulsar positions from published VLBI measurements with those obtained from pulsar timing data from the Nanshan and Parkes radio telescopes in order to relate the two reference frames. We find that the timing positions differ significantly from the VLBI positions (and also differ between different ephemerides). A statistically significant change in the obliquity of the ecliptic of $2.16\pm0.33$\,mas is found for the JPL ephemeris DE405, but no significant rotation is found in subsequent JPL ephemerides. The accuracy with which we can relate the two frames is limited by the current uncertainties in the VLBI reference source positions and in matching the pulsars to their reference source. Not only do the timing positions depend on the ephemeris used in computing them, but also different segments of the timing data lead to varying position estimates. These variations are mostly common to all ephemerides, but slight changes are seen at the 10$\mu$as level between ephemerides.

755 new CCD observations during the years 2014-2016 have been reduced to derive the precise positions of Triton, the first satellite of Neptune. The observations were made by the 1 m telescope at Yunnan Observatory over fifteen nights. The positions of Triton are measured with respect to the stars in Gaia DR1 star catalogue. The theoretical position of Triton was retrieved from the Jet Propulsion Laboratory (JPL) ephemeris nep081xl, while the position of Neptune was obtained from the JPL ephemeris DE431mx. Our results show that the mean O-Cs (observed minus computed) are 0.042 and -0.006 arcsec in right ascension and declination respectively. The dispersions of our observations are estimated at about 0.012 arcsec in each direction.

In order to obtain high precision CCD positions of Himalia, the sixth Jovian satellite, a total of 598 CCD observations have been obtained during the years 2015-2016. The observations were made by using the 2.4 m and 1 m telescopes administered by Yunnan Observatories over 27 nights. Several factors which would influence the positional precision of Himalia were analyzed, including the reference star catalogue used, the geometric distortion and the phase effect. By taking advantage of its unprecedented positional precision, the recently released catalogue Gaia DR1 was chosen to match reference stars in the CCD frames of both Himalia and open clusters which were observed for deriving the geometric distortion. The latest version of SOFA library was used to calculate the positions of reference stars. The theoretical positions of Himalia were retrieved from the Jet Propulsion Laboratory Horizons System which includes the satellite ephemeris JUP300, while the positions of Jupiter were based on the planetary ephemeris DE431. Our results showed that the means of observed minus computed (O-C) residuals are 0.071 and -0.001 arcsec in right ascension and declination, respectively. Their standard deviations are estimated at about 0.03 arcsec in each direction.

The geometric distortion of CCD field of view has direct influence on the positional measurements of CCD observations. In order to obtain high precision astrometric results, the geometric distortion should be derived and corrected precisely. As presented in our previous work Peng et al. (2012), a convenient solution has been carried out and also been made with successful application to Phoebe's observations. In order to further improve the solution, an orthogonal method based on the Zernike polynomials is used in this work. Four nights of CCD observations including Himalia, the sixth satellite of Jupiter, and open clusters (NGC1664 or NGC2324) on each night have been processed to make an application. The observations were obtained from the 2.4 m telescope administered by Yunnan Observatories. The catalog UCAC4 was used to match reference stars in all of the CCD frames. The ephemeris of Himalia is retrieved from the (IMCCE). Our results show that the means of observed minus computed (O-C) positional residuals are -0.034 and -0.026 arcsec in right ascension and declination, respectively. The corresponding standard deviations are 0.031 and 0.028 arcsec. The measurement dispersion is significantly improved than that by using our previous solution.

Pulse arrival time measurements over fourteen years with the Nanshan 25-m and Parkes 64-m radio telescopes have been used to determine the average profile and timing properties for PSR J1757$-$2421 (B1754$-$24). Analysis of the radio profile data shows a large variation of spectral index across the profile and an unusual increase in component separation with increasing frequency. The timing observations show that PSR B1754$-$24 underwent a large glitch with a fractional increase in spin frequency of $\Delta\nu_{\rm g}/\nu\sim 7.8\times10^{-6}$ in May 2011. At the time of the glitch there was a large permanent jump in the rotational frequency accompanied by two exponential recovery terms with timescales of 15(6) and 97(15) days, respectively.

We present a timing analysis of two Rossi X-ray Timing Explorer observations of the microquasar GRS 1915+105 during the heartbeat state. The phase-frequency-power maps show that the intermediate-frequency aperiodic X-ray variability weakens as the source softens in the slow rise phase, and when the quasi-periodic oscillation disappears in the rise phase of the pulse of the double-peaked class its sub-harmonic is still present with a hard phase lag. In the slow rise phase, the energy-frequency-power maps show that most of the aperiodic variability is produced in the corona, and may also induce the aperiodic variability observed at low energies from an accretion disk, which is further supported by the soft phase lag especially in the intermediate-frequency range (with a time delay up to 20 ms). In the rise phase of the pulse, the low-frequency aperiodic variability is enhanced significantly and there is a prominent hard lag (with a time delay up to 50 ms), indicating that the variability is induced by extension of the disk toward small radii as implied by the increase in flux and propagates into the corona. However, during the hard pulse of the double-peaked class, the variability shows no significant lag, which may be attributed to an optically thick corona. These timing results are generally consistent with the spectral results presented by Neilsen et al. (2011, 2012) which indicated that the slow rise phase corresponds to a local Eddington limit and the rise phase of the pulse corresponds to a radiation pressure instability in the disk.

We present a new model for the distribution of free electrons in the Galaxy, the Magellanic Clouds and the intergalactic medium (IGM) that can be used to estimate distances to real or simulated pulsars and fast radio bursts (FRBs) based on their dispersion measure (DM). The Galactic model has an extended thick disk representing the so-called warm interstellar medium, a thin disk representing the Galactic molecular ring, spiral arms based on a recent fit to Galactic HII regions, a Galactic Center disk and seven local features including the Gum Nebula, Galactic Loop I and the Local Bubble. An offset of the Sun from the Galactic plane and a warp of the outer Galactic disk are included in the model. Parameters of the Galactic model are determined by fitting to 189 pulsars with independently determined distances and DMs. Simple models are used for the Magellanic Clouds and the IGM. Galactic model distances are within the uncertainty range for 86 of the 189 independently determined distances and within 20\% of the nearest limit for a further 38 pulsars. We estimate that 95\% of predicted Galactic pulsar distances will have a relative error of less than a factor of 0.9. The predictions of YMW16 are compared to those of the TC93 and NE2001 models showing that YMW16 performs significantly better on all measures. Timescales for pulse broadening due to interstellar scattering are estimated for (real or simulated) Galactic and Magellanic Cloud pulsars and FRBs.

In order to study the potential in high precision CCD astrometry of irregular satellites, we have made experimental observations for Himalia, the sixth and irregular satellite of Jupiter. A total of 185 CCD observations were obtained by using the 2.4 m telescope and 1 m telescope at Yunnan Observatories over ten nights. Preliminary analysis for the observations were made, including the geometric distortion, the atmospheric refraction, and also the phase effect. All positions of Himalia are measured relative to the reference stars from the catalogue UCAC4 in each CCD field of view. The theoretical positions of Himalia were retrieved from the IMCCE, while the positions of Jupiter were obtained based on the planetary theory INPOP13C. The results show that the mean (O-C) (observed minus computed) residuals are -0.004 and -0.002 arcsec in right ascension and declination, respectively. The standard deviation of (O-C) residuals are estimated at about 0.04 arcsec in each direction.

We formulate the lunar cratering distribution and verify the cratering asymmetries generated by the main-belt asteroids (MBAs) as well as the near-Earth objects (NEOs). Based on a planar model that excludes the terrestrial and lunar gravitations on the impactors and assuming the impactor encounter speed with Earth $v_{\rm{enc}}$ is higher than the lunar orbital speed $v_{\rm{M}}$, we rigorously integrated the lunar cratering distribution, and derived its approximation to the first order of $v_{\rm{M}}/v_{\rm{enc}}$. Numerical simulations of lunar bombardment by the MBAs during the late heavy bombardment were performed with an Earth-Moon distance $a_{\rm{M}}$ = 20--60 Earth radii in five cases. The analytical model directly proves the existence of a leading/trailing asymmetry and the absence of near/far asymmetry. The approximate form of the leading/trailing asymmetry is $(1 + A_1 \cos\beta)$, which decreases as the apex distance $\beta$ increases. The numerical simulations show evidence of a pole/equator asymmetry as well as the leading/trailing asymmetry, and the former is empirically described as $(1 + A_2 \cos2\varphi)$, which decreases as the latitude modulus $|\varphi|$ increases. The amplitudes $A_{1,2}$ are reliable measurements of asymmetries. Our analysis explicitly indicates the quantitative relations between cratering distribution and bombardment conditions (impactor properties and the lunar orbital status) like $A_1 \propto v_{\rm{M}}/v_{\rm{enc}}$, resulting in a method for reproducing the bombardment conditions through measuring the asymmetry. Mutual confirmation between analytical model and numerical simulations is found in terms of the cratering distribution and its variation with $a_{\rm{M}}$. Estimates of $A_1$ for crater density distributions generated by the MBAs and the NEOs are 0.101--0.159 and 0.117, respectively.

We determine the proportions of two mixed crater populations distinguishable by size distributions on the Moon. A "multiple power-law" model is built to formulate crater size distribution $N(D) \propto D^{-\alpha}$ whose slope $\alpha$ varies with crater diameter $D$. Fitted size distribution of lunar highland craters characterized by $\alpha = 1.17 \pm 0.04$, $1.88 \pm 0.07$, $3.17 \pm 0.10$ and $1.40 \pm 0.15$ for consecutive $D$ intervals divided by 49, 120 and 251 km and that of lunar Class 1 craters with a single slope $\alpha = 1.96 \pm 0.14$, are taken as Population 1 and 2 crater size distribution respectively, whose sum is then fitted to the size distribution of global lunar craters with $D$ between 10 and 100 km. Estimated crater densities of Population 1 and 2 are $44 \times 10^{-5}$ and $5 \times 10^{-5}$ km$^{-2}$ respectively, leading to the proportion of the latter $10 \%$. The results underlines the need for considering the Population 1 craters and the relevant impactors, the primordial main-belt asteroids which dominated the late heavy bombardment.

We present proper motions for 15 pulsars which are observed regularly by the Nanshan 25-m radio telescope. Two methods, the frequentist method (Coles et al.2011) and the Bayesian (Lentati et al. 2014) method, are used and the results are compared. We demonstrate that the two methods can be applied to young pulsar data sets that exhibit large amounts of timing noise with steep spectral exponents and give consistent results. The measured positions also agree with very-long-baseline interferometric positions. Proper motions for four pulsars are obtained for the first time, and improved values are obtained for five pulsars.

Aims. To make a detailed study of the nulling and subpulse drifting in PSR J1727$-$2739 for investigation of its radiation properties. Methods. The observations were carried out on 20 March, 2004 using the Parkes 64-m radio telescope, with a central frequency of 1518 MHz. A total of 5568 single pulses were analysed. Results. This pulsar shows well defined nulls with lengths lasting from 6 to 281 pulses and separated by burst phases ranging from 2 to 133 pulses. We estimate a nulling fraction of around 68\%. No emission in the average pulse profile integrated over all null pulses is detected with significance above $3\sigma$. Most transitions from nulls to bursts are within a few pulses, whereas the transitions from bursts to nulls exhibit two patterns of decay: decrease gradually or rapidly. In the burst phase, we find that there are two distinct subpulse drift modes with vertical spacing between the drift bands of $9.7 \pm 1.6$ and $5.2 \pm 0.9$ pulse periods, while sometimes there is a third mode with no subpulse drifting. Some mode transitions occur within a single burst, while others are separated by nulls. Different modes have different average pulse profiles. Possible physical mechanisms are discussed.

This paper reports on polarimetric radiation properties based on the switching modes of normal PSR B2020+28 by analysing the data acquired from the Nanshan 25-m radio telescope at 1556 MHz. With nearly 8 hours quasi-continuous observation, the data presented some striking and updated phenomena. The change of relative intensity between the leading and trailing components is the predominant feature of mode switching. The intensity ratio between the leading and trailing components are measured for the individual profiles averaged over 30 seconds. It is found that there is an excess of high ratios over the normal distribution, which indicates that two modes exist in the pulsar. The distribution of abnormal mode has a narrower width indicating that the abnormal mode is more stable than the normal mode. A total of 76 mode switching events are detected in our data. It spends 89\% in the normal mode and 11\% in the abnormal mode. The intrinsic distributions of mode timescales are constrained with power-law distributions. The significant difference in the index of the duration distribution between normal and abnormal modes possibly indicates that the timescale for the abnormal mode to get stable is shorter than that for the normal mode. The frequent switching between both modes may indicate that the oscillations between different magnetospheric states are rapid.

Based on a multiwavelength study, the interstellar medium and young stellar objects (YSOs) around the HII region Sh2-82 have been analyzed. Two molecular clumps were found from the archival data of the Galactic Ring Survey, and using the Two Micron All-Sky Survey catalog, we found two corresponding young clusters embedded in the molecular clumps. The very good relations between CO emission, infrared shells and YSOs suggest that it is probably a triggered star formation region from the expansion of Sh2-82. We further used the data from the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire from Spitzer to study the YSOs within the two clumps, confirming star formation in this region. By spectral energy distribution fits to each YSO candidate with infrared excess, we derived the slope of the initial mass function. Finally, comparing the HII region's dynamical age and the fragmentation time of the molecular shell, we discard the "collect and collapse" process as being the triggering mechanism for YSO formation. Sh2-82 can be a mixture of other processes such as radiative-driven implosion and/or collisions with pre-existing clumps.

In this paper, we investigate the electron Landau-level stability and its influence on the electron Fermi energy, $E_{\rm F}(e)$, in the circumstance of magnetars, which are powered by magnetic field energy. In a magnetar, the Landau levels of degenerate and relativistic electrons are strongly quantized. A new quantity $g_{n}$, the electron Landau-level stability coefficient is introduced. According to the requirement that $g_{n}$ decreases with increasing the magnetic field intensity $B$, the magnetic-field index $\beta$ in the expression of $E_{\rm F}(e)$ must be positive. By introducing the Dirac$-\delta$ function, we deduce a general formulae for the Fermi energy of degenerate and relativistic electrons, and obtain a particular solution to $E_{\rm F}(e)$ in a superhigh magnetic field (SMF). This solution has a low magnetic-field index of $\beta=1/6$, compared with the previous one, and works when $\rho\geq 10^{7}$~g cm$^{-3}$ and $B_{\rm cr}\ll B\leq 10^{17}$~Gauss. By modifying the phase space of relativistic electrons, a SMF can enhance the electron number density $n_e$, and decrease the maximum of electron Landau level number, which results in a redistribution of electrons. According to Pauli exclusion principle, the degenerate electrons will fill quantum states from the lowest Landau level to the highest Landau level. As $B$ increases, more and more electrons will occupy higher Landau levels, though $g_{n}$ decreases with the Landau level number $n$. The enhanced $n_{e}$ in a SMF means an increase in the electron Fermi energy and an increase in the electron degeneracy pressure. The results are expected to facilitate the study of the weak-interaction processes inside neutron stars and the magnetic-thermal evolution mechanism for megnetars.

Based on the basic definition of Fermi energy of degenerate and relativistic electrons, we obtain a special solution to electron Fermi energy, $E_{\rm F}(e)$, and express $E_{\rm F}(e)$ as a function of electron fraction, $Y_{e}$, and matter density, $\rho$. Several useful analytical formulae for $Y_{e}$ and $\rho$ within classical models and the work of Dutra et al. 2014 (Type-2) in relativistic mean field theory are obtained using numerically fitting. When describing the mean-field Lagrangian, density, we adopt the TMA parameter set, which is remarkably consistent with with the updated astrophysical observations of neutron stars. Due to the importance of the density dependence of the symmetry energy, $S$, in nuclear astrophysics, a brief discussion on the symmetry parameters $S_v$ and $L$ (the slope of $S$) is presented. Combining these fit formulae with boundary conditions for different density regions, we can evaluate the value of $E_{\rm F}(e)$ in any given matter density, and obtain a schematic diagram of $E_{\rm F}(e)$ as a continuous function of $\rho$. Compared with previous study on the electron Fermi energy in other models, our methods of calculating $E_{\rm F}(e)$ are more simple and convenient, and can be universally suitable for the relativistic electron regions in the circumstances of common neutron stars. We have deduced a general expression of $E_{\rm F}(e)$ and $n_{e}$, which could be used to indirectly test whether one EoS of a NS is correct in our future studies on neutron star matter properties. Since URCA reactions are expected in the center of a massive star due to high-value electron Fermi energy and electron fraction, this study could be useful in the future studies on the NS thermal evolution.

By introducing Dirac ?-function in superhigh magnetic field, we deduce a general formula for pressure of degenerate and relativistic electrons, Pe, which is suitable for superhigh magnetic fields, discuss the quantization of Landau levels of electrons, and consider the quantum electrodynamic(QED) effects on the equations of states (EOSs) for different matter systems. The main conclusions are as follows: the stronger the magnetic field strength, the higher the electron pressure becomes; compared with a common radio pulsar, a magnetar could be a more compact oblate spheroid-like deformed neutron star due to the anisotropic total pressure; and an increase in the maximum mass of a magnetar is expected because of the positive contribution of the magnetic field energy to the EOS of the star. Since this is an original work in which some uncertainties could exist, to further modify and perfect our theory model should be considered in our future studies.

Due to the lack of long term pulsed emission in quiescence and the strong timing noise, it is impossible to directly measure the braking index $n$ of a magnetar. Based on the estimated ages of their potentially associated supernova remnants (SNRs), we estimate the values of the mean braking indices of eight magnetars with SNRs, and find that they cluster in a range of $1\sim$42. Five magnetars have smaller mean braking indices of $1<n<3$, and we interpret them within a combination of magneto-dipole radiation and wind aided braking, while the larger mean braking indices of $n>3$ for other three magnetars are attributed to the decay of external braking torque, which might be caused by magnetic field decay. We estimate the possible wind luminosities for the magnetars with $1<n<3$, and the dipolar magnetic field decay rates for the magnetars with $n>3$ within the updated magneto-thermal evolution models. Although the constrained range of the magnetars' braking indices is tentative, due to the uncertainties in the SNR ages, which come from distance uncertainties and the unknown conditions of the expanding shells, our method provides an effective way to constrain the magnetars' braking indices if the measurements of the SNRs' ages are reliable, which can be improved by future observations.

Observations show that there is a proton spectral "break" with E$_{break}$ at 1-10MeV in some large CME-driven shocks. Theoretical model usually attribute this phenomenon to a diffusive shock acceleration. However, the underlying physics of the shock acceleration still remains uncertain. Although previous numerical models can hardly predict this "break" due to either high computational expense or shortcomings of current models, the present paper focuses on simulating this energy spectrum in converged two shocks by Monte Carlo numerical method. Considering the Dec 13 2006 CME-driven shock interaction with an Earth bow shock, we examine whether the energy spectral "break" could occur on an interaction between two shocks. As result, we indeed obtain the maximum proton energy up to 10MeV, which is the premise to investigate the existence of the energy spectral "break". Unexpectedly, we further find a proton spectral "break" appears distinctly at the energy $\sim$5MeV.

Observations from multiple spacecraft show that there are energy spectral "breaks" at 1-10MeV in some large CME-driven shocks. However, numerical models can hardly simulate this property due to high computational expense. The present paper focuses on analyzing these energy spectral "breaks" by Monte Carlo particle simulations of an isolated CME-driven shock. Taking the Dec 14 2006 CME-driven shock as an example, we investigate the formation of this energy spectral property. For this purpose, we apply different values for the scattering time in our isolated shock model to obtain the highest energy "tails", which can potentially exceed the "break" energy range. However, we have not found the highest energy "tails" beyond the "break" energy range, but instead find that the highest energy "tails" reach saturation near the range of energy at 5MeV. So, we believe that there exists an energy spectral "cut off" in an isolated shock. If there is no interaction with another shock, there would not be formation of the energy spectral "break" property.

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.

Anisotropic bursts of gravitational radiation produced by events such as super-massive black hole mergers leave permanent imprints on space. Such gravitational wave "memory" (GWM) signals are, in principle, detectable through pulsar timing as sudden changes in the apparent pulse frequency of a pulsar. If an array of pulsars is monitored as a GWM signal passes over the Earth, the pulsars would simultaneously appear to change pulse frequency by an amount that varies with their sky position in a quadrupolar fashion. Here we describe a search algorithm for such events and apply the algorithm to approximately six years of data from the Parkes Pulsar Timing Array. We find no GWM events and set an upper bound on the rate for events which could have been detected. We show, using simple models of black hole coalescence rates, that this non-detection is not unexpected.

By taking into account the surface diffuseness correction for unstable nuclei, the accuracy of the macroscopic-microscopic mass formula is further improved. The rms deviation with respect to essentially all the available mass data falls to 298 keV, crossing the 0.3 MeV accuracy threshold for the first time within the mean-field framework. Considering the surface effect of the symmetry potential which plays an important role in the evolution of the "neutron skin" toward the "neutron halo" of nuclei approaching the neutron drip line, we obtain an optimal value of the symmetry energy coefficient J=30.16 MeV. With an accuracy of 258 keV for all the available neutron separation energies and of 237 keV for the alpha-decay Q-values of super-heavy nuclei, the proposed mass formula is particularly important not only for the reliable description of the r-process of nucleosynthesis but also for the study of the synthesis of super-heavy nuclei.

In a superhigh magnetic field, direct Urca reactions can proceed for an arbitrary proton concentration. Since only the electrons with high energy $E$ ($E > Q$, $Q$ is the threshold energy of inverse $\beta-$decay) at large Landau levels can be captured, we introduce the Landau level effect coefficient $q$ and the effective electron capture rate $\Gamma_{\rm eff}$. By using $\Gamma_{\rm eff}$, the values of $L_{\rm X}$ and $L_{\rm \nu}$ are calculated, where and $L_{\rm \nu}$, $L_{\rm X}$ are the average neutrino luminosity of Anomalous X-ray Pulsars (AXPs) and the average X-ray luminosity of AXPs, respectively. The complete process of electron capture inside a magnetar is simulated numerically.

In this paper, we modify our previous research carefully, and derive a new expression of electron energy density in superhigh magnetic fields. Based on our improved model, we re-compute the electron capture rates and the magnetic fields' evolutionary timescales $t$ of magnetars. According to the calculated results, the superhigh magnetic fields may evolve on timescales $\sim (10^{6}-10^{7})$ yrs for common magnetars, and the maximum timescale of the field decay, $t\approx 2.9507 \times 10^{6}$ yrs, corresponding to an initial internal magnetic field $B_{\rm 0}= 3.0 \times 10^ {15}$ G and an initial inner temperature $T_{\rm 0}= 2.6 \times 10^ {8}$ K. Motivated by the results of the neutron star-supernova remnant(SNR) association of Zhang $\&$ Xie(2011), we calculate the maximum $B_{\rm 0}$ of magnetar progenitors, $B_{\rm max}\sim (2.0\times 10^{14}-2.93 \times 10^{15})$ G when $T_{\rm 0}= 2.6 \times 10^ {8}$ K. When $T_{\rm 0}\sim 2.75 \times 10^ {8}-~1.75 \times 10^ {8}$ K, the maximum $B_{\rm 0}$ will also be in the range of $\sim 10^{14}-10^{15}$ G, not exceeding the upper limit of magnetic field of a magnetar under our magnetar model. We also investigate the relationship between the spin-down ages of magnetars and the ages of their SNRs, and explain why all AXPs associated with SNRs look older than their real ages, whereas all SGRs associated with SNRs appear younger than they are.

As soon as the energy of electrons near the Fermi surface are higher than $Q$, the threshold energy of inverse $\beta-$ decay, the electron capture process will dominate. The resulting high-energy neutrons will destroy anisotropic ${}^3P_2$ neutron superfluid Cooper pairs. By colliding with the neutrons produced in the process $n+ (n\uparrow n\downarrow)\longrightarrow n+ n+ n$, the kinetic energy of the outgoing neutrons will be transformed into thermal energy. The transformed thermal energy would transported from the star interior to the star surface by conduction, then would be transformed into radiation energy as soft X-rays and gamma-rays. After a highly efficient modulation within the pulsar magnetosphere, the surface thermal emission (mainly soft X/$\gamma$-ray emission) has been shaped into a spectrum with the observed characteristics of magnetars. By introducing two important parameters: Landau level-superfluid modified factor and effective X/$\gamma$-ray coefficient, we numerically simulate the process of magnetar cooling and magnetic field decay, and then compute magnetars' soft X/$\gamma$-ray luminosities $L_{X}$. Further, we obtain aschematic diagrams of $L_{X}$ as a function of magnetic field strength $B$. The observations are compared with the calculations.

In this paper, we consider the effect of Landau levels on the decay of superhigh magnetic fields of magnetars. Applying ${}^3P_2$ anisotropic neutron superfluid theory yield a second-order differential equation for a superhigh magnetic field $B$ and its evolutionary timescale $t$. The superhigh magnetic fields may evolve on timescales $\sim (10^{6}-10^{7})$ yrs for common magnetars. According to our model, the activity of a magnetar may originate from instability caused by the high electron Fermi energy.

Assuming that the timescale of the magnetic field decay is approximately equal to that of the stellar cooling via neutrino emission, we obtain a one-to-one relationship between the effective surface thermal temperature and the inner temperature. The ratio of the effective neutrino luminosity to the effective X-ray luminosity decreases with decaying magnetic field.

Pulsar nulling is a phenomenon of sudden cessation of pulse emission for a number of periods. The nulling fraction was often used to characterize the phenomenon. We propose a new method to analyse pulsar nulling phenomenon, by involving two key parameters, the nulling degree, $\chi$, which is defined as the angle in a rectangular coordinates for the numbers of nulling periods and bursting periods, and the nulling scale, $ N $, which is defined as the effective length of the consecutive nulling periods and bursting periods. The nulling degree $\chi$ can be calculated by $\tan \chi = N_{\rm nulling} / N_{\rm bursting} $ and the mean is related to the nulling fraction, while the nulling scale, $ N $, is also a newly defined fundamental parameter which indicates how often the nulling occurs. We determined the distributions of $\chi$ and $ N $ for 10 pulsars by using the data in literature. We found that the nulling degree $\chi$ indicates the relative length of nulling to that of bursting, and the nulling scale $ N $ is found to be related to the derivative of rotation frequency and hence the loss rate of rotational energy of pulsars. Their deviations reflect the randomness of the nulling process.

In this paper, we propose a new method to use the strong lensing data sets to constrain a cosmological model. By taking the ratio $\mathcal{D}^{obs}_{ij}=\theta_{\mathrm{E_{\mathrm{i}}}}\sigma_{\mathrm{0_{\mathrm{j}}}}^2/\theta_{\mathrm{E_{\mathrm{j}}}}\sigma_{\mathrm{0_{\mathrm{i}}}}^2$ as cosmic observations, one can {\it completely} eliminate the uncertainty caused by the relation $\sigma_{\mathrm{SIS}}=f_{\mathrm{E}}\sigma_0$ which characterizes the relation between the stellar velocity dispersion $\sigma_0$ and the velocity dispersion $\sigma_{SIS}$. Via our method, a relative tight constraint to the cosmological model space can be obtained, for the spatially flat $\Lambda$CDM model as an example $\Omega_m=0.143_{- 0.143-0.143-0.143}^{+ 0.000769+0.143+0.489}$ in $3\sigma$ regions. And by using this method, one can also probe the nature of dark energy and the spatial curvature of our Universe.

In this paper, an approximate method of calculating the Fermi energy of electrons ($E_{F}(e)$) in a high-intensity magnetic field, based on the analysis of the distribution of a neutron star magnetic field, has been proposed. In the interior of a Neutron star, different forms of intense magnetic field could exist simultaneously and a high electron Fermi energy could be generated by the release of magnetic field energy. The calculation results show that: $E_{F}(e)$ is related to density $\rho$, the mean electron number per baryon $Y_{e}$ and magnetic field strength $B$.

A statistical study was made on both the energy dependence of the low-frequency quasi-periodic oscillation (QPO) amplitude (LFQPO amplitude spectrum) and the LFQPO amplitude from all the RXTE observations of GRS 1915+105. Based on the two-branch correlation of the LFQPO frequency and the hardness ratio, the observations that were suitable for evaluating the LFQPO amplitude spectrum were divided into two groups. According to a comparison between the radio and X-ray emissions, we deduced that the jets during the two groups of observations are very different. A negative correlation between the LFQPO frequency and the radio flux was found for one group. The LFQPO amplitude spectrum was fitted by a power-law with an exponential cutoff in order to describe it quantitatively. It reveals that as the LFQPO frequency increases, the power-law hardens. And the cutoff energy firstly decreases, and then smoothly levels off. The fit also shows that the LFQPO amplitude spectra of the two groups are essentially the same, suggesting that the LFQPO seems not to originate from the jet. The LFQPO amplitude spectra are hard, indicating a possible origin of the LFQPO in the corona. As the LFQPO frequency increases, the LFQPO amplitude firstly increases and then decreases. The effects of the low pass filter and the jet on the LFQPO amplitude are discussed.

Using data from the ATNF pulsar catalogue, 186 binary pulsars are shown in the magnetic field versus spin period (B-P) diagram, and their relationship to the spin-up line is investigated. Generally speaking, pulsars in binary systems should be below the spin-up line when they get enough accretion mass from their companions. It is found that there are seven binary pulsars above the spin-up line. Based on the parameters of these seven binary systems, we describe possible reasons why they are above the spin-up line.

186 binary pulsars are shown in the magnetic field versus spin period (B-P) dia- gram, and their relations to the millisecond pulsars can be clearly seen. We declaim a minimum accretion rate for the millisecond pulsar formation both from the observation and theory. If the accretion rate is lower than the minimum accretion rate, the pulsar in binary system will not become a millisecond pulsar after the evolution.

We made a phase-resolved timing analysis of GRS 1915+105 in its {\rho} state and obtained detailed {\rho} cycle evolutions of the frequency, the amplitude and the coherence of low-frequency quasi-periodic oscillation (LFQPO). We combined our timing results with the spectral study by Neilsen et al. to perform an elaborate comparison analysis. Our analyses show that the LFQPO frequency does not scale with the inner disk radius, but it is related to the spectral index, indicating a possible correlation between the LFQPOs and the corona. The LFQPO amplitude spectrum and other results are naturally explained by tying the LFQPO to the corona. The similarities of the spectra of variability parameters between the LFQPO from {\rho} state and those from more steady states indicate that the LFQPOs of GRS 1915+105 in very different states seem to share the same origin.

Properties of X-ray luminosities in low-mass X-ray binaries (LMXBs) mainly depend on donors. We have carried out a detailed study of donors in persistent neutron-star LMXBs (PLMXBs) by means of a population synthesis code. PLMXBs with different donors have different formation channels. Our numerical simulations show that more than 90% of PLMXBs have main sequence (MS) donors, and PLMXBs with red giant (RG) donors via stellar wind (Wind) are negligible. In our model, most of neutron stars (NSs) in PLMXBs with hydrogen-rich donors form via core-collapse supernovae, while more than 90% of NSs in PLMXBs with naked helium star (He) donors or white dwarf (WD) donors form via an evolution-induced collapse via helium star ($1.4 \leq M_{\rm He}/M_\odot \leq 2.5$) or an accretion-induced collapses for an accreting ONeMg WD.

We present astrometric results on two millisecond pulsars, PSR B1257+12 and PSR J1022+1001, as carried out through VLBI. For PSR B1257+12, a model-independent distance of $710_{\rm -38}^{\rm +43}$ pc and proper motion of ($\mu_{\rm \alpha}=46.44\pm0.08$ mas/yr, $\mu_{\rm \delta}=-84.87\pm0.32$ mas/yr) were obtained from 5 epochs of VLBA and 4 epochs of EVN observations, spanning about 2 years. The two dimensional proper motion of PSR J1022+1001 ($\mu_{\rm \alpha} \sim -10.13$ mas/yr, $\mu_{\delta} \sim 16.89$ mas/yr) was also estimated, using 3 epochs of EVN observations. Based on our results, the X-ray efficiency of PSR B1257+12 should be in the same range as other millisecond pulsars, and not as low as previously thought.

Gravitational wave bursts produced by supermassive binary black hole mergers will leave a persistent imprint on the space-time metric. Such gravitational wave memory signals are detectable by pulsar timing arrays as a glitch event that would seem to occur simultaneously for all pulsars. In this paper, we describe an initial algorithm which can be used to search for gravitational wave memory signals. We apply this algorithm to the Parkes Pulsar Timing Array data set. No significant gravitational wave memory signal is founded in the data set.

We report on the properties of strong pulses from PSR B0656+14 by analyzing the data obtained using Urumqi 25-m radio telescope at 1540 MHz from August 2007 to September 2010. In 44 hrs of observational data, a total of 67 pulses with signal-to-noise ratios above a 5-{\sigma} threshold were detected. The peak flux densities of these pulses are 58 to 194 times that of the average profile, and the pulse energies of them are 3 to 68 times that of the average pulse. These pulses are clustered around phases about 5 degrees ahead of the peak of the average profile. Comparing with the width of the average profile, they are relatively narrow, with the full widths at half-maximum range from 0.28 to 1.78 degrees. The distribution of pulse-energies of the pulses follows a lognormal distribution. These sporadic strong pulses detected from PSR B0656+14 are different in character from the typical giant pulses, and from its regular pulses.

Symbiotic X-ray binaries (SyXBs) comprise a rare class of low-mass X-ray binaries. We study the Galactic SyXBs, which we consider as detached binaries composed of low-mass giants and wind-fed neutron star companions, by simulation of the interaction of a magnetized neutron star (NS) with its environment and utilizing a population synthesis code. We focus mainly on the parameters that influence observational appearance of the SyXB: the donor wind velocity (vw) and the angular momentum distribution in the shell of matter settling onto NS. We estimate the birthrate of SyXB as $\sim 4.1\times 10^{-5}$ yr$^{-1}$ to $ \sim 6.6\times 10^{-6}$ yr$^{-1}$ and their number in the Galaxy as $\sim$(100 -- 1000). Assumed stellar wind velocity from cool giants is the input parameter that influences the model SyXBs population most. Among known SyXBs or candidate systems, 4U 1954+31 and IGR J16358-4724 in which NS have very long spin periods may host quasi-spherically accreting NSs. GX 1+4 has a peculiar long-term spin behaviour and it may also be a quasi-spherical wind-accreting source. We cannot identify whether there are wind-fed accretion disks in 4U 1700+24, Sct X-1, IRXS J180431.1-273932 and 2XMM J174016.0-290337.

Context. The early phase of the coalescence of supermassive black hole (SMBH) binaries from their host galaxies provides a guaranteed source of low-frequency (nHz-$\mu$Hz) gravitational wave (GW) radiation by pulsar timing observations. These types of GW sources would survive the coalescing and be potentially identifiable. Aims. We aim to provide an outline of a new method for detecting GW radiation from individual SMBH systems based on the Sloan Digital Sky Survey (SDSS) observational results, which can be verified by future observations. Methods. Combining the sensitivity of the international Pulsar Timing Array (PTA) and the Square Kilometer Array (SKA) detectors, we used a binary population synthesis (BPS) approach to determine GW radiation from close galaxy pairs under the assumption that SMBHs formed at the core of merged galaxies. We also performed second post-Newtonian approximation methods to estimate the variation of the strain amplitude with time. Results. We find that the value of the strain amplitude \emph{h} varies from about $10^{-14}$ to $10^{-17}$ using the observations of 20 years, and we estimate that about 100 SMBH sources can be detected with the SKA detector.

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.

In order to depict the quantization of Landau levels, we introduce Dirac $\delta$ function, and gain a concise expression for the electron Fermi energy, $E_{F}(e) \propto B^{1/4}$. The high soft X-ray luminosities of magnetars may be naturally explained by our theory.

The mode switching phenomenon of PSR B0329+54 is investigated based on the long-term monitoring from September 2003 to April 2009 made with the Urumqi 25m radio telescope at 1540 MHz. At that frequency, the change of relative intensity between the leading and trailing components is the predominant feature of mode switching. The intensity ratios between the leading and trailing components are measured for the individual profiles averaged over a few minutes. It is found that the ratios follow normal distributions, where the abnormal mode has a wider typical width than the normal mode, indicating that the abnormal mode is less stable than the normal mode. Our data show that 84.9% of the time for PSR B0329+54 was in the normal mode and 15.1% was in the abnormal mode. From the two passages of eight-day quasi-continuous observations in 2004, and supplemented by the daily data observed with 15 m telescope at 610 MHz at Jodrell Bank Observatory, the intrinsic distributions of mode timescales are constrained with the Bayesian inference method. It is found that the gamma distribution with the shape parameter slightly smaller than 1 is favored over the normal, lognormal and Pareto distributions. The optimal scale parameters of the gamma distribution is 31.5 minutes for the abnormal mode and 154 minutes for the normal mode. The shape parameters have very similar values, i.e. 0.75^{+0.22}_{-0.17} for the normal mode and 0.84^{+0.28}_{-0.22} for the abnormal mode, indicating the physical mechanisms in both modes may be the same. No long-term modulation of the relative intensity ratios was found for both the modes, suggesting that the mode switching was stable. The intrinsic timescale distributions, for the first time constrained for this pulsar, provide valuable information to understand the physics of mode switching.

Using the data from the Rossi X-Ray Timing Explorer satellite, we investigate the spectral evolution along a "Z" track and a "v" track on the hardness-intensity diagrams of the first transient Z source XTE J1701-462. The spectral analyses suggest that the inner disk radius depends on the mass accretion rate, in agreement with the model prediction, R_in \propto ((dM/dt)_disk)^{2/7}, for a radiation pressure dominated accretion disk interacting with the magnetosphere of a neutron star (NS). The changes in the disk mass accretion rate (dM/dt)_disk are responsible for the evolution of the "Z" or "v" track. The radiation pressure thickens the disk considerably, and also produces significant outflows. The NS surface magnetic field strength, derived from the interaction between the magnetosphere and the radiation pressure dominated accretion disk, is ~(1--3)X10^9 G, which is possibly between normal atoll and Z sources. A significant hard tail is detected in the horizontal branches and we discuss several possible origins of the hard tail.

We report on variations in the mean position angle of the 20 millisecond pulsars being observed as part of the Parkes Pulsar Timing Array (PPTA) project. It is found that the observed variations are dominated by changes in the Faraday rotation occurring in the Earth's ionosphere. Two ionospheric models are used to correct for the ionospheric contribution and it is found that one based on the International Reference Ionosphere gave the best results. Little or no significant long-term variation in interstellar RM was found with limits typically about 0.1 rad m$^{-2}$ yr$^{-1}$ in absolute value. In a few cases, apparently significant RM variations over timescales of a few 100 days or more were seen. These are unlikely to be due to localised magnetised regions crossing the line of sight since the implied magnetic fields are too high. Most probably they are statistical fluctuations due to random spatial and temporal variations in the interstellar electron density and magnetic field along the line of sight.

We analyze individual pulses detected from RRAT J1819$-$1458. From April 2007 to April 2010, we carried out observations using the Nanshan 25-m radio telescope of Urumqi Observatory at a central frequency of 1541.25 MHz. We obtain a dispersion measure $DM=195.7\pm0.3$ pc cm^{-3} by analyzing all the 423 detected bursts. The tri-band pattern of arrival time residuals is confirmed by a single pulse timing analysis. Twenty-seven bimodal bursts located in the middle residual band are detected, and, profiles of two typical bimodal bursts and two individual single-peak pulses are presented. We determine the statistical properties of SNR and W$_{50}$ of bursts in different residual bands. The W$_{50}$ variation with SNR shows that the shapes of bursts are quite different from each other. The cumulative probability distribution of intensity for a possible power law with index $\alpha=1.6\pm0.2$ is inferred from the number of those bursts with $SNR\ge6$ and high intensities.

Polarization profiles are presented for 20 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array project. The observations used the Parkes multibeam receiver with a central frequency of 1369 MHz and the Parkes digital filterbank pulsar signal-processing system PDFB2. Because of the large total observing time, the summed polarization profiles have very high signal/noise ratios and show many previously undetected profile features. Thirteen of the 20 pulsars show emission over more than half of the pulse period. Polarization variations across the profiles are complex and the observed position angle variations are generally not in accord with the rotating-vector model for pulsar polarization. Never-the-less, the polarization properties are broadly similar to those of normal (non-millisecond) pulsars, suggesting that the basic radio emission mechanism is the same in both classes of pulsar. The results support the idea that radio emission from millisecond pulsars originates high in the pulsar magnetosphere, probably close to the emission regions for high-energy X-ray and gamma-ray emission. Rotation measures were obtained for all 20 pulsars, eight of which had no previously published measurements.

In this paper, we report a detailed investigation of the multiwavelength properties of a newly detected gamma-ray pulsar, PSR J2021+4026, in both observational and theoretical aspects. We firstly identify an X-ray source in the XMM-Newton serendipitous source catalogue, 2XMM J202131.0+402645, located within the 95% confidence circle of PSR J2021+4026. With an archival Chandra observation, this identification provides an X-ray position with arcsecond accuracy which is helpful in facilitating further investigations. Searching for the pulsed radio emission at the position of 2XMM J202131.0+402645 with a 25-m telescope at Urumqi Astronomical Observatory resulted in null detection and places an upper-limit of 0.1~mJy for any pulsed signal at 18~cm. Together with the emission properties in X-ray and gamma-ray, the radio quietness suggests PSR J2021+4026 to be another member of Geminga-like pulsars. In the radio sky survey data, extended emission features have been identified in the gamma-ray error circle of PSR J2021+4026. We have also re-analyzed the gamma-ray data collected by FERMI's Large Area Telescope. We found that the X-ray position of 2XMM J202131.0+402645 is consistent with that of the optimal gamma-ray timing solution. We have further modeled the results in the context of outer gap model which provides us with constraints for the pulsar emission geometry such as magnetic inclination angle and the viewing angle. We have also discussed the possibility of whether PSR J2021+4026 has any physical association with the supernova remnant G78.2+2.1 (gamma-Cygni).

Seven years of pulse time-of-arrival measurements have been collected from observations of the young pulsar PSR B2334+61 using the Nanshan radio telescope of Urumqi Observatory. A phase-connected solution has been obtained over the whole data span, 2002 August to 2009 August. This includes a very large glitch that occurred between 2005 August 26 and September 8 (MJDs 53608 and 53621). The relative increase in rotational frequency for this glitch, $\Delta\nu_{g}/\nu~\sim~20.5\times10^{-6}$, is the largest ever seen. Although accounting for less than 1\% of the glitch, there were two well-defined exponential decay terms with time constants of 21 and 147 days respectively. There was also a large long-term increase in the spindown rate with $\Delta\dot\nu_p/\dot\nu \sim 0.011$ at the time of the glitch. A highly significant oscillation with a period of close to one year is seen in the post-glitch residuals. It is very unlikely that this can be accounted for by a pulsar position error or proper motion -- it appears to result from effects interior to the neutron star. Implications of these results for pulsar glitch models are discussed.

Glitches detected in pulsar timing observations at the Nanshan radio telescope of Urumqi Observatory between 2002 July and 2008 December are presented. In total, 29 glitches were detected in 19 young pulsars, with this being the first detection of a glitch in 12 of these pulsars. Fractional glitch amplitudes range from a few parts in 10^{-11} to 3.9 x 10^{-6}. Three "slow" glitches are identified in PSRs J0631+1036, B1822-09 and B1907+10. Post-glitch recoveries differ greatly from pulsar to pulsar and for different glitches in the same pulsar. Most large glitches show some evidence for exponential post-glitch recovery on timescales of 100 -- 1000 days, but in some cases, e.g., PSR B1758-23, there is little or no recovery. Beside exponential recoveries, permanent increases in slowdown rate are seen for the two large glitches in PSRs B1800-21 and B1823-13. These and several other pulsars also show a linear increase in nudot following the partial exponential recovery, which is similar to the Vela pulsar post-glitch behaviour. Analysis of the whole sample of known glitches show that fractional glitch amplitudes are correlated with characteristic age with a peak at about 10^5 years, but there is a spread of two or three orders of magnitude at all ages. Glitch activity is positively correlated with spin-down rate, again with a wide spread of values.

We present an analysis of timing residual (noise) of 54 pulsars obtained from 25-m radio telescope at Urumqi Observatory with a time span of 5~8 years, dealing with statistics of the Hurst parameter. The majority of these pulsars were selected to have timing noise that look like white noise rather than smooth curves. The results are compared with artificial series of different constant pairwise covariances. Despite the noise like appearance, many timing residual series showed Hurst parameters significantly deviated from that of independent series. We concluded that Hurst parameter may be capable of detecting dependence in timing residual and of distinguishing chaotic behavior from random processes.

By assuming an aspherical stellar wind with an equatorial disk from a red giant, we investigate the production of Type Ia supernovae (SNe Ia) via symbiotic channel. We estimate that the Galactic birthrate of SNe Ia via symbiotic channel is between $1.03\times 10^{-3}$ and $2.27\times 10^{-5}$ yr$^{-1}$, the delay time of SNe Ia has wide range from $\sim$ 0.07 to 5 Gyr. The results are greatly affected by the outflow velocity and mass-loss rate of the equatorial disk. Using our model, we discuss the progenitors of SN 2002ic and SN 2006X.

We describe a comprehensive pulsar monitoring campaign for the Large Area Telescope (LAT) on the {\em Fermi Gamma-ray Space Telescope} (formerly GLAST). The detection and study of pulsars in gamma rays give insights into the populations of neutron stars and supernova rates in the Galaxy, into particle acceleration mechanisms in neutron star magnetospheres, and into the "engines" driving pulsar wind nebulae. LAT's unprecedented sensitivity between 20 MeV and 300 GeV together with its 2.4 sr field-of-view makes detection of many gamma-ray pulsars likely, justifying the monitoring of over two hundred pulsars with large spin-down powers. To search for gamma-ray pulsations from most of these pulsars requires a set of phase-connected timing solutions spanning a year or more to properly align the sparse photon arrival times. We describe the choice of pulsars and the instruments involved in the campaign. Attention is paid to verifications of the LAT pulsar software, using for example giant radio pulses from the Crab and from PSR B1937+21 recorded at Nan\c{c}ay, and using X-ray data on PSR J0218+4232 from XMM-Newton. We demonstrate accuracy of the pulsar phase calculations at the microsecond level. Data Table 1 is only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ .

Since mid-2005, a pulsar searching system has been operating at 18 cm on the 25-m radio telescope of Urumqi Observatory. Test observations on known pulsars show that the system can perform the intended task. The prospect of using this system to observe 3EG sources and other target searching tasks is discussed.

This Letter is to investigate the physics of a newly discovered phenomenon -- contracting flare loops in the early phase of solar flares. In classical flare models, which were constructed based on the phenomenon of expansion of flare loops, an energy releasing site is put above flare loops. These models can predict that there is a vertical temperature gradient in the top of flare loops due to heat conduction and cooling effects. Therefore, the centroid of an X-ray looptop source at higher energy bands will be higher in altitude, for which we can define as normal temperature distribution. With observations made by {\it RHESSI}, we analyzed 10 M- or X-class flares (9 limb flares). For all these flares, the movement of looptop sources shows an obvious U-shaped trajectory, which we take as the signature of contraction-to-expansion of flare loops. We find that, for all these flares, normal temperature distribution does exist, but only along the path of expansion. The temperature distribution along the path of contraction is abnormal, showing no spatial order at all. The result suggests that magnetic reconnection processes in the contraction and expansion phases of these solar flares are different.

We report the timing-analysis results obtained for RRAT J1819-1458 from regular timing observations at 1.54 GHz using the Urumqi 25 m radio telescope between 2007 April to 2008 March. RRAT J1819-1458 is a relatively young and highly magnetized neutron star discovered by its sporadic short bursts in the Parkes Multibeam Pulsar Survey data. In 94 hrs of observation data, we detected a total of 162 dispersed bursts of RRAT J1819-1458 with the signal-to-noise ratios (S/Ns) above 5-$\sigma$ threshold. Among them, 5 bursts clearly show two-component structure. The S/N of the strongest burst is 13.3. The source's DM measured through our data is 196.0 $\pm$0.4 pc cm$^{-3}$. The timing position, frequency and its first derivative were determined using standard pulsar timing techniques on the arrival times of these individual bursts. The accuracy of the solved rotating parameters are improved comparing with that in previous publication. Our timing position with 2-$\sigma$ error is consistent with the position of its X-ray counterpart CXOU J181934.1-145804. The effect of timing noise and the phase fluctuation of the individual short bursts on the timing residuals is briefly discussed. The distribution of the timing residuals is bimodal, which cannot be explained readily by timing irregularity.

Quasi-continuous observations of PSR B03239+54 over 20 days using the Nanshan 25-m telescope at 1540 MHz have been used to study the effects of refractive scintillation on the pulsar flux density and diffractive scintillation properties. Dynamic spectra were obtained from datasets of 90 min duration and diffractive parameters derived from a two-dimensional auto-correlation analysis. Secondary spectra were also computed but these showed no significant evidence for arc structure. Cross correlations between variations in the derived parameters were much lower than predicted by thin screen models and in one case was of opposite sign to the prediction. Observed modulation indices were larger than predicted by thin screen models with a Kolmogorov fluctuation spectrum. Structure functions were computed for the flux density, diffractive timescale and decorrelation bandwidth. These indicated a refractive timescale of $8\pm 2$ h, much shorter than predicted by the thin screen model. The measured structure-function slope of $0.4\pm 0.2$ is also inconsistent with scattering by a single thin screen for which a slope of 2.0 is expected. All observations are consistent with scattering by an extended medium having a Kolmogorov fluctuation spectrum which is concentrated towards the pulsar. This interpretation is also consistent with recent observations of multiple diffuse scintillation arcs for this pulsar.

Six glitches have been recently observed in the rotational frequency of the young pulsar PSR B1737-30 (J1740-3015) using the 25-m Nanshan telescope of Urumqi Observatory. With a total of 20 glitches in 20 years, it is one of the most frequently glitching pulsars of the about 1750 known pulsars. Glitch amplitudes are very variable with fractional increases in rotation rate ranging from 10^{-9} to 10^{-6}. Inter-glitch intervals are also very variable, but no relationship is observed between interval and the size of the preceding glitch. There is a persistent increase in |\dot\nu|, opposite in sign to that expected from slowdown with a positive braking index, which may result from changes in the effective magnetic dipole moment of the star during the glitch.

With its relatively long observation time per pointing, the Parkes multibeam survey was effective in detecting nulling pulsars. We have made 2-hour observations of 23 pulsars which showed evidence for pulse nulling or mode changing in the survey data. Because of the low flux density of these pulsars, in most cases averaging times of between 10 s and 60 s were necessary and so this analysis is insensitive to very short nulls. Seven of the pulsars had null fractions of more than 40% with the largest having a lower limit of 95%. Mode changes were observed in six pulsars with clear relationships between nulling and mode changing in some cases. Combined with earlier results, the data suggest that large null fractions are more related to large characteristic age than to long pulse period. The observations suggest that nulling and mode changing are different manifestations of the same phenomenon.

We have measured the positions of 74 pulsars from regular timing observations using the Nanshan radio telescope at Urumqi Observatory between 2000 January and 2004 August (MJD 51500 -- 53240). Proper motions were determined for these pulsars by comparing their current positions with positions given in pulsar catalogues. We compare our results to earlier measurements in the literature and show that, in general, the values agree. New or improved proper motions are obtained for 16 pulsars. The effect of period fluctuations and other timing noise on the determination of pulsar positions is investigated. For our sample, the mean and rms transverse velocities are 443 and 224 km/s respectively, agreeing with previous work even though we determine distances using the new NE2001 electron density model.

We report here on extensive radio observations of the binary system containing PSR B1259-63 and the Be star SS 2883, made around the time of the 2004 periastron. This is the fourth periastron to have been observed in detail. As in previous observations, changes in the pulsar's dispersion and rotation measures are detected over a period spanning 200 days. We show that the eclipse of the pulsed emission lasts from 16 days prior to periastron to 15 days after periastron and is consistent from one periastron to the next. We demonstrate that the timing solution proposed by Wang et al. (2004) provides a good fit through the 2004 periastron. The light curve of the transient unpulsed radio emission is broadly similar from one periastron to the next. For this periastron, however, the light curve is strongly peaked post-periastron with rather little enhancement prior to periastron, in contrast to the 2000 periastron where the peak flux densities were more equal. These observations remain consistent with the interpretation that the pulsar passes through the dense circumstellar disk of the Be star just before and just after periastron. The observed differences from one periastron to the next can be ascribed to variations in the local disk density and magnetic field structure at the time the pulsar enters the disk.

From 2001 January to 2002 June, we monitored PSRs B0329+54, B0823+26, B1929+10, B2020+28 and B2021+51 using the Nanshan 25-m radio telescope of Urumqi Observatory to study their diffractive interstellar scintillation (DISS). The average interval between observations was about 9 days and the observation duration ranged between 2 and 6 hours depending on the pulsar. Wide variations in the DISS parameters were observed over the 18-month data span. Despite this, the average scintillation velocities are in excellent agreement with the proper motion velocities. The average two-dimensional autocorrelation function for PSR B0329+54 is well described by a thin-screen Kolmogorov model, at least along the time and frequency axes. Observed modulation indices for the DISS time scale and bandwidth and the pulsar flux density are greater than values predicted for a Kolmogorov spectrum of electron density fluctuations. Correlated variations over times long compared to the nominal refractive scintillation time are observed, suggesting that larger-scale density fluctuations are important. For these pulsars, the scintillation bandwidth as a function of frequency has a power-law index (~3.6) much less than expected for Kolmogorov turbulence (~4.4). Sloping fringes are commonly observed in the dynamic spectra, especially for PSR B0329+54. The detected range of fringe slopes are limited by our observing resolution. Our observations are sensitive to larger-scale fringes and hence smaller refractive angles, corresponding to the central part of the scattering disk.

In this paper we report unusual glitches in two young pulsars, PSR J1825-0935 (B1822-09) and PSR J1835-1106. For PSR J1825-0935, a slow glitch characterised by a temporary decrease in the slowdown rate occurred between 2000 December 31 to 2001 December 6. This event resulted in a permanent increase in frequency with fractional size $\Delta\nu/\nu\sim31.2(2)\times10^{-9}$, however little effect remained in slowdown rate. The glitch in PSR J1835-1106 occurred abruptly in November 2001 (MJD 52220\pm3) with $\Delta\nu/\nu\sim14.6(4)\times10^{-9}$ and little or no change in the slow-down rate. A significant change in $\ddot\nu$ apparently occurred at the glitch with $\ddot\nu$ having opposite sign for the pre- and post-glitch data.

A starquake mechanism for pulsar glitches is developed in the solid quark star model. It is found that the general glitch natures (i.e., the glitch amplitudes and the time intervals) could be reproduced if solid quark matter, with high baryon density but low temperature, has properties of shear modulus \mu = 10^{30~34} erg/cm^3 and critical stress \sigma_c = 10^{18~24} erg/cm^3. The post-glitch behavior may represent a kind of damped oscillations.

This paper summarizes the results of 13 years of timing observations of a unique binary pulsar, PSR B1259$-$63, which has a massive B2e star companion. The data span encompasses four complete orbits and includes the periastron passages in 1990, 1994, 1997 and 2000. Changes in dispersion measure occurring around the 1994, 1997 and 2000 periastrons are measured and accounted for in the timing analysis. There is good evidence for a small glitch in the pulsar period in 1997 August, not long after the 1997 periastron, and a significant frequency second derivative indicating timing noise. We find that spin-orbit coupling with secular changes in periastron longitude and projected semi-major axis ($x$) cannot account for the observed period variations over the whole data set. While fitting the data fairly well, changes in pulsar period parameters at each periastron seem ruled out both by X-ray observations and by the large apparent changes in pulsar frequency derivative. Essentially all of the systematic period variations are accounted for by a model consisting of the 1997 August glitch and step changes in $x$ at each periastron. These changes must be due to changes in the orbit inclination, but we can find no plausible mechanism to account for them. It is possible that timing noise may mask the actual changes in orbital parameters at each periastron, but the good fit to the data of the $x$ step-change model suggests that short-term timing noise is not significant.

A pulsar timing system has been operating in the 18-cm band at the Urumqi Astronomical Observatory 25-m telescope since mid-1999. Frequency resolution allowing dedispersion of the pulsar signals is provided by a 2X128X2.5 MHz filterbank/digitiser system. Observations of 74 pulsars over more than 12 months have resulted in updated pulsar periods and period derivatives, as well as improved positions. Comparison with previous measurements showed that the changes in period and period derivative tend to have the same sign and to be correlated in amplitude. A model based on unseen glitches gives a good explanation of the observed changes, suggesting that long-term fluctuations in period and period derivatives are dominated by glitches. In 2000 July, we detected a glitch of relative amplitude Delta_nu/nu~24X1e-9 in the Crab pulsar. The post-glitch decay appears similar to other large Crab glitches.

Timing observations of 40 mostly young pulsars using the ATNF Parkes radio telescope between 1990 January and 1998 December are reported. In total, 20 previously unreported glitches and ten other glitches were detected in 11 pulsars. These included 12 glitches in PSR J1341$- $6220, corresponding to a glitch rate of 1.5 glitches per year. We also detected the largest known glitch, in PSR J1614$-$5047, with $\Delta\nu_g/\nu \approx 6.5 \times 10^{-6}$ where $\nu = 1/P$ is the pulse frequency. Glitch parameters were determined both by extrapolating timing solutions to inter-glitch intervals and by phase-coherent timing fits across the glitch(es). Analysis of glitch parameters, both from this work and from previously published results, shows that most glitches have a fractional amplitude $\Delta\nu_g/\nu$ of between $10^{-8}$ and $10^{-6}$. There is no consistent relationship between glitch amplitude and the time since the previous glitch or the time to the following glitch, either for the ensemble or for individual pulsars. As previously recognised, the largest glitch activity is seen in pulsars with ages of order 10$^4$ years, but for about 30 per cent of such pulsars, no glitches were detected in the 8-year data span. There is some evidence for a new type of timing irregularity in which there is a significant increase in pulse frequency over a few days, accompanied by a decrease in the magnitude of the slowdown rate. Fits of an exponential recovery to post-glitch data show that for most older pulsars, only a small fraction of the glitch decays. In some younger pulsars, a large fraction of the glitch decays, but in others, there is very little decay.