High resolution, pulse to pulse observation of the 2016 Vela glitch and its relaxation provided an opportunity to probe the neutron star internal structure and dynamics with unprecedented detail. We use the observations of this glitch to infer superfluid characteristics in the framework of the vortex creep model. The glitch rise time constraint of 12.6 seconds put stringent limits on the angular momentum exchange between the crustal superfluid and the observed crust. Together with the observed excess acceleration in the rotation rate as compared to the post-glitch equilibrium value this discriminates crustal superfluid-crust lattice and core superfluid-crustal normal matter coupling time-scales. An evident decrease in the crustal rotation rate immediately before the glitch is consistent with the formation of a new vortex trap zone that initiates the large scale vortex unpinning avalanche. Formation of vortex trap by a crust breaking quake induces short-lived magnetospheric changes. The long term post-glitch spin-down rate evolution reveals the moments of inertia and recoupling time-scales of the superfluid layers participating in the glitch and leads to an estimation of the time to the next glitch which agrees with the time interval between the 2016 and 2019 glitches. Our results are consistent with theoretical estimates of effective neutron and proton masses in the superfluid. We also constrain the vortex line-flux tube pinning energy per intersection as 2 MeV.

Nonequilibrium conditions imposed by neutrino cooling through the liquid-solid transition lead to disorder in the solid crust of neutron stars. Disorder reduces the superfluid fraction, $\rho_s/\rho$, at densities above that of neutron drip, $\rho_d \approx 4\times 10^{11}\,g/cm^3$. For an amorphous solid crust the suppression of $\rho_s$ is small, except in the highest density regions of the crust. In contrast to the strong reduction in neutron conduction predicted for coherent Bragg scattering in a crystalline crust, the disordered solid crust supports sufficient neutron superfluid density to account for pulsar glitches.

The Crab pulsar displayed its largest glitch on 2017 November. An extended initial spin-up phase of this largest glitch was resolved, for the first time with high cadence of observations both in radio and X-rays on a time-scale of 2 days. A combination of crustquake and vortex unpinning models is invoked to account for the extended spin-up, magnitude and post-glitch relaxation characteristics of this glitch. We evaluate the extended spin-up followed by the familiar spin-down as due to the creep response to the initial induced inward motion of some vortex lines pinned to broken crustal plates moving inward towards the rotation axis, together with the common and familiar post-glitch creep response to the sudden outward motion of vortices unpinned at the glitch. Our analysis confirms that the number of unpinned vortices participating in glitches are similar in all Crab glitches, and within an order of magnitude in all glitches from all pulsars. This typical number of unpinned vortices is related to the broken plate size in quakes as triggers for vortex unpinning avalanches. The physical determinant of this universal broken plate size is in turn the critical strain angle in the neutron star crust. Occurrence of this largest Crab glitch after a relatively long inactive period is consistent with accumulation of the pinned vorticity to be tapped.

We study the evolution of newborn neutron stars in high-mass X-ray binaries interacting with a wind-fed super-Eddington disk. The inner disk is regularized to a radiation-dominated quasi-spherical configuration for which we calculate the inner radius of the disk, the total luminosity of the system and the torque acting on the neutron star accordingly, following the evolution of the system through the ejector and early propeller stages. We find that the systems with $B \gtrsim 10^{13}$ G pass through a short ($\sim 20\,{\rm yr}$) ejector stage appearing as supernova impostors followed by a propeller stage lasting $\sim 10^3\,{\rm yr}$. In the super-critical propeller stage the system is still bright ($L\sim 10^{40}\,{\rm erg\, s^{-1}}$) due to the spindown power and therefore appears as an ultra-luminous X-ray source (ULX). The system evolves into pulsating ULX (PULX) when the neutron star spins down to a period ($P\sim 1$ s) allowing for accretion onto its surface to commence. Systems with lower magnetic fields, $B \sim 10^{11}$ G, pass through a long ($10^5\,{\rm yr}$) super-critical propeller stage with luminosities similar to those of the ultra-luminous super-soft sources (ULS), $L \lesssim 10^{40}\,{\rm erg\, s^{-1}}$. The equilibrium periods of these systems in the accretion stage is about $10\,{\rm ms}$, which is much smaller than the typical period range of PULX observed to date. Such systems could have a larger population, but their pulsations would be elusive due to the smaller size of the magnetosphere. Our results suggest that the ULS and some nonpulsating ULX are rapidly spinning and highly magnetized young neutron stars at the super-critical propeller stage.

The two component model of neutron star dynamics describing the behaviour of the observed crust coupled to the superfluid interior has so far been applied to radio pulsars for which the external torques are constant on dynamical timescales. We recently solved this problem under arbitrary time dependent external torques. Our solutions pertain to internal torques that are linear in the rotation rates, as well as to the extremely non-linear internal torques of the vortex creep model. Two-component models with linear or nonlinear internal torques can now be applied to magnetars and to neutron stars in binary systems, with strong variability and timing noise. Time dependent external torques can be obtained from the observed spin-down (or spin-up) time series, $\dot{\Omega}(t)$.

We discuss the minimum glitch size of Crab observed by Espinoza et al. (2014). Modelling the crustquake as a trigger mechanism we estimate the size of the broken plate. The plate size obtained, D $\sim 100 m$ is comparable to plate size estimates for PSR J1119--6127. The plate size naturally leads to an estimate of the number of unpinned vortices involved in the glitch, $N \sim 10^{13} $. This number is of the same order in all Crab and Vela pulsar glitches. The mimimum glitch relates the constancy of all these numbers among different pulsars to the basic plate size involved in crust breaking. This in turn depends on the critical strain angle $\theta_{cr}$ of the Coulomb lattice in the neutron star crust. The minimum glitch size implies $\theta_{cr} \sim 10^{-1}$ in agreement with theoretical and computational estimates.

