Source author record

Frederic A. Rasio

Frederic A. Rasio appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.HE astro-ph.SR astro-ph.EP astro-ph gr-qc astro-ph.IM astro-ph.CO physics.comp-ph

Catalog footprint

What is connected

70works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Compact Object Modeling in the Globular Cluster 47 Tucanae

The globular cluster 47~Tucanae (47~Tuc) is one of the most massive star clusters in the Milky Way and is exceptionally rich in exotic stellar populations. For several decades it has been a favorite target of observers, and yet it is computationally very challenging to model because of its large number of stars ($N\gtrsim 10^6$) and high density. Here we present detailed and self-consistent 47~Tuc models computed with the \texttt{Cluster Monte Carlo} code (\texttt{CMC}). The models include all relevant dynamical interactions coupled to stellar and binary evolution, and reproduce various observations, including the surface brightness and velocity dispersion profiles, pulsar accelerations, and numbers of compact objects. We show that the present properties of 47~Tuc are best reproduced by adopting an initial stellar mass function that is both bottom-heavy and top-light relative to standard assumptions \citep[as in, e.g.,][]{Kroupa2001}, and an initial Elson profile \citep{Elson1987} that is overfilling the cluster's tidal radius. We include new prescriptions in \texttt{CMC} for the formation of binaries through giant star collisions and tidal captures, and we show that these mechanisms play a crucial role in the formation of neutron star binaries and millisecond pulsars in 47~Tuc; our best-fit model contains $\sim 50$ millisecond pulsars, $70\%$ of which are formed through giant collisions and tidal captures. Our models also suggest that 47~Tuc presently contains up to $\sim 200$ stellar-mass black holes, $\sim 5$ binary black holes, $\sim 15$ low-mass X-ray binaries, and $\sim 300$ cataclysmic variables.

preprint2022arXiv

Eccentric Mergers of Intermediate-Mass Black Holes from Evection Resonances in AGN Disks

We apply the theory of nonlinear resonance capture to the problem of a black hole binary (BHB) orbiting a supermassive black hole (SMBH) while embedded in the accretion disk of an active galactic nucleus (AGN). If successful, resonance capture can trigger dramatic growth in the BHB eccentricity, with important consequences for the BHB merger timescale, as well as for the gravitational wave (GW) signature of such an eccentric merger. This resonance capture may occur when the orbital period around the SMBH (the "outer binary") and the apsidal precession of the BHB (the "inner binary") are in a 1:1 commensurability. This effect is analogous to the phenomenon of lunar evection resonance in the early Sun-Earth-Moon system, with the distinction that in the present case, the BHB apsidal precession is due to general relativity, rather than rotationally-induced distortion. In contrast to the case of lunar evection, however, the inner binary undergoes orbital decay driven by GW emission, rather than orbital expansion driven by tidal dissipation. This distinction fundamentally alters the three-body dynamics, forbidding resonance capture, and limiting eccentricity growth. However, if the BHB migrates through of a gaseous AGN disk, the change in the outer binary can counterbalance the suppressing effect of BHB decay, permitting evection resonance capture and the production of eccentric BHB mergers. We compute the likelihood of resonance capture assuming an agnostic distribution of parameters for the three bodies involved and for the properties of the AGN disk. We find that intermediate-mass ratio BHBs (involving an intermediate-mass black hole and a stellar-mass black hole) are the most likely to be captured into evection resonance and thus undergo an eccentric merger. We also compute the GW signature of these mergers, showing that they can enter the LISA band while eccentric.

preprint2022arXiv

Formation of Low-mass Black Holes and Single Millisecond Pulsars in Globular Clusters

Close encounters between neutron stars and main-sequence stars occur in globular clusters and may lead to various outcomes. Here we study encounters resulting in tidal disruption of the star. Using $N$-body models, we predict the typical stellar masses in these disruptions and the dependence of the event rate on host cluster properties. We find that tidal disruption events occur most frequently in core-collapsed globular clusters and that roughly $25\%$ of the disrupted stars are merger products (i.e., blue straggler stars). Using hydrodynamic simulations, we model the tidal disruptions themselves (over timescales of days) to determine the mass bound to the neutron star and the properties of the accretion disks formed. In general, we find that roughly $80-90\%$ of the initial stellar mass becomes bound to the neutron star following disruption. Additionally, we find that neutron stars receive impulsive kicks of up to about $20\,$km/s as a result of the asymmetry of unbound ejecta; these kicks place these neutron stars on elongated orbits within their host cluster, with apocenter distances well outside the cluster core. Finally, we model the evolution of the (hypercritical) accretion disks on longer timescales (days to years after disruption) to estimate the accretion rate onto the neutron stars and accompanying spin-up. As long as $\gtrsim1\%$ of the bound mass accretes onto the neutron star, millisecond spin periods can be attained. We argue the growing numbers of isolated millisecond pulsars observed in globular clusters may have formed, at least in part, through this mechanism. In the case of significant mass growth, some of these neutron stars may collapse to form low-mass ($\lesssim3\,M_{\odot}$) black holes.

preprint2022arXiv

Gravitational Microlensing Rates in Milky Way Globular Clusters

Many recent observational and theoretical studies suggest that globular clusters (GCs) host compact object populations large enough to play dominant roles in their overall dynamical evolution. Yet direct detection, particularly of black holes and neutron stars, remains rare and limited to special cases, such as when these objects reside in close binaries with bright companions. Here we examine the potential of microlensing detections to further constrain these dark populations. Based on state-of-the-art GC models from the CMC Cluster Catalog, we estimate the microlensing event rates for black holes, neutron stars, white dwarfs, and, for comparison, also for M dwarfs in Milky Way GCs, as well as the effects of different initial conditions on these rates. Among compact objects, we find that white dwarfs dominate the microlensing rates, simply because they largely dominate by numbers. We show that microlensing detections are in general more likely in GCs with higher initial densities, especially in clusters that undergo core collapse. We also estimate microlensing rates in the specific cases of M22 and 47 Tuc using our best-fitting models for these GCs. Because their positions on the sky lie near the rich stellar backgrounds of the Galactic bulge and the Small Magellanic Cloud, respectively, these clusters are among the Galactic GCs best-suited for dedicated microlensing surveys. The upcoming 10-year survey with the Rubin Observatory may be ideal for detecting lensing events in GCs.

preprint2022arXiv

Hydrodynamics of Collisions and Close Encounters between Stellar Black Holes and Main-sequence Stars

Recent analyses have shown that close encounters between stars and stellar black holes occur frequently in dense star clusters. Depending upon the distance at closest approach, these interactions can lead to dissipating encounters such as tidal captures and disruptions, or direct physical collisions, all of which may be accompanied by bright electromagnetic transients. In this study, we perform a wide range of hydrodynamic simulations of close encounters between black holes and main-sequence stars that collectively cover the parameter space of interest, and we identify and classify the various possible outcomes. In the case of nearly head-on collisions, the star is completely disrupted with roughly half of the stellar material becoming bound to the black hole. For more distant encounters near the classical tidal-disruption radius, the star is only partially disrupted on the first pericenter passage. Depending upon the interaction details, the partially disrupted stellar remnant may be tidally captured by the black hole or become unbound (in some cases, receiving a sufficiently large impulsive kick from asymmetric mass loss to be ejected from its host cluster). In the former case, the star will undergo additional pericenter passages before ultimately being disrupted fully. Based on the properties of the material bound to the black hole at the end of our simulations (in particular, the total bound mass and angular momentum), we comment upon the expected accretion process and associated electromagnetic signatures that are likely to result.

preprint2022arXiv

Merger rates of intermediate-mass black hole binaries in nuclear star clusters

Repeated mergers of stellar-mass black holes (BHs) in dense star clusters can produce intermediate-mass black holes (IMBHs). In particular, nuclear star clusters at the centers of galaxies have deep enough potential wells to retain most of the BH merger products, in spite of the significant recoil kicks due to anisotropic emission of gravitational radiation. These events can be detected in gravitational waves (GWs), which represent an unprecedented opportunity to reveal IMBHs. In this paper, we analyze the statistical results of a wide range of numerical simulations, which encompass different cluster metallicities, initial BH seed masses, and initial BH spins, and we compute the merger rate of IMBH binaries. We find that merger rates are in the range $0.01$-$10$\,Gpc$^{-3}$\,yr$^{-1}$ depending on IMBH masses. We also compute the number of multi-band detections in ground-based and space-based observatories. Our model predicts that a few merger events per year should be detectable with LISA, DECIGO, ET, and LIGO for IMBHs with masses $\lesssim 1000\msun$, and a few tens of merger events per year with DECIGO, ET, and LIGO only.

preprint2022arXiv

Repeated mergers, mass-gap black holes, and formation of intermediate-mass black holes in nuclear star clusters

Current theoretical models predict a mass gap with a dearth of stellar black holes (BHs) between roughly $50\,M_\odot$ and $100\,M_\odot$, while, above the range accessible through massive star evolution, intermediate-mass BHs (IMBHs) still remain elusive. Repeated mergers of binary BHs, detectable via gravitational wave emission with the current LIGO/Virgo/Kagra interferometers and future detectors such as LISA or the Einstein Telescope, can form both mass-gap BHs and IMBHs. Here we explore the possibility that mass-gap BHs and IMBHs are born as a result of successive BH mergers in dense star clusters. In particular, nuclear star clusters at the centers of galaxies have deep enough potential wells to retain most of the BH merger products after they receive significant recoil kicks due to anisotropic emission of gravitational radiation. We show that a massive stellar BH seed can easily grow to $\sim 10^3 - 10^4\,M_\odot$ as a result of repeated mergers with other smaller BHs. We find that lowering the cluster metallicity leads to larger final BH masses. We also show that the growing BH spin tends to decrease in magnitude with the number of mergers, so that a negative correlation exists between final mass and spin of the resulting IMBHs. Assumptions about the birth spins of stellar BHs affect our results significantly, with low birth spins leading to the production of a larger population of massive BHs.

preprint2021arXiv

Intermediate-mass Black Holes from High Massive-star Binary Fractions in Young Star Clusters

Black holes formed in dense star clusters, where dynamical interactions are frequent, may have fundamentally different properties than those formed through isolated stellar evolution. Theoretical models for single star evolution predict a gap in the black hole mass spectrum from roughly $40-120\,M_{\odot}$ caused by (pulsational) pair-instability supernovae. Motivated by the recent LIGO/Virgo event GW190521, we investigate whether black holes with masses within or in excess of this "upper-mass gap" can be formed dynamically in young star clusters through strong interactions of massive stars in binaries. We perform a set of $N$-body simulations using the CMC cluster-dynamics code to study the effects of the high-mass binary fraction on the formation and collision histories of the most massive stars and their remnants. We find that typical young star clusters with low metallicities and high binary fractions in massive stars can form several black holes in the upper-mass gap and often form at least one intermediate-mass black hole. These results provide strong evidence that dynamical interactions in young star clusters naturally lead to the formation of more massive black hole remnants.

preprint2021arXiv

Modeling Dense Star Clusters in the Milky Way and Beyond with the Cluster Monte Carlo Code

