Source author record

Sambaran Banerjee

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

20 published item(s)

preprint2022arXiv

Binary black hole mergers from young massive clusters in the pair-instability supernova mass gap

The recent discovery of the binary black hole (BBH) merger event GW190521, between two black holes (BHs) of $\approx100M_\odot$, and as well as other massive BBH merger events involving BHs within the pair-instability supernova (PSN) mass gap have sparked widespread debate on the origin of such extreme gravitational-wave (GW) events. In this study, I investigate whether dynamical interactions in young massive clusters (YMCs) serves as a viable scenario for assembling PSN-gap BBH mergers. To that end, I explore a grid of 40 new evolutionary models of a representative YMC of initial mass $M_{\rm cl}=7.5\times10^4M_\odot$ ($N\approx1.28\times10^5$) and size $r_h=2$ pc, with all BH progenitor stars being initially in primordial binaries. All cluster models are evolved with the direct, relativistic N-body code NBODY7 incorporating up to date remnant formation, BH natal spin, and general-relativistic (GR) merger recoil schemes. The BBH mergers from these model cluster computations agree well with the masses and effective spin parameters of the GW events in the latest GW transient catalogue (GWTC). In particular, GW190521-like, i.e., $\approx200M_\odot$, low aligned spin events are produced via dynamical merger among BHs derived from star-star merger products. GW190403-like, i.e., PSN-gap, highly asymmetric and aligned events result from mergers involving BHs that are spun up via matter accretion or binary interaction. The present YMC models yield a present day, intrinsic merger rate density of $0-3.8\times10^{-2}{{\rm~yr}^{-1}{\rm Gpc}^{-3}}$ for GW190521-type events. They produce GW190403-like events at a rate within $0-1.6\times10^{-1}{{\rm~yr}^{-1}{\rm Gpc}^{-3}}$ and their total BBH-merger yield within the PSN gap is $0-8.4\times10^{-1}{{\rm~yr}^{-1}{\rm Gpc}^{-3}}$.

preprint2021arXiv

Black hole mergers in compact star clusters and massive black hole formation beyond the mass-gap

We present direct N-body simulations, carried out with Nbody6++GPU, of young and compact low metallicity star clusters with $1.1\times 10^5$ stars, a velocity dispersion of $\sim$ 10 $\mathrm{km\,s^{-1}}$, a half mass radius $R_h=0.6$ pc, and a binary fraction of $10\%$ including updated evolution models for stellar winds and pair-instability supernovae (PISNe). Within the first tens of megayears of evolution, each cluster hosts several black hole (BH) merger events which nearly cover the complete mass range of primary and secondary BH masses for current LIGO/Virgo/Kagra gravitational wave detections. The importance of gravitational recoil is estimated statistically. We present several possible formation paths of massive BHs above the assumed lower PISNe mass-gap limit ($45 M_\odot$) into the intermediate-mass BH (IMBH) regime ($> 100 M_\odot$) which include collisions of stars and BHs as well as the direct collapse of stellar merger remnants with low mass cores. The stellar evolution updates result in the early formation of higher mass stellar BHs than for the previous model. The resulting higher collision rates with massive stars support the rapid formation of massive BHs. For models assuming a high accretion efficiency for star-BH mergers, we present a first-generation formation scenario for GW190521-like events, a merger of two BHs in the PISN mass-gap, which is dominated by star-BH mergers. This IMBH formation path is independent of gravitational recoil and therefore conceivable in dense stellar systems with low escape velocities. One simulated cluster even forms an IMBH binary (153$M_\odot$,173$M_\odot$) which is expected to merge within a Hubble time.

preprint2021arXiv

Merger rate density of stellar-mass binary black holes from young massive clusters, open clusters, and isolated binaries: comparisons with LIGO-Virgo-KAGRA results