G. Srinivasan et al. (1990) proposed a simple and elegant explanation for the reduction of the neutron star magnetic dipole moment during binary evolution leading to low mass X-ray binaries and eventually to millisecond pulsars: Quantized vortex lines in the neutron star core superfluid will pin against the quantized flux lines of the proton superconductor. As the neutron star spins down in the wind accretion phase of binary evolution, outward motion of vortex lines will reduce the dipole magnetic moment in proportion to the rotation rate. The presence of a toroidal array of flux lines makes this mechanism inevitable and independent of the angle between the rotation and magnetic axes. The incompressibility of the flux-line array (Abrikosov lattice) determines the epoch when the mechanism will be effective throughout the neutron star. Flux vortex pinning will not be effective during the initial young radio pulsar phase. It will, however, be effective and reduce the dipole moment in proportion with the rotation rate during the epoch of spindown by wind accretion as proposed by Srinivasan et al. The mechanism operates also in the presence of vortex creep.

The two component model describes neutron star dynamics incorporating the response of the superfluid interior. Conventional solutions and applications involve constant external torques, as appropriate for radio pulsars on dynamical timescales. We present the general solution of two component dynamics under arbitrary time dependent external torques, with internal torques that are linear in the rotation rates, or with the extremely non-linear internal torques due to vortex creep. The two-component model incorporating the response of linear or nonlinear internal torques can now be applied not only to radio pulsars but also to magnetars and to neutron stars in binary systems, with strong observed variability and noise in the spin-down or spin-up rates. Our results allow the extraction of the time dependent external torques from the observed spin-down (or spin-up) time series, $\dot{\Omega}(t)$. Applications are discussed.

The inter-glitch timing of the Vela pulsar is characterized by a constant second derivative of the rotation rate. This takes over after the post-glitch exponential relaxation, and is completed at about the time of the next glitch. The vortex creep model explains the second derivatives in terms of non-linear response to the glitch. We present inter-glitch timing fits to the present sample covering 16 large glitches, taking into account the possibility that in some glitches part of the step in spin-down rate may involve a "persistent shift", as observed in the Crab pulsar. Modifying the expression for the time between glitches with this hypothesis leads to better agreement with the observed inter-glitch time intervals. We extrapolate the inter-glitch model fits to obtain spin-down rates just prior to each glitch, and use these to calculate the braking index n = 2.81 +/- 0.12. The next glitch should occur around Dec. 22, 2017 +/- 197 days if no persistent shift is involved, but could occur as early as July 27, 2016 +/- 152 days if the 2013 glitch gave rise to a typical Vela persistent shift. Note added: Literally while we were submitting the first version of this paper, on Dec. 12, 2016, we saw ATel # 9847 announcing a Vela pulsar glitch which has arrived 138 days after our prediction with a persistent shift, within the 1 sigma uncertainty of 152 days.

The optical excess in the spectra of dim isolated neutron stars (XDINs) is a significant fraction of their rotational energy loss-rate. This is strikingly different from the situation in isolated radio pulsars. We investigate this problem in the framework of the fallback disc model. The optical spectra can be powered by magnetic stresses on the innermost disc matter, as the energy dissipated is emitted as blackbody radiation mainly from the inner rim of the disc. In the fallback disc model, XDINs are the sources evolving in the propeller phase with similar torque mechanisms. In this this model, the ratio of the total magnetic work that heats up the inner disc matter is expected to be similar for different XDINs. Optical luminosities that are calculated consistently with the the optical spectra and the theoretical constraints on the inner disc radii give very similar ratios of the optical luminosity to the rotational energy loss rate for all these sources. These ratios indicate that a significant fraction of the magnetic torque heats up the disc matter while the remaining fraction expels disc matter from the system. For XDINs, the contribution of heating by X-ray irradiation to the optical luminosity is negligible in comparison with the magnetic heating. The correlation we expect between the optical luminosities and the rotational energy loss-rates of XDINs can be a property of the systems with low X-ray luminosities, in particular those in the propeller phase.

Two energetic hard X-ray bursts have recently triggered the Fermi and Swift space observatories from the rotation powered pulsar, PSR J1119-6127. We have performed in depth spectral and temporal analyses of these two events. Our extensive searches in both observatory data for lower luminosity bursts uncovered 10 additional events from the source. We report here on the timing and energetics of the 12 bursts from PSR J1119-6127 during its burst active phase of 2016 July 26 and 28. We also found a spectral softer X-ray flux enhancement in a post burst episode, which shows evidence of cooling. We discuss here the implications of these results on the nature of this unusual high-field radio pulsar, which firmly place it within the typical magnetar population.

Treatment of the vortex motion in the superfluids of the inner crust and the outer core of neutron stars is a key ingredient in modeling a number of pulsar phenomena, including glitches and magnetic field evolution. After recalculating the microscopic vortex velocity in the inner crust, we evaluate the velocity for the vortices in the outer core for the first time. The vortex motion between pinning sites is found to be substantially faster in the inner crust than in the outer core, $v_0^{\rm crust} \sim 10^{7}\mbox{\cms} \gg v_0^{\rm core} \sim 1\mbox{\cms}$. One immediate result is that vortex creep is always in the nonlinear regime in the outer core in contrast to the inner crust, where both nonlinear and linear regimes of vortex creep are possible. Other implications for pulsar glitches and magnetic field evolution are also presented.

The correlation between the frequency and the absolute value of the frequency derivative of the kilohertz QPOs observed for the first time from 4U 1636-53 is a simple consequence and indicator of the existence of a non-Keplerian rotation rate in the accretion disk boundary layer. This paper interprets the observed correlation, showing that the observations provide strong evidence in support of the fundamental assumption of disk accretion models around slow rotators, that the boundary layer matches the Keplerian disk to the neutron star magnetosphere.

