Isolated pulsars are rotating neutron stars with accurately measured angular velocities $\Omega$, and their time derivatives that show unambiguously that the pulsars are slowing down. The commonly accepted view is that it arises through emission of magnetic dipole radiation (MDR) from a rotating magnetized body. The calculated energy loss by a rotating pulsar with a constant moment of inertia is assumed proportional to a model dependent power of $\Omega$. This relation leads to the power law $\dot{\Omega}$ = -K $\Omega^{\rm n}$ where $n$ is called the braking index. The MDR model predicts $n$ exactly equal to 3. Selected observations of isolated pulsars provide rather precise values of $n$, individually accurate to a few percent or better, in the range 1$ <$ n $ < $ 2.8, which is consistently less than the predictions of the MDR model. In spite of an extensive investigation of various modifications of the MDR model, no satisfactory explanation of observation has been found yet. The aim of this work is to determine the deviation of the value of $n$ from the canonical $n = 3$ for a star with a frequency dependent moment of inertia in the region of frequencies from zero to the Kepler velocity (onset of mass shedding by a rotating deformed star), in the macroscopic MDR model. For the first time, we use microscopic realistic equations of state (EoS) of the star to determine its behavior and structure. In addition, we examine the effects of the baryonic mass M$_{\rm B}$ of the star, and possible core superfluidity, on the value of the braking index within the MDR model.

Isolated pulsars are rotating neutron stars with accurately measured angular velocities $\Omega$, and their time derivatives which show unambiguously that the pulsars are slowing down. Although the exact mechanism of the spin-down is a question of debate in detail, the commonly accepted view is that it arises through emission of magnetic dipole radiation (MDR) from a rotating magnetized body. Other processes, including the emission of gravitational radiation, and of relativistic particles (pulsar wind), are also being considered. The calculated energy loss by a rotating pulsar with a constant moment of inertia is assumed proportional to a model dependent power of $\Omega$. This relation leads to the power law $\dot{\Omega}$ = -K $\Omega^{\rm n}$ where $n$ is called the braking index. The MDR model predicts $n$ exactly equal to 3. Selected observations of isolated pulsars provide rather precise values of $n$, individually accurate to a few percent or better, in the range 1$ <$ n $ < $ 2.8, which is consistently less than the predictions of the MDR model. In spite of an extensive investigation of various modifications of the MDR model, no satisfactory explanation of observation has been found as yet. We employ four physically realistic equations of state, and two computational codes, to model the dynamical effects of rotation on the braking index in the MDR model. In addition to this we simulate an effect on moment of inertia where we assume a certain amount of superfluid matter has manifest between the crust and core region of the star thus essentially eliminating momentum transfer between the two regions. We find that the effects of rotation on braking index are significant at high frequencies, but have little effect at frequencies consistent with the most accurately measured pulsars to-date.

The innermost stable cicular orbit (ISCO) of an accretion disc orbiting a neutron star (NS) is often assumed a unique prediction of general relativity. However, it has been argued that ISCO also appears around highly elliptic bodies described by Newtonian theory. In this sense, the behaviour of an ISCO around a rotating oblate neutron star is formed by the interplay between relativistic and Newtonian effects. Here we briefly explore the consequences of this interplay using a straightforward analytic approach as well as numerical models that involve modern NS equations of state. We examine the ratio K between the ISCO radius and the radius of the neutron star. We find that, with growing NS spin, the ratio K first decreases, but then starts to increase. This non-monotonic behaviour of K can give rise to a neutron star spin interval in which ISCO appears for two very different ranges of NS mass. This may strongly affect the distribution of neutron stars that have an ISCO (ISCO-NS). When (all) neutron stars are distributed around a high mass M0, the ISCO-NS spin distribution is roughly the same as the spin distribution corresponding to all neutron stars. In contrast, if M0 is low, the ISCO-NS distribution can only have a peak around a high value of spin. Finally, an intermediate value of M0 can imply an ISCO-NS distribution divided into two distinct groups of slow and fast rotators. Our findings have immediate astrophysical applications. They can be used for example to distinguish between different models of high-frequency quasiperiodic oscillations observed in low-mass NS X-ray binaries.