We describe the public release of the Cluster Monte Carlo Code (CMC) a parallel, star-by-star $N$-body code for modeling dense star clusters. CMC treats collisional stellar dynamics using Hénon's method, where the cumulative effect of many two-body encounters is statistically reproduced as a single effective encounter between nearest-neighbor particles on a relaxation timescale. The star-by-star approach allows for the inclusion of additional physics, including strong gravitational three- and four-body encounters, two-body tidal and gravitational-wave captures, mass loss in arbitrary galactic tidal fields, and stellar evolution for both single and binary stars. The public release of CMC is pinned directly to the COSMIC population synthesis code, allowing dynamical star cluster simulations and population synthesis studies to be performed using identical assumptions about the stellar physics and initial conditions. As a demonstration, we present two examples of star cluster modeling: first, we perform the largest ($N = 10^8$) star-by-star $N$-body simulation of a Plummer sphere evolving to core collapse, reproducing the expected self-similar density profile over more than 15 orders of magnitude; second, we generate realistic models for typical globular clusters, and we show that their dynamical evolution can produce significant numbers of black hole mergers with masses greater than those produced from isolated binary evolution (such as GW190521, a recently reported merger with component masses in the pulsational pair-instability mass gap).

preprint2021arXiv

The Observed Rate of Binary Black Hole Mergers can be Entirely Explained by Globular Clusters

Since the first signal in 2015, the gravitational-wave detections of merging binary black holes (BBHs) by the LIGO and Virgo collaborations (LVC) have completely transformed our understanding of the lives and deaths of compact object binaries, and have motivated an enormous amount of theoretical work on the astrophysical origin of these objects. We show that the phenomenological fit to the redshift-dependent merger rate of BBHs from Abbott et al. (2020) is consistent with a purely dynamical origin for these objects, and that the current merger rate of BBHs from the LVC could be explained entirely with globular clusters alone. While this does not prove that globular clusters are the dominant formation channel, we emphasize that many formation scenarios could contribute a significant fraction of the current LVC rate, and that any analysis that assumes a single (or dominant) mechanism for producing BBH mergers is implicitly using a specious astrophysical prior.

preprint2020arXiv

A Dynamical Survey of Stellar-Mass Black Holes in 50 Milky Way Globular Clusters

Recent numerical simulations of globular clusters (GCs) have shown that stellar-mass black holes (BHs) play a fundamental role in driving cluster evolution and shaping their present-day structure. Rapidly mass-segregating to the center of GCs, BHs act as a dynamical energy source via repeated super-elastic scattering, delaying onset of core collapse and limiting mass segregation for visible stars. While recent discoveries of BH candidates in Galactic and extragalactic GCs have further piqued interest in BH-mediated cluster dynamics, numerical models show that even if significant BH populations remain in today's GCs, they are typically in configurations that are not directly detectable. We demonstrated in Weatherford et al. (2018) that an anti-correlation between a suitable measure of mass segregation ($Δ$) in observable stellar populations and the number of retained BHs in GC models can be applied to indirectly probe BH populations in real GCs. Here, we estimate the number and total mass of BHs presently retained in 50 Milky Way GCs from the ACS Globular Cluster Survey by measuring $Δ$ between populations of main sequence stars, using correlations found between $Δ$ and BH retention in the CMC Cluster Catalog models. We demonstrate that the range in $Δ$'s distribution from our models matches that for observed GCs to a remarkable degree. Our results further provide the narrowest constraints to-date on the retained BH populations in the GCs analyzed. Of these 50 GCs, we identify NGCs 2808, 5927, 5986, 6101, and 6205 to presently contain especially large BH populations, each with total BH mass exceeding $10^3\,\rm{M_{\odot}}$.

preprint2020arXiv

COSMIC Variance in Binary Population Synthesis

The formation and evolution of binary stars is a critical component of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling and/or merging compact binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observations of electromagnetic and gravitational wave sources. Binary population synthesis is a tool that balances physical modeling with simulation speed to produce large binary populations on timescales of days. We present a community-developed binary population synthesis suite: COSMIC which is designed to simulate compact-object binary populations and their progenitors. As a proof of concept, we simulate the Galactic population of compact binaries and their gravitational wave signal observable by the Laser Interferometer Space Antenna (LISA). We find that $\sim10^8$ compact binaries reside in the Milky Way today, while $\sim10^4$ of them may be resolvable by LISA.

preprint2020arXiv

Demographics of triple systems in dense star clusters

Depending on the stellar type, more than $\sim 50$\% and $\sim 15\%$ of stars in the field have at least one and two stellar companions, respectively. Hierarchical systems can be assembled dynamically in dense star clusters, as a result of few-body encounters among stars and/or compact remnants in the cluster core. In this paper, we present the demographics of stellar and compact-object triples formed via binary--binary encounters in the \texttt{CMC Cluster Catalog}, a suite of cluster simulations with present-day properties representative of the globular clusters (GCs) observed in the Milky Way. We show how the initial properties of the host cluster set the typical orbital parameters and formation times of the formed triples. We find that a cluster typically assembles hundreds of triples with at least one black hole (BH) in the inner binary, while only clusters with sufficiently small virial radii are efficient in producing triples with no BHs, as a result of the BH-burning process. We show that a typical GC is expected to host tens of triples with at least one luminous component at present day. We discuss how the Lidov-Kozai mechanism can drive the inner binary of the formed triples to high eccentricities, whenever it takes place before the triple is dynamically reprocessed by encountering another cluster member. Some of these systems can reach sufficiently large eccentricities to form a variety of transients and sources, such as blue stragglers, X-ray binaries, Type Ia Supernovae, Thorne-Zytkow objects, and LIGO/Virgo sources.

preprint2020arXiv

Gravitational-wave captures by intermediate-mass black holes in galactic nuclei

Intermediate-mass black holes (IMBHs) have not been detected beyond any reasonable doubt, despite their potential role as massive seeds for quasars and sources of tidal disruption events, ultra-luminous X-ray sources, dwarf galaxy feedback, and hypervelocity stars. Gravitational wave (GW) observations can help to find and confirm the existence of IMBHs. Current and upcoming detectors, such as LIGO, Virgo, KAGRA, LISA, ET, and DECIGO promise to identify the full range from stellar-mass to supermassive black holes (SMBHs). In this paper, we address the question of whether IMBHs can produce GWs in galactic nuclei. We consider the possibility that stellar black holes (SBHs) form bound systems and later coalesce with an IMBH through gravitational captures in the dense nucleus. We show that this mechanism is efficient for IMBH masses in the range $\sim 3\times 10^3$ M$_\odot$--$2\times 10^4$ M$_\odot$. We find that the typical distributions of peak frequencies and merger timescales depend mainly on the IMBH mass. In particular, the typical peak frequency is about $0.2\,$Hz, $0.1\,$Hz, $0.09\,$Hz, and $0.05\,$Hz for $M_{\rm IMBH}=5\times 10^3$ M$_\odot$, $8\times 10^3$ M$_\odot$, $1\times 10^4$ M$_\odot$, and $2\times 10^4$ M$_\odot$, respectively. Our results show that, at design sensitivity, both DECIGO and ET should be able to detect these IMBH--SBH mergers. Furthermore, most of the mergers will appear eccentric ($e \gtrsim 0.1$), providing an indication of their dynamical origin.

preprint2020arXiv

GW190412 as a Third-Generation Black Hole Merger from a Super Star Cluster

We explore the possibility that GW190412, a binary black hole merger with a non-equal-mass ratio and significantly spinning primary, was formed through repeated black hole mergers in a dense super star cluster. Using a combination of semi-analytic prescriptions for the remnant spin and recoil kick of black hole mergers, we show that the mass ratio and spin of GW190412 are consistent with a binary black hole whose primary component has undergone two successive mergers from a population of $\sim 10M_{\odot}$ black holes in a high-metallicity environment. We then explore the production of GW190412-like analogs in the CMC Cluster Catalog, a grid of 148 $N$-body star cluster models, as well as a new model, behemoth, with nearly $10^7$ particles and initial conditions taken from a cosmological MHD simulation of galaxy formation. We show that the production of binaries with GW190412-like masses and spins is dominated by massive super star clusters with high metallicities and large central escape speeds. While many are observed in the local universe, our results suggest that a careful treatment of these massive clusters, many of which may have been disrupted before the present day, is necessary to characterize the production of unique gravitational-wave events produced through dynamics.

preprint2020arXiv

Merging Black Holes in the Low-mass and High-mass Gaps from 2+2 Quadruple Systems

The origin of the black hole (BH) binary mergers observed by LIGO-Virgo is still uncertain, as are the boundaries of the stellar BH mass function. Stellar evolution models predict a dearth of BHs both at masses $\gtrsim 50$ M$_\odot$ and $\lesssim 5$ M$_\odot$, thus leaving low- and high-mass gaps in the BH mass function. A natural way to form BHs of these masses is through mergers of neutron stars (NSs; for the low-mass gap) or lower-mass BHs (for the high-mass gap); the low- or high-mass-gap BH produced as a merger product can then be detected by LIGO-Virgo if it merges again with a new companion. We show that the evolution of a 2+2 quadruple system can naturally lead to BH mergers with component masses in the low- or high-mass gaps. In our scenario, the BH in the mass gap originates from the merger of two NSs, or two BHs, in one of the two binaries and the merger product is imparted a recoil kick (from anisotropic gravitational wave emission), which triggers its interaction with the other binary component of the quadruple system. The outcome of this three-body interaction is usually a new eccentric compact binary containing the BH in the mass gap, which can then merge again. The merger rate is $\sim 10^{-7} - 10^{-2}$ Gpc$^{-3}$ yr$^{-1}$ and $\sim 10^{-3} - 10^{-2}$ Gpc$^{-3}$ yr$^{-1}$ for BHs in the low-mass and high-mass gap, respectively. As the sensitivity of gravitational wave detectors improves, tighter constraints will soon be placed on the stellar BH mass function.

preprint2020arXiv

Modeling Dense Star Clusters in the Milky Way and Beyond with the $\texttt{CMC}$ Cluster Catalog

We present a set of 148 independent $N$-body simulations of globular clusters (GCs) computed using the code $\texttt{CMC}$ ($\texttt{Cluster Monte Carlo}$). At an age of $\sim10-13\,$Gyr, the resulting models cover nearly the full range of cluster properties exhibited by the Milky Way GCs, including total mass, core and half-light radii, metallicity, and galactocentric distance. We use our models to investigate the role that stellar-mass black holes play in the process of core collapse. Furthermore, we study how dynamical interactions affect the formation and evolution of several important types of sources in GCs, including low-mass X-ray binaries, millisecond pulsars, blue stragglers, cataclysmic variables, Type Ia supernovae, calcium-rich transients, and merging compact binaries. While our focus here is on old, low-metallicity GCs, our $\texttt{CMC}$ simulations follow the evolution of clusters over a Hubble time, and they include a wide range of metallicities (up to solar), so that our results can also be used to study younger and higher-metallicity star clusters.

preprint2019arXiv

On the Rate of Neutron Star Binary Mergers from Globular Clusters

The first detection of gravitational waves from a neutron star - neutron star (NS-NS) merger, GW170817, and the increasing number of observations of short gamma-ray bursts (SGRBs) have greatly motivated studies of the origins of NS-NS and neutron star - black hole (NS-BH) binaries. We calculate the merger rates of NS-NS and NS-BH binaries from globular clusters (GCs) using realistic GC simulations with the \texttt{CMC} cluster catalog. We use a large sample of models with a range of initial numbers of stars, metallicities, virial radii and galactocentric distances, representative of the present-day Milky Way GCs, to quantify the inspiral times and volumetric merger rates as a function of redshift, both inside and ejected from clusters. We find that over the complete lifetime of most GCs, stellar BHs dominate the cluster cores and prevent the mass segregation of NSs, thereby reducing the dynamical interaction rates of NSs so that at most a few NS binary mergers are ever produced. We estimate the merger rate in the local universe to be $\sim\rm{0.02\,Gpc^{-3}\,yr^{-1}}$ for both NS-NS and NS-BH binaries, or a total of $\sim 0.04$~Gpc$^{-3}$~yr$^{-1}$ for both populations. These rates are about 5 orders of magnitude below the current empirical merger rate from LIGO/Virgo. We conclude that dynamical interactions in GCs do not play a significant role in enhancing the NS-NS and NS-BH merger rates.