Glitches are sudden changes in rotation frequency and spin-down rate, observed from pulsars of all ages. Standard glitches are characterized by a positive step in angular velocity ($\Delta\Omega$ $ > $ $0$) and a negative step in the spin-down rate ($\Delta \dot \Omega$ $ < $ $0$) of the pulsar. There are no glitch-associated changes in the electromagnetic signature of rotation-powered pulsars in all cases so far. For the first time, in the last glitch of PSR J1119-6127, there is clear evidence for changing emission properties coincident with the glitch. This glitch is also unusual in its signature. Further, the absolute value of the spin-down rate actually decreases in the long term. This is in contrast to usual glitch behaviour. In this paper we extend the vortex creep model in order to take into account these peculiarities. We propose that a starquake with crustal plate movement towards the rotational poles of the star induces inward vortex motion which causes the unusual glitch signature. The component of the magnetic field perpendicular to the rotation axis will decrease, giving rise to a permanent change in the pulsar external torque.

The dim isolated neutron stars (XDINs) have periods in the same range as the anomalous X-ray pulsars (AXPs) and the soft gamma-ray repeaters (SGRs). We apply the fallback disk model, which explains the period clustering and other properties of AXP/SGRs, to the six XDINs with measured periods and period derivatives. Present properties of XDINs are obtained in evolutionary scenarios with surface dipole magnetic fields $B_0 \sim 10^{12}$ G. The XDINs have gone through an accretion epoch with rapid spin-down earlier, and have emerged in their current state, with the X-ray luminosity provided by neutron star cooling and no longer by accretion. Our results indicate that the known XDINs are not likely to be active radio pulsars, as the low $B_0$, together with their long periods place these sources clearly below the "death valley".

A region of toroidally oriented quantized flux lines must exist in the proton superconductor in the core of the neutron star. This region will be a site of vortex pinning and creep. Entrainment of the neutron superfluid with the crustal lattice leads to a requirement of superfluid moment of inertia associated with vortex creep in excess of the available crustal moment of inertia. This will effect constraints on the equation of state. The toroidal flux region provides the moment of inertia necessary to complement the crust superfluid with postglitch relaxation behavior fitting the observations.

We investigate the X-ray enhancement and the long-term evolution of the recently discovered, second "low-B magnetar" Swift J1822.3-1606 in the frame of the fallback disk model. During a soft gamma burst episode, the inner disk matter is pushed back to larger radii forming a density gradient at the inner disk. Subsequent relaxation of the inner disk could account for the observed X-ray enhancement light curve of Swift J1822.3-1606. We obtain model fits to the X-ray data with basic disk parameters similar to those employed to explain the X-ray outburst light curves of other AXPs and SGRs. The long period (8.4 s) of the neutron star can be reached by the effect of the disk torques in the long-term accretion phase (1-3 x 10^5 yrs). The currently ongoing X-ray enhancement could be due to a transient accretion epoch, or the source could still be in the accretion phase in quiescence. Considering these different possibilities, we determine the model curves that could represent the long-term rotational and the X-ray luminosity evolution of Swift J1822.3-1606, which constrain the strength of the magnetic dipole field to the range of 1-2 x 10^12 G on the surface of the neutron star.

Recent measurements showed that the period derivative of the 'high-B' radio pulsar PSR J1734-3333 is increasing with time. For neutron stars evolving with fallback disks, this rotational behavior is expected in certain phases of the long-term evolution. Using the same model as employed earlier to explain the evolution of anomalous X-ray pulsars and soft gamma-ray repeaters, we show that the period, the first and second period derivatives and the X-ray luminosity of this source can simultaneously acquire the observed values for a neutron star evolving with a fallback disk. We find that the required strength of the dipole field that can produce the source properties is in the range of 10^{12} - 10^{13} G on the pole of the neutron star. When the model source reaches the current state properties of PSR J1734-3333, accretion onto the star has not started yet, allowing the source to operate as a regular radio pulsar. Our results imply that PSR J1734-3333 is at an age of ~ 3 x 10^4 - 2 x 10^5 years. Such sources will have properties like the X-ray dim isolated neutron stars or transient AXPs at a later epoch of weak accretion from the diminished fallback disk.

The presence of matter with angular momentum, in the form of a fallback disk around a young isolated neutron star will determine its evolution. This leads to an understanding of many properties of different classes of young neutron stars, in particular a natural explanation for the period clustering of AXPs, SGRs and XDINs. The spindown or spinup properties of a neutron star are determined by the dipole component of the magnetic field. The natural possibility that magnetars and other neutron stars may have different strengths of the dipole and higher multipole components of the magnetic field is now actually required by observations on the spindown rates of some magnetars. This talk gives a broad overview and some applications of the fallback disk model to particular neutron stars. Salient points are: (i) A fallback disk has already been observed around the AXP 4U 0142+61 some years ago. (ii) The low observed spindown rate of the SGR 0418+5729 provides direct evidence that the dipole component of the field is in the $10^{12} G$ range. All properties of the SGR 0418+5729 at its present age can be explained by spindown under torques from a fallback disk. (iii) The anomalous braking index of PSR J1734-3333 can also be explained by the fallback disk model which gives the luminosity, period, period derivative and the period second derivative at the present age. (iv) These and all applications to a variety of other sources employ the same disk physics and evolution, differing only in the initial conditions of the disk.