I investigate the roles of cluster dynamics and massive binary evolution in producing stellar-remnant binary black hole (BBH) mergers over the cosmic time. To that end, dynamical BBH mergers are obtained from long-term direct N-body evolutionary models of $\sim10^4M_\odot$, pc-scale young massive clusters (YMC) evolving into moderate-mass open clusters (OC). Fast evolutionary models of massive isolated binaries (IB) yield BBHs from binary evolution. Population synthesis in a Model Universe is then performed, taking into account observed cosmic star-formation and enrichment histories, to obtain BBH-merger yields from these two channels observable at the present day and over cosmic time. The merging BBH populations from the two channels are combined by applying a proof-of-concept Bayesian regression chain, taking into account observed differential intrinsic BBH merger rate densities from the second gravitational-wave transient catalogue (GWTC-2). The analysis estimates an OB-star binary fraction of $f_{\rm Obin}\gtrsim90$% and a YMC formation efficiency of $f_{\rm YMC}\sim10^{-2}$, being consistent with recent optical observations and large scale structure formation simulations. The corresponding combined Model Universe present-day, differential intrinsic BBH merger rate density and the cosmic evolution of BBH merger rate density both agree well with those from GWTC-2. The analysis also suggests that despite significant 'dynamical mixing' at low redshifts, BBH mergers at high redshifts ($z_{\rm event}\gtrsim1$) could still be predominantly determined by binary-evolution physics. Caveats in the present approach and future improvements are discussed.

preprint2020arXiv

BSE versus StarTrack: implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

The masses of stellar-remnant black holes (BH), as a result of their formation via massive single- and binary-stellar evolution, is of high interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star (BNS) mergers. Here we present new developments in the N-body evolution program NBODY7 in regards to its stellar-remnant formation and related schemes. We demonstrate that the newly-implemented stellar-wind and remnant-formation schemes in the NBODY7 code's BSE sector, such as the 'rapid' and the 'delayed' supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produces neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main sequence (ZAMS) mass and metallicity, with those from the StarTrack population-synthesis program. We also demonstrate the new implementations of various natal-kick mechanisms on NSs and BHs such as the 'convection-asymmetry-driven', 'collapse-asymmetry-driven', and 'neutrino-emission-driven' kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. We find that the SN material fallback causes the convection-asymmetry kick to effectively retain similar number and mass of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry kick would cause nearly all BHs to retain in clusters irrespective of remnant formation model and metallicity, whereas the inference of a large number of BHs in GCs would potentially rule out the neutrino-driven kick mechanism. Pre-SN mergers of massive primordial binaries would cause BH masses to deviate from the single-star ZAMS mass-remnant mass relation. Such mergers, at low metallicities, can produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster.

preprint2020arXiv

Formation of low-spinning 100 Msun black holes

It is speculated that a merger of two massive stellar-origin BHs in a dense stellar environment may lead to the formation of a massive BH in the pair-instability mass gap (50-135 Msun). Such a merger-formed BH is expected to typically have a high spin (a=0.7). If such a massive BH acquires another BH it may lead to another merger detectable by LIGO/Virgo in gravitational waves. Acquiring a companion may be hindered by gravitational-wave kick/recoil, which accompanies the first merger and may quickly remove the massive BH from its parent globular or nuclear cluster. We test whether it is possible for a massive merger-formed BH in the pair-instability gap to be retained in its parent cluster and have low spin. Such a BH would be indistinguishable from a primordial BH. We employed results from numerical relativity calculations of black hole mergers to explore the range of gravitational-wave recoil velocities for various combinations of merging BH masses and spins. We compared merger-formed massive BH speeds with typical escape velocities from globular and nuclear clusters. We show that a globular cluster is highly unlikely to form and retain a 100 Msun BH if the spin of the BH is low (a<0.3) as such BHs acquire high recoil speeds (>200 km/s) that exceed typical escape speeds from globular clusters (50 km/s). However, a very low-spinning (a=0.1) and massive (100 Msun) BH could be formed and retained in a galactic nuclear star cluster. Even though such massive merger-formed BHs with such low spins acquire high speeds during formation (400 km/s), they may avoid ejection since massive nuclear clusters have high escape velocities (300-500 km/s). A future detection of a massive BH in the pair-instability mass gap with low spin would therefore not be proof of the existence of primordial BHs, which are sometimes claimed to have low spins and arbitrarily high masses.

preprint2020arXiv