preprint2019arXiv

The Role of "Black Hole Burning" in the Evolution of Dense Star Clusters

As self-gravitating systems, dense star clusters exhibit a natural diffusion of energy from their innermost to outermost regions, which leads to a slow and steady contraction of the core until it ultimately collapses under gravity. However, in spite of the natural tendency toward so-called "core collapse," the globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the core-collapsed and non-core-collapsed clusters. This suggests an internal energy source is at work, delaying the onset of core collapse in many clusters. Primordial binary stars have been thought for a long time to provide this energy source, but recent analyses have cast doubt upon the corresponding "binary-burning" mechanism as a viable explanation. Over the past decade, a large amount of both observational and theoretical work has suggested that many stellar-mass black holes (BHs) are retained in typical clusters today and that they play a dynamically-significant role in these clusters throughout their entire lifetimes. Here we review our latest understanding of the formation and evolution of BH populations in GCs and demonstrate that, through their dynamical interaction with their host cluster, BHs can naturally explain the distinction between core-collapsed and non-core-collapsed clusters through a process we call "black hole burning."

preprint2016arXiv

Binary Black Hole Mergers from Globular Clusters: Masses, Merger Rates, and the Impact of Stellar Evolution

The recent discovery of GW150914, the binary black hole merger detected by Advanced LIGO, has the potential to revolutionize observational astrophysics. But to fully utilize this new window into the universe, we must compare these new observations to detailed models of binary black hole formation throughout cosmic time. Expanding upon our previous work (Rodriguez et al., 2015), we study merging binary black holes formed in globular clusters using our Monte Carlo approach to stellar dynamics. We have created a new set of 52 cluster models with different masses, metallicities, and radii to fully characterize the binary black hole merger rate. These models include all the relevant dynamical processes (such as two-body relaxation, strong encounters, and three-body binary formation) and agree well with detailed direct N-body simulations. In addition, we have enhanced our stellar evolution algorithms with updated metallicity-dependent stellar wind and supernova prescriptions, allowing us to compare our results directly to the most recent population synthesis predictions for merger rates from isolated binary evolution. We explore the relationship between a cluster's global properties and the population of binary black holes that it produces. In particular, we derive a numerically calibrated relationship between the merger times of ejected black hole binaries and a cluster's mass and radius. We explore the masses and mass ratios of these binaries as a function of redshift, and find a merger rate of ~5 Gpc$^{-3}$ yr$^{-1}$ in the local universe, with 80% of sources having total masses from $32M_{\odot}$ to $64M_{\odot}$. Under standard assumptions, approximately 1 out of every 7 binary black hole mergers in the local universe will have originated in a globular cluster, but we also explore the sensitivity of this result to different assumptions for binary stellar evolution. [Abridged]

preprint2016arXiv

Distinguishing Between Formation Channels for Binary Black Holes with LISA

The recent detections of GW150914 and GW151226 imply an abundance of stellar-mass binary-black-hole mergers in the local universe. While ground-based gravitational-wave detectors are limited to observing the final moments before a binary merges, space-based detectors, such as the Laser Interferometer Space Antenna (LISA), can observe binaries at lower orbital frequencies where such systems may still encode information about their formation histories. In particular, the orbital eccentricity and mass of binary black holes in the LISA frequency band can be used together to discriminate between binaries formed in isolation in galactic fields and those formed in dense stellar environments such as globular clusters. In this letter, we explore the orbital eccentricity and mass of binary-black-hole populations as they evolve through the LISA frequency band. Overall we find that there are two distinct populations discernible by LISA. We show that up to ~90% of binaries formed either dynamically or in isolation have eccentricities measurable by LISA. Finally, we note how measured eccentricities of low-mass binary black holes evolved in isolation could provide detailed constraints on the physics of black-hole natal kicks and common-envelope evolution.

preprint2016arXiv

Dynamical Formation of the GW150914 Binary Black Hole

We explore the possibility that GW150914, the binary black hole merger recently detected by Advanced LIGO, was formed by gravitational interactions in the core of a dense star cluster. Using models of globular clusters with detailed $N$-body dynamics and stellar evolution, we show that a typical cluster with a mass of $3\times10^5M_{\odot}$ to $6\times10^5M_{\odot}$ is optimal for forming GW150914-like binary black holes that will merge in the local universe. We identify the most likely dynamical processes for forming GW150914 in such a cluster, and we show that the detection of GW150914 is consistent with the masses and merger rates expected for binary black hole mergers from globular clusters. Our results show that dynamical processes provide a significant and well-understood pathway for forming binary black hole mergers in the local universe. Understanding the contribution of dynamics to the binary black hole merger problem is a critical step in unlocking the full potential of gravitational-wave astronomy.

preprint2016arXiv

Illuminating Black Hole Binary Formation Channels with Spins in Advanced LIGO

The recent detections of the binary black hole mergers GW150914 and GW151226 have inaugurated the field of gravitational-wave astronomy. For the two main formation channels that have been proposed for these sources, isolated binary evolution in galactic fields and dynamical formation in dense star clusters, the predicted masses and merger rates overlap significantly, complicating any astrophysical claims that rely on measured masses alone. Here, we examine the distribution of spin- orbit misalignments expected for binaries from the field and from dense star clusters. Under standard assumptions for black-hole natal kicks, we find that black-hole binaries similar to GW150914 could be formed with significant spin-orbit misalignment only through dynamical processes. In particular, these heavy-black-hole binaries can only form with a significant spin-orbit anti-alignment in the dynamical channel. Our results suggest that future detections of merging black hole binaries with measurable spins will allow us to identify the main formation channel for these systems.

preprint2016arXiv

Merging black hole binaries in galactic nuclei: implications for advanced-LIGO detections

Motivated by the recent detection of gravitational waves from the black hole binary merger GW150914, we study the dynamical evolution of black holes in galactic nuclei where massive star clusters reside. With masses of ~10^7M_Sun and sizes of only a few parsecs, nuclear star clusters are the densest stellar systems observed in the local universe and represent a robust environment where (stellar mass) black hole binaries can dynamically form, harden and merge. We show that due to their large escape speeds, nuclear star clusters can keep a large fraction of their merger remnants while also evolving rapidly enough that the holes can sink back to the central regions where they can swap in new binaries that can subsequently harden and merge. This process can repeat several times and produce black hole mergers of several tens of solar masses similar to GW150914 and up to a few hundreds of solar masses, without the need of invoking extremely low metallicity environments or implausible initial conditions. We use a semi-analytical approach to describe the formation and dynamics of black holes in massive star clusters. We find a black hole binary merger rate per volume from nuclear star clusters of ~1.5 Gpc^-3 yr^-1, implying up to a few tens of possible detections per year with Advanced LIGO. Our models suggest a local merger rate of 0.3- 1 Gpc^-3 yr^-1 for high mass black hole binaries similar to GW150914 (total mass >~ 50 M_Sun, redshift z< 0.3); a merger rate comparable to that of high mass black hole binaries that are dynamically assembled in globular clusters. Finally, we show that if all black holes receive high natal kicks, >~50km s^-1, then nuclear star clusters could dominate the local merger rate of binary black holes compared to the merger rate of similar binaries produced in either globular clusters or through isolated binary evolution.

preprint2016arXiv

Million-Body Star Cluster Simulations: Comparisons between Monte Carlo and Direct $N$-body

We present the first detailed comparison between million-body globular cluster simulations computed with a Hénon-type Monte Carlo code, CMC, and a direct $N$-body code, NBODY6++GPU. Both simulations start from an identical cluster model with $10^6$ particles, and include all of the relevant physics needed to treat the system in a highly realistic way. With the two codes "frozen" (no fine-tuning of any free parameters or internal algorithms of the codes) we find excellent agreement in the overall evolution of the two models. Furthermore, we find that in both models, large numbers of stellar-mass black holes (> 1000) are retained for 12 Gyr. Thus, the very accurate direct $N$-body approach confirms recent predictions that black holes can be retained in present-day, old globular clusters. We find only minor disagreements between the two models and attribute these to the small-$N$ dynamics driving the evolution of the cluster core for which the Monte Carlo assumptions are less ideal. Based on the overwhelming general agreement between the two models computed using these vastly different techniques, we conclude that our Monte Carlo approach, which is more approximate, but dramatically faster compared to the direct $N$-body, is capable of producing a very accurate description of the long-term evolution of massive globular clusters even when the clusters contain large populations of stellar-mass black holes.

preprint2015arXiv

Binary Black Hole Mergers from Globular Clusters: Implications for Advanced LIGO

The predicted rate of binary black hole mergers from galactic fields can vary over several orders of magnitude and is extremely sensitive to the assumptions of stellar evolution. But in dense stellar environments such as globular clusters, binary black holes form by well-understood gravitational interactions. In this letter, we study the formation of black hole binaries in an extensive collection of realistic globular cluster models. By comparing these models to observed Milky Way and extragalactic globular clusters, we find that the mergers of dynamically-formed binaries could be detected at a rate of ~100 per year, potentially dominating the binary black hole merger rate. We also find that a majority of cluster-formed binaries are more massive than their field-formed counterparts, suggesting that Advanced LIGO could identify certain binaries as originating from dense stellar environments.

preprint2015arXiv

Black hole mergers and blue stragglers from hierarchical triples formed in globular clusters

Hierarchical triple-star systems are expected to form frequently via close binary-binary encounters in the dense cores of globular clusters. In a sufficiently inclined triple, gravitational interactions between the inner and outer binary can cause large-amplitude oscillations in the eccentricity of the inner orbit ("Lidov-Kozai cycles"), which can lead to a collision and merger of the two inner components. In this paper we use Monte Carlo models of dense star clusters to identify all triple systems formed dynamically and we compute their evolution using a highly accurate three-body integrator which incorporates relativistic and tidal effects. We find that a large fraction of these triples evolve through a non-secular dynamical phase which can drive the inner binary to higher eccentricities than predicted by the standard secular perturbation theory (even including octupole-order terms). We place constraints on the importance of Lidov-Kozai-induced mergers for producing: (i) gravitational wave sources detectable by Advanced LIGO (aLIGO), for triples with an inner pair of stellar black holes; and (ii) blue straggler stars, for triples with main-sequence-star components. We find a realistic aLIGO detection rate of black hole mergers due to the Lidov-Kozai mechanism of 1yr^-1, with about 20% of these having a finite eccentricity when they first chirp into the aLIGO frequency band. While rare, these events are likely to dominate among eccentric compact object inspirals that are potentially detectable by aLIGO. For blue stragglers, we find that the Lidov-Kozai mechanism can contribute only up to ~10% of their total numbers in globular clusters.

preprint2015arXiv

Formation of Black Hole Low-Mass X-ray Binaries in Hierarchical Triple Systems