When an accreting star is close to rotational equilibrium between the dipole component of the stellar magnetic field and the accretion disk, the star's rotation rate is roughly of the order of the Keplerean rotation rate at the inner boundary of the disk, estimated as the conventional Alfven radius. A range of frequencies higher than this equilibrium rotation frequency can naturally arise if the accretion flow is channeled by higher multipoles of the star's magnetic field. The higher multipole components of the magnetic field will balance the material stresses of the accretion flow at radii closer to the star. The Kepler frequencies associated with these generalized Alfven radii increase with the order of the multipole. Other frequency bands, like the epicyclic frequencies associated with the accretion flow, may in turn be higher than the Kepler frequencies. We present expressions for the spectrum of higher frequencies arising due to these effects. Kilohertz quasi-periodic oscillation frequencies that are much higher than the rotation rate of the neutron star, as observed from the recently discovered 11 Hz (P = 90 ms) X-ray pulsar IGR J17480-2446 in the globular cluster Terzan 5, may be due to modulation of the accretion rate by the excitation of these modes in the accretion flow. The very high QPO frequencies observed from the soft gamma repeaters SGR 1806-20 (P = 5.2 s) and SGR 1900+14 (P = 7.5 s) may also correspond to these characteristic frequencies if SGRs accrete from fallback disks around them.

The low mass X-ray binary (LMXB) IGR J17480-2446 in the globular cluster Terzan 5 harbors an 11 Hz accreting pulsar. This is the first object discovered in a globular cluster with a pulsar spinning at such low rate. The accreting pulsar is anomalous because its characteristics are very different from the other five known slow accreting pulsars in galactic LMXBs. Many features of the 11 Hz pulsar are instead very similar to those of accreting millisecond pulsars, spinning at frequencies >100 Hz. Understanding this anomaly is valuable because IGR J17480-2446 can be the only accreting pulsar discovered so far which is in the process of becoming an accreting millisecond pulsar. We first verify that the neutron star (NS) in IGR J17480-2446 is indeed spinning up by carefully analyzing X-ray data with coherent timing techniques that account for the presence of timing noise. We then study the present Roche lobe overflow epoch and the two previous spin-down epochs dominated by magneto dipole radiation and stellar wind accretion. We find that IGR J17480-2446 is very likely a mildly recycled pulsar and suggest that it has started a spin-up phase in an exceptionally recent time, that has lasted less than a few 10^7 yr. We also find that the total age of the binary is surprisingly low (<10^8 yr) when considering typical parameters for the newborn NS and propose different scenarios to explain this anomaly.

The period derivative bound for SGR 0418+5729 (Rea et al. 2010) establishes the magnetic dipole moment to be distinctly lower than the magnetar range, placing the source beyond the regime of isolated pulsar activity in the P - dP/dt diagram and giving a characteristic age > 2 \times 10^{7} years, much older than the 10^5 year age range of SGRs and AXPs. So the spindown must be produced by a mechanism other than dipole radiation in vacuum. A fallback disk will spin down a neutron star with surface dipole magnetic field in the 10^{12} G range and initial rotation period P_0 ~ 100 ms to the 9.1 s period of SGR 0418+5729 in a few 10^4 to ~10^5 years. The current upper limit to the period derivative gives a lower limit of \sim 10^5 years to the age that is not sensitive to the neutron star's initial conditions. The total magnetic field on the surface of SGR 0418+5729 could be significantly larger than its 10^{12} G dipole component.

We show that the period clustering of anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs), their X-ray luminosities, ages and statistics can be explained with fallback disks with large initial specific angular momentum. The disk evolution models are developed by comparison to self-similar analytical models. The initial disk mass and angular momentum set the viscous timescale. An efficient torque, with (1 - w^2) dependence on the fastness parameter w leads to period clustering in the observed AXP-SGR period range under a wide range of initial conditions. The timescale t_0 for the early evolution of the fallback disk, and the final stages of fallback disk evolution, when the disk becomes passive, are the crucial determinants of the evolution. The disk becomes passive at temperatures around 100 K, which provides a natural cutoff for the X-ray luminosity and defines the end of evolution in the observable AXP and SGR phase. This low value for the minimum temperature for active disk turbulence indicates that the fallback disks are active up to a large radius greater than ~10^{12} cm. We find that transient AXPs and SGRs are likely to be older than their persistent cousins. A fallback disk with mass transfer rates corresponding to the low quiescent X-ray luminosities of the transient sources in early evolutionary phases would have a relatively lower initial mass, such that the mass-flow rate in the disk is not sufficient for the inner disk to penetrate into the light cylinder of the young neutron star, making mass accretion onto the neutron star impossible. The transient AXP phase therefore must start later. The model results imply that the transient AXP/SGRs, although older, are likely to be similar in number to persistent sources (abridged).

We investigate dynamical coupling timescales of a neutron star's superfluid core, taking into account the interactions of quantized neutron vortices with quantized flux lines of the proton superconductor in addition to the previously considered scattering of the charged components against the spontaneous magnetization of the neutron vortex line. We compare the cases where vortex motion is constrained in different ways by the array of magnetic flux tubes associated with superconducting protons. This includes absolute pinning to and creep across a uniform array of flux lines. The effect of a toroidal arrangement of flux lines is also considered. The inclusion of a uniform array of flux tubes in the neutron star core significantly decreases the timescale of coupling between the neutron and proton fluid constituents in all cases. For the toroidal component, creep response similar to that of the inner crust superfluid is possible.

Basic ideas about the torques on the neutron star and the existence of an equilibrium rotation period followed from the recognition that most X-ray binaries contain accretion powered neutron stars. The evolution of binaries through a phase of accretion onto the neutron star, eventually leading to a post-accretion radio pulsar phase, was initially discussed as a way to understand the scarcity of binaries among the radio pulsars and the relatively short rotation periods of the first discovered binary radio pulsars in terms of magnetic fields that would be smaller than the familiar $10^{12} G$ range. The discovery of the millisecond pulsars made us realize that the fields can be much lower in a new class of radio pulsars that have been spun up by accretion in LMXBs. The predicted spin-down rates of the millisecond pulsar was soon confirmed. The observers' search for millisecond X-ray periods was on, leading first to the discovery of QPOs, and eventually to the discovery of the X-ray millisecond pulsars. The theorists' quest for explanations of why X-ray millisecond pulsations are not observed from LMXBs also started right away.

