Distribution Of Natal Kicks Research Articles

Cannonball or Bowling Ball: Proper Motion and Parallax for PSR J0002+6216

We report the results of careful astrometric measurements of the cannonball pulsar J0002+6216 carried out over 3 yr using the High Sensitivity Array. We significantly refine the proper motion to μ = 35.3 ± 0.6 mas yr−1 and place new constraints on the distance, with the overall effect of lowering the velocity and increasing the inferred age to 47.60 ± 0.80 kyr. Although the pulsar is brought more in line with the standard natal kick distribution, this new velocity has implications for the morphology of the pulsar wind nebula that surrounds it, the density of the interstellar medium through which it travels, and the age of the supernova remnant (CTB 1) from which it originates.

Evolution of eccentric high-mass X-ray binaries

Context. The formation of neutron stars is associated with powerful astrophysical transients such as supernovae. In many cases, asymmetries in the supernova explosions are thought to be responsible for the large observed velocities of neutron stars. Aims. We aim to study the complete evolutionary history of one particular eccentric high-mass X-ray binary containing a neutron star, GX 301–2, and characterize the natal kick at the time of neutron star formation. Methods. We used the publicly available stellar-evolution code MESA to evolve binaries from their initial stages until the core-collapse scenario. We incorporated a natal kick distribution based on observations to continue the evolution during the X-ray binary phase and search for candidates matching current observations of GX 301–2. Results. We find that the range of initial masses is constrained to be less than around 30 M⊙ depending on the initial mass ratio, as higher initial masses will most likely end up producing a black hole. In the completely conservative mass-transfer scenario under study, only is an interaction between the stars when the donor is still burning Hydrogen in its core, the so-called Case A of mass transfer, able to produce progenitors for GX 301–2. The natal kick study favours kicks of variable strength, which in turn increases the tilt angle between the orbital angular momentum and the spin of the neutron star. Conclusions. We conclude that only a narrow initial progenitor parameter space is able to produce a binary such as GX 301–2 when assuming a completely conservative mass transfer. Additionally, the strength of the natal kick can span a wide range of values, but it can be constrained when considering new data concerning the systemic velocity of the binary. Finally, we derive the fraction of the expected number of binaries such as GX 301–2 in the Galaxy to be ∼6 × 10−5, implying a really low chance of finding a binary similar to GX 301–2.

Binary Natal Kicks in the Galactic Center: X-Ray Binaries, Hypervelocity Stars, and Gravitational Waves

Theoretical and observational studies suggest that stellar binaries exist in large numbers in galactic nuclei like our own Galactic Center. Neutron stars (NSs), and debatedly, black holes and white dwarfs, receive natal kicks at birth. In this work, we study the effect of two successive natal kicks on a population of stellar binaries orbiting the massive black hole (MBH) in our Galactic Center. These natal kicks can significantly alter the binary orbit in a variety of ways, and also the orbit of the binary around the MBH. We found a variety of dynamical outcomes resulting from these kicks, including a steeper cusp of single NSs relative to the initial binary distribution. Furthermore, hypervelocity star and binary candidates, including hypervelocity X-ray binaries, are a common outcome of natal kicks. In addition, we show that the population of X-ray binaries in the Galactic Center can be used as a diagnostic for the BH natal kick distribution. Finally, we estimate the rate of gravitational wave events triggered by natal kicks, including binary mergers and EMRIs.

The observed velocity distribution of young pulsars – II. Analysis of complete PSRπ

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.

Tidal disruption events on to stellar black holes in triples

ABSTRACT Stars passing too close to a black hole can produce tidal disruption events (TDEs), when the tidal force across the star exceeds the gravitational force that binds it. TDEs have usually been discussed in relation to massive black holes that reside in the centres of galaxies or lurk in star clusters. We investigate the possibility that triple stars hosting a stellar black hole (SBH) may be sources of TDEs. We start from a triple system made up of three main-sequence stars and model the supernova (SN) kick event that led to the production of an inner binary comprised of an SBH. We evolve these triples with a high-precision N-body code and study their TDEs as a result of Kozai–Lidov oscillations. We explore a variety of distributions of natal kicks imparted during the SN event, various maximum initial separations for the triples, and different distributions of eccentricities. We show that the main parameter that governs the properties of the SBH–MS binaries that produce a TDE in triples is the mean velocity of the natal kick distribution. Smaller σ’s lead to larger inner and outer semimajor axes of the systems that undergo a TDE, smaller SBH masses, and longer time-scales. We find that the fraction of systems that produce a TDE is roughly independent of the initial conditions, while estimate a TDE rate of $2.1\times 10^{-4}{\!-\!}4.7 \, \mathrm{yr}^{-1}$, depending on the prescriptions for the SBH natal kicks. This rate is almost comparable to the expected TDE rate for massive black holes.

