The inspiral and merger of eccentric binaries leads to gravitational waveforms distinct from those generated by circularly merging binaries. Dynamical environments can assemble binaries with high ...eccentricity and peak frequencies within the LIGO band. In this paper, we study binary-single stellar scatterings occurring in dense stellar systems as a source of eccentrically inspiraling binaries. Many interactions between compact binaries and single objects are characterized by chaotic resonances in which the binary-single system undergoes many exchanges before reaching a final state. During these chaotic resonances, a pair of objects has a non-negligible probability of experiencing a very close passage. Significant orbital energy and angular momentum are carried away from the system by gravitational wave (GW) radiation in these close passages, and in some cases this implies an inspiral time shorter than the orbital period of the bound third body. We derive the cross section for such dynamical inspiral outcomes through analytical arguments and through numerical scattering experiments including GW losses. We show that the cross section for dynamical inspirals grows with increasing target binary semi-major axis a and that for equal-mass binaries it scales as a super(2/7). Thus, we expect wide target binaries to predominantly contribute to the production of these relativistic outcomes. We estimate that eccentric inspirals account for approximately 1% of dynamically assembled non-eccentric merging binaries. While these events are rare, we show that binary-single scatterings are a more effective formation channel than single-single captures for the production of eccentrically inspiraling binaries, even given modest binary fractions.
Binary stars are common. While only those with small separations may exchange gas with one another, even the widest binaries interact with their gaseous surroundings. Drag forces and accretion rates ...dictate how these systems are transformed by these interactions. We perform three-dimensional hydrodynamic simulations of Bondi-Hoyle-Lyttleton flows, in which a binary moves supersonically relative to a homogeneous medium, using the adaptive mesh refinement code FLASH. We simulate a range of values of the initial semimajor axis of the orbit relative to the gravitational focusing impact parameter of the pair. When the binary separation is less than the gravitational focusing impact parameter, the pair orbits within a shared bow shock. When the pair is wider, each object has an individual bow shock structure. The long-term evolution of the binary is determined by the timescales for accretion, slowing of the center of mass, and orbital inspiraling. We find a clear hierarchy of these timescales; a binary's center-of-mass motion is slowed over a shorter timescale than the pair inspirals or accretes. In contrast to previous analytic predictions, which assume an unperturbed background medium, we find that the timescale for orbital inspiraling is proportional to the semimajor axis to the 0.19 0.01 power. This positive scaling indicates that gaseous drag forces can drive binaries either to coalescence or to the critical separation at which gravitational radiation dominates their further evolution. We discuss the implications of our results for binaries embedded in the interstellar medium, active galactic nuclei disks, and common envelope phases.
Galactic center black holes appear to be nearly universally surrounded by dense stellar clusters. When these black holes go through an active accretion phase, the multiple components of the accretion ...disk, stellar cluster, and black hole system all coexist. We analyze the effect of drag forces on highly eccentric stellar orbits incurred as stars puncture through the disk plane. Disk crossings dissipate orbital energy, drawing eccentric stars into more circular orbits. For high surface density disks, such as those found around black holes accreting near the Eddington mass accretion limit, the magnitude of this energy dissipation can be larger than the mean scatterings that stars receive by two-body relaxation. One implication of this is the presence of a disk "loss cone" for highly eccentric stellar orbits where the dissipation from disk interaction outweighs scatter via two-body relaxation. The disk loss cone is larger than the tidal disruption loss cone for near-Eddington black hole accretion rates. Stellar orbits within the disk loss cone are lost from the overall cluster as stellar orbits are circularized and stars are potentially ablated by their high-velocity impacts with the disk. We find, however, that the presence of the disk loss cone has a minimal effect on the overall rate of stellar tidal disruptions. Stars are still efficiently fed to the black hole from more-distant stellar orbits that receive large-enough per-orbit scatter to jump over the disk loss cone and end up tidally disrupted.
We present a study of binary-single interactions with energy-loss terms such as tidal dissipation and gravitational-wave (GW) emission added to the equation of motion. The inclusion of such terms ...leads to the formation of compact binaries that form during the three-body interaction through two-body captures. These binaries predominantly merge relatively promptly at high eccentricity, with several observable and dynamical consequences to follow. Despite their possibility for being observed in both present and upcoming transient surveys, their outcomes are not firmly constrained. In this paper, we present an analytical framework that allows to estimate the cross section of such two-body captures, which permits us to study how the corresponding rates depend on the initial orbital parameters, the mass hierarchy, the type of interacting object, and the energy dissipation mechanism. This formalism is applied here to study the formation of two-body GW captures, for which we estimate absolute and relative rates relevant to Advanced LIGO detections. It is shown that two-body GW captures should have compelling observational implications if a sizable fraction of detected compact binaries are formed via dynamical interactions.
ABSTRACT
We analyse how drag forces modify the orbits of objects moving through extended gaseous distributions. We consider how hydrodynamic (surface area) drag forces and dynamical friction ...(gravitational) drag forces drive the evolution of orbital eccentricity. While hydrodynamic drag forces cause eccentric orbits to become more circular, dynamical friction drag can cause orbits to become more eccentric. We develop a semi-analytic model that accurately predicts these changes by comparing the total work and torque applied to the orbit at periapse and apoapse. We use a toy model of a radial power-law density profile, ρ ∝ r−γ, to determine that there is a critical γ = 3 power index, which separates the eccentricity evolution in dynamical friction: orbits become more eccentric for γ < 3 and circularize for γ > 3. We apply these findings to the infall of a Jupiter-like planet into the envelope of its host star. The hydrostatic envelopes of stars are defined by steep density gradients near the limb and shallower gradients in the interior. Under the influence of gaseous dynamical friction, an infalling object’s orbit will first decrease in eccentricity and then increase. The critical separation that delineates these regimes is predicted by the local density slope and is linearly dependent on polytropic index. More broadly, our findings indicate that binary systems may routinely emerge from common envelope phases with non-zero eccentricities that were excited by the dynamical friction forces that drove their orbital tightening.
ABSTRACT When embedded in dense cluster cores, intermediate-mass black holes (IMBHs) acquire close stellar or stellar-remnant companions. These companions are not only gravitationally bound, but also ...tend to hierarchically isolate from other cluster stars through series of multibody encounters. In this paper we study the demographics of IMBH companions in compact star clusters through direct N-body simulations. We study clusters initially composed of 105 or 2 × 105 stars with IMBHs of 75 and 150 solar masses, and we follow their evolution for 6-10 Gyr. A tight, innermost binary pair of IMBH and stellar object rapidly forms. The IMBH has a companion with an orbital semimajor axis at least three times tighter than the second-most-bound object over 90% of the time. These companionships have typical periods on the order of years and are subject to cycles of exchange and destruction. The most frequently observed, long-lived pairings persist for ∼107 years. The demographics of IMBH companions in clusters are diverse: they include both main-sequence, giant stars and stellar remnants. Companion objects may reveal the presence of an IMBH in a cluster in one of several ways. The most-bound companion stars routinely suffer grazing tidal interactions with the IMBH, offering a dynamical mechanism to produce repeated flaring episodes like those seen in the IMBH candidate HLX-1. The stellar winds of companion stars provide a minimum quiescent accretion rate for IMBHs, with implications for radio searches for IMBH accretion in globular clusters. Finally, gravitational wave inspirals of compact objects occur with promising frequency.
Recent observations have revealed that the remnants of stellar-coalescence transients are bipolar. This raises the questions of how these bipolar morphologies arise and what they teach us about the ...mechanisms of mass ejection during stellar mergers and common-envelope phases. In this paper, we analyze hydrodynamic simulations of the lead-in to binary coalescence, a phase of unstable Roche lobe overflow that takes the binary from the Roche limit separation to the engulfment of the more compact accretor within the envelope of the extended donor. As mass transfer runs away at increasing rates, gas trails away from the binary. Contrary to previous expectations, early mass loss from the system remains bound to the binary and forms a circumbinary torus. Later ejecta, generated as the accretor grazes the surface of the donor, have very different morphologies and are unbound. These two components of mass loss from the binary interact as later, higher-velocity ejecta collide with the circumbinary torus formed by earlier mass loss. Unbound ejecta are redirected toward the poles, and escaping material creates a bipolar outflow. Our findings show that the transition from bound to unbound ejecta from coalescing binaries can explain the bipolar nature of their remnants, with implications for our understanding of the origin of bipolar remnants of stellar-coalescence transients and, perhaps, some preplanetary nebulae.
We perform the first systematic study of how dynamical stellar tides and general relativistic (GR) effects affect the dynamics and outcomes of binary-single interactions. For this, we have ...constructed an N-body code that includes tides in the affine approximation, where stars are modeled as self-similar ellipsoidal polytropes, and GR corrections using the commonly used post-Newtonian formalism. Using this numerical formalism, we are able resolve the leading effect from tides and GR across several orders of magnitude in both stellar radius and initial target binary separation. We find that the main effect from tides is the formation of two-body tidal captures that form during the chaotic and resonant evolution of the triple system. The two stars undergoing the capture spiral in and merge. The inclusion of tides can thus lead to an increase in the stellar coalescence rate. We also develop an analytical framework for calculating the cross section of tidal inspirals between any pair of objects with similar mass. From our analytical and numerical estimates, we find that the rate of tidal inspirals relative to collisions increases as the initial semimajor axis of the target binary increases and the radius of the interacting tidal objects decreases. The largest effect is therefore found for triple systems hosting white dwarfs and neutron stars (NSs). In this case, we find the rate of highly eccentric white dwarf-NS mergers to likely be dominated by tidal inspirals. While tidal inspirals occur rarely, we note that they can give rise to a plethora of thermonuclear transients, such as Ca-rich transients.
This paper examines the properties of flows around objects embedded within common envelopes in the simplified context of a "wind tunnel." We establish characteristic relationships between key common ...envelope flow parameters like the Mach number and density scale height. Our wind tunnel is a three-dimensional, Cartesian geometry hydrodynamic simulation setup that includes the gravity of the primary and secondary stars and allows us to study the coefficients of drag and accretion experienced by the embedded object. Accretion and drag lead to a transformation of an embedded object and its orbit during a common envelope phase. We present two suites of simulations spanning a range of density gradients and Mach numbers-relevant for flow near the limb of a stellar envelope to the deep interior. In one suite, we adopt an ideal gas adiabatic exponent of , in the other, . We find that coefficients of drag rise in flows with steeper density gradients and that coefficients of drag and accretion are consistently higher in the more compressible, flow. We illustrate the impact of these newly derived coefficients by integrating the inspiral of a secondary object through the envelopes of ( ) and ( ) giants. In these examples, we find a relatively rapid initial inspiral because, near the stellar limb, dynamical friction drag is generated mainly from dense gas focused from deeper within the primary-star's envelope. This rapid initial inspiral timescale carries potential implications for the timescale of transients from early common envelope interaction.
This paper investigates the recent stellar-merger transient M31LRN 2015 in the Andromeda galaxy. We analyze published optical photometry and spectroscopy along with a Hubble Space Telescope detection ...of the color and magnitude of the pre-outburst source. The transient outburst is consistent with dynamically driven ejecta at the onset of a common envelope episode, which eventually leads to the complete merger of a binary system. The light curve appears to contain two components: first of fast ejecta driven by shocks at the onset of common envelope, and later, ∼0.3 M of further ejecta as the secondary becomes more deeply engulfed within the primary. Just prior to merger, we find that the primary star is a 3-5.5 M subgiant-branch primary star with a radius of 30-40 R . Its position in the color-magnitude diagram shows that it is growing in radius, consistent with a picture where it engulfs its companion. By matching the properties of the primary star to the transient outburst, we show that the optical transient lasts less than 10 orbits of the original binary, which had a pre-merger period of ∼10 days. We consider the possible orbital dynamics leading up to the merger, and argue that if the system merged due to the Darwin tidal instability it implies a lower mass main-sequence companion of 0.1-0.6 M . This analysis represents a promising step toward a more detailed understanding of flows in common envelope episodes through direct observational constraints.