The observational characteristics of quasi-periodic oscillations (QPOs) from accreting neutron stars strongly indicate the oscillatory modes in the innermost regions of accretion disks as a likely source of the QPOs. The inner regions of accretion disks around neutron stars can harbor very high frequency modes related to the radial epicyclic frequency $\kappa $. The degeneracy of $\kappa $ with the orbital frequency $\Omega $ is removed in a non-Keplerian boundary or transition zone near the magnetopause between the disk and the compact object. We show, by analyzing the global hydrodynamic modes of long wavelength in the boundary layers of viscous accretion disks, that the fastest growing mode frequencies are associated with frequency bands around $\kappa $ and $\kappa \pm \Omega $. The maximum growth rates are achieved near the radius where the orbital frequency $\Omega $ is maximum. The global hydrodynamic parameters such as the surface density profile and the radial drift velocity determine which modes of free oscillations will grow at a given particular radius in the boundary layer. In accordance with the peak separation between kHz QPOs observed in neutron-star sources, the difference frequency between two consecutive bands of the fastest growing modes is always related to the spin frequency of the neutron star. This is a natural outcome of the boundary condition imposed by the rotating magnetosphere on the boundary region of the inner disk.

The recent detection of the anomalous X-ray pulsar (AXP) 4U 0142+61 in the mid-infrared with the {\it Spitzer} observatory by Z.Wang and coworkers constitutes the first instance of a disk around an AXP. We show, by analyzing earlier optical and near-IR data together with the recent data, that the overall broadband data set can be reproduced by a single model of an irradiated and viscously heated disk.

The spin periods of the neutron stars in most Low Mass X-ray Binary (LMXB) systems still remain undetected. One of the models to explain the absence of coherent pulsations has been the suppression of the beamed signal by Compton scattering of X-ray photons by electrons in a surrounding corona. We point out that simultaneously with wiping out the pulsation signal, such a corona will upscatter (pulsating or not) X-ray emission originating at and/or near the surface of the neutron star leading to appearance of a hard tail of Comptonized radiation in the source spectrum. We analyze the hard X-ray spectra of a selected set of LMXBs and demonstrate that the optical depth of the corona is not likely to be large enough to cause the pulsations to disappear.

The anomalous X-ray pulsar 4U 0142+61 was recently detected in the mid infrared bands with the SPITZER Observatory (Wang, Chakrabarty & Kaplan 2006). This observation is the first instance for a disk around an AXP. From a reanalysis of optical and infrared data, we show that the observations indicate that the disk is likely to be an active disk rather than a passive dust disk beyond the light cylinder, as proposed in the discovering paper. Furthermore, we show that the irradiated accretion disk model can also account for all the optical and infrared observations of the anomalous X-ray pulsars in the persistent state.

The cooling and reheating histories of dim isolated neutron stars(DINs) are discussed. Energy dissipation due to dipole spindown with ordinary and magnetar fields, and due to torques from a fallback disk are considered as alternative sources of reheating which would set the temperature of the neutron star after the initial cooling era. Cooling or thermal ages are related to the numbers and formation rates of the DINs and therefore to their relations with other isolated neutron star populations. Interaction with a fallback disk, higher multipole fields and activity of the neutron star are briefly discussed.

Almost all pulsars with anomalous positive $\ddot \Omega $ measurements (corresponding to anomalous braking indices in the range 5$<n<$100), including all the pulsars with observed large glitches ($\Delta\Omega/\Omega$ $>$ 10$^{-7}$) as well as post glitch or interglitch $\ddot \Omega$ measurements obey the scaling between $\ddot \Omega$ and glitch parameters originally noted in the Vela pulsar. Negative second derivative values can be understood in terms of glitches that were missed or remained unresolved. We discuss the glitch rates and a priori probabilities of positive and negative braking indices according to the model developed for the Vela pulsar. This behavior supports the universal occurrence of a nonlinear dynamical coupling between the neutron star crust and an interior superfluid component. The implied lower limit to dynamical energy dissipation in a neutron star with spindown rate $\dot \Omega$ is $\dot E_{diss}> 1.7 \times 10 ^{-6} \dot E_{rot}$. Thermal luminosities and surface temperatures due to dynamical energy dissipation are estimated for old neutron stars which are spinning down as rotating magnetic dipoles beyond the pulsar death line.

The long term (~1.5 years) X-ray enhancement and the accompanying infrared enhancement light curves of the anomalous X-ray pulsar 1E 2259+58 following the major bursting epoch can be accounted for by the relaxation of a fall back disk that has been pushed back by a gamma-ray flare. The required burst energy estimated from the results of our model fits is low enough for such a burst to have remained below the detection limits. We find that an irradiated disk model with a low irradiation efficiency is in good agreement with both X-ray and infrared data. Non-irradiated disk models also give a good fit to the X-ray light curve, but are not consistent with the infrared data for the first week of the enhancement.

In the inner regions of accretion disks around compact objects, the orbital frequency of the gas deviates from the local Keplerian value. For long-wavelength modes in this region, the radial epicyclic frequency kappa is higher than the azimuthal frequency Omega. This has significant implications for models of the twin kHz QPOs observed in many neutron-star sources that traditionally identify the frequencies of the two kHz QPOs with dynamical frequencies in the accretion disk. The recognition that the highest frequency in the transition or boundary region of the disk is actually the epicyclic frequency also modifies significantly the constraints imposed by the observation of high-frequency QPOs on the mass and radius of the compact objects.