Stellar-mass black holes in young massive and open stellar clusters IV: updated stellar-evolutionary and black hole spin models and comparisons with the LIGO-Virgo O1/O2 merger-event data

I present a set of long-term, direct, relativistic many-body computations of model dense stellar clusters with up-to-date stellar-evolutionary, supernova (SN), and remnant natal-kick models, including pair instability and pulsation pair instability supernova (PSN and PPSN), using an updated version of NBODY7 N-body simulation program. The N-body model also includes stellar evolution-based natal spins of BHs and treatments of binary black hole (BBH) mergers based on numerical relativity. These, for the first time in a direct N-body simulation, allow for second-generation BBH mergers. The set of 65 evolutionary models have initial masses $10^4M_\odot-10^5M_\odot$, sizes 1 pc-3 pc, metallicity $0.0001-0.02$, with the massive stars in primordial binaries and they represent young massive clusters (YMC) and moderately massive open clusters (OC). Such models produce dynamically-paired BBH mergers that agree well with the observed masses, mass ratios, effective spin parameters, and final spins of the LVC O1/O2 merger events, provided BHs are born with low or no spin but spin up after undergoing a BBH merger or matter accretion onto it. In particular, the distinctly higher mass, effective spin parameter, and final spin of GW170729 merger event is naturally reproduced, as also the mass asymmetry of the O3 event GW190412. The computed models produce massive, $\sim100M_\odot$ BBH mergers with primary mass within the 'PSN gap' and also yield mergers involving remnants in the 'mass gap'. They also suggest that YMCs and OCs produce persistent, Local-Universe GW sources detectable by LISA. Such clusters are also capable of producing eccentric LIGO-Virgo mergers.

preprint2016arXiv

Blue straggler formation at core collapse

Among the most striking feature of blue straggler stars (BSS) is the presence of multiple sequences of BSSs in the colour-magnitude diagrams (CMDs) of several globular clusters. It is often envisaged that such a multiple BSS sequence would arise due a recent core collapse of the host cluster, triggering a number of stellar collisions and binary mass transfers simultaneously over a brief episode of time. Here we examine this scenario using direct N-body computations of moderately massive star clusters (of order 10^4 Msun ). As a preliminary attempt, these models are initiated with approx. 8-10 Gyr old stellar population and King profiles of high concentrations, being "tuned" to undergo core collapse quickly. BSSs are indeed found to form in a "burst" at the onset of the core collapse and several of such BS-bursts occur during the post-core-collapse phase. In those models that include a few percent primordial binaries, both collisional and binary BSSs form after the onset of the (near) core-collapse. However, there is as such no clear discrimination between the two types of BSSs in the corresponding computed CMDs. We note that this may be due to the less number of BSSs formed in these less massive models than that in actual globular clusters.

preprint2015arXiv

How can young massive clusters reach their present-day sizes?

The classic question that how young massive star clusters attain their shapes and sizes, as we find them today, remains to be a challenge. Both observational and computational studies of star-forming massive molecular gas clouds infer that massive cluster formation is primarily triggered along the small-scale ($\lesssim0.3$ pc) filamentary substructures within the clouds. The present study is intended to investigate the possible ways in which a filament-like-compact, massive star cluster (effective radius 0.1-0.3 pc) can expand $\gtrsim10$ times, still remaining massive enough ($\gtrsim10^4 M_\odot$), to become a young massive star cluster, as we observe today. To that end, model massive clusters (of initially $10^4 M_\odot-10^5 M_\odot$) are evolved using Sverre Aarseth's state-of-the-art N-body code NBODY7. All the computed clusters expand with time, whose sizes (effective radii) are compared with those observed for young massive clusters, of age $\lesssim100$ Myr, in the Milky Way and other nearby galaxies. It is found that beginning from the above compact sizes, a star cluster cannot expand by its own, i.e., due to two-body relaxation, stellar-evolutionary mass loss, dynamical heating by primordial binaries and stellar-mass black holes, up to the observed sizes of young massive clusters; they always remain much more compact compared to the observed ones. This calls for additional mechanisms that can boost the expansion of a massive cluster after its assembly. Using further N-body calculations, it is shown that a substantial residual gas expulsion, with $\approx30$% star formation efficiency, can indeed swell the newborn embedded cluster adequately. The limitations of the present calculations and their consequences are discussed.