On the formation history of Galactic double neutron stars

Double neutron stars (DNSs) have been observed as Galactic radio pulsars, and the recent discovery of gravitational waves from the DNS merger GW170817 adds to the known DNS population. We perform rapid population synthesis of massive binary stars and discuss model predictions, including formation rates, mass distributions, and delay time distributions. We vary assumptions and parameters of physical processes such as mass transfer stability criteria, supernova kick distributions, remnant mass distributions and common-envelope energetics. We compute the likelihood of observing the orbital period-eccentricity distribution of the Galactic DNS population under each of our population synthesis models, allowing us to quantitatively compare the models. We find that mass transfer from a stripped post-helium-burning secondary (case BB) onto a neutron star is most likely dynamically stable. We also find that a natal kick distribution composed of both low (Maxwellian $\sigma=30\rm~km~s^{-1}$) and high ($\sigma=265\rm~km~s^{-1}$) components is preferred over a single high-kick component. We find that the observed DNS mass distribution can place strong constraints on model assumptions.

The black hole retention fraction in star clusters

Context. Recent research has been constraining the retention fraction of black holes (BHs) in globular clusters by comparing the degree of mass segregation with N-body simulations. They are consistent with an upper limit of the retention fraction being 50% or less. Aims. In this work, we focus on direct simulations of the dynamics of BHs in star clusters. We aim to constrain the effective distribution of natal kicks that BHs receive during supernova (SN) explosions and to estimate the BH retention fraction.Methods. We used the collisional N-body code nbody6 to measure the retention fraction of BHs for a given set of parameters, which are: the initial mass of a star cluster, the initial half-mass radius, and σBH, which sets the effective Maxwellian BH velocity kick distribution. We compare these direct N-body models with our analytic estimates and newest observational constraints. Results. The numerical simulations show that for the one-dimensional velocity kick dispersion σBH < 50 km s−1, clusters with radii of 2 pc and that are initially more massive than 5 × 103 M⊙ retain more than 20% of BHs within their half-mass radii. Our simple analytic model yields a number of retained BHs that is in good agreement with the N-body models. Furthermore, the analytic estimates show that ultra-compact dwarf galaxies should have retained more than 80% of their BHs for σBH ≤ 190 km s−1. Although our models do not contain primordial binaries, in the most compact clusters with 103 stars, we have found evidence of delayed SN explosions producing a surplus of BHs compared to the IMF due to dynamically formed binary stars. These cases do not occur in the more populous or expanded clusters.

Explaining LIGO’s observations via isolated binary evolution with natal kicks

We compare binary evolution models with different assumptions about black-hole natal kicks to the first gravitational-wave observations performed by the LIGO detectors. Our comparisons attempt to reconcile merger rate, masses, spins, and spin-orbit misalignments of all current observations with state-of-the-art formation scenarios of binary black holes formed in isolation. We estimate that black holes (BHs) should receive natal kicks at birth of the order of $\sigma\simeq 200$ (50) km/s if tidal processes do (not) realign stellar spins. Our estimate is driven by two simple factors. The natal kick dispersion $\sigma$ is bounded from above because large kicks disrupt too many binaries (reducing the merger rate below the observed value). Conversely, the natal kick distribution is bounded from below because modest kicks are needed to produce a range of spin-orbit misalignments. A distribution of misalignments increases our models' compatibility with LIGO's observations, if all BHs are likely to have natal spins. Unlike related work which adopts a concrete BH natal spin prescription, we explore a range of possible BH natal spin distributions. Within the context of our models, for all of the choices of $\sigma$ used here and within the context of one simple fiducial parameterized spin distribution, observations favor low BH natal spin.

The cosmic merger rate of stellar black hole binaries from the Illustris simulation

The cosmic merger rate density of black hole binaries (BHBs) can give us an essential clue to constraining the formation channels of BHBs, in light of current and forthcoming gravitational wave detections. Following a Monte Carlo approach, we couple new population-synthesis models of BHBs with the Illustris cosmological simulation, to study the cosmic history of BHB mergers. We explore six population-synthesis models, varying the prescriptions for supernovae, common envelope, and natal kicks. In most considered models, the cosmic BHB merger rate follows the same trend as the cosmic star formation rate. The normalization of the cosmic BHB merger rate strongly depends on the treatment of common envelope and on the distribution of natal kicks. We find that most BHBs merging within LIGO's instrumental horizon come from relatively metal-poor progenitors (<0.2 Zsun). The total masses of merging BHBs span a large range of values, from ~6 to ~82 Msun. In our fiducial model, merging BHBs consistent with GW150914, GW151226 and GW170104 represent ~6, 3, and 12 per cent of all BHBs merging within the LIGO horizon, respectively. The heavy systems, like GW150914, come from metal-poor progenitors (<0.15 Zsun). Most GW150914-like systems merging in the local Universe appear to have formed at high redshift, with a long delay time. In contrast, GW151226-like systems form and merge all the way through the cosmic history, from progenitors with a broad range of metallicities. Future detections will be crucial to put constraints on common envelope, on natal kicks, and on the BHB mass function.