Basic rotational and magnetic properties of neutron superfluids and proton superconductors in neutron stars are reviewed. The modes of precession of the neutron superfluid are discussed in detail. We emphasize that at finite temperature, pinning of superfluid vortices does not offer any constraint on the precession. Any pinning energies can be surmounted by thermal activation and there exists a dynamical steady state in which the superfluid follows the precession of the crust at a small lag angle between the crust and superfluid rotation velocity vectors. At this small lag the system is far from the critical conditions for unpinning, even if the observed precession of the crust may entail a large angle between the figure axis and the crust's rotation velocity vector. We conclude that if long period modulations of pulse arrival times and pulse shapes observed in a pulsar like the PSR B1828-11 are due to the precession of the neutron star, this does not have any binding implications about the existence of pinning by flux lines or the existence of Type II superconductivity in the neutron star.

We investigate the rotational dynamics of disk accretion around a strongly magnetized neutron star with an aligned dipole field. The magnetospheric field is assumed to thread the disk plasma both inside and outside the corotation radius. As a result of disk-star interaction, the magnetic torque on the disk affects the structure of accretion flow to yield the observed spin-up or spin-down rates for a source of given fastness, magnetic field strength, and mass accretion rate. Within the model we obtain a prescription for the dynamical viscosity of such magnetically modified solutions for a Keplerian disk. We then use this prescription to find a model solution for the rotation rate profile throughout the entire disk including the non-Keplerian inner disk. We find that the non-Keplerian angular velocity transition region is not necessarily narrow for a source of given spin state. The boundary layer approximation, as in the standard magnetically threaded disk model, holds only in the case of dynamical viscosity decreasing all the way to the innermost edge of the disk. These results are applied to several observed disk-fed X-ray pulsars that have exhibited quasi-periodic oscillations (QPOs). The QPO frequencies provide a constraint on the fastness parameter and enable one to determine uniquely the width of the angular velocity transition zone for each source within model assumptions. We discuss the implications of these results on the value of the critical fastness parameter for a magnetized star in spin equilibrium. Applications of our model are also made with relevant parameters from recent numerical simulations of quasi-stationary disk-magnetized star interactions.

The radiation pressure of a radio pulsar does not necessarily disrupt a surrounding disk. The position of the inner radius of a thin disk around a neutron star can be estimated by comparing the electromagnetic energy density generated by the neutron star with the kinetic energy density of the disk. Inside the light cylinder, the near zone electromagnetic field is essentially the dipole magnetic field, and the inner radius is the conventional Alfven radius. Far outside the light cylinder, in the radiation zone, $E=B$ and the electromagnetic energy density is $<S >/c \propto 1/r^2$ where $S$ is the Poynting vector. Shvartsman (1970) argued that a stable equilibrium can not be found in the radiative zone because the electromagnetic energy density dominates over the kinetic energy density, with the relative strength of the electromagnetic stresses increasing with radius. In order to check whether this is true also near the light cylinder, we employ global electromagnetic field solutions for rotating oblique magnetic dipoles (Deutsch 1955). Near the light cylinder the electromagnetic energy density increases steeply enough with decreasing $r$ to balance the kinetic energy density at a stable equilibrium. The transition from the near zone to the radiation zone is broad. The radiation pressure of the pulsar can not disrupt the disk for values of the inner radius up to about twice the light cylinder radius if the rotation axis and the magnetic axis are orthogonal. This allowed range beyond the light cylinder extends much further for small inclination angles. We discuss implications of this result for accretion driven millisecond pulsars and young neutron stars with fallback disks.

We assume that all the energy loss of the putative pulsar in SN 1987A would contribute to the luminosity of the remnant, which acts like a bolometer. The bolometric luminosity of SN 1987A provides an upper bound on the pulsar's rate of energy loss. An isolated pulsar spinning down by magnetic dipole radiation alone, with initial rotation periods of 10-30 ms, as extrapolated for galactic young pulsars, can have a luminosity below the bolometric bound if either the magnetic field is weak, B $\sim$ 10$^{9}-10^{10}$ G, or if it is so strong that the pulsar luminosity decays rapidly, B $\sim$ 10$^{16}$ G.

We discuss the options of the fall-back disk model of Anomalous X-Ray Pulsars (AXPs). We argue that the power-law index of the mass inflow rate during the propeller stage can be lower than those employed in earlier models. We take into account the effect of the super-critical mass inflow at the earliest stages on the inner radius of the disk and argue that the system starts as a propeller. Our results show that, assuming a fraction of the mass inflow is accreted onto the neutron star, the fall-back disk scenario can produce AXPs for acceptable parameters.

We investigate whether young neutron stars with fall-back disks can produce Anomalous X-Ray Pulsars (AXPs) within timescales indicated by the ages of associated supernova remnants. The system passes through a propeller stage before emerging as an AXP or a radio pulsar. The evolution of the disk is described by a diffusion equation which has self-similar solutions with either angular momentum or total mass of the disk conserved. We associate these two types of solutions with accretor and propeller regimes, respectively. Our numerical calculations of thin disk models with changing inner radius take into account the super-critical accretion at the early stages, and electron scattering and bound-free opacities with rich metal content. Our results show that, assuming a fraction of the mass inflow is accreted onto the neutron star, the fall-back disk scenario can produce AXPs for acceptable parameters.

