Highly magnetized, fast rotating neutron stars called pulsars are one of the most stable natural rotators in the universe. Their clock-like behavior has a widespread range of applications from detections of the stochastic background of gravitational waves   to testing theories of gravity  . However, the clocks are far from perfect as they slow down with time leading to increase in spin frequency v and its first time derivative according to the spin down law  given as
or its equivalent in terms of the spin period and its first time derivative given as
where the exponent n is the braking index and K is a positive constant based on standard spin-down model of a dipole rotator in a vacuum. The surface magnetic field can be explicitly expressed  as
where I and R are respectively, the moment of inertia and the radius; α is the angle between the magnetic and the spin axes of pulsar and c is the speed of light. On the other hand, direct measurements of pulsar braking index n can be accomplished through measurements of pulse spin period and its first and second time derivatives respectively and is given   as
where is period second time derivative. Conversely, the radiations emitted by pulsars are at the expense of another physical quantity, which can be inherent to the pulsar or external. Those pulsars whose emissions are at the expense of the rotational kinetic energy are known as rotation-powered pulsars . For these pulsars (radius of 10 km and moment of inertia ~ 1045 g∙cm−2), the estimate of the surface magnetic field is given  as
The detailed long-term timing observations of these pulsars have revealed that the otherwise smooth spindown of most radio pulsars is prone to a wide variety of disruptions called timing noise, which have proven extremely useful in the use of pulsars as physical tools. At present, the observed timing activities are broadly grouped into glitches and timing noise. Glitches are sudden tiny jumps in pulsar spin rate which are usually accompanied by increase in the magnitude of frequency derivative . Glitches are uncommon spectacular events characterized by recovery of the jumps in frequency and its first time derivatives on a large range of timescales  . Currently, it is generally agreed that glitches come from various complex dynamical changes within the interior of neutron star; hence, the study could offer valuable information about the internal structure and dynamics of neutron stars   . On the contrary, timing noise, is the random, sustained fluctuations in the pulsars’ spin parameters, that lasts from days to years and results in excess timing phase residuals . It is believed to occur often in rotation-powered pulsars in which the pulse arrival times wander continuously about the fitted ephemeris  . Even though the predominance of timing noise (as unmodeled structures in pulsar timing residuals in excess of what is expected from measurement error) within the pulsar samples seems not fully established   , our understanding of the phenomenon is still largely poor  . Meanwhile, timing noise has been attributed to some form of variations in the coupling between neutron star crust and its superfluid core   and/or torque fluctuations in the magnetosphere    .
Hitherto, a number of parameters were introduced to illustrate and measure the amount of timing noise in pulsars. Theoretically, a qualitative measurement of the level of timing noise activity in pulsars can be obtained from the surplus phase residuals remaining after accounting for its deterministic spin-down  .  in that effect, fitted the observed rotational phase of each pulsar with a function form
which is a Taylor series expansion around a reference time t0, where the terms , , and are the spin frequency, its first and second time derivatives respectively, while is the phase at a time t. The activity parameter is the logarithm of the ratio of the root-mean-square deviation over a 12 ms observation for a given pulsar to the same quantity in the crab pulsar , given as
where m is the order of the polynomial fit to the data, T is the time span of the fit in days, is the timing noise contribution to the gross root mean square phase residuals obtained from the quadrature difference between the mean squares of the observed phase residuals and the measurement errors  while is the timing noise of the crab pulsar (in ms).  introduced the stability parameter expressly used to characterize the intensity of the fluctuations observed in pulsars, which is based on the amplitude of the observed spin frequency second time derivative and can be written as
where the is frequency second time derivative over a period of 108. The cubic term dominates the variance for a timing noise in the group of average of any red process with a monotonically decaying spectrum. Pulsar clock stability parameter introduced by  is a dimensionless Allan variance-like parameter used to quantify the degree of rotational stability in rotation-powered pulsars and is defined as
where is the root-mean-square of , measured over observing time span T . Conventionally, T = 10 yrs. Meanwhile, a direct quantification of the contamination in from the standard timing activity by all forms of timing fluctuation is 
where and are, respectively, the variances of the timing residuals over T after a second and third-order polynomial fits. The developments of timing noise parameters with various spin-down properties have been used for detailed statistical study of timing noise in pulsar populations     . Specifically, some statistical analyses of timing noise phenomenon in rotation-powered pulsars involving spin-down properties have shown that timing noise is common among young pulsars      . However, no result of timing observations has revealed that timing irregularities are strongly dominated by the observed surface magnetic field of a sizeable number of pulsars over long timescales (>10 years). Thus, the main idea behind this paper is to explore the contributions of the pulsar’s surface magnetic field strength on timing noise irregularities and the intrinsic evolution of pulsars.
2. Data Analysis and Results
The data used in the current analysis is a large homogeneous sample of Jodrell Bank Observatory (JBO) radio pulsars with significant measurements of timing noise parameters estimated from radio timing data with time span length of over four decades. The activity parameter was computed using Equation (7) and taking parameters from  on the assumption that the method used effectively whitened the timing residuals of the pulsars in the sample. A sample of 366 objects with complete information on rotation-powered pulsars’ stability parameter ( ) and pulsar clock stability parameter were compiled from . Also, the timing noise statistic σ23 was calculated from the information in  using Equation (10). The corresponding surface magnetic field was extracted from the Australian Telescope National Facility (ATNF) Pulsar website available at the website http://www.atnf.csiro.au/research/pulsar/psrcat. We used both analytic and statistical methods for this study. The statistical method was used to estimate the measures of dispersion of our sample while the correlation analysis was used to check for any form of relationship between the parameters. The Pearson product moment enabled us to calculate the correlation coeﬃcient r, which provides the strength of correlation between the parameters. The distributions of the activity parameter and pulsar clock stability parameter are shown in Figure 1(a) and Figure 1(b). For the activity parameter, it ranges from −2.24 to 1.83 with mean value value of ~ −0.83 ± 0.02, and median value of ~ −0.46. The figure shows that about ~210 (~81%) of the pulsars in the sample have activity parameter less than zero. However, in the case of Pulsar clock stability parameter, the values are from −13.13 to −7.04 with mean value of ~ −10.81 ± 0.02 and median value of −10.09 with few ~112 (constituting 23% of pulsars) have clock stability parameter values less than the mean value. Figure 1(c) and Figure 1(d) are the distributions of stability parameter and timing noise statistic σ23. The stability parameter ranges from −2.50 to −0.92 with mean and median values of ~ −1.38 ± 0.03 and ~ −1.02 respectively in a multimodal configuration. For the timing noise statistic, it ranges from −1.08 to 3.50 with mean value of
Figure 1. Histograms showing distributions of (a) Activity parameter (b) pulsar clock stability parameter (c) stability parameter (d) Timing Noise statistic of our sample.
~1.27 ± 0.02 and median value of 1.01. The distribution appears to be bimodal, which could suggest two kinds of timing noise statistics, and skews to the right (with skewness (μ) in the range: 0.03 ≤ μ ≤ 0.02) with many pulsars (~78%) having values less than the mean value. In order to evaluate for any possible relationship between timing noise parameters and surface magnetic field, we made plots of timing noise parameters as a function of surface magnetic field.
In doing this, a logarithmic value of surface magnetic field was used to reduce the nearness in the value. Figure 2(a) & Figure 2(b) show the bivariate scatter plots of the activity parameter and pulsar clock stability parameter against the surface magnetic field for the pulsars in our sample. In the figures, the magnetic field of pulsar population in our sample is in the range 1 × 107 ≤ B ≤ 1.58 × 1014 Gauss. The figures clearly show an upward sloping trend with closely clustered scatter points depicting that pulsars with large values of both activity parameter and pulsar clock stability parameter have large values of surface magnetic field. A linear regression analysis of our data yields a positive correlation coefficient (r ~0.88) for the two variables being studied. The exact form of their relationship is
Figure 2. Scatter plots of (a) activity parameter (b) pulsar clock stability parameter (c) stability parameter (d) timing noise statistic against the logarithm of surface magnetic field.
However, Figure 2(c) & Figure 2(d) show the scatter plots of the stability parameter and timing noise statistic σ23 against the surface magnetic field. The plots, which show an upward sloping tendency indicates that pulsars with large surface magnetic field exhibit high timing noise irregularities. A regression analyses of our data show that the stability parameter and timing noise statistic σ23 have significant positive correlation with the surface magnetic field (r ~ > 0.80) in each case with the linear relationship of the form
implying that surface magnetic is a key promoter of timing noise irregularity in these pulsars.
The result of the statistical analysis of a large data base of rotation powered pulsars with significant measurements of timing noise parameters and their effects on the intrinsic surface magnetic field has been presented. The surface magnetic field is a derived and rotational parameter and equally an intrinsic property of radio pulsar. The radio pulsars’ spin period and its first time derivative is a key property of radio pulsars . The surface magnetic field is dependent on the pulsar spin period and its first time derivative, thus, the existence of a positive relationship between the pulsar spin periods and its first time derivative on the timing noise of pulsars is remarkable . Hence, the surface magnetic field could be observed to show positive relationship with the timing noise parameters. This implies that timing noise is widespread and prevalent among radio pulsars that have large surface magnetic field. We observed that pulsars with high values of timing noise parameters are associated with large surface magnetic field. Previous results show that young pulsars (104 years) appear to show more timing noise activity than their old, more stable counterparts      . Consequently, the existence of this dependence of timing noise activity on the surface magnetic field is prominent among young pulsars with large surface magnetic field. A significant feature in the plots of the timing noise parameters against the surface magnetic field with the values of the surface magnetic field in the of 1011 - 1012 Gauss is remarkable in that this range corresponds to the triple pulsar system PSR 1714-34 in the globular system. These high values of surface magnetic field with the timing noise parameters could be attributed to the acceleration of pulsars in the globular cluster due to their strong gravitational field. It is believed that timing noise is a direct consequence of the source of the intrinsic strong gravitation and fluctuations in the globular system. Moreover, the very strong (r ≥ 0.80) correlations which characterize the relationships between timing noise parameters and the surface magnetic field is a clear suggestion that pulsars with high surface magnetic fields on the average display more timing noise than their counterparts.
The observed surface magnetic field of the radio pulsars has been carefully studied for characterization with the pulsar timing noise. Our results showed that the strength of radio pulsar surface magnetic field have obvious significant effects on the timing noise parameters. There is existence of positive relationship between the parameters that characterize timing noise and the radio pulsar surface magnetic field. This is a clear indication that pulsars with large surface magnetic field are highly susceptible to timing noise activity.
 Shannon, R.M., Ravi, V., Lentati, L.T., Lasky, P.D., Hobbs, G., Kerr, M., et al. (2015) Gravitational Waves from Binary Supermassive Blackhole Missing in Pulsar Observations. Science, 349, 1522-1525.
 Arzoumanian, Z., Baker, P.T., Brazier, A., Burke-Spolaor, S., Chamberlin, S.J., Chatterjee, S., et al. (2018) The NANOGrav 11 Year Dataset: Pulsar Timing Constraints on the Stochastic Gravitational Wave Background. Astrophysical Journal, 859, 47.
 Kramer, M., Stairs, I.H., Manchester, R.N., McLaughlin, M.A., Lyne, A.G., Ferdman, R.D., et al. (2006) Tests of General Relativity from Timing the Double Pulsar. Science, 314, 97-102.
 Cameron, A.D., Champion, D.J., Kramer, M., Bailes, M., Barr, E.D., Bassa, C.G., et al. (2018) The High Timing Resolution Universe Pulsar Survey-XIII. PSR J1757-1854, the Most Accelerated Binary Pulsar. Monthly Notices of the Royal Astronomical Society: Letters, 475, L57-L61.
 Allen, M.P. and Horvath, J.E. (1997) Glitches, Torque Evolution and the Dynamics of Young Pulsars. Monthly Notices of the Royal Astronomical Society, 287, 615-621.
 Yu, M., Manchester, R.N., Hobbs, G., Johnston, S., Kaspi, V.M., Keith, M., et al. (2013) Detection of 107 Glitches in 36 Southern Pulsars. Monthly Notices of the Royal Astronomical Society, 429, 688-724.
 Eya, I.O., Urama, J.O. and Chukwude, A.E. (2017) Angulsar Momentum Transfer and Fractional Moment of Inertia in Pulsar Glitches. Astrophysical Journal, 840, 56.
 Hobbs, G., Lyne, A.G. and Kramar, M. (2010) An Analysis of the Timing Irregularities for 366 Pulsars. Monthly Notices of the Royal Astronomical Society, 402, 1027-1048.
 Eya, I.O., Urama, J.O. and Chukwude, A.E. (2019) On the Distributions of Pulsar Glitch Sizes and Inter-Glitch Time Intervals. Research in Astronomy and Astrophysics, 19, 89.
 Kramer, M., Lyne, A.G., O’Brien, J.T., Jordan, C.A. and Lorimer, D.R. (2006) Periodically Active Pulsar Giving Insight into Magnetospheric Physics. Science, 312, 549-551.
 Urama, J.O., Link, B. and Weisberg, J.M. (2006) A Strong-Correlation in Radio Pulsars with Implications for Torque Variations. Monthly Notices of the Royal Astronomical Society: Letters, 370, L76-L79.
 Chukwude, A.E. and Odo, F.C. (2016) Assessing the Effects of Timing Irregularities on Radio PulsarsAnomalous Braking Indices. Research in Astronomy and Astrophysics, 16, 150.
 Iyida, E.U., Eze, C.I., Eya, I.O. and Chukwude, A.E. (2018) On the Effects of Spin Properties on Timing Noise Parameters of Rotation-Powered Pulsars. International Journal of Astrophysics and Space Science, 6, 101-107.
 D’Alessandro, F., McCulloch, P.M., King, E.A., Hamilton, P.A. and McConnell, D. (1993) Timing Observations of Southern Pulsars—1987 to 1991. Monthly Notices of the Royal Astronomical Society, 261, 883.