Central Compact Objects (CCOs) are X-ray sources with luminosity ranging between $10^{32}$-$10^{34}$ erg~s$^{-1}$, located at the centres of supernova remnants. Some of them have been confirmed to be neutron stars. Timing observations have allowed the estimation of their dipole magnetic field, placing them in the range $\sim10^{10}$-$10^{11}$ G. The decay of their weak dipole fields, mediated by the Hall effect and Ohmic dissipation, cannot provide sufficient thermal energy to power their X-ray luminosity, as opposed to magnetars whose X-ray luminosities are comparable. Motivated by the question of producing high X-ray power through magnetic field decay while maintaining a weak dipole field, we explore the evolution of a crustal magnetic field that does not consist of an ordered axisymmetric structure, but rather comprises a tangled configuration. This can be the outcome of a non-self-excited dynamo, buried inside the crust by fallback material following the supernova explosion. We find that such initial conditions lead to the emergence of the magnetic field from the surface of the star and the formation of a dipolar magnetic field component. An internal tangled magnetic field of the order of $10^{14}$ G can provide sufficient Ohmic heating to the crust and power CCOs, while the dipole field it forms is approximately $10^{10}$ G, as observed in CCOs.

Understanding the natal kicks, or birth velocities, of neutron stars are essential for understanding the evolution of massive binaries as well as double neutron star formation. We use maximum likelihood methods as published in Verbunt et al. to analyse a new large dataset of parallaxes and proper motions measured by Deller et al. This sample is roughly three times larger than number of measurements available before. For both the complete sample and its younger part (spin-down ages $\tau < 3$ Myr), we find that a bimodal Maxwellian distribution describes the measured parallaxes and proper motions better than a single Maxwellian with probability of 99.3 and 95.0 per cent respectively. The bimodal Maxwellian distribution has three parameters: fraction of low-velocity pulsars and distribution parameters $\sigma_1$ and $\sigma_2$ for low and high-velocity modes. For a complete sample, these parameters are as follows: $42_{-15}^{+17}$ per cent, $\sigma_1=128_{-18}^{+22}$ km s$^{-1}$ and $\sigma_2 = 298\pm 28$ km s$^{-1}$. For younger pulsars, which are assumed to represent the natal kick, these parameters are as follows: $20_{-10}^{+11}$ per cent, $\sigma_1=56_{-15}^{+25}$ km s$^{-1}$ and $\sigma_2=336\pm 45$ km s$^{-1}$. In the young population, $5\pm 3$ per cent of pulsars has velocities less than 60 km s$^{-1}$. We perform multiple Monte Carlo tests for the method taking into account realistic observational selection. We find that the method reliably estimates all parameters of the natal kick distribution. Results of the velocity analysis are weakly sensitive to the exact values of scale-lengths of the Galactic pulsar distribution.

We investigate the extent to which resonances between an oscillating background of ultra-light axion and a binary Keplerian system can affect the motion of the latter. These resonances lead to perturbations in the instantaneous time-of-arrivals, and to secular variations in the period of the binary. While the secular changes at exact resonance have recently been explored, the instantaneous effects have been overlooked. In this paper, we examine the latter using N-body simulations including the external oscillatory forcing induced by the axion background. While the secular effects are restricted to a narrow width near the resonance, the instantaneous changes, albeit strongest close to resonances, are apparent for wide range of configurations. We compute the signal-to-noise ratio (SNR) as a function of semi-major axis for a detection of axion oscillations through the R\{o} mer delay. The latter can be extracted from the time-of-arrivals of binary pulsars. The SNR broadly increases with increasing binary eccentricity in agreement with the secular expectation. However, we find that it differs significantly from the scaling a^{5/2} around the lowest orders of resonance. Future observations could probe these effects away from resonances and, therefore, constrain a much broader range of axion masses provided that binary pulsar systems are found near the central region of our Galaxy, and that the time-or-arrival measurement accuracy reaches < 10 ns