We present the first results of the observations of the extremely bright optical afterglow of gamma-ray burst (GRB) 030329 with the 1.5m Russian-Turkish telescope RTT150 (TUBITAK National Observatory, Bakyrlytepe, Turkey). RTT150 was one of the first 1.5m-class telescopes pointed to the afterglow. Observations were started approximately 6 hours after the burst. During the first 5 hours of our observations the afterglow faded exactly as a power law with index -1.19+-0.01 in each of the BVRI Bessel filters. After that, in all BVRI filters simultaneously we observe a steepening of the power law light curve. The power law decay index smoothly approaches the value ~= -1.9, observed by other observatories later. This power law break occurs at t-t_0 =0.57 days and lasts for +-0.1 days. We observe no variability above the gradual fading with the upper limits 10--1% on time scales 0.1--1000s. Spectral flux distribution in four BVRI filters corresponds to the power law spectrum with spectral index \alpha=0.66+-0.01. The change of the power law decay index in the end of our observations can be interpreted as a signature of collimated ultrarelativistic jet. The afterglow flux distribution in radio, optical and x-rays is consistent with synchrotron spectrum. We continue our observations of this unique object with RTT150.

We show that the giant flares of soft-gamma ray repeaters (E~10^44 ergs) can push the inner regions of a fall-back disk out to larger radii by radiation pressure, while matter remains bound to the system for plausible parameters. The subsequent relaxation of this pushed-back matter can account for the observed enhanced X-ray emission after the August 27 giant flare of SGR 1900+14. Based on the results of our models, we estimate that the ratio of the fluences of the enhanced X-ray emissions to that of the preceding bursts remains constant for a particular SGR with similar preburst inner disk conditions, which is consistent with the four different burst observations of SGR 1900+14.

Interaction with possible fallback material, along with the magnetic fields and rotation rates at birth should determine the fates and categories of young neutron stars. This paper addresses some issues related to pure or hybrid accretion models for explaining the properties of young neutron stars.

We suggest a new scenario to explain the outburst light curves of black hole soft X-ray transients together with the secondary maximum and the bump seen on their decay phases. Our explanations are based on the disk instability models considering the effect of X-ray irradiation. The scenario is consistent with the observed X-ray delays by a few days with respect to the optical rise for both the main outburst and the later maxima. We test our scenario by numerically solving the disk diffusion equation. The obtained model curve fits well to the observed X-ray outburst photon flux curve of the black hole soft X-ray transient GS/GRS 1124-68, a typical representative of the four BH SXTs with very similar light curves.

Large pulsar glitches seem to be common to all radio pulsars and to exhibit a universal behaviour connecting the rate of occurrence, event size and interglitch relaxation that can be explained if the glitches are due to angular momentum exchange in the neutron star. This has implications for the energy dissipation rate. The large timing excursions observed in SGRs are not similar to the pulsar glitches. The one glitch observed in an AXP is very similar to pulsar glitches and demonstrates that typical neutron star glitch response occurs when the external torque is constant for an extended period of time. In the long run the timing variations of both the AXPs and the SGRs are much stronger than the timing noise of radio pulsars but comparable with the behaviour of accreting sources. For the DTNs, the thermal luminosities do not seem compatible with cooling ages and statistics if these sources are magnetars. The energy dissipation rates implied by the dynamical behaviour of glitching pulsars will provide the luminosities of the DTNs under propeller spindown torques supplied by fossil disks on neutron stars with conventional 10$^{12}$ G dipole magnetic fields.

Genetic connections between anomalous X-ray pulsars (AXPs) or soft gamma ray repeaters (SGRs)and supernova remnants (SNRs) can hold important clues to the nature of these objects through the properties of the associated SNRs. We examine various criteria jointly to assess the reliability of possible connections. 3 AXPs and 2 SGRs are most probably genetically connected with SNRs. Distances and ages are discussed and their likely ranges are indicated.

Neutron stars may be surrounded by fall-back disks formed from supernova core-collapse. If the disk circumscribes the light-cylinder, the neutron star will be an active radio pulsar spinning down under the propeller spin-down torque applied by the disk as well as the usual magnetic dipole radiation torque. Evolution across the P-P_dot diagram is very rapid when pulsar spin-down is dominated by the propeller torque. This explains the distribution of pulsars in the P-P_dot diagram.

Using X-ray data from the Rossi X-Ray Timing Explorer (RXTE), we report the pulse timing results of the accretion powered high mass X-ray binary (HMXRB) pulsar 4U 1907+09 covering a time span of almost two years. We measured three new pulse periods in addition to the previously measured four pulse periods. We are able to connect pulse arrival times in phase for more than a year. The source has been spinning down almost at a constant rate with a spin down rate of \dot \nu = (-3.54 \pm 0.02) \times 10^{-14} Hz s ^{-1} for more than 15 years. Residuals of pulse arrival times yield a very low level of random walk noise strengths \sim 2 \times 10^{-20} rad ^{2} sec ^{-3} on a time scale of 383 days, which is four decades lower than that of the HMXRB pulsar Vela X-1. The noise strength is only a factor of 5 greater than that of the low mass X-ray binary pulsar (LMXRB) 4U 1626-67. The low level of the timing noise and the very stable spin down rate of 4U 1907+09 makes this source unique among the HMXRBs, providing another example, in addition to 4U 1626-67, of long term quiet spin down from an accreting source. These examples show that the extended quiet spin down episodes observed in the anomalous X-ray pulsars (AXPs) pulsars 1RXS J170849.0-400910 and 1E 2259+586 do not necessarly imply that these sources are not accreting pulsars.