The formation of Black Hole (BH) Low-Mass X-ray Binaries (LMXB) poses a theoretical challenge, as low-mass companions are not expected to survive the common-envelope scenario with the BH progenitor. Here we propose a formation mechanism that skips the common-envelope scenario and relies on triple-body dynamics. We study the evolution of hierarchical triples, following the secular dynamical evolution up to the octupole-level of approximation, including general relativity, tidal effects and post-main-sequence evolution, such as mass loss, changes to stellar radii and supernovae. During the dynamical evolution of the triple system, the "eccentric Kozai-Lidov" mechanism can cause large eccentricity excitations in the LMXB progenitor, resulting in three main BH-LMXB formation channels. Here we define BH-LMXB candidates as systems where the inner BH companion star crosses its Roche limit. In the "eccentric" channel (~ 81% of the LMXBs in our simulations), the donor star crosses its Roche limit during an extreme eccentricity excitation, while still on a wide orbit. Second, we find a "giant" LMXB channel (~ 11%), where a system undergoes only moderate eccentricity excitations, but the donor star fills its Roche lobe after evolving toward the giant branch. Third, we identify a "classical" channel (~8%), where tidal forces and magnetic braking shrink and circularize the orbit to short periods, triggering mass transfer. Finally, for the giant channel, we predict an eccentric ($\sim 0.3-0.6$), preferably inclined (~40, ~140 degreed) tertiary, typically on a wide enough orbit (~10^4AU), to potentially become unbound later in the triple evolution.

preprint2015arXiv

Tidally-driven Roche-Lobe Overflow of Hot Jupiters with MESA

Many exoplanets have now been detected in orbits with ultra-short periods, very close to the Roche limit. Building upon our previous work, we study the possibility that mass loss through Roche lobe overflow (RLO) may affect the evolution of these planets, and could possibly transform a hot Jupiter into a lower-mass planet (hot Neptune or super-Earth). We focus here on systems in which the mass loss occurs slowly ("stable mass transfer" in the language of binary star evolution) and we compute their evolution in detail with the binary evolution code MESA. We include the effects of tides, RLO, irradiation and photo-evaporation of the planet, as well as the stellar wind and magnetic braking. Our calculations all start with a hot Jupiter close to its Roche limit, in orbit around a sun-like star. The initial orbital decay and onset of RLO are driven by tidal dissipation in the star. We confirm that such a system can indeed evolve to produce lower-mass planets in orbits of a few days. The RLO phase eventually ends and, depending on the details of the mass transfer and on the planetary core mass, the orbital period can remain around a few days for several Gyr. The remnant planets have a rocky core and some amount of envelope material, which is slowly removed via photo-evaporation at nearly constant orbital period; these have properties resembling many of the observed super-Earths and sub-Neptunes. For these remnant planets we also predict an anti-correlation between mass and orbital period; very low-mass planets ($M_{\rm pl}\,\lesssim\,5\,M_{\oplus}$) in ultra-short periods ($P_{\rm orb}$<1d) cannot be produced through this type of evolution.

preprint2014arXiv

From Hot Jupiters to Super-Earths via Roche Lobe Overflow

Through tidal dissipation in a slowly spinning host star the orbits of many hot Jupiters may decay down to the Roche limit. We expect that in most cases the ensuing mass transfer will be stable. Using detailed numerical calculations we find that this evolution is quite rapid, potentially leading to complete removal of the gaseous envelope in a few Gyr, and leaving behind an exposed rocky core ("hot super-Earth"). Final orbital periods are quite sensitive to the details of the planet's mass-radius relation, and to the effects of irradiation and photo-evaporation, but could be as short as a few hours, or as long as several days. Our scenario predicts the existence of planets with intermediate masses ("hot Neptunes") that should be found precisely at their Roche limit and in the process of losing mass through Roche lobe overflow. The observed excess of small single-planet candidate systems observed by Kepler may also be the result of this process. If so, the properties of their host stars should track those of the hot Jupiters. Moreover, the number of systems that produced hot Jupiters could be 2-3 times larger than one would infer from contemporary observations.

preprint2014arXiv

Planets on the Edge

Hot Jupiters formed through circularization of high-eccentricity orbits should be found at orbital separations $a$ exceeding $twice$ that of their Roche limit $a_{\rm R}$. Nevertheless, about a dozen giant planets have now been found well within this limit ($a_{\rm R}<a<2a_{\rm R}$), with one coming as close as 1.2$a_{\rm R}$. In this Letter, we show that orbital decay (starting beyond 2$a_{\rm R}$) driven by tidal dissipation in the star can naturally explain these objects. For a few systems (WASP-4 and 19), this explanation requires the linear reduction in convective tidal dissipation proposed originally by Zahn (1966, 1989) and verified by recent numerical simulations Penev et al. (2007), but rules out the quadratic prescription proposed by Goldreich & Nicholson (1977). Additionally, we find that WASP-19-like systems could potentially provide direct empirical constraints on tidal dissipation, as we could soon be able to measure their orbital decay through high precision transit timing measurements.

preprint2014arXiv

Sowing the seeds of massive black holes in small galaxies: Young clusters as the building blocks of Ultra-Compact-Dwarf Galaxies

Interacting galaxies often have complexes of hundreds of young stellar clusters of individual masses $\sim 10^{4-6}~M_\odot$ in regions that are a few hundred parsecs across. These cluster complexes interact dynamically, and their coalescence is a candidate for the origin of some ultracompact dwarf galaxies (UCDs). Individual clusters with short relaxation times are candidates for the production of intermediate-mass black holes of a few hundred solar masses, via runaway stellar collisions prior to the first supernovae in a cluster. It is therefore possible that a cluster complex hosts multiple intermediate-mass black holes that may be ejected from their individual clusters due to mergers or binary processes, but bound to the complex as a whole. Here we explore the dynamical interaction between initially free-flying massive black holes and clusters in an evolving cluster complex. We find that, after hitting some clusters, it is plausible that the massive black hole will be captured in an ultracompact dwarf forming near the center of the complex. In the process, the hole typically triggers electromagnetic flares via stellar disruptions, and is also likely to be a prominent source of gravitational radiation for the advanced ground-based detectors LIGO and VIRGO. We also discuss other implications of this scenario, notably that the central black hole could be considerably larger than expected in other formation scenarios for ultracompact dwarfs.

preprint2014arXiv

The Dynamical Evolution of Stellar Black Holes in Globular Clusters

Our current understanding of the stellar initial mass function and massive star evolution suggests that young globular clusters may have formed hundreds to thousands of stellar-mass black holes, the remnants of stars with initial masses from $\sim 20 - 100\, M_\odot$. Birth kicks from supernova explosions may eject some black holes from their birth clusters, but most should be retained. Using a Monte Carlo method we investigate the long-term dynamical evolution of globular clusters containing large numbers of stellar black holes. We describe numerical results for 42 models, covering a range of realistic initial conditions, including up to $1.6\times10^6$ stars. In almost all models we find that significant numbers of black holes (up to $\sim10^3$) are retained all the way to the present. This is in contrast to previous theoretical expectations that most black holes should be ejected dynamically within a few Gyr. The main reason for this difference is that core collapse driven by black holes (through the Spitzer "mass segregation instability") is easily reverted through three-body processes, and involves only a small number of the most massive black holes, while lower-mass black holes remain well-mixed with ordinary stars far from the central cusp. Thus the rapid segregation of stellar black holes does not lead to a long-term physical separation of most black holes into a dynamically decoupled inner core, as often assumed previously. Combined with the recent detections of several black hole X-ray binary candidates in Galactic globular clusters, our results suggest that stellar black holes could still be present in large numbers in many globular clusters today, and that they may play a significant role in shaping the long-term dynamical evolution and the present-day dynamical structure of many clusters.

preprint2014arXiv

The Dynamics of the Multi-planet System Orbiting Kepler-56

Kepler-56 is a multi-planet system containing two coplanar inner planets that are in orbits misaligned with respect to the spin axis of the host star, and an outer planet. Various mechanisms have been proposed to explain the broad distribution of spin-orbit angles among exoplanets, and these theories fall under two broad categories. The first is based on dynamical interactions in a multi-body system, while the other assumes that disk migration is the driving mechanism in planetary configuration and that the star (or disk) is titled with respect to the planetary plane. Here we show that the large observed obliquity of Kepler-56 system is consistent with a dynamical origin. In addition, we use observations by Huber et al. (2013) to derive the obliquity's probability distribution function, thus improving the constrained lower limit. The outer planet may be the cause of the inner planets' large obliquities, and we give the probability distribution function of its inclination, which depends on the initial orbital configuration of the planetary system. We show that even in the presence of precise measurement of the true obliquity, one cannot distinguish the initial configurations. Finally we consider the fate of the system as the star continues to evolve beyond the main sequence, and we find that the obliquity of the system will not undergo major variations as the star climbs the red giant branch. We follow the evolution of the system and find that the innermost planet will be engulfed in ~129 Myr. Furthermore we put an upper limit of ~155 Myr for the engulfment of the second planet. This corresponds to ~ 3% of the current age of the star.

preprint2014arXiv

Tidal Dissipation and Obliquity Evolution in Hot Jupiter Systems

Two formation scenarios have been proposed to explain the tight orbits of hot Jupiters. They could be formed in orbits with a small inclination (with respect to the stellar spin) via disk migration, or in more highly inclined orbits via high-eccentricity migration, where gravitational interactions with a companion and tidal dissipation are at play. Here we target hot Jupiter systems where the misalignment $λ$ has been inferred observationally and we investigate whether their properties are consistent with high-eccentricity migration. Specifically, we study whether stellar tides can be responsible for the observed distribution of $λ$ and orbital separations. Improving on previous studies, we use detailed models for each star, thus accounting for how convection (and tidal dissipation) depends on stellar properties. In line with observations suggesting that hotter stars have higher $λ$, we find that $λ$ increases as the amount of stellar surface convection decreases. This trend supports the hypothesis that tides are the mechanism shaping the observed distribution of $λ$. Furthermore, we study the past orbital evolution of five representative systems, chosen to cover a variety of temperatures and misalignments. We consider various initial orbital configurations and integrate the equations describing the coupled evolution of the orbital separation, stellar spin, and misalignment. We account for stellar tides and wind mass loss, stellar evolution, and magnetic braking. We show that the current properties of these five representative systems can be explained naturally, given our current understanding of tidal dissipation and with physically motivated assumptions for the effects driving the orbital evolution.

preprint2013arXiv

Extreme orbital evolution from hierarchical secular coupling of two giant planets

Observations of exoplanets over the last two decades have revealed a new class of Jupiter-size planets with orbital periods of a few days, the so-called "hot Jupiters". Recent measurements using the Rossiter-McLaughlin effect have shown that many (~ 50%) of these planets are misaligned; furthermore, some (~ 15%) are even retrograde with respect to the stellar spin axis. Motivated by these observations, we explore the possibility of forming retrograde orbits in hierarchical triple configurations consisting of a star-planet inner pair with another giant planet, or brown dwarf, in a much wider orbit. Recently Naoz et al. (2011) showed that in such a system, the inner planet's orbit can flip back and forth from prograde to retrograde, and can also reach extremely high eccentricities. Here we map a significant part of the parameter space of dynamical outcomes for these systems. We derive strong constraints on the orbital configurations for the outer perturber that could lead to the formation of hot Jupiters with misaligned or retrograde orbits. We focus only on the secular evolution, neglecting other dynamical effects such as mean-motion resonances, as well as all dissipative forces. For example, with an inner Jupiter-like planet initially on a nearly circular orbit at 5 AU, we show that a misaligned hot Jupiter is likely to be formed in the presence of a more massive planetary companion (> 2 MJ) within 140 AU of the inner system, with mutual inclination 50 degrees and eccentricity above 0.25. This is in striking contrast to the test-particle approximation, where an almost perpendicular configuration can still cause large eccentricity excitations, but flips of an inner Jupiter-like planet are much less likely to occur. The constraints we derive can be used to guide future observations, and, in particular, searches for more distant companions in systems containing a hot Jupiter.