preprint2015arXiv

Initial conditions of formation of starburst clusters: constraints from stellar dynamics

How starburst clusters form out of molecular clouds is still an open question. In this article, I highlight some of the key constraints in this regard, that one can get from the dynamical evolutionary properties of dense stellar systems. I particularly focus on secular expansion of massive star clusters and hierarchical merging of sub-clusters, and discuss their implications vis-a-vis the observed properties of young massive clusters. The analysis suggests that residual gas expulsion is necessary for shaping these clusters as we see them today, irrespective of their monolithic or hierarchical mode of formation.

preprint2015arXiv

Possible smoking-gun evidence for initial mass segregation in re-virialized post-gas expulsion globular clusters

We perform a series of direct $N$-body calculations to investigate the effect of residual gas expulsion from the gas-embedded progenitors of present-day globular clusters (GCs) on the stellar mass function (MF). Our models start either tidally filling or underfilling, and either with or without primordial mass segregation. We cover 100 Myr of the evolution of modeled clusters and show that the expulsion of residual gas from initially mass-segregated clusters leads to a significantly shallower slope of the stellar MF in the low- ($m\leq 0.50 M_\odot$) and intermediate-mass ($\simeq 0.50-0.85 M_\odot$) regime. Therefore, the imprint of residual gas expulsion and primordial mass segregation might be visible in the present-day MF. We find that the strength of the external tidal field, as an essential parameter, influences the degree of flattening, such that a primordially mass-segregated tidally-filling cluster with $r_h/r_t\geq 0.1$ shows a strongly depleted MF in the intermediate stellar mass range. Therefore, the shape of the present-day stellar MF in this mass range probes the birth place of clusters in the Galactic environment. We furthermore find that this flattening agrees with the observed correlation between the concentration of a cluster and its MF slope, as found by de Marchi et al.. We show that if the expansion through the residual gas expulsion in primordial mass segregated clusters is the reason for this correlation then GCs most probably formed in strongly fluctuating local tidal fields in the early proto-Milky Way potential, supporting the recent conclusion by Marks \& Kroupa.

preprint2014arXiv

A perfect starburst cluster made in one go: the NGC 3603 young cluster

Understanding how distinct, near-spherical gas-free clusters of very young, massive stars shape out of vast, complex clouds of molecular hydrogen is one of the biggest challenges in astrophysics. A popular thought dictates that a single gas cloud fragments into many new-born stars which, in turn, energize and rapidly expel the residual gas to form a gas-free cluster. This study demonstrates that the above classical paradigm remarkably reproduces the well-observed central, young cluster (HD 97950) of the Galactic NGC 3603 star-forming region, in particular, its shape, internal motion and the mass distribution of stars, naturally and consistently follow from a single model calculation. Remarkably, the same parameters (star formation efficiency, gas expulsion time scale and delay) reproduce HD 97950 as were found to reproduce the Orion Nebula Cluster, Pleiades and R136. The present results thereby provide intriguing evidences of formation of star clusters through single-starburst events followed by significant residual gas expulsion.

preprint2014arXiv

The formation of NGC 3603 young starburst cluster: "prompt" hierarchical assembly or monolithic starburst?

