Kick Distribution Research Articles

Brownian Particle in a Poisson‐Shot‐Noise Active Bath: Exact Statistics, Effective Temperature, and Inference

AbstractThe dynamics of an overdamped Brownian particle in a thermal bath that contains a dilute solution of active particles is studied. The particle moves in a harmonic potential and experiences Poisson shot‐noise kicks with specified amplitude distribution due to moving active particles in the bath. From the Fokker–Planck equation for the particle dynamics, the stationary solution for the displacement distribution is derived along with the moments characterizing mean, variance, skewness, and kurtosis, as well as finite‐time first and second moments. An effective temperature is also computed through the fluctuation–dissipation theorem and show that equipartition theorem holds for all zero‐mean kick distributions, including those leading to non‐Gaussian stationary statistics. For the case of Gaussian‐distributed active kicks, a re‐entrant behavior from non‐Gaussian to Gaussian stationary states and a heavy‐tailed leptokurtic distribution across a wide range of parameters are found as seen in recent experimental studies. Further analysis reveals statistical signatures of the irreversible dynamics of the particle displacement in terms of the time asymmetry of cross‐correlation functions. Fruits of the work is the development of an compact inference scheme that may allow experimentalists to extract the rate and moments of underlying shot‐noise solely from the statistics the particle position.

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.

An observationally derived kick distribution for neutron stars in binary systems

ABSTRACT Understanding the natal kicks received by neutron stars (NSs) during formation is a critical component of modelling the evolution of massive binaries. Natal kicks are an integral input parameter for population synthesis codes, and have implications for the formation of double NS systems and their subsequent merger rates. However, many of the standard observational kick distributions that are used are obtained from samples created only from isolated NSs. Kick distributions derived in this way overestimate the intrinsic NS kick distribution. For NSs in binaries, we can only directly estimate the effect of the natal kick on the binary system, instead of the natal kick received by the NS itself. Here, for the first time, we present a binary kick distribution for NSs with low-mass companions. We compile a catalogue of 145 NSs in low-mass binaries with the best available constraints on proper motion, distance, and systemic radial velocity. For each binary, we use a three-dimensional approach to estimate its binary kick. We discuss the implications of these kicks on system formation, and provide a parametric model for the overall binary kick distribution, for use in future theoretical modelling work. We compare our results with other work on isolated NSs and NSs in binaries, finding that the NS kick distributions fit using only isolated pulsars underestimate the fraction of NSs that receive low kicks. We discuss the implications of our results on modelling double NS systems, and provide suggestions on how to use our results in future theoretical works.

VFTS 243 as predicted by the BPASS fiducial models

ABSTRACT The recent discovery of an unambiguous quiescent black hole (BH) and main-sequence O star companion in VFTS 243 opens the door to new constraints on theoretical stellar evolution and population models looking to reproduce the progenitors of BH–BH binaries. Here, we show that the binary population and spectral synthesis fiducial models natively predict VFTS 243-like systems: We find that VFTS 243 likely originated from a binary system in an ∼15 d orbit with primary mass ranging from 40 to 50 $\mathrm{\, M}_\odot$ and secondary star with initial mass of 24–25 $\mathrm{\, M}_\odot$. Additionally, we find that the death of the primary star must have resulted in a low-energy explosion E < 1050 erg. With a uniform prior, we find that the kick velocity of the newborn BH was ≤10 $\, \text{km s}^{-1}$. The very low eccentricity reported for VFTS 243 and the subsequent conclusion by the authors that the supernova kick must have been very small are in line with the peak in the posterior distribution between 0 and 5 $\, \text{km s}^{-1}$ found from our numerical simulations performed with a uniform prior. Finally, the reduced Hobbs kick distribution commonly used in BH population synthesis is strongly disfavoured.

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.

Constraints to neutron-star kicks in high-mass X-ray binaries withGaiaEDR3