Evolution of close binaries often proceeds through the common envelope stage. The physics of the envelope ejection (CEE) is not yet understood, and several mechanisms were suggested to be involved. These could give rise to different timescales for the CEE mass-loss. In order to probe the CEE-timescales we study wide companions to post-CE binaries. Faster mass-loss timescales give rise to higher disruption rates of wide binaries and result in larger average separations. We make use of data from Gaia DR2 to search for ultra-wide companions (projected separations $10^3$-$2\times 10^5$ a.u. and $M_2 > 0.4$ M$_\odot$) to several types of post-CEE systems, including sdBs, white-dwarf post-common binaries, and cataclysmic variables. We find a (wide-orbit) multiplicity fraction of $1.4\pm 0.2$ per cent for sdBs to be compared with a multiplicity fraction of $5.0\pm 0.2$ per cent for late-B/A/F stars which are possible sdB progenitors. The distribution of projected separations of ultra-wide pairs to main sequence stars and sdBs differs significantly and is compatible with prompt mass loss (upper limit on common envelope ejection timescale of $10^2$ years). The smaller statistics of ultra-wide companions to cataclysmic variables and post-CEE binaries provide weaker constraints. Nevertheless, the survival rate of ultra-wide pairs to the cataclysmic variables suggest much longer, $\sim10^4$ years timescales for the CEE in these systems, possibly suggesting non-dynamical CEE in this regime.

The origin of ultra-wide massive binaries (orbital separations $10^3-2\times 10^5$~AU) and their properties are not well characterized nor understood. Here we use the second Gaia data release to search for wide astrometric companions to Galactic O-B5 stars which share similar parallax and proper motion with the primaries. Using the data we characterize the frequency and properties of such binaries. We find an ultra-wide multiplicity fraction of $4.4\pm0.5$ per cent, to our completeness limit (up to $\approx 17$~mag; down to G-stars at distances of 0.3-2~kpc, excluding stars in clusters). The secondary mass-function is generally consistent with a Kroupa initial stellar function; if extrapolated to lower mass companion stars we then might expect a wide-binary fraction of $\sim 27\pm5\%$. In addition we use these data as a verification sample to test the existence of ultra-wide binaries among neutron stars (NSs) and black holes (BHs). We propose that the discovery of such binary can provide unique constraints on the weakest natal kicks possible for NSs/BHs. If a compact object is formed in an ultra-wide binary and receives a very-low natal kick, such a binary should survive as a common proper motion pair. We therefore use Gaia data to search for ultra-wide companions to pulsars (normal and millisecond ones) and X-ray binaries. We find no reliable pairs. Future data could potentially provide stringent constraints through this method.

In this note we propose that recently discovered radio pulsar J0250+5854 with 23.5 sec spin period is presently at the Hall attractor stage. This can explain low temperature and absence of magnetar-like activity of this source together with its spin period and period derivative. We present results of calculations of the evolution of this source in a simple model of magnetic field decay. The neutron star could start its evolution as a magnetar with initial field $\sim 10^{14}-10^{15}$ G for realistic range of parameter $Q$ describing crust imperfections. Future measurements of surface temperature and age of this neutron star might help to probe this hypothesis.

We use the Bayesian approach to write the posterior probability density for the three-dimensional velocity of a pulsar and for its kinematic age. As a prior, we use the bimodal velocity distribution found in a recent article by Verbunt, Igoshev & Cator (2017). When we compare the kinematic ages with spin-down ages, we find that in general, they agree with each other. In particular, maximum likelihood analysis sets the lower limit for the exponential magnetic field decay timescale at $8$ Myr with a slight preference of $t_\mathrm{dec} \approx 12$ Myr and compatible with no decay at all. One of the objects in the study, pulsar B0950+08 has kinematic and cooling ages $\approx 2$ Myr which is in strong contradiction with its spin-down age $\tau\approx 17$ Myr. The 68 per cent credible range for the kinematic age is 1.2--8.0 Myr. We conclude that the most probable explanation for this contradiction is a combination of magnetic field decay and long initial period. Further timing, UV and X-ray observations of B0950+08 are required to constrain its origin and evolution better.