The formation of very young massive clusters or "starburst" clusters is currently one of the most widely debated topic in astronomy. The classical notion dictates that a star cluster is formed in-situ in a dense molecular gas clump followed by a substantial residual gas expulsion. On the other hand, based on the observed morphologies of many young stellar associations, a hierarchical formation scenario is alternatively suggested. A very young (age $\approx$ 1 Myr), massive ($>10^4M_\odot$) star cluster like the Galactic NGC 3603 young cluster (HD 97950) is an appropriate testbed for distinguishing between such "monolithic" and "hierarchical" formation scenarios. A recent study by Banerjee and Kroupa (2014) demonstrates that the monolithic scenario remarkably reproduces the HD 97950 cluster. In the present work, we explore the possibility of the formation of the above cluster via hierarchical assembly of subclusters. These subclusters are initially distributed over a wide range of spatial volumes and have various modes of sub-clustering. Unlike the monolithic model that reproduces HD 97950 very well, the same is found to be prohibitive with hierarchical assembly alone. Only those systems which assemble promptly into a single cluster (in $\lesssim1$ Myr) from a close separation (all within $\lesssim2$ pc) could match the observed density profile of HD 97950 after $\approx70$% gas dispersal. These results therefore suggest that NGC 3603 young cluster has formed essentially monolithically, i.e., either in-situ or via a prompt assembly, followed by a substantial residual gas expulsion. Both scenarios are consistent with the inferred young age and the small age spread of this cluster. Future observations of molecular cloud filaments with ALMA and proper motion measurements of young clusters with Gaia will provide more direct tests of such birth environments. [Abridged]

preprint2013arXiv

Did the infant R136 and NGC 3603 clusters undergo residual gas expulsion?

Based on kinematic data observed for very young, massive clusters that appear to be in dynamical equilibrium, it has recently been argued that such young systems set examples where the early residual gas-expulsion did not happen or had no dynamical effect. The intriguing scenario of a star cluster forming through a single starburst has thereby been challenged. Choosing the case of the R136 cluster of the Large Magellanic Cloud, the most cited one in this context, we perform direct N-body computations that mimic the early evolution of this cluster including the gas-removal phase (on a thermal timescale). Our calculations show that under plausible initial conditions as consistent from observational data, a large fraction (> 60%) of a gas-expelled, expanding R136-like cluster is bound to regain dynamical equilibrium by its current age. Therefore, the recent measurements of velocity dispersion in the inner regions of R136, that indicate that the cluster is in dynamical equilibrium, are consistent with an earlier substantial gas expulsion of R136 followed by a rapid re-virialization (in approx. 1 Myr). Additionally, we find that the less massive Galactic NGC 3603 Young Cluster (NYC), with a substantially longer re-virialization time, is likely to be found deviated from dynamical equilibrium at its present age (approx. 1 Myr). The recently obtained stellar proper motions in the central part of the NYC indeed suggest this and are consistent with the computed models. This work significantly extends previous models of the Orion Nebula Cluster which already demonstrated that the re-virialization time of young post-gas-expulsion clusters decreases with increasing pre-expulsion density.

preprint2013arXiv

R144 : a very massive binary likely ejected from R136 through a binary-binary encounter

R144 is a recently confirmed very massive, spectroscopic binary which appears isolated from the core of the massive young star cluster R136. The dynamical ejection hypothesis as an origin for its location is claimed improbable by Sana et al. due to its binary nature and high mass. We demonstrate here by means of direct N-body calculations that a very massive binary system can be readily dynamically ejected from a R136-like cluster, through a close encounter with a very massive system. One out of four N-body cluster models produces a dynamically ejected very massive binary system with a mass comparable to R144. The system has a system mass of $\approx$ 355 Msun and is located at 36.8 pc from the centre of its parent cluster, moving away from the cluster with a velocity of 57 km/s at 2 Myr as a result of a binary-binary interaction. This implies that R144 could have been ejected from R136 through a strong encounter with an other massive binary or single star. In addition, we discuss all massive binaries and single stars which are ejected dynamically from their parent cluster in the N-body models.

preprint2012arXiv

On the true shape of the upper end of the stellar initial mass function: The case of R136