preprint2013arXiv

Introducing CAFein, a New Computational Tool for Stellar Pulsations and Dynamic Tides

Here we present CAFein, a new computational tool for investigating radiative dissipation of dynamic tides in close binaries and of non-adiabatic, non-radial stellar oscillations in isolated stars in the linear regime. For the latter, CAFein computes the non-adiabatic eigenfrequencies and eigenfunctions of detailed stellar models. The code is based on the so-called Riccati method, a numerical algorithm that has been successfully applied to a variety of stellar pulsators, and which doesn't suffer of the major drawbacks of commonly used shooting and relaxation schemes. Here we present an extension of the Riccati method to investigate dynamic tides in close binaries. We demonstrate CAFein's capabilities as a stellar pulsation code both in the adiabatic and non-adiabatic regime, by reproducing previously published eigenfrequencies of a polytrope, and by successfully identifying the unstable modes of a stellar model in the $β$ Cephei/SPB region of the Hertzsprung-Russell diagram. Finally, we verify CAFein's behavior in the dynamic tides regime by investigating the effects of dynamic tides on the eigenfunctions and orbital and spin evolution of massive Main Sequence stars in eccentric binaries, and of hot Jupiter host stars. The plethora of asteroseismic data provided by the NASA's Kepler satellite, some of which include the direct detection of tidally excited stellar oscillations, make CAFein quite timely. Furthermore, the increasing number of observed short-period detached double white dwarfs (WD) and the observed orbital decay in the tightest of such binaries open up a new possibility of investigating WD interiors through the effects of tides on their orbital evolution

preprint2013arXiv

Mock Observations of Blue Stragglers in Globular Cluster Models

We created artificial color-magnitude diagrams of Monte Carlo dynamical models of globular clusters, and then used observational methods to determine the number of blue stragglers in those clusters. We compared these blue stragglers to various cluster properties, mimicking work that has been done for blue stragglers in Milky Way globular clusters to determine the dominant formation mechanism(s) of this unusual stellar population. We find that a mass-based prescription for selecting blue stragglers will choose approximately twice as many blue stragglers than a selection criterion that was developed for observations of real clusters. However, the two numbers of blue stragglers are well-correlated, so either selection criterion can be used to characterize the blue straggler population of a cluster. We confirm previous results that the simplified prescription for the evolution of a collision or merger product in the BSE code overestimates their lifetimes. We show that our model blue stragglers follow similar trends with cluster properties (core mass, binary fraction, total mass, collision rate) as the true Milky Way blue stragglers, as long as we restrict ourselves to model clusters with an initial binary fraction higher than 5%. We also show that, in contrast to earlier work, the number of blue stragglers in the cluster core does have a weak dependence on the collisional parameter Gamma in both our models and in Milky Way globular clusters.

preprint2013arXiv

Retention of Stellar-Mass Black Holes in Globular Clusters

Globular clusters should be born with significant numbers of stellar-mass black holes (BHs). It has been thought for two decades that very few of these BHs could be retained through the cluster lifetime. With masses ~10 MSun, BHs are ~20 times more massive than an average cluster star. They segregate into the cluster core, where they may eventually decouple from the remainder of the cluster. The small-N core then evaporates on a short timescale. This is the so-called Spitzer instability. Here we present the results of a full dynamical simulation of a globular cluster containing many stellar-mass BHs with a realistic mass spectrum. Our Monte Carlo simulation code includes detailed treatments of all relevant stellar evolution and dynamical processes. Our main finding is that old globular clusters could still contain many BHs at present. In our simulation, we find no evidence for the Spitzer instability. Instead, most of the BHs remain well-mixed with the rest of the cluster, with only the innermost few tens of BHs segregating significantly. Over the 12 Gyr evolution, fewer than half of the BHs are dynamically ejected through strong binary interactions in the cluster core. The presence of BHs leads to long-term heating of the cluster, ultimately producing a core radius on the high end of the distribution for Milky Way globular clusters (and those of other galaxies). A crude extrapolation from our model suggests that the BH--BH merger rate from globular clusters could be comparable to the rate in the field.

preprint2013arXiv

Secular Dynamics in Hierarchical Three-Body Systems

The secular approximation for the evolution of hierarchical triple configurations has proven to be very useful in many astrophysical contexts, from planetary to triple-star systems. In this approximation the orbits may change shape and orientation, on time scales longer than the orbital time scales, but the semi major axes are constant. For example, for highly inclined triple systems, the Kozai-Lidov mechanism can produce large-amplitude oscillations of the eccentricities and inclinations. Here we revisit the secular dynamics of hierarchical triple systems. We derive the secular evolution equations to octupole order in Hamiltonian perturbation theory. Our derivation corrects an error in some previous treatments of the problem that implicitly assumed a conservation of the z-component of the angular momentum of the inner orbit (i.e., parallel to the total angular momentum of the system). Already to quadrupole order, our results show new behaviors including the possibility for a system to oscillate from prograde to retrograde orbits. At the octupole order, for an eccentric outer orbit, the inner orbit can reach extremely high eccentricities and undergo chaotic flips in its orientation. We discuss applications to a variety of astrophysical systems, from stellar triples to merging compact binaries and planetary systems. Our results agree with those of previous studies done to quadrupole order only in the limit in which one of the inner two bodies is a massless test particle and the outer orbit is circular;our results agree with previous studies at octupole order for the eccentricity evolution, but not for the inclination evolution.

preprint2013arXiv

Stellar Collisions and Blue Straggler Stars in Dense Globular Clusters

Blue straggler stars (BSS) are abundantly observed in all Galactic globular clusters (GGC) where data exist. However, observations alone cannot reveal the relative importance of various formation channels or the typical formation times for this well studied population of anomalous stars. Using a state-of-the-art Hénon-type Monte Carlo code that includes all relevant physical processes, we create 128 models with properties typical of the observed GGCs. These models include realistic numbers of single and binary stars, use observationally motivated initial conditions, and span large ranges in central density, concentration, binary fraction, and mass. Their properties can be directly compared with those of observed GGCs. We can easily identify the BSSs in our models and determine their formation channels and birth times. We find that for central densities above ~10^3 Msun/pc^3 the dominant formation channel is stellar collisions while for lower density clusters, mass transfer in binaries provides a significant contribution (up to ~ 60% in our models). The majority of these collisions are binary-mediated, occurring during 3-body and 4-body interactions. As a result a strong correlation between the specific frequency of BSSs and the binary fraction in a cluster can be seen in our models. We find that the number of BSSs in the core shows only a weak correlation with the collision rate estimator Γtraditionally used by observers, in agreement with the latest Hubble Space Telescope (ACS) data. Using an idealized "full mixing" prescription for collision products, our models indicate that the BSSs observed today may have formed several Gyrs ago. However, denser clusters tend to have younger (~1 Gyr) BSSs.

preprint2012arXiv

A Parallel Monte Carlo Code for Simulating Collisional N-body Systems

We present a new parallel code for computing the dynamical evolution of collisional N-body systems with up to N~10^7 particles. Our code is based on the the Henon Monte Carlo method for solving the Fokker-Planck equation, and makes assumptions of spherical symmetry and dynamical equilibrium. The principal algorithmic developments involve optimizing data structures, and the introduction of a parallel random number generation scheme, as well as a parallel sorting algorithm, required to find nearest neighbors for interactions and to compute the gravitational potential. The new algorithms we introduce along with our choice of decomposition scheme minimize communication costs and ensure optimal distribution of data and workload among the processing units. The implementation uses the Message Passing Interface (MPI) library for communication, which makes it portable to many different supercomputing architectures. We validate the code by calculating the evolution of clusters with initial Plummer distribution functions up to core collapse with the number of stars, N, spanning three orders of magnitude, from 10^5 to 10^7. We find that our results are in good agreement with self-similar core-collapse solutions, and the core collapse times generally agree with expectations from the literature. Also, we observe good total energy conservation, within less than 0.04% throughout all simulations. We analyze the performance of the code, and demonstrate near-linear scaling of the runtime with the number of processors up to 64 processors for N=10^5, 128 for N=10^6 and 256 for N=10^7. The runtime reaches a saturation with the addition of more processors beyond these limits which is a characteristic of the parallel sorting algorithm. The resulting maximum speedups we achieve are approximately 60x, 100x, and 220x, respectively.

preprint2012arXiv

Binary Neutron Star Mergers

We review the current status of studies of the coalescence of binary neutron star systems. We begin with a discussion of the formation channels of merging binaries and we discuss the most recent theoretical predictions for merger rates. Next, we turn to the quasi-equilibrium formalisms that are used to study binaries prior to the merger phase and to generate initial data for fully dynamical simulations. The quasi-equilibrium approximation has played a key role in developing our understanding of the physics of binary coalescence and, in particular, of the orbital instability processes that can drive binaries to merger at the end of their lifetimes. We then turn to the numerical techniques used in dynamical simulations, including relativistic formalisms, (magneto-)hydrodynamics, gravitational-wave extraction techniques, and nuclear microphysics treatments. This is followed by a summary of the simulations performed across the field to date, including the most recent results from both fully relativistic and microphysically detailed simulations. Finally, we discuss the likely directions for the field as we transition from the first to the second generation of gravitational-wave interferometers and while supercomputers reach the petascale frontier.

preprint2012arXiv

Constraining Intermediate-Mass Black Holes in Globular Clusters

Decades after the first predictions of intermediate-mass black holes (IMBHs) in globular clusters (GCs) there is still no unambiguous observational evidence for their existence. The most promising signatures for IMBHs are found in the cores of GCs, where the evidence now comes from the stellar velocity distribution, the surface density profile, and, for very deep observations, the mass-segregation profile near the cluster center. However, interpretation of the data, and, in particular, constraints on central IMBH masses, require the use of detailed cluster dynamical models. Here we present results from Monte Carlo cluster simulations of GCs that harbor IMBHs. As an example of application, we compare velocity dispersion, surface brightness and mass-segregation profiles with observations of the GC M10, and constrain the mass of a possible central IMBH in this cluster. We find that, although M10 does not seem to possess a cuspy surface density profile, the presence of an IMBH with a mass up to 0.75% of the total cluster mass, corresponding to about 600 Msun, cannot be excluded. This is also in agreement with the surface brightness profile, although we find it to be less constraining, as it is dominated by the light of giants, causing it to fluctuate significantly. We also find that the mass-segregation profile cannot be used to discriminate between models with and without IMBH. The reason is that M10 is not yet dynamically evolved enough for the quenching of mass segregation to take effect. Finally, detecting a velocity dispersion cusp in clusters with central densities as low as in M10 is extremely challenging, and has to rely on only 20-40 bright stars. It is only when stars with masses down to 0.3 Msun are included that the velocity cusp is sampled close enough to the IMBH for a significant increase above the core velocity dispersion to become detectable.

preprint2012arXiv

Formation of Massive Black Holes in Dense Star Clusters. II. IMF and Primordial Mass Segregation

