skip to main content


Title: The dominant mechanism(s) for populating the outskirts of star clusters with neutron star binaries
ABSTRACT

It has been argued that heavy binaries composed of neutron stars (NSs) and millisecond pulsars (MSPs) can end up in the outskirts of star clusters via an interaction with a massive black hole (BH) binary expelling them from the core. We argue here, however, that this mechanism will rarely account for such observed objects. Only for primary masses ≲100 M⊙ and a narrow range of orbital separations should a BH–BH binary be both dynamically hard and produce a sufficiently low recoil velocity to retain the NS binary in the cluster. Hence, BH binaries are in general likely to eject NSs from clusters. We explore several alternative mechanisms that would cause NS/MSP binaries to be observed in the outskirts of their host clusters after a Hubble time. The most likely mechanism is a three-body interaction involving the NS/MSP binary and a normal star. We compare to Monte Carlo simulations of cluster evolution for the globular clusters NGC 6752 and 47 Tuc, and show that the models not only confirm that normal three-body interactions involving all stellar-mass objects are the dominant mechanism for putting NS/MSP binaries into the cluster outskirts, but also reproduce the observed NS/MSP binary radial distributions without needing to invoke the presence of a massive BH binary. Higher central densities and an episode of core collapse can broaden the radial distributions of NSs/MSPs and NS/MSP binaries due to three-body interactions, making these clusters more likely to host NSs in the cluster outskirts.

 
more » « less
NSF-PAR ID:
10478649
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
527
Issue:
3
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 6913-6925
Size(s):
["p. 6913-6925"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Observations of X-ray binaries indicate a dearth of compact objects in the mass range from ∼2 − 5  M ⊙ . The existence of this (first mass) gap has been used to discriminate between proposed engines behind core-collapse supernovae. From LIGO/Virgo observations of binary compact remnant masses, several candidate first mass gap objects, either neutron stars (NSs) or black holes (BHs), were identified during the O3 science run. Motivated by these new observations, we study the formation of BH-NS mergers in the framework of isolated classical binary evolution, using population synthesis methods to evolve large populations of binary stars (Population I and II) across cosmic time. We present results on the NS to BH mass ratios ( q  =  M NS / M BH ) in merging systems, showing that although systems with a mass ratio as low as q  = 0.02 can exist, typically BH-NS systems form with moderate mass ratios q  = 0.1 − 0.2. If we adopt a delayed supernova engine, we conclude that ∼30% of BH-NS mergers may host at least one compact object in the first mass gap (FMG). Even allowing for uncertainties in the processes behind compact object formation, we expect the fraction of BH-NS systems ejecting mass during the merger to be small (from ∼0.6 − 9%). In our reference model, we assume: (i) the formation of compact objects within the FMG, (ii) natal NS/BH kicks decreased by fallback, (iii) low BH spins due to Tayler-Spruit angular momentum transport in massive stars. We find that ≲1% of BH-NS mergers will have any mass ejection and about the same percentage will produce kilonova bright enough to have a chance of being detected with a large (Subaru-class) 8 m telescope. Interestingly, all these mergers will have both a BH and an NS in the FMG. 
    more » « less
  2. ABSTRACT

    We present multi-epoch spectroscopic follow-up of a sample of ellipsoidal variables selected from Gaia Data Release 3 (DR3) as candidates for hosting quiescent black holes (BHs) and neutron stars (NSs). Our targets were identified as BH/NS candidates because their optical light curves – when interpreted with models that attribute variability to tidal distortion of a star by a companion that contributes negligible light – suggest that the companions are compact objects. From the likely BH/NS candidates identified in recent work accompanying Gaia DR3, we select 14 of the most promising targets for follow-up. We obtained spectra for each object at 2–10 epochs, strategically observing near conjunction to best constrain the radial velocity semi-amplitude. From the measured semi-amplitudes of the radial velocity curves, we derive minimum companion masses of $M_{2,\, \rm min} \le 0.5 \, {\rm M}_{\odot }$ in all cases. Assuming random inclinations, the typical inferred companion mass is $M_2 \sim 0.15\, {\rm M}_{\odot }$. This makes it unlikely that any of these systems contain a BH or NS, and we consider alternative explanations for the observed variability. We can best reproduce the observed light curves and radial velocities with models for unequal-mass contact binaries with star-spots. Some of the objects in our sample may also be detached main-sequence binaries, or even single stars with pulsations or star-spot variability masquerading as ellipsoidal variation. We provide recommendations for future spectroscopic efforts to further characterize this sample and more generally to search for compact object companions in close binaries.

     
    more » « less
  3. ABSTRACT

    Dynamical interactions involving binaries play a crucial role in the evolution of star clusters and galaxies. We continue our investigation of the hydrodynamics of three-body encounters, focusing on binary black hole (BBH) formation, stellar disruption, and electromagnetic (EM) emission in dynamical interactions between a BH-star binary and a stellar-mass BH, using the moving-mesh hydrodynamics code AREPO. This type of encounters can be divided into two classes depending on whether the final outcome includes BBHs. This outcome is primarily determined by which two objects meet at the first closest approach. BBHs are more likely to form when the star and the incoming BH encounter first with an impact parameter smaller than the binary’s semimajor axis. In this case, the star is frequently disrupted. On the other hand, when the two BHs encounter first, frequent consequences are an orbit perturbation of the original binary or a binary member exchange. For the parameters chosen in this study, BBH formation, accompanied by stellar disruption, happens in roughly one out of four encounters. The close correlation between BBH formation and stellar disruption has possible implications for EM counterparts at the binary’s merger. The BH that disrupts the star is promptly surrounded by an optically and geometrically thick disc with accretion rates exceeding the Eddington limit. If the debris disc cools fast enough to become long-lived, EM counterparts can be produced at the time of the BBH merger.

     
    more » « less
  4. ABSTRACT

    Multibody dynamical interactions of binaries with other objects are one of the main driving mechanisms for the evolution of star clusters. It is thus important to bring our understanding of three-body interactions beyond the commonly employed point-particle approximation. To this end, we here investigate the hydrodynamics of three-body encounters between star–black hole (BH) binaries and single stars, focusing on the identification of final outcomes and their long-term evolution and observational properties, using the moving-mesh hydrodynamics code AREPO. This type of encounter produces five types of outcomes: stellar disruption, stellar collision, weak perturbation of the original binary, binary member exchange, and triple formation. The two decisive parameters are the binary phase angle, which determines which two objects meet at the first closest approach, and the impact parameter, which sets the boundary between violent and non-violent interactions. When the impact parameter is smaller than the semimajor axis of the binary, tidal disruptions and star-BH collisions frequently occur when the BH and the incoming star first meet, while the two stars mostly merge when the two stars meet first instead. In both cases, the BHs accrete from an accretion disc at super-Eddington rates, possibly generating flares luminous enough to be observed. The stellar collision products either form a binary with the BH or remain unbound to the BH. Upon collision, the merged stars are hotter and larger than the main sequence stars of the same mass at similar age. Even after recovering their thermal equilibrium state, stellar collision products, if isolated, would remain hotter and brighter than main sequence stars until becoming giants.

     
    more » « less
  5. Abstract

    Gravitational-wave (GW) detections of merging neutron star–black hole (NSBH) systems probe astrophysical neutron star (NS) and black hole (BH) mass distributions, especially at the transition between NS and BH masses. Of particular interest are the maximum NS mass, minimum BH mass, and potential mass gap between them. While previous GW population analyses assumed all NSs obey the same maximum mass, if rapidly spinning NSs exist, they can extend to larger maximum masses than nonspinning NSs. In fact, several authors have proposed that the ∼2.6Mobject in the event GW190814—either the most massive NS or least massive BH observed to date—is a rapidly spinning NS. We therefore infer the NSBH mass distribution jointly with the NS spin distribution, modeling the NS maximum mass as a function of spin. Using four LIGO–Virgo NSBH events including GW190814, if we assume that the NS spin distribution is uniformly distributed up to the maximum (breakup) spin, we infer the maximum nonspinning NS mass is2.70.4+0.5M(90% credibility), while assuming only nonspinning NSs, the NS maximum mass must be >2.53M(90% credibility). The data support the mass gap’s existence, with a minimum BH mass at5.41.0+0.7M. With future observations, under simplified assumptions, 150 NSBH events may constrain the maximum nonspinning NS mass to ±0.02M, and we may even measure the relation between the NS spin and maximum mass entirely from GW data. If rapidly rotating NSs exist, their spins and masses must be modeled simultaneously to avoid biasing the NS maximum mass.

     
    more » « less