Context.All neutron star progenitors in neutron-star high-mass X-ray binaries (NS HMXBs) undergo a supernova event that may lead to a significant natal kick impacting the motion of the whole binary system. The space observatoryGaiaperforms a deep optical survey with exquisite astrometric accuracy, for both position and proper motions, that can be used to study natal kicks in NS HMXBs.Aims.Our aim is to survey the observed Galactic NS HMXB population and to quantify the magnitude of the kick imparted onto their NSs, and to highlight any possible differences arising between the various HMXB types.Methods.We performed a census of Galactic NS HMXBs and cross-matched existing detections in X-rays, optical, and infrared with theGaiaEarly Data Release 3 database. We retrieved their parallaxes, proper motions, and radial velocities (when available), and performed a selection based on the quality of the parallax measurement. We then computed their peculiar velocities with respect to the rotating reference frame of the Milky Way, and including their respective masses and periods, we estimated their kick velocities through Markov chain Monte Carlo simulations of the orbit undergoing a supernova event.Results.We infer the posterior kick distributions of 35 NS HMXBs. After an inconclusive attempt at characterising the kick distributions with Maxwellian statistics, we find that the observed NS kicks are best reproduced by a gamma distribution of mean 116−15+18km s−1. We note that supergiant systems tend to have higher kick velocities than Be HMXBs. The peculiar velocity versus non-degenerate companion mass plane hints at a similar trend, supergiant systems having a higher peculiar velocity independently of their companion mass.

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.

Detecting patterns in the behaviour of goalkeepers and kickers in the penalty shootout: a between-gender comparison among score situations

ABSTRACT There is plenty of research on penalty kicking in men's soccer, with a focus on either the goalkeeper or the penalty taker. Yet women's soccer and their playing behaviour are under-represented in research. The current study was designed to examine gender differences in the choice patterns of expert kickers and goalkeepers during penalty shooting in relation to the previously documented right-oriented bias. Using videos of penalties from shootouts at the highest level of men's and women's international soccer, we recorded the goalkeeper's dive direction and accurately measured the location of the ball as it crossed the goal line. We created a map of all kicks and their outcomes (goal, no goal, off-target). Our mapping procedure allowed us to use alternative definitions for the goal centre width and to extract the corresponding kick distribution between right, centre and left. In addition to analysing the potential right-oriented bias in women goalies’ behaviour, we also analysed the joint distribution of kickers’ and goalkeepers’ choices for each score situation (behind, tied, or ahead). Our findings indicate that the goalkeepers’ general tendency was to dive more often to the right, while the kickers’ tendency was to shoot to the right of the goalies. Moreover, this latter tendency of kicking to the goalies’ right was found to be stronger among the female kickers. Finally, our analysis refutes the claim that goalkeepers exhibit a detrimental right-oriented bias, with this conclusion being even stronger among female goalies.

The Fast Radio Burst-emitting Magnetar SGR 1935+2154—Proper Motion and Variability from Long-term Hubble Space Telescope Monitoring

Abstract We present deep Hubble Space Telescope (HST) near-infrared (NIR) observations of the magnetar SGR 1935+2154 from 2021 June, approximately 6 yr after the first HST observations, a year after the discovery of fast-radio-burst-like emission from the source, and in a period of exceptional high-frequency activity. Although not directly taken during a bursting period the counterpart is a factor of ∼1.5–2.5 brighter than seen at previous epochs with F140W(AB) = 24.65 ± 0.02 mag. We do not detect significant variations of the NIR counterpart within the course of any one orbit (i.e., on minutes to hour timescales), and contemporaneous X-ray observations show SGR 1935+2154 to be at the quiescent level. With a time baseline of 6 yr from the first identification of the counterpart we place stringent limits on the proper motion (PM) of the source, with a measured PM of μ = 3.1 ± 1.5 mas yr−1. The direction of PM indicates an origin of SGR 1935+2154 very close to the geometric center of SNR G57.2+08, further strengthening their association. At an adopted distance of 6.6 ± 0.7 kpc, the corresponding tangential space velocity is ν T = 97 ± 48 km s−1 (corrected for differential Galactic rotation and peculiar solar motion), although its formal statistical determination may be compromised owing to few epochs of observation. The current velocity estimate places it at the low end of the kick distribution for pulsars, and makes it among the lowest known magnetar kicks. When collating the few-magnetar kick constraints available, we find full consistency between the magnetar kick distribution and the much larger pulsar kick sample.

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.

