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.
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- NSF-PAR ID:
- 10474312
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 527
- Issue:
- 2
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 2734-2749
- Size(s):
- ["p. 2734-2749"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
ABSTRACT Strong dynamical interactions among stars and compact objects are expected in a variety of astrophysical settings, such as star clusters and the disks of active galactic nuclei. Via a suite of three-dimensional hydrodynamics simulations using the moving-mesh code arepo, we investigate the formation of transient phenomena and their properties in close encounters between an $2\, {\rm M}_{\odot }$ or $20\, {\rm M}_{\odot }$ equal-mass circular binary star and single $20\, {\rm M}_{\odot }$ black hole (BH). Stars can be disrupted by the BH during dynamical interactions, naturally producing electromagnetic transient phenomena. Encounters with impact parameters smaller than the semimajor axis of the initial binary frequently lead to a variety of transients whose electromagnetic signatures are qualitatively different from those of ordinary disruption events involving just two bodies. These include the simultaneous or successive disruptions of both stars and one full disruption of one star accompanied by successive partial disruptions of the other star. On the contrary, when the impact parameter is larger than the semimajor axis of the initial binary, the binary is either simply tidally perturbed or dissociated into bound and unbound single stars (‘micro-Hills’ mechanism). The dissociation of $20\, {\rm M}_{\odot }$ binaries can produce a runaway star and an active BH moving away from one another. Also, the binary dissociation can either produce an interacting binary with the BH, or a non-interacting, hard binary; both could be candidates of BH high- and low-mass X-ray binaries. Hence, our simulations especially confirm that strong encounters can lead to the formation of the (generally difficult to form) BH low-mass X-ray binaries.
-
more » « less
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.
-
Context. The origin of the observed population of Wolf-Rayet (WR) stars in low-metallicity galaxies, such as the Small Magellanic Cloud (SMC), is not yet understood. Standard, single-star evolutionary models predict that WR stars should stem from very massive O-type star progenitors, but these are very rare. On the other hand, binary evolutionary models predict that WR stars could originate from primary stars in close binaries. Aims. We conduct an analysis of the massive O star, AzV 14, to spectroscopically determine its fundamental and stellar wind parameters, which are then used to investigate evolutionary paths from the O-type to the WR stage with stellar evolutionary models. Methods. Multi-epoch UV and optical spectra of AzV 14 are analyzed using the non-local thermodynamic equilibrium (LTE) stellar atmosphere code PoWR. An optical TESS light curve was extracted and analyzed using the PHOEBE code. The obtained parameters are put into an evolutionary context, using the MESA code. Results. AzV 14 is a close binary system with a period of P = 3.7058 ± 0.0013 d. The binary consists of two similar main sequence stars with masses of M 1, 2 ≈ 32 M ⊙ . Both stars have weak stellar winds with mass-loss rates of log Ṁ /( M ⊙ yr −1 ) = −7.7 ± 0.2. Binary evolutionary models can explain the empirically derived stellar and orbital parameters, including the position of the AzV 14 components on the Hertzsprung-Russell diagram, revealing its current age of 3.3 Myr. The model predicts that the primary will evolve into a WR star with T eff ≈ 100 kK, while the secondary, which will accrete significant amounts of mass during the first mass transfer phase, will become a cooler WR star with T eff ≈ 50 kK. Furthermore, WR stars that descend from binary components that have accreted significant amount of mass are predicted to have increased oxygen abundances compared to other WR stars. This model prediction is supported by a spectroscopic analysis of a WR star in the SMC. Conclusions. Inspired by the binary evolutionary models, we hypothesize that the populations of WR stars in low-metallicity galaxies may have bimodal temperature distributions. Hotter WR stars might originate from primary stars, while cooler WR stars are the evolutionary descendants of the secondary stars if they accreted a significant amount of mass. These results may have wide-ranging implications for our understanding of massive star feedback and binary evolution channels at low metallicity.more » « less
-
more » « less
Abstract We present analyses of improved photometric and spectroscopic observations for two detached eclipsing binaries at the turnoff of the open cluster NGC 752: the 1.01 days binary DS And and the 15.53 days BD +37 410. For DS And, we find
M 1= 1.692 ± 0.004 ± 0.010M ⊙,R 1= 2.185 ± 0.004 ± 0.008R ⊙,M 2= 1.184 ± 0.001 ± 0.003M ⊙, andR 2= 1.200 ± 0.003 ± 0.005R ⊙. We either confirm or newly identify unusual characteristics of both stars in the binary: the primary star is found to be slightly hotter than the main-sequence turnoff and there is a more substantial discrepancy in its luminosity compared to models (model luminosities are too large by about 40%), while the secondary star is oversized and cooler compared to other main-sequence stars in the same cluster. The evidence points to nonstandard evolution for both stars, but most plausible paths cannot explain the low luminosity of the primary star. BD +37 410 only has one eclipse per cycle, but extensive spectroscopic observations and the Transiting Exoplanet Survey Satellite light curve constrain the stellar masses well:M 1= 1.717 ± 0.011M ⊙andM 2= 1.175 ± 0.005M ⊙. The radius of the main-sequence primary star near 2.9R ⊙definitively requires large convective core overshooting (>0.2 pressure scale heights) in models for its mass, and multiple lines of evidence point toward an age of 1.61 ± 0.03 ± 0.05 Gyr (statistical and systematic uncertainties). Because NGC 752 is currently undergoing the transition from nondegenerate to degenerate He ignition of its red clump stars, BD +37 410 A directly constrains the star mass where this transition occurs.