A promising mechanism to form intermediate-mass black holes (IMBHs) is the runaway merger in dense star clusters, where main-sequence stars collide and form a very massive star (VMS), which then collapses to a black hole. In this paper we study the effects of primordial mass segregation and the importance of the stellar initial mass function (IMF) on the runaway growth of VMSs using a dynamical Monte Carlo code for N-body systems with N as high as 10^6 stars. Our code now includes an explicit treatment of all stellar collisions. We place special emphasis on the possibility of top-heavy IMFs, as observed in some very young massive clusters. We find that both primordial mass segregation and the shape of the IMF affect the rate of core collapse of star clusters and thus the time of the runaway. When we include primordial mass segregation we generally see a decrease in core collapse time (tcc). Moreover, primordial mass segregation increases the average mass in the core, thus reducing the central relaxation time, which also decreases tcc. The final mass of the VMS formed is always close to \sim 10^-3 of the total cluster mass, in agreement with the previous studies and is reminiscent of the observed correlation between the central black hole mass and the bulge mass of the galaxies. As the degree of primordial mass segregation is increased, the mass of the VMS increases at most by a factor of 3. Flatter IMFs generally increase the average mass in the whole cluster, which increases tcc. For the range of IMFs investigated in this paper, this increase in tcc is to some degree balanced by stellar collisions, which accelerate core collapse. Thus there is no significant change in tcc for the somewhat flatter global IMFs observed in very young massive clusters.

preprint2012arXiv

Monte Carlo Simulations of Globular Cluster Evolution. VI. The Influence of an Intermediate Mass Black Hole

We present results of a series of Monte Carlo simulations investigating the imprint of a central intermediate-mass black hole (IMBH) on the structure of a globular cluster. We investigate the three-dimensional and projected density profiles, and stellar disruption rates for idealized as well as realistic cluster models, taking into account a stellar mass spectrum and stellar evolution, and allowing for a larger, more realistic, number of stars than was previously possible with direct N-body methods. We compare our results to other N-body and Fokker-Planck simulations published previously. We find, in general, very good agreement for the overall cluster structure and dynamical evolution between direct N-body simulations and our Monte Carlo simulations. Significant differences exist in the number of stars that are tidally disrupted by the IMBH, which is most likely an effect of the wandering motion of the IMBH, not included in the Monte Carlo scheme. These differences, however, are negligible for the final IMBH masses in realistic cluster models as the disruption rates are generally much lower than for single-mass clusters. As a direct comparison to observations we construct a detailed model for the cluster NGC 5694, which is known to possess a central surface brightness cusp consistent with the presence of an IMBH. We find that not only the inner slope but also the outer part of the surface brightness profile agree well with observations. However, there is only a slight preference for models harboring an IMBH compared to models without.

preprint2012arXiv

On the Formation of Hot Jupiters in Stellar Binaries

We study the production of Hot Jupiters (HJs) in stellar binaries. We show that the "eccentric Kozai-Lidov" (EKL) mechanism can play a key role in the dynamical evolution of a star-planet-star triple system. We run a large set of Monte Carlo simulations including the secular evolution of the orbits, general relativistic precession, and tides, and we determine the semi-major axis, eccentricity, inclination and spin-orbit angle distributions of the HJs that are produced. We explore the effect of different tidal friction parameters on the results. We find that the efficiency of forming HJs when taking the EKL mechanism into account is higher then previously estimated. Accounting for the frequency of stellar binaries, we find that this production mechanism can account for about 30% of the observed HJ population. Current observations of spin-orbit angles are consistent with this mechanism producing \sim 30% of all HJs, and up to 100% of the misaligned systems. Based on the properties of binaries without a HJ in our simulations, we predict the existence of many Jupiter-like planets with moderately eccentric and inclined orbits and semi-major axes of several AU.

preprint2012arXiv

Planets in Open Clusters Detectable by Kepler

While hundreds of planets have been discovered around field stars, only a few are known in star clusters. To explain the lack of short-period giant planets in globular clusters (GC), such as 47 Tucane and ω Centauri, it has been suggested that their low metallicities may have prevented planet formation. Alternatively, the high rates of close stellar encounters in these clusters may have influenced the formation and subsequent evolution of planetary systems. How common are planets in clusters around normal main-sequence stars? Here we consider whether this question can be addressed using data from the Kepler mission. The Kepler field of view contains 4 low-density (relative to GCs) open clusters where the metallicities are about solar (or even higher) and stellar encounters are much less frequent than in typical GCs. We provide detailed $N$-body models and show that most planets in Kepler-detectable orbits are not significantly perturbed by stellar encounters in these open clusters. We focus on the most massive cluster, NGC 6791, which has super-solar metallicity, and find that if planets formed in this cluster at the same frequency as observed in the field, Kepler could detect 1 -- 20 transiting planets depending on the planet-size distribution and the duration of data collection. Due to the large distance to NGC 6791 Kepler will have to search relatively faint ($K_p<20$) stars for the full extended mission to achieve such a yield.

preprint2012arXiv

Understanding the Dynamical State of Globular Clusters: Core-Collapsed vs Non Core-Collapsed

We study the dynamical evolution of globular clusters using our Hénon-type Monte Carlo code for stellar dynamics including all relevant physics such as two-body relaxation, single and binary stellar evolution, Galactic tidal stripping, and strong interactions such as physical collisions and binary mediated scattering. We compute a large database of several hundred models starting from broad ranges of initial conditions guided by observations of young and massive star clusters. We show that these initial conditions very naturally lead to present day clusters with properties including the central density, core radius, half-light radius, half-mass relaxation time, and cluster mass, that match well with those of the old Galactic globular clusters. In particular, we can naturally reproduce the bimodal distribution in observed core radii separating the "core-collapsed" vs the "non core-collapsed" clusters. We see that the core-collapsed clusters are those that have reached or are about to reach the equilibrium "binary burning" phase. The non core-collapsed clusters are still undergoing gravo-thermal contraction.

preprint2011arXiv

Angular Momentum Changes Due to Direct Impact Accretion in a Simplified Binary System

We model a circular mass-transferring binary system to calculate the exchange of angular momentum between stellar spins and the orbit due to direct impact of the mass transfer stream onto the surface of the accretor. We simulate mass transfer by calculating the ballistic motion of a point mass ejected from the $L_1$ point of the donor star, conserving the total linear and angular momentum of the system, and treating the stars as uniform density spheres with main sequence radii determined by their masses. We show that, contrary to previous assumptions in the literature, direct impact does not always act as a sink of orbital angular momentum and may in fact increase it by facilitating the transfer of angular momentum from the spin of the donor to the orbit. Here, we show an example of the exchange of angular momentum, as well as a measure of the orbital angular momentum changes for a variety of binary star systems with main sequence components.

preprint2011arXiv

Hot Jupiters from Secular Planet--Planet Interactions

