Dynamical constraints on the origin of the young B-stars in the Galactic center

Regular star formation is thought to be inhibited close to the massive black hole (MBH) in the Galactic center. Nevertheless, tens of young main sequence B stars have been observed in an isotropic distribution close to it. Various models have been suggested for the formation of the B-stars closest to the MBH (<0.05 pc; the S-stars), typically involving the migration of these stars from their original birthplace to their currently observed position. Here we explore the orbital phase space distribution of the B-stars throughout the central pc expected from the various suggested models for the origin of the B-stars. We find that most of these models have difficulties in explaining, by themselves, both the population of the S-stars (<0.05 pc), and the population of the young B-stars further away (up to 0.5 pc). Most models grossly over-predict the number of B-stars up to 0.5 pc, given the observed number of S-stars. Such models include the intermediate-mass black hole assisted cluster inspiral scenario, Kozai-like perturbations by two disks, spiral density waves migration in a gaseous disk, and some of the eccentric disk instability models. We focus on one of the other models, the massive perturber induced binary disruption, which is consistent with both the S-stars and the extended population of B-stars further away. For this model we use analytical arguments and N-body simulations to provide further observational predictions. These could be compared with future observations to further support this model, constrain it or refute it. These predictions include the radial distribution of the young B-stars, their eccentricity distribution and its dependence on distance from the MBH (higher eccentricities at larger distances from the MBH), as well as less specific expectations regarding their mass function.Comment: Comments are welcome

The properties of dynamically ejected runaway and hyper-runaway stars

Runaway stars are stars observed to have large peculiar velocities. Two mechanisms are thought to contribute to the ejection of runaway stars, both involve binarity (or higher multiplicity). In the binary supernova scenario a runaway star receives its velocity when its binary massive companion explodes as a supernova (SN). In the alternative dynamical ejection scenario, runaway stars are formed through gravitational interactions between stars and binaries in dense, compact clusters or cluster cores. Here we study the ejection scenario. We make use of extensive N-body simulations of massive clusters, as well as analytic arguments, in order to to characterize the expected ejection velocity distribution of runaways stars. We find the ejection velocity distribution of the fastest runaways (>~80 km s^-1) depends on the binary distribution in the cluster, consistent with our analytic toy model, whereas the distribution of lower velocity runaways appears independent of the binaries properties. For a realistic log constant distribution of binary separations, we find the velocity distribution to follow a simple power law; Gamma(v) goes like v^(-8/3) for the high velocity runaways and v^(-3/2) for the low velocity ones. We calculate the total expected ejection rates of runaway stars from our simulated massive clusters and explore their mass function and their binarity. The mass function of runaway stars is biased towards high masses, and depends strongly on their velocity. The binarity of runaways is a decreasing function of their ejection velocity, with no binaries expected to be ejected with v>150 km s^-1. We also find that hyper-runaways with velocities of hundreds of km s^-1 can be dynamically ejected from stellar clusters, but only at very low rates, which cannot account for a significant fraction of the observed population of hyper-velocity stars in the Galactic halo.Comment: Now matching published ApJ versio

Dynamical evolution of the young stars in the Galactic center: N-body simulations of the S-stars

We use N-body simulations to study the evolution of the orbital eccentricities of stars deposited near (<0.05 pc) the Milky Way massive black hole (MBH), starting from initial conditions motivated by two competing models for their origin: formation in a disk followed by inward migration; and exchange interactions involving a binary star. The first model predicts modest eccentricities, lower than those observed in the S-star cluster, while the second model predicts higher eccentricities than observed. The N-body simulations include a dense cluster of 10 M_sun stellar black holes (SBHs), expected to accumulate near the MBH by mass segregation. Perturbations from the SBHs tend to randomize the stellar orbits, partially erasing the dynamical signatures of their origin. The eccentricities of the initially highly eccentric stars evolve, in 20 Myr (the S-star lifespan), to a distribution that is consistent at the ~95 % level with the observed eccentricity distribution. In contrast, the eccentricities of the initially more circular orbits fail to evolve to the observed values in 20 Myr, arguing against the disk migration scenario. We find that 20 % - 30 % of the S-stars are tidally disrupted by the MBH over their lifetimes, and that the S-stars are not likely to be ejected as hypervelocity stars outside the central 0.05 pc by close encounters with stellar black holes.Comment: 6 pages, 2 figures. Minor corrections, Sumitted to Ap

Wind-shearing in gaseous protoplanetary disks and the evolution of binary planetesimals

One of the first stages of planet formation is the growth of small planetesimals. This early stage occurs much before the dispersal of most of the gas from the protoplanetary disk. Due to their different aerodynamic properties, planetesimals of different sizes and shapes experience different drag forces from the gas during this time. Such differential forces produce a wind-shearing (WISH) effect between close by, different size planetesimals. For any two planetesimals, a WISH radius can be considered, at which the differential acceleration due to the wind becomes greater than the mutual gravitational pull between the planetesimals. We find that the WISH radius could be much smaller than the Hill radius, i.e. WISH could play a more important role than tidal perturbations by the star. Here we study the WISH radii for planetesimal pairs of different sizes and compare the effects of wind and gravitational shearing (drag force vs. gravitational tidal force). We then discuss the role of WISH for the stability and survival of binary planetesimals. Binaries are sheared apart by the wind if they are wider than their WISH radius. WISH-stable binaries can inspiral and possibly coalesce due to gas drag. Here, we calculate the WISH radius and the gas drag-induced merger timescale, providing stability and survival criteria for gas-embedded binary planetesimals. Our results suggest that even WISH-stable binaries may merge in times shorter than the lifetime of the gaseous disk. This may constrain currently observed binary planetesimals to have formed far from the star or at a late stage after the dispersal of most of the disk gas. We note that the WISH radius may also be important for other processes such as planetesimal erosion and planetesimal encounters and collisions in a gaseous environment.Comment: ApJ, in pres