On the impact of neutron star binaries’ natal-kick distribution on the Galactic r-process enrichment

We study the impact of the neutron star binaries' (NSBs) natal kick distribution on the Galactic r-process enrichment. We model the growth of a Milky Way type halo based on N-body simulation results and its star formation history based on multi epoch abundance matching techniques. We consider the NSBs that merge well beyond the galaxy's effective radius ($>2\times R_\mathrm{eff}$) do not contribute to Galactic r-process enrichment. Assuming a power-law delay-time distribution (DTD) function ($\propto t^{-1}$) with $t_\mathrm{min}=30$ Myr for binaries' coalescence timescales, and an exponential profile for their natal kick distribution with an average value of 180 km s$^{-1}$, we show that up to $\sim$ 40% of all formed NSBs do not contribute to r-process enrichment by $z=0$, either because they merge far from the galaxy at a given redshift (up to $\sim$ 25%) or have not yet merged by today ($\sim$15%). Our result is largely insensitive to the details of the DTD function. Assuming a constant coalescence timescale of 100 Myr well approximates the adopted DTD with 30% of the NSBs not contributing to r-process enrichment. Our results, although rather dependent on the adopted natal kick distribution, represent a first step towards estimating the impact of natal kicks and DTD functions on r-process enrichment of galaxies that would need to be incorporated in the hydrodynamical simulations.

Lidov–Kozai Cycles with Gravitational Radiation: Merging Black Holes in Isolated Triple Systems

Abstract We show that a black-hole binary with an external companion can undergo Lidov–Kozai cycles that cause a close pericenter passage, leading to a rapid merger due to gravitational-wave emission. This scenario occurs most often for systems in which the companion has a mass comparable to the reduced mass of the binary and the companion orbit has a semimajor axis within a factor of ∼10 of the binary semimajor axis. Using a simple population-synthesis model and three-body simulations, we estimate the rate of mergers in triple black-hole systems in the field to be about six per Gpc3 per year in the absence of natal kicks during black-hole formation. This value is within the low end of the 90% credible interval for the total black hole–black hole merger rate inferred from the current LIGO results. There are many uncertainties in these calculations, the largest of which is the unknown distribution of natal kicks. Even modest natal kicks of 40 km s−1 will reduce the merger rate by a factor of 40. A few percent of these systems will have eccentricity greater than 0.999 when they first enter the frequency band detectable by aLIGO (above 10 Hz).

Effects of Stellar-Mass Black Holes on Massive Star Cluster Evolution

AbstractRecent observations have revealed the existence of stellar mass black hole (BH) candidates in some globular clusters (GC) in the Milky Way and in other galaxies. Given that the detection of BHs is challenging, these detections likely indicate the existence of large populations of BHs in these clusters. This is in direct contrast to the past understanding that at most a handful of BHs may remain in old GCs due to quick mass segregation and rapid mutual dynamical ejection. Modern realistic star-by-star numerical simulations suggest that the retention fraction of BHs is typically much higher than what was previously thought. The BH dynamics near the cluster center leads to dynamical formation of new binaries and dynamical ejections, and acts as a persistent and significant energy source for these clusters. We have started exploring effects of BHs on the global evolution and survival of star clusters. We find that the evolution as well as survival of massive star clusters can critically depend on the details of the initial assumptions related to BH formation physics, such as natal kick distribution, and the initial stellar mass function (IMF). In this article we will present our latest results.

Investigating stellar-mass black hole kicks

We investigate whether stellar-mass black holes have to receive natal kicks in order to explain the observed distribution of low-mass X-ray binaries containing black holes within our Galaxy. Such binaries are the product of binary evolution, where the massive primary has exploded forming a stellar-mass black hole, probably after a common envelope phase where the system contracted down to separations of order 10-30 Rsun. We perform population synthesis calculations of these binaries, applying both kicks due to supernova mass-loss and natal kicks to the newly-formed black hole. We then integrate the trajectories of the binary systems within the Galactic potential. We find that natal kicks are in fact necessary to reach the large distances above the Galactic plane achieved by some binaries. Further, we find that the distribution of natal kicks would seem to be similar to that of neutron stars, rather than one where the kick velocities are reduced by the ratio of black hole to neutron-star mass (i.e. where the kicks have the same momentum). This result is somewhat surprising; in many pictures of stellar-mass black-hole formation, one might have expected black holes to receive kicks having the same momentum (rather than the same speed) as those given to neutron stars.