About 25 per cent of `hot Jupiters' (extrasolar Jovian-mass planets with close-in orbits) are actually orbiting counter to the spin direction of the star. Perturbations from a distant binary star companion can produce high inclinations, but cannot explain orbits that are retrograde with respect to the total angular momentum of the system. Such orbits in a stellar context can be produced through secular (that is, long term) perturbations in hierarchical triple-star systems. Here we report a similar analysis of planetary bodies, including both octupole-order effects and tidal friction, and find that we can produce hot Jupiters in orbits that are retrograde with respect to the total angular momentum. With distant stellar mass perturbers, such an outcome is not possible. With planetary perturbers, the inner orbit's angular momentum component parallel to the total angular momentum need not be constant. In fact, as we show here, it can even change sign, leading to a retrograde orbit. A brief excursion to very high eccentricity during the chaotic evolution of the inner orbit allows planet-star tidal interactions to rapidly circularize that orbit, decoupling the planets and forming a retrograde hot Jupiter.

preprint2011arXiv

Twin Binaries: Studies of Stability, Mass Transfer, and Coalescence

Motivated by suggestions that binaries with almost equal-mass components ("twins") play an important role in the formation of double neutron stars and may be rather abundant among binaries, we study the stability of synchronized close and contact binaries with identical components in circular orbits. In particular, we investigate the dependency of the innermost stable circular orbit on the core mass, and we study the coalescence of the binary that occurs at smaller separations. For twin binaries composed of convective main-sequence stars, subgiants, or giants with low mass cores (M_c <~0.15M, where M is the mass of a component), a secular instability is reached during the contact phase, accompanied by a dynamical mass transfer instability at the same or at a slightly smaller orbital separation. Binaries that come inside this instability limit transfer mass gradually from one component to the other and then coalesce quickly as mass is lost through the outer Lagrangian points. For twin giant binaries with moderate to massive cores (M_c >~0.15M), we find that stable contact configurations exist at all separations down to the Roche limit, when mass shedding through the outer Lagrangian points triggers a coalescence of the envelopes and leaves the cores orbiting in a central tight binary. In addition to the formation of binary neutron stars, we also discuss the implications of our results for the production of planetary nebulae with double degenerate central binaries.

preprint2010arXiv

How planet-planet scattering can create high-inclination as well as long-period orbits

Recent observations have revealed two new classes of planetary orbits. Rossiter- Mclaughlin (RM) measurements have revealed hot Jupiters in high-obliquity orbits. In addition, direct-imaging has discovered giant planets at large (~ 100 AU) separations via direct-imaging technique. Simple-minded disk-migration scenarios are inconsistent with the high-inclination (and even retrograde) orbits as seen in recent RM measurements. Furthermore, forming giant planets at large semi-major axis (a) may be challenging in the core-accretion paradigm. We perform many N-body simulations to explore the two above-mentioned orbital architectures. Planet-planet scattering in a multi-planet system can naturally excite orbital inclinations. Planets can also get scattered to large distances. Large-a planetary orbits created from planet-planet scattering are expected to have high eccentricities (e). Theoretical models predict that the observed long-period planets, such as Fomalhaut-b have moderate e \approx 0.3. Interestingly, these are also in systems with disks. We find that if a massive-enough outer disk is present, a scattered planet may be circularized at large a via dynamical friction from the disk and repeated scattering of the disk particles.

preprint2010arXiv

Pervasive Orbital Eccentricities Dictate the Habitability of Extrasolar Earths

The long-term habitability of Earth-like planets requires low orbital eccentricities. A secular perturbation from a distant stellar companion is a very important mechanism in exciting planetary eccentricities, as many of the extrasolar planetary systems are associated with stellar companions. Although the orbital evolution of an Earth-like planet in a stellar binary is well understood, the effect of a binary perturbation to a more realistic system containing additional gas giant planets has been very little studied. Here we provide analytic criteria confirmed by a large ensemble of numerical integrations that identify the initial orbital parameters leading to eccentric orbits. We show that an extra-solar earth is likely to experience a broad range of orbital evolution dictated by the location of a gas-giant planet, necessitating more focused studies on the effect of eccentricity on the potential for life.

preprint2010arXiv

Tidal Evolution of Close-in Planets

Recent discoveries of several transiting planets with clearly non-zero eccentricities and some large inclinations started changing the simple picture of close-in planets having circular and well-aligned orbits. Two major scenarios to form such planets are planet migration in a disk, and planet--planet interactions combined with tidal dissipation. The former scenario can naturally produce a circular and low-obliquity orbit, while the latter implicitly assumes an initially highly eccentric and possibly high-obliquity orbit, which are then circularized and aligned via tidal dissipation. We investigate the tidal evolution of transiting planets on eccentric orbits. We show that the current and future orbital evolution of these systems is likely dominated by tidal dissipation, and not by a more distant companion. Although most of these close-in planets experience orbital decay all the way to the Roche limit, there are two characteristic evolution paths for them, depending on the relative efficiency of tidal dissipation inside the star and the planet. We point out that the current observations may be consistent with one of them. Our results suggest that at least some of the close-in planets with non-zero orbital eccentricity may have been formed by tidally circularizing an initially eccentric orbit. We also find that even when the stellar spin-orbit misalignment is observed to be small at present, some systems could have had a highly misaligned orbit in the past. Finally, we also re-examine the recent claim by Levrard et. al., who found that all orbital and spin parameters evolve on a similar timescale to orbital decay.

preprint2010arXiv

Unstable Planetary Systems Emerging Out Of Gas Disks

The discovery of over 400 extrasolar planets allows us to statistically test our understanding of formation and dynamics of planetary systems via numerical simulations. Traditional N-body simulations of multiple-planet systems without gas disks have successfully reproduced the eccentricity (e) distribution of the observed systems, by assuming that the planetary systems are relatively closely packed when the gas disk dissipates, so that they become dynamically unstable within the stellar lifetime. However, such studies cannot explain the small semi-major axes (a) of extrasolar planetary systems, if planets are formed, as the standard planet formation theory suggests, beyond the ice line. In this paper, we numerically study the evolution of three-planet systems in dissipating gas disks, and constrain the initial conditions that reproduce the observed semi-major axis and eccentricity distributions simultaneously. We adopt the initial conditions that are motivated by the standard planet formation theory, and self-consistently simulate the disk evolution, and planet migration by using a hybrid N-body and 1D gas disk code. We also take account of eccentricity damping, and investigate the effect of saturation of corotation resonances on the evolution of planetary systems. We find that the semi-major axis distribution is largely determined in a gas disk, while the eccentricity distribution is determined after the disk dissipation. We also find that there may be an optimum disk mass which leads to the observed a-e distribution. Our simulations generate a larger fraction of planetary systems trapped in mean-motion resonances (MMRs) than the observations, indicating that the disk's perturbation to the planetary orbits may be important to explain the observed rate of MMRs. We also find much lower occurrence of planets on retrograde orbits than the current observations of close-in planets suggest.

preprint2009arXiv

Monte Carlo Simulations of Globular Cluster Evolution. V. Binary Stellar Evolution

We study the dynamical evolution of globular clusters containing primordial binaries, including full single and binary stellar evolution using our Monte Carlo cluster evolution code updated with an adaptation of the single and binary stellar evolution codes SSE/BSE from Hurley et. al (2000, 2002). We describe the modifications we have made to the code. We present several test calculations and comparisons with existing studies to illustrate the validity of the code. We show that our code finds very good agreement with direct N-body simulations including primordial binaries and stellar evolution. We find significant differences in the evolution of the global properties of the simulated clusters using stellar evolution compared to simulations without any stellar evolution. In particular, we find that the mass loss from stellar evolution acts as a significant energy production channel simply by reducing the total gravitational binding energy and can significantly prolong the initial core contraction phase before reaching the binary-burning quasi steady state of the cluster evolution as noticed in Paper IV. We simulate a large grid of clusters varying the initial cluster mass, binary fraction, and concentration and compare properties of the simulated clusters with those of the observed Galactic globular clusters (GGCs). We find that our simulated cluster properties agree well with the observed GGC properties. We explore in some detail qualitatively different clusters in different phases of their evolution, and construct synthetic Hertzprung-Russell diagrams for these clusters.

preprint2007arXiv

On the rarity of double black hole binaries: consequences for gravitational-wave detection

Double black hole binaries are among the most important sources of gravitational radiation for ground-based detectors such as LIGO or VIRGO. Even if formed with lower efficiency than double neutron star binaries, they could dominate the predicted detection rates, since black holes are more massive than neutron stars and therefore could be detected at greater distances. Here we discuss an evolutionary process that can very significantly limit the formation of close double black hole binaries: the vast majority of their potential progenitors undergo a common envelope (CE) phase while the donor, one of the massive binary components, is evolving through the Hertzsprung gap. Our latest theoretical understanding of the CE process suggests that this will probably lead to a merger, inhibiting double black hole formation. Barring uncertainties in the physics of CE evolution, we use population synthesis calculations, and find that the corresponding reduction in the merger rate of double black holes formed in galactic fields is so great (by ~500) that their contribution to inspiral detection rates for ground-based detectors could become relatively small (1 in 10) compared to double neutron star binaries. A similar process also reduces the merger rates for double neutron stars, by factor of ~5, eliminating most of the previously predicted ultracompact NS-NS systems. Our predicted detection rates for Advanced LIGO are now much lower for double black holes (~2/yr), but are still quite high for double neutron stars (~20/yr). If double black holes were found to be dominant in the detected inspiral signals, this could indicate that they mainly originate from dense star clusters (not included here) or that our theoretical understanding of the CE phase requires significant revision.

preprint2006arXiv

Monte Carlo Simulations of Globular Cluster Evolution. IV. Direct Integration of Strong Interactions

We study the dynamical evolution of globular clusters containing populations of primordial binaries, using our newly updated Monte Carlo cluster evolution code with the inclusion of direct integration of binary scattering interactions. We describe the modifications we have made to the code, as well as improvements we have made to the core Monte Carlo method. We present several test calculations to verify the validity of the new code, and perform many comparisons with previous analytical and numerical work in the literature. We simulate the evolution of a large grid of models, with a wide range of initial cluster profiles, and with binary fractions ranging from 0 to 1, and compare with observations of Galactic globular clusters. We find that our code yields very good agreement with direct N-body simulations of clusters with primordial binaries, but yields some results that differ significantly from other approximate methods. Notably, the direct integration of binary interactions reduces their energy generation rate relative to the simple recipes used in Paper III, and yields smaller core radii. Our results for the structural parameters of clusters during the binary-burning phase are now in the tail of the range of parameters for observed clusters, implying that either clusters are born significantly more or less centrally concentrated than has been previously considered, or that there are additional physical processes beyond two-body relaxation and binary interactions that affect the structural characteristics of clusters.

preprint2005arXiv

Binary Mergers and Growth of Black Holes in Dense Star Clusters

We model the dynamical evolution of primordial black holes (BHs) in dense star clusters using a simplified treatment of stellar dynamics in which the BHs are assumed to remain concentrated in an inner core, completely decoupled from the background stars. Dynamical interactions involving BH binaries are computed exactly and are generated according to a Monte Carlo prescription. Recoil and ejections lead to complete evaporation of the BH core on a timescale ~10^9 yr for typical globular cluster parameters. Orbital decay driven by gravitational radiation can make binaries merge and, in some cases, successive mergers can lead to significant BH growth. Our highly simplified treatment of the cluster dynamics allows us to study a large number of models and to compute statistical distributions of outcomes, such as the probability of massive BH growth and retention in a cluster. We find that, in most models, there is a significant probability (~20-80%) of BH growth with final masses > 100 M_{\sun}. In one case, a final BH formed with mass ~ 620 M_{\sun}. However, if the typical merger recoil speed (due to asymmetric emission of gravitational radiation) significantly exceeds the cluster escape speed, no growth ever occurs. Independent of the recoil speed, we find that BH-BH mergers enhanced by dynamical interactions in cluster cores present an important source of gravitational waves for ground--based laser interferometers. Under optimistic conditions, the total rate of detections by Advanced LIGO, could be as high as a few tens of events per year from inspiraling BHs from clusters.

preprint2004arXiv

Collisions Between Single Stars in Dense Clusters: Runaway Formation of a Massive Object

Using Monte Carlo codes, we follow the collisional evolution of clusters in a variety of scenarios. We consider the conditions under which a cluster of main sequence stars may undergo rapid core collapse due to mass segregation, thus entering a phase of runaway collisions, forming a very massive star (VMS, M >= 1000 Msun) through repeated collisions between single stars. Although collisional mass loss is accounted for realistically, we find that a VMS forms even in proto-galactic nuclei models with a high velocity dispersion (many 100 km/s). Such a VMS may be a progenitor for an intermediate-mass black hole (M >= 100 Msun). In contrast, in galactic nuclei hosting a central massive black hole, collisions are found to be disruptive. The stars which are subject to collisions are progressively ground down by high-velocity collisions and a merger sequence appears impossible.

preprint2004arXiv

Formation of Massive Black Holes in Dense Star Clusters. I. Mass Segregation and Core Collapse

We study the early dynamical evolution of young, dense star clusters using Monte Carlo simulations for systems with up to N~10^7 stars. Rapid mass segregation of massive main-sequence stars and the development of the Spitzer instability can drive these systems to core collapse in a small fraction of the initial half-mass relaxation time. If the core collapse time is less than the lifetime of the massive stars, all stars in the collapsing core may then undergo a runaway collision process leading to the formation of a massive black hole. Here we study in detail the first step in this process, up to the occurrence of core collapse. We have performed about 100 simulations for clusters with a wide variety of initial conditions, varying systematically the cluster density profile, stellar IMF, and number of stars. We also considered the effects of initial mass segregation and stellar evolution mass loss. Our results show that, for clusters with a moderate initial central concentration and any realistic IMF, the ratio of core collapse time to initial half-mass relaxation time is typically ~0.1, in agreement with the value previously found by direct N-body simulations for much smaller systems. Models with even higher central concentration initially, or with initial mass segregation (from star formation) have even shorter core-collapse times. Remarkably, we find that, for all realistic initial conditions, the mass of the collapsing core is always close to ~10^-3 of the total cluster mass, very similar to the observed correlation between central black hole mass and total cluster mass in a variety of environments. We discuss the implications of our results for the formation of intermediate-mass black holes in globular clusters and super star clusters, ultraluminous X-ray sources, and seed black holes in proto-galactic nuclei.

preprint2001arXiv

Binaries and Globular Cluster Dynamics

We summarize the results of recent theoretical work on the dynamical evolution of globular clusters containing primordial binaries. Even a very small initial binary fraction (e.g., 10%) can play a key role in supporting a cluster against gravothermal collapse for many relaxation times. Inelastic encounters between binaries and single stars or other binaries provide a very significant energy source for the cluster. These dynamical interactions also lead to the production of large numbers of exotic systems such as ultracompact X-ray binaries, recycled radio pulsars, double degenerate systems, and blue stragglers. Our work is based on a new parallel supercomputer code implementing Henon's Monte Carlo method for simulating the dynamical evolution of dense stellar systems in the Fokker-Planck approximation. This new code allows us to calculate very accurately the evolution of a cluster containing a realistic number of stars (N ~ 10^5 - 10^6) in typically a few hours to a few days of computing time. The discrete, star-by-star representation of the cluster in the simulation makes it possible to treat naturally a number of important processes, including single and binary star evolution, all dynamical interactions of single stars and binaries, and tidal interactions with the Galaxy.

preprint2001arXiv

Monte Carlo Simulations of Globular Cluster Evolution - II. Mass Spectra, Stellar Evolution and Lifetimes in the Galaxy

We study the dynamical evolution of globular clusters using our new 2-D Monte Carlo code, and we calculate the lifetimes of clusters in the Galactic environment. We include the effects of a mass spectrum, mass loss in the Galactic tidal field, and stellar evolution. We consider initial King models containing N = 10^5 - 3x10^5 stars, and follow the evolution up to core collapse, or disruption, whichever occurs first. We find that the lifetimes of our models are significantly longer than those obtained using 1-D Fokker-Planck (F-P) methods. We also find that our results are in very good agreement with recent 2-D F-P calculations, for a wide range of initial conditions. Our results show that the direct mass loss due to stellar evolution can significantly accelerate the mass loss through the tidal boundary, causing most clusters with a low initial central concentration (Wo <~ 3) to disrupt quickly in the Galactic tidal field. Only clusters born with high initial central concentrations (Wo >~ 7) or steep initial mass functions are likely to survive to the present and undergo core collapse. We also study the orbital characteristics of escaping stars, and find that the velocity distribution of escaping stars in collapsing clusters looks significantly different from the distribution in disrupting clusters. We calculate the lifetime of a cluster on an eccentric orbit in the Galaxy, such that it fills its Roche lobe only at perigalacticon. We find that such an orbit can extend the lifetime by at most a factor of a few compared to a circular orbit in which the cluster fills its Roche lobe at all times.

preprint2000arXiv

Evolution of Stellar Collision Products in Globular Clusters - II. Off-axis Collision

We continue our exploration of collisionally merged stars in the blue straggler region of the color-magnitude diagram. We report the results of new SPH calculations of parabolic collisions between two main-sequence stars, with the initial structure and composition profiles of the parent stars having been determined from stellar evolution calculations. Parallelization of the SPH code has permitted much higher numerical resolution of the hydrodynamics. We also present evolutionary tracks for the resulting collision products, which emerge as rapidly rotating blue stragglers. The rotating collision products are brighter, bluer and remain on the main sequence longer than their non-rotating counterparts. In addition, they retain their rapid rotation rates throughout their main sequence lifetime. Rotationally-induced mixing strongly affects the evolution of the collision products, although it is not sufficient to mix the entire star. We discuss the implications of these results for studies of blue straggler populations in clusters. This work shows that off-axis collision products cannot become blue stragglers unless they lose a large fraction of their initial angular momentum. The mechanism for this loss is not apparent, although some possibilities are discussed.

preprint1999arXiv

Coalescing Binary Neutron Stars

Coalescing compact binaries with neutron star or black hole components provide the most promising sources of gravitational radiation for detection by the LIGO/VIRGO/GEO/TAMA laser interferometers now under construction. This fact has motivated several different theoretical studies of the inspiral and hydrodynamic merging of compact binaries. Analytic analyses of the inspiral waveforms have been performed in the Post-Newtonian approximation. Analytic and numerical treatments of the coalescence waveforms from binary neutron stars have been performed using Newtonian hydrodynamics and the quadrupole radiation approximation. Numerical simulations of coalescing black hole and neutron star binaries are also underway in full general relativity. Recent results from each of these approaches will be described and their virtues and limitations summarized.

preprint1996arXiv

The Effect of Encounters on the Eccentricity of Binaries in Clusters

We derive analytical expressions for the change in the orbital eccentricity of a binary following a distant encounter with a third star on a hyperbolic or parabolic orbit. To establish the accuracy of these expressions, we present detailed comparisons with the results of direct numerical integrations of the equations of motion for the three bodies. We treat with particular care the difficult case of a binary with zero initial eccentricity. In this case, we show that the eccentricity $δe$ induced by the encounter declines in general as a power-law, $δe\propto (a/\rp)^{5/2}$, where $a$ is the binary semi-major axis and $\rp$ is the periastron distance of the encounter. This power-law arises from the octupole-level secular perturbation of the binary. In contrast, non-secular quadrupole-level perturbations induce an eccentricity change that declines exponentially with $\rp$. These non-secular effects can become dominant at sufficiently small $\rp$, for a sufficiently high relative velocity, or for a sufficiently massive perturber. We also derive cross sections for eccentricity change and compare our results with those of previous studies based on numerical scattering experiments. Our results have important implications for a number of astrophysical problems including, in particular, the evolution of binary millisecond pulsars in globular clusters.

preprint1995arXiv

ON BLUE STRAGGLER FORMATION BY DIRECT COLLISIONS OF MAIN-SEQUENCE STARS

We report the results of new SPH calculations of parabolic collisions between main-sequence (MS) stars. The stars are assumed to be close to the MS turn-off point in a globular cluster and are therefore modeled as $n=3$, $Γ=5/3$ polytropes. We find that the high degree of central mass concentration in these stars has a profound effect on the hydrodynamics. In particular, very little hydrodynamic mixing occurs between the dense, helium-rich inner cores and the outer envelopes. As a result, and in contrast to what has been assumed in previous studies, blue stragglers formed by direct stellar collisions are not necessarily expected to have anomalously high helium abundances in their envelopes or to have their cores replenished with fresh hydrogen fuel.

preprint1995arXiv

THE MINIMUM MASS RATIO OF W URSAE MAJORIS BINARIES

The minimum mass ratio for tidal stability of a contact binary containing two unevolved main-sequence stars is calculated to be $\qmin\simeq0.09$ in the case of a mostly radiative primary, and is higher if an appreciable fraction of the mass lies in a convective envelope. At least one observed system, AW UMa, has a mass ratio just below this value ($q=0.075$), implying, if the system is stable, that the primary must be slightly evolved and must have a very shallow convective envelope. Contact binaries with mass ratios significantly below that of AW UMa should not be observed, since they are tidally unstable and quickly merge into a single, rapidly rotating object, on a timescale $\sim10^3$--$10^4\,$yr.

preprint1995arXiv

THE ORBITAL ECCENTRICITIES OF BINARY MILLISECOND PULSARS IN GLOBULAR CLUSTERS

Low-mass binary millisecond pulsars (LMBPs) are born with very small orbital eccentricities, typically of order $e_i\sim10^{-6}$--$10^{-3}$. In globular clusters, however, higher eccentricities $e_f\gg e_i$ can be induced by dynamical interactions with passing stars. Here we show that the cross section for this process is much larger than previously estimated. This is because, even for initially circular binaries, the induced eccentricity $e_f$ for an encounter with pericenter separation $r_{\rm p}$ beyond a few times the binary semi-major axis $a$ declines only as a power-law, $e_f\propto(r_{\rm p}/a)^ {-5/2}$, and {\it not\/} as an exponential. We find that {\it all\/} currently known \lmbps in clusters were probably affected by interactions, with their current eccentricities typically an order of magnitude or more greater than at birth.

Frederic A. Rasio

What is connected

Connect this record

See the researcher in context

Building this map preview

70 published item(s)

Compact Object Modeling in the Globular Cluster 47 Tucanae

Eccentric Mergers of Intermediate-Mass Black Holes from Evection Resonances in AGN Disks

Formation of Low-mass Black Holes and Single Millisecond Pulsars in Globular Clusters

Gravitational Microlensing Rates in Milky Way Globular Clusters

Hydrodynamics of Collisions and Close Encounters between Stellar Black Holes and Main-sequence Stars

Merger rates of intermediate-mass black hole binaries in nuclear star clusters

Repeated mergers, mass-gap black holes, and formation of intermediate-mass black holes in nuclear star clusters

Intermediate-mass Black Holes from High Massive-star Binary Fractions in Young Star Clusters

Modeling Dense Star Clusters in the Milky Way and Beyond with the Cluster Monte Carlo Code

The Observed Rate of Binary Black Hole Mergers can be Entirely Explained by Globular Clusters

A Dynamical Survey of Stellar-Mass Black Holes in 50 Milky Way Globular Clusters

COSMIC Variance in Binary Population Synthesis

Demographics of triple systems in dense star clusters

Gravitational-wave captures by intermediate-mass black holes in galactic nuclei

GW190412 as a Third-Generation Black Hole Merger from a Super Star Cluster

Merging Black Holes in the Low-mass and High-mass Gaps from 2+2 Quadruple Systems

Modeling Dense Star Clusters in the Milky Way and Beyond with the $\texttt{CMC}$ Cluster Catalog

On the Rate of Neutron Star Binary Mergers from Globular Clusters

The Role of "Black Hole Burning" in the Evolution of Dense Star Clusters

Binary Black Hole Mergers from Globular Clusters: Masses, Merger Rates, and the Impact of Stellar Evolution

Distinguishing Between Formation Channels for Binary Black Holes with LISA

Dynamical Formation of the GW150914 Binary Black Hole

Illuminating Black Hole Binary Formation Channels with Spins in Advanced LIGO

Merging black hole binaries in galactic nuclei: implications for advanced-LIGO detections

Million-Body Star Cluster Simulations: Comparisons between Monte Carlo and Direct $N$-body

Binary Black Hole Mergers from Globular Clusters: Implications for Advanced LIGO

Black hole mergers and blue stragglers from hierarchical triples formed in globular clusters

Formation of Black Hole Low-Mass X-ray Binaries in Hierarchical Triple Systems

Tidally-driven Roche-Lobe Overflow of Hot Jupiters with MESA

From Hot Jupiters to Super-Earths via Roche Lobe Overflow

Planets on the Edge

Sowing the seeds of massive black holes in small galaxies: Young clusters as the building blocks of Ultra-Compact-Dwarf Galaxies

The Dynamical Evolution of Stellar Black Holes in Globular Clusters

The Dynamics of the Multi-planet System Orbiting Kepler-56

Tidal Dissipation and Obliquity Evolution in Hot Jupiter Systems

Extreme orbital evolution from hierarchical secular coupling of two giant planets

Introducing CAFein, a New Computational Tool for Stellar Pulsations and Dynamic Tides

Mock Observations of Blue Stragglers in Globular Cluster Models

Retention of Stellar-Mass Black Holes in Globular Clusters

Secular Dynamics in Hierarchical Three-Body Systems

Stellar Collisions and Blue Straggler Stars in Dense Globular Clusters

A Parallel Monte Carlo Code for Simulating Collisional N-body Systems

Binary Neutron Star Mergers

Constraining Intermediate-Mass Black Holes in Globular Clusters

Formation of Massive Black Holes in Dense Star Clusters. II. IMF and Primordial Mass Segregation

Monte Carlo Simulations of Globular Cluster Evolution. VI. The Influence of an Intermediate Mass Black Hole

On the Formation of Hot Jupiters in Stellar Binaries

Planets in Open Clusters Detectable by Kepler

Understanding the Dynamical State of Globular Clusters: Core-Collapsed vs Non Core-Collapsed

Angular Momentum Changes Due to Direct Impact Accretion in a Simplified Binary System

Hot Jupiters from Secular Planet--Planet Interactions

Twin Binaries: Studies of Stability, Mass Transfer, and Coalescence

How planet-planet scattering can create high-inclination as well as long-period orbits

Pervasive Orbital Eccentricities Dictate the Habitability of Extrasolar Earths

Tidal Evolution of Close-in Planets

Unstable Planetary Systems Emerging Out Of Gas Disks

Monte Carlo Simulations of Globular Cluster Evolution. V. Binary Stellar Evolution

On the rarity of double black hole binaries: consequences for gravitational-wave detection

Monte Carlo Simulations of Globular Cluster Evolution. IV. Direct Integration of Strong Interactions

Binary Mergers and Growth of Black Holes in Dense Star Clusters

Collisions Between Single Stars in Dense Clusters: Runaway Formation of a Massive Object

Formation of Massive Black Holes in Dense Star Clusters. I. Mass Segregation and Core Collapse

Binaries and Globular Cluster Dynamics

Monte Carlo Simulations of Globular Cluster Evolution - II. Mass Spectra, Stellar Evolution and Lifetimes in the Galaxy

Evolution of Stellar Collision Products in Globular Clusters - II. Off-axis Collision

Coalescing Binary Neutron Stars

The Effect of Encounters on the Eccentricity of Binaries in Clusters

ON BLUE STRAGGLER FORMATION BY DIRECT COLLISIONS OF MAIN-SEQUENCE STARS

THE MINIMUM MASS RATIO OF W URSAE MAJORIS BINARIES

THE ORBITAL ECCENTRICITIES OF BINARY MILLISECOND PULSARS IN GLOBULAR CLUSTERS