The demographics of neutron star – white dwarf mergers

Context. 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. Aims. Here we characterise the rates of such events, their delay-time distributions, their progenitors, and the distribution of their properties. Methods. We use binary population 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. Results. 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 time-integrated 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 is slightly biased towards earlier times. They typically explode in young 100 Myr < τ < 1 Gyr environments, but have a tail distribution extending to long, gigayear-timescales. Models including significant kicks give rise to relatively wide offset distribution extending to hundreds of kiloparsecs. Conclusions. 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 recently observed rapidly evolving SNe, while they are less likely to be related to the class of Ca-rich SNe.

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.

Neutron star kicks – II. Revision and further testing of the conservation of momentum ‘kick’ model

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.

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.

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.

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.

Gravitational wave recoils in non-axisymmetric Robinson–Trautman spacetimes

We examine the gravitational wave recoil waves and the associated net kick velocities in non-axisymmetric Robinson-Trautman spacetimes. We use characteristic initial data for the dynamics corresponding to non-head-on collisions of black holes. We make a parameter study of the kick distributions, corresponding to an extended range of the incidence angle $\rho_0$ in the initial data. For the range of $\rho_0$ examined ($3^{\circ} \leq \rho_0 \leq 110^{\circ}$) the kick distributions as a function of the symmetric mass parameter $\eta$ satisfy a law obtained from an empirical modification of the Fitchett law, with a parameter $C$ that accounts for the non-zero net gravitational momentum wave fluxes for the equal mass case. The law fits accurately the kick distributions for the range of $\rho_0$ examined, with a rms normalized error of the order of $5 \%$. For the equal mass case the nonzero net gravitational wave momentum flux increases as $\rho_0$ increases, up to $\rho_0 \simeq 55^{\circ}$ beyond which it decreases. The maximum net kick velocity is about $190 {\rm km/s}$ for for the boost parameter considered. For $\rho_0 \geq 50^{\circ}$ the distribution is a monotonous function of $\eta$. The angular patterns of the gravitational waves emitted are examined. Our analysis includes the two polarization modes present in wave zone curvature.

Low-mass X-ray binary populations in galaxy outskirts: Globular clusters and supernova kicks

For the first time, we have systematically explored the population of discrete X-ray sources in the outskirts of early-type galaxies. Based on a broad sample of 20 galaxies observed with Chandra we detected overdensity of X-ray sources in their outskirts. The overdensity appears as halos of resolved sources around the galaxies. These halos are broader than the stellar light, extending out to at least ~ 10 Re (Re is the effective radius). These halos are composed of sources fainter than ~5E38 erg/s, whereas the more luminous sources appear to follow the distribution of the stellar light, suggesting that the excess source population consists of neutron star binaries. Dividing the galaxy sample into four groups according to their stellar mass and specific frequency of globular clusters, we find that the extended halos are present in all groups except for the low-mass galaxies with low globular cluster content. We propose that the extended halos may be comprised of two independent components, low-mass X-ray binaries (LMXBs) located in globular clusters (GCs), which are known to have a wider distribution than the stellar light, and neutron star (NS) LMXBs kicked out of the main body of the parent galaxy by supernova explosions. The available deep optical and X-ray data of NGC 4365 support this conclusion. For this galaxy we identified 60.1+/-10.8 excess sources in the 4-10 Re region of which ~ 40% are located in GCs, whereas ~ 60% are field LMXBs. We interpret the latter as kicked NS LMXBs. We discuss the implications of these results for the natal kick distributions of black holes and neutron stars.