The mergers of neutron stars (NSs) and white dwarfs (WDs) could give rise to explosive transients, potentially observable with current and future transient surveys. However, the expected properties and distribution of such events is not well understood. Here we characterize the rates of such events, their delay time distribution, their progenitors and the distribution of their properties. We use binary populations synthesis models and consider a wide range of initial conditions and physical processes. In particular we consider different common-envelope evolution models and different NS natal kick distributions. We provide detailed predictions arising from each of the models considered. We find that the majority of NS-WD mergers are born in systems in which mass-transfer played an important role, and the WD formed before the NS. For the majority of the mergers the WDs have a carbon-oxygen composition (60-80%) and most of the rest are with oxygen-neon WDs. The rates of NS-WD mergers are in the range of 3-15% of the type Ia supernovae (SNe) rate. Their delay time distribution is very similar to that of type Ia SNe, but slightly biased towards earlier times. They typically explode in young 0.1-1Gyr environments, but have a tail distribution extending to long, Gyrs-timescales. Models including significant kicks give rise to relatively wide offset distribution extending to hundreds of kpcs. The demographic and physical properties of NS-WD mergers suggest they are likely to be peculiar type Ic-like SNe, mostly exploding in late type galaxies. Their overall properties could be related to a class of rapidly evolving SNe recently observed, while they are less likely to be related to the class of Ca-rich SNe.

We report on the first X-ray observation of the 0.28 s isolated radio pulsar PSR J1154--6250 obtained with the XMM-Newton observatory in February 2018. A point-like source is firmly detected at a position consistent with that of PSR J1154--6250. The two closest stars are outside the 3$\sigma$ confidence limits of the source position and thus unlikely to be responsible for the observed X-ray emission. The energy spectrum of the source can be fitted equally well either with an absorbed power-law with a steep photon index $\Gamma\approx 3.3$ or with an absorbed blackbody with temperature $kT=0.21\pm 0.04$~keV and emitting radius $R_\mathrm{BB} \approx 80$ m (assuming a distance of 1.36~kpc). The X-ray luminosity of $4.4\times 10^{30}$ erg s$^{-1}$ derived with the power-law fit corresponds to an efficiency of $\eta_X = L^\mathrm{unabs}_X/\dot E = 4.5\times 10^{-3}$, similar to those of other old pulsars. The X-ray properties of PSR J1154--6250 are consistent with an old age and suggest that the spatial coincidence of this pulsar with the OB association Cru OB1 is due to a chance alignment.

Recently, numerical calculations of the magnetic field evolution in neutron stars demonstrated the possible existence of a Hall attractor, a stage at which the evolution of the field driven by the Hall cascade ends. The existence of such a stage in neutron star magnetic evolution is very important, and can be potentially probed by observations. Here we discuss three types of objects which could have reached this stage. First, we briefly describe the evolution of normal radio pulsars with ages about a few hundred thousand years. Then we analyse in more detail observations of RX J1856.5-3754, one of the Magnificent Seven, focusing on the surface temperature distribution and comparing model predictions with the temperature map inferred from X-ray observations. Finally, we discuss the necessity of the Hall attractor stage to explain the hypothetical existence of accreting magnetars. We conclude that at the moment there is no direct confirmation of the Hall attractor stage in known sources. However, more detailed observations in the near future can demonstrate existence (or absence) of this stage of the crustal magnetic field evolution.

Several candidates for accreting magnetars have been proposed recently by different authors. Existence of such systems contradicts the standard magnetic field decay scenario where a large magnetic field of a neutron star reaches $\lesssim$ few$\times 10^{13}$G at ages $\gtrsim 1$ Myr. Among other sources, the high mass X-ray binary 4U0114+65 seems to have a strong magnetic field around $10^{14}$ G. We develop a new Bayesian estimate for the kinematic age and demonstrate that 4U0114+65 has kinematic age 2.4-5 Myr ($95\%$ credential interval) since the formation of the neutron star. We discuss which conditions are necessary to explain the potential existence of magnetars in accreting high-mass binaries with ages about few Myrs and larger. Three necessary ingredients are: the Hall attractor to prevent rapid decay of dipolar field, relatively rapid cooling of the crust in order to avoid Ohmic decay due to phonons, and finally, low values of the parameter $Q$ to obtain long Ohmic time scale due to impurities. If age and magnetic field estimates for proposed accreting magnetars are correct, then these systems set the strongest limit on the crust impurity for a selected sample of neutron stars and provide evidence in favour of the Hall attractor.

