ABSTRACT
Understanding the natal kicks, or birth velocities, of neutron stars is essential for understanding the evolution of massive binaries and double neutron star formation. We use maximum ...likelihood methods as published in Verbunt et al. to analyse a new large data set 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 τ < 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 σ1 and σ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 σ2 = 298 ± 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 σ2 = 336 ± 45 km s−1. In the young population, 5 ± 3 per cent of pulsars have 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.
Abstract
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) that ...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 natal kick 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.
Abstract
The speeds of young isolated pulsars are generally inferred from their observed 2D velocities on the plane of the sky under the assumption that the unobserved radial velocity is not special, ...i.e., that the measured 2D velocity is an isotropic projection of the full 3D velocity. However, if pulsar spins are preferentially aligned with kicks, then the observer’s viewing angle relative to the pulsar velocity vector is in fact special because the direction of the spin impacts the detectability of the pulsar. This means that the measured 2D velocity of observable pulsars is not an isotropic projection, which affects inference on 3D velocities. We estimate this effect and conclude that it could lead to a ∼15% systematic overestimate of neutron-star natal kicks if young pulsars have high obliquity angles and narrow beams, but the exact correction factor depends on the distribution of beam-spin and spin-kick misalignment angles and beam widths.
ABSTRACT
Supernova explosion and the associated neutron star (NS) natal kicks are important events on a pathway of a binary to become a gravitational wave source, an X-ray binary, or a millisecond ...radio pulsar. Weak natal kicks often lead to binary survival, while strong kicks frequently disrupt the binary. In this article, we aim to further constrain NS natal kicks in binaries. We explore binary population synthesis models by varying prescription for natal kick, remnant mass, and mass accretion efficiency. We introduce a robust statistical technique to analyse combined observations of different nature. Using this technique, we further test different models using parallax and proper motion measurements for young isolated radio pulsars and similar measurements for Galactic Be X-ray binaries (BeXs). Our best model for natal kicks is consistent with both measurements and contains a fraction of w = 0.2 ± 0.1 weak natal kicks with $\sigma _1 = 45^{+25}_{-15}$ km s−1, the remaining natal kicks are drawn from the high-velocity component, same as in previous works: σ2 = 336 km s−1. We found that currently used models for natal kicks of NSs produced by electron capture supernova (ecSN; combination of Maxwellian σ = 265 km s−1 and σ = 30 km s−1 for electron capture) are inconsistent or marginally consistent with parallaxes and proper motions measured for isolated radio pulsars. We suggest a new model for natal kicks of ecSN, which satisfy both observations of isolated radio pulsars and BeXs.
ABSTRACT
We numerically model evolution of magnetic fields inside a neutron star under the influence of ambipolar diffusion in the weak-coupling mode in the one-fluid MHD approximation. Our ...simulations are 3D and performed in spherical coordinates. Our model covers the neutron star core and includes crust where the magnetic field decay is due to Ohmic decay. We discover an instability of poloidal magnetic field under the influence of ambipolar diffusion. This instability develops in the neutron star core and grows on a time-scale of 0.2 dimensionless times, reaching saturation by 2 dimensionless times. The instability leads to formation of azimuthal magnetic field with azimuthal wavenumber m = 14 (at the moment of saturation) which keeps merging and reaches m = 4 by 16 dimensionless times. Over the course of our simulations (16 dimensionless times) the surface dipolar magnetic field decays, reaching 20 per cent of its original value and keeps decaying. The decay time-scale for the total magnetic energy is six dimensionless times. The ambipolar diffusion induces electric currents in the crust where these currents dissipate efficiently. Strong electric currents in the crust lead to heating, which could correspond to luminosities of ≈1029 erg s−1 during hundreds of Myrs for an initial magnetic field of 1014 G. Ambipolar diffusion leads to formation of small-scale magnetic fields at the neutron star surface.
ABSTRACT
Off-centred dipole configurations have been suggested to explain different phenomena in neutron stars, such as natal kicks, irregularities in polarization of radio pulsars and properties of ...X-ray emission from millisecond pulsars. Here, for the first time, we model magnetothermal evolution of neutron stars with crust-confined magnetic fields and off-centred dipole moments. We find that the dipole shift decays with time if the initial configuration has no toroidal magnetic field. The decay time-scale is inversely proportional to magnetic field. The octupole moment decreases much faster than the quadrupole. Alternatively, if the initial condition includes strong dipolar toroidal magnetic field, the external poloidal magnetic field evolves from centred dipole to off-centred dipole. The surface thermal maps are very different for configurations with weak B = 1013 G and strong B = 1014 G magnetic fields. In the former case, the magnetic equator is cold while in the latter case, it is hot. We model light curves and spectra of our magnetothermal configurations. We found that in the case of cold equator, the pulsed fraction is small (below a few per cent in most cases) and spectra are well described with a single blackbody. Under the same conditions, models with stronger magnetic fields produce light curves with pulsed fraction of tens of per cent. Their spectra are significantly better described with two blackbodies. Overall, the magnetic field strength has a more significant effect on bulk thermal emission of neutron stars than does the field geometry.
ABSTRACT
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 time-scales for the CEE mass-loss. In order to probe the CEE-time-scales we study wide companions to post-CE binaries. Faster mass-loss time-scales 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 ultrawide companions (projected separations 103–2 × 105 au and M2 > 0.4 M⊙) 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 ± 0.2 per cent for sdBs to be compared with a multiplicity fraction of 5.0 ± 0.2 per cent for late-B/A/F stars which are possible sdB progenitors. The distribution of projected separations of ultrawide pairs to main sequence stars and sdBs differs significantly and is compatible with prompt mass-loss (upper limit on common envelope ejection time-scales of 102 yr). The smaller statistics of ultrawide companions to cataclysmic variables and post-CEE binaries provide weaker constraints. Nevertheless, the survival rate of ultrawide pairs to the cataclysmic variables suggest much longer, ∼104 yr time-scales for the CEE in these systems, possibly suggesting non-dynamical CEE in this regime.
Neutron stars are natural physical laboratories allowing us to study a plethora of phenomena in extreme conditions. In particular, these compact objects can have very strong magnetic fields with ...non-trivial origin and evolution. In many respects, its magnetic field determines the appearance of a neutron star. Thus, understanding the field properties is important for the interpretation of observational data. Complementing this, observations of diverse kinds of neutron stars enable us to probe parameters of electro-dynamical processes at scales unavailable in terrestrial laboratories. In this review, we first briefly describe theoretical models of the formation and evolution of the magnetic field of neutron stars, paying special attention to field decay processes. Then, we present important observational results related to the field properties of different types of compact objects: magnetars, cooling neutron stars, radio pulsars, and sources in binary systems. After that, we discuss which observations can shed light on the obscure characteristics of neutron star magnetic fields and their behaviour. We end the review with a subjective list of open problems.
Abstract
Central compact objects (CCOs) are young neutron stars emitting thermal X-rays with bolometric luminosities
L
X
in the range of 10
32
–10
34
erg s
−1
. Gourgouliatos, Hollerbach, and Igoshev ...recently suggested that peculiar emission properties of CCOs can be explained by tangled magnetic field configurations formed in a stochastic dynamo during the proto–neutron star stage. In this case the magnetic field consists of multiple small-scale components with negligible contribution of global dipolar field. We study numerically three-dimensional magnetothermal evolution of tangled crustal magnetic fields in neutron stars. We find that all configurations produce complicated surface thermal patterns that consist of multiple small hot regions located at significant separations from each other. The configurations with initial magnetic energy of (2.5–10) × 10
47
erg have temperatures of hot regions that reach ≈ 0.2 keV, to be compared with the bulk temperature of ≈ 0.1 keV in our simulations with no cooling. A factor of two in temperature is also seen in observations of CCOs. The hot spots produce periodic modulations in light curve with typical amplitudes of ≤9%–11%. Therefore, the tangled magnetic field configuration can explain thermal emission properties of some CCOs.