No Arabic abstract
At least $70%$ of massive OBA-type stars reside in binary or higher-order systems. The dynamical evolution of these systems can lend insight into the origins of extreme phenomena such as X-ray binaries and gravitational wave sources. In one such dynamical process, the Eccentric Kozai-Lidov (EKL) Mechanism, a third companion star alters the secular evolution of a binary system. For dynamical stability, these triple systems must have a hierarchical configuration. We explore the effects of a distant third companions gravitational perturbations on a massive binarys orbital configuration before significant stellar evolution has taken place ($leq 10$ Myr). We include tidal dissipation and general relativistic precession. With large ($38,000$ total) Monte-Carlo realizations of massive hierarchical triples, we characterize imprints of the birth conditions on the final orbital distributions. Specifically, we find that the final eccentricity distribution over the range $0.1-0.7$ is an excellent indicator of its birth distribution. Furthermore, we find that the period distributions have a similar mapping for wide orbits. Finally, we demonstrate that the observed period distribution for approximately $10$ Myr-old massive stars is consistent with EKL evolution.
Interacting binaries are of general interest as laboratories for investigating the physics of accretion, which gives rise to the bulk of high-energy radiation in the Galaxy. They allow us to probe stellar evolution processes that cannot be studied in single stars. Understanding the orbital evolution of binaries is essential in order to model the formation of compact binaries. Here we focus our attention on studying orbital evolution driven by angular momentum loss through stellar winds in massive binaries. We run a suite of hydrodynamical simulations of binary stars hosting one mass losing star with varying wind velocity, mass ratio, wind velocity profile and adiabatic index, and compare our results to analytic estimates for drag and angular momentum loss. We find that, at leading order, orbital evolution is determined by the wind velocity and the binary mass ratio. Small ratios of wind to orbital velocities and large accreting companion masses result in high angular momentum loss and a shrinking of the orbit. For wider binaries and binaries hosting lighter mass-capturing companions, the wind mass-loss becomes more symmetric, which results in a widening of the orbit. We present a simple analytic formula that can accurately account for angular momentum losses and changes in the orbit, which depends on the wind velocity and mass ratio. As an example of our formalism, we compare the effects of tides and winds in driving the orbital evolution of high mass X-ray binaries, focusing on Vela X-1 and Cygnus X-1 as examples.
We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrapolation from higher mass ratios, we find that for $qlesssim 10^{-3}$ the eccentricity growth rate becomes negative, i.e., the binary {it circularises} as it shrinks. This behaviour is due to the subset of interacting stars captured in metastable counter-rotating orbits; those stars tend to extract angular momentum from the binary, promoting eccentricity growth for $q>10^{-3}$, but tend to inject angular momentum into the binary driving it towards circularisation for $q<10^{-3}$. The physical origin of this behaviour requires a detailed study of the orbits of this subset of stars and is currently under investigation. Our findings might have important consequences for intermediate MBHs (IMBHs) inspiralling onto MBHs (e.g. a putative $10^3rm M_{odot}$ black hole inspiralling onto SgrA$^*$).
The number of binaries containing black holes or neutron stars depends critically on the fraction of binaries that survive supernova explosions. We searched for surviving star plus remnant binaries in a sample of 49 supernova remnants (SNR) containing 23 previously identified compact remnants and three high mass X-ray binaries (HMXB), finding no new interacting or non-interacting binaries. The upper limits on any main sequence stellar companion are typically <0.2Msun and are at worst <3Msun. This implies that f<0.1 of core collapse SNRs contain a non-interacting binary, and f=0.083 (0.032<f<0.17) contain an interacting binary at 90% confidence. We also find that the transverse velocities of HMXBs are low, with a median of only 12~km/s for field HMXBs, so surviving binaries will generally be found very close to the explosion center. We compare the results to a standard StarTrack binary population synthesis (BPS) model, finding reasonable agreement with the observations. In particular, the BPS models predict that 5% of SNe should leave a star plus remnant binary.
Dense stellar clusters are natural sites for the origin and evolution of exotic objects such as relativistic binaries (potential gravitational wave sources), blue stragglers, etc. We investigate the secular dynamics of a binary system driven by the global tidal field of an axisymmetric stellar cluster in which the binary orbits. In a companion paper (Hamilton & Rafikov 2019a) we developed a general Hamiltonian framework describing such systems. The effective (doubly-averaged) Hamiltonian derived there encapsulates all information about the tidal potential experienced by the binary in its orbit around the cluster in a single parameter $Gamma$. Here we provide a thorough exploration of the phase-space of the corresponding secular problem as $Gamma$ is varied. We find that for $Gamma > 1/5$ the phase-space structure and the evolution of binary orbital element are qualitatively similar to the Lidov-Kozai problem. However, this is only one of four possible regimes, because the dynamics are qualitatively changed by bifurcations at $Gamma = 1/5,0,-1/5$. We show how the dynamics are altered in each regime and calculate characteristics such as secular evolution timescale, maximum possible eccentricity, etc. We verify the predictions of our doubly-averaged formalism numerically and find it to be very accurate when its underlying assumptions are fulfilled, typically meaning that the secular timescale should exceed the period of the binary around the cluster by $gtrsim 10-10^2$ (depending on the cluster potential and binary orbit). Our results may be relevant for understanding the nature of a variety of exotic systems harboured by stellar clusters.
The nature of the progenitors and explosion mechanism of Type Iax supernovae (SNe Iax) remain a mystery. The single-degenerate (SD) systems that involve the incomplete pure deflagration explosions of near-Chandrasekhar-mass white dwarfs (WDs) have recently been proposed for producing SNe Iax, in which non-degenerate companions are expected to survive from SN explosions. In this work we concentrate on the main-sequence (MS) donor SD progenitor systems. By mapping the computed companion models from three-dimensional hydrodynamical simulations of ejecta-companion interaction into a one-dimensional stellar evolution code MESA, we investigate the long-term appearance and observational signatures of surviving MS companions of SNe Iax by tracing their post-impact evolution. Depending on different MS companion models, it is found that the shocked surviving companion stars can significantly expand and evolve to be more luminous (5-500 Lsun) for a time-scale of 10-1e4 yr. Comparing with the late-time light curve of an observed SN Iax (SN 2005hk), it is suggested that surviving MS companions of SNe Iax would expect to be visible about 1000 days after the explosion when SN itself has been faded.