We argue that comparison with observations of theoretical models for the velocity distribution of pulsars must be done directly with the observed quantities, i.e. parallax and the two components of proper motion. We develop a formalism to do so, and apply it to pulsars with accurate VLBI measurements. We find that a distribution with two maxwellians improves significantly on a single maxwellian. The `mixed' model takes into account that pulsars move away from their place of birth, a narrow region around the galactic plane. The best model has 42% of the pulsars in a maxwellian with average velocity sigma sqrt{8/pi}=120 km/s, and 58% in a maxwellian with average velocity 540 km/s. About 5% of the pulsars has a velocity at birth less than 60\,km/s. For the youngest pulsars (tau_c<10 Myr), these numbers are 32% with 130 km/s, 68% with 520 km/s, and 3%, with appreciable uncertainties.

The aim of this work is to study the imprints that different models for black hole (BH) and neutron star (NS) formation have on the Galactic distribution of X-ray binaries (XRBs) which contain these objects. We find that the root mean square of the height above the Galactic plane of BH- and NS-XRBs is a powerful proxy to discriminate among different formation scenarios, and that binary evolution following the BH/NS formation does not significantly affect the Galactic distributions of the binaries. We find that a population model in which at least some BHs receive a (relatively) high natal kick fits the observed BH-XRBs best. For the NS case, we find that a high NK distribution, consistent with the one derived from the measurement of pulsar proper motion, is the most preferable. We also analyse the simple method we previously used to estimate the minimal peculiar velocity of an individual BH-XRB at birth. We find that this method may be less reliable in the bulge of the Galaxy for certain models of the Galactic potential, but that our estimate is excellent for most of the BH-XRBs.

It has long been unclear if the small-scale magnetic structures on the neutron star (NS) surface could survive the fall-back episode. The study of the Hall cascade (Cumming, Arras and Zweibel 2004; Wareing and Hollerbach 2009) hinted that energy in small scales structures should dissipate on short timescales. Our new 2D magneto-thermal simulations suggest the opposite. For the first $\sim$10 kyrs after the fall-back episode with accreted mass $10^{-3} M_\odot$, the observed NS magnetic field appears dipolar, which is insensitive to the initial magnetic topology. In framework of the Ruderman & Sutherland (1975) vacuum gap model during this interval, non-thermal radiation is strongly suppressed. After this time the initial (i.e. multipolar) structure begins to re-emerge through the NS crust. We distinguish three evolutionary epochs for the re-emergence process: the growth of internal toroidal field, the advection of buried poloidal field, and slow Ohmic diffusion. The efficiency of the first two stages can be enhanced when small-scale magnetic structure is present. The efficient re-emergence of high order harmonics might significantly affect the curvature of the magnetospheric field lines in the emission zone. So, only after few $10^4$ yrs would the NS starts shining as a pulsar again, which is in correspondence with radio silence of central compact objects (CCOs). In addition, these results can explain the absence of good candidates for thermally emitting NSs with freshly re-emerged field among radio pulsars (Bogdanov, Ng and Kaspi 2014), as NSs have time to cool down, and supernova remnants can already dissipate.

We present a new method for probing the thermal electron content of the Galaxy by spectral analysis of background point sources in the absorption-only limit to the radiative transfer equation. In this limit, calculating the spectral index, $\alpha$, of these sources using a natural logarithm results in an additive factor, which we denote $\alpha_\mathrm{EM}$, resulting from the absorption of radiation due to the Galactic thermal electron population. We find that this effect is important at very low frequencies ($\nu\lesssim200$ MHz), and that the frequency spacing is critical. We model this effect by calculating the emission measure across the sky. A (smooth) thermal electron model for the Galaxy does not fit the observed emission measure distribution, but a simple, cloud-based model to represent the clumpy nature of the warm interstellar medium does. This model statistically reproduces the Galactic emission measure distribution as obtained independently from $H_\alpha$ data well. We find that at the lowest frequencies ($\sim10-50$ MHz), the observed spectral index for a large segment of the Galaxy below Galactic latitudes of $\lesssim15^\circ$ could be changed significantly (i.e., $\alpha_\mathrm{EM}\gtrsim0.1$). This method therefore provides a correction to low-frequency spectral index measurements of extragalactic sources, and provides a sensitive probe of the thermal electron distribution of the Galaxy using current and next-generation low-frequency radio telescopes. We show that this effect should be robustly detectable individually in the strongest sources, and statistically in source samples at a level of $\alpha_\mathrm{EM}\gtrsim0.18,0.06$, and 0.02 for source densities of 10, 100 and 1,000 sources per square degree.

We use a Bayesian approach to derive the distance probability distribution for one object from its parallax with measurement uncertainty for two spatial distribution priors, viz. a homogeneous spherical distribution and a galactocentric distribution - applicable for radio pulsars - observed from Earth. We investigate the dependence on measurement uncertainty, and show that a parallax measurement can underestimate or overestimate the actual distance, depending on the spatial distribution prior. We derive the probability distributions for distance and luminosity combined, and for each separately, when a flux with measurement error for the object is also available, and demonstrate the necessity of and dependence on the luminosity function prior. We apply this to estimate the distance and the radio and gamma-ray luminosities of PSR J0218+4232. The use of realistic priors improves the quality of the estimates for distance and luminosity, compared to those based on measurement only. Use of a wrong prior, for example a homogeneous spatial distribution without upper bound, may lead to very wrong results.

We analyse the origin of the magnetic field decay in normal radio pulsars found by us in a recent study. This decay has a typical time scale $\sim 4 \times 10^5$~yrs, and operates in the range $\sim 10^5$~--~few$\times 10^5$~yrs. We demonstrate that this field evolution may be either due to the Ohmic decay related to the scattering from phonons, or due to the Hall cascade which reaches the Hall attractor. According to our analysis the first possibility seems to be more reliable. So, we attribute the discovered field decay mainly to the Ohmic decay on phonons which is saturated at the age few$\times 10^5$~yrs, when a NS cools down to the critical temperature below which the phonon scattering does not contribute much to the resistivity of the crust. Some role of the Hall effect and attractor is not excluded, and will be analysed in our further studies.

We use a modified pulsar current analysis to study magnetic field decay in radio pulsars. In our approach we analyse the flow, not along the spin period axis as has been performed in previous studies, but study the flow along the direction of growing characteristic age, $\tau=P/(2\dot P)$. We perform extensive tests of the method and find that in most of the cases it is able to uncover non-negligible magnetic field decay (more than a few tens of per cent during the studied range of ages) in normal radio pulsars for realistic initial properties of neutron stars. However, precise determination of the magnetic field decay timescale is not possible at present. The estimated timescale may differ by a factor of few for different sets of initial distributions of neutron star parameters. In addition, some combinations of initial distributions and/or selection effects can also mimic enhanced field decay. We apply our method to the observed sample of radio pulsars at distances $<10$ kpc in the range of characteristic ages $8 \times 10^4 < \tau < 10^6$ years where, according to our study, selection effects are minimized. By analysing pulsars in the Parkes Multibeam and Swinburne surveys we find that, in this range, the field decays roughly by a factor of two. With an exponential fit this corresponds to the decay time scale $\sim 4 \times 10^5$ yrs. With larger statistics and better knowledge of the initial distribution of spin periods and magnetic field strength, this method can be a powerful tool to probe magnetic field decay in neutron stars.

The variety of the observational appearance of young isolated neutron stars must find an explanation in the framework of some unifying approach. Nowadays it is believed that such scenario must include magnetic field decay, the possibility of magnetic field emergence on a time scale $\lesssim 10^4$--$10^5$ yrs, significant contribution of non-dipolar fields, and appropriate initial parameter distributions. We present our results on the initial spin period distribution, and suggest that inconsistences between distributions derived by different methods for samples with different average ages can uncover field decay or/and emerging field. We describe a new method to probe the magnetic field decay in normal pulsars. The method is a modified pulsar current approach, where we study pulsar flow along the line of increasing characteristic age for constant field. Our calculations, performed with this method, can be fitted with an exponential decay for ages in the range $8\times 10^4$--$3.5 \times 10^5$ yrs with a time scale $\sim 5 \times 10^5$ yrs. We discuss several issues related to the unifying scenario. At first, we note that the dichotomy, among local thermally emitting neutron stars, between normal pulsars and the Magnificent Seven remains unexplained. Then we discuss the role of high-mass X-ray binaries in the unification of neutron star evolution. We note, that such systems allow to check evolutionary effects on a time scale longer than what can be probed with normal pulsars alone. We conclude with a brief discussion of importance of discovering old neutron stars accreting from the interstellar medium.

Recently, Lee et al. used Gaussian mixture models (GMM) to study the radio pulsar population. In the distribution of normal pulsars in the P-dotP plane, they found four clusters. We develop this approach further and apply it to different synthetic pulsar populations in order to determine whether the method can effectively select groups of sources that are physically different. We check several combinations of initial conditions as well as the models of pulsar evolution and the selection effects. We find that in the case of rapidly evolving objects, the GMM is oversensitive to parameter variations and does not produce stable results. We conclude that the method does not help much to identify the sub-populations with different initial parameters or/and evolutionary paths. For the same reason, the GMM does not discriminate effectively between theoretical population synthesis models of normal radio pulsars.

We analyze different possibilities to explain the wide initial spin period distribution of radio pulsars presented by Noutsos et al. (2013). With a population synthesis modeling we demonstrate that magnetic field decay can be used to interpret the difference between the recent results by Noutsos et al. (2013) and those by Popov & Turolla (2012), where a much younger population of NSs associated with supernova remnants with known ages has been studied. In particular, an exponential field decay with tau_mag = 5Myrs can produce a "tail" in the reconstructed initial spin period distribution (as obtained by Noutsos et al. 2013) up to P_0 > 1 s starting with a standard gaussian with P_0 = 0.3 s and sigma_P = 0.15 s. Another option to explain the difference between initial spin period distributions from Noutsos et al. (2013) and Popov & Turolla (2012) -- the emerging magnetic field -- is also briefly discussed.

Using the data of the ATNF pulsar catalog we study the relation connected the real age t of young neutron stars (NS) and their spin-down age \tau. We suppose that this relation is independent from both initial period of the NS and its initial surface magnetic field, and that the laws of the surface magnetic field decay are similar for all NSs in the Milky Way. We further assume that the birth-rate of pulsars was constant during at least last 200 million years. With these assumptions we were able to restore the history of the magnetic field decay for the galactic NSs. We reconstruct the universal function f(t)=B(t)/B0, where B0 is the initial magnetic field and B(t) is the magnetic field of NS at the age t. The function f(t) can be fitted by a power law with power index \alpha = -1.17.

The properties of the spin-down age are investigated. Based on assumption about a uniform magnetic field decay law we suggest a new method which allows us to shed light on magnetic field decay. This method is applied for following selection: isolated non-millisecond pulsars from the ATNF catalog are chosen. Pulsars in the selection are with the spin-down ages from 4 \cdot 10^4 to 2 \cdot 10^6 years. In order to avoid observational selection we take into account only pulsars which are closer to the Sun than 10 kpc. For this selection we restore the uniform magnetic field decay law. It appears that the magnetic field decays three times from 4 \cdot 10^4 to 3.5 \cdot 10^5 years. This function is approximated by modified power-law. We also estimate the birthrate of pulsars in our Galaxy and find that it should be about 2.9 pulsars per century.