Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe Hydrodynamic Origin of Neutron Star Kicks
105
0
0.0
(
0
)
Authors
J. Nordhaus
Ask ChatGPT about the research
No Arabic abstract
We present results from a suite of axisymmetric, core-collapse supernova simulations in which hydrodynamic recoil from an asymmetric explosion produces large proto-neutron star (PNS) velocities. We use the adaptive-mesh refinement code CASTRO to self-consistently follow core-collapse, the formation of the PNS and its subsequent acceleration. We obtain recoil velocities of up to 620 km/s at ~1 s after bounce. These velocities are consistent with the observed distribution of pulsar kicks and with PNS velocities obtained in other theoretical calculations. Our PNSs are still accelerating at several hundred km/s at the end of our calculations, suggesting that even the highest velocity pulsars may be explained by hydrodynamic recoil in generic, core-collapse supernovae.
rate research
Read More
In Bray and Eldridge (2017), we proposed a simple neutron star kick formula, v kick = alpha (M ejecta / M remnant) + beta to explain the observed 2D velocities of young single neutron stars. Using this kick we found that there is no statistically significant preference for a kick orientation nor for any of the three initial mass function (IMF) slopes tested, and that populations including binary stars reproduced the kick distribution better than single star only populations. However, recent analysis by Janka (2017), prompted us to revisit our basic assumptions and our new analysis has led to revised best-fit kick values of alpha=100 km per second and beta = -170 km per second. The reduction of beta to a negative value is due to using the 2D observed kick velocity distribution rather than the modelled 3D velocity distribution for neutron stars (NS). To further test the validity of the new kick, we have created synthetic populations of runaway star and double neutron star (DNS) binaries at solar metallicity (Z=0.02) using our best-fit kick. We find our new kick values create runaway star velocities and DNS period distributions in agreement with the comparable observational distributions with only the DNS eccentricities in tension with the observations. From our DNS and BH-BH datasets we estimate a predicted DNS merger rate at solar metallicity of 3,864 (+1,570/-2,371) per cubic Gpc per yr and a BH-BH merger rate of 5 (+40/-1) per cubic Gpc per yr.
The first neutron star-neutron star (NS-NS) merger was discovered on August 17, 2017 through gravitational waves (GW170817) and followed with electromagnetic observations. This merger was detected in an old elliptical galaxy with no recent star formation. We perform a suite of numerical calculations to understand the formation mechanism of this merger. We probe three leading formation mechanisms of double compact objects: classical isolated binary star evolution, dynamical evolution in globular clusters and nuclear cluster formation to test whether they are likely to produce NS-NS mergers in old host galaxies. Our simulations with optimistic assumptions show current NS-NS merger rates at the level of 10^-2 yr^-1 from binary stars, 5 x 10^-5 yr^-1 from globular clusters and 10^-5 yr^-1 from nuclear clusters for all local elliptical galaxies (within 100 Mpc^3). These models are thus in tension with the detection of GW170817 with an observed rate 1.5 yr^-1 (per 100 Mpc^3; LIGO/Virgo estimate). Our results imply that either (i) the detection of GW170817 by LIGO/Virgo at their current sensitivity in an elliptical galaxy is a statistical coincidence; or that (ii) physics in at least one of our three models is incomplete in the context of the evolution of stars that can form NS-NS mergers; or that (iii) another very efficient (unknown) formation channel with a long delay time between star formation and merger is at play.
The idea that gamma-ray bursts might be a kind of phenomena associated with neutron star kicks was first proposed by Dar & Plaga (1999). Here we study this mechanism in more detail and point out that the neutron star should be a high speed one (with proper motion larger than $sim 1000$ km/s). It is shown that the model agrees well with observations in many aspects, such as the energetics, the event rate, the collimation, the bimodal distribution of durations, the narrowly clustered intrinsic energy, and the association of gamma-ray bursts with supernovae and star forming regions. We also discuss the implications of this model on the neutron star kick mechanism, and suggest that the high kick speed were probably acquired due to the electromagnetic rocket effect of a millisecond magnetar with an off-centered magnetic dipole.
The collapse of a massive stars core, followed by a neutrino-driven, asymmetric supernova explosion, can naturally lead to pulsar recoils and neutron star kicks. Here, we present a two-dimensional, radiation-hydrodynamic simulation in which core collapse leads to significant acceleration of a fully-formed, nascent neutron star (NS) via an induced, neutrino-driven explosion. During the explosion, a ~10% anisotropy in the low-mass, high-velocity ejecta lead to recoil of the high-mass neutron star. At the end of our simulation, the NS has achieved a velocity of ~150 km s$^{-1}$ and is accelerating at ~350 km s$^{-2}$, but has yet to reach the ballistic regime. The recoil is due almost entirely to hydrodynamical processes, with anisotropic neutrino emission contributing less than 2% to the overall kick magnitude. Since the observed distribution of neutron star kick velocities peaks at ~300-400 km s$^{-1}$, recoil due to anisotropic core-collapse supernovae provides a natural, non-exotic mechanism with which to obtain neutron star kicks.
We propose a simple model to explain the velocity of young neutron stars. We attempt to confirm a relationship between the amount of mass ejected in the formation of the neutron star and the `kick velocity imparted to the compact remnant resulting from the process. We assume the velocity is given by $v_{rm kick}=alpha,(M_{rm ejecta} / M_{rm remnant}) + beta,$. To test this simple relationship we use the BPASS (Binary Population and Spectral Synthesis) code to create stellar population models from both single and binary star evolutionary pathways. We then use our Remnant Ejecta and Progenitor Explosion Relationship (REAPER) code to apply different $alpha$ and $beta$ values and three different `kick orientations then record the resulting velocity probability distributions. We find that while a single star population provides a poor fit to the observational data, the binary population provides an excellent fit. Values of $alpha=70, {rm km,s^{-1}}$ and $beta=110,{rm km,s^{-1}}$ reproduce the cite{RN165} observed 2-dimensional velocities and $alpha=70, {rm km,s^{-1}}$ and $beta=120,{rm km,s^{-1}}$ reproduce their inferred 3-dimensional velocity distribution for nearby single neutron stars with ages less than 3 Myrs. After testing isotropic, spin-axis aligned and orthogonal to spin-axis `kick orientations, we find no statistical preference for a `kick orientation. While ejecta mass cannot be the only factor that determines the velocity of supernovae compact remnants, we suggest it is a significant contributor and that the ejecta based `kick should replace the Maxwell-Boltzmann velocity distribution currently used in many population synthesis codes.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