Context: The shape of the stellar initial mass function (IMF) of a star cluster near its upper mass limit is a focal topic of investigation as it determines the high mass stellar content and hence the dynamics of the cluster at its embedded phase as well as during its young gas-free phase. The massive stellar content of a young cluster, however, can be substantially modified due to the dynamical ejections of the massive stars so that the present-day high-mass stellar mass function (hereafter MF) can be different than that with which the cluster is born. Aims: In the present study, we provide a preliminary estimate of this evolution of the high-mass IMF of a young cluster due to early ejections of massive stars, using the Large Magellanic Cloud massive, young cluster R136 as an example. Methods: To that end, we utilize the results of the state-of-the-art calculations by Banerjee, Kroupa & Oh (2012) comprising direct N-body computations of realistic, binary-rich, mass-segregated models of R136. In particular, these calculations provide the ejection fraction of stars as a function of stellar mass. Results: We find that if the measured IMF of R136 is granted to be canonical, as observations indicate, then the "true" high-mass IMF of R136 at its birth must be at least moderately top-heavy when corrected for the dynamical escape of massive stars. Conclusions: The top-heaviness of the true high-mass IMF over the observationally determined one is a general feature of massive, young clusters where the dynamical ejection of massive stars is efficient. We discuss its implications and possible improvements over our current estimate.

preprint2012arXiv

The emergence of super-canonical stars in R136-type star-burst clusters

[abridged] Among the most remarkable features of the stellar population of R136, the central, young, massive star cluster in the 30 Doradus complex of the Large Magellanic Cloud, are the single stars whose masses substantially exceed the canonical stellar upper mass limit of 150 M_sun. A recent study by us, viz., that of Banerjee, Kroupa & Oh (2012; Paper I), which involves realistic N-body computations of star clusters mimicking R136, indicates that such "super-canonical" (SC) stars can be formed out of a dense stellar population with a canonical initial mass function (IMF) through dynamically induced mergers of the most massive binaries. Here we study the formation of SC stars in the R136 models of Paper I in detail. To avoid forming extraneous SC stars from initially highly eccentric primordial binaries as in Paper I, we compute additional models with only initially circular primordial binaries. We also take into account the mass-evolution of the SC stars using detailed stellar evolutionary models that incorporate updated treatments of stellar winds. We find that SC stars begin to form via dynamical mergers of massive binaries from approx. 1 Myr cluster age. We obtain SC stars with initial masses up to approx. 250 M_sun from these computations. Multiple SC stars are found to remain bound to the cluster simultaneously within a SC-lifetime. These properties of the dynamically formed SC stars are consistent with those observed in R136. In fact, the stellar evolutionary models of SC stars imply that had they formed primordially along with the rest of the R136 cluster, i.e., violating the canonical upper limit, they would have evolved below the canonical 150 M_sun limit by approx. 3 Myr, the likely age of R136, and would not have been observable as SC stars at the present time in R136. This strongly supports the dynamical formation scenario of the observed SC stars in R136.

preprint2011arXiv

A new type of compact stellar population: dark star clusters

Among the most explored directions in the study of dense stellar systems is the investigation of the effects of the retention of supernova remnants, especially that of the massive stellar remnant black holes (BHs), in star clusters. By virtue of their eventual high central concentration, these stellar mass BHs potentially invoke a wide variety of physical phenomena, the most important ones being emission of gravitational waves (GWs), formation of X-ray binaries, and modification of the dynamical evolution of the cluster. Here we propose, for the first time, that rapid removal of stars from the outer parts of a cluster by the strong tidal field in the inner region of our Galaxy can unveil its BH sub-cluster, which appears as a star cluster that is gravitationally bound by an invisible mass. We study the formation and properties of such systems through direct N-body computations and estimate that they can be present in significant numbers in the inner region of the Milky Way. We call such objects "dark star clusters" (DSCs) as they appear dimmer than normal star clusters of similar mass and they comprise a predicted, new class of entities. The finding of DSCs will robustly cross-check BH retention; they will not only constrain the uncertain natal kicks of BHs, thereby the widely debated theoretical models of BH formation, but will also pinpoint star clusters as potential sites for GW emission for forthcoming ground-based detectors such as the Advanced LIGO. Finally, we also discuss the relevance of DSCs for the nature of IRS 13E.

preprint2011arXiv

Runaway massive stars from R136: VFTS 682 is very likely a "slow runaway"