The similarity of rotation periods of, the anomalous X-ray pulsars (AXPs), the soft gamma ray repeaters (SGRs) and the dim thermal neutron stars (DTNs) suggests a common mechanism with an asymptotic spindown phase through the propeller and early accretion stages. The DTNs are in the propeller stage. Their luminosities arise from frictional heating in the neutron star. If the 8.4 s rotation period of the DTN RXJ 0720.4-3125 is close to its rotational equilibrium period, the propeller torque indicates a magnetic field in the 10$^{12}$ Gauss range. The mass inflow rate onto the propeller is of the order of the AXP accretion rates. The limited range of rotation periods, taken to be close to equilibrium periods, and magnetic fields in the range 5 E11- 5 E12 Gauss correspond to mass inflow rates 3.2 E14 gm/s < \dot{M} < 4.2 E17 gm/s. Observed spindown rates of the AXPs and SGRs also fit in with these fields rather than magnetar fields periods. The source of the mass inflow is a remnant accretion disk formed as part of the fallback during the supernova explosion. These classes of sources thus represent the alternative pathways for those neutron stars that do not become radio pulsars. For the highest mass inflow rates the propeller action may support enough circumstellar material so that the optical thickness to electron scattering destroys the X-ray beaming, and the rotation period is not observable. These are the radio quiet neutron stars (RQNSs) at the centers of supernova remnants Cas A, Puppis A, RCW 103 and 296.5+10.

The ROSAT HRI was used to monitor X-ray emission from the Vela Pulsar. Six observations span 2-1/2 years and 3 glitches. The summed data yield a determination of the pulse shape, and X-ray emission from the pulsar is found to be 12 % pulsed with one broad and two narrow peaks. One observation occurred 15 days after a large glitch. No change in pulse structure was observed and any change in X-ray luminosity, if present, was less than 3 %. Implications for neutron star structure are discussed.

We report the results of analysis of observations of the Vela Pulsar by PCA on RXTE. Our data consists of two parts. The first part contains observations at 1, 4, and 9 days after the glitch in 1996 and has 27000 sec. total exposure time. The second part of observations were performed three months after this glitch and have a total exposure time of 93000 sec. We found pulsations in both sets. The observed spectrum is a power-law with no apparent change in flux or count rate. The theoretical expectations of increase in flux due to internal heating after a glitch are smaller than the uncertainty of the observations.

The similarity of rotation periods from three interesting classes of neutron stars, the anomalous X-ray pulsars (AXPs), the soft gamma ray repeaters (SGRs) and the dim isolated thermal neutron stars (DTNs) suggests a common mechanism with an asymptotic spindownphase, extending through the propeller and early accretion stages. The DTNs are interpreted as sources in the propeller stage. Their low luminosities arise from frictional heating in the neutron star. SGRs and AXPs are accreting at $\dot{M} \sim 10^{15} gm/s $. The limited range of near equilibrium periods corresponds to a limited range of mass inflow rates $\dot{M}$. For lower rates the source of mass inflow may be depleted before the asymptotic stage is reached, while sources with higher $\dot{M}$ or later ages possess circumstellar material that is optically thick to electron scattering, destroying the X-ray beaming and the modulation at the rotation period. The model works with conventional magnetic fields of 10$^{11}-10^{12}$ G, obviating the need to postulate magnetars. Frequently sampled timing observations of AXPs, SGRs and DTNs can distinguish between this explanation and the magnetar model.

We have investigated the stability of the pulse frequency second derivatives ($\ddot \nu $) of PSR 0823+26, PSR 1706-16, PSR 1749-28, PSR 2021+51 which show significant quadratic trends in their pulse frequency histories in order to determine whether the observed second derivatives are secular or they arise as part of noise processes. We have used TOA data extending to more than three decades which are the longest time spans ever taken into account in pulse timing analyses. We investigated the stability of pulse frequency second derivative in the framework of low resolution noise power spectra (Deeter 1984) estimated from the residuals of pulse frequency and TOA data. We have found that the $\ddot \nu $ terms of these sources arise from the red torque noise in the fluctuations of pulse frequency derivatives which may originate from the external torques from the magnetosphere of pulsar.

We note that the solution of accretion disk dynamics for an initial delta-function mass distribution gives a light curve that fits both the rise and the decay pattern of the outburst light curves of black-hole soft X-ray transients (BSXTs) until the onset of the first mini outburst quite well. The Green's function solution of Lynden-Bell & Pringle (1974) is employed for two differenttime-independent viscosity laws to calculate the expected count rates of X-ray photons in the Ginga energy bands as a function of time. For both models basic characteristics of the outburst light curves of two typical sources GS 2000+25 and GS/GRS 1124-68 are reproduced together with plausible values of the thin disk parameter $\alpha$ and the recurrence times. This agreement with the outburst light curves and the source properties during quiescence support the idea of mass accumulation and the sporadic release of accumulated mass at the outer disk.

Using the archival ROSAT observation of TT Ari, X-ray energy spectra in different orbital phases and power spectra of the intensity time series are presented. Spectral fits show that the source gets brighter during the observation. The orbital modulation of the X-ray counting rate and bremsstrahlung temperature suggests that soft X-ray emission peaks in the orbital phase interval 0.75-0.90, when an outer disk hot spot is near the line of sight. This correlates with the orbital modulation of C IV($\lambda$1549) absorption. Timing analysis indicates that while the source gets brighter, the frequency of the 1mHz oscillation is not correlated with X-ray intensity. This implies that in the X-rays from TT Ari, the beat frequency model is not appropriate for explaining the changes in the 1mHz oscillations. \keywords{Accretion discs - stars: TT Ari - stars: cataclysmic variables - X-ray: binaries - X-rays

We report observations of PSR's 0740-28, 1737-30, 1822-09, 1915+13 and 1916+14 with ROSAT. In the 0.1-2.1 keV range upper limits to luminosity are derived for power law and blackbody spectra, using a range of $N_{H}$ estimates. The upper limit to the blackbody luminosity from PSR1822-09 turns out to be consistent with standard cooling curves. For the other pulsars the upper limits are not restrictive as they are much larger than the luminosities predicted by the models.