We conduct a theoretical study on the ejection of runaway massive stars from R136 --- the central massive, star-burst cluster in the 30 Doradus complex of the Large Magellanic Cloud. Specifically, we investigate the possibility of the very massive star (VMS) VFTS 682 being a runaway member of R136. Recent observations of the above VMS, by virtue of its isolated location and its moderate peculiar motion, have raised the fundamental question whether isolated massive star formation is indeed possible. We perform the first realistic N-body computations of fully mass-segregated R136-type star clusters in which all the massive stars are in primordial binary systems. These calculations confirm that the dynamical ejection of a VMS from a R136-like cluster, with kinematic properties similar to those of VFTS 682, is common. Hence the conjecture of isolated massive star formation is unnecessary to account for this VMS. Our results are also quite consistent with the ejection of 30 Dor 016, another suspected runaway VMS from R136. We further note that during the clusters' evolution, mergers of massive binaries produce a few single stars per cluster with masses significantly exceeding the canonical upper-limit of 150 solar mass. The observations of such single super-canonical stars in R136, therefore, do not imply an IMF with an upper limit greatly exceeding the accepted canonical 150 solar mass limit, as has been suggested recently, and they are consistent with the canonical upper limit.

preprint2011arXiv

Stellar mass black holes in star clusters: gravitational wave emission and detection rates

We investigate the dynamics of stellar-mass black holes (BH) in star clusters focusing on the dynamical formation of BH-BH binaries, which are very important sources of gravitational waves (GW). We examine the properties of these BH-BH binaries through direct N-body computations of Plummer clusters, having initially N(0) <= 10^5 low mass stars and a population of stellar mass BHs, using the state-of-the-art N-body integrator "NBODY6". We find that the stellar mass BHs segregate rapidly into the cluster core to form a central dense sub-cluster of BHs in which BH-BH binaries form via 3-body encounters. While most of the BH binaries finally escape from the cluster by recoils due to super-elastic encounters with the single BHs, we find that for clusters with N(0) >= 5 X 10^4, typically a few of them dynamically harden to the extent that they can merge via GW emission within the cluster. Also, for each of such clusters, there are a few escaped BH binaries that merge within a Hubble time, most of the mergers happening within a few Gyr of cluster evolution. These results imply that the intermediate-aged massive clusters constitute the most important class of star clusters that can produce dynamical BH-BH mergers at the present epoch. The BH-BH merger rates obtained from our computations imply a significant detection rate (~ 30/yr) for the "Advanced LIGO" GW detector that will become operative in the near future. Finally, we briefly discuss our ongoing development on this work incorporating the formation of BHs in star clusters from stellar evolution. In particular, we highlight the effect of stellar metallicity on the BH sub-cluster driven expansion of a star cluster's core.

preprint2009arXiv

Stellar-mass black holes in star clusters: implications for gravitational wave radiation

We study the dynamics of stellar-mass black holes (BH) in star clusters with particular attention to the formation of BH-BH binaries, which are interesting as sources of gravitational waves (GW). We examine the properties of these BH-BH binaries through direct N-body simulations of star clusters using the GPU-enabled NBODY6 code. We perform simulations of N <= 10^5 Plummer clusters of low-mass stars with an initial population of BHs. Additionally, we do several calculations of star clusters confined within a reflective boundary mimicking only the core of a massive cluster. We find that stellar-mass BHs with masses ~ 10 solar mass segregate rapidly into the cluster core and form a sub-cluster of BHs within typically 0.2 - 0.5 pc radius, which is dense enough to form BH-BH binaries through 3-body encounters. While most BH binaries are ejected from the cluster by recoils received during super-elastic encounters with the single BHs, few of them harden sufficiently so that they can merge via GW emission within the cluster. We find that for clusters with $N \ga 5\times 10^4$, typically 1 - 2 BH-BH mergers occur within them during the first ~ 4 Gyr of evolution. Also for each of these clusters, there are a few escaping BH binaries that can merge within a Hubble time, most of the merger times being within a few Gyr. These results indicate that intermediate-age massive clusters constitute the most important class of candidates for producing dynamical BH-BH mergers. Old globular clusters cannot contribute significantly to the present-day BH-BH merger rate since most of the mergers from them would have occurred earlier. In contrast, young massive clusters are too young to produce significant number of BH-BH mergers. Our results imply significant BH-BH merger detection rates for the proposed "Advanced LIGO" GW